Search NASA⌕ Search

SEARCH · Search NASA

Results for “CO 2 injection”

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 73 records · Page 4
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

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 OSTI2023

Assessing the potential of composite confining systems for secure and long-term CO 2 retention in geosequestration

A potential geologic target for CO 2 storage should ensure secure containment of injected CO 2 . Traditionally, this objective has been achieved by targeting reservoirs with overlying seals-regionally extensive, low permeability units that have been proven capable of retaining buoyant fluid accumulations over geologic time. However, considering that the amount of CO 2 is limited by a decadal injection period, vertical migration of CO 2 can be effectively halted by a composite system of discontinuous shale/silt/mudstone barriers in bedded sedimentary rocks. Here, we studied the impact of depositional architectures in a composite confining system on CO 2 migration and confinement at reservoir scale. We stochastically generated lithologically heterogeneous reservoir models containing discontinuous barriers consistent with statistical distributions of net-sand-to-gross-shale ratio (NTG) and horizontal correlation lengths derived from well log data and observations of producing hydrocarbon fields in Southern Louisiana. We then performed an extensive suite of reservoir simulations of CO 2 injection and post-injection to evaluate the sensitivity of CO 2 plume migration and pressure response of the composite system to a series of geologic and fluid parameters including the lateral continuity of barriers, NTG, permeability anisotropy within the sand body, and capillary pressure contrast between the sand and shale facies. The results indicate that lateral continuity of barriers and NTG are the dominant controls on CO 2 plume geometry and pressure build-up in the reservoir, while the impact of NTG is particularly pronounced. The significance of intraformational barriers becomes apparent as they facilitate the local capillary trapping of CO 2 . Those barriers improve the pore space occupancy by promoting a more dispersed shape of the plume and ultimately retard the buoyancy-driven upward migration of the plume post injection.

58 GEOSCIENCES↗
DOE OSTI2021

Injection data analysis using material balance time for CO 2 storage capacity estimation in deep closed saline aquifers

Estimating the ultimate storage capacity of deep saline aquifers is important to address the formation potential to store the envisioned large volumes of CO 2 . Injection data (i.e. injection rate, bottomhole pressure, and cumulative injected volume of CO 2 ) are routinely recorded during storage operations. These data contain valuable information on the subsurface (e.g. the reservoir pore volume and the formation storage capacity) that can be extracted. In this paper, we present a two-step graphical technique to infer the pore volume and the ultimate storage capacity of closed saline aquifers by analyzing the available injection data. First, the pore volume is inferred through adapting the concept of the material balance time. Material balance time is an approximate superposition time function developed to interpret production data from oil and gas wells operating at variable pressure/rate conditions during the boundary-dominated flow period. Using material balance techniques, the ultimate storage capacity is then estimated through linear extrapolation of the average pressure trend to the maximum allowable pressure the formation can withstand. The average pressure is not available in practice, but is can be obtained from the injection data. Two approaches are presented in this study to calculate the average pressure; namely the rigorous and the approximate approaches. Unlike the rigorous approach, the approximate approach does not require a prior knowledge of some reservoir properties (e.g. relative permeability, absolute permeability, formation porosity and thickness) to calculate the average pressure. To investigate its potential and reliability in analyzing CO 2 injection data, the proposed technique is applied to four synthetic cases representing different well operating conditions. Results indicate that the approximate approach consistently overestimates the actual (simulated) storage capacity as compared to the rigorous approach. The agreement - between the inferred and the simulated reservoir pore volume, and between the analytical and numerical estimates of storage capacity - validates the potential application of the technique to CO 2 storage in closed saline aquifers. The technique is further substantiated through application to a field data set utilized from a commercial-scale geological storage (CGS) project. Finally, field data interpretation shows that the proposed technique can be utilized to identify the degree of hydraulic continuity and reservoir compartmentalization within a target formation by interpreting the corresponding pressure and rate responses.

02 PETROLEUM↗
DOE OSTI2025

Factors controlling injection-induced rupture of intersecting faults during geological sequestration of CO 2

This study addresses coupled multiphase fluid flow and geomechanics effects on potential fault activation associated with subsurface CO 2 injection around intersecting faults. An enhanced fault-representation model is used to capture geomechanical responses of two intersecting faults with finite length during CO 2 injection. The faults are embedded in a strike-slip stress regime of a caprock-reservoir-basement system with the faults represented by zero-thickness interfaces with adjacent finite-thickness damage zones. A sensitivity analysis is conducted to study the effect of fault permeability, slip-weakening behavior, well location relative to the orientation of faults, and well placement (the number and location of injection wells). Five metrics (pressure, CO 2 plume, shear state on the fault, as well as shear displacement and stress path at selected fault monitoring points) are selected to assess CO 2 migration and reactivation of intersecting faults. The results show that induced ruptures are favored by low permeability faults due to high pressure buildup and by slip-weakening behavior resulting from fault strength reduction. The location of one injection well relative to fault orientation determines the magnitude of changes in effective normal stress and shear stress, affecting the location of induced ruptures. Well placement (two injection wells used in the paper) dominates pressure diffusion around the intersection and tips of faults. This redistributes changes in effective normal stress caused by each injection well, influencing the spatial distribution of ruptures along faults. A larger injection volume induces far-field ruptures that are controlled by stress transfer within the injection layer. The findings presented here can provide valuable insights into engineering operations for a long-term, safe, and reliable geologic CO 2 storage.

Fault permeability↗
DOE OSTI2021

Assessing the sequestration time scales of some ocean-based carbon dioxide reduction strategies

Abstract Ocean-based carbon dioxide (CO 2 ) removal (CDR) strategies are an important part of the portfolio of approaches needed to achieve negative greenhouse gas emissions. Many ocean-based CDR strategies rely on injecting CO 2 or organic carbon (that will eventually become CO 2 ) into the ocean interior, or enhancing the ocean’s biological pump. These approaches will not result in permanent sequestration, because ocean currents will eventually return the injected CO 2 back to the surface, where it will be brought into equilibrium with the atmosphere. Here, a model of steady state global ocean circulation and mixing is used to assess the time scales over which CO 2 injected in the ocean interior remains sequestered from the atmosphere. There will be a distribution of sequestration times for any single discharge location due to the infinite number of pathways connecting a location at depth with the sea surface. The resulting probability distribution is highly skewed with a long tail of very long transit times, making mean sequestration times much longer than typical time scales. Deeper discharge locations will sequester purposefully injected CO 2 much longer than shallower ones and median sequestration times are typically decades to centuries, and approach 1000 years in the deep North Pacific. Large differences in sequestration times occur both within and between the major ocean basins, with the Pacific and Indian basins generally having longer sequestration times than the Atlantic and Southern Oceans. Assessments made over a 50 year time horizon illustrates that most of the injected carbon will be retained for injection depths greater than 1000 m, with several geographic exceptions such as the Western North Atlantic. Ocean CDR strategies that increase upper ocean ecosystem productivity with the goal of exporting more carbon to depth will have mainly a short-term influence on atmospheric CO 2 levels because ∼70% will be transported back to the surface ocean within 50 years. The results presented here will help plan appropriate ocean CDR strategies that can help limit climate damage caused by fossil fuel CO 2 emissions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Assessing Impacts on Pressure Stabilization and Leasing Acreage for CO 2 Storage Utilizing Oil Migration Concepts

Favorable geological storage for CO 2 has long been pictured as large anticlines with thick sandstones, similar to oil reservoirs in the petroleum system. Unlike oil, however, stored CO 2 does not need to be recoverable, which raises the possibility of using dissolution and residual trapping to augment the capacity of buoyant traps and tap more of the bulk rock volume. The work presented builds on that idea, asking the following question: If we inject CO 2 down to a syncline – analogous to the carrier bed in the petroleum system – how would this injection mechanism impact storage capacity and plume shape, migration, and stabilization? To address this question, we built a reservoir model, based on seismic interpretation of Middle Miocene strata, offshore Galveston, Texas. 3-D seismic and well logs were used to characterize key intervals. Reservoirs chosen for modeling are progradational-aggradational sands with mud intercalation. They have a higher degree of heterogeneity than the more conventional reservoirs commonly targeted for CO 2 storage. Modeling investigated how far the CO 2 plume would migrate under two scenarios: (1) injecting CO 2 at the base of the salt withdrawal basin (syncline scenario) and (2) injecting CO 2 at the base of the structural closure, similar to a common injection well location for EOR purposes (base scenario). For each scenario, we separately simulated injection of 30 MT of CO 2 and 60 MT of CO 2 continuously for 30 years and observe the plume and pressure evolution for 100 years after the injection stops. The simulation shows that injecting the CO 2 into a syncline limits the vertical migration of CO 2 , thus making synclinal injection more secure. In the syncline scenario, the geological layer around the injection point is more heterogeneous than the layer in the base scenario; thus, the CO 2 tends to migrate laterally. Additionally, in the syncline scenario, the plume does not even reach the upper part of the anticline, allowing us to safely store an additional amount of CO 2 into the reservoir. Furthermore, the simulation also shows that in the syncline scenario, the times needed for the reservoir to reach its stabilized pressure after the end of injections are faster. To summarize, CO 2 injection at the base of a syncline could provide additional storage, increase the safety of the project from the limited vertical plume migration, and expedite plume stabilization, which could result in the decrease of monitoring frequency as the project runs, thus lowering the operating cost of the project in the long run.

58 GEOSCIENCES↗
DOE OSTI2021

Using Controlled-Source Electromagnetic (CSEM) for CO 2 Storage Monitoring in the North Dakota CarbonSAFE Project

Conference presentation for International Meeting for Applied Geophysics & Energy (IMAGE), September 26 – October 1, 2021, Denver, CO. Geophysical methods are key for characterizing the geologic formations to store CO 2 and monitor the injected CO 2 over time to ensure containment. In the integrated multimeasurement geophysical approach considered for this project, it is expected that the controlled-source electromagnetic (CSEM) method is a strong contributor to mapping the CO 2 movement. A feasibility study of the CSEM method, including 1D and 3D modeling and a field noise test, was conducted to determine its effectiveness in monitoring CO 2 in the Broom Creek and Deadwood Formations. The study results demonstrate that the CSEM method can be used for CO 2 storage monitoring in the study area.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2026

Modeling supercritical CO 2 flow and mineralization in reactive host rocks with PFLOTRAN v7.0

Understanding the flow and reactivity of CO 2 injected into geological reservoirs is important for many subsurface applications including secure geologic carbon storage (GCS), critical mineral extraction, enhanced geothermal systems (EGS), and enhanced oil recovery (EOR). Traditionally, subsurface CO 2 injection for GCS applications has focused on geologic formations with favorable subsurface configurations for CO 2 migration and trapping through non-reactive mechanisms such as structural, solubility, and petrophysical trapping. Recently, CO 2 -reactive rocks such as mafic and ultramafic basalts have been investigated for their potential to react with injected CO 2 in situ to simultaneously dissolve host rock minerals and mineralize CO 2 as carbonates. Engineering rapid CO 2 mineralization in the subsurface is attractive because of the increased density of stored CO 2 , the additional safety factors associated with solidification, and the potential to extract valuable critical minerals. Here we present recent developments in the parallel flow and reactive transport simulator PFLOTRAN to model coupled CO 2 -brine flow and reactive transport for a wide range of injection and production applications involving reactive CO 2 -brine systems. These developments are based on the well established and trusted CO 2 flow capabilities in the STOMP-CO 2 simulator. New capabilities added to PFLOTRAN include new CO 2 -brine equations of state with optional thermal coupling, several new constitutive relationships like capillary pressure smoothing and scanning path hysteresis, a fully implicit well model, and native linkage with PFLOTRAN's well-established reactive transport libraries. A series of benchmarks between PFLOTRAN and STOMP-CO 2 verify the newly developed CO 2 -brine flow capabilities. Demonstrations of coupled CO 2 -brine flow modeling and reactive transport show how CO 2 mineralization can be engineered in reactive host rocks. Finally, an example use case involving copper leaching by CO 2 and critical mineral extraction is presented to showcase the strengths of this new implementation. Several limitations still remain, including limited availability of field data to parameterize models. Future work should constrain the evolution of mineral surface area during mineralization and the temperature and/or pH dependence of geochemical reactions for specific systems of interest.

Critical Minerals↗
DOE OSTI2020

Quantification of CO 2 Mineralization at the Wallula Basalt Pilot Project

In 2013, the Pacific Northwest National Laboratory led a geologic carbon sequestration field demonstration where ~1000 tonnes of CO 2 was injected into several deep Columbia River Basalt zones near Wallula, Washington. Rock core samples extracted from the injection zone two years after CO 2 injection revealed nascent carbonate mineralization that was qualitatively consistent with expectations from laboratory experiments and reactive transport modeling. In this study, we report on a new detailed analysis of the 2012 pre-injection and 2015 post-injection hydrologic tests that capitalizes on the difference in fluid properties between scCO 2 and water to assess changes in near-field, wellbore, and reservoir conditions that are apparent approximately two years following the end of injection. This comparative hydrologic test analysis method provides a new way to quantify the amount of injected CO 2 that was mineralized in the field test. Modeling results indicate that approximately 60% of the injected CO 2 was sequestered via mineralization within two years, with the resulting carbonates occupying ~4% of the available reservoir pore space. The method presented here provides a new monitoring tool to assess the fate of CO 2 injected into chemically reactive basalt formations but could also be adapted for long-term monitoring and verification within more traditional subsurface carbon storage reservoirs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Assessment of CO 2 trapping mechanisms in partially depleted oil-bearing sands

The goal of this work is to evaluate various CO 2 sequestration mechanisms occurring in the Morrow B Sandstone in the Farnsworth Unit. A history-matched numerical simulation model was created using extensive geological, petrophysical, and operational data collected from the field. The numerical model is competent to investigate the impact of residual, structural-stratigraphic, solubility, and mineral trapping mechanisms on the fluid transportation dynamics and petrophysical property variations. The model forecasts the field response of 20 years of WAG injections. Afterward, all wells were shut-in, and the reservoir was allowed to evolve for 1000 years to investigate the fate of injected CO 2 . In this paper, we assess the impacts of various trapping mechanisms on oil recovery and CO 2 storage efficacy. By analyzing the results reported from the numerical simulation model, the in-situ fluid composition and mineralogy changes are also investigated. More importantly, we seek to confirm the petrophysical property variations due to the CO 2 injection with observations from laboratory measurements. The experiences gained from this study provide valuable insights regarding physiochemical storage induced by the CO 2 injection activities and serve as a benchmark case for future CO 2 enhanced oil recovery (EOR) projects involving reactive solute transport.

42 ENGINEERING↗
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 OSTI2021

Challenges of designing, operating, and managing brine extraction at a carbon dioxide storage site

Developing and deploying commercial-scale carbon dioxide (CO 2 ) storage at a site requires operational strategies to manage pressure and CO 2 plume distribution. In the natural gas storage and oil and gas industries, fluid extraction has been successfully implemented as a strategy to manage pressure and control fluid distribution (i.e., hydrocarbon recovery and natural gas storage efficiency, respectively). The technical challenges of integrating brine extraction at CO 2 storage sites are expected to be similar, if not identical, to those in the natural gas storage and petroleum industries. Major technical challenges of using brine extraction to manage storage in a storage unit include (1) designing and completing the brine extraction wells, (2) locating the brine extraction well, (3) assessing the efficacy of brine extraction, and (4) handling the extracted brine. The design and completion of the brine extraction well include the well type (i.e., vertical or horizontal) and perforated interval. The location of the well is influenced by the timing of extraction operations (i.e., prior or after CO 2 injection starts) and project objectives (i.e., desired directional movement of the CO 2 plume and avoid or minimize CO 2 extraction). The efficacy of brine extraction requires a baseline projection of the pressure and CO 2 plume distribution. Options for handling extracted brine include subsurface disposal, surface disposal, and industrial use. Here this paper addresses means of overcoming technical challenges that may likely be encountered before or after CO 2 injection starts.

25 ENERGY STORAGE↗
DOE OSTI2023

Legacy Well Leakage Risk Analysis at the Farnsworth Unit Site

This paper summarizes the results of the risk analysis and characterization of the CO 2 and brine leakage potential of Farnsworth Unit (FWU) site wells. The study is part of the U.S. DOE’s National Risk Assessment Partnership (NRAP) program, which aims to quantitatively evaluate long-term environmental risks under conditions of significant geologic uncertainty and variability. To achieve this, NRAP utilizes risk assessment and computational tools specifically designed to quantify uncertainties and calculate the risk associated with geologic carbon dioxide (CO 2 ) sequestration. For this study, we have developed a workflow that utilizes physics-based reservoir simulation results as input to perform leakage calculations using NRAP Tools, specifically NRAP-IAM-CS and RROM-Gen. These tools enable us to conduct leakage risk analysis based on ECLIPSE reservoir simulation results and to characterize wellbore leakage at the Farnsworth Unit Site. We analyze the risk of leakage from both individual wells and the entire field under various wellbore integrity distribution scenarios. The results of the risk analysis for the leakage potential of FWU wells indicate that, when compared to the total amount of CO 2 injected, the highest cemented well integrity distribution scenario (FutureGen high flow rate) exhibits approximately 0.01% cumulative CO 2 leakage for a 25-year CO 2 injection duration at the end of a 50-year post-injection monitoring period. In contrast, the highest possible leakage scenario (open well) shows approximately 0.1% cumulative CO 2 leakage over the same time frame.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Machine learning and deep learning for mineralogy interpretation and CO 2 saturation estimation in geological carbon Storage: A case study in the Illinois Basin

Carbon capture and storage (CCS) is a promising approach to simultaneously maintaining energy security and reducing carbon dioxide (CO 2 ) emissions under the current energy portfolio that is dominated by fossil fuel energy. Pre-injection formation characterization and post-injection CO 2 monitoring are two critical tasks to guarantee storage efficiency in CCS. The CCS projects in the Illinois Basin, the first large-scale CO 2 injection into saline aquifers in the United States, employed conventional and the latest pulsed neutron logging (PNL) tools for mineralogy interpretation and CO 2 saturation estimation, which provide valuable references for future CCS projects. Because of the inherent fuzziness of petrophysical measurements and complex subsurface heterogeneity, interpreting well-logging data is time-consuming, and its accuracy can be user-biased. In recent years, data-driven methods have been widely used to capture the non-linear patterns between input features and interpretation results. This work applied and evaluated four commonly used machine learning (ML) models, including ridge regression (RR), random forest (RF), gradient boosting regression (GBR), support vector regression (SVR), and one deep learning (DL) model, the artificial neural network (ANN). We optimized the hyperparameters of the four ML models and the DL model using the simulated annealing algorithm and the grid search strategy, respectively. The input features of the mineralogy interpretation models were eleven conventional well-logging parameters, and the label data (i.e., ground truth) were the porosity and volumetric fractions of six minerals, including quartz, feldspar, dolomite, calcite, clay, and iron minerals. The results demonstrated that the GBR and RF models were superior in predicting volumetric fractions of minerals and porosity; label data with low coefficient of variation (CV) values tended to yield better performance. For CO 2 saturation estimation, the RF was the best-performing model, followed by SVR, ANN, GBR, and RR. Furthermore, we conducted feature importance ranking using the permutation importance algorithm and found that the formation sigma and well pressure were the most important features in this study. In conclusion, the study of CCS projects in the Illinois Basin bridges the gap between the limited knowledge and understanding of geological carbon storage and the increasing demand for reliable, cost-effective, and sustainable energy solutions.

58 GEOSCIENCES↗
DOE OSTI2021

NRAP-Open-IAM Multisegmented Wellbore Reduced-Order Model

Geologic carbon storage is one of the promising strategies to mitigate climate change by reducing the emission of carbon dioxide to the atmosphere. As part of the National Risk Assessment Partnership (NRAP), a systems-level stochastic analysis tool called the open source integrated assessment model, NRAP-Open-IAM, has been developed to estimate and manage the risk of containment loss at a geological carbon sequestration site. NRAP-Open-IAM contains several wellbore leakage model components that estimate the fluid leak rate that may occur through compromised legacy wells due to the increase in pressure resulting from CO 2 injection activities. Coupled to a reservoir component model, these components estimate the leakage of CO 2 and/or brine from a storage reservoir to overlying aquifer layers and the atmosphere through legacy wells. This report presents the theoretical framework and quality testing of the multisegmented wellbore reduced-order model. The model allows for segmenting of the legacy wells passing through the overlying stratigraphy into several intervals to simulate a site’s specific stratigraphic and hydrogeologic properties. For quality assurance, the analytical model is validated against numerical reservoir flow simulations for single and multiple aquifer(s) models. The results indicate that the model accurately predicts the transport of two-phase fluids (brine and injected CO 2 ) through the well over time. A detailed description of the model helps users to understand the model and provides a basis for future improvements.

58 GEOSCIENCES↗
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↗
DOE OSTI2023

Selecting representative geological realizations to model subsurface CO 2 storage under uncertainty

Carbon capture and storage (CCS) is one of the quickest and most effective solutions for reducing carbon emissions. The majority of subsurface storage occurs in saline aquifers, for which geological information is lacking which in turn results in geological uncertainty. To evaluate uncertainty in CO 2 injection projections, the use of multiple geological realizations (GRs) has been practiced very commonly. In this approach, hundreds or thousands of high-resolution GRs is used that quickly becomes computationally expensive. This issue can be addressed with representative geological realizations (RGRs) that preserve the uncertainty domain of the ensemble GRs. Here, in this study, we propose the use of unsupervised machine learning (UML) frameworks, including dissimilarity measurement, dimensionality reduction, clustering and sampling algorithms ta select a predetermined number of RGRs. We compare the simulation outputs of the RGR sets and the ensemble using the Kolmogorov–Smirnov (KS) test to select the best UML. The UML frameworks and their associated selection processes are evaluated using a saline aquifer with a single CO 2 injection well and 200 GRs with varying uncertain petrophysical characteristics. The best UML framework is selected to use only 5% of the GRs while maintaining the uncertainty domain of the ensemble GRs. In addition, the best UML framework is tested using a saline aquifer with three CO 2 injection wells and varied GRs. The results show that our proposed UML framework can be used to choose RGRs, capturing the whole uncertainty domain. Our approach leads to a significant reduction in the computational cost associated with scenario testing, decision-making, and development planning for CO 2 storage sites under geological uncertainty.

58 GEOSCIENCES↗