Three-dimensional numerical investigation of the interaction between multiple hydraulic fractures in horizontal wells
Introduction
Hydraulic fracturing has been extensively applied to the oil and gas industry [1], Enhanced Geothermal Systems [2], and the mining industry [3]. Complex, non-planar hydraulic fracture (HF) geometry is often observed in multi-stage hydraulic fracturing in horizontal wells. The interaction of multiple HFs is of major importance in the routine design of well completions and treatment strategies. The opening of an HF causes an alteration in the neighboring stress field, which in turn affects the propagation of the adjacent HFs. This phenomenon is referred to as the “stress shadow effect” [1]. Green and Sneddon [4] presented an analytical solution for a pressurized HF opening and the induced stress in an elastic solid. Based on Green and Sneddon's solution, an analysis of the fracture-induced stresses has been given to calculate the stability of reservoirs [5]. A three-dimensional stress solution produced by an elliptical crack has also been developed [6]. However, calculating stress variations induced by the interaction between multiple curved HFs using analytical solutions remains a challenging task.
Stress shadow has been documented for decades from fracture mapping data. Field observation has shown that the maximum increase in compressive stress can be expected to occur at the HF plane, with the stress disturbance radiating into the surrounding rock for hundreds of feet [1]. Laboratory experiments showed that the stress shadow effect is evident with the presence of closely spaced notches [7]. The stress shadow effect may either inhibit the propagation of neighboring HFs or induce a deviation of neighboring HFs. An older HF impact the shape and path of newer HF. The newer HF path curved toward and coalesced with the older HF [8]. Because of experimental limitations, it is difficult to perform a sensitivity study of the stress interference process for multiple fracture propagation under laboratory conditions.
Numerical simulation has been extensively used to study the interaction between multiple HFs. A comprehensive overview of numerical methods for hydraulic fracturing can be found in previous publications [9], [10], [11]. Only a few numerical methods are introduced in this study, including the finite element method (FEM), boundary element method (BEM), displacement discontinuity method (DDM), distinct element method (DEM), and lattice spring method. A finite-element model has been employed to study the interaction between multiple HFs [12]. The results showed that the stress induced by the opening of the HF or the fluid leak-off could alter the magnitudes and orientations of the principal stresses, thereby influencing the pathway of the HF. The 3D multiple-cluster fracturing problems have been investigated using the XFEM-based CZM (cohesive zone model based on the extended FEM) [13]. Results showed that the HF propagation pathway depends on a complex stress distribution (or stress shadowing effect). The interaction between HFs causes HF to coalesce, grow parallel, or diverge depending on cluster spacing. A pseudo-3D model based on the BEM has been utilized to investigate the simultaneous growth of multiple HFs [14]. The simulations showed that the HFs do not always propagate perpendicular to the far-field, minimum principal stress direction. The HF pattern complexity is influenced by the magnitude of the HF net pressure relative to the differential stress. The interaction between multiple neighboring cracks has been investigated using the BEM model [15]. Modeling results showed that the creation of outer cracks could act to either promote or inhibit adjacent cracks propagation depending on the spacing and the loading conditions. A plane strain simulator based on the DDM has been developed for simulating multiple HFs growth [16]. The simulations showed that the HFs can curve toward or away from one another, potentially intersecting. The curving of HFs is associated with a combination of opening and sliding along the previous HF, as well as the disturbance of the local stress field. An HF-propagation model based on the DDM has been presented [17]. Model results indicated that non-planar geometry reduced propagation occurs for interior HFs, resulting in a substantial restriction in HF aperture and a reduction in length. A pseudo 3D-based HF model based on the 3D-DDM has been developed [18]. Model results showed that the height growth may be promoted or suppressed for parallel HFs. The accurate prediction on the HF height and width profile require the incorporation of stress shadow effect.
The DEM has been extensively applied in different practices in rock engineering [19], [20], [21]. Simulations of stress interference among multiple HFs based on the DEM has been presented [22]. Simulation results showed that stress interference might impact the extent of the induced HF and treatment effectiveness. Numerical investigation of the stress shadow effect using the DEM code (both mechanical-only and hydro-mechanical coupled modes) has been presented [23]. The results showed that multiple HFs enhance the increase in minimum principal stress in the inter-fracture zone. Additionally, the increase in minimum principal stress reduces the HF apertures.
In the conventional DEM, discontinuities are regarded as distinct boundaries between blocks. The failure of intact blocks and the propagation of HFs may be simulated by implementing a Voronoi or Trigon tessellation within the intact blocks [24], [25]. Therefore, the HF trajectory is constrained by the geometry of the contacts. Based on Synthetic Rock Mass and Lattice methods, a fully-coupled code, XSite, has been presented, which allows simulating asymmetrical non-planar HF propagation in a full 3D setting [26]. The simulation of non-planar propagation of multiple HFs using XSite has been presented [27]. Simulation results revealed that the propagation of a middle HF is enhanced with an increase in the cluster spacing but decreases for higher in-situ stress difference. The mechanical interference among multiple HFs has been investigated using XSite [28]. The simulation captured the HF front segmentation and asymmetric growth due to the combined influence of the stress shadowing and perforation/in-situ stress misalignment.
In the present study, the lattice simulator, XSite, is employed to study the interference between multiple HFs. The simulator allows sensitivity analyses of the key parameters influencing multiple HFs propagation to be undertaken. Several simulations of simultaneous propagation of multiple HFs are presented, in which the effects of the various parameters (in-situ stresses, Young’s modulus (E), the viscosity of the fluid (μ), and injection rate (Q)) are investigated. Multiple HFs propagation in a layered formation is also simulated.
Section snippets
Modeling methodology
The lattice method is developed based on the DEM, with particles and contacts replaced by nodes and springs, respectively. As shown in Fig. 1, the lattice is a quasi-random array of nodes connected by normal and shear springs, which can fail in a brittle manner. The springs represent the elasticity of the rock mass. The fracturing of the intact material can be modeled by spring breakage. The joints can be inserted into the lattice spring network using a smooth joint model approach [26]. In
Model description
Fig. 3 depicts the setup of the base model, the volume of the rock block was assumed as 5 m × 4 m × 4 m. A vertical starter crack (2) is positioned at the center of the block, normal to the x-axis. The starter cracks (1), (3) are placed on the two sides of starter crack (2), and the distance between two adjacent cracks is 1.5 m. All starter crack has a radius of 0.2 m. Base model input parameters are listed in Table 1.
Effect of in-situ stress
Three types of stress regime were considered, i.e. σx = σy = σz = 5 MPa; σy = σz = 5 MPa, σx = 7.5 MPa; σx = σy = 5 MPa, σz = 7.5 MPa. Other parameters were assumed as: E = 20 GPa, μ = 1 mPa·s. Simulations were conducted under a constant injection rate of 0.001 m3/s for 2 s.
For the case of σx = σy = σz = 5 MPa (i.e., an isotropic stress field), in the very early stage of injection, three HFs perpendicular to the wellbore was initiated from the vertical starter cracks. The further propagation of
Discussion
A propagating HF alters the stress field near the HF, causing an increase in the local minimum principal stress. The increment of minimum horizontal stress may be greater than that of maximum horizontal stress, leading to a change in the direction of the maximum horizontal stress [17]. The induced stress field can have an impact on the further propagation of HFs. For the case of σx = σy = σz = 5 MPa (isotropic stress field), the propagating HF increased the compressive stress perpendicular to
Conclusion
In this paper, a lattice-spring code, XSite, was adopted to study the interference between multiple HFs. Simulations were conducted to account for the effects of in-situ stresses, Young’s modulus, the viscosity of the fluid, and the injection rate. The propagation of multiple HFs in a layered formation is also simulated. The following conclusions can be drawn.
- (1)
The interaction of multiple fractures can alter the local stress field, which in turn affects fracture propagation. Simultaneous
Declaration of Competing Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Acknowledgments
We thank Dr. Branko Damjanac at Itasca Consulting Group, Inc. for providing XSite used for the modeling and the guidance of the HF modeling. The first author wishes to thank the China Scholarship Council for their support for visiting Simon Fraser University. This work has been supported by the National Key Research and Development Program of China (grant no. 2017YFC0603003).
References (49)
- et al.
Natural fractures in soft coal seams and their effect on hydraulic fracture propagation: A field study
J Pet Sci Eng
(2020) - et al.
Laboratory scale research on the impact of stress shadow and natural fractures on fracture geometry during horizontal multi-staged fracturing in shale
Int J Rock Mech Min Sci
(2018) - et al.
Applying complex fracture model and integrated workflow in unconventional reservoirs
J Pet Sci Eng
(2014) - et al.
Overview of numerical models for interactions between hydraulic fractures and natural fractures: Challenges and limitations
Comput Geotech
(2016) - et al.
Numerical methods for hydraulic fracture propagation: A review of recent trends
J Nat Gas Sci Eng
(2018) - et al.
Finite element simulations of interactions between multiple hydraulic fractures in a poroelastic rock
Int J Rock Mech Min Sci
(2017) - et al.
Developments in the characterization of complex rock slope deformation and failure using numerical modelling techniques
Eng Geol
(2006) - et al.
Investigation of the influence of stress shadows on horizontal hydraulic fractures from adjacent lateral wells
J Unconv Oil Gas Resour
(2015) - et al.
Numerical simulation of strainbursts using a novel initiation method
Comput Geotech
(2019) - et al.
The application of a modified Voronoi logic to brittle fracture modelling at the laboratory and field scale
Int J Rock Mech Min Sci
(2014)
Numerical simulation of non-planar fracture propagation in multi-cluster fracturing with natural fractures based on Lattice methods
Eng Fract Mech
Initiation and growth of a hydraulic fracture from a circular wellbore
Int J Rock Mech Min Sci
Investigating the interaction of hydraulic fracture with pre-existing joints based on lattice spring modeling
Comput Geotech
Fractography Basics
Fractography Fail Anal Polym Elsevier
Stress shadow size and aperture of hydraulic fractures in unconventional shales
J Pet Sci Eng
The effect of stress distribution on the shape and direction of hydraulic fractures in layered media
Eng Fract Mech
Effect of discontinuity stress shadows on hydraulic fracture re-orientation
Int J Rock Mech Min Sci
Optimizing horizontal completion techniques in the Barnett shale using microseismic fracture mapping
SPE Annu Tech Conf Exhib Soc Petroleum Eng
Enhanced Geothermal Systems (EGS): Hydraulic fracturing in a thermo–poroelastic framework
J Pet Sci Eng
The distribution of stress in the neighbourhood of a flat elliptical crack in an elastic solid
Math Proc Cambridge Philos Soc
Analysis and prediction of microseismicity induced by hydraulic fracturing
SPE J
Experimental investigation of the interaction among closely spaced hydraulic fractures
45th US Rock Mech/Geomech Symp
XFEM-based CZM for the simulation of 3D multiple-cluster hydraulic fracturing in quasi-brittle shale formations
Rock Mech Rock Eng
Cited by (15)
Effects of heat treatment on mechanical properties and degradation behavior of hollow glass microsphere reinforced Mg-15Al-5Zn-1.5Cu composites
2023, Transactions of Nonferrous Metals Society of China (English Edition)Combined finite-discrete element modelling of hydraulic fracturing in reservoirs with filled joints
2023, Geoenergy Science and EngineeringStudy of presplit blasting under high in-situ stress
2023, Engineering Fracture MechanicsCombined finite-discrete element modelling of hydraulic fracturing in deep geologically complex reservoirs
2023, International Journal of Rock Mechanics and Mining SciencesHydraulic fracture development in conglomerate reservoirs simulated using combined finite-discrete element method
2023, Engineering Fracture MechanicsMechanical behaviors of coal measures and ground control technologies for China's deep coal mines – A review
2023, Journal of Rock Mechanics and Geotechnical EngineeringCitation Excerpt :Simultaneously propagating fractures can move toward or away from each other, which results in complicated non-planar geometry and branching of HFs, as shown in Fig. 33. An adjacent layer with a low modulus restricts the fracture height growth, whereas the contrast between lateral growth and stress interference will be relatively enhanced (Zhao et al., 2021). These findings further provide insight into the design and optimization of multi-stage HF in horizontal boreholes.