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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
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↗
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.↗
NASA NTRS2019

Mars Atmospheric Chemistry Simulations with the GEM-Mars General Circulation Model

General Circulation Models with interactive physical and chemistry processes are the state-of-the-art tools for an integrated view and understanding of the Martian atmosphere and climate system. The GEM-Mars model currently includes 16 tracers for chemical composition and applies a fully online, interactive calculation of the photo- and gas phase chemistry of carbon dioxide (CO2) and water vapor (H2O). These species largely control the chemical composition of the neutral Mars atmosphere through their photolysis products and their subsequent interactions. Water vapor undergoes a complex cycle on Mars as it is transported and interacts with ice reservoirs both at the planet's surface and at water ice clouds, which in turn provide radiative feedbacks. In the photochemical cycles involving CO2 and H2O, the abundances of 5 species have been reported by previous investigations with significant spatio-temporal coverage: CO2, H2O, CO, O3 and H2O2. This paper presents the current status of the atmospheric chemistry simulations in GEM-Mars by comparing them to a selection of these observational datasets as well as to oxygen dayglow emission from O2(a1∆g). The results are consistent with previous model-data comparisons and illustrate that the water cycle and the photochemistry are well implemented in the model. In particular, the simulation of the key reservoir species H2O2 provides a good match to the available data. Model-data biases for ozone columns and oxygen airglow are related to the simulated water vapor vertical profile, as these species have important column contributions from vertical layers at the top of the hygropause.

F Daerden↗
NASA NTRS2009

Sodium Variable Conductance Heat Pipe for Radioisotope Stirling Systems

In a Stirling radioisotope system, heat must continually be removed from the General Purpose Heat Source (GPHS) modules to maintain the modules and surrounding insulation at acceptable temperatures. Normally, the Stirling convertor provides this cooling. If the converter stops in the current system, the insulation is designed to spoil, preventing damage to the GPHS, and also ending the mission. An alkali-metal Variable Conductance Heat Pipe (VCHP) has been designed to allow multiple stops and restarts of the Stirling convertor in an Advanced Stirling Radioisotope Generator (ASRG). When the Stirling convertor is turned off, the VCHP will activate when the temperatures rises 30 C above the setpoint temperature. A prototype VCHP with sodium as the working fluid was fabricated and tested in both gravity aided and against gravity conditions for a nominal heater head temperature of 790 C. The results show very good agreement with the predictions and validate the model. The gas front was located at the exit of the reservoir when heater head temperature was 790 C while cooling was ON, simulating an operating Advanced Stirling Converter (ASC). When cooling stopped, the temperature increased by 30 C, allowing the gas front to move past the radiator, which transferred the heat to the case. After resuming the cooling flow, the front returned at the initial location turning OFF the VCHP. The against gravity working conditions showed a colder reservoir and faster transients.

Tarau, Calin↗
DOE OSTI2024

Unveiling the Potential of MeshGraphNets for Predicting Subsurface Evolution in Carbon Storage Projects

This is the conference paper accompanying an oral presentation “Unveiling the Potential of MeshGraphNets for Predicting Subsurface Evolution in Carbon Storage Projects” at the 17th International Conference on Greenhouse Gas Control Technologies GHGT-17 held in Calgary, Canada, October 20-24 , 2024. Carbon capture and storage (CCS) technology is critical for mitigating climate change but requires effective subsurface reservoir management to ensure safe containment of injected CO2. Accurate predictions of reservoir pressure and saturation are essential for assessing long-term CCS performance. Traditional numerical simulations, while effective, are computationally intensive, time-consuming, and constrained by data discretization. Previous work has shown the effectiveness of MeshGraphNets (MGN), a graph-based machine learning framework, as an innovative alternative for predicting reservoir behavior. MGN leverages graph neural networks (GNNs) and mesh representations to model complex geological formations, offering superior adaptability across different discretizations and reservoir configurations. Classic MGN implementations utilize an autoregressive technique to predict future behavior based on current predictions, but this technique is hampered by error accumulation over time. To enhance the model accuracy in time-series predictions, this study implemented a multi-step rollout strategy that integrates autoregressive predictions during training to stabilize prediction of saturation over time. Using the Illinois Basin – Decatur Project (IBDP) dataset, comprising 100 simulations of CO2 injection, pressure, and saturation changes, the framework demonstrated its ability to learn spatial dependencies and temporal dynamics. With inputs including permeabilities, porosities, and injection rates, MGN accurately predicted CO2 plume evolution over time, even with limited training data. Moreover, the addition of a multi-step rollout procedure during training improved the ability of MGN to predict stably over time by ~15%. This research positions MGN, enhanced with multi-step rollout capabilities, as a robust and efficient tool for CCS applications. It advances the field by enabling precise, computationally efficient predictions of reservoir behavior, providing a foundation for the broader adoption of machine learning frameworks in CCS and other geoscience domains.

Holcomb, Paul↗
DOE OSTI2024

Geologic Characterization of the South Georgia Rift Basin for Source Proximal CO2 Storage

The project Geologic Characterization of the South Georgia Rift Basin for Source Proximal CO2 Storage is one of 9 site characterization projects that were implemented as part of ARRA (American Recovery and Reinvestment Act). Data from this project was used to improve resolution of data in NATCARB in the area of study. Data related to this study has already been incorporated in NATCARB Atlas. The South Carolina Research Foundation and partners evaluated the feasibility of CCS in the Jurassic/ Triassic (J / TR) saline formations of the buried Mesozoic South Georgia Rift (SGR) Basin that extends from South Carolina into Georgia. The J / TR sequence, based on preliminary assessment of limited geologic and geophysical data, appears to have both the appropriate areal extent and multiple horizons to permanently and safely store CO2 The presence of several igneous rock layers within the sequence may potentially provide adequate seals to prevent upward CO2 migration into the Coastal Plain aquifer systems. Approximately 81 kilometers of 2-D seismic reflection data were collected by Bay Geophysical, Inc. to explore a portion of the SGR located in southern Georgia. The 81 kilometers were divided into two lines approximately 40.5 kilometers each, with Line 1 intersecting Georgia well GGS 3457. Line 2 intersects Line 1 at the southern portion of Line 1 to maximize the extent of coverage away from GGS-3457 (a deep well drilled in the 1980s for oil and gas exploration). This well had a set of usable logs, including gamma and neutron logs that provided promising results related to CO2 storage. Results showed sandstone with porosity values greater than 10 percent and a thickness of 120 meters. The design of the seismic shot was to extrapolate information away from the well and to better define the extent of the SGR and the necessary reservoir and caprock for a successful CO2 injection. A numerical simulation model of CO2 Injection and migration was developed based on the geology log for the GGS-3457 well. The simulation model was used to investigate the feasibility of injecting 30 million metric tons of CO2 into SGR J / TA sediments and integrity of the diabase layers as seals to prevent CO2 migration.

2-D seismic↗
NASA NTRS1992

Thermal and dynamical regimes of single- and two-phase magmatic flow in dikes

The coupling between thermal and dynamical regimes of single- and two-phase magmatic flow in dikes, due to temperature-dependent viscosity and dissipation, was investigated using finite element calculations of magma flow in dikelike channels with length-to-width ratios of 1000:1 or more. Solutions of the steady state equations governing magma flow are obtained for a variety of conditions ranging from idealized plane-parallel models to cases involving nonparallel geometry and two-phase flows. The implications of the numerical simulations for the dynamics of flow in a dike-reservoir system and the consequences of dike entrance conditions on magmatic storage are discussed. Consideration is also given to an unmixing/self-lubrication mechanism which may be important for the lubrication of silicic magmas rising to the earth's surface in mixed magma ascent scenarios, which naturally segregates magma mixtures of two components with differing viscosities to minimize the driving pressure gradient.

Carrigan, Charles R.↗
DOE OSTI2024

A Novel 'Smart Microchip Proppants' Technology for Precision Diagnostics of Hydraulic Fracture Networks (Edited Final Report)

This project introduces innovative technology to improve subsurface characterization, visualization, and diagnostics of unconventional reservoirs (fossil resources). Through a collaborative effort involving the University of Kansas, UCLA, MicroSilicon Inc., and EOG Resources, the project aims to deliver precision diagnostics for hydraulic fractures using novel high-resolution imaging technology based on smart microchip proppants. Additionally, it seeks to enhance the accuracy and predictability of integrated numerical, and machine-learning modeling techniques for hydraulic fracture characterization and simulation. This groundbreaking technology addresses significant gaps in understanding unconventional and tight reservoir behavior and optimizing well-completion strategies, enabling more cost-efficient recovery of unconventional resources.

02 PETROLEUM↗
DOE OSTI2026

H I Depletion Begins Well Beyond the Virial Radius: A FAST Stacking Study of 36 Galaxy Clusters to 5 × R 200

Abstract We present a stacking study of the neutral atomic hydrogen (H i ) content in and around 36 local galaxy clusters at z < 0.07, using a combination of the FAST All Sky H i survey (FASHI) and the extensive spectroscopic catalog mainly from the Dark Energy Spectroscopic Instrument (DESI). We employ spectral stacking techniques to probe the average H i mass and HI-to-stellar mass ratio ( M HI / M * ) for member galaxies down to stellar masses of M * ∼ 10 9 M ⊙ , spanning a projected cluster-centric distance of up to 5 R 200 . Our analysis reveals a pronounced environmental effect; both M HI and M HI / M * decrease steadily toward the cluster center, dropping by ∼0.5 dex on average from the outskirts to the core. Crucially, we find that M HI / M * of galaxies remain lower than the field galaxies even at the 5 R 200 . This provides direct, statistical evidence for substantial gas stripping and preprocessing in the cluster outskirts, likely occurring in infalling groups and large-scale filaments. By further splitting the sample by g − r color, we show that the H i deficiency persists at fixed galaxy color; even the bluest cluster members exhibit ∼0.5 dex lower M HI / M * than field galaxies of similar color, reflecting environmental effects on the cold gas reservoir prior to full optical transformation. The total H i mass within clusters and their outskirts agrees broadly with predictions from cosmological simulations. Our results underscore the critical role of the extended cluster environment in quenching galaxies by depleting their cold gas reservoirs well before they enter the dense cluster core.

Cheng, Cheng [Chinese Academy of Sciences South Am↗
NASA NTRS1997

A blast of gas in the latest Paleocene: simulating first-order effects of massive dissociation of oceanic methane hydrate

Carbonate and organic matter deposited during the latest Paleocene thermal maximum is characterized by a remarkable -2.5% excursion in delta 13C that occurred over approximately 10(4) yr and returned to near initial values in an exponential pattern over approximately 2 x 10(5) yr. It has been hypothesized that this excursion signifies transfer of 1.4 to 2.8 x 10(18) g of CH4 from oceanic hydrates to the combined ocean-atmosphere inorganic carbon reservoir. A scenario with 1.12 x 10(18) g of CH4 is numerically simulated here within the framework of the present-day global carbon cycle to test the plausibility of the hypothesis. We find that (1) the delta 13C of the deep ocean, shallow ocean, and atmosphere decreases by -2.3% over 10(4) yr and returns to initial values in an exponential pattern over approximately 2 x 10(5) yr; (2) the depth of the lysocline shoals by up to 400 m over 10(4) yr, and this rise is most pronounced in one ocean region; and (3) global surface temperature increases by approximately 2 degrees C over 10(4) yr and returns to initial values over approximately 2 x 10(6) yr. The first effect is quantitatively consistent with the geologic record; the latter two effects are qualitatively consistent with observations. Thus, significant CH4 release from oceanic hydrates is a plausible explanation for observed carbon cycle perturbations during the thermal maximum. This conclusion is of broad interest because the flux of CH4 invoked during the maximum is of similar magnitude to that released to the atmosphere from present-day anthropogenic CH4 sources.

Non-NASA Center↗
NASA NTRS1984

Spatial analysis in recreation resource management for the Berlin Lake Reservoir Project

Spatial analysis of geographic information systems and the acquisition and use of remotely-sensed data within the U.S. Army Corps of Engineers is an emerging Technology Work units have been developed under te Remote Sensing Research and Development Program, which are most relevant to the productive needs of the Corps in both the military and civil works missions. Corps participation in the SPOT simulation champaign is one such example of this research. This paper describes the application of spatial analysis and remote sensing in recreation resource managmaster planning at the Berlin Lake Reservoir Project within the Pittsburgh District. SPOT simulator data was acquired over Berlin Lake, Site No. 10, on July 8, 1983. The first part of this paper describes the background of the U.S. Army Corps of Engineers and the Berlin Lake project, the geographic information system being developed, and the planned use of SPOT and similar data. The remainder of the paper describes the results on an analysis of the simulated SPOT data conducted at the NASA Goddard Institut for Space Studies.

Edwardo, H. A.↗
NASA NTRS1990

CELSS simulations for a lunar outpost

The paper describes a five-year simulation of two mission scenarios consisting of 14 different 'events' that could take place at a lunar outpost. The time-dependent status of the life support consumables was calculated in response to the two selected mission scenarios. This application demonstrates that complex sequences of events are reproducible for understanding of integrated mission operations. A Controlled Ecological Life Support System (CELSS) can regenerate air, water, and food from the wastes generated in the habitat and do so with safety and reliability. The CELSS Emulator is under development at the NASA's Johnson Space Center with a purpose of investigating computer simulations of integrated CELSS operations involving humans, plants, process machinery, and reservoirs. The Emulator Version 2.0 has been implemented to provide a mission-scenario-analysis capability. Thus, the future space exploration missions, lunar or Mars, can be analyzed by using 'events' to build mission timelines.

Cullingford, H. S.↗
NASA NTRS1993

Models for direct Monte Carlo simulation of coupled vibration-dissociation

A new model for reactive collisions is developed within the framework of a particle method, which simulates coupled vibration-dissociation (CVD) behavior in high-temperature gases. The fundamental principles of particle simulation methods are introduced with particular attention given to the probability functions employed to select thermal and reactive collisions. Reaction probability functions are derived which favor vibrationally excited molecules as reaction candidates. The new models derived here are used to simulate CVD behavior during thermochemical relaxation of constant-volume O2 reservoirs, as well as the dissociation incubation behavior of postshock N2 flows for comparisons with previous models and experimental data.

Haas, Brian L.↗
DOE OSTI2025

Southwest Regional Partnership on Carbon Sequestration: Phase III (Final Scientific/Technical Report)

The Southwest Regional Partnership on Carbon Sequestration (SWP) is one of 7 regional partnerships formed in 2003 under the U.S. Department of Energy’s (DOE) Regional Carbon Sequestration Partnerships (RCSPs) initiative. The overall purpose of the initiative was to help determine and implement the technology, infrastructure, and regulations most appropriate to promote carbon storage in different regions of the country. Covering Arizona, Colorado, New Mexico, Oklahoma, Utah, and parts of Texas, Wyoming, and Kansas, the SWP evaluated regional carbon storage and utilization potential and focused on technologies and sites that could complement the region’s strong position in energy production. The project progressed through three phases: • Phase I (2003–2005): Characterized regional geologic formations and CO 2 sources, assessed sequestration potential, and identified pilot test sites. • Phase II (2005–2013): Conducted small-scale field tests to validate sequestration methods, including geologic and terrestrial projects. • Phase III (2008–2022): Demonstrated large-scale CO 2 injection at a commercial oil field to test monitoring, verification, and long-term storage strategies. This report covers Phase III. The final project site, the Farnsworth Unit (FWU) in Texas, provided real-world testing of reservoir characterization, monitoring, and risk evaluation tools and processes that could be used in any commercial scale carbon capture, utilization, and storage (CCUS) project. Extensive data collection and analysis helped refine best practices for reservoir characterization, injection monitoring, and storage verification. The SWP contributed to national databases, DOE best practice manuals, and regional geological assessments to support future sequestration efforts. Key lessons learned include the importance of robust data management, strategic site selection, regulatory navigation, and effective industry collaboration. The project’s findings will inform ongoing and future carbon storage initiatives. Task 1 (Regional Characterization) • The SWP continued to participate in national outreach efforts and NATCARB. • The SWP evaluated multiple potential sites before selecting the FWU as the primary field test location. Task 2 (Public Outreach and Education) • The SWP contributed to national databases, DOE best practice manuals, and regional geological assessments to support future sequestration efforts. Task 3 (Permitting and Regulatory Compliance) • The SWP ensured compliance with federal and state regulations, including National Environmental Policy Act (NEPA) requirements. • The SWP obtained all necessary permits for drilling, injection, and monitoring activities. Task 4 (Site Characterization and Planning) • The SWP developed work plans for four key activities: characterization, simulation, monitoring and verification, and risk evaluation. • The SWP collected and synthesized legacy data from multiple sources to build initial static geological models and dynamic reservoir models demonstrating project feasibility. • The SWP conducted an initial risk evaluation and developed mitigation plans. Task 5 (Field Operations and Data Collection) • The SWP drilled, logged, and cored three characterization wells to gather critical subsurface data. • The SWP conducted multiple geophysical surveys, including 3D seismic, crosswell seismic, and vertical seismic profiling, to improve reservoir characterization. Task 6 (Monitoring and Verification) • The SWP performed extensive geological characterization using data from characterization wells and seismic surveys. • The SWP established a surface monitoring network to track CO 2 flux in soil gas, groundwater chemistry, and near-surface atmospheric CO 2 levels. • The SWP built and refined reservoir models to study the effects of relative permeability on simulation behavior and improve calibration with experimental data. Task 7 (Risk Assessment and Model Refinement) • The SWP conducted multiple studies to evaluate reservoir integrity, predict CO 2 plume behavior and improve predictive modeling capabilities. • The SWP refined geological models and used them to enhance the accuracy of simulation models. • The SWP continued quantitative risk assessment of top-ranked risks and strengthened the link between qualitative and quantitative risk methodologies.

02 PETROLEUM↗
DOE OSTI2024

3D modeling of deep borehole electromagnetic measurements with energized casing source for fracture mapping at the Utah Frontier Observatory for Research in Geothermal Energy

Here, we present a 3D numerical modelling analysis evaluating the deployment of a borehole electromagnetic measurement tool to detect and image a stimulated zone at the Utah Frontier Observatory for Research in Geothermal Energy geothermal site. As the depth to the geothermal reservoir is several kilometres and the size of the stimulated zone is limited to several 100 m, surface-based controlled-source electromagnetic measurements lack the sensitivity for detecting changes in electrical resistivity caused by the stimulation. To overcome the limitation, the study evaluates the feasibility of using a three-component borehole magnetic receiver system at the Frontier Observatory for Research in Geothermal Energy site. To provide sufficient currents inside and around the enhanced geothermal reservoir, we use an injection well as an energized casing source. To efficiently simulate energizing the injection well in a realistic 3D resistivity model, we introduce a novel modelling workflow that leverages the strengths of both 3D cylindrical-mesh-based electromagnetic modelling code and 3D tetrahedral-mesh-based electromagnetic modelling code. The former is particularly well-suited for modelling hollow cylindrical objects like casings, whereas the latter excels at representing more complex 3D geological structures. In this workflow, our initial step involves computing current densities along a vertical steel-cased well using a 3D cylindrical electromagnetic modelling code. Subsequently, we distribute a series of equivalent current sources along the well's trajectory within a complex 3D resistivity model. We then discretize this model using a tetrahedral mesh and simulate the borehole electromagnetic responses excited by the casing source using a 3D finite-element electromagnetic code. This multi-step approach enables us to simulate 3D casing source electromagnetic responses within a complex 3D resistivity model, without the need for explicit discretization of the well using an excessive number of fine cells. We discuss the applicability and limitations of this proposed workflow within an electromagnetic modelling scenario where an energized well is deviated, such as at the Frontier Observatory for Research in Geothermal Energy site. Using the workflow, we demonstrate that the combined use of the energized casing source and the borehole electromagnetic receiver system offer measurable magnetic field amplitudes and sensitivity to the deep localized stimulated zone. The measurements can also distinguish between parallel-fracture anisotropic reservoirs and isotropic cases, providing valuable insights into the fracture system of the stimulated zone. Besides the magnetic field measurements, vertical electric field measurements in the open well sections are also highly sensitive to the stimulated zone and can be used as additional data for detecting and imaging the target. We can also acquire additional multiple-source data by grounding the surface electrode at various locations and repeating borehole electromagnetic measurements. This approach can increase the number of monitoring data by several factors, providing a more comprehensive dataset for analysing the deep-localized stimulated zone. The numerical analysis indicates that it is feasible to use the combination of the energized casing and downhole electromagnetic measurements in monitoring localized stimulated zone at large depths.

58 GEOSCIENCES↗
DOE OSTI2026

Optimization of well design and CO 2 injection strategy for risk reduction in Class VI geological carbon sequestration wells

The safety and durability of Class VI wells are critical for geological carbon sequestration (GCS). However, current GCS operations face unique challenges: unlike traditional Class II wells, Class VI CO 2 injection wells operate at rates up to 100 times higher, dramatically increasing the risk of wellbore leakage and structural compromise due to severe temperature drops and associated mechanical stresses. Despite existing guidelines on material selection, there remains a substantial gap in understanding how rapid CO 2 injection rates, low surface temperatures, and variable reservoir conditions interact to threaten long-term well integrity. This study presents a comprehensive, original workflow integrating advanced analytical and numerical models for both well flow and well integrity analysis. By systematically simulating a wide range of field-relevant scenarios—including variations in injection rate, CO 2 temperature, and reservoir pressure—this work provides the first cross-validated assessment of cooling effects on wellbore. The results reveal that extreme temperature drops, up to 60 °C, can occur under high injection rates, particularly in depleted reservoirs, significantly increasing the risk of cement failure. Building on these insights, the study proposes innovative, practical well design and operational strategies, including ductile cement formulations, pre-stressing techniques, advanced insulation coatings, and proactive management of injection rates. The safety of Class VI well extends beyond simply using CO 2 resistant materials. Cement materials should possess optimal thermo-hydraulic-mechanical-chemical properties for effective performance. This work provides a scientific basis for optimizing Class VI well designs, with direct benefits for minimizing environmental risk, lowering operational costs, and enhancing the long-term reliability of GCS.

25 ENERGY STORAGE↗
DOE OSTI2024

Simulation of Multiphase Flow and Poromechanical Effects Around Injection Wells in CO 2 Storage Sites

In geological CO 2 storage operations, wellbore deformations and leakage pathways formations can occur around injection and abandoned wells subjected to high rates and long-term CO 2 injection. To guide engineering design and prevent CO 2 leakage risks, a full understanding of the underlying physics and robust numerical models is necessary to evaluate the response of underground formations in the near wellbore region and in the reservoir. In this study, a multi-scale and multi-physics open-source simulator (GEOS) is used to simulate multiphase flow and poromechanical deformations over time in three dimensions. The governing equations for mechanical deformations of the rock body and multiphase compositional fluid flow within the rock matrix are solved with a fully coupled finite element and finite volume approach. The Drucker–Prager model with friction hardening is applied to simulate elastoplastic deformation and a multiphase fluid model with power-law correlations for relative permeability is used to model the migration of CO 2 plume, which are coupled with numerical implicit scheme. Simulation results are verified against multiple analytical solutions for multiphase flow and wellbore problems, thus demonstrating the accuracy of this advanced simulator. In two engineering applications, here we highlight the impact of elastoplastic deformation and coupled modeling for assessing induced displacements and stress perturbations, which are more pronounced in the near wellbore regions. This work focuses on short-term processes in the vicinity of injection wells where stress evolutions, rock deformations and multiphase compositional flow and transport are simulated jointly to ensure wellbore stability and prevent damage. This fully coupled geomechanical model can simulate multiphase flow and any associated poromechanical effects within the CO 2 storage site and in the surrounding formations. Such a large-scale, long-term, multi-physics simulation model is useful in many ways: it can guide operational decisions for CO 2 injection, assess the containment potential and risks of a site, and analyze the wellbore stability and integrity during and after CO 2 injection.

58 GEOSCIENCES↗