SEARCH · Search NASA

Results for “seismic data processing”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 541 records · Page 30

All sources NASA NTRS DOE OSTI

DOE OSTI2019

EGS Collab Project: Status, Tests, and Data

The EGS (Enhanced Geothermal Systems) Collab project is performing stimulation and flow experiments in highly-monitored and well-characterized intermediate-scale (approximately10 to 20 meter) field test beds at a depth of approximately 1,500 meters in the Sanford Underground Research Facility (SURF) in the Black Hills of South Dakota. Our fracture stimulation and interwell flow tests are performed to better understand processes that control formation of effective subsurface heat exchangers that are critical to the development and success of EGS. Different EGS Collab stimulations will be performed under dissimilar stress conditions produce data for model comparisons that better differentiate stimulation mechanisms and the evolution of permeability enhancement in crystalline rock. EGS Collab experiments provide a means of testing tools, concepts and strategies that could later be employed under geothermal reservoir conditions at DOE’s Frontier Observatory for Research in Geothermal Energy (FORGE) and other enhanced geothermal systems. Key to the project is using numerical simulations in the experiment design and interpretation of results, and the collection of high-quality stimulation/flow, thermal, and geologic/geophysical data. A wide range of geologic and geophysical data have been acquired and interpreted for pre-stimulation testbed characterization. Multiple stimulation and flow tests have been performed at SURF, and multiple data sets have been collected and analyzed during these tests including electrical resistance tomography, microseismic, continuous active source seismic monitoring, distributed temperature, and the SIMFIP tool has been used to identify fracture formation and characteristics. An overview of some of the numerical simulations and experiments performed at SURF are presented.

15 GEOTHERMAL ENERGY↗

DOE OSTI2015

Preliminary results of local earthquake tomography around Bali, Lombok, and Sumbawa regions

Bali, Sumbawa, and Lombok regions are located in active tectonic influence by Indo-Australia plate subducts beneath Sunda plate in southern part and local back-arc thrust in northern part the region. Some active volcanoes also lie from eastern part of Java, Bali, Lombok and Sumbawa regions. Previous studies have conducted subsurface seismic velocity imaging using regional and global earthquake data around the region. In this study, we used P-arrival time from local earthquake networks compiled by MCGA, Indonesia within time periods of 2009 up to 2013 to determine seismic velocity structure and simultaneously hypocenter adjustment by applying seismic tomography inversion method. For the tomographic inversion procedure, we started from 1-D initial velocity structure. We evaluated the resolution of tomography inversion results through checkerboard test and calculating derivative weigh sum. The preliminary results of tomography inversion show fairly clearly high seismic velocity subducting Indo-Australian and low velocity anomaly around volcano regions. The relocated hypocenters seem to cluster around the local fault system such as back-arc thrust fault in northern part of the region and around local fault in Sumbawa regions. Our local earthquake tomography results demonstrated consistent with previous studies and improved the resolution. For future works, we will determine S-wave velocity structure using S-wave arrival time to enhance our understanding of geological processes and for much better interpretation.

58 GEOSCIENCES↗

DOE OSTI2018

Integrated pre-feasibility study for CO 2 geological storage in the Cascadia Basin, offshore Washington State, British Columbia

The Cascadia CarbonSAFE project conducted a pre-feasibility study to evaluate technical and nontechnical aspects of storing 50 million metric tons (MMT) of carbon dioxide in a safe, ocean basalt reservoir offshore Washington State and British Columbia. Sub-seafloor basalts are very common on Earth and enable geological mineralization as a long-term storage mechanism, permanently sequestering the carbon in solid rock form. This project evaluated the offshore storage complex, developed potential industrial source/transport scenarios, built an inventory of existing geophysical/geological data and environmental monitoring capabilities, conducted laboratory studies and reservoir modeling to determine storage capacity, and analyzed economic factors, regulatory requirements, and project management risks. Our team included researchers at Columbia University in New York, Pacific Northwest National Laboratory (PNNL) in Washington, University of Victoria in British Columbia, University of California at Santa Cruz, and GHG Underground in Maine, as well as collaborators at University of Iceland in Reykjavik and University of Washington in Saint Louis. In the study region, large emitters (>100,000 MT/year) generate a total of approximately 40 MMT annually from stationary sources, such as power plants, refineries, ammonia production, and mineral processing plants. Five potential industry source/transport scenarios were identified (three in U.SA S. and two in Canada) to potentially provide 50 MMT of carbon dioxide to the offshore reservoir. One scenario included net carbon-negative sources that would reduce atmospheric levels. The inventory of existing data illustrated that basalt properties in the region are potentially beneficial for long-term storage. Permeability is on the order of 0.1 to 1 Darcies within the uppermost 600 m of ocean crust and density-derived porosity reaches 10–20% in thin layers that provide act as flow channels and provide access to porous and permeable basalt. Fine-grained sediments overlying basalt provide a low-permeability barrier separating the storage reservoir from the overlying ocean. Subsurface monitoring of the storage complex benefits from a unique opportunity close proximity of to the nearby NEPTUNE cabled seafloor observatory, which actively transmits data to shore and could be used to monitor reservoir properties. Injection pressure changes, reservoir chemistry, seismicity, well head corrosion, and potential leakage monitoring would be feasible in real time. Using existing data and published models, the STOMP-CO 2 modeling code developed at PNNL was used to simulate carbon dioxide injection, storage, and mineralization in the reservoir. Preliminary simulations indicate that a 50 MMT carbon dioxide plume injected over a 20-year period remains well within the reservoir area, both laterally and vertically, for at least a 50-year period of post-injection monitoring. Laboratory studies on core samples from Cascadia measured dissolution rates in CO 2 –basalt–seawater mixtures and indicated large variability in rates with the potential for complete mineralization. Compared to published data, the dissolution rate measured in these samples is faster, decreasing the time until full conversion to carbonate minerals would occur. Adding geochemical reactivity in the STOMP-CO 2 model indicates that carbon dioxide would be fully converted to carbonate minerals in 135 years or less after injection ceases. The project assessed the viability and costs for the five scenarios using either dedicated offshore pipelines or shipping vessels. Carbon dioxide capture costs were not considered. Shipping incurred lower upfront CAPEX costs, and because most viable sources were located near the shore and waterways, shipping is preferred; no dedicated carbon dioxide pipelines exist in the region. The project also considered potential risks from mixed-use impacts on surrounding lands and ocean areas, such as active use by commercial and recreational fisherman, and other community considerations in planning for public engagement in the future. Other project risks were identified, including project management, permitting and compliance, and financing, as well as technical issues such as long-term instrumentation, plume migration, and model development. Geological concerns like reservoir permeability, seal, and induced seismicity held relatively low risk. The regulatory framework for offshore storage in the U.S. and Canada may be covered by general legal/regulatory programs, but is not yet fully defined. Existing programs appear to restrict the transportation of carbon dioxide across national boundaries (e.g., Canada to U.S.) and its sub-seabed disposal in some circumstances. Our assessment suggests that developing specific and separate pilot projects will help to advance the technical feasibility and regulatory requirements in each country, and lessons learned from such activities may be transferable elsewhere around the globe.

42 ENGINEERING↗

NASA NTRS1985

Heat and extension at mid- and lower crustal levels of the Rio Grande rift

The process by which large amounts (50 to 200 percent) of crustal extension are produced was concisely described by W. Hamilton in 1982 and 1983. More recently, England, Sawyer, P. Morgan and others have moved toward quantifying models of lithospheric thinning by incorporating laboratory and theoretical data on rock rheology as a function of composition, temperature, and strain rate. Hamilton's description identifies three main crustal layers, each with a distinctive mechanical behavior; brittle fracturing and rotation in the upper crust, discontinuous ductile flow in the middle crust and laminar ductile flow in the lower crust. The temperature and composition dependent brittle-ductile transition essentially defines the diffuse boundary between upper and middle crust. It was concluded that the heat responsible for the highly ductile nature of the lower crust and the lensoidal and magma body structures at mid-crustal depths in the rift was infused into the crust by relatively modest ( 10 percent by mass) magmatic upwelling (feeder dikes) from Moho levels. Seismic velocity-versus-depth data, supported by gravity modeling and the fact that volumes of rift related volcanics are relatively modest ( 6000 cubic km) for the Rio Grande system, all imply velocities and densities too small to be consistent with a massive, composite, mafic intrusion in the lower crust.

Olsen, K. H.↗

DOE OSTI2015

Poromechanics of stick-slip frictional sliding and strength recovery on tectonic faults

Pore fluids influence many aspects of tectonic faulting including frictional strength aseismic creep and effective stress during the seismic cycle. But, the role of pore fluid pressure during earthquake nucleation and dynamic rupture remains poorly understood. Here we report on the evolution of pore fluid pressure and porosity during laboratory stick-slip events as an analog for the seismic cycle. We sheared layers of simulated fault gouge consisting of glass beads in a double-direct shear configuration under true triaxial stresses using drained and undrained fluid conditions and effective normal stress of 5–10 MPa. Shear stress was applied via a constant displacement rate, which we varied in velocity step tests from 0.1 to 30 µm/s. Here, we observe net pore pressure increases, or compaction, during dynamic failure and pore pressure decreases, or dilation, during the interseismic period, depending on fluid boundary conditions. In some cases, a brief period of dilation is attendant with the onset of dynamic stick slip. Our data show that time-dependent strengthening and dynamic stress drop increase with effective normal stress and vary with fluid conditions. For undrained conditions, dilation and preseismic slip are directly related to pore fluid depressurization; they increase with effective normal stress and recurrence time. Microstructural observations confirm the role of water-activated contact growth and shear-driven elastoplastic processes at grain junctions. These results indicate that physicochemical processes acting at grain junctions together with fluid pressure changes dictate stick-slip stress drop and interseismic creep rates and thus play a key role in earthquake nucleation and rupture propagation.

58 GEOSCIENCES↗

DOE OSTI2022

Unsupervised learning from three-component accelerometer data to monitor the spatiotemporal evolution of meso-scale hydraulic fractures

Enhanced geothermal systems can provide a substantial share of the global energy demand. There exist several hurdles in the engineering implementations of such geothermal systems. One such hurdle is the accurate monitoring of the fracture networks created in subsurface through hydraulic stimulation of these systems. Micro seismicity associated with the stimulation is the primary means to locate the event hypocenters for estimating the stimulated rock volume. Existing methods for location the hypocenters are restricted to only the highest amplitude impulsive signals that are simultaneously detected on several sensors. Consequently, a large portion (usually ~99%) of the measurements are left unused. In this paper, an unsupervised manifold-approximation followed by clustering of 3-component accelerometer data is used to analyze the seismicity recorded on a monitoring well. With this method, a larger portion of the measured signal is used for the monitoring of the hydraulic fracture network. We analyze the EGS Collab experiment 1 microseismic data, recorded at the Sanford Underground Research Facility, South Dakota. Using the data from a single three-component accelerometer, the polarization features viz. Azimuth, incidence, rectilinearity, and planarity are used as inputs for the unsupervised manifold approximation followed by clustering. Our study shows that density-based clusters in the projected 3D space correspond to distinct types of hydraulically fractured zones around the injection point. Finally, we show that the temporal evolution of these clusters can be used to track fracture creation and propagation.

58 GEOSCIENCES↗

DOE OSTI2022

Quantitative interpretation of time-lapse seismic data at Farnsworth field unit: Rock physics modeling, and calibration of simulated time-lapse velocity responses

Here, this study investigates the contribution of fluid saturation variation to the time-lapse velocity response by performing fluid substitution modeling. The methodology is exemplified by the time-lapse seismic monitoring of carbon dioxide at Farnsworth field unit (FWU). In order to evaluate the fluid distribution in a matured oil reservoir, the Southwest Regional Partnership (SWP) acquired multiple vertical seismic profile (VSP) surveys at different times during the CO 2 –water alternatinggas (WAG) injection period. In this work, we present a thorough methodology for computing the elastic response of the saturated rock for different fluid saturations using a site-specific petro-elastic model (PEM). The output from the PEM was combined with results from a fluid compositional model to compute the seismic velocities at times corresponding to each VSP survey. To produce a calibrated simulated response, the measured time-lapse seismic velocities were integrated into the numerical simulation model. The mismatches between the predicted and measured time-lapse velocities were minimized through an iterative calibration process using a trained artificial neural network proxy (ANN) coupled with a particle swarm optimizer (PSO). Our study indicates that the hybrid optimization workflow can effectively perform the history matching. With an accurate prediction of the hydrodynamic properties, the migration of CO 2 within the subsurface was modeled by predicting the spatial velocity distribution for a radius of 305 m around the injection well. The technology demonstrated and the expertise gained from this study can guide similar CO 2 -WAG projects.

58 GEOSCIENCES↗

DOE OSTI2020

The Local Seismoacoustic Wavefield of a Research Nuclear Reactor and Its Response to Reactor Power Level

In this article, we describe the seismoacoustic wavefield recorded outdoors but inside the facility fence of the High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory (Tennessee). HFIR is a research nuclear reactor that generates neutrons for scattering, irradiation research, and isotope production. This reactor operates at a nominal power of 85 MW, with a full-power period between 24 and 26 days. This study uses data from a single seismoacoustic station that operated for 60 days and sampled a full operating reactor cycle, that is, full-power operation and end-of-cycle outage. The analysis presented here is based on identifying signals that characterize the steady, that is, full-power operation and end-of-cycle outage, and transitional, that is, start-up and shutdown, states of the reactor. We found that the overall seismoacoustic energy closely follows the main power cycle of the reactor and identified spectral regions excited by specific reactor operational conditions. In particular, we identified a tonal noise sequence with a fundamental frequency around 21.4 Hz and multiple harmonics that emerge as the reactor reaches 90% of nominal power in both seismic and acoustic channels. We also utilized temperature measurements from the monitoring system of the reactor to suggest links between the operation of reactor’s subsystems and seismoacoustic signals. We demonstrate that seismoacoustic monitoring of an industrial facility can identify and track some industrial processes and detect events related to operations that involve energy transport.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2016

Frequency domain analysis of errors in cross-correlations of ambient seismic noise

We analyse random errors (variances) in cross-correlations of ambient seismic noise in the frequency domain, which differ from previous time domain methods. Extending previous theoretical results on ensemble averaged cross-spectrum, we estimate confidence interval of stacked cross-spectrum of finite amount of data at each frequency using non-overlapping windows with fixed length. The extended theory also connects amplitude and phase variances with the variance of each complex spectrum value. Analysis of synthetic stationary ambient noise is used to estimate the confidence interval of stacked cross-spectrum obtained with different length of noise data corresponding to different number of evenly spaced windows of the same duration. This method allows estimating Signal/Noise Ratio (SNR) of noise cross-correlation in the frequency domain, without specifying filter bandwidth or signal/noise windows that are needed for time domain SNR estimations. Based on synthetic ambient noise data, we also compare the probability distributions, causal part amplitude and SNR of stacked cross-spectrum function using one-bit normalization or pre-whitening with those obtained without these pre-processing steps. Natural continuous noise records contain both ambient noise and small earthquakes that are inseparable from the noise with the existing pre-processing steps. Using probability distributions of random cross-spectrum values based on the theoretical results provides an effective way to exclude such small earthquakes, and additional data segments (outliers) contaminated by signals of different statistics (e.g. rain, cultural noise), from continuous noise waveforms. This technique is applied to constrain values and uncertainties of amplitude and phase velocity of stacked noise cross-spectrum at different frequencies, using data from southern California at both regional scale (~35 km) and dense linear array (~20 m) across the plate-boundary faults. A block bootstrap resampling method is used to account for temporal correlation of noise cross-spectrum at low frequencies (0.05–0.2 Hz) near the ocean microseismic peaks.

58 GEOSCIENCES↗

DOE OSTI2018

Central and Eastern North America Ground-Motion Characterization (NGA-East)

This document is the final project report of the Next Generation Attenuation for Central and Eastern North America (CENA) project (NGA-East). The NGA-East objective was to develop a new ground-motion characterization (GMC) model for the CENA region. The GMC model consists of a set of new ground-motion models (GMMs) for median and standard deviation of ground motions and their associated weights to be used with logic-trees in probabilistic seismic hazard analyses (PSHA). NGA-East is a large multidisciplinary project coordinated by the Pacific Earthquake Engineering Research Center (PEER), at the University of California. The project has two components: (1) a set of scientific research tasks, and (2) a model-building component following the framework of the “Seismic Senior Hazard Analysis Committee (SSHAC) Level 3” (Budnitz et al. 1997; NRC 2012). Component (2) is built on the scientific results of component (1) of the NGA-East project. This report documents the tasks under component (2) of the project. Under component (1) of NGA-East, several scientific issues were addressed, including: (a) development of a new database of ground motion data recorded in CENA; (b) development of a regionalized ground-motion map for CENA, (c) definition of the reference site condition; (d) simulations of ground motions based on different methodologies; and (e) development of numerous GMMs for CENA. The scientific tasks of NGA-East were all documented as a series of PEER reports. The scope of component (2) of NGA-East was to develop the complete GMC. This component was designed as a SSHAC Level 3 study with the goal of capturing the ground motions’ center, body, and range of the technically defensible interpretations in light of the available data and models. The SSHAC process involves four key tasks: evaluation, integration, formal review by the Participatory Peer Review Panel (PPRP), and documentation (this report). Key tasks documented in this report include review and evaluation of the empirical groundmotion database, the regionalization of ground motions, and screening sets of candidate GMMs. These are followed by the development of new median and standard deviation GMMs, the development of new analyses tools for quantifying the epistemic uncertainty in ground motions, and the documentation of implementation guidelines of the complete GMC for PSHA computations. Appendices include further documentation of the relevant SSHAC process and additional supporting technical documentation of numerous sensitivity analyses results. The PEER reports documenting component (1) of NGA-East are also considered “attachments” to the current report and are all available online on the PEER website (https://peer.berkeley.edu/). The final NGA-East GMC model includes a set of 17 GMMs defined for 24 ground-motion intensity measures, applicable to CENA in the moment magnitude range of 4.0 to 8.2 and covering distances up to 1500 km. Standard deviation models are also provided for site-specific analysis (single-station standard deviation) and for general PSHA applications (ergodic standard deviation). Adjustment factors are provided for consideration of source-depth effects and hanging-wall effects, as well as for hazard computations at sites in the Gulf Coast region.

58 GEOSCIENCES↗

NASA NTRS1997

Structure and Kinematics of a Complex Crater: Upheaval Dome, Southeast, Utah

Two vastly different phenomena, extraterrestrial impact and salt diapirism, have been proposed for the origin of Upheaval Dome. Upheaval Dome is a about 2.5-km-diameter structural dome surrounded by a 5-km-diameter ring structural depression, which is in turn flanked by extensive, nearly flat-lying Colorado Plateau strata. Seismic refraction data and geologic mapping indicate that the dome originated by the collapse of a transient cavity formed by impact; data also show that rising salt has had a negligible influence on dome development. Evidence for this includes several factors: (1) a rare lag deposit of impactite is present; (2) fan-tailed fracture surfaces (shatter surfaces) and a few shattercones are present; (3) the top of the underlying salt horizon is at least 500 m below the center of the dome, with no exposures of salt in the dome to support the possibility that a salt diapir has ascended through it; (4) sedimentary strata in the center are significantly imbricated by top-to-the-center thrust faulting and are complexly folded; (5) top-to-the-center low-angle normal faults are found at the perimeter of the structure; and (6) clastic dikes are widespread. The scarcity of melt rocks and shock fabrics is attributed to approximately 0.5 km of erosion; the structures of the dome reflect processes of complex crater development at a depth of about 0.5 km below the crater floor. Based on mapping and kinematic analysis, we infer that the dome formed mainly by centerward motion of rock units along listric faults. Outcrop-scale folding and upturning of beds, especially common in the center, largely resulted from this motion. In addition, we have detected some centerward motion of fault-bounded wedges resulting from displacements on subhorizontal faults that conjoin and die out within horizontal bedding in the perimeter of the structure. Collectively, the observed deformation accounts for the creation of both the central uplift and the encircling ring syncline.

Kriens, B. J.↗

DOE OSTI2019

Accelerating full waveform inversion via source stacking and cross-correlations

SUMMARY Accurate synthetic seismic wavefields can now be computed in 3-D earth models using the spectral element method (SEM), which helps improve resolution in full waveform global tomography. However, computational costs are still a challenge. These costs can be reduced by implementing a source stacking method, in which multiple earthquake sources are simultaneously triggered in only one teleseismic SEM simulation. One drawback of this approach is the perceived loss of resolution at depth, in particular because high-amplitude fundamental mode surface waves dominate the summed waveforms, without the possibility of windowing and weighting as in conventional waveform tomography. This can be addressed by redefining the cost-function and computing the cross-correlation wavefield between pairs of stations before each inversion iteration. While the Green’s function between the two stations is not reconstructed as well as in the case of ambient noise tomography, where sources are distributed more uniformly around the globe, this is not a drawback, since the same processing is applied to the 3-D synthetics and to the data, and the source parameters are known to a good approximation. By doing so, we can separate time windows with large energy arrivals corresponding to fundamental mode surface waves. This opens the possibility of designing a weighting scheme to bring out the contribution of overtones and body waves. It also makes it possible to balance the contributions of frequently sampled paths versus rarely sampled ones, as in more conventional tomography. Here we present the results of proof of concept testing of such an approach for a synthetic 3-component long period waveform data set (periods longer than 60 s), computed for 273 globally distributed events in a simple toy 3-D radially anisotropic upper mantle model which contains shear wave anomalies at different scales. We compare the results of inversion of 10 000 s long stacked time-series, starting from a 1-D model, using source stacked waveforms and station-pair cross-correlations of these stacked waveforms in the definition of the cost function. We compute the gradient and the Hessian using normal mode perturbation theory, which avoids the problem of cross-talk encountered when forming the gradient using an adjoint approach. We perform inversions with and without realistic noise added and show that the model can be recovered equally well using one or the other cost function. The proposed approach is computationally very efficient. While application to more realistic synthetic data sets is beyond the scope of this paper, as well as to real data, since that requires additional steps to account for such issues as missing data, we illustrate how this methodology can help inform first order questions such as model resolution in the presence of noise, and trade-offs between different physical parameters (anisotropy, attenuation, crustal structure, etc.) that would be computationally very costly to address adequately, when using conventional full waveform tomography based on single-event wavefield computations.

Romanowicz, Barbara↗

NASA NTRS2012

Seismic Imager Space Telescope

A concept has been developed for a geostationary seismic imager (GSI), a space telescope in geostationary orbit above the Pacific coast of the Americas that would provide movies of many large earthquakes occurring in the area from Southern Chile to Southern Alaska. The GSI movies would cover a field of view as long as 300 km, at a spatial resolution of 3 to 15 m and a temporal resolution of 1 to 2 Hz, which is sufficient for accurate measurement of surface displacements and photometric changes induced by seismic waves. Computer processing of the movie images would exploit these dynamic changes to accurately measure the rapidly evolving surface waves and surface ruptures as they happen. These measurements would provide key information to advance the understanding of the mechanisms governing earthquake ruptures, and the propagation and arrest of damaging seismic waves. GSI operational strategy is to react to earthquakes detected by ground seismometers, slewing the satellite to point at the epicenters of earthquakes above a certain magnitude. Some of these earthquakes will be foreshocks of larger earthquakes; these will be observed, as the spacecraft would have been pointed in the right direction. This strategy was tested against the historical record for the Pacific coast of the Americas, from 1973 until the present. Based on the seismicity recorded during this time period, a GSI mission with a lifetime of 10 years could have been in position to observe at least 13 (22 on average) earthquakes of magnitude larger than 6, and at least one (2 on average) earthquake of magnitude larger than 7. A GSI would provide data unprecedented in its extent and temporal and spatial resolution. It would provide this data for some of the world's most seismically active regions, and do so better and at a lower cost than could be done with ground-based instrumentation. A GSI would revolutionize the understanding of earthquake dynamics, perhaps leading ultimately to effective warning capabilities, to improved management of earthquake risk, and to improved public safety policies. The position of the spacecraft, its high optical quality, large field of view, and large field of regard will make it an ideal platform for other scientific studies. The same data could be simply reused for other studies. If different data, such as multi-spectral data, is required, additional instruments could share the telescope.

Sidick, Erkin↗

DOE OSTI2019

EGS Collab Earth Modeling: Integrated 3D Model of the Testbed

The EGS Collab project is conducting a series of stimulation and interwell flow tests in an intermediate scale (~10-20 m) testbed located on the 4850 level in the Sanford Underground Research Facility (SURF) in Lead, South Dakota. The testbed consists of eight ~200 ft (~60 m) HQ-diameter (9.6 cm) boreholes that are drilled into the Poorman Formation from the western rib of the West Access Drift. Of the eight boreholes, one borehole is used as an injection/stimulation well, another sub-parallel borehole located about 10 m away from the injection well is used as a production well, and rest of the other boreholes are used as monitoring wells. For 3D visualization of the testbed as well as to provide spatially accurate and consistent parameter data for geomechanical and geophysical process modeling tasks, we are using Leapfrog Software (Seequent Limited). For analysis of fractures and creation of conditioned stochastic discrete fracture network (CS-DFN) model of the testbed, we are using FracMan (Golder Associates Inc.). As of mid-May 2019, the EGS Collab Leapfrog database contains the general layout of the drifts/shafts at the 4850 and 4100 levels; generalized geologic framework model of a domain that spans both the 4850 and 4100 levels; the CS-DFN of the testbed; and layouts of the EGS Collab testbed boreholes, kISMET testbed boreholes, and other existing boreholes at the 4850 level. This information has been obtained from the SURF Vulcan database, detailed laser scan and point surveys of borehole collar locations, directional surveys of borehole trajectories, and mapping of features in the drifts, boreholes, and cores. Similarly, other types of data that have been incorporated into the Leapfrog database include static features of the testbed such as natural discontinuities (e.g., fractures, weeps), location of monitoring instruments, and stimulation notches. Baseline geophysical characterization results (e.g., seismic, ERT, temperature, etc.) are also imported into the Leapfrog database. The location and size of packers, various point-source signals detected during stimulation and flow tests (e.g., temperature anomalies, micro-earthquake events, etc.) are also compiled in the database. All feature locations are reported using the local Homestake Mine coordinate system. In this paper, we present some of the geostructural characteristics as well as the status of the testbed as represented by Leapfrog/FracMan database/visualizations.

15 GEOTHERMAL ENERGY↗

DOE OSTI2023

The Thermal Plumbing System of Stromboli Volcano, Aeolian Islands (Italy) Inferred From Electrical Conductivity and Induced Polarization Tomography

Abstract We performed the first 3D island‐scale tomography of the electrical conductivity of Stromboli volcano (Aeolian Islands, Italy) using 2D acquisition lines (37.2 km) and a total of 18,880 measurements and 2,402 unique electrode locations. This 3D data set was inverted using a Gauss‐Newton algorithm, parallel‐processing on an unstructured tetrahedral mesh containing 678,420 finite‐element nodes and 3,580,145 elements to account for the topography of the volcanic island. The tomogram exhibits a conductive body (10 −2 –1.0 S m −1 ) consistent with the location of CO 2 and temperature anomalies observed at the ground surface. It corresponds to the hydrothermal system with high electrical conductivity associated with alteration. In order to confirm this interpretation, a 2.5D large‐scale induced polarization tomography was performed crossing the volcano. The joint interpretation of the conductivity and normalized chargeability is done with a petrophysical model previously tested and verified at both shield‐ and strato‐volcanoes. This model implies that alteration (through the effect of the cation exchange capacity associated with clay minerals and zeolites) plays a strong role in both controlling the electrical conductivity and normalized chargeability at Stromboli volcano. A temperature tomogram, derived from the geoelectrical measurements, is consistent with surface temperature anomalies and the Very Long Period (VLP) seismicity related to the mild‐explosive activity. This survey displays at 600 m a.s.l. a lateral shift in the highest temperature location, also corresponding to the source of VLP seismicity. Structural boundaries have a major role in the hottest hydrothermal fluids rising below the active crater terrace of Stromboli volcano.

58 GEOSCIENCES↗

DOE OSTI2019

Using Machine Learning to Discern Eruption in Noisy Environments: A Case Study Using CO2-Driven Cold-Water Geyser in Chimayó, New Mexico

In this paper, we present an approach based on machine learning (ML) to distinguish eruption and precursory signals of Chimayó geyser (New Mexico, U.S.A.) under noisy environmental conditions. This geyser can be considered a natural analog of CO 2 intrusion into shallow water aquifers. By studying this geyser, we can understand upwelling of CO 2 -rich fluids from depth, which has relevance to leak monitoring in a CO 2 sequestration project. ML methods such as random forests (RFs) are known to be robust multiclass classifiers and perform well under unfavorable, noisy conditions. However, the extent of the RF method’s accuracy is poorly understood for this CO 2 -driven geysering application. The current study aims to quantify the performance of RF classifiers to discern the geyser state. Toward this goal, we first present the data collected from the seismometer that is installed near the Chimayó geyser. The seismic signals collected at this site contain different types of noises such as daily temperature variations, animal movement near the geyser, and human activity. First, we filter the signals from these noises by combining the Butterworth high-pass (BH) filter and an autoregressive (AR) method in a multilevel fashion. We show that by combining these filtering techniques in a hierarchical fashion leads to a reduction in noise in the seismic data without removing the precursors and eruption event signals. We then use RF on the filtered data to classify the state of geyser into three classes: remnant noise, precursor, and eruption states. RF classifier is constructed based on the comprehensive features extracted using the Tsfresh Python package. We show that the classification accuracy using RF on the filtered data is greater than 90%. We also evaluate the accuracy of other classical time-series methods such as dynamic time warping (DTW) on filtered data along with RF on partially filtered data in which we remove the daily trends. Classification accuracy shows that DTW performs poorly (44%) and RF on partially filtered data performs decently (87%). Denoising seismic signals from both daily trends and human activity enhances RF classifier performance by 7%. Finally, these aspects make the proposed ML framework attractive for event discrimination and signal enhancement under noisy conditions, with strong potential for application to monitor leaks in CO 2 sequestration.

58 GEOSCIENCES↗

DOE OSTI2018

Integrated characterization of CO 2 storage reservoirs on the Rock Springs Uplift combining geomechanics, geochemistry, and flow modeling

This project improves understanding of the effects CO 2 injection and storage on geomechanical, petrophysical, and other properties on two potential reservoirs, the Madison Limestone and Weber Sandstone, at the Rock Springs Uplift, southwest Wyoming. Rock physics studies, laboratory experiments, and computer simulations were employed. Two groups of rock samples were prepared, one aged with formation brine and the second aged with CO 2 -saturated formation brine and subsequently used for geomechanical, geochemical, and petrophysical testing. Rock physics studies provided a facies classification to discriminate among different lithologies and to investigate the feasibility of predicting petrophysical properties using seismic inversion. By applying a zone-wise regression relation established in this facies analysis, a 3D reservoir model was developed. It was determined that time-lapse seismic data are not feasible for monitoring CO 2 displacement at the Rock Springs Uplift or to mitigate leakage-risks situations because pressure and fluid saturation have small effects on elastic velocities and synthetic seismograms. The small variations that do occur are due to the low porosity of the sandstone and the high stiffness of the dolomite. In the Madison Limestone, mineral dissolution due to reaction with CO 2 -saturated brine was observed in highly permeable pathways while mineral precipitation happened in intergranular or micropores. Changes in pore size were quantified via Time-Domain Nuclear Magnetic Resonance transverse relaxation time (T 2 ) and diffusion coefficient distributions. As opposed to continuous injection of CO 2 -saturated brine, a subtle increase in porosity after static ‘soaking’ with CO 2 -saturated brine increases permeability due to the heterogeneity of changes in pore and pore throat sizes. Porosity-permeability relationships are fitted using a power-law function that can be adopted to model time-evolved porosity and permeability in a transport-limited CO 2 injection process. The effect of CO 2 -saturated brine on the linear elastic properties of Weber Sandstone and Madison Limestone cannot be generally concluded. Considering horizontal and vertical orientations and different confining pressures, no consistent relationship between stress (or strain) data and reaction with CO 2 -saturated brine was observed in the nonlinear plastic and post-failure regime. The change of elastic constants due to reaction with CO 2 -saturated brine is more significant in Madison Limestone than in Weber Sandstone; however, no consistent trend was observed. Reservoir simulations suggest that injectivity for Weber Sandstone is adversely impacted by low permeability. Changes in confining pressure do not impact injection rate and can be ignored while the effect of changes to pore pressure cannot. The effect of confining pressure and pore pressure on the size of the CO 2 plume in the Madison Limestone is negligible because of higher permeability. Geomechanical and petrophysical properties induced by CO 2 -brine-rock reactions were time-dependent. Simulations considered these time-dependent properties and indicate that injectivity, porosity distribution, pressure distribution and CO 2 saturation were affected.

58 GEOSCIENCES↗

NASA NTRSDate not supplied

Characterizing Martian Volcanic Provinces’ Magmatic Evolution and Chemistry through Equations of State Modeling Initial Study

Here we discuss a novel interdisciplinary approach to investigating igneous compositions of large volcanic provinces on Mars using remote sensing data sets (Mars Odyssey Gamma Ray and neutron Spectrometer suite--GRS, and gravity) to inform petrologic and thermoelastic modeling. Martian volcanic provinces, starting from Noachian to Amazonian age, including Elysium (EVP), are locations of great geologic interest which have been active over long time scales from Hesperian to Amazonian. Regional scale change in eruptive processes are poorly understood. Compared to large igneous provinces on Earth, the martian volcanic activity has persisted for orders of magnitude longer. Therefore, changes in mantle chemistry, pressure, and temperature during lithospheric cooling are expected to produce significant changes in the conditions of magma production, storage and ascent, affecting the degree of fractional crystallization and crustal contamination. In this study, we perform a detailed modeling to test the hypothesis that compositional variability within volcanic provinces resulted from spatiotemporal changes in the depth of magma formation and present initial results. Our methods include constraining the pressure and temperature conditions of EVP as a case study of geologically recent magmatic evolution on Mars using GRS informed surface chemistry constraints and pMELTS modeling. Second, we place constraints on the density and seismic velocities of the EVP melt through thermoelastic modeling with a local gravity analysis. This analysis is also extended to Noachian aged volcanic provinces on Mars initial results estimates mantle pressure for each sub-region is 16 kbar pressure, whereas, degree of partial melting is low and varies between 10 to 12. This study aims to validate our theoretical model and develop a perspective view of spatiotemporal changes at the interior of Mars throughout the time.

Martian science, Mars, geology, geophysics↗