Search NASA⌕ Search

SEARCH · Search NASA

Results for “Fracturing”

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 109 records · Page 6

Seismic Analysis of Spatio-Temporal Fracture Generation During EGS Resource Development - Deviatoric MT, Fracture Network, and Final Report

This submission contains 167 deviatoric moment tensor (MT) solutions for the seismicity observed two years prior and three years post start of injection activities at The Geysers Prati 32 EGS Demonstration. Also included is a statistical representation of the properties of 751 fractures derived from the analysis of seismicity observed two years prior and three years post start of injection activities at The Geysers Prati 32 EGS Demonstration Project. The locations of the fractures are taken from microseismic hypocenters, the fracture areas are derived from moment magnitudes via scaling relationships, and the azimuths (sigma 1) and dips (sigma 3) are derived from the results of stress analyses.

15 GEOTHERMAL ENERGY↗

Post-Fracturing Experiment Simulation of Hydraulic Fracture Propagation in a Deep Mine Using a Fully Coupled 3D Network-Flow and Quasi-static Discrete Element Model

Aimed at supporting the hydraulic fracturing experiments in a vertical borehole in the Poorman phyllite at approximately 1500 m depth at the Sanford Underground Research Facility (SURF), we perform post-hydraulic fracturing experiment simulations and compare the simulated borehole pressure changes and flow rates with actual measurements using a newly developed fully coupled three-dimensional (3D) network flow and quasi-static discrete element model (DEM). The quasi-static DEM model, which is constructed by Delaunay tessellation of the rock volume, considers rock fabric heterogeneities by using the “disordered” DEM mesh and adding random perturbations to the stiffness and tensile/shear strengths of individual DEM elements and the elastic beam between them. A conjugate 3D flow network based on DEM lattice is constructed to calculate the fluid flow in both the fracture and porous matrix. One distinctive advantage of the DEM model is that fracturing is naturally described by the breakage of elastic beams between DEM elements. It is also extremely convenient to introduce mechanical anisotropy into the model by simply assigning orientation-dependent tensile/shear strengths to the elastic beams.

15 GEOTHERMAL ENERGY↗

Strain-dependent partial slip on rock fractures under seismic-frequency torsion: Seismic-Frequency Fracture Partial Slip

Measurements of nonlinear modulus and attenuation of fractures provide the opportunity to probe their mechanical state. We have adapted a low-frequency torsional apparatus to explore the seismic signature of fractures under low normal stress, simulating low effective stress environments such as shallow or high pore pressure reservoirs. We report strain-dependent modulus and attenuation for fractured samples of Duperow dolomite (a carbon sequestration target reservoir in Montana), Blue Canyon Dome rhyolite (a geothermal analog reservoir in New Mexico), and Montello granite (a deep basement disposal analog from Wisconsin). We use a simple single effective asperity partial slip model to fit our measured stress-strain curves, and solve for the friction coefficient, contact radius, and full slip condition. These observations have the potential to develop into new field techniques for measuring differences in frictional properties during reservoir engineering manipulations and estimate the stress conditions where reservoir fractures and faults begin to fully slip.

58 GEOSCIENCES↗

A novel targeted plugging and fracture‐adaptable gel used as a diverting agent in fracturing

Abstract Diverting agents that could plug a targeted location and fill different shapes of fractures have attracted an increasing attention in the field of fracturing. In this study, a series of novel diverting agents were prepared by mixing aldehyde sodium alginate (ASA) and amino gelatin (AG) through a Schiff base reaction to form a primary network, polylactic acid (PLA) fibers were added to produce a 3D network, and Ca 2+ was mixed in as an ionic cross‐linking agent to strengthen the cross‐linking points. It was found that gelation time, swelling behavior, and mechanical performance could be controlled by varying the content of aldehyde groups and amino groups. Gelation time could be controlled in 1‐10 minutes. Plugging behavior and degradation tests demonstrated that the hydrogel could suit fractures with various shapes and plug effectively under high temperatures and pressures, and it could eventually become acidic fluid to help clean the wellbore. This novel ASA/AG/PLA hydrogel constitutes a promising environmentally friendly diverting agent in the application of diversion fracturing.

Zhang, Yuyuan↗

Acid Erosion of Carbonate Fractures and Accessibility of Arsenic-Bearing Minerals: In Operando Synchrotron-Based Microfluidic Experiment

Underground flows of acidic fluids through fractured rock can create new porosity and increase accessibility to hazardous trace elements such as arsenic. In this study, we developed a custom microfluidic cell for an in operando synchrotron experiment using X-ray attenuation. The experiment mimics reactive fracture flow by passing an acidic fluid over a surface of mineralogically heterogeneous rock from the Eagle Ford shale. Over 48 h, calcite was preferentially dissolved, forming an altered layer 200–500 μm thick with a porosity of 63–68% and surface area >10$\times$ higher than that in the unreacted shale as shown by xCT analyses. Calcite dissolution rate quantified from the attenuation data was 3 $\times$ 10 –4 mol/m 2 s and decreased to 3 $\times$ 10 –5 mol/m 2 s after 24 h because of increasing diffusion limitations. Erosion of the fracture surface increased access to iron-rich minerals, thereby increasing access to toxic metals such as arsenic. Quantification using XRF and XANES microspectroscopy indicated up to 0.5 wt % of As(-I) in arsenopyrite and 1.2 wt % of As(V) associated with ferrihydrite. This study provides valuable contributions for understanding and predicting fracture alteration and changes to the mobilization potential of hazardous metals and metalloids.

58 GEOSCIENCES↗

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↗

Patterns of fracture and tidal stresses due to nonsynchronous rotation - Implications for fracturing on Europa

This study considers the global patterns of fracture that would result from nonsynchronous rotation of a tidally distorted planetary body. The incremental horizontal stresses in a thin elastic or viscous shell due to a small displacement of the axis of maximum tidal elongation are derived, and the resulting stress distributions are applied to interpret the observed pattern of fracture lineaments on Europa. The observed pattern of lineaments can be explained by nonsynchronous rotation if these features formed by tension fracturing and dike emplacement. Tension fracturing can occur for a small displacement of the tidal axis, so that the resulting lineaments may be consistent with other evidence suggesting a young age for the surface.

Helfenstein, P.↗

Composite interlaminar fracture toughness: Three-dimensional finite element modeling for mixed mode 1, 2 and 3 fracture

A computational method/procedure is described which can be used to simulate individual and mixed mode interlaminar fracture progression in fiber composite laminates. Different combinations of Modes 1, 2, and 3 fracture are simulated by varying the crack location through the specimen thickness and by selecting appropriate unsymmetric laminate configurations. The contribution of each fracture mode to strain energy release rate is determined by the local crack closure methods while the mixed mode is determined by global variables. The strain energy release rates are plotted versus extending crack length, where slow crack growth, stable crack growth, and rapid crack growth regions are easily identified. Graphical results are presented to illustrate the effectiveness and versatility of the computational simulation for: (1) evaluating mixed-mode interlaminar fracture, (2) for identifying respective dominant parameters, and (3) for selecting possible simple test methods.

Murthy, P. L. N.↗

Composite interlaminar fracture toughness - Three-dimensional finite-element modeling for mixed mode I, II, and fracture

A computational method/procedure is described which can be used to simulate individual and mixed mode interlaminar fracture progression in fiber composite laminates. Different combinations of Modes 1, 2, and 3 fracture are simulated by varying the crack location through the specimen thickness and by selecting appropriate unsymmetric laminate configurations. The contribution of each fracture mode to strain energy release rate is determined by the local crack closure methods while the mixed mode is determined by global variables. The strain energy release rates are plotted versus extending crack length, where slow crack growth, stable crack growth, and rapid crack growth regions are easily identified. Graphical results are presented to illustrate the effectiveness and versatility of the computational simulation for: (1) evaluating mixed-mode interlaminar fracture, (2) for identifying respective dominant parameters, and (3) for selecting possible simple test methods.

Murthy, Pappu L. N.↗

Effect of initial delamination on Mode 1 and Mode 2 interlaminar fracture toughness and fatigue fracture threshold

Static and fatigue double-cantilever beam (DCB) and end-notch flexure (ENF) tests were conducted to determine the effect of the simulated initial delamination in interlaminar fracture toughness, G(sub c), and fatigue fracture threshold, G(sub th). Unidirectional, 24-ply specimens of S2/SP250 glass/epoxy were tested using Kapton inserts of four different thickness - 13, 25, 75, and 130 microns, at the midplane at one end, or with tension or shear precracks, to simulate an initial delamination. To determine G(sub c), the fatigue fracture threshold below which no delamination growth would occur in less than 1 x 10(exp 6) cycles, fatigue tests were conducted by cyclically loading specimens until delamination growth was detected. Consistent values of model 1 fracture toughness, G(sub Ic), were measured from DCB specimens with inserts of thickness 75 microns or thinner, or with shear precracks. The fatigue DCB tests gave similar values of G(sub Ith) for the 13, 25, and 75 microns specimens. Results for the shear precracked specimens were significantly lower that for specimens without precracks. Results for both the static and fatigue ENF tests showed that measured G(IIc) and G(IIth) values decreased with decreasing insert thickness, so that no limiting thickness could be determined. Results for specimens with inserts of 75 microns or thicker were significantly higher than the results for precracked specimens or specimens with 13 or 25 microns inserts.

Murri, Gretchen Bostaph↗

Prediction of Composite Laminate Fracture: Micromechanics and Progressive Fracture

This report describes an investigation to predict first-ply failure and final fracture in selected composite laminates subjected to inplane loads. The laminates were composed of glass fiber and graphite fibers in epoxy matrices. Failure envelopes based on first-ply failure and laminate fracture were generated for combined loading of these laminates. Predictions were evaluated by micromechanics-based theory and progressive fracture. The results show that, for most cases, combined tensile loading significantly enhanced the laminate fracture stress in comparison to the uniaxial loading.

Gotsis, P. K.↗

A fully coupled method for massively parallel simulation of hydraulically driven fractures in 3-dimensions: FULLY COUPLED PARALLEL SIMULATION OF HYDRAULIC FRACTURES IN 3-D

This study describes a fully coupled finite element/finite volume approach for simulating field-scale hydraulically driven fractures in three dimensions, using massively parallel computing platforms. The proposed method is capable of capturing realistic representations of local heterogeneities, layering and natural fracture networks in a reservoir. A detailed description of the numerical implementation is provided, along with numerical studies comparing the model with both analytical solutions and experimental results. The results demonstrate the effectiveness of the proposed method for modeling large-scale problems involving hydraulically driven fractures in three dimensions.

03 NATURAL GAS↗

Risk Prediction of New Adjacent Vertebral Fractures After PVP for Patients with Vertebral Compression Fractures: Development of a Prediction Model

PurposeWe aim to determine the predictors of new adjacent vertebral fractures (AVCFs) after percutaneous vertebroplasty (PVP) in patients with osteoporotic vertebral compression fractures (OVCFs) and to construct a risk prediction score to estimate a 2-year new AVCF risk-by-risk factor condition.Materials and MethodsPatients with OVCFs who underwent their first PVP between December 2006 and December 2013 at Hospital A (training cohort) and Hospital B (validation cohort) were included in this study. In training cohort, we assessed the independent risk predictors and developed the probability of new adjacent OVCFs (PNAV) score system using the Cox proportional hazard regression analysis. The accuracy of this system was then validated in both training and validation cohorts by concordance (c) statistic.Results421 patients (training cohort: n = 256; validation cohort: n = 165) were included in this study. In training cohort, new AVCFs after the first PVP treatment occurred in 33 (12.9%) patients. The independent risk factors were intradiscal cement leakage and preexisting old vertebral compression fracture(s). The estimated 2-year absolute risk of new AVCFs ranged from less than 4% in patients with neither independent risk factors to more than 45% in individuals with both factors.ConclusionsThe PNAV score is an objective and easy approach to predict the risk of new AVCFs.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Investigation of fracture in porous materials: a phase-field fracture study informed by ReaxFF

Microscopic features (e.g, pore shapes, sizes, and distribution) in porous material substantially affect the overall mechanical properties such as stiffness and strength. In turn, these material properties determine the macroscopic behaviors of fracture in the porous material. In certain cases, macroscopic properties can be derived from the porous skeleton and void ratio (i.e., porosity), but in many other cases, such derivation is a challenging task. This paper presents a numerical investigation of microporosity and micropore shapes effect on the macrofracture behavior in porous amorphous silica. For this study, we extend the recently-proposed combined molecular dynamic (MD) and phase-field (PF) fracture modeling approach by including different pore shapes in the atomistic domain. In the MD simulations, we adopt ReaxFF to evaluate the material properties, where four different micropore cases are considered. Based on the material properties derived from MD simulations, the macrofracture propagation of porous media is studied using hybrid PF simulation. In the characterization of the pore structure, the concept of pore ligament is proposed to relate the pore shape and the critical energy release rate. Two classical fracture problems were used to evaluate the effect of pore shape on the macrofracture behavior. The results of the case study show that although the micropore shapes change the macrofracture behaviors, these effects vary with the geometry and loading conditions of macroscopic boundary value problems. Furthermore, the case study also shows that the influence of micropore structure can be captured at the macroscopic level through the material properties derived from the MD simulations.

42 ENGINEERING↗

Application of Koopman operator for model-based control of fracture propagation and proppant transport in hydraulic fracturing operation

This work explores the application of the recently developed Koopman operator approach for model identification and feedback control of a hydraulic fracturing process. Controlling fracture propagation and proppant transport with precision is a challenge due in large part to the difficulty of constructing approximate models that accurately capture the characteristic moving boundary and highly-coupled dynamics exhibited by the process. Koopman operator theory is particularly attractive here as it offers a way to explicitly construct linear representations for even highly nonlinear dynamics. The method is data-driven and relies on lifting the states to an infinite-dimensional space of functions called observables where the dynamics are governed by a linear Koopman operator. Here this work considers two problems: (a) fracture geometry control, and (b) proppant concentration control. In both cases, an approximate linear model of the corresponding dynamics is constructed and used to design a model predictive controller (MPC). The manuscript shows that in the case of highly nonlinear dynamics, as observed in the proppant concentration, use of canonical functions in the observable basis fails. In such cases, a priori system knowledge can be leveraged to choose the required basis. The numerical experiments demonstrate that the Koopman linear model shows excellent agreement with the real system and successfully achieves the desired target values maximizing the oil and gas productivity. Additionally, due to its linear structure, the Koopman models allow convex MPC formulations that avoid any issues associated with nonlinear optimization.

42 ENGINEERING↗

Micromechanics of Fracture Propagation During Multistage Stress Relaxation and Creep in Brittle Rocks

Time-dependent rock deformation caused by the initiation and growth of fractures leads to the weakening of the rock mass. Understanding the fracturing mechanisms involved in the time-dependent behavior in brittle rocks is very important and to achieve this goal, a systematic series of three types of experiments was performed on double-flawed prismatic Barre granite specimens under unconfined compression. The first series aimed to identify the failure mechanism in the short-term failure mode under monotonic loading, the, second series involved multistage relaxation (constant strain) experiments to analyze the damage at different strain levels, and the third series explored the fracture propagation under multistage creep (constant load) experiments. The spatial and temporal evolution of cracking mechanisms were evaluated using the acoustic emission (AE) and two-dimensional digital image correlation (2D-DIC) techniques to observe the whole crack growth process as well as the accumulated inelastic strain at the specified region of interest. Results suggest that in the case of multistage creep experiments, the time to failure was less compared to the multistage relaxation, when loaded above the crack damage threshold (CD) estimated from the monotonic testing. Furthermore, the frequency magnitude distribution of the AE events generated in the three loading conditions followed the Gutenberg Richter model. A relatively lower b-value was obtained for the creep experiments, indicative of high energy AE events and faster crack growth. In addition, the AE and DIC results also revealed high evolution of tensile cracks at-different stages of creep and relaxation compared to shear and mixed-mode cracks.

58 GEOSCIENCES↗

An Investigation of CO2-Responsive Preformed Particle Gel for Conformance Control of CO2 Flooding in Reservoirs With Fractures or Fracture-Like Channels

Summary Gel treatment is an important technique to solve early CO2 breakthrough and excess–CO2–production problems, caused by the low viscosity and low density of CO2, as well as the heterogeneity of reservoirs with fractures or fracture–like channels. However, there is no reported work on gel that increases its volume after reacting with CO2 (termed CO2–responsive gel) for the conformance control of CO2 flooding. In this paper, the intrinsic properties of a CO2–responsive preformed particle gel (CR–PPG) were evaluated in a water/supercritical–CO2 (scCO2) environment in high–pressure vessels. Continuous scCO2 injection and CR–PPG treatment were conducted in fractured sandstone cores, to probe their plugging performance to scCO2 flow in a high–permeability–contrast system. The volumetric swelling ratio (VSR) of the CR–PPG increased by approximately two times in the presence of scCO2, compared with a sample under similar conditions in the absence of scCO2. The CR–PPG swelling ratio decreased with increasing NaCl concentration. Under the same conditions, the temperature did not have an apparent effect on the swelling ratio after 31 days of swelling. In coreflooding experiments, the placed CR–PPG resisted a considerable pressure up to 617.0 psi before breakthrough. After a shut–in process, CO2–breakthrough pressure was detected at 437.2 psi. It is observed that the shut–in process improved the plugging performance of CR–PPG to CO2 as revealed by the increase in the residual resistance factor. Controlling the shut–in time was found to be effective in augmenting the increase in the residual resistance factor, by increasing the VSR of placed CR–PPG. Resistance of CR–PPG to some real field challenges, including a high pressure gradient and long–term exposure to CO2, was also reported for field–applicability concerns.

Engineering↗