Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geothermal Systems”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Monitoring Fracture Hydromechanical Evolution in the Lab and Field Using Unsupervised Metric Learning

Fractures evolve in time through thermal‐hydraulic‐mechanical‐chemical (THMC) processes that alter their long‐range hydraulic transport properties and modify subsurface behavior and activities. The location of subsurface fractures makes it necessary to use remote sensing techniques such as passive or active seismic monitoring for fracture characterization. In this paper, we develop a machine learning approach to monitor the evolution of fracture properties using passive seismic sources in a laboratory setting and using active seismic monitoring from the Sanford Underground Research Facility in Lead, South Dakota, at a depth of 1.25 km in amphibolite rock during stimulation of natural fractures as well as during induced fracturing. The unsupervised metric learning technique applies tandem neural networks (twin (Siamese) or triplet) with contrastive loss and adaptive margins to track slowly varying systems for which class or similarity labels are not available. The approach adopts locality‐sensitive hashing to divide time‐ordered contiguous data into an arbitrary number of pseudo‐classes. Contrastive‐loss training with many hash bins generates an evolving latent‐space trajectory. This approach enables unsupervised metric learning for seismic data stacks under the condition of contiguous state sampling and slowly varying fracture properties. The displacement discontinuity theory provides a mechanistic foundation for the fracture‐dependent trajectories that are related to relaxation of fractures with time‐dependent specific stiffness responding to changes in stress or fluid saturation.

02 PETROLEUM↗
DOE OSTI2019

Fracture Tracer Injection Response to Pressure Perturbations at an Injection Well

The EGS Collab project constructed an approximately ten-meter scale field site where fracture stimulation and flow/transport models can be validated against controlled, in-situ experiments. The first multi-well experimental site was established at the 4850 level of the Stanford Underground Research Facility (SURF) in the Homestake Mine located in South Dakota. Hydraulic fractures were created at an injection well drilled sub-horizontal from the drift. A flow system was established in one set of fractures by injection water at approximately 400 ml/min between a set of packers 164 feet from the drift wall in the injection well through a hydraulically stimulated fracture. Injected water was recovered from five locations in 4 nearby wells. From the end of October to the middle of November of 2018, a series of fracture characterization tests were conducted using a series of 10 tracer tests (7 which used C-dot and chloride as conservative solute tracers) to assess the flow pathway in the stimulated fracture. The injected tracers were detected in three of the five water production locations where the total water recovery ranged from approximately 50 to 80% of the injected water depending on the day the test was being conducted. Analysis of a series of tracers during this two-week period suggest a large change in the flow fracture field occurred during this testing period. A comparison of the tracer breakthrough curves at the production well showed a marked decrease in the initial and peak concentration over time, whereas the OB well exhibited an increase in the initial and peak concentration arrival during this fracture characterization testing period. These changes are believed to be in response to a number of higher-pressure short-term injections at the injection well in early November (2nd to the 6th). Results of this testing suggest that the fracture flow pathways can be altered as a result of the pressure perturbations in the injection well on the integrity of adjacent monitoring wells. Results from the EGS Collab project will support the DOE Geothermal Technology Office FORGE and other EGS development efforts.

15 GEOTHERMAL ENERGY↗
DOE OSTI2015

Dynamics within geyser conduits, and sensitivity to environmental perturbations: Insights from a periodic geyser in the El Tatio geyser field, Atacama Desert, Chile

Despite more than 200. years of scientific study, the internal dynamics of geyser systems remain poorly characterized. As a consequence, there remain fundamental questions about what processes initiate and terminate eruptions, and where eruptions begin. Over a one-week period in October 2012, we collected down-hole measurements of pressure and temperature in the conduit of an exceptionally regular geyser (132. s/cycle) located in the Chilean desert. Here, we identified four stages in the geyser cycle: (1) recharge of water into the conduit after an eruption, driven by the pressure difference between water in the conduit and in a deeper reservoir; (2) a pre-eruptive stage that follows the recharge and is dominated by addition of steam from below; (3) the eruption, which occurs by rapid boiling of a large mass of water at the top of the water column, and decompression that propagates boiling conditions downward; and (4) a relaxation stage during which pressure and temperature decrease until conditions preceding the recharge stage are restored. Eruptions are triggered by the episodic addition of steam coming from depth, suggesting that the dynamics of the eruptions are dominated by geometrical and thermodynamic complexities in the conduit and reservoir. Further evidence favoring the dominance of internal processes in controlling periodicity is also provided by the absence of responses of the geyser to environmental perturbations (air pressure, temperature and probably also Earth tides).

Bubble trap↗
DOE OSTI2015

Geomechanical simulation of the stress tensor rotation caused by injection of cold water in a deep geothermal reservoir

We present a three-dimensional thermohydromechanical numerical study of the evolution and distribution of the stress tensor within the northwest part of The Geysers geothermal reservoir (in California), including a detailed study of the region around one injection well from 2003 to 2012. Initially, after imposing a normal faulting stress regime, we calculated local changes in the stress regime around injection wells. Our results were compared with previously published studies in which the stress state was inferred from inverting the focal plane mechanism of seismic events. Our main finding is that changes in stress tensor orientation are caused by injection-induced progressive cooling of the reservoir, as well as by the seasonal variations in injection rate. Because of the gravity flow and cooling around a liquid zone formed by the injection, the vertical stress reduction is larger and propagates far below the injection well. At the same time, the horizontal stress increases, mostly because of stress redistribution below and above the cooling area. These two phenomena cause the rotation of the stress tensor and the appearance of a strike-slip regime above, inside, and below the cooling area. The cooling and the associated rotation of the stress regime can play a significant role in the observed long-term deepening of the microseismicity below active injection wells.

58 GEOSCIENCES↗
DOE OSTI2023

Enhancing flexibility for climate change using seasonal energy storage (aquifer thermal energy storage) in distributed energy systems

Long-term energy storage is expected to play a vital role in the deep decarbonization of building energy sectors, while enhancing the flexibility of buildings to withstand future climate variations. However, it is challenging to design distributed multi-energy systems (DMES) while taking into account the uncertainties introduced by climate change, since stochastic optimization of such systems is difficult. The present study introduces a stochastic optimization model to address this bottleneck, taking into account DMES including aquifer thermal energy storage (ATES) as the long-term thermal storage. For the first time, a novel optimization algorithm links ATES with the DMES optimization model with the support of a simplified geotechnical model. Subsequently, a case study was conducted, focusing on a residential district in Chicago where the impact of future climate condition, energy demand, and solar and wind energy potentials were evaluated using Weather Research and Forecasting (WRF) data (up to 2080) and the EnergyPlus model. The study revealed that ATES is an attractive way to improve the renewable energy penetration level and minimize the dependence on fossil fuels with reasonable support from the grid to assist the fluctuations in both demand and generation. Furthermore, ATES notably reduces fuel consumption and dependence while greatly enhancing the flexibility of the energy system to withstand fluctuations in demand and renewable energy generation brought by future climate variations. These qualities will make ATES an important part of distributed energy systems, even though it is not currently the lowest cost alternative due to lack of technology maturity. Furthermore, the design platform introduced in the present study can be used to design DMES enhancing flexibility to accommodate future climate variations.

25 ENERGY STORAGE↗
DOE OSTI2016

Effects of Asperity Distribution on Fluid Flow and Induced Seismicity During Deep Geothermal Exploitation

This work investigates the injection-induced seismic response of a heterogeneous fault plane, featuring low-permeability asperities embedded into a high-permeability damage zone. We simulate the pressure evolution with a hydrogeological simulator, accounting for the heterogeneous fault plane. Seismicity occurs then on the asperities, represented as unstable patches reactivating by means of the Mohr-Coulomb criterion. The hydrological and seismic modules are implicitly coupled to account for effects of asperity reactivation on the permeability. Results show that permeability changes may cause at a later time a change in seismicity propagation. Finally, we also investigated such effects by varying the density of asperities.

58 GEOSCIENCES↗
DOE OSTI2015

The Northwest Geysers EGS Demonstration Project, California

Our project goal is to demonstrate the feasibility of stimulating a deep high-temperature reservoir (HTR) (up to 400 °C, 750 °F). There were two previously abandoned wells, Prati State 31 (PS-31) and Prati 32 (P-32), reopened and deepened to be used as an injection and production doublet to stimulate the HTR. The deepened portions of both wells have conductive temperature gradients of 10 °F/100 ft (182 °C/km), produce connate native fluids and magmatic gas, and the rocks were isotopically unexchanged by meteoric water. The ambient temperature meteoric water injected into these hot dry rocks has evidently created a permeability volume of several cubic kilometers as determined by seismic monitoring. Preliminary isotopic analyses of the injected and produced water indicate that 50–75% of the steam from the created EGS reservoir is injection-derived.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Induced seismicity red-light thresholds for enhanced geothermal prospects in the Netherlands

Prospects for geothermal energy in the Netherlands have renewed concerns around induced earthquakes. Risks from induced earthquakes are managed by traffic light protocols (TLPs), where the red-light is chosen as the stop-point before exceeding a tolerance to risk. Here, we simulate post-shut-in earthquake scenarios based on realistic information for the Netherlands. We focus on three risk metrics: aggregates like nuisance and damage impacts and also local personal risk (LPR) – a likelihood of building collapse fatality for an individual. Our results show that the severity of these risks varies spatially by orders of magnitude. Prior induced seismicity (e.g., the 2012 Huizinge event) provides a reference baseline to calibrate the Dutch earthquake risk tolerances. We find that these calibrated risk tolerances are similar to those observed in North America, suggesting an underlying sociological ‘license to operate.’ Furthermore, the use of calibrated risk tolerances results in nuisance concerns completely eclipsing the other two metrics. We compare our results to a hypothetical Groningen geothermal operation and find that our approach sets red-light thresholds approximately one magnitude unit below the ML 3.6 Huizinge event. Overall, our results provide a first-order recommendation for red-light thresholds and proactive management of Dutch enhanced geothermal induced seismicity.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Monitoring spatiotemporal evolution of fractures during hydraulic stimulations at the first EGS collab testbed using anisotropic elastic-waveform inversion

The EGS Collab project acquired continuous active-source seismic monitoring (CASSM) data before, during, and after hydraulic stimulations at the first testbed at the depth of 4850 ft (1478 m) at the Sanford Underground Research Facility in Lead, South Dakota, for monitoring fracture creation and evolution. CASSM acquisition was conducted using 24 hydrophones, 18 accelerometers, and 17 piezoelectric sources within four fracture-parallel wells and two orthogonal wells. 3D anisotropic traveltime tomography and anisotropic elastic-waveform inversion of the campaign cross-borehole seismic data show that the rock within the stimulation region is a heterogeneous horizontal transverse isotropic medium. Here we use these inversion results as the initial models and apply 3D anisotropic first-arrival traveltime tomography and 3D anisotropic elastic-waveform inversion to the CASSM data acquired after each stimulation in May, 2018 and December, 2018. We observe the spatiotemporal evolution of seismic velocities and anisotropic parameters caused by hydraulic fracture stimulations, showing the regions of rock alternation caused by hydraulic fracture stimulation.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Numerical simulation of flow and mixing in fracture intersections

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

58 GEOSCIENCES↗
DOE OSTI2024

High resolution identification and quantification of diffuse deep groundwater discharge in mountain rivers using continuous boat-mounted helium measurements

Discharge of deeply sourced groundwater to streams is difficult to locate and quantify, particularly where both discrete and diffuse discharge points exist, but diffuse discharge is one of the primary controls on solute budgets in mountainous watersheds. The noble gas helium is a unique identifier of deep groundwater discharge because groundwater with long residence times is commonly enriched in helium. In this study, a portable mass spectrometer was used to measure longitudinal variation in dissolved helium concentrations in two mountainous rivers at high spatial resolution not feasible with traditional sampling techniques. Helium profiles were then simulated using a mass-balance model to quantify longitudinal variation in groundwater discharge to the receiving rivers. Results indicate helium concentrations were enriched by multiple orders of magnitude above atmospheric equilibrium in both rivers and that this persisted for up to 18 km below observed pulse inputs in the Colorado River. Helium mass-balance models match observed longitudinal patterns with the exception of sharp initial increases in helium observed in the rivers. Increased longitudinal groundwater discharge rates correspond to mapped geologic structures in both watersheds that likely transport deep geothermal water. Models show variable sensitivity to spatial assignment of input variables representing the groundwater source, illustrating the importance of collecting data from discrete groundwater discharges where possible. The methodology shows promise for field experiments designed to assess air–water exchange rates and to quantify total groundwater discharge from a combination of discrete and diffuse sources.

Deep groundwater↗
DOE OSTI2015

Improve earthquake hypocenter using adaptive simulated annealing inversion in regional tectonic, volcano tectonic, and geothermal observation

Observation of earthquakes is routinely used widely in tectonic activity observation, and also in local scale such as volcano tectonic and geothermal activity observation. It is necessary for determining the location of precise hypocenter which the process involves finding a hypocenter location that has minimum error between the observed and the calculated travel times. When solving this nonlinear inverse problem, simulated annealing inversion method can be applied to such global optimization problems, which the convergence of its solution is independent of the initial model. In this study, we developed own program codeby applying adaptive simulated annealing inversion in Matlab environment. We applied this method to determine earthquake hypocenter using several data cases which are regional tectonic, volcano tectonic, and geothermal field. The travel times were calculated using ray tracing shooting method. We then compared its results with the results using Geiger’s method to analyze its reliability. Our results show hypocenter location has smaller RMS error compared to the Geiger’s result that can be statistically associated with better solution. The hypocenter of earthquakes also well correlated with geological structure in the study area. Werecommend using adaptive simulated annealing inversion to relocate hypocenter location in purpose to get precise and accurate earthquake location.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

TReactMech v4.217

TReactMech couples geomechanical processes (poroelasticity, failure, and inelastic strain) with multiphase nonisothermal flow (derived from TOUGH2) and reactive geochemical transport. At its core is the reactive-transport code TOUGHREACT v4.13. TReactMech is an efficient hybrid parallel simulator, solving the geomechanics using finite elements and MPI/PetSc, the multiphase flow using integrated finite difference and MPI/PETSc, and the reactive chemistry using OpenMP. The advantages of TReactMech are in its multiphase flow capabilities (e.g., supercritical CO2, supercritical water, air) and parallel geomechanics including full 3-D stress tensor, shear and tensile failure, coupled to porosity and permeability changes. It is backwardly compatible with TOUGH2 and TOUGHREACT v4.13, allowing for easier transitions between the codes. TReactMech can be used to simulate many natural and engineered subsurface systems, including geothermal reservoirs, borehole heat exchangers, geologic carbon sequestration, geologic storage of nuclear waste, groundwater resources, weathering, sediment diagenesis, seafloor hydrothermal circulation, hydrofracturing in unconventional reservoirs, and injection/production-induced surface deformation.

Sonnenthal, Eric↗
DOE OSTI2013

Geologic Map and GIS Data for the Wabuska Geothermal Area

Wabuska—ESRI geodatabase (ArcGeology v1.3): - Contains all the geologic map data, including faults, contacts, folds, veins, dikes, unit polygons, and attitudes of strata. - List of stratigraphic units and stratigraphic correlation diagram. - One cross‐section.

15 GEOTHERMAL ENERGY↗