Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “subsurface system”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 379 records · Page 21

4D electrical resistivity tomography for assessing the influence of vegetation and subsurface moisture on railway cutting condition

Instability of slopes, embankments, and cuttings on the railway network is increasingly prevalent globally. Monitoring vulnerable infrastructure aids in geotechnical asset management, and improvements to transport safety and efficiency. Here, we examine the use of a novel, near-real-time Electrical Resistivity Tomography (ERT) monitoring system for assessing the stability of a railway cutting in Leicestershire, United Kingdom. In 2015, an ERT monitoring system was installed across a relict landslide (grassed) and an area of more stable ground on either side (wooded), to monitor changes in electrical resistivity through time and space, and to assess the influence of different types of vegetation on the stability of transportation infrastructure. Two years of 4-Dimensional ERT monitoring results are presented here, and petrophysical relationships developed in the laboratory are applied to calibrate the resistivity models in order to provide an insight into hydrogeological pathways within a railway cutting. The influence of vegetation type on subsurface moisture pathways and on slope stability is also assessed – here we find that seasonal subsurface changes in moisture content and soil suction are exacerbated by the presence of trees (wooded area). This results in shrink-swell behaviour of the clays comprising the railway cutting, resulting in fissuring and a reduction in shear strength, leading to instability. As such, it is proposed that on slopes comprised of expansive soils, grassed slopes are beneficial for stability. Insights into the use of 4-D ERT for monitoring railway infrastructure gained from this study may be applied to the monitoring of critical geotechnical assets elsewhere.

58 GEOSCIENCES↗

Gas diffusion through variably-water-saturated zeolitic tuff: Implications for transport following a subsurface nuclear event

Noble gas transport through geologic media has important applications in the characterization of underground nuclear explosions (UNEs). Without accurate transport models, it is nearly impossible to distinguish between xenon signatures originating from civilian nuclear facilities and UNEs. Understanding xenon transport time through the earth is a key parameter for interpreting measured xenon isotopic ratios. One of the most challenging aspects of modeling gas transport time is accounting for the effect of variable water saturation of geological media. In this study, we utilize bench-scale laboratory experiments to characterize the diffusion of krypton, xenon, and sulfur hexafluoride (SF6) through intact zeolitic tuff under different saturations. Here, we demonstrate that the water in rock cores with low partial saturation dramatically affects xenon transport time compared to that of krypton and SF6 by blocking sites in zeolitic tuff that preferentially adsorb xenon. This leads to breakthrough trends that are strongly influenced by the degree of the rock saturation. Xenon is especially susceptible to this phenomenon, a finding that is crucial to incorporate in subsurface gas transport models used for nuclear event identification. We also find that the breakthrough of SF6 diverges significantly from that of noble gases within our system. When developing field scale models, it is important to understand how the behavior of xenon deviates from chemical tracers used in the field, such as SF6 (Carrigan et al., 1996). These new insights demonstrate the critical need to consider the interplay between rock saturation and fission product sorption during transport modeling, and the importance of evaluating specific interactions between geomedia and gases of interest, which may differ from geomedia interactions with chemical tracers.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Development of Land‐River Two‐Way Hydrologic Coupling for Floodplain Inundation in the Energy Exascale Earth System Model

Abstract Floodplain inundation links river and land systems through significant water, sediment, and nutrient exchanges. However, these two‐way interactions between land and river are currently missing in most Earth System Models. In this study, we introduced the two‐way hydrological coupling between the land component, E3SM Land Model, and the river component, Model for Scale Adaptive River Transport, in Energy Exascale Earth System Model (E3SM) to study the impacts of floodplain inundation on land and river processes. We calibrated the river channel geometry and developed a new data‐driven inundation scheme to improve the simulation of inundation dynamics in E3SM. The new inundation scheme captures 96% of the spatial variation of inundation area in a satellite inundation product at global scale, in contrast with 7% when the default inundation scheme of E3SM was used. Global simulations including the new inundation scheme performed at resolution with and without land‐river two‐way coupling was used to quantify the impact of coupling. Comparisons show that two‐way coupling modifies the water and energy cycle in 20% of the global land cells. Specifically, riverine inundation is reduced by two‐way coupling, but inland inundation is intensified. Wetter periods are more impacted by the two‐way coupling at the global scale, while regions with different climates exhibit different sensitivities. The two‐way exchange of water between the land and river components of E3SM provides the foundation for enabling two‐way coupling of land‐river sediment and biogeochemical fluxes. These capabilities will be used to improve understanding of the interactions between water and biogeochemical cycles and their response to human perturbations.

54 ENVIRONMENTAL SCIENCES↗

Joint Optimization of Well Completions and Controls for CO 2 Enhanced Oil Recovery and Storage

CO 2 storage through CO 2 enhanced oil recovery (EOR) is considered as one of the technologies to help promote larger scale deployment of CO 2 storage because of associated economic benefits through oil recovery, 45Q tax credits and the utilization of existing infrastructure. The objective of this study is to demonstrate how optimal reservoir management and operation strategies (including well completions and controls) can be used to optimize both CO 2 storage and oil recovery. The optimization problem was focused on jointly estimating the well completions (i.e., fraction of injection/production well perforations in each reservoir layer) and CO 2 injection/oil production controls that maximize the net present value (NPV) in a CO 2 EOR and storage operation. We utilized the newly developed StoSAG algorithm, one of the most efficient optimization algorithms in the reservoir management community, to solve the optimization problem. The performance of joint optimization approach was compared with the performance of well control only optimization approach. In addition, the performance of co-optimization of CO 2 storage and oil recovery approach was compared with the performances of maximization of only CO 2 storage and maximization of only oil recovery approaches. The optimization results showed that a joint optimization of well completions and well controls can achieve an 8.84% higher final NPV than the one obtained from the optimization of only well controls. It was observed that the NPV incremental for joint optimization is mainly due to the fact that the optimal well completions and controls approach results in efficient CO 2 storage and oil production from different reservoir layers depending on the differences in individual layer properties. Comparison of co-optimization (i.e., maximization of NPV) and maximization of only CO 2 storage or only oil recovery showed that the co-optimization and maximization of only oil recovery result in significantly higher final NPV than that obtained through maximization of only CO 2 storage approach while maximization of only CO 2 storage can achieve significantly higher CO 2 storage in the reservoir compared to the other two scenarios. The similar results for co-optimization and maximization of oil production are obtained because of the difference in oil revenue compared to CO 2 storage tax credit. To the best of our knowledge, this is the first study in oil/gas industry and CO 2 storage community to perform joint optimization of well completions and well controls in the fields. We expect that the proposed optimization framework will be a useful and efficient tool for field engineers to optimally manage CO 2 EOR projects to maximize revenue through oil recovery as well as CO 2 storage by taking advantage of the new 45Q tax law.

Artificial intelligence↗

Mixed heavy metal stress induces global iron starvation response

Abstract Multiple heavy metal contamination is an increasingly common global problem. Heavy metals have the potential to disrupt microbially mediated biogeochemical cycling. However, systems-level studies on the effects of combinations of heavy metals on bacteria are lacking. For this study, we focused on the Oak Ridge Reservation (ORR; Oak Ridge, TN, USA) subsurface which is contaminated with several heavy metals and high concentrations of nitrate. Using a native Bacillus cereus isolate that represents a dominant species at this site, we assessed the combined impact of eight metal contaminants, all at site-relevant concentrations, on cell processes through an integrated multi-omics approach that included discovery proteomics, targeted metabolomics, and targeted gene-expression profiling. The combination of eight metals impacted cell physiology in a manner that could not have been predicted from summing phenotypic responses to the individual metals. Exposure to the metal mixture elicited a global iron starvation response not observed during individual metal exposures. This disruption of iron homeostasis resulted in decreased activity of the iron-cofactor-containing nitrate and nitrite reductases, both of which are important in biological nitrate removal at the site. We propose that the combinatorial effects of simultaneous exposure to multiple heavy metals is an underappreciated yet significant form of cell stress in the environment with the potential to disrupt global nutrient cycles and to impede bioremediation efforts at mixed waste sites. Our work underscores the need to shift from single- to multi-metal studies for assessing and predicting the impacts of complex contaminants on microbial systems.

59 BASIC BIOLOGICAL SCIENCES↗

3-D Geologic Controls of Hydrothermal Fluid Flow at Brady geothermal field, Nevada, USA

In many hydrothermal systems, fracture permeability along faults provides pathways for groundwater to transport heat from depth. Faulting generates a range of deformation styles that cross-cut heterogeneous geology, resulting in complex patterns of permeability, porosity, and hydraulic conductivity. Vertical connectivity (a throughgoing network of permeable areas that allows advection of heat from depth to the shallow subsurface) is rare and is confined to relatively small volumes that have highly variable spatial distribution. This local compartmentalization of connectivity represents a significant challenge to understanding hydrothermal circulation and for exploring, developing, and managing hydrothermal resources. Here, we present an evaluation of the geologic characteristics that control this compartmentalization in hydrothermal systems through 3-D analysis of the Brady geothermal field in western Nevada. A published 3-D geologic map of the Brady area is used as a basis to develop structural and geological variables that are hypothesized to control or effect permeability or connectivity. The 3-D distribution of these variables is compared to the distribution of productive and non-productive fluid flow intervals along production wells and non-productive wells via principal component analysis (PCA). This comparison elucidates which geologic and structural variables are most closely associated with productive fluid flow intervals. Results indicate that production intervals at Brady are located: (1) within or near to known and stress-loaded macro-scale faults, and (2) in areas of high fault and fracture density.

58 GEOSCIENCES↗

Cellular remains in a ~3.42-billion-year-old subseafloor hydrothermal environment

Subsurface habitats on Earth host an extensive extant biosphere and likely provided one of Earth’s earliest microbial habitats. Though the site of life’s emergence continues to be debated, evidence of early life provides insights into its early evolution1. Here we present our discovery of exceptionally well-preserved ~3.40-billion-years-old carbonaceous filamentous microfossils that inhabited a paleo-subseafloor hydrothermal vein system of the Barberton greenstone belt (BGB) in South Africa. The distribution of the microfossils in growth position on the walls of the initial conduits created by low- temperature hydrothermal fluid and their morphological and chemical characteristics as investigated over a range of scales, are consistent with an assemblage of methanogens that thrived in an ultramafic volcanic substrate.

Cavalazzi, Barbara↗

3-D Geologic Controls of Hydrothermal Fluid Flow at Brady Geothermal Field, Nevada using PCA

In many hydrothermal systems, fracture permeability along faults provides pathways for groundwater to transport heat from depth. Faulting generates a range of deformation styles that cross-cut heterogeneous geology, resulting in complex patterns of permeability, porosity, and hydraulic conductivity. Vertical connectivity (a through going network of permeable areas that allows advection of heat from depth to the shallow subsurface) is rare and is confined to relatively small volumes that have highly variable spatial distribution. This local compartmentalization of connectivity represents a significant challenge to understanding hydrothermal circulation and for exploring, developing, and managing hydrothermal resources. Here, we present an evaluation of the geologic characteristics that control this compartmentalization in hydrothermal systems through 3-D analysis of the Brady geothermal field in western Nevada. A published 3-D geologic map of the Brady area is used as a basis to develop structural and geological variables that are hypothesized to control or effect permeability or connectivity. The 3-D distribution of these variables is compared to the distribution of productive and non-productive fluid flow intervals along production wells and non-productive wells via principal component analysis (PCA). This comparison elucidates which geologic and structural variables are most closely associated with productive fluid flow intervals. Results indicate that production intervals at Brady are located: (1) within or near to known and stress-loaded macro-scale faults, and (2) in areas of high fault and fracture density. This submission includes the published journal article detailing this work, the published 3-D geologic map of the Brady Geothermal Area used as a basis to develop structural and geological variables that are hypothesized to control or effect permeability or connectivity, 3-D well data, along which geologic data were sampled for PCA analyses, and associated metadata file. This work was done using existing R programs.

15 GEOTHERMAL ENERGY↗

Stability of Pt Skin Intermetallic Core Catalysts and Adsorption Properties for the Oxygen Reduction Reaction

Density functional theory calculations were used to determine the stability of metal slabs consisting of a Pt surface monolayer and intermetallic supporting layers made from combinations of six transition metal elements (Pt, Fe, Co, Ni, Cu, and Ag), as a model system for Pt skin intermetallic core nanoparticle catalysts. The stability of the slabs is largely determined by strain at the interface of the Pt skin and the subsurface intermetallic, which was described by a lattice matching parameter (r). The surface charge on the Pt skin was found to be correlated with the average electronegativity (EN) of the intermetallic core, so this average EN was used as a descriptor for how the electronic coupling (or ligand effect) affects adsorption energies. A total of 46 slabs were investigated in terms of their stability, from which 10 stable slabs were selected for further studies of adsorbate binding (OOH*, O*, and OH*) that are intermediates in the oxygen reduction reaction (ORR). The correlation between all three adsorption energies and descriptors r and EN was found. Using a linear fit between our descriptors and the calculated adsorption energies, the overpotential for the ORR was obtained as a function of r and EN, from which a volcano plot was produced. The volcano peak was found at r = 0.96 or at EN = 2.025. Interestingly, neither r nor EN was a sufficient single reactivity descriptor as the data points were well off the general trend in both linear fits; this implies that both the strain effect and the ligand effect influence the adsorption energies, although they are partly correlated. The (r, EN) target peak parameters were used to screen over 241 intermetallic combinations of transition metal elements as active ORR activity. In conclusion, this analysis identified 11 intermetallic compounds which can support a Pt skin to have a high predicted ORR activity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Predicting Thermal Performance of an Enhanced Geothermal System From Tracer Tests in a Data Assimilation Framework

Abstract Predicting the thermal performance of an enhanced geothermal system (EGS) requires a comprehensive characterization of the underlying fracture flow patterns from practically available data such as tracer data. However, due to the inherent complexities of subsurface fractures and the generally insufficient geological/geophysical data, interpreting tracer data for fracture flow characterization and thermal prediction remains a challenging task. The present study aims to tackle the challenge by leveraging a data assimilation method to maximize the utilization of information inherently contained in tracer data, and meanwhile maintain the flexibility to handle various uncertainties. A tracer data interpretation framework was proposed with the following three components integrated: (a) We use principal component analysis (PCA) to reduce the dimensionality of model parameter space. (b) We use ES‐MDA (ensemble smoother with multiple data assimilation) to invert for fracture aperture/flow fields and obtain posterior model ensembles for uncertainty quantification. Various data types are assimilated jointly to improve the predictive ability of the posterior ensemble. (c) The inverted fracture aperture fields are then incorporated into reservoir models to predict thermal performance. We developed a field‐scale EGS model to verify the ability of the framework to characterize highly heterogeneous fracture aperture/flow fields and predicting thermal performance. We also applied the framework to a mesoscale field experiment to demonstrate its potential application in real‐world geothermal reservoirs. The results indicate that the proposed framework can effectively retrieve fracture flow information from tracer data for thermal prediction and uncertainty quantification, and thus provide informative guidance for EGS optimization and risk management.

15 GEOTHERMAL ENERGY↗

Pulsed Thermal Tomography Nondestructive Examination of Additively Manufactured Reactor Materials and Components. Third Annual Progress Report

Additive manufacturing (AM) of high-strength corrosion resistance alloys for nuclear energy applications, such as stainless steel and Inconel, is currently based on laser powder bed fusion (LPBF) process. Some of the challenges with using LPBF method for nuclear manufacturing include the possibility of introducing pores into metallic structures. Probability of crack initiation at the pore depends on size, shape, and orientation of the defect. Pulsed Infrared Thermography Imaging (PIT) provides a capability for non-destructive evaluation (NDE) of sub-surface defects in arbitrary size structures. The PIT method is based on recording material surface temperature transients with infrared (IR) camera following thermal pulse delivered on material surface with flash light. The PIT method has advantages for NDE of actual AM structures because the method involves one-sided non-contact measurements and fast processing of large sample areas captured in one image. Following initial qualification of an AM component for deployment in a nuclear reactor, a PIT system can also be used for in-service nondestructive evaluation (NDE) applications. In this report, we describe recent progress in enhancing PIT capabilities in detecting microscopic subsurface defects in metals, and classifying shapes and orientation of pores in thermal images. For detection of microscopic defects in PIT imaging data, we have developed Spatial Temporal Denoised Thermal Source Separation (STDTSS) unsupervised machine learning (ML) image processing algorithm. We show that flat bottom hole (FBH) defects as small as 200µm in SS316 and IN718 specimens, can be detected with STDTSS algorithm. To the best of our knowledge, these are the smallest detected defects which are reported in literature. For classification of defects shapes, we have previously developed thermal tomography (TT) algorithm to obtain depth reconstructions of material defects from data cube of sequentially recorded surface temperatures. However, interpretation of TT images is non-trivial because of blurring with increasing depth. To address this challenge, we have developed a deep learning convolutional neural network (CNN) to classify size and orientation subsurface defects in simulated TT images.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Time-lapse electrical resistivity and induced polarization monitoring of a simulated ecosystem-scale coastal flooding experiment

Accurate estimation of the spatio-temporal variation in soil saturation and salt concentration is valuable for calibrating and evaluating Earth system models as well as detecting important eco-hydrological interactions. Such representative models play a vital role in investigating the stability of ecosystems facing changing hydrologic disturbance regimes. Ecosystem-scale field experiments are a useful way to capture hydrological disturbance and gain a deeper understanding of the complex mechanisms driving environmental changes. However, previous studies have not explored the use of quantitative imaging at a large spatial scale to reproduce breakthrough curves and inform system response and recovery following hydrologic disturbance events. Therefore, in this study, non-invasive geophysical methods were used to capture the spatio-temporal variation in subsurface saturation and salt concentrations during an ecosystem-scale coastal forest flooding experiment.

coastal↗

Quantifying Subsurface Biogeochemical Variability in a High Altitude Watershed During Winter Isolation

Shallow subsurface microbial and geochemical processes in watersheds are dynamic and responsive to seasonal and long-term environmental change. At high elevation or latitude, such changes may occur during winter months when normal sampling is impossible, impractical, unsafe, or may actually alter the processes being studied. Yet accurate modeling of biogeochemical parameters in these environments requires sampling that can capture events and persistent cold climate trends. Gaps in current models exist where data have not been collected during extended snow- or ice-covered periods. The goals of our exploratory study were to examine the biogeochemical processes occurring in a well-studied high altitude watershed during winter when the system is largely inaccessible. We tested the hypothesis that during snow cover in the East River (ER) watershed, episodic excursions of microbial community structure and biogeochemical processes and concentrations fluctuated from values extrapolated from pre- and post-snow time periods. Through our research we identified previously unresolved subsurface-surface transport phenomena and a heretofore under-appreciated methane cycle occurring during wintertime months, thus providing valuable insight into watershed processes. Our research showed that autonomous sampling techniques, first used in marine systems, is extremely versatile and can be adapted for use in harsh continental settings to examine temporal changes in freshwater chemistry and microbial community assemblage. Further, we found that by sampling during periods when systems are hard to access, will lend insight and better accountability for global budgets for methane, an important greenhouse gas. By strengthening our knowledge of wintertime biogeochemical relationships in high altitude watersheds our research leads to more complete descriptions of cryptic processes and ultimately will inform and improve the reactive transport models that simulate processes in these systems.

54 ENVIRONMENTAL SCIENCES↗

4D Proxy Imaging of Fracture Dilation and Stress Shadowing Using Electrical Resistivity Tomography During High Pressure Injections into a Dense Rock Formation

Fluid flow through fractured rock systems is governed in large part by the distribution, interconnectivity, and size of fracture apertures. In-situ stress is one of the primary factors controlling fracture aperture, and one that is altered significantly during high-pressure fluid injections or extractions. Interactions between stress, pore pressure, aperture, and fluid flow can result in complex and evolving poroelastic behavior with significant implications regarding the predictability and risk of developing and managing deep subsurface reservoirs (geothermal, fossil energy, and geologic carbon sequestration). In saturated crystalline rocks, bulk electrical conductivity is highly sensitive to porosity, and therefore to fracture aperture size and distribution. We demonstrate the use of time-lapse 3D electrical resistivity tomography for remotely monitoring stress induced changes in aperture distribution during high pressure injections into a fractured rock system at a scale of tens of meters. Results reveal a complex and continuously evolving stress field involving aperture dilations in the natural fracture system and aperture contractions in adjacent zones of shadow stress. Results provide information about the spatiotemporal changes in the system behavior and point to the potential of electrical imaging for autonomously and remotely monitoring evolving stress conditions by proxy through changes in bulk electrical conductivity.

electrical resisitivity tomography, stress, 4D Tim↗

Evaluation of distributed process-based hydrologic model performance using only a priori information to define model inputs

Fully distributed, integrated surface–subsurface hydrological models (ISSHMs) have seen renewed interest due to availability of better software, high performance computing facilities, and high-resolution, spatially extensive data products. ISSHMs are valuable as tools for advancing system understanding as they can resolve multiple processes defined on the plot scale including three-dimensional interaction of surface water and groundwater. Here, we evaluated the performance of an ISSHM, the Advanced Terrestrial Simulator (ATS), on seven diverse catchments across the continental US using widely available data products to define model inputs without calibration. We compare the ATS-simulated streamflow and evapotranspiration with gauge observations and MODIS-derived evapotranspiration, respectively. Using the Kling-Gupta Efficiency (KGE) as metric, ATS with default data products performed reasonably well at 6 of 7 catchments for streamflow. However, in one of those 6 catchments ATS had poor performance on baseflow and ATS’s overall performance was thus judged to be inadequate despite the acceptable KGE. ATS performance for evapotranspiration was good in all 7 catchments using default data products. In the two catchments where ATS streamflow performance using default data products was not acceptable, the performance was significantly improved by using local information on subsurface properties below the soil. We also compare the model-simulated streamflow and evapotranspiration with the Sacramento soil moisture accounting (SAC-SMA) model, a semi-distributed model that was calibrated on a catchment-by-catchment basis. Uncalibrated ATS performance is comparable to the calibrated SAC-SMA model in terms of streamflow while ATS performance is similar to or better (much better in certain catchments) in reproducing MODIS-derived evapotranspiration. Reasonably good performance of ATS without catchment-specific calibration provides new confidence in the ISSHM class of models and community data products as tools for advancing understanding of watershed function in a changing environment.

54 ENVIRONMENTAL SCIENCES↗

Short-Wave Infrared Upconverting Nanoparticles

Optical technologies enable real-time, noninvasive analysis of complex systems but are limited to discrete regions of the optical spectrum. While wavelengths in the short-wave infrared (SWIR) window (typically, 1700-3000 nm) should enable deep subsurface penetration and reduced photodamage, there are few luminescent probes that can be excited in this region. Here, we report the discovery of lanthanide-based upconverting nanoparticles (UCNPs) that efficiently convert 1740 or 1950 nm excitation to wavelengths compatible with conventional silicon detectors. Screening of Ln3+ ion combinations by differential rate equation modeling identifies Ho3+/Tm3+ or Tm3+ dopants with strong visible or NIR-I emission following SWIR excitation. Experimental upconverted photoluminescence excitation (U-PLE) spectra find that 10% Tm3+-doped NaYF4 core/shell UCNPs have the strongest 800 nm emission from SWIR wavelengths, while UCNPs with an added 2% or 10% Ho3+ show the strongest red emission when excited at 1740 or 1950 nm. Mechanistic modeling shows that addition of a low percentage of Ho3+ to Tm3+-doped UCNPs shifts their emission from 800 to 652 nm by acting as a hub of efficient SWIR energy acceptance and redistribution up to visible emission manifolds. Parallel experimental and computational analysis shows rate equation models are able to predict compositions for specific wavelengths of both excitation and emission. These SWIR-responsive probes open a new IR bioimaging window, and are responsive at wavelengths important for vision technologies.

Qi, Xiao↗

Electrical Resistivity Tomography based monitoring of stress perturbations to optimize placement of high-precision strain meters

The Center for Understanding Subsurface Signals and Permeability is a new U.S. Department of Energy Earthshot Center focused on understanding and predicting the long-term evolution of permeability in enhanced geothermal systems. The center will use a highly instrumented testbed within the Sanford Underground Research Facility to conduct field scale experiments that elucidate and test capabilities to simulate geochemical-geomechanical interactions and permeability evolution. Here we demonstrate initial developments using previously collected electrical resistivity tomography (ERT) monitoring data with high-performance multi-physics modelling advancements to inform the optimal location of two new monitoring boreholes. Specifically, ERT monitoring data collected during shear stimulation testing shows marked responses to changes in stress during borehole pressurization. We demonstrate how the same response is being simulated, ultimately to train a machine-learning algorithm to estimate rock properties and enable enhanced prediction of stress and strain responses anticipated during future testing campaigns.

Stress, EGS, CUSSP, 3D Electrical Imaging↗

Techno-economic feasibility of geothermal energy production using inactive oil and gas wells for district heating and cooling systems in Tuttle, Oklahoma

Oil and gas wells have been drilled in the United States and often encounter hot co-produced water possibly suitable for geothermal direct-use applications. Here this study evaluated geothermal resources, heating and cooling demand, and techno-economic potential of four oil and gas wells repurposed for geothermal energy production to serve two public schools and 250 nearby houses in Tuttle, Oklahoma, USA. Subsurface geology in the study area primarily consists of Permian to Mississippian sandstone and limestone formations approximately from 2 km to 3.5 km depth. With a geothermal gradient of 25.3 °C/km, the temperature of geothermal energy production was estimated as 68.2 °C at 2.1 km. The potential of calcite scaling within boreholes and pipes was minimal. Using the characterized reservoir properties, techno-economic analyses were conducted for evaluating levelized costs of geothermal energy production for direct-use heating and cooling and district heating systems with five different production scenarios. Heating and cooling demand in two schools and 250 houses and additional 0.3 MW and 0.6 MW base heating loads for potential geothermal direct-use applications in the study area were also incorporated in the modeling. Results indicated that the levelized cost of heating of geothermal district heating systems utilizing existing wells were significantly lower than those of conventional geothermal energy systems that include well drillings. The geothermal energy production with doublet and quartet configurations was generally sufficient to supply total annual heating demand in the district, while a peaking boiler was used to supply peak loads. Doublet configurations showed higher production temperature with a lower thermal drawdown than the quartet configurations. The doublet system at 2.1 km production depth with 9.3 kg/s flow rate represented the optimal levelized cost of $71/MWh and 91.2% utilization factor. These results imply that the Tuttle geothermal district heating system is techno-economically feasible and comparable to geothermal district heating systems in the United States and the regional natural gas rate for residentials by recycling inactive oil and gas wells with no drillings.

15 GEOTHERMAL ENERGY↗