Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fracture control”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Understanding Water Controls on Shale Gas/Oil Mobilization into Fractures (FY 2019-2021(FY2022, No Cost Extension), Final Report)

The large volumes of injected and produced water, and costs associated with both hydraulic fracturing fluids and produced water treatment motivate development of hydraulic fracturing strategies that use water more efficiently. The overall objective of this project is to improve understanding on how water is distributed within unconventional shales during and after hydraulic fracturing in order to provide a more rational basis for reservoir stimulation strategies that increase hydrocarbon recovery per unit of water used. The strong capillary imbibition of water into shale matrix pores and immiscibility with hydrocarbons makes it highly desirable to reduce water use during reservoir development. Given the importance, and competing influences of capillarity and gravity in controlling water in fractured porous media, this project also seeks to understand how these factors interact to obstruct or allow gas flow.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2020

Experimental and theoretical study of the restrained shrinkage cracking of early age well cement

The shrinkage cracking occurring in oil & gas well cement material is a major factor that reduces the cement sheath integrity and impairs the durability of an underground well. Shrinkage mainly takes place in the early hydration period of the cement, and the resulting restrained shrinkage cracking can form a leaking network for oil & gas or carbon dioxide stored in an underground reservoir. It is thus important to understand the early age behavior of well cement and the leaking pathway formation in the well. This study describes the application of the digital image correlation (DIC) method as a tool to measure the strain development and crack distributions of well cement in the lab during early phase changes. A theoretical model is presented to calculate the stress distribution built up in the cement sheath during early ages of well cement hydration, which helps to interpret the experimental observations and determine the major factors influencing the cracking sequences. This study presents the whole hydration – drying shrinkage – and cracking mechanism of early age well cement under restrained ring conditions experimentally and theoretically. Finally, the results offer useful insights on the key material and geometric parameters which control the fracture of a well cementing structure.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Nanomechanical investigation of the interplay between pore morphology and crack orientation of amorphous silica

Porous amorphous silica (a-SiO 2 ) is of both fundamental and practical interests, as they exhibit a large specific surface area and tunable porous network. However, the brittle nature of a-SiO 2 and the presence of pre-existing cracks at both micro- and nano-scales lead to complex mechanical behavior. Here, we systematically investigate the effects of pre-existing crack and its orientation on the mechanical properties of a-SiO 2 with varying pore shapes using reactive molecular dynamics simulations. We demonstrate that pore shape will primarily influence the Young’s modulus (E) and critical energy release rate (G IC ). We further investigate the impact of pore shape and crack orientation by local characterization of the structural parameters. By defining the high stress and inter-mediate regions, the overall mechanical properties are found to be greatly influenced by the pore shape which can be reflected through the spatial distribution of von Mises stress. Overall, G IC is found to increase with the increase of ligament length (also known as pore wall thickness). Meanwhile, the effect of the pre-existing crack on the crack propagation process is confirmed by analyzing the density distribution evolution. These results highlight the interplay between pore morphology and crack orientation in controlling the fracture behaviors in brittle porous materials.

42 ENGINEERING↗
DOE OSTI2022

Stress-controlled fatigue of HfNbTaTiZr high-entropy alloy and associated deformation and fracture mechanisms

The stress-controlled fatigue tests are carried out at a stress ratio of 0.1 and a frequency of 10 Hz, and span both low-cycle and high-cycle regimes by varying the applied stress amplitudes. The high-cycle fatigue regime gives a fatigue strength of 497 MPa and a fatigue ratio of 0.44. At equivalent conditions, the alloy's fatigue strength is greater than all other high-entropy alloys (HEAs) with reported high-cycle fatigue data, dilute body-centered cubic alloys, and many structural alloys such as steels, titanium alloys, and aluminum alloys. Through in-depth analyses of crack-propagation trajectories, fracture-surface morphologies and deformation plasticity by means of various microstructural analysis techniques and theoretical frameworks, the alloy's remarkable fatigue resistance is attributed to delayed crack initiation in the high-cycle regime, which is achieved by retarding the formation of localized persistent slip bands, and its good resistance to crack propagation in the low-cycle regime, which is accomplished by intrinsic toughening backed up by extrinsic toughening. Additionally, the stochastic nature of the fatigue data is neatly captured with a 2-parameter Weibull model.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Learning the factors controlling mineral dissolution in three-dimensional fracture networks: applications in geologic carbon sequestration

We perform a set of high-fidelity simulations of geochemical reactions within three-dimensional discrete fracture networks (DFN) and use various machine learning techniques to determine the primary factors controlling mineral dissolution. The DFN are partially filled with quartz that gradually dissolves until quasi-steady state conditions are reached. At this point, we measure the quartz remaining in each fracture within the domain as our primary quantity of interest. We observe that a primary sub-network of fractures exists, where the quartz has been fully dissolved out. This reduction in resistance to flow leads to increased flow channelization and reduced solute travel times. However, depending on the DFN topology and the rate of dissolution, we observe substantial variability in the volume of quartz remaining within fractures outside of the primary subnetwork. This variability indicates an interplay between the fracture network structure and geochemical reactions. We characterize the features controlling these processes by developing a machine learning framework to extract their relevant impact. Specifically, we use a combination of high-fidelity simulations with a graph-based approach to study geochemical reactive transport in a complex fracture network to determine the key features that control dissolution. We consider topological, geometric and hydrological features of the fracture network to predict the remaining quartz in quasi-steady state. We found that the dissolution reaction rate constant of quartz and the distance to the primary sub-network in the fracture network are the two most important features controlling the amount of quartz remaining. This study is a first step towards characterizing the parameters that control carbon mineralization using an approach with integrates computational physics and machine learning.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Stratigraphically Controlled Stress Variations at the Hydraulic Fracture Test Site-1 in the Midland Basin, TX

We investigated the relationship between stratigraphy, stress, and microseismicity at the Hydraulic Fracture Test Site-1. The site comprises two sets of horizontal wells in the Wolfcamp shale and a deviated well drilled after hydraulic fracturing. Regional stresses indicate normal/strike-slip faulting with E-W compression. Stress measurements in vertical and horizontal wells show that the minimum principal stress varies with depth. Strata with high clay and organic content show high values of the least compressive stress, consistent with the theory of viscous stress relaxation. By integrating data from core, logs, and the hydraulic fracturing stages, we constructed a stress profile for the Wolfcamp sequence, which predicts how much pressure is required for hydraulic fracture growth. We applied the results to fracture orientation data from image logs to determine the population of pre-existing faults that are expected to slip during stimulation. We also determined microseismic focal plane mechanisms and found slip on steeply dipping planes striking NW, consistent with the orientations of potentially active faults predicted by the stress model. This case study represents a general approach for integrating stress measurements and rock properties to predict hydraulic fracture growth and the characteristics of injection-induced microseismicity.

microseismicity↗
DOE OSTI2021

The influence of preexisting host rock damage on fault network localization

The transition from stable to unstable fracture propagation occurs when fractures begin to interact and link. Thus, fracture network coalescence controls how rocks and engineered structures fail. Here, to constrain the factors that influence localization in shear zones under brittle conditions, we build discrete element method models with a rough fault embedded in a shear zone. We add varying numbers of diffuse, randomly-placed weaknesses to examine the influence of diffuse damage on fracture network localization. The number of weaknesses controls the localization behavior of the fault network and the final fault geometry. We quantify localization using the Gini coefficient of the fracture volume, which measures the nonuniformity in a population. Each model generally increases in localization toward failure. However, models with more diffuse damage experience delocalization phases that are superimposed on the overall trend of increasing localization. The observed link between delocalization and host rock damage may help explain the varying localization of low magnitude seismicity in southern California. Models with more diffuse damage produce more complex fault geometries comprised of several parallel strands of wing cracks. The propagation of these wing cracks reduces the shear stress acting on the model boundaries, indicating that this fracture development increases the mechanical efficiency of the system.

58 GEOSCIENCES↗
DOE OSTI2020

Microstructure and fracture toughness characterization of three 9Cr ODS EUROFER steels with different thermo-mechanical treatments

Ferritic martensitic ODS steels are one of the candidate structural materials for future Gen-IV nuclear fission and fusion reactors. Here, the dependence of fracture toughness on microstructure was investigated by comparing three 9Cr ODS EUROFER steels manufactured through different thermo-mechanical processing routes. Quasi-static fracture toughness testing was performed with sub-sized C(T) specimens and microstructural characterization was carried out using scanning electron microscopy, electron backscatter diffraction, and transmission electron microscopy. It was found that at lower test temperatures (-100 – 22 °C), the fracture toughness was primarily controlled by crack initiation at sub-micron particles and by production of secondary cracks during fracture. At higher temperatures (above 100 °C), fracture toughness was predominantly controlled by the matrix ductility and the grain boundary strength with a relatively ductile coarse-grained alloy demonstrating higher fracture toughness compared to high-strength fine-grained alloys. These results and discussion show that variations in thermomechanical treatments can produce significant differences in microstructure and fracture toughness behavior of ferritic martensitic ODS steels.

36 MATERIALS SCIENCE↗
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 OSTI2022

Investigation of Coupled Processes in Fractures and the Bordering Matrix via a Micro‐Continuum Reactive Transport Model

Abstract In multi‐mineral fractured rocks, the altered porous layer on the fracture surface resulting from preferential dissolution of the fast‐reacting minerals can have profound impacts on subsequent chemical‐physical alteration of the fractures. This study adopts the micro‐continuum approach to provide further understanding of reactive transport processes in the altered layer (AL), and mass exchanges with the bordering matrix and fracture. The modeling framework couples the Darcy‐Brinkman‐Stokes (DBS) solver in COMSOL Multiphysics and the geochemical modeling capability of CrunchFlow. Three‐dimensional steady state simulations with systematically varied chemical‐physical parameters of the AL were performed to examine the impacts of individual factors and processes. Our simulation results confirm previous observations that dissolution of the fast‐reacting mineral (i.e., calcite) is largely controlled by diffusion across the AL. We also show that dissolution of the slow‐reacting mineral (i.e., dolomite), which controls AL development and fracture enlargement, increases with surface area and has a complex dependence on different local rate‐limiting processes. In particular, advection can result in evident spatial variations in the local dissolution rates of dolomite, although it does not affect the bulk chemistry significantly. The difference in the spatial patterns between simulations with and without advection in the AL is more noticeable in the locations with smaller apertures, with up to 20% difference in local reaction rates. Therefore, it is important to include a full depiction of advection, diffusion, and reactions for accurately capturing local dynamics that control long‐term fracture evolution.

58 GEOSCIENCES↗
DOE OSTI2021

Machine-learning predictions of the shale wells’ performance

The ultra-low permeability nature of shale reservoirs leads to an extended linear flow and necessitates horizontal wells with multi-stage engineered fractures to efficiently extract hydrocarbons resources. These artificially-generated and naturally-occurring fractures form complex networks that create complex flow regimes which control oil production. These fractures are neither identical nor equally-spaced, which leads to a production profile with a masked onset of the boundary-dominated flow. The combination of the extended linear flow with the indeterminate onset of the boundary-dominated flow challenges the current deterministic analytic approaches to forecast the estimated ultimate recovery (EUR). In this work, we propose a novel machine-learning approach which overcomes these challenges and provides reliable EUR estimates based on field-wide analyses. We implement a novel unsupervised machine learning (ML) methodology, which allows for automatic identification of the optimal number of features (signals) present in the data based on non-negative matrix/tensor factorization coupled with k-means clustering incorporating regularization and physics constraints. In the presented analyses, the input data to the ML algorithm is the available (public) production history from the field collected at existing unconventional reservoirs. We validate our approach through hindcasting of the production data, where we achieved an excellent agreement. In addition, our approach is able to identify the poorly-performing wells, which could benefit from early refracing. Our approach provides fast and accurate estimations of the well performance without presumptions about the state of the well or the flow regime.

03 NATURAL GAS↗
DOE OSTI2023

Chemical controls on the propagation and healing of subcritical fractures

Human activities involving subsurface reservoirs—resource extraction, carbon and nuclear waste storage—alter thermal, mechanical, and chemical steady-state conditions in these systems. Because these systems exist at lithostatic pressures, even minor chemical changes can cause chemically assisted deformation. Therefore, understanding how chemical effects control geomechanical properties is critical to optimizing engineering activities. The grand challenge in predicting the effect of chemical processes on mechanical properties lays in the fact that these phenomena take place at molecular scales, while they manifest all the way to reservoir scales. To address this fundamental challenge, we investigated chemical effects on deformation in model and real systems spanning molecular- to centimeter scales. We used theory, experiment, molecular dynamics simulation, and statistical analysis to (1) identify the effect of simple reactions, such as hydrolysis, on molecular structures in interfacial regions of stressed geomaterials; (2) quantify chemical effects on the bulk mechanical properties, fracture and displacement for granular rocks and single crystals; (3) develop initial understanding of universal scaling for individual displacement events in layered geomaterials; and (4) develop analytic approximations for the single-chain mechanical response utilizing asymptotically correct statistical thermodynamic theory. Taken together, these findings advance the challenging field of chemo-mechanics.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Inference of Induced Fracture Geometries Using Fiber-Optic Distributed Strain Sensing in Hydraulic Fracture Test Site 2

Fiber-optic Distributed Strain Sensor observations in a vertical well allow rapid inference of the time evolution of nearby hydraulic fracture heights and their vertical aperture distributions during stimulation and depressurization. We describe a computationally efficient deconvolutional model based on an elastostatic kernel for fracture opening within an isotropic, homogeneous, linear elastic medium. Furthermore, an application of this model to observations from the Hydraulic Fracture Test Site 2 reveals that ultimate fracture heights are typically greater than those indicated by microseismic monitoring, and some hydraulic fractures continue their upward growth for several hours after high-pressure injections ceased.

58 GEOSCIENCES↗