SEARCH · Engineering Papers

Results for “inverse modeling”

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 505 records · Page 28

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Multi‐Scale Seismic Imaging of the Ridgecrest, CA, Region With Waveform Inversion of Regional and Dense Array Data

Abstract We develop an inversion procedure for deriving multi‐scale velocity models with waveform inversions of earthquake and ambient noise data at multi‐frequency bands recorded by regional and dense sensor configurations. The method is applied for the area around the 2019 Ridgecrest earthquake rupture zones, utilizing data recorded by regional stations and dense 2D and 1D arrays with station spacings of ∼5 km and ∼100 m, respectively. Starting with regional Vp, Vs models and locations of Ridgecrest aftershocks, the velocity models and event locations are improved iteratively by inversions of waveforms recorded by regional stations and the 2D array, using a minimum spectral element size of ∼600 m. Waveforms from local events recorded by dense 1D arrays across the M7.1 rupture zone with frequencies of up to 10 Hz are used to resolve small‐scale features of the rupture zone and shallow crust with a local spectral element size of 80 m. The refined models provide self‐consistent descriptions of the rupture zone and the shallow crust embedded in the regional structures. The results reveal pronounced low Vs and high Vp/Vs in the M6.4 and M7.1 rupture zones coinciding with concentrations of seismicity, and also around the Garlock fault and in several local basins. We also observe clear velocity contrasts across the Garlock fault with polarity reversals along strike and with depth. The obtained multi‐scale velocity models can be used to improve derivations of earthquake source properties, simulations of dynamic ruptures and ground motions, and the understanding of fault and tectonic processes in the region.

Geochemistry & Geophysics↗

DOE OSTI2021

Improving Subsurface Stress Characterization for Carbon Dioxide Storage Projects

Researchers at New Mexico Tech, Los Alamos National Laboratory (LANL), Sandia National Laboratory (SNL), and other collaborators are developing a methodology to improve characterization of stress in the subsurface by means of a model-based inversion of reservoir engineering data, time lapse seismic measurements, and microseismicity. Historically, various geophysical techniques have been used in efforts to understand the state of stress in the subsurface through direct and indirect imaging of stress sensitive features (faults and fractures) and observations of transient stress related observations (time variant elastic moduli and microseismicity). Standard direct and indirect techniques for seismic fault and fracture imaging suffer from detectability limits due to reliance on high data quality and multiplicity. Variations in effective elastic properties from inversion of high quality seismic (Vertical Seismic Profile (VSP)) data may be used to infer stress near the wellbore through integration with independent experimental characterization of the stress-velocity relationship. Given a sufficiently robust observation network, microseismic emissions may be inverted to characterize focal mechanisms which, together with supporting assumptions and constraints, inform estimates of in-situ stress. While each of these techniques contributes in part to the characterization of stress within limited spatial and temporal domains, no one method provides an unambiguous stress measurement or a predictive capability over a site scale spatial extent. Challenges associated with solution non-uniqueness, measurement ambiguity, and irregular sampling may be greatly minimized through combination of one or more of independent measurements within a common framework in which realistic geological, hydrodynamic, geomechanical, and seismological constituent process models may act as constraints.

02 PETROLEUM↗

DOE OSTI2021

Seismic Source Modeling Software Enhancements (FY21)

Seismic source modeling allows researchers both to simulate how a source that induces seismic waves interacts with the Earth to produce observed seismograms and, inversely, to infer what the time histories, sizes, and force distributions were for a seismic source given observed seismograms. In this report, we discuss improvements made in FY21 to our software as applies to both the forward and inverse seismic source modeling problems. For the forward portion of the problem, we have added the ability to use full 3-D nonlinear simulations by implementing 3-D time varying boundary conditions within Sandia’s linear seismic code Parelasti. Secondly, on the inverse source modeling side, we have developed software that allows us to invert seismic gradiometer-derived observations in conjunction with standard translational motion seismic data to infer properties of the source that may improve characterization in certain circumstances. First, we describe the basic theory behind each software enhancement and then demonstrate the software in action with some simple examples.

58 GEOSCIENCES↗

NASA NTRS2014

A Decade of Volcanic Observations from Aura and the A-Train

Aura observations have made many seminal contributions to volcanology. Prior to the Aura launch, satellite observations of volcanic degassing (e.g., from TOMS) were mostly restricted to large eruptions. However, the vast majority of volcanic gases are released during quiescent 'passive' degassing between eruptions. The improved sensitivity of Aura OMI permitted the first daily, space-borne measurements of passive volcanic SO2 degassing, providing improved constraints on the source locations and magnitude of global SO2 emissions for input to atmospheric chemistry and climate models. As a result of this unique sensitivity to volcanic activity, OMI data were also the first satellite SO2 measurements to be routinely used for volcano monitoring at several volcano observatories worldwide. Furthermore, the Aura OMI SO2 data also offer unprecedented sensitivity to volcanic clouds in the UTLS, elucidating the transport, fate and lifetime of volcanic SO2 and providing critical input to aviation hazard mitigation efforts. Another major advance has been the improved vertical resolution of volcanic clouds made possible by synergy between Aura and other A-Train instruments (e.g., AIRS, CALIPSO, CloudSat), advanced UV SO2 altitude retrievals, and inverse trajectory modeling of detailed SO2 cloud maps. This altitude information is crucial for climate models and aviation hazards. We will review some of the highlights of a decade of Aura observations of volcanic activity and look ahead to the future of volcanic observations from space.

Aura↗

DOE OSTI2023

Model data for numerical evaluation of photosensitive tracers as a strategy for separating surface and subsurface transient storage in streams

This model-data archive pertains to a study aimed at the numerical evaluation of the photosensitive tracers as a potential strategy for separating the effects of surface and hyporheic storage zones (SSZs and HSZs, respectively). Separating effects of SSZ and HSZ are important for accurately representing stream function as HSZs and SSZs expose solutes to significantly different biogeochemical conditions. We perform numerical experiments using a multiscale reactive transport model for stream corridors implemented in ATS code, which allows for representing multiple storage zones with their respective travel time distributions and biogeochemistry. For each of the numerical experiment, we provide in this model-data archive python wrapper to drive ATS (forward_model.py), synthetic observation (synthetic_btc.py, BTC_observed.csv), forward model files including ATS input (multiscale_transport.tpl), PFLOTRAN inputs (reactions_channel.tpl, reactions_hz.tpl, reactions_sz.tpl), MCMC (mcmc_run.py), predictive uncertainty (pred_uncert.py, BTCs_simulated_day.csv, BTCs_simulated_night.csv), MCMC outputs (tracer_test_GR.npy, tracer_test_logps.npy, tracer_test_parameters.npy) and Jupyter notebook for post-processing and visualization (post-processing.ipynb). For the denitrification application, ATS input file (denitrification_multisubgrid.tpl), PFLOTRAN input files (denitrification_channel.in, denitrification_hz.in), predictive uncertainty (pred_uncert_denitrification.py, BTCs_simulated_DO.csv, BTCs_simulated_DOC.csv, BTCs_simulated_Nitrate.csv).

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Learning viscoelasticity models from indirect data using deep neural networks

In this study, we propose a novel approach to model viscoelasticity materials, where rate-dependent and non-linear constitutive relationships are approximated with deep neural networks. We assume that inputs and outputs of the neural networks are not directly observable, and therefore common training techniques with input–output pairs for the neural networks are inapplicable. To that end, we develop a novel computational approach to both calibrate parametric and learn neural-network-based constitutive relations of viscoelasticity materials from indirect displacement data in the context of multiple-physics systems. We show that limited displacement data holds sufficient information to quantify the viscoelasticity behavior. We formulate the inverse computation – modeling viscoelasticity properties from observed displacement data – as a PDE-constrained optimization problem and minimize the error functional using a gradient-based optimization method. The gradients are computed by a combination of automatic differentiation and implicit function differentiation rules. The effectiveness of our method is demonstrated through numerous benchmark problems in geomechanics and porous media transport.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1982

High-resolution X-ray observations of M87 - Nucleus, jet and radio halo

X-ray observations of M87 made using the Einstein Observatory HRI show the presence of four components. A broad, diffuse component associated with hot gas in the galaxy follows an approximate inverse power law; the pressure of the gas could be sufficient to confine the inner and outer radio lobes. A diffuse, asymmetrical X-ray distribution to the east and southwest appears similar to the radio morphology and is consistent with an inverse Compton model. An X-ray source centered on the nucleus of the galaxy is extended by 2-3 arcsec and may be due to thermal bremsstrahlung; a possible compact source of less than 10 to the 41st ergs/sec is weak compared with other active galaxies, and its implications for the presence of a massive black hole are discussed. A structure correlated with the radio-optical jet is also seen; the X-ray flux is consistent with an extrapolation of the optical synchrotron spectrum. The four components suggest that all three canonical astrophysical X-ray emission mechanisms are present in M87: thermal bremsstrahlung, inverse Compton scattering, and synchrotron radiation.

Schreier, E. J.↗

NASA NTRS2013

Analyzing Carbon Dioxide and Methane Emissions in California Using Airborne Measurements and Model Simulations

Greenhouse gas (GHG) concentrations have increased over the past decades and are linked to global temperature increases and climate change. These changes in climate have been suggested to have varying effects, and uncertain consequences, on agriculture, water supply, weather, sea-level rise, the economy, and energy. To counteract the trend of increasing atmospheric concentrations of GHGs, the state of California has passed the California Global Warming Act of 2006 (AB-32). This requires that by the year 2020, GHG (e.g., carbon dioxide (CO2) and methane (CH4)) emissions will be reduced to 1990 levels. To quantify GHG fluxes, emission inventories are routinely compiled for the State of California (e.g., CH4 emissions from the California Greenhouse Gas Emissions Measurement (CALGEM) Project). The major sources of CO2 and CH4 in the state of California are: transportation, electricity production, oil and gas extraction, cement plants, agriculture, landfills/waste, livestock, and wetlands. However, uncertainties remain in these emission inventories because many factors contributing to these processes are poorly quantified. To alleviate these uncertainties, a synergistic approach of applying air-borne measurements and chemical transport modeling (CTM) efforts to provide a method of quantifying local and regional GHG emissions will be performed during this study. Additionally, in order to further understand the temporal and spatial distributions of GHG fluxes in California and the impact these species have on regional climate, CTM simulations of daily variations and seasonality of total column CO2 and CH4 will be analyzed. To assess the magnitude and spatial variation of GHG emissions and to identify local “hot spots”, airborne measurements of CH4 and CO2 were made by the Alpha Jet Atmospheric eXperiment (AJAX) over the San Francisco Bay Area (SFBA) and San Joaquin Valley (SJV) in January and February 2013 during the Discover-AQ-CA study. High mixing ratios of GHGs were observed in-flight with a high degree of spatial variability. To provide an additional method to quantify GHG emissions, and analyze AJAX measurement data, the GEOS-Chem CTM is used to simulate SFBA/SJV GHG measurements. A nested-grid version of GEOS-Chem will be applied and utilizes varying emission inventories and model parameterizations to simulate GHG fluxes/emissions. The model considers CO2 fluxes from fossil fuel use, biomass/biofuel burning, terrestrial and oceanic biosphere exchanges, shipping and aviation, and production from the oxidation of carbon monoxide, CH4, and non-methane volatile organic carbons. The major sources of CH4 simulated in GEOS-Chem are domesticated animals, rice fields, natural gas leakage, natural gas venting/flaring (oil production), coal mining, wetlands, and biomass burning. Preliminary results from the comparison between available observations (e.g., AJAX and CALGEM CH4 emission maps) and GEOS-Chem results will be presented, along with a discussion of CO2 and CH4 source apportionment and the use of the GEOS-Chem-adjoint to perform inverse GHG modeling.

Johnson, Matthew S.↗

NASA NTRS2011

Using High-Resolution Forward Model Simulations of Ideal Atmospheric Tracers to Assess the Spatial Information Content of Inverse CO2 Flux Estimates

Attribution of observed atmospheric carbon concentrations to emissions on the country, state or city level is often inferred using "inversion" techniques. Such computations are often performed using advanced mathematical techniques, such as synthesis inversion or four-dimensional variational analysis, that invoke tracing observed atmospheric concentrations backwards through a transport model to a source region. It is, to date, not well understood how well such techniques can represent fine spatial (and temporal) structure in the inverted flux fields. This question is addressed using forward-model computations with idealized tracers emitted at the surface in a large number of grid boxes over selected regions and examining how distinctly these emitted tracers can be detected downstream. Initial results show that tracers emitted in half-degree grid boxes over a large region of the Eastern USA cannot be distinguished from each other, even at short distances over the Atlantic Ocean, when they are emitted in grid boxes separated by less than five degrees of latitude - especially when only total-column observations are available. A large number of forward model simulations, with varying meteorological conditions, are used to assess how distinctly three types observations (total column, upper tropospheric column, and surface mixing ratio) can separate emissions from different sources. Inferences inverse modeling and source attribution will be drawn.

Pawson, Steven↗

DOE OSTI2020

Spectroscopic Probe Molecule Selection Using Quantum Theory, First-Principles Calculations, and Machine Learning

Probe molecule vibrational spectra have a long history of being used to characterize materials including metals, oxides, metal-organic frameworks, and even human proteins. Furthermore, recent advances in machine learning have enabled computationally generated spectra to aid in detailed characterization of complex surfaces with probe molecules. Despite widespread use of probe molecules, the science of probe molecule selection is underdeveloped. Here, we develop physical concepts, including orbital interaction energy and the energy overlap integral, to explain and predict the ability of probe molecules to discriminate structural descriptors. We resolve the crystal orbital overlap population (COOP) to specific molecular orbitals and quantify their bonding character, which directly influences vibrational frequencies. Using only a single adsorbate calculation from density function theory (DFT), we compute the interaction energy of individual adsorbate molecular orbitals with adsorption site atomic orbitals across many different sites. Combining the molecular orbital resolved COOP and changes in orbital interaction energy enables probe molecule selection for improved discrimination of various sites. We demonstrate these concepts by comparing the predicted effectiveness of carbon monoxide (CO), nitric oxide (NO), and ethylene (C 2 H 4 ) to probe Pt adsorption sites. Finally, using a previously developed machine learning framework, we show that models trained on hundreds of thousand C 2 H 4 spectra, computed from DFT, which regress surface binding-type and generalized coordination number (GCN), outperform those trained using CO and NO spectra. Lastly, a python package, pDOS_overlap, for implementing the electron density based analysis on any combination of adsorbates and materials, is also made available.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

NASA NTRS2004

Strategy Developed for Selecting Optimal Sensors for Monitoring Engine Health

Sensor indications during rocket engine operation are the primary means of assessing engine performance and health. Effective selection and location of sensors in the operating engine environment enables accurate real-time condition monitoring and rapid engine controller response to mitigate critical fault conditions. These capabilities are crucial to ensure crew safety and mission success. Effective sensor selection also facilitates postflight condition assessment, which contributes to efficient engine maintenance and reduced operating costs. Under the Next Generation Launch Technology program, the NASA Glenn Research Center, in partnership with Rocketdyne Propulsion and Power, has developed a model-based procedure for systematically selecting an optimal sensor suite for assessing rocket engine system health. This optimization process is termed the systematic sensor selection strategy. Engine health management (EHM) systems generally employ multiple diagnostic procedures including data validation, anomaly detection, fault-isolation, and information fusion. The effectiveness of each diagnostic component is affected by the quality, availability, and compatibility of sensor data. Therefore systematic sensor selection is an enabling technology for EHM. Information in three categories is required by the systematic sensor selection strategy. The first category consists of targeted engine fault information; including the description and estimated risk-reduction factor for each identified fault. Risk-reduction factors are used to define and rank the potential merit of timely fault diagnoses. The second category is composed of candidate sensor information; including type, location, and estimated variance in normal operation. The final category includes the definition of fault scenarios characteristic of each targeted engine fault. These scenarios are defined in terms of engine model hardware parameters. Values of these parameters define engine simulations that generate expected sensor values for targeted fault scenarios. Taken together, this information provides an efficient condensation of the engineering experience and engine flow physics needed for sensor selection. The systematic sensor selection strategy is composed of three primary algorithms. The core of the selection process is a genetic algorithm that iteratively improves a defined quality measure of selected sensor suites. A merit algorithm is employed to compute the quality measure for each test sensor suite presented by the selection process. The quality measure is based on the fidelity of fault detection and the level of fault source discrimination provided by the test sensor suite. An inverse engine model, whose function is to derive hardware performance parameters from sensor data, is an integral part of the merit algorithm. The final component is a statistical evaluation algorithm that characterizes the impact of interference effects, such as control-induced sensor variation and sensor noise, on the probability of fault detection and isolation for optimal and near-optimal sensor suites.

Source record↗

DOE OSTI2023

Experimental demonstration of >20 kJ laser energy coupling in 1-cm hydrocarbon-filled gas pipe targets via inverse Bremsstrahlung absorption with applications to MagLIF

Laser propagation experiments using four beams of the National Ignition Facility to deliver up to 35 kJ of laser energy at 351 nm laser wavelength to heat magnetized liner inertial fusion-scale (1 cm-long), hydrocarbon-filled gas pipe targets to ~keV electron temperatures have demonstrated energy coupling >20 kJ with essentially no backscatter in 15% critical electron density gas fills with 0–19 T applied axial magnetic fields. The energy coupling is also investigated for an electron density of 11.5% critical and for applied field strengths up to 24 T at both densities. This spans a range of Hall parameters 0 < ω ce τ ei ≲ 2, where a Hall parameter of 0.5 is expected to reduce electron thermal conduction across the field lines by a factor of 4–5 for the conditions of these experiments. At sufficiently high applied field strength (and therefore Hall parameter), the measured laser propagation speed through the targets increases in the measurements, consistent with reduced perpendicular electron thermal transport; this reduces the coupled energy to the target once the laser burns through the gas pipe. In conclusion, the results compare well with a 1D analytic propagation model for inverse Bremsstrahlung absorption.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Latitudinal gradient in the respiration quotient and the implications for ocean oxygen availability

Significance The loss of ocean oxygen caused by climate warming is a serious environmental issue and could lead to major declines in animal habitats. However, current Earth system models struggle in explaining observed trends in global ocean oxygen. Using a model, we demonstrate that marine oxygen concentrations depend sensitively on a rarely considered quantity known as the respiration quotient to prescribe the oxygen consumed per mole of organic carbon respired. Using a combination of direct chemical measurements across major biomes in the eastern Pacific Ocean plus changes in the global distribution of oxygen and carbon, we show how the respiration quotient increases systematically with temperature. The findings imply that we can experience more intense declines in ocean oxygen with warming.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1994

Visible-infrared remote-sensing model and applications for ocean waters

Remote sensing has become important in the ocean sciences, especially for research involving large spatial scales. To estimate the in-water constituents through remote sensing, whether carried out by satellite or airplane, the signal emitted from beneath the sea surface, the so called water-leaving radiance (L(w)), is of prime importance. The magnitude of L(w) depends on two terms: one is the intensity of the solar input, and the other is the reflectance of the in-water constituents. The ratio of the water-leaving radiance to the downwelling irradiance (E(d)) above the sear surface (remote-sensing reflectance, R(sub rs)) is independent of the intensity of the irradiance input, and is largely a function of the optical properties of the in-water constituents. In this work, a model is developed to interpret r(sub rs) for ocean water in the visible-infrared range. In addition to terms for the radiance scattered from molecules and particles, the model includes terms that describe contributions from bottom reflectance, fluorescence of gelbstoff or colored dissolved organic matter (CDOM), and water Raman scattering. By using this model, the measured R(sub rs) of waters from the West Florida Shelf to the Mississippi River plume, which covered a (concentration of chlorophyll a) range of 0.07 - 50 mg/cu m, were well interpreted. The average percentage difference (a.p.d.) between the measured and modeled R(sub rs) is 3.4%, and, for the shallow waters, the model-required water depth is within 10% of the chart depth. Simple mathematical simulations for the phytoplankton pigment absorption coefficient (a(sub theta)) are suggested for using the R(sub rs) model. The inverse problem of R(sub rs), which is to analytically derive the in-water constituents from R(sub rs) data alone, can be solved using the a(sub theta) functions without prior knowledge of the in-water optical properties. More importantly, this method avoids problems associated with a need for knowledge of the shape and value of the chlorophyll-specific absorption coefficient. The simulation was tested for a wide range of water types, including waters from Monterey Bay, the West Florida Shelf, and the Mississippi River plume. Using the simulation, the R(sub rs)-derived in-water absorption coefficients were consistent with the values from in-water measurements (r(exp 2) greater than 0.94, slope approximately 1.0). In the remote-sensing applications, a new approach is suggested for the estimation of primary production based on remote sensing. Using this approach, the calculated primary production (PP) values based upon remotely sensed data were very close to the measured values for the euphotic zone (r(exp 2) = 0.95, slope 1.26, and 32% average difference), while traditional, pigment-based PP model provided values only one-third the size of the measured data. This indicates a potential to significantly improve the accuracy of the estimation of primary production based upon remote sensing.

Lee, Zhongping↗

NASA NTRS1981

Solar wind electron densities from Viking dual-frequency radio measurements

Simultaneous phase coherent, two-frequency measurements of the time delay between the earth station and the Viking spacecraft have been analyzed in terms of the electron density profiles from 4 solar radii to 200 solar radii. The measurements were made during a period of solar activity minimum (1976-1977) and show a strong solar latitude effect. The data were analyzed with both a model independent, direct numerical inversion technique and with model fitting, yielding essentially the same results. It is shown that the solar wind density can be represented by two power laws near the solar equator proportional to r exp -2.7 and r exp -2.04. However, the more rapidly falling term quickly disappears at moderate latitudes (approximately 20 deg) leaving only the inverse-square behavior.

Muhleman, D. O.↗

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 OSTI2025

Phase II Field Demonstration at Lansing Smith Generating Plant, Southport, Florida (Final Report)

The Final Technical Report: Field Test Design and Pressure Management Strategies for Phase II Field Demonstration of Optimal Pressure Control, Plume Management, and Produced Water Strategies presents the culmination of multi-year efforts under the U.S. Department of Energy’s Brine Extraction and Storage Test (BEST) program, led by EPRI in partnership with Gulf Power at the Plant Smith site near Panama City, Florida. The project was designed to evaluate and demonstrate the technical feasibility of managing subsurface pressures and fluid movement associated with large-scale CO₂ injection, using low-salinity water as a proxy fluid. Through a combination of field injection testing, reservoir modeling, and optimization studies, the research team developed and refined active and passive brine extraction strategies aimed at controlling injection-induced pressure buildup, mitigating risks of fault activation, and managing plume migration. The field demonstration incorporated a new injection well (TIW-2), a new monitoring/extraction well (TEMW-A), and an existing well (TIW-1) repurposed for passive pressure relief. Complementary geophysical monitoring was designed to track plume development and evaluate the effectiveness of pressure management strategies. The report details the integrated workflow encompassing wellfield development, regulatory permitting, model calibration, and survey design. It includes a comprehensive electromagnetic (EM) modeling and inversion study used to develop a cost-effective, time-lapse geophysical monitoring plan capable of imaging the evolving low-salinity plume within the highly saline Lower Tuscaloosa formation. Reservoir simulation results guided the design of the 17-month injection program and the timing of active extraction to maintain formation pressures below a hypothetical fault reactivation threshold. Supporting analyses evaluated potential injectivity risks related to clay fines migration and geochemical incompatibility, identifying practical mitigation measures such as salinity and pH control. The final design recommends a series of crosswell EM surveys—conducted before, during, and after injection—as the most effective and economical approach for plume imaging, supplemented by continuous downhole pressure and flow monitoring. Collectively, the study provides a field-ready framework for cost-effective pressure management and monitoring in support of future CO₂ storage projects.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2022

A Bayesian Approach for In-Situ Stress Prediction and Uncertainty Quantification for Subsurface Engineering

Many subsurface engineering applications require accurate knowledge of the in-situ state of stress for their safe design and operation. Existing methods to meet this need primarily include field measurements for estimating one or more of the principal stresses from a borehole, or optimization methods for constructing a 3D geomechanical model in terms of geophysical measurements. These methods, however, often contain considerable uncertainty in estimating the state of stress. Here, in this paper, we build on a Bayesian approach to quantify uncertainty in stress estimations for subsurface engineering applications. This approach can provide an estimate of the 3D distribution of stress throughout the volume of interest and provide an estimate of the uncertainty arising from the stress measurement, the rheology parameters, and a paucity of measurements. The value of this approach is demonstrated using stress measurements from the In Salah carbon storage site, which was one of the world’s first industrial carbon capture and storage projects. This demonstration shows the application of this Bayesian approach for estimating the initial state of stress for In Salah and quantifying the uncertainty in the estimated stress. Also, an assessment of a maximum injection pressure to prevent geomechanical risks from CO 2 injection pressures is provided in terms of the probability distribution of the minimum principal stress quantified by the approach. With the In Salah case study, this paper demonstrates that using the Bayesian approach can provide additional insights for site explorations and/or project operations to make informed-site decisions for subsurface engineering applications.

58 GEOSCIENCES↗