TY - JOUR
T1 - Propagation mechanisms and parametric influence in multiple interacting hydraulic fractures
T2 - A 3-D G/XFEM hydro-mechanical modeling
AU - Mukhtar, Faisal M.
AU - Shauer, Nathan
AU - Duarte, C. Armando
N1 - F. M. Mukhtar and C. A. Duarte would like to acknowledge the financial support provided by the Deanship of Research Oversight and Coordination at King Fahd University of Petroleum & Minerals (KFUPM) under Research Grant SB191024. The support of the College of Petroleum Engineering & Geosciences at KFUPM in facilitating the use of computing resources is also gratefully acknowledged.
PY - 2022/8/10
Y1 - 2022/8/10
N2 - Using a robust three-dimensional generalized/eXtended finite element method (G/XFEM) algorithm, this paper presents a comprehensive study on multiple hydraulic fracture propagation and their interactions under different treatment conditions. Aimed at capturing the complex multiphysics behavior, the resulting nonplanar fracture footprints under mixed-mode conditions are simulated using a formulation, which couples the solid/rock domain equations and the fluid flow within the fractures/crevices. Due to the small length scale of the fracture process zone relative to the surrounding formation size, as typically encountered in tight shale reservoirs, linear elastic fracture mechanics based on a regularized Irwin criterion is adopted to describe the fractured solid/rock response while a power law is used for the fluid flow by assuming a Newtonian fluid behavior. Equipped with the capability for mesh adaptivity and automatic time step search algorithms, the G/XFEM utilized to discretize the resulting system of coupled nonlinear equations allows the adoption of independent meshes for the background solid domain and the fracture surfaces without any matching/compatibility requirement between the two. Fracture propagation directions are decided based on the modes I, II, and III stress intensity factors that are extracted using the displacement correlation method. The presented model is first applied to a pair of misaligned fractures and then to an array of en échelon fractures for qualitative verification against a literature prediction by the boundary element method (BEM) and an observed field behavior, respectively. Next, simulation of several sets and configurations of multiple hydraulic fractures, resulting in a total of 23 parametric studies, are carried out to investigate the influences of fracture spacing, injection fluid viscosity, number of fracture clusters, and the type of remote stress conditions.
AB - Using a robust three-dimensional generalized/eXtended finite element method (G/XFEM) algorithm, this paper presents a comprehensive study on multiple hydraulic fracture propagation and their interactions under different treatment conditions. Aimed at capturing the complex multiphysics behavior, the resulting nonplanar fracture footprints under mixed-mode conditions are simulated using a formulation, which couples the solid/rock domain equations and the fluid flow within the fractures/crevices. Due to the small length scale of the fracture process zone relative to the surrounding formation size, as typically encountered in tight shale reservoirs, linear elastic fracture mechanics based on a regularized Irwin criterion is adopted to describe the fractured solid/rock response while a power law is used for the fluid flow by assuming a Newtonian fluid behavior. Equipped with the capability for mesh adaptivity and automatic time step search algorithms, the G/XFEM utilized to discretize the resulting system of coupled nonlinear equations allows the adoption of independent meshes for the background solid domain and the fracture surfaces without any matching/compatibility requirement between the two. Fracture propagation directions are decided based on the modes I, II, and III stress intensity factors that are extracted using the displacement correlation method. The presented model is first applied to a pair of misaligned fractures and then to an array of en échelon fractures for qualitative verification against a literature prediction by the boundary element method (BEM) and an observed field behavior, respectively. Next, simulation of several sets and configurations of multiple hydraulic fractures, resulting in a total of 23 parametric studies, are carried out to investigate the influences of fracture spacing, injection fluid viscosity, number of fracture clusters, and the type of remote stress conditions.
KW - generalized/eXtended finite element method
KW - hydraulic fractures
KW - hydro-mechanical simulation
KW - mixed-mode fracture propagation
KW - multiple fracture interactions
UR - http://www.scopus.com/inward/record.url?scp=85130411696&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85130411696&partnerID=8YFLogxK
U2 - 10.1002/nag.3378
DO - 10.1002/nag.3378
M3 - Article
AN - SCOPUS:85130411696
SN - 0363-9061
VL - 46
SP - 2033
EP - 2059
JO - International Journal for Numerical and Analytical Methods in Geomechanics
JF - International Journal for Numerical and Analytical Methods in Geomechanics
IS - 11
ER -