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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Isotopic fractionation as in-situ sensor of subsurface reactive flow and precursor for rock failure

Greater utilization of subsurface reservoirs perturbs in-situ chemical-mechanical conditions with wide ranging consequences from decreased performance to project failure. Understanding the chemical precursors to rock deformation is critical to reducing the risks of these activities. To address this need, we investigated the coupled flow-dissolution- precipitation-adsorption reactions involving calcite and environmentally-relevant solid phases. Experimentally, we quantified (1) stable isotope fractionation processes for strontium during calcite nucleation and growth, and during reactive fluid flow; (2) consolidation behavior of calcite assemblages in the common brines. Numerically, we quantified water weakening of calcite using molecular dynamics simulations; and quantified the impact of calcite dissolution rate on macroscopic fracturing using finite element models. With microfluidic experiments and modeling, we show the effect of local flow fields on the dissolution kinetics of calcite. Taken together across a wide range of scales and methods, our studies allow us to separate the effects of reaction, flow, and transport, on calcite fracturing and the evolution of strontium isotopic signatures in the reactive fluids.

42 ENGINEERING↗
DOE OSTI2023

Using Data-Driven Prediction of Downstream 1D River Flow to Overcome the Challenges of Hydrologic River Modeling

Methods for downstream river flow prediction can be categorized into physics-based and empirical approaches. Although based on well-studied physical relationships, physics-based models rely on numerous hydrologic variables characteristic of the specific river system that can be costly to acquire. Moreover, simulation is often computationally intensive. Conversely, empirical models require less information about the system being modeled and can capture a system’s interactions based on a smaller set of observed data. This article introduces two empirical methods to predict downstream hydraulic variables based on observed stream data: a linear programming (LP) model, and a convolutional neural network (CNN). We apply both empirical models within the Colorado River system to a site located on the Green River, downstream of the Yampa River confluence and Flaming Gorge Dam, and compare it to the physics-based model Streamflow Synthesis and Reservoir Regulation (SSARR) currently used by federal agencies. Results show that both proposed models significantly outperform the SSARR model. Moreover, the CNN model outperforms the LP model for hourly predictions whereas both perform similarly for daily predictions. Although less accurate than the CNN model at finer temporal resolution, the LP model is ideal for linear water scheduling tools.

13 HYDRO ENERGY↗
DOE OSTI2020

Theoretical power-law relationship between permeability and formation factor

Permeability k and formation factor F are important petrophysical properties that have broad applications to reservoir engineering. Although well-log analysis has advanced reservoir characterizations, direct measurement of k is still not available. Accordingly, resistivity logs and more particularly formation factor have been frequently used to estimate k. In this study, we present a theoretical framework to relate k to F. We use power-law scaling from continuum percolation theory to develop a theoretic relationship between k and F, which is similar in form to a published power-law equation with empirical coefficients. Our result gives a physical interpretation to this empirical relationship. We also compare power-law scaling laws from percolation theory with data derived from experimental measurements and numerical simulations in the literature. Finally, using the Curve Fitting toolbox of MATLAB, the parameters of the power-law equations are determined, and the significance of the non-universal power laws are discussed.

42 ENGINEERING↗
DOE OSTI2022

Survival estimates across five life stages of redfin ( Perca fluviatilis ) exposed to simulated pumped-storage hydropower stressors

The global prevalence of pumped-storage hydropower (PSH) is expected to grow exponentially as countries transition to renewable energy sources. Compared to conventional hydropower, little is currently known regarding PSH impacts on aquatic biota. This study estimated the survival of five life stages (egg, two larval stages, juvenile and adult) of redfin (European) perch (Perca fluviatilis) following passage through a PSH facility during the pumping phase. This was achieved by simulating the individual stressors expected to occur during passage through a 2000-MW PSH facility using laboratory-simulated (shear strain and extreme compression) and modelling (blade strike, BS) approaches. Our results indicate that redfin could survive the shear, pressure and BS stressors expected within the PSH facility, but impacts varied among life stages. Juvenile survival was >70% across all shear strain rates, while the survival of eggs and larvae declined markedly as strain rate increased. All life stages had high survival when exposed to rapid compression and BS. The high survival of redfin to the stressors tested suggests the PSH facility could facilitate the passage of redfin during the pumping phase from the lower to the higher elevation reservoir. This outcome would be welcomed in situations where the species is native, but could have adverse implications for the conservation of native biota where the species is considered a pest.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Are we modeling the properties of unconventional shales correctly?

Technological advancements have led to impressive growth in hydrocarbon (HC) production from unconventional shale reservoirs. However, major limiting factors in shale gas development are i) low recovery of gas in place (<20%), (ii) a rapid decline in well productivity, iii) release of contaminants, despite the use of advanced hydraulic fracturing fluids and multiple fracturing stages. This is mainly due to a lack of understanding of the nanoscale physicochemical properties of shales, especially that of kerogen, the macromolecule that is not only the source of the majority of the HC’s in shales but also holds most of the HC’s in adsorbed state. Over the last few years, a large number of studies have been published on simulating the physicochemical properties of shale utilizing molecular models of kerogen. However, the molecular models of kerogen input into these simulations are based on the kerogen “type” derived from a very limited number of shale samples. In this paper, we examine the variations that can exist in kerogen structure within a particular kerogen “type” across different shale basins and single shale basin at similar thermal maturity levels. We propose that using kerogen “type” based structural models for molecular simulations could lead to inaccurate estimation of HC reserves, HC recovery, HC production, frackability, and quantity and quality of produced water. Our study highlights the need for developing a better classification of kerogen based on its “molecular structure” instead of “type” for a more accurate prediction of physicochemical properties of shales.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2024

Numerical Simulations in Support of a Long-Term Test of Gas Production From Hydrate Accumulations on the Alaska North Slope: Water Production and Associated Design and Management Issues

Here, we investigated numerical simulation strategies for a long-term test of depressurization-induced gas production from the B1 Sand of Unit B at the Hydrate-01 Stratigraphic Test Well. The main objective of this study was to estimate fluid production rates (with emphasis on water production) under a variety of conditions and production scenarios and contribute new insights to the design and management of the field test. In the first part of the study, we investigated the system response to a three-step depressurization process using two limiting sets of flow properties─the expected maximum and minimum intrinsic and effective permeabilities─for the very heterogeneous reservoir. In the second part, we investigated the effect of the production interval length and placement within the formation relative to the boundaries of the hydrate-bearing unit. The best performing well configuration was used in the third part of the study, which used the most representative subsurface flow properties to investigate the effect of the depressurization strategy on the production performance. The best overall performance (largest gas production with modest water production and a strong response at the observation wells) was obtained with a 10 m-long well situated 3 m below the top of the formation and a three-step depressurization scheme at 15-day intervals to a terminal bottomhole pressure of 2.8 MPa. The overall production performance was enhanced by a faster rate of depressurization. Estimated water production rates in all cases were limited and easily manageable. None of the tested well configurations or depressurization strategies significantly reduced water production without also severely reducing gas production. In all the investigated cases, 95% of the long-term fraction of produced water was replenished by inflows from the boundaries and could not be reduced. These substantial water inflows are an unavoidable feature of HU-B and cannot be easily mitigated by a hydraulic control.

02 PETROLEUM↗
DOE OSTI2024

Characterizing Hydraulic Fracture Propagation Before Fracture Hits

Estimating the distance from hydraulic fracture tip to monitor well can be very useful for fracture characterization, well spacing optimization, and preventing parent-child well interference. Heart-shape signal is referred to as the extensional precursor of fracture hit recorded by cross-well strain measurements and can be served as a vital tool to make such estimation. This study incorporates the 3D Displacement Discontinuity Method to understand the impact of fracture geometry and monitor well offset on the heart-shape signal’s characteristics. Results from numerical simulation and analytical solutions reveal a strong linear correlation between the spatial extent of heart shape signal and fracture tip distance. This relationship was further developed to predict tip distance using field data from the Hydraulic Fracture Test Site 2. A reasonable approximation result from field data further validates the methodology. In addition, it is worth noting that the estimation accuracy depends on the ratio between fracture dimension and tip distance. The findings of this study offer a novel approach for real-time monitoring and characterizing hydraulic fracture propagation, which can be further used for well spacing optimization in unconventional and Enhanced Geothermal System reservoir development, as well as cap rock integrity monitoring for carbon sequestration projects.

Jin, Ge↗
DOE OSTI2023

The stability of FeH x and hydrogen transport at Earth’s core mantle boundary

Iron hydride in Earth’s interior can be formed by the reaction between hydrous minerals (water) and iron. Studying iron hydride improves our understanding of hydrogen transportation in Earth’s interior. Our high-pressure experiments found that face-centered cubic (fcc) FeH x (x≤1) is stable up to 165 GPa, and our ab initio molecular dynamics simulations predicted that fcc FeH x transforms to a superionic state under lower mantle conditions. In the superionic state, H-ions in fcc FeH become highly diffusive-like fluids with a high diffusion coefficient of ~3.7 × 10 -4 cm 2 s -1 , which is comparable to that in the liquid Fe-H phase. The densities and melting temperatures of fcc FeH x were systematically calculated. Similar to superionic ice, the extra entropy of diffusive H-ions increases the melting temperature of fcc FeH. Further, the wide stability field of fcc FeH enables hydrogen transport into the outer core to create a potential hydrogen reservoir in Earth’s interior, leaving oxygen-rich patches (ORP) above the core mantle boundary (CMB).

58 GEOSCIENCES↗
DOE OSTI2023

Enabling site-specific well leakage risk estimation during geologic carbon sequestration using a modular deep-learning-based wellbore leakage model

Geologic carbon sequestration (GCS) is a promising technology for mitigating net carbon emissions and growing climate concern by storing CO 2 in reservoirs. Oil and gas brownfields are an attractive option for CO 2 storage, but these sites have many historical wellbores from petroleum production and can be a potential leakage pathway for CO 2 or formation brine. Therefore, risk management of GCS operations requires an assessment of potential well leakage. Due to the high uncertainty of the system, stochastic approaches are ideal for quantifying the range of risk behaviors, but they must be computationally efficient in the face of complex physics. Here, we develop a new physics-centric deep learning wellbore model to predict the leakage of CO 2 and brine through leaky wellbores. Multi-physics numerical simulations were used to generate data sets, and physics-informed features were introduced. Neural networks were optimized with an automated searching algorithm. Feature analysis quantifies the impact of each feature on model prediction and confirms the role of physics-inspired parameters. The model shows high predictive performance across a wide range of geologic and injection conditions and well attributes. In conclusion, a case study illustrates how the model is applied to assess well leakage in GCS operations.

58 GEOSCIENCES↗
DOE OSTI2024

The Role and Lifetime of Dissociative Heterogeneous Processes in Improving Simulated Ozone on Mars

Abstract Ozone simulated in Mars Global Climate Models (MGCMs) is used to assess the underlying chemistry occurring in the atmosphere. Currently, ozone total column abundance (TCA) is under‐predicted in MGCMs by up to 120%, implying missing or inaccurate chemistry in models. Heterogeneous reactions of hydroxyl radicals (HO X ) have been offered as an explanation for some of this bias, because they cause ozone to increase at locations where it's currently under‐predicted. We use four simulations to compare modeled ozone TCA with observations from the UVIS spectrometer aboard the ExoMars Trace Gas Orbiter to improve the representation of heterogeneous processes and their impact on ozone. We use a gas‐phase only run, a dissociative scheme, an adsorbed HO X retention scheme, and a hybrid scheme that combines the dissociative mechanism with the retention of HO X on water ice. We find retention of HO X is dependent on water ice sublimation, and ozone abundance increases when water ice persists for longer periods (1–20 sols). Over time, the loss of HO X causes a depletion in H 2 O 2 concentration (HO X reservoir), and thus allows ozone concentration to increase. When adsorbed HO X are desorbed and dissociate into other by‐products, HO X are not immediately available to destroy ozone. This results in larger ozone concentrations than if desorbed HO X are released directly back into their gaseous states. When using the hybrid scheme, ozone TCA is increased up to 50% where the ozone deficit is greatest, demonstrating the best agreement with observations, and implying that HO X radicals are both retained when adsorbed and dissociate.

Brown, M. A. J.↗
DOE OSTI2024

Heart Shape to Fracture Distance: Characterizing Hydraulic Fracture Propagation before Hits

Estimating the distance from the hydraulic fracture tip to the monitor well can be useful for fracture characterization, well spacing optimization, and preventing parent-child well interference. A heart-shaped signal is referred to as the extensional precursor of a fracture hit recorded by crosswell strain measurements and can serve as a vital tool for such estimation. This study incorporates the 3D displacement discontinuity method (DDM) to understand the impact of fracture geometry and monitor well offset on the heart-shaped signal’s characteristics. Results from numerical simulation and analytical solutions reveal a strong linear correlation between the spatial extent of the heart-shaped signal and the fracture tip distance. This relationship was further developed to predict tip distance using field data from the Hydraulic Fracture Test Site 2 (HFTS2). A reasonable approximation result from field data further validates the methodology. In addition, it is worth noting that the estimation accuracy depends on the ratio between fracture dimension and tip distance. The findings of this study offer a novel approach for real-time monitoring and characterizing hydraulic fracture propagation, which can be further used for well spacing optimization in unconventional and enhanced geothermal system reservoir development, as well as caprock integrity monitoring for carbon sequestration projects.

58 GEOSCIENCES↗
DOE OSTI2022

Observations of phase changes in monoolein during high viscous injection

Serial crystallography of membrane proteins often employs high-viscosity injectors (HVIs) to deliver micrometre-sized crystals to the X-ray beam. Typically, the carrier medium is a lipidic cubic phase (LCP) media, which can also be used to nucleate and grow the crystals. However, despite the fact that the LCP is widely used with HVIs, the potential impact of the injection process on the LCP structure has not been reported and hence is not yet well understood. The self-assembled structure of the LCP can be affected by pressure, dehydration and temperature changes, all of which occur during continuous flow injection. These changes to the LCP structure may in turn impact the results of X-ray diffraction measurements from membrane protein crystals. To investigate the influence of HVIs on the structure of the LCP we conducted a study of the phase changes in monoolein/water and monoolein/buffer mixtures during continuous flow injection, at both atmospheric pressure and under vacuum. The reservoir pressure in the HVI was tracked to determine if there is any correlation with the phase behaviour of the LCP. The results indicated that, even though the reservoir pressure underwent (at times) significant variation, this did not appear to correlate with observed phase changes in the sample stream or correspond to shifts in the LCP lattice parameter. During vacuum injection, there was a three-way coexistence of the gyroid cubic phase, diamond cubic phase and lamellar phase. During injection at atmospheric pressure, the coexistence of a cubic phase and lamellar phase in the monoolein/water mixtures was also observed. The degree to which the lamellar phase is formed was found to be strongly dependent on the co-flowing gas conditions used to stabilize the LCP stream. A combination of laboratory-based optical polarization microscopy and simulation studies was used to investigate these observations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Effects of Charged Solute-Solvent Interaction on Reservoir Temperature during Subsurface CO2 Injection

A short-term side-effect of CO2 injection is a developing low-pH front that forms ahead of the bulk water injectant, due to differences in solute diffusivity. Observations of downhole well temperature show a reduction in aqueous-phase temperature with the arrival of a low-pH front, followed by a gradual rise in temperature upon the arrival of a high concentration of bicarbonate ion. In this work, we model aqueous-phase transient heat advection and diffusion, with the volumetric energy generation rate computed from solute-solvent interaction using the Helgeson–Kirkham–Flowers (HKF) model, which is based on the Born Solvation model, for computing specific molar heat capacity and the enthalpy of charged electrolytes. A computed injectant water temperature profile is shown to agree with the actual bottom hole sampled temperature acquired from sensors. The modeling of aqueous-phase temperature during subsurface injection simulation is important for the accurate modeling of mineral dissolution and precipitation because forward dissolution rates are governed by a temperature-dependent Arrhenius model.

58 GEOSCIENCES↗
DOE OSTI2021

Pericentric passage-driven star formation in satellite galaxies and their hosts: CLUES from local group simulations

Local Group satellite galaxies show a wide diversity of star formation histories (SFHs) whose origin is yet to be fully understood. Using hydrodynamical simulations from the Constrained Local UniversE project, we study the SFHs of satellites of Milky Way-like galaxies in a cosmological context: while in the majority of the cases the accretion on to their host galaxy causes the satellites to lose their gas, with a subsequent suppression in SF, in about 25 per cent of our sample we observe a clear enhancement of SF after infall. Peaks in SF clearly correlate with the satellite pericentric passage around its host and, in one case, with a satellite–satellite interaction. We identify two key ingredients that result in enhanced SF after infall: galaxies must enter the host’s virial radius with a reservoir of cold gas M gas , inf/M vir , inf ≳ 10 -2 and with a minimum pericentric distance ≳10 kpc (mean distance ~50 kpc for the full sample), in order to form new stars due to compression of cold gas at pericentric passage. On the other hand, satellites that infall with little gas or whose pericentric distance is too small, have their gas ram-pressure stripped and subsequent SF quenched. The pericentric passage of satellites likewise correlates with SF peaks in their hosts, suggesting that this mechanism induces bursts of SF in satellites and central galaxies alike, in agreement with recent studies of our Galaxy’s SFH. Our research findings can explain the recently reported multiple stellar populations observed in dwarf galaxies such as Carina and Fornax, and should be taken into account in semi-analytic models of galaxy formation and satellite quenching.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Thermal modeling of three lakes within the continuous permafrost zone in Alaska using the LAKE 2.0 model

Abstract. Lakes in the Arctic are important reservoirs of heat with much lower albedo in summer and greater absorption of solar radiation than surrounding tundra vegetation. In the winter, lakes that do not freeze to their bed have a mean annual bed temperature >0 ∘C in an otherwise frozen landscape. Under climate warming scenarios, we expect Arctic lakes to accelerate thawing of underlying permafrost due to warming water temperatures in the summer and winter. Previous studies of Arctic lakes have focused on ice cover and thickness, the ice decay process, catchment hydrology, lake water balance, and eddy covariance measurements, but little work has been done in the Arctic to model lake heat balance. We applied the LAKE 2.0 model to simulate water temperatures in three Arctic lakes in northern Alaska over several years and tested the sensitivity of the model to several perturbations of input meteorological variables (precipitation, shortwave radiation, and air temperature) and several model parameters (water vertical resolution, sediment vertical resolution, depth of soil column, and temporal resolution). The LAKE 2.0 model is a one-dimensional model that explicitly solves vertical profiles of water state variables on a grid. We used a combination of meteorological data from local and remote weather stations, as well as data derived from remote sensing, to drive the model. We validated modeled water temperatures with data of observed lake water temperatures at several depths over several years for each lake. Our validation of the LAKE 2.0 model is a necessary step toward modeling changes in Arctic lake ice regimes, lake heat balance, and thermal interactions with permafrost. The sensitivity analysis shows us that lake water temperature is not highly sensitive to small changes in air temperature or precipitation, while changes in shortwave radiation and large changes in precipitation produced larger effects. Snow depth and lake ice strongly affect water temperatures during the frozen season, which dominates the annual thermal regime of Arctic lakes. These findings suggest that reductions in lake ice thickness and duration could lead to more heat storage by lakes and enhanced permafrost degradation.

Clark, Jason A. (ORCID:0000000248105261)↗
DOE OSTI2020

Transient flow modeling in fractured media using graphs

In this work, we describe a method to simulate transient fluid flows in fractured media using an approach based on graph theory. Our approach builds on past work where the graph-based approach was successfully used to simulate steady-state fluid flows in fractured media. We find a mean computational speedup of the order of 1400 from an ensemble of a 100 discrete fracture networks in contrast to the O(10 4 ) speedup that was obtained for steady-state flows earlier. However, the transient flows considered here involve an additional degree of complexity that was not present in the steady-state flows considered previously with a graph-based approach, that of time marching and solution of the flow equations within a time-stepping scheme. We verify our method with an analytical test case and demonstrate its use on a practical problem related to fluid flows in hydraulically fractured reservoirs. By enabling the study of transient flows, we create an opportunity for a wide set of possibilities where a steady-state approximation is not sufficient, such as the example motivated by hydraulic fracturing that we present here. This work validates the concept that graphs are able to reliably capture the topological properties of the fracture network and serve as effective surrogates in an uncertainty-quantification framework.

58 GEOSCIENCES↗
DOE OSTI2023

Quantifying the Contributions of Aerosol- and Snow-Produced ClNO 2 through Observations and 1D Modeling

Nitryl chloride (ClNO 2 ) is a radical reservoir that forms and accumulates in the nocturnal atmospheric boundary layer influenced by combustion emissions and chloride (e.g., sea salt and road salt). Upon sunrise, ClNO 2 rapidly photolyzes to generate highly reactive chlorine radicals (Cl • ) that affect the air quality by generating secondary air pollutants. Recent studies have shown road salt aerosols and saline snowpack to be sources of ClNO 2 in the wintertime urban environment; however, the quantitative contributions of each chloride source are not known. In this study, we examine the vertically resolved contributions of aerosol particles and saline snowpack as sources of ClNO 2 by using an observationally constrained snow–atmosphere coupled one-dimensional model applied to wintertime Kalamazoo, Michigan, U.S. Model simulations show that ClNO 2 emitted from urban snowpack can be vertically transported throughout the entire atmospheric boundary layer and can be a significant source of ClNO 2 , contributing up to ~60% of the ClNO 2 budget near the surface. Here, modeled snowpack ClNO 2 emission rates were 6 (±7) times higher than the observationally derived emission rates, suggesting that not all snow chloride is available for reaction. ClNO 2 production from both aerosol particles and snow emissions are required to best simulate the observed surface-level ClNO 2 . Using the traditional bulk parameterization for ClNO 2 produced from particles significantly overestimated ClNO 2 due to the assumption of having equivalent dinitrogen pentoxide (N 2 O 5 ) uptake and chloride availability for the entire particle population. In comparison, the chemically resolved surface area-based parameterization slightly underestimated the observations and had uncertainties deriving from ClNO 2 production from residential wood burning particles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Pressure stability in explicitly coupled simulations of poromechanics with application to CO 2 sequestration

We study in detail the pressure stabilizing effects of the non-iterated fixed-stress splitting in poromechanical problems which are nearly undrained and incompressible. When applied in conjunction with a spatial discretization which does not satisfy the discrete inf–sup condition, namely a mixed piecewise linear–piecewise constant spatial discretization, the explicit fixed-stress scheme can have a pressure stabilizing effect in transient problems. This effect disappears, however, upon time step refinement or the attainment of steady state. The interpretation of the scheme as an Augmented Lagrangian method similar to Uzawa iteration for incompressible flow helps explain these results. Moreover, due to the slowly evolving solution within undrained seal regions, we show that the explicit fixed-stress scheme requires very large time steps to reveal its pressure stabilizing effect in examples of geologic CO 2 sequestration. We note that large time steps can result in large errors in drained regions, such as the aquifer or reservoir regions of these examples, and can prevent convergence of nonlinear solvers in the case of multiphase flows, which can make the explicit scheme an unreliable source of pressure stabilization. We conclude by demonstrating that pressure jump stabilization is as effective in the explicit fixed-stress setting as in the fully implicit setting for undrained problems, while maintaining the stability and convergence of the fixed-stress split for drained problems.

58 GEOSCIENCES↗