Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “subsurface system”

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

Tidal frequencies and quasiperiodic subsurface water level variations dominate redox dynamics in a salt marsh system

Salt marshes are hotspots of nutrient processing en route to sensitive coastal environments. Whilst our understanding of these systems has improved over the years, we still have limited knowledge of the spatiotemporal variability of critical biogeochemical drivers within salt marshes. Sea-level rise will continue to force change on salt marsh functioning, highlighting the urgency of filling this knowledge gap. Our study was conducted in a central California estuary experiencing extensive marsh drowning and relative sea-level rise, making it a model system for such an investigation. Here we instrumented three marsh positions subjected to different degrees of tidal inundation (6.7%, 8.9%, and 11.2% of the time for the upper, middle, and lower marsh positions, respectively), providing locations with varied biogeochemical characteristics and hydrological interactions at the site. We continuously monitored redox potential (Eh) at depths of 0.1, 0.3, and 0.5 m, subsurface water levels (WL), and temperature at 0.7 m depth at each marsh position. To understand how drivers of subsurface biogeochemical processes fluctuate across tidal cycles, we used wavelet analyses to explain the interactions between Eh and WL. We found that tidal forcing significantly affects key drivers of biogeochemical processes by imparting controls on Eh variability, likely driving subsurface hydro-biogeochemistry of the salt marsh. Wavelet coherence showed that the Eh-WL relationship is nonlinear, and their lead–lag relationship is variable. We found that precipitation events perturb Eh at depth over timescales of hours, even though WL shows relatively minimal change during events. Furthermore, this work highlights the importance of high frequency in situ measurements, such as Eh, to help explain factors that govern subsurface biogeochemistry and hydrological processes in salt marshes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Hydropower Geotechnical Foundations: Current Practice and Innovation Opportunities for Low-Head Applications

Hydropower is a renewable energy resource that produces electricity from flowing water under pressure. Engineered hydropower structures, such as dams, are used to create a hydraulic head, enabling a turbinegenerator unit to convert pressurized flow into electricity. While hydropower has been a source of renewable energy since antiquity, new development in the United States has slowed in recent decades. Based on recent resource assessments, the largest opportunity to expand hydropower in the United States is from new stream-reach development (i.e., new hydropower development along stream-reaches that do not currently have hydroelectric facilities or other forms of infrastructure, such as dams). Roughly 75% of identified new stream-reach development potential is from low-head sites (less than 30 ft of head), which typically suffer from smaller power densities and higher normalized costs, given economies of scale. Hydropower developers and other stakeholders are thus interested in strategies to reduce initial capital costs while practicing sustainable development to maximize environmental compatibility with minimal disruption to natural aquatic life, sediment, and water flows. Historically, civil works have represented a significant cost driver for new hydropower development, with the foundation system representing a major cost component and source of uncertainty. The foundation system is the collection of engineered structural features (e.g. cutoff trenches, walls, grouting, anchors) constructed at or below the preconstruction ground surface that interfaces between the overlying structures (superstructures) and the bed material below (subsurface). Development of a hydropower foundation system must consider the various characteristics of the surrounding stream environment and subsurface while adhering to the engineering requirements of the superstructure that it supports. The care of water, excavation, and other construction activities are important features of foundation design and construction. The design and construction cost of the foundation system is largely dependent on the site geology and riverbed composition and is influenced by the level of geotechnical assessment required and conducted. Thus, a hydropower facility’s geotechnical foundation is often highly site-specific, with proper site selection and assessment being important to project success. The foundation system is designed to provide structural stability (of the foundation and dam), limit seepage, ensure public safety, and maintain functionality for the project life, during both construction and facility operations. Inadequate foundation or dam design can result in dam failure and the uncontrolled release of significant volumes of water, which could cause a high number of casualties and extensive property damage downstream of the failure. According to the Association of State Dam Safety Officials, approximately 30% of all historical dam failures in the United States are attributed to foundation or abutment defects, and another 20% are attributed to piping or seepage through the embankment, foundation, or abutment. To ameliorate these safety considerations, foundations often require massive amounts of construction material (e.g., grout, concrete, engineered dam fill) and long construction times. Foundation design also requires significant analysis prior to construction because the initial in-stream and abutment subsurface conditions are site-specific, and sufficient data for them often are lacking. Current practice requires on-site assessment, using expensive drilling and invasive and non-invasive investigation methods, to determine the expected cost of foundation material and treatment. Additionally, foundation construction often requires site dewatering (and other care of water activities), which involves constructing temporary diversion structures upstream and often downstream, called cofferdams, and water diversion systems that route water around the construction site. Cofferdams and water diversion systems can drastically increase construction costs and contribute to environmental disruption, including modification of flow patterns and benthic habitats. Given the technical, economic, and environmental challenges associated with hydropower foundations, opportunities exist to improve the current state of practice and to develop new and innovative solutions to Hydropower is a renewable energy resource that produces electricity from flowing water under pressure. Engineered hydropower structures, such as dams, are used to create a hydraulic head, enabling a turbinegenerator unit to convert pressurized flow into electricity. While hydropower has been a source of renewable energy since antiquity, new development in the United States has slowed in recent decades. Based on recent resource assessments, the largest opportunity to expand hydropower in the United States is from new stream-reach development (i.e., new hydropower development along stream-reaches that do not currently have hydroelectric facilities or other forms of infrastructure, such as dams). Roughly 75% of identified new stream-reach development potential is from low-head sites (less than 30 ft of head), which typically suffer from smaller power densities and higher normalized costs, given economies of scale. Hydropower developers and other stakeholders are thus interested in strategies to reduce initial capital costs while practicing sustainable development to maximize environmental compatibility with minimal disruption to natural aquatic life, sediment, and water flows. Historically, civil works have represented a significant cost driver for new hydropower development, with the foundation system representing a major cost component and source of uncertainty. The foundation system is the collection of engineered structural features (e.g. cutoff trenches, walls, grouting, anchors) constructed at or below the preconstruction ground surface that interfaces between the overlying structures (superstructures) and the bed material below (subsurface). Development of a hydropower foundation system must consider the various characteristics of the surrounding stream environment and subsurface while adhering to the engineering requirements of the superstructure that it supports. The care of water, excavation, and other construction activities are important features of foundation design and construction. The design and construction cost of the foundation system is largely dependent on the site geology and riverbed composition and is influenced by the level of geotechnical assessment required and conducted. Thus, a hydropower facility’s geotechnical foundation is often highly site-specific, with proper site selection and assessment being important to project success. The foundation system is designed to provide structural stability (of the foundation and dam), limit seepage, ensure public safety, and maintain functionality for the project life, during both construction and facility operations. Inadequate foundation or dam design can result in dam failure and the uncontrolled release of significant volumes of water, which could cause a high number of casualties and extensive property damage downstream of the failure. According to the Association of State Dam Safety Officials, approximately 30% of all historical dam failures in the United States are attributed to foundation or abutment defects, and another 20% are attributed to piping or seepage through the embankment, foundation, or abutment. To ameliorate these safety considerations, foundations often require massive amounts of construction material (e.g., grout, concrete, engineered dam fill) and long construction times. Foundation design also requires significant analysis prior to construction because the initial in-stream and abutment subsurface conditions are site-specific, and sufficient data for them often are lacking. Current practice requires on-site assessment, using expensive drilling and invasive and non-invasive investigation methods, to determine the expected cost of foundation material and treatment. Additionally, foundation construction often requires site dewatering (and other care of water activities), which involves constructing temporary diversion structures upstream and often downstream, called cofferdams, and water diversion systems that route water around the construction site. Cofferdams and water diversion systems can drastically increase construction costs and contribute to environmental disruption, including modification of flow patterns and benthic habitats. Given the technical, economic, and environmental challenges associated with hydropower foundations, opportunities exist to improve the current state of practice and to develop new and innovative solutions to challenges frequently encountered with traditional approaches. With this understanding, it is critically important to understand and document the current state of practice for hydropower geotechnical foundations, identify key challenges, and define opportunities for innovative solutions. To this end, this report documents the current state of practice across the three main phases of geotechnical foundation development: (1) geotechnical site assessment, (2) design, and (3) construction for hydropower systems. It also describes the major challenges with conventional approaches and identifies opportunities for innovation to reduce hydropower foundations costs, timelines, and risks. Key takeaways from this report include the following: Approximately 80% of available low-head sites are expected to have foundations on soil beds rather than rock beds, suggesting that rockfill and earthfill dams may be the most cost-effective conventional dam type for new projects.; Geotechnical and geologic investigation activities are time-consuming and expensive but are essential to define the parameters and criteria needed for foundation design.; Certain riverbed soil and bedrock types present significant technical challenges or require expensive foundation construction, which can prove financially prohibitive for low-head project development.; Modular hydropower design and prefabricated modular foundations represent a promising but unproven paradigm for new hydropower development. Design and construction approaches using optimized and highly repeatable, reliable components would benefit project cost, time, and risk but require additional research and development.; Temporary construction features for foundations, including cofferdams, water diversion, and water control systems, can prove costly and have inherent construction risk.; For economically viable development, hydropower geotechnical foundations should be limited to 4 to 15% of the project’s total initial capital costs. Many proposed projects have experienced cost overruns attributable to foundation difficulties or surprises during construction. These overruns may have been due to inadequate investigations, lack of adequate engineering effort to tailor the structures to site geology and topography, and/or contractual terms, among other considerations.; Challenges for hydropower foundations and opportunities for innovative technology solutions are identified in the following areas (consistent with the three main phases of foundation development): Geotechnical site assessment, Foundation design and materials, Construction methods and technology. Ultimately, this report aims to provide information about geotechnical foundations for low-head hydropower and to motivate transformative technologies to support hydropower growth.

13 HYDRO ENERGY↗
DOE OSTI2023

Low-order aquifer thermal energy storage model for geothermal system simulation

This paper presents a low-order aquifer thermal energy storage (ATES) model for simulation of combined subsurface and above-surface energy systems. The model is included in the Modelica IBPSA Library, which is a free open-source library with basic models for building and district energy and control systems. The model uses a lumped-component method, in which the transient conductive-convective heat and mass transfer equation is radially discretized. To verify the accuracy of the model, we present an inter-model comparison from a simulation test suite. Results show that the Modelica ATES model is in good agreement, with a normalized mean bias error for yearly variation of aquifer temperatures of 1.6×10−2 and 9×10−5 at 1 m and 10 m distance from the well.

Maccarini, Alessandro↗
DOE OSTI2022

Unsupervised learning from three-component accelerometer data to monitor the spatiotemporal evolution of meso-scale hydraulic fractures

Enhanced geothermal systems can provide a substantial share of the global energy demand. There exist several hurdles in the engineering implementations of such geothermal systems. One such hurdle is the accurate monitoring of the fracture networks created in subsurface through hydraulic stimulation of these systems. Micro seismicity associated with the stimulation is the primary means to locate the event hypocenters for estimating the stimulated rock volume. Existing methods for location the hypocenters are restricted to only the highest amplitude impulsive signals that are simultaneously detected on several sensors. Consequently, a large portion (usually ~99%) of the measurements are left unused. In this paper, an unsupervised manifold-approximation followed by clustering of 3-component accelerometer data is used to analyze the seismicity recorded on a monitoring well. With this method, a larger portion of the measured signal is used for the monitoring of the hydraulic fracture network. We analyze the EGS Collab experiment 1 microseismic data, recorded at the Sanford Underground Research Facility, South Dakota. Using the data from a single three-component accelerometer, the polarization features viz. Azimuth, incidence, rectilinearity, and planarity are used as inputs for the unsupervised manifold approximation followed by clustering. Our study shows that density-based clusters in the projected 3D space correspond to distinct types of hydraulically fractured zones around the injection point. Finally, we show that the temporal evolution of these clusters can be used to track fracture creation and propagation.

58 GEOSCIENCES↗
DOE OSTI2025

The rigorous upscaling of advection-dominated transport in heterogeneous porous media via the Method of Finite Averages

Systems involving advection-dominated transport through heterogeneous porous and fractured media are ubiquitous in subsurface engineering applications. However, upscaling such systems continues to challenge rigorous modeling efforts, particularly when advection is stronger than diffusion at fine spatial scales (i.e., when the Péclet number is greater than one at length scales that characterize a system’s unit-cells, representative elementary volumes, or averaging regions). Here, in this work, we propose and validate a strategy for extending the Method of Finite Averages (MoFA), a rigorous upscaling methodology for heterogeneous porous media, to upscale transport systems experiencing stronger advection than diffusion at fine scales (i.e., fine-scale Péclet numbers greater than one). We detail the strategy, the physical conditions under which it can be applied while retaining a priori modeling error guarantees, and implement the strategy to obtain a MoFA model for advective-diffusive transport that accommodates advective physics at fine spatial scales. We then perform two numerical experiments considering systems with system-scale Péclet numbers of 300 and 1000 — which correspond to fine-scale Péclet numbers of 30 and 100, respectively — to verify that the error guarantees are met under the strategy. After, we conduct a numerical study to demonstrate the strategy’s advantages over the original MoFA methodology. The results suggest that rigorously-upscaled transport models for heterogeneous porous media experiencing advective physics at finer spatial scales can be derived through MoFA and resolved orders of magnitude faster than their pore-scale counterparts. The results also suggest that the presented strategy is limited to modeling shallow concentration gradients when there are large differences between the time scales related to advection and a system’s temporally-varying boundary conditions. This limitation hinders the strategy’s practicality in modeling more advective systems, and as such, opportunity exists for developing additional strategies that accommodate rapidly-varying boundary conditions — and consequentially, steeper concentration gradients — while modeling advective systems with MoFA.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Versatile soil gas concentration and isotope monitoring: optimization and integration of novel soil gas probes with online trace gas detection

Abstract. Gas concentrations and isotopic signatures can unveil microbial metabolisms and their responses to environmental changes in soil. Currently, few methods measure in situ soil trace gases such as the products of nitrogen and carbon cycling or volatile organic compounds (VOCs) that constrain microbial biochemical processes like nitrification, methanogenesis, respiration, and microbial communication. Versatile trace gas sampling systems that integrate soil probes with sensitive trace gas analyzers could fill this gap with in situ soil gas measurements that resolve spatial (centimeters) and temporal (minutes) patterns. We developed a system that integrates new porous and hydrophobic sintered polytetrafluoroethylene (sPTFE) diffusive soil gas probes that non-disruptively collect soil gas samples with a transfer system to direct gas from multiple probes to one or more central gas analyzer(s) such as laser and mass spectrometers. Here, we demonstrate the feasibility and versatility of this automated multiprobe system for soil gas measurements of isotopic ratios of nitrous oxide (δ18O, δ15N, and the 15N site preference of N2O), methane, carbon dioxide (δ13C), and VOCs. First, we used an inert silica matrix to challenge probe measurements under controlled gas conditions. By changing and controlling system flow parameters, including the probe flow rate, we optimized recovery of representative soil gas samples while reducing sampling artifacts on subsurface concentrations. Second, we used this system to provide a real-time window into the impact of environmental manipulation of irrigation and soil redox conditions on in situ N2O and VOC concentrations. Moreover, to reveal the dynamics in the stable isotope ratios of N2O (i.e., 14N14N16O, 14N15N16O, 15N14N16O, and 14N14N18O), we developed a new high-precision laser spectrometer with a reduced sample volume demand. Our integrated system – a tunable infrared laser direct absorption spectrometry (TILDAS) in parallel with Vocus proton transfer reaction mass spectrometry (PTR-MS), in line with sPTFE soil gas probes – successfully quantified isotopic signatures for N2O, CO2, and VOCs in real time as responses to changes in the dry–wetting cycle and redox conditions. Broadening the collection of trace gases that can be monitored in the subsurface is critical for monitoring biogeochemical cycles, ecosystem health, and management practices at scales relevant to the soil system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Study of an Unrefined Humate Solution as a Possible Attenuation-based Remedy for Uranium Contamination in Acidic Groundwater - 20351

The Savannah River Site (SRS), joined the environmental cleanup program in 1981 after the Cold War, and was designated as a hazardous waste management facility. About 1.8 billion gallons of acid wastes were disposed into the F/H Area Seepage Basin that included many radionuclides and dissolved metals, resulting in highly contaminated groundwater plumes with pH of 3 - 5.5. The acidity of the plumes contributes to the mobility of several constituents of concern (COC) such as H-3, U-238, I-129, and Sr-90 for the F-Area plume and H-3, Sr-90 and mercury for the H-Area plume. An active treatment (pump-and-treat and re-injection) system was constructed and put in place in 1997 to address the removal of metals and radionuclides from the contaminated groundwater. As remediation projects advanced, active treatments were transitioned to more passive or enhanced-passive approaches such as the subsurface barrier with base injection system, which has been controlling the mobility of contaminants. However, new and more efficient attenuation-based remedies are always sought. Several studies have proposed that humic substances can be used to remediate sites contaminated with heavy metals by creating a permeable reactive barrier. Humic substances (HS) are major components of soil organic matter. HS are polyfunctional organic macromolecules formed by the chemo-microbiological decomposition of biomass or dead organic matter. These substances are usually divided into three main fractions: humin (insoluble at all pH values), humic acid (soluble at pH greater than 3.5), and fulvic acid (soluble at all pH values). Humic substances are major components of the soil at SRS and, in addition, are also helpful in the removal process of uranium. This research focuses on uranium (VI), which is a key contaminant of concern in the F-Area groundwater plume. The interaction of uranium with sediment in the presence or absence of humic substances involves complex mechanisms that are not yet well understood. The interaction of U(VI) with humic substances can affect the adsorption of U on sediment, altering its mobility in the subsurface. The objective of this study is to study the mobility of uranium in the presence of humic acid and to determine if sediments amended by humic acid can enhance the sorption of heavy metals onto sediments. In this research, experiments were conducted using background (clean) F-Area aquifer sediments to understand and predict uranium mobility. A blended material, humate, containing both humic and fulvic acids and a chemically modified humate (KW15 modified humic or mod-HA) was studied as a possible amendment for uranium remediation in SRS groundwater. Batch experiments were conducted in triplicates using 20 ml of aqueous suspension in DI water containing SRS sediment and mod-HA and studied the effect of pH, kinetics of humate sorption onto SRS sediment and also studied the sorption of uranium onto mod-HA amended SRS sediment. Results indicate that mod-HA amended sediment increases uranium removal significantly. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Predicting Flow in Fracture Networks With Quantum Algorithms

Uncertainty quantification plays a crucial role in the modeling of subsurface flow. For instance, uncertainties in the properties of geologic fracture networks significantly impact flow, requiring numerous simulations to accurately estimate quantities of interest. However, each simulation is computationally expensive because it requires solving a large linear system to capture features that involve both small and large fractures. An example is in percolation, where the interaction of many small fractures (which cumulatively can have a large surface area) with the rock matrix must be modeled precisely. Quantum computing is an emerging tool with the potential to address this issue. Quantum algorithms offer a significant speedup in solving linear systems, achieving efficiencies that are challenging to match with classical approaches. These classical approaches include direct solvers, such as LU decomposition, and iterative methods, notably preconditioned conjugate gradient, commonly used in subsurface modeling to solve large sparse systems. However, applying quantum algorithms to geologic fracture flow requires careful attention to algorithmic and problem-specific constraints to fully realize this quantum advantage. In this work we describe a quantum algorithm for generalized Monte Carlo applications with a quadratic speedup over the classical approaches which can be combined with the quantum speedup, currently under investigation, for solving quantum linear systems for subsurface flow. We show that for quantum algorithms the computational cost of estimating a quantity of interest for a statistical ensemble of networks is roughly the same as that of a single realization, essentially implying that one can get uncertainty quantification for free.

58 GEOSCIENCES↗
DOE OSTI2024

Aromatic amino acid metabolism and active transport regulation are implicated in microbial persistence in fractured shale reservoirs

Abstract Hydraulic fracturing has unlocked vast amounts of hydrocarbons trapped within unconventional shale formations. This large-scale engineering approach inadvertently introduces microorganisms into the hydrocarbon reservoir, allowing them to inhabit a new physical space and thrive in the unique biogeochemical resources present in the environment. Advancing our fundamental understanding of microbial growth and physiology in this extreme subsurface environment is critical to improving biofouling control efficacy and maximizing opportunities for beneficial natural resource exploitation. Here, we used metaproteomics and exometabolomics to investigate the biochemical mechanisms underpinning the adaptation of model bacterium Halanaerobium congolense WG10 and mixed microbial consortia enriched from shale-produced fluids to hypersalinity and very low reservoir flow rates (metabolic stress). We also queried the metabolic foundation for biofilm formation in this system, a major impediment to subsurface energy exploration. For the first time, we report that H. congolense WG10 accumulates tyrosine for osmoprotection, an indication of the flexible robustness of stress tolerance that enables its long-term persistence in fractured shale environments. We also identified aromatic amino acid synthesis and cell wall maintenance as critical to biofilm formation. Finally, regulation of transmembrane transport is key to metabolic stress adaptation in shale bacteria under very low well flow rates. These results provide unique insights that enable better management of hydraulically fractured shale systems, for more efficient and sustainable energy extraction.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2026

Comparative physiological and genomic characterization of a novel Nitrobacter vulgaris strain from a nitrate-contaminated subsurface

Nitrite-oxidizing bacteria (NOB) represent a crucial node in the global nitrogen cycle. By catalyzing the second step of nitrification—the oxidation of nitrite to nitrate to generate energy for growth—NOB activity controls the fate of nitrite (NO 2 - ) in aerobic environments. Despite thriving in diverse environments, including soils, freshwater, marine ecosystems, subsurface habitats, and water treatment systems, organisms capable of nitrite oxidation are confined to Nitrobacter, Nitrospira, Nitrospina, Nitrotoga, and a few other specific lineages. The genus Nitrobacter, recognized for its facultative heterotrophic metabolism, is often associated with high-nitrogen environments. Here, we report the physiological characterization of a novel strain, Nitrobacter vulgaris strain MLSD-S22, isolated from a nitrate- and heavy-metal-contaminated subsurface. Growth inhibition experiments revealed that strain MLSD-S22 and the N. vulgaris type strain Z exhibited similar sensitivities to nitrite and nitrate, with nitrite being the most inhibitory. Microrespirometry demonstrated that the two N. vulgaris strains and Nitrobacter winogradskyi Nb-255 possessed higher affinities for nitrite and oxygen than previously reported for Nitrobacter, suggesting potential to compete in low-substrate environments. Long-read DNA sequencing provided a complete genome for strain MLSD-S22, revealing two plasmids and an intact nitrous oxide (N 2 O) reduction operon—an unexpected feature for Nitrobacter. While N 2 O reduction activity was not observed under the tested conditions, this discovery raises questions about the contribution of Nitrobacter NOB to the N 2 O sink. These findings broaden the physiological and genomic diversity of Nitrobacter, offering new insights into their adaptation strategies and providing a framework for future evaluation of their potential roles in nitrogen loss.

Nitrobacter↗
DOE OSTI2020

Decovalex-2019 (Executive Summary)

The DECOVALEX Project is an on-going international research collaboration, established in 1992, to advance the understanding and modeling of coupled Thermal (T), Hydrological (H), Mechanical (M) and Chemical (C) processes in geological in geological systems. DECOVALEX was initially motivated by the recognition that prediction of these coupled effects is an essential part of the performance and safety assessment of geologic disposal systems for radioactive waste and spent nuclear fuel. Later it was realized that these processes also play a critical role in other subsurface engineering activities, such as subsurface CO 2 storage, enhanced geothermal systems, and unconventional oil and gas production through hydraulic fracturing. Research teams from many countries (e.g., Canada, China, Czech Republic, Finland, France, Germany, Japan, Republic of Korea, Spain, Sweden, Switzerland, Taiwan, United Kingdom, and the United States) various institutions have participated in the DECOVALEX Project over the years, providing a wide range of perspectives and solutions to these complex problems. These institutions represent radioactive waste management organizations, national research institutes, regulatory agencies, universities, as well as industry and consulting groups. The overall aim of DECOVALEX-2019 was to increase the understanding of various thermo-hydro-mechanical-chemical processes of importance for radionuclide release and transport from a repository to the biosphere and how they can be described and modelled using mathematical models. The scientific and technical objectives are: to increase the basic understanding of T-H-M-C coupled processes in fractured rocks (crystalline, sedimentary, argillaceous) and buffer materials; to investigate the predictive capabilities of different codes to field experiments and to perform verification of codes; to exchange experimental data, and improve the understanding of the constitutive behavior of crystalline and argillaceous rock masses and buffer materials; and to perform THMC calculations in a performance/safety assessment context.

58 GEOSCIENCES↗
DOE OSTI2020

DECOVALEX-2019 (Task E Final Report)

The DECOVALEX Project is an on-going international research collaboration, established in 1992, to advance the understanding and modeling of coupled Thermal (T), Hydrological (H), Mechanical (M) and Chemical (C) processes in geological in geological systems. DECOVALEX was initially motivated by the recognition that prediction of these coupled effects is an essential part of the performance and safety assessment of geologic disposal systems for radioactive waste and spent nuclear fuel. Later it was realized that these processes also play a critical role in other subsurface engineering activities, such as subsurface CO 2 storage, enhanced geothermal systems, and unconventional oil and gas production through hydraulic fracturing. Research teams from many countries (e.g., Canada, China, Czech Republic, Finland, France, Germany, Japan, Republic of Korea, Spain, Sweden, Switzerland, Taiwan, United Kingdom, and the United States) various institutions have participated in the DECOVALEX Project over the years, providing a wide range of perspectives and solutions to these complex problems. These institutions represent radioactive waste management organizations, national research institutes, regulatory agencies, universities, as well as industry and consulting groups. This document is the final report of Task E which was proposed and coordinated by Andra, the National Radioactive Waste Management Agency in France, presenting the technical definitions of the problems studied, approaches applied, achievements made and outstanding issues for future research. The purpose of Task E of the DECOVALEX-2019 project is to investigate upscaling THM modelling from small-scale experiments (some cubic meters) to full-scale experiments (some ten cubic meters) and finally to the scale of the waste repository (cubic kilometers). To achieve this aim, the data of two in-situ heating experiments performed by Andra (the French National Radioactive Waste Management Agency) in the Meuse/Haute-Marne Underground Research Laboratory (MHM URL) have formed the basis for the understanding of the THM behavior of the COx at different scales. The first experiment provided the reference values of the THM parameters by means of a calibration exercise and they were used for a blind prediction and an interpretative analysis of the second one.

58 GEOSCIENCES↗
DOE OSTI2021

Using Complex Probability Amplitudes to Simulate Solute Transport in Composite Porous Media

Probability amplitudes are fundamental to quantum mechanics and offer robust descriptions of complicated systems, which have allowed physicists to explain behaviors inaccessible to classical physics. This article ponders how some of the same conceptual underpinnings of the mathematics used for modeling quantum systems might be applied to subsurface water resources problems and speculates how these tools could facilitate applications on quantum computers. A probability amplitude-based model for describing advective-dispersive transport in porous media using linear operators is investigated. The proposed complex valued model decomposes spreading into two “sub-continuum partial dispersion” coefficients, and this recovers classical spreading when the sum of these coefficients is the Fickian dispersion coefficient. However, the probability amplitudes have a manyto-one relationship to a probability distribution, so it embeds a level of heterogeneity into seemingly equivalent functions. Two propagators with different sub-continuum coefficients may have the same macroscopic behavior when either is considered in isolation, but when they act on the other the system’s behavior changes. Additionally, differences in the amplitudes cause a reduction in spreading as velocity correlations are disrupted, despite both propagators having identical dispersion coefficients, and this cannot be achieved using classical methods without changing the dispersion coefficient. The main point is that these amplitude-based models offer a way to embed information about the system into the propagators, instead of just “averaging it out” when making an upscaled model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Methods and systems for obtaining reconstructed low-frequency seismic data for determining a subsurface feature

A computer-implemented method for obtaining reconstructed seismic data for determining a subsurface feature, includes: determining an initial training velocity model, training a machine learning model based on first training seismic data and second training seismic data generated from the training velocity model, the first training seismic data corresponding to one or more first frequencies, the second training seismic data corresponding to one or more second frequencies lower than the one or more first frequencies, obtaining, based on measured seismic data and the machine learning model, reconstructed seismic data corresponding to the one or more second frequencies, generating a velocity model based on the measured seismic data, the reconstructed seismic data, and a full waveform inversion (FWI), and when the generated velocity model does not satisfy a preset condition, updating the training velocity model based on the generated velocity model, to obtain updated reconstructed seismic data for determining a subsurface feature.

58 GEOSCIENCES↗
DOE OSTI2026

Modeling supercritical CO 2 flow and mineralization in reactive host rocks with PFLOTRAN v7.0

Understanding the flow and reactivity of CO 2 injected into geological reservoirs is important for many subsurface applications including secure geologic carbon storage (GCS), critical mineral extraction, enhanced geothermal systems (EGS), and enhanced oil recovery (EOR). Traditionally, subsurface CO 2 injection for GCS applications has focused on geologic formations with favorable subsurface configurations for CO 2 migration and trapping through non-reactive mechanisms such as structural, solubility, and petrophysical trapping. Recently, CO 2 -reactive rocks such as mafic and ultramafic basalts have been investigated for their potential to react with injected CO 2 in situ to simultaneously dissolve host rock minerals and mineralize CO 2 as carbonates. Engineering rapid CO 2 mineralization in the subsurface is attractive because of the increased density of stored CO 2 , the additional safety factors associated with solidification, and the potential to extract valuable critical minerals. Here we present recent developments in the parallel flow and reactive transport simulator PFLOTRAN to model coupled CO 2 -brine flow and reactive transport for a wide range of injection and production applications involving reactive CO 2 -brine systems. These developments are based on the well established and trusted CO 2 flow capabilities in the STOMP-CO 2 simulator. New capabilities added to PFLOTRAN include new CO 2 -brine equations of state with optional thermal coupling, several new constitutive relationships like capillary pressure smoothing and scanning path hysteresis, a fully implicit well model, and native linkage with PFLOTRAN's well-established reactive transport libraries. A series of benchmarks between PFLOTRAN and STOMP-CO 2 verify the newly developed CO 2 -brine flow capabilities. Demonstrations of coupled CO 2 -brine flow modeling and reactive transport show how CO 2 mineralization can be engineered in reactive host rocks. Finally, an example use case involving copper leaching by CO 2 and critical mineral extraction is presented to showcase the strengths of this new implementation. Several limitations still remain, including limited availability of field data to parameterize models. Future work should constrain the evolution of mineral surface area during mineralization and the temperature and/or pH dependence of geochemical reactions for specific systems of interest.

Critical Minerals↗
DOE OSTI2022

Coupling Subsurface and Above-Surface Models for Optimizing the Design of Borefields and District Heating and Cooling Systems

Accurate dynamic energy simulation is important for the design and sizing of district heating and cooling systems with geothermal heat exchange for seasonal energy storage. Current modeling approaches in building and district energy simulation tools typically consider heat conduction through the ground between boreholes without flowing groundwater. While detailed simulation tools for subsurface heat and mass transfer exist, these fall short in simulating above-surface energy systems. To support the design and operation of such systems, the study developed a coupled model including a software package for building and district energy simulation, and software for detailed heat and mass transfer in the subsurface. For the first, it uses the open-source Modelica Buildings Library, which includes dynamic simulation models for building and district energy and control systems. For the heat and mass transfer in the soil, it uses the TOUGH simulator. The TOUGH family of codes can model heat and multi-phase, multi-component mass transport for a variety of fluid systems, as well as chemical reactions, in fractured porous media. The study validated the coupled modeling approach by comparing the simulation results with one from the g-function based ground response model. It then looked into effects when the water table and the regional groundwater flow are considered in the ground, from the perspective of heat exchange between borehole and ground, and the electrical consumption of the district heating and cooling systems. To access the simulation models, please find the links in the submission: -- For coupled approach validation: see model Buildings.Fluid.Geothermal.Borefields.Examples.BorefieldsWithTough and Buildings.Examples.DistrictReservoirNetworks.Examples.Reservoir3Variable_TOUGH from the "Modelica Building Library" resource, branch issue1495_tough_interface, commit a2667c0. -- For the study of the effect of water table: see model Buildings.Examples.DistrictReservoirNetworks.Examples.Reservoir3Variable_TOUGH from he "Modelica Building Library" resource, branch issue1495_tough_interface_moreIO, commit 760de49. -- For the study of the effect of regional groundwater flow: see Buildings.Examples.DistrictReservoirNetworks.Examples.Reservoir3Variable_TOUGH from he "Modelica Building Library" resource, branch issue1495_tough_interface_moreIO_3D, commit c2a2d2a. The coupling interface script "GrounResponse.py" can be found from the above links in the folder Buildings/Resources/Python-Sources. Also, the needed files for TOUGH simulation are in the folder Buildings/Resources/Python-Sources/ToughFiles that can be accessed through the above links. A brief description of these files is given below; detailed specifications for the first three files may be found in the TOUGH3 Users Guide (Jung et al., 2018) https://tough.lbl.gov/documentation/tough-manuals/. (1) INCON - initial conditions for each grid block (2) INFILE - main input file with material properties and control parameters (3) MESH - description of the computational grid (4) readsave - Modelica/TOUGH interface program: read the final output of TOUGH simulation after TOUGH time step and prepare for transfer to Modelica for next Modelica time step (5) readsave.inp - input parameters for program readsave (6) writeincon - Modelica/TOUGH interface program: write the output of Modelica after Modelica time step and prepare for transfer to TOUGH as initial conditions for the next TOUGH step (7) writeincon.inp - input parameters for program writeincon

15 GEOTHERMAL ENERGY↗