Results for “Legacy wells” — Page 2

DOE OSTI2026

Hydrogeological assessment of CO2 containment assurance and wellbore integrity at a Gulf Coast storage site

Abstract A large-scale carbon capture and storage (CCS) initiative on the Texas Gulf Coast serves as a premier demonstration of the U.S. Department of Energy’s CarbonSAFE program. Targeting deep saline formations, specifically Oligo-Miocene deltaic sequences, the project aims to establish technical and commercial viability for geologic CO2 storage within a major industrial corridor. This study provides a rigorous hydrogeological assessment to support Class VI permitting by quantifying the high degree of containment security. Utilizing a compositional reservoir simulator, we developed a suite of 27 distinct simulation cases to evaluate vertical plume dynamics near both planned injection wells and proximal legacy infrastructure. To ensure numerical accuracy near wellbores, we implemented a refined mesh strategy, determining that a 5.6 ft × 5.6 ft grid refinement offered the optimal balance between computational efficiency and descriptive precision. The modeling framework utilized a systematic sensitivity-based approach to evaluate the mechanical redundancy of the subsurface system by performing a bounding analysis of wellbore interfaces against hypothetical high-permeability microannuli. By systematically isolating competing physical drivers, including permeability, porosity, gas hysteresis, thermal gradients, salinity, and solubility trapping (quantified via Henry’s law with dynamically adjusted coefficients), this work moves beyond binary assessments to establish a nuanced hierarchy of containment factors. The results confirm that primary trapping mechanisms (e.g., gas hysteresis and solubility), combined with the site's unique geomechanical stratigraphy, significantly restrict vertical mobility and reinforce the robust containment security of the reservoir. Baseline results demonstrate substantial vertical separation between the CO2 plume and the upper confining system, ensuring robust containment. Sensitivity analysis reveals that even under highly conservative bounding scenarios—assuming theoretical 10-Darcy pathways at specific wellbore locations—the 2,900-ft thick multi-layered confining zone remains a reliable barrier. In these hypothetical upper-bound cases, peak upward fluxes of CO2 and saltwater after 15 years of injection remain localized and dissipate rapidly within the lower sections of the confining interval, leaving the integrity of the seal uncompromised. Furthermore, the study identifies that while localized wellbore pathways define theoretical upper bounds of vertical migration, the Area of Review (AoR) is primarily sensitive to regional thermal gradients and hysteresis, which can influence the AoR by over 3,000 acres in pessimistic configurations. Also, primary trapping mechanisms, specifically gas hysteresis and solubility, work in tandem with the Gulf Coast’s unique geomechanical stratigraphy to significantly restrict vertical mobility. Ductile, smectite-rich mudstones facilitate natural borehole convergence and the self-healing of potential conduits, creating a natural geomechanical bridge that effectively mitigates migration potential at both current injection points and legacy-well locations. This comprehensive modeling effort demonstrates that the integration of high-resolution wellbore simulations and regional geomechanical observations confirms the long-term storage security of the studied site, providing a physics-based foundation for industrial-scale CCS deployments. This modeling framework establishes a baseline for future research into coupled geomechanical effects, such as time-dependent borehole convergence, to further refine long-term containment projections. Acknowledgements We thank the Gulf Coast Carbon Center (GCCC) at the Bureau of Economic Geology for foundational research support. We appreciate Alex Bump for technical guidance and David Hoffman for model mesh generation. This work used TACC’s Frontera cluster for simulations and CMG Ltd. software licenses provided to UT-Austin. This material is based upon work supported by the Department of Energy under Award Number DE-FE0032338. Disclaimer This material is based upon work supported by the U.S. Department of Energy’s Fossil Energy and Carbon Management Office under the CarbonSAFE program, award Number DE-FE0032338. The views expressed herein do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

58 GEOSCIENCES↗

DOE OSTI2026

Hydrogeological assessment of CO2 containment assurance and wellbore integrity at a Gulf Coast storage site

Abstract A large-scale carbon capture and storage (CCS) initiative on the Texas Gulf Coast serves as a premier demonstration of the U.S. Department of Energy’s CarbonSAFE program. Targeting deep saline formations, specifically Oligo-Miocene deltaic sequences, the project aims to establish technical and commercial viability for geologic CO2 storage within a major industrial corridor. This study provides a rigorous hydrogeological assessment to support Class VI permitting by quantifying the high degree of containment security. Utilizing a compositional reservoir simulator, we developed a suite of 27 distinct simulation cases to evaluate vertical plume dynamics near both planned injection wells and proximal legacy infrastructure. To ensure numerical accuracy near wellbores, we implemented a refined mesh strategy, determining that a 5.6 ft × 5.6 ft grid refinement offered the optimal balance between computational efficiency and descriptive precision. The modeling framework utilized a systematic sensitivity-based approach to evaluate the mechanical redundancy of the subsurface system by performing a bounding analysis of wellbore interfaces against hypothetical high-permeability microannuli. By systematically isolating competing physical drivers, including permeability, porosity, gas hysteresis, thermal gradients, salinity, and solubility trapping (quantified via Henry’s law with dynamically adjusted coefficients), this work moves beyond binary assessments to establish a nuanced hierarchy of containment factors. The results confirm that primary trapping mechanisms (e.g., gas hysteresis and solubility), combined with the site's unique geomechanical stratigraphy, significantly restrict vertical mobility and reinforce the robust containment security of the reservoir. Baseline results demonstrate substantial vertical separation between the CO2 plume and the upper confining system, ensuring robust containment. Sensitivity analysis reveals that even under highly conservative bounding scenarios—assuming theoretical 10-Darcy pathways at specific wellbore locations—the 2,900-ft thick multi-layered confining zone remains a reliable barrier. In these hypothetical upper-bound cases, peak upward fluxes of CO2 and saltwater after 15 years of injection remain localized and dissipate rapidly within the lower sections of the confining interval, leaving the integrity of the seal uncompromised. Furthermore, the study identifies that while localized wellbore pathways define theoretical upper bounds of vertical migration, the Area of Review (AoR) is primarily sensitive to regional thermal gradients and hysteresis, which can influence the AoR by over 3,000 acres in pessimistic configurations. Also, primary trapping mechanisms, specifically gas hysteresis and solubility, work in tandem with the Gulf Coast’s unique geomechanical stratigraphy to significantly restrict vertical mobility. Ductile, smectite-rich mudstones facilitate natural borehole convergence and the self-healing of potential conduits, creating a natural geomechanical bridge that effectively mitigates migration potential at both current injection points and legacy-well locations. This comprehensive modeling effort demonstrates that the integration of high-resolution wellbore simulations and regional geomechanical observations confirms the long-term storage security of the studied site, providing a physics-based foundation for industrial-scale CCS deployments. This modeling framework establishes a baseline for future research into coupled geomechanical effects, such as time-dependent borehole convergence, to further refine long-term containment projections. Acknowledgements We thank the Gulf Coast Carbon Center (GCCC) at the Bureau of Economic Geology for foundational research support. We appreciate Alex Bump for technical guidance and David Hoffman for model mesh generation. This work used TACC’s Frontera cluster for simulations and CMG Ltd. software licenses provided to UT-Austin. This material is based upon work supported by the Department of Energy under Award Number DE-FE0032338. Disclaimer This material is based upon work supported by the U.S. Department of Energy’s Fossil Energy and Carbon Management Office under the CarbonSAFE program, award Number DE-FE0032338. The views expressed herein do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

58 GEOSCIENCES↗

DOE OSTI2023

Integrating risk assessment methods for carbon storage: A case study for the quest carbon capture and storage facility

This paper documents how tools developed by the National Risk Assessment Partnership (NRAP) complement the Bowtie risk assessment approach through the demonstration of a risk-based Area of Review (AoR) determination for the Quest Carbon Capture and Storage (CCS) facility near Edmonton, Alberta with the NRAP Open-Source Integrated Assessment Model (NRAP Open-IAM). The Bowtie risk assessment for the Quest CCS site was developed by Shell Canada Ltd., who operates the asset on behalf of the Athabasca Oil Sand Project — a joint venture, owned by Canadian Natural Resources Ltd, Chevron Canada and Shell Canada Ltd. The Quest Bowtie treats site risk comprehensively – considering a large set of technical and project risks (e.g., operational risks, above ground infrastructure, public perception, etc.). NRAP tools, methods, and analyses are focused on quantitative assessment of subsurface technical risks. NRAP-Open-IAM was used to evaluate potential groundwater impacts at the Quest site, i.e., a sufficient change above background groundwater quality levels to be measurable, should CO 2 and/or brine leak from the reservoir via a faulty and/or compromised wellbore. Our results show there is a low potential for groundwater impact and support a further reduction (from ~460 to 100 km 2 ) to the established AoR for the Quest site. Integrating both approaches early in the site characterization period could significantly benefit commercial-scale projects by drastically reducing the region over which leakage risks (e.g., legacy wells) need to be evaluated, i.e., the AoR. This is especially true for brownfield sites, which could contain hundreds to thousands of legacy wells. Further, basing monitoring network design on those regions where impacts are most likely to occur has the potential to facilitate more efficient deployment of monitoring technologies and improve the likelihood of detecting a leak should one occur.

58 GEOSCIENCES↗

DOE OSTI2021

A Review of Well Integrity Based on Field Experience at Carbon Utilization and Storage Sites

Maintaining well integrity is critical to the success of geologic carbon storage (GCS) and carbon dioxide enhanced oil recovery (CO 2 -EOR) operations. Wells that experience leakage because of integrity issues can potentially become a risk to the environment or human health if they release previously captured CO 2 back into the atmosphere or into freshwater aquifers. There are many GCS and CO 2 -EOR sites in operation around the world. However, well integrity experiences at these sites are not widely documented in the public domain. This study details findings from a survey of well integrity experiences elicited from operators of GCS and CO 2 -EOR sites. The goal of the survey was to obtain information about site characteristics and operator experiences with well integrity, monitoring methods, and risk assessment of legacy wells. Literature relevant to the survey questions was also reviewed and summarized to provide context for survey responses and identify areas where field experiences with well integrity do and do not align with the current state of research. We highlight the current state-of-practice, identify research needs, and provide context for future interactions between researchers, operators, and regulators on issues related to well integrity.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

A Review of International Field Experience with Well Integrity at Carbon Utilization and Storage Sites

This report details findings from a survey of well integrity experiences elicited from operators of geologic carbon storage (GCS) and carbon dioxide enhanced oil recovery (CO2-EOR) sites around the world. The survey consisted of 41 questions organized in four sections and its goal was to obtain information about site characteristics and operator experiences with well integrity, monitoring methods, and risk assessment of legacy wells. Current literature relevant to the survey questions was also reviewed and summarized to provide context for survey responses and identify areas where field experiences with well integrity do and do not align with the current state of research.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

NETL Well Integrity Workshop: Identifying Well Integrity Research Needs for Subsurface Energy Infrastructure

Wells are a critical component of subsurface energy infrastructure. Ensuring the integrity of wells as engineered pathways for the safe extraction, injection, and storage of fluids in the subsurface is key to maximizing the effectiveness and resilience of that infrastructure. Addressing well integrity issues in a technically robust manner that promotes environmental sustainability and social equity is also an important focus of the United States (U.S.) Department of Energy’s (DOE) Office of Fossil Energy and Carbon Management. Industry best practices, regulatory standards, modern monitoring data acquisition and control systems, and decades of research and development have dramatically improved the performance and reliability of wells for hydrocarbon extraction and underground injection in the oil and gas industry. Yet, important innovation is required to improve and ensure well integrity performance in engineered geologic systems where operational environments (fluid composition, temperature, pressure, and/or stress conditions) and long functional life cycles of well systems present unique challenges. Additionally, work is needed to understand and manage the long-term integrity and risks associated with legacy wells—especially those located adjacent to and presenting hazards for new subsurface activity.

02 PETROLEUM↗

DOE OSTI2020

Effect of Computational Schemes on Coupled Flow and Geo-Mechanical Modeling of CO 2 Leakage through a Compromised Well

Carbon capture, utilization, and storage (CCUS) describes a set of technically viable processes to separate carbon dioxide (CO 2 ) from industrial byproduct streams and inject it into deep geologic formations for long-term storage. Legacy wells located within the spatial domain of new injection and production activities represent potential pathways for fluids (i.e., CO 2 and aqueous phase) to leak through compromised components (e.g., through fractures or micro-annulus pathways). The finite element (FE) method is a well-established numerical approach to simulate the coupling between multi-phase fluid flow and solid phase deformation interactions that occur in a compromised well system. We assumed the spatial domain consists of a three-phases system: a solid, liquid, and gas phase. For flow in the two fluids phases, we considered two sets of primary variables: the first considering capillary pressure and gas pressure (PP) scheme, and the second considering liquid pressure and gas saturation (PS) scheme. Fluid phases were coupled with the solid phase using the full coupling (i.e., monolithic coupling) and iterative coupling (i.e., sequential coupling) approaches. The challenge of achieving numerical stability in the coupled formulation in heterogeneous media was addressed using the mass lumping and the upwinding techniques. Numerical results were compared with three benchmark problems to assess the performance of coupled FE solutions: 1D Terzaghi’s consolidation, Liakopoulos experiments, and the Kueper and Frind experiments. We found good agreement between our results and the three benchmark problems. For the Kueper and Frind test, the PP scheme successfully captured the observed experimental response of the non-aqueous phase infiltration, in contrast to the PS scheme. These exercises demonstrate the importance of fluid phase primary variable selection for heterogeneous porous media. We then applied the developed model to the hypothetical case of leakage along a compromised well representing a heterogeneous media. Considering the mass lumping and the upwinding techniques, both the monotonic and the sequential coupling provided identical results, but mass lumping was needed to avoid numerical instabilities in the sequential coupling. Additionally, in the monolithic coupling, the magnitude of primary variables in the coupled solution without mass lumping and the upwinding is higher, which is essential for the risk-based analyses.

deformation flow↗

DOE OSTI2020

Identifying Abandoned Well Sites Using Database Records and Aeromagnetic Surveys

Oil and natural gas are primary sources of energy in the United States. Improved drilling and extracting techniques have led to a renewed interest in historic oil and gas fields, but limited records of legacy wells make new drilling efforts more difficult, as abandoned wells may provide conduits for liquids and gases to migrate into groundwater reservoirs or the atmosphere. Well finding using aeromagnetic surveys pinpoints the location of steel-cased wells, detecting both active and abandoned wells, including buried casings lacking aboveground markers. Here, we present six aeromagnetic surveys conducted in Pennsylvania and Wyoming as case studies, comparing the magnetic points to locations known in databases. In all study sites, more magnetic points were detected than recorded in databases. Based on differences between theoretical database well counts and the actual number of wells detected in surveys, we estimated the total number of wells in Pennsylvania to be 395 000–466 000 and 181 000–182 000 in Wyoming. Extrapolating to the national level, in this work, we estimate the average number of wells in the continental United States is 6.04 ± 19.97 million wells with 1.16 ± 3.84 million of those designated as abandoned wells, within the range of previous abandoned well count estimations. Although aeromagnetic surveys are limited to detecting steel-cased wells and do not differentiate sites based on well status, this study nevertheless demonstrates the utility of aeromagnetic surveys in well finding efforts across the country and shows limitations in database records of oil and natural gas wells.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Assessing Suitable Geologic Carbon Storage Sites Across Utah

Utah has a wealth of potential geological reservoirs for carbon dioxide storage (CS) and a long history of geologic research resulting in an abundance of available subsurface data to evaluate CS potential. Reservoirs may include sandstone, carbonate, and basalt; these rock types are plentiful in Utah’s subsurface and the complex Phanerozoic history throughout the state requires evaluating each geologic region individually for promising reservoir-seal pairs for CO2 storage. Classifying Utah by geologic provinces (or “geo-regions”) allows for customized thinking about suitable CS reservoir and seal distribution, CO2 point sources, land use, and existing infrastructure. Preliminary results from this study highlight the geologic CS potential across 15 geo-regions. Four regions stand out as having high CS potential: the Uinta Basin, San Rafael Swell, Paradox Basin, and the southern Basin and Range Province. The Uinta Basin and San Rafael Swell geo-regions are well suited for CS and have several projects ongoing to evaluate Cretaceous Frontier and Naturita Formations, Jurassic Navajo Sandstone and Entrada Sandstone, and Permian Weber Sandstone reservoir units that lie beneath robust sealing units like the ~5000-ft-thick Mancos Shale and Carmel Formation. Reservoirs such as the Navajo and Weber Sandstones have been demonstrated to be suitable reservoirs through a long history of oil and gas exploration in Utah. New areas of interest include the southern Basin and Range in southwest Utah, where the Jurassic Navajo Sandstone is overlain by the sealing Carmel Formation at suitable depths (>3000 ft), and has good porosities based on outcrop analogue data. Just to the north (e.g., central Basin and Range), legacy wells and 2D seismic data show possible salt and subsurface basalt flows that may provide additional possible CS reservoirs and seals. This geo-region also has the advantage of being coupled with geothermal energy resources that may be used to power burgeoning direct air capture technologies. In the Paradox Basin of southeastern Utah, the Leadville Limestone is a potential storage reservoir beneath the thick (4000–5000 ft), salt-bearing Pennsylvanian Paradox Formation. Although the northern and western parts of Utah offer CS potential, these areas typically contain less infrastructure and subsurface penetrations, creating geologic uncertainty associated with subsurface seals and reservoirs due to a lack of data. Overall, this statewide assessment and ranking is the first step to aid in evaluating CS potential across Utah and provides a foundation for future research in the most favorable locations.

58 GEOSCIENCES↗

DOE OSTI2022

The Cypress Sandstone Seal System

The Cypress Sandstone is the youngest and shallowest unit in the Illinois Basin that was featured in the United States Carbon Utilization and Storage Atlas IV as a target for saline carbon storage with an estimated 0.2 to 2.3 GT of storage potential. Additional research on a residual oil zone (ROZ) developed within the Cypress Sandstone has delineated 27 prospects with approximately 290.8 million m3 (1.8 billion barrels) of oil in place. 21 to 31 million m 3 (144 to 196 million barrels) of oil is estimated to be recoverable using carbon dioxide enhanced oil recovery (CO 2 -EOR). Storage of CO 2 associated with EOR in these ROZ prospects alone, not accounting for associated main pay zones (MPZs), underlying brine formation, or intervals adjacent to or between prospects, is estimated to be up to 10.4 billion tonnes. The Cypress Sandstone is thus well understood to be a CO 2 injection target, both for EOR and associated storage. However, the seal system overlying the Cypress Sandstone is poorly understood. Unlike deeper formations such as the Mt. Simon Sandstone or the St. Peter Sandstone which are either in use as a CO 2 sink or being characterized for prospective storage, respectively, the Cypress is not overlain by hundreds of feet of impermeable shale. Rather, the Cypress is overlain by a lithologically variable interval that is composed generally of shales and limestones with some sandstone in the part of the Basin where the Cypress is deep enough to facilitate CO 2 storage. Also, due to its status as one of the shallowest and most prolific oil reservoirs in the Basin, the seal system overlying the Cypress Sandstone has a relatively high number of legacy well penetrations. The purpose of this report is to characterize the Cypress Sandstone seal system using well logs and available core. Gross thickness, lithology (facies), and mineralogy of seals is described and mapped across the Basin. The column height of CO 2 that can be held is calculated using capillary pressure data from a representative core.

02 PETROLEUM↗

DOE OSTI2026

A Review and Outlook on Experimental Advances and Innovations in Geological CO 2 Storage: Insights from Depleted Gas Reservoirs and Saline Aquifers

Geological storage of carbon dioxide (CO 2 ) in depleted gas reservoirs and deep saline aquifers is a key part of global decarbonization efforts. As carbon capture and storage advances toward commercial-scale deployment, the credibility and scalability of laboratory experiments are increasingly vital for guiding safe and effective field implementation. This review offers a comprehensive, cross-scale evaluation of experimental methodologies, including core flooding, high-pressure, high-temperature systems, microfluidic visualization, and emerging systems such as multilayer commingled/compartmentalized core flooding, 3D-printed micromodels, and AI-powered digital twins. These innovations are demonstrated to enhance representativeness, reproducibility, and real-time insight, thereby addressing the limitations of conventional workflows. A critical analysis of methodological gaps, such as inconsistent pressure–temperature conditions, oversimplified brine chemistry, and a lack of standardization, reveals experimental sources of scale translation errors and performance uncertainty. By comparing the unique challenges of depleted gas reservoirs (such as low water saturation and legacy well leakage) to those of saline aquifers (including pressure buildup and caprock integrity), this review identifies formation-specific priorities for experimental design. Novel contributions include a synthesis of best practices, integration strategies for model calibration, and recommendations for standardizing core handling, saturation procedures, and reporting protocols. Furthermore, this work serves as a guide for developing robust, field-relevant experimental strategies that can increase the deployment and regulatory acceptance of CO 2 storage technologies at scale.

58 GEOSCIENCES↗

DOE OSTI2020

LLNL Kimberlina 1.2 NUFT Simulations June 2018 (v2)

This dataset contains the output 6,000, 3-dimensional reactive multi-phase flow and transport aquifer simulations of brine and CO2 leakage into a protective aquiver in California’s San Joaquin Valley and input data files detailing the geologic mesh, aquifer physical properties and CO2 and brine injection rates. This data set was generated as an ongoing effort with the US DOE National Risk Assessment Partnership (NRAP) to evaluate the effectiveness of monitoring techniques to detect brine and CO2 leakage from legacy wells into underground sources of drinking water overlaying a CO2 storage reservoir. Each simulation contains a unique set of input parameters, generated stochastically. The outputs consist of these upper three geologic layers (from top): the Etchegoin, Macoma-Chanac, Santa Margarita-McLure formations. These simulations span the several distances (1, 3 and 6 km or wells W31-0.2, W31-0.5 and W31-1.0, respectively) from the CO2 injector, initiated from bottom hole pressure and saturation to calculate wellbore leakage from the storage reservoir, with low and high regional groundwater gradients and wellbore leakage into 5 leaky nodes. The dataset includes 1,000 unique simulations for each distance, which each contain a unique aquifer heterogeneity, aquifer and caprock permeability, and two model generations are included with a high permeability (prod07) and hybrid permeability (prod09). The range of permeability distributions is listed in Table 1. Each model generation consists of 3,000 simulations. Included in the dataset are the leakage rates determined from 2D wellbore models which utilize the pressure and CO2 saturation from LBL's reservoir simulations, NUFT mesh files with distributed lithology, NUFT rocktab files which describe the material properties for the geologic layers and the NUFT input files and post-processed output 'ntab' files. Each ntab file contains spatial (rows) and temporal (columns) model output tables for each model cell, the locations (x,y,z) and dimensions for each cells (dx, dy, dz). Table 1. Permeability distribution ranges for prod07 and prod09 model generations Geologic Layer: Permeability Range (log10 m^2) prod07 prod09 Etchegoin -12.92 to -10.92 -13.70 to -11.44 Macoma-Chanac -12.72 to -10.72 -13.50 to -11.24 Santa Margarita-McLure -12.70 to -10.70 -13.48 to -11.22 The input files used to generate the model include which are included in the dataset are: Time series of CO2 leakage input into the model (ex: Q_brn.W31-0.2.sim1000.layers123.tab) Time series of CO2 leakage input into the model (ex: Q_CO2.W31-0.2.sim1000.layers123.tab) Physical properties of the aquifer materials detailing the aquifer porosity, solid density, partitioning coefficients, permeabilities and van-Genuchten parameters detailed in a NUFT rocktab file: (ex: sim1000.usnt.rocktab) Numerical mesh and geologic data assigned to each model cell detailed in a NUFT genmsh format (ex: sim1000.mesh_k16.prod07.trans.genmsh) The primary output parameters are: pH (use absolute value) Change in TDS (mg/kg) Change in Pressure (Pa) Change CO2 gas saturation (fraction range 0.0-1.0) for example, the directory /p/lscratchh/mansoor1/nrap/kimberlina/prod09/mainfiles/sim1000/W31- 0.2 contains: sim1000.W31-0.2.trans.pH.red.ntab sim1000.W31-0.2.no_bg.trans.TDS.red.ntab sim1000.W31-0.2.usnt.P.deltabg.red.ntab sim1000.W31-0.2.usnt.CO2_sat.deltabg.red.ntab Each row in the NTAB files consist of model output per numerical grid cell. Each output file contains 33 columns (variables), including the information of numerical records, geologic location and sizes and the simulated parameter values over time. The first 13 variables are about numerical records and relative geologic information for a simulation grid: 1. index: simulation index 2. i: the ith grid of x-axis 3. j: the ith grid of y-axis 4. k: the ith grid of z-axis 5. element_ref: element reference 6. nuft_ind: nuft index 7. x: grid location in the x axis direction 8. y: grid location in the y axis direction 9. z: grid location in the z axis direction 10. dx: grid length in the x axis direction 11. dy: grid length in the y axis direction 12. dz: grid length in the z axis direction 13. volume: volume of the simulation grid The remainder (14, 15, 16...) variables are the simulated parameter values over time, take Pressure as an example, are: 14. 0.0y: initial pressure per cell. 15. 10.0y: simulated pressure at the end of the 10th year. 16. 20.0y: simulated pressure at the end of the 20th year. ... (repeated for every 10 years until 200 years)... The model extends 10,000 m, 5,000 m and 1,411 m in the x,y and z dimensions, respectively. The mesh consists of 164,832 cells with mesh dimensions of 101 x 51 x 32 (nx, ny, nz), with cell dimensions ranging from 100 m laterally (along x and y-axis) and model layers are as designated in the z-axis: Layer 1: atmosphere (1e-30 m thick) Layer 2: upper caprock (10 m thick) Layers 3-13: Etchegoin (536.23 m thck) Layers 14-27: Macoma-Chanac (679.04 m thick) Layers 28-32: Santa Margarita-McLure (185.94 m thick) The wellbore is placed along node i=51, j=26, and extends vertically along 5 nodes from the top to the bottom of the model. Special instructions when extracting files: Each Gzip archive (ex: prod07.sim1000-sim00099.tar.gz) contains 100 simulations. Gzip archives should be transferred into base directories (ie. In Linux: mkdir prod07; mv prod07.*.tar.gz prod07/.) before extracting, or files will be overwritten. Each sub-simulation tree should have the following file structure pattern (using the linux 'tree' command): |-- prod07 | |-- sim0001 | |-- W31-0.2 | | |-- Q_brn.W31-0.2.sim0001.layers123.tab | | |-- Q_co2.W31-0.2.sim0001.layers123.tab | | |-- sim0001.W31-0.2.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.2.trans.pH.red.ntab | | |-- sim0001.W31-0.2.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.2.usnt.P.deltabg.red.ntab | |-- W31-0.5 | | |-- Q_brn.W31-0.5.sim0001.layers123.tab | | |-- Q_co2.W31-0.5.sim0001.layers123.tab | | |-- sim0001.W31-0.5.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.5.trans.pH.red.ntab | | |-- sim0001.W31-0.5.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.5.usnt.P.deltabg.red.ntab | |-- W31-1.0 | | |-- Q_brn.W31-1.0.sim0001.layers123.tab | | |-- Q_co2.W31-1.0.sim0001.layers123.tab | | |-- sim0001.W31-1.0.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-1.0.trans.pH.red.ntab | | |-- sim0001.W31-1.0.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-1.0.usnt.P.deltabg.red.ntab | |-- sim0001.mesh_k16.prod07.trans.genmsh Disclaimer This document was prepared as an account of work sponsored by an agency of the United States government. Neither the United States government nor Lawrence Livermore National Security, LLC, nor any of their employees makes any warranty, expressed or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States government or Lawrence Livermore National Security, LLC. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States government or Lawrence Livermore National Security, LLC, and shall not be used for advertising or product endorsement purposes. Lawrence Livermore National Laboratory is operated by Lawrence Livermore National Security, LLC, for the U.S. Department of Energy, National Nuclear Security Administration under Contract DE-AC52-07NA27344. This report was reviewed and released as LLNL-MI-753464.

aquifer↗

DOE OSTI2022

NRAP-Open-IAM: Generic Aquifer Component Development and Testing

The Generic Aquifer Model calculates the concentrations of dissolved salt and dissolved CO 2 surrounding a leaking legacy well. The Generic Aquifer model can also estimate the size of an “impact plume” where concentration changes exceed user-specified thresholds. The model is a component of NRAP-Open-IAM, an open-source Integrated Assessment Model (IAM) developed by the National Risk Assessment Partnership (NRAP) to perform risk assessment for geologic CO 2 storage. The input parameters were selected to cover a wide range of groundwater aquifers and leakage rates. The generic aquifer model was developed using a generative adversarial deep learning network, trained using a large synthetic dataset of STOMP multiphase flow simulations. The deep learning model predictions of dissolved salt and dissolved CO 2 in the aquifer compare well to the original STOMP simulation results. The extent of aquifer impacted by leaking CO 2 or brine is calculated using a user-defined mass fraction threshold. The aquifer impact volumes calculated based on STOMP simulation results compare well to those calculated based on the deep learning model. In a provided python script, gridded observation results from the generic aquifer component of NRAP-Open-IAM are converted to HDF5 format files for monitoring design with the DREAM code.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

NRAP-Open-IAM Analytical Reservoir Model: Development and Testing

Geological carbon sequestration (GCS) is a key technology for reducing global carbon dioxide (CO 2 ) emissions. Over the last decade, the U.S. Department of Energy has invested in understanding the science base, developing practical implementation methods, and demonstrating secure GCS technologies to mitigate the environmental impacts associated with the atmospheric release of CO 2 . As part of the National Risk Assessment Partnership, a systems-level risk assessment tool, called the NRAP-Open-IAM, has been developed to conduct risk assessment and enable safe operations at a GCS site. The current NRAP-Open-IAM contains a simple reservoir model component that calculates the evolution of CO 2 saturation and fluid pressure in a storage reservoir during CO 2 injection operations. This report presents the development and testing of a new analytical reservoir reduced-order model (ROM), which is extended from an existing semi-analytical model for estimation of CO 2 and brine leakage along legacy wells, and enhances the capability of the NRAP-Open-IAM to simulate more types of reservoir conditions. The developed model is validated against three reference studies, and the results indicate that the new ROM predicts the behavior of the two-phase fluids (brine and injected CO 2 ) well and is applicable to different reservoir simulation boundary conditions (i.e., constant pressure boundary and infinite-acting boundary) without a priori user specification of the boundary type. Sensitivity analysis for a set of model parameters is performed using 4,000 synthetic cases prepared via a fully automated process and using machine-learning-based feature selection. The stochastic analysis identifies gravitational number (i.e., ratio of gravitational forces to viscous force) and distance between the injection well and observation location as the most impactful parameters for matching the pressure and CO 2 saturation, respectively, between the numerical simulations and the ROM. This report details the possible ROM uncertainties and serves as a guide for users to understand the use and limitations of this ROM. The code implementation of the model will be released as a module within the NRAP-Open-IAM.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

NRAP-Open-IAM Multisegmented Wellbore Reduced-Order Model: Improvement and Quality Assurance

The multisegmented wellbore model (MSW) semi-analytically estimates the amount of CO 2 and brine leakage from a leaking legacy well by segmenting it into intervals to simulate site-specific stratigraphic and hydrogeologic properties. The model is a component of the National Risk Assessment Partnership Open-Source Integrated Assessment Model (NRAP-Open-IAM), which was developed to perform risk assessment for geologic CO 2 storage. The new wellbore leakage model, which uses deep learning networks for a caprock segment, was developed to enhance the analytical MSW. The model was trained and validated using a synthetic data set of Subsurface Transport Over Multiple Phases (STOMP) multiphase flow simulations from various geological, well attribute, and operational conditions to ensure its quality. The results demonstrate that the model is more accurate than the existing model in predicting the transport of two-phase fluids (brine and injected CO 2 ) through the well. This report provides a detailed explanation of the model development and quality assurance.

58 GEOSCIENCES↗

DOE OSTI2022

H-2-SALT: Storing Fossil Energy as Hydrogen in Salt Caverns

Hydrogen storage in subsurface caverns is known as perhaps largest and longest-term energy storage system. While few such systems are currently in service (e.g., Texas; Teesside, UK), they are envisioned as being a pillar of the current energy transition that will allow intermittent power production from renewable sources (e.g., wind, solar) to be balanced with demand at grid-scale. In addition, hydrogen energy storage can allow existing fossil power plants to run more economically either by minimizing startups/shutdowns or by allowing them to take advantage of arbitrage. Finally, stored hydrogen can feed a variety of non-power users such as pipeline gas, heating, transportation, and manufacturing. The H-2-SALT paper study assessed the feasibility of a power-to-hydrogen system utilizing salt cavern storage of hydrogen in bedded salt in central Kansas, where over 750 such caverns have been constructed to date, of which approximately 350 are still in service. Legacy well and salt cavern data were collected to develop a regional geological database for two sites in Central Kansas and adjacent areas for use in more detailed geological and geomechanical site characterization efforts. This study found commercial viability and competitiveness of large-scale, electrolytic H 2 production and storage ($1.78 per kg H 2 ) that can be used for both electrical power supply to the grid during high-priced electricity periods (including co-firing with NG and fuel cell power production) and sale of H 2 for various industries, such as petrochemicals or transportation. The commercial potential of the H-2-SALT system is also relatively high because each of its components operate commercially today. It is recommended that a stratigraphic test well be drilled to collect log data and core samples that can reduce uncertainty surrounding geological and geomechanical properties of the salt beds that will be required for further cavern and energy storage system design studies.

03 NATURAL GAS↗

DOE OSTI2024

Jurassic deserts and CO2 sequestration: the Glen Canyon Group in Utah, western USA

The Early Jurassic in the Western Interior of USA consisted of an extensive desert environment, including one of the largest ergs in geologic history. The Navajo erg has been estimated to extend as much as 2.2 million square kilometers, though the preserved extent is somewhat smaller. The Mesozoic was a global greenhouse phase, and during the Jurassic this region experienced fluctuations in climate aridity, reflected in the depositional environment. In addition to extensive dunefields, interdune deposits (lakes and oases) and fluvial systems are also documented within formations of the Glen Canyon Group. The Navajo Sandstone has received much attention as a potential CO2 injection target in recent years. It consists of thick, aeolian sandstones with high porosity and permeability, occurs in both outcrop and subcrop, and has industry data. Mapping of stratigraphic and sedimentologic changes within the Navajo Sandstone has been undertaken in localized areas, however piecing together these studies and developing regional models for CO2 potential has not received as much attention. There is significant industry data across Utah, and utilizing this data to pivot from a focus on hydrocarbon extraction to CO2 injection is an effective way to move forward with green energy, and to meet carbon neutral emission goals. Using legacy well data, we are developing a more comprehensive study of the Navajo Sandstone as a potential CO2 reservoir, in addition to identifying other zones of interest within the Glen Canyon Group, and improving understanding of stratigraphic complexity within one of the most significant aeolian systems in the world.

Mahon, Elizabeth↗

DOE OSTI2025

TCCSP Project Implementation Plan

The TCCSP CarbonSAFE Phase II efforts focused on developing the technical, community, and economic foundation for a commercial-scale regional geologic storage complex for CO 2 captured. TCCSP has conducted a thorough evaluation of the regional, local, and site-specific geology, legacy well infrastructure, injection site design, and project planning as part of the CarbonSAFE Phase II TCCSP. Specially, the work performed under the CarbonSAFE Phase II has significantly furthered the advancement towards the commerciality of the TCCSP storage complex by delivering on the following objectives: Objective 1: Acquire site-specific characterization data to validate carbon dioxide (CO 2 ) storage estimates of reservoir units and the competency of confining layers in the region to prevent CO 2 leakage to develop a commercial-scale CO 2 project. Objective 2: Assess technical, economic, regulatory, and stakeholder aspects of the TCCSP.

54 ENVIRONMENTAL SCIENCES↗