Search NASA⌕ Search

SEARCH · Search NASA

Results for “fracture geometry”

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

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

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Results of Multiple Tracer Injections into Fractures in the EGS Collab Testbed-1

The EGS Collab project constructed an intermediate scale (~10-20 m) testbed at the 4850 level of the Stanford Underground Research Facility (SURF) in South Dakota for testing and validating fracture stimulation and flow/transport models. This testbed consists of eight ~200 ft (~60 m) HQ-diameter (9.6 cm) boreholes that are drilled into the crystalline rocks of the Poorman Formation from the West Access Drift tunnel. Of the eight boreholes, one borehole is used as an injection/stimulation well, while another sub-parallel borehole located about 10 m away from the injection well is used as a production well, and rest of the other boreholes are used as geophysical/fluid sampling monitoring wells. Hydraulic stimulation activities were conducted at three locations along the injection hole in an attempt to create direct fracture connections to the production hole. A flow system has been established between injection and production boreholes through a set of hydraulically stimulated fractures propagated from a notch located at 164 ft in the injection hole. Although we planned for a single production well, a set of natural fractures in the testbed is believed to have intersected the stimulated hydraulic fractures and provided additional flow paths for water to be transported to the drift through multiple monitoring boreholes and weep zones. As the flow tests continued after stimulation activities, single or a combination of two or three producers become dominant producers at different times, mostly as a response to intersection of hydraulic fracture and testbed wells, activation of natural fractures and making new leak points to the monitoring wells, resealing of leaky wells, and so on. Since late October 2018, multiple tracers were injected into the fracture system at the 164 ft location that involves both stimulated and natural fractures, and tracers were recovered from multiple locations in nearby wells and weep. The cumulative water recoveries over the time have ranged from 50 to 90%; however, the injected tracer recoveries from various production sources are much less (ranging from a few percentages to 38%). Changes in water/tracer recoveries, shifting of major producing wells from one well to the others, and other observations (e.g., microseismic, electric resistivity tomography, etc.) indicate a testbed that has undergone several changes since October 2018. In this paper, we present tracer recovery data accumulated during several tracer campaigns, provide conceptual testbed flow pathways, and provide observations that suggest a evolutionary nature of fracture volume and fracture geometry in the testbed.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

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.

Jin, Ge↗
DOE OSTI2022

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.

58 GEOSCIENCES↗
DOE OSTI2021

Observations and Modeling of Fiber-Optics Strain on Hydraulic Fracture Height Growth in HFTS-2

Understanding fracture height growth can be of great significance to optimizing field development and improving recovery. The Hydraulic Fracturing Test Site 2 (HFTS-2) has provided a unique opportunity and an advanced dataset to allow us to observe and understand fracture geometries rigorously. Low frequency distributed acoustic sensing (LF-DAS) data from a vertical well in HFTS-2 showed three key observations: (i) excessive upward height growth (>1000 ft) and limited downward growth of the hydraulic fractures during pumping, (ii) considerable additional upward fracture height growth (~300 ft) after well shut in, and (iii) very complex LF-DAS strain rate patterns for a small fiber-to-stage offset. Advanced geomechanical modeling was performed to simulate the hydraulic fracture propagation and the resulting strain responses in the vertical direction. The modeling results demonstrated asymmetric upward and downward fracture height growths as observed in HFTS-2 with a similar upward height growth rate. Simulated waterfall plots of vertical strain rate showed distinct patterns for different fiber-to-fracture distances. The upward-growing fracture tip can be clearly identified by the interfaces between compressing and extending zones. It was also found that the complex strain rate patterns observed in HFTS-2 for small fiber-stage offsets were not caused by the mechanical layering but possibly result from the simultaneous propagation of multiple hydraulic fractures at different rates. Furthermore, the simulation results improved the understating of the HFTS-2 LF-DAS data and the simulated strain rate patterns could also serve as templates for fracture height interpretation from LF-DAS data in future.

58 GEOSCIENCES↗
DOE OSTI2021

Microseismic at HFTS2: A Story of Three Stimulated Wells

A microseismic data set was acquired at Hydraulic Fracturing Test Site 2 (HFTS2) in the Delaware Basin (Texas, USA) by the HFTS2 consortium as a part of data-collection program designed to help maximize the efficiency of hydraulic fracturing, understand the impacts of nearby depleted zones, and increase the production of hydrocarbons. Advanced methods of microseismic monitoring applied at HFTS2 include 1) native synchronization of five monitor arrays, 2) real-time location of microseismic events, 3) joint event hypocenter/ velocity-model inversion that locates microseismic events simultaneously with constructing azimuthally anisotropic velocity models, 4) high-resolution reservoir imaging, utilizing microseismic events as downhole sources. HFTS2 microseismic data set illustrated the influence of geologic complexity, variations in treatment designs, and the presence of parent wells on the efficiency of hydraulic fracturing. Three Wolfcamp wells, Boxwood 1H, 2H, and 4H, were stimulated sequentially, exhibiting very different patterns of microseismicity. When the Boxwood 1H well is stimulated, its toe stages quickly reacted to the presence of depleted Bitterroot section to the east of the Boxwood pad, extending microseismic events over 3,000 ft to the east. This behavior changes dramatically after the treatment passes the Bitterroot section at the toe, and the heel stages in 1H showed tight microseismic grouping around the stimulated well, indicative of treatment of virgin reservoir rock. The Boxwood 2H demonstrates this same pattern of behavior between the toe and heel stages of the well, suggesting that the depleted Bitterroot section remains easier to refracture than adjacent non-depleted rock. Microseismic activity associated with Boxwood 4H treatment is distinctly different from those in Boxwood wells 1H and 2H. The events generated during stimulation of 4H are hardly influenced by the depleted Bitterroot section and fill the space between the previously stimulated 1H and 2H wells. Altogether, microseismic patterns suggest that treatment fluids tend to find the weakest passages through geologic formations, and the relative well positions and treatment sequence have a significant impact on fracture geometry, as indicated by microseismic event locations.

58 GEOSCIENCES↗
DOE OSTI2024

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.

58 GEOSCIENCES↗
DOE OSTI2022

Observations and Modeling of Fiber Optic Strain on Hydraulic Fracture Height Growth in Hydraulic Fracturing Test Site 2 (HFTS-2)

Summary Understanding fracture height growth can be of great significance to optimizing field development and improving recovery. The Hydraulic Fracturing Test Site 2 (HFTS-2) has provided a unique opportunity and an advanced data set to allow us to observe and understand fracture geometries rigorously. Low-frequency distributed acoustic sensing (LF-DAS) data from a vertical well in HFTS-2 showed three key observations: (i) excessive upward height growth (>1,000 ft) and limited downward growth of the hydraulic fractures during pumping, (ii) considerable additional upward fracture height growth (~300 ft) after well shut in, and (iii) very complex LF-DAS strain rate patterns for a small fiber-to-stage offset. Advanced geomechanical modeling was performed to simulate the hydraulic fracture propagation and the resulting strain responses in the vertical direction. The modeling results demonstrated asymmetric upward and downward fracture height growths as observed in HFTS-2 with a similar upward height growth rate. Simulated waterfall plots of vertical strain rate showed distinct patterns for different fiber-to-fracture distances. The upward-growing fracture tip can be clearly identified by the interfaces between compressing and extending zones. It was also found that the complex strain rate patterns observed in HFTS-2 for small fiber-stage offsets were not caused by the mechanical layering but possibly result from the simultaneous propagation of multiple hydraulic fractures at different rates. The simulation results improved the understating of the HFTS-2 LF-DAS data, and the simulated strain rate patterns could also serve as templates for fracture height interpretation from LF-DAS data in future.

Engineering↗
DOE OSTI2023

HFTS-1 (Final Report)

Despite the long history of hydraulic fracturing, the optimal number of fracturing stages during multi-stage fracture stimulation in horizontal wells can be a complex subject. In addition to the increased expense of multistage fracturing in horizontal wells, increasing the number of fracturing stages does not always correlate with a commensurate increase in production. The problem is the application of a uniform fracture stimulation design to all stages with no consideration for geological variations along the wellbore. The result is an inefficient use and costly waste of energy and water. Optimization of the fracturing process requires an understanding of the cause-and-effect relationship between fracturing parameters and local geological properties at a given location along the wellbore. Realizing that the generalized rock mechanics theories and hypotheses are not truly applicable to fractured and laminated shales, quantifiable impacts of a shale’s geomechanical and depositional features are a prerequisite for design and implementation of optimized hydraulic fractures. The overarching goal of this project is to understand and define the relationships of shale geology and fracture dynamics using detailed field data that includes coring of the fracture domain. Finally, resource recovery from shale formations is estimated to be less than fifteen percent. Research proposed in this project will establish the foundations for investigating enhanced recovery techniques for increased resource recovery in existing fracture treated wells. Natural gas as an EOR fluid is in the initial stages of being used broadly in the oil and gas industry and as such, many aspects of the process need to be researched and addressed prior to widespread acceptance. In addition, re-fracturing techniques, particularly using cemented liners also need to be further investigated. This project explores several of these EOR techniques in greater detail, not just at the theoretical level, but through actual field test sites. This report delves into the research program, identified tasks, several studies, observations, results, and conclusions under the Hydraulic Fracturing Test Site (HFTS) project. This test project is subdivided into three distinct phases and each phase is discussed in detail under separate sections. The original test site for Phase 1 as well as the Phase 2 site are in Midland Basin in Texas and operated by Laredo Petroleum. In the first two phases, the program targeted the Wolfcamp formation in Midland Basin. Phase 3 of the study targeted the EagleFord play in South Texas with Devon Energy as the operator. Significant conclusions from this study include: 1) lateral and vertical fracture growth mapping as well as proppant transport behavior, 2) Use of field gas and its applicability for huff n puff EOR injection in the Midland Basin, 3) Development and validation of proppant log, both at core scale as well as in-fill well implementation by analyzing drilling mud return samples, 4) Applicability of liner refracturing including best practices for enhanced unit level productivity. This report will detail all these studies as well as technical results and conclusions from said studies. These include general subsurface characterization work, completion designs as implemented, studies looking at fracture geometry, high resolution microseismic study, advanced crosswell strain mapping studies including drainage characterization, core characterization results, proppant analysis (core and observation wells), geochemistry, etc. to name a few.

02 PETROLEUM↗
DOE OSTI2025

Surrogate models for development of unconventional shale reservoirs by an integrated numerical approach of hydraulic fracturing, flow and geomechanics, and machine learning

We develop well-completion surrogate models by taking an integrated workflow of hydraulic fracturing, flow, geomechanics, and machine learning simulation. There are three steps in the proposed workflow. First, history-matching processes are conducted with the field data including pumping and production data for characterization. Second, full-physics simulation is performed with various parameters of the field development (e.g., cluster spacing, clusters per stage, pumping rates and times, amount of proppant, and well spacing) to generate multiple simulation results by changing the parameters of the completion design with well-known hydraulic fracturing, reservoir, geomechanics simulators to calculate fracture geometry, reservoir depressurization, induced stress changes. The workflow is demonstrated over a field in the Southern Midland Basin. Here, we take two completion scenarios: a single well case followed by a multi-well case. Finally, a Long Short-Term Memory (LSTM) machine learning algorithm is employed to create surrogate models that can replicate the full-physics simulation results. Furthermore, results show that the trained models applied in the single well and multi-well cases for a particular geological system can provide good accuracy close to those provided by full-physics simulations. Specifically, the site-specific surrogate models can predict fracture parameters (length, height, and surface area) and cumulative production accurately with computational efficiency, suggesting our proposed workflow can be used as a pragmatic tool for expediting the well completion optimization process.

Geomechanics↗
DOE OSTI2020

Three-dimensional fracture continuum characterization aided by surface time-domain electromagnetics and hydrogeophysical joint inversion—proof-of-concept

Efficient and safe production of hydraulically fractured reservoirs benefits from the prediction of their geometrical attributes. Geophysical methods have the potential to provide data that are sensitive to fracture geometries, alleviating the typically sparse nature of in situ reservoir observations. Moreover, surface-based methods can be logistically and economically attractive since they avoid operational interference with the injection well infrastructure. Here, we investigate the potential of the surface-based time-domain electromagnetic (EM) method. EM methods can play an important role owing to their sensitivity to injection-induced fluid property changes. Two other advantageous factors are the EM signal-enhancing effect of vertical steel-cased wells and the fact that injected proppants can be enhanced to produce a stronger electrical conductivity contrast with the reservoir’s connate fluid. Nevertheless, an optimal fracture characterization will no doubt require the integration of EM and reservoir injection and production data. We hence carry out our investigations within a hydrogeophysical parameter estimation framework where EM data and injection flow rates are combined in a fully coupled way. Given the interdisciplinary nature of coupled hydrogeophysical inverse modeling, we dedicate one section to laying out key aspects in a didactic manner.

58 GEOSCIENCES↗
DOE OSTI2025

Joint state-parameter estimation for the reduced fracture model via the united filter

Here, in this paper, we introduce an effective United Filter method for jointly estimating the solution state and physical parameters in flow and transport problems within fractured porous media. Fluid flow and transport in fractured porous media are critical in subsurface hydrology, geophysics, and reservoir geomechanics. Reduced fracture models, which represent fractures as lower-dimensional interfaces, enable efficient multi-scale simulations. However, reduced fracture models also face accuracy challenges due to modeling errors and uncertainties in physical parameters such as permeability and fracture geometry. To address these challenges, we propose a United Filter method, which integrates the Ensemble Score Filter (EnSF) for state estimation with the Direct Filter for parameter estimation. EnSF, based on a score-based diffusion model framework, produces ensemble representations of the state distribution without deep learning. Meanwhile, the Direct Filter, a recursive Bayesian inference method, estimates parameters directly from state observations. The United Filter combines these methods iteratively: EnSF estimates are used to refine parameter values, which are then fed back to improve state estimation. Numerical experiments demonstrate that the United Filter method surpasses the state-of-the-art Augmented Ensemble Kalman Filter, delivering more accurate state and parameter estimation for reduced fracture models. This framework also provides a robust and efficient solution for PDE-constrained inverse problems with uncertainties and sparse observations.

Bayesian inference↗
DOE OSTI2020

Application of Koopman operator for model-based control of fracture propagation and proppant transport in hydraulic fracturing operation

This work explores the application of the recently developed Koopman operator approach for model identification and feedback control of a hydraulic fracturing process. Controlling fracture propagation and proppant transport with precision is a challenge due in large part to the difficulty of constructing approximate models that accurately capture the characteristic moving boundary and highly-coupled dynamics exhibited by the process. Koopman operator theory is particularly attractive here as it offers a way to explicitly construct linear representations for even highly nonlinear dynamics. The method is data-driven and relies on lifting the states to an infinite-dimensional space of functions called observables where the dynamics are governed by a linear Koopman operator. Here this work considers two problems: (a) fracture geometry control, and (b) proppant concentration control. In both cases, an approximate linear model of the corresponding dynamics is constructed and used to design a model predictive controller (MPC). The manuscript shows that in the case of highly nonlinear dynamics, as observed in the proppant concentration, use of canonical functions in the observable basis fails. In such cases, a priori system knowledge can be leveraged to choose the required basis. The numerical experiments demonstrate that the Koopman linear model shows excellent agreement with the real system and successfully achieves the desired target values maximizing the oil and gas productivity. Additionally, due to its linear structure, the Koopman models allow convex MPC formulations that avoid any issues associated with nonlinear optimization.

42 ENGINEERING↗
DOE OSTI2024

Integrating Experiments and Well Logs to Predict Caney Shale Static Mechanical Properties During Production with Supervised Machine Learning

Caney shale is one of the emerging oil reservoirs in Oklahoma. Understanding the impact of effective stress on its mechanical properties is critical for predicting hydraulic fracture geometry and overall hydrocarbon production. The objective of our study is to evaluate the impact of effective stress on the dynamic Young’s modulus using ultrasonic velocity measurements for Caney shale samples. A triaxial cell was utilized to measure ultrasonic (P-wave and S-wave) velocities for ten downhole Caney shale samples under various effective stresses to indirectly assess the impact of pore pressure change. The dynamic Young’s moduli estimated from these measurements were integrated with available conventional well logs (excluding sonic logs) and triaxial test results from Benge et al. (2021) to predict the static Young’s modulus using Random Forest (RF) and Extreme Gradient Boosting (XGBoost) models. The results showed that the estimated dynamic Young’s moduli from ultrasonic measurements were higher than the corresponding static Young’s modulus of cores from the same vertical well at similar depths. With increasing effective stress, the dynamic Young’s modulus increased for all samples. The estimated dynamic-to-static correction factor tended to be higher in zones of high neutron porosity (PHIN) and low density compared to other zones. Finally, SHapley Additive exPlanations (SHAP) for RF and XGBoost models identified depth, gamma ray (GR), and PHIN as key features for predicting the static Young’s modulus. This study enhances our understanding of the dynamic and static Young’s moduli for the Caney shale interval, as a function of effective stress and conventional well logs. The findings from this study can improve predictions of production throughout the well's lifespan by offering insights into the mechanical property degradation resulting from pore pressure depletion.

Kholy, Sherif M.↗
DOE OSTI2022

HFTS-2 (Final Report)

This report delves into the research program, identified tasks, various studies, observations, results, and conclusions under the Hydraulic Fracturing Test Site (HFTS-2) project. This test site is located in Loving County close to the Texas, New Mexico border. The program targeted the Wolfcamp formation in Delaware Basin. The aim of the project was to improve our understanding of hydraulic fracturing processes to not only optimize unconventional resource development in US Basins but to minimize environmental footprint as well. Various new diagnostic technologies were demonstrated at this site, including distributed acoustic microseismic, high resolution strain using Rayleigh frequency band (DSS-RFS), and in-fill proppant logging tool. Significant conclusions from this study include: 1) lateral and vertical fracture growth mapping as well as proppant transport behavior, 2) Use of DSS-RFS as a diagnostic tool to understand cluster level behavior, 3) Validation of proppant log, both at core scale as well as in-fill well implementation, and 4) New modeling tools to map far-field strain response using distributed strain sensing data. Specifically, far-field fracture distribution is influenced by linear fracture corridors defined by the stress-state as well as clustering behavior influenced by spatio-temporal effects. Depletion impacts on new stimulation was demonstrated and technically quantified. Significant vertical fracture height growth was observed. However, our understanding is limited by available gauge data availability. Numerical strain modeling techniques have been developed and demonstrated to help understand observed strain response behavior. Finally, proppant log has been further validated as a useful tool to map propped zones within the SRV. It has also shown very strong correlation with observed drawdown behavior at this test site. This report will detail many of these studies as well as technical results and conclusions from said studies. These include general subsurface characterization work, completion designs as implemented, studies looking at fracture geometry, high resolution microseismic study including source mechanisms, advanced DSS-RFS studies including drainage characterization, core characterization results, proppant analysis (core and in-fill child well), geochemistry, etc.

02 PETROLEUM↗
DOE OSTI2024

Phase Field Modeling of Near-Wellbore Hydraulic Fracture Nucleation and Propagation

The goal of Task 2 of the project titled "Closing the loop between in situ stress complexity and EGS fracture complexity" is to model fracture nucleation and propagation in the near wellbore region and, in combination with the experimental study conducted as part of Task 4, investigate the relationship between the in situ stress and the near-wellbore fracture patterns. To this end, We have identified phase field as a promising modeling approach to model the nucleation and propagation process in the near wellbore region. In phase field, fractures are represented by a diffuse variable, i.e the damage. Compared to approaches that model fractures as sharp interfaces, phase field does not require any re-meshing or element insertion algorithm, which makes the treatment of complex fracture geometries straightforward, even with simple meshes.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Closing the Loop Between In Situ Stress Complexity and EGS Fracture Complexity (Quarterly Research Performance Progress Report: Q3Y3)

As part of Task 1, we have started by testing our modeling capabilities by reproducing isothermal DFIT simulations presented in the literature. To start, we have been working towards adjusting numerical parameters to reproduce DFIT tests with different fracture geometries. Then, we have focused modeling the DFIT conducted at well 58-32 (cycle 4 in zone 2). Specifically, we have carefully analyzed the influence of the aperture-contact stress relationship and investigated how two different laws affect model results. We have also introduced thermal effects in the model and rerun the simulations to estimate the effect of cooling on the DFIT response. Due to the small amount of fluid injected during a DFIT, thermal effects seem to be moderate.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Characterizing Stress Roughness at Utah FORGE Through Simulation of Hydraulic Fracture Growth During 16A Stage 3 Stimulation

In most geologic formations, the vertical gradient of the minimum horizontal stress (𝜎 h ) exceeds the hydrostatic gradient, often driving hydraulic fractures to propagate upward. However, observations indicate that the upward growth of real hydraulic fractures is less pronounced than theoretical predictions based on smoothly varying 𝜎 h fields. In fact, the layered structure of sedimentary rocks hinders fracture propagation across layers, partially explaining the limited height growth observed in practice (Zoback et al., 2022). While crystalline rocks lack the pervasive layered fabric of sedimentary formations, they still possess structural fabric formed over their geologic history, resulting in inherently "rough" in situ stress fields. Recent studies have identified another key factor: the "roughness" of in situ stress, characterized by temporal, relatively short-wavelength fluctuations superimposed on the overall stress gradient. This stress roughness leads to apparent toughness anisotropy, where the vertical toughness appears significantly larger than the horizontal toughness (P. Fu et al., 2019; Dontsov & Suarez-Rivera, 2021). Neglecting the effects of rock fabric and stress roughness can yield inaccurate predictions of hydraulic fracture geometry and growth rates at larger length scales.

58 GEOSCIENCES↗
DOE OSTI2025

Closing the Loop Between In Situ Stress Complexity and EGS Fracture Complexity: Quarterly Research Performance Progress Report

As part of Task 1, we have started by testing our modeling capabilities by reproducing isothermal DFIT simulations presented in the literature. To start, we have been working towards adjusting numerical parameters to reproduce DFIT tests with different fracture geometries. Then, we have focused modeling the DFIT conducted at well 58-32 (cycle 4 in zone 2). Specifically, we have carefully analyzed the influence of the aperture-contact stress relationship and investigated how two different laws affect model results.

58 GEOSCIENCES↗