Simultaneous Multifracture Treatments: Fully Coupled Fluid Flow and Fracture Mechanics for Horizontal Wells

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.

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.


Geofluids ◽  
2018 ◽  
Vol 2018 ◽  
pp. 1-20 ◽  
Author(s):  
Quansheng Liu ◽  
Lei Sun ◽  
Pingli Liu ◽  
Lei Chen

Simultaneous multiple fracturing is a key technology to facilitate the production of shale oil/gas. When multiple hydraulic fractures propagate simultaneously, there is an interaction effect among these propagating hydraulic fractures, known as the stress-shadow effect, which has a significant impact on the fracture geometry. Understanding and controlling the propagation of simultaneous multiple hydraulic fractures and the interaction effects between multiple fractures are critical to optimizing oil/gas production. In this paper, the FDEM simulator and a fluid simulator are linked, named FDEM-Fluid, to handle hydromechanical-fracture coupling problems and investigate the simultaneous multiple hydraulic fracturing mechanism. The fractures propagation and the deformation of solid phase are solved by FDEM; meanwhile the fluid flow in the fractures is modeled using the principle of parallel-plate flow model. Several tests are carried out to validate the application of FDEM-Fluid in hydraulic fracturing simulation. Then, this FDEM-Fluid is used to investigate simultaneous multiple fractures treatment. Fractures repel each other when multiple fractures propagate from a single horizontal well, while the nearby fractures in different horizontal wells attract each other when multiple fractures propagate from multiple parallel horizontal wells. The in situ stress also has a significant impact on the fracture geometry.


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.


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.


2015 ◽  
Author(s):  
B.. Lecampion ◽  
J.. Desroches ◽  
X.. Weng ◽  
J.. Burghardt ◽  
J.E.. E. Brown

Abstract There is accepted evidence that multistage fracturing of horizontal wells in shale reservoirs results in significant production variation from perforation cluster to perforation cluster. Typically, between 30 and 40% of the clusters do not significantly contribute to production while the majority of the production comes from only 20 to 30% of the clusters. Based on numerical modeling, laboratory and field experiments, we investigate the process of simultaneously initiating and propagating several hydraulic fractures. In particular, we clarify the interplay between the impact of perforation friction and stress shadow on the stability of the propagation of multiple fractures. We show that a sufficiently large perforation pressure drop (limited entry) can counteract the stress interference between different growing fractures. We also discuss the robustness of the current design practices (cluster location, limited entry) in the presence of characterized stress heterogeneities. Laboratory experiments highlight the complexity of the fracture geometry in the near-wellbore region. Such complex fracture path results from local stress perturbations around the well and the perforations, as well as the rock fabric. The fracture complexity (i.e., the merging of multiple fractures and the reorientation towards the preferred far-field fracture plane) induces a strong nonlinear pressure drop on a scale of a few meters. Single entry field experiments in horizontal wells show that this near-wellbore effect is larger in magnitude than perforation friction and is highly variable between clusters, without being predictable. Through a combination of field measurements and modeling, we show that such variability results in a very heterogeneous slurry rate distribution; and therefore, proppant intake between clusters during a stage, even in the presence of limited entry techniques. We also note that the estimated distribution of proppant intake between clusters appears similar to published production log data. We conclude that understanding and accounting for the complex fracture geometry in the near-wellbore is an important missing link to better engineer horizontal well multistage completions.


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.


2005 ◽  
Vol 127 (3) ◽  
pp. 248-256 ◽  
Author(s):  
Hossein Jahediesfanjani ◽  
Faruk Civan

Coalbed methane (CBM) reservoirs are characterized as naturally fractured, dual porosity, low permeability, and water saturated gas reservoirs. Initially, the gas, water, and coal are at thermodynamic equilibrium under prevailing reservoir conditions. Dewatering is essential to promote gas production. This can be accomplished by suitable completion and stimulation techniques. This paper investigates the efficiency and performance of the openhole cavity, hydraulic fractures, frack and packs, and horizontal wells as potential completion methods which may reduce formation damage and increase the productivity in coalbed methane reservoirs. Considering the dual porosity nature of CBM reservoirs, numerical simulations have been carried out to determine the formation damage tolerance of each completion and stimulation approach. A new comparison parameter, named as the normalized productivity index Jnp(t) is defined as the ratio of the productivity index of a stimulated well to that of a nondamaged vertical well as a function of time. Typical scenarios have been considered to evaluate the CBM properties, including reservoir heterogeneity, anisotropy, and formation damage, for their effects on Jnp(t) over the production time. The results for each stimulation technique show that the value of Jnp(t) declines over the time of production with a rate which depends upon the applied technique and the prevailing reservoir conditions. The results also show that horizontal wells have the best performance if drilled orthogonal to the butt cleats. Long horizontal fractures improve reservoir productivity more than short vertical ones. Open-hole cavity completions outperform vertical fractures if the fracture conductivity is reduced by any damage process. When vertical permeability is much lower than horizontal permeability, production of vertical wells will improve while productivity of horizontal wells will decrease. Finally, pressure distribution of the reservoir under each scenario is strongly dependent upon the reservoir characteristics, including the hydraulic diffusivity of methane, and the porosity and permeability distributions in the reservoir.


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.


Sign in / Sign up

Export Citation Format

Share Document