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
-
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
-
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
-
Abstract Molecular dynamics (MD) has served as a powerful tool for designing materials with reduced reliance on laboratory testing. However, the use of MD directly to treat the deformation and failure of materials at the mesoscale is still largely beyond reach. In this work, we propose a learning framework to extract a peridynamics model as a mesoscale continuum surrogate from MD simulated material fracture data sets. Firstly, we develop a novel coarse-graining method, to automatically handle the material fracture and its corresponding discontinuities in the MD displacement data sets. Inspired by the weighted essentially non-oscillatory (WENO) scheme, the key idea lies at an adaptive procedure to automatically choose the locally smoothest stencil, then reconstruct the coarse-grained material displacement field as the piecewise smooth solutions containing discontinuities. Then, based on the coarse-grained MD data, a two-phase optimization-based learning approach is proposed to infer the optimal peridynamics model with damage criterion. In the first phase, we identify the optimal nonlocal kernel function from the data sets without material damage to capture the material stiffness properties. Then, in the second phase, the material damage criterion is learnt as a smoothed step function from the data with fractures. As a result, a peridynamics surrogate is obtained. As a continuum model, our peridynamics surrogate model can be employed in further prediction tasks with different grid resolutions from training, and hence allows for substantial reductions in computational cost compared with MD. We illustrate the efficacy of the proposed approach with several numerical tests for the dynamic crack propagation problem in a single-layer graphene. Our tests show that the proposed data-driven model is robust and generalizable, in the sense that it is capable of modeling the initialization and growth of fractures under discretization and loading settings that are different from the ones used during training.more » « less
