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

A coupled thermo-hydro-mechanical model for simulating leakoff-dominated hydraulic fracturing with application to geologic carbon storage

A potential risk of injecting CO2 into storage reservoirs with marginal permeability (≲ 10 mD (1 mD = 10 –15 m 2 )) is that commercial injection rates could induce fracturing of the reservoir and/or the caprock. Such fracturing is essentially fluid-driven fracturing in the leakoff-dominated regime. Recent studies suggested that fracturing, if contained within the lower portion of the caprock complex, could substantially improve the injectivity without compromising the overall seal integrity. Modeling this phenomenon entails complex coupled interactions among the fluids, the fracture, the reservoir, and the caprock. Here, we develop a simple method to capture all these interplays in high fidelity by sequentially coupling a hydraulic fracturing module with a coupled thermal-hydrological-mechanical (THM) model for nonisothermal multiphase flow. The model was made numerically tractable by taking advantage of self-stabilizing features of leakoff-dominated fracturing. The model is validated against the PKN solution in the leakoff-dominated regime. Moreover, we employ the model to study thermo-poromechanical responses of a fluid-driven fracture in a field-scale carbon storage reservoir that is loosely based on the In Salah project's Krechba reservoir. The model reveals complex yet intriguing behaviors of the reservoir-caprock-fluid system with fracturing induced by cold CO 2 injection. We also study the effects of the in situ stress contrast between the reservoir and caprock and thermal contraction on the vertical containment of the fracture. The proposed model proves effective in simulating practical problems on length and time scales relevant to geological carbon storage.

58 GEOSCIENCES↗
DOE OSTI2020

Assessing mechanical response of CO 2 storage into a depleted carbonate reef using a site-scale geomechanical model calibrated with field tests and InSAR monitoring data

Geomechanical risks of injection have raised concerns regarding secure CO 2 storage. In this work, a combined monitoring and modeling approach is used to assess the stress changes and surface uplift associated with CO 2 injection into a depleted carbonate reef of the Michigan basin. A site-scale geomechanical model is built by assigning mechanical properties of formations using well-log and experimental data. Gravity load is applied to the model to estimate the vertical component of stress as well as different lateral boundary displacement scenarios to estimate horizontal stresses. We used a poroelastic pressure-dependent model (instead of a linear elastic mechanical earth model) to calibrate initial stresses using hydraulic fracture test data measured at depleted reservoir status. Multi-phase fluid flow-geomechanical simulations are performed to estimate the poroelastic response during (1) primary depletion (2) field-scale CO 2 injection phase (3) a hypothetical forecast scenario in which well bottom hole pressure (BHP) reach 45000 KPa. The predicted surface uplift is less than 1 mm at the end of the field-scale CO 2 injection phase which is in good agreement with Interferometric Synthetic Aperture Radar (InSAR) uplift measurement. Although the InSAR data shows an insignificant uplift, hydromechanical modeling of injection shows that CO 2 injection still causes reservoir deformation emphasizing the role of carbonate overburden and reservoir formation mechanical properties and limited size of reef on diminishing the surface deformation. Modeling indicates poroelastic response of caprock matters to estimate uplift. The lower permeability of the top two layers provides additional barrier to large uplift. Also, history of subsidence due to production should be accounted to predict uplift due to a follow up injection correctly. This report shows the significance of combining a calibrated geomechanical model with field measured stresses and monitoring data to be used as a tool to ensure the safety of CO 2 storage.

42 ENGINEERING↗
DOE OSTI2024

Dynamic Binary Complexes (DBC) as Super-Adjustable Viscosity Modifiers for Hydraulic Fracturing Fluids

In the preceding project year two, we refined three DBC formulations from a selection of over 50 different chemistries. The optimization study primarily encompassed testing for (i) reversibility extent, (ii) performance in the presence of chemical additives, (iii) adhesion and friction behavior during displacement in wellbores and pipelines, (iv) corrosion protection performance, and (v) injection performance with model fracture systems at the laboratory scale. Highly promising results obtained from all these tests signify the significant potential of DBCs in enhancing hydrocarbon recovery from unconventional reservoirs. The primary activities in the third project year included publishing experimental findings across multiple articles and conducting outreach initiatives. Throughout the year, we undertook tasks such as replicating experimental results, further optimizing various formulations and their associated experimental sets, and conducting additional tests to address missing components based on reviewer feedback and suggestions. We also explored the surfactant and friction-reduction aspects of selected formulations through drag reduction tests. In addition, we constructed an improved fracturing performance setup and performed flow injection tests. The specific DBC formulations focused on during this project period were A8/B1, A12/B5, and A10/B12. We also obtained results for additional DBC formulations and a select few commercial fracturing fluids for the purpose of comparison. Within the project's scope, we aim to enhance the experimental findings with the development of various models. The first two years focused on two key aspects: (i) the creation of a high-fidelity hydraulic fracturing model for non-Newtonian fluids to gain insights into the implementation of DBC fluids in fracking environments, and (ii) the development of a multiphase flow simulator for estimating total production, fluid saturation in the reservoir, and the creation of a fracture propagation model and kinetic Monte Carlo (kMC) models for diverse applications. In the third year, we delved into the fundamental nanostructural properties of DBCs, exploring aspects such as material chemistry, pH tunability, and control of DBC formation and stability. Subsequently, in the extension year, we conducted a systematic investigation of various building blocks containing primary, secondary, and tertiary amine functional groups to understand their impact on rheological and viscoelastic properties. Furthermore, we explored a Dissipative Particle Dynamics (DPD) model to simulate self-assembly processes with precision, creating a high-fidelity representation of relevant nanostructures. The tasks performed this year with the significant results obtained have been discussed in Section 2.

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

Tuning the Legacy Survey of Space and Time (LSST) Observing Strategy for Solar System Science

Abstract The Vera C. Rubin Observatory is expected to start the Legacy Survey of Space and Time (LSST) in early to mid-2025. This multiband wide-field synoptic survey will transform our view of the solar system, with the discovery and monitoring of over five million small bodies. The final survey strategy chosen for LSST has direct implications on the discoverability and characterization of solar system minor planets and passing interstellar objects. Creating an inventory of the solar system is one of the four main LSST science drivers. The LSST observing cadence is a complex optimization problem that must balance the priorities and needs of all the key LSST science areas. To design the best LSST survey strategy, a series of operation simulations using the Rubin Observatory scheduler have been generated to explore the various options for tuning observing parameters and prioritizations. We explore the impact of the various simulated LSST observing strategies on studying the solar system’s small body reservoirs. We examine what are the best observing scenarios and review what are the important considerations for maximizing LSST solar system science. In general, most of the LSST cadence simulations produce ±5% or less variations in our chosen key metrics, but a subset of the simulations significantly hinder science returns with much larger losses in the discovery and light-curve metrics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Explaining Neural Spike Activity for Simulated Bio-plausible Network through Deep Sequence Learning

With significant improvements in large-scale simulations of brain models, there is a growing need to develop tools for rapid analysis and interpreting the simulation results. In this work, we explore the potential of sequential deep learning models to understand and explain the network dynamics among the neurons extracted from a large-scale neural simulation in STACS (Simulation Tool for Asynchronous Cortical Stream). Our method employs a representative neuroscience model that abstracts the cortical dynamics with a reservoir of randomly connected spiking neurons with a low stable spike firing rate throughout the simulation duration. We subsequently analyze the spike dynamics of the simulated spiking neural network through an autoencoder model and an attention-based mechanism.

Kulkarni, Shruti↗
DOE OSTI2021

Thermal-Hydrological-Mechanical Modelling of Stockton University Reservoir Cooling System, Fine Scale Stress Test Modelling

Mesh, properties, initial conditions, injection/withdrawal rates for modelling thermal, hydrological, and mechanical effects of fluid injection to and withdrawal from ground for Stockton University reservoir cooling system (aquifer storage cooling system), Galloway, New Jersey, for unscheduled two hour injection at 133 % designed capacity, on fine scale grid, with some results. Second simulation of J.T. Smith, E. Sonnenthal, P. Dobson, P. Nico, and M. Worthington, 2021. Thermal-hydrological-mechanical modeling of Stockton University reservoir cooling system, Proceedings of the 46th Workshop on Geothermal Reservoir Engineering, Stanford University, SGP-TR-218, from which Figures 6-9, pertain.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Multi‐scale simulation of rock compaction through breakage models with microstructure evolution

Abstract. Regional subsidence due to fluid depletion includes the interaction among multiple physical processes. Specifically, rock compaction is governed by coupled hydro-mechanical feedbacks involving fluid flow, effective stress change and pore collapse. Although poroelastic models are often used to explain the delay between depletion and subsidence, recent evidence indicates that inelastic effects could alter the rock microstructure, thus exacerbating coupling effects. Here, a constitutive law built within the framework of Breakage Mechanics is proposed to account for the inherent connection between rock microstructure, hydraulic conductivity, and pore compaction. Furthermore, it is embedded into a 1-D hydromechanical coupled finite element analysis (FEA) to explore the effects of micro-structure rearrangement on the development of reservoir compaction. Numerical examples with the proposed model are compared with simulations under constant hydraulic conductivity to illustrate the model capability to capture the non-linear processes of reservoir compaction induced by fluid depletion.

Buscarnera, Giuseppe↗
NASA NTRS2016

Considerations in Managing the Fill Rate of the Grand Ethiopian Renaissance Dam Reservoir Using a System Dynamics Approach.

The purpose of this study is to evaluate simulated fill rate scenarios for the Grand Ethiopian Renaissance Dam while taking into account plausible climate change outcomes for the Nile River Basin. The region lacks a comprehensive equitable water resource management strategy, which creates regional security concerns and future possible conflicts. We employ climate estimates from 33 general circulation models within a system dynamics model as a step in moving toward a feasible regional water resource management strategy. We find that annual reservoir fill rates of 8-15% are capable of building hydroelectric capacity in Ethiopia while concurrently ensuring a minimum level of stream flow disruption into Egypt before 2039. Insofar as climate change estimates suggest a modest average increase in stream flow into the Aswan, climate changes through 2039 are unlikely to affect the fill rate policies. However, larger fill rates will have a more detrimental effect on stream flow into the Aswan, particularly beyond a policy of 15%. While this study demonstrates that a technical solution for reservoir fill rates is feasible, the corresponding policy challenge is political. Implementation of water resource management strategies in the Nile River Basin specifically and Africa generally will necessitate a national and regional willingness to cooperate.

system dynamics↗
DOE OSTI2025

Characterization of kerogen nanopores using 2D NMR relaxation and MD simulations

The characterization of kerogen nanopores is crucial for predicting the geostorage capacity and recoverability of natural gas in unconventional gas shale reservoirs. Towards this end, a powerful technique is presented which integrates 2D NMR T 1 -T 2 relaxation measurements with molecular dynamics (MD) simulations of hydrocarbons confined in the nanopores of kerogen. The integrated NMR-MD technique is demonstrated using T 1 -T 2 measurements of kerogen isolates and organic-rich chalks saturated with heptane, together with MD simulations of heptane completely dissolved in a realistic kerogen model. The NMR-MD results are used to extract the swelling ratio and nanopore size distribution of kerogen as a function of depth in the reservoir. The effects of organic nanoconfinement on the T 1 relaxation dispersion and T 2 residual dipolar coupling of heptane are investigated, as well as the effect of downhole effective stress on the kerogen nanopore size as a function of depth and compaction. Potential applications in partially depleted gas shale reservoirs are discussed, including CO 2 utilization/geostorage, geostorage of green H 2 , and integration of the NMR-MD technique with thermodynamic models for predicting the competitive sorption of gas mixtures in kerogen.

Compaction↗
DOE OSTI2026

Chemically Enabled CO 2 -Enhanced Oil Recovery in Multi-Porosity, Hydrothermally Altered Carbonates in the Southern Michigan Basin (Final Technical Report)

This is the Final Technical Report for the project "Chemically Enabled CO 2 -Enhanced Oil Recovery in Multi-Porosity, Hydrothermally Altered Carbonates in the Southern Michigan Basin." Over the course of six years of collaboration between Battelle and project partners, all stated objectives of the program have been completed, including full geological characterization of the TBR trend (See companion report for Task2), laboratory and modeling experiments to determine the optimum composition and design of CO 2 -EOR operations in the TBR trend, execution of a field test of chemically-enhanced CO 2 in a TBR well, and integration of the data and learnings gathered during these efforts into a full-trend development plan. Detailed reporting on these activities, their outcomes, and implications for trend-wide development is provided in the report. This report and encompassed data will provide TBR field operators with detailed information on what worked, what did not work, and how to proceed with production optimization of their TBR assets using chemically-enhanced CO 2 -EOR. CO 2 -EOR is a relatively well understood and broadly implemented strategy for increasing incremental production across the oil and gas industry, but its application has been primarily focused on reservoirs with limited heterogeneity. The intention of this project was show first that the same physical mechanisms that improve recovery factors in homogeneous reservoirs (namely wettability alteration, viscosity alteration, oil swelling, and mobility control) are at play in heterogeneous reservoirs. This was proven by the project’s laboratory studies and dynamic simulations, with the potential exception of mobility control, which needs further study. The second intention was to demonstrate via direct field testing that CO 2 -EOR can work in a strongly heterogeneous reservoir. While the field test strategy implemented during this project did not succeed in producing oil, data gathered during the test sheds light on what may work for field operators who try chemically-enhanced CO 2 -EOR within their own reservoirs, significantly reducing the level of uncertainty carried by first-of-a-kind commercial efforts that could (and should) follow this test. Simultaneously, the project has identified several large-volume ethanol plants and other sources of CO 2 emissions in the region and provided a handrail that CO 2 emitters and field operators can leverage to capture, transport, and inject that CO 2 into their fields. This project has also shown that, in many cases, the economics of CO 2 -EOR in the TBR are attractive. And finally, by completing a project of this scope in the southern Michigan Basin, the project has contributed to the knowledge base and operational experience of field operators, state regulatory agencies, local service companies, and state universities, with CO 2 -EOR projects which should allow follow-on projects to proceed safely and efficiently.

02 PETROLEUM↗
DOE OSTI2020

Gas Production from Hot Water Circulation through Hydraulic Fractures in Methane Hydrate-Bearing Sediments: THC-Coupled Simulation of Production Mechanisms

Methane hydrates, widely found in permafrost and deep marine sediments, have great potential as a future energy 6 source. Conventional production schemes perform poorly for challenging hydrate reservoirs with low permeability. We propose an 7 efficient production scheme by combining hydraulic fracturing from horizontal wells and hot water circulation through fractures. A 8 fully coupled thermo-hydro-chemical (THC) model is developed to simulate the key physical processes during gas production from 9 a hydrate reservoir representative of typical geological settings in Shenhu, South China Sea. We found that the gas production 10 process has two distinct stages divided by thermal breakthrough: a relatively short prebreakthrough stage and a postbreakthrough 11 stage yielding stable gas production. Heat advection along and near the hydraulic fracture dominates the prebreakthrough stage, 12 whereas conduction-driven thermal recovery in the volume around fractures dominates the postbreakthrough stage. We identified 13 that the steady-state injection temperature has a strong effect on the performance of the proposed scheme while the fluid mass 14 circulation rate has a moderate impact beyond a threshold. The proposed scheme proves to be efficient and robust over a range of 15 reservoir conditions with respect to initial hydrate saturation and intrinsic permeability, including their spatial heterogeneities, 16 thereby offering a promising solution for challenging reservoir conditions.

gas hydrate, fracture flow, gas production, numeri↗
NASA NTRS1982

Snowmelt runoff modeling in simulation and forecasting modes with the Martinec-Mango model

The Martinec-Rango snowmelt runoff model was applied to two watersheds in the Rio Grande basin, Colorado-the South Fork Rio Grande, a drainage encompassing 216 sq mi without reservoirs or diversions and the Rio Grande above Del Norte, a drainage encompassing 1,320 sq mi without major reservoirs. The model was successfully applied to both watersheds when run in a simulation mode for the period 1973-79. This period included both high and low runoff seasons. Central to the adaptation of the model to run in a forecast mode was the need to develop a technique to forecast the shape of the snow cover depletion curves between satellite data points. Four separate approaches were investigated-simple linear estimation, multiple regression, parabolic exponential, and type curve. Only the parabolic exponential and type curve methods were run on the South Fork and Rio Grande watersheds for the 1980 runoff season using satellite snow cover updates when available. Although reasonable forecasts were obtained in certain situations, neither method seemed ready for truly operational forecasts, possibly due to a large amount of estimated climatic data for one or two primary base stations during the 1980 season.

Shafer, B.↗
DOE OSTI2023

Projecting Future Colorado River Basin Water and Hydropower Operations

Drought, variable renewable energy (VRE) resources (primarily wind and photovoltaic technologies), and power markets can all impact hydropower plant operations. We analyze how these three factors collectively may potentially change the future operation and economics of federal dams and hydropower plants located in the Colorado River Basin (CRB). Federal CRB energy and capacity resources, including 4,200 MW of hydropower capacity, are marketed by Western Area Power Administration (WAPA) and sold through long-term Firm Electric Service (FES) contracts to entities located in 15 western and central States that serve millions of electricity consumers. The parameters of FES contracts allow a wide range of scheduling flexibility for FES customers. The FES customers will adjust the timing of energy delivery according to the changing water, power market, and VRE contribution level to the power grid. We simulate and project water and hydropower operations with two models: 1) reservoir operation and 2) power system operation (production cost) with various temporal granularity. We consider two sets of grid scenarios (one representing today's infrastructure, and one representing higher VRE contributions in 2036) and two sets of water scenarios (drought and normal hydrology). We modeled scenarios to understand the power market-energy price levels and patterns, FES customer scheduling, hydropower economics, WAPA finances and hydropower operation for various hydrologic conditions, and VRE contribution level to the power grid. The detailed modeling framework includes: (1) creating ensembles of CRB hydrologic and hydropower futures, (2) selecting representative futures, (3) simulation of FES customer contracts, (4) simulating western U.S. power grid economic scheduling and dispatch under two VRE contribution levels, utilizing information from modeling steps (1)-(4), (5) assessment of reservoir water and hydropower operation, and (6) comparing scenario operations and economics. Our study team comprising WAPA, and an interdisciplinary team of national labs will gain insight into how WAPA can prepare for and adapt its practices in response to an evolving U.S. power grid.

Colorado River↗
DOE OSTI2026

Basin-scale analysis of Mokelumne River Formation for multi-well CO 2 injection

Large-scale carbon sequestration will likely require multiple projects injecting CO 2 into the same subsurface formation, raising concerns about safe operation and efficient use of storage capacity. This study evaluates the long-term response of the Mokelumne River Formation in California’s Sacramento Basin to multi-megaton CO 2 injection using three geologic models of the formation and the open-source simulator GEOS. The analysis focuses on three aspects of reservoir performance: (1) average pressure increase and dissolved CO 2 mass after 30 years for varying well counts and injection rates, (2) pressure interference in a multi-well configuration, and (3) dynamic storage capacity with identification of overpressure-prone regions. The results show that average formation pressure increases linearly with injected mass, while CO 2 dissolution exhibits mixed scaling: approximately linear with the number of wells but sublinear with injection rate, indicating that distributing injection across more wells enhances dissolution more effectively than increasing per-well rates. Pressure-interference effects are significant, with lower-permeability conditions delaying their onset but amplifying their magnitude at later times. Dynamic capacity, defined by the first occurrence of pressure exceeding the local overburden-based limit anywhere in the formation, varies across geologic models and assumed overburden pressure gradients. A lower fidelity geologic model predicts nearly twice the storage capacity of the two higher fidelity models, which consistently estimate approximately 1 Gt under the upper-bound overburden pressure gradient considered for the Sacramento Basin. In all model scenarios, overpressure develops away from injection wells, particularly in higher-elevation regions, highlighting the importance of basin-scale modelling for identifying risks beyond the immediate well vicinity.

Basin-scale↗
DOE OSTI2025

Three-Dimensional Pore Networks in Miocene Stevens Sandstone of California: Implications for CO 2 Geologic Storage

The Miocene Stevens Sandstone in the San Joaquin Basin of California is increasingly recognized as a promising candidate for CO 2 geological storage due to the enormous storage capacity, proven sealing, and existing infrastructure. In this study, computed microtomography imaging and pore network modeling were employed to investigate the influence of pore geometry and wettability on the CO 2 injectivity and residual trapping. Image analysis revealed that a significant fraction of the cement and matrix consists of microporous regions. The microporosity can substantially increase the overall pore space, yet its contribution to permeability remains modest, particularly in samples with low permeability. The intrinsic heterogeneity of turbidite reservoirs further complicates the reservoir properties among different layers. Two-phase flow simulations under varying wettability conditions (water-wet, weak water-wet, and neutral-wet) demonstrated that the CO 2 injection is predominantly controlled by macropores. CO 2 invades microporous regions only after these larger pores are filled. The presence of microporosity leads to a decrease in both initial and residual CO 2 saturations, with the magnitude of the reduction being influenced by wettability. Neutral-wet scenarios exhibit higher CO 2 mobility and thus lower residual trapping than water-wet scenarios. The results imply that heterogeneity in pore geometry and cement distribution across different layers can result in stratified CO 2 flow pathways, complicating efforts to predict injection performance. Overall, the Stevens Sandstone shows considerable promise for CO 2 geologic storage, but effective implementation will require detailed characterization of the pore structure as well as the integration of reactive fluid flow to account for potential mineral dissolution and fines migration.

fluids↗
DOE OSTI2023

TReactMech v4.217

TReactMech couples geomechanical processes (poroelasticity, failure, and inelastic strain) with multiphase nonisothermal flow (derived from TOUGH2) and reactive geochemical transport. At its core is the reactive-transport code TOUGHREACT v4.13. TReactMech is an efficient hybrid parallel simulator, solving the geomechanics using finite elements and MPI/PetSc, the multiphase flow using integrated finite difference and MPI/PETSc, and the reactive chemistry using OpenMP. The advantages of TReactMech are in its multiphase flow capabilities (e.g., supercritical CO2, supercritical water, air) and parallel geomechanics including full 3-D stress tensor, shear and tensile failure, coupled to porosity and permeability changes. It is backwardly compatible with TOUGH2 and TOUGHREACT v4.13, allowing for easier transitions between the codes. TReactMech can be used to simulate many natural and engineered subsurface systems, including geothermal reservoirs, borehole heat exchangers, geologic carbon sequestration, geologic storage of nuclear waste, groundwater resources, weathering, sediment diagenesis, seafloor hydrothermal circulation, hydrofracturing in unconventional reservoirs, and injection/production-induced surface deformation.

Sonnenthal, Eric↗
NASA NTRS1984

An investigation of the phototropic effect on seedling orientation in a microgravity environment: A student involvement project

A microgravity growth chamber was designed to investigate the phototropic response of radish seedlings. Enclosed in a one fourth inch thick, hexagonal, fiberglass-foam spacepak nineteen inches across corners, the experiment consists of a growth chamber and germination tray, a water reservoir and solenoid valve, a fluorescent light for photo simulation, a Minolta X700 camera with programmable back, a 50 mm macro lens and flash, a battery pack, and a computer controller. Two temperature sensors and one light sensor located in the walls of the growth chamber provide temperature and illumination data. A computer provides 8 K command and 34 K data storage capability. The experiment was not activated during the STS flight because a malfunctioning latching relay stuck and reduced the battery power level.

Barainca, J. W.↗