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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Identifying Potential Geochemical and Microbial Impacts of Hydrogen Storage in a Deep Saline Aquifer

Hydrogen is valuable commodity and a promising energy carrier for variable energy production. Storage of hydrogen may occur through injection of hydrogen or a hydrogen/methane gas blend in subsurface reservoirs. However, the geochemical and biological reactions that may impact the stored hydrogen are not yet understood. Therefore, we collected samples from a deep storage aquifer located in the St. Peter Formation in southern Illinois. The reservoir material was primarily quartz with sulphur and iron deposits, while the major constituents of the fluid were chloride and sulphate. 16S rRNA gene amplicon sequencing revealed a low biomass microbial community that contained no obvious hydrogen-consuming bacteria. Next, we enriched a field sample to increase the biomass and completed a metagenomic analysis, finding a low number of genes present that are associated with hydrogen consumption. Then, we completed a series of reactor experiments under reservoir conditions with 15% H2/85% CH4 gas simulating a short-term hydrogen storage, high withdrawal scenario. We found minimal changes in the geochemistry or microbiology for the reactor experiments. This work suggests that short-term storage may be highly successful, although significant additional work needs to be completed in order to accurately evaluate the risks associated with long-term hydrogen storage scenarios. It is essential we continue to expand our understanding of the dynamics present in saline aquifers and provide new insights into how hydrogen storage may impact underground geological storage environments.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Nationwide cost and capacity estimates for sedimentary basin geothermal power and implications for geologic CO 2 storage

Sedimentary basins are naturally porous and permeable subsurface formations that underlie approximately half of the United States. In addition to being targets for geologic CO 2 storage, these resources could supply geothermal power: sedimentary basin geothermal heat can be extracted with water or CO 2 and used to generate electricity. The geothermal power potential of these basins and the accompanying implication for geologic CO 2 storage are, however, understudied. Here, we use the Sequestration of CO 2 Tool (SCO2T PRO ) and the generalizable GEOthermal techno-economic simulator (genGEO) to address this gap by a) estimating the cost and capacity of sedimentary basin geothermal power plants across the United States and b) comparing those results to nationwide CO 2 sequestration cost and storage potential estimates. We find that across the United States, using CO 2 as a geothermal heat extraction fluid reduces the cost of sedimentary basin power compared to using water, and some of the lowest cost capacity occurs in locations not typically considered for their geothermal resources (e.g., Louisiana, South Dakota). Additionally, using CO 2 effectively doubles the sedimentary basin geothermal resource base, equating to hundreds of gigawatts of new capacity, by enabling electricity generation in geologies that are otherwise (with water) too impermeable, too thin, too cold, or not deep enough. We find there is competition for the best sedimentary basin resources between water- and CO 2 -based power, but no overlap between the lowest-cost resources for CO 2 storage and CO 2 -based power. In this way, our results suggest that deploying CO 2 -based power may increase the cost of water based systems (by using the best resources) and the cost of CO 2 storage (by storing CO 2 in locations that otherwise may not be targeted). As such, our findings demonstrate that determining the best role for sedimentary basins within the energy transition may require balancing tradeoffs between competing priorities.

CPG↗

DOE OSTI2024

Investigation of acoustic waves under subsurface conditions to improve the predictions of rock mechanical properties and natural fracture characteristics

Mechanical properties and natural fracture characteristics are critical to investigate for subsurface engineering applications, including carbon storage, well drilling, and stimulation, as they govern rock stability, fluid flow, and mechanical behavior under stress. This dissertation integrates experimental and machine learning approaches to enhance the prediction and understanding of these properties by analyzing acoustic wave behavior under varied subsurface conditions. First, the influence of temperature, pore pressure, and supercritical CO2 (scCO2) saturation on poroelastic properties is examined using Gray Berea sandstone samples. The results show that temperature and pore pressure significantly affect the bulk modulus and Biot’s coefficient, while scCO2 saturation impacts rock compressibility, informing strategies for effective geological carbon storage. The study extends this understanding by experimentally evaluating the impact of reservoir depletion on the dynamic mechanical properties of the emerging Caney shale in South Oklahoma with the employment of unsupervised machine learning to predict static mechanical properties across the Caney shale. Integrating petrophysical data and chemostratigraphy, the workflow—featuring K-means clustering, principal component analysis (PCA), and inverse distance weighting (IDW)—improves stratigraphic characterization and the estimation of static-to-dynamic modulus ratios, which is vital for optimizing drilling and stimulation strategies. Finally, the work explores how natural fracture characteristics in shale influence acoustic waveforms and shear wave splitting (SWS) analysis. Experimental data on fractured samples under different stress and temperature conditions, combined with machine learning models such as K-nearest neighbors (KNN) and extreme gradient boosting (XGBoost), reveal key fracture properties impacting SWS and wave propagation. Together, these studies provide a comprehensive framework for linking acoustic wave behavior with rock properties, advancing the methods for monitoring and predicting geomechanical changes. The insights offered valuable implications for safer, more efficient CO2 injection, hydrocarbon extraction, and subsurface management.

Elkholy, Sherif↗

DOE OSTI2024

CarbonSAFE Phase II: Optimizing Alabama’s CO 2 Storage in Shelby County, Alabama (Project OASIS), Milestone M3 - Site Specific Drilling Report

The Phase II Storage Complex Feasibility project, entitled "Optimizing Alabama's CO 2 Storage in Shelby County, Alabama (Project OASIS)," is a part of the DOE/NETL's CarbonSAFE initiative. The Project is managed by the Southern States Energy Board (SSEB), and includes participation from Advanced Resources International, Inc. (ARI), Crescent Resource Innovation, Southern Company, Alabama A&M University, Auburn University, and Oklahoma State University. Project OASIS is working to establish the foundation for a commercial-scale geologic storage complex for CO 2 captured from Plant Gaston (home of the National Carbon Capture Center) and surrounding industrial sources of CO 2 located in Shelby County, Alabama. The Project objectives are: • Demonstrate that the subsurface saline formations at the storage complex can store commercial volumes of CO 2 safely and permanently. • Develop a comprehensive Community Benefits Plan. • Develop the infrastructure framework for a CO 2 storage hub. • Develop a rigorous risk registry and to conduct a comprehensive risk assessment. • Develop a monitoring plan. • Develop a comprehensive site characterization plan to support an Underground Injection Control Class VI Permit in a future Phase III program. • Evaluate the commercial viability of the project. Project OASIS is about 30 miles southeast of Birmingham, Alabama within a geologic province called the Valley and Ridge (Figure 1.1). The Valley and Ridge Province comprises a sequence of Paleozoic carbonate and clastic rocks that underwent structural deformation during the Alleghanian Orogeny. Storage prospects occur in relatively flat lying structural panels located between thrust faults. Available geologic studies related to hydrocarbon exploration suggest that Cambro-Ordovician carbonates and Cambrian clastic units offer multiple potential storage intervals, and that regional confining systems are present, such as the tectonically thickened Floyd-Parkwood Shale. The surface property is owned by a timber and land stewardship company, The Westervelt Company, Inc., who worked with the Project Team to select and prepare adequate sites for geologic assessment. The purpose of drilling the Westover Stratigraphic Test Well #2 was to collect geologic data to model the feasibility of commercial scale CO 2 injection and storage. This includes geological and geophysical evaluations, reservoir engineering analyses, and risk assessments.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2024

Managing Subsurface Pressure Buildup and Interference in Commercial-Scale CO 2 Storage Project with Proximal Injection Wells

Large-scale decarbonization using carbon capture and storage (CCS) is likely to involve many commercial-scale CO 2 storage projects located in close proximity to each other. This close proximity raises concerns over pressure interference among the storage projects. Pressure interference between injection and storage efforts can reduce the practicable CO 2 storage resource and force wells to inject CO 2 at a lower rate to avoid the fracture pressure thresholds per United States Environmental Protection Agency (EPA) Class VI well regulations to preserve injection and confining zone integrity and potentially mitigate against inducing seismic activity. These analyses employ numerical full-physics reservoir modeling to evaluate how pressure buildup fronts and CO 2 plumes evolve under commercial-scale injection volumes of CO 2 in which multiple storage sites located in close proximity occur in tandem. The simulation models mimic injection at pseudo basin-scale and assume homogeneous saline formation(s) as storage targets with a pair of upper and lower homogeneous seal layer/s. These analyses specifically investigate the efficacy of two basin-wide reservoir pressure management strategies in addressing the technical challenges associated with pressure buildup and CO 2 plume commingling. The strategies explored include: 1) enlarging the area of injection well spacing (WS) and 2) storing CO 2 in a stacked sequence (SSS) of saline formations compared to a single formation. The storage and confining zones properties assumed were held common across the scenarios, unless specified otherwise. Analyses results show that after injecting 4 million tons per year for 30 years using 4 separate wells (each injecting 1 million metric tons per year), the radius of CO 2 plume extends to a mere 3 km or less from injection wells. Meanwhile, the radius of pressure buildup ranges on the order of tens to a few hundreds of kilometers, depending on the magnitude of pressure buildup threshold that one would use to define the front. CO 2 plume commingling from different injection wells appears to occur 50 years post-injection, especially under scenarios with narrowly spaced (i.e., < 5 km apart) injection well locations. Findings from sensitivity cases on the well spacing suggest that storage formations modeled would require different well spacing to avoid fracture pressure thresholds. For instance, modeled storage formations with high fracture gradients (i.e., 0.8 psi/ft) would need less than 5–km well spacing, whereas those with lower fracture gradients (i.e., 0.7 psi/ft) would need approximately 20–km well spacing, based on assumed modeling parameters. Under stacked injection, the pressure challenges (described above) still exist but are more alleviated due to distributing the same injection volume across more available reservoir volume. These analyses demonstrate that stacked-sequence storage can effectively address the challenges, while still providing the same target CO 2 storage volumes and allowing a large number of storage projects to be deployed in the same basin by better utilizing the available storage resource across different reservoir depths. Among cases modeled, the resulting pressure buildup front is most suppressed when each storage project distributes injection volumes over several wells, each of which injects a portion of the total CO 2 across the stacked sequence. This strategy results in the smallest CO 2 aerial footprint amongst scenarios evaluated but also shows the largest reduction in the pressure buildup at the top of perforation at the injection wells (upwards of approximately 42 percent compared to the commercial-scale single-formation storage), the result of which is crucial to maintain caprock integrity. The findings presented by this research draw attention to the importance of greater coordination among storage operators and regulatory stakeholders to foster the upscaling and deployment of CCS. These analyses provide insights into required decision-making when considering multi-project deployment in a shared basin. Because these analyses evaluate a very specific geologic situation, they bear further investigations across other geologic situations.

42 ENGINEERING↗

DOE OSTI2024

Enhancing predictive understanding and accuracy in geological carbon dioxide storage monitoring: Simulation and history matching of tracer transport dynamics

Co-injection of conservative tracers with carbon dioxide (CO 2 ) is a viable tool for monitoring subsurface processes during geological CO 2 storage (GCS). This research investigates the simulation and history-matching of a gas tracer (sulfur hexafluoride, SF 6 ) during CO 2 flooding, employing a core flooding result in Berea sandstone. Four extensively used saturation functions are assessed for their efficacy in history matching of CO 2 /SF 6 injection at the core scale. The history-matching process incorporates particle swarm optimization (PSO) to fine-tune constitutive relationships parameters. Next, employing transport models at the aquifer scale, we interrogate the impact on tracer transport and mixing of saturation function uncertainties, arising from the non-uniqueness of constitutive relationships parameters and saturation function types. To assess the effects of geological heterogeneity on behavior of tracer breakthrough curves (BTCs), we employ two normalized parameters assessing the degree of mixing and SF 6 breakthrough time. The aquifer-scale investigation encompasses both homogeneous and heterogeneous systems with and without capillary heterogeneity effects. Our findings underscore the critical importance of addressing saturation function uncertainties, emphasizing the significance of auxiliary experiments and innovative methodologies to enhance predictive accuracy. The findings highlight significant disparities in arrival times, BTC peaks, tails, and mixing levels, even under optimal conditions. Heterogeneity, with or without capillary heterogeneity, plays a crucial role in shaping BTC variations, resulting in accelerated SF 6 breakthrough times and reduced BTC peaks. Evaluation of monitoring points distant from the injector reveals a dampening effect on the SF6 BTC peak, particularly in heterogeneous systems with capillary heterogeneity, where the peak is halved. These insights underscore the challenges associated with tracer monitoring and the necessity for enhanced methodologies to improve predictive accuracy in subsurface processes during GCS.

58 GEOSCIENCES↗

DOE OSTI2025

High-Temperature Aquifer Thermal Energy Storage (HT-ATES) Projects in Germany and the Netherlands—Review and Lessons Learned

Aquifer thermal energy storage (ATES) is a concept that can help to address heating and cooling needs through the use of the subsurface as a seasonal thermal energy storage (STES) system. Over 2800 ATES systems have been deployed with storage temperatures typically below 25 °C and only a few with higher temperatures (>40 °C), which would increase the energy density and utility of the stored thermal fluids. Until now, only a few high-temperature aquifer thermal energy storage (HT-ATES) projects have been initiated and are still in operation. These HT-ATES projects have encountered a range of technical and non-technical challenges. This study reviews ten such projects: four in Germany and six in the Netherlands. The non-technical issues include public acceptance, a lack of regulatory framework for these systems, managing overlapping uses of the subsurface, managing changes with the providers and off-takers of thermal energy, and obtaining financing to implement these projects. Common technical issues include geological factors such as incomplete characterization of the subsurface and reservoir heterogeneity; geochemical issues such as mineral scaling, corrosion, and biofouling; lower than expected thermal recovery; and issues with system design and reliability. This review highlights benefits and challenges faced by HT-ATES projects with the goal to use the lessons learned to improve the siting, design, development, and operation of such systems. Recommendations include improved initial subsurface site characterization, use of coupled process models to optimize system design and predict system performance, cascaded uses of stored thermal energy to better utilize the stored heat, monitoring networks to provide feedback on system performance, and expanded system scale to allow for continued operation even when maintenance of some system components is required. Techno-economic modeling and risk analysis could be used to optimize such HT-ATES project design and identify key factors that will affect sustained economic viability. In addition, design flexibility is important for these systems to allow for changing conditions regarding the supply and demand of thermal energy. Adopting these findings should improve the performance and reduce the risks for future HT-ATES projects worldwide.

15 - GEOTHERMAL ENERGY↗

DOE OSTI2021

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

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

02 PETROLEUM↗

DOE OSTI2024

Optimizing Carbon Capture, Transport, and Storage: Overcoming Challenges with Machine Learning and Cost-Benefit Analysis

Carbon Capture, Utilization, and Storage (CCUS) is a critical strategy for reducing CO₂ emissions and mitigating climate change. However, its widespread deployment faces numerous challenges across the capture, transport, and storage phases. These challenges include the technical complexity of predicting subsurface behaviors during CO₂ injection, ensuring long-term storage integrity, optimizing transportation networks, and balancing the economic and environmental trade-offs. Addressing these issues requires an integrated approach combining advanced subsurface modeling with system-level analyses to assess costs, risks, and benefits. This presentation provides an overview of studies conducted by the National Energy Technology Laboratory (NETL) to tackle these challenges. NETL’s efforts encompass cutting-edge research in subsurface fluid behavior machine learning predictions, alongside the development of innovative tools for system optimization and economic evaluation. By bridging technical expertise and strategic analysis, NETL aims to advance the deployment of CCUS technologies to support global decarbonization efforts. Presented at the Carnegie Mellon University CEE IESS Student Seminar October 4, 2024.

Shih, Chung Yan↗

DOE OSTI2025

Mechanistic understanding of carbon mineralization in fracture systems using microfluidics

Carbon mineralization in mafic and ultramafic rocks presents an opportunity for permanent carbon storage in the Earth's subsurface. However, due to their lower permeability, pre-existing fracture networks are key for mineralization to occur. Therefore, to fully develop this technology, a mechanistic understanding of the mineralization behavior in fractures with the consideration of hydrodynamic components is required. We use high-pressure microfluidics to investigate key mechanisms influencing dissolution–precipitation in a fracture network. The experiments were conducted in micromodels made of natural rocks with a comb-shaped flow channel to mimic a fracture network. This enabled studying the effect of injection rate on coupled dissolution–precipitation in advection and diffusion-dominated flow paths. We used gypsum carbonation as an analog reaction to allow for realistic experimental time frames due to its rapid reaction kinetics. The experimental work is coupled with high-fidelity numerical simulations to enhance our understanding of the parameters affecting the mineralization reaction. Our results demonstrate the importance of flow rate on the rate and nature of the gypsum carbonation reaction revealing that higher flow rates enable deeper penetration of the mineral precipitation front into the dead-end channels. This is an important finding since for sustained mineralization in a fracture network, precipitation in dead-ends while still allowing for flowing fractures is critical. Detailed characterization of the precipitates showed that lower flow rates led to porous and loose precipitates in the form of aragonite while higher flow rates mimicked supersaturation behavior leading to the formation of calcite. The reactive transport simulations further demonstrated the significance of flow velocity in advection-dominated channels to influence the efficiency of carbon mineralization in diffusion-dominated channels, potentially clogging of dead-end channels. These findings highlight the need for coupling chemical, mechanical, and hydrodynamic processes to evaluate the nature and extent of carbon mineralization in fractured media critical for permanent storage in mafic and ultramafic formations. This research further highlights the need for more investigation in potential subsurface fracture generation techniques to aid carbon mineralization.

25 ENERGY STORAGE↗

DOE OSTI2024

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↗

DOE OSTI2025

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]↗

DOE OSTI2023

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↗

DOE OSTI2024

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.↗

DOE OSTI2022

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.↗

DOE OSTI2024

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↗

DOE OSTI2023

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↗

DOE OSTI2022

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↗