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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Foot-scale Evaluation of CO2-Responsive Polymer and CO2 Binary Fluid as an Alternative Fracturing Fluid for Enhanced Geothermal Systems

A foot-scale high-temperature true-triaxial fracturing apparatus was used to conduct hydraulic fracture tests on rock samples. The fracturing tests used different fluids including water, CO2, CO2 with water, and CO2 with a 1% polyallylamine (PAA)/water solution. The tests were conducted with both constant pressure injection and constant flow rate injection modes. CO2-based fracturing fluids were found to produce higher breakdown pressures, high transient flow rates, and produce higher-conductivity fractures as compared to water-based fracturing fluids. Additionally, faster pressurization rates with CO2-based fracturing fluids (obtained when fracturing at constant flow rate mode) are found to be associated with higher fracture conductivities. When fracturing with CO2-based fluids in the presence of a PAA aqueous solution, for example if the rock is saturated with 1wt% PAA in water, the volume expansion caused by CO2-induced crosslinking of PAA leads to a faster pressure increase due to the associated volume expansion and increase in viscosity. It was also found that CO2 as a fracturing fluid injected in hot dry rock (HDR) attain the highest fracture conductivity only when injected at very high flow rates, followed very closely by the CO2/PAA fracturing fluid system that generates fractures with, on average, similarly high conductivity values though independently of injection flow rate and using 1/6 of the mass of CO2 as compared to CO2 in HDR. Breakdown pressures were also similar for CO2 stimulation in HDR and CO2/PAA fluid system under identical injection flow rates. Finally, the well-known low viscosity of CO2 phase prevents the efficient transport of proppants while the reacted CO2/PAA fluid system is known to form a high viscosity binary fracturing fluid system with the potential to carry proppant in addition of being a fine fracturing fluid.

Enhanced Geothermal Systems, polyallylamine, Co2-r↗
DOE OSTI2019

Spatial Simulation and Characterization of Three-Dimensional Fractures in Gejiu tin District, Southwest China, Using GEOFRAC

Fracture is an important factor controlling mineralization and ore distribution in metallic deposits. Fracture also affects mechanical stability of drifts. A plausible three-dimensional discrete fracture network (DFN) was constructed using GEOFRAC, a geostatistical method of conditioning directions (strikes and dips) and locations of sample fractures. The northern part of the Gejiu tin district, southwest China, was selected, using 10,212 fractures sampled on drift roofs at five different levels and along a 5463 m length. Key parameters of this DFN are fracture density, direction, and the connective condition of disks that form a fracture plane. The simulated fractures were verified by good correspondence of their directions with those of the sample fractures; furthermore, continuous fractures longer than 1 km matched the faults observed in outcrops. The most notable result was clarification of fracture control by the orebody shape: Simulated gentle fractures (30° dips or less) formed layered orebodies and created boundaries between each layer, while steep fractures (60° dips or more) displaced the layers by acting as faults. Such fault-type fractures and the densely fractured portions require caution with respect to mining safety. It was concluded that GEOFRAC is useful for three-dimensional DFN modeling in underground mines.

58 GEOSCIENCES↗
DOE OSTI2016

A practical model for fluid flow in discrete-fracture porous media by using the numerical manifold method

Here in this study, a numerical manifold method (NMM) model is developed to analyze flow in porous media with discrete fractures in a non-conforming mesh. This new model is based on a two-cover-mesh system with a uniform triangular mathematical mesh and boundary/fracture-divided physical covers, where local independent cover functions are defined. The overlapping parts of the physical covers are elements where the global approximation is defined by the weighted average of the physical cover functions. The mesh is generated by a tree-cutting algorithm. A new model that does not introduce additional degrees of freedom (DOF) for fractures was developed for fluid flow in fractures. The fracture surfaces that belong to different physical covers are used to represent fracture flow in the direction of the fractures. In the direction normal to the fractures, the fracture surfaces are regarded as Dirichlet boundaries to exchange fluxes with the rock matrix. Furthermore, fractures that intersect with Dirichlet or Neumann boundaries are considered. Simulation examples are designed to verify the efficiency of the tree-cutting algorithm, the calculation's independency from the mesh orientation, and accuracy when modeling porous media that contain fractures with multiple intersections and different orientations. The simulation results show good agreement with available analytical solutions. Finally, the model is applied to cases that involve nine intersecting fractures and a complex network of 100 fractures, both of which achieve reasonable results. The new model is very practical for modeling flow in fractured porous media, even for a geometrically complex fracture network with large hydraulic conductivity contrasts between fractures and the matrix.

58 GEOSCIENCES↗
DOE OSTI2017

A numerical manifold method model for analyzing fully coupled hydro-mechanical processes in porous rock masses with discrete fractures

A numerical manifold method (NMM) model was developed for fully coupled analysis of hydro-mechanical (HM) processes in porous rock masses with discrete fractures. Using an NMM two-cover-mesh system of mathematical and physical covers, fractures are conveniently discretized by dividing the mathematical cover along fracture traces to physical cover, resulting in a discontinuous model on a non-conforming mesh. In this model, discrete fracture deformation (e.g. open and slip) and fracture fluid flow within a permeable and deformable porous rock matrix are rigorously considered. For porous rock, direct pore-volume coupling was modeled based on an energy-work scheme. For mechanical analysis of fractures, a fracture constitutive model for mechanically open states was introduced. For fluid flow in fractures, both along-fracture and normal-to-fracture fluid flow are modeled without introducing additional degrees of freedom. When the mechanical aperture of a fracture is changing, its hydraulic aperture and hydraulic conductivity is updated. At the same time, under the effect of coupled deformation and fluid flow, the contact state may dynamically change, and the corresponding contact constraint is updated each time step. Therefore, indirect coupling is realized under stringent considerations of coupled HM effects and fracture constitutive behavior transfer dynamically. To verify the new model, examples involving deformable porous media containing a single and two sets of fractures were designed, showing good accuracy. Finally, the model was applied to analyze coupled HM behavior of fractured porous rock domains with complex fracture networks under effects of loading and injection.

58 GEOSCIENCES↗
DOE OSTI2013

Modifications of Carbonate Fracture Hydrodynamic Properties by CO 2 -Acidified Brine Flow

Acidic reactive flow in fractures is relevant in subsurface activities such as CO{sub 2} geological storage and hydraulic fracturing. Understanding reaction-induced changes in fracture hydrodynamic properties is essential for predicting subsurface flows such as leakage, injectability, and fluid production. In this study, x-ray computed tomography scans of a fractured carbonate caprock were used to create three dimensional reconstructions of the fracture before and after reaction with CO{sub 2}-acidified brine (Ellis et al., 2011, Greenhouse Gases: Sci. Technol., 1:248-260). As expected, mechanical apertures were found to increase substantially, doubling and even tripling in some places. However, the surface geometry evolved in complex ways including ‘comb-tooth’ structures created from preferential dissolution of calcite in transverse sedimentary bands, and the creation of degraded zones, i.e. porous calcite-depleted areas on reacted fracture surfaces. These geometric alterations resulted in increased fracture roughness, as measured by surface Z{sub 2} parameters and fractal dimensions D{sub f}. Computational fluid dynamics (CFD) simulations were conducted to quantify the changes in hydraulic aperture, fracture transmissivity and permeability. The results show that the effective hydraulic apertures are smaller than the mechanical apertures, and the changes in hydraulic apertures are nonlinear. Overestimation of flow rate by a factor of two or more would be introduced if fracture hydrodynamic properties were based on mechanical apertures, or if hydraulic aperture is assumed to change proportionally with mechanical aperture. The differences can be attributed, in part, to the increase in roughness after reaction, and is likely affected by contiguous transverse sedimentary features. Hydraulic apertures estimated by the 1D statistical model and 2D local cubic law (LCL) model are consistently larger than those calculated from the CFD simulations. In addition, a novel ternary segmentation method was devised to handle the degraded zones, allowing for a bounding analysis of the effects on hydraulic properties. We found that the degraded zones account for less than 15% of the fracture volume, but cover 70% to 80% of the fracture surface. When the degraded zones are treated as part of the fracture, the fracture transmissivities are two to four times larger because the fracture surfaces after reaction are not as rough as they would be if one considers the degraded zone as part of the rock. Therefore, while degraded zones created during geochemical reactions may not significantly increase mechanical aperture, this type of feature cannot be ignored and should be treated with prudence when predicting fracture hydrodynamic properties.

roughness↗
DOE OSTI2021

Propagation of Cryogenic Thermal Fractures from Unconfined PMMA Boreholes

In cryogenic fracturing, a rock surface exposed to cryogenic fluids undergoes a large thermal gradient, and the resultant local tensile stress overcomes rock strength and initiates fractures. This study investigates the development of cracks generated from the cryogenic treatment of a borehole under no external confining stress on specimens. The experiments were performed on transparent PMMA specimens to observe fracture proliferation around boreholes. Liquid nitrogen was flowed through the boreholes to cool the borehole surface. The results show that initial fracture growth is characterized by abrupt starts and stops, and as the fracture propagates outward, the growth appears more continuous. In an early stage, horizontal/radial fractures and vertical fractures are the defining patterns. Horizontal fractures tend to be separated by a specific exclusion distance (i.e., spacing between cracks). While distinct horizontal/vertical fractures and exclusion distance manifest themselves at an early stage, fractures resulting from fracture interactions and curvatures can develop into complex shapes at later stages. Cryogenic thermal loading induces distinctively curved fractures. The tendency of curvature may prevent greater penetration. An increase in the borehole pressure during liquid nitrogen flow, however, can lessen fracture tortuosity and facilitate radial propagation. A high flow pressure and rate are also advantageous in that they accelerate cooling and fracture propagation.

42 ENGINEERING↗
DOE OSTI2022

Scaling Behavior of Thermally Driven Fractures in Deep Low‐Permeability Formations: A Plane Strain Model With 1‐D Heat Conduction

Abstract Injection of cold fluids through/into deep formations may cause significant cooling, thermal stress, and possible thermal fracturing. In this study, the thermal fracturing of low‐permeability formations under one‐dimensional heat conduction was investigated using a plane strain model. Dimensionless governing equations, with dimensionless fracture length , aperture , spacing , time , and effective confining stress , were derived. Solution of single thermal fracture was derived analytically, while solution of multiple fractures with constant (or dynamic) spacing were obtained using the displacement discontinuity method (and stability analysis). For single fracture, increases nonlinearly with and then transitions to scaling law , indicating that late‐time fracture length increases linearly with the square root of cooling time. For constantly spaced fractures, deviates from the single‐fracture solution at a later for a larger , showing slower propagation under inter‐fracture stress interaction . For dynamically spaced fractures, fracture arrest induced by stress interaction was determined by the stability analysis; the fully transient solution provides evolution of dimensionless fracture length, spacing, aperture, and pattern; a similar scaling law, with , obtained shows the effect of both stress interaction and fracture arrest. The solution and scaling law provide fast predictions for all reservoir and cooling conditions using (single) model parameter . Application to a geothermal site with demonstrates that thermal fractures reach 0.67, 6.25, and 78.00 m in length, 0.49, 2.30, and 13.00 m in spacing, and 0.43, 2.09, and 12.19 mm in aperture at 1, 100, and 10,000 days.

Chen, Bin↗
DOE OSTI2025

Continuous surface-to-distributed acoustic sensor snapshots explain reactivation of individual natural fractures during an unconventional reservoir stimulation

ABSTRACT Fiber-optic sensing technologies allow petroleum engineering teams to detect hydraulic fracture interaction with boreholes during unconventional reservoir stimulation. In combination with high-repeatability seismic sources, the same distributed acoustic sensors (DASs) enable vertical seismic profiling (VSP) of the fracture evolution away from the boreholes. We discovered clear signatures of seismic scattering on activated fractures during nine days of continuous seismic monitoring of the fracturing stages at the Austin Chalk/Eagle Ford Field Laboratory. The present study applies a novel approach for quantitative analysis of the scattering events in terms of the evolution of the geometry and elastic stiffness of individual fractures. Our characterization strategy sequentially refines the fracture models: from a stack of 1D soft layers to 3D rectangular inclusions. First, we estimate the number of fracture locations and reflectivity using a modified sparse-spike deconvolution of the stacked VSP traces. The fracture set consists of five fractures spaced by 15–30 m with a reflectivity of approximately 1%. Then, we develop a scattering integral method to refine these estimates along with an inversion of the fracture top and bottom for each monitoring vintage. We find that, initially, some of the fractures are located above the monitoring fiber with the height of approximately 100 m. Then we integrate the seismic interpretation with the low-frequency DAS and pressure and microseismic monitoring to reconstruct the activation process of the fractures. Most likely, some of the natural fractures slowly grew downward to the monitoring fiber as a result of fluid injections in the stimulated well. This led to bright strain anomalies but did not trigger seismicity. The top of the fractures remained almost constant and were limited by a lithologic boundary/stress barrier. To our knowledge, this is the first time VSP data enabled tracking of the fracture evolution with such high spatial and temporal resolution, which was previously only available for crosswell surveys and at a much smaller scale.

Glubokovskikh, Stanislav↗
DOE OSTI2019

Relationships between lineal fracture intensity and chemical composition in the Marcellus Shale, Appalachian Basin

Within mudrock reservoirs, brittle zones undergo failure during hydraulic stimulation, creating numerous artificial fractures which enable hydrocarbons to be liberated from the reservoir. Natural fractures in mudrock reduce the tensile strength of the host rock, creating planes of weaknesses that are hypothesized to be reactivated during hydraulic stimulation. Combined, brittleness and natural fractures contribute to creating more abundant and complex fracture networks during hydraulic stimulation. Research efforts toward quantifying rock brittleness have resulted in numerous mineral-/compositional-based indices, which are used during petrophysical analysis to predict zones most conducive to hydraulic stimulation. In contrast, investigations on the relationship between chemical composition and core-scale natural fractures are limited. For this study, we collected high-resolution energy-dispersive X-ray fluorescence (XRF) data, calibrated with a wave-dispersive XRF, from a Marcellus Shale core. Additionally, we characterized corescale natural fractures in terms of length, width, in-filling material or lack thereof, and orientation. Following the characterization, we transformed the natural fracture data into a continuous P10 (lineal fracture intensity) curve, expressed as the number of fractures per a one-half foot window. Using these data sets, we investigated the relationship between rock composition and natural fracture intensity. Regression analyses recorded positive relationships between natural fracture intensity and calcium, silicon/aluminum, and total organic carbon (TOC), and negative relationships with silicon and aluminum. Aluminum recorded the strongest (negative) relationship ([Formula: see text]) with natural fracture intensity. To access the degree to which natural fractures can be predicted based on chemical composition, we applied a partial least-squares analysis, a multivariate method, and recorded an [Formula: see text]. Our study illustrates that although numerous factors are responsible for natural fracture genesis, such fractures predictively concentrate in areas of similar chemical composition, largely in zones with low aluminum concentrations.

Geochemistry & Geophysics↗
DOE OSTI2021

Fracture toughness evaluations for spent nuclear fuel dry storage canister welds and spent nuclear fuel clad-pellet structures

Integrity of spent nuclear fuel (SNF) interim storage canisters is very important to the safety of the back-end nuclear fuel cycle. Stress corrosion cracking (SCC) potential of interim storage canister has been considered as a high priority. Because no post-weld heat treatment was required for forming these canisters, the high tensile residual stress existed within these canister welds. This can change the fracture resistance capacity significantly as well as increase SCC potential. Due to relative thin shell thickness of a canister weldment, the spiral notch torsion test (SNTT) method was used to estimate the canister weldment fracture toughness. SNTT was developed to measure the intrinsic fracture toughness of structural materials using small specimens. The SNTT method has been applied to a wide variety of structural materials, such as low-alloy steels, stainless steel, aluminum alloy, ceramics, concrete, and composites. The SNTT system operates by applying pure torsion to cylindrical specimens with a notch line that spirals around the specimen at a 45° pitch. In order to carry out pure torsion load mode, biaxial tension/torsion tester was developed accordingly to perform SNTT protocol. Moreover, applying fracture mechanics approach to SNF system reliability investigation is warranted due to the inherent flaws and hydride structures existed in a SNF system after nuclear reactor operation. However, none of the existing fracture toughness data deal with fuel cladding specific geometry or spent fuel material conditions, such as cladding structure with the pellet-inserts and the associated pellet clad mechanical interactions induced mixed-mode damage mechanisms. Thus, the development of an intrinsic fracture mechanics approach that is suitable for SNF materials is needed. Furthermore, due to thin wall and small dimension of clad tubing structure, the SNTT method was used to estimate the clad tubing structure fracture toughness. Fracture testing were performed on the received stainless steel canister weldment, most SNTT weld samples fracture initiation sites are at heat-affected zone (HAZ) regions. The estimated fracture toughness J Q ’ for the baseline SS304 steel is at 283 kJ/m². The estimated JQ’ for the SS304/308 weld and baseline metals are 148 kJ/m 2 and 459 kJ/m 2 , respectively. Out of cell fracture testing for spent fuel structure were carried out on the surrogate rods made of Zr-4 clad and alumina inserts, the estimated fracture toughness values for baseline Zr-4 cladding with alumina-pellet inserts are: (1) For SNTT samples with a short or medium crack length, between 5.4-mm and 8-mm, the estimate J IQ upon fracture initiation for the baseline Zr-4 cladding is at 50 kJ/m 2 with 2-sigma uncertainty of 3.26 kJ/m 2 , and the associated K IQ is at 67.46 MPa$\sqrt{m}$; and (2) For SNTT samples with a long crack length, around 13-mm, the crack initialization is deviated from that of the Mode-I tensile fracture and appears to be a mixed-mode fracture of Mode I - tensile stress and Mode III - out of plane shear stress; the estimated J MQ is at 18.9 kJ/m 2 , the associated K MQ is at 41.4 MPa$\sqrt{m}$.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Experimental Investigation of Low-Frequency Distributed Acoustic Sensor Responses to Two Parallel Propagating Fractures

Low-frequency distributed acoustic sensing (LF-DAS) is a diagnostic tool for hydraulic fracture propagation with far-field monitoring using fiber optic sensors. LF-DAS senses strain rate variation caused by stress field change due to fracture propagation. Fiber optic sensors are installed in the monitoring wells in the vicinity of a fractured well. From the strain responses, fracture propagation can be evaluated. To understand subsurface conditions with multiple propagating fractures, a laboratory-scale hydraulic fracture experiment was performed simulating the LF-DAS response to fracture propagation with embedded distributed optical fiber strain sensors under these conditions. The experiment was performed using a transparent cube of epoxy with two parallel radial initial flaws centered in the cube. Fluid was injected into the sample to generate fractures along the initial flaws. The experiment used distributed high-definition fiber optic strain sensors with tight spatial resolutions. The sensors were embedded at two different locations on opposite sides of the initial flaws, serving as observation/monitoring locations. We also employed finite element modeling to numerically solve the linear elastic equations of equilibrium continuity and stress–strain relationships. The measured strains from the experiment were compared to simulation results from the finite element model. The experimentally derived strain and strain-rate waterfall plots from this study show the responses to both fractures propagating, while the fracture at the lower position took most of the fluid during the experiment. Interestingly, a fracture first began propagating from the upper flaw of the two flaws, but once the lower fracture was initiated, it grew much faster than the upper fracture. Both fibers were intercepted by the lower fracture, further verifying the strain signature as a fracture is approaching and intersecting an offset fiber.

Chemistry↗
DOE OSTI2023

Multiscale spatial analysis of fracture nodes in two dimensions

Spatial arrangement of fractures as a function of scale is an important component of fracture quantification for inferential and predictive modeling. Available methods that analyze fracture spatial arrangement are based on one-dimensional spacing data; therefore, they are limited to semi-parallel fractures. Such methods cannot be applied to fracture networks in higher dimensions, particularly when fractures have different orientations. Here to characterize fracture arrangements in two dimensions, we propose using Ripley’s K-function, as a method of point pattern analysis, to quantify spatial arrangement of fracture nodes. Fracture nodes, such as barycenters, intersection points, and tips, are point-based representations of fracture locations and connectivity within the fracture network. We introduce formulations for isotropic as well as directional analyses of spatial arrangement. In addition, we derive formulations for edge correction in circular and rectangular study domains. Finally, we demonstrate applications of Ripley’s K-function on two natural fracture datasets. Our proposed method supports quantification and characterization of fracture spatial arrangements that allow practitioners to build representative models of fractures in the subsurface.

02 PETROLEUM↗
DOE OSTI2021

Predicting Fracture Network Development in Crystalline Rocks

The geometric properties of fractures influence whether they propagate, arrest, or coalesce with other fractures. Thus, quantifying the relationship between fracture network characteristics may help predict fracture network development, and perhaps precursors to catastrophic failure. To constrain the relationship and predictability of fracture characteristics, we deform eight one centimeter tall rock cores under triaxial compression while acquiring in situ X-ray tomograms. The tomograms reveal precise measurements of the fracture network characteristics above the spatial resolution of 6.5 µm. We develop machine learning models to predict the value of each characteristic using the other characteristics, and excluding the macroscopic stress or strain imposed on the rock. The models predict fracture development more accurately in the experiments performed on granite and monzonite, than the experiments on marble. Fracture network development may be more predictable in these igneous rocks because their microstructure is more mechanically homogeneous than the marble, producing more systematic fracture development that is not strongly impeded by grain contacts and cleavage planes. The varying performance of the models suggest that fracture volume, length, and aperture are the most predictable of the characteristics, while fracture orientation is the least predictable. Orientation does not correlate with length, as suggested by the idea that the orientation evolves with increasing differential stress and thus fracture length. This difference between the observed and expected relationship between orientation and length highlights the influence of mechanical heterogeneities and local stress perturbations on fracture growth as fractures propagate, link, and coalesce.

58 GEOSCIENCES↗
DOE OSTI2015

Stimuli-Responsive/Rheoreversible Hydraulic Fracturing Fluids as a Greener Alternative to Support Geothermal and Fossil Energy Production

Cost-effective yet safe creation of high-permeability reservoirs within deep bedrock is the primary challenge for the viability of enhanced geothermal systems (EGS) and unconventional oil/gas recovery. Although fracturing fluids are commonly used for oil/gas, standard fracturing methods are not developed or proven for EGS temperatures and pressures. Furthermore, the environmental impacts of currently used fracturing methods are only recently being determined. Widespread concerns about the environmental contamination have resulted in a number of regulations for fracturing fluids advocating for greener fracturing processes. To enable EGS feasibility and lessen environmental impact of reservoir stimulation, an environmentally benign, CO2-activated, rheoreversible fracturing fluid that enhances permeability through fracturing (at significantly lower effective stress than standard fracturing fluids) due to in situ volume expansion and gel formation is investigated herein. The chemical mechanism, stability, phase-change behavior, and rheology for a novel polyallylamine (PAA)-CO2 fracturing fluid was characterized at EGS temperatures and pressures. Hydrogel is formed upon reaction with CO2 and this process is reversible (via CO2 depressurization or solubilizing with a mild acid) allowing removal from the formation and recycling, decreasing environmental impact. Rock obtained from the Coso geothermal field was fractured in laboratory experiments under various EGS temperatures and pressures with comparison to standard fracturing fluids, and the fractures were characterized with imaging, permeability measurement, and flow modeling. This novel fracturing fluid and process may vastly reduce water usage and the environmental impact of fracturing practices and effectively make EGS production and unconventional oil/gas exploitation cost-effective and cleaner.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2014

Fractured rock stress‐permeability relationships from in situ data and effects of temperature and chemical‐mechanical couplings

Abstract The purpose of this paper is to (i) review field data on stress‐induced permeability changes in fractured rock; (ii) describe estimation of fractured rock stress‐permeability relationships through model calibration against such field data; and (iii) discuss observations of temperature and chemically mediated fracture closure and its effect on fractured rock permeability. The field data that are reviewed include in situ block experiments, excavation‐induced changes in permeability around tunnels, borehole injection experiments, depth (and stress) dependent permeability, and permeability changes associated with a large‐scale rock‐mass heating experiment. Data show how the stress‐permeability relationship of fractured rock very much depends on local in situ conditions, such as fracture shear offset and fracture infilling by mineral precipitation. Field and laboratory experiments involving temperature have shown significant temperature‐driven fracture closure even under constant stress. Such temperature‐driven fracture closure has been described as thermal overclosure and relates to better fitting of opposing fracture surfaces at high temperatures, or is attributed to chemically mediated fracture closure related to pressure solution (and compaction) of stressed fracture surface asperities. Back‐calculated stress‐permeability relationships from field data may implicitly account for such effects, but the relative contribution of purely thermal‐mechanical and chemically mediated changes is difficult to isolate. Therefore, it is concluded that further laboratory and in situ experiments are needed to increase the knowledge of the true mechanisms behind thermally driven fracture closure, and to further assess the importance of chemical‐mechanical coupling for the long‐term evolution of fractured rock permeability.

58 GEOSCIENCES↗
DOE OSTI2013

Evaluation of known-boundary and resistivity constraints for improving cross-borehole DC electrical resistivity imaging of discrete fractures

There is a need to better characterize discrete fractures in contaminated hard rock aquifers to determine the fate of remediation injections away from boreholes and also to evaluate hydraulic fracturing performance. A synthetic cross borehole electrical resistivity study was conducted assuming a discrete fracture model of an existing contaminated site with known fracture locations. Four boreholes and two discrete fracture zones, assumed to be the dominant electrical and hydraulically conductive pathways, were explicitly modeled within an unstructured tetrahedral mesh. We first evaluated different regularization constraints starting with an uninformed smoothness-constrained inversion, to which a priori information was incrementally added. We found major improvements when (1) smoothness regularization constraints were relaxed (or disconnected) along boreholes and fractures, (2) a homogeneous conductivity was assumed along boreholes, and (3) borehole conductivity constraints that could be determined from a specific conductance log were applied. We also evaluated the effect of including borehole packers on fracture zone model recovery. We found the fracture zone conductivities with the inclusion of packers were comparable to similar trials excluding the use of packers regardless of electrical potential changes. The misplacement of fracture regularization disconnects (FRDs) can easily be misinterpreted as actual fracture locations. Conductivities within these misplaced disconnects were near the starting model value and removing smoothing between boreholes and assumed fracture locations helped in identifying incorrectly located FRDs. Analysis confirmed improved model sensitivity when regularization disconnects were (1) applied along the boreholes and fracture zones, and (2) FRDs were placed where actual fractures existed. Our findings demonstrate that structural constraints used after careful evaluation of a priori information are critical to improve imaging of fracture electrical conductivities, locations and orientations.

geophysics, inversion, resistivity, fracture, cont↗
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

Deterministic Discrete Fracture Network (DFN) Model for the EGS Collab Project on the 4850 Level of the Sanford Underground Research Facility (SURF)

The EGS Collab is conducting hydraulic fracture stimulation and fluid circulation experiments in the Sanford Underground Research Facility (SURF) located in Lead, South Dakota. A total of eight ~60m-long subhorizontal boreholes were drilled from the 4850 Level (~1.5 km below the ground surface) into the crystalline rock of this former mine. Six of these holes are used for geophysical monitoring, one is used for hydraulic fracture stimulation, and the remaining hole was designed as a production borehole that receives water from the injection well via the induced and natural fracture system. The primary goal of creating the discrete fracture network model is to show that these modeling methods are critical for the development of enhanced geothermal systems (EGS). This includes the prediction of rock behavior during fracturing and during an extended period of water flow between the parallel injection and production boreholes. Understanding the results from the induced fracturing and flow is complicated by the presence of significant natural fractures that interact with the stimulation and/or flow pathways. The delineation and characterization of natural fractures is thus an important part of the project, and therefore a model of the Discrete Fracture Network (DFN) was developed on a deterministic basis. The DFN was populated using observations and interpretations integrated from drift (horizontal passageways that allow access in the underground) fracture mapping, analysis of core recovered from the eight boreholes, borehole televiewer logs and videos, and observations of flow between and within boreholes and in the drift. The natural fracture system is dominated by a pervasive northwest-trending, steeply dipping shear system that is identifiable in the drifts and the core. Hydraulic fracture stimulation, flow/tracer circulation tests, and geophysical monitoring revealed that the behavior of the injected water, and perhaps the growth of induced fractures, has been significantly influenced by the existing fractures identified in the DFN.

40 EE - Geothermal Technologies Office (EE-4G)↗