SEARCH · Engineering Papers

Results for “FRACTURES”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 379 records · Page 21

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Distinct Fracture Mineralogy That is Out of Equilibrium With Modern Groundwaters Provides Important Context for Subsurface Life

Rock fracture surfaces in the crust are essential habitat for microorganisms. Fracture‐groundwater interfaces provide physical substrates for biofilm growth and are sources of carbon, nutrients, and electron donors and acceptors. To better understand geochemical processes impacting fracture surfaces and the subsurface microbiome, we identified fractures in archived rock cores from the Soudan formation, which is known to host saline groundwaters and isolated microbial communities dependent on rock‐water interactions. Cores with open fractures were thin sectioned and studied via electron microprobe and synchrotron X‐ray fluorescence microprobe. Most fracture surfaces had mineralogy distinct from that of the bulk rock. Chlorite minerals were abundant on fracture surfaces and had elemental compositions suggesting deposition during late‐stage hydrothermal alteration. Fracture‐lining chlorites likely limit access to iron oxide and sulfide minerals that are active in subsurface biogeochemical cycles. Calcium‐rich rinds were also observed along fracture edges. These rinds were too thin and poorly ordered to be identified via light microscopy or X‐ray diffraction; however, Ca K‐edge micro‐X‐ray absorption near‐edge structure spectroscopy identified them as carbonates, minerals not observed in the bulk rock. Thermodynamic modeling shows that carbonate precipitation is largely unfavorable in Soudan groundwaters, indicating that fracture edge conditions differed from those in modern water samples. Because of the low carbon concentrations in Soudan groundwaters, carbonate rinds likely play an important role in subsurface carbon cycling and may mark fracture surfaces that once hosted biofilms. Overall, this study suggests that fracture alteration can both play an active role in and suppress rock‐water interactions essential to subsurface life.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Classification and Localization of Fracture-Hit Events in Low-Frequency Distributed Acoustic Sensing Strain Rate with Convolutional Neural Networks

Summary Distributed acoustic sensing (DAS) has been used in the oil and gas industry as an advanced technology for surveillance and diagnostics. Operators use DAS to monitor hydraulic fracturing activities, examine well stimulation efficacy, and estimate complex fracture system geometries. Particularly, low-frequency DAS can detect geomechanical events such as fracture hits because hydraulic fractures propagate and create strain rate variations in the rock. Analysis of DAS data today is mostly done post-job and subject to interpretation methods. However, the continuous and dense data stream generated live by DAS poses the opportunity for more efficient and accurate real-time data-driven analysis. The objective of this study is to develop a machine learning-based workflow that can identify and locate fracture-hit events in simulated strain rate responses correlated with low-frequency DAS data. In this paper, “fracture hit” refers to a hydraulic fracture originating from a stimulated well intersecting an offset well. We start with building a single fracture propagation model to produce strain rate patterns observed at a hypothetical monitoring well. This model is used to generate two sets of strain rate responses with one set containing fracture-hit events. The labeled synthetic data are then used to train a custom convolutional neural network (CNN) model for identifying the presence of fracture-hit events. The same model is trained again for locating the event with the output layer of the model replaced with linear units. We achieved near-perfect predictions for both event classification and localization. These promising results prove the feasibility of using CNN for real-time event detection from fiber-optic sensing data. Additionally, we use edge detection techniques to recognize fracture-hit event patterns in strain rate images. The fracture-hit location can be identified using recognized pixels in the image. The accuracy of edge detection-based location identification is also plausible, but edge detection is dependent on the assumption of pattern shape and image quality, hence it is less robust compared to CNN models. This comparison further supports the need for CNN applications in image-based real-time fiber-optic sensing event detection.

Engineering↗

DOE OSTI2023

Interpretation of Fracture Initiation Points by In-Well Low-Frequency Distributed Acoustic Sensing in Horizontal Wells

Summary Low-frequency distributed acoustic sensing (LF-DAS) exploits the optical phase shift of Rayleigh backscatter in fiber-optic cables to obtain distributed measurements of changes in strain and temperature. Fiber-optic cables are often installed for multistage hydraulic fracture diagnostics in horizontal wells. LF-DAS in an untreated well provides far-field strain measurements, while offset wells are hydraulically fractured. Such a configuration is called crosswell LF-DAS sensing. Crosswell LF-DAS measurements have proved useful in diagnosing fracture hits, fracture azimuth, planarity, cluster efficiency, fracture propagation rates, and the dynamic distance to the fracture front. In contrast, in-well LF-DAS is conducted on the actively fractured well. Due to cool fracture fluid being injected at high injection rates, the strain component of the LF-DAS response is largely obscured by temperature changes. In permanent fiber-optic cable installations, distributed temperature sensing (DTS) is often conducted simultaneously with LF-DAS. An opportunity exists to decouple the temperature and strain components of LF-DAS sensors to observe strain changes on in-well LF-DAS. The LF-DAS response is modeled as linearly dependent on strain and temperature changes. Theoretical LF-DAS temperature and strain sensitivity coefficients are derived based on the changes to the index of refraction and length of the fiber. Using the DTS measurements, temperature changes are computed, smoothed, filtered, and compared to the LF-DAS response. Crossplots of the in-well LF-DAS measurements and temperature changes from DTS measurements far from the actively fractured region are used to validate the theoretical sensitivity coefficients. Uncertainty in the temperature component of the LF-DAS response is quantified. The difficulty in corresponding the different spatial and temporal resolutions of the DTS and LF-DAS measurements is overcome by comparing the responses over a moving temporal and spatial window. If the LF-DAS response at the center of the window agrees with the DTS response within uncertainty, the measurement is filtered out. After filtering, the remaining nonzero in-well LF-DAS measurements are due to changes in strain. The data are then visualized in waterfall plots. The results indicate that the theoretical and observed strain and temperature coefficients agree within 10%. After the temperature component of the in-well LF-DAS response is extracted, the remaining nonzero measurements are located primarily within the actively treated region. Locations with peaks in the strain response are interpreted to indicate fracture initiation points. These fracture initiation points are compared with in-well DTS and high-frequency DAS noise measurements across multiple stages to better understand fracture initiation along the horizontal well.

Engineering↗

DOE OSTI2014

Scale-Dependent Fracture-Matrix Interactions and Their Impact on Radionuclide Transport: Development of efficient particle-tracking methods

Matrix Diffusion and Adsorption within a rock matrix are important mechanisms for retarding transport of radionuclides in fractured rock. Due to computational limitations and difficulties in characterizing complex subsurface systems, diffusive exchange between a fracture network and surrounding rock matrix is often modeled using simplified conceptual representations. There is significant uncertainty in “effective” parameters used in these models, such as the “effective matrix diffusivity”. Often, these parameters are estimated by fitting sparse breakthrough data, and estimated values fall outside meaningful ranges, because simplified interpretive models do not consider complex three-dimensional flow. There is limited understanding of the relationship between the effective parameters and rock mass characteristics including network structure and matrix properties. There is also evidence for an apparent scale-dependence in “effective matrix diffusion” coefficients. These observations raise questions on whether fracture-matrix interaction parameters estimated from small-scale tracer tests can be used for predicting radionuclide fate and transport at the scale of DOE field sites. High-resolution three-dimensional Discrete-Fracture-Network-Matrix (DFNM) models based on well-defined local scale transport equations can help to address some of these questions. Due to tremendous advances in computational technology over the last 10 years, DFNM modeling in relatively large domains is now feasible. The overarching objective of our research is to use DFNM modeling to improve fundamental understanding of how effective parameters in conceptual models are related to fracture network structure and matrix properties. An advanced three-dimensional DFNM model is being developed, which combines upscaled particle-tracking algorithms for fracture-matrix interaction and a parallel fracture-network flow simulator. The particle-tracking algorithms allow complexity in flow fields at different scales, and track transport across fracture-matrix interfaces based on rigorous local approximations to the transport equations. This modeling approach can incorporate aperture variability, multi-scale preferential flow and matrix heterogeneity. We developed efficient particle-tracking methods for handling matrix diffusion and adsorption on fracture walls and demonstrated their efficiency for use within the context of large-scale complex fracture network models with variability in apertures across a network of fractures and within individual fractures.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2016

2015 Accomplishments-Tritium aging studies on stainless steel. Effects of hydrogen isotopes, crack orientation, and specimen geometry on fracture toughness

This study reports on the effects of hydrogen isotopes, crack orientation, and specimen geometry on the fracture toughness of stainless steels. Fracture toughness variability was investigated for Type 21-6-9 stainless steel using the 7K0004 forging. Fracture toughness specimens were cut from the forging in two different geometric configurations: arc shape and disc shape. The fracture toughness properties were measured at ambient temperature before and after exposure to hydrogen gas and compared to prior studies. There are three main conclusions that can be drawn from the results. First, the fracture toughness properties of actual reservoir forgings and contemporary heats of steel are much higher than those measured in earlier studies that used heats of steel from the 1980s and 1990s and forward extruded forgings which were designed to simulate reservoir microstructures. This is true for as-forged heats as well as forged heats exposed to hydrogen gas. Secondly, the study confirms the well-known observation that cracks oriented parallel to the forging grain flow will propagate easier than those oriented perpendicular to the grain flow. However, what was not known, but is shown here, is that this effect is more pronounced, particularly after hydrogen exposures, when the forging is given a larger upset. In brick forgings, which have a relatively low amount of upset, the fracture toughness variation with specimen orientation is less than 5%; whereas, in cup forgings, the fracture toughness is about 20% lower than that forging to show how specimen geometry affects fracture toughness values. The American Society for Testing Materials (ASTM) specifies minimum specimen section sizes for valid fracture toughness values. However, sub-size specimens have long been used to study tritium effects because of the physical limitation of diffusing hydrogen isotopes into stainless steel at mild temperatures so as to not disturb the underlying forged microstructure. This study shows that fracture toughness values of larger specimens are higher and more representative of the material’s fracture behavior in a fully constrained tritium reservoir. The toughness properties measured for sub-size specimens were about 65-75% of the values for larger specimens. While the data from sub-size specimens are conservative, they may be overly so. The fracture toughness properties from sub-size specimens are valuable in that they can be used for tritium effects studies and show the same trends and alloy differences as those seen from larger specimen data. Additional work is planned, including finite element modeling, to see if sub-size specimen data could be adjusted in some way to be more closely aligned with the actual material behavior in a fully constrained pressure vessel.

08 HYDROGEN↗

DOE OSTI2016

Characterization of the Oriskany and Berea Sandstones: Evaluating Biogeochemical Reactions of Potential Sandstone–Hydraulic Fracturing Fluid Interaction

The Marcellus shale, located in the mid-Atlantic Appalachian Basin, has been identified as a source for natural gas and targeted for hydraulic fracturing recovery methods. Hydraulic fracturing is a technique used by the oil and gas industry to access petroleum reserves in geologic formations that cannot be accessed with conventional drilling techniques (Capo et al., 2014). This unconventional technique fractures rock formations that have low permeability by pumping pressurized hydraulic fracturing fluids into the subsurface. Although the major components of hydraulic fracturing fluid are water and sand, chemicals, such as recalcitrant biocides and polyacrylamide, are also used (Frac Focus, 2015). There is domestic concern that the chemicals could reach groundwater or surface water during transport, storage, or the fracturing process (Chapman et al., 2012). In the event of a surface spill, understanding the natural attenuation of the chemicals in hydraulic fracturing fluid, as well as the physical and chemical properties of the aquifers surrounding the spill site, will help mitigate potential dangers to drinking water. However, reports on the degradation pathways of these chemicals are limited in existing literature. The Appalachian Basin Marcellus shale and its surrounding sandstones host diverse mineralogical suites. During the hydraulic fracturing process, the hydraulic fracturing fluids come into contact with variable mineral compositions. The reactions between the fracturing fluid chemicals and the minerals are very diverse. This report: 1) describes common minerals (e.g. quartz, clay, pyrite, and carbonates) present in the Marcellus shale, as well as the Oriskany and Berea sandstones, which are located stratigraphically below and above the Marcellus shale; 2) summarizes the existing literature of the degradation pathways for common hydraulic fracturing fluid chemicals [polyacrylamide, ethylene glycol, poly(diallyldimethylammonium chloride), glutaraldehyde, guar gum, and isopropanol]; 3) reviews the known research about the interactions between several hydraulic fracturing chemicals [e.g. polyacrylamide, ethylene glycol, poly(diallyldimethylammonium chloride), and glutaraldehyde] with the minerals (quartz, clay, pyrite, and carbonates) common to the lithologies of the Marcellus shale and its surrounding sandstones; and 4) characterizes the Berea sandstone and analyzes the physical and chemical effects of flowing guar gum through a Berea sandstone core.

58 GEOSCIENCES↗

DOE OSTI2017

Imaging Fracture Networks Using Joint Seismic and Electrical Change Detection Techniques (Final Report)

During the initial phase of this SubTER project, we conducted a series of high resolution seismic imaging campaigns designed to characterize induced fractures. Fractures were emplaced using a novel explosive source, designed at Sandia National Laboratories, that limits damage to the borehole. This work provided evidence that fracture locations could be imaged at inch scales using high-frequency seismic tomography but left many fracture properties (i.e. permeability) unresolved. We present here the results of the second phase of the project, where we developed and demonstrated emerging seismic and electrical geophysical imaging technologies that characterize 1) the 3D extent and distribution of fractures stimulated from the explosive source, 2) 3D fluid transport within the stimulated fracture network through use of a contrasting tracer, and 3) fracture attributes through advanced data analysis. Focus was placed upon advancing these technologies toward near real-time acquisition and processing in order to help provide the feedback mechanism necessary to understand and control fracture stimulation and fluid flow. Results from this study include a comprehensive set of 4D cross-hole seismic and electrical data that take advantage of change detection methodologies allowing for perturbations associated with the fracture emplacement and particulate tracer to be isolated. During the testing the team also demonstrated near real-time 4D electrical resistivity tomography imaging and 4D seismic tomography using the CASSM approach with a temporal resolution approaching 1 minute. All of the data collected were used to develop methods of estimating fracture attributes from seismic data, develop methods of assimilating disparate and transient data sets to improve fracture network imaging resolution, and advance capabilities for near real-time inversion of cross-hole tomographic data. These results are illustrated here. Advancements in these areas are relevant to all situations where fracture emplacement is used for reservoir stimulation (e.g. Enhanced Geothermal Systems (EGS) and tight shale gases).

42 ENGINEERING↗

DOE OSTI2025

Role of Intersections in Fracture Connectivity

Networks of intersecting fractures often provide the flow paths through subsurface reservoirs. Assessing network connectivity is challenging because fracture intersections compose a vanishingly small fraction of the network void volume. In this paper, motivated by 3D X-ray imaging of the simplest element of fracture network, that is, two orthogonal fractures, we perform a percolation and finite-size scaling analysis to study the connectivity provided by fracture intersections. The conditions when an intersection enhances connectivity across a sample depend on spatial correlations in the fracture aperture distributions, on the stress state, and on the direction of flow. Here we consider three flow directions: (a) across intersections, (b) parallel to intersections and (c) around corners. For (a), intersections provide minimal enhancement of connectivity because they contribute little additional void area. For (b), intersections increase the probability of a connected path near threshold by enabling 3D connected pathways that are not possible in parallel fractures. Flow around corners, (c), is fundamentally the result of the intersection connecting two fractures in series and spatial correlations are broken around corners, suppressing the connectivity relative to (a). When the connected fractures are stressed equally, a joint percolation threshold emerges that continues to have scale invariance. However, when the fractures are stressed unequally, the system has mixed percolation without clearly defined percolation thresholds. In all cases, percolation probabilities are found to be scale dependent which has important consequences for the connectivity of larger fracture networks composed of the fundamental element studied here.

02 PETROLEUM↗

DOE OSTI2016

Simulation of Hydraulic and Natural Fracture Interaction Using a Coupled DFN-DEM Model

The presence of natural fractures will usually result in a complex fracture network due to the interactions between hydraulic and natural fracture. The reactivation of natural fractures can generally provide additional flow paths from formation to wellbore which play a crucial role in improving the hydrocarbon recovery in these ultra-low permeability reservoir. Thus, accurate description of the geometry of discrete fractures and bedding is highly desired for accurate flow and production predictions. Compared to conventional continuum models that implicitly represent the discrete feature, Discrete Fracture Network (DFN) models could realistically model the connectivity of discontinuities at both reservoir scale and well scale. In this work, a new hybrid numerical model that couples Discrete Fracture Network (DFN) and Dual-Lattice Discrete Element Method (DL-DEM) is proposed to investigate the interaction between hydraulic fracture and natural fractures. Based on the proposed model, the effects of natural fracture orientation, density and injection properties on hydraulic-natural fractures interaction are investigated.

Discrete Fracture Network↗

DOE OSTI2020

Results of Multiple Tracer Injections into Fractures in the EGS Collab Testbed-1

The EGS Collab project constructed an intermediate scale (~10-20 m) testbed at the 4850 level of the Stanford Underground Research Facility (SURF) in South Dakota for testing and validating fracture stimulation and flow/transport models. This testbed consists of eight ~200 ft (~60 m) HQ-diameter (9.6 cm) boreholes that are drilled into the crystalline rocks of the Poorman Formation from the West Access Drift tunnel. Of the eight boreholes, one borehole is used as an injection/stimulation well, while another sub-parallel borehole located about 10 m away from the injection well is used as a production well, and rest of the other boreholes are used as geophysical/fluid sampling monitoring wells. Hydraulic stimulation activities were conducted at three locations along the injection hole in an attempt to create direct fracture connections to the production hole. A flow system has been established between injection and production boreholes through a set of hydraulically stimulated fractures propagated from a notch located at 164 ft in the injection hole. Although we planned for a single production well, a set of natural fractures in the testbed is believed to have intersected the stimulated hydraulic fractures and provided additional flow paths for water to be transported to the drift through multiple monitoring boreholes and weep zones. As the flow tests continued after stimulation activities, single or a combination of two or three producers become dominant producers at different times, mostly as a response to intersection of hydraulic fracture and testbed wells, activation of natural fractures and making new leak points to the monitoring wells, resealing of leaky wells, and so on. Since late October 2018, multiple tracers were injected into the fracture system at the 164 ft location that involves both stimulated and natural fractures, and tracers were recovered from multiple locations in nearby wells and weep. The cumulative water recoveries over the time have ranged from 50 to 90%; however, the injected tracer recoveries from various production sources are much less (ranging from a few percentages to 38%). Changes in water/tracer recoveries, shifting of major producing wells from one well to the others, and other observations (e.g., microseismic, electric resistivity tomography, etc.) indicate a testbed that has undergone several changes since October 2018. In this paper, we present tracer recovery data accumulated during several tracer campaigns, provide conceptual testbed flow pathways, and provide observations that suggest a evolutionary nature of fracture volume and fracture geometry in the testbed.

15 GEOTHERMAL ENERGY↗

DOE OSTI2021

Sealing of fractures in a representative CO 2 reservoir caprock by migration of fines

The impact of fines migration on fracture transmissivity reduction was investigated by injecting a brine solution containing a suspension of 0.1 wt. % kaolinite particles with a mean particle size distribution of 9.6 µm through fractured shale core samples. The fractures had apertures estimated to be approximately 100 µm. A mass balance approach was used to determine the quantity of kaolinite that was deposited within the fractures (influent – effluent = amount deposited in fractures). Large fractions (44 to 90%) of the suspended kaolinite pumped through the fractures were deposited within the fractures. Based on fracture volumes estimated with X-ray computed tomography, it was determined that approximately 10 to 17% of the fracture volume was filled with kaolinite at the point when flow was completely restricted. These results indicate that 100 µm fractures in CO 2 reservoir caprocks could be sealed within hours if the brines passing through the fractures contain a proportional volume of particulates to the tests performed in this laboratory study.

04 OIL SHALES AND TAR SANDS↗

DOE OSTI2017

Osteophytes and fracture calluses share developmental milestones and are diminished by unloading

Osteophytes are a typical radiographic finding during osteoarthritis (OA), but the mechanisms leading to their formation are not well known. Comparatively, fracture calluses have been studied extensively; therefore, drawing comparisons between osteophytes and fracture calluses may lead to a deeper understanding of osteophyte formation. Here, we compared the time courses of osteophyte and fracture callus formation, and investigated mechanisms contributing to development of these structure. Additionally, we investigated the effect of mechanical unloading on the formation of both fracture calluses and osteophytes. Mice underwent either transverse femoral fracture or non-invasive anterior cruciate ligament rupture. Fracture callus and osteophyte size and ossification were evaluated after 3, 5, 7, 14, 21, or 28 days. Additional mice were subjected to hindlimb unloading after injury for 3, 7, or 14 days. Protease activity and gene expression profiles after injury were evaluated after 3 or 7 days of normal ambulation or hindlimb unloading using in vivo fluorescence reflectance imaging (FRI) and quantitative PCR. We found that fracture callus and osteophyte growth achieved similar developmental milestones, but fracture calluses formed and ossified at earlier time points. Hindlimb unloading ultimately led to a threefold decrease in chondro/osteophyte area, and a twofold decrease in fracture callus area. Unloading was also associated with decreased inflammation and protease activity in injured limbs detected with FRI, particularly following ACL rupture. qPCR analysis revealed disparate cellular responses in fractured femurs and injured joints, suggesting that fracture calluses and osteophytes may form via different inflammatory, anabolic, and catabolic pathways. © 2017 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 36:699–710, 2018.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2019

Modeling Three-Dimensional Fluid-Driven Propagation of Multiple Fractures using TOUGH-FEMM

Here, a novel numerical simulation tool named TOUGH-FEMM is presented and applied to model three-dimensional (3D) hydraulic fracturing in porous rock. The fluid flow in both fractures and porous rock is modeled using TOUGH2, which is a well-established code for analysis of multiphase and multi-component fluid flow. Rock deformations associated with fracture propagation are modeled using finite element-meshfree method (FEMM). FEMM is an approach to simulate fracture propagation without remeshing, in which the fracture path does not need to be predetermined. Fracture mechanics with mixed-mode stress intensity factors are employed to detect fracture instability and determine the direction of fracture propagation. TOUGH-FEMM is confirmed for modeling fluid-driven fracture propagation in 3D through a number of simulation examples, including modeling of hydraulic fracturing laboratory experiments and by comparison to independent numerical simulation results for multiple interacting hydraulic fractures at ten to hundred meter scale.

58 GEOSCIENCES↗

DOE OSTI2021

Experimental Investigation of Elastodynamic Nonlinear Response of Dry Intact, Fractured and Saturated Rock

Nonlinear elastodynamic response of fractured rocks carries crucial information on fracture features that can be exploited to forecast flow properties, friction constitutive behavior and poromechanical response. Well-controlled laboratory experiments are designed to measure the nonlinear elastodynamic response of Westerly granite in three states: dry intact, dry fractured and saturated fractured. We study the effect of fracturing and saturation in modifying the elastodynamic response of the rock. Each sample is tested at a normal stress level of 15 MPa. We measure the elastodynamic response of an intact L-shaped sample of Westerly Granite subjected to normal stress oscillations of prescribed amplitudes (0.2–1 MPa) and frequencies (0.1, 1, 10 Hz). Ultrasonic waves transmitted across the sample are used to monitor the evolution of wave velocity before, during and after dynamic stressing. Furthermore, the nonlinearity of the elastodynamic response is measured in terms of: (1) the offset in normalized wave velocity; (2) the amplitude of wave velocity fluctuation during the oscillations; and (3) recovery rate of the wave velocity post-oscillation. We observe that the three nonlinearity parameters show a similar trend. Irrespective of the parameter, the nonlinearity measures higher for sample in dry-intact condition than that for dry-fractured and the saturated-fractured sample exhibits smaller nonlinearity than the dry-fractured sample. As expected, the saturated sample exhibits less nonlinearity than the dry intact and fractured samples due to the presence of interstitial fluid and the resulting increased interface stiffness. Conversely, the dry intact rock shows a higher nonlinearity than the dry fractured. We use numerical simulations to show that the presence of fracture significantly alters the strain distribution across the bulk of the sample and only the contacting asperities are highly strained, thus resulting in a decrease in the measured elastodynamic nonlinearity.

58 GEOSCIENCES↗

DOE OSTI2019

Investigation of Stress Field and Fracture Development During Shale Maturation Using Analog Rock Systems

The emergence of hydrocarbons within shale as a major recoverable resource has sparked interest in fluid transport through these tight mudstones. Recent studies suggest the importance to recovery of microfracture networks that connect localized zones with large organic content to the inorganic matrix. This paper presents a joint modeling and experimental study to examine the onset, formation, and evolution of microfracture networks as shale matures. Both the stress field and fractures are simulated and imaged. A novel laboratory-scale, phase-field fracture propagation model was developed to characterize the material failure mechanisms that play a significant role during the shale maturation process. The numerical model developed consists of coupled solid deformation, pore pressure, and fracture propagation mechanisms. Benchmark tests were conducted to validate model accuracy. Laboratory-grade gelatins with varying Young’s modulus were used as scaled-rock analogs in two dimensional Hele-Shaw cell apparatus. Yeast within the gelatin generates gas in a fashion analogous to hydrocarbon formation as shale matures. These setups allow study and visualization of host rock elastic-brittle fracture and fracture network propagation mechanisms. The experimental setup was fitted to utilize photoelasticity principles coupled with birefringence properties of gelatin to explore visually the stress field of the gelatin as the fracture network developed. Stress optics image analysis and Linear Elastic Fracture Mechanics (LEFM) principles for crack propagation were used to monitor fracture growth for each gelatin type. Observed and simulated responses suggest gas diffusion within and deformation of the gelatin matrix as predominant mechanisms for energy dissipation depending on gelatin strength. LEFM, an experimental estimation of principal stress development with fracture growth, at different stages was determined for each gelatin rheology. Here, the interplay of gas diffusion and material deformation determines the resulting frequency and pattern of fractures. Results correlate with Young’s modulus. Experimental and computed stress fields reveal that fractures resulting from internal gas generation are similar to, but not identical to, type 1 opening mode.

58 GEOSCIENCES↗

DOE OSTI2024

Numerical investigation of multiphase flow through self-affine rough fractures

Multiphase flow through fractures has great significance in subsurface energy recovery and gas storage applications. Different fracture and flow properties affect flow through a fracture which is difficult to control in laboratory experiments. Here, we perform lattice Boltzmann simulations in an ensemble of synthetically generated fractures. Drainage simulations are performed at different capillary numbers, wettability, and viscosity ratios. We track the invading front and quantify breakthrough saturations and show that roughness and wettability have a strong effect on fluid invasion through a complex fracture. Invading a more viscous fluid results in more stable displacement regardless of the capillary number while at very low capillary numbers, fluid migration is dependent on the inherent structure of the fracture. We develop a fluid displacement phase diagram in a single rough fracture and compare our results from that in the literature. Finally, we extend the phase diagrams across multiple fractures and demonstrate the importance of natural fracture features of roughness and wettability in identifying stable versus unstable displacement regimes during multiphase flow through rough fractures. Our work presents an end-to-end numerical pathway for testing on experimental data and expanding numerical data sets for testing combinations of different physical phenomenon and make valuable predictions on fluid flow through rough fractures.

02 PETROLEUM↗

DOE OSTI2018

Effectiveness of supercritical-CO 2 and N 2 huff-and-puff methods of enhanced oil recovery in shale fracture networks using microfluidic experiments

Current oil recovery methods in hydraulically fractured shale reservoirs have a low recovery efficiency of about 10%. The objective of this work is to investigate the effectiveness of nitrogen and supercritical carbon dioxide in the huff-and-puff method for enhanced oil recovery as a means to re-energize reservoirs and improve recovery rates. Here, we conduct direct visualization experiments with a microfluidic system to reveal the mechanisms and to quantify the recovery rates of oil from fracture networks. We compared the effectiveness of water, nitrogen and supercritical carbon dioxide at reservoir conditions in a process mimicking the huff-and-puff method in both dead-end and connected fracture systems. The microfluidic chips were made of glass and placed in a confining pressure system pressurized to 10 MPa, 50 °C. The system was allowed to equilibrate, and then depressurized to simulate huff-and-puff oil recovery. Fluorescence microscope images were continuously taken to visualize and calculate residual oil saturation as a function of pressure drawdown. As the system was depressurized from 10 MPa, gas exsolution from the oil liquid phase, including bubble nucleation, growth, and coalescence, appeared to be the main driver for mobilizing oil from the fracture networks. Injection of supercritical CO 2 resulted in the highest recovery rate with an average end-point recovery of about 90% in the connected fracture network and 60% in the dead-end fracture network. N 2 has lower solubility in oil and hence showed a lower recovery rate of 40% in the connected fracture network and 25% in the dead-end fracture network. Injection of water had no effect on oil mobilization since water is insoluble, immiscible and incompressible. The main mechanism of enhanced recovery was gas exsolution from the liquid phase as pressure was decreased below the bubble point pressure. Because the gas was distributed throughout the oil phase, bubble nucleation, growth, coalescence, and elongation occurred throughout the fracture network. Expansion of the gas forced oil out of the network through piston displacement in continuous oil areas and film flow in the dispersed oil areas. Bubbles began as spheres and grew until they touched the fracture walls where they elongated along the fracture length. The bubble growth rate depended on local mass transfer from liquid to gas phase and gas volume expansion due to pressure drop. The efficiency of the huff-and-puff process is dependent on the solubility and miscibility of the injection fluid with oil. High gas solubility allows for more bubble nucleation, growth and expansion during the depressurization cycle.

58 GEOSCIENCES↗

DOE OSTI2014

Geochemical and Geomechanical Effects on Wellbore Cement Fractures

Experimental studies were conducted using batch reactors, X-ray microtomograpy (XMT), and computational fluid dynamics (CFD) simulation to determine changes in cement fracture surfaces, fluid flow pathways, and permeability with geochemical and geomechanical processes. Composite Portland cement-basalt caprock core with artificial fractures was prepared and reacted with CO2-saturated groundwater at 50°C and 10 MPa for 3 to 3.5 months under static conditions to understand the geochemical and geomechanical effects on the integrity of wellbores containing defects. Cement-basalt interface samples were subjected to mechanical stress at 2.7 MPa before the CO2 reaction. XMT provided three-dimensional (3-D) visualization of the opening and interconnection of cement fractures due to mechanical stress. After the CO2 reaction, XMT images revealed that calcium carbonate precipitation occurred extensively within the fractures in the cement matrix, but only partially along fractures located at the cement-basalt interface. The permeability calculated based on CFD simulation was in agreement with the experimentally measured permeability. The experimental results imply that the wellbore cement with fractures is likely to be healed during exposure to CO2-saturated groundwater under static conditions, whereas fractures along the cement-caprock interface are still likely to remain vulnerable to the leakage of CO2. CFD simulation for the flow of different fluids (CO2-saturated brine and supercritical CO2) using a pressure difference of 20 kPa and 200 kPa along ~2 cm-long cement fractures showed that a pressure gradient increase resulted in an increase of CO2 fluids flux by a factor of only ~3-9 because the friction of CO2 fluids on cement fracture surfaces increased with higher flow rate as well. At the same pressure gradient, the simulated flow rate was higher for supercritical CO2 than CO2-saturated brine by a factor of only ~2-3, because the viscosity of supercritical CO2 is much lower than that of CO2-saturated brine. The study suggests that in deep geological reservoirs the geochemical and geomechanical processes have coupled effects on the wellbore cement fracture evolution and fluid flow along the fracture surfaces.

36 MATERIALS SCIENCE↗