Mechanisms of Simultaneous Hydraulic-Fracture Propagation From Multiple Perforation Clusters in Horizontal Wells

SPE Journal ◽  
2016 ◽  
Vol 21 (03) ◽  
pp. 1000-1008 ◽  
Author(s):  
Kan Wu ◽  
Jon E. Olson

Summary Simultaneous multiple-fracture treatments in horizontal wellbores are becoming a prevalent approach to economically develop unconventional resources in shale reservoirs. One challenge to efficiently use the technique is the generation of effective hydraulic fractures from all perforation clusters. In this work, we conducted a fundamental study of physical mechanisms controlling simultaneous multiple-fracture propagation and discussed the potential approaches to improve nonuniform development of multiple fractures. This study was investigated by our recently developed 3D fracture-propagation model that captures the coupled elastic deformation of the rock with fluid flow in the horizontal wellbore and within the fractures. The model demonstrated that fracture geometry was controlled by both the stress-shadow effects and dynamic partitioning of flow rate. The analysis results indicated that the nonuniform development of a multiple-fracture array, for example, a three-fracture array in this study, was induced by the uneven partitioning of flow rate into each fracture, which was dependent on the flow resistance from wellbore friction, perforation friction, and fracture propagation. Furthermore, the stress shadowing from the exterior fractures exerted additional stress on the interior fractures and increased the resistance of fracture propagation, resulting in the interior fractures receiving much less fluid. To minimize the negative effects of stress shadowing and favor more-uniform fracture growth, we investigated potential approaches to promote uniform partitioning of flow rate through adjusting the flow resistance between multiple fractures. The results showed that adjusting perforation friction can provide an effective way to modify the partitioning of flow rate and mitigate the negative effects of stress shadowing. The mechanisms investigated in this study are consistent with field observations. Our approach can help field operators to improve the effectiveness of multiple fracturing treatments and maximize the production.

SPE Journal ◽  
2014 ◽  
Vol 20 (02) ◽  
pp. 337-346 ◽  
Author(s):  
Kan Wu ◽  
Jon E. Olson

Summary Successfully creating multiple hydraulic fractures in horizontal wells is critical for unconventional gas production economically. Optimizing the stimulation of these wells will require models that can account for the simultaneous propagation of multiple, potentially nonplanar, fractures. In this paper, a novel fracture-propagation model (FPM) is described that can simulate multiple-hydraulic-fracture propagation from a horizontal wellbore. The model couples fracture deformation with fluid flow in the fractures and the horizontal wellbore. The displacement discontinuity method (DDM) is used to represent the mechanics of the fractures and their opening, including interaction effects between closely spaced fractures. Fluid flow in the fractures is determined by the lubrication theory. Frictional pressure drop in the wellbore and perforation zones is taken into account by applying Kirchoff's first and second laws. The fluid-flow rates and pressure compatibility are maintained between the wellbore and the multiple fractures with Newton's numerical method. The model generates physically realistic multiple-fracture geometries and nonplanar-fracture trajectories that are consistent with physical-laboratory results and inferences drawn from microseismic diagnostic interpretations. One can use the simulation results of the FPM for sensitivity analysis of in-situ and fracture treatment parameters for shale-gas stimulation design. They provide a physics-based complex fracture network that one can import into reservoir-simulation models for production analysis. Furthermore, the results from the model can highlight conditions under which restricted width occurs that could lead to proppant screenout.


Geofluids ◽  
2021 ◽  
Vol 2021 ◽  
pp. 1-18
Author(s):  
Zixi Jiao ◽  
Anlin Zhang ◽  
Longhuan Du ◽  
Yang Yang ◽  
Hua Fan

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 investigated, and the effects of engineering parameters and in situ conditions on the nonuniform development of multiple fractures are discussed. A 3D fracture propagation model was established with ABAQUS to show the influence of the stress shadow effects and dynamic partitioning of the flow rate by simulating the propagation of multiple competing fractures generated in the perforation clusters. Based on the results of these simulations, simultaneous flow in multiple fractures can propagate evenly. Through adjusting the number of perforations in each cluster or the perforation diameter, the effect of the stress shadow can be significantly reduced by increasing the perforation friction, and the factors that affect the development of multiple fractures are changed, from the stress shadow effect to the dynamic partitioning of the flow rate. When the stress shadow effect is dominant, increasing the fracturing fluid viscosity promotes the uniform development of multiple fractures and increases the fracture width. When the dynamic partitioning of the flow rate is dominant, increasing the injection rate greatly affects the uniform development of multiple fractures.


SPE Journal ◽  
2019 ◽  
Vol 24 (05) ◽  
pp. 2292-2307 ◽  
Author(s):  
Jizhou Tang ◽  
Kan Wu ◽  
Lihua Zuo ◽  
Lizhi Xiao ◽  
Sijie Sun ◽  
...  

Summary Weak bedding planes (BPs) that exist in many tight oil formations and shale–gas formations might strongly affect fracture–height growth during hydraulic–fracturing treatment. Few of the hydraulic–fracture–propagation models developed for unconventional reservoirs are capable of quantitatively estimating the fracture–height containment or predicting the fracture geometry under the influence of multiple BPs. In this paper, we introduce a coupled 3D hydraulic–fracture–propagation model considering the effects of BPs. In this model, a fully 3D displacement–discontinuity method (3D DDM) is used to model the rock deformation. The advantage of this approach is that it addresses both the mechanical interaction between hydraulic fractures and weak BPs in 3D space and the physical mechanism of slippage along weak BPs. Fluid flow governed by a finite–difference methodology considers the flow in both vertical fractures and opening BPs. An iterative algorithm is used to couple fluid flow and rock deformation. Comparison between the developed model and the Perkins–Kern–Nordgren (PKN) model showed good agreement. I–shaped fracture geometry and crossing–shaped fracture geometry were analyzed in this paper. From numerical investigations, we found that BPs cannot be opened if the difference between overburden stress and minimum horizontal stress is large and only shear displacements exist along the BPs, which damage the planes and thus greatly amplify their hydraulic conductivity. Moreover, sensitivity studies investigate the impact on fracture propagation of parameters such as pumping rate (PR), fluid viscosity, and Young's modulus (YM). We investigated the fracture width near the junction between a vertical fracture and the BPs, the latter including the tensile opening of BPs and shear–displacement discontinuities (SDDs) along them. SDDs along BPs increase at the beginning and then decrease at a distance from the junction. The width near the junctions, the opening of BPs, and SDDs along the planes are directly proportional to PR. Because viscosity increases, the width at a junction increases as do the SDDs. YM greatly influences the opening of BPs at a junction and the SDDs along the BPs. This model estimates the fracture–width distribution and the SDDs along the BPs near junctions between the fracture tip and BPs and enables the assessment of the PR required to ensure that the fracture width at junctions and along intersected BPs is sufficient for proppant transport.


2021 ◽  
Vol 18 (6) ◽  
pp. 954-969
Author(s):  
Yunlin Gao ◽  
Huiqing Liu ◽  
Chao Pu ◽  
Huiying Tang ◽  
Kun Yang ◽  
...  

Abstract To extract more gas from shale gas reservoirs, the spacing among hydraulic fractures should be made smaller, resulting in a significant stress shadow effect. Most studies regarding the stress shadow effect are based on the assumption of homogeneity in rock properties. However, strong heterogeneity has been observed in shale reservoirs, and the results obtained with homogeneous models can be different from practical situations. A series of case studies have been conducted in this work to understand the effects of mechanical heterogeneity on multiple fracture propagation. Fracture propagation was simulated using the extended finite element method. A sequential Gaussian simulation was performed to generate a heterogeneous distribution of geomechanical properties. According to the simulation results, the difficulty of fracture propagation is negatively correlated with the Young's modulus and Poisson's ratio, and positively correlated with tensile strength. When each of the multiple fractures propagates in a homogeneous area with different mechanical properties, the final geometry of the fracture is similar to homogeneous conditions. When the rock parameter is a random field or heterogeneity perpendicular to the propagation direction of fracture, the fracture will no longer take the wellbore as the center of symmetry. Based on the analysis of fracture propagation in random fields, a small variance of elastic parameters can result in asymmetrical propagation of multiple fractures. Moreover, the asymmetrical propagation of hydraulic fractures is more sensitive to the heterogeneity of Poisson's ratio than Young's modulus. This study emphasises the importance of considering geomechanical heterogeneity and provides some meaningful suggestions regarding hydraulic fracturing designs.


SPE Journal ◽  
2019 ◽  
Vol 24 (03) ◽  
pp. 1364-1377 ◽  
Author(s):  
Vyacheslav Guk ◽  
Mikhail Tuzovskiy ◽  
Don Wolcott ◽  
Joe Mach

Summary Horizontal wells with multiple hydraulic fractures have become a standard completion for the development of tight oil and gas reservoirs. Successful optimization of multiple-fracture design on horizontal wells began empirically in the Barnett Shale in the late 1990s (Steward 2013; Gertner 2013). More recently, research has focused on further improving fracturing performance by developing a model-derived optimum. Some researchers have focused on an economic optimum on the basis of multiple runs of an analytical or numerical model (Zhang et al. 2012; Saputelli et al. 2014). With such an approach, a new set of model runs is necessary to optimize the design each time the input parameters change significantly. Running multiple simulations for every optimization case might not always be practical. An alternative approach is to develop well-performance curves with dimensionless variables on the basis of the performance model. Such an approach was the basis for unified fracture design (UFD) for a single fracture in a vertical well (Economides et al. 2002). However, a similar systemized method to calculate the optimum for a horizontal well with multiple hydraulic fractures was missing. The objective of this study was to develop a rigorous and unified dimensionless optimization technique with type curves for the case of multiple transverse fractures in a horizontal well—an extension of UFD. The mathematical problem was solved in dimensionless variables. Multiple fractures include the proppant number (NP), penetration ratio (Ix), dimensionless conductivity (CfD), and aspect ratio (yeD) for each fracture, which is inversely proportional to the number of fractures. The direct boundary element (DBE) method was used to generate the dimensionless productivity index (JD) for a given range of these parameters (28,000 runs) for the pseudosteady-state case. Finally, total well JD was plotted as a function of the number of fractures for various NP. The effect of minimum fracture width was studied, and the optimization curves were adjusted for three cases of minimum fracture width. The provided dimensionless type curves can be used to identify the optimized number of fractures and their geometry for a given set of parameters, without running a more complicated numerical model multiple times. First, the proppant mass (and hence, NP) used for the fracture design can be selected on the basis of economic or other considerations. For this purpose, a relationship between total JD and NP, which accounts for the minimum fracture width requirement, was provided. Then, the optimal number of fractures can be calculated for a given NP using the generated type curves with minimum width constraints. The following observations were made during the study on the basis of the performed runs: For a given volume or proppant, NP, total JD for multiple fractures increases to an asymptote as the number of fractures increases. This asymptote represents a technical potential for multiple fractures and for high proppant numbers (NP≥100), with a technical potential of 3πNP. Below this asymptote, the more fractures that are created for a fixed NP, the larger the JD. In practice, minimum fracture width constrains the fracture geometry, and therefore maximum JD. For the case when 20/40 sand is used for multiple hydraulic fracturing of a 0.01-md formation with square total area, the optimal number of factures is approximately NP25. Application of horizontal drilling technology with multiple fractures assumes the availability of high proppant numbers. It was shown mathematically that the alternative low proppant numbers (NP≤20 for the previous case) are impractical for multiple fractures, because total JD cannot be significantly higher than JD for an optimized single fracture in the same area. This means that low formation permeability and/or high proppant volumes are needed for multiple fracture treatments.


Geofluids ◽  
2020 ◽  
Vol 2020 ◽  
pp. 1-16
Author(s):  
Chaoneng Zhao ◽  
Yongquan Hu ◽  
Jinzhou Zhao ◽  
Qiang Wang ◽  
Pei He ◽  
...  

The fracture propagation in hydraulic fracturing is described as a nonlinear problem dynamic boundary. Due to the limitation of mesh refinement, it is difficult to obtain the real crack propagation path using conventional numerical methods. Meshless methods (MMs) are an effective method to eliminate the dependence on the computational grid in the simulation of fracture propagation. In this paper, a hydraulic fracture propagation model is established based on the element-free Galerkin (EFG) method by introducing jump and branch enrichment functions. Based on the proposed method, three types of fracturing technology are investigated. The results reveal that the stress interference between fractures has an important impact on the propagation path. For the codirectional fracturing simultaneously, fractures propagate in a repel direction. However, the new fracture is attracted and eventually trapped by the adjacent fracture in the sequential fracturing case. For the opposite simultaneous fracturing in multiwells, two fractures with a certain lateral spacing will deflect toward each other. The effect of stress shadow should be used rationally in the optimization of construction parameters; for the single well multistage fracturing, the stage spacing should be out of stress inversion area, while for the simultaneous fracturing of multiple wells, stress inversion zones should be used to maximize communication between natural fractures. Overall, this study establishes a novel and effective approach of using MM to simulate the propagation of hydraulic fractures, which can serve as a useful reference for understanding the mechanism of hydraulic fracture propagation under various conditions.


Energies ◽  
2021 ◽  
Vol 14 (4) ◽  
pp. 894
Author(s):  
Tianwei Sun ◽  
Qingdong Zeng ◽  
Huilin Xing

In this study, a coupled thermo-hydro-mechanical model to simulate multiple hydraulic fracture propagation is presented. Fracture propagation with elastic deformation is described by using a displacement discontinuity method. The temperature distribution and induced thermal stress are calculated via a semi-analytical method in an explicit way. An iterative scheme is proposed to solve the coupling between fracture propagation with fluid flow and induced thermal stress. The numerical model is validated against related analytical solutions. Several numerical cases are modeled to investigate the controlling factors for uniform growth of multiple fractures. Results show that using non-uniform fracture spacings and proper increasing the spacing for fractures away from the heel of wellbore promote the uniform growth of multiple fractures by comparison with using uniform fracture spacings. Increasing the perforation diameter for the middle cluster also works. Besides, single-wing fracturing could greatly improve the uniform growth of multiple hydraulic fractures. Finally, it shows that the thermal stress has a significant influence on fracture geometrical size but has limited effect on fracture propagation path. In addition, the thermal effect promotes the uniform growth of multiple fractures.


Energies ◽  
2019 ◽  
Vol 12 (7) ◽  
pp. 1335 ◽  
Author(s):  
Jun Xie ◽  
Haoyong Huang ◽  
Yu Sang ◽  
Yu Fan ◽  
Juan Chen ◽  
...  

Recently, the Changning shale gas field has been one of the most outstanding shale plays in China for unconventional gas exploitation. Based on the more practical experience of hydraulic fracturing, the economic gas production from this field can be optimized and gradually improved. However, further optimization of the fracture design requires a deeper understanding of the effects of engineering parameters on simultaneous multiple fracture propagation. It can increase the effective fracture number and the well performance. In this paper, based on the Changning field data, a complex fracture propagation model was established. A series of case studies were investigated to analyze the effects of engineering parameters on simultaneous multiple fracture propagation. The fracture spacing, perforating number, injection rate, fluid viscosity and number of fractures within one stage were considered. The simulation results show that smaller fracture spacing implies stronger stress shadow effects, which significantly reduces the perforating efficiency. The perforating number is a critical parameter that has a big impact on the cluster efficiency. In addition, one cluster with a smaller perforating number can more easily generate a uniform fracture geometry. A higher injection rate is better for promoting uniform fluid volume distribution, with each cluster growing more evenly. An increasing fluid viscosity increases the variation of fluid distribution between perforation clusters, resulting in the increasing gap between the interior fracture and outer fractures. An increasing number of fractures within the stage increases the stress shadow among fractures, resulting in a larger total fracture length and a smaller average fracture width. This work provides key guidelines for improving the effectiveness of hydraulic fracture treatments.


SPE Journal ◽  
2020 ◽  
Vol 25 (02) ◽  
pp. 662-691 ◽  
Author(s):  
Ripudaman Manchanda ◽  
Shuang Zheng ◽  
Sho Hirose ◽  
Mukul M. Sharma

Summary This paper presents the formulation and results from a coupled finite-volume (FV)/finite-area (FA) model for simulating the propagation of multiple hydraulically driven fractures in two and three dimensions at the wellbore and pad scale. The proposed method captures realistic representations of local heterogeneities, layering, fracture turning, poroelasticity, interactions with other fractures, and proppant transport. We account for competitive fluid and proppant distribution between multiple fractures from the wellbore. Details of the model formulation and its efficient numerical implementation are provided, along with numerical studies comparing the model with both analytical solutions and field results. The results demonstrate the effectiveness of the proposed method for the comprehensive modeling of hydraulically driven fractures in three dimensions at a pad scale.


Sign in / Sign up

Export Citation Format

Share Document