Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “seismic data”

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.

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

At least 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Quantifying modeling uncertainty in simplified beam models for building response prediction

The use of simple models for response prediction of building structures is preferred in earthquake engineering for risk evaluations at regional scales, as they make computational studies more feasible. The primary impediment in their gainful use presently is the lack of viable methods for quantifying (and reducing upon) the modeling errors/uncertainties they bear. This study presents a Bayesian calibration method wherein the modeling error is embedded into the parameters of the model. Here, the method is specifically described for coupled shear-flexural beam models here, but it can be applied to any parametric surrogate model. The major benefit the method offers is the ability to consider the modeling uncertainty in the forward prediction of any degree-of-freedom or composite response regardless of the data used in calibration. The method is extensively verified using two synthetic examples. In the first example, the beam model is calibrated to represent a similar beam model but with enforced modeling errors. In the second example, the beam model is used to represent the detailed finite element model of a 52-story building. Both examples show the capability of the proposed solution to provide realistic uncertainty estimation around the mean prediction.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

The future of subsurface monitoring: AEC’s breakthroughs in CCS technology

Carbon capture and storage (CCS) has emerged as a key solution in the fight against climate change. However, for CCS to succeed, it is crucial to ensure that the sequestered CO2 stays safely trapped underground. The U.S. Department of Energy (DOE) has emphasized the need for advancements in subsurface monitoring, measurement, reporting, and verification. Aside from caprock integrity failure, the other primary failure points usually involve defective cement in the casing annulus of wellbores or plugged and abandoned wells. In addition, many energy producers (e.g., oil and gas, geothermal) and storage and disposal operators (e.g., H2 and water) must deal with the same issue. Poorly placed or degraded cement can create pathways for gas or fluid to escape from casing annuli and in plugged and abandoned or orphan wells, posing environmental risks. Yet, a reliable and cost-effective way to monitor cement and well integrity over multiple decades is still unavailable. Traditional geophysical methods like 4D seismic imaging and surface-based electromagnetic monitoring lack the resolution and accuracy for detecting these types of failures (Vasco et al., 2022; Fawad and Mondol, 2021). Wireline logging is expensive to run continuously and is obtrusive to the operation. While fiber optics can potentially be a solution, its bulkiness can significantly compromise the cement's integrity. To address these challenges, the Advanced Energy Consortium (AEC) at The University of Texas at Austin’s Bureau of Economic Geology (the Bureau) has been pioneering research in subsurface monitoring using its portfolio of distributed autonomous microfabricated sensors for harsh subsurface environments since 2008. A class of these microsensors [System on a Chip (SoC)] can be mixed in cement and permanently placed without compromising the cement column; the sensors would then communicate with each other or a data acquisition (DAQ) master node. Another class of the AEC microsensors can be fully autonomous, with rechargeable micro-batteries capable of exceeding 100°C, flash memory, and, currently, a pressure and temperature sensor. They are designed to circulate in mud, geothermal fluids, U-loops, or pipelines. They can log data into memory and are unobtrusive to operations. Our team has been working on a multi-year DOE-funded project (DE-FE0031856)—supported by $2.95M in federal funding and $0.75M in cost-matching from the AEC—to demonstrate SoC sensor utility for CO2 leakage monitoring in CCS applications. This multi-institutional collaboration developed a novel sensing architecture utilizing radiofrequency (RF) microsensors embedded within the cement sheath. These sensors detect CO2 migration and are interrogated via a Smart Casing Collar (SCC).

58 GEOSCIENCES↗
DOE OSTI2020

Evidence of shock-compressed stishovite above 300 GPa

SiO 2 is one of the most fundamental constituents in planetary bodies, being an essential building block of major mineral phases in the crust and mantle of terrestrial planets (1–10 M E ). Silica at depths greater than 300 km may be present in the form of the rutile-type, high pressure polymorph stishovite (P4 2 /mnm) and its thermodynamic stability is of great interest for understanding the seismic and dynamic structure of planetary interiors. Previous studies on stishovite via static and dynamic (shock) compression techniques are contradictory and the observed differences in the lattice-level response is still not clearly understood. Here, laser-induced shock compression experiments at the LCLS- and SACLA XFEL light-sources elucidate the high-pressure behavior of stishovite on the lattice-level under in situ conditions on the Hugoniot to pressures above 300 GPa. We find stishovite is still (meta-)stable at these conditions, and does not undergo any phase transitions. This contradicts static experiments showing structural transformations to the CaCl 2 , α-PbO 2 and pyrite-type structures. However, rate-limited kinetic hindrance may explain our observations. These results are important to our understanding into the validity of EOS data from nanosecond experiments for geophysical applications.

58 GEOSCIENCES↗
DOE OSTI2026

Electromagnetic Analysis of ITER Electron Cyclotron Emission Model

The ITER electron cyclotron emission (ECE) diagnostic system is located at diagnostics shielding module 2 (DSM2), equatorial port 9 (EP9), to measure electron temperature profile and electron temperature fluctuations, and also assess nonthermal electron distributions via the oblique view. Therefore, ECE has both radial and oblique views with two hot sources for calibration and two couples of mirrors for two different optical views. It shall receive that the large electromagnetic (EM) loads up to 100-MN/m3 force density due to the eddy currents on the copper mirrors generated by the short transient plasma disruption. The simplified EM analysis model using magnetic field ( B ) data and flux variations (dB/dt) method based on the worst case of plasma disruption, MD_DW_EXP16MS_CATIII, has been developed to calculate the force and moments for each ECE component and support structure for EP9 DSM2 Bay2 and Bay3. The B and dB/dt methods use the constant B and dB/dt assumed for each local bay space, but the real B and dB/dt are basically varied along the radial direction of each bay. The complex global model with the local ECE components and support structure using Maxwell transient has been developed to do EM analysis to compare the results of the B and dB/dt models. The volumetric EM force density of the whole ECE components and support structure can be used for the subsequent structural analysis. In conclusion, combined with the thermal and nuclear loads and seismic and initial loads, the ECE integration analysis can be finalized for the operation case.

Fang, J. [Princeton Plasma Physics Laboratory (PPP↗
DOE OSTI2020

Applying Waveform Correlation to Mining Blasts Using a Global Sparse Network

Agencies that monitor for underground nuclear tests are interested in techniques that automatically characterize mining blasts to reduce the human analyst effort required to produce high-quality event bulletins. Waveform correlation is effective in finding similar waveforms from repeating seismic events, including mining blasts. We report the results of an experiment that uses waveform templates recorded by multiple International Monitoring System stations of the Comprehensive Nuclear-Test-Ban Treaty for up to 10 years prior to detect and identify mining blasts that occur during single weeks of study. We discuss approaches for template selection, threshold setting, and event detection that are specialized for mining blasts and a sparse, global network. We apply the approaches to two different weeks of study for each of two geographic regions, Wyoming and Scandinavia, to evaluate the potential for establishing a set of standards for waveform correlation processing of mining blasts that can be effective for operational monitoring systems with a sparse network. We compare candidate events detected with our processing methods to the Reviewed Event Bulletin of the International Data Centre to develop an intuition about potential reduction in analyst workload.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2024

High-Pressure Melting Experiments of Fe 3 C and a Thermodynamic Model of Fe-C Liquids for the Earth's Core

Melting experiments of Fe 3 C were conducted to 85 GPa in laser-heated diamond anvil cells with in situ X-ray diffraction and post-experiment textural observation. From the determined pressure-temperature conditions of the melting curve for Fe 3 C, together with literature data on the melting point of diamond and eutectic point of the system Fe-Fe 3 C/Fe 7 C 3 under high pressures, we established a self-consistent thermodynamic model for high-pressure melting of the system Fe-C including the mixing parameters for liquids. The results show that mixing of Fe and C liquids is negatively nonideal from 1 bar to the pressure at the center of the Earth. The departure from ideal mixing becomes progressively larger with increasing pressure, which leads to greatly stabilized liquids under core pressures. The modeled carbon content in eutectic melts under core pressures is 3.3–4.4 wt%. From the Gibbs free energy, we derived an internally consistent parameters for Fe-C outer cores which included the crystallizing points at their bottoms, isentropic thermal profiles, and densities and longitudinal seismic wave speeds (Vp). While the addition of carbon in excess of the eutectic melt composition effectively reduces the density of iron liquid, the Vp of iron liquid is not greatly changed. Therefore, the low density and high Vp of PREM relative to pure iron cannot be reconciled by an Fe-C liquid. Therefore, the Earth's core cannot be approximated by the system Fe-C and should include another light element.

58 GEOSCIENCES↗
DOE OSTI2020

Effects of composition and pressure on electronic states of iron in bridgmanite

Electronic states of iron in the lower mantle's dominant mineral, (Mg,Fe,Al)(Fe,Al,Si)O 3 bridgmanite, control physical properties of the mantle including density, elasticity, and electrical and thermal conductivity. However, the determination of electronic states of iron has been controversial, in part due to different interpretations of Mössbauer spectroscopy results used to identify spin state, valence state, and site occupancy of iron. We applied energy-domain Mössbauer spectroscopy to a set of four bridgmanite samples spanning a wide range of compositions: 10–50% Fe/total cations, 0–25% Al/total cations, 12–100% Fe 3+ /total Fe. Measurements performed in the diamond-anvil cell at pressures up to 76 GPa below and above the high to low spin transition in Fe 3+ provide a Mössbauer reference library for bridgmanite and demonstrate the effects of pressure and composition on electronic states of iron. Results indicate that although the spin transition in Fe 3+ in the bridgmanite B-site occurs as predicted, it does not strongly affect the observed quadrupole splitting of 1.4 mm/s, and only decreases center shift for this site to 0 mm/s at ~70 GPa. Thus center shift can easily distinguish Fe 3+ from Fe 2+ at high pressure, which exhibits two distinct Mössbauer sites with center shift ~1 mm/s and quadrupole splitting 2.4–3.1 and 3.9 mm/s at ~70 GPa. Correct quantification of Fe 3+ /total Fe in bridgmanite is required to constrain the effects of composition and redox states in experimental measurements of seismic properties of bridgmanite. Furthermore, in Fe-rich, mixed-valence bridgmanite at deep-mantle-relevant pressures, up to ~20% of the Fe may be a Fe 2.5+ charge transfer component, which should enhance electrical and thermal conductivity in Fe-rich heterogeneities at the base of Earth's mantle.

58 GEOSCIENCES↗
DOE OSTI2021

Stratigraphically Controlled Stress Variations at the Hydraulic Fracture Test Site-1 in the Midland Basin, TX

We investigated the relationship between stratigraphy, stress, and microseismicity at the Hydraulic Fracture Test Site-1. The site comprises two sets of horizontal wells in the Wolfcamp shale and a deviated well drilled after hydraulic fracturing. Regional stresses indicate normal/strike-slip faulting with E-W compression. Stress measurements in vertical and horizontal wells show that the minimum principal stress varies with depth. Strata with high clay and organic content show high values of the least compressive stress, consistent with the theory of viscous stress relaxation. By integrating data from core, logs, and the hydraulic fracturing stages, we constructed a stress profile for the Wolfcamp sequence, which predicts how much pressure is required for hydraulic fracture growth. We applied the results to fracture orientation data from image logs to determine the population of pre-existing faults that are expected to slip during stimulation. We also determined microseismic focal plane mechanisms and found slip on steeply dipping planes striking NW, consistent with the orientations of potentially active faults predicted by the stress model. This case study represents a general approach for integrating stress measurements and rock properties to predict hydraulic fracture growth and the characteristics of injection-induced microseismicity.

microseismicity↗
DOE OSTI2020

Applying Waveform Correlation to Mining Blasts Using a Global Sparse Network

Agencies that monitor for underground nuclear tests are interested in techniques that automatically characterize mining blasts to reduce the human analyst effort required to produce high-quality event bulletins. Waveform correlation is effective in finding similar waveforms from repeating seismic events, including mining blasts. We report the results of an experiment that uses waveform templates recorded by multiple International Monitoring System stations of the Preparatory Commission for the Comprehensive Nuclear-Test-Ban Treaty Organization for up to 10 years prior to the time period of interest to detect and identify mining blasts that occur during single weeks of study. We discuss approaches for template selection, threshold setting, and event detection that are specialized for mining blasts and a sparse, global network. We apply the approaches to two different weeks of study for each of two geographic regions, Wyoming and Scandinavia, to evaluate the potential for establishing a set of standards for waveform correlation processing of mining blasts that can be effective for operational monitoring systems with a sparse network. We compare candidate events detected with our processing methods to the Reviewed Event Bulletin of the International Data Centre to develop an intuition about potential reduction in analyst workload.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2024

The Collaborative Seismic Earth Model: Generation 2

Geological interpretations, earthquake source inversions and ground motion modeling, among other applications, require models that jointly resolve crustal and mantle structure. With the second generation of the Collaborative Seismic Earth Model (CSEM2), we present a global multi-resolution tomographic Earth model that serves this purpose. The model evolves through successive regional- and global-scale refinements. While the first generation aggregated regional models, with this study, we ensure consistency between all individual submodels, resulting in a model that accurately explains wave propagation across scales. Recent regional tomographic models were incorporated, comprising continental-scale inversions for Asia and Africa, as well as regional inversions for the Western US, Central Andes, Iran, and Southeast Asia. Across all regional refinements, over 793,000 source-receiver pairs contributed. Moreover, the long-wavelength Earth model (LOWE) introduces large-scale structures outside of pre-existing local refinements. A full-waveform inversion for global anisotropic P-and S-wave speed structure over a total of 194 iterations with a minimum period of 50 s on a large data set of 1 hr of waveform data from 2,423 earthquakes and over 6 million source-receiver pairs ensures that regional updates in the crust and uppermost mantle translate into updates of deeper, global-scale structure. To test the performance of CSEM2, we evaluate waveform fits between observed and synthetic seismograms at 50 s for an independent data set on the global scale, and on the regional scale for lower periods. We accurately simulate waveforms within and across regional refinements, maintaining the original resolution of the submodels embedded in the global framework.

58 GEOSCIENCES↗
DOE OSTI2021

Deep convolutional autoencoders as generic feature extractors in seismological applications

The idea of using a deep autoencoder to encode seismic waveform features and then use them in different seismological applications is appealing. In this paper, we designed tests to evaluate this idea of using autoencoders as feature extractors for different seismological applications, such as event discrimination (i.e., earthquake vs. noise waveforms, earthquake vs. explosion waveforms), and phase picking. These tests involve training an autoencoder, either undercomplete or overcomplete, on a large amount of earthquake waveforms, and then using the trained encoder as a feature extractor with subsequent application layers (either a fully connected layer, or a convolutional layer plus a fully connected layer) to make the decision. By comparing the performance of these newly designed models against the baseline models trained from scratch, we conclude that the autoencoder feature extractor approach may only outperform the baseline under certain conditions, such as when the target problems require features that are similar to the autoencoder encoded features, when a relatively small amount of training data is available, and when certain model structures and training strategies are utilized. The model structure that works best in all these tests is an overcomplete autoencoder with a convolutional layer and a fully connected layer to make the estimation.

58 GEOSCIENCES↗
DOE OSTI2022

Advancing Geophysical Techniques to Image a Stratigraphic Hydrothermal Resource

Sedimentary-hosted geothermal energy systems are permeable structural, structural-stratigraphic, and/or stratigraphic horizons with sufficient temperature for direct use and/or electricity generation. Sedimentary-hosted (i.e., stratigraphic) geothermal reservoirs may be present in multiple locations across the central and eastern Great Basin of the USA, thereby constituting a potentially large base of untapped, economically accessible energy resources. Sandia National Laboratories has partnered with a multi disciplinary group of collaborators to evaluate a stratigraphic system in Steptoe Valley, Nevada using both established and novel geophysical imaging techniques. The goal of this study is to inform an optimized strategy for subsequent exploration and development of this resource and analogous ones. Building from prior Nevada Play Fairway Analysis (PFA), this team is primarily 1) collecting additional geophysical data, 2) employing novel joint geophysical inversion/modeling techniques to update existing 3D geologic models, and 3) integrating the geophysical results to produce a working, geologically constrained thermo-hydrological reservoir model. Prior PFA work highlights Steptoe Valley as a favorable resource basin that likely has both sedimentary and hydrothermal characteristics. However, there remains significant uncertainty on the nature and architecture of the resource(s) at depth, which increases the risk in exploratory drilling. Newly acquired gravity, magnetic, magnetotelluric, and controlled-source electromagnetic data products, in conjunction with new and preexisting geoscientific measurements and observations, are being integrated and evaluated for efficacy in understanding stratigraphic geothermal resources and mitigating exploration risk. Furthermore, the influence of hydrothermal activity on sedimentary-hosted reservoirs in favorable structural settings, and whether fault-controlled systems may locally enhance temperature and permeability in some deep stratigraphic reservoirs, will also be evaluated.

Geothermal, Sedimentary Heat, Geophysics, Seismic,↗
DOE OSTI2022

Foam Fracturing Study for Stimulation Development of Enhanced Geothermal Systems

This is a final technical report for the project: Foam Fracturing Study for Stimulation Development of Enhanced Geothermal Systems (EGS). The goal is to demonstrate the feasibility of foam fracturing in EGS applications. The project, led by Oak Ridge National Laboratory (ORNL), was conducted in collaboration with Temple University. The report describes the research activities with Task 1 at ORNL: foam fracturing testing system development and experimental study on foam fracturing, and Task 2 at Temple University: foam testing and foam characterization. Main findings are: 1. A foam fracturing test system has been developed at ORNL, which can be used to perform foam fracturing under pressure up to 6,000 psi. The system monitors foam density during fracturing online and is capable of testing materials in both monotonic and cyclic (up to 50 Hz) injections. 2. Foam fracturing tests were carried out on Charcoal black granite specimens with a blind borehole to the middle length. Two diameters of blind borehole were tested; G2 series: 9.53 mm and G3 series: 4.76 mm. N2-in-water foam was used with AOS as a surfactant. 3. There was a hole-size effect on fracture initiation pressure. The effect is smaller in the case of foam, which was influenced by the high penetrability of gas in foam. Breakdown pressure showed a behavior just as that of fracture pressure; namely an increased value for small hole samples, while the effect in water fracture was more impressive than in foam fracture. 4. Water mass was reduced in foam fracturing within similar range of breakdown pressures. In G2 series, it was decreased from 10.44 g for water fracturing to 5.17 g, representing more than 50% water reduction. Therefore, there is the potential to reduce water use in EGS stimulation through foam fracturing. 5. Use of cyclic injection has the potential to reduce the breakdown pressure and seismicity in EGS application. Experiments using 4-s cycle period found that specimens can be fractured with a low number of cycles. The fatigue pressure was approximately 64 - 77% of monotonic breakdown pressure for water fracturing and 58 - 94% of the breakdown pressure for foam fracturing. 6. A foam stability testing system has been developed that can test foam at 220 Deg C to 2,000 psi. Tested components of candidate foams included two gases: N2 and CO2; 4 surfactants: AOS, SDS, NP-40 and CTAC; 5 stabilizing agents: guar, bentonite clay, borate salt, silica NPs, and GO. 7. N2 and AOS provided the most stable performance over the tested ranges. Furthermore, the AOS foam with stabilizing agents of guar and borate salt (crosslinker) offered the highest half-life of 20 minutes at 200 Deg C and 1,000 psi. 8. Arrhenius equation and modified power law have been demonstrated to fit well the half-time vs. temperature and pressure data, respectively. These relations can be useful to provide the suggestion for future foam stability study. This submission contains the supporting data developed during the project: 1) A final technical report 2) Granite fracturing data in monotonic and cyclic injections with water and N2 foam Foam performance data in various temperatures and pressures, including half-time, is submitted separately.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Joint Inversion of Body Wave Arrival Times, Surface Wave Dispersion Data and Receiver Functions: Method and Application to South China

In this work, we have developed a new joint inversion method that incorporates body wave arrival times, surface wave dispersion and receiver functions to simultaneously update earthquake locations and constrain three-dimensional P-wave (Vp) and S-wave velocity (Vs) models. Due to complementary sensitivities of the three types of data, the proposed joint inversion algorithm can reduce the intrinsic non-uniqueness of inversions using fewer types of data and better determine smooth velocity variations and velocity discontinuities. Synthetic tests demonstrate the advantages of this new joint inversion algorithm in resolving velocity structures, especially in constraining velocity gradients across the Moho interface. We have applied the proposed joint inversion algorithm to image the lithosphere velocity structure of south China. The inverted Vp and Vs models fit body wave arrival times, surface wave dispersion and receiver functions well. We further analyzed the distribution of Vs gradients across the Moho interface in detail, which helps us better understand tectonics in south China.

58 GEOSCIENCES↗
DOE OSTI2026

Survey of Deep Learning and Physics-Based Approaches in Computational Wave Imaging

Computational wave imaging (CWI) extracts hidden structure and physical properties of a volume of material by analyzing wave signals that traverse that volume. Applications include seismic exploration of the Earth’s subsurface, acoustic imaging and nondestructive testing (NDT) in material science, and ultrasound computed tomography (USCT) in medicine. Current approaches for solving CWI problems can be divided into two categories: those rooted in traditional physics and those based on deep learning. Physics-based methods stand out for their ability to provide high-resolution and quantitatively accurate estimates of acoustic properties within the medium. However, they can be computationally intensive and are susceptible to ill-posedness and nonconvexity typical of CWI problems. Machine learning (ML)-based computational methods have recently emerged, offering a different perspective to address these challenges. Diverse scientific communities have independently pursued the integration of deep learning in CWI. This review discusses how contemporary scientific ML techniques, and deep neural networks in particular, have been developed to enhance and integrate with traditional physics-based methods for solving CWI problems. We present a structured framework that consolidates existing research spanning multiple domains, including computational imaging, wave physics, and data science. This study concludes with important lessons learned from existing ML-based methods and identifies technical hurdles and emerging trends through a systematic analysis of the extensive literature on this topic.

42 ENGINEERING↗
DOE OSTI2022

Thermal equation of state for zoisite: Implications for the transportation of water into the upper mantle and the high-velocity anomaly in the Farallon plate

Seismic tomography studies have revealed a high-velocity anomaly at depths between 100 km and 300 km in the Farallon plate. However, the reasons for the high-velocity anomaly continue to be debated. An analysis of the mineral proportions of eclogites exhumed from the Farallon plate shows that the average amount of zoisite in eclogite is ∼16.0 vol%. Therefore, the presence of zoisite eclogite needs to be considered to explain the high-velocity anomaly of the Farallon plate. However, the thermal equation of state and stability of zoisite have not been fully investigated under high pressure−temperature (P-T) conditions. We investigated the high-pressure and high-temperature behavior of natural zoisite utilizing synchrotron single crystal-X-ray diffraction (XRD). The results indicate that zoisite is metastable up to 24.8 GPa and 700 K. We obtained the ambient unit-cell volume V0 = 901.26(3) Å3 by synchrotron single crystal-XRD measurement. We also fitted the pressure-volume-temperature data to a high-temperature Birch-Murnaghan equation of state and obtained the zero-pressure bulk modulus K0 = 134.7(8) GPa, the temperature derivative of the bulk modulus (∂K/∂T)P = −0.011(4) GPa/K, and the thermal expansion coefficients α0V = 1.9(7) × 10−5 K−1 and α1V = 3(2) × 10−8 K−2 when the pressure derivative of bulk modulus (K′0) is fixed at 4. By incorporating the results from previous studies, we calculated the density and bulk sound velocity profiles of zoisite eclogite along the Farallon plate geotherm. Finally, we infer that zoisite could carry water to depths of ∼300 km within cold subducting slabs and that the bulk sound velocity of typical zoisite eclogite with ∼61.0 vol% omphacite, ∼23.0 vol% garnet, and ∼16.0 vol% zoisite could cause the high-velocity anomalies at depths of 100−300 km in the Farallon plate.

Geology↗
DOE OSTI2023

The Kimberlina synthetic multiphysics dataset for CO 2 monitoring investigations

Abstract We present a synthetic multi‐scale, multi‐physics dataset constructed from the Kimberlina 1.2 CO 2 reservoir model based on a potential CO 2 storage site in the Southern San Joaquin Basin of California. Among 300 models, one selected reservoir‐simulation scenario produces hydrologic‐state models at the onset and after 20 years of CO 2 injection. Subsequently, these models were transformed into geophysical properties, including P‐ and S‐wave seismic velocities, saturated density where the saturating fluid can be a combination of brine and supercritical CO 2 , and electrical resistivity using established empirical petrophysical relationships. From these 3D distributions of geophysical properties, we have generated synthetic time‐lapse seismic, gravity and electromagnetic responses with acquisition geometries that mimic realistic monitoring surveys and are achievable in actual field situations. We have also created a series of synthetic well logs of CO 2 saturation, acoustic velocity, density and induction resistivity in the injection well and three monitoring wells. These were constructed by combining the low‐frequency trend of the geophysical models with the high‐frequency variations of actual well logs collected at the potential storage site. In addition, to better calibrate our datasets, measurements of permeability and pore connectivity have been made on cores of Vedder Sandstone, which forms the primary reservoir unit. These measurements provide the range of scales in the otherwise synthetic dataset to be as close to a real‐world situation as possible. This dataset consisting of the reservoir models, geophysical models, simulated time‐lapse geophysical responses and well logs forms a multi‐scale, multi‐physics testbed for designing and testing geophysical CO 2 monitoring systems as well as for imaging and characterization algorithms. The suite of numerical models and data have been made publicly available for downloading on the National Energy Technology Laboratory's (NETL) Energy Data Exchange (EDX) website.

58 GEOSCIENCES↗
DOE OSTI2021

Hydraulic fracture mapping using wavelet-based fusion of wave transmission and emission measurements

Fracture characterization is essential to hydrocarbon and geothermal exploration and production as well as carbon geo-sequestration. Active and passive measurements have been extensively applied to map fractures across various scales. Active and passive methods have inherent advantages and limitations that complement each other. Here, in this study, to improve the mapping of embedded fractures and the surrounding geomechanically altered regions, we integrate active shear-wave transmission measurement with passive acoustic-emission measurement collected during a lab-scale hydraulic-fracturing experiment. The proposed data fusion leverages the advantages of the two modalities of measurements, while minimizing their individual limitations. Polarized shear-wave transmission (active) measurement was collected before and after the hydraulic fracturing, and the acoustic emission (passive) measurement was collected during hydraulic fracturing. Two sets of two-dimensional maps of fracture and fracture-induced damage in axial, median and frontal planes were obtained by separately processing active and passive measurements. These 2D maps are then fed to wavelet-based image fusion technique to integrate the two sources of information for the reliable mapping/imaging of the embedded fractures and the surrounding geomechanically altered regions.

58 GEOSCIENCES↗