SEARCH · Engineering Papers
Results for “geothermal fields”
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.
From Fluid Flow to Coupled Processes in Fractured Rock: Recent Advances and New Frontiers
Abstract Quantitative predictions of natural and induced phenomena in fractured rock is one of the great challenges in the Earth and Energy Sciences with far‐reaching economic and environmental impacts. Fractures occupy a very small volume of a subsurface formation but often dominate fluid flow, solute transport and mechanical deformation behavior. They play a central role in CO 2 sequestration, nuclear waste disposal, hydrogen storage, geothermal energy production, nuclear nonproliferation, and hydrocarbon extraction. These applications require predictions of fracture‐dependent quantities of interest such as CO 2 leakage rate, hydrocarbon production, radionuclide plume migration, and seismicity; to be useful, these predictions must account for uncertainty inherent in subsurface systems. Here, we review recent advances in fractured rock research covering field‐ and laboratory‐scale experimentation, numerical simulations, and uncertainty quantification. We discuss how these have greatly improved the fundamental understanding of fractures and one's ability to predict flow and transport in fractured systems. Dedicated field sites provide quantitative measurements of fracture flow that can be used to identify dominant coupled processes and to validate models. Laboratory‐scale experiments fill critical knowledge gaps by providing direct observations and measurements of fracture geometry and flow under controlled conditions that cannot be obtained in the field. Physics‐based simulation of flow and transport provide a bridge in understanding between controlled simple laboratory experiments and the massively complex field‐scale fracture systems. Finally, we review the use of machine learning‐based emulators to rapidly investigate different fracture property scenarios and accelerate physics‐based models by orders of magnitude to enable uncertainty quantification and near real‐time analysis.
Lessons Learned from AskGDR: Usage and Impact Analysis of the Geothermal Data Repository's AI Research Assistant: Preprint
In October of 2024, the Department of Energy's (DOE) Geothermal Data Repository (GDR) team officially launched AskGDR, an AI research assistant resulting from the integration of a Large Language Model (LLM) with the metadata and supporting documents associated with GDR datasets. AskGDR allows GDR users to ask deeper questions about the origin of datasets, the methods used to collect them, and the findings they help support. Using Retrieval Augmented Generation (RAG), AskGDR can be used to summarize findings spread across dozens of papers and technical reports or to extract relevant information describing a single data field. However, generative AI is experimental. The National Renewable Energy Laboratory (NREL) has been collecting metrics on AskGDR and documenting lessons learned during its deployment. This paper will outline the efficacy and impact of AskGDR through analysis of its use, operating costs, number and types of questions asked, and the quality of answers provided.
Spontaneous crack healing in calcite reveals the influence of dynamic strain evolution and surface chemistry
The mechanics of fracture healing in calcite remain poorly constrained yet are fundamental to managing fluid transport in geothermal reservoirs and hydrocarbon systems. Here, we apply microfocused synchrotron Laue X-ray diffraction and infrared spectroscopy to investigate subcritical crack healing in a 1 mm-thick calcite crystal subjected to controlled loading in a double-torsion device. Over a 44-hour period following load removal, we map the evolution of residual strain fields surrounding the crack tip and observe a progressive increase in compressive strain perpendicular to the crack plane accompanied by infrared spectroscopic signatures that reveal enhanced accumulation of water at the healed interface. The correlation between strain evolution and surface chemistry suggests that spontaneous crack healing in calcite is driven by dynamic anelastic relaxation coupled with irreversible fluid-mineral interactions. These findings offer insight into time-dependent crack closure processes in carbonates and highlight the role of chemically-mediated plasticity in subsurface fracture evolution.
Well-Controlled Hydro-thermo-chemical Stimulation of a Coso Geothermal Well to Enhance Productivity
The goal of this project was to demonstrate controlled fracturing and permeability increase to enhance resource production and power generation at Coso. It was planned to stimulate at least two intervals in the Well 73-18RD with zonal isolation using advanced high-temperature packer technology. The anticipation was that a hydraulic fracture would form in one zone and in the second zone natural fractures would reopen, shear and extend, thereby increasing fluid communication in the field. In both cases, the intent was to carefully design the stimulation to take advantage of the relative roles of hydraulic pressure, thermal stress and chemical (“hydro-thermo-chemical”) dissolution. The project had planned to use the Welltec Annular Barrier (WAB) system developed as part of another DOE-funded project which was on schedule to be ready for field deployment. If successful, the Coso Operating Company would incorporate this approach into their ongoing well maintenance in other wells increasing the total MWs produced. The project had multiple tasks and the objective of Task 1 was to develop a workover plan, to build a geologic model, to design a stimulation program, and to work with the operating and services companies to develop a detailed budget and obtain the necessary permits and sundry notices. In Phase 2 the workover plan would be implemented, the designed stimulation would be carried out, and subsequently the well would logged and tested to measure the impact of stimulation. Per sponsor request, a Pre-Task 1 activity with a Go/No Go decision point was scheduled in Task 1.
Heart Shape to Fracture Distance: Characterizing Hydraulic Fracture Propagation before Hits
Estimating the distance from the hydraulic fracture tip to the monitor well can be useful for fracture characterization, well spacing optimization, and preventing parent-child well interference. A heart-shaped signal is referred to as the extensional precursor of a fracture hit recorded by crosswell strain measurements and can serve as a vital tool for such estimation. This study incorporates the 3D displacement discontinuity method (DDM) to understand the impact of fracture geometry and monitor well offset on the heart-shaped signal’s characteristics. Results from numerical simulation and analytical solutions reveal a strong linear correlation between the spatial extent of the heart-shaped signal and the fracture tip distance. This relationship was further developed to predict tip distance using field data from the Hydraulic Fracture Test Site 2 (HFTS2). A reasonable approximation result from field data further validates the methodology. In addition, it is worth noting that the estimation accuracy depends on the ratio between fracture dimension and tip distance. The findings of this study offer a novel approach for real-time monitoring and characterizing hydraulic fracture propagation, which can be further used for well spacing optimization in unconventional and enhanced geothermal system reservoir development, as well as caprock integrity monitoring for carbon sequestration projects.
Characterizing Hydraulic Fracture Propagation Before Fracture Hits
Estimating the distance from hydraulic fracture tip to monitor well can be very useful for fracture characterization, well spacing optimization, and preventing parent-child well interference. Heart-shape signal is referred to as the extensional precursor of fracture hit recorded by cross-well strain measurements and can be served as a vital tool to make such estimation. This study incorporates the 3D Displacement Discontinuity Method to understand the impact of fracture geometry and monitor well offset on the heart-shape signal’s characteristics. Results from numerical simulation and analytical solutions reveal a strong linear correlation between the spatial extent of heart shape signal and fracture tip distance. This relationship was further developed to predict tip distance using field data from the Hydraulic Fracture Test Site 2. A reasonable approximation result from field data further validates the methodology. In addition, it is worth noting that the estimation accuracy depends on the ratio between fracture dimension and tip distance. The findings of this study offer a novel approach for real-time monitoring and characterizing hydraulic fracture propagation, which can be further used for well spacing optimization in unconventional and Enhanced Geothermal System reservoir development, as well as cap rock integrity monitoring for carbon sequestration projects.
Geothermal Deep Direct-Use for Low-Carbon Heating: A Case Study at Cornell University
Cornell University is developing a geothermal deep direct-use system to provide reliable, low carbon, and cost-competitive heating for its main campus in Ithaca, New York. This geothermal system, also referred to as Earth Source Heat, would provide the majority of the campus heating load and represents a key pillar for the university to become carbon neutral by 2035. Over the last 10 years, numerous desktop, seismic, and geophysical studies have been undertaken at Cornell to characterize the local subsurface, predict long-term reservoir performance, design the integration into the existing district heating network, and estimate capital and levelized costs. Recently, the university received funding for drilling an observatory well to obtain in-situ subsurface data (e.g., reservoir temperatures, pressures, fluid chemistry, stress field), significantly lowering subsurface uncertainty and de-risking the project. Drilling of the test well is planned for fall 2021. This paper provides an overview and status of the project, including a summary of past studies and a description of the proposed well and planned testing activities.
Prediction and Analysis of Utah FORGE Injection Activities using a Coupled Thermo-hydro-mechanical and Earthquake (THM+E) Modeling Workflow
A coupled thermo-hydro-mechanical (THM) numerical workflow that is capable of modeling seismic slip is critical for the successful development of enhanced geothermal systems (EGS). By integrating key physical processes, this workflow enables accurate simulation of temperature and pressure diffusions, stress changes, and induced seismicity. As a result, it serves as a vital tool for predicting induced seismicity and optimizing reservoir stimulation strategies. The Utah FORGE (Frontier Observatory for Research in Geothermal Energy) project, located near Milford, Utah, is a U.S. Department of Energy initiative aimed at advancing EGS technology. In April 2024, eight new stimulation stages (Stages 3R-10) were conducted in well 16A (injection well) subsequent to the first series of stimulation (Stages 1-3) performed in April, 2022. To monitor the induced seismicity, geophones were deployed in wells 58-32, 56-32, and 78B-32, while fiber optic cables were also installed in wells 16B, 78-32, and 78B-32 to collect microseismic data and detect frac hits Preliminary analyses of microseismic catalogs and fiber optic data suggest that the stimulated fractures in Stages 3R–6 closely align with that generated during Stage 3, indicating that the new stimulations were likely reactivating the previously stimulated fracture. To better understand the underlying process, a comprehensive modeling approach that can accurately capture thermal, hydrological, mechanical, and seismic responses is essential. In this work, we propose and utilize a coupled thermo-hydro-mechanical and earthquake (THM+E) simulation workflow to numerically investigate the stimulation activities on well 16A. The specific objective is to confirm whether the new stimulation stages (Stages 3R–6) reactivated fractures previously stimulated during Stage 3. For this purpose, we perform THM+E simulations individually for Stages 3, 3R, 4, and 5, incorporating the discrete fracture networks (DFNs) created by the plane-fitting technique based on the microseismic catalogs. The simulation workflow consists of two separate models: a THM model and an earthquake model, coupled in a one-way manner. Detailed descriptions of the workflow are provided in Section 3. Simulation results are presented in terms of injection pressure, permeability evolution, and predicted seismic catalogs, which are then compared with field data for further analyses. This report is structured as follows. In Section 2, we present detailed analyses of the field data and propose the hypothesis that the new stimulation stages (Stages 3R–6) were probably reactivating the previously stimulated fractures in Stage 3. In Section 3, we introduce the coupled THM+E workflow and the problem setup to validate our hypothesis, followed by the simulation results for each stage in Section 4. Meanwhile, discussions are included to analyze the model predictions and their comparison with field data. Lastly, we conclude the report and outline future plans in Section 5.
WHOLESCALE Catalog of Rock Samples at San Emidio Nevada collected in January 2021 Version 2.0
This submission includes an update to the WHOLESCALESamples2021January.xls file where the strikes and dips initially reported have been corrected to comply with the right hand rule. The updated excel file and a link to the original submission are included in this report. The original submission "WHOLESCALE Catalog of Rock Samples at San Emidio Nevada collected in January 2021" contains information on thirty-six rock samples collected from San Emidio, Nevada during January, 2021 for Subtask 2.3 of the WHOLESCALE project. The following resources include a .zip of rock sample photos taken in the field, a .zip of rock sample photos taken in the laboratory at UW-Madison, and an excel catalog of rock samples with information on sample name, rock type, coordinates of sample location, structural measurements, field notes, observations for plug preparation (e.g., weathering, ability to be cut and cored), and rock descriptions. It should be noted that not every sample was photographed in the field. Names and descriptions of rock formation units are taken from Rhodes et al. (2011). The README.txt file is a description of this submission.
Interpretation of hydrothermal conditions, production-injection induced effects, and evidence for enhanced geothermal system-type heat exchange in response to >30 years of production at Roosevelt Hot Springs, Utah, USA
The Roosevelt Hot Springs hydrothermal system is located at the base of the Mineral Mountains in southwestern Utah on the eastern side of the Basin and Range. Hydrothermal activity is related to relatively recent bimodal magmatism, and the system is hosted in coarsely crystalline rock made of Oligocene–Miocene granitoids and Precambrian gneiss. The hydrothermal plume covers ~5 km 2 , with a maximum temperature of 268 °C at ~750 m depth, and a vertically extensive fault-fracture mesh east of the Opal Mound fault controls the upflow of hydrothermal fluids. Power generation (currently 38 MWe gross) began in 1984, and up through 2016, four wells were used for fluid production, and three wells were used for edge-field injection. Chemical analyses of produced fluids show that modern reservoir fluid compositions are similar to but more concentrated than those at the start of production, having near-neutral pH, total dissolved solids of 7000–10,000 mg/kg, and ionic ratios of Cl/HCO 3 ~50–100, Cl/SO4 ~50–100, and Na/K ~4–5. Chemical geothermometers indicate equilibration temperatures that mainly range between 240° and 300 °C. Early production induced a steep drop in pressure (~3.0–3.5 MPa), which was accompanied by a 250–300 m lowering of piezometric levels in wells and development of a shallow steam zone across the system. Hydrothermal fluid compositions evolved continuously in response to production-related steam-loss and injection breakthrough, which is reflected by gradual increases in chloride of up to 35% and stable isotope ratios of up to ~2‰δ 18 O and ~10‰δD. Simple mixing model calculations suggest that there has been a significant amount, ~10–20 MW th , of sustained multi-decadal heat mining and enhanced geothermal system (EGS)–type heat transfer by the injectate as it returns to the production zone. Overall, the two factors that have sustained long-term power production (currently 38 MWe gross) are the increased upflow of deep chloride water and, to a lesser extent, the mining of heat at <1 km depth.
A Novel Approach to Map Permeability Using Passive Seismic Emission Tomography
Newly acquired magnetotelluric data and passive seismic data collected with tightly spaced geophone arrays are combined with historic drilling, active seismic, and potential fields data to generate 3-D permeability maps. A cooperative inversion methodology has been developed using active seismic, magnetotelluric, and gravity data in order to produce more robust velocity models for passive seismic data processing without requiring expensive 3-D active seismic surveys. The cooperative inversion estimates velocities from other geophysical data where no prior seismic velocity information is available at two geothermal sites in Nevada: San Emidio and Crescent Valley.
A Modular Mechanism for Downhole Weight-on-Bit and Torque Reaction in Small Diameter Boreholes
Drilling systems that use downhole rotation must react torque either through the drill-string or near the motor to achieve effective drilling performance. Problems with drill-string loading such as buckling, friction, and twist become more severe as hole diameter decreases. Therefore, for small holes, reacting torque downhole without interfering with the application of weight-on-bit, is preferred. In this paper, we present a novel mechanism that enables effective and controllable downhole weight on bit transmission and torque reaction. This scalable design achieves its unique performance through four key features: (1) mechanical advantage based on geometry, (2) direction dependent behavior using rolling and sliding contact, (3) modular scalability by combining modules in series, and (4) torque reaction and weight on bit that are proportional to applied axial force. As a result, simple mechanical devices can be used to react large torques while allowing controlled force to be transmitted to the drill bit. We outline our design, provide theoretical predictions of performance, and validate the results using full-scale testing. The experimental results include laboratory studies as well as limited field testing using a percussive hammer. These results demonstrate effective torque reaction, axial force transmission, favorable scaling with multiple modules, and predictable performance that is proportional to applied force.
A combined experimental and modelling study of granite hydrothermal alteration
Geochemical reactions can induce significant changes of rock reservoir porosity and permeability via mineral dissolution and precipitation processes, affecting the long-term fluid behaviour within various geological systems. Here, the understanding and quantification of these reactions rely on field and experimental studies and on the predictions of reactive transport models. The present study was aimed at assessing the extent to which current geochemical models integrating available mineral dissolution/precipitation rate equations can reproduce the experimental data obtained from 4 to 17-day long hydrothermal alteration experiments of a muscovite-biotite granite and, thus, help provide an accurate description of the evolution of geothermal systems within granitic reservoirs. The experiments were conducted at a constant temperature of 180 °C and over an aqueous fluid pH range of 2 to 8.5, using both mixed-flow and static batch reactors. Modelled major element (K, Al, Si, Ca, and Mg) concentrations were generally in satisfactory agreement with the corresponding measured elemental fluxes – the differences between modelled and experimental values were generally within a factor of 5 – and the predicted identity and mass of formed secondary phases were consistent with the microscopic observations of the reacted solids. However, larger differences between measured and modelled element concentrations were observed when significant amounts of secondary phases formed, notably at pH 2 to 3, and for longer-term batch experiments. Much of this concentration difference stems from the underestimation of the amounts of Al-phases formed at acid to near-neutral pH. Although an idealized rock composition was considered, the observed mismatch between model calculations and experimental data can be attributed to inadequate mineral precipitation reaction rates and a poor description of reactive surface areas in existing geochemical modelling codes. More accurate quantification of precipitation kinetics, including nucleation and growth, and improved descriptions of the temporal change of mineral surface area would enhance the predictive capabilities of reactive transport models and benefit, particularly, the efforts aimed at increasing the sustainability of EGS reservoirs.
A Non-Invasive Approach for Elucidating the Spatial Distribution of In-Situ Stress in Deep Subsurface Geologic Formations Considered for CO 2 Storage (Final Report)
The primary objectives of this project were to: develop a new, non-invasive method to define the orientation and magnitude of in-situ stresses in deep geologic formations (depths that exceed 1500 meters) considered for carbon dioxide (CO 2 ) sequestration; to demonstrate the method at a real field site; and, to extend current technical capabilities and reduce costs of acquiring information about deep-stress fields. This investigation was segmented into five tasks, including Project Management (Task 1) and four technical tasks (Tasks 2 through 5). The results of the technical tasks are described in this document.
Control development and sizing analysis for a 5th generation district heating and cooling network using Modelica
5th generation district heating and cooling system (5GDHC) are a relatively new concept. They use a single district loop near ambient temperature to provide heating and cooling. This paper improves on the modelling and control of a 5GDHC system called the reservoir network. The study updates the sewage heat exchanger plant model to more realistically represent seasonal changes, uses refined pump models with variable efficiency, introduces a ground coupled district pipe model to consider the inertia of the district network and implements a new control strategy for geothermal storage and sewage heat exchanger. The new approach reduced operating costs, mainly due to pumping cost for storage, sewage heat exchanger plant and distribution pump, while increasing the overall robustness of the approach in different sizing conditions. Thanks to the new controller, the pumping consumption was reduced by 21% with respect to the original baseline. Furthermore, the new control makes the system take better advantage of design changes, when reducing borehole field size and increasing the sewage heat exchanger size, the pumping energy savings become 29% with respect to the original baseline. Lastly, borehole field temperature stability was analyzed through 40 years of simulation.
Utah FORGE: Sanvean Technologies Drilling Data from Well 16B(78)-32
Included here is Sanvean Technologies bit sensor data amalgamated with data from National Oilwell Varco's (NOV) BlackBox tool for Reed Hycalog bits used during drilling of Well 16B(78)-32. The dataset contains information collected at the bit while drilling including rate of penetration (ROP), top drive torque, and bit box temperature. The data was recorded at the bit box and top sub of the motor. RPM was measured by onboard gyro recording continuously in each sensor, and shock levels were also recorded on X, Y and Z axis. This data was merged with EDR in time format and saved in file sets (the zipped files) then output into CSV files. Please note: fields in the CSV files, such as the date field, may need to be formatted to display properly. There is an additional zipped folder in each dataset that is password protected. Sanvean GameChanger Viewer software must be used to view this password protected data. Information on how to use and download this free software is also included here.
Utah FORGE 2-2439v2: Characterizing In-Situ Stress with Laboratory Modelling and Field Measurements - 2024 Annual Workshop Presentation
This is a presentation on A Multi-Component Approach to Characterizing In-Situ Stress at the Utah FORGE Site: Laboratory Modelling and Field Measurements project by The University of Pittsburgh, presented by Andrew Bunger. The project characterizes the stress in the Utah FORGE EGS reservoir using three methods: Method 1: Demonstrate complimentary laboratory rock-core stress estimation combined with Machine Learning approach for measuring in-situ stress from field sonic log data; Method 2: Complete field based in-situ measurement (mini-frac); and Method 3: Develop a mechanics-based method for connection near wellbore stress measurements to stresses away from the well-bore. This presentation was featured in the Utah FORGE R&D Annual Workshop on August 14, 2024.