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 577 records · Page 32

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

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

Elkholy, Sherif↗

Mobility of Source Zone Heavy Metals and Radionuclides: The Mixed Roles of Fermentative Activity on Fate and Transport of U and Cr. Final Report

Various U. S. Department of Energy (DOE) low and medium-level radioactive waste sites contain mixtures of heavy metals, radionuclides and assorted organic materials. In addition, there are numerous sites around the world that are contaminated with a mixture of organic and inorganic contaminants. In most sites, over time, water infiltrates the wastes, and releases metals, radionuclides and other contaminants causing transport into the surrounding environment. We investigated the role of fermentative microorganisms in such sites that may control metal, radionuclide and organics migration from source zones. The project was initiated based on the following overarching hypothesis: Metals, radionuclides and other contaminants can be mobilized by infiltration of water into waste storage sites. Microbial communities of lignocellulose degrading and fermenting microorganisms present in the subsurface of contaminated DOE sites can significantly impact migration by directly reducing and immobilizing metals and radionuclides while degrading complex organic matter to low molecular weight organic compounds. These low molecular weight organic acids and alcohols can increase metal and radionuclide mobility by chelation (i.e., certain organic acids) or decrease mobility by stimulating respiratory metal reducing microorganisms. We demonstrated that fermentative organisms capable of affecting the fate of Cr6+, U6+ and trinitrotoluene can be isolated from organic-rich low level waste sites as well as from less organic rich subsurface environments. The mechanisms, pathways and extent of contaminant transformation depend on a variety of factors related to the type of organisms present, the aqueous chemistry as well as the geochemistry and mineralogy. This work provides observations and quantitative data across multiple scales that identify and predict the coupled effects of fermentative carbon and electron flow on the transport of radionuclides, heavy metals and organic contaminants in the subsurface; a primary concern of the DOE Environmental Remediation Science Division (ERSD) and Subsurface Geochemical Research (SBR) Program.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

EGS Collab Project Experiment 1 Overview and Progress

Enhanced Geothermal Systems (EGS) offer the potential to extract and use large quantities of clean energy, but questions remain on reservoir creation and sustainability. The EGS Collab project, supported by the US Department of Energy’s Geothermal Technologies Office, is establishing a suite of highly monitored and well-characterized intermediate-scale (~10-20 m) field test beds along with fracture stimulation and interwell flow tests to better understand processes that control formation of effective subsurface heat exchangers. EGS Collab tests will provide a means of testing tools and concepts that could later be employed under geothermal reservoir conditions at DOE’s Frontier Observatory for Research in Geothermal Energy (FORGE) or enhanced geothermal systems. Key to the project is using numerical simulations in the experiment design and interpretation of results. Our first set of experiments is underway at the Sanford Underground Research Facility (SURF) in South Dakota. To date, stepwise stimulations have been performed at two locations, with the final stimulation connecting our injection and production wells. Numerous data have been collected and are currently being analyzed.

Enhanced Geothermal Systems, EGS Collab, Stimulati↗

Fourier-MIONet: Fourier-enhanced multiple-input neural operators for multiphase modeling of geological carbon sequestration

Geologic carbon sequestration (GCS) is a safety-critical technology that aims to reduce the amount of carbon dioxide in the atmosphere, which also places high demands on reliability. Multiphase flow in porous media is essential to understand CO 2 migration and pressure fields in the subsurface associated with GCS. However, numerical simulation for such problems in 4D is computationally challenging and expensive, due to the multiphysics and multiscale nature of the highly nonlinear governing partial differential equations (PDEs). It prevents us from considering multiple subsurface scenarios and conducting real-time optimization. Here, we develop a Fourier-enhanced multiple-input neural operator (Fourier-MIONet) to learn the solution operator of the problem of multiphase flow in porous media. Fourier-MIONet utilizes the recently developed framework of the multiple-input deep neural operators (MIONet) and incorporates the Fourier neural operator (FNO) in the network architecture. Once Fourier-MIONet is trained, it can predict the evolution of saturation and pressure of the multiphase flow under various reservoir conditions, such as permeability and porosity heterogeneity, anisotropy, injection configurations, and multiphase flow properties. Compared to the enhanced FNO (U-FNO), the proposed Fourier-MIONet has 90% fewer unknown parameters, and it can be trained in significantly less time (about 3.5 times faster) with much lower CPU memory (<15%) and GPU memory (<35%) requirements, to achieve similar prediction accuracy. In addition to the lower computational cost, Fourier-MIONet can be trained with only 6 snapshots of time to predict the PDE solutions for 30 years. Furthermore, we observed that Fourier-MIONet can maintain good accuracy when predicting out-of-distribution (OOD) data. The excellent generalizability of Fourier-MIONet is enabled by its adherence to the physical principle that the solution to a PDE is continuous over time. Furthermore, the developed Fourier-MIONet makes it possible to solve the long-time evolution of geological carbon sequestration in a large-scale three-dimensional space accurately and efficiently.

97 MATHEMATICS AND COMPUTING↗

Pore-Scale Characterization of Biogeochemical Controls on Iron and Uranium Speciation under Flow Conditions

Etched silicon microfluidic pore network models (micromodels) with controlled chemical and redox gradients, mineralogy, and microbiology under continuous flow conditions are used for the incremental development of complex microenvironments that simulate subsurface conditions. We demonstrate the colonization of micromodel pore spaces by an anaerobic Fe(III)-reducing bacterial species (Geobacter sulfurreducens) and the enzymatic reduction of a bioavailable Fe(III) phase within this environment. Using both X-ray Microprobe and X-ray Absorption Spectroscopy, we investigate the combined effects of the precipitated Fe(III) phases and the microbial population on uranium biogeochemistry under flow conditions. Precipitated Fe(III) phases within the micromodel were most effectively reduced in the presence of an electron shuttle (AQDS), and Fe(II) ions adsorbed onto the precipitated mineral surface without inducing any structural change. In the absence of Fe(III), U(VI) was effectively reduced by the microbial population to insoluble U(IV), which was precipitated in discrete regions associated with biomass. In the presence of Fe(III) phases, however, both U(IV) and U(VI) could be detected associated with biomass, suggesting re-oxidation of U(IV) by localized Fe(III) phases. These results demonstrate the importance of the spatial localization of biomass and redox active metals, and illustrate the key effects of pore-scale processes on contaminant fate and reactive transport.

Environmental Molecular Sciences Laboratory↗

Phosphate-Induced Immobilization of Uranium in Hanford Sediments

Phosphate can be added to subsurface environments to immobilize U(VI) contamination. The efficacy of immobilization depends on the site-specific groundwater chemistry and aquifer sediment properties. Batch and column experiments were performed with sediments from the Hanford 300 Area in Washington State and artificial groundwater prepared to emulate the conditions at the site. Batch experiments revealed enhanced U(VI) sorption with increasing phosphate addition. X-ray absorption spectroscopy measurements of samples from the batch experiments found that U(VI) was predominantly adsorbed at conditions relevant to the column experiments and most field sites (low U(VI) loadings, <25 μM), and U(VI) phosphate precipitation occurred only at high initial U(VI) (>25μM) and phosphate loadings. While batch experiments showed the transition of U(VI) uptake from adsorption to precipitation, the column study was more directly relevant to the subsurface environment because of the high solid:water ratio in the column and the advective flow of water. In column experiments, nearly six times more U(VI) was retained in sediments when phosphate-containing groundwater was introduced to U(VI)-loaded sediments than when the groundwater did not contain phosphate. This enhanced retention persisted for at least one month after cessation of phosphate addition to the influent fluid. Sequential extractions and laser-induced fluorescence spectroscopy of sediments from the columns suggested that the retained U(VI) was primarily in adsorbed forms. These results indicate that in situ remediation of groundwater by phosphate addition provides lasting benefit beyond the treatment period via enhanced U(VI) adsorption to sediments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Monitoring Preferential Flow of Water in Sand Using Thermoacoustics Wave Imaging

Abstract Accurate predictions of fluid flow, mass transport, and reaction rates critically impact the efficiency and reliability of subsurface exploration and sustainable use of subsurface resources. Quantitative dynamical sensing and imaging can play a pivotal role in the ability to make such predictions. Geophysical thermoacoustic technology has the potential to provide the aforementioned capabilities, since it builds upon the principle that electromagnetic and mechanical wave fields can be coupled through a thermodynamic process. In this letter, we present laboratory experiments featuring the efficacy of thermoacoustic imaging in the monitoring of preferential flow of water in porous media. Our laboratory experimental equipment can be readily packaged in a form factor that fits in a borehole, and the use of multiple acoustic transducers—which can be combined with volumetric coding techniques—has the potential to provide quasi‐real‐time imaging (0.5 Hz video rate) of regions in close proximity (a few meters) of an open field well.

Geology↗

Thermo-Hydrological Modeling of Thermal Energy Storage in a Depleted Oil Reservoir

Thermal energy storage in oil and gas reservoirs leverages the existing surface and subsurface infrastructure, which can pave the way for economic production of geothermal energy. Existing studies on geothermal energy storage are focused mostly on the use of aquifers with more homogeneous rock and fluid properties. Coupling of heat and fluid flow in a multiphase-multicomponent system, such as an oil reservoir, is imperative especially if existing oil field assets need to be repurposed as required for a sustainable energy transition. The objective is to model the subsurface thermo-hydrological processes associated with reservoir performance and operational sustainability. The model evaluates formation pressure and temperature within the reservoir and at the injection/production wells during multiple charge and discharge cycles. Hot water (approximately 200 degrees C) heated by Concentrating Solar Power (CSP) at high pressure is injected into the existing oil reservoir for storage and produced as thermal energy for power generation, which will be accompanied by enhanced oil recovery. To demonstrate the coupled fluid and heat flow during the injection/production cycle in the subsurface reservoir, TOUGH3 (developed by Berkeley Lab) is used to simulate the thermo-hydrological (TH) processes in a multiphase, multicomponent system. Two well geometries are considered within the reservoir grid: 1) a single-well huff-n-puff system (same well is used for injection and production), and 2) an isolated injection-production well doublet. Seasonal charge and discharge cycling are implemented based on the scheduling specified in the model input file. The model reports pressure, temperature, enthalpy, liquid fluxes, heat fluxes, pore velocities, and changes in porosity & permeability due to temperature and pressure variations during the cyclic Reservoir Thermal Energy Storage (RTES) operations. The results from the simulations can be used to optimize the operational parameters (such as well spacing and injection/production rates) and round-trip efficiency for surface power-plants coupled with thermal energy storage over time. They can also serve as important inputs for levelized cost of storage estimations. The research will help to design and integrate surface renewable energy sources, such as concentrating solar power (CSP), with RTES to help balance out power supply and demand on the grid.

CSP↗

Wetness index based on landscape position and topography ($\mathrm{WILT}$): Modifying TWI to reflect landscape position

Water and land resource management planning benefits greatly from accurate prediction and understanding of the spatial distribution of wetness. The topographic wetness index (TWI) was conceived to predict relative surface wetness, and thus hydrologic responsiveness, across a watershed based on the assumption that shallow slope-parallel flow is a major driver of the movement and distribution of soil water. The index has been extensively used in modeling of landscape characteristics responsive to wetness, and some studies have shown the TWI performs well in landscapes where interflow is a dominant process. However, groundwater flow dominates the hydrology of low-slope landscapes with high subsurface conductivities, and the TWI assumptions are not likely to perform well in such environments. For groundwater dominated systems, we propose a hybrid wetness index (Wetness Index based on Landscape position and Topography, WILT) that inversely weights the upslope contributing area by the distance to the nearest surface water feature and the depth to groundwater. When explicit depth to groundwater data are not available, height above and separation from surface water features can act as surrogates for proximity to groundwater. Here, the resulting WILT map provides a more realistic spatial distribution of relative wetness across a low-slope Coastal Plain landscape as demonstrated by improved prediction of hydric soils, depth to groundwater, nitrogen and carbon concentrations in the A horizon of the soil profile, and sensitivity to DEM scale.

60 APPLIED LIFE SCIENCES↗

Role of Outer Membrane C-Type Cytochromes MtrC and OmcA in Shewanella Oneidensis MR-1 Cell Production, Accumulation, and Detachment During Respiration on Hematite

Solid phase iron oxides are considered to be important terminal electron acceptors for microbial respiration in many anoxic environments. Besides the knowledge that cells attach to and reduce these substrates, other aspects of surface-associated cell behavior and the related cell surface components that influence cell-mineral interactions are not well understood. In the present study, wild-type cells of the dissimilatory iron-reducing bacterium Shewanella oneidensis MR-1 formed thin biofilms one-to-two cell layers in thickness when respiring on natural specular hematite under flow conditions similar to those which exist in aquatic sediments and subsurface environments. The distribution of cells within the biofilm indicated that direct contact was not required for electron transfer from cells to the mineral surface. Detached biomass in the form of single cells represented >99% of the surface-associated wild-type cell production from respiration on hematite over the biofilm life cycle. A mutant deficient in the outer membrane c35 type cytochrome OmcA, while still able to respire and replicate on hematite, established a lower steady-state cell density on the mineral surface than that of the wild-type strain. A mutant deficient in MtrC, another outer membrane c-type cytochrome, and a mutant deficient in both cytochromes were unable to reduce sufficient amounts of hematite to support detectable growth on the mineral surface. When considered in the context of previous work, the results support a growing body of evidence that the relative importance of OmcA and MtrC to cell respiration and replication depends on the form of iron oxide available as terminal electron acceptor.

Mitchell, Andrew C↗

Composite Analysis Vadose Zone Transport of Selected Radionuclides Released from DOE O 435.1 Sources.

The main objectives of the composite analysis (CA) compliance case vadose zone (VZ) facet are to simulate transport of radionuclide releases from the surface to the water table, and to provide radionuclide transfer rates to groundwater for inclusion into the Plateau-to-River (P2R) saturated zone (SZ) model (CP-57037, Model Package Report: Plateau to River Groundwater Model, Version 8.3). The parallel exascale version of the Subsurface Transport Over Multiple Phases (eSTOMP1) simulator is used to simulate flow and transport for the CA VZ models. The CA compliance case approach to address the scale and distribution of contaminant sources in the VZ facet is to subdivide the Hanford Site Central Plateau into 26 models that contain contaminant sources and liquid discharges likely to commingle during migration through the VZ to the water table.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Composite Analysis Vadose Zone Transport of Selected Hanford Site Composite Analysis: Vadose Zone Transport of Selected Radionuclides Released from DOE O 435.1 Waste Sites

The main objectives of the composite analysis (CA) compliance case vadose zone (VZ) facet are to simulate transport of radionuclide releases from the surface to the water table, and to provide radionuclide transfer rates to groundwater for inclusion into the Plateau-to-River (P2R) saturated zone (SZ) model (CP-57037, Model Package Report: Plateau to River Groundwater Model, Version 8.3). The parallel exascale version of the Subsurface Transport Over Multiple Phases (eSTOMP) simulator is used to simulate flow and transport for the CA VZ models. The CA compliance case approach to address the scale and distribution of contaminant sources in the VZ facet is to subdivide the Hanford Site Central Plateau into 26 models that contain contaminant sources and liquid discharges likely to commingle during migration through the VZ to the water table.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Unveiling the Essential Parameters Driving Mineral Reactions during CO2 Storage in Carbonate Aquifers through Proxy Models

Numerical simulation is a commonly employed technique for studying carbon dioxide (CO2) storage processes in porous media, particularly saline aquifers. It enables the representation of diverse trapping mechanisms and the assessment of CO2 retention capacity within the subsurface. The intricate physicochemical phenomena involved necessitate the incorporation of multiphase flow, accurate depiction of fluid and rock properties, and their interactions. Among these factors, geochemical reaction rates and mechanisms are pivotal for successful CO2 trapping in carbonate reactive rocks. However, research on kinetic parameters and the influence of lithology on CO2 storage remains limited. This limitation is partly due to the challenges faced in laboratory experiments, where the time scale of the reactions and the lack of in situ conditions hinder accurate measurement of mineral reaction rates. This study employs proxy models constructed using response surfaces calibrated with simulation results to address uncertainties associated with geochemical reactions. Monte Carlo simulation is utilized to explore a broader range of parameters and identify influential factors affecting CO2 mineralization. The findings indicate that an open database containing kinetic parameters can support uncertainty assessment. Additionally, the proxy models effectively represent objective functions related to CO2 injectivity and mineralization, with calcite dissolution playing a predominant role. pH, calcite concentration, and CO2 injection rate significantly impact dolomite precipitation, while quartz content remains unaffected.

58 GEOSCIENCES↗

Implementation of Enhanced Attenuation at the DOE Mound Site OU-1 Landfill: Accelerating Progress and Reducing Costs - 16270

A field demonstration of enhanced attenuation (EA) was initiated in 2014 to expedite the remediation of industrial solvents in the groundwater impacted by the former landfill Operable Unit 1 (OU-1) at the US Department of Energy (DOE) Mound, Ohio, Site. EA uses active engineering solutions to alter the target site so that the contaminant plume will passively stabilize and shrink. For OU-1, the EA was implemented by targeted injection of an electron donor (vegetable oil) to create 'structured geochemical zones'. Volatile organic compounds (VOCs) in the altered subsurface system encounter alternating anaerobic and aerobic environments along the plume flow path. The effect of these structured geochemical zones is effective and rapid degradation of parent compounds such as trichloroethene (TCE) and tetrachloroethene (PCE), and daughter compounds such as cis-1,2-dichloroethene and vinyl chloride. EA at Mound OU-1 provides a transition to natural attenuation and is an alternative to the baseline pump-and-treat (P and T) system. Industrial solvents (primarily TCE) and other VOCs that originated from the former solid-waste landfill contaminated the groundwater in the Buried Valley Aquifer beneath the Mound site. A groundwater P and T system was initiated in 1996 to control the plume, and a soil vapor extraction system was installed and operated from 1997 to 2003 to remove the bulk of the VOC sources from the vadose zone. Physical removal of the landfill waste and contaminated soil was performed between 2006 and 2010. Approximately 76,100 cubic meters (99,500 cubic yards) of material were removed from the OU-1 landfill area. Operation of the groundwater P and T system continued throughout this period. The principles and technical basis of the ongoing EA field demonstration were developed by the Interstate Technology and Regulatory Council. The goal of the EA implementation is to provide a bridge between active groundwater treatments and monitored natural attenuation. EA relies on understanding and manipulating the contaminant mass balance. Technologies are applied to sustainably reduce the contaminant loading from the source and/or increase the rate of natural attenuation processes within the plume. In 2014, a core team that included DOE, the US Environmental Protection Agency (EPA), the Ohio EPA, and technical support organizations evaluated remedial options for OU-1. During the 2014 evaluation, the team projected a cleanup time frame for the baseline P and T system of 26 additional years and a cleanup time frame for the EA remedy of 13 years. On the basis of that evaluation, EPA approved the implementation of the EA remedy that required the shutdown of the P and T system. The combination of technologies implemented for OU-1 EA included (a) neat vegetable oil in the vadose zone in the former landfill and (b) emulsified vegetable oil within the footprint of the groundwater plume. In the first part of the deployment, neat oil spreads laterally and forms a thin layer on the water table beneath residual soil sources to intercept and reduce future VOC loading and to reduce oxygen inputs to the local groundwater. In the second part of the deployment, emulsified oil forms active bioremediation anaerobic treatment zones within the plume footprint to degrade existing groundwater contaminants (via reductive dechlorination, co-metabolism, or both) and stimulates long-term attenuation capacity in the down-gradient aerobic plume (via co-metabolism). The results from the first year of the field demonstration indicate that the dissolved TCE and PCE plumes have decreased in size and mass. Statistical tests indicate that the VOC concentrations in the majority of the source-area wells are decreasing. The VOC concentrations throughout the plume are projected to be below the target maximum contaminant level of 2 to 5 μg/L within 5 to 10 years-sooner than the design projection of 13 years. The biogeochemistry and microbial community within the structured geochemical zones are consistent with the design and appear to be stable. At Mound OU-1, EA has accelerated progress toward the remedial action objectives and reduced costs. Ongoing monitoring will determine whether the EA remedy has transitioned the site into monitored natural attenuation status and whether the EA will continue to be effective, timely, and sustainable. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Play-fairway analysis for geothermal exploration: Examples from the Great Basin, western USA

Elevated permeability within fault systems provides pathways for circulation of geothermal fluids. Future geothermal development depends on precise and accurate location of such fluid flow pathways in order to both accurately assess geothermal resource potential and increase drilling success rates. The collocation of geologic characteristics that promote permeability in a given geothermal system define the geothermal ‘fairway’, the location(s) where upflow zones are probable and where exploration efforts including drilling should be focused. We define the geothermal fairway as the collocation of 1) fault zones that are ideally oriented for slip or dilation under ambient stress conditions, 2) areas with a high spatial density of fault intersections, and 3) lithologies capable of supporting dense interconnected fracture networks. Areas in which these characteristics are concomitant with both elevated temperature and fluids are probable upflow zones where economic-scale, sustainable temperatures and flow rates are most likely to occur. Employing a variety of surface and subsurface data sets, we test this ‘play-fairway’ exploration methodology on two Great Basin geothermal systems, the actively producing Brady’s geothermal system and a ‘greenfield’ geothermal prospect at Astor Pass, NV. These analyses, based on 3D structural and stratigraphic framework models, reveal subsurface characteristics about each system, well beyond the scope of standard exploration methods. At Brady’s, the geothermal fairways we define correlate well with successful production wells and pinpoint several drilling targets for maintaining or expanding production in the field. In addition, hot-dry wells within the Brady’s geothermal field lie outside our defined geothermal fairways. At Astor Pass, our play-fairway analysis provides for a data-based conceptual model of fluid flow within the geothermal system and indicates several targets for exploration drilling.

15 GEOTHERMAL ENERGY↗

Subsurface Characterization and Machine Learning Predictions at Brady Hot Springs Results

Geothermal power plants typically show decreasing heat and power production rates over time. Mitigation strategies include optimizing the management of existing wells - increasing or decreasing the fluid flow rates across the wells - and drilling new wells at appropriate locations. The latter is expensive, time-consuming, and subject to many engineering constraints, but the former is a viable mechanism for periodic adjustment of the available fluid allocations. Data and supporting literature from a study describing a new approach combining reservoir modeling and machine learning to produce models that enable strategies for the mitigation of decreased heat and power production rates over time for geothermal power plants. The computational approach used enables translation of sets of potential flow rates for the active wells into reservoir-wide estimates of produced energy and discovery of optimal flow allocations among the studied sets. In our computational experiments, we utilize collections of simulations for a specific reservoir (which capture subsurface characterization and realize history matching) along with machine learning models that predict temperature and pressure timeseries for production wells. We evaluate this approach using an "open-source" reservoir we have constructed that captures many of the characteristics of Brady Hot Springs, a commercially operational geothermal field in Nevada, USA. Selected results from a reservoir model of Brady Hot Springs itself are presented to show successful application to an existing system. In both cases, energy predictions prove to be highly accurate: all observed prediction errors do not exceed 3.68% for temperatures and 4.75% for pressures. In a cumulative energy estimation, we observe prediction errors that are less than 4.04%. A typical reservoir simulation for Brady Hot Springs completes in approximately 4 hours, whereas our machine learning models yield accurate 20-year predictions for temperatures, pressures, and produced energy in 0.9 seconds. This paper aims to demonstrate how the models and techniques from our study can be applied to achieve rapid exploration of controlled parameters and optimization of other geothermal reservoirs. Includes a synthetic, yet realistic, model of a geothermal reservoir, referred to as open-source reservoir (OSR). OSR is a 10-well (4 injection wells and 6 production wells) system that resembles Brady Hot Springs (a commercially operational geothermal field in Nevada, USA) at a high level but has a number of sufficiently modified characteristics (which renders any possible similarity between specific characteristics like temperatures and pressures as purely random). We study OSR through CMG simulations with a wide range of flow allocation scenarios. Includes a dataset with 101 simulated scenarios that cover the period of time between 2020 and 2040 and a link to the published paper about this project, where we focus on the Machine Learning work for predicting OSR's energy production based on the simulation data, as well as a link to the GitHub repository where we have published the code we have developed (please refer to the repository's readme file to see instructions on how to run the code). Additional links are included to associated work led by the USGS to identify geologic factors associated with well productivity in geothermal fields. Below are the high-level steps for applying the same modeling + ML process to other geothermal reservoirs: 1. Develop a geologic model of the geothermal field. The location of faults, upflow zones, aquifers, etc. need to be accounted for as accurately as possible 2. The geologic model needs to be converted to a reservoir model that can be used in a reservoir simulator, such as, for instance, CMG STARS, TETRAD, or FALCON 3. Using native state modeling, the initial temperature and pressure distributions are evaluated, and they become the initial conditions for dynamic reservoir simulations 4....

15 GEOTHERMAL ENERGY↗

Collapse of Reacted Fracture Surface Decreases Permeability and Frictional Strength

Abstract Geochemical and geomechanical perturbations of the subsurface caused by the injection of fluids present risks of leakage and seismicity. This study investigated how acidic fluid flow affects hydraulic and frictional properties of fractures using experiments with 3.8‐cm‐long specimens of Eagle Ford shale, a laminated shale with carbonate‐rich strata. In low‐pressure flow cells, one set of samples was exposed to acidic brine and another set was exposed to neutral brine. X‐ray computed tomography and energy‐dispersive X‐ray spectroscopy revealed that samples exposed to acidic brine were calcite‐depleted and had developed a porous altered layer, while the other set showed no evidence of alteration. After reaction, samples were compressed and sheared in a triaxial cell that supplied normal stress and differential pore pressure at prescribed sliding velocities, independently measuring friction and permeability. During the initial compression, the porous altered layer collapsed into fine particles that filled the fracture. This effectively impeded flow and sealed the fracture, resulting in fracture permeability to decrease 1 to 2 orders of magnitude relative to the unaltered fractures. This is a favorable outcome in subsurface applications where the goal is to reduce leakage risks. However, during shear the reacted fracture had lower frictional strength because the fine‐grained particles in the collapsed layer prevented the formation of interlocking microasperities. Therefore, coupled geochemical and geomechanical processes that could favorably seal fractures could also increase the likelihood of induced seismicity. These findings have important implications for geological carbon sequestration, pressurized fluid energy storage, geothermal energy, and other subsurface technologies.

58 GEOSCIENCES↗

Zero-Valent Iron Permeable Reactive Barrier to Remediate Volatile Organic Compounds in Groundwater - 20209

The US DOE-SRS, the US EPA, and the South Carolina Department of Health and Environmental Control determined it was appropriate to perform a non-time critical removal action at the P-Area Groundwater Operable Unit at the SRS to reduce the mass and downgradient transport of trichloroethylene in the P-Area groundwater plume. Contaminated groundwater discharges to a nearby stream, Steel Creek, within the SRS boundaries, resulting in trichloroethylene concentrations above the maximum contaminant level. Impact to surface water is limited in areal extent and supported by recently collected characterization data. The P-Area Groundwater Operable Unit encompasses the groundwater beneath an industrial area within SRS, P Area, where the P-Reactor once operated. The boundaries of the P-Area Groundwater Operable Unit extend northwest to Steel Creek, northeast toward PAR Pond, and southeast to Meyers Branch. Groundwater in the Upper Three Runs Aquifer of the P-Area Groundwater Operable Unit has been impacted by reactor and facility operations between 1954 and 1991, including tritium and volatile organic compounds. The P-Area surface units contributing to groundwater contamination were remediated as part of the P-Area Operable Unit in 2011. The P-Reactor closure is one of the first of its kind in the DOE Complex and is one of only a few full-sized production reactors in the US to undergo completion of final closure activities. The nature and extent of groundwater contamination was determined using a variety of investigative approaches such as groundwater monitoring wells, direct-push technology, and surface water samples. Groundwater contamination associated with trichloroethylene is primarily exhibited in a narrow plume that extends from the source area at P-Reactor and west to Steel Creek. Maximum contaminant level exceedances in groundwater occur over an area of ∼6.9 hectares for trichloroethylene with concentrations as high as 7.7 milligrams per liter. To the west of the P-Area facility area, the trichloroethylene groundwater plume is controlled by a buried geologic feature, assumed to be an old stream bed, that further narrows the groundwater plume in what has been designated as the 'neck area.' This narrowing of the groundwater plume provides an ideal location for a treatment barrier. The non-time critical removal action alternative chosen is to install a zero-valent iron permeable reactive barrier within the neck area of the trichloroethylene groundwater plume, perpendicular to groundwater flow direction. This technology will provide a treatment barrier that will reduce trichloroethylene groundwater concentrations by 90% and has an anticipated useful life of at least 25 years. A pre-design investigation was performed in the neck area to confirm site lithology, hydrogeology, geochemistry, and extent of trichloroethylene contamination prior to a final design. A treatability study, conducted as part of the pre-design investigation, indicated that the subsurface and groundwater in the PArea Groundwater Operable Unit is compatible with the zero-valent iron and will not lead to excessive buildup from mineralization/precipitation or biofouling. Probabilistic modeling was conducted using field and laboratory data to determine the expected performance of the zero-valent iron permeable reactive barrier. The model simulations indicated that a 3.81-centimeter thick barrier would provide greater than 90% reduction of trichloroethylene groundwater concentration. The final design of the zero-valent iron permeable reactive barrier is a barrier that will extend 80.5 linear meters in a 'zigzag' orientation to best transect the trichloroethylene plume and account for varying groundwater flow. The barrier will be installed from 13.7 meters below ground surface to 41.1 meters below ground surface for 65.8 linear meters and from 13.7 meters below ground surface to 36.6 meters below ground surface for 14.6 linear meters, the base of which is 'keyed' into a low permeability zone. The barrier is designed to a thickness of 10.2 centimeters, which was determined to reduce the trichloroethylene groundwater concentrations by greater than 90% with a safety factor of 2.67. A total of approximately 689 metric tons of zero-valent iron will be injected through 22 injection wells spaced 3.66 meters apart, using guar to suspend the zero-valent iron. Zero-valent iron permeable reactive barrier construction will be monitored through 23 installed resistivity receivers offset 7.32 meters from the zerovalent iron permeable reactive barrier. The zero-valent iron will be energized with a low-voltage 100 Hertz signal during injection and will be monitored using the resistivity receivers to ensure complete coalescence of the zero-valent iron permeable reactive barrier. The zero-valent iron was sized to have a hydraulic conductivity greater than the natural subsurface, thus promoting groundwater flow through the barrier. As contaminated groundwater contacts the zero-valent iron, volatile organic compounds, including trichloroethylene, are immediately degraded to harmless compounds such as ethylene. The performance of the zero-valent iron permeable reactive barrier will be monitored using three upgradient monitoring well clusters, six downgradient monitoring well clusters, and four in-wall monitoring wells. The in-wall monitoring wells will indicate immediate reduction of trichloroethylene mass in the groundwater and allow analyses of the zero-valent iron permeable reactive barrier health. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗