An official website of the United States government
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
A phase‐field model for hydraulic fracture nucleation and propagation in porous media
Abstract 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):
- 1933367
- PAR ID:
- 10482248
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- International Journal for Numerical and Analytical Methods in Geomechanics
- Volume:
- 47
- Issue:
- 16
- ISSN:
- 0363-9061
- Format(s):
- Medium: X Size: p. 3065-3089
- Size(s):
- p. 3065-3089
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
ABSTRACT This work proposes a novel approach for coupling non‐isothermal fluid dynamics with fracture mechanics to capture thermal effects within fluid‐filled fractures accurately. This method addresses critical aspects of calculating fracture width in enhanced geothermal systems, where the temperature effects of fractures are crucial. The proposed algorithm features an iterative coupling between an interface‐capturing phase‐field fracture method and interface‐tracking thermo‐fluid‐structure interaction using arbitrary Lagrangian–Eulerian coordinates. We use a phase‐field approach to represent fractures and reconstruct the geometry to frame a thermo‐fluid‐structure interaction problem, resulting in pressure and temperature fields that drive fracture propagation. We developed a novel phase‐field interface model accounting for thermal effects, enabling the coupling of quantities specific to the fluid‐filled fracture with the phase‐field model through the interface between the fracture and the intact solid domain. We provide several numerical examples to demonstrate the capabilities of the proposed algorithm. In particular, we analyze mesh convergence of our phase‐field interface model, investigate the effects of temperature on crack width and volume in a static regime, and highlight the method's potential for modeling slowly propagating fractures.more » « less
-
ABSTRACT:Creation of a fracture network in a hydraulic fracturing process is essential for subsurface energy extraction and CO2 sequestration. It is facilitated by reactivation of pre-existing intersecting weak layers and cemented cracks in the rock. In this study, a poromechanical model is developed for the hydraulic fracturing process in rocks containing such pre-existing weak layers. Based on the mixture theory, the crack band model is used to simulate the growth of a crack system. The governing equations with the parameters for hydromechanical coupling are derived, to describe the evolution of the opening and branching of cracks caused by water injection. Microplane model M7 is adopted to characterize the deformation and fracturing of the solid skeleton of the rock, and the Poiseuille law is used to characterize fluid flow through the hydraulic fractures. Numerical simulations are performed to reproduce and interpret recently published laboratory-scale hydraulic fracturing experiments conducted at Los Alamos National Laboratory (LANL). In these experiments, the rock was represented by confined plaster slabs containing orthogonal intersecting weak layers of higher porosity. Numerical simulations reveal how poromechanical characteristics such as the Biot coefficient and the fluid injection rate lead to various typical fracture modes observed in the experiments. These modes include formation of one dominant planar crack or various orthogonal fracture networks.more » « less
-
Abstract Krauklis waves are generated by pressure disturbances in fluid‐filled cavities and travel along the solid‐fluid interface. Their far‐field radiation, observed in seismic data from volcanoes or hydraulic fracturing, is known as long‐period events. Characterized by low velocity and resonance, Krauklis waves help estimate fracture size and discern fluids in saturated fractures. Despite numerous theoretical models analyzing Krauklis waves, the existing paradigms are founded on static flow conditions. However, in geological contexts, the assumption of static flow may not be valid. We developed an experimental apparatus using a tri‐layer model consisting of a pair of aluminum plates to examine the effect of fluid flow on Krauklis waves. We employed an infusion syringe pump to inject fluids into the fracture under different flow rates. We used water, oil, and an aqueous solution of Polyethylene glycol as fracture fluids. We calculated resonant frequency, phase velocity, and quality factor to characterize the Krauklis waves. Our findings reveal that an increase in flow rate leads to a higher phase velocity, higher quality factor, and a shift to higher resonant frequency when the flow is in the direction of initial wave propagation while decreasing amplitude. Additionally, when the flow is in the opposite direction of initial wave propagation, we note higher wave absorption and distortion of the Krauklis waves. Our observations unequivocally affirm that fluid flow leaves strong signatures on the Krauklis waves, providing a robust basis for characterizing fluid dynamics within geological settings through the analysis of Krauklis wave.more » « less
-
Davidson, Cam; Wirth, Karl (Ed.)Fault-tip damage zones develop in response to fault propagation and displacement and are caused by the local amplification of stresses at the fault tip. Understanding the geometry and intensity of damage zones is crucial for evaluating earthquake hazards and assessing the potentials of oil and gas production, geothermal energy, and groundwater resources. Fractures initiate as a result of stresses exceeding rock strength and propagate based on the stress field at the fault tip. We investigate the damage zone of a fault segment within the Sevier normal fault zone near Orderville, Utah, focusing on fractures that developed within the Jurassic Navajo Sandstone, the Temple Cap Formation, and the oldest beds of the Carmel Formation. Because normal faults grow laterally as slip and displacement increase, we focus on the tip zone of a fault segment where fracturing is well-exposed. We executed a series of unmanned aerial vehicle (UAV) flights to capture high-resolution imagery of inaccessible rock exposures. We use these images to construct structure-from-motion (SfM) virtual outcrop models (VOMs) that we georeference and analyze using Agisoft Metashape. We collected and analyzed fracture orientation and intensity data in the field and with VOMs. Both types of data reveal a distribution of fracture intensity that is consistent with inner and outer damage zones similar to previous studies of other fault systems. Adjacent to the tip, the inner damage zone has a higher fracture intensity on the hanging wall compared to the footwall. This high fracture intensity on the hanging wall ends 30 meters over from the fault core where the intensity of the outer damage zone of the hanging wall becomes similar to that within the inner damage zone of the footwall. Laterally, along strike of the fault tip, intense fracturing ends 60 meters to the south and all fracturing ends 350 meters from the fault tip. Our results have implications for the spatial distribution of fracturing and related permeability in similar normal fault systems.more » « less
