SEARCH · Search NASA
Results for “Fracture analysis”
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.
Interpretation of Hydraulic Fractures Based on Microseismic Response in the Marcellus Shale, Monongalia County, West Virginia, USA: Implications for Shale Gas Production
Summary Hydraulic fracturing is critical for extracting shale gas in the subsurface. The treatment technique of multistage hydraulic fracturing is widely used to maximize production. In multistage hydraulic fracturing, not all pumping stages make the same contribution for production, although the designed stimulation process is almost same in every pumping stage. In this study, we characterize the microseismic responses of hydraulic fracturing. Significant variations of microseismic characteristics are observed among the different pumping stages. Pre-existing natural fractures are examined along the horizontal well in selected stages. Combined with the Mohr’s circle analysis, results of the fracture study show that induced hydraulic fractures can be captured by pre-existing natural fractures. Induced hydraulic fractures are simulated by the unconventional fracture model (UFM), and the result reveals that the stress-shadow effect diverts the direction of hydraulic fractures. The diverted hydraulic fracture affects the development of the hydraulic fracture network, which has influence on production. The treatment of two-step pumping is investigated by application of hydraulic diffusivity. Hydraulic fracturing performance and production could be optimized by the treatment of two-step pumping in a single stage. The second pumping creates new fractures and fills the fractures with additional proppants to maintain production for a long duration.
Computed Tomography Scanning and Geophysical Measurements of the Enhanced Geothermal Systems (EGS) Collab SURF Core
The computed tomography (CT) facilities and the Multi-Sensor Core Logger (MSCL) at the National Energy Technology Laboratory (NETL) in Morgantown, West Virginia were used to characterize core from 8 wells from the Enhanced Geothermal Systems (EGS) Collab project located in the Sanford Underground Research Facility’s (SURF) 4850 level. The work in this report falls in Experiment 1 in a three-phase experiment, with the primary purpose to increase understanding of hydraulic fracturing. The primary impetus of this work is a collaboration between U.S. Department of Energy (DOE) and the EGS Collab. The resultant datasets are presented in this report and can be accessed from NETL's Energy Data eXchange (EDX) online system using the following link: https://edx.netl.doe.gov/dataset/surf-collab-core. All equipment and techniques used were non-destructive, enabling future examinations and analyses to be performed on these cores. None of the equipment used was suitable for direct visualization of the pore space in the core, although fractures and discontinuities were detectable with the methods tested. Low resolution CT imagery with the NETL medical CT scanner was performed on the entire core. Qualitative analysis of the medical CT images, coupled with X-ray fluorescence (XRF), P-wave, and magnetic susceptibility measurements from the MSCL were useful in identifying zones of interest for more detailed analysis as well as fractured zones. The ability to quickly identify key areas for more detailed study with higher resolution will save time and resources in future studies. The combination of methods used provides a multi-scale analysis of the core; the resulting macro and micro descriptions are relevant to many subsurface energy related examinations traditionally performed at NETL.
Effect of Shear Displacement and Stress Changes on Fracture Hydraulic Aperture and Flow Anisotropy
Fluid flow in fractures has the potential to drastically change the economic, environmental, and safety risks associated with a subsurface operation. This work focuses on the analysis of experimental shear fracturing and subsequent shear displacement of a shale specimen that contains multiple pre-existing natural fractures. A triaxial direct shear apparatus with concurrent µCT imaging allowed continuous monitoring of sample permeability and changes in fracture geometry during fracture creation, displacement, and changes in stress. Steady-state, low Reynolds number fluid flow simulations were used to determine the changes of each individual fracture’s hydraulic aperture due to displacement and effective stress. The development of anisotropy in the hydraulic aperture of a generated shear fracture was determined for shear displacements of 0.38, 0.74, and 2.04 mm. At high confining stress (30 MPa), shear displacement significantly reduced the hydraulic aperture of the shear fracture to below flow detectable limits and X-ray imaging resolution (23 × 23 × 23 µm 3 ). The shear generated fracture was found to be more sensitive to changes in effective stress compared to the pre-existing natural fractures. Existing models derived from fluid flow equations accounting for stochastic fracture roughness were able to partially fit the simulated flow properties of the pre-existing fractures, but no model was found to fit the flow properties of the newly formed shear fracture. These results are the first report of simultaneous changes in the geometry and flow properties of surrounding, parallel natural fractures as a shear fracture is generated with varying stress and shear displacement conditions.
Fracture mechanisms of sodium silicate glasses
Reactive classical molecular dynamics simulations of sodium silicate glasses, xNa 2 O–(100 – x)SiO 2 (x = 10–30), under quasi-static loading, were performed for the analysis of molecular scale fracture mechanisms. Mechanical properties of the sodium silicate glasses were consistent with experimentally reported values, and the amount of crack propagation varied with reported fracture toughness values. The most crack propagation occurred in NS20 systems (20-mol% Na 2 O) compared with the other simulated compositions. Dissipation via two mechanisms, the first through sodium migration as a lower activation energy process and the second through structural rearrangement as a higher activation energy process, was calculated and accounted for the energy that was not stored elastically or associated with the formation of new fracture surfaces. Here, a correlation between crack propagation and energy dissipation was identified, with systems with higher crack propagation exhibiting less energy dissipation. Sodium silicate glass compositions with lower energy dissipation also exhibited the most sodium movement and structural rearrangement within 10 Å of the crack tip during loading. Therefore, high sodium mobility near the crack tip may enable energy dissipation without requiring formation of structural defects. Therefore, the varying mobilities of the network modifiers near crack tips influence the brittleness and the crack growth rate of modified amorphous oxide systems.
Grain Boundary and Lattice Fracture Toughness of UO 2 Measured Using Small-Scale Mechanics
Small-scale mechanical testing has been used to isolate the properties of individual grain boundaries and regions of the lattice in undoped UO 2 polycrystals. Notched microcantilever specimens were prepared by focused ion beam milling and tested in bending. Fracture toughness values were calculated from finite element analysis of the sample geometry at fracture load using the J-integral method. Lattice and grain boundary fracture toughness values of 1.17 ± 0.16 MPa m 1/2 and 1.08 ± 0.08 MPa m 1/2 , respectively, were obtained. The similarity between the lattice and grain boundary fracture toughness values is consistent with the preference for transgranular room-temperature fracture in UO 2 . The grain boundary and lattice fracture toughness values also agree well with the toughness of dense polycrystals, which suggests that microstructural-scale toughening mechanisms are limited in such samples.
Extended Rice–Thomson analysis and atomistic simulations revealing grain boundary effects on fracture in refractory high-entropy alloys
Significance This work serves to extend the fundamental ductile vs. brittle fracture theory, specifically the Rice–Thomson criterion, by introducing a grain boundary ahead of an initiating crack which propagates at an oblique angle to impinge the boundary. Atomistic fracture simulations on two refractory complex concentrated alloys, the brittle NbMoTaW and the ductile Nb 45 Ta 25 Ti 15 Hf 15 , demonstrate qualitative correspondence with the extended Rice–Thomson criterion and experimental observations. Abstract Understanding how grain boundaries mediate fracture remains a critical challenge in designing ductile, high-performance refractory alloys. Here, we extend the Rice–Thomson criterion to account for the angle between cracks and the impinging grain boundaries (GBs), capturing the competition between intergranular fracture and dislocation-mediated plasticity. Using machine learning interatomic potentials, we performed molecular statics simulations to probe fracture mechanisms in nanocrystalline NbMoTaW and Nb 45 Ta 25 Ti 15 Hf 15 , each with two different grain sizes, revealing trends consistent with experimental observations and the extended Rice model. Comparison with averaged R-curves for bulk samples demonstrates that GBs enhance ductility in Nb 45 Ta 25 Ti 15 Hf 15 in both grain sizes investigated. In contrast, GBs only locally improve fracture resistance in NbMoTaW when cracks are temporarily pinned at GBs inclined at high angles from the crack, but generally promote brittle intergranular fracture. These contrasting behaviors are attributed to differences in GB cohesion, reflecting clear alloying trends that align with ab-initio calculations and trends observed experimentally. Our results bridge classical fracture theory, atomistic simulations, and experimental observations, providing a comprehensive understanding of the fracture mechanisms in nanocrystalline refractory complex concentrated alloys.
Code for Experiment in Publication “Sensitivity Analysis in the Presence of Intrinsic Stochasticity for Discrete Fracture Network Simulations”
Following the Open Research requirements for AGU journals, we must release the code used to perform the experiment described in our recent publication, posted at (https://arxiv.org/abs/2312.04722). This code fits a joint emulator to data from a Discrete Fracture Network (DFN) simulation, performed using the open-source software DFNworks (https://dfnworks.lanl.gov/). All code to be released implements existing methods; there are no novel algorithms nor any major innovations to existing software.
In-situ 4D visualization and analysis of temperature-driven creep in an oil shale propped fracture
Here, the proppant embedment due to creep in shales is a known issue affecting the useable lifetime of wells in unconventional oil and gas recovery. One of the factors influencing creep is the presence of organics, whose properties can be very sensitive to temperature. In this work we investigated for the first time the role of temperature-induced creep increase in proppant embedment in an organics-rich Green River oil shale sample via in-situ synchrotron X-ray micro-computed tomography. We observed that temperatures as low as 75 °C already induce fast creep, with a fracture aperture closing rate of 13 μm/h and a loss of fracture conductivity rate of 8.7%/h, due only to proppant embedment, in the measured interval at the first heating stage. Local displacement data analysis provided evidence for markedly plastic deformation around the proppant-shale contacts, in contrast with the brittle proppant embedment observed on more cemented and less organics-rich shales at room temperature. The results highlight how the problem of temperature-dependent mechanical behavior might be more important than previously thought, in shales with a high content in organics, and that in-situ micro-imaging techniques can play a key role in understanding the underlying mechanisms, contributing to solve creep-related problems associated with hydraulic fracturing in complex scenarios.
A layered 2D computational framework: Theory and applications to nuclear fuel behavior
Nuclear fuel performance computer codes have been developed over the last 50 years to analyze fuel behavior under various operating conditions. Traditionally, these codes used quasi-two-dimensional (also commonly known as 1.5D) representations of the fuel rod, which model the rod using a set of one-dimensional axisymmetric models that represent the behavior at specific axial positions on the rod. Modern fuel performance codes have the ability to investigate full three-dimensional (3D) effects and couple to other physics-based codes for true multiphysics simulations. However, with increasing complexity comes increasing computational costs. Many phenomena of interest involve azimuthally-varying behavior that cannot be represented using the aforementioned quasi-two-dimensional approach, but do not require the use of a full 3D model. To efficiently address these problems, there is a need for a computational framework that provides a compromise between the quasi-two-dimensional and full 3D models. In this paper, we present a new quasi-three-dimensional approach that represents the fuel as a set of 2D planar models that represent the behavior of the fuel cross-section at various axial positions. Presented here are the theory behind the methodology, test cases to illustrate proper implementation, and practical applications of its use in the BISON fuel performance code for a variety of cases in nuclear fuel analysis, such as fuel fracture, axial fuel relocation, and cladding distension and oxidation.
High temperature fracture toughness of single-layer CrAlN and CrAlSiN hard coatings
Hardness and toughness are among of the most important physical properties of hard coatings and it is important to analyze them at temperatures present during coating applications. CrAlN and CrAlSiN are hard coating materials frequently used in high temperature applications. Here, while both possess high oxidation resistance and thermal stability, they differ in their microstructure. In this work, the fracture toughness of magnetron-sputtered single-layer CrAlN and CrAlSiN hard coatings was determined up to 700 °C by in-situ scanning electron microscopy micro-cantilever deflection tests. The results show a decreasing trend of the fracture toughness with increasing temperature. While both coatings' fracture toughness trends are similar up to 400 °C, they diverge at higher temperatures. In general, the evolution of the fracture toughness with increasing temperature shows similarities to the observed changes in hardness in the same temperature range.
Global Sensitivity Analysis of a Reactive Transport Model for Mineral Scale Formation During Hydraulic Fracturing
Injection of water-based hydraulic fracturing fluid (HFF) into tight shale gas/oil formations can increase formation permeability and enhance production rates, but this process frequently causes mineral scale formation that can occlude pore space and hinder flow. To identify the most important factors that control the formation of mineral scales, we applied a novel global sensitivity analysis method—distance-based generalized sensitivity analysis (DGSA)—to a reactive transport model (RTM) that was previously built and calibrated to simulate precipitation of barite [BaSO4] and iron (hydr)oxide [Fe(OH) 3 ] in shale matrices and on fracture surfaces. Reactive transport simulations were run with model parameters randomly sampled based on assigned uncertainties. Modeling results for barite and Fe(OH)3 formation were clustered using machine-learning algorithms. A list of ranked critical input parameters was obtained after statistical quantification of cumulative distribution functions of input parameters. We found that barite formation is most sensitive to the rate of sulfate ion generation, which is determined by the pyrite dissolution rate coefficient and oxidant availability. In addition, barite formation is sensitive to the initial amounts of barite in HFF and shale, followed by barite thermodynamics/kinetics. For Fe(OH) 3 formation, the ranked factors are Fe(OH)3 precipitation rate coefficients, initial HFF pH, initial Fe(OH) 3 amount in HFF, and oxidant availability. Overall, our results provide insights into managing mineral scale formation during hydraulic fracturing to enhance production. Meanwhile, this study serves as an example of global sensitivity analysis of RTMs using the efficient, straightforward, and open-source DGSA method.
Utah FORGE: Seismic DAS and Geophone Borehole Data Processing and 3D Imaging of Vp/Vs Ratio in the 2024 Stimulated Reservoir
This dataset includes a final report and a 3D velocity model derived from seismic DAS and geophone borehole data collected during the April 2024 stimulation of the reservoir at Utah FORGE. The report details the processing of over 50,000 P- and S-wave travel times used in a tomographic inversion to estimate the Vp/Vs ratio, revealing anomalies adjacent to wells 16A and 16B that may be associated with injected fracturing fluids. Additionally, Wadati analysis of more than 27,000 pairs of differential P- and S-wave travel times supports this interpretation with a high Vp/Vs estimate of 1.86, averaged over the dimensions of the earthquake cluster. The accompanying data file provides a 3D model of P- and S-wave velocities and Vp/Vs ratios, structured with spatial coordinates and velocity values.
Revisiting the first Sandia Fracture Challenge with transient deformation heating and strain localization considerations
In this study, we present a systematic experimental study and accompanying theoretical analysis of the dependence of ductile fracture on strain localization, strain hardening rates, deformation-induced thermal softening, and transient heat conduction under spatially uniform as well as spatially heterogeneous deformation. Spatially uniform cases are studied via standard dogbone shaped specimens subject to uniaxial tension. Spatially heterogeneous deformation is studied via the so-called Sandia Fracture Challenge (SFC) specimen, which is a standard compact tension specimen modified with three machined holes in front of a blunt notch. We utilize the same precipitation hardened martensitic stainless steel used in the first SFC experiments, i.e. 15-5 PH with an H1075 heat treatment. We also study 15-5 PH in Condition A and with the H900 heat treatment, each of which has a different hardening behavior and ductility. We find that the ductility does not correlate with the deformation to first crack initiation in the SFC specimen. Instead, hardening rates are better correlated. Moreover, a re-examination finds that an accurate calibration of the hardening rates is strongly correlated with accurate blind predictions of ductile fracture in the SFC specimens. Given that thermal softening can greatly affect hardening rates, we provided a thermographic analysis of deformation-induced heating and transient heat transfer in both the dogbone shaped samples and the SFC specimens. Transient heating in stainless steels is found to have first-order effects on the strain at the onset of necking, strain to fracture, as well as strain localization and ductile fracture in complex geometries even under so-called quasi-static loading rates.
An adaptive space-time phase field formulation for dynamic fracture of brittle shells based on LR NURBS
Here, we present an adaptive space-time phase field formulation for dynamic fracture of brittle shells. Their deformation is characterized by the Kirchhoff–Love thin shell theory using a curvilinear surface description. All kinematical objects are defined on the shell’s mid-plane. The evolution equation for the phase field is determined by the minimization of an energy functional based on Griffith’s theory of brittle fracture. Membrane and bending contributions to the fracture process are modeled separately and a thickness integration is established for the latter. The coupled system consists of two nonlinear fourth-order PDEs and all quantities are defined on an evolving two-dimensional manifold. Since the weak form requires C 1 -continuity, isogeometric shape functions are used. The mesh is adaptively refined based on the phase field using Locally Refinable (LR) NURBS. Time is discretized based on a generalized-α method using adaptive time-stepping, and the discretized coupled system is solved with a monolithic Newton–Raphson scheme. Finally, the interaction between surface deformation and crack evolution is demonstrated by several numerical examples showing dynamic crack propagation and branching.
Beyond the Surface: Non-Invasive Low-Field NMR Analysis of Microbially-Induced Calcium Carbonate Precipitation in Shale Fractures
Microbially-induced calcium carbonate precipitation (MICP) is a biological process in which microbially-produced urease enzymes convert urea and calcium into solid calcium carbonate (CaCO 3 ) deposits. MICP has been demonstrated to reduce permeability in shale fractures under elevated pressures, raising the possibility of applying this technology to enhance shale reservoir storage safety. For this and other applications to become a reality, non-invasive tools are needed to determine how effectively MICP seals shale fractures at subsurface temperatures. In this study, two different MICP strategies were tested on 2.54 cm diameter and 5.08 cm long shale cores with a single fracture at 60 °C. Flow-through, pulsed-flow MICP-treatment was repeatedly applied to Marcellus shale fractures with and without sand (“proppant”) until reaching approximately four orders of magnitude reduction in apparent permeability, while a single application of polymer-based “immersion” MICP-treatment was applied to an Eagle Ford shale fracture with proppant. Low-field nuclear magnetic resonance (LF-NMR) and X-Ray computed microtomography (micro-CT) techniques were used to assess the degree of biomineralization. With the flow-through approach, these tools revealed that while CaCO 3 precipitation occurred throughout the fracture, there was preferential precipitation around proppant. Without proppant, the same approach led to premature sealing at the inlet side of the core. In contrast, immersion MICP-treatment sealed off the fracture edges and showed less mineral precipitation overall. This study highlights the use of LF-NMR relaxometry in characterizing fracture sealing and can help guide NMR logging tools in subsurface remediation efforts.