Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “subsurface storage”

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 235 records · Page 13

Pore-Scale Study on the Positive Feedback Between Stress and Porosity Caused by Pressure Solution in Porous Media

Pressure solution is an important process in the evolution of sedimentary rocks, which provide storage space for most of our petroleum resources. It directly influences the generation, migration, and storage of petroleum fluids in subsurface sedimentary rocks. Here, in this paper, we develop a pore-scale, mechanochemical model to demonstrate a possible positive feedback between the local porosity and pore surface stress, in which a higher local porosity causes a higher local pore surface stress, thus enhancing pressure solution and consequently further increasing the local porosity. Pore surface stress represents stress on a solid grain adjacent to a pore. Specifically, the pore-scale, mechanochemical model directly simulates the stress distribution over solid and pore surfaces using a finite element model. The dissolution of solids at the solid-pore interfaces under a far-from-equilibrium condition is simulated using a first-order kinetics model that accounts for the local stress distribution. The updated pore geometry, caused by pore surface dissolution, is then used in the stress simulation in the next numerical iteration. Two types of porous media, the Oriskany sandstone and an artificial porous medium with spherical pores, were tested in the mechanochemical simulation. The positive stress-porosity feedback during pressure solution was observed in both samples. In addition, the model quantitatively illustrated the distribution of local mineral dissolution rates on all pore surfaces, as well as its relation to the effective mineral dissolution rate of the entire sample. Based on the comparison between the two porous media, the local mineral dissolution was regulated by pore space distribution, geometry, and coalescence during pressure solution. This work is the first that uses direct, pore-scale numerical simulation to demonstrate the positive stress-porosity feedback during pressure solution, which has the potential to advance the understanding of the mechanical-chemical (MC) coupling in many geological processes that are relevant to subsurface energy systems, such as the recovery of petroleum hydrocarbons and geothermal energy.

CT scanning↗

High-Resolution 3D Acoustic Borehole Integrity Monitoring (Research Performance Progress Report Q2 2021)

Real-time, in-situ, high spatial resolution (sub-cm) imaging of the near-borehole environment would revolutionize wellbore diagnostics and integrity assessment by direct observation of defects. It is becoming increasingly apparent, that better understanding of the near-wellbore environment is required to meet the safety and operational needs in challenging environments such as those present in subsurface energy extraction (geothermal) and storage (CO2 sequestration) applications. Therefore, it is important to have a more robust ability to image the near-borehole and reliably detect defects. It was proposed to further develop and improve our advanced 3D imaging system to evaluate casing defects (e.g. corrosion) and cement quality in either open- or cased-borehole with the ultimate goal to develop a commercially deployable technology. The system consists of a unique acoustic source (LANL) and advanced inversion techniques for image processing (LANL, ORNL). This system will provide comprehensive borehole integrity monitoring with improved resolution over existing techniques. As an application of this imaging system, we will characterize the effectiveness of next-generation wellbore completion technology (NETL, SNL), and will demonstrate that, unlike current technology, the proposed approach can successfully characterize foamed cements.

15 GEOTHERMAL ENERGY↗

Report Q3FY21 - High-Resolution 3D Acoustic Borehole Integrity Monitoring

Real-time, in-situ, high spatial resolution (sub-cm) imaging of the near-borehole environment would revolutionize wellbore diagnostics and integrity assessment by direct observation of defects. It is becoming increasingly apparent, that better understanding of the near-wellbore environment is required to meet the safety and operational needs in challenging environments such as those present in subsurface energy extraction (geothermal) and storage (CO 2 sequestration) applications. Therefore, it is important to have a more robust ability to image the near-borehole and reliably detect defects. It was proposed to further develop and improve our advanced 3D imaging system to evaluate casing defects (e.g. corrosion) and cement quality in either open- or cased-borehole with the ultimate goal to develop a commercially deployable technology. The system consists of a unique acoustic source (LANL) and advanced inversion techniques for image processing (LANL, ORNL). This system will provide comprehensive borehole integrity monitoring with improved resolution over existing techniques. As an application of this imaging system, we will characterize the effectiveness of next-generation wellbore completion technology (NETL, SNL), and will demonstrate that, unlike current technology, the proposed approach can successfully characterize foamed cements.

42 ENGINEERING↗

High-Resolution 3D Acoustic Borehole Integrity Monitoring (Final Report)

Real-time, in-situ, high spatial resolution (sub-cm) imaging of the near-borehole environment would revolutionize wellbore diagnostics and integrity assessment by direct observation of defects. It is becoming increasingly apparent, that better understanding of the near-wellbore environment is required to meet the safety and operational needs in challenging environments such as those present in subsurface energy extraction (geothermal) and storage (CO2 sequestration) applications. Therefore, it is important to have a more robust ability to image the near-borehole and reliably detect defects. It was proposed to further develop and improve our advanced 3D imaging system to evaluate casing defects (e.g. corrosion) and cement quality in either open- or cased-borehole with the ultimate goal to develop a commercially deployable technology. The system consists of a unique acoustic source (LANL) and advanced inversion techniques for image processing (LANL, ORNL). This system will provide comprehensive borehole integrity monitoring with improved resolution over existing techniques. As an application of this imaging system, we will characterize the effectiveness of next-generation wellbore completion technology (NETL, SNL), and will demonstrate that, unlike current technology, the proposed approach can successfully characterize foamed cements.

42 ENGINEERING↗

A microfluidic study of transient flow states in permeable media using fluorescent particle image velocimetry

Velocity fields in flow in permeable media are of great importance to many subsurface processes such as geologic storage of CO 2 , oil and gas extraction, and geothermal systems. Steady-state flow is characterized by velocity fields that do not change significantly over time. The flow field transitions to a new steady state once it experiences a disturbance such as a change in flow rate or in pressure gradient. This transition is often assumed to be instantaneous, which justifies the expression of constitutive relations as functions of instantaneous phase saturations. This work examines the evolution of velocity fields in a surrogate quasi-2D permeable medium using a microfluidic device, a microscopy system, and a high-speed camera. Tracer particles are injected into the medium along with DI water. The evolution of the velocity field is examined by tracing these particles in the captured images using the standard high-density particle image velocimetry algorithm founded on cross-correlation. The results suggest that the transition between steady states for an incompressible fluid takes a finite and non-negligible amount of time that is independent of the magnitude of the change in pressure gradient. The existence of transient states and the nature of the response during these states are readily interpreted by the principle of least action where flow gradually establishes an optimal configuration such that energy dissipation is 2 minimized. The findings provide evidence against the applicability of the assumption that flowing phases relax instantaneously to their steady states and, hence, against the accuracy of the classical multiphase extension of Darcy's law.

42 ENGINEERING↗

Technoeconomic Analysis of Infrastructure Buildout Scenarios

The success of CCUS deployment in the southeast will depend on the optimized pipeline network from CO 2 point sources to geologic sinks for CO 2 storage. An optimal pipeline framework will reduce both environmental impacts and pipeline building cost. However, finding an optimal pipeline/transport infrastructure design is a non-trivial task and requires simultaneous usage of large volumes of data, computational resources, and a state-of-the-art simulator. Such a transport infrastructure model must consider point sources for CO 2 capture and associated volumes, sinks for CO 2 storage, and transportation from source to sink via pipeline networks. These considerations must be addressed simultaneously and systemically in an optimization process, in which a defined objective function (i.e., capital, variable CO 2 capture, transport, and storage cost function) is required to be minimized with the consideration of practical constraints (e.g., CO 2 flow through pipelines is less than the maximum capacity). Although optimization has been widely used in subsurface resources production and CO 2 storage, its application in large scale CCUS infrastructure design is rarely reported in the literature. SimCCS, developed by Los Alamos National Laboratory, is an open-source CCUS pipeline infrastructure design toolset that facilitates the optimization of pipeline infrastructure networks.

99 GENERAL AND MISCELLANEOUS↗

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

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

54 ENVIRONMENTAL SCIENCES↗

Modeling the Cost of Onshore CO2 Pipeline Transport and Onshore CO2 Saline Storage

This paper describes the FECM/NETL CO2 Transport Cost Model (CO2_T_COM), a technoeconomic model of CO2 transport by pipeline, and the FECM/NETL CO2 Saline Storage Cost Model (CO2_S_COM), a technoeconomic model of storage of CO2 in a deep, subsurface saline formation. The results of applying CO2_T_COM to calculate break-even CO2 prices to transport CO2 at different mass flow rates and different distances are presented. Similarly, CO2_S_COM is used to calculate the break-even CO2 price for storing CO2 in 314 potential storage formations across the US. These break-even prices are used to construct cost-supply curves for CO2 storage which are presented on a national and regional basis. CO2_T_COM and CO2_S_COM are the most comprehensive open source technoeconomic models available for analyzing CO2 pipeline transport and CO2 saline storage. To be presented at the SPE?AAPG/SEG Carbon Capture Utilization and Storage Conference in Houston, TX, March 11-13, 2024.

Morgan, David↗

Predicting Geologic Behavior in Carbon Storage Projects Using Graph Neural Network

This study was invited to presented at NVIDIA's GTC conference to highlight the potential of Graph Neural Network as a novel and promising methodology for predicting pressure and saturation evolution in carbon storage projects. Carbon capture and storage (CCS) technology plays a pivotal role in mitigating greenhouse gas emissions, facilitating the transition to a low-carbon future. Effective management of subsurface reservoirs is essential to ensure the safe and efficient storage of captured carbon dioxide (CO₂). Accurate predictions of pressure and saturation over time are critical for evaluating the long-term performance and integrity of CCS projects. In recent years, Graph Neural Network (GNN) has emerged as a powerful framework for analyzing complex data in graph-structured domains. This abstract explores the application of GNN to forecast pressure and saturation evolution in carbon storage projects. Traditional numerical simulations of subsurface reservoirs have proven successful in providing pressure and saturation forecasts. However, these simulations involve massive amounts of computational effort and require extensive domain expertise for proper model calibration and validation. Graph Neural Operator offers an alternative approach that harnesses the inherent graph structure of reservoirs, where nodes represent reservoir grid cells and edges represent the geological connectivity between them.

Shih, Chung Yan↗

Seismic Elastic Double-Beam Characterization of Faults and Fractures for CO₂ Storage Site Selection

Site characterization for underground injection and storage of gigatonne-scale CO₂ requires reliable and cost-effective methods to detect and characterize faults and fractures and to assess their stress state and fault activation potential. This is critical, as wastewater injection and disposal have been shown to activate faults and induce earthquakes, and CO₂ leakage remains a key concern for long-term storage. In this project, we developed seismic methods to detect and characterize large-scale sedimentary and crystalline basement faults and associated small-scale fractures below conventional seismic imaging resolution using multicomponent (9C) surface seismic data. Machine learning was used to automatically interpret large-scale faults, providing key information for estimating the maximum magnitude of potential induced earthquakes. High-fidelity imaging was achieved by exploiting redundancy across multiple elastic wave modes, where independent images from different modes and frequencies cross-validate each other. We also used our nonlinear signal comparison (NLSC) method for ground roll removal, improving data quality in complex near-surface conditions. The methods were validated using field data acquired in central Montana. Results show that basement faults extend into the sedimentary section and that small-scale fractures are widespread above the basement. The inferred stress orientation is consistent with regional stress data, and the estimated maximum induced earthquake magnitude is small (Mw ~2.3). The developed workflow provides a practical approach for fault and fracture characterization and for assessing induced seismicity and leakage risk. It is directly applicable to CO₂ storage site selection and to other subsurface systems.

02 PETROLEUM↗

Design, construction, and validation of an in-situ groundwater trace element analyzer with applications in carbon storage

Abstract It is estimated that carbon emissions should reach net-zero by 2050 to meet important climate targets. Carbon capture is likely necessary to reach these targets, requiring a long-term storage solution such as geological carbon sequestration. However, as with any subsurface activity, leakage can occur, potentially impacting groundwater quality near the storage site. Rapid detection is essential to mitigate damage to this resource. Since CO 2 will acidify groundwater, the concentrations of acid soluble minerals and associated cations will increase. Thus, an in-situ, real-time element analysis system based on laser-induced breakdown spectroscopy (LIBS) is under development to monitor these elements. The system splits the traditional LIBS system into a miniature, all-optical sensor head built around a passively Q-switch laser fiber coupled to a control unit. Previous work has validated the LIBS technique for use at high pressure as well as the split system design. In this work, a fieldable prototype sensor is developed and tested in an onsite monitoring well where trace elements concentrations (approx. 0–3 ppm) were tracked over 20 days. These concentrations varied in response to local rainfall, diluting with increased rain, demonstrating the ability of a LIBS-based sensor to track trace elements under real-world conditions.

54 ENVIRONMENTAL SCIENCES↗

Basic Energy Sciences Roundtable: Foundational Science for Carbon Dioxide Removal Technologies (Brochure)

Basic research to identify and understand the fundamental principles governing carbon dioxide removal processes of capture, conversion, and storage is essential for achieving zero carbon emissions by 2050. In March 2022, the Office of Basic Energy Sciences (BES)—in coordination with the US Department of Energy (DOE) technology Offices of Energy Efficiency and Renewable Energy, Fossil Energy and Carbon Management, and Nuclear Energy—held a roundtable titled, “Foundational Science for Carbon Dioxide Removal Technologies,” to discuss the scientific and technical barriers for CO 2 capture, conversion, and storage. Five priority research opportunities were identified to address these scientific and technical challenges and accelerate progress toward the realization of zero carbon emissions. The PROs are as follows: 1) Master Interfacial Processes of CO2 Transport and Reactivity Across Multiple Length and Time Scales; 2) Create Materials that Simultaneously Exhibit Multiple Properties for CO 2 Capture and Release or Conversion; 3) Discover Unconventional Pathways for Energy-Efficient CO 2 Capture, Release, and Conversion; 4) Control Multiphase Interactions Required for CO 2 Conversion into Minerals and Materials; and 5) Achieve Predictive Understanding of Coupled Processes in Complex Subsurface Geologic Systems for Secure Carbon Storage.

36 MATERIALS SCIENCE↗

A Project Lifetime Approach to the Management of Induced Seismicity Risk at Geologic Carbon Storage Sites

The geologic storage of carbon dioxide (CO 2 ) is one method that can help reduce atmospheric CO 2 by sequestering it into the subsurface. Large-scale deployment of geologic carbon storage, however, may be accompanied by induced seismicity. We present a project lifetime approach to address the induced seismicity risk at these geologic storage sites. This approach encompasses both technical and nontechnical stakeholder issues related to induced seismicity and spans the time period from the initial consideration phase to postclosure. These recommendations are envisioned to serve as general guidelines, setting expectations for operators, regulators, and the public. They contain a set of seven actionable focus areas, the purpose of which are to deal proactively with induced seismicity issues. Although each geologic carbon storage site will be unique and will require a custom approach, these general best practice recommendations can be used as a starting point to any site-specific plan for how to systematically evaluate, communicate about, and mitigate induced seismicity at a particular reservoir.

58 GEOSCIENCES↗

Impacts of Topography‐Driven Water Redistribution on Terrestrial Water Storage Change in California Through Ecosystem Responses

Abstract Lateral subsurface flow plays an essential role in sustaining the terrestrial ecosystem, but it is not explicitly represented in most Earth System Models. In this study, we implemented an explicit lateral saturated flow model into the E3SM land model (ELM). The model explicitly describes lateral flow in the saturated zone by representing, for each model grid, an idealized hillslope consisting of five hydrologically connected soil columns. We conducted three model experiments driven by 0.125° atmospheric forcing data during 1980–2015 over California using models of the default ELM, a modified version of ELM to enhance infiltration, and the model with the lateral saturated flow model. The simulated runoff, evapotranspiration, and terrestrial water storage anomaly (TWSA) from the three simulations were evaluated against available observations, and the model explicitly representing lateral flow performs best. The new model produces greater gridcell‐averaged evapotranspiration especially over the mountainous regions with moderate relief and seasonally dry climates. Most importantly, it improves the modeled seasonal variations, interannual variabilities, and the recent decadal decline of TWSA. Many of these improvements can be attributed to the enhanced ecosystem resilience to droughts as demonstrated by transpiration increases caused by lateral flow. Model sensitivity experiments suggest that subsurface runoff is most sensitive to the ratio between horizontal and vertical saturated hydraulic conductivity, followed by hillslope planforms (convergent, divergent, and uniform), number of columns, and lower boundary conditions. Future work should effectively characterize hillslopes in global models and explore the long‐term influences of lateral water movement on modeled biogeochemical cycle.

54 ENVIRONMENTAL SCIENCES↗

The Kimberlina synthetic multiphysics dataset for CO 2 monitoring investigations

Abstract We present a synthetic multi‐scale, multi‐physics dataset constructed from the Kimberlina 1.2 CO 2 reservoir model based on a potential CO 2 storage site in the Southern San Joaquin Basin of California. Among 300 models, one selected reservoir‐simulation scenario produces hydrologic‐state models at the onset and after 20 years of CO 2 injection. Subsequently, these models were transformed into geophysical properties, including P‐ and S‐wave seismic velocities, saturated density where the saturating fluid can be a combination of brine and supercritical CO 2 , and electrical resistivity using established empirical petrophysical relationships. From these 3D distributions of geophysical properties, we have generated synthetic time‐lapse seismic, gravity and electromagnetic responses with acquisition geometries that mimic realistic monitoring surveys and are achievable in actual field situations. We have also created a series of synthetic well logs of CO 2 saturation, acoustic velocity, density and induction resistivity in the injection well and three monitoring wells. These were constructed by combining the low‐frequency trend of the geophysical models with the high‐frequency variations of actual well logs collected at the potential storage site. In addition, to better calibrate our datasets, measurements of permeability and pore connectivity have been made on cores of Vedder Sandstone, which forms the primary reservoir unit. These measurements provide the range of scales in the otherwise synthetic dataset to be as close to a real‐world situation as possible. This dataset consisting of the reservoir models, geophysical models, simulated time‐lapse geophysical responses and well logs forms a multi‐scale, multi‐physics testbed for designing and testing geophysical CO 2 monitoring systems as well as for imaging and characterization algorithms. The suite of numerical models and data have been made publicly available for downloading on the National Energy Technology Laboratory's (NETL) Energy Data Exchange (EDX) website.

58 GEOSCIENCES↗

Impacts on mechanical strength of chemical reactions induced by hydrous supercritical CO2 in Boise Sandstone

Geomechanics experiments were used to assess mechanical alteration of Boise Sandstone promoted by reactions with supercritical carbon dioxide (scCO 2 ) and water vapor. During geologic carbon storage, scCO 2 is injected into subsurface reservoirs, forming buoyant plumes. At brine-plume interfaces, scCO 2 can dissolve into native brines, and water from brines can partition into scCO 2 , forming hydrous scCO 2 . This study investigates the effect of hydrous scCO 2 on the strength of Boise Sandstone. Samples are first exposed to recirculating hydrous scCO 2 for 24 h at 70 °C and 13.8 MPa scCO 2 pressure. Samples are reacted with scCO 2 with added water contents up to 500 mL. After scCO 2 exposure, samples are deformed at room temperature under confining pressures of 3.4, 6.9, and 10.3 MPa. The results demonstrate that hydrous scCO 2 induces chemical reactions in Boise Sandstone, with ions migrating from the solid into the hydrous scCO 2 phase. At the longer time-scales, these reactions could lead to mechanical weakening in the samples; however, on the scale of our experiments, the strength changes are within sample variability. Because the solubility of water in scCO 2 is extremely low (0.008 mol H 2 O per 1 mol CO 2 ), the mineral dissolution of Boise Sandstone was under 0.002 wt.%. Additionally, mineral grains and pore throats in Boise Sandstone are cemented with quartz, which is not susceptible to dissolution at these conditions. Our results indicate that humidity in scCO 2 plumes is unlikely to sustain chemical reactions and induce long term strength changes in quartz cemented sandstones due to resistant mineralogies and low water solubility.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A robust framework for frictional fault contact in geological formations using a stabilized augmented Lagrangian approach

Numerical simulations are essential to evaluate the performance and safety of engineered subsurface systems such as geological carbon storage sites, enhanced geothermal fields, and oil and gas reservoirs. A key challenge lies in accurately modeling the frictional contact behavior along fault surfaces. This problem involves inequality constraints that arise from the physics of frictional slip, requiring specialized numerical methods to handle the resulting highly nonlinear and path-dependent behavior. Here, in this work, we address this challenge using an Augmented Lagrangian Method (ALM) implemented via the Uzawa algorithm. The formulation employs mixed finite element spaces, combining low-order piecewise linear displacements within the 3D domain cells with piecewise constant tractions defined on the fault surfaces. Furthermore, to ensure stability and satisfy the inf-sup condition, the discrete displacement space is enriched with face bubble functions on both sides of the contact interfaces. This approach offers several advantages over other stabilization techniques that rely on additional terms, and it integrates naturally in the Uzawa framework.

58 GEOSCIENCES↗

Time-lapse seismic data inversion for estimating reservoir parameters using deep learning

Geologic carbon sequestration involves the injection of captured carbon dioxide ([Formula: see text]) into subsurface formations for long-term storage. The movement and fate of the injected [Formula: see text] plume is of great concern to regulators because monitoring helps to identify potential leakage zones and determines the possibility of safe long-term storage. To address this concern, we design a deep-learning framework for [Formula: see text] saturation monitoring to determine the geologic controls on the storage of the injected [Formula: see text]. We use different combinations of porosities and permeabilities for a given reservoir to generate saturation and velocity models. We train the deep-learning model with a few time-lapse seismic images and their corresponding changes in saturation values for a particular [Formula: see text] injection site. The deep-learning model learns the mapping from the change in the time-lapse seismic response to the change in [Formula: see text] saturation during the training phase. We then apply the trained model to data sets comprising different time-lapse seismic image slices (corresponding to different time instances) generated using different porosity and permeability distributions that are not part of the training to estimate the [Formula: see text] saturation values along with the plume extent. Our algorithm provides a deep-learning assisted framework for the direct estimation of [Formula: see text] saturation values and plume migration in heterogeneous formations using the time-lapse seismic data. Our method improves the efficiency of time-lapse inversion by streamlining the large number of intermediate steps in the conventional time-lapse inversion workflow. This method also helps to incorporate the geologic uncertainty for a given reservoir by accounting for the statistical distribution of porosity and permeability during the training phase. Tests on different examples verify the effectiveness of our approach.

Geochemistry & Geophysics↗