Search NASA⌕ Search

SEARCH · Search NASA

Results for “subsurface flow”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Interactions between molecular-scale processes and hyporheic exchange for understanding Fe-S-C cycling in riparian wetlands (Final Report)

Wetlands represent some of the most productive ecosystems on the planet and critically influence global environmental health. Specifically, wetlands promote water quality by transforming nutrients and organic compounds and sequestering metals and contaminants. Riparian wetland hyporheic zones, where toxic surface water and anoxic groundwater mix, exhibit dynamic conditions that drive steep redox gradients and promote hotspots of diverse and fluctuating microbial activity. Changes in climate, water quality, and water quantity can disturb hydrologic flow and biogeochemical processing in these environments. Understanding how sulfate loading, impacted by hydrologic flux and anthropogenic inputs, influences iron and carbon cycling in wetlands will be crucial for predicting water quality issues driven by iron mineral precipitation and sorption, such as the release of heavy metals and other toxic elements. We used a fully integrated multi-scale and multi-method approach to develop a mechanistic understanding of how hydrologic flow influences coupled iron and sulfur cycles in riparian wetlands. This entailed hydrological, geochemical, and microbial observations at two locations: an anthropogenic sulfate-impacted riparian wetland in northern Minnesota, and a low-sulfate Fe-rich riparian wetland in Tims Branch at the Savannah River Site (SRS). These sites were characterized by hydrologically dynamic conditions where the stream and wetland systems oscillated between gaining (upward flow) and losing (downward flow) conditions that recharged the system occasionally with oxidants that fueled a variety of biogeochemical reactions. Aqueous geochemical measurements of surface water, groundwater, and porewater samples were made alongside solid-phase geochemical analyses of sediment gravity cores. Bulk X-ray absorption spectroscopy at the Advanced Photon Source (APS; Argonne) interrogated the speciation and distribution of Fe and S mineral phases of the sediments. Interestingly, it was discovered that, despite strongly reducing conditions, Fe(III) compounds and a variety of intermediate valence S compounds were stable in the subsurface. Indeed, compounds like thiosulfate, S(0), and intermediate valence organosulfur compounds were more prevalent than FeS and pyrite. The composition of the sediments did change with changes in hydrologic flow, showing their reactivity in changing redox conditions. The abundance of these intermediate S compounds was likely formed as a result of anaerobic oxidation by aqueous and solid-phase Fe(III) compounds, fueling a cryptic S cycle that is driving the breakdown of organic matter in the hyporheic zone. Microbiome surveys showed a core community that seemed stable across the landscape, but changed with increasing depth into the sediment. The community composition did not seem to change dramatically with changes in season or hydrologic flow, except for some organisms that have the potential to contribute to S cycling. More work is needed to confirm their functional activity. These fine process-scale analyses were placed within a dynamic field context using physical flow parameters from surface water and groundwater level measurements. These data helped shed light on sulfur-driven biogeochemical processes in hydrologically dynamic riparian wetlands, addressing a gap in our understanding about the impacts of pollution and other anthropogenic changes on ecologically sensitive environments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

High Frequency Monitoring and Nitrate Sourcing Reveals Baseflow and Stormflow Controls on Total Dissolved Nitrogen and Carbon Export Along a Rural‐Urban Gradient

Abstract Efforts to reduce nitrogen and carbon loading from developed watersheds typically target specific flows or sources, but across gradients in development intensity there is no consensus on the contribution of different flows to total loading or sources of nitrogen export. This information is vital to optimize management strategies leveraging source reductions, stormwater controls, and restorations. We investigate how solute loading and sources vary across flows and land‐use using high frequency monitoring and stable nitrate isotope analysis from five catchments with different sanitary infrastructure, along a gradient in development intensity. High frequency monitoring allowed estimation of annual loading and attribution to storm versus baseflows. Nitrate loads were 16 kg/km 2 /yr. from the forested catchment and ranged from 68 to 119 kg/km 2 /yr., across developed catchments, highest for the septic served site. Across developed catchments, baseflow contributions ranged from 40% of N loading to 75% from the septic served catchment, and the contribution from high stormflows increased with development intensity. Stormflows mobilized and mixed many surface and subsurface nitrate sources while baseflow nitrate was dominated by fewer sources which varied by catchment (soil, wastewater, or fertilizer). To help inform future sampling designs, we demonstrate that grab sampling and targeted storm sampling would likely fail to accurately predict annual loadings within the study period. The dominant baseflow loads and subsurface stormflows are not treated by surface water management practices primarily targeted to surface stormflows. Using a balance of green and gray infrastructure and stream/riparian restoration may target specific flow paths and improve management.

Delesantro, Joseph M.↗
DOE OSTI2023

Use of Stable Mercury Isotopes to Assess Mercury and Methylmercury Transformation and Transport across Critical Interfaces from the Molecular to the Watershed Scale (Final Report)

This project titled “Use of Stable Mercury Isotopes to Assess Mercury and Methylmercury Transformation and Transport across Critical Interfaces from the Molecular to the Watershed Scale” represents a collaborative effort between the University of Michigan (Jason Demers, PI) and Oak Ridge National Laboratory (Scott Brooks, co-I). Much has been learned about mercury (Hg) cycling in stream ecosystems, and East Fork Poplar Creek (EFPC) in particular, through decades of previous research. Nevertheless, some of the most fundamental questions regarding the sources of bioavailable Hg and its transformation to toxic methylmercury (MeHg) have remained unanswered. These fundamental questions include: (1) what are the sources and biogeochemical processes that lead to the input of dissolved Hg to stream water across critical subsurface interfaces within stream ecosystems, and EFPC in particular? and (2) what are the sources and biogeochemical processes that control the production and fate of bioaccumulative MeHg within stream ecosystems, and in EFPC in particular? To address these fundamental questions, our project aimed to couple laboratory experiments and field observations, both utilizing natural abundance Hg stable isotope techniques, to identify the processes responsible for generating mobile, bioavailable dissolved Hg from recalcitrant legacy sources within critical subsurface zones (e.g., streambed hyporheic zone, riparian floodplain subsurface). We used the isotopic signature of this bioavailable dissolved Hg to track its mobilization across these critical interfaces in order to link diffuse subsurface sources of dissolved Hg with increases in surface water dissolved Hg flux measured at the watershed scale. Additionally, our research aimed to determine the isotopic composition of MeHg within these same critical subsurface zones. We directly assessed the isotopic composition of MeHg within biota in order to gain insight into which subsurface sources of inorganic Hg and toxic MeHg are available for bioaccumulation within the EFPC ecosystem. Net fluxes of dissolved Hg along the flow path of EFPC were shown to vary spatially and temporally. In the Upper EFPC, within the Y12 boundary, stream water flux of dissolved Hg consistently decreased between the outfall and the downstream boundary of Y12 (57% ± 29%, 1SD). Within the Upper EFPC, an assessment of Hg isotopic composition suggested that losses were strongly reaction-driven, although isotopic diagnostics did not conform to any known processes. Downstream of Y12, in the upper reach of the Lower EFPC, dissolved Hg fluxes tended to increase during the dormant season (net gain of 11-120%), and decrease during the growing season (net loss of 23% +/- 18%, 1SD). In the downstream-most reach of Lower EFPC, dissolved Hg fluxes increased by 12-108% in 9 out of 10 monthly assessments. Overall, diffuse fluxes from the non-Y12 watershed accounted for 34% (+/- 17%, 1SD) of all dissolved Hg exported during base flow. Within Lower EFPC, an assessment of Hg isotopic composition was consistent with the contribution of diffuse Hg inputs from high-concentration hotspots within riparian floodplains and streambed hyporheic pore water. To investigate remobilization of recalcitrant Hg from legacy sediment sources, we developed procedures that coupled isotopic analysis with sequential extractions of streambed sediment. We found that the proportion of weakly-bound Hg within EFPC streambed sediment was relatively small, but could still account for a large proportion of the annual flux of dissolved Hg from EFPC. These sequential extractions also showed that this weakly-bound Hg fraction could be replenished from the much larger fraction of recalcitrant Hg in sediment. The isotopic composition of these weakly-bound and remobilized recalcitrant Hg fractions within the sediment was consistent with high-concentration dissolved Hg hotspots within hyporheic pore water. Thus, this research provided novel evidence that legacy mercury sources within streambed sediment could provide an ongoing contribution of dissolved Hg to surface waters. Finally, we developed new methods for the direct determination of the MeHg isotopic composition of organisms, which allowed a more direct evaluation of inorganic Hg and MeHg sources accumulating in the food web. We found that fish and aquatic invertebrates in both EFPC and a regional background site obtained inorganic Hg and MeHg from multiple isotopically distinct sources, including sediment, suspended particulates, and periphyton. Photodemethylation was found to be an important reaction influencing MeHg dynamics at both sites. However, the balance of microbial methylation and demethylation processes differed between the two streams, with fractionation resulting from methylation and demethylation processes being relatively in balance within the regional background site, whereas microbial methylation appeared to be dominant over microbial demethylation within the EFPC ecosystem. Broadly, the application of Hg isotopic analysis in this study led to numerous novel insights regarding the biogeochemical cycling of Hg in stream ecosystems. This research project promoted the development of two new approaches, including the coupling of sequential extractions with Hg isotopic analysis to assess remobilization of recalcitrant Hg within sediments, and a new method for the isotopic analysis of MeHg isolated from environmental samples. Both of these efforts represent advances in capacity for the field of mercury isotopic analysis and environmental assessment. Overall, this study demonstrates that the application of Hg stable isotope techniques continues to provide new insights into the biogeochemical cycling of Hg in complex aquatic environments, both within the EFPC and beyond.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Geophysical Imaging of Geothermal Systems Spanning Various Geologic Settings

Despite the wide range of geologic settings in which geothermal systems reside, many systems share characteristics that can be targeted by geophysical techniques. Over the last few years, the authors have used potential field and magnetotelluric (MT) methods to jointly image and characterize numerous geothermal systems in various geologic and tectonic settings. This includes systems in the Basin and Range, Washington Cascade Range, Northeastern Oregon, and the Snake River Plain. The principal geothermal targets are areas of elevated temperatures near fault interactions that promote permeability and fluid flow. Gravity, magnetic, and MT data are all sensitive to different physical parameters that together form a complementary set of tools for characterizing subsurface structure and identifying favorable geothermal targets. In general, potential fields delineate important geologic structures and MT locates fluids and zones of hydrothermal alteration. Joint interpretation of these complementary methods, along with other geologic information, provides a robust framework for geothermal prospecting and understanding a hydrothermal system. Specifically, common geothermal targets lie along potential field gradients and are associated with an electrically conductive anomaly that resides in the less-dense, less-magnetic material.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

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↗
DOE OSTI2014

Multiphase Fluid Flow in Deformable Variable-Aperture Fractures (Final Report, February 1, 2009 – January 31, 2014)

Fractures provide flow paths that can potentially lead to fast migration of fluids or contaminants. A number of energy-­related applications involve fluid injections that significantly perturb both the pressures and chemical composition of subsurface fluids. These perturbations can cause both mechanical deformation and chemical alteration of host rocks with potential for significant changes in permeability. In fractured rock subjected to coupled chemical and mechanical stresses, it can be difficult to predict the sign of permeability changes, let alone the magnitude. This project integrated experimental and computational studies to improve mechanistic understanding of these coupled processes and develop and test predictive models and monitoring techniques. The project involved three major components: (1) study of two-­phase flow processes involving mass transfer between phases and dissolution of minerals along fracture surfaces (Detwiler et al., 2009; Detwiler, 2010); (2) study of fracture dissolution in fractures subjected to normal stresses using experimental techniques (Ameli, et al., 2013; Elkhoury et al., 2013; Elkhoury et al., 2014) and newly developed computational models (Ameli, et al., 2014); (3) evaluation of electrical resistivity tomography (ERT) as a method to detect and quantify gas leakage through a fractured caprock (Breen et al., 2012; Lochbuhler et al., 2014). The project provided support for one PhD student (Dr. Pasha Ameli; 2009-­2013) and partially supported a post-­doctoral scholar (Dr. Jean Elkhoury; 2010-­2013). In addition, the project provided supplemental funding to support collaboration with Dr. Charles Carrigan at Lawrence Livermore National Laboratory in connection with (3) and supported one MS student (Stephen Breen; 2011-­2013). Major results from each component of the project include the following: (1) Mineral dissolution in fractures occupied by two fluid phases (e.g., oil-­water or water-­CO 2 ) causes changes in local capillary forces and redistribution of fluids. These coupled processes enhance channel formation and the potential for development of fast flow paths through fractures. (2) Dissolution in fractures subjected to normal stress can result in behaviors ranging from development of dissolution channels and rapid permeability increases to fracture healing and significant permeability decreases. The timescales associated with advective transport of dissolved ions in the fracture, mineral dissolution rates, and diffusion within the adjacent porous matrix dictate the sign and magnitude of the resulting permeability changes. Furthermore, a high--resolution mechanistic model that couples elastic deformation of contacts and aperture-­dependent dissolution rates predicts the range of observed behaviors reasonably well. (3) ERT has potential as a tool for monitoring gas leakage in deep formations. Using probabilistic inversion methods further enhances the results by providing uncertainty estimates of inverted parameters.

02 PETROLEUM↗
DOE OSTI2021

Subsurface Characterization and Machine Learning Predictions at Brady Hot Springs: Preprint

Subsurface data analysis, reservoir modeling, and machine learning (ML) techniques have been applied to the Brady Hot Springs (BHS) geothermal field in Nevada, USA to further characterize the subsurface and assist with optimizing reservoir management. Hundreds of reservoir simulations have been conducted in TETRAD-G and CMG STARS to explore different injection and production fluid flow rates and allocations and to develop a training data set for ML. This process included simulating the historical injection and production since 1979 and prediction of future performance through 2040. ML networks were created and trained using TensorFlow based on multilayer perceptron (MLP), long short-term memory (LSTM), and convolutional neural network (CNN) architectures. These networks took as input selected flow rates, injection temperatures, and historical field operation data and produced estimates of future production temperatures. This approach was first successfully tested on a simplified single fracture doublet system, followed by the application to the BHS reservoir. Using an initial BHS dataset with 37 simulated scenarios, the trained and validated network predicted the production temperature for 6 production wells with the mean absolute percentage error of less than 8%. In a complementary analysis effort, the principal component analysis applied to 13 BHS geological parameters revealed that vertical fracture permeability shows the strongest correlation with fault density and fault intersection density. A new BHS reservoir model was developed considering the fault intersection density as proxy for permeability. This new reservoir model helps to explore under-exploited zones in the reservoir. A data gathering plan to obtain additional subsurface data was developed; it includes temperature surveying for three idle injection wells, at which the reservoir simulations indicate high bottom-hole temperatures. The collected data assist with calibrating the reservoir model and may lead to converting these wells to producers to access under-exploited zones in the reservoir. Data gathering activities are planned for the first quarter of 2021.

40 EE - Geothermal Technologies Office (EE-4G)↗
DOE OSTI2021

2019 Groundwater Monitoring Report Project Shoal Area: Subsurface Corrective Action Unit 447

The Project Shoal Area is a site in Nevada where an underground nuclear test was conducted in 1963. It later came to be known as the Shoal, Nevada, Site. Surface contamination at the site has been remediated, but investigation of groundwater contamination resulting from the test is still in the corrective action process. Annual sampling and water-level monitoring are conducted as part of the subsurface corrective action strategy, which has focused on revising the site conceptual model and evaluating the adequacy of the monitoring well network. It has also included enhancements to the monitoring well network to address uncertainties in the groundwater flow direction and the cause of rising water levels in site wells west of the shear zone since the first hydrologic characterization (HC) wells were installed in 1996. Revisions to the site conceptual model and enhancements to the monitoring strategy were provided to Nevada Division of Environmental Protection (NDEP) in the Addendum to: Corrective Action Decision Document/Corrective Action Plan (CADD/CAP) for the Subsurface Corrective Action Unit 447 Shoal, Nevada, Site. NDEP approved the addendum to the CADD/CAP, which included an expanded contaminant boundary and compliance boundary for the site.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Reversible Fe(II) uptake/release by magnetite nanoparticles

Magnetite commonly coexists with aqueous Fe2+ (Fe2+(aq)) in anoxic subsurface environments. Complex interactions between magnetite and Fe2+(aq) profoundly impact redox potential fluctuations in surrounding environment and biogeochemical cycles of important elements and contaminants. However, the ability of magnetite to act as a source/sink of electron equivalents through fluctuations in solution pH or the activity of Fe2+(aq) remains poorly quantified. We systematically studied the interrelationships between equilibrium Fe2+(aq) concentrations and structural versus surface-localized Fe?II)/Fe?III) ratios in magnetite using micro X-ray diffraction and synchrotron-based X-ray magnetic circular dichroism, respectively, under different controlled experimental conditions. Relative to pH 7, at pH 6 proton-promoted dissolution yields Fe2+(aq) release from magnetite nanoparticles, coupled to a decrease in the structural Fe?II)/Fe?III) ratio by electron hopping along the octahedral sublattice from the particle interior to the surface. At pH 8, magnetite sorbs Fe2+(aq), increasing both the structural and surface-localized Fe?II)/Fe?III) ratio. Amendments of Fe2+(aq) inhibit acidic Fe2+(aq) release and promote Fe2+(aq) uptake at more basic conditions, whereas increasing magnetite loading facilitates Fe2+(aq)–magnetite interaction at the same respective pH extremes. The reversible flow of Fe?II) across the magnetite–solution interface under different conditions implies that the redox reactivity of magnetite nanoparticles is quickly responsive to changes in environmental conditions, such as an increase in pH due to groundwater passing through carbonate-rich rocks, via a dynamic redistribution of electron equivalents between particle interiors and the solid/water interface.

Peng, Huan↗
DOE OSTI2012

A Microsoft Excel Program for Two-Well Tracer Test Data Analysis

Two-well tracer tests are often conducted to investigate subsurface solute transport in the field. Analyzing breakthrough curves in the extraction and monitoring wells using numerical methods is nontrivial due to highly nonuniform flow conditions. We extended, and implemented analytical solutions for the convection-dispersion equation for an injection-extraction well-duplet in a homogeneous confined aquifer under steady state conditions. Functions were provided to calculate the concentrations in the extraction and monitoring wells at any location due to a step or pulse injection. Type curves for a step injection were compared with those calculated by numerically integrating the solution for a pulse injection. The results from the two approaches are similar when the dispersivity is small. As the dispersivity increases, the latter was found to be more accurate but requires more computing time. The code was verified by comparing the results with published type curves and applied to analyze data from the literature. It can be used as a first approximation for two-well tracer test data analysis, and to check accuracy of numerical solutions. The code and example files are publically-available.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2013

Analysis of the Monitoring Network at the Salmon, Mississippi, Site

The Salmon site in southern Mississippi was the location of two underground nuclear tests and two methane-oxygen gas explosion tests conducted in the Tatum Salt Dome at a depth of 2,715 feet below ground surface. The U.S. Atomic Energy Commission (a predecessor agency of the U.S. Department of Energy [DOE]) and the U.S. Department of Defense jointly conducted the tests between 1964 and 1970. The testing operations resulted in surface contamination at multiple locations on the site and contamination of shallow aquifers. No radionuclides from the nuclear tests were released to the surface or to groundwater, although radionuclide-contaminated drill cuttings were brought to the surface during re-entry drilling. Drilling operations generated the largest single volume of waste materials, including radionuclide-contaminated drill cuttings and drilling fluids. Nonradioactive wastes were also generated as part of the testing operations. Site cleanup and decommissioning began in 1971 and officially ended in 1972. DOE conducted additional site characterization between 1992 and 1999. The historical investigations have provided a reasonable understanding of current surface and shallow subsurface conditions at the site, although some additional investigation is desirable. For example, additional hydrologic data would improve confidence in assigning groundwater gradients and flow directions in the aquifers. The U.S. Environmental Protection Agency monitored groundwater at the site as part of its Long-Term Hydrologic Monitoring Program from 1972 through 2007, when DOE's Office of Legacy Management (LM) assumed responsibility for site monitoring. The current monitoring network consists of 28 monitoring wells and 11 surface water locations. Multiple aquifers which underlie the site are monitored. The current analyte list includes metals, radionuclides, and volatile organic compounds (VOCs).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Environmental ISR (Quarterly Progress Report Q3 FY19 April-June 2019)

Progress have been good on both projects. On the subsurface radionuclide project, graduate student, Emily Tran, submitted our first joint manuscript resulting from this collaboration to Water Resources. The manuscript covers batch and flow through experiments using radionuclides with and without bentonite colloids through a naturally fractured chalk rock. Solutions used were collected from a field site in Israel with high ionic strength and filtered prior to using in the experiment. Results show that the bentonite colloids are unstable under these high ionic strengths in carbonate rock in Israel. The U is not significantly eluted with or without bentonite colloids and cerium precipitates as a mobile colloid in the absence of colloids.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2015

Developing a 3-D Seismic-Attribute Framework Model of Yucca Flat, Nevada National Security Site

Using data collected during underground nuclear testing in Yucca Flat at the Nevada National Security Site, an existing 3-D hydrostratigraphic framework model (HFM) that was developed for groundwater modeling is being converted to a 3-D seismic-attribute framework model (SFM). The resulting SFM can be used to better understand seismic wave propagation from underground chemical explosive tests planned for Yucca Flat as part of the Source Physics Experiment. The Yucca Flat HFM covers 1,250 square kilometers and depicts the distribution of subsurface geologic units, called hydrostratigraphic units, according to their ability to transmit groundwater. Many of the rock properties and geologic structures that control groundwater flow also influence seismic wave transmission, so the team reasoned that the HFM could be modified to create an SFM. Existing geologic and geophysical data, including interval velocity measurements from down-hole seismic surveys conducted in 188 holes across Yucca Flat, were used to evaluate the HFM and develop a 3-D velocity property model. The velocity model was integrated with the HFM so that unit volumes within the HFM could be evaluated according to their seismic character. Results to date show, as anticipated, that most of the 3-D hydrostratigraphic unit volumes in the HFM can function as seismic-stratigraphic unit volumes. Thus, with relatively minor modifications, the HFM can be converted to an SFM.

seismic, Yucca Flat, model↗
DOE OSTI2019

Three-dimensional geologic mapping to assess geothermal potential: examples from Nevada and Oregon

Geologic structure plays an important role in controlling fluid flow in geothermal systems. In particular, very complex structural settings, consisting of many closely spaced and intersecting faults, host many geothermal systems. To elucidate the key geologic factors that affect fault-controlled geothermal circulation, it is critical to precisely characterize the structural and stratigraphic geometries in these complex settings. Here, we present a methodology and the results of 3D geologic analyses of two geothermal systems in the Basin and Range, USA. This methodology is a quantitative and geologically focused technique that can be used to precisely characterize geothermal areas, in a time when future geothermal growth demands increased exploration precision and efficiency. Surficial and subsurface geologic and geophysical data are synthesized in the construction of detailed 3D geologic maps of geothermal areas. Based on these 3D geologic maps, we examine several geologic attributes that control permeability development and geothermal fluid flow along faults. We use the stress state of faults and the distribution of structural discontinuities (i.e., fault intersections and fault terminations) to identify locations of upflow along faults in these geothermal systems. These results and the methodology presented herein are directly applicable to structurally controlled geothermal fields in the Basin and Range and worldwide. As development focus shifts toward blind geothermal resources, integration of precisely characterized subsurface structural information into exploration methods will be increasingly critical to continued growth in geothermal exploration and development.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

High-Resolution Simulations of Geological CO 2 Injection: Application to the SPE11 Benchmark

Geological carbon sequestration (GCS) will play a critical role in decarbonization and in facilitating the transition to clean energy systems. Because CO 2 is highly mobile, ensuring its safe and permanent injection into subsurface geological formations involves monitoring over larger spatial domains and longer time periods than is typical for hydrocarbon reservoirs. This can benefit from simulation tools capable of modeling key CO 2 trapping mechanisms, particularly those optimized for speed and scalability on high-performance computing systems. Using isothermal versions of the SPE11B and SPE11C benchmark cases, we conduct a mesh refinement study simulating CO 2 injection into kilometer-scale rock formations at centimeter resolution with the GEOS open-source simulation framework. We focus on how mesh refinement improves the accuracy of convective mixing in both 2D and 3D simulations. The computational costs associated with achieving a converged solution highlight the need for predictive upscaling techniques. A systematic performance scaling analysis—including both central processing unit (CPU) and graphics processing unit (GPU) architectures—complements the “Results” section.

Geosciences↗
DOE OSTI2011

Hydrological parameter estimations from a conservative tracer test with variable-density effects at the Boise Hydrogeophysical Research Site

Reliable predictions of groundwater flow and solute transport require an estimation of the detailed distribution of the parameters (e.g., hydraulic conductivity, effective porosity) controlling these processes. However, such parameters are difficult to estimate because of the inaccessibility and complexity of the subsurface. In this regard, developments in parameter estimation techniques and investigations of field experiments are still challenging and necessary to improve our understanding and the prediction of hydrological processes. Here we analyze a conservative tracer test conducted at the Boise Hydrogeophysical Research Site in 2001 in a heterogeneous unconfined fluvial aquifer. Some relevant characteristics of this test include: variable-density (sinking) effects because of the injection concentration of the bromide tracer, the relatively small size of the experiment, and the availability of various sources of geophysical and hydrological information. The information contained in this experiment is evaluated through several parameter estimation approaches, including a grid-search-based strategy, stochastic simulation of hydrological property distributions, and deterministic inversion using regularization and pilot-point techniques. Doing this allows us to investigate hydraulic conductivity and effective porosity distributions and to compare the effects of assumptions from several methods and parameterizations. Our results provide new insights into the understanding of variabledensity transport processes and the hydrological relevance of incorporating various sources of information in parameter estimation approaches. Among others, the variable-density effect and the effective porosity distribution, as well as their coupling with the hydraulic conductivity structure, are seen to be significant in the transport process. The results also show that assumed prior information can strongly influence the estimated distributions of hydrological properties.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Computational Analysis of Coupled Geoscience Processes in Fractured and Deformable Media

Prediction of flow, transport, and deformation in fractured and porous media is critical to improving our scientific understanding of coupled thermal-hydrological-mechanical processes related to subsurface energy storage and recovery, nonproliferation, and nuclear waste storage. Especially, earth rock response to changes in pressure and stress has remained a critically challenging task. In this work, we advance computational capabilities for coupled processes in fractured and porous media using Sandia Sierra Multiphysics software through verification and validation problems such as poro-elasticity, elasto-plasticity and thermo-poroelasticity. We apply Sierra software for geologic carbon storage, fluid injection/extraction, and enhanced geothermal systems. We also significantly improve machine learning approaches through latent space and self-supervised learning. Additionally, we develop new experimental technique for evaluating dynamics of compacted soils at an intermediate scale. Overall, this project will enable us to systematically measure and control the earth system response to changes in stress and pressure due to subsurface energy activities.

58 GEOSCIENCES↗