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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2024

On the Feasibility of Deep Geothermal Wells Using Numerical Reservoir Simulation

This study examines the geothermal energy extraction potential from the basement rock within the Denver–Julesburg Basin, focusing on the flow performance and heat extraction efficiency of different geothermal well configurations. It specifically compares U-shaped, V-shaped, inclined V-shaped, and pipe-in-pipe configurations against enhanced geothermal system setups. Through numerical modeling, we evaluated the thermal behavior of these systems under various operational scenarios and fracture conditions. The results suggest that while closed-loop systems offer moderate temperature increases, Enhanced geothermal system configurations show substantial potential for high-temperature extraction. This underscores the importance of evaluating well configurations in complex geological settings. The insights from this study aid in strategic geothermal energy planning and development, marking significant advancements in geothermal technology and setting a foundation for future explorations and optimizations.

15 GEOTHERMAL ENERGY↗

DOE OSTI2022

Oil Recovery Prediction for Polymer Flood Field Test of Heavy Oil on Alaska North Slope Via Machine Assisted Reservoir Simulation

The first ever polymer flood field pilot to enhance the recovery of heavy oils on the Alaska North Slope is ongoing. This study constructs and calibrates a reservoir simulation model to predict the oil recovery performance of the pilot through machine-assisted reservoir simulation techniques. To replicate the early water breakthrough observed during waterflooding, transmissibility contrasts are introduced into the simulation model, forcing viscous fingering effects. In the ensuing polymer flood, these transmissibility contrasts are reduced to replicate the restoration of injection conformance during polymer flooding, as indicated by a significant decrease in water cut. Later, transmissibility contrasts are reinstated to replicate a water surge event observed in one of the producing wells during polymer flooding. This event may represent decreased injection conformance from fracture overextension; its anticipated occurrence in the other production well is included in the final forecast. The definition of polymer retention in the simulator incorporates the tailing effect reported in laboratory studies; this tailing effect is useful to the simultaneous history match of producing water cut and produced polymer concentration. The top 24 best-matched simulation models produced at each stage of the history matching process are used to forecast oil recovery. The final forecast clearly demonstrates that polymer flooding significantly increases the heavy oil production for this field pilot compared to waterflooding alone. This exercise displays that a simulation model is only valid for prediction if flow behavior in the reservoir remains consistent with that observed during the history matched period. Critically, this means that a simulation model calibrated for waterflooding may not fully capture the benefits of an enhanced oil recovery process such as polymer flooding. Therefore, caution is recommended in using basic waterflood simulation models to scope potential enhanced oil recovery projects.

Keith, Cody Douglas↗

DOE OSTI2021

A Graphics Processing Unit–Based, Industrial Grade Compositional Reservoir Simulator

Summary Recently, graphics processing units (GPUs) have been demonstrated to provide a significant performance benefit for black-oil reservoir simulation, as well as flash calculations that serve an important role in compositional simulation. A comprehensive approach to compositional simulation based on GPUs has yet to emerge, and the question remains as to whether the benefits observed in black-oil simulation persist with a more complex fluid description. We present a positive answer to this question through the extension of a commercial GPU-based black-oil simulator to include a compositional description based on standard cubic equations of state (EOSs). We describe the motivations for the selected nonlinear formulation, including the choice of primary variables and iteration scheme, and support for both fully implicit methods (FIMs) and adaptive implicit methods (AIMs). We then present performance results on an example sector model and simplified synthetic case designed to allow a detailed examination of runtime and memory scaling with respect to the number of hydrocarbon components and model size, as well as the number of processors. We finally show results from two complex asset models (synthetic and real) and examine performance scaling with respect to GPU generation, demonstrating that performance correlates strongly with GPU memory bandwidth. NOTE: This paper is also published as part of the 2021 SPE Reservoir Simulation Conference Special Issue.

Engineering↗

DOE OSTI2025

Geothermal Reservoir Simulation Analysis in Support of Electricity Co-Production Feasibility Study at the Blackburn Oil Field, Nevada

Geothermal electricity co-production is a viable option for oil reservoirs producing large water cuts with elevated wellhead temperatures. Repurposing existing oil wells significantly reduces initial investment costs historically associated with geothermal resource utilization. The National Renewable Energy Laboratory (NREL), partnering with Gradient Geothermal, Inc. (formerly known as Transitional Energy) and Grant Canyon Oil and Gas, has been tasked to evaluate the feasibility of geothermal electricity co-production at the Blackburn Oil Field with Organic Rankine Cycle (ORC) generators. The Devonian steady-state reservoir has historically been producing high water cuts of 240 degrees F (115.6 degrees C) observed at the wellhead without documented pressure drawdown or thermal breakthrough. An estimated initial reservoir temperature of approx. 260 degrees F (126.7 degrees C) has been observed in the field and history-matched in a wellbore production analysis and reservoir simulation. Our objective was to develop a conceptual geological model of the subsurface, simulate a natural-state reservoir, model production scenarios, and complete a technical feasibility analysis to accomplish this task. Through extensive modeling and the use of available proprietary and public data, it was possible simulate three scenarios that indicated minimal thermal decline over the duration of a simulated ten-year production and re-injection scheme.

Blackburn Nevada↗

DOE OSTI2024

Geothermal Reservoir Simulation Analysis in Support of Electricity Co-Production Feasibility Study at the Blackburn Oil Field, Nevada: Preprint

Geothermal electricity co-production is a viable option for oil reservoirs producing large water cuts with elevated wellhead-observed temperatures. Repurposing existing oil wells significantly reduces initial investment costs historically associated with geothermal resource utilization. The National Renewable Energy Laboratory (NREL), partnering with Gradient Geothermal, Inc. (formerly known as Transitional Energy) and Grant Canyon Oil and Gas, has been tasked to evaluate the feasibility of geothermal electricity co-production at the Blackburn Oil Field with Organic Rankine Cycle (ORC) generators. The Devonian steady-state reservoir has historically been producing high water cuts of 240 degrees F (115.6 degrees C) observed at the wellhead without documented pressure drawdown or thermal breakthrough. An estimated initial reservoir temperature of approx. 260degreesF (126.7 degrees C) has been observed in the field and history-matched in a wellbore production analysis and reservoir simulation. Our objective was to develop a conceptual geological model of the subsurface, simulate a natural-state reservoir, model production scenarios, and complete a technical feasibility analysis to accomplish this task. Through extensive modeling and the use of available proprietary and public data, it was possible simulate three scenarios that indicated minimal thermal decline over the duration of a simulated ten-year production and re-injection scheme.

Blackburn Nevada↗

DOE OSTI2022

Investigating Performance of Geothermal Reservoirs With Bi-Directional Heat Transfer for Seasonal Thermal Storage Using Thermal-Hydraulic Reservoir Simulations

This paper presents initial simulation results of deep subsurface seasonal thermal storage of waste heat from a natural gas-fired power plant during the summer months for providing heating during winter months at the Cornell University main campus in Ithaca, New York. An injection temperature profile is calculated with the MEnU tool and numerical reservoir simulations are conducted with the simulator FALCON. Both a porous medium and large single fracture geothermal reservoir are considered at about 3 km depth with initial reservoir temperature of 75 degrees C. Simulations indicate that a much smaller reservoir - which is technically easier and less risky to develop - is sufficient when operating the reservoir for seasonal thermal storage instead of for continuous heat production. For example, storing 13 GWh of heat seasonally in a 100 m thick porous medium based on the Galway formation requires about 65 m of well spacing to limit the yearly variation in production temperature to 10 degrees C. Operating this reservoir for continuous heat production would require a well spacing of 365 m to limit the thermal decline to 10 degrees C after 10 years. Neglecting wellbore heat losses and assuming uniform and homogeneous properties, a round-trip thermal efficiency of 100% can be obtained if the yearly averaged injection temperature matches the initial reservoir temperature.

district heating↗

DOE OSTI2021

Multilevel Spectral Coarsening for Graph Laplacian Problems with Application to Reservoir Simulation

We extend previously developed two-level coarsening procedures for graph Laplacian problems written in a mixed saddle point form to the fully recursive multilevel case. The resulting hierarchy of discretizations gives rise to a hierarchy of upscaled models, in the sense that they provide approximation in the natural norms (in the mixed setting). This property enables us to utilize them in three applications: (i) as an accurate reduced model, (ii) as a tool in multilevel Monte Carlo simulations (in application to finite volume discretizations), and (iii) for providing a sequence of nonlinear operators in a full approximation scheme for solving nonlinear pressure equations discretized by the conservative two-point flux approximation. Finally, we illustrate the potential of the proposed multilevel technique in all three applications on a number of popular benchmark problems used in reservoir simulation.

multilevel Monte Carlo↗

DOE OSTI2021

A large-scale comparison of Artificial Intelligence and Data Mining (AI&DM) techniques in simulating reservoir releases over the Upper Colorado Region

In recent years, the Artificial Intelligence and Data Mining (AI&DM) models have become popular tools in assisting various aspects of reservoir operation. However, the practical uses are still rarely reported. Comparison experiment of many AI&DM models over a large number of reservoir cases is particularly valuable to help reservoir operators first examine the usefulness and transferability of different AI&DM models, and then identify the most stable and reliable AI&DM model in assist of various decision-making processes. In this study, a total of 12 AI&DM models with different parameterizations and simulation scenarios are comprehensively tested out and compared in simulating the controlled reservoir outflows of 33 reservoir cases over the Upper Colorado Region, United States. Results show that the Random Forecast and the Long-Short-Term-Memory model could consistently derive the best statistical performance than other models under the baseline simulation scenario. The employed AI&DM models could obtain satisfactory statistical interquartile ranges (25–75%) between [0.6–0.9], [0.3–0.8], and [0.2–0.8], for CORR, NSE, and KGE measurements, respectively, and [1.5–6.5], [–15 to 20], and [0.5–8.5] for the normalized RMSE, PBIAS and RSR measurements, respectively. Results also show Multi-Layer Perceptron model and Extreme Gradient Boosting Tree Algorithm produced more stable and superior performance than other models under more complex input scenarios. We also found that the performance of different AI&DM models are closely relevant to the reservoir elevations, sizes, and functionalities. Discussions were made about the sensitivity of AI&DM models’ parameterizations and the key advantages of AI&DM models over the rule-based reservoir models. We further identify that the main advantage of AI&DM models is the flexibility in designing input structures, whereas the rule-based simulation model is rather limited. Future studies were suggested regarding the best way reservoir operators and researchers could use, select, and apply different AI&DM models in simulating reservoir releases under different natural and modeling environments. Finally, this comparison study also serves as a reference and a piece of groundwork for further promoting the practical uses of AI&DM models in assisting reservoir operation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Validation of strongly coupled geomechanics and gas hydrate reservoir simulation with multiscale laboratory tests

In this work, we validate a coupled flow-geomechanics simulator for gas hydrate deposits, named T+M AM , performing two meter-scale laboratory experiments of gas hydrates for production by depressurization, replicating the gas hydrate deposit in the Ulleung Basin, East Sea, South Korea. The first experiment with a sand-only specimen is a 1D 1 m-scale depressurization test based on the excess gas method, which represents the grain coating hydrate growth. On the other hand, the second is a 3D 1.5 m-scale test with the excess water method for a sand-mud alternating layer system, representing the pore filling hydrate growth. We measure production and displacement at the top with different depressurization levels. In particular, the 3D test exhibits high coupling strength of substantial deformation induced by incompressibility of water and high deformability of the specimen. For validation, we match pressure, flow rate, and displacement between the experimental data and numerical results. Thus, we identify that T+M AM is a reliable simulator, which can be applied to fields in both permafrost and deep oceanic hydrate deposits of strongly coupled flow and geomechanics systems. This validation also implies that other coupled simulators based on the same coupling formulation as T+M AM can be validated when individual flow and geomechanics simulators are stable and reliable.

02 PETROLEUM↗

DOE OSTI2024

Reservoir Simulations and Accuracy Measures for CO 2 and Pressure Plumes for the Broom Creek Formation in Central North Dakota

Conference paper presented at 17th International Conference on Greenhouse Gas Control Technologies (GHGT-17), Calgary, Alberta, Canada, October 20–24, 2024. This study estimated the accuracy of the predicted 5-year carbon dioxide (CO 2 ) and pressure plumes for a hypothetical geological storage project injecting CO 2 into the Broom Creek Formation near central North Dakota.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2025

Deep Learning-based Surrogate Model for Efficient Reservoir Simulation in Large-scale Geological Carbon Storage: Application in IBDP Dataset

This project introduces an advanced deep learning (DL)-based surrogate modeling approach to enhance the efficiency and accuracy of large-scale geological carbon storage (GCS) simulations. Using the Illinois Basin Decatur Project (IBDP) dataset as training data, the study employs a residual U-Net architecture to predict critical state variables such as pressure and CO₂ saturation, as well as CO₂ plume migration. By incorporating key geological parameters (e.g., porosity, permeability, and rock facies) and physics-informed inputs like the diffusive time of flight and time step, the DL model effectively reduces computational complexity while maintaining robust physical constraints. Compared to traditional simulators like Eclipse, the DL model achieves remarkable accuracy, with a root mean square error (RMSE) of 1.57 psi for pressure and 0.007 for saturation, and dramatically reduces computational time from hours to just 69.9 seconds for 50-step simulations. These results demonstrate the potential of innovative DL methodologies to improve the predictivity and operational efficiency of GCS simulations, providing a reliable foundation for decision-making in CCS operations. Supported by the SMART initiative, this project underscores the success of leveraging computational innovations to advance CCS technologies.

advanced deep learning↗

DOE OSTI2025

Influence of Reservoir Convection on Heat Extraction with Closed-Loop Geothermal Systems

Over the past few years, advancements in closed-loop geothermal systems (CLGS), also called advanced geothermal systems (AGS), have sparked a renewed interest in these types of designs. CLGS have certain advantages over traditional and enhanced geothermal systems (EGS), including not requiring in-situ reservoir permeability, conservation of the circulating fluid, and allowing for different fluids, including working fluids directly driving a turbine at the surface. CLGS may be attractive in environments where water resources are limited, rock contaminants must be avoided, and stimulation treatments are not available (e.g., due to regulatory or technical reasons). Despite these advantages, CLGS have some challenges, including limited surface area for heat transfer and requiring long wellbores and laterals to obtain multi-MW output in conduction-only reservoirs. CLGS have been investigated in conduction-only systems. In this paper, we explore the impact of both forced and natural convection on the levels of heat extraction with a CLGS deployed in a hot wet rock reservoir. We bound potential benefits of convection by investigating liquid reservoirs over a range of natural and forced convective coefficients. Additionally, we investigate the effects of permeability, porosity, and geothermal temperature gradient in the reservoir on CLGS outputs. Reservoir simulations indicate that reservoir permeabilities of at least ~100 mD are required for natural convection to increase the heat output with respect to a conduction-only scenario. The impact increases with increasing reservoir temperature. When subject to a forced convection flow field, Darcy velocities of at least 10-7 m/s are required to obtain an increase in heat output.

advanced geothermal systems↗

DOE OSTI2025

Seismic H2: Version 1.0

Seismic-H2 is an integrated software package for geological hydrogen reservoir simulation, optimization, and leakage monitoring. The package includes multiple components: (1) code used for modeling seismic wave propagation in 3D heterogeneous elastic media based on finite-difference method to support detection of geological hydrogen storage reservoir leakage; (2) 3D reservoir simulations of leaks from an underground reservoir and 3D simulations of saline aquifers and depleted gas reservoirs; (3) seismic monitoring costs of passive and active seismic monitoring required for UHS; (4) rock physics calculations and interpolations for converting the reservoir simulations from part (2) into the elastic media models in part (1); (5) pre-processing seismic data; and lastly (6), a GUI interface that combines these different components.

Creasy, Neala↗

DOE OSTI2021

History Matching and Performance Prediction of a Polymer Flood Pilot in Heavy Oil Reservoir on Alaska North Slope

The first-ever polymer flood pilot to enhance heavy oil recovery on Alaska North Slope (ANS) is ongoing. After more than 2.5 years of polymer injection, significant benefit has been observed from the decrease in water cut from 65% to less than 15% in the project producers. The primary objective of this study is to develop a robust history-matched reservoir simulation model capable of predicting future polymer flood performance. In this work, the reservoir simulation model has been developed based on the geological model and available reservoir and fluid data. In particular, four high transmissibility strips were introduced to connect the injector-producer well pairs, simulating short-circuiting flow behavior that can be explained by viscous fingering and reproducing the water cut history. The strip transmissibilities were manually tuned to improve the history matching results during the waterflooding and polymer flooding periods, respectively. It has been found that higher strip transmissibilities match the sharp water cut increase very well in the waterflooding period. Then the strip transmissibilities need to be reduced with time to match the significant water cut reduction. The viscous fingering effect in the reservoir during waterflooding and the restoration of injection conformance during polymer flooding have been effectively represented. Based on the validated simulation model, numerical simulation tests have been conducted to investigate the oil recovery performance under different development strategies, with consideration for sensitivity to polymer parameter uncertainties. The oil recovery factor with polymer flooding can reach about 39% in 30 years, twice as much as forecasted with continued waterflooding. Besides, the updated reservoir model has been successfully employed to forecast polymer utilization, a valuable parameter to evaluate the pilot test’s economic efficiency. All the investigated development strategies indicate polymer utilization lower than 3.5 lbs/bbl in 30 years, which is economically attractive.

Wang, Xindan↗

DOE OSTI2023

Techno-Economic Analysis for a Potential Geothermal District Heating System in Tuttle, Oklahoma: Preprint

Geothermal deep direct use (DDU) has potential across a wide swath of the United States but is underutilized due to challenging project economics associated with developing a deep geothermal resource for a large-scale and variable heat demand. The National Renewable Energy Laboratory (NREL) and University of Oklahoma evaluated the feasibility of a geothermal district heating (GDH) and cooling system in two schools and 250 houses by utilizing existing oil and gas (O&G) wells in Tuttle, Oklahoma. Heating and cooling demand in the two schools and a typical single-family house were modeled using EnergyPlus building energy simulation software. The modeling results indicated that annual heating demand in two schools and 250 houses is approximately 2.61 GWhth, and cooling demand in the two schools is approximately 2.65 GWhth. In this scope, the techno-economic analysis (TEA) was conducted using the GEOPHIRES tool combined with the TOUGH2 reservoir simulator. The reservoir performance, including geothermal heat production capacity, was modeled by the reservoir simulator TOUGH2. Then, levelized cost of heat (LCOH) was calculated using GEOPHIRES version 3.0, which includes new features such as hourly heat load optimization and peak performance evaluation. Geothermal reservoir temperature was estimated as 90.5 degrees C at a total depth of 3.3 km by the regional average temperature gradient of 22.8 degrees C/km and validated by cation geothermometer calculations. Four production scenarios with two different well configurations and two different heat load profiles have been developed for well flow rates ranging between 3.1 kg/s and 9.3 kg/s. The LCOH of the district heating and cooling system was calculated between $95 and $210/MWh ($28/MMBtu to $62/MMBtu) for four different production scenarios. Typical natural gas prices for residential customers in Oklahoma have ranged from 9 to 19 $/MMBtu over the past decade, which indicates a challenge for deployment of such a GDH system.

deep direct-use↗

DOE OSTI2023

Techno-Economic Analysis for a Potential Geothermal District Heating System in Tuttle, Oklahoma

deep direct-use↗

DOE OSTI2021

Predictive Data-driven Platform for Subsurface Energy Production

Subsurface energy activities such as unconventional resource recovery, enhanced geothermal energy systems, and geologic carbon storage require fast and reliable methods to account for complex, multiphysical processes in heterogeneous fractured and porous media. Although reservoir simulation is considered the industry standard for simulating these subsurface systems with injection and/or extraction operations, reservoir simulation requires spatio-temporal “Big Data” into the simulation model, which is typically a major challenge during model development and computational phase. In this work, we developed and applied various deep neural network-based approaches to (1) process multiscale image segmentation, (2) generate ensemble members of drainage networks, flow channels, and porous media using deep convolutional generative adversarial network, (3) construct multiple hybrid neural networks such as convolutional LSTM and convolutional neural network-LSTM to develop fast and accurate reduced order models for shale gas extraction, and (4) physics-informed neural network and deep Q-learning for flow and energy production. We hypothesized that physicsbased machine learning/deep learning can overcome the shortcomings of traditional machine learning methods where data-driven models have faltered beyond the data and physical conditions used for training and validation. We improved and developed novel approaches to demonstrate that physics-based ML can allow us to incorporate physical constraints (e.g., scientific domain knowledge) into ML framework. Outcomes of this project will be readily applicable for many energy and national security problems that are particularly defined by multiscale features and network systems.

58 GEOSCIENCES↗

DOE OSTI2022

History Matching and Prediction of a Polymer Flood Pilot in Heavy Oil Reservoir on Alaska North Slope

The first-ever polymer flood pilot to enhance heavy oil recovery on Alaska North Slope is ongoing. After more than 3 years of polymer injection, significant benefit has been observed from the decrease in water cut from 65% to less than 15% in the project producers. The primary objective of this study is to develop a robust history-matched reservoir simulation model capable of predicting future polymer flood performance. In this work, the reservoir simulation model has been developed based on the geological model and available reservoir and fluid data. In particular, four high transmissibility strips were introduced to connect the injector-producer well pairs, simulating short-circuiting flow behavior that can be explained by viscous fingering and reproducing the water cut history. The strip transmissibilities were manually tuned to improve the history matching results during the waterflooding and polymer flooding periods, respectively. It has been found that higher strip transmissibilities match the sharp water cut increase very well in the waterflooding period. Then the strip transmissibilities need to be reduced with time to match the significant water cut reduction. The viscous fingering effect in the reservoir during waterflooding and the restoration of injection conformance during polymer flooding have been effectively represented. Based on the validated simulation model, numerical simulation tests have been conducted to investigate the oil recovery performance under different development strategies, with consideration for sensitivity to polymer parameter uncertainties. The oil recovery factor with polymer flooding can reach about 39% in 30 years, twice as much as forecasted with continued waterflooding. Besides, the updated reservoir model has been successfully employed to forecast polymer utilization, a valuable parameter to evaluate the pilot test’s economic efficiency. All the investigated development strategies indicate polymer utilization lower than 3.5 lbs/bbl in 30 years, which is less than that of the same polymer used in a polymer pilot in Argentina.

Wang, Xindan↗