Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “3D 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.

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

At least 1,441 records · Page 80
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Development and validation of non-axisymmetric heat flux simulations with 3D fields using the HEAT code

A new comprehensive module to simulate heat fluxes from three-dimensional (3D) magnetic fields has been implemented in the HEAT code. Especially compact tokamaks like SPARC require tools to predict and manage large heat fluxes. Existing release versions of HEAT can only simulate axisymmetric heat flux on 3D plasma facing components. The new module uses an M3D-C1 perturbed equilibrium and the MAFOT code to trace field lines of the perturbed 3D magnetic field. Heat flux is then assigned to the resulting footprints via a 3D layer model. The model distinguishes between the scrape-off layer, the magnetic lobes and the private flux region, and employs only 0D parameters like the layer width, diffusive spread and the last closed flux surface position in the perturbed edge to generate a heat flux profile. The magnitude is normalized to the total input power. Resulting heat flux simulations are compared and validated against infrared measurements in the DIII-D tokamak with applied 3D fields; good agreement is found for several cases. The new module can now be applied to the SPARC tokamak; a preliminary result for applied rotating 3D fields is shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRSDate not supplied

Modeling Stellar Jitter and its Effects on the Stellar Spectrum for More Precise Radial Velocity Measurements

Detection of Earth-size exoplanets is a technological and data analysis challenge. Future progress in Earth-mass exoplanet detection is expected from the development of extreme precision radial velocity (EPRV) measurements. The required precision of these measurements necessitates developing a new physics-based data analysis methodology to discriminate planetary signals from host-star-related effects, taking stellar variability and instrumental uncertainties into account. In this work, we investigate and quantify stellar disturbances of the planet-hosting solar-type star HD121504 from 3D radiative modeling obtained with the StellarBox code. The model was used to determine the statistical properties of the turbulent plasma and to obtain synthetic spectroscopic observations for several Fe I lines at different locations on the stellar disk to mimic high-resolution spectroscopic observations. Presented are the initial results of our effort to characterize stellar jitter properties as well as future objectives.

SMD↗
DOE OSTI2014

Evaluation of Cloud-resolving and Limited Area Model Intercomparison Simulations using TWP-ICE Observations. Part 2: Rain Microphysics

Ten 3D cloud-resolving model (CRM) simulations and four 3D limited area model (LAM) simulations of an intense mesoscale convective system observed on January 23-24, 2006 during the Tropical Warm Pool – International Cloud Experiment (TWP-ICE) are compared with each other and with observations and retrievals from a scanning polarimetric radar, co-located UHF and VHF vertical profilers, and a Joss-Waldvogel disdrometer in an attempt to explain published results showing a low bias in simulated stratiform rainfall. Despite different forcing methodologies, similar precipitation microphysics errors appear in CRMs and LAMs with differences that depend on the details of the bulk microphysics scheme used. One-moment schemes produce too many small raindrops, which biases Doppler velocities low, but produces rain water contents (RWCs) that are similar to observed. Two-moment rain schemes with a gamma shape parameter (μ) of 0 produce excessive size sorting, which leads to larger Doppler velocities than those produced in one-moment schemes, but lower RWCs than observed. Two moment schemes also produce a convective median volume diameter distribution that is too broad relative to observations and thus, may have issues balancing raindrop formation, collision coalescence, and raindrop breakup. Assuming a μ of 2.5 rather than 0 for the raindrop size distribution improves one-moment scheme biases, and allowing μ to have values greater than 0 may improve two-moment schemes. Under-predicted stratiform rain rates are associated with under-predicted ice water contents at the melting level rather than excessive rain evaporation, in turn likely associated with convective detrainment that is too high in the troposphere and mesoscale circulations that are too weak. In addition to stronger convective updrafts than observed, limited domain size prevents a large, well-developed stratiform region from developing in CRMs, while a dry bias in ECMWF analyses does the same to the LAMs.

cloud-resolving↗
DOE OSTI2016

ON THE OBSERVATION AND SIMULATION OF SOLAR CORONAL TWIN JETS

We present the first observation, analysis, and modeling of solar coronal twin jets, which occurred after a preceding jet. Detailed analysis on the kinetics of the preceding jet reveals its blowout-jet nature, which resembles the one studied in Liu et al. However, the erupting process and kinetics of the twin jets appear to be different from the preceding one. Lacking detailed information on the magnetic fields in the twin jet region, we instead use a numerical simulation using a three-dimensional (3D) MHD model as described in Fang et al., and find that in the simulation a pair of twin jets form due to reconnection between the ambient open fields and a highly twisted sigmoidal magnetic flux, which is the outcome of the further evolution of the magnetic fields following the preceding blowout jet. Based on the similarity between the synthesized and observed emission, we propose this mechanism as a possible explanation for the observed twin jets. Combining our observation and simulation, we suggest that with continuous energy transport from the subsurface convection zone into the corona, solar coronal twin jets could be generated in the same fashion addressed above.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2013

Flight Deck Weather Avoidance Decision Support: Implementation and Evaluation

Weather related disruptions account for seventy percent of the delays in the National Airspace System (NAS). A key component in the weather plan of the Next Generation of Air Transportation System (NextGen) is to assimilate observed weather information and probabilistic forecasts into the decision process of flight crews and air traffic controllers. In this research we explore supporting flight crew weather decision making through the development of a flight deck predicted weather display system that utilizes weather predictions generated by ground-based radar. This system integrates and presents this weather information, together with in-flight trajectory modification tools, within a cockpit display of traffic information (CDTI) prototype. that the CDTI features 2D and perspective 3D visualization models of weather. The weather forecast products that we implemented were the Corridor Integrated Weather System (CIWS) and the Convective Weather Avoidance Model (CWAM), both developed by MIT Lincoln Lab. We evaluated the use of CIWS and CWAM for flight deck weather avoidance in two part-task experiments. Experiment 1 compared pilots' en route weather avoidance performance in four weather information conditions that differed in the type and amount of predicted forecast (CIWS current weather only, CIWS current and historical weather, CIWS current and forecast weather, CIWS current and forecast weather and CWAM predictions). Experiment 2 compared the use of perspective 3D and 21/2D presentations of weather for flight deck weather avoidance. Results showed that pilots could take advantage of longer range predicted weather forecasts in performing en route weather avoidance but more research will be needed to determine what combinations of information are optimal and how best to present them.

CWAM↗
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 NTRS1993

Some aspects of the interaction between chemical and dynamic processes relating to the Antarctic ozone hole

Observational and modeling studies have been conducted to examine the interaction between the chemical and dynamical processes that occur during springtime in the lower stratosphere of the Southern Hemisphere. The temporal evolution of the ozone distribution and the circulation during 1987 is contrasted with that for 1988 as an illustrative example of how dynamical processes and the resulting meteorological conditions modulate the ozone depletion. Concurrently with the observational analysis, an effort was initiated to simulate the ozone depletion during austral spring using a 3D chemical/transport model. The model includes a parameterized representation of the heterogeneous processes thought to be important in this region. The simulation indicates that the inclusion of this additional chemistry, which results in the release of free chlorine and the redistribution of odd nitrogen into reservoir species, reproduces many aspects of the observations.

Eckman, R. S.↗
DOE OSTI2023

Temperature Uncertainty Modeling with Proxy Structural Data as Geostatistical Constraints for Well Siting: An Example Applied to Granite Springs Valley, NV, USA

Utilizing existing temperature and structural information around Granite Springs Valley, Nevada, we build 3D stochastic temperature models with the aim of evaluating the 3D uncertainty of temperature and choosing between candidate exploration well locations . The data used to support the modeling are measured temperatures and structural proxies from 3D geologic modeling, the latter considered "secondary" data. Two stochastic geostatistical techniques are explored for incorporating the structural proxies: cosimulation and local varying mean. With both the cosimulation and local varying mean methods, many equally likely temperature models (i.e., realizations) are produced, from which temperature probability profiles are calculated at candidate well locations. To aid in choosing between the candidate locations, two quantities summarize the temperature probabilities: Vprior and entropy. Vprior quantifies the likelihood for economic temperatures at each candidate location, whereas entropy identifies where new information has the most potential to reduce uncertainty. In general, the cosimulation realizations have smoother spatial structure, and extrapolate high temperatures at candidate locations that are located along the direction of the longest spatial correlation, which are down dip from existing temperature logs. The smooth realizations result in tight temperature probability profiles that are easier to interpret, but they have unrealistic temperature reversals in some locations because the cosimulation technique does not enforce a conductive geothermal gradient as a baseline (i.e., linearly increasing temperature with depth). The local varying mean results produce realizations with more realistic geothermal gradients, with temperatures increasing downward since a depth-temperature relationship is included. However, because they have much noisier spatial nature compared to cosimulation, it is harder to interpret the temperature probability profiles. The different local varying mean results allow the geologist to determine which proxy (e.g., dilation versus distance to fault termination) should be used given the specific geothermal system. In general, Vprior from local varying mean results identify locations that are close to high values for the structural proxies: areas with highe r probabilities for higher temperatures. The entropy results identify where uncertainty is greatest and therefore new drilling information could be most useful. Though these techniques provide useful information, even when applied to areas of sparse data, our comp arison of these two techniques demonstrates the need for new geothermal geostatistics techniques that combine the advantages of these two methods and that are tailored to the spatial uncertainty issues inherent in geothermal exploration.

3D temperature modeling↗
DOE OSTI2011

EFFECTS OF ROTATION ON STOCHASTICITY OF GRAVITATIONAL WAVES IN THE NONLINEAR PHASE OF CORE-COLLAPSE SUPERNOVAE

By performing three-dimensional (3D) simulations that demonstrate the neutrino-driven core-collapse supernovae aided by the standing accretion shock instability (SASI), we study how the spiral modes of the SASI can impact the properties of the gravitational-wave (GW) emission. To see the effects of rotation in the nonlinear postbounce phase, we give a uniform rotation on the flow advecting from the outer boundary of the iron core, the specific angular momentum of which is assumed to agree with recent stellar evolution models. We compute fifteen 3D models in which the initial angular momentum and the input neutrino luminosities from the protoneutron star are changed in a systematic manner. By performing a ray-tracing analysis, we accurately estimate the GW amplitudes generated by anisotropic neutrino emission. Our results show that the gravitational waveforms from neutrinos in models that include rotation exhibit a common feature; otherwise, they vary much more stochastically in the absence of rotation. The breaking of the stochasticity stems from the excess of the neutrino emission parallel to the spin axis. This is because the compression of matter is more enhanced in the vicinity of the equatorial plane due to the growth of the spiral SASI modes, leading to the formation of the spiral flows circulating around the spin axis with higher temperatures. We point out that recently proposed future space interferometers like Fabry-Perot-type DECIGO would permit the detection of these signals for a Galactic supernova.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

De novo atomic protein structure modeling for cryoEM density maps using 3D transformer and HMM

Accurately building 3D atomic structures from cryo-EM density maps is a crucial step in cryo-EM-based protein structure determination. Converting density maps into 3D atomic structures for proteins lacking accurate homologous or predicted structures as templates remains a significant challenge. Here, we introduce Cryo2Struct, a fully automated de novo cryo-EM structure modeling method. Cryo2Struct utilizes a 3D transformer to identify atoms and amino acid types in cryo-EM density maps, followed by an innovative Hidden Markov Model (HMM) to connect predicted atoms and build protein backbone structures. Cryo2Struct produces substantially more accurate and complete protein structural models than the widely used ab initio method Phenix. Additionally, its performance in building atomic structural models is robust against changes in the resolution of density maps and the size of protein structures.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Fuel properties to enable lifted-flame combustion

The Fuel Properties to Enable Lifted-Flame Combustion project responded directly to solicitation DE-FOA-0000239 AOI 1A, Fuels and Lubricants for Advanced Combustion Regimes. This subtopic was intended to encompass clean and highly-efficient, liquid-fueled combustion engines to achieve extremely low engine-out nitrogen oxides (NOx) and particulate matter (PM) as a target and similar efficiency as state-of-the-art direct injection diesel engines. The intent of this project was to identify how fuel properties can be used to achieve controllable Leaner Lifted Flame Combustion (LLFC) with low NOx and PM emissions. Specifically, this project was expected to identify and test key fuel properties to enable LLFC and their compatibility with current fuel systems and to enhance combustion models to capture the effect of fuel properties on advanced combustion. Successful demonstration of LLFC may reduce the need for after treatment devices, thereby reducing costs and improving thermal efficiency. The project team consisted of key technical personnel from Ford Motor Company (FMC), the University of Wisconsin-Madison (UW), Sandia National Laboratories (SNL) and Lawrence Livermore National Laboratories (LLNL). Each partner had key roles in achieving project objectives. FMC investigated fuel properties relating to LLFC and sooting tendency. Together, FMC and UW developed and integrated 3D combustion models to capture fuel property combustion effects. FMC used these modeling results to develop a combustion system and define fuel properties to support a single-cylinder demonstration of fuel-enabled LLFC. UW investigated modeling the flame characteristics and emissions behavior of different fuels, including those with different cetane number and oxygen content. SNL led spray combustion experiments to quantify the effect of key fuel properties on combustion characteristics critical for LLFC, as well as single cylinder optical engine experiments to improve fundamental understanding of flame lift-off, generate model validation data, and demonstrate LLFC concurrent with FMC efforts. Additionally, LLNL was added to the project during the second year to develop a detailed kinetic mechanism for a key oxygenate to support CFD modeling. Successful completion of this project allowed the team to enhance fundamental understanding of LLFC, improve the state of current combustion models and increase understanding of desired fuel properties. This knowledge also improves our knowledge of how cost effective and environmentally friendly renewable fuels can assist in helping meet future emission and greenhouse gas regulations.

Kurtz, Eric↗
DOE OSTI2021

Exosphere Modeling of Proxima b: A Case Study of Photochemical Escape with a Venus-like Atmosphere

Exoplanets orbiting M dwarfs within habitable zones are exposed to stellar environments more extreme than that terrestrial planets experience in our solar system, which can significantly impact the atmospheres of the exoplanets and affect their habitability and sustainability. This study provides the first prediction of hot oxygen corona structure and the associated photochemical loss from a 1 bar CO 2 -dominated atmosphere of a Venus-like rocky exoplanet, where dissociative recombination of O 2 + ions is assumed to be the major source reaction for the escape of neutral O atoms and formation of the hot O corona (or exospheres) as on Mars and Venus. We employ a 3D Monte Carlo code to simulate the exosphere of Proxima Centauri b (PCb) based on the ionosphere simulated by a 3D magnetohydrodynamic model. Our simulation results show that variability of the stellar wind dynamic pressure over one orbital period of PCb does not affect the overall spatial structure of the hot O corona but contributes to the change in the global hot O escape rate that varies by an order of magnitude. The escape increases dramatically when the planet possesses its intrinsic magnetic fields as the ionosphere becomes more extended with the presence of a global magnetic field. Furthermore, the extended hot O corona may lead to a more extended H exosphere through collisions between thermal H and hot O, which exemplifies the importance of considering nonthermal populations in exospheres to interpret future observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Computing the Central Charge of the 3D Ising CFT Using Quantum Finite Elements

The 3D Ising conformal field theory (CFT) describes different physical systems, such as uniaxial magnets or fluids, at their critical points. In absence of an analytical solution for the 3D Ising model, the scaling dimensions and operator product expansion (OPE) coefficients characterizing this CFT must be determined numerically. The currently most-cited values for these quantities have been obtained from the conformal bootstrap, while lattice calculations have so far only produced reliable results for the scaling dimensions involved in calculating the critical exponents. Using Quantum Finite Elements to investigate critical \(\phi^4\)-theory on \(\mathbb{R}\times\mathbb{S}^2\), we show in this work that it is possible to extract scaling dimensions and OPE coefficients of the 3D Ising CFT by fitting the lattice four-point function with expectations from the operator product expansion for the radially quantized CFT and extrapolating to the continuum limit. This way, we have for the first time been able to use Monte Carlo simulations to compute the central charge of the theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Warpage-Resistant, Under-Extrusion-Free, High-Surface-Quality Additive Manufacturing Process for Polyethylene-Based Composite Radiation Shielding Material

Polyethylene (PE) is one of the best shielding materials for primary space radiation due to its high hydrogen content. For effective secondary neutron shielding, boron-rich fillers are incorporated to enhance performance. The semicrystalline nature and high thermal expansion coefficient of PE impede its adoption for in situ additive manufacture in space via the fused deposition modeling (FDM) 3D printing. Here, we developed an optimized PE blend to mitigate the effects of under-extrusion and warpage. Guided by studies on extrusion and warpage, we developed an optimal set of printing parameters for the proposed PE blend. The optimum PE blend─both in its pure form and when doped with fillers─has been tested on different FDM printers. The printed structures exhibit high and uniform density, smooth surfaces, no warpage, and competitive mechanical properties. The FDM-printed plates demonstrate efficient shielding from thermal neutrons, predicted via modeling and confirmed experimentally using extended Q-range small-angle neutron scattering.

additive manufacturing↗
DOE OSTI2013

IS THE SUN LIGHTER THAN THE EARTH? ISOTOPIC CO IN THE PHOTOSPHERE, VIEWED THROUGH THE LENS OF THREE-DIMENSIONAL SPECTRUM SYNTHESIS

We consider the formation of solar infrared (2-6 {mu}m) rovibrational bands of carbon monoxide (CO) in CO5BOLD 3D convection models, with the aim of refining abundances of the heavy isotopes of carbon ({sup 13}C) and oxygen ({sup 18}O, {sup 17}O), to compare with direct capture measurements of solar wind light ions by the Genesis Discovery Mission. We find that previous, mainly 1D, analyses were systematically biased toward lower isotopic ratios (e.g., R {sub 23} {identical_to} {sup 12}C/{sup 13}C), suggesting an isotopically 'heavy' Sun contrary to accepted fractionation processes that were thought to have operated in the primitive solar nebula. The new 3D ratios for {sup 13}C and {sup 18}O are R {sub 23} = 91.4 {+-} 1.3 (R {sub Circled-Plus} = 89.2) and R {sub 68} = 511 {+-} 10 (R {sub Circled-Plus} = 499), where the uncertainties are 1{sigma} and 'optimistic'. We also obtained R {sub 67} = 2738 {+-} 118 (R {sub Circled-Plus} = 2632), but we caution that the observed {sup 12}C{sup 17}O features are extremely weak. The new solar ratios for the oxygen isotopes fall between the terrestrial values and those reported by Genesis (R {sub 68} = 530, R {sub 67} = 2798), although including both within 2{sigma} error flags, and go in the direction favoring recent theories for the oxygen isotope composition of Ca-Al inclusions in primitive meteorites. While not a major focus of this work, we derive an oxygen abundance, {epsilon}{sub O} {approx} 603 {+-} 9 ppm (relative to hydrogen; log {epsilon} {approx} 8.78 on the H = 12 scale). The fact that the Sun is likely lighter than the Earth, isotopically speaking, removes the necessity of invoking exotic fractionation processes during the early construction of the inner solar system.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2006

"One-Stop Shopping" for Ocean Remote-Sensing and Model Data

OurOcean Portal 2.0 (http:// ourocean.jpl.nasa.gov) is a software system designed to enable users to easily gain access to ocean observation data, both remote-sensing and in-situ, configure and run an Ocean Model with observation data assimilated on a remote computer, and visualize both the observation data and the model outputs. At present, the observation data and models focus on the California coastal regions and Prince William Sound in Alaska. This system can be used to perform both real-time and retrospective analyses of remote-sensing data and model outputs. OurOcean Portal 2.0 incorporates state-of-the-art information technologies (IT) such as MySQL database, Java Web Server (Apache/Tomcat), Live Access Server (LAS), interactive graphics with Java Applet at the Client site and MatLab/GMT at the server site, and distributed computing. OurOcean currently serves over 20 real-time or historical ocean data products. The data are served in pre-generated plots or their native data format. For some of the datasets, users can choose different plotting parameters and produce customized graphics. OurOcean also serves 3D Ocean Model outputs generated by ROMS (Regional Ocean Model System) using LAS. The Live Access Server (LAS) software, developed by the Pacific Marine Environmental Laboratory (PMEL) of the National Oceanic and Atmospheric Administration (NOAA), is a configurable Web-server program designed to provide flexible access to geo-referenced scientific data. The model output can be views as plots in horizontal slices, depth profiles or time sequences, or can be downloaded as raw data in different data formats, such as NetCDF, ASCII, Binary, etc. The interactive visualization is provided by graphic software, Ferret, also developed by PMEL. In addition, OurOcean allows users with minimal computing resources to configure and run an Ocean Model with data assimilation on a remote computer. Users may select the forcing input, the data to be assimilated, the simulation period, and the output variables and submit the model to run on a backend parallel computer. When the run is complete, the output will be added to the LAS server for

Li, P. Peggy↗
DOE OSTI2020

Charged Wellbore Casing Controlled Source Electromagnetics (CWC-CSEM) for Reservoir Imaging and Monitoring (Final report)

This project addresses the needs of the U.S. Department of Energy (DOE) to develop advanced monitoring technologies and protocols to track the fate of subsurface carbon dioxide (CO2) plumes for carbon storage. Specifically, the project seeks to develop and test a unique and novel system of technologies consisting of electromagnetic data acquisition, coupled multiphysics imaging, and reservoir model enhancement to understand the migration and long-term distribution of CO2 in the subsurface. The overarching objective is to develop an integrated approach for long term monitoring of carbon storage. The two main components of the project include the methodology development and the test of the method at a field site. The methodology component consists of 1) developing the field procedure and protocol for collecting time-lapse controlled-source electromagnetic (CSEM) data with source electric current injected into the subsurface through wellbore casings; 2) building of background 3D electrical conductivity utilizing multiple sources of data such as supplemental surface transient EM (TEM) surveys, well-logs, and seismic structural information, for enhancing CSEM signal from reservoir depths; and 3) coupled multiphysics simulations and inversion of CSEM data constrained by production data and by structural information from seismic imaging of the reservoir and overlying formations. The testing component used the field site of Bell Creek Oil Field, which served both as a field laboratory for the method development as well as a test site to evaluate the CSEM signal strengths and the methodology developed in this research project. We have accomplished all the proposed tasks and developed the methodology as planned. These include the procedure for time-lapse CSEM data acquisition, data processing techniques, integration with 3D conductivity model building, fast reservoir simulation for history matching using machine learning, and interpreting CSEM data with coupling to the reservoir modeling. Collectively, the outcome of these tasks form a coherent workflow that can be applied to monitor dedicated carbon storage in saline reservoirs. The testing component evaluated the applicability and limitations of the method, and concluded that the method would be ideal for monitoring dedicated carbon storage sites utilizing saline reservoirs.

54 ENVIRONMENTAL SCIENCES↗