skip to main content

This content will become publicly available on August 25, 2024

Title: A phase‐field model for hydraulic fracture nucleation and propagation in porous media

Many geo‐engineering applications, for example, enhanced geothermal systems, rely on hydraulic fracturing to enhance the permeability of natural formations and allow for sufficient fluid circulation. Over the past few decades, the phase‐field method has grown in popularity as a valid approach to modeling hydraulic fracturing because of the ease of handling complex fracture propagation geometries. However, existing phase‐field methods cannot appropriately capture nucleation of hydraulic fractures because their formulations are solely energy‐based and do not explicitly take into account the strength of the material. Thus, in this work, we propose a novel phase‐field formulation for hydraulic fracturing with the main goal of modeling fracture nucleation in porous media, for example, rocks. Built on the variational formulation of previous phase‐field methods, the proposed model incorporates the material strength envelope for hydraulic fracture nucleation through two important steps: (i) an external driving force term, included in the damage evolution equation, that accounts for the material strength; (ii) a properly designed damage function that defines the fluid pressure contribution on the crack driving force. The comparison of numerical results for two‐dimensional test cases with existing analytical solutions demonstrates that the proposed phase‐field model can accurately model both nucleation and propagation of hydraulic fractures. Additionally, we present the simulation of hydraulic fracturing in a three‐dimensional domain with various stress conditions to demonstrate the applicability of the method to realistic scenarios.

more » « less
Award ID(s):
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
International Journal for Numerical and Analytical Methods in Geomechanics
Medium: X Size: p. 3065-3089
["p. 3065-3089"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Despite its critical role in the study of earthquake processes, numerical simulation of the entire stages of fault rupture remains a formidable task. The main challenges in simulating a fault rupture process include the complex evolution of fault geometry, frictional contact, and off‐fault damage over a wide range of spatial and temporal scales. Here, we develop a phase‐field model for quasi‐dynamic fault nucleation, growth, and propagation, which features two standout advantages: (i) it does not require any sophisticated algorithms to represent fault geometry and its evolution; and (ii) it allows for modeling fault nucleation, propagation, and off‐fault damage processes with a single formulation. Built on a recently developed phase‐field framework for shear fractures with frictional contact, the proposed formulation incorporates rate‐ and state‐dependent friction, radiation damping, and their impacts on fault mechanics and off‐fault damage. We show that the numerical results of the phase‐field model are consistent with those obtained from well‐verified approaches that model the fault as a surface of discontinuity, without suffering from the mesh convergence issue in the existing continuous approaches to fault rupture (e.g., the stress glut method). Further, through numerical examples of fault propagation in various settings, we demonstrate that the phase‐field approach may open new opportunities for investigating complex earthquake processes that have remained overly challenging for the existing numerical methods.

    more » « less
  2. Abstract

    Fluid injection into rock formations can either produce complex branched hydraulic fractures, create simple planar fractures, or be dominated by porous diffusion. Currently, the optimum injection parameters to create branched fractures are unknown. We conducted repeatable hydraulic fracturing experiments using analog‐rock samples with controlled heterogeneity to quantify the fluid parameters that promote fracture branching. A large range of injection rates and fluid viscosities were used to investigate their effects on induced fracture patterns. Paired with a simple analytical model, our results identify the threshold at which fracture transitions from an isolated planar crack to branched cracks when closed natural fractures exist. These results demonstrate that this transition can be controlled by injection rate and fluid viscosity. In relation to the field practices, the present model predicts slickwater and lower viscosity fluid injections promote fracture branching, with the Marcellus shale used as an example.

    more » « less
  3. null (Ed.)
    Abstract. The continuum of behavior that emerges during fracturenetwork development in crystalline rock may be categorized into threeend-member modes: fracture nucleation, isolated fracture propagation, andfracture coalescence. These different modes of fracture growth producefracture networks with distinctive geometric attributes, such as clusteringand connectivity, that exert important controls on permeability and theextent of fluid–rock interactions. To track how these modes of fracturedevelopment vary in dominance throughout loading toward failure and thushow the geometric attributes of fracture networks may vary under theseconditions, we perform in situ X-ray tomography triaxial compressionexperiments on low-porosity crystalline rock (monzonite) under upper-crustalstress conditions. To examine the influence of pore fluid on the varyingdominance of the three modes of growth, we perform two experiments undernominally dry conditions and one under water-saturated conditions with 5 MPa ofpore fluid pressure. We impose a confining pressure of 20–35 MPa and thenincrease the differential stress in steps until the rock failsmacroscopically. After each stress step of 1–5 MPa we acquire athree-dimensional (3D) X-ray adsorption coefficient field from which weextract the 3D fracture network. We develop a novel method of trackingindividual fractures between subsequent tomographic scans that identifieswhether fractures grow from the coalescence and linkage of several fracturesor from the propagation of a single fracture. Throughout loading in all ofthe experiments, the volume of preexisting fractures is larger than that ofnucleating fractures, indicating that the growth of preexisting fracturesdominates the nucleation of new fractures. Throughout loading until close tomacroscopic failure in all of the experiments, the volume of coalescingfractures is smaller than the volume of propagating fractures, indicatingthat fracture propagation dominates coalescence. Immediately precedingfailure, however, the volume of coalescing fractures is at least double thevolume of propagating fractures in the experiments performed at nominallydry conditions. In the water-saturated sample, in contrast, although thevolume of coalescing fractures increases during the stage preceding failure,the volume of propagating fractures remains dominant. The influence ofstress corrosion cracking associated with hydration reactions at fracturetips and/or dilatant hardening may explain the observed difference infracture development under dry and water-saturated conditions. 
    more » « less
  4. Abstract

    Hydraulic fractures that grow in close proximity to one an other interact and compete for fluid that is injected to the wellbore, leading to dominance of some fractures and suppression of others. This phenomenon is ubiquitously encountered in stimulation of horizontal wells in the petroleum industry and it also bears possible relevance to emplacement of multiple laterally propagating swarms of magma‐driven dykes. Motivated by a need to validate mechanical models, this paper focuses on laboratory experiments and their comparison to simulation results for the behavior of multiple, simultaneously growing hydraulic fractures. The experiments entail the propagation of both uniformly and nonuniformly spaced hydraulic fractures by injection of glucose or glycerin‐based solutions into transparent (polymethyl methacrylate) blocks. Observed fracture growth is then compared to predictions of a fully coupled, parallel‐planar 3D hydraulic fracturing simulator. Results from experiments and simulations confirm the suppression of inner fractures when the spacing between the fractures is uniform. For certain non‐uniform spacing, both experiments and simulations show mitigated suppression of the central fractures. Specifically, the middle fracture in a 5‐fracture array grows nearly equally to the outer fractures from the beginning of injection. Furthermore, with some delay, the other two fractures that are suppressed with uniformly spaced configurations grow, and eventually achieve a velocity exceeding the other three fractures in the array. Hence, these experiments give the first laboratory evidence of a model‐predicted behavior wherein certain nonuniform fracture spacings result in drastic increases in the growth of all fractures within the array.

    more » « less
  5. Solids dispersion is an important part of hydraulic fracturing, both in helping to understand phenomena such as tip screen-out and spreading of the pad, and in new process variations such as cyclic pumping of proppant. Whereas many frac fluids have low viscosity, e.g. slickwater, others transport proppant through increased viscosity. In this context, one method for influencing both dispersion and solids-carrying capacity is to use a yield stress fluid as the frac fluid. We propose a model framework for this scenario and analyse one of the simplifications. A key effect of including a yield stress is to focus high shear rates near the fracture walls. In typical fracturing flows this results in a large variation in shear rates across the fracture. In using shear-thinning viscous frac fluids, flows may vary significantly on the particle scale, from Stokesian behaviour to inertial behaviour across the width of the fracture. Equally, according to the flow rates, Hele-Shaw style models give way at higher Reynolds number to those in which inertia must be considered. We develop a model framework able to include this range of flows, while still representing a significant simplification over fully three-dimensional computations. In relatively straight fractures and for fluids of moderate rheology, this simplifies into a one-dimensional model that predicts the solids concentration along a streamline within the fracture. We use this model to make estimates of the streamwise dispersion in various relevant scenarios. This model framework also predicts the transverse distributions of the solid volume fraction and velocity profiles as well as their evolutions along the flow part. 
    more » « less