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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Post-forming Room Temperature Brittle Fracture in a High-Strength Low-Alloy Steel Sheet After Various Forming Modes

The effects of various sheet forming modes on the post-forming room temperature fracture behavior of a 550-MPa yield strength high-strength low-alloy (HSLA) steel are investigated. In-plane biaxial stretching, plane strain (cold rolling), uniaxial tension, cylindrical cup drawing, and in-plane compression (IPC) are examined up to a von Mises effective prestrain ε{sub eff} = 0.7. Most of the results pertain to sub-size Charpy-type impact specimens prepared from the prestrained material. After sufficiently large prestrains, the fracture behavior is highly anisotropic. Low-energy fracture modes, namely cleavage and intergranular fracture, occur after all modes of prestraining along plane directions that are perpendicular to the principal compressive prestrain. These planes are parallel to the sheet surface after biaxial stretching and cold rolling and correspond to brittle splits extending into the primary, ductile fracture surface. After cup drawing, the brittle planes are oriented perpendicular to the circumferential compressive prestrain. This results in very low energy fractures propagating in the length direction of a fully drawn cup. After in-plane compression, brittle fracture occurs along planes that are perpendicular to the compression direction. The fracture toughness for a crack propagating along such a plane after IPC to ε{sub eff} = 0.38 was very low (12.7 MPa √m), much lower than in the undeformed condition (247 MPa√m). Changes in grain shape and in crystallographic texture caused by the various prestraining modes, as well as the microstructural damage at non-metallic inclusions and carbides, are examined to try to understand the fracture behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Evolution of planar fractures in limestone: The role of flow rate, mineral heterogeneity and local transport processes

In this study, the processes that affect the evolution of fracture geometry and reactive transport in limestone fractures were evaluated using an integrated experimental and modeling approach. Flow-through experiments were performed to investigate the alteration of three artificially made planar limestone fractures, which were exposed to HCl-acidified water flow (pH 3.8 ± 0.1) at three different rates. Reactive infiltration instability developed in the fractures because of small fracture aperture variations and naturally occurring heterogeneity in the rock, i.e., the presence of clay spots of higher porosity. The dissolution regime changed from conical to dominant and ramified wormhole with increasing flow rate, i.e., Péclet number. Maps of fracture aperture obtained from X-ray micro-tomography imaging, together with measurements of permeability and analyses of fluid chemistry, provided experimental constraints for modeling. Numerical simulations were performed using a 2.5D continuum reactive transport model to examine the development of flow patterns and the reactivity and permeability evolution over time, and to interrogate the impacts of fluid reactivity, mineral heterogeneity and diffusive boundary layer. The modeling results showed that fracture alteration is very sensitive to initial fluid reactivity and mineralogical heterogeneity, in addition to flow rate. Finally, the changes of local transport processes during fracture opening also affected the reactivity of the fractured samples.

58 GEOSCIENCES↗
DOE OSTI2024

Relative permeabilities for two-phase flow through wellbore cement fractures

Multiple fluids are likely to exist in fractures and flow paths associated with leaky wellbores, including liquids (e.g., crude oil) and gases (e.g., gas exsolved from liquid). These fluids occupy and move through different portions of the pore spaces within the fractures depending on many factors, including fluid properties, fracture size, and the amount of the different fluids. Upward leakage of any phase, through the fracture, can contaminate water-bearing formations, create hazardous surface conditions, and compromise the functionality of the wellbore. Early signs of wellbore leaks may be expressed by anomalous pressure behavior at surface monitoring points on cavern storage wells. These pressure anomalies are difficult to interpret, necessitating knowledge of the factors that affect the multiphase flow in fractures and porous media. These parameters are critical to modeling multiphase flow in fractures. This insight can guide further diagnosis and maximize leak remediation. Here, our study focuses on the relationship of the liquid–gas relative permeabilities for representative variable-aperture wellbore cement fracture. To obtain the relative permeability of each phase, two-phase flow tests were conducted where both fluids were flowing simultaneously through a fractured wellbore cement specimen under a range of factors, namely (1) aperture size, (2) capillary numbers, and (3) viscosity ratio. The flow experiments were conducted under a range of confining stresses and flow velocities, using nitrogen gas and silicone oils (of different viscosities) in a specially designed pressure vessel. The sum of gas and oil relative permeabilities were found to be less than one under all conditions, which indicates that the presence of one phase affects the permeability of the other phase, and vice versa. Since the gas phase flow conditions include a significant inertial flow component in addition to viscous flow, the inertial flow coefficients at different saturation states are presented. The factors affecting the relationship between the relative permeabilities are discussed in detail. A new mathematical model for estimating the relative permeability of wellbore cement fracture is presented and experimentally validated.

58 GEOSCIENCES↗
DOE OSTI2014

The development of in situ fracture toughness evaluation techniques in hydrogen environment

Reliability of hydrogen pipelines and storage tanks is significantly influenced by the mechanical performance of the structural materials exposed in the hydrogen environment. Fracture behavior and fracture toughness are of specific interest since they are relevant to many catastrophic failures. However, many conventional fracture testing techniques are difficult to be realized under the presence of hydrogen. Thus it is desired to develop novel in situ techniques to study the fracture behavior of structural materials in hydrogen environments. In this study, special testing apparatus were designed to facilitate in situ fracture testing in H 2 . A torsional fixture was developed to utilize an emerging fracture testing technique, Spiral Notch Torsion Test (SNTT). The design concepts will be discussed. Preliminary in situ testing results indicated that the exposure to H 2 significantly reduces the fracture toughness of 4340 high strength steels by up to 50 percent. Furthermore, SNTT tests conducted in air demonstrated a significant fracture toughness reduction in samples subject to simulated welding heat treatment using Gleeble, which illustrated the effect of welding on the fracture toughness of this material.

08 HYDROGEN↗
DOE OSTI2012

Vertebral Compression Fracture (VCF) After Spine Stereotactic Body Radiation Therapy (SBRT): Analysis of Predictive Factors

Purpose: Vertebral compression fractures (VCFs) are increasingly observed after spine stereotactic body radiation therapy (SBRT). The aim of this study was to determine the risk of VCF after spine SBRT and identify clinical and dosimetric factors predictive for VCF. The analysis incorporated the recently described Spinal Instability Neoplastic Score (SINS) criteria. Methods and Materials: The primary endpoint of this study was the development of a de novo VCF (ie, new endplate fracture or collapse deformity) or fracture progression based on an existing fracture at the site of treatment after SBRT. We retrospectively scored 167 spinal segments in 90 patients treated with spine SBRT according to each of the 6 SINS criteria. We also evaluated the presence of paraspinal extension, prior radiation, various dosimetric parameters including dose per fraction ({>=}20 Gy vs <20 Gy), age, and histology. Results: The median follow-up was 7.4 months. We identified 19 fractures (11%): 12 de novo fractures (63%) and 7 cases of fracture progression (37%). The mean time to fracture after SBRT was 3.3 months (range, 0.5-21.6 months). The 1-year fracture-free probability was 87.3%. Multivariate analysis confirmed that alignment (P=.0003), lytic lesions (P=.007), lung (P=.03) and hepatocellular (P<.0001) primary histologies, and dose per fraction of 20 Gy or greater (P=.004) were significant predictors of VCF. Conclusions: The presence of kyphotic/scoliotic deformity and the presence of lytic tumor were the only predictive factors of VCF based on the original 6 SINS criteria. We also report that patients with lung and hepatocellular tumors and treatment with SBRT of 20 Gy or greater in a single fraction are at a higher risk of VCF.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2022

Experimentally validated multiphysics modeling of fracture induced by thermal shocks in sintered UO 2 pellets

Uranium Dioxide (UO 2 ) fuel powers almost all commercial Nuclear Power Plants (NPPs) worldwide, generating carbon-free energy and contributing to the fight against climate change. UO 2 fuel incurs damage and fractures due to large thermal gradients that develop across the fuel pellet during normal and transient operating conditions. A comprehensive understanding of the underlying mechanisms by which these processes take place is still lacking. A combined experimental and computational approach is utilized here to quantify the behavior of UO 2 fuel fracture induced by thermal shock. Here, this work introduces both (1) an experimental study to understand the fuel fracturing behavior of sintered UO 2 pellets when exposed to thermal shock, and (2) a Multiphysics phase-field fracture model capable of simulating this process. Parametric studies were conducted to evaluate the effects of uncertainties in fracture properties on the fracture behavior of UO 2 due to thermal shocking. A set of energy release rate (or equivalently fracture toughness) and contract area (the part of the fuel pellet in direct contact with the cold bath) were able to capture the overall fracture trends of the corresponding experimental data. Our combined approach presents a new method for accounting for the effects of microstructure and sample size on the energy release rate/fracture toughness. The experimental data were collected from multiple experiments that exposed UO 2 pellets to high-temperature conditions (589–676 °C) followed by a quench in sub-zero water. This work demonstrates that joint experimental and computational efforts are able to advance the understanding of thermal fracture in the primary fuel source for existing and future NPPs.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Multiscale spatial analysis of fracture arrangement and pattern reconstruction using Ripley's K-function

This work presents novel multiscale spatial data analytics using Ripley's K-function, as a measure of spatial interaction, to study one-dimensional arrangement of fractures. Fracture spatial arrangements are classified into clustered, anticlustered, or indistinguishable from random by testing statistical significance of the calculated Ripley's K-function. Characterizations of fracture arrangements are performed as a function of length scale and position. Analysis of the K-function along the study interval identifies where fracture clustering and anticlustering occur. Here, a simulation technique is also introduced here to statistically reconstruct spatial arrangements and to generate fracture realizations that are spatially similar to the fractures observed in the field. With this simulation technique, one can also fill spatial gaps in fracture measurements where data are absent, unreliable, or unused. Synthetic as well as field-measured 1D fracture datasets are used for testing and demonstration. Methods introduced in this work can be readily applied to fracture datasets observed in outcrops, borehole image logs, and cores.

58 GEOSCIENCES↗
DOE OSTI2013

Fracture modes in tubular LSFCO ceramic membranes under graded reducing conditions

Highlights: ► Microstructural evolution in LSFCO membranes under graded environment is reported. ► The role of chemically induced stresses and oxygen deficiency is evaluated. ► The stress distribution is modeled by a point defect model. - Abstract: Chromium (III) oxide (Cr{sub 2}O{sub 3})-doped LaSrFeO{sub 3} perovskite, La{sub 0.2}Sr{sub 0.8}Fe{sub 0.8}Cr{sub 0.2}O{sub 3−δ} (LSFCO), is being considered as a potential material for applications in solid oxide fuel cells, gas separation membranes, and electrochemical reactors because of its high electro-catalytic activity. Similar to other perovskites, the performance and mechanical strength of LSFCO materials are significantly affected by environment and temperature. Here, we report a fracture gradient phenomenon in tubular C-ring-shaped LSFCO ceramic membranes under graded reducing conditions. The graded reducing condition was produced by flushing N{sub 2} on the outer side of the C-ring membranes at 1000 °C while keeping the inner side untreated. The rings were then diametrically compressed to fracture, and the resultant fracture morphology was analyzed with a scanning electron microscope (SEM). A fracture gradient with three distinct regions across the thickness of the membranes was identified on the split surfaces. In the outer region of the C-ring specimen exposed to N{sub 2}, a mixed inter/transgranular fracture with a predominant intergranular pattern was observed. In the middle section of the fracture surface, a characteristic transgranular fracture of the perovskite grains was found. At the inner region of the ring, a mixed inter/transgranular fracture with a predominant transgranular pattern occurred. The mechanism of gradient fractures was attributed both to chemically induced stresses caused by oxygen diffusion and to the formation of a separate phase of oxygen-deficient perovskite in the parent perovskite. The stresses generated were modeled by a point defect model. This work provides significant information on microstructure evolutions of tubular LSFCO membranes under graded reducing atmospheres.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Fatigue and fracture behavior of bulk metallic glasses and their composites

A fundamental understanding of the fatigue and fracture behavior of bulk metallic glasses (BMGs) and their composites is of critical significance for designing new BMG systems and developing new manufacturing and processing techniques so as to broaden the scope of applications of BMGs and their composites. However, the fatigue and fracture studies on BMGs are limited so far, compared to other mechanical properties. The present work reviews the fatigue and fracture behavior of BMGs and their composites, as well as that of metallic-glass films, ribbons, and wires. The grand challenge for the fatigue and fracture performance of BMGs is: What produces a large difference among the fatigue and fracture results of BMGs? According to the fatigue and fracture investigations of crystalline alloys including recently invented high entropy alloys, many factors could be involved, such as the composition, material quality, specimen geometry, chemical environment, surface condition, temperature, cyclic frequency, mean stress, and residual stress, etc. Based on this challenge, the present work will review and address the factors affecting the fatigue and fracture behavior of BMGs and their composites. Here, the mechanisms of fatigue-crack initiation, propagation, and fracture of BMGs and their composites in different loading conditions and environments will be outlined, analyzed, and discussed. Future research directions of fatigue and fracture of BMGs and their composites are provided for reference.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Chemical effects on subcritical fracture in silica from molecular dynamics simulations

Fracture toughness of silicates is reduced in aqueous environments due to water-silica interactions at the crack tip. In this paper, to investigate this effect, classical molecular dynamics (MD) simulations using the bond-order based reactive force field (ReaxFF) were used to simulate silica fracture. The chemical and mechanical aspects were separated by simulating fracture in: (a) vacuum with dynamic loading, (b) an aqueous environment with dynamic loading, and (c) an aqueous environment with static subcritical mechanical loading to track silica dissolution. The addition of water to silica fracture reduced the silica fracture toughness by ~25%, a trend consistent with experimentally reported results. Analysis of Si-O bonds in the process zone and calculations of dissipation energy associated with fracture indicated that water relaxes the entire process zone, and not just the surface. Additionally, the crack tip sharpens during fracture in water and an increased number of microscopic propagation events occur. This results in earlier fracture in systems with increasing mechanical loading in aqueous conditions, despite the lack of significant silica dissolution. Finally and therefore, the threshold for Si-O bond breakage has been lowered in the presence of water and the reduction in fracture toughness is due to structural and energetic changes in the silica, rather than specific dissolution events.

58 GEOSCIENCES↗
DOE OSTI2025

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

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

02 PETROLEUM↗
DOE OSTI2026

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.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Transforming microseismic clouds into near real-time visualization of the growing hydraulic fracture

SUMMARY Microseismic observations during unconventional reservoir stimulation are typically seen as a proxy for clusters of hydraulic fractures and the extent of the stimulated reservoir. Such straightforward interpretation is often misleading and fails to provide a physically reasonable image of the fracturing process. This paper demonstrates the application of a physics-based machine learning algorithm which enables a rapid and accurate fracture mapping from the microseismic data. Our training and validation data set relies on a history-matched geomechanical modelling workflow implemented in GEOS software for the Hydraulic Fracturing Test Site 1 (HFTS-1) project. For this study we augmented the simulated fracture growth through geostatistical modelling of induced seismicity, so that the synthetic microseismic catalogue matches the main statistical properties of the field observations. We formulated the problem of mapping the actual fracture in the clutter of events to parallel common video segmentation workflows: several past video frames (microseismic density snapshots) are passed through a deep convolutional network to classify whether a given voxel is associated with a fracture or intact rock. We found that for accurate fracture mapping, the network’s input and architecture must be augmented to incorporate the fluid injection parameters (pressure, rate, concentration of proppant, and location of the perforation within the cluster). The error rate for the network reached as little as 10 per cent of the fracture area, while a conventional microseismic interpretation approach yielded ∼300 per cent. Our approach also yields must faster predictions than conventional methods (minutes instead of weeks), and could enable engineers to make rapid decisions regarding engineering parameters (pumping rate, viscosity) in real time during stimulation.

58 GEOSCIENCES↗
DOE OSTI2015

Challenges and opportunities for fractured rock imaging using 3D cross-borehole electrical resistivity

There is an increasing need to characterize discrete fractures away from boreholes to better define fracture distributions and monitor solute transport. We performed a 3D evaluation of static and time-lapse cross-borehole electrical resistivity tomography (ERT) data sets from a limestone quarry in which flow and transport are controlled by a bedding-plane feature. Ten boreholes were discretized using an unstructured tetrahedral mesh, and 2D panel measurements were inverted for a 3D distribution of conductivity. We evaluated the benefits of 3D versus 2.5D inversion of ERT data in fractured rock while including the use of borehole regularization disconnects (BRDs) and borehole conductivity constraints. High-conductivity halos (inversion artifacts) surrounding boreholes were removed in static images when BRDs and borehole conductivity constraints were implemented. Furthermore, applying these constraints focused transient changes in conductivity resulting from solute transport on the bedding plane, providing a more physically reasonable model for conductivity changes associated with solute transport at this fractured rock site. Assuming bedding-plane continuity between fractures identified in borehole televiewer data, we discretized a planar region between six boreholes and applied a fracture regularization disconnect (FRD). Although the FRD appropriately focused conductivity changes on the bedding plane, the conductivity distribution within the discretized fracture was nonunique and dependent on the starting homogeneous model conductivity. Synthetic studies performed to better explain field observations showed that inaccurate electrode locations in boreholes resulted in low-conductivity halos surrounding borehole locations. These synthetic studies also showed that the recovery of the true conductivity within an FRD depended on the conductivity contrast between the host rock and fractures. Our findings revealed that the potential exists to improve imaging of fractured rock through 3D inversion and accurate modeling of boreholes. However, deregularization of localized features can result in significant electrical conductivity artifacts, especially when representing features with a high degree of spatial uncertainty.

resistivity, fractured rock, inversion, imaging, g↗
DOE OSTI2021

Impact Fracture and Fragmentation of Glass via the 3D Combined Finite-Discrete Element Method

A driving technical concern for the automobile industry is their assurance that developed windshield products meet Federal safety standards. Besides conducting innumerable glass breakage experiments, product developers also have the option of utilizing numerical approaches that can provide further insight into glass impact breakage, fracture, and fragmentation. The combined finite-discrete element method (FDEM) is one such tool and was used in this study to investigate 3D impact glass fracture processes. To enable this analysis, a generalized traction-separation model, which defines the constitutive relationship between the traction and separation in FDEM cohesive zone models, was introduced. The mechanical responses of a laminated glass and a glass plate under impact were then analyzed. For laminated glass, an impact fracture process was investigated and results were compared against corresponding experiments. Correspondingly, two glass plate impact fracture patterns, i.e., concentric fractures and radial fractures, were simulated. The results show that for both cases, FDEM simulated fracture processes and fracture patterns are in good agreement with the experimental observations. The work demonstrates that FDEM is an effective tool for modeling of fracture and fragmentation in glass.

42 ENGINEERING↗
DOE OSTI2022

Postmortem and Antemortem Forensic Assessment of Pediatric Fracture Healing from Radiographs and Machine Learning Classification

A timeline of pediatric bone healing using fracture healing characteristics that can be assessed solely using radiographs would be practical for forensic casework, where the fracture event may precede death by days, months, or years. However, the dating of fractures from radiographs is difficult, imprecise, and lacks consensus, as only a few aspects of the healing process are visible on radiographs. Multiple studies in both the clinical and forensic literature have attempted to develop a usable scale to assess pediatric bone healing on radiographs using various healing characteristics. In contrast to the orthopedic definition, a fracture in forensic casework is only considered to be healed when the area around the fracture has been remodeled to the point that the fracture is difficult to detect on a radiograph or on the surface of the bone itself, a process that can take several years. We subjectively assessed visible characteristics of healing in radiograms of fractures occurring in 942 living children and adolescents. By dividing these assessments into learning and test (validation) sets, the accuracy of a newly proposed fracture healing scale was compared to a previous study. Two machine learning models were used to test predictions of the new scale. All three models produced similar estimates with substantial imprecision. Results corroborate the Malone model with an independent dataset and support the efficacy of using less complex models to estimate fracture age in children.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Fractal Characterization of Multimodal, Multiscale Images of Shale Rock Fracture Networks

An array of multimodal and multiscale images of fractured shale rock samples and analogs was collected with the aim of improving the numerical representation of fracture networks. 2D and 3D reconstructions of fractures and matrices span from 10 –6 to 100 m. The origin of the fracture networks ranged from natural to thermal maturation to hydraulically induced maturation. Images were segmented to improve fracture identification. Then, the fractal dimension and length distribution of the fracture networks were calculated for each image dataset. The resulting network connectivity and scaling associations are discussed at length on the basis of scale, sample and order of magnitude. Fracture network origin plays an important role in the resulting fracture systems and their scaling. Rock analogs are also evaluated using these descriptive tools and are found to be faithful depictions of maturation-induced fracture networks.

58 GEOSCIENCES↗
DOE OSTI2024

Fracture Network Influence on Rock Damage and Gas Transport following an Underground Explosion

Simulations of rock damage and gas transport following underground explosions that omit preexisting fracture networks in the subsurface cannot fully characterize the influence of geo-structural variability on gas transport. Previous studies do not consider the impact that fracture network structure and variability have on gas seepage. In this study, we develop a sequentially coupled, axi-symmetric model to look at the damage pattern and resulting gas breakthrough curves following an underground explosion given different fracture network realizations. We simulate 0.327 and 0.164 kT chemical explosives with burial depths of 100 m for 90 stochastically generated fracture networks. Gases quickly reach the surface in 30% of the higher yield simulations and 5% of the lower yield simulations. The fast breakthrough can be attributed to the formation of connected pathways between fractures to the surface. The formation of a connected damage pathway to the surface is not clearly correlated with the fracture intensity (P32) in our simulations. Breakthrough curves with slower transport are highly variable depending on the fracture network sample. The variability in the breakthrough behavior indicates that ignoring the influence of fracture networks on rock damage, which strongly influences the hydraulic properties following an underground explosion, will likely lead to a large underestimation of the uncertainty in the gas transport to the surface. This work highlights the need for incorporation of fracture networks into models for accurately predicting gas seepage following underground explosions.

58 GEOSCIENCES↗