Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fracturing fluid”

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

Coupled Thermo–Hydro–Mechanical Modeling of Hydro-Shearing Stimulation in an Enhanced Geothermal System in the Raft River Geothermal Field, USA

The Raft River Enhanced Geothermal System (EGS) demonstration project aims to improve the geothermal production by enhancing the reservoir permeability via shear stimulation in well RRG-9. In this paper, we performed a series of 3D thermo– hydro–mechanical (THM) simulations to investigate the influence of hydraulic and thermal effects on the development of this EGS. The model includes synthetic fracture populations based on borehole televiewer images and in-situ stress measurements from well RRG-9. Fracture permeability evolution is determined using an empirical permeability law developed from laboratory experiments. The model was calibrated by comparing the hydraulic response to field observations, including wellhead pressure, injection rate, and well injectivity. Particularly, we analyzed the enhancement of reservoir permeability and the spatial extent of the stimulation zone for the given injection schedule. Our results indicate that the permeability enhancement of fractured geothermal reservoir is caused by the combined effects of injection-induced cooling and fluid pressure increase. The decrease of temperature plays a dominant role in reactivation of natural fractures under the hydro-shearing mechanism, while the higher injection pressure promotes shear failure and enlarges the stimulation zone. Finally, for the specified extensional stress state, the model favors greatest permeability enhancement along the maximum horizontal principal stress, a moderate vertical enhancement, and a smaller gain along the minimum horizontal principal stress.

58 GEOSCIENCES↗
DOE OSTI2021

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

Continuous surface-to-distributed acoustic sensor snapshots explain reactivation of individual natural fractures during an unconventional reservoir stimulation

ABSTRACT Fiber-optic sensing technologies allow petroleum engineering teams to detect hydraulic fracture interaction with boreholes during unconventional reservoir stimulation. In combination with high-repeatability seismic sources, the same distributed acoustic sensors (DASs) enable vertical seismic profiling (VSP) of the fracture evolution away from the boreholes. We discovered clear signatures of seismic scattering on activated fractures during nine days of continuous seismic monitoring of the fracturing stages at the Austin Chalk/Eagle Ford Field Laboratory. The present study applies a novel approach for quantitative analysis of the scattering events in terms of the evolution of the geometry and elastic stiffness of individual fractures. Our characterization strategy sequentially refines the fracture models: from a stack of 1D soft layers to 3D rectangular inclusions. First, we estimate the number of fracture locations and reflectivity using a modified sparse-spike deconvolution of the stacked VSP traces. The fracture set consists of five fractures spaced by 15–30 m with a reflectivity of approximately 1%. Then, we develop a scattering integral method to refine these estimates along with an inversion of the fracture top and bottom for each monitoring vintage. We find that, initially, some of the fractures are located above the monitoring fiber with the height of approximately 100 m. Then we integrate the seismic interpretation with the low-frequency DAS and pressure and microseismic monitoring to reconstruct the activation process of the fractures. Most likely, some of the natural fractures slowly grew downward to the monitoring fiber as a result of fluid injections in the stimulated well. This led to bright strain anomalies but did not trigger seismicity. The top of the fractures remained almost constant and were limited by a lithologic boundary/stress barrier. To our knowledge, this is the first time VSP data enabled tracking of the fracture evolution with such high spatial and temporal resolution, which was previously only available for crosswell surveys and at a much smaller scale.

Glubokovskikh, Stanislav↗
DOE OSTI2023

Influence of extreme fracture flow channels on the thermal performance of open-loop geothermal systems at commercial scale

Abstract Adequate stewardship of geothermal resources requires accurate forecasting of long-term thermal performance. In enhanced geothermal systems and other fracture-dominated reservoirs, predictive models commonly assume constant-aperture fractures, although spatial variations in aperture can greatly affect reservoir permeability, fluid flow distribution, and heat transport. Whereas previous authors have investigated the effects of theoretical random aperture distributions on thermal performance, here we further explore the influence of permeability heterogeneity considering field-constrained aperture distributions from a meso-scale field site in northern New York, USA. Using numerical models of coupled fluid flow and heat transport, we conduct thermal–hydraulic simulations for a hypothetical reservoir consisting of a relatively impervious porous matrix and a single, horizontal fracture. Our results indicate that in highly channelized fields, most well design configurations and operating conditions result in extreme rates of thermal drawdown (e.g., 50% drop in production well temperatures in under 2 years). However, some other scenarios that account for the risks of short-circuiting can potentially enhance heat extraction when mass flow rate is not excessively high, and the direction of geothermal extraction is not aligned with the most permeable features in the reservoir. Through a parametric approach, we illustrate that well separation distance and relative positioning play a major role in the long-term performance of highly channelized fields, and both can be used to help mitigate premature thermal breakthrough.

Rangel-Jurado, Nicolás↗
DOE OSTI2024

Revealing complex subsurface dynamics with continuous seismic monitoring: Observations using distributed acoustic sensing and surface orbital vibrators during hydraulic fracturing

Understanding hydraulic fracturing is crucial to improving the stimulation of unconventional reservoirs and increasing fluid production. This study develops a novel seismic monitoring technology using distributed acoustic sensing (DAS) and surface orbital vibrators (SOV) to capture fracture seismic response and mechanical properties at high temporal intervals. We analyze continuous time-lapse vertical seismic profiling (VSP) data acquired every hour during the first nine days of treatment of an unconventional reservoir in the Austin Chalk/Eagle Field Laboratory. The VSP data contain clear seismic signals scattered from the activated fractures. The spatiotemporal changes of the fracture reflectivity revealed by the SOV/DAS data correlate well with the observations of fracture locations inferred from low-frequency DAS data. These results capture the fracture opening and closure processes, as well as highlight potential prestage activations of the fractures due to hydraulic connectivity with preexisting fracture systems. Therefore, analysis of the presented data set provides a unique opportunity to understand fracture initiation and subsequent evolution, not only in the context of unconventional resources but also in enhanced geothermal systems.

Correa, Julia↗
DOE OSTI2026

Detecting and Characterizing Fracture Zones Using a Convolutional Neural Network

This project directly supports the Geothermal Technologies Office (GTO) objectives outlined in the Multi-Year Program Plan (MYPP) by advancing two key research areas: “Exploration and Characterization” and “Data, Modeling, and Analysis.” This project has successfully demonstrated a pre-drilling ability to image and characterize the distribution and connectivity of subsurface faults and fractures, key parameters for identifying permeable pathways that enable geothermal fluids to circulate and produce energy. Specifically, we developed and implemented innovative machine learning methodologies to enhance geothermal exploration. Large-scale faults were detected using a Convolutional Neural Network (CNN), while small-scale fractures were characterized using a novel Double-Beam Neural Network (DBNN). These tools have proven both technically effective and cost-efficient by reducing reliance on expensive exploratory drilling. Through collaboration with our geothermal industry partner, this research has significantly advanced techniques for identifying hidden geothermal systems and extending the productive lifespan of existing geothermal fields. We applied our methods to two geothermal fields—Soda Lake (Nevada) and Lightning Dock (New Mexico)—to identify shallow steam-charged fracture zones and characterize deep faults at depths of 1.5-2 km. The steam zone identified at the Soda Lake geothermal field showed excellent agreement with prior drilling data, validating the effectiveness of our approaches. In addition, the analysis revealed three new prospective drilling targets for further development and verification. The outcomes of this project improve our scientific understanding of geothermal reservoir behavior, enhance exploration efficiency, extend the economic life of existing geothermal plants. Ultimately, these advancements contribute to GTO’s goal of achieving more sustainable, affordable, and data-driven geothermal energy development across the United States.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

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.

02 PETROLEUM↗
DOE OSTI2021

Accessing the Subsurface Biosphere Within Rocks Undergoing Active Low-Temperature Serpentinization in the Samail Ophiolite (Oman Drilling Project)

The Oman Drilling Project established an “Active Alteration” multi-borehole observatory in peridotites undergoing low-temperature serpentinization in the Samail Ophiolite. The highly serpentinized rocks are in contact with strongly reducing fluids. Distinct hydrological regimes, governed by differences in rock porosity and fracture density, give rise to steep redox (Eh +200 to –750 mV) and pH (pH range 8.5–11.2) gradients within the 300–400 m deep boreholes. The serpentinites and fluids host an active subsurface ecosystem. Microbial cell abundances in serpentinite vary at least six orders of magnitude, from ≤3.5 × 10 1 to 2.9 × 10 7 cells/g. Low levels of biological sulfate reduction (2–1,000 fmol/cm 3 /day) can be detected in rock cores, particularly in rocks in contact with reduced groundwaters with pH < 10.5. Thermodesulfovibrio is the predominant sulfate reducer identified via metagenomic sequencing of adjacent groundwater communities. Here, we infer that transport and reaction of microbially generated sulfide with the serpentine and brucite assemblages gives rise to optical darkening and sulfide overprinting, including the formation of tochilinite-vallerite group minerals, potentially serving as an indicator that this system is inhabited by microbial life. Olivine mesh-cores replaced with ferroan brucite and minor awaruite, abundant veins containing hydroandradite garnet and polyhedral serpentine, and late-stage carbonate veins are suggested as targets for future spatially resolved life-detection investigations. The high-quality whole-round core samples that have been preserved can be further probed to define how life distributes itself and functions within a system where chemical disequilibria are sustained by low-temperature water/rock interaction, and how biosignatures of in situ microbial activity are generated.

58 GEOSCIENCES↗
DOE OSTI2022

Long-term permeability evolution of shale seal rocks with argon and scCO 2

In this work, we conducted a suite of experiments to evaluate the long-term permeability evolution of shale rocks under constant effective stress, before, during, and after interaction with supercritical carbon dioxide (scCO 2 ). To do so, we measured the time-dependent evolution of argon permeabilities for the pre- and post-reacted samples (i.e., before and after long-term interaction with scCO 2 ). In addition, we obtained permeability evolution during long-term interaction with scCO 2 . The samples showed either relatively constant permeabilities or a moderate decrease during pre reaction long-term argon tests. The permeability evolution during long-term CO 2 tests showed continuous increase, continuous decrease, or cycles of increase/decrease in permeability. The long term response of the samples to CO 2 included phenomena such as 20 (i) salt precipitation, (ii) swelling-induced cracks, and (iii) carbonate dissolution. While it is obvious that salt precipitation and swelling-induced cracks decrease and increase the permeability, respectively, the sample response to carbonate dissolution proved to be more complex and may increase or decrease the permeability. The permeability evolution during post-reaction long-term argon injection is also affected by the contribution of each of these three phenomena, during long-term interaction with CO 2 . We observe increase, decrease, and constant permeability evolution during post-reaction argon tests. Our experiments reveal that the initial permeability of the samples plays a significant role on the long-term permeability response of shales in the presence argon and scCO 2 fluids. This study shows, when shales are hydraulically fractured with CO 2 , their initial permeability has a more significant role than their permeability evolution over time.

58 GEOSCIENCES↗
DOE OSTI2022

Coupling Subsurface and Above-Surface Models for Optimizing the Design of Borefields and District Heating and Cooling Systems

Accurate dynamic energy simulation is important for the design and sizing of district heating and cooling systems with geothermal heat exchange for seasonal energy storage. Current modeling approaches in building and district energy simulation tools typically consider heat conduction through the ground between boreholes without flowing groundwater. While detailed simulation tools for subsurface heat and mass transfer exist, these fall short in simulating above-surface energy systems. To support the design and operation of such systems, the study developed a coupled model including a software package for building and district energy simulation, and software for detailed heat and mass transfer in the subsurface. For the first, it uses the open-source Modelica Buildings Library, which includes dynamic simulation models for building and district energy and control systems. For the heat and mass transfer in the soil, it uses the TOUGH simulator. The TOUGH family of codes can model heat and multi-phase, multi-component mass transport for a variety of fluid systems, as well as chemical reactions, in fractured porous media. The study validated the coupled modeling approach by comparing the simulation results with one from the g-function based ground response model. It then looked into effects when the water table and the regional groundwater flow are considered in the ground, from the perspective of heat exchange between borehole and ground, and the electrical consumption of the district heating and cooling systems. To access the simulation models, please find the links in the submission: -- For coupled approach validation: see model Buildings.Fluid.Geothermal.Borefields.Examples.BorefieldsWithTough and Buildings.Examples.DistrictReservoirNetworks.Examples.Reservoir3Variable_TOUGH from the "Modelica Building Library" resource, branch issue1495_tough_interface, commit a2667c0. -- For the study of the effect of water table: see model Buildings.Examples.DistrictReservoirNetworks.Examples.Reservoir3Variable_TOUGH from he "Modelica Building Library" resource, branch issue1495_tough_interface_moreIO, commit 760de49. -- For the study of the effect of regional groundwater flow: see Buildings.Examples.DistrictReservoirNetworks.Examples.Reservoir3Variable_TOUGH from he "Modelica Building Library" resource, branch issue1495_tough_interface_moreIO_3D, commit c2a2d2a. The coupling interface script "GrounResponse.py" can be found from the above links in the folder Buildings/Resources/Python-Sources. Also, the needed files for TOUGH simulation are in the folder Buildings/Resources/Python-Sources/ToughFiles that can be accessed through the above links. A brief description of these files is given below; detailed specifications for the first three files may be found in the TOUGH3 Users Guide (Jung et al., 2018) https://tough.lbl.gov/documentation/tough-manuals/. (1) INCON - initial conditions for each grid block (2) INFILE - main input file with material properties and control parameters (3) MESH - description of the computational grid (4) readsave - Modelica/TOUGH interface program: read the final output of TOUGH simulation after TOUGH time step and prepare for transfer to Modelica for next Modelica time step (5) readsave.inp - input parameters for program readsave (6) writeincon - Modelica/TOUGH interface program: write the output of Modelica after Modelica time step and prepare for transfer to TOUGH as initial conditions for the next TOUGH step (7) writeincon.inp - input parameters for program writeincon

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Deep Learning–Assisted Multiobjective Optimization of Geological CO 2 Storage Performance under Geomechanical Risks

In geological CO 2 storage, designing the optimal well control strategy for CO 2 injection to maximize CO 2 storage while minimizing the associated geomechanical risks is not trivial. This challenge arises due to pressure buildup, CO 2 plume migration, the highly nonlinear nature of geomechanical responses to rock-fluid interaction, and the high computational cost associated with coupled flow and geomechanics simulations. In this paper, we introduce a novel optimization framework to address these challenges. The optimization problem is formulated as follows: maximize total CO 2 storage while minimizing geomechanical risks by adjusting the injection schedules within bounded constraints. The geomechanical risks are primarily driven by injection-induced pressure build-up, which is characterized by ground displacement and the induced microseismicity. We used the Fourier neural operator (FNO)-based deep learning model to construct surrogate models, replacing the time-consuming coupled flow and geomechanics simulations for evaluating the aforementioned objective functions. The developed surrogate models have been incorporated into a multiobjective optimization framework through a genetic algorithm to reduce the computational burden. The proposed optimization framework reduces the computational cost from approximately 2,400 hours, when using objective function evaluations based on physics-based simulations, to around 20 minutes. A set of Pareto-optimal solutions of the proposed workflow yields nontrivial optimal decisions, reducing the microseismicity potential and the vertical displacement. This Pareto front highlights the optimal trade-offs between CO 2 storage amount, safety, and ground displacement, emphasizing the need for careful optimization and management of injection strategies to achieve a balanced outcome. The novelty of this work is twofold. First, we demonstrate the importance of incorporating the minimization of the geomechanical risks as objective functions into the CO 2 storage optimization workflow to mitigate the potential risk of induced microseismicity and ground displacement. Second, we leverage the FNO-based surrogate models to optimize a real-field CO 2 storage operation.

42 ENGINEERING↗
DOE OSTI2025

High-Resolution Simulations of Geological CO 2 Injection: Application to the SPE11 Benchmark

Geological carbon sequestration (GCS) will play a critical role in decarbonization and in facilitating the transition to clean energy systems. Because CO 2 is highly mobile, ensuring its safe and permanent injection into subsurface geological formations involves monitoring over larger spatial domains and longer time periods than is typical for hydrocarbon reservoirs. This can benefit from simulation tools capable of modeling key CO 2 trapping mechanisms, particularly those optimized for speed and scalability on high-performance computing systems. Using isothermal versions of the SPE11B and SPE11C benchmark cases, we conduct a mesh refinement study simulating CO 2 injection into kilometer-scale rock formations at centimeter resolution with the GEOS open-source simulation framework. We focus on how mesh refinement improves the accuracy of convective mixing in both 2D and 3D simulations. The computational costs associated with achieving a converged solution highlight the need for predictive upscaling techniques. A systematic performance scaling analysis—including both central processing unit (CPU) and graphics processing unit (GPU) architectures—complements the “Results” section.

Geosciences↗
DOE OSTI2021

Recent Thermal and Hydraulic Stimulation Results at Raft River, ID EGS Site

The Raft River geothermal system, located 160 km northwest of Salt Lake City, Utah is the site of a Department of Energy Enhanced Geothermal System (EGS) stimulation demonstration project. The target well RRG-9 ST-1, is roughly 1.6 km (1 mile) south of the power plant. This well has undergone a series of thermal and hydraulic stimulation treatments. Furthermore, the well passes through approximately 1,524 m MD (5,000 ft. MD) of laterally discontinuous Quaternary and Tertiary volcanoclastic and volcanic rocks that unconformably overlie the Precambrian basement. The Elba Quartzite formation located in the Precambrian basement is the target zone for the stimulation program. Fluid injected into the well exits through a fracture zone in the Elba Quartzite at 1,567 m TVD (5,140 ft. TVD). In June 2013 low-rate stimulation began at RRG-9 ST-1 using injection water from the geothermal plant. The average injection rate, temperature, and well head pressure of the low-rate stimulation was 163 Lpm (43 gpm), 34°C, and 1,931kPa (280 psig), respectively. This injection continued through late August of 2013 when wellhead pressures were intentionally increased to approximately 5,516 kPa (800 psig) at approximately 984 Lpm (260 gpm). From September 12 to September 24, 2013 cold well water was injected. The average surface injection temperature was approximately 12°C. On September 25, 2013 injection of the warmer plant water was resumed. This was continued through the rest of 2013 and into 2014. Between April 1, 2014 and April 4, 2014, a higher rate hydraulic stimulation was conducted. Colder water (~10°C) was injected at rates of 20 bpm to 30 bpm (3,178 Lpm to 4,770 Lpm) at wellhead pressures of 5,861 kPa (850 psig) over approximately 24 hours. Microseismic activity was recorded near the Tertiary-Precambrian contact. Subsequently, injection of water from the plant was resumed. As of mid-January, 2015, over 666 million liters (176 million gallons) have been injected into the well. The fate of this water is the subject of ongoing numerical analyses and continuous site monitoring. Various analyses have been implemented to monitor and assess the ongoing stimulation program. Modified Hall plot analysis shows increasing conductivity and/or a decreasing skin factor around the wellbore. Since April, 2014, the injectivity index has increased significantly from 0.28 Lpm/kPa (0.5 gpm/psig) to 0.93 Lpm/kPa (1.7 gpm/psig).

enhanced geothermal system↗
DOE OSTI2023

Corrosive Influence of Carbon Dioxide on Crack Initiation in Quartz: Comparison With Liquid Water and Vacuum Environments

The stimulation of crack growth in quartz and siliceous materials by injecting carbon dioxide (CO 2 ) represents a key technology in long-term carbon storage and in the development of natural gas wells. While this technology is widely used, the molecular impact of CO 2 interactions on the solid matrix is only incompletely understood. In this work, we employ reactive molecular dynamics simulations to study how the CO 2 fluid environment affects the mechanical properties of pre-cracked single-crystal quartz. The thermodynamic conditions of interest are those relevant to subsurface reservoirs. We report how structural properties of quartz—bond length distribution and crack tip shape—evolve upon introduction of a fluid. These properties are directly related to macroscopic quantities of the global stress-strain curves, thus reaffirming the inherent coupling across multiple scales for fluid-solid interactions in the subsurface. Here, we find that CO 2 reduces the fracture toughness of quartz by 12.1% compared to that of quartz in vacuum, thereby promoting crack growth and enhancing fluid transport in the subsurface.

58 GEOSCIENCES↗
DOE OSTI2021

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

Isotopic fractionation as in-situ sensor of subsurface reactive flow and precursor for rock failure

Greater utilization of subsurface reservoirs perturbs in-situ chemical-mechanical conditions with wide ranging consequences from decreased performance to project failure. Understanding the chemical precursors to rock deformation is critical to reducing the risks of these activities. To address this need, we investigated the coupled flow-dissolution- precipitation-adsorption reactions involving calcite and environmentally-relevant solid phases. Experimentally, we quantified (1) stable isotope fractionation processes for strontium during calcite nucleation and growth, and during reactive fluid flow; (2) consolidation behavior of calcite assemblages in the common brines. Numerically, we quantified water weakening of calcite using molecular dynamics simulations; and quantified the impact of calcite dissolution rate on macroscopic fracturing using finite element models. With microfluidic experiments and modeling, we show the effect of local flow fields on the dissolution kinetics of calcite. Taken together across a wide range of scales and methods, our studies allow us to separate the effects of reaction, flow, and transport, on calcite fracturing and the evolution of strontium isotopic signatures in the reactive fluids.

42 ENGINEERING↗
DOE OSTI2026

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.

materials science↗