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 19 records

Geomechanical and Hydrogeological Evaluation of a Shallow Hydraulic Fracture at the Devine Fracture Pilot Site, Medina County, Texas

UT-Austin’s Devine Fracture Pilot Site, 50 miles southwest of San Antonio, Texas, has been targeted for a comprehensive, multidisciplinary development of fracture diagnostic techniques that are cross-validated by ground-truth data acquisition near a recently created, 175-ft-deep, horizontal hydraulic fracture (Ahmadian et al. 2018 Demonstration of proof of concept of electromagnetic geophysical methods for high resolution illumination of induced fracture networks. In Proceedings of the SPE Hydraulic Fracturing Technology Conference and Exhibition, The Woodlands, Texas, USA, 23–25 January 2018. SPE-189858-MS.). To evaluate the fracture diagnostic methods at this site, we conducted injection tests with a predefined volumetric flow-rate profile, resembling a diagnostic fracture injection test on September 2020. Subsequently, we developed hydrogeological and geomechanical models based on flow-rate and bottomhole-pressure measurements. History-matching efforts using a simplified layer-cake hydrogeological model resulted in the field-scale formation permeability of 9.87 × 10 –15 m 2 (10 mD) and Darcy-scale fracture permeability. The analysis of the bottomhole pressure and injection-rate history showed that (1) the newly created horizontal fracture was closed adjacent to the injection well pre-injection and (2) the initial pump-pressure increase at a nominal volumetric injection rate led to near-well fracture reopening, fluid conductivity increase, and abrupt injection-rate increase. To overcome hydrogeological-model limitations of predicting fracture reopening throughout injection, we extended the modeling to a finite-element, poroelastic analysis of horizontal-fracture growth using a cohesive-zone model. Using this fracture-reopening model, we improved the history match of the transient-pressure response during the experiment by adjusting the hydromechanical properties. Furthermore, post-injection pressure transient analyses helped reduce uncertainty in the overburden-stress gradient, and the initial hydraulic-fracturing simulation verified the plausibility of achieving the surveyed propped fracture area.

02 PETROLEUM↗

Diagnostic assessment of reservoir response to fracturing: a case study from Hydraulic Fracturing Test Site (HFTS) in Midland Basin

Abstract This paper outlines a data collection and diagnostics case study involving multiple horizontal shale wells. We look at well production profiles using rate transient analysis, differences in near wellbore complexity, geologic variations within the area of interest, as well as compositional differences in the rocks based on cores obtained from within the stimulated reservoir. The Hydraulic Fracturing Test Site is a multi-well experiment involving stimulation of unconventional shale wells in the southeastern Midland portion of the Permian Basin. The targeted formations include both the upper as well as the middle Wolfcamp formations, also referred alternatively as Wolfcamp A and Wolfcamp B. Data integration and analysis shared in this paper help us understand the various geologic controls impacting well productivity, particularly the wide variance observed between the Wolfcamp A and Wolfcamp B formations. Rate transient analysis indicates similar system permeabilities for stimulated wells. However, we observe higher effective fracture half-lengths for upper Wolfcamp wells. Using observations from 3D seismic interpretations (such as pad scale faults) as well as petrophysical and image log data, we highlight the substantial differences in stimulation as we move along the well laterals from the heel toward the toe sections. These differences are further reconciled with observations from zones with high data density at the core locations through stimulated rock, as well as independent data such as microseismic emissions. At the test site, Wolfcamp A was found to be relatively quartz rich with significant heterogeneity whereas Wolfcamp B is richer in clay and organic content. This impacts the geomechanical characteristics of the rock mass with much higher natural fracture density in the shallower interval. Thus, the fracture growth is more uniform in the deeper interval and more heterogeneous with branching likely in upper interval. Increased complexity also leads to consistently better productivity from the wells in the shallower interval as demonstrated from RTA results. This case study is unique because it provides valuable insights from actual sampling of the stimulated zones in hydraulically fractured wells and helps understand impact of various factors that contribute toward variability in well production. The findings from this study provides insights into need for optimization of completion designs in the various Wolfcamp landing zones, such as optimization of cluster or fracture spacing in various Wolfcamp intervals. In addition, it provides a useful template for data collection and research direction in future field test sites of similar nature in unconventional reservoirs.

Energy & Fuels↗

Sequentially coupled flow and geomechanical simulation with a discrete fracture model for analyzing fracturing fluid recovery and distribution in fractured ultra-low permeability gas reservoirs

More accurate characterization and prediction of the in-situ distribution of fracturing fluid in fractured reservoirs are needed for enhancing well productivity. In this study, an implicit-sequentially coupled flow/geomechanics simulator incorporating an efficient discrete fracture model is developed to model fluid distribution and recovery performance of ultra-low permeability gas reservoirs. The finite-volume and finite-element methods are used for space discretization of the flow and geomechanics equations, respectively, while the backward Euler method is employed for time discretization. The flow and geomechanics equations are solved sequentially based on fixed-stress splitting. An efficient discrete-fracture model is used to explicitly model the fractured system. Flexible unstructured gridding is employed to model arbitrarily-oriented fractures. The interrelations among pore volume, permeability and geomechanical conditions are considered dynamically using two-way coupled flow and geomechanics computations. The geometry of fracture (networks) due to hydraulic fracturing has significant impacts on the fracturing fluid recovery efficiency and ensuing fluid distribution. Under the same injection volume, the fracturing fluid recovery is higher when the fracture geometry is planar. Fluid recovery is relatively lower whenever natural fractures are activated during fracturing treatments; flowback time is also shortened when complex fracture network with enlarged fracture interface is present. Fracturing fluid in hydraulic fractures may leak off into the natural fractures and subsequently imbibes into the surrounding matrix due to capillarity effects. The fracturing fluid recovery and in-situ fluid distribution are sensitive to the shut-in duration and fracture closure behavior. This study analyzes the coupled flow-geomechanical responses of fractured gas reservoirs during the post-fracturing periods. Understanding the fate of the fracturing fluid can provide insights on, to some extent, the stimulated fracture volume, size of the water invasion zone, and efficiency of the fracturing design. The simulation predictions can also provide more accurate initial reservoir conditions (e.g. distributions of different phases and pressure) for long-term well performance estimation.

42 ENGINEERING↗

Close Observation of Hydraulic Fracturing at EGS Collab Experiment 1: Fracture Trajectory, Microseismic Interpretations, and the Role of Natural Fractures

Abstract Despite the wide application of hydraulic fracturing in enhanced geothermal system (EGS) development and unconventional hydrocarbon production, our understanding of the many factors affecting the propagation of hydraulic fractures has relied on circumstantial evidence. This is partly due to the scarcity of direct observations in the subsurface. The EGS Collab project attempts to address these issues in the context of EGS research by performing intermediate‐scale (∼10 m) hydraulic stimulation experiments in a thoroughly characterized and heavily instrumented underground testbed. This paper analyzes the data collected from the first suite of hydraulic fracturing tests in this testbed, consisting of seven stimulation episodes. High‐quality microseismic data delineate five planar features very clearly. Combining fracture‐wellbore intersection observations from distributed temperature sensing (DTS) and visual observations from an open‐hole well, as well as prior in‐situ stress measurements, we conclude with high certainty that the four larger planes were hydraulic fractures. The growth of the hydraulic fractures was temporarily halted by a prominent, open natural fracture in the testbed but they eventually crossed and slightly reoriented under continued stimulation. Mineral‐filled (i.e., healed) natural fractures, though prevalent in the testbed, did not have perceivable effects on hydraulic fracture propagation. The high‐quality, mutually corroborating data sets allowed conclusions to be drawn with high confidence and attests to the advantage of intermediate‐scale experiments in subsurface research.

58 GEOSCIENCES↗

Fracture and Stress Evolution on Europa: New Insights Into Fracture Interpretation and Ice Thickness Estimates Using Fracture Mechanics Analyses

The work completed during the funding period has provided many important insights into fracturing behavior in Europa's ice shell. It has been determined that fracturing through time is likely to have been controlled by the effects of nonsynchronous rotation stresses and that as much as 720 deg of said rotation may have occurred during the visible geologic history. It has been determined that there are at least two distinct styles of strike-slip faulting and that their mutual evolutionary styles are likely to have been different, with one involving a significant dilational component during shear motion. It has been determined that secondary fracturing in perturbed stress fields adjacent to older structures such as faults is a prevalent process on Europa. It has been determined that cycloidal ridges are likely to experience shear stresses along the existing segment portions as they propagate, which affects propagation direction and ultimately induces tailcracking at the segment tip than then initiates a new cycle of cycloid segment growth. Finally, it has been established that mechanical methods (e.g., flexure analysis) can be used to determine the elastic thickness of the ice shell, which, although probably only several km thick, is likely to be spatially variable, being thinner under bands but thicker under ridged plains terrain.

Kattenhorn, Simon↗

The Role of Chemistry in Fracture Pattern Development and Opportunities to Advance Interpretations of Geological Materials

Abstract Fracture pattern development has been a challenging area of research in the Earth sciences for more than 100 years. Much has been learned about the spatial and temporal complexity inherent to these systems, but severe challenges remain. Future advances will require new approaches. Chemical processes play a larger role in opening‐mode fracture pattern development than has hitherto been appreciated. This review examines relationships between mechanical and geochemical processes that influence the fracture patterns recorded in natural settings. For fractures formed in diagenetic settings (~50 to 200 °C), we review evidence of chemical reactions in fractures and show how a chemical perspective helps solve problems in fracture analysis. We also outline impediments to subsurface pattern measurement and interpretation, assess implications of discoveries in fracture history reconstruction for process‐based models, review models of fracture cementation and chemically assisted fracture growth, and discuss promising paths for future work. To accurately predict the mechanical and fluid flow properties of fracture systems, a processes‐based approach is needed. Progress is possible using observational, experimental, and modeling approaches that view fracture patterns and properties as the result of coupled mechanical and chemical processes. A critical area is reconstructing patterns through time. Such data sets are essential for developing and testing predictive models. Other topics that need work include models of crystal growth and dissolution rates under geological conditions, cement mechanical effects, and subcritical crack propagation. Advances in machine learning and 3‐D imaging present opportunities for a mechanistic understanding of fracture formation and development, enabling prediction of spatial and temporal complexity over geologic timescales. Geophysical research with a chemical perspective is needed to correctly identify and interpret fractures from geophysical measurements during site characterization and monitoring of subsurface engineering activities.

58 GEOSCIENCES↗

Laboratory investigation of hydraulic fracturing in granitic rocks using active and passive seismic monitoring

SUMMARY Knowledge of the fracturing processes can be important for the optimization of pressurized fluid injection operations in the deep underground rock mass. Active and passive seismic monitoring techniques have been used in the field for tracking or mapping the propagating hydraulic fracture. Although both these monitoring techniques provide valuable information about the generated fracture network, it is difficult for either technique to comprehensibly identify the different processes associated with hydraulic fracturing. The combined active and passive monitoring has the potential for better characterization of the complex hydraulic fracturing phenomena. In this study, laboratory hydraulic fracturing experiments with combined active and passive seismic monitoring were conducted on true triaxially loaded Barre granite cubes with different fluid injection rates. The seismic inelastic fracturing was detected by 16 passive acoustic emission sensors, where 3678 and 2370 seismic source events were detected for the high and low injection rate experiments, respectively. For active monitoring, strong variations in the attributes of signals were observed which were transmitted through four source–receiver pairs, placed both perpendicular and parallel to the generated hydraulic fracture. Positive velocity changes were observed for active sensor pairs with ray paths passing through the generated hydraulic fracture indicating fluid permeation, whereas isolated dry deformation was characterized by a slight but permanent velocity decrease. Compared to velocity, the energy of the active signals was 1–2 orders of magnitude more sensitive to different hydraulic fracturing processes. However, the sensitivity and signatures of the active signal attributes were found to be dependent on the frequency range and direction of ray path with respect to the location of the generated fracture network. Using the coupled evaluation of the active and passive signals we were able to systematically identify various hydraulic fracturing processes including: (1) aseismic deformation, (2) fracture initiation and fluid permeation, (3) pressure build-up, (4) fracture propagation and (5) pressure release and leak-off. The results of this study showed that combining the respective advantages of active and passive seismic techniques and using both of them to monitor the failure processes can facilitate a more comprehensive understanding and better control of the hydraulic stimulations in subsurface operations.

Geochemistry & Geophysics↗

Developing Upscaling Approach for Swarming Hydraulic Fractures Observed at Hydraulic Fracturing Test Site through Multiscale Simulations

We show that this work aims to address a challenge posed by recent observations of tightly spaced hydraulic fractures in core samples from the hydraulic fracturing test site (HFTS) in the Middle Wolfcamp Formation. Many fractures in retrieved cores have subfoot spacing, which is at odds with conventional models in which usually one hydraulic fracture is initiated per cluster. Models assuming a single fracture at each cluster, although a common practice, often predict excessive fracture propagation that is inconsistent with microseismic observation. Here, we aim to develop a numerical approach to effectively account for densely spaced hydraulic fractures in field-scale simulations. Because it is impractical to explicitly model all aforementioned fractures, we develop a new upscaling law that enables existing simulation tools to predict reservoir response to fracture swarms. The upscaling law is derived based on an energy equivalence argument and validated through multiscale simulations using a high-fidelity code, GEOS. The swarming fractures are first modeled with a spacing that is much smaller than the cluster spacing; these fractures are then approximated by an upscaled, single fracture based on the proposed upscaling law. The upscaled fracture is shown to successfully match the energy input rate and produce the total fracture aperture and average propagation length of the explicitly simulated swarm. Afterward, the upscaling approach is further implemented in 3D field-scale simulations and validated against the HFTS microseismic data of a horizontal well. Our results show that hydraulic fracture swarming can significantly affect fracture propagation behaviors compared with the propagation of single fractures as assumed by conventional modeling approaches. Under the considered situations, the conventional treatment yields fast propagation speed that far exceeds that indicated by the microseismic data. We also illustrate that this discrepancy can be reduced readily through the implementation of the upscaling law. Our results demonstrate the importance of accounting for the fracture swarming effect in field-scale simulations and the efficacy of this approach to enable realistic predictions of reservoir responses to fracture swarms, without the need to model tightly spaced fractures individually.

02 PETROLEUM↗

Experimental study of three-dimensional CO 2 -water drainage and fracture-matrix interactions in fractured porous media

Estimation of CO 2 storage capacity in fractured porous reservoirs requires a better understanding of the CO 2 -water flow fundamentals in fracture-matrix systems. There are few core-flood experiments on three-dimensional CO 2 -water drainage in a fracture-matrix system, yet none have examined the impacts of the capillary continuity of fractures and fracture-matrix interactions. In this study, twelve drainage experiments were conducted in four fracture-matrix columns, each comprising a vertical stack of a cylindrical rock core, a filter paper serving as an analogue to a horizontal fracture, and a ceramic plate. Here, the core sample was surrounded by open space at the top and circumferential sides to model fractures, allowing for 3-D CO 2 -water drainage in the rock core and displaced water draining across the horizontal fracture. X-ray computed tomography was conducted to visualize the dynamic invasion/drainage processes in four rock core samples showing contrasts in anisotropy, permeability, and heterogeneity. Experimental results show (1) the equilibrium CO 2 saturations in the rock matrix vary from 0.10 to 0.60 at controlled capillary pressures up to 200 kPa, (2) the CO 2 saturations in the matrix increase with water saturation in the horizontal fracture resulting in better capillary continuity for water to drain across the fracture, and (3) the fracture water saturation can be enhanced by the non-uniform fracture capillary pressure and countercurrent flow of CO 2 and water across the fracture-matrix interface. In the case of high fracture water saturation, matrix CO 2 saturation largely depends on matrix anisotropy and heterogeneity. The core-scale experimental results contribute to understand the fracture-matrix interactions and CO 2 storage efficiency in fractured porous media.

54 ENVIRONMENTAL SCIENCES↗

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

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

15 GEOTHERMAL ENERGY↗

Pressure Transient Analyses and Poroelastic Modeling of Hydraulic Fracture Dilation for Multiple Injections at the Devine Fracture Pilot Site

Our team has conducted electromagnetic (EM) surveys for the past six years to monitor hydraulic-fracture behavior at the Devine Fracture Pilot Site (DFPS). The sub-horizontal orientation of a shallow hydraulic fracture at the DFPS provides uniform access to the fracture area for interrogation and data collection. Ahmadian et al. (2023) suggested a possible correlation between spatiotemporal changes in the flow rate, bottomhole pressure (BHP), and the observed surface recorded electric field at the DFPS. In this paper, we present the development of poroelastic forward models and pressure transient analyses (PTAs) to support the development of a multiphysics inverse model for these EM surveys. First, we conducted PTAs of the shut-in periods after six injections out of 10 to determine the fracture closure pressure (FCP) or the overburden pressure used in a poroelastic fracture reopening model. Second, we developed a finite-element poroelastic model throughout five injection cycles to include the effect of the cumulative injected volumes due to the previous injections on current fracture dilation in the presence of highly permeable unpropped and propped zones adjacent to the cohesive layer that models fracture reopening. Fracture reopening in this poroelastic model is based on a calibrated traction-separation response using the bottomhole pressure collected in two injection campaigns in 2020 and 2022. We used the outcomes of a previous simulation study of the primary hydraulic-fracturing stimulation to define the dimension of an unpropped fracture zone ahead of the propped fracture area. The PTAs led to FCPs consistent with those obtained using the injection data collected at the DFPS in 2020. Further, these analyses showed that at later injections, the fracture closure occurred at a later time with respect to the shut-in time, inferring the effect of cumulative injected volumes in previous injections. The simulation results show that considering the propped and unpropped fracture zones improves our poroelastic model in predicting the injection-well BHP. The numerical simulation results demonstrate a significant excess pore pressure near the fracture because of the preceding formation loadings by the previous injections. The obtained fracture dilation area and fluid pressure distribution provide a basis to improve the development of a multiphysics inverse model. Furthermore, in an iteratively coupled scheme, this pressure distribution can be introduced into EM models to render a holistic view of the causative mechanisms for the surface signal anomalies.

02 PETROLEUM↗

Hydraulic-Fracture-Height Growth Under the Combined Influence of Stress Barriers and Natural Fractures

A fully coupled finite-element/finite-volume code is used here to model 3D hydraulically driven fractures under the influence of strong vertical variations in closure stress interacting with natural fractures. Previously unknown 3D interaction mechanisms on fracture-height growth are revealed. Slipping of a natural fracture, triggered by elevated fluid pressure from an intersecting hydraulic fracture, can induce both increases and decreases of normal stress in the minimum-horizontal-stress direction, toward the center and tip of the natural fracture, respectively. Consequently, natural fractures are expected to be able to both encourage and inhibit the progress of hydraulic fractures propagating through stress barriers, depending on the relative locations between the intersecting fractures. Once the hydraulic fracture propagates above the stress barrier through the weakened segment near a favorably located natural fracture, a configuration consisting of two opposing fractures cuts the stress barrier from above and below. The fluid pressure required to break the stress barrier under such opposing-fracture configurations is substantially lower than that required by a fracture penetrating the same barrier from one side. Sensitivity studies of geologic conditions and operational parameters have also been performed to explore the feasibility of controlled fracture height. The interactions between hydraulic fractures, natural fractures, and geologic factors such as stress barriers in three dimensions are shown to be much more complex than in two dimensions. Although it is impossible to exhaust all the possible configurations, the ability of a 3D, fully coupled numerical model to naturally capture these processes is well-demonstrated.

42 ENGINEERING↗

Fracture Toughness Evaluation for Sandia Mock-up Stainless Steel Canister Weldment Using Spiral Notch Torsion Fracture Toughness Test

Many radioactive materials within the nuclear fuel cycle present a significant hazard. Such materials include spent nuclear fuel, high-level waste, legacy waste and other nuclear materials that have no current outlet than storage. These materials are often held in long-term storage as an interim stage within their lifecycle. Lifecycles can include reuse or disposal. Safety and security of spent nuclear fuel (SNF) interim storage installations are very important, due to a great concentration of fission products, actinides and activation products. Stress corrosion cracking of interim storage containers has been identified as a high priority data gap by the Department of Energy. However, little has been done with regards to the canister weld material properties, including the fracture toughness, and their impact on stress corrosion, until recently from Sandia National Laboratory canister mockup program. Furthermore, because no post-weld heat treatment was required, the associated high tensile residual stress within these canister welds can drive the initiation and growth of stress corrosion cracking (SCC) cracks. Thus, the fracture mechanical property, such as fracture toughness, is an essential information to understand the canister weldment fracture resistance behavior and to support mitigation protocol development to counter SCC. In this work, we carried out fracture toughness evaluation on the received mock-up canister weldments from Sandia, using spiral notch torsion fracture test technology. The details of the weldment fracture toughness evaluation program development at ORNL are presented in this progress report. The objective of this research is to collect fracture toughness data on the mock-up canister weldment from Sandia National Laboratory. Data will be collected under quasi-static fracture loading using the Oak Ridge National Laboratory (ORNL) developed Spiral Notch Torsion Test (SNTT) technology carried out on a biaxial tension/torsion tester. These data will be used to support the ongoing SNF dry storage canister modeling activities and to address regulatory issues associated with canister container integrity and the follow-on SNF system transport reliability. SNTT has been a recent breakthrough in measuring fracture toughness for different materials, including metals, ceramics, concrete and polymers composites. Due to its high geometry constraint and unique loading condition, SNTT can be used to measure the fracture toughness with smaller specimens without concern of size effects. The application of SNTT to brittle materials has been proved to be successful. The micro-cracks induced by original notches in brittle materials could ensure crack growth in SNTT samples. Therefore, no fatigue pre-cracks are needed. The application of SNTT to the ductile material to generate valid toughness data will require a test sample with sufficient crack length to increase the sample fracture tip localized constraint. Fatigue pre-crack growth techniques are employed to introduce sharp crack front into the sample. Previously, only rough calculations were applied to estimate the compliance evolution in the SNTT crack growth process; while accurate quantitative descriptions have never been attempted. This generates an urgent need to understand the crack evolution during the SNTT fracture testing process of ductile materials. The newly developed governing equations for SNTT crack growth estimate are used to control the fatigue crack growth to effectively reaching the target crack length. The detailed SNTT approach and its estimated fracture toughness for the as-received Sandia canister weldment and the associated baseline 304 stainless steel are presented in this report. In the past SNTT technology has also been successfully applied to investigate the fracture behavior of X52 and X80 steels and the X52 welded materials used for hydrogen infrastructures. The current study in supporting Sandia SNF Dry Storage Mock-up Canister Program is the first attempt of applying SNTT to the highly ductile stainless steel and the associated weld materials. The SNTT test results indicate that SNTT method is a reliable test approach with good repeatability in applying to SS304 steel and SS304/308 weld material. The estimate JQ upon fracture for the baseline SS304 steel is at 1618 lb/in. The estimate JQ upon fracture for the SS304/308 weld from HAZ and weld specimens are 720 lb/in and 850 lb/in, respectively.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Effects of Dead‐End Fractures on Non‐Fickian Transport in Three‐Dimensional Discrete Fracture Networks

Abstract Understanding mechanistic causes of non‐Fickian transport in fractured media is important for many hydrogeologic processes and subsurface applications. This study elucidates the effects of dead‐end fractures on non‐Fickian transport in three‐dimensional (3D) fracture networks. Although dead‐end fractures have been identified as low‐velocity regions that could delay solute transport, the direct relation between dead‐end fractures and non‐Fickian transport has been elusive. We systematically generate a large number of 3D discrete fracture networks with different fracture length distributions and fracture densities. We then identify dead‐end fractures using a novel graph‐based method. The effect of dead‐end fractures on solute residence time maximizes at the critical fracture density of the percolation threshold, leading to strong late‐time tailing. As fracture density increases beyond the percolation threshold, the network connectivity increases, and dead‐end fractures diminish. Consequently, the increase in network connectivity leads to a reduction in the degree of late‐time tailing. We also show that dead‐end fractures can inform about main transport paths, such as the mean tortuosity of particle trajectories. This study advances our mechanistic understanding of solute transport in 3D fracture networks.

58 GEOSCIENCES↗

Impact of Concurrent Solubilization and Fines Migration on Fracture Aperture Growth in Shales during Acidized Brine Injection

Acidic hydraulic fracturing fluid chemically and physically alters shale rock fabric during injection and shut-in, creating a “reaction-altered zone” along the fracture faces. To better characterize the variable thickness and composition of this reaction-altered zone under advective flow, we take a coupled experimental and modeling approach. A fluidic cell, with six fiducial markers, is first fabricated to keep the rock sample in place during the core floods and to allow image alignment of acquired images. Then, we conduct a series of reactive core floods in a clay-rich siliceous Wolfcamp shale sample with 10 wt % carbonate, using a synthetic fracturing fluid under no confining stress and at room temperature. High-resolution computed tomography (CT) scans are periodically conducted to observe the spatial alteration of the fracture network. We then perform scanning electron microscopy (SEM-EDS) on the two orthogonal surfaces (fracture surface and freshly cut profile face) to generate high-resolution elemental maps that show the change in mineralogy, both with distance along a given flow path along a fracture surface and with depth from the fracture surface into the shale matrix. These results are contrasted against a two-dimensional (2D) advection-diffusion reaction model developed previously for batch reactions between shale and synthetic fracturing fluids. The model simulates the geochemical interaction occurring at the fracture/matrix interface and penetrating into the shale matrix during the reactive core flood. Both model and experimental results show that the acidic brine is neutralized during the core flood, corresponding to an increase in fracture aperture as a function of fluid volume injected with the greatest change near the inlet. SEM-EDS scans reveal significant dissolution of carbonates on the fracture surface without pyrite oxidation. The reactive transport model indicates that carbonate depletion into the shale interior should be observable, yet SEM-EDS shows no discernible loss of carbonate in the orthogonal profile face. The combination of these observations suggests an additional fracture evolution mechanism in the reactive system, i.e., fines migration. We show that fines migration enhances the access of fracturing fluid to the matrix resulting in a more pronounced fracture widening. Finally, we conclude that coupled mineral dissolution and fines migration govern fracture aperture growth during acidized brine injection. In this work, we effectively show the underlying risk of relying solely on models that do not include an important (transport) process that can alter the system significantly and propose a combined chemomechanical mechanism for fracture evolution appropriate for this shale mineralogy.

58 GEOSCIENCES↗