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 217 records · Page 12

CO2 Saline Storage Talk at the University of Wyoming

These are slides presented virtually to a class on Carbon Capture and Storage (CCS) at the University of Wyoming in Laramie WY. The talk was on the basics of carbon storage in saline formations and reviewed basic geology, important physical properties of CO2 and brine, how fluids flow in the subsurface, environmental risks associated with CO2 storage, the regulations governing CO2 saline storage, and high-level design considerations. The FECM/NETL CO2 Saline Storage Cost Model was used to illustrate how geologic properties vary in different regions of the US (PA, IL and WY) and how the performance of CO2 storage and the cost of CO2 CO2 saline storage depends on geology (October 23, 2024).

Morgan, David↗

An analysis of Potential Benefits of Adding Diurnal Thermal Storage in Geothermal Heat Pump Systems

Geothermal heat pumps (GHP) are highly efficient for space heating and cooling, utilizing the subsurface of the ground as seasonal energy storage through borehole heat exchangers (i.e., heat stored in summer is used to provide heating in winter, and cold stored in winter is used to provide cooling in summer). On the other hand, diurnal thermal energy storage can actively be charged and discharged for load shifting. Conventional thermal energy storage uses tanks filled with thermal storage media (such as phase change materials) and heat exchangers. It occupies building floor space or land areas outside the building, which may not be available especially in residential buildings. This paper will introduce a novel configuration that can utilize existing GHP system for both diurnal and seasonal thermal storage. This system can actively charge the ground on daily basis by injecting heat or cold in borehole heat exchangers. A simulation model was created for the novel configuration. Computer simulations were conducted using the model for the new GHP system serving a typical single-family house in three U.S. cities representing hot, warm, and cold climates. Based on simulation results, the potential benefits of load shifting and downsize of the expensive ground heat exchangers used by the GHP system were assessed.

Anees, Fady [ORNL]↗

Minimizing exposure to legacy wells and avoiding conflict between storage projects: Exploring area of review as a screening tool

Elevated pressure from large-volume injection is a key driver of risk and project cost. If transmissive features (e.g., non-isolating wells or fracture systems) are present, increased injection-zone pressure can drive fluids from depth toward protected freshwater resources. In US Carbon Capture and Storage (CCS) law, the area at risk is known as the Area of Review (AoR). The size and number of potentially transmissive features to be evaluated and possibly remediated or managed is a function of the size and location of the AoR. The size of the AoR depends on several variables, including properties of the injection zone, properties of protected resources, and injection rate and duration. Evaluation of the intersection of these variables across a portfolio of sites highlights the injection zone depth and boundary conditions as top-level controls. Deep injection, use of multiple stacked injection zones, reduced injection rate and choice of injection well location can all be used to minimize AoR and the number of potentially transmissive features within it. Here, we introduce the concept of pressure space (defined as connected pore volume times pressure) as the key subsurface commodity for CO 2 storage and we suggest that it forms a more robust basis for leasing and regulation than pore space alone.

58 GEOSCIENCES↗

Nitrogen Vacancy Centers in Diamond for Stress Sensing Applications: Theory and Experiments

Sensing of stress under high elevated environmental conditions with high resolution has critical importance for range of applications including earth’s subsurface scanning and geological CO2 storage monitoring. Using first principles density functional theory (DFT) approach combined with the theoretical modelling of low energy Hamiltonian, here we investigate a novel approach to detect unprecedented level of pressure by taking advantage of solid-state electron’s spin of Nitrogen vacancy (NV) centers in nanodiamond. We computationally explore the effect of strain on the defect band edges and band gaps by varying the lattice parameters of diamond supercell hosting a single NV center. A low energy Hamiltonian is developed that includes effect of stress on the energy level of ±1 spin manifold at the ground state. By quantifying the energy level shift and split, we predict the pressure sensing of up to 0.3 MPa/√Hz. We show superiority of the quantum sensing over traditional optical sensing techniques by discussing our results from DFT and theoretical modelling for the frequency shift per unit pressure. The proposed quantum stress sensing could be useful to measure earth’s subsurface vibrations. Our results open avenues for development of sensing technology with a high sensitivity and resolution under extreme pressure conditions that potentially has a wider applicability than existing pressure sensing technologies.

Paudel, Hari P.↗

Causes of Missing Snowmelt Following Drought

Abstract Water management in snowy mountainous regions hinges on forecasting snowmelt runoff. However, droughts are altering snowpack‐runoff relationships with ongoing debate about the driving mechanisms. For example, in 2021 in California, less than half of predicted streamflow arrived. Mechanisms proposed for this “missing” streamflow included changes in evapotranspiration ( ET ), rainfall, and subsurface moisture conditions. Here, we demonstrate that ET in drought years generates dry subsurface conditions that reduce runoff in subsequent years. A model including this legacy of depleted moisture storage reduced median error in 2021 forecasts from 60% to 20% at 15 minimally disturbed basins and from 18% to 2% at 6 water supply basins in the Sierra Nevada (basins range in area from 5 to 23,051 km 2 and mean annual precipitation from 814 to 1,549 mm). Our findings indicate that the relationship between snowpack and runoff will evolve as plant ecosystems respond to climate change and alter subsurface water storage dynamics.

Lapides, Dana A.↗

Numerical investigation of multiphase flow through self-affine rough fractures

Multiphase flow through fractures has great significance in subsurface energy recovery and gas storage applications. Different fracture and flow properties affect flow through a fracture which is difficult to control in laboratory experiments. Here, we perform lattice Boltzmann simulations in an ensemble of synthetically generated fractures. Drainage simulations are performed at different capillary numbers, wettability, and viscosity ratios. We track the invading front and quantify breakthrough saturations and show that roughness and wettability have a strong effect on fluid invasion through a complex fracture. Invading a more viscous fluid results in more stable displacement regardless of the capillary number while at very low capillary numbers, fluid migration is dependent on the inherent structure of the fracture. We develop a fluid displacement phase diagram in a single rough fracture and compare our results from that in the literature. Finally, we extend the phase diagrams across multiple fractures and demonstrate the importance of natural fracture features of roughness and wettability in identifying stable versus unstable displacement regimes during multiphase flow through rough fractures. Our work presents an end-to-end numerical pathway for testing on experimental data and expanding numerical data sets for testing combinations of different physical phenomenon and make valuable predictions on fluid flow through rough fractures.

02 PETROLEUM↗

Proxy-based Bayesian inversion of strain tensor data measured during well tests

Recent instrument developments have made it possible to measure the strain tensor caused by injecting or pumping fluid from aquifers or reservoirs, but the full value of these data is limited because the long runtimes of poroelastic forward models makes it impractical to use many inversion schemes. This limits the interpretation of strain data for managing the recovery of resources or storage of wastes in the subsurface. This paper describes a method of inverting deformation data using a poroelastic numerical simulator so the results can be used to manage reservoirs or aquifers. We developed a workflow designed to reduce the number of simulations sufficiently to make it feasible to use DREAMzs, an advanced Bayesian inversion method that translates the uncertainties from different sources into unbiased posterior parameter distributions and uncertainty envelopes around the field data. Using a KNN proxy model for the poroelastic simulator is key to reducing the overall computations, and the workflow includes a strategy for ensuring the proxy model results converge on the results from the simulator. The workflow is tested using an idealized example that verifies the ability to correctly identify parameters and characterize noise used to perturb the data. Field data from an injection test at an oil reservoir near Tulsa, Oklahoma, are also used to evaluate the efficacy of the workflow with a real dataset. The workflow identified 265 history matching solutions out of 1240 total simulation runs (21% acceptance ratio), where the results were used to characterize posterior parameter distribution and evaluate the prediction uncertainty. Furthermore, this workflow is significant because it enables strain tensor, or other geomechanical measurements to be interpreted to guide decision-making during energy and environmental processes in the subsurface.

42 ENGINEERING↗

Insights into Controls of Mineralogy and Pore Structure on the Density of Methane Adsorption Phase in Shales under Supercritical Conditions

With gradually increasing enhancing population and industrialization, exploitation of shale gas is enhancing in the world to satisfy the growing demand for energy worldwide. To exploit shale gas efficiently, accurate assessment of shale gas in place (GIP) is necessary for determining production strategies. Unlike free and dissolved gas, adsorbed gas contributes to the shale GIP up to 85% due to the well-developed surface area and micropores. Measured excess adsorption of methane (up to over 95% of subsurface shale gas) is often corrected to absolute adsorption to obtain the actual amount of adsorption under supercritical (or geological) conditions. During the correction, adsorption phase density (APD) is critical. However, the APD of subsurface shale gas and the effects of shale properties (e.g., mineralogy and pore structure) on APD remain poorly understood. A series of high-pressure methane adsorption isotherms (HPMAI) on Caney Shales were collected and analyzed in conjunction with other United States and Chinese shales from the literature at 35–125 °C and up to 15 MPa. Here, a three-layer Ono–Kondo (OK3) model is utilized to derive the temperature- and pressure-dependent APD coupling low-pressure nitrogen adsorption isotherms (LPNAI) and HPMAI. X-ray diffraction and organic geochemistry are combined to reveal the mineralogy. Brunauer–Emmett–Teller, Barrett–Joyner–Halenda, and Horvath–Kawazoe analyses via LPNAI are used to investigate the pore structures. Results show that APD increases with organic matter (OM) proxied by total organic carbon and decreases with clay minerals and the sum of quartz and feldspar. OM dramatically contributes to the APD as multiple-layer adsorption exists, and the APD for OM could be 1.4–8.5 times that for clay minerals. Other inorganic minerals contribute less to APD. The properties and constitution of the surface area instead of the volume fraction contribute to the APD in shales. APD does not show an obvious correlation with micropore volume, likely related to the ratio of micropore volume to the total pore volume. Here, we provide a significant and comprehensive study of petrological factors that impact the APD of subsurface shale gas, which will improve the estimation of supercritical adsorption and shale GIP under reservoir conditions. Also, the findings in this work can provide applications for subsurface carbon dioxide adsorption and storage.

58 GEOSCIENCES↗

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↗