Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Subsurface geologic modeling”

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 343 records · Page 19

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↗

Characterization of Most Promising Sequestration Formations in the Rocky Mountain Region

The project Characterization of Most Promising Sequestration Formations in the Rocky Mountain Region is one of 9 site characterization projects that were implemented as part of ARRA (American Recovery and Reinvestment Act). Data from this project was used to improve resolution of data in NATCARB in the area of study. Data related to this study has already been incorporated in NATCARB Atlas. The Rocky Mountain Carbon Capture and Storage (RMCCS) project investigated multiple geologic formations and characterized a local site on the Colorado Plateau for future CCS opportunities. The RMCCS project focused on the Cretaceous Dakota, Jurassic Entrada, and Pennsylvanian Weber Sandstones, the three largest regional formations. All formations in this project are potential CO2 storage resources for future power plants, natural gas processing plants, cement plants, and oil shale development projects. The area adjacent to Craig, Colorado, (Sand Wash Basin) was the area selected for detailed geologic characterization on the RMCCS project. The basin was selected in part because the geology can be extrapolated to other sites on the Colorado Plateau. Field mapping and seismic surveys were conducted to identify and evaluate the basin's structural configuration. A 9,745-foot deep characterization well was drilled to collect 131 feet of core and a suite of geophysical well log data. Petrophysical tests on samples of core were used to calibrate geophysical log data, which can be used to obtain storage resource estimates and evaluate associated uncertainty as well as simulate the hydrologic behavior of injected CO2. A detailed analysis of the primary formations (Dakota, Entrada and Weber sandstones) yielded a more accurate CO2 storage resource assessment for these formations within the Colorado Plateau; RMCCS estimates indicate a total CO2 storage resource of more than 38,000 million metric tons. The characterization of the Sand Wash Basin (2-D seismic surveys, multiple well logs and lithological, petrophysical and geochemical analyses) allowed for a detailed 3-D model to be constructed. The model served as the framework for analyses ranging from CO2 storage resource, injectivity, and subsurface flow to uncertainty estimates to evaluation of risk.

2-D seismic↗

Understanding and Predictability of Integrated Mountain Hydroclimate (Workshop Report)

Mountainous systems cover approximately 23% of Earth’s land and are distributed across all continents. They can capture and store atmospheric moisture that is then cycled through the terrestrial surface and subsurface system, released to downstream communities, and cycled back to the atmosphere. Mountain hydroclimate—characterized by steep gradients, geological, ecological, and biogeochemical diversity—is influenced by topographic forcing and elevated warming and susceptible to large subseasonal to multidecadal variability and rapid changes. Terrestrial hydrological and biogeochemical cycles also experience cascading effects from global warming impacts, such as multidecadal declines in mountain snowpack, longer growing seasons, and increased frequency and severity of extreme events like droughts and wildfires. However, little is known about the effects of these impacts and their feedbacks on climate systems and surface-subsurface compartments. Also unknown are the full implications of changing hydroclimate and extreme events on hydro biogeochemical cycles across atmosphere, terrestrial, and human systems in mountain regions and beyond. This knowledge gap is critical, given human reliance on mountain systems for stable water supply and quality. Mountain systems’ increasing vulnerability to climate change and human perturbations motivates the need to improve understanding of integrated mountain hydroclimate (IMHC) systems and their feedbacks and impacts on humans across scales. However, due to large heterogeneity and strong gradients, coupled natural-human processes in mountain regions present significant challenges for observations, modeling, predictions, and projections. Motivated by gaps in mountain hydroclimate understanding, observations, and modeling and the need for credible projections of future changes, the U.S. Department of Energy’s (DOE) Biological and Environmental Research (BER) program organized a virtual workshop on “Understanding and Predictability of Integrated Mountain Hydroclimate.” Sponsored by BER’s Earth and Environmental Systems Sciences Division (EESSD), the workshop aimed to inform and catalyze EESSD’s growing interests in enhancing predictive understanding of IMHC. Organizers structured the workshop to identify (1) knowledge gaps, (2) observational and modeling challenges, (3) short-term (1 to 3 years) and long-term (10 years and beyond) research opportunities, and (4) strategies for fostering collaboration and coordination. To address the outstanding challenges of IMHC, the workshop included two sessions organized by disciplinary, cross-disciplinary, and crosscutting science topics. The disciplinary and cross-disciplinary topics focused on essential IMHC elements: atmosphere, terrestrial, and human systems and their interactions. Breakout sessions on disciplinary and cross-disciplinary topics facilitated identification of crosscutting topics and central emerging themes. Session 1 focused on connecting existing DOE investments to accelerate progress related to scientific challenges in understanding mountain hydroclimate. In Session 2, participants further explored key Session 1 takeaways through the lens of multiagency collaborations and coordination.

54 ENVIRONMENTAL SCIENCES↗

Determining the 3D Subsurface Density Structure of Taurus Littrow Valley Using Apollo 17 Gravity Data

Surface gravity surveys can detect subsurface density variations that can reveal subsurface geologic features. In 1972, the Apollo 17 (A17) mission conducted the Traverse Gravimeter Experiment (TGE) using a gravimeter that measured the local gravity field near Taurus Littrow Valley (TLV), located on the south-eastern rim of the Serenitatis basin. TLV is hypothesized to be a basaltfilled radial graben resulting from the impact that formed Mare Serenitatis. It is bounded by both the North and South Massifs (NM and SM) as well as other smaller mountains to the East that are thought to be mainly composed of brecciated highland material. The TGE is the first and only successful gravity survey on the surface of the Moon. Other more recent satellite surveys, such as NASA's Gravity Recovery and Interior Laboratory (GRAIL) mission (2011- 2012), have produced the best global gravity field to date (approx. 13km resolution). However, these satellite surveys are not sensitive enough to detect fine-scale (<1km) lunar subsurface structures. This underscores the value of the data collected at the surface by A17. In the original analysis of the data a 2D forward-modelling approach was used to derive a thickness of the subsurface basalt layer of 1.0 km by assuming a simple flat-faced rectangular geometry and using densities derived from Apollo lunar samples. We are investigating whether modern 3D modelling techniques in combination with high-resolution topographical and image datasets can reveal additional fine-scale subsurface structure in TLV.

TGE↗

Coupling subsurface and above-surface models for optimizing the design of borefields and district heating and cooling systems in the presence of varying water-table depth. In: Proceedings, 46th Workshop on Geothermal Reservoir Engineering

Dynamic energy simulation is important for the design and sizing of district heating and cooling systems with geothermal heat exchange. Current modeling approaches in building and district energy simulation tools typically consider heat conduction through the ground between boreholes, without flow of groundwater. On the other hand, detailed simulation tools for subsurface heat and mass transfer exist, but these fall short in simulating above-surface energy systems. To support the design and operation of such systems, we have developed a coupled model including a software package for building and district energy simulation, and software for detailed heat and mass transfer in the ground. For the first, we use 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, we use the TOUGH simulator. TOUGH 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. In previous work, we described the coupling of these software packages, including how time-dependent boundary conditions for the borehole walls are synchronized for use in Modelica and TOUGH. We verified that the coupled Modelica/TOUGH code produced consistent results with the original Modelica code for an idealized problem in which heat transfer was purely by conduction in a uniform geologic medium. Here, we examine less idealized problems for which TOUGH’s advanced capabilities for modeling fluid flow are required. The first problem has a shallow vadose zone, and the second problem has a thicker vadose zone with a water-table depth that varies in time, which requires a fine vertical grid discretization for the TOUGH model.

Doughtry, Christine↗

Integration of Large-Scale Electrical Imaging into Geological Framework Development and Refinement

Geologic framework models (GFMs) are critical to the construction of reliable simulation models of groundwater flow and contaminant transport. To support GFM development, direct information (e.g., core samples, fluid samples, hydraulic testing) tends to be sparse and separated by large distances relative to the spatial scales of aquifer heterogeneity. There are additional challenges associated with highly contaminated legacy waste sites, where drilling is particularly costly, and invasive sampling requires specialized handling and disposal of hazardous materials. At these sites in particular, non-invasive geophysical imaging can play an important role in filling spatial gaps between boreholes and reducing characterization costs by optimizing and minimizing the number of necessary boreholes. Here this paper presents a case study demonstrating the use of large-scale (> 30 km 2 ) electrical mapping to identify hydrostratigraphy and potential paleochannels at the Hanford Site, located in Washington State, USA. In two field campaigns, over 36 line-kilometers of electrical resistivity tomography (ERT) data were collected along 14 transects. ERT surveys were sited and performed to image critical aspects (e.g., paleochannels, stratigraphic contacts) of the subsurface, demonstrating a general workflow for integrating ERT with GFM development. Inconsistencies between the GFM and ERT were catalogued to provide a basis for future site characterization using complementary geophysical methods and (or) direct sampling.

58 GEOSCIENCES↗

Influence of geochemical reactions on the creep behavior of Mt. Simon sandstone

The capture and subsurface storage of carbon dioxide is a sustainable option that is currently pursued worldwide to mitigate greenhouse gas effect. However, predicting the long-term mechanical integrity of CO 2 underground storage systems remains a challenge. To address that question, it is essential to understand the influence of fluid-rock chemo-mechanical interactions on the long-term and on the time-dependent mechanical properties. In turn, the long-term mechanical response and the time-dependent mechanical behavior can be represented by the creep response. We investigate the impact of CO 2 -induced geochemical reaction on the creep response of Mt. Simon sandstone with a 50–400 micron grain size. We perform static and dynamic flow experiments on Mt. Simon sandstone specimens under geological conditions, at a temperature of T = 50 - 53 °C and for CO 2 pressures of P = 8:62; 17:2 MPa under both static flow and dynamic flow-through conditions. After aging, we employ creep indentation testing, high-resolution SEM-EDS, computer vision, machine learning, and micromechanics modeling to probe changes on the microstructure and mechanical properties. Following both static and dynamic flow-through experiments, we observe a 10–22% decrease in quartz volume fraction and an increase in both the microporosity (7–28%) and nanoporosity (60–65%). Additional CO 2 -induced microstructural changes include an enlargement of pore throats and the formation of channels. These observations point to the presence of CO 2 -induced K-feldspar dissolution and clay dissolution reactions. Here, the macroscopic creep behavior is logarithmic and the macroscopic creep modulus varies depending on the microporosity and the relative quartz and feldspar content. As a result of these geochemical reactions and of the related microstructural changes, a 55-60% decrease in the macroscopic logarithmic creep modulus is predicted.

58 GEOSCIENCES↗

Deep Learning for Modeling Enhanced Geothermal Systems

Enhanced Geothermal Systems (EGS) offer a vast potential to expand the use of geothermal energy. Heat is extracted from this engineered system by injecting cold water into a subsurface fractures, which are in contact with the hot dry rock, and pulled through the production wells. Creating EGS requires improving the natural permeability of hot crystalline rocks. To develop economically viable EGS reservoirs, significant technical barriers (e.g., better stimulation technologies without adequate water and/or permeability) and non-technical barriers (e.g., land access and permitting) must be overcome. In this short conference paper, we present a workflow to address a part of this challenge – “How to develop economically viable EGS using existing technologies?”. Our workflow called the GeoThermalCloud (GTC) for EGS, leverages recent advances in machine learning, deep learning, and cloud computing. This GTC framework is open-source and available at https://github.com/SmartTensors/GeoThermalCloud.jl. The GTC framework provides trained deep learning (DL) models to estimate the net present value of a given EGS design scenario. The Geothermal Design Tool (https://github.com/GeoDesignTool/GeoDT.git), a fast and simplified multi-physics solver, is used to develop a database for training DL models. The database consists of EGS design parameters (inputs to DL model) and their net present value (output of DL model) in uncertain geologic systems. The EGS design parameters for constructing this training database are based on Utah FORGE but include the options of more wells and deeper depths. The DL models are trained by ingesting the EGS design parameters and estimating the corresponding net present value. Such an emulation allows us to screen various EGS designs quickly and identify good development strategies by coupling them with optimization techniques. Our preliminary results show promise in DL emulation of net present value. However, a lot more work is needed to improve the predictive capability of DL models (i.e., extensive hyperparameter tuning is necessary). This will be the primary focus of our future work.

artificial neural networks, geothermal↗

Mitigating Injection‐Induced Seismicity Along Basement Faults by Extraction: Application to 2016–2018 Pohang Earthquakes

Abstract The risk of induced earthquake is considered as the major obstacle for sustainable subsurface energy activities. As a mitigation approach, mechanistic risk assessment can be essential to minimize the potential of seismic hazard. Coupled hydro‐mechanical modeling approach is integrated with geomechanical analysis to better assess the occurrence of seismicity induced by injection and/or extraction, which will provide a basis to update operational conditions (e.g., number and location of mitigating wells, rate and duration of operations) for mitigating earthquake nucleation. The Coulomb stability analysis and seismicity rate estimates from this generic study indicate that simultaneous extraction can relieve injection‐induced pore‐pressure buildup, but also transfer and accumulate elastic energy along the faults. Field‐scale simulation with variation in well design for the Pohang EGS site, South Korea, show that multiple wells aligned across a preexisting fault form gradients of pressure and stress fields normal to the fault, which can enhance the direct impact of diffusion process and poroelastic stressing on the fault stability perturbation. Our mechanistic study suggests multiphysics‐oriented mitigation strategies that consider (1) possible mechanical deformation, caused by hydro‐mechanical coupling effects, and (2) proper configuration of well operations based on site‐specific geological characteristics.

Chang, K. W.↗

Final Report DE-FE0031785 Mohsen Ahmadian, Ph.D. Demonstration of Proof of Concept of a Multiphysics Approach for Real-Time Remote Monitoring of Dynamic Changes in Pressure and Salinity in Hydraulically Fractured Networks

Hydraulic fracturing has evolved into a multistep process with varying flow rates, carrier fluids (e.g., gel or slickwater), proppant loadings, and proppant grain sizes. As a result, primary recovery from a hydraulically fractured tight-oil reservoir is often a tiny fraction of the original oil in place, ranging between 5 and 10%. As stated in the FOA1990, “part of this problem is due to the inability of current well completion processes to effectively stimulate the entire reservoir volume in contact with the wellbore. Innovative technologies are needed that can help improve the effectiveness of reservoir completion methods, maximize stimulated reservoir volumes, and optimize recovery over the entire producing life span of a well”. We first need to enhance the current fracture diagnostic techniques to improve a well-completion design. However, detecting and delineating a subsurface hydraulic fracture is extremely difficult because the induced fracture network is only fractionally propped, and these propped fractures are generally very thin. Microseismic and tiltmeter monitoring techniques can provide information on the fracture extent but provide little or no information on the movement and final distribution of proppant or production fluids. On the other hand, electromagnetic (EM) imaging has shown the capability to monitor proppant distribution throughout the fracture area, especially in the presence of Electrically Active Proppants (EAPs). A previous EM survey of hydraulic fracturing at the Devine Fracture Pilot Site (DFPS) and subsequent EM code developments demonstrated this survey as a robust technique to remotely interrogate the extent of the EAP-filled hydraulic fracture during its propagation. The objectives of the project were threefold: (1) to capitalize on the material properties of an EAP to demonstrate remote monitoring of relative changes in pressure, pressure, and flow that are commonly encountered during production from a hydraulically fractured reservoir; (2) to evaluate EM imaging tools, to achieve Objective 1 in near real-time; and (3) to develop a multi-physics joint inversion approach to precisely predict flow patterns and physiochemical changes within an EAP-filled fracture network. This research project was built upon our previous work at the Devine Test Site managed by the Bureau of Economic Geology (BEG) at The University of Texas at Austin (UT-Austin). It also leveraged a significant investment from the Advanced Energy Consortium (AEC) to address the DOE's interest in subsurface flow, containment, and characterization by multiple signals. This three-year and three-month project succeeded in demonstrating the feasibility of a real-time dynamic fluid flow mapping technique at Technology Readiness Level 5 (TRL5) by utilizing a commercially available surface-based Controlled-Source Electromagnetic (CSEM) method (Objectives 1, 2). We demonstrated that injections into an EAP-filled fracture could be successfully coupled with real-time electric field measurements on the surface, leading to remote monitoring of dynamic changes within the EAP-filled fracture. Furthermore, the observed electric field in our study is influenced by bottomhole pressure, flow rate, and salinity, which is demonstrated by comparing these parameters with the electrical field potentials. EM simulations solely based on assumptions of fracture conductivity changes during injection did not reproduce the whole measured electric field magnitudes. Preliminary estimates showed that including Streaming Potential (SP) in our geophysical model is likely needed to reduce the simulation misfit.

02 PETROLEUM↗

Numerical simulation of flow and mixing in fracture intersections

Fluid transport through fractured geological formations is strongly influenced by the redistribution of solutes at fracture intersections. In this study, we perform detailed numerical simulations of flow and scalar transport within the intersection of two smooth, planar fractures. The analysis focuses on the mixing ratio, the proportion of solute flux exiting along the outlet branch aligned with the primary inlet flow direction, relative to the total solute flux at the outlets. We systematically investigate how the mixing ratio varies with four key parameters: Peclet number, Reynolds number, flow rate ratio between outlet branches, and fracture intersection angle. Results show that the mixing ratio decreases with increasing Peclet number and outlet flow rate ratio, consistent with reduced diffusive spreading and enhanced streamline routing. While low Reynolds numbers have minimal impact, inertial effects at higher Reynolds numbers significantly increase the mixing ratio. Additionally, acute and obtuse intersection angles alter flow partitioning and modify the solute distribution at the outlets. These findings provide a quantitative basis for incorporating physically realistic mixing behavior—intermediate between complete mixing and streamline-following assumptions—into network-scale transport models. The results have direct relevance to subsurface energy systems, including geothermal energy production, carbon sequestration, and contaminant remediation.

58 GEOSCIENCES↗

Skillful Multi‐Month Predictions of Ecosystem Stressors in the Surface and Subsurface Ocean

Abstract Anthropogenic carbon emissions and associated climate change are driving rapid warming, acidification, and deoxygenation in the ocean, which increasingly stress marine ecosystems. On top of long‐term trends, short term variability of marine stressors can have major implications for marine ecosystems and their management. As such, there is a growing need for predictions of marine ecosystem stressors on monthly, seasonal, and multi‐month timescales. Previous studies have demonstrated the ability to make reliable predictions of the surface ocean physical and biogeochemical state months to years in advance, but few studies have investigated forecast skill of multiple stressors simultaneously or assessed the forecast skill below the surface. Here, we use the Community Earth System Model (CESM) Seasonal to Multiyear Large Ensemble (SMYLE) along with novel observation‐based biogeochemical and physical products to quantify the predictive skill of dissolved inorganic carbon (DIC), dissolved oxygen, and temperature in the surface and subsurface ocean. CESM SMYLE demonstrates high physical and biogeochemical predictive skill multiple months in advance in key oceanic regions and frequently outperforms persistence forecasts. We find up to 10 months of skillful forecasts, with particularly high skill in the Northeast Pacific (Gulf of Alaska and California Current Large Marine Ecosystems) for temperature, surface DIC, and subsurface oxygen. Our findings suggest that dynamical marine ecosystem prediction could support actionable advice for decision making.

Environmental Sciences & Ecology↗

Geostatistical Mapping of Salinity Conditioned on Borehole Logs, Montebello Oil Field, California

We present a geostatistics-based stochastic salinity estimation framework for the Montebello Oil Field that capitalizes on available total dissolved solids (TDS) data from groundwater samples as well as electrical resistivity (ER) data from borehole logging. Data from TDS samples (n = 4924) was coded into an indicator framework based on falling below four selected thresholds (500, 1000, 3000, and 10,000 mg/L). Collocated TDS-ER data from the surrounding groundwater basin were then employed to produce a kernel density estimator to establish conditional probabilities for ER data (n = 8 boreholes) falling below the selected TDS thresholds within the Montebello Oil Field area. Directional variograms were estimated from these indicator coded data, and 500 TDS realizations from conditional indicator simulation were generated for the subsurface region above the Montebello Oil Field reservoir. Simulations were summarized as 3D maps of median TDS, most likely salinity class, and probability for exceeding each of the specified TDS thresholds. Results suggested TDS was below 500 mg/L in most of the study area, with a trend toward higher values (500 to 1000 mg/L) to the southwest; consistent with the average regional groundwater flow direction. Discrete localized zones of TDS greater than 1000 mg/L were observed, with one of these zones in the greater than 10,000 mg/L range; however, these areas were not prevalent. The probabilistic approach used here is adaptable and is readily modified to include additional data and types and can be employed in time-lapse salinity modeling through Bayesian updating.

54 ENVIRONMENTAL SCIENCES↗

The effect of original and initial saturation on residual nonwetting phase capillary trapping efficiency

Injection of supercritical carbon dioxide (CO 2 ) into geological formations is a strategy for both atmospheric greenhouse gas reduction (climate change mitigation) and enhanced oil recovery. To understand CO 2 trapping efficiency, the capillary trapping behaviors that immobilize subsurface fluids were analyzed at the pore-scale using pairs of proxy fluids representing the range of in situ nonwetting and wetting fluid properties encountered in geologic storage reservoirs. The pairs of fluids were cycled through imbibition and drainage processes using a flow cell apparatus containing a sintered glass bead column. Computed x-ray microtomography (microCT) was used to identify immobilized nonwetting fluid after imbibition and drainage events. From microCT images, the trapped residual (post-secondary imbibition) nonwetting phase was spatially correlated to both the original (post-primary imbibition) and the initial (post-primary drainage) nonwetting phase; this relationship is referred to here as the original saturation dependence (SO-dependence) and initial saturation dependence (SI-dependence), respectively. Significant trends of decreasing SO- and SI-dependence with increasing wetting and nonwetting fluid phase viscosities were observed. This finding implies that the amount of CO 2 injected and ultimately trapped is dependent on the nonwetting phase (e.g. oil or gas) already present in the formation, as well as on the manner in which supercritical CO 2 is initially injected, which are factors not considered in current trapping models. To explore the potential effect that varying viscosity and interfacial tension (IFT) might have on this process, we report on a variety of fluid pairs with different viscosities and IFT.

42 ENGINEERING↗

Thermopile Energy Harvesting for Subsurface Wellbore Sensors (Final Report)

Robust in situ power harvesting underlies all efforts to enable downhole autonomous sensors for real-time and long-term monitoring of CO 2 plume movement and permeance, wellbore health, and induced seismicity. This project evaluated the potential use of downhole thermopile arrays, known as thermoelectric generators (TEGs), as power sources to charge sensors for in situ real-time, long-term data capture and transmission. Real-time downhole monitoring will enable “Big Data” techniques and machine learning, using massive amounts of continuous data from embedded sensors, to quantify short- and long-term stability and safety of enhanced oil recovery and/or commercial-scale geologic CO 2 storage. This project evaluated possible placement of the TEGs at two different wellbore locations: on the outside of the casing; or on the production tubing. TEGs convert heat flux to electrical power, and in the borehole environment, would convert heat flux into or out of the borehole into power for downhole sensors. Such heat flux would be driven by pumping of cold or hot fluids into the borehole—for instance, injecting supercritical CO 2 —creating a thermal pulse that could power the downhole sensors. Hence, wireless power generation could be accomplished with in situ TEG energy harvesting. This final report summarizes the project’s efforts that accomplished the creation of a fully operational thermopile field unit, including selection of materials, laboratory benchtop experiments and thermal-hydrologic modeling for design and optimization of the field-scale power generation test unit. Finally, the report describes the field unit that has been built and presents results of performance and survivability testing. The performance and survivability testing evaluated the following: 1) downhole power generation in response to a thermal gradient produced by pumping a heated fluid down a borehole and through the field unit; and 2) component survivability and operation at elevated temperature and pressure conditions representative of field conditions. The performance and survivability testing show that TEG arrays are viable for generating ample energy to power downhole sensors, although it is important to note that developing or connecting to sensors was beyond the scope of this project. This project’s accomplishments thus traversed from a low Technical Readiness Level (TRL) on fundamental concepts of the application and modeling to TRL-5 via testing of the fully integrated field unit for power generation in relevant environments. A fully issued United States Patent covers the wellbore power harvesting technology and applications developed by this project.

47 OTHER INSTRUMENTATION↗

Simulator for Hydrologic Unstructured Domains (SHUD v1.0): numerical modeling of watershed hydrology with the finite volume method

Abstract. Hydrologic modeling is an essential strategy for understanding and predicting natural flows, particularly where observations are lacking in either space or time or where complex terrain leads to a disconnect in the characteristic time and space scales of overland and groundwater flow. However, significant difficulties remain for the development of efficient and extensible modeling systems that operate robustly across complex regions. This paper introduces the Simulator for Hydrologic Unstructured Domains (SHUD), an integrated, multiprocess, multiscale, flexible-time-step model, in which hydrologic processes are fully coupled using the finite volume method. SHUD integrates overland flow, snow accumulation/melt, evapotranspiration, subsurface flow, groundwater flow, and river routing, thus allowing physical processes in general watersheds to be realistically captured. SHUD incorporates one-dimensional unsaturated flow, two-dimensional groundwater flow, and a fully connected river channel network with hillslopes supporting overland flow and baseflow. The paper introduces the design of SHUD, from the conceptual and mathematical description of hydrologic processes in a watershed to the model's computational structures. To demonstrate and validate the model performance, we employ three hydrologic experiments: the V-catchment experiment, Vauclin's experiment, and a model study of the Cache Creek Watershed in northern California. Ongoing applications of the SHUD model include hydrologic analyses of hillslope to regional scales (1 m2 to 106 km2), water resource and stormwater management, and interdisciplinary research for questions in limnology, agriculture, geochemistry, geomorphology, water quality, ecology, climate and land-use change. The strength of SHUD is its flexibility as a scientific and resource evaluation tool where modeling and simulation are required.

58 GEOSCIENCES↗

Conceptual Model Testing Related to SDU 6 Drainwell Observations

From its inception in the early 1950s through the end of the Cold War in the early 1990s, the Savannah River Site (SRS) produced nuclear materials for national defense in five reactors. Additionally, irradiated reactor fuel and target tubes were dissolved in nitric acid to recover plutonium and uranium using the PUREX (Plutonium Uranium Reduction EXtraction) process. Liquid waste from these chemical separations processes was then stored onsite in 51 underground tanks. Eight waste storage tanks have been operationally closed (i.e. cleaned and grouted) and the remaining tanks hold a mixture of liquids, insoluble solids, and precipitated salts (SRMC-LWP-2022-00001), the latter generated by evaporating water from the liquid waste. Waste is currently being retrieved from tanks and separated into 1) high-radioactivity, low-volume, and 2) low-radioactivity, high-volume components, principally through the Salt Waste Processing Facility (SWPF) (SRMC-LWP-2023-00001). The former waste stream is vitrified in the Defense Waste Processing Facility (DWPF), stored onsite, and destined for offsite disposal in a deep geologic repository. The latter stream is mixed with dry cementitious materials in the Saltstone Production Facility (SPF) and the wet slurry placed in onsite Saltstone Disposal Units (SDUs) within the Saltstone Disposal Facility (SDF), where it hardens into a cement waste form termed saltstone. A low-infiltration surface cover system will be placed over the SDF at closure, where SDUs will then be in the subsurface post-closure (SRR-CWDA-2019-00001).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Prediction of Distributed River Sediment Respiration Rates Using Community-Generated Data and Machine Learning

River sediment microbial respiration is a key indicator of ecosystem functioning and the biogeochemical fluxes across this critical zone link surface and subsurface waters. As such, there is tremendous interest in measuring and mapping these respiration rates. Respiration observations are expensive and labor intensive; there is limited data available to the community. An open science, collaborative initiative is collecting samples for respiration rate analysis and multi-scale metadata; this evolving data set is being used for making machine learning (ML) predictions at unsampled sites to help inform continued community engagement. However, it is a challenge to find an optimum configuration for ML models to work with this feature-rich (i.e., 100+ possible input variables) data set. Here, we present results from a two-tiered approach to managing the analysis of this complex data set: (a) a stacked ensemble of models that automatically optimizes hyperparameters and manages the training of many models and (b) feature permutation importance to detect the most important features in the models. The major elements of this workflow are modular, portable, open, and cloud-based thus making this implementation a potential template for other applications. The models developed here predict that sediment organic matter chemistry is one of the most important features for predicting sediment respiration rate. Other larger-scale, important features fall into the categories of climatic, ecological, geological, and fluvial settings. Leveraging these larger-scale features to generate data-driven estimates of river sediment respiration rates reveals spatially consistent but heterogeneous patterns across the river network of the Columbia River Basin.

54 ENVIRONMENTAL SCIENCES↗