Search NASA⌕ Search

SEARCH · Search NASA

Results for “subsurface permeability”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Estimating Watershed Subsurface Permeability From Stream Discharge Data Using Deep Neural Networks

Subsurface permeability is a key parameter in watershed models that controls the contribution from the subsurface flow to stream flows. Since the permeability is difficult and expensive to measure directly at the spatial extent and resolution required by fully distributed watershed models, estimation through inverse modeling has had a long history in subsurface hydrology. The wide availability of stream surface flow data, compared to groundwater monitoring data, provides a new data source to infer soil and geologic properties using integrated surface and subsurface hydrologic models. As most of the existing methods have shown difficulty in dealing with highly nonlinear inverse problems, we explore the use of deep neural networks for inversion owing to their successes in mapping complex, highly nonlinear relationships. We train various deep neural network (DNN) models with different architectures to predict subsurface permeability from stream discharge hydrograph at the watershed outlet. The training data are obtained from ensemble simulations of hydrographs corresponding to an permeability ensemble using a fully-distributed, integrated surface-subsurface hydrologic model. The trained model is then applied to estimate the permeability of the real watershed using its observed hydrograph at the outlet. Our study demonstrates that the permeabilities of the soil and geologic facies that make significant contributions to the outlet discharge can be more accurately estimated from the discharge data. Their estimations are also more robust with observation errors. Compared to the traditional ensemble smoother method, DNNs show stronger performance in capturing the nonlinear relationship between permeability and stream hydrograph to accurately estimate permeability. Our study sheds new light on the value of the emerging deep learning methods in assisting integrated watershed modeling by improving parameter estimation, which will eventually reduce the uncertainty in predictive watershed models.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1998

Magmatic Intrusions and a Hydrothermal Origin for Fluvial Valleys on Mars

Numerical models of Martian hydrothermal systems demonstrate that systems associated with magmatic intrusions greater than several hundred cubic kilometers can provide sufficient groundwater outflow to form the observed fluvial valleys, if subsurface permeability exceeds about 1.0 darcy. Groundwater outflow increases with increasing intrusion volume and subsurface permeability and is relatively insensitive to intrusion depth and subsurface porosity within the range considered here. Hydrothermally-derived fluids can melt through 1 to 2 km thick ice-rich permafrost layers in several thousand years. Hydrothermal systems thus provide a viable alternative to rainfall for providing surface water for valley formation. This mechanism can form fluvial valleys not only during the postulated early warm, wet climatic epoch, but also during more recent epochs when atmospheric conditions did not favor atmospheric cycling of water. The clustered distribution of the valley networks on a given geologic surface or terrain unit of Mars may also be more compatible with localized, hydrothermally-driven groundwater outflow than regional rainfall. Hydrothermal centers on Mars may have provided appropriate environments for the initiation of life or final oases for the long-term persistence of life.

Gulick, Virginia C↗
DOE OSTI2025

The Role of Snowmelt and Subsurface Heterogeneity in Headwater Hydrology of a Mountainous Catchment in Colorado: A Model‐Data Integration Approach

Mountainous headwater streams are sustained by both snowmelt‐driven streamflow and groundwater discharge in the Upper Colorado River Basin. However, predicting headwater stream discharge magnitude and peak flow timing is challenging in mountainous terrains, where snowmelt rates vary with vegetation type and elevation, and heterogeneous subsurface physical properties influence groundwater storage and its release. We used a model‐data integration approach to investigate the roles of snowmelt and subsurface structure in stream discharge and groundwater level. We ran an ensemble of 100 integrated surface‐subsurface hydrologic models for a mountainous headwater catchment near Crested Butte, Colorado, USA. We also evaluated and calibrated these models against observed data sets, including snow depth measurements using distributed temperature probes, stream discharge, and groundwater levels. Calibration with multiple data sources using neural density estimators has further constrained uncertainty in subsurface properties and snowmelt rates. Results indicated that observed slower snowmelt rates in evergreen forests delayed the peak flow and baseflow onset. In upstream areas with lower subsurface permeability, water was stored within the subsurface but was not released as interflow or shallow groundwater flow, and thereby not contributing to downstream streamflow during recession limb periods. Double peaks in groundwater occurred in areas with spatial subsurface heterogeneity, in our case due to the contrast between granodiorite and Mancos shale. These process‐based insights into groundwater and snowmelt dynamics in mountainous headwaters will help improve predictions of headwater hydrology.

Wang, Lijing [University of Connecticut, Storrs, C↗
DOE OSTI2022

Deep learning to estimate permeability using geophysical data

Time-lapse electrical resistivity tomography (ERT) is a popular geophysical method to estimate three-dimensional (3D) permeability fields from electrical potential difference measurements. Traditional inversion and data assimilation methods are used to ingest this ERT data into hydrogeophysical models to estimate permeability. Due to ill-posedness and the curse of dimensionality, existing inversion strategies provide poor estimates and low resolution of the 3D permeability field. Recent advances in deep learning provide us with powerful algorithms to overcome this challenge. This paper presents a deep learning (DL) framework to estimate the 3D subsurface permeability from time-lapse ERT data. To test the feasibility of the proposed framework, we train DL-enabled inverse models on simulation data. Each measurement in both synthetic and field data is standardized by removing the mean and scaling the time-series to unit variance. This pre-processing step is necessary to bring simulation data closer to field observations. Subsurface process models based on hydrogeophysics are used to generate this synthetic data. Training performed on limited simulation data resulted in the DL model over-fitting. An advanced data augmentation based on mixup is implemented to generate additional training samples to overcome this issue. This mixup technique creates weakly labeled (low-fidelity) samples from strongly labeled (high-fidelity) data. The weakly labeled training data is then used to develop DL-enabled inverse models and reduce over-fitting. As both time-lapse ERT (1133048 features/realization) and 3D permeability (585453 features/realization) data samples are from a high-dimensional space, principal component analysis (PCA) is employed to reduce dimensionality. Encoded ERT and encoded permeability are generated using the trained PCA estimators. A deep neural network is then trained to map the encoded ERT to encoded permeability. This mixup training and unsupervised learning allowed us to build a fast and reasonably accurate DL-based inverse model under limited simulation data. Results show that proposed weak supervised learning can capture salient spatial features in the 3D permeability field. Quantitatively, the average mean squared error (in terms of the natural log) on the strongly labeled training, validation, and test datasets is less than 0.5. The R 2 -score (global metric) is greater than 0.75, and the percent error in each cell (local metric) is less than 10%. Finally, an added benefit in terms of computational cost is that the proposed DL-based inverse model is at least O(10 4 ) times faster than running a forward model once it is trained. Data generation, DL model training, and hyperparameter tuning to identify optimal neural network architectures utilized high-performance computing resources while the DL inference is performed on a standard laptop. Approximately, O(10 5 ) processor hours are used for generating data and DL tuning and training. We acknowledge that the data generation and DL model development are expensive. But once a DL model is trained, it can be re-used for inversion rapidly for the given system, with set physics and domain. Note that traditional inversion may require multiple forward model simulations (e.g., in the order of 10 to 1000), which are very expensive. This computational savings ≈ O(10 5 ) – O(10 7 )) makes the proposed DL-based inverse model attractive for subsurface imaging and real-time ERT monitoring applications due to fast and yet reasonably accurate estimations of permeability field.

58 GEOSCIENCES↗
DOE OSTI2026

Coupled Investigation of Fracture Permeability Impact on Reservoir Stress and Seismic Slip Behavior (Final Technical Report)

Enhanced Geothermal Systems (EGS) produce clean energy by circulating fluid through hot rock deep underground and bringing that heat to the surface to generate electricity. For this process to work reliably, fluids must be able to move efficiently through networks of natural or engineered fractures in the rock. Enhancing and maintaining subsurface permeability over time is essential for sustainable energy production. However, fluid injection changes the underground temperature, pressure, rock stress, and chemistry, which can alter permeability and sometimes trigger earthquakes. Predicting these interconnected processes remains a key challenge. To address this, we combined high-temperature laboratory experiments with high-fidelity simulations to better understand how fractures in geothermal reservoirs evolve over time. Our experiments measured how fractures respond to stress, slip, slip rate, and chemical reactions under geothermal conditions. These data were integrated into coupled thermal-hydrological-mechanical-chemical and earthquake (THMC+E) models tailored to the Utah FORGE site. The validated modeling framework improves predictions of reservoir performance and seismic response and helps guide operational decisions. This work reduces technical risk and strengthens the scientific foundation needed to make geothermal energy a reliable and scalable clean energy resource.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

An Integrated Feasibility Study of Reservoir Thermal Energy Storage (RTES) in Portland, OR, USA

In regions with long cold overcast winters and sunny summers, Deep Direct-Use (DDU) can be coupled with Reservoir Thermal Energy Storage (RTES) technology to take advantage of pre-existing subsurface permeability to save summer heat for later use during cold seasons. Many aquifers worldwide are underlain by permeable regions (reservoirs) containing brackish or saline groundwater that has limited beneficial use due to poor water quality. We investigate the utility of these relatively deep, slow flowing reservoirs for RTES by conducting an integrated feasibility study in the Portland Basin, Oregon, USA, developing methods and obtaining results that can be widely applied to RTES systems elsewhere. As a case study, we have conducted an economic and social cost-benefit analysis for the Oregon Health and Science University (OHSU), a teaching hospital that is recognized as critical infrastructure in the Portland Metropolitan Area. Our investigation covers key factors that influence feasibility including 1) the geologic framework, 2) heat and fluid flow modeling, 3) capital and maintenance costs, 4) the regulatory framework, and 5) operational risks. By pairing a model of building seasonal heat demand with an integrated model of RTES resource supply, we determine that the most important factors that influence RTES efficacy in the study area are operational schedule, well spacing, the amount of summer heat stored (in our model, a function of solar array size), and longevity of the system. Generally, heat recovery efficiency increases as the reservoir and surrounding rocks warm, making RTES more economical with time. Selecting a base-case scenario, we estimate a levelized cost of heat (LCOH) to compare with other sources of heating available to OHSU and find that it is comparable to unsubsidized solar and nuclear, but more expensive than natural gas. Additional benefits of RTES include energy resiliency in the event that conventional energy supplies are disrupted (e.g., natural disaster) and a reduction in fossil fuel consumption resulting in a smaller carbon footprint. Key risks include reservoir heterogeneity and a possible reduction in permeability through time due to scaling (mineral precipitation). Lastly, a map of thermal energy storage capacity for the Portland Basin yields a total of 43,400 GWh, suggesting tremendous potential for RTES in the Portland Metropolitan Area.

14 SOLAR ENERGY↗
DOE OSTI2024

A Workflow for Characterizing Legacy Wells as Potential Leakage Pathways for Integration to NRAP-Open-IAM

Carbon capture and storage is a crucial component of climate change mitigation strategies, involving the capture of carbon dioxide (CO2) from point sources and its injection into permeable subsurface formation. Many suitable CO2 storage sites coincide with legacy wells since the conditions that kept hydrocarbons in-situ for thousands of years are also ideal for storage of carbon dioxide. To protect underground sources of drinking water (USDW) during greenhouse gas injection, the Environmental Protection Agency (EPA) mandates area of review evaluations. These evaluations ensure that drinking water sources would not be contaminated by injected fluids. They include identification of legacy wellbores, integrity assessments, and implementing any necessary corrective action. Previous assessment approaches of legacy wells include high-level scoring of regional data and well construction and abandonment evaluation. This work describes a novel methodology that evaluates well construction and abandonment, ranks them based on complexity, and performs a risk assessment with NRAP-Open-IAM. A workflow of the methodology is presented, highlighting its capabilities and limitations.

Wise, Jarrett↗
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↗
NASA NTRS1990

Subsurface energy storage and transport for solar-powered geysers on Triton

The location of active geyser-like eruptions and related features close to the current subsolar latitude on Triton suggests a solar energy source for these phenomena. Solid-state greenhouse calculations have shown that sunlight can generate substantially elevated subsurface temperatures. A variety of models for the storage of solar energy in a subgreenhouse layer and for the supply of gas and energy to a geyser are examined. 'Leaky greenhouse' models with only vertical gas transport are inconsistent with the observed upper limit on geyser radius of about 1.5 km. However, lateral transport of energy by gas flow in a porous N2 layer with a block size on the order of a meter can supply the required amount of gas to a source region about 1 km in radius. The decline of gas output to steady state may occur over a period comparable with the inferred active geyser lifetime of 5 earth years. The required subsurface permeability may be maintained by thermal fracturing of the residual N2 polar cap. A lower limit on geyser source radius of about 50 to 100 m predicted by a theory of negatively buoyant jets is not readily attained.

Kirk, Randolph L.↗
NASA NTRS1996

Hydrogeological Interpretation of Candidate Origin Sites for Martian Meteorite ALH84001

Barlow (this meeting) has identified two potential source craters for the martian meteorite ALH84001. The craters are at 11.7 deg S, 243.3 deg W (Mare Tyrrhenum site) and 14.0 deg S, 343.5 deg W (Sinus Sabaeus site). As noted by Barlow, both craters lie in the heavily cratered terrain (HCT) and are adjacent to fluvial valleys, Here I explore the fluvial history of these areas based upon the surrounding valley morphology. The most prominent valley network at the Sabaeus site is Evros Vallis. This wide, flat-floored valley is approximately 600 km long with an average width of 2.5 km and a depth of 220 m. The eroded volume of the entire Evros network is approximately 6 x 10(exp 11) cc. This is typical for networks located in the heavily cratered terrain (e.g. Warrego and Parana Valles). Evros is also an isolated valley system. No similar networks are found in the surrounding terrain. Thus it is unlikely that Evros formed as a result of widespread rainfall. A localized water source, such as discharge of a hydrothermal system or localized melting of snowfall, seems more consistent Previous modeling has demonstrated that only hydrothermal systems associated with high permeability subsurfaces can discharge sufficient water to form a valley network. The bulk of the discharge from such systems is consequently low temperature, slightly heated water Precipitation of calcium carbonate by low temperature fluids is consistent with most interpretations of the geochemistry of ALH84001. Available imagery at the Tyrrhenum site is of lesser quality. While eroded units of the HCT are nearby, there are no comparable well developed valley networks at this site. Erosion is instead manifested predominantly as gullies on slopes. This style of erosion suggests that water was not present at this site for the length of time as at the more integrated Sabaeus site. The superposition of fluidized ejecta blankets suggests however that ground water or ground ice was still present at this locality fairly recently, if the chosen crater is indeed correct.

Gulick, Virginia C.↗
DOE OSTI2024

Computationally efficient and error aware surrogate construction for numerical solutions of subsurface flow through porous media

Limiting the injection rate to restrict the pressure below a threshold at a critical location can be an important goal of simulations that model the subsurface pressure between injection and extraction wells. The pressure is approximated by the solution of Darcy’s partial differential equation for a given permeability field. The subsurface permeability is modeled as a random field since it is known only up to statistical properties. This induces uncertainty in the computed pressure. Solving the partial differential equation for an ensemble of random permeability simulations enables estimating a probability distribution for the pressure at the critical location. These simulations are computationally expensive, and practitioners often need rapid online guidance for real-time pressure management. An ensemble of numerical partial differential equation solutions is used to construct a Gaussian process regression model that can quickly predict the pressure at the critical location as a function of the extraction rate and permeability realization. The Gaussian process surrogate analyzes the ensemble of numerical pressure solutions at the critical location as noisy observations of the true pressure solution, enabling robust inference using the conditional Gaussian process distribution. Our first novel contribution is to identify a sampling methodology for the random environment and matching kernel technology for which fitting the Gaussian process regression model scales as O ( n log n ) instead of the typical O ( n 3 ) rate in the number of samples n used to fit the surrogate. The surrogate model allows almost instantaneous predictions for the pressure at the critical location as a function of the extraction rate and permeability realization. Our second contribution is a novel algorithm to calibrate the uncertainty in the surrogate model to the discrepancy between the true pressure solution of Darcy’s equation and the numerical solution. Finally, although our method is derived for building a surrogate for the solution of Darcy’s equation with a random permeability field, the framework broadly applies to solutions of other partial differential equations with random coefficients.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Hawai‘i Play Fairway (Final Report)

Most of Hawai'i's geothermal resources are blind—their manifestations, such as hot springs and steam vents, do not appear on the ground surface because the heated water flows far below. With the exception of $K\bar{i}lauea East Rift Zone$, in most areas of Hawai'i, high lateral permeability in the first kilometer below ground surface prevents surface thermal features from developing. As a methodology for discovering these blind resources, Play Fairway Analysis (PFA) involves finding potential locations of blind hydrothermal systems and describing potential geothermal sources in rift-zone settings. Using the PFA to find Hawai'i's geothermal resources, the University of Hawai'i (UH) conducted the Hawai‘i Play Fairway Project, Hawai'i's first statewide geothermal resource assessment since 1985. Sponsored by the U.S. Department of Energy, the Hawai'i Play Fairway Project provided an updated resource assessment, a roadmap for additional exploration activities, and the identification of areas for further exploration. Benefitting from UH's core competency in earth sciences and experienced geothermal researchers, the project comprised three phases. During the first phase, the team identified, compiled, and ranked existing geologic, groundwater, and geophysical datasets relevant to subsurface heat, fluid and permeability. Using a Bayesian statistical approach, the team developed a statistical methodology to integrate these data into a resource probability map. The team evaluated the confidence in the probability value and considered development viability of areas with geothermal resources. With these analyses, the team identified 10 locations in the Hawaiian Islands for exploration activities. For the second phase, the team collected new groundwater data in 10 locations across the state and new geophysical data on $L\bar{a}na‘i, Maui$, and central Hawai'i Island and modeled topographically induced stress to better characterize subsurface permeability. Analyzing the subsurface stresses, the team evaluated the potential for fracture-induced permeability. The team inverted the MT and gravity data to produce 3D models of resistivity and density, respectively, on $L\bar{a}na‘i$, across $Haleakal\bar{a}'s$ SW rift (Maui), and surrounding Mauna Kea (Hawai‘i Island). The team developed and applied a new method for incorporating depth information about resistivity, density, and potential for fracture-induced permeability into the statistical method for computing resource probability in these three focus areas. The team incorporated the new groundwater results with the new geophysical results and the calculations of potential for fracture-induced permeability to produce updated maps of resource probability and confidence. Through combining data from the first and second phases, the team determined locations for further exploration during the third phase. For MT and gravity surveys, the team recommended $Kaua'i's$ $L\bar{i}hu'e$ $Basin$, the east rift of $Maui's$ $Haleakal\bar{a}$ volcano, and the southwest rift of Hawai'i Island's Mauna Loa volcano. The MT and gravity surveys aimed to enable improved confidence in the resource potential in these locations. For drilling deep groundwater well(s), the team recommended Southeast Mauna Kea and $L\bar{a}na's$ $P\bar{a}l\bar{a}wai$ $Basin$. During the third phase, further exploration involved drilling a groundwater well in $L\bar{a}na's$ $P\bar{a}l\bar{a}wai$ $Basin$ and performing more geophysical surveys. We deepened an existing water well proximal to our target area on $L\bar{a}na'i$ due to funding constraints that precluded us from spudding a new well that would exceed 1km depth. Drilling was preceded by a number of substantial elements including: writing an Environmental Assessment and the subsequent legal process, performance of deviation logging, lowering a camera down the well, coordinating site preparation with $P\bar{u}lama$ $L\bar{a}na'i$, shipping the UH-owned rig interisland, procuring supplies, and leading 3 community meetings on $L\bar{a}na'i$. Drilling occurred 24/7 the entire month of June 2019 over which time $L\bar{a}na'i$ $Well$ 10 was deepened from 427 m to 1057 m, with continuous core collected. We measured a roughly linear temperature gradient averaging 42°C/km and a maximum bottom hole temperature of 66°C. This gradient is more than twice the background for Hawai'i and within a range of gradients measured in this depth range for some exploration wells within KERZ. We consider these results encouraging for $L\bar{a}na'i's$ resource potential and recommend following with a slim hole within $L\bar{a}na'i's$ caldera (our target zone) to ~ 2 km. Further, the positive implications such results have for the island of O‘ahu are substantial - the shield stage of O'ahu's volcanoes ended 1-2 My earlier. However, O'ahu uses more electricity than the rest of the islands combined, and the utility recently called for 500-700MW of firm, dispatchable renewable electricity on O'ahu by 2033. In Phase 3, we also collected limited new encouraging groundwater data, and updated our thoughts on the probabilities of fluid and permeability at resource depths (Pr F = 1; Pr P = mostly unconstrained). Ultimately, we advocate for using our final probability of heat, and confidence in this probability, to drive the next phase of exploration. We contend further development of geothermal in Hawai‘i will enable the state to achieve its 100% renewable policy objective and Hawai'i to transition off of fossil fuels through geothermal discovery and development. The project not only produced a large amount of data and expanded the existing knowledge of Hawai'i's geothermal resources, but also produced publications, theses, presentations, core photos, datasets, media reports, television interviews, community events, and a blog. Students and new professionals benefitted from the project's hands-on research experiences and educational opportunities and earned awards and recognition.

15 GEOTHERMAL ENERGY↗
NASA NTRS2022

On the Stability and Phase Behavior of Titan's Subsurface Liquid Columns

On Titan, liquid hydrocarbon may stay in the subsurface porous permeable crust known as the alkanofer, analogous to water in Earth’s aquifer. In addition to pressure gradient, the subsurface liquid in alkanofers is subject to vertical compositional grading due to the gravity and temperature gradient. The common wisdom is that the liquid would normally stay underground in a stability established by the pressure that increases with depth as observed in aquifers on Earth. However, Titan’s liquids consist of nitrogen and hydrocarbons, mainly methane and ethane, the behavior of which is very sensitive to temperature and pressure. Consequently, the liquid density does not always increase with depth, thus may introduce a reverse density profile that leads to vertical convective instability of the liquid column. If reverse density profiles are present, capillary pressures arising from liquid trapped within small pores in the crust can help with the column stability. The liquid held in the capillaries can seal the space below it thus helping with the stability, unless the overpressure built from underneath becomes larger than the capillary pressure, which causes leakage to allow the liquids to seep upward from the deep. This situation is analogous to hydrocarbon seeps on Earth, where oil and natural gas escape the reservoir and flow slowly through network of cracks to the surface. An algorithm based on an extended Gibbs equation commonly used in petroleum reservoir engineering is employed in this work to produce pressure, density, and compositional profiles for evaluating the stability and phase behavior of Titan’s subsurface fluids.

Sugata Tan↗
DOE OSTI2020

Experimental Investigation of Hydraulic Fracturing and Stress Sensitivity of Fracture Permeability Under Changing Polyaxial Stress Conditions

Understanding and predicting fracture propagation and subsequent fluid flow characteristics is critical to geoenergy technologies that engineer and/or utilize favorable geological conditions to store or extract fluids from the subsurface. Fracture permeability decreases nonlinearly with increasing normal stress, but the relationship between shear displacement and fracture permeability is less well understood. We utilize the new Geo–Reservoir Experimental Analogue Technology (GREAT cell), which can apply polyaxial stress states and realistic reservoir temperatures and pressures to cylindrical samples and has the unique capability to alter both the magnitude and orientation of the radial stress field by increments of 11.25° during an experiment. We load synthetic analogue materials and real rock samples to stress conditions representative of 500–1,000 m depth, investigate the hydraulic stimulation process, and then conduct flow experiments while changing the fluid pressure and the orientation of the intermediate and minimum principal stresses. High–resolution circumferential strain measurements combined with fluid pressure data indicate fracture propagation can be both stable (no fluid pressure drop) and unstable (fluid pressure drop). The induced fractures exhibit both opening and shear displacements during their creation and/or during fluid flow with changing radial stress states. Flow tests during radial stress field rotation reveal that fracture normal effective stress has first–order control on fracture permeability but increasing fracture offset can lead to elevated permeabilities at maximum shear stress. Here, the results have implications for our conceptual understanding of fracture propagation as well as fluid flow and deformation around fractures.

58 GEOSCIENCES↗
DOE OSTI2021

Interfacial and Confinement-Mediated Organization of Gas Hydrates, Water, Organic Fluids, and Nanoparticles for the Utilization of Subsurface Energy and Geological Resources

Harnessing the subsurface geologic environments in an efficient and environmentally sustainable manner is challenged by uncertainties associated with predicting the fate of fluids and sustaining porosity and permeability in subsurface geologic environments. Some of these uncertainties arise from confined and interfacially induced structures of fluids in subsurface geologic environments. The formation of gas hydrates, phase transitions of confined fluids, assembly and deposition of heavy hydrocarbons, and the agglomeration and fate of nanoparticles in confined environments are summarized in this review. Nanoscale confinement contributes to anisotropic structures and dynamics of fluids, which is the basis for anomalous phase transition thermodynamics, reactivity, transport, and geomechanical behavior. In this review, we discuss the structures of confined fluids and deviation in observed properties from bulk fluids. The factors influencing the structures of confined fluids can be generally divided into two groups: (a) pore characteristics including pore size, pore surface chemistry, and pore geometry and (b) confined fluid/solid characteristics such as molecular structure, concentrations, charges, pore filling, and presence of additives. Scientific advancements and knowledge gaps in our understanding of the structures of confined fluids and the associated differences in observed properties compared to bulk fluids are discussed. Here, the phenomena discussed in this review are of particular relevance to our efforts in harnessing the subsurface environments for a low carbon future by increasing the utilization of geothermal energy, using CO 2 as a working fluid, and storing CO 2 in subsurface geologic environments.

42 ENGINEERING↗
DOE OSTI2025

Hidden Features: How Subsurface and Landscape Heterogeneity Govern Hydrologic Connectivity and Stream Chemistry in a Montane Watershed

ABSTRACT Hydrologic connectivity is defined as the connection among stores of water within a watershed and controls the flux of water and solutes from the subsurface to the stream. Hydrologic connectivity is difficult to quantify because it is goverened by heterogeniety in subsurface storage and permeability and responds to seasonal changes in precipitation inputs and subsurface moisture conditions. How interannual climate variability impacts hydrologic connectivity, and thus stream flow generation and chemistry, remains unclear. Using a rare, four‐year synoptic stream chemistry dataset, we evaluated shifts in stream chemistry and stream flow source of Coal Creek, a montane, headwater tributary of the Upper Colorado River. We leveraged compositional principal component analysis and end‐member mixing to evaluate how seasonal and interannual variation in subsurface moisture conditions impacts stream chemistry. Overall, three main findings emerged from this work. First, three geochemically distinct end members were identified that constrained stream flow chemistry: reach inflows, and quick and slow flow groundwater contributions. Reach inflows were impacted by historic base and precious metal mine inputs. Bedrock fractures facilitated much of the transport of quick flow groundwater and higher‐storage subsurface features (e.g., alluvial fans) facilitated the transport of slow flow groundwater. Second, the contributions of different end members to the stream changed over the summer. In early summer, stream flow was composed of all three end members, while in late summer, it was composed predominantly of reach inflows and slow flow groundwater. Finally, we observed minimal differences in proportional composition in stream chemistry across all four years, indicating seasonal variability in subsurface moisture and spatial heterogeneity in landscape and geologic features had a greater influence than interannual climate fluctuation on hydrologic connectivity and stream water chemistry. These findings indicate that mechanisms controlling solute transport (e.g., hydrologic connectivity and flow path activation) may be resilient (i.e., able to rebound after perturbations) to predicted increases in climate variability. By establishing a framework for assessing compositional stream chemistry across variable hydrologic and subsurface moisture conditions, our study offers a method to evaluate watershed biogeochemical resilience to variations in hydrometeorological conditions.

Johnson, Keira [College of Earth, Ocean, and Atmos↗
DOE OSTI2023

Physics-embedded inverse analysis with algorithmic differentiation for the earth’s subsurface

Abstract Inverse analysis has been utilized to understand unknown underground geological properties by matching the observational data with simulators. To overcome the underconstrained nature of inverse problems and achieve good performance, an approach is presented with embedded physics and a technique known as algorithmic differentiation. We use a physics-embedded generative model, which takes statistically simple parameters as input and outputs subsurface properties (e.g., permeability or P-wave velocity), that embeds physical knowledge of the subsurface properties into inverse analysis and improves its performance. We tested the application of this approach on four geologic problems: two heterogeneous hydraulic conductivity fields, a hydraulic fracture network, and a seismic inversion for P-wave velocity. This physics-embedded inverse analysis approach consistently characterizes these geological problems accurately. Furthermore, the excellent performance in matching the observational data demonstrates the reliability of the proposed method. Moreover, the application of algorithmic differentiation makes this an easy and fast approach to inverse analysis when dealing with complicated geological structures.

54 ENVIRONMENTAL SCIENCES↗