Investigation of Rupture and Slip Mechanisms of Hydraulic Fractures in Multiple-Layered Formations

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 ◽  
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.


2021 ◽  
Author(s):  
Jiacheng Wang ◽  
Jon Olson

Abstract We propose an adaptive Eulerian-Lagrangian (E-L) proppant module and couple it with our simplified three-dimensional displacement discontinuity method (S3D DDM) hydraulic fracture model. The integrated model efficiently calculates proppant transport during three-dimensional (3D) hydraulic fracture propagation in multi-layer formations. The results demonstrate that hydraulic fracture height growth mitigates the form of proppant bed, so the proppant placement is more uniform in the hydraulic fracture under a smaller stress contrast. A higher fracturing fluid viscosity improves the suspension of proppant particles and generates a fracture larger in height and width but shorter in length. Lower proppant density and particle size reduce the proppant settling and create more uniform proppant placements, while they do not affect the hydraulic fracture geometry. Moreover, a larger proppant particle size limits the accessibility of the hydraulic fracture to the proppant, so the larger proppant particles do not fill the fracture tip and edge where the fracture width is small.


2021 ◽  
pp. 1-14
Author(s):  
Qian Gao ◽  
Ahmad Ghassemi

Summary The impacts of formation layering on hydraulic fracture containment and on pumping energy are critical factors in a successful stimulation treatment. Conventionally, it is considered that the in-situ stress is the dominant factor controlling the fracture height. The influence of mechanical properties on fracture height growth is often ignored or is limited to consideration of different Young’s moduli. Also, it is commonly assumed that the interfaces between different layers are perfectly bounded without slippage, and interface permeability is not considered. In-situ experiments have demonstrated that variation of modulus and in-situ stress alone cannot explain the containment of hydraulic fractures observed in the field (Warpinski et al. 1998). Enhanced toughness, in-situ stress, interface slip, and energy dissipation in the layered rocks should be combined to contribute to the fracture containment analysis. In this study, we consider these factors in a fully coupled 3D hydraulic fracture simulator developed based on the finite element method. We use laboratory and numerical simulations to investigate these factors and how they affect hydraulic fracture propagation, height growth, and injection pressure. The 3D fully coupled hydromechanical model uses a special zero-thickness interface element and the cohesive zone model (CZM) to simulate fracture propagation, interface slippage, and fluid flow in fractures. The nonlinear mechanical behavior of frictional sliding along interface surfaces is considered. The hydromechanical model has been verified successfully through benchmarked analytical solutions. The influence of layered Young’s modulus on fracture height growth in layered formations is analyzed. The formation interfaces between different layers are simulated explicitly through the use of the hydromechanical interface element. The impacts of mechanical and hydraulic properties of the formation interfaces on hydraulic fracture propagation are studied. Hydraulic fractures tend to propagate in the layer with lower Young’s modulus so that soft layers could potentially act as barriers to limit the height growth of hydraulic fractures. Contrary to the conventional view, the location of hydraulic fracturing (in softer vs. stiffer layers) does affect fracture geometry evolution. In addition, depending on the mechanical properties and the conductivity of the interfaces, the shear slippage and/or opening along the formation interfaces could result in flow along the interface surfaces and terminate the fracture growth. The frictional slippage along the interfaces can serve as an effective mechanism of containment of hydraulic fractures in layered formations. It is suggested that whether a hydraulic fracture would cross a discontinuity depends not only on the layers’ mechanical properties but also on the hydraulic properties of the discontinuity; both the frictional slippage and fluid pressure along horizontal formation interfaces contribute to the reinitiation of a hydraulic fracture from a pre-existing flaw along the interfaces, producing an offset from the interception point to the reinitiation point.


Energies ◽  
2021 ◽  
Vol 14 (9) ◽  
pp. 2610
Author(s):  
Wenzheng Liu ◽  
Qingdong Zeng ◽  
Jun Yao ◽  
Ziyou Liu ◽  
Tianliang Li ◽  
...  

Rock yielding may well take place during hydraulic fracturing in deep reservoirs. The prevailing models based on the linear elastic fracture mechanics (LEFM) are incapable of describing the evolution process of hydraulic fractures accurately. In this paper, a hydro-elasto-plastic model is proposed to investigate the hydraulic fracture propagation in deep reservoirs. The Drucker–Prager plasticity model, Darcy’s law, cubic law and cohesive zone model are employed to describe the plastic deformation, matrix flow, fracture flow and evolution of hydraulic fractures, respectively. Combining the embedded discrete fracture model (EDFM), extended finite element method (XFEM) and finite volume method, a hybrid numerical scheme is presented to carry out simulations. A dual-layer iterative procedure is developed based on the fixed-stress split method, Picard iterative method and Newton–Raphson iterative method. The iterative procedure is used to deal with the coupling between nonlinear deformation with fracture extension and fluid flow. The proposed model is verified against analytical solutions and other numerical simulation results. A series of numerical cases are performed to investigate the influences of rock plasticity, internal friction angle, dilatancy angle and permeability on hydraulic fracture propagation. Finally, the proposed model is extended to simulate multiple hydraulic fracture propagation. The result shows that plastic deformation can enhance the stress-shadowing effect.


2019 ◽  
Author(s):  
Jizhou Tang ◽  
Lihua Zuo ◽  
Lizhi Xiao ◽  
Kan Wu ◽  
Bin Qian ◽  
...  

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 ◽  
2020 ◽  
Vol 2020 ◽  
pp. 1-12 ◽  
Author(s):  
Song Yaobin ◽  
Lu Weiyong ◽  
He Changchun ◽  
Bai Erhu

According to the theory of plane mechanics involving the interaction of hydraulic and natural fractures, the law of hydraulic fracture propagation under the influence of natural fractures is verified using theoretical analysis and RFPA2D-Flow numerical simulation approaches. The shear and tensile failure mechanisms of rock are simultaneously considered. Furthermore, the effects of the approach angle, principal stress difference, tensile strength and length of the natural fracture, and elastic modulus and Poisson’s ratio of the reservoir on the propagation law of a hydraulic fracture are investigated. The following results are obtained: (1) The numerical results agree with the experimental data, indicating that the RFPA2D-Flow software can be used to examine the hydraulic fracture propagation process under the action of natural fractures. (2) In the case of a low principal stress difference and low approach angle, the hydraulic fracture likely causes shear failure along the tip of the natural fracture. However, under a high stress difference and high approach angle, the hydraulic fracture spreads directly through the natural fracture along the original direction. (3) When natural fractures with a low tensile strength encounter hydraulic fractures, the hydraulic fractures likely deviate and expand along the natural fractures. However, in the case of natural fractures with a high tensile strength, the natural fracture surface is closed, and the hydraulic fracture directly passes through the natural fracture, propagating along the direction of the maximum principal stress. (4) Under the same principal stress difference, a longer natural fracture corresponds to the easier initiation and expansion of a hydraulic fracture from the tip of the natural fracture. However, when the size of the natural fracture is small, the hydraulic fracture tends to propagate directly through the natural fracture. (5) A smaller elastic modulus and larger Poisson’s ratio of the reservoir result in a larger fracture initiation pressure. The presented findings can provide theoretical guidance regarding the hydraulic fracturing of reservoirs with natural fractures.


Sign in / Sign up

Export Citation Format

Share Document