Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reservoir Simulation”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

Imaging a Free-Geometry Magma Reservoir Using FEM-Based Linear Inverse Modeling of InSAR Data: Application to Rabaul Caldera, Papua New Guinea

We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluidfilled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.

Inverse theory↗
DOE OSTI2024

Carbonated Brine Injection for a Pilot Site as a Low-Risk Geologic Carbon Storage Strategy: A Simulation Study

This technical report presents a preliminary study for CBI implementation in the Pena Creek pilot site in Dimmit County, Texas, targeting the Edwards Limestone Formation for SWD, in collaboration with NGL Energy Partners. The objective of this study is to find the conservative CO 2 molality for safe storage and potential risks for CBI at this pilot site. Geomodel was performed to get site-specific reservoir geology information. Reservoir modeling was performed to investigate dissolved CO 2 molality and pressure changes at relevant injection conditions. Wellbore corrosion and compatibility was discussed with corrosion modeling. Dissolved CO 2 long-term fate in the reservoir was examined with reactive transport modeling. This plan would set the scientific foundation as the initial step for future implementation in the field.

58 GEOSCIENCES↗
DOE OSTI2025

Carbon Storage in Fold‐and‐Thrust Belts: An Overlooked Gigatonne Storage Opportunity

This study presents numerical investigations of the trapping characteristics of fold-and-thrust belt structures, defining three carbon capture and storage (CCS) play types that could be used to store commercial volumes (millions of tonnes) of CO 2 . Specifically, we present simulations of CO 2 storage in three fold-and-thrust belt models comprising a thrust-ramp, duplex, and thrust-fold geometry. To constrain these play types in realistic geology, each model is based on a study site, including a novel investigation of a greenfield saline reservoir in Virginia, USA, being considered for commercial carbon storage and two well-characterized petroleum fields: the Wilburton field in Oklahoma, USA, and the Incahuasi field in Bolivia. Our results provide insight into several key parameters, such as the long-term security of injected CO 2 in these geologies and injection strategies for maximizing storage efficiency while reducing pressure-related risk. These results improve the understanding of CCS in fold-and-thrust belt storage sites globally by describing general storage parameters that may be applied to site-specific projects. We find that thrust-ramp geometries may securely trap CO 2 through solubility and hydrodynamic trapping under suitable reservoir conditions, duplex structures may store some quantities of CO 2 but are pressure-constrained, and that thrust-ramp structures may store large quantities of CO 2 by maximizing fetch volume, which simultaneously lowers geomechanical risk by reducing pressure buildup along zones of weakness.

42 ENGINEERING↗
DOE OSTI2025

Technical Report on Subsurface Monitoring of the Brady Hot Spring Geothermal Site, Nevada, based upon Full Waveform Inversion

Abilities to accurately characterize the subsurface in a geothermal setting is key to assess and support production. An important element of geothermal reservoir monitoring is also the ability to investigate fluid transport within fracture network. This report focuses on improving subsurface imaging and monitoring in geothermal settings using full waveform inversion based on the adjoint method and time-lapse imaging. To assess our method, we rely on a dense seismic dataset collected in 2016 at the Brady Hot Springs geothermal site in Nevada for the DOE-funded project Poroelastic Tomography by Adjoint Inverse Modeling of Data from Seismology, Geodesy, and Hydrology. This dataset captures subsurface changes across four stages of geothermal power plant operations, which involve varying rates of fluid injection and extraction. Two velocity models were previously derived from this dataset using different methods: one based on travel times and another on sweep interferometry. Our first step is to refine these models using adjoint tomography, which has been applied successfully at global and regional-scales but is less common at the reservoir-scale. Two approaches are then explored for time-lapse analysis: directly comparing refined tomographic models from different stages or backpropagating waveform differences relative to a baseline tomographic model. The main take away is that both approaches highlight similar reservoir behaviors, but the latter approach is more computationally effective in capturing small-scale changes in subsurface properties. For this work, we leverage the use of Salvus (www.mondaic.com), an end-to-end seismic imaging solution, relying on the spectral element method to compute forward and adjoint simulations, and developed by Mondaic Ltd. It includes integrated workflow management that handles waveform and metadata, launches simulations, computes waveform misfits and adjoint sources, and iterates for model updates by nonlinear optimization.

15 GEOTHERMAL ENERGY↗
NASA NTRS1972

Pressure ramp programmer; IMBLMS Phase B4 Additional Tasks: Task 3.0 pressure ramp programmer

A pressure ramp programmer model was designed, fabricated and tested. This model, in conjunction with an automatic blood pressure monitor, automatically controls the pressure in the blood pressure monitor arterial cuff. The cuff pressurization cycle is designed to maximize accuracy and repeatability of blood pressure measurements. The key feature of this automatic cycle is rapid blood pressure cuff bleed down from an initial setting until systolic (diastolic) pressure is encountered followed by a short repressurization and slow bleed, long enough to permit accurate systolic (diastolic) pressure determination. The system includes a pressure reservoir which bleeds the cuff through a precision needle valve; a solenoid valve which permits rapid pressurization from the reservoir; and a pressure sensor which provides information for bleed rate and set point controls. Korotkoff sound signals from a microphone in the blood pressure cuff (not part of the system) provide decision information to the digital control system. The system completed a series of engineering tests using simulated Korotkoff sound inputs. The system performed successfully in all cases and was stable over an extended period of time.

Fogal, G. L.↗
NASA NTRS2017

Imaging the Complex Geometry of a Magma Reservoir Using Fem-Based Linear Inverse Modeling of InSAR Data: Application to Rabaul Caldera, Papua New Guinea

We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluid-filled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.

Inverse theory↗
DOE OSTI2024

5-2428: Coupled Investigation of Fracture Permeability Impact on Reservoir Stress and Seismic Slip Behavior [Slides]

Our goal is to develop, apply and validate a holistic thermal, hydrologic, mechanical, and chemical (THMC) workflow that also includes evaluation of induced seismic slip in EGS reservoirs. We will integrate experimental and modelling approaches to reduce parameter uncertainty and better integrate friction constitutive laws for faults to improve our ability to predict and mitigate seismic hazard at Utah FORGE and future EGS sites. We propose a novel approach that incorporates 3D physics-based Earthquake simulations in THMC models, using lab measurements of friction parameters, herein referred to as “THMC+E” models. These simulations will enable exploration of various circumstances that may hinder EGS success and develop mitigation strategies.

58 GEOSCIENCES↗
NASA NTRS1985

Transient Analysis of a Magnetic Heat Pump

An experimental heat pump that uses a rare earth element as the refrigerant is modeled using NASTRAN. The refrigerant is a ferromagnetic metal whose temperature rises when a magnetic field is applied and falls when the magnetic field is removed. The heat pump is used as a refrigerator to remove heat from a reservoir and discharge it through a heat exchanger. In the NASTRAN model the components modeled are represented by one-dimensional ROD elements. Heat flow in the solids and fluid are analyzed. The problem is mildly nonlinear since the heat capacity of the refrigerant is temperature-dependent. One simulation run consists of a series of transient analyses, each representing one stroke of the heat pump. An auxiliary program was written that uses the results of one NASTRAN analysis to generate data for the next NASTRAN analysis.

Schroeder, E. A.↗
NASA NTRS1985

Laboratory studies on the reactions between chlorine, sulfur dioxide, and oxygen - Implications for the Venus stratosphere

Fourier transform IR spectrophotometry is used to monitor the reactants and products in a Venus stratosphere simulation study involving the photolysis of mixtures of Cl2 and SO2, with and without O2 present in an atmosphere of N2. When several speculative reactions inferred from these experiments are incorporated by the Yung and DeMore (1982) model of Venus stratospheric chemistry, it emerges that SO2Cl2 is a key reservoir species for chlorine, and that the reaction between Cl and SO2 furnishes an important cycle for the destruction of O2 and the conversion of SO2 to H2SO4, thereby providing a possible solution to the photochemistry of the Venus stratosphere.

Demore, W. B.↗
DOE OSTI2026

Complex Electrical Conductivity of a Single‐Fractured Rock: Fracture‐and‐Matrix Coupling Mechanism and Aperture Size Predictions

Fractured rocks play a crucial role in myriad natural and engineered systems, including Earth's critical zone, oil/gas/geothermal reservoirs, and geological CO 2 /H 2 /waste storage systems. While complex electrical conductivity is extensively used to estimate the pore and grain sizes of conventional porous rocks and soils, it is rarely used to predict the aperture size of fractured rocks and this remains poorly understood. Here, integrating theory, simulations, and experiments, we show that under external fields, fractured rocks follow the fracture-and-matrix coupling to make the bulk complex conductivity non-linear with respect to water conductivity. We find that the relaxation time and quadrature conductivity for porous media do not apply to fractured rocks, but, instead, reasonably accurate predictions of aperture size can be made based on the true formation factor. This study unravels the fundamental mechanism governing conduction and polarization of fractured rocks and paves the way for the non-invasive investigation of global fractured rocks.

58 GEOSCIENCES↗
NASA NTRS1997

Timescales of Land Surface Evapotranspiration Response

Soil and vegetation exert strong control over the evapotranspiration rate, which couples the land surface water and energy balances. A method is presented to quantify the timescale of this surface control using daily general circulation model (GCM) simulation values of evapotranspiration and precipitation. By equating the time history of evaporation efficiency (ratio of actual to potential evapotranspiration) to the convolution of precipitation and a unit kernel (temporal weighting function), response functions are generated that can be used to characterize the timescales of evapotranspiration response for the land surface model (LSM) component of GCMS. The technique is applied to the output of two multiyear simulations of a GCM, one using a Surface-Vegetation-Atmosphere-Transfer (SVAT) scheme and the other a Bucket LSM. The derived response functions show that the Bucket LSM's response is significantly slower than that of the SVAT across the globe. The analysis also shows how the timescales of interception reservoir evaporation, bare soil evaporation, and vegetation transpiration differ within the SVAT LSM.

Scott, Russell↗
NASA NTRS2015

Core Analysis Combining MT (TIPPER) and Dielectric Sensors (Sans EC) in Earth and Space

On terrestrial planets and moons of our solar system cores reveal details about a geological structure's formation, content, and history. The strategy for the search for life is focused first on finding water which serves as a universal solvent, and identifying the rocks which such solvent act upon to release the constituent salts, minerals, ferrites, and organic compounds and chemicals necessary for life. Dielectric spectroscopy measures the dielectric properties of a medium as a function of frequency. Reflection measurements in the frequency range from 300 kHz to 300 MHz were carried out using RF and microwave network analyzers interrogating SansEC Sensors placed on clean geological core samples. These were conducted to prove the concept feasibility of a new geology instrument useful in the field and laboratory. The results show that unique complex frequency spectra can be acquired for a variety of rock core samples. Using a combination of dielectric spectroscopy and computer simulation techniques the magnitude and phase information of the frequency spectra can be converted to dielectric spectra. These low-frequency dielectric properties of natural rock are unique, easily determined, and useful in characterizing geology. TIPPER is an Electro-Magnetic Passive-Source Geophysical Method for Detecting and Mapping Geothermal Reservoirs and Mineral Resources. This geophysical method uses distant lightning and solar wind activity as its energy source. The most interesting deflections are caused by the funneling of electrons into more electrically conductive areas like mineralized faults, water or geothermal reservoirs. We propose TIPPER to be used with SansEC for determining terrain/ocean chemistry, ocean depth, geomorphology of fracture structures, and other subsurface topography characteristics below the ice crust of Jovian moons. NASA envisions lander concepts for exploration of these extraterrestrial icy surfaces and the oceans beneath. One such concept would use a nuclear powered heated tip for melting through the ice sheath of Europa and inserting a down hole SansEC with TIPPER interface. NASA's Juno space probe already on the way to Jupiter as part of the Exploration New Frontiers Program and the planned Europa mission will conduct detailed reconnaissance of Jupiter's moon Europa and investigate whether the icy moon could harbor conditions suitable for life. It has already been observed that Jovian moons have auroras that may serve as naturally occurring active energy sources for a TIPPER instrument.

Mound, Michael C.↗
NASA NTRS2015

Inconsistent Strategies to Spin up Models in CMIP5: Implications for Ocean Biogeochemical Model Performance Assessment

During the fifth phase of the Coupled Model Intercomparison Project (CMIP5) substantial efforts were made to systematically assess the skill of Earth system models. One goal was to check how realistically representative marine biogeochemical tracer distributions could be reproduced by models. In routine assessments model historical hindcasts were compared with available modern biogeochemical observations. However, these assessments considered neither how close modeled biogeochemical reservoirs were to equilibrium nor the sensitivity of model performance to initial conditions or to the spin-up protocols. Here, we explore how the large diversity in spin-up protocols used for marine biogeochemistry in CMIP5 Earth system models (ESMs) contributes to model-to-model differences in the simulated fields. We take advantage of a 500-year spin-up simulation of IPSL-CM5A-LR to quantify the influence of the spin-up protocol on model ability to reproduce relevant data fields. Amplification of biases in selected biogeochemical fields (O2, NO3, Alk-DIC) is assessed as a function of spin-up duration. We demonstrate that a relationship between spin-up duration and assessment metrics emerges from our model results and holds when confronted with a larger ensemble of CMIP5 models. This shows that drift has implications for performance assessment in addition to possibly aliasing estimates of climate change impact. Our study suggests that differences in spin-up protocols could explain a substantial part of model disparities, constituting a source of model-to- model uncertainty. This requires more attention in future model intercomparison exercises in order to provide quantitatively more correct ESM results on marine biogeochemistry and carbon cycle feedbacks.

histories↗
DOE OSTI2024

Cyclic Injection Leads to Larger and More Frequent Induced Earthquakes under Volume-Controlled Conditions

As carbon storage technologies advance globally, methods to understand and mitigate induced earthquakes become increasingly important. Although the physical processes that relate increased subsurface pore pressure changes to induced earthquakes have long been known, reliable methods to forecast and control induced seismic sequences remain elusive. Suggested reservoir engineering scenarios for mitigating induced earthquakes typically involve modulation of the injection rate. Some operators have implemented periodic shutdowns (i.e., effective cycling of injection rates) to allow reservoir pressures to equilibrate (e.g., Paradox Valley) or shut-in wells after the occurrence of an event of concern (e.g., Basel, Switzerland). Other proposed scenarios include altering injection rates, actively managing pressures through coproduction of fluids, and preinjection brine extraction. Here, in this work, we use 3D physics-based earthquake simulations to understand the effects of different injection scenarios on induced earthquake rates, maximum event magnitudes, and postinjection seismicity. For comparability, the modeled injection considers the same cumulative volume over the project’s operational life but varies the schedule and rates of fluid injected. Simulation results show that cyclic injection leads to more frequent and larger events than constant injection. Furthermore, with intermittent injection scenario, a significant number of events are shown to occur during pauses in injection, and the seismicity rate remains elevated for longer into the postinjection phase compared to the constant injection scenario.

58 GEOSCIENCES↗
NASA NTRS2003

Numerical Simulations of the Evolution of the CO2 Atmosphere of Mars: 4.53 Ga to the Present

In this paper we approach four key questions that are central to our understanding of the nature of climate change on Mars: 1. By what mechanism has the evolution of the CO2 atmosphere, dominated at early times by violent process, been guided to the moderate quantities present today? 2.How do obliquity cycles affect the migration of CO2 between the various reservoirs of CO2, and what effect does this have on the bulk atmosphere of the planet in time? 3. Is the geophysical evidence that the last few Myr experienced periods that are substantially wetter than the present attributable to a substantially more massive greenhouse atmosphere? 4. Is the current atmospheric pressure determined primarily by the partial pressure of CO2 in cold (approx. 148 K) ice caps, or is it because disequilibrium water, and hence weathering, ceases at pressures below the triple point of water?

C V Manning↗
DOE OSTI2026

Reservoir-scale model of geologic hydrogen production from serpentinization: Cyclic injection in a dual-permeability fracture-matrix system

Geologic hydrogen (GeoH 2 ) from serpentinization is a promising low-carbon resource, but its reservoir-scale behavior remains poorly understood. We develop a dual-permeability reactive transport model for an injector–producer well pair in ultramafic rock that couples multiphase flow, heat transfer, geochemistry, and porosity–permeability evolution. Here, the model is calibrated to olivine flow-through experiments and upscaled to two-year long simulations with continuous injection and cyclic injection with shut-in-to-injection ratios (SIR = 1, 0.5, 0.1). Olivine reacts along high-flux pathways to form lizardite and magnetite, increasing pH; H 2 (aq) and H 2 (g) peak early and then decline as exsolution and advective export outpace local generation. Continuous injection yields the highest cumulative H 2 but the lowest water-use efficiency (1.138 x 10 –6 mol/kgw). Cyclic injection increases this ratio to 1.35 x 10 –6 , 1.377 x 10 –6 , and 1.182 x 10 –6 mol/kgw for SIR = 1, 0.5, and 0.1, respectively; SIR = 0.5 provides a ~20% improvement over continuous injection and the best compromise for pilot design.

08 HYDROGEN↗
NASA NTRS1988

Attitude motion of a non-attitude-controlled cylindrical satellite

In 1985, two non-attitude-controlled satellites were each placed in a low earth orbit by the Scout Launch Vehicle. The satellites were cylindrical in shape and contained reservoirs of hydrazine fuel. Three-axis magnetometer measurements, telemetered in real time, were used to derive the attitude motion of each satellite. Algorithms are generated to deduce possible orientations (and magnitudes) of each vehicle's angular momentum for each telemetry contact. To resolve ambiguities at each contact, a force model was derived to simulate the significant long-term effects of magnetic, gravity gradient, and aerodynamic torques on the angular momentum of the vehicles. The histories of the orientation and magnitude of the angular momentum are illustrated.

Wilkinson, C. K.↗
NASA NTRS1992

Laser-driven hypersonic air-breathing propulsion simulator

A feasibility study is presented of simulating airbreathing propulsion on small scale hypersonic models using laser energy. The laser heat addition scheme allows simultaneous inlet and exhaust flows during wind tunnel testing of models with scramjet models. The proposed propulsion simulation concept has extended the Kantrowitz (1974) idea to propulsive wind tunnel models of hypersonic aircraft. Critical issues in aeropropulsive testing of models based on a ramjet power plant are addressed which include transfer of the correct amount of energy to the flowing gas, efficient absorption of laser energy into the gas, and test performance under tunnel reservoir conditions and at reasonable Reynolds numbers.

Joshi, Prakash B.↗