Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FRACTURES”

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

Geochemical Impact of Acid Spearhead and Slickwater Stimulation on Wolfcamp Shale from the Hydraulic Fracturing Test Site

The Hydraulic Fracturing Test Site 1 (HFTS-1) was a field study performed in the Wolfcamp Formation in the West Texas Permian (Midland) Basin, USA, with a focus on improving the efficiency of hydraulic fracturing. Investigating site-specific rock-fluid geochemical interactions during hydraulic fracturing is an important step to understanding the impact on formation shale porosity, permeability, and long-term shale gas production. During field operations in this region, hydraulic fracturing fluid (HFF) injection usually starts with a concentrated acid spearhead for rapid rock dissolution, followed by the injection of near-neutral pH slickwater containing chemicals and proppants. A multistep sequential injection approach was used to investigate different stages of rock-fluid interactions. The carbonate content in the host rock is important when acid spearhead is considered, as carbonate mineral dissolution is rapid and can result in porosity and permeability changes in the shale matrix. Here, in this study, we designed flow-through experiments using fractured carbonate-rich and clay-rich Wolfcamp shale cores with (1) a short-time acid soaking step and (2) a long-term slickwater flow-through step to simulate the injection method used at HFTS-1. The fluid chemistry was analyzed. A thorough mineralogical progression [e.g., Calcium (Ca) dissolution and iron (Fe) redox progression] in the cores during HFF injection was also characterized and imaged by synchrotron microprobe. Reactive transport modeling was performed based on the experimental setup. The results showed that the acid spearhead is a crucial step in creating a reaction front by mineral dissolution, especially in carbonate-rich shales. A slight layer of ferrihydrite precipitated during the slickwater flow-through period. This study provides insights into potential geochemical impact due to hydraulic fracturing operations in the Permian Basin.

clastic rock↗
DOE OSTI2013

Fracture-Permeability Relationships of Utica Shale: Early Results

In October 2013, we will begin a Los Alamos Lab-funded project on the hydraulic fracturing of shale. The goal is to use unique LANL microfluidic and triaxial core flood experiments integrated with state-of-the-art numerical simulation to reveal fundamental dynamics of fracture-fluid interactions. Our particular focus is on the development of CO 2 -based fracturing fluids and fracturing techniques to enhance production, greatly reduce wastewater, while simultaneously sequestering CO 2 . Our more general goal is to determine fundamental fracture aperture-permeability relationships and the efficiency of hydrocarbon extraction as a function of fracture characteristics.

02 PETROLEUM↗
DOE OSTI2015

Anisotropic Elastic-Waveform Modeling for Fracture Characterization in EGS Reservoirs

Enhanced geothermal systems (EGS) contain newly created fractures in addition to possible existing fractures. Accurate characterization and monitoring of EGS reservoirs are crucial for optimal placement of new wells and effective extraction of geothermal heat. The fractured reservoirs behave as anisotropic media where seismic waves propagate with different velocities along different directions. In addition, the anisotropic properties of fluid-filled fracture zones could be different from those of dry fracture zones. We develop an optimized rotated staggered-grid elastic-wave finite-difference method for simulating seismic-wave propagation in heterogeneous, anisotropic media. Our new method uses a few extra grid points and optimized finite-difference coefficients based on the space-time dispersion relation, and reduce numerical dispersion of the conventional rotated staggered-grid finite-difference scheme. We validate our new method using synthetic vertical-seismic-profiling (VSP) data for an anisotropic geophysical model built with geologic features found at the Raft River EGS reservoir. This improved and optimized rotated staggered-grid finite-difference method provides an essential tool for analyzing VSP data, reverse-time migration, and elastic-waveform inversion in anisotropic, fractured reservoirs.

58 GEOSCIENCES↗
DOE OSTI2015

Multi-scale modeling of inter-granular fracture in UO2

A hierarchical multi-scale approach is pursued in this work to investigate the influence of porosity, pore and grain size on the intergranular brittle fracture in UO2. In this approach, molecular dynamics simulations are performed to obtain the fracture properties for different grain boundary types. A phase-field model is then utilized to perform intergranular fracture simulations of representative microstructures with different porosities, pore and grain sizes. In these simulations the grain boundary fracture properties obtained from molecular dynamics simulations are used. The responses from the phase-field fracture simulations are then fitted with a stress-based brittle fracture model usable at the engineering scale. This approach encapsulates three different length and time scales, and allows the development of microstructurally informed engineering scale model from properties evaluated at the atomistic scale.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2014

Tracer Methods for Characterizing Fracture Creation in Engineered Geothermal Systems

The aim of this proposal is to develop, through novel high-temperature-tracing approaches, three technologies for characterizing fracture creation within Engineered Geothermal Systems (EGS). The objective of a first task is to identify, develop and demonstrate adsorbing tracers for characterizing interwell reservoir-rock surface areas and fracture spacing. The objective of a second task is to develop and demonstrate a methodology for measuring fracture surface areas adjacent to single wells. The objective of a third task is to design, fabricate and test an instrument that makes use of tracers for measuring fluid flow between newly created fractures and wellbores. In one method of deployment, it will be used to identify qualitatively which fractures were activated during a hydraulic stimulation experiment. In a second method of deployment, it will serve to measure quantitatively the rate of fluid flowing from one or more activated fracture during a production test following a hydraulic stimulation.

15 GEOTHERMAL ENERGY↗
DOE OSTI2019

Use of Neutron Flux Calculated by Shift in a Grizzly Reactor Pressure Vessel Fracture Simulation

The reactor pressure vessel (RPV) plays the critical role of containing the reactor and coolant in a light water reactor nuclear power plant, and must be able to safely perform this function under a variety of normal and off-normal transient loading conditions. The RPV is subjected to harsh environmental conditions (elevated temperature and high neutron flux) whenever the plant is operating, which leads to embrittlement of the steel over time. RPVs contain populations of flaws introduced during the manufacturing process, and the primary safety concern is that during a transient event, a fracture could initiate at the site of one of these flaws and lead to rupture of the RPV. RPVs are very rugged structures, and the probability of this occurring is very low, but as the RPV's steel becomes embrittled over time, its resistance to fracture decreases and this becomes more likely. Assessing the likelihood of RPV fracture during a transient event involves performing a probabilistic fracture mechanics analysis, in which random sampling procedures are used to assess the probability of fracture of a population of flaws. Grizzly offers unique capabilities for modeling the RPV using 1D, 2D, or 3D representations, taking advantage of parallel computing for this type of analysis. One of the important inputs to Grizzly is the spatial distribution of the neutron fluence, which has a large impact on the material embrittlement, and should thus be characterized as accurately as possible. The VERA core simulator developed by the CASL program takes advantage of high performance computing resources to perform high-fidelity multiphysics simulations of reactor physics, including the effects of neutron transport, fuel performance, and thermal hydraulics. Recently a capability was developed to also evaluate neutron fluxes away from the reactor core through the Shift code. This allows for a detailed 3D map of the spatial variation of the fluence in components such as the RPV to be computed, which is ideal for use with Grizzly. This report documents a procedure developed to use the fluences computed in Shift in a Grizzly simulation, which allows for improved realisim in the prediction of the effects of environmental exposure on propensity for reactor pressure vessel fracture. This capability demonstrated on RPVs with multiple types of configurations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Fracture Network Quantification during CO2 Injection

This is the presentation prepared for the ARMA 2025 (59th US Rock Mechanics/Geomechanics Symposium) Conference held in Santa Fe, New Mexico, June 8-11, 2025. Accurate mapping and quantification of these networks are essential to ensure the integrity of CO2 storage reservoirs, understand and reduce potential leakage, and maintain long-term environmental safety. This study presents a novel machine learning-driven approach, integrated with geomechanical analysis, to quantify fracture networks and assess their spatial distribution during CO2 injection. This paper combines microseismic monitoring data with principles of hydraulic diffusivity and geomechanical analysis to characterize reservoir scale fracture network. The novelty of our approach lies in its capacity to assimilate time-dependent pressure data and microseismicity into a cohesive framework, which not only identifies microseismic triggering fronts but also tracks fracture distribution during active injection. Besides, leveraging image log data and analysis our approach also provides another angle of the insights to solidate the fracture networks understanding and geomechanical impacts. Key results from our study include the detection of over 100 distinct fracture clusters across the injection site, with fracture orientations strongly correlated with the prevailing in-situ stress field.

CO2 storage and sequestration↗
DOE OSTI2024

A multiphysics coupling framework for exascale simulation of fracture evolution in subsurface energy applications

Predicting the evolution of fractured media is challenging due to coupled thermal, hydrological, chemical and mechanical processes that occur over a broad range of spatial scales, from the microscopic pore scale to field scale. We present a software framework and scientific workflow that couples the pore scale flow and reactive transport simulator Chombo-Crunch with the field scale geomechanics solver in GEOS to simulate fracture evolution in subsurface fluid-rock systems. This new multiphysics coupling capability comprises several novel features. An HDF5 data schema for coupling fracture positions between the two codes is employed and leverages the coarse resolution of the GEOS mechanics solver which limits the size of data coupled, and is, thus, not taxed by data resulting from the high resolution pore scale Chombo-Crunch solver. The coupling framework requires tracking of both before and after coarse nodal positions in GEOS as well as the resolved embedded boundary in Chombo-Crunch. We accomplished this by developing an approach to geometry generation that tracks the fracture interface between the two different methodologies. The GEOS quadrilateral mesh is converted to triangles which are organized into bins and an accessible tree structure; the nodes are then mapped to the Chombo representation using a continuous signed distance function that determines locations inside, on and outside of the fracture boundary. The GEOS positions are retained in memory on the Chombo-Crunch side of the coupling. The time stepping cadence for coupled multiphysics processes of flow, transport, reactions and mechanics is stable and demonstrates temporal reach to experimental time scales. The approach is validated by demonstration of 9 days of simulated time of a core flood experiment with fracture aperture evolution due to invasion of carbonated brine in wellbore-cement and sandstone. We also demonstrate usage of exascale computing resources by simulating a high resolution version of the validation problem on OLCF Frontier.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Simulating Fracture Sealing by Granular LCM Particles in Geothermal Drilling

Lost circulation occurs when the returned fluid is less than what is pumped into the well due to loss of fluid to pores or fractures. A lost-circulation event is a common occurrence in a geothermal well. Typical geothermal reservoirs are often under-pressured and have larger fracture apertures. A severe lost-circulation event is costly and may lead to stuck pipe, well instability, and well abandonment. One typical treatment is adding lost-circulation materials (LCMs) to seal fractures. Conventional LCMs fail to properly seal fractures because their mechanical limit is exceeded at elevated temperatures. In this paper, parametric studies in numerical simulations are conducted to better understand different thermal effects on the sealing mechanisms of LCMs. The computational fluid dynamics (CFDs) and the discrete element method (DEM) are coupled to accurately capture the true physics of sealing by granular materials. Due to computational limits, the traditional Eulerian–Eulerian approach treats solid particles as a group of continuum matter. With the advance of modern computational power, particle bridging is achievable with DEM to track individual particles by modeling their interactive forces between each other. Particle–fluid interactions can be modeled by coupling CFD algorithms. Fracture sealing capability is investigated by studying the effect of four individual properties including fluid viscosity, particle size, friction coefficient, and Young’s modulus. It is found that thermally degraded properties lead to inefficient fracture sealing.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

A New Modeling Framework for Multi-Scale Simulation of Hydraulic Fracturing and Production from Unconventional Reservoirs

This paper describes a new modeling framework for microscopic to reservoir-scale simulations of hydraulic fracturing and production. The approach builds upon a fusion of two existing high-performance simulators for reservoir-scale behavior: the GEOS code for hydromechanical evolution during stimulation and the TOUGH+ code for multi-phase flow during production. The reservoir-scale simulations are informed by experimental and modeling studies at the laboratory scale to incorporate important micro-scale mechanical processes and chemical reactions occurring within the fractures, the shale matrix, and at the fracture-fluid interfaces. These processes include, among others, changes in stimulated fracture permeability as a result of proppant behavior rearrangement or embedment, or mineral scale precipitation within pores and microfractures, at µm to cm scales. In our new modeling framework, such micro-scale testing and modeling provides upscaled hydromechanical parameters for the reservoir scale models. We are currently testing the new modeling framework using field data and core samples from the Hydraulic Fracturing Field Test (HFTS), a recent field-based joint research experiment with intense monitoring of hydraulic fracturing and shale production in the Wolfcamp Formation in the Permian Basin (USA). Below, we present our approach coupling the reservoir simulators GEOS and TOUGH+ informed by upscaled parameters from micro-scale experiments and modeling. We provide a brief overview of the HFTS and the available field data, and then discuss the ongoing application of our new workflow to the HFTS data set.

58 GEOSCIENCES↗
DOE OSTI2013

Intergranular fracture in UO2: derivation of traction-separation law from atomistic simulations

In this study, the intergranular fracture behavior of UO2 was studied by molecular dynamics simulations using the Basak potential. In addition, the constitutive traction-separation law was derived from atomistic data using the cohesive-zone model. In the simulations a bicrystal model with the (100) symmetric tilt E5 grain boundaries was utilized. Uniaxial tension along the grain boundary normal was applied to simulate Mode-I fracture. The fracture was observed to propagate along the grain boundary by micro-pore nucleation and coalescence, giving an overall intergranular fracture behavior. Phase transformations from the Fluorite to the Rutile and Scrutinyite phases were identified at the propagating crack tips. These new phases are metastable and they transformed back to the Fluorite phase at the wake of crack tips as the local stress concentration was relieved by complete cracking. Such transient behavior observed at atomistic scale was found to substantially increase the energy release rate for fracture. Insertion of Xe gas into the initial notch showed minor effect on the overall fracture behavior.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

Using well casing as an electrical source to monitor hydraulic fracture fluid injection

The depth to surface resistivity (DSR) method transmits current from a source located in a cased or openhole well to a distant surface return electrode while electric field measurements are made at the surface over the target of interest. This paper presents both numerical modelling results and measured data from a hydraulic fracturing field test where conductive water was injected into a resistive shale reservoir during a hydraulic fracturing operation. Modelling experiments show that anomalies due to hydraulic fracturing are small but measureable with highly sensitive sensor technology. The field measurements confirm the model results,showing that measured differences in the surface fields due to hydraulic fracturing have been detected above the noise floor. Our results show that the DSR method is sensitive to the injection of frac fluids; they are detectable above the noise floor in a commercially active hydraulic fracturing operation, and therefore this method can be used for monitoring fracture fluid movement.

02 PETROLEUM↗
DOE OSTI2016

Numerical Study of Critical Role of Rock Heterogeneity in Hydraulic Fracture Propagation

Log and seismic data indicate that most shale formations have strong heterogeneity. Conventional analytical and semi-analytical fracture models are not enough to simulate the complex fracture propagation in these highly heterogeneous formation. Without considering the intrinsic heterogeneity, predicted morphology of hydraulic fracture may be biased and misleading in optimizing the completion strategy. In this paper, a fully coupling fluid flow and geomechanics hydraulic fracture simulator based on dual-lattice Discrete Element Method (DEM) is used to predict the hydraulic fracture propagation in heterogeneous reservoir. The heterogeneity of rock is simulated by assigning different material force constant and critical strain to different particles and is adjusted by conditioning to the measured data and observed geological features. Based on proposed model, the effects of heterogeneity at different scale on micromechanical behavior and induced macroscopic fractures are examined. From the numerical results, the microcrack will be more inclined to form at the grain weaker interface. The conventional simulator with homogeneous assumption is not applicable for highly heterogeneous shale formation.

Dual Lattice Discrete Element Method↗
DOE OSTI2018

The distribution, orientation, and characteristics of natural fractures for Experiment 1 of the EGS Collab Project, Sanford Underground Research Facility

The EGS Collab project is focused on understanding and predicting permeability enhancement and evolution in crystalline rocks. To accomplish this, the project is creating a suite of intermediate-scale (~10-20 m) field test beds coupled with stimulation and interwell flow tests that will provide a basis to better understand the fracture geometries and processes that control heat transfer between rock and stimulated fractures. As part of the site characterization effort for the first experimental test bed, our team has worked on mapping the distribution, orientation, and nature of open and healed fractures exposed along the drift wall and within the eight bore holes drilled for this test bed. The fractures have been characterized through detailed description of continuous cores obtained from these boreholes, evaluation of televiewer logs, and mapping of fractures and seeps exposed along the drift wall. The fracture data are being compiled and interpreted for slip and dilation tendencies, and will be incorporated into coupled-process geomechanical flow and transport models to better constrain the planned flow and tracer tests.

15 GEOTHERMAL ENERGY↗
DOE OSTI2013

Intergranular fracture in UO{sub 2}: derivation of traction-separation law from atomistic simulations

In this study, the intergranular fracture behavior of UO{sub 2} was studied by molecular dynamics simulations using the Basak potential. In addition, the constitutive traction-separation law was derived from atomistic data using the cohesive-zone model. In the simulations a bicrystal model with the (100) symmetric tilt Σ5 grain boundaries was utilized. Uniaxial tension along the grain boundary normal was applied to simulate Mode-I fracture. The fracture was observed to propagate along the grain boundary by micro-pore nucleation and coalescence, giving an overall intergranular fracture behavior. Phase transformations from the Fluorite to the Rutile and Scrutinyite phases were identified at the propagating crack tips. These new phases are metastable and they transformed back to the Fluorite phase at the wake of crack tips as the local stress concentration was relieved by complete cracking. Such transient behavior observed at atomistic scale was found to substantially increase the energy release rate for fracture. Insertion of Xe gas into the initial notch showed minor effect on the overall fracture behavior. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Fracture Network Quantification during CO2 Injection

This is the conference paper accompanying an oral presentation at the ARMA 2025 (59th US Rock Mechanics/Geomechanics Symposium) Conference held in Santa Fe, New Mexico, June 8-11, 2025. Accurate mapping and quantification of these networks are essential to ensure the integrity of CO2 storage reservoirs, understand and reduce potential leakage, and maintain long-term environmental safety. This study presents a novel machine learning-driven approach, integrated with geomechanical analysis, to quantify fracture networks and assess their spatial distribution during CO2 injection. This paper combines microseismic monitoring data with principles of hydraulic diffusivity and geomechanical analysis to characterize reservoir scale fracture network. The novelty of our approach lies in its capacity to assimilate time-dependent pressure data and microseismicity into a cohesive framework, which not only identifies microseismic triggering fronts but also tracks fracture distribution during active injection. Besides, leveraging image log data and analysis our approach also provides another angle of the insights to solidate the fracture networks understanding and geomechanical impacts. Key results from our study include the detection of over 100 distinct fracture clusters across the injection site, with fracture orientations strongly correlated with the prevailing in-situ stress field.

CO2 storage and sequestration↗
DOE OSTI2015

Shale gas and non-aqueous fracturing fluids: Opportunities and challenges for supercritical CO₂

Hydraulic fracturing of shale formations in the United States has led to a domestic energy boom. Currently, water is the only fracturing fluid regularly used in commercial shale oil and gas production. Industry and researchers are interested in non-aqueous working fluids due to their potential to increase production, reduce water requirements, and to minimize environmental impacts. Using a combination of new experimental and modeling data at multiple scales, we analyze the benefits and drawbacks of using CO₂ as a working fluid for shale gas production. We theorize and outline potential advantages of CO₂ including enhanced fracturing and fracture propagation, reduction of flow-blocking mechanisms, increased desorption of methane adsorbed in organic-rich parts of the shale, and a reduction or elimination of the deep re-injection of flow-back water that has been linked to induced seismicity and other environmental concerns. We also examine likely disadvantages including costs and safety issues associated with handling large volumes of supercritical CO₂. The advantages could have a significant impact over time leading to substantially increased gas production. In addition, if CO₂ proves to be an effective fracturing fluid, then shale gas formations could become a major utilization option for carbon sequestration.

03 NATURAL GAS↗
DOE OSTI2019

EGS Collab Experiment 1: Common Discrete Fracture Network

This package includes data and models that support hydraulic fracture stimulation and fluid circulation experiments in the Sanford Underground Research Facility (SURF). A paper by Schwering et al. (2020) describes the deterministic basis for developing a "common" discrete fracture network (CDFN) model of significant natural fractures in EGS Collab Testbed 1 on the 4850-Level of SURF. The ReadMe for this model shows drift, wells, scanlines, fracture data, interpreted fractures, and geophysical visualizations. There is also a summary of the data that was used in this experiment and includes results from reviewing core, televiewer (TV) logs, core-TV depth/feature registration, and from mapping weeps in the 4850-Level drift. The CDFN is intended to be a baseline model of the pre-stimulated testbed (though some observations from stimulation helped inform the model).

15 GEOTHERMAL ENERGY↗