Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 595 records · Page 33

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.↗

Investigation of enhanced CO 2 storage in deep saline aquifers by WAG and brine extraction in the Minnelusa sandstone, Wyoming

Geological CO 2 sequestration in deep saline aquifers has been extendedly investigated to reach the goal of carbon neutral as large amount of CO 2 can be reduced in a short time. Here, this study investigates the feasibility of water-alternating-gas (WAG) injection and brine extraction on enhancing CO 2 storage in a target deep saline aquifer, Minnelusa Sandstone in the Powder River Basin (PRB) of Wyoming. An integrated numerical model is developed to account for four scenarios: (1) CO 2 continuous injection (CI) through one injection well, (2) WAG injection through one injection well, (3) CO 2 CI through one injection well and brine extraction through one producer, (4) WAG injection through one injection well and brine extraction through one producer. Simulation results suggest that WAG injection and brine extraction corporately make it possible to enhance CO 2 injectivity while securing CO 2 storage safety. For instance, WAG injection considerably reduces structural trapping contribution while enhances dissolved and residual trapping contributions. As for brine extraction, it can decrease the maximum averaged reservoir pressure by 37% and 22% for CI and WAG injection, respectively. Besides, sensitive analyses of the operational parameters for the fourth scenarios are performed. Results reveal that, when the amount of total CO 2 injection has been predetermined, CO 2 injection time and rate within each period of WAG should be small whereas it is preferable to have large water injection time and rate. This is to secure desirable dissolved and residual trapping contributions (store CO 2 with safety). Sensitive degree results confirm that CO 2 injection rate, water injection time, and CO 2 injection time of WAG injection exhibit significant effects on CO 2 storage whereas the impacts of producer bottom-hole pressure (BHP), water injection rate, and CO 2 -water injection time are relatively weaker. This study not only sheds light on achieving double-win goal: enhance CO 2 injectivity and store CO 2 with safety, but also can be a critical reference for other CO 2 storage in deep saline aquifers.

58 GEOSCIENCES↗

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↗

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↗

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.↗

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↗

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↗

Water and Climate Impacts on ERCOT Long-Term Systems Assessment

Water and climate change pose many potential challenges to the electric power system. Substantial water is withdrawn every day to support thermoelectric power generating unit operations, and changes to water supply have the potential to affect generation dispatch. Climate change can accelerate growing demand for electricity, which can necessitate additional generating capacity, often in locations with limited water supply. Drought conditions also threaten thermoelectric power plant operations due to streamflow and reservoir levels dropping below intake structures, or water temperatures exceeding a power plants' permitted operating conditions. Here we explore how future climate change might influence decisions related to electricity capacity expansion planning in Electric Reliability Council of Texas (ERCOT) using a multi-model framework. Specifically, water resource modeling is used to simulate climate impacts on the future water supply for thermoelectric and hydropower generation for four future climate projections. Separately, temperature impacts on electricity load are evaluated for these scenarios. These climate impacts are applied to five alternative electricity futures in an electricity capacity expansion model that projects future generation and transmission capacity additions in ERCOT. Results indicate that climate has a measurable influence on future generation and transmission capacity needs, with temperature-driven increases in peak and average load resulting in 5-15 GW additional generating capacity and up to 1 GW additional transmission capacity. Additional capacity is a diverse mix of PV, natural gas, and wind, depending on the makeup of economic, policy and technology assumptions. Climate impacts increase total system costs 2-5%, while the marginal cost of energy and emissions are not affected substantially by climate change effects.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Observations and Modeling of Fiber-Optics Strain on Hydraulic Fracture Height Growth in HFTS-2

Understanding fracture height growth can be of great significance to optimizing field development and improving recovery. The Hydraulic Fracturing Test Site 2 (HFTS-2) has provided a unique opportunity and an advanced dataset to allow us to observe and understand fracture geometries rigorously. Low frequency distributed acoustic sensing (LF-DAS) data from a vertical well in HFTS-2 showed three key observations: (i) excessive upward height growth (>1000 ft) and limited downward growth of the hydraulic fractures during pumping, (ii) considerable additional upward fracture height growth (~300 ft) after well shut in, and (iii) very complex LF-DAS strain rate patterns for a small fiber-to-stage offset. Advanced geomechanical modeling was performed to simulate the hydraulic fracture propagation and the resulting strain responses in the vertical direction. The modeling results demonstrated asymmetric upward and downward fracture height growths as observed in HFTS-2 with a similar upward height growth rate. Simulated waterfall plots of vertical strain rate showed distinct patterns for different fiber-to-fracture distances. The upward-growing fracture tip can be clearly identified by the interfaces between compressing and extending zones. It was also found that the complex strain rate patterns observed in HFTS-2 for small fiber-stage offsets were not caused by the mechanical layering but possibly result from the simultaneous propagation of multiple hydraulic fractures at different rates. Furthermore, the simulation results improved the understating of the HFTS-2 LF-DAS data and the simulated strain rate patterns could also serve as templates for fracture height interpretation from LF-DAS data in future.

58 GEOSCIENCES↗

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↗

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↗

CO 2 geological storage: Critical insights on plume dynamics and storage efficiency during long-term injection and post-injection periods

In the quest for net-zero greenhouse gas emission, CO 2 capture and storage technologies play a prominent role. Estimates suggest a potential of sequestering up to several thousand metric gigatons of CO 2 . However, some of the main challenges in its development are an assessment of reservoir-specific storage capacities due to variabilities and heterogeneities in the underlying properties (structure, formations, facies, petrophysics, reservoir architecture, and others) and understanding the migration of CO 2 in the subsurface. How much CO 2 can be stored in a reservoir, and how quickly can it be injected? What will be the dominant storing (trapping) mechanisms? These questions, among others, remain as challenges to the commercial deployment of this technology. This article critically investigates CO 2 plume characteristics and determines the evolving contributions of different trapping mechanisms during a 100-year injection period. Here, simulation results on the CO 2 plume dynamics indicate that lateral propagation is much larger than vertical propagation. The ratio of average lateral to vertical dispersion ranged between 6.3 and 18.8 for the cases investigated. On average, for all the scenarios investigated, the free phase (mobile supercritical CO 2 ) is the dominant trapping mechanism (in terms of the amount of CO 2 stored) during the injection period, accounting up to 45% of the CO 2 injected. Both residual (30%) and solubility (25%) trapping storage ratios decline with time during injection, and no significant mineral trapping appeared during the 100-year injection period. However, the behavior of CaCO 3 suggests a predominant dissolution of CaCO 3 , especially in the vicinity of the injection well.

03 NATURAL GAS↗

Operation of Natural Gas Pipeline Networks With Storage Under Transient Flow Conditions

Here, we formulate a nonlinear optimal control problem for intraday operation of a natural gas pipeline network that includes storage reservoirs. The dynamics of compressible gas flow through pipes, compressors, reservoirs, and wells are considered. In particular, a reservoir is modeled as a rigid, hollow container that stores gas under isothermal conditions and uniform density, and a well is modeled as a vertical pipe. For each pipe, flow dynamics are described by a coupled partial differential equation (PDE) system in density and mass flux variables, with momentum dissipation modeled using the Darcy–Wiesbach friction approximation. Compressors are modeled as scaling up the pressure of gas between the inlet and outlet. The governing equations for all network components are spatially discretized and assembled into a nonlinear differential-algebraic equation (DAE) system, which synthesizes above-ground pipeline and subsurface reservoir dynamics into a single reduced-order model. We seek to maximize an objective function that quantifies economic profit and network efficiency subject to the flow equations and inequalities that represent operating limitations. The problem is solved using a primal–dual interior point solver, and the solutions are validated in computational experiments and simulations on several pipeline test networks to demonstrate the effectiveness of the proposed methodology.

03 NATURAL GAS↗

Utah FORGE 5-2428: Fracture Permeability Impact on Reservoir Stress and Seismic Slip Behavior - Workshop Presentation

This is a presentation on the Fracture Permeability Impact on Seismic Slip Behavior project by Lawrence Livermore National Laboratory, presented by Dr. Kayla A. Kroll. The project's objective is to develop, apply and validate a holistic thermal, hydrologic, mechanical, and chemical (THMC) workflow that includes evaluation of induced seismic slip in EGS reservoirs. This presentation was featured in the Utah FORGE R&D Annual Workshop on September 8, 2023. The workshop provided a valuable opportunity to explore the progress made in each of the 17 Research and Development projects funded under Solicitation 2020-1 which aim to enhance our understanding of the crucial factors influencing the development of Enhanced Geothermal Systems (EGS) reservoirs and resources.

15 GEOTHERMAL ENERGY↗

Utah FORGE 5-2557: Fluid and Temperature in Fracture Mechanics and Coupled THMC Processes - Workshop Presentation

This is a presentation on the Role of Fluid and Temperature in Fracture Mechanics and Coupled Thermo-Hydro-Mechanical-Chemical (THMC) Processes for Enhanced Geothermal Systems project by Purdue University, presented by Distinguished Professor of Physics & Astronomy, Laura J. Pyrak-Nolte. The project's objective was to develop and validate a macroscopic model that accounts for local deformation/frictional behavior, seismic/aseismic behavior, chemical reactions, and determine the adequacy of classic Coulomb failure vs. rate-and-state friction. This presentation was featured in the Utah FORGE R&D Annual Workshop on September 8, 2023. The workshop provided a valuable opportunity to explore the progress made in each of the 17 Research and Development projects funded under Solicitation 2020-1 which aim to enhance our understanding of the crucial factors influencing the development of Enhanced Geothermal Systems (EGS) reservoirs and resources.

15 GEOTHERMAL ENERGY↗

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.↗

Reactive chemical transport simulations of geologic carbon sequestration: Methods and applications

Chemical reaction simulations are considerably used to quantitatively assess the long-term geologic carbon sequestration (GCS), such as CO 2 sequestration capacity estimations, leakage pathway analyses, enhanced oil recovery (EOR) efficiency studies, and risk assessments of sealing formations (caprocks), wellbores, and overlying underground water resources. All these require a deep understanding of the CO 2 - associated chemical reactions. To ensure long-term, safe CO 2 sequestration in the intended formations, modeling is the only way to plausibly assess the CO 2 flow, reaction, and transport over thousands of years. This review summarizes the multiple methodologies for describing homogeneous and heterogeneous chemical reaction patterns and multiscale application examples, the recent progress and current status of chemical reaction simulations for GCS, and the impact of such simulations on geological CO2 sequestration performance. Technical gaps and future challenges are also discussed for further study. The trends and challenges of such studies include: (1) the combination of coupled chemical, mechanical, and transport processes with calibrated experiments and associated uncertainty/risk assessments; (2) enhancement of the ability to simulate detailed geophysical and geochemical equations to mimic in situ conditions; and (3) characterization of multiscale subsurface systems with detailed conceptual models and assignment of suitable boundary conditions for field-scale sequestration fields. One major gap remaining is the current lack of accurate (and scale-justified) kinetic and equilibrium chemical reaction parameters under reservoir conditions. Advanced models that couple chemical, mechanical, and transport processes with scale-justified parameters, from lab to field-scale experiments, are required for quantitative assessments of sequestration capacity and the long-term safety of GCS projects.

58 GEOSCIENCES↗

Radio AGN in nearby dwarf galaxies: the important role of AGN in dwarf galaxy evolution

ABSTRACT We combine deep optical and radio data, from the Hyper Suprime-Cam and the Low-Frequency Array (LOFAR), respectively, to study 78 radio active galactic nuclei (AGN) in nearby (z < 0.5) dwarf galaxies. Comparison to a control sample, matched in stellar mass and redshift, indicates that the AGN and controls reside in similar environments, show similar star formation rates (which trace gas availability) and exhibit a comparable incidence of tidal features (which indicate recent interactions). We explore the AGN properties by combining the predicted gas conditions in dwarfs from a cosmological hydrodynamical simulation with a Monte Carlo suite of simulated radio sources, based on a semi-analytical model for radio-galaxy evolution. In the subset of LOFAR-detectable simulated sources, which have a similar distribution of radio luminosities as our observed AGN, the median jet powers, ages, and accretion rates are ∼1035 W, ∼5 Myr, and ∼10−3.4 M⊙ yr−1, respectively. The median mechanical energy output of these sources is ∼100 times larger than the median binding energy expected in dwarf gas reservoirs, making AGN feedback plausible. Since special circumstances (in terms of environment, gas availability, and interactions) are not necessary for the presence of AGN, and the central gas masses are predicted to be an order of magnitude larger than that required to fuel the AGN, AGN triggering in dwarfs is likely to be stochastic and a common phenomenon. Together with the plausibility of energetic feedback, this suggests that AGN could be important drivers of dwarf galaxy evolution, as is the case in massive galaxies.

Davis, F. (ORCID:0000000285224980)↗