XFEM-Based CZM for the Simulation of 3D Multiple-Cluster Hydraulic Fracturing in Quasi-Brittle Shale Formations

Published on Aug 2, 2016in Rock Mechanics and Rock Engineering4.14
· DOI :10.1007/S00603-016-1057-2
Mahdi Haddad7
Estimated H-index: 7
(University of Texas at Austin),
Kamy Sepehrnoori67
Estimated H-index: 67
(University of Texas at Austin)
Sources
Abstract
The cohesive zone model (CZM) honors the softening effects and plastic zone at the fracture tip in a quasi-brittle rock, e.g., shale, which results in a more precise fracture geometry and pumping pressure compared to those from linear elastic fracture mechanics. Nevertheless, this model, namely the planar CZM, assumes a predefined surface on which the fractures propagate and therefore restricts the fracture propagation direction. Notably, this direction depends on the stress interactions between closely spaced fractures and can be acquired by integrating CZM as the segmental contact interaction model with a fully coupled pore pressure–displacement model based on extended finite element method (XFEM). This integrated model, called XFEM-based CZM, simulates the fracture initiation and propagation along an arbitrary, solution-dependent path. In this work, we modeled a single stage of 3D hydraulic fracturing initiating from three perforation clusters in a single-layer, quasi-brittle shale formation using planar CZM and XFEM-based CZM including slit flow and poroelasticity for fracture and matrix spaces, respectively, in Abaqus. We restricted the XFEM enrichment zones to the stimulation regions as enriching the whole domain leads to extremely high computational expenses and unrealistic fracture growths around sharp edges. Moreover, we validated our numerical technique by comparing the solution for a single fracture with KGD solution and demonstrated several precautionary measures in using XFEM in Abaqus for faster solution convergence, for instance the initial fracture length and mesh refinement. We demonstrated the significance of the injection rate and stress contrast in fracture aperture, injection pressure, and the propagation direction. Moreover, we showed the effect of the stress distribution on fracture propagation direction comparing the triple-cluster fracturing results from planar CZM with those from XFEM-based CZM. We found that the stress shadowing effect of hydraulic fractures on each other can cause these fractures to coalesce, grow parallel, or diverge depending on cluster spacing. We investigated the effect of this arbitrary propagation direction on not only the fractures’ length, aperture, and the required injection pressure, but also the fractures’ connection to the wellbore. This connection can be disrupted due to the near-wellbore fracture closure which may embed proppant grains on the fracture wall or screen out the fracture at early times. Our results verified that the near-wellbore fracture closure strongly depends on the following: (1) the implemented model, planar or XFEM-based CZM; and (2) fracture cluster spacing. Ultimately, we proposed the best fracturing scenario and cluster spacing to maintain the fractures connected to the wellbore.
📖 Papers frequently viewed together
1,475 Citations
65 Citations
124 Citations
References31
Newest
#1Mahdi Haddad (University of Texas at Austin)H-Index: 7
#2Jing DuH-Index: 2
Last. Sandrine Vidal-GilbertH-Index: 2
view all 3 authors...
10 CitationsSource
#1Mahdi Haddad (University of Texas at Austin)H-Index: 7
#2Kamy Sepehrnoori (University of Texas at Austin)H-Index: 67
Abstract Economic production from shale gas cannot be achieved by natural mechanisms alone; it requires technologies such as hydraulic fracturing in multiple stages along a horizontal wellbore. Developing numerical models for hydraulic fracturing is essential since a successful fracturing job in a shale formation cannot be generalized to another due to different shale characteristics, and restricted access to the field data acquisition. The cohesive zone model (CZM) identifies the plastic zone a...
65 CitationsSource
#1María González (LSU: Louisiana State University)H-Index: 5
#2A. Dahi Taleghani (LSU: Louisiana State University)H-Index: 11
Last. Jon E. Olson (University of Texas at Austin)H-Index: 45
view all 3 authors...
36 CitationsSource
#1Ali S. Ziarani (Trican Well Service)H-Index: 6
#2Cheney Chen (Trican Well Service)H-Index: 2
Last. Dana RoneyH-Index: 2
view all 5 authors...
6 CitationsSource
#1Mahdi Haddad (University of Texas at Austin)H-Index: 7
#2Kamy Sepehrnoori (University of Texas at Austin)H-Index: 67
18 CitationsSource
#1Pengcheng Fu (LLNL: Lawrence Livermore National Laboratory)H-Index: 19
#2Scott M. Johnson (LLNL: Lawrence Livermore National Laboratory)H-Index: 12
Last. Charles R. Carrigan (LLNL: Lawrence Livermore National Laboratory)H-Index: 24
view all 3 authors...
184 CitationsSource
#1T. Mohammadnejad (Sharif University of Technology)H-Index: 7
#2Amir R. Khoei (Sharif University of Technology)H-Index: 34
In this paper, a fully coupled numerical model is developed for the modeling of the hydraulic fracture propagation in porous media using the extended finite element method in conjunction with the cohesive crack model. The governing equations, which account for the coupling between various physical phenomena, are derived within the framework of the generalized Biot theory. The fluid flow within the fracture is modeled using the Darcy law, in which the fracture permeability is assumed according to...
212 CitationsSource
#1Kai Huang (A&M: Texas A&M University)H-Index: 10
#2Z. Zhang (A&M: Texas A&M University)H-Index: 2
Last. Ahmad Ghassemi (A&M: Texas A&M University)H-Index: 31
view all 3 authors...
SUMMARY Propagation of fractures, especially those emanating from wellbores and closed natural fractures, often involves Mode I and Mode II, and at times Mode III, posing significant challenges to its numerical simulation. When an embedded inclined fracture is subjected to compression, the fracture edge is constrained by the surrounding materials so that its true propagation pattern cannot be simulated by 2D models. In this article, a virtual multidimensional internal bond (VMIB) model is presen...
34 CitationsSource
#1Sau-Wai Wong (UCR: University of California, Riverside)H-Index: 10
#2Mikhail Geilikman (UCR: University of California, Riverside)H-Index: 1
Last. Guanshui Xu (UCR: University of California, Riverside)H-Index: 23
view all 3 authors...
The technology of multiple hydraulic fracture stimulation in horizontal wells has trans‐ formed the business of oil and gas exploitation from extremely tight, unconventional hydro‐ carbon bearing rock formations. The fracture stimulation process typically involves placing multiple fractures stage-by-stage along the horizontal well using diverse well completion technologies. The effective design of such massive fracture stimulation requires an under‐ standing of how multiple hydraulic fractures w...
18 CitationsSource
We present a method for the simulation of 3-D hydraulic fracturing in fully saturated porous media. The discrete fracture(s) is driven by the fluid pressure. A cohesive fracture model is adopted where the fracture follows the face of the elements around the fracture tip which is closest to the normal direction of the maximum principal stress at the fracture tip. No predetermined fracture path is needed. This requires continuous updating of the mesh around the crack tip to take into account the e...
89 CitationsSource
Cited By54
Newest
#1Bo Wang (China University of Petroleum)H-Index: 5
#2Fujian Zhou (China University of Petroleum)H-Index: 17
Last. Li Lizhe (China University of Petroleum)
view all 5 authors...
Abstract Near-wellbore diverting fracturing (NWDF) can successfully generate new diverting fractures and dramatically enhance the stimulation effects for vertical wells. Investigating the fracture geometry and the injection pressure response is necessary for the candidate well selection and the stimulation effect evaluation. This work establishes a fully fluid–solid coupling model to simulate the process of NWDF under various parameter combinations, including the horizontal stress contrast, Youn...
1 CitationsSource
This paper presents the first time implementation of the eXtended Finite Element Method (XFEM) in the general purpose commercial software COMSOL Multiphysics. An enrichment strategy is proposed, consistent with the structure of the software. To this end, for each set of enrichment functions, an additional Solid Mechanics module is incorporated into the numerical framework, coupled with compatible modifications to the internal variables. The Linear Elastic Fracture Mechanics (LEFM) is exclusively...
#1Sai Liu (SJTU: Shanghai Jiao Tong University)
#2Zhennan Zhang (SJTU: Shanghai Jiao Tong University)H-Index: 2
Last. Zhongwei Huang (China University of Petroleum)
view all 3 authors...
Abstract The numerical simulation of a 3D complex reservoir has been an important but difficult problem. One of the biggest challenges is how to deal with the large number of natural fractures with...
Source
#1Zixi Jiao (Sichuan University)
#2Anlin Zhang (Sichuan University)
Last. Hua Fan
view all 5 authors...
Simultaneous multiple-fracture treatments in horizontal wellbores have become one of the key methods for economically and efficiently developing oil and gas resources in unconventional reservoirs. However, field data show that some perforation clusters have difficulty propagating fractures due to the internal mechanism of competing initiation and propagation among the fractures. In this paper, the physical mechanisms that influence simultaneous multiple-fracture initiation and propagation are in...
Source
#1Bin Sun (SEU: Southeast University)H-Index: 8
Last. Zhaodong Xu (SEU: Southeast University)H-Index: 1
view all 3 authors...
Source
Source
#1Peng Tan (China University of Petroleum)H-Index: 10
#2Yan Jin (China University of Petroleum)H-Index: 49
Last. Huiwen Pang (China University of Petroleum)H-Index: 1
view all 3 authors...
Abstract The transition zone was defined as an interface with certain thickness, strong heterogeneity and anisotropy in the vertical profile for layered formations, which had great effect on hydraulic fracture (HF) height propagation. In this paper, a numerical model for transversely isotropic layered shale with transition zone was established by utilizing the extended finite element method (XFEM) based on cohesive zone model (CZM). The effects of in-situ stress, dip angle, anisotropy and tensil...
4 CitationsSource
#1Xian Shi (China University of Petroleum)
#2Yong Qin (PetroChina)H-Index: 2
Last. Songcai Han (China University of Petroleum)H-Index: 7
view all 7 authors...
Abstract Hydraulic fracturing is a necessary production enhancement procedure for low permeable conglomerate reservoirs. To study the hydraulic fracture (HF) propagation mechanism in strong heterogeneous conglomerate reservoirs, a numerical simulation based on the modified global cohesive zone method (CZM) was proposed. Two types of geological conglomerate models were built, and natural fractures were also considered in each model. A comprehensive discussion of the treatment parameters for the f...
Source
#1Chang Huang (LSU: Louisiana State University)H-Index: 1
#2Shengli Chen (LSU: Louisiana State University)H-Index: 15
1 CitationsSource
#1Fang ShiH-Index: 6
#2Jishan Liu (UWA: University of Western Australia)H-Index: 50
Abstract Field data suggest that hydraulic fractures under in situ conditions are three-dimensional and non-planar in shale gas reservoirs. In this study, a fully coupled deformation, fracture growth, and fluid flow model is developed to simulate this complex phenomenon with a focus on the reduction of the total computational effort. An extended finite element method with high-order tip enrichment functions is used to simulate the rock deformation while the finite element method for the laminar ...
1 CitationsSource