Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “3D Groundwater Model”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Tracking Radioactive Isotopes in HVAC and Application for Hot Cell Analyses - 20155

In recent years, the surge in the number of isotope production facilities under design has increased the need to analyze isotope migration in facility Heating, Ventilation, and Air Conditioning (HVAC) during normal operations and accident scenarios. These facilities are used to produce isotopes for medical, security, and industrial applications and as such are subject to license and regulatory requirements 10CFR20, 10CFR30, 10CFR50, and 10CFR70. GOTHIC, a general-purpose thermal-hydraulics software package, includes the ability to model isotopic tracers, radioactive decay and isotope migration as well as HEPA and charcoal filters for isotope retention. A GOTHIC model was developed by Zachry Nuclear Engineering (ZNE) to examine the effect of negative room pressure, HEPA filtration, and HVAC fluctuations on radiation areas, hot cells and gloveboxes. Radioactive tracers were used to simulate the concentration of spills within contaminated areas, track the migration of isotopes of interest, and determine the isotopic retention and buildup on facility HEPA filters. A variety of isotopes with concern to dose (e.g., Kr-85, Sr-90, I-131, etc.), including their decay and progeny, are included in the analysis. Negative pressures are maintained in the regions of interest by a representative central HVAC system equipped with a volumetric fan that exhausts to the environment after a series of isolation valves and HEPA filters. GOTHIC is an industry trusted tool for providing engineering solutions for a variety of applications, including fission product tracking, aerosol and particulate transport and ventilation assessments. The software provides an integrated analysis environment that includes a graphical user interface (GUI) for constructing analysis models, a numerical solver that includes parallel processing capabilities and a post-processor for evaluating simulation results. It solves the conservation equations for mass, momentum and energy for multicomponent, multi-phase flow in lumped parameter and multi-dimensional geometries (1, 2, or full 3D), including the effects of turbulence, diffusion and buoyancy. It has been developed and maintained under a Quality Assurance program in compliance with the requirements of 10CFR50 Appendix B and applicable portions of ASME NQA-1 since 1995. GOTHIC has been used for assessing both forced and natural convection conditions for a wide range of applications, including: - Tracking concentration of hazardous gases and chemicals for habitability and safety assessments - Determining ventilation and filtration requirements and optimizing location and arrangement of these systems - Room heat-up, including diverse and Flexible coping strategies for Extended Loss of AC Power (FLEX/ELAP) - Equipment Qualification (EQ) A distinctive feature of GOTHIC is the ability to track many different fields/substances in a simulation, including user defined tracer elements, in the liquid, vapor and droplet fields as well as surfaces and filters. This capability allows GOTHIC to model fission product transport and release or the removal of particulates or harmful toxins from exhaust gases using a spray scrubber or other types of filtration systems. GOTHIC also includes models for engineered equipment, such as fans, filters, charcoal filters, dryers/demisters, dampers, etc. The aerosols and other filtered material are removed or accumulated in these components. The range of aerosol and radiological applications that GOTHIC has been used for includes: - Source Term: Primary Coolant (Equilibrium) Activity; Non-Water Coolant Source/Leakage. - Conditions for Iodine Re-evolution: Sump/Suppression Pool Conditions and pH; - RWST Conditions and pH. - Isotope Removal Mechanisms: Containment Sprayed and Unsprayed Region Mixing; Charcoal Filter Heating due to Iodine decay. - Radionuclide Transport and Decay: Post-LOCA Release in containment; Transport between connected Compartments and vent systems; Groundwater transport of radionuclides. - Non-Newtonian Fluid modeling for sludge, waste tanks, etc. ADAMS ML071581053 (titled 'Best Practice Guidelines for the use of CFD in Nuclear Reactor Safety Applications') poses guidelines for applying single phase CFD codes in nuclear reactor safety problems and GOTHIC is listed as a 'tool for 3D flows' and 'dispersal and deposition of radionuclides.' The Nuclear Quality Assurance (NQA) pedigree of GOTHIC is an important aspect for applications in the nuclear industry. The fundamental tracer models (convective transport, molecular and turbulent diffusion, removal mechanisms, etc.) have been verified using analytical solutions and validated against applicable separate effects tests. Also, GOTHIC has been benchmarked to many integrated effects tests, including Phebus FP (Fission Product). GOTHIC gives good agreement for the buildup and decay of fission products in Phebus Test 3. The model developed by ZNE demonstrates GOTHIC's applicability and acceptability for use in analyzing the migration and retention of radioactive isotopes and their progeny in normal operation and accident scenario analyses for isotope production facilities, hot cells, and gloveboxes. The tracer activities calculated by GOTHIC can then be used in downstream radiation transport and shielding codes like RADTRAD-NAI{sup C}, MCNP{sup R}, and MicroShield{sup R} to determine on-site and of-site doses. (authors)

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

Development of inter-grid-cell lateral unsaturated and saturated flow model in the E3SM Land Model (v2.0)

Abstract. The lateral transport of water in the subsurface is important in modulating terrestrial water energy distribution. Although a few land surface models have recently included lateral saturated flow within and across grid cells, it is not a default configuration in the Climate Model Intercomparison Project version 6 experiments. In this work, we developed the lateral subsurface flow model within both unsaturated and saturated zones in the Energy Exascale Earth System Model (E3SM) Land Model version 2 (ELMv2.0). The new model, called ELMlat, was benchmarked against PFLOTRAN, a 3D subsurface flow and transport model, for three idealized hillslopes that included a convergent hillslope, divergent hillslope, and tilted V-shaped hillslope with variably saturated initial conditions. ELMlat showed comparable performance against PFLOTRAN in terms of capturing the dynamics of soil moisture and groundwater table for the three benchmark hillslope problems. Specifically, the mean absolute errors (MAEs) of the soil moisture in the top 10 layers between ELMlat and PFLOTRAN were within 1 %±3 %, and the MAEs of water table depth were within ±0.2 m. Next, ELMlat was applied to the Little Washita experimental watershed to assess its prediction of groundwater table, soil moisture, and soil temperature. The spatial pattern of simulated groundwater table depth agreed well with the global groundwater table benchmark dataset generated from a global model calibrated with long-term observations. The effects of lateral groundwater flow on the energy flux partitioning were more prominent in lowland areas with shallower groundwater tables, where the difference in simulated annual surface soil temperature could reach 0.3–0.4 ∘C between ELMv2.0 and ELMlat. Incorporating lateral subsurface flow in ELM improves the representation of the subsurface hydrology, which will provide a good basis for future large-scale applications.

58 GEOSCIENCES↗
DOE OSTI2022

Hanford Site Composite Analysis: Comparison of One- and Three-Dimensional Simulation of Vadose Zone Transport and Transfer to Groundwater for Selected Waste Sites

The purpose of this document is to compare 1D and 3D technetium-99 (Tc-99) VZ transport and transfer rates to groundwater for five selected waste sites (two tanks, two cribs, and one trench) to demonstrate the need to consider 3D flow and transport mechanisms and commingling of contaminant plumes originating from adjacent waste sites. The waste sites selected for the comparison evaluation are the 241-A-105 Tank in the A Farms Area model, the 216-B-18 Crib and the 216-B-25 Trench in the BC Cribs and Trenches model, the 216-B-49 Crib in the B Complex model, and the 241-SX-115 Tank in the S Farms Area model. The selected waste sites all represent different disposal or leak conditions, with and without potential impacts due to release-fluid commingling, occurring at the Hanford Site Central Plateau.

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

The role of the spatial heterogeneity and correlation length of surface wettability on two-phase flow in a CO 2 -water-rock system

This study characterized and modeled heterogeneous surface wettability in sandstone and investigated the role of spatial heterogeneity and correlation length of surface wettability on relative permeability in a supercritical CO 2 (scCO 2 )-brine-rock system. Understanding the role of wettability heterogeneity on relative permeability is essential to geological CO 2 sequestration, oil and gas recovery, and contaminated groundwater remediation. Although numerous studies have attempted to understand the influences of surface wettability, capillary number (Ca), and viscosity ratio, the role of the spatial variation and correlation length of surface wettability on two-phase flow in three-dimensional (3D) porous media has not been unraveled due to the challenges in the measurement and representation of realistic rock surface wettability. In this work, we conducted in-situ measurements of surface contact angle (CA) in a Bentheimer sandstone after CO 2 flooding using micro-computed tomography (micro-CT), and found that the pore-scale CA distribution on rock surfaces followed a log-normal distribution associated with a spatial correlation length. Based on the statistical information from CT scanning, a Gaussian random field was used to model CA distributions that had desired standard deviations and spatial correlation lengths, which were then adjusted within a certain range of values for sensitivity analyses to study their combined effects on the two-phase flow in the porous medium using the lattice Boltzmann (LB) method. The LB two-phase flow simulation was accelerated using hybrid, multicore parallel computing to overcome the challenges in simulating multiphase flow in a large 3D domain having 800 × 800 × 600 nodes. The simulation results showed that the surface wettability heterogeneity (i.e., standard deviation of CA) had a lesser effect on the relative permeability of the wetting fluid (water) but a more significant impact on the relative permeability of the non-wetting fluid (scCO 2 ). The Corey model was used to fit the LB-simulated relative permeability curves of water and scCO 2 and showed that the variations in the relative permeability curves for both water and scCO 2 increased as the standard deviation and spatial correlation length of CA increased. This study illustrated that the assumption of homogeneous surface wettability may cause errors in multiphase flow simulations. Furthermore, the impacts of both the standard deviation and spatial correlation length of CAs should be accounted for. This is the first study that explored the spatial correlation lengths associated with CA distributions on sandstone surfaces and comprehensively investigated the roles of both spatial variation and correlation length of CA on two-phase flow properties in 3D porous media. The optimized LB multiphase flow model was proved a powerful tool to study the interplays and combined effects of these statistical parameters, which had critical applications in numerous natural and engineering processes that involved multiphase flow in porous media.

58 GEOSCIENCES↗
DOE OSTI2026

PFLOTRAN modeling data and scripts associated with “Refining the Hydrogeologic Framework of a Large River Corridor Model Using Waterborne Transient Electromagnetics”

NOTE: The manuscript associated with this data package is currently in review. The data may be revised based on reviewer feedback. Upon manuscript acceptance, this data package will be updated with the final dataset and additional metadata. This data package is associated with the publication “Refining the Hydrogeologic Framework of a Large River Corridor Model Using Waterborne Transient Electromagnetics” submitted to Water Resources Research (Terry et al. 2025). The data package contains the groundwater modeling dataset from PFLOTRAN software. It includes the python script for mesh generation, boundary condition setting, PFLOTRAN input deck formation and postprocessing. It couples groundwater flow and species transport for Hanford Reach river corridor and pipelines the model generation and processing. This model can be used to easily generate the model and analysis for Hanford site. It can also be adjusted to other hydrologic area with ease. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. The data package consists of 6 folders: (1) “data” contains all necessary data as input and intermediate data for processing; (2) “mesh” contains all mesh related files to generate mesh in Hanford Reach river corridor; (3) “model_run” contains the generated script for PFLOTRAN modeling; (4) “notebooks” contains all the Python script to generate the model; (5) “output” contains all the output from the computation; (6) “postprocessing” contains the Python script to generate scientific figure for manuscript. All files are .csv (comma-separated values), .h5 (HDF5 format), .in (input files), .ipynb (Jupyter notebooks), .p (Python pickle), .png (images), .PNG (images), .py (Python scripts), .pyc (Python bytecode), .r (R scripts), .sh (shell scripts), .txt (text files), .vtu (3D mesh/visualization format), .xz (compressed archive), or .zip (compressed archive).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Incorporating Heterogeneity into 3D Geologic Models for the Alluvial Basin Reference Case

The unsaturated zone alluvium reference case is one of several geologic systems under consideration by the U.S. Department of Energy Office of Nuclear Energy for hosting repositories for spent nuclear fuel and associated waste (Sevougian et al., 2019). As noted by Mariner et al. (2018), the generic alluvial basin offers positive attributes that merit its consideration as a reference case by the Spent Fuel and Waste Science and Technology (SFWST) campaign. Hundreds of alluvial basins and sub-basins are scattered across the arid western United States (Figure 1-1). Precipitation and infiltration rates are relatively low with high evapotranspiration, resulting in deep regional water tables. Deep regional water tables allow for significant vertical separation between proposed repository depths and groundwater and thus longer transport paths from a hypothetical repository to an aquifer. Accumulations of alluvial sediments within these basins are typically on the order of hundreds of meters, and locally may exceed 1,000 m in thickness, as is the case for the Deming sub-basin in southern New Mexico. In alluvial basins, alluvium creates a thick geologic host medium that serves as the natural barrier system (NBS) in the conceptual model framework of a geologic disposal system and isolates the waste packages from receptors in the biosphere. Further, alluvial basin fill is typically comprised of stacked playa and lacustrine deposits along the basin axis (Perry et al., 2018). Characterized by low permeability, these layers protect the biosphere above the repository and the groundwater resources below the repository by limiting flow and transport of radionuclides.

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

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

Dynamic risk assessment for geologic CO 2 sequestration

At a geologic CO 2 sequestration (GCS) site, geologic uncertainty usually leads to large uncertainty in the predictions of properties that influence metrics for leakage risk assessment, such as CO 2 saturations and pressures in potentially leaky wellbores, CO 2 /brine leakage rates, and leakage consequences such as changes in drinking water quality in groundwater aquifers. The large uncertainty in these risk-related system properties and risk metrics can lead to over-conservative risk management decisions to ensure safe operations of GCS sites. The objective of this work is to develop a novel approach based on dynamic risk assessment to effectively reduce the uncertainty in the predicted risk-related system properties and risk metrics. We demonstrate our framework for dynamic risk assessment on two case studies: a 3D synthetic example and a synthetic field example based on the Rock Springs Uplift (RSU) storage site in Wyoming, USA. Results show that the U.S. National Risk Assessment Partnership’s Open Source Integrated Assessment Model (NRAP-Open-IAM) coupled with a conformance evaluation can be used to effectively quantify and reduce the uncertainty in the predictions of risk-related system properties and risk metrics in GCS.

58 GEOSCIENCES↗
DOE OSTI2021

The effects of well damage and completion designs on geoelectrical responses in mature wellbore environments

Well integrity is one of the major concerns in long-term geologic storage sites due to the potential risk of well leakage and groundwater contamination. Evaluating changes in electrical responses due to energized steel-cased wells has the potential to quantify and predict possible wellbore failures because any kind of breakage or corrosion along highly conductive well casings will have an impact on the distribution of the subsurface electrical potential. However, realistic wellbore-geoelectrical models that can fully capture fine-scale details of well completion design and the state of well damage at the field scale require extensive computational effort, or they can even be intractable to simulate. To overcome this computational burden while still keeping the model realistic, we have used the hierarchical finite-element method that represents electrical conductivity at each dimensional component (1D edges, 2D planes, and 3D cells) of a tetrahedral mesh. This allows well completion designs with real-life geometric scales and well systems with realistic, detailed, progressive corrosion and damage in our models. We have developed a comparison of possible discretization approaches of a multicasing completion design in the finite-element model. The effects of the surface casing and the coupling between concentric well casings as well as the effects of the degree and the location of well damage on the electrical responses are also examined. As a result, we analyze real surface electric field data to detect wellbore integrity failure associated with damage.

58 GEOSCIENCES↗
DOE OSTI2023

Building an EPA Class VI Permit Application

Summary To accelerate the commercialization of carbon capture and storage (CCS), the US Department of Energy (US DOE) is building on decades of characterization efforts and pilot-scale projects through their CarbonSAFE program. Administered through their National Energy Technology Laboratory, this program seeks to bring fully integrated projects to the sector that can store more than 50 million tonnes of CO2 over a 30-year period. The program, which was enacted before the enhancement of Internal Revenue Code Section 45Q, is in the capture assessment, characterization, and permitting phase. The objectives of this paper are to discuss (a) the injection permitting requirements of the CarbonSAFE projects; (b) information gathering in support of the permit; (c) the timelines of field development and permit-related activities; (d) the major technical components of the field development plan; and (e) early feedback from the regulators toward acceptance of the permit. In Mississippi, more than 30,000 acres have been characterized by six deep characterization wells, a deep groundwater well, and 92 line miles of 2D seismic as part of the CarbonSAFE Project ECO2S. During the acquisition of seismic data, all receiver lines were live, which resulted in the generation of a pseudo-3D seismic design. The incorporation of a 3D seismic survey was not included as part of this project due to logistical difficulties presented by the undulating, wooded surface terrain. A suite of openhole geophysical logs was taken from each well, allowing for a detailed interpretation of prospective storage reservoirs and confining intervals to complement the analysis carried out on the 290 ft of a whole core that was cut through the prospective confining zone and storage reservoir. The detailed geologic and reservoir data were assembled and entered into a 3D model to assess the injection capacity and the area of review (AoR). This information fed into the detailed corrective action, monitoring, testing, and postinjection site care (PISC) modeling. The results have been exceptional. The geologic assessment has revealed three primary storage targets, ranging in depth from 3,500 ft to 6,000 ft. These storage reservoirs net 1,300 ft of sandstone, with mean porosity and permeability of 29% and 3.6 darcies, respectively. Together, these reservoirs have storage capacities that may exceed 20 million tonnes per square mile, making this a gigatonne prospect. Forward modeling of the project resulted in an AoR of 16 sq miles, injecting about 8000 t/d, for 30 years, via two deep injection wells. The excellent confining characteristics of the caprock, relatively simple geologic structure, and lack of historical well drilling activity in this area provide excellent containment of the injected CO2. Based on this work, the project has proposed 20 years of PISC. To date, only two US CO2 injection permits have been granted. These projects relied on a singular capture point feeding a singular sequestration point (source to sink), and considerations have not been made to garner CO2 emissions from other industrial sources. The Kemper County Storage Complex is a first-of-its-kind storage hub concept that looks to develop an area capable of storing significant quantities of CO2 from the region. Also, this work will show how characterization efforts, geological and numerical modeling efforts, and plan development were constructed in support of permit and incentives acceptance.

Energy & Fuels↗
DOE OSTI2021

Stratigraphic Identification with Airborne Electromagnetic Methods at the Hanford Site, Washington

Stratigraphic units can influence the fate and transport of subsurface contaminants within groundwater. Units having coarse-grained sediments act as preferential flow pathways, and therefore can accelerate the transport of contaminants to reach human and ecological receptors. At legacy waste sites, detailed knowledge of subsurface stratigraphy can be used for effective monitoring and remediation planning to help minimize risk to human health and the environment. Airborne electromagnetic (AEM) methods can non-invasively provide information on kilometer-scale or larger subsurface stratigraphic features and fill informational gaps in directly sampled data from sparsely located boreholes. In this paper, we present inversion results of a 412 line-km frequency-domain AEM survey to delineate subsurface stratigraphic features at the Hanford Site, located in southeastern Washington State. The inversion was performed using a massively parallel 3D electromagnetic modeling and inversion code, where the modeling is based on solving frequency-domain Maxwell’s equations using an unstructured-mesh finite-element method and the inversion employs a Gauss-Newton optimization scheme. The results are compared to an underlying geologic framework model (GFM), built by interpolating contact depths of stratigraphic units interpreted from site borehole datasets. In areas with good borehole coverage, the inversion results show a good match with the GFM to a depth of about 60 m. Outside of these areas, the inversion results exhibit inconsistencies from the assumptions made to create the GFM, demonstrating that the AEM survey results can be used to improve the understanding of the geological conceptual model.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Effects of Dissolution Regimes on Flow Channelization and Solute Transport in 3D Fracture Networks: Insights From Graph‐Based Reactive Transport Modeling

We investigate how mineral dissolution reshapes flow pathways and solute transport in three‐dimensional discrete fracture networks using a computationally efficient graph‐based reactive transport model. The DFNs are inspired by field‐site observations of fractured carbonate and represent realistic connectivity and structural heterogeneity. Flow is simulated with the Reynolds equation, and dissolution follows first‐order kinetics with diffusive limitations captured through an effective mass‐transfer coefficient. By systematically varying two key dimensionless parameters, the effective Damköhler number (Da), governing reaction versus advection rates, and a transport parameter (Da), analogous to the Thiele modulus, distinct flow channelization regimes emerge: mildly channelized at low G, highly channelized at intermediate Da, and extreme wormhole formation at high Da and low G. Eulerian and Lagrangian analyses, including breakthrough curves, particle tortuosity, dispersivity, and flow channeling indicators quantitatively characterize the progression of dissolution‐driven network restructuring. Across all regimes, initial fracture heterogeneity persists. The results underscore how the interplay between this initial structure, advection, reaction, and diffusion critically shapes subsurface flow pathways, with implications for applications ranging from groundwater remediation to enhanced geothermal systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Subsurface Characterization for Evaluating Geothermal Resource Potential from Existing Oil and Gas Wells in Tuttle, Oklahoma: Preprint

Oil and gas (O&G) wells often encounter co-produced hot water, possibly suitable for geothermal direct-use applications. The City of Tuttle is located on the eastern part of the Anadarko sedimentary basin in Oklahoma with high heat-in-place potential and recovery capability at depth. This study aims at demonstrating the potential of geothermal energy production for direct-use applications in two public schools and 250 nearby houses in Tuttle via repurposing existing O&G wells. In this scope, geochemistry, geology, and borehole log data were collected and incorporated into a 3D conceptual subsurface model. A digital elevation model (DEM) was used to represent the study area topography with four O&G wells. In addition, hydrogeochemical characteristics of the geothermal fluid and scaling potential were analyzed using ternary diagrams and chemical ratios to develop mixing models. The subsurface geology model indicated that the study area primarily consists of Permian to Mississippian Sandstone and Limestone formations, implying a porosity ranging between 12% and 22%, and a permeability up to 3.90E-14 m2 in certain reservoir levels. The reservoir temperature is expected to be ranging between 80 degrees C to 95 degrees C around 3 km depth with an average temperature gradient of 22.8 degrees C/km. Chemical geothermometers also estimated the reservoir temperature as 90 degrees C. Findings of the chemical model demonstrated that the geothermal fluid is Sodium-Potassium-Chloride-Sulfate type and possibly mixed with shallow groundwater resulting in higher Ca and Mg concentrations and lower Na/K ratio implying lower calcite scaling. These results comprehensively characterize the potential of geothermal resources in the study area and imply that geothermal energy production by repurposing existing O&G wells is suitable for low-temperature direct-use applications.

GEOTHERMAL ENERGY↗
DOE OSTI2023

Subsurface Characterization for Evaluating Geothermal Resource Potential from Existing Oil and Gas Wells in Tuttle, Oklahoma

Oil and gas (O&G) wells often encounter co-produced hot water, possibly suitable for geothermal direct-use applications. The City of Tuttle is located on the eastern part of the Anadarko sedimentary basin in Oklahoma with high heat-in-place potential and recovery capability at depth. This study aims at demonstrating the potential of geothermal energy production for direct-use applications in two public schools and 250 nearby houses in Tuttle via repurposing existing O&G wells. In this scope, geochemistry, geology, and borehole log data were collected and incorporated into a 3D conceptual subsurface model. A digital elevation model (DEM) was used to represent the study area topography with four O&G wells. In addition, hydrogeochemical characteristics of the geothermal fluid and scaling potential were analyzed using ternary diagrams and chemical ratios to develop mixing models. The subsurface geology model indicated that the study area primarily consists of Permian to Mississippian Sandstone and Limestone formations, implying a porosity ranging between 12% and 22%, and a permeability up to 3.90E-14 m2 in certain reservoir levels. The reservoir temperature is expected to be ranging between 80 degrees C to 95 degrees C around 3 km depth with an average temperature gradient of 22.8 degrees C/km. Chemical geothermometers also estimated the reservoir temperature as 90 degrees C. Findings of the chemical model demonstrated that the geothermal fluid is Sodium-Potassium-Chloride-Sulfate type and possibly mixed with shallow groundwater resulting in higher Ca and Mg concentrations and lower Na/K ratio implying lower calcite scaling. These results comprehensively characterize the potential of geothermal resources in the study area and imply that geothermal energy production by repurposing existing O&G wells is suitable for low-temperature direct-use applications.

gas wells↗
DOE OSTI2023

Advances in 3D Geologic Modeling of Alluvial Basins with a Focus on Facies and Property Modeling

The unsaturated zone alluvium reference case is one of several geologic systems under consideration by the U.S. Department of Energy Office of Nuclear Energy for hosting repositories for spent nuclear fuel and associated waste (Sevougian et al., 2019). As noted by Mariner et al. (2018), the generic alluvial basin offers positive attributes that merit its consideration as a reference case by the Spent Fuel and Waste Science and Technology (SFWST) campaign. There are hundreds of alluvial basins and sub-basins scattered across the arid western United States (Figure 1-1). Precipitation and infiltration rates are relatively low with high evapotranspiration, resulting in vertical separation between repository and water table and thus longer transport paths to an aquifer. Accumulations of alluvial sediments within these basins are typically on the order of hundreds of meters, and locally may exceed 1,000 m in thickness, as is the case for the Deming sub-basin in southern New Mexico. A thick geologic host medium, which serves as the natural barrier system (NBS) in the conceptual model framework of a geologic disposal system, isolates the waste packages from receptors in the biosphere. Further, alluvial basin fill is typically comprised of stacked playa and lacustrine deposits along the basin axis (Perry et al., 2018). Characterized by low permeability, these layers protect the biosphere above the repository and the groundwater resources below the repository.

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

River Dynamics Control Transit Time Distributions and Biogeochemical Reactions in a Dam-Regulated River Corridor

Residence Time Distributions (RTDs) exert an important control on biogeochemical translation in the hydrological systems. Previous tracer studies have revealed that RTDs often followed time-invariant exponential, lognormal, or heavy-tailed RTDs that have power-law behavior for long tails in headwater or low-order streams. However, there is a recognition that RTDs can be more complicated and time-variable in response to dynamic hydrological forcing. Here, we use particle tracking to estimate RTDs along the Hanford Reach of the Columbia River and to quantify the influences of river stage fluctuations. Particle tracking is conducted using the velocity field from high-resolution 3D groundwater flow simulations. The effects of hydrological forcing on the residence time distribution were evaluated by varying river flow boundary conditions and releasing particles in different time windows. Our results revealed that dynamic stage fluctuations created rapidly changing losing-gaining conditions in the river and led to highly transient RTDs, which contributes to multiple modes of RTDs. Dam-induced high-frequency (sub-daily) flow variation contributes to the short-time (sub-daily) component of the RTDs. Deviation of the reactant consumption under the single-mode assumption compared to the multimodal RTDs is relatively small (~5%) and appears when the Damköhler number was close to one.More specifically, high-frequency (daily) dam-induced stage variations have higher impacts on biogeochemical reactions with faster reaction rate, and can potentially increase the reactant consumption by 26.68%. Based on these findings, we suggest that current river basin models could be improved by including bank storage and more complex RTDs influenced by both short and long term river stage fluctuations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

3D Seismic Anatomy of a Watershed Reveals Climate–Topography Coupling That Drives Water Flowpaths and Bedrock Weathering

To investigate how bedrock transforms to soil, we mapped the topography of the interface demarcating onset of weathering under an east-west trending shale watershed in the Valley and Ridge province in the U.S.A. Using wave-equation travel-time tomography from a seismic array of >4000 geophones, we obtained a 3D P-wave velocity (Vp) model that resolves structures ~20 m below land surface (mbls). The depth of mobile soil and the onset of dissolution of chlorite roughly match Vp = 600 m/s and Vp = 2700 m/s, respectively. Chlorite dissolution initiates porosity growth in the shale matrix. Depth to the 2700 m/s contour is greater under the N- as compared to S-facing hillslopes and under sub-planar as compared to concave-up land surfaces. Broadly, the geometries of the ‘soil’ and ‘chlorite’ Vp contours are consistent with the calculated potential for shear fracture opening under weak regional compression. However, this calculated fracture potential does not consistently explain observations related to N- vs. S-facing aspect nor fracture density observed by borehole televiewer. Apparently, regional compression is only a secondary influence on Vp: the primary driver of P-wave slowing in the upper layers of this catchment is topographic control of reactive water flowpaths and their integrated effects on weathering. Here, the Vp result is best explained as the long-term integrated effect of groundwater flow-induced geochemical weathering of shale in response to climate-driven patterns of micro- and macro-topography.

54 ENVIRONMENTAL SCIENCES↗