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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Geologic stress modulates fluid mixing at fracture intersections

Fracture intersections are critical links that enable flow and transport in subsurface fracture networks, and their behavior strongly influences fluid mixing in a network. Although all subsurface fractures are subjected to geological stress, we lack a fundamental understanding of how fracture intersection geometry evolves under stress and how these changes influence fluid mixing. Here, we combine 3D printing, 3D X-ray tomographic imaging, and 3D pore-scale numerical simulations to reveal stress-induced changes in intersection geometry and their impact on mixing. Mixing is found to be strongly affected by partial closure of an intersection under stress. As an intersection closes, the void area for fluid flow and diffusion decreases leading to substantial deviations between conventional mixing models and full pore-scale modeling. To address this, we propose a modified mixing model that accounts for intersection deformation, which is essential for accurate modeling of solute transport and mixing through fracture networks.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

A New Workflow of X-ray CT Image Processing and Data Analysis of Structural Features in Rock Using Open-Source Software

X-ray computed tomography (CT) images of rock specimens often contain artifacts which must be corrected before scientific analyses are performed. Here, we present a new workflow of automated image processing to utilize poor-quality X-ray CT scan images. The workflow runs on the open-source image analysis software and efficiently separates desired features from low-contrast scanned images. The new workflow is a two-step technique using contrast enhancement and automated feature segmentation to generate noise-free binary images. The results of binary images using the proposed workflow and using a conventional thresholding technique are analyzed to show the quality of the proposed method. The paper also presents a workflow of estimating the structural geometries of features in two and three dimensions. The results of the structural feature analyses and computational time were compared between the open-source (ImageJ) and commercial image analysis software (Bruker Computed Tomography Analyzer). The commercial software was more computationally efficient, but the task-specific macros in open-source software enabled the user-desired automation in image processing and data extraction of desired structural features of comparable quality.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Numerical and Analytical Modeling of Flow Partitioning in Partially Saturated Fracture Networks

Infiltration processes in fractured-porous media remain a crucial, yet not very well understood component of recharge and vulnerability assessment. Under partially-saturated conditions flows in fractures, percolating fracture networks, and fault zones contribute to the fastest spectrum of infiltration velocities via preferential pathways. Specifically, the partitioning dynamics at fracture intersections determine the magnitude of flow fragmentation into vertical and horizontal components and hence the bulk flow velocity and dispersion of fracture networks. In this work, we derive an approximate analytical solution for the partitioning process and validate it using smoothed particle hydrodynamics simulations. The developed transfer function is conceptually based on simulation results and laboratory experiments carried out in previous works. It allows to efficiently simulate flow through fracture networks with simple cubic structure and arbitrary number of fractures and aperture sizes via linear response theory and convolution of a given input signal. We derive a non-dimensional bulk flow velocity ($\widetilde{v}$) and dispersion coefficient ($\widetilde{D}$) to characterize fracture networks in terms of dimensionless horizontal and vertical time scales $\tau_m$ and $\tau_0$. The dispersion coefficient is shown to strongly depend on the horizontal time scale and converges towards a constant value of $0.08$ within reasonable fluid and geometrical parameter ranges, while the non-dimensional velocity exhibits a characteristic $\widetilde{v} \sim \tau_m^{-1/2}$ scaling. Given that hydraulic information is often only available at limited places within (fractured-porous) aquifer system, such as boreholes or springs, our study intends to provide a rudimentary analytical concept to potentially reconstruct internal fracture network geometries from external boundary information, e.g., the dispersive properties of discharge (groundwater level fluctuations).

42 ENGINEERING↗
DOE OSTI2022

Opening-mode fracturing and cementation during hydrocarbon generation in shale: An example from the Barnett Shale, Delaware Basin, West Texas

Relative timing of fracturing is a key input for predictive fracture models, but timing information for fractures is commonly difficult to obtain. In this study, we used crosscutting relations and fluid inclusion assemblage temperatures from fracture cements from a few well-documented sampled fractures, combined with a one-dimensional burial history model, to establish timing for three generations of opening-mode fractures in a Barnett Shale core from the southern part of the Delaware Basin, Pecos County, West Texas. A burial history model is presented for the cored well and matched to measured vitrinite reflectance in samples from the core, and bottomhole temperature in the well. Here, the earliest fractures (group 1) likely formed due to early fluid-expulsion events (ca. 300 Ma) and were folded during host-rock compaction. Later group 2 fractures are sealed with fibrous barite containing primary, liquid hydrocarbon inclusions (mean homogenization temperature [T h ] = –9°C) and aqueous fluid inclusions (mean T h = 108.1°C). Group 2 fractures likely formed in response to fluid overpressure associated with cracking of type II kerogen to oil. Group 3 vertical fractures are up to 2 m in height with kinematic apertures ranging from less than 0.05 to 1.4 mm, partly open, and strike dominantly 010°–020°. Sequentially trapped aqueous fluid inclusions in fracture-spanning quartz cement bridges (mean T h = 110°C in crack-seal texture and 128°C in post-crack-seal fracture cement) record fracture opening under increasing temperature, inferred to reflect increasing burial, with continued overpressuring during the Triassic to Late Cretaceous. Some group 3 fractures may have continued to fill during Cenozoic uplift.

02 PETROLEUM↗
DOE OSTI2020

An embedded 3D fracture modeling approach for simulating fracture-dominated fluid flow and heat transfer in geothermal reservoirs

In this paper, we describe an efficient modeling approach, named embedded discrete fracture method (EDFM), for incorporating arbitrary 3D, discrete fractures, such as hydraulic fractures or faults, into modeling fracture-dominated fluid flow and heat transfer in fractured geothermal reservoirs. This technique allows 3D discrete fractures to be discretized independently from surrounding rock volume and inserted explicitly into a primary fracture/matrix grid, generated without including 3D discrete fractures in prior. An effective computational algorithm is developed to discretize these 3D discrete fractures and construct local connections between 3D fractures and fracture/matrix grid blocks representing the surrounding rock volume. The constructed gridding information on 3D fractures is then added to the primary grid. This embedded fracture modeling approach can be directly implemented into a developed geothermal reservoir simulator via the integral finite difference (IFD) method or with TOUGH2 technology. This embedded fracture modeling approach is very promising and computationally efficient to handle realistic 3D discrete fractures with complicated geometries, connections, and spatial distributions. Compared with other fracture modeling approaches, it avoids cumbersome 3D unstructured, local refining procedures, and increases computational efficiency by simplifying Jacobian matrix size and sparsity, while maintaining enough accuracy. Several numeral simulations are presented to demonstrate the utility and robustness of the proposed technique. Our numerical experiments show that this approach captures all the key patterns about fluid flow and heat transfer dominated by fractures in these cases. Thus, this approach is readily available to the simulation of fractured geothermal reservoirs with both artificial and natural fractures.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Specimen Size and Geometry Effects on the Master Curve Fracture Toughness Measurements of EUROFER97 and F82H Steels

EUROFER97 and F82H are two leading reduced-activation ferritic-martensitic (RAFM) steels for fusion blanket applications. Exposure to the harsh environment of fusion reactors can result in severe degradation of materials fracture toughness (FT). Thus, the post-irradiation evaluation of FT is critical to understanding the material behavior. Due to the space constraint of irradiation facilities, the development of small specimen test techniques (SSTT) is necessary to evaluate the performance of irradiated materials. In this study, we evaluated the specimen size and geometry effects on the ductile-to-brittle transition FT of EUROFER97 batch-3 and F82H-BA12 steels. The specimen thicknesses ranged from 1.65 to 12.7 mm and the geometries included 1.65 mm bend bar, 4 mm mini-compact tension (miniCT), and 0.5T compact tension (CT) specimens. Fracture toughness testing and evaluations were performed using the Master Curve method in the ASTM E1921-19 standard. After size correction to 1T size using the Master Curve method, no specimen size effect was observed between the 4 mm miniCT and 0.5T CT specimens for the Master Curve reference temperature T0Q, while the bend bars yielded a higher T0Q. A strong effect of fatigue precrack front straightness on T0Q for 0.5T CT specimens was observed. The minimum number of specimens needed for each specimen geometry has been determined.

Chen, Xiang↗
DOE OSTI2025

Fracture toughness measurement of ceramics with fuel pellet geometry

The crack initiation of nuclear fuel pellets, driven by extreme temperature gradients, is a critical factor impacting reactor performance. Although this fact is widely acknowledged, there is limited fracture toughness data for specimens with representative geometry and fabricated using relevant processing methods. Here, to address this, a new chevron-notched short-rod specimen configuration with diameter and length both equal to 10 mm, representative of a conventional fuel pellet and tested in tension, was developed in an attempt to validly measure Mode I fracture toughness, K Ic . Test calibration was carried out by testing polycrystalline ZrO 2 , SiC, and Si 3 N 4 ceramics fabricated with conventional fuel pellet geometry and known fracture toughnesses. The results suggest that the method described in this paper produces valid K Ic measurements using conventional ceramic fuel pellets.

Ceramics↗
DOE OSTI2022

Fracture toughness of schist, amphibolite, and rhyolite from the Sanford Underground Research Facility (SURF), Lead, South Dakota

The Cracked Chevron Notched Brazilian Disc (CCNBD) method was selected for Mode I fracture toughness tests on Poorman schist, Yates amphibolite, and rhyolite dikes from the EGS Collab site at the SURF in Lead, South Dakota. The effects of lithology, anisotropy, and loading rate were investigated. Fracture toughness was greatest in amphibolite, with schist and rhyolite having similar toughness values ($K$ amphibolite > $K$ rhyolite ≈ $K$ schist ). The effects of anisotropy on fracture toughness were investigated in the foliated schist samples. Schist samples were prepared in three geometries (divider, arrester, and short transverse) which controlled how the fracture would propagate relative to foliations. The divider geometry was strongest and short transverse geometry was the weakest ($K$ divider > $K$ arrester > $K$ shorttransverse ). Fracture toughness was observed to decrease with decreasing loading rate. Optical and SEM microscopy revealed that for the short transverse geometry, fractures tended to propagate along grain boundaries, whereas in arrester and divider geometries fractures tended to propagate through grains. In foliated samples, the tortuosity of the fracture observed in thin section was greater in arrester and divider geometries than in short transverse geometries.

58 GEOSCIENCES↗
DOE OSTI2024

Evaluation of Lost Circulation Material Sealing for Geothermal Drilling

Lost circulation is a pervasive problem in geothermal wells that can create prohibitive costs during drilling. The main issue with treatment is that the mechanism of plug formation is poorly understood. Here we applied two experimental approaches to characterize the clogging effectiveness of different materials. Fracture flow tests with different geometries were conducted with various individual materials and mixtures at relevant conditions. A high-temperature flow loop system was also developed to inject single- and mixed-material plugs into a gravel pack with a non-uniform geometry to compare with the fracture tests. The fracture tests revealed that single materials tended to form no plug or an unstable plug, while mixtures of materials were uniformly better at sealing fractures. Gravel pack tests at high temperatures show most of the materials are intact but degraded. The fibrous materials can create partial or unstable plugs in the gravel pack, but mixed-material plugs are far more effective at clogging. Both test types suggest that (1) mixed materials are more effective at blocking fluid flow and (2) fibrous materials seal fracture openings better, while granular materials seal inside fractures or pore throats better. Further research is needed to study the long-term stability of different plug configurations.

clogging tests↗
DOE OSTI2020

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.

42 ENGINEERING↗
DOE OSTI2021

Computationally Efficient Multiscale Neural Networks Applied to Fluid Flow in Complex 3D Porous Media

Abstract The permeability of complex porous materials is of interest to many engineering disciplines. This quantity can be obtained via direct flow simulation, which provides the most accurate results, but is very computationally expensive. In particular, the simulation convergence time scales poorly as the simulation domains become less porous or more heterogeneous. Semi-analytical models that rely on averaged structural properties (i.e., porosity and tortuosity) have been proposed, but these features only partly summarize the domain, resulting in limited applicability. On the other hand, data-driven machine learning approaches have shown great promise for building more general models by virtue of accounting for the spatial arrangement of the domains’ solid boundaries. However, prior approaches building on the convolutional neural network (ConvNet) literature concerning 2D image recognition problems do not scale well to the large 3D domains required to obtain a representative elementary volume (REV). As such, most prior work focused on homogeneous samples, where a small REV entails that the global nature of fluid flow could be mostly neglected, and accordingly, the memory bottleneck of addressing 3D domains with ConvNets was side-stepped. Therefore, important geometries such as fractures and vuggy domains could not be modeled properly. In this work, we address this limitation with a general multiscale deep learning model that is able to learn from porous media simulation data. By using a coupled set of neural networks that view the domain on different scales, we enable the evaluation of large ( $$>512^3$$ > 512 3 ) images in approximately one second on a single graphics processing unit. This model architecture opens up the possibility of modeling domain sizes that would not be feasible using traditional direct simulation tools on a desktop computer. We validate our method with a laminar fluid flow case using vuggy samples and fractures. As a result of viewing the entire domain at once, our model is able to perform accurate prediction on domains exhibiting a large degree of heterogeneity. We expect the methodology to be applicable to many other transport problems where complex geometries play a central role.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Flaw-insensitive fracture of a micrometer-sized brittle metallic glass

Brittle materials, such as oxide glasses, are usually very sensitive to flaws, giving rise to a macroscopic fracture strength that is much lower than that predicted by theory. The same applies to metallic glasses (MGs), with the important difference that these glasses can exhibit certain plastic strain prior to catastrophic failure. Here we consider the strongest metallic alloy known, a ternary Co 55 Ta10B 35 MG. In this work, we show that this macroscopically brittle glass is flaw-insensitive at the micrometer scale. This discovery emerges when testing pre-cracked specimens with self-similar geometries, where the fracture stress does not decrease with increasing pre-crack size. The fracture toughness of this ultra-strong glassy alloy is further shown to increase with increasing sample size. Both these findings deviate from our classical understanding of fracture mechanics, and are attributed to a transition from toughness-controlled to strength-controlled fracture below a critical sample size.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Fracture Description of the HFTS-2 Slant Core, Delaware Basin, West Texas

Previous work in the Midland Basin Hydraulic Fracture Test Site project (HFTS1) has shown the value of core recovered from a stimulated volume in providing information on the geometry and extent of hydraulic fractures. This project (HFTS2) is in the Wolfcamp Formation of the Delaware Basin. A slant core through the stimulated volume was acquired above and below a stimulated well and the results of the fracture description are presented here. Here, nine, 3-inch-diameter whole cores from the slant well were examined for fractures prior to slabbing, together with a CT scan of the core. We used criteria developed during the HFTS1 project to identify hydraulic, natural, drilling-induced and core-handling fractures. Material collected from a few fracture faces was cleaned, sieved, and examined under a microscope to look for proppant. A total of 1261 fractures, including 500 hydraulic fractures, were described in 948 feet of core. Hydraulic fractures are not evenly distributed through the cores; they tend to occur in clusters. Many are remarkably planar, smooth, and featureless, while others have twist hackles, steps, plumose, and kinks. Natural, subvertical fractures, mostly less than 1 mm in width, sealed with blocky calcite, are oriented NE-SW (Set 1) and WNW-ESE (Set 2). Two Set 1 fractures have drusy cement on parted faces and were possibly open prior to drilling. A few subvertical sealed fractures with different orientations are also present. Bed-parallel natural fractures, most having fibrous cement (beef), and some showing top-to-north bed-parallel shear are present. Faults with normal/oblique movement were found in one part of the core. Drilling induced and core-handling breaks were identified in all cores, and were also noted to have reactivated some natural fractures. The results from this work will allow hydraulic fracture distribution relative to the stimulated well to be assessed. The number of fractures is generally greater than the number of perforations, but these can be compared for different stages and perforation clusters. Surface features provide information about the direction of hydraulic fracture propagation and about segmentation and bifurcation that might increase the number of fracture strands. Reactivation of natural fractures, including bedding-parallel fractures, occurs although the mechanism of parting is mostly unknown. The wider HFTS2 project will make use of this fracture description to improve hydraulic fracture modeling (geomechanics and engineering), to help explain pressure depletion and fiber optics observations (reservoir depletion and frac geometry), to help verify indirect fracture diagnostic techniques such as microseismic monitoring (geophysics), and to help calibrate horizontal well image logs.

58 GEOSCIENCES↗
DOE OSTI2021

New Fracture Diagnostic Tool for Unconventionals: High-Resolution Distributed Strain Sensing via Rayleigh Frequency Shift during Production in Hydraulic Fracture Test 2

Fiber Optic monitoring in unconventional reservoirs has proven to be an invaluable diagnostic tool for assessing both near-wellbore stimulation effectiveness and to help describe the far-field frac geometries created by hydraulic fracture stimulation. Unfortunately, gaining any detailed qualitative and quantitative understanding of the near-wellbore frac geometry or cluster/stage productivity during production via Fiber Optic (FO) has proven to be more difficult, particularly in wells producing liquids. A new FO diagnostic method, Distributed Strain Sensing based on Rayleigh Frequency Shift (DSS-RFS), first demonstrated for oil and gas applications in the Hydraulic Test Site 2 (HFTS2) provides new insights about the characteristics of near-wellbore-region (NWR) during production. DSS-RFS is different from other FO strain measurements because it relies on accurate measurement of frequency shifts of Rayleigh backscattered spectrum obtained by scanning the fiber with a coherent optical time-domain reflectometer with a range of laser frequencies using a tunable-wavelength laser system. Changes in strain are measured with an extremely high spatial resolution of 20 cm and with high signal-to-noise ratios over long distances. In HFTS2, strain changes for the entire wellbore have been measured twice during scheduled shut-in and reopening operations (February 2020 and September 2020). After removing temperature effects, consistent strain changes have been observed at the location of most perforation clusters. These are caused by near wellbore fracture aperture changes due to pressure increases during shut-in within the near-wellbore fracture network. The strain-change patterns from the DSS-RFS during shut-in correlate very well with the location of clusters and allow for the definition of extending intervals with positive strain signals at each cluster and slightly compressing intervals with negative strain signals between the clusters and in the non stimulated intervals. The locations of the measured positive strain peaks also show good correspondence to DAS acoustic intensity measurements acquired during the stimulation. The geometry and magnitude of the strain changes differ significantly between the two tested completion designs in the same well. During shut in and reopening each cluster exhibit its own strain-change / pressure path. In addition, the September 2020 dataset also revealed the existence of small but measurable strain changes as consequence of pressure decline during production. These strain changes also correlate well with the presence of producing clusters, but the strain-rate signals are opposite to that obtained during shut-in and reopening operations. Although Downloaded from http://onepetro.org/URTECONF/proceedings-pdf/21URTC/2-21URTC/D021S031R002/2477551/urtec-2021-5408-ms.pdf/1 by Carol Worster on 28 February 2022 URTeC 5408 2 we are still in the early stages of exploring the potential of this novel FO technique, we believe that the highly detailed information contained in the measurement of strain changes using DSS-RFS during production can significantly improve our understanding of near-wellbore hydraulic fracture characteristics and the relationships between stimulation and production from unconventional oil and gas wells.

04 OIL SHALES AND TAR SANDS↗