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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Optimal CO 2 storage management considering safety constraints in multi-stakeholder multi-site GCS projects: A Markov game perspective

Geological carbon storage (GCS) projects could involve a diverse array of stakeholders or players from public, private, and regulatory sectors, each with different objectives and responsibilities. Given the complexity, scale, and long-term nature of GCS operations, determining whether individual stakeholders can independently optimize their interests — or whether collaborative coalition agreements are needed — remains a central question for effective GCS project planning and management. To access large, high-quality storage resources, future GCS deployment may increasingly occur in geologically connected sites, where shared geological features such as pressure space and reservoir pore capacity can lead to competitive behavior among stakeholders. In this work, we propose a paradigm based on Markov games to quantitatively investigate how different coalition structures affect the goals of stakeholders. We frame this multi-stakeholder multi-site problem as a multi-agent reinforcement learning problem with safety constraints. Our approach enables agents to learn optimal strategies while complying with safety regulations. We present an example where multiple operators are injecting CO 2 into their respective project areas in a geologically connected basin. To address the high computational cost of repeated simulations of high fidelity models, a previously developed surrogate model based on the Embed-to-Control (E2C) framework is employed. Our results demonstrate the effectiveness of the proposed framework in addressing optimal management of CO 2 storage when multiple stakeholders with different objectives and goals are involved.

58 GEOSCIENCES↗
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↗
NASA NTRS1978

Performance of a thermionic converter module utilizing emitter and collector heat pipes

A thermionic converter module simulating a configuration for an out-of-core thermionic nuclear reactor was designed, fabricated, and tested. The module consists of three cylindrical thermionic converters. The tungsten emitter of the converter is heated by a tungsten, lithium heat pipe. The emitter heat pipes are immersed in a furnace, insulated by MULTI-FOIL thermal insulation, and heated by tungsten radiation filaments. The performance of each thermionic converter was characterized before assembly into the module. Dynamic voltage, current curves were taken using a 60 Hz sweep and computerized data acquisition over a range of emitter, collector, and cesium-reservoir temperatures. An output power of 215 W was observed at an emitter temperature of 1750 K and a collector temperature of 855 K for a two diode module. With a three diode module, an output power of 270 W was observed at an average emitter temperature of 1800 K and a Collector temperature of 875 K.

Kroeger, E. W.↗
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↗
NASA NTRS2016

Production of Peroxy Nitrates in Boreal Biomass Burning Plumes over Canada During the BORTAS Campaign

The observations collected during the BOReal forest fires on Tropospheric oxidants over the Atlantic using Aircraft and Satellites (BORTAS) campaign in summer 2011 over Canada are analysed to study the impact of forest fire emissions on the formation of ozone (O3 and total peroxy nitrates (sigma)PNs, (sigma)ROONO2. The suite of measurements on board the BAe-146 aircraft, deployed in this campaign, allows us to calculate the production of O3 and of (sigma)PNs, a long-lived NOx reservoir whose concentration is supposed to be impacted by biomass burning emissions.In fire plumes, profiles of carbon monoxide (CO), which is a well-established tracer of pyrogenic emission, show concentration enhancements that are in strong correspondence with a significant increase of concentrations of (sigma)PNs, where as minimal increase of the concentrations of O3 and NO2 is observed. The (sigma)PN and O3 productions have been calculated using the rate constants of the first- and second-order react Pions of volatile organic compound (VOC) oxidation. The (sigma)PN and O3 productions have also been quantified by 0-D model simulation based on the Master Chemical Mechanism. Both methods show that in fire plumes the average production of (sigma)PNs and O3 are greater than in the background plumes, but the increase of (sigma)PN production is more pronounced than the O3 production. The average (sigma)PN production in fire plumes is from 7 to 12 times greater than in the background, whereas the average O3 production in fire plumes is from 2 to 5 times greater than in the background. These results suggest that, at least for boreal forest fires and for the measurements recorded during the BORTAS campaign,fire emissions impact both the oxidized NOy and O3;but (1)(sigma)PN production is amplified significantly more thanO3 production and (2) in the forest fire plumes the ratio between the O3 production and the (sigma)PN production is lower than the ratio evaluated in the background air masses, thus confirming that the role played by the (sigma)PNs produced during biomass burning is significant in the O3 budget. The implication of these observations is that fire emissions in some cases, for example boreal forest fires and in the conditions reported here, may influence more long-lived precursors of O3than short-lived pollutants, which in turn can be transported and eventually diluted in a wide area.

biomass burning↗
NASA NTRS2019

Integration of Airborne and Ground Observations of Nitryl Chloride in the Seoul Metropolitan Area and the Implications on Regional Oxidation Capacity During KORUS-AQ 2016

Nitryl chloride (ClNO2) is a radical reservoir species that releases chlorine radicals upon photolysis. An integrated analysis of the impact of ClNO2on regional photochemistry in the Seoul metropolitan area (SMA) during the Korea–United States Air Quality Study (KORUS-AQ) 2016 field campaign is presented. Comprehensive multiplatform observations were conducted aboard the NASA DC-8 and at two ground sites (Olympic Park, OP; Taehwa Research Forest, TRF), representing an urbanized area and a forested suburban region, respectively. Positive correlations between daytime Cl2 and ClNO2 were observed at both sites, the slope of which was dependent on O3 levels. The possible mechanisms are explored through box model simulations constrained with observations. The overall diurnal variations in ClNO2 at both sites appeared similar but the night-time variations were systematically different. For about half of the observation days at the OP site the level of ClNO2 increased at sunset but rapidly decreased at around midnight. On the other hand, high levels were observed throughout the night at the TRF site. Significant levels of ClNO2 were observed at both sites for 4–5 h after sunrise. Airborne observations, box model calculations, and back-trajectory analysis consistently show that these high levels of ClNO2 in the morning are likely from vertical or horizontal transport of air masses from the west. Box model results show that chlorine-radical-initiated chemistry can impact the regional photochemistry by elevating net chemical production rates of ozone by ∼25 % in the morning.

Nitryl chloride (ClNO2)↗
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↗
NASA NTRS2022

Mapping Firn Saturation over Greenland using NASA’s Soil Moisture Active Passive Satellite

Mapping the spatial extent of recently identified englacial hydrological features (i.e., ice slabs and perennial firn aquifers) formed by meters-thick water-saturated firn layers over the percolation facies of the Greenland Ice Sheet using L-band microwave radiometry has recently been demonstrated. However, these initial maps are binary, and do not provide a parameter to estimate the spatial variability in the thickness and volumetric fraction of meltwater stored within the firn pore space. Here, we exploit enhanced-resolution vertical-polarization L-band brightness temperature (T^(B)_(V)) imagery (2015–2019) generated using observations collected over Greenland by NASA’s Soil Moisture Active Passive (SMAP) satellite and a simple two-layer L-band brightness temperature model. We map water-saturated firn layers via a “firn saturation” parameter, and interpret our results together with ice slab and perennial firn aquifer spatial extents, estimates of snow accumulation simulated via the Regional Atmospheric Climate Model (RACMOp2.3), and airborne radar surveys collected via NASA’s Operation IceBridge (OIB) campaigns. We find that variable firn saturation parameter values are mapped in lower snow accumulation ice slab areas in western, northern, and northeastern Greenland, where firn is colder and water-saturated firn layers seasonally refreeze as solid-ice. Higher firn saturation parameter values are mapped in higher snow accumulation perennial firn aquifer areas in southeastern, southern, and northwestern Greenland, where firn is near the melting point, and meters-thick water-saturated firn layers exist. Our results have implications for identifying expansive englacial reservoirs that store significant volumes of meltwater in locations that are vulnerable to meltwater-induced hydrofracturing and accelerated outlet glacier flow year-round.

Greenland ice sheet↗
NASA NTRS1992

Relationship of peroxyacetyl nitrate to active and total odd nitrogen at northern high latitudes: influence of reservoir species on NOx and O3

Measurements of peroxyacetyl nitrate (PAN), NO, NO2, HNO3, NOy (total odd nitrogen), and O3 were made in the high-latitude troposphere over North America and Greenland (35 degrees to 82 degrees N) during the Arctic Boundary Layer Expedition (ABLE 3A) (July-August 1988) throughout 0-to 6-km altitudes. These data are analyzed to quantitatively describe the relationships between various odd nitrogen species and assess their significance to global tropospheric chemistry. In the free troposphere, PAN was as much as 25 times more abundant than NOx. PAN to NOx ratio increased with increasing altitude and latitude. PAN was found to be the single most abundant reactive nitrogen species in the free troposphere and constituted a major fraction of NOy, PAN to NOy ratios were about 0.1 in the boundary layer and increased to 0.4 in the free troposphere. A 2-D global photochemical model with C1-C3 hydrocarbon chemistry is used to compare model predictions with measured results. A sizable portion (approximately 50%) of the gaseous reactive nitrogen budget is unaccounted for, and unknown organic nitrates and pernitrates are expected to be present. Model calculations (August 1, 70 degrees N) show that a major fraction of the observed NOx (50 to 70% of median) may find its source in the available PAN reservoir. PAN and the unknown reservoir species may have the potential to control virtually the entire NOx availability of the high latitude troposphere. It is predicted that the summer NOx and O3 mixing ratios in the Arctic/sub-Arctic troposphere would be considerably lower in the absence of the ubiquitous PAN reservoir. Conversely, this PAN reservoir may be responsible for the observed temporal increase in tropospheric O3 at high latitudes.

NASA Center ARC↗
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↗
NASA NTRS2009

Using Volcanic Ash to Remove Dissolved Uranium and Lead

Experiments have shown that significant fractions of uranium, lead, and possibly other toxic and/or radioactive substances can be removed from an aqueous solution by simply exposing the solution, at ambient temperature, to a treatment medium that includes weathered volcanic ash from Pu'u Nene, which is a cinder cone on the Island of Hawaii. Heretofore, this specific volcanic ash has been used for an entirely different purpose: simulating the spectral properties of Martian soil. The treatment medium can consist of the volcanic ash alone or in combination with chitosan, which is a natural polymer that can be produced from seafood waste or easily extracted from fungi, some bacteria, and some algae. The medium is harmless to plants and animals and, because of the abundance and natural origin of its ingredient( s), is inexpensive. The medium can be used in a variety of ways and settings: it can be incorporated into water-filtration systems; placed in contact or mixed with water-containing solids (e.g., soils and sludges); immersed in bodies of water (e.g., reservoirs, lakes, rivers, or wells); or placed in and around nuclear power plants, mines, and farm fields.

McKay, David S.↗
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↗