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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Effect of initial water saturation on the performance of fracturing fluids with and without polyallylamine under simulated EGS conditions

Objectives/Scope: StimuFrac (US Patents 9,873,828 B2 and 9,447,315 B2), a CO 2 - reactive polymer aqueous solution [polyallylamine (PAA) 1wt% in water] combined with CO 2 , can be used as a potentially less water-intensive fracturing fluid for enhanced geothermal systems (EGS). Our previous results show that in hot dry rock (HDR), PAA/CO 2 fracturing fluids outperformed other fluids such as water, CO 2 , and CO 2 /water in generating large fractures with less fluid consumed. The objective of this work is to investigate the effect of initial water saturation on the performance of StimuFrac by conducting hydraulic fracturing tests with ½ foot cubic rock samples held under representative EGS stress/temperature conditions and by using cyclic injection strategies (under constant injection rate). The resulting fracture hydraulic conductivities, breakdown pressures, and volumes of fluids required are compared. Methods/Procedures/Process: To simulate geothermal reservoir conditions, in all tests, the rock sample was held under triaxial confinement and at 200 °C, and different volumes of water were initially injected into the rock sample before any fracturing processes were initiated. For the single-cycle PAA (or water) alternating CO 2 (PAG or WAG) injection fracturing experiments, one complete cycle consisted of two steps: (1) injecting a PAA slug (or water slug) followed by (2) injecting CO 2 to initiate and propagate the fracture. For experiments involving multiple injection cycles, the CO 2 injection pressure is increased until it peaks and begins to decline (indicating fracture initiation at this moment), and then continued being injected for another 30 seconds to propagate the fracture. Then, these two-step cycles [injection of PAA (or water) followed by CO 2 injection (up to 2-4 mL/min)] are repeated. Applications/Significance/Novelty: The results of this study suggest that water saturation significantly affects the fracturing fluid transmission into the rock pore space, thus affecting the fracture initiation and propagation. In this study, fracturing tests via a single injection cycle or multiple injection cycles were performed. Splitting the rock samples in half after testing reveals that fracture propagation is significantly limited under high water saturation conditions (three-day initial water injection) compared to stimulation experiments performed in hot dry rock. The fractures propagate less than 1/3 of the distance from the wellbore to the outer rock surface, and in some cases, no fracture is generated. This may be caused by leak-off dominating the fracturing process and the fluid injection rate is insufficient to overcome leak-off, even under high injection rate conditions. Additionally, CO 2 could be leaking off into the wellbore annulus and this may be making it more difficult to generate sufficiently high-pressure gradients away from the near-wellbore region. Under low water saturation conditions (dry rock or after 1-day initial water injection), PAA/CO 2 consistently generated significantly larger fractures compared with the other fluids. CO 2 generated large fractures only in the hot dry rock and only when using high injection rates, though data variability is high.

58 GEOSCIENCES↗
DOE OSTI2021

Joint impedance and facies inversion of time-lapse seismic data for improving monitoring of CO 2 incidentally stored from CO 2 EOR

Time-lapse seismic monitoring is an effective and proven technology for mapping the distribution of CO 2 in a subsurface reservoir. When injected CO 2 displaces other reservoir fluids, porous-medium properties are changed and thus the seismic impedance changes, causing time-lapse seismic amplitude differences in the injection zones. The analysis and interpretation of images created from these amplitude differences can provide information about reservoir architecture and the CO 2 migration within the reservoir. Incorporating seismic inversion and rock physics into the interpretation of time-lapse seismic data can considerably improve the modeling and monitoring to detect and assess the location of CO 2 over time. The joint inversion method presented in this paper has an integral representation of the geology in the inversion algorithm using elastic facies, which provides information about the spatial distribution of the geologic heterogeneities controlling the movement of fluids in the reservoir. The method was successfully applied to time-lapse seismic data from a mature oil field undergoing CO 2 enhanced oil recovery. The estimated seismic acoustic impedances and facies reflect the characteristics of individual geologic facies and fluid conditions of the reservoir subject to CO 2 injection. Finally, the probabilities estimated by the joint impedance and facies inversion for the reservoir's litho-fluid facies can be used for forecasting CO 2 saturation and pressure changes within the target reservoir.

4D seismic↗
DOE OSTI2024

Time-lapse seismic inversion for CO 2 saturation with SeisCO2Net: An application to Frio-II site

Seismic monitoring of geological CO 2 storage (GCS) involves highly nonlinear seismic inversion and petrophysical inversion, making it challenging to estimate CO 2 volume efficiently and detect possible early CO 2 leakages. Deep learning (DL) using convolutional neural networks (CNNs) has shown promise in solving highly nonlinear seismic inversion problems. However, direct estimation of CO 2 plume extent/saturation from time-lapse seismic gathers using DL is still underexplored, with no reported field applications to date. The investigation of field data is primarily hindered by scarcity of field data for neural network training. Other obstacles include highly nonlinear seismic-petrophysics inverse relationship, and presence of noise in field seismic data. We introduce SeisCO2Net, a deep CNN that predicts CO 2 saturation maps directly from time-lapse full waveform shot gathers. For training, we use site-specific geological information, fluid flow physics, rock physics, and seismic modeling to generate synthetic datasets that closely resemble the CO 2 storage site. Synthetic tests show promising results, inspiring us to apply SeisCO2Net's trained weights on field data collected at Frio-II GCS site by leveraging transfer learning principles. As reference, we compare SeisCO2Net's predicted CO 2 saturation maps with results obtained from physics-based inversion. Our analyses show both methods display similar CO 2 plume shapes, reasonable CO 2 plume characteristics, and comparable saturation values. Our results suggest pre-training CNNs on physics-informed synthetic datasets and then applying the learned weights to field data is a viable approach to estimating field CO 2 saturation. This method effectively addresses the scarcity of field training data, thus encouraging the feasibility of long-term GCS monitoring.

58 GEOSCIENCES↗
DOE OSTI2026

Impact of Rock Wettability and Mineralization on CO 2 Storage Efficiency in Basalt Reservoirs

Basalt formations have emerged as highly promising targets for geological CO 2 storage due to their abundance of reactive silicate minerals that can rapidly convert dissolved CO 2 into stable carbonate minerals. The efficiency of CO 2 trapping in basalts, however, depends not only on their geochemical reactivity but also on reservoir wettability, which governs how CO 2 partitions among structural, residual, solubility, and mineral trapping mechanisms. However, wettability and mineral kinetics have largely been examined separately, wettability for plume migration and kinetics in terms of mineralization, leaving their coupled impact unquantified. Here, we systematically integrate th ese processes in a reservoir-scale reactive transport model to quantify the impact of three distinct wettability conditions (water-wet, mixed-wet, and CO 2 -wet) on plume migration, mineral precipitation, pore-scale blockage, and long-term immobilization. Results show that under ideal conditions (no mineralization), water-wet systems exhibited plume aspect ratios up to 28% higher than CO 2 -wet systems and promoted at least 60% greater residual trapping. When mineralization was incorporated, mineral trapping accounted for 57%, 53%, and 46% of the total immobilized CO 2 in water-wet, mixed-wet, and CO 2 -wet systems, respectively. Dissolved CO 2 contributed 43-46%, while residual trapping remained below 5%, and mobile CO 2 was substantially reduced due to mineralization, demonstrating that as the system becomes more water-wet, CO 2 is more effectively converted into stable carbonate minerals, although extensive mineral precipitation may lead to pore clogging. A global Sobol-Morris sensitivity analysis further reveals that the kinetic parameters of reactive minerals are the most influential factor, followed by wettability, while salinity has only a minor influence. Furthermore, these findings highlight that wettability is not merely a flow property but a fundamental control on CO 2 mineralization in basalts, with direct implications for site selection for storage projects.

Basalts↗
DOE OSTI2026

WBM-Hydro-ReEDS

WBM-hydro-ReEDS is a tool for converting daily reservoir release simulations from the WMB model to monthly, plant-level hydropower generation estimates, and then to inputs for the ReEDS model.

Turner, Sean [Oak Ridge National Laboratory (ORNL)↗
DOE OSTI2021

A Non-Invasive Approach for Elucidating the Spatial Distribution of In-Situ Stress in Deep Subsurface Geologic Formations Considered for CO 2 Storage (Final Report)

The primary objectives of this project were to: develop a new, non-invasive method to define the orientation and magnitude of in-situ stresses in deep geologic formations (depths that exceed 1500 meters) considered for carbon dioxide (CO 2 ) sequestration; to demonstrate the method at a real field site; and, to extend current technical capabilities and reduce costs of acquiring information about deep-stress fields. This investigation was segmented into five tasks, including Project Management (Task 1) and four technical tasks (Tasks 2 through 5). The results of the technical tasks are described in this document.

02 PETROLEUM↗
DOE OSTI2024

Impacts of Alternative Operations and Renewable Energy Deployment on Columbia River Hydropower

The Columbia River Treaty Tribes in the Pacific Northwest - the Nez Perce, Umatilla, Warm Springs, and Yakama - hold treaty-reserved fishing rights for the Columbia River, one of the world's most productive salmon rivers and a critical resource for these tribes. However, the tribes have expressed that current operating regimes of hydropower dams throughout the Columbia River Basin (CRB) do not fully account for tribal fishing rights and have negatively impacted fish populations. Four tribes acting together through the Columbia River Inter-Tribal Fish Commission (CRITFC)'s 2022 Energy Vision stated that current power system models often do not fully account for the many constraints faced by hydropower facilities in the CRB. As the deployment of variable renewable energy (VRE) like solar and wind continues to increase, power system flexibility will become increasingly important. As a result, there is a growing need to better understand the true capabilities of the hydropower generation fleet while accurately accounting for ecological constraints (Northwest Power and Conservation Council 2022). Understanding hydropower's role in a future grid with a higher share of VRE can inform water resources planning to address ecological needs. The goal of this study is to examine the impacts of alternative hydropower operation rules and weather variability on hydropower generation and grid operations in VRE and transmission infrastructure deployment scenarios. The study evaluates how hydropower scheduling practices can reduce ecological impacts to the region's salmon populations. More specifically, the study examines the impact of today's dam water release rules (called "adjusted base water rules" in this study) and new, ecologically informed water release rules (called "ecosystem water rules" in this study) on two VRE scenarios (current renewable energy levels and higher renewable energy levels) using six weather years for each scenario (2008-2013). CRITFC developed the ecosystem water rules, which capture a portion of the changes they recommend for CRB hydropower operations. These analyses use a water resource planning model (OASIS) and a production cost model (PLEXOS) to simulate CRB reservoir cascade and Western Interconnection power grid operation.

13 HYDRO ENERGY↗
DOE OSTI2026

Oil in Alaska North Slope Gas Hydrate Reservoir: Micro-CT and Flow Simulation Insights into Permeability

Gas hydrate-bearing sands on the Alaska North Slope (ANS) host minor volumes of crude oil whose impact on formation permeability has never been quantified relative to a gas hydrate reservoir system. Here, we combine in situ pressure-core microcomputed-tomography (μ-CT), thermogravimetric analysis (TGA), gas-chromatography–mass-spectrometry (GC–MS), scanning-electron microscopy coupled with energy-dispersive spectroscopy (SEM-EDS), and Stokes-flow simulation to (i) quantify oil saturation, (ii) infer its source, and (iii) evaluate its influence on permeability after gas-hydrate dissociation. μ-CT and SEM-EDS imaging identify the host sediment as silty, with some grain-coating clay present. Eight representative 3003-voxel subvolumes extracted from a preserved pressure core (Core 15P-3c, from the B1 sand (unit B), HYDRATE 02 Geo Data Well (GDW), 886.96–887.29 m measured depth, MD) exhibit porosities of 0.34–0.37 and an average oil saturation of 0.09 ± 0.03, in comparison to an independent TGA-based estimate of 0.16. GC-MS fingerprints obtained from an adjacent depressurized core (Core 17P-1, 891.37–891.44 m MD), together with oil pore habit revealed by μ-CT, suggest that the oil is partly native but also includes contributions from mineral oil-based drilling fluid contamination. Further analysis of oil–porewater interactions and flushing tests confirms that native oil saturation in the analyzed samples ranges from 0.04 to 0.08 with an average of approximately 0.06, and the oil is highly immobile. Flow simulations demonstrate that the native oil saturation of ≈0.06 reduces permeability by approximately 50% due to both pore blockage and increased flow-path tortuosity. These findings reveal the potential presence of native oil within the ANS gas hydrate reservoir and suggest that future simulation models may need to account for its impact on permeability to improve long-term performance predictions of gas and water production.

02 PETROLEUM↗
DOE OSTI2020

Charged Wellbore Casing Controlled Source Electromagnetics (CWC-CSEM) for Reservoir Imaging and Monitoring (Final report)

This project addresses the needs of the U.S. Department of Energy (DOE) to develop advanced monitoring technologies and protocols to track the fate of subsurface carbon dioxide (CO2) plumes for carbon storage. Specifically, the project seeks to develop and test a unique and novel system of technologies consisting of electromagnetic data acquisition, coupled multiphysics imaging, and reservoir model enhancement to understand the migration and long-term distribution of CO2 in the subsurface. The overarching objective is to develop an integrated approach for long term monitoring of carbon storage. The two main components of the project include the methodology development and the test of the method at a field site. The methodology component consists of 1) developing the field procedure and protocol for collecting time-lapse controlled-source electromagnetic (CSEM) data with source electric current injected into the subsurface through wellbore casings; 2) building of background 3D electrical conductivity utilizing multiple sources of data such as supplemental surface transient EM (TEM) surveys, well-logs, and seismic structural information, for enhancing CSEM signal from reservoir depths; and 3) coupled multiphysics simulations and inversion of CSEM data constrained by production data and by structural information from seismic imaging of the reservoir and overlying formations. The testing component used the field site of Bell Creek Oil Field, which served both as a field laboratory for the method development as well as a test site to evaluate the CSEM signal strengths and the methodology developed in this research project. We have accomplished all the proposed tasks and developed the methodology as planned. These include the procedure for time-lapse CSEM data acquisition, data processing techniques, integration with 3D conductivity model building, fast reservoir simulation for history matching using machine learning, and interpreting CSEM data with coupling to the reservoir modeling. Collectively, the outcome of these tasks form a coherent workflow that can be applied to monitor dedicated carbon storage in saline reservoirs. The testing component evaluated the applicability and limitations of the method, and concluded that the method would be ideal for monitoring dedicated carbon storage sites utilizing saline reservoirs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

High-Resolution Simulations of Geological CO 2 Injection: Application to the SPE11 Benchmark

Geological carbon sequestration (GCS) will play a critical role in decarbonization and in facilitating the transition to clean energy systems. Because CO 2 is highly mobile, ensuring its safe and permanent injection into subsurface geological formations involves monitoring over larger spatial domains and longer time periods than is typical for hydrocarbon reservoirs. This can benefit from simulation tools capable of modeling key CO 2 trapping mechanisms, particularly those optimized for speed and scalability on high-performance computing systems. Using isothermal versions of the SPE11B and SPE11C benchmark cases, we conduct a mesh refinement study simulating CO 2 injection into kilometer-scale rock formations at centimeter resolution with the GEOS open-source simulation framework. We focus on how mesh refinement improves the accuracy of convective mixing in both 2D and 3D simulations. The computational costs associated with achieving a converged solution highlight the need for predictive upscaling techniques. A systematic performance scaling analysis—including both central processing unit (CPU) and graphics processing unit (GPU) architectures—complements the “Results” section.

Geosciences↗
DOE OSTI2026

ELM‐MOSART‐DOC: A Large‐Scale Riverine Dissolved Organic Carbon Model and Its Application Over the United States

Riverine dissolved organic carbon (DOC), primarily sourced from soil organic carbon (SOC), plays a crucial role in regional and global carbon cycles. However, the complexities of the underlying mechanisms and limited observations present significant challenges for predictive understanding of DOC at regional or larger scales. Recently, we developed a machine learning‐based (ML) map of DOC transformation rates, bridging the gap between SOC and DOC leaching flux and simplifying terrestrial DOC representation. Building on this advancement, we introduce ELM‐MOSART‐DOC, a DOC module integrated into the riverine component of the Energy Exascale Earth System Model (E3SM)—the Model for Scale Adaptive River Transport (MOSART). ELM‐MOSART‐DOC simulates DOC transport and transformation across both headwater streams and river networks, including those managed. Model validation demonstrates the ability of ELM‐MOSART‐DOC to accurately capture long‐term average DOC concentrations, with Kling‐Gupta Efficiency (KGE) scores of 0.58 and 0.76 at large and local stations, respectively. We further assess the impact of reservoirs through different simulation schemes, revealing that reservoirs significantly alter DOC fluxes by regulating streamflow patterns and promoting DOC mineralization. Model simulations indicate that reservoirs reduce total DOC flux from the Mississippi River into the ocean by 7.5%, with the long‐term average annual export decreasing from 3.34 to 3.14 teragrams (Tg) per year. ELM‐MOSART‐DOC integrates process‐based modeling with ML parameterization to enhance the predictive understanding of riverine biogeochemical processes. This approach reduces uncertainties in modeling regional and global carbon cycle ESMs and provides new insights into carbon cycling and its implications for global environmental change.

Li, Lingbo [Univ. of Houston, TX (United States); ↗
DOE OSTI2023

Enhancing the representation of water management in global hydrological models

Abstract. This study enhances an existing global hydrological model (GHM), Xanthos, by adding a new water management module that distinguishes between the operational characteristics of irrigation, hydropower, and flood control reservoirs. We remapped reservoirs in the Global Reservoir and Dam (GRanD) database to the 0.5∘ spatial resolution in Xanthos so that a single lumped reservoir exists per grid cell, which yielded 3790 large reservoirs. We implemented unique operation rules for each reservoir type, based on their primary purposes. In particular, hydropower reservoirs have been treated as flood control reservoirs in previous GHM studies, while here, we determined the operation rules for hydropower reservoirs via optimization that maximizes long-term hydropower production. We conducted global simulations using the enhanced Xanthos and validated monthly streamflow for 91 large river basins, where high-quality observed streamflow data were available. A total of 1878 (296 hydropower, 486 irrigation, and 1096 flood control and others) out of the 3790 reservoirs are located in the 91 basins and are part of our reported results. The Kling–Gupta efficiency (KGE) value (after adding the new water management) is ≥ 0.5 and ≥ 0.0 in 39 and 81 basins, respectively. After adding the new water management module, model performance improved for 75 out of 91 basins and worsened for only 7. To measure the relative difference between explicitly representing hydropower reservoirs and representing hydropower reservoirs as flood control reservoirs (as is commonly done in other GHMs), we use the normalized root mean square error (NRMSE) and the coefficient of determination (R2). Out of the 296 hydropower reservoirs, the NRMSE is > 0.25 (i.e., considering 0.25 to represent a moderate difference) for over 44 % of the 296 reservoirs when comparing both the simulated reservoir releases and storage time series between the two simulations. We suggest that correctly representing hydropower reservoirs in GHMs could have important implications for our understanding and management of freshwater resource challenges at regional-to-global scales. This enhanced global water management modeling framework will allow the analysis of future global reservoir development and management from a coupled human–earth system perspective.

13 HYDRO ENERGY↗
DOE OSTI2021

A simulation study of carbon storage with active reservoir management

As part of the Integrated Midcontinent Stacked Carbon Storage Hub (IMSCS-Hub) project led by Battelle Memorial Institute, a study was conducted to determine the feasibility of storing carbon dioxide (CO 2 ) in the stacked saline rock formations of the Sleepy Hollow Field (SHF), located in Red Willow County, southern Nebraska. A series of CO 2 injection simulation scenarios, with and without active reservoir management (ARM; brine extraction), were evaluated to investigate the feasibility of storing 50+ million tonnes (Mt) of CO 2 . The results indicated CO 2 injection combined with ARM may enable permanent storage of 50+ Mt of CO 2 . The area of review (AOR), the area in which underground sources of drinking water (USDWs) might be endangered during CO 2 injection, was assessed for the simulation scenarios. In comparison to a case without ARM, brine extraction resulted in a much smaller AOR, covering an area of 42 square miles (108.8 km 2 ), roughly one-fourth the size of an AOR resulting from CO 2 injection without ARM. Our findings, presented in this paper, indicate that CO 2 injection with ARM can improve the CO 2 storage capacity of a geologic storage complex up to 100% while also reduce the rate of pressure buildup in the subsurface, resulting in a 75% reduction in AOR. This may help in lowering carbon capture and storage project costs, risks, and effort needed to meet monitoring requirements for a storage project.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Reservoir Dynamics in Proposed Operational Scenarios Where CO 2 -EOR Fields are Transitioned From CO 2 -Flood Enhanced Oil Recovery to Dedicated Carbon Storage: A Field Case Study

This study models the transition of CO 2 -enhanced oil recovery (CO 2 -EOR) fields to dedicated carbon storage using a generalized reservoir model informed by real field data, SACROC. Through scenario-based simulations, the work evaluates how reservoir depletion levels, boundary conditions, fluid properties, and domain size affect pressure buildup, CO 2 plume migration, and long-term containment performance. The results inform practices for repurposing oil fields into reliable CO 2 storage sites.

02 PETROLEUM↗
DOE OSTI2021

Developing Upscaling Approach for Swarming Hydraulic Fractures Observed at Hydraulic Fracturing Test Site through Multiscale Simulations

We show that this work aims to address a challenge posed by recent observations of tightly spaced hydraulic fractures in core samples from the hydraulic fracturing test site (HFTS) in the Middle Wolfcamp Formation. Many fractures in retrieved cores have subfoot spacing, which is at odds with conventional models in which usually one hydraulic fracture is initiated per cluster. Models assuming a single fracture at each cluster, although a common practice, often predict excessive fracture propagation that is inconsistent with microseismic observation. Here, we aim to develop a numerical approach to effectively account for densely spaced hydraulic fractures in field-scale simulations. Because it is impractical to explicitly model all aforementioned fractures, we develop a new upscaling law that enables existing simulation tools to predict reservoir response to fracture swarms. The upscaling law is derived based on an energy equivalence argument and validated through multiscale simulations using a high-fidelity code, GEOS. The swarming fractures are first modeled with a spacing that is much smaller than the cluster spacing; these fractures are then approximated by an upscaled, single fracture based on the proposed upscaling law. The upscaled fracture is shown to successfully match the energy input rate and produce the total fracture aperture and average propagation length of the explicitly simulated swarm. Afterward, the upscaling approach is further implemented in 3D field-scale simulations and validated against the HFTS microseismic data of a horizontal well. Our results show that hydraulic fracture swarming can significantly affect fracture propagation behaviors compared with the propagation of single fractures as assumed by conventional modeling approaches. Under the considered situations, the conventional treatment yields fast propagation speed that far exceeds that indicated by the microseismic data. We also illustrate that this discrepancy can be reduced readily through the implementation of the upscaling law. Our results demonstrate the importance of accounting for the fracture swarming effect in field-scale simulations and the efficacy of this approach to enable realistic predictions of reservoir responses to fracture swarms, without the need to model tightly spaced fractures individually.

02 PETROLEUM↗
DOE OSTI2022

Sand Migration Simulation during Gas Production from Gas Hydrate Reservoir at Kuparuk 7–11–12 site in the Prodhoe Bay Unit, Alaska

Uncontrolled sand production impedes continuous gas production from a hydrate reservoir as observed in the past field-scale gas production tests. Sand mobilization is strongly linked with sediment deformation and high pressure gradient. Sand production occurs when the mobilized sands reach the well. Throughout gas production from a hydrate reservoir, because of hydrate dissociation, deformation and pressure gradient change both in time and space. Therefore, it is necessary to consider the entire process to identify where sand is likely mobilized and how much mobilized sand could reach the well. Here, this study utilizes a coupled thermal, hydrological, chemical, and mechanical (thermo-hydro-chemo-mechanical) sand migration model to simulate a year long gas production from the hydrate reservoir at the Kuparuk 7–11–12 site in the Prodhoe Bay Unit, Alaska. It is found that sand would mainly come from the lower portion of the production zone where faster hydrate dissociation occurs. The relatively faster hydrate dissociation coupled with the fact that hydrate-bearing sediments and hydrate-free sediments have similar stiffness but different strengths causes complex stress transfer, which results in excessive shear deformation when the upper portion of the production zone starts to dissociate. This is evident not only near the well but also away from the well, leading to a large amount of sand mobilization. Another focus of the modeling study is evaluation of the effect of sand migration on gas production. The comparison of two extreme cases–with sandscreen and without sandscreen–suggests that migrated and settled solids around a sand control device can prevent pressure drop across the near-wellbore zone and lead to reduction of the gas production rate.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Simulation study of chemo-mechanical impacts of CO 2 injection in morrow b sandstone reservoir

Here, this study investigates the impacts of geomechanical and geochemical changes on carbon storage in a partially depleted oil reservoir, using results from four different coupled simulation models. Models were used to examine the relative importance of storage mechanisms, and how changing reservoir parameters might affect these mechanisms through time. The study uses data from a Morrowan sandstone reservoir in the Farnsworth Unit (FWU), Ochiltree County, Texas which is currently undergoing CO 2 enhanced oil recovery (EOR). Partially depleted oil reservoirs such as the FWU offer attractive carbon utilization and/or storage targets because of existing infrastructure and economic benefits from incremental oil recovery as well as tax credits. However, prediction of storage capacity or long-term fluid migration in these fields can be difficult because of the wide variation in formation fluids and operational histories that may have undergone. CO 2 injection can cause complex geomechanical and geochemical responses in a reservoir as a result of interplay between dynamic changes in pore pressure, reservoir temperature, fluid composition, and interactions between formation fluids, CO 2 , and reservoir rock. Thus, multiple coupled numerical simulation models must be developed and used to more precisely understand what CO 2 storage mechanisms are most significant, as well as the long-term fate of the stored CO 2 . Our study used results from hydrodynamic, coupled hydro-geomechanical, coupled hydro-geochemical, and coupled hydro-geomechanical-geochemical models to examine how changes in geomechanical and geochemical properties can impact the injectivity or storage capacity of CO 2 . Models simulated historical field operations and then forward-modeled a water-alternate-gas (WAG) operation for 20 years, followed by a 1000-year post-injection monitoring. The work demonstrates that in this specific reservoir, geomechanical impacts appear to be more significant than any geochemical processes, which showed minimal impact on reservoir properties. Furthermore, these simulation scenarios corroborated other studies showing more than 90% of CO 2 injected are being stored and will remain stored in the reservoir. The study was able to draw on a rich set of data from this currently operational project for the work. Because the Morrow is a widespread oil-bearing reservoir in the midcontinent of the United States and is fairly typical of many sandstone reservoirs elsewhere, the study has significance beyond the boundaries of our field area.

58 GEOSCIENCES↗