scholarly journals Hydraulic Fracture Propagation Near the Cavity in a Poroelastic Media

2021 ◽  
Vol 11 (22) ◽  
pp. 11004
Author(s):  
Anton Azarov ◽  
Andrey Patutin ◽  
Sergey Serdyukov

In this paper, we investigate the problem of the propagation of hydraulic fractures in a poroelastic medium that has a circular cavity. The research was conducted using the extended finite element method (XFEM) implemented in the ABAQUS software package. The problem was considered in a plane formulation. The initial crack was oriented parallel to the surface of the cavity. It was shown that the path of the hydraulic fracture depends strongly on the hydrostatic stress in the medium and the distance between the crack and the cavity. We studied the influences of the poroelastic parameters, such as permeability and the Biot coefficient, on the propagation of cracks. It was shown that the cracks were less curved when the coupled problem of poroelasticity was considered. The features of fluid pressure changes inside the fracture and at the opening of the mouth were studied. It was shown that the fluid pressure in the fracture during injection was minimally sensitive to the state of the stress in the medium, to the position of the initial crack, and to the poroelastic parameters. The solution to the problem in this setting can be used to simulate hydraulic fracturing close to mine workings during a controlled roof’s collapse to prevent it from hanging, and the formation of impervious screens to reduce airflow from the mine to degassing boreholes through the rock, for example.

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.


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.


SPE Journal ◽  
2018 ◽  
Vol 24 (01) ◽  
pp. 302-318 ◽  
Author(s):  
Jixiang Huang ◽  
Joseph P. Morris ◽  
Pengcheng Fu ◽  
Randolph R. Settgast ◽  
Christopher S. Sherman ◽  
...  

Summary A fully coupled finite-element/finite-volume code is used to model 3D hydraulically driven fractures under the influence of strong vertical variations in closure stress interacting with natural fractures. Previously unknown 3D interaction mechanisms on fracture-height growth are revealed. Slipping of a natural fracture, triggered by elevated fluid pressure from an intersecting hydraulic fracture, can induce both increases and decreases of normal stress in the minimum-horizontal-stress direction, toward the center and tip of the natural fracture, respectively. Consequently, natural fractures are expected to be able to both encourage and inhibit the progress of hydraulic fractures propagating through stress barriers, depending on the relative locations between the intersecting fractures. Once the hydraulic fracture propagates above the stress barrier through the weakened segment near a favorably located natural fracture, a configuration consisting of two opposing fractures cuts the stress barrier from above and below. The fluid pressure required to break the stress barrier under such opposing-fracture configurations is substantially lower than that required by a fracture penetrating the same barrier from one side. Sensitivity studies of geologic conditions and operational parameters have also been performed to explore the feasibility of controlled fracture height. The interactions between hydraulic fractures, natural fractures, and geologic factors such as stress barriers in three dimensions are shown to be much more complex than in two dimensions. Although it is impossible to exhaust all the possible configurations, the ability of a 3D, fully coupled numerical model to naturally capture these processes is well-demonstrated.


SPE Journal ◽  
2016 ◽  
Vol 21 (04) ◽  
pp. 1302-1320 ◽  
Author(s):  
Mark W. McClure ◽  
Mohsen Babazadeh ◽  
Sogo Shiozawa ◽  
Jian Huang

Summary We developed a hydraulic-fracturing simulator that implicitly couples fluid flow with the stresses induced by fracture deformation in large, complex, 3D discrete-fracture networks (DFNs). The code is efficient enough to perform field-scale simulations of hydraulic fracturing in DFNs containing thousands of fractures, without relying on distributed-memory parallelization. The simulator can describe propagation of hydraulic fractures and opening and shear stimulation of natural fractures. Fracture elements can open or slide, depending on their stress state, fluid pressure, and mechanical properties. Fracture sliding occurs in the direction of maximum resolved shear stress. Nonlinear empirical equations are used to relate normal stress, fracture opening, and fracture sliding to fracture aperture and transmissivity. Fluid leakoff is treated with a semianalytical 1D leakoff model that accounts for changing pressure in the fracture over time. Fracture propagation is modeled with linear-elastic fracture mechanics. The Forchheimer equation (Forchheimer 1901) is used to simulate non-Darcy pressure drop in the fractures because of high flow rate. A crossing criterion is implemented that predicts whether propagating hydraulic fractures will cross natural fractures or terminate against them, depending on orientation and stress anisotropy. Height containment of propagating hydraulic fractures between bedding layers can be modeled with a vertically heterogeneous stress field or by explicitly imposing hydraulic-fracture-height containment as a model assumption. Limitations of the model are that all fractures must be vertical; the mechanical calculations assume a linearly elastic and homogeneous medium; proppant transport is not included; and the locations of potentially forming hydraulic fractures must be specified in advance. Simulations were performed of a single propagating hydraulic fracture with and without leakoff to validate the code against classical analytical solutions. Field-scale simulations were performed of hydraulic fracturing in a densely naturally fractured formation. The simulations demonstrate how interaction with natural fractures in the formation can help explain the high net pressures, relatively short fracture lengths, and broad regions of microseismicity that are often observed in the field during stimulation in low-permeability formations, and that are not predicted by classical hydraulic-fracturing models. Depending on input parameters, our simulations predicted a variety of stimulation behaviors, from long hydraulic fractures with minimal leakoff into surrounding fractures to broad regions of dense fracturing with a branching network of many natural and newly formed fractures.


Energies ◽  
2019 ◽  
Vol 12 (9) ◽  
pp. 1634 ◽  
Author(s):  
Juhyun Kim ◽  
Youngjin Seo ◽  
Jihoon Wang ◽  
Youngsoo Lee

Most shale gas reservoirs have extremely low permeability. Predicting their fluid transport characteristics is extremely difficult due to complex flow mechanisms between hydraulic fractures and the adjacent rock matrix. Recently, studies adopting the dynamic modeling approach have been proposed to investigate the shape of the flow regime between induced and natural fractures. In this study, a production history matching was performed on a shale gas reservoir in Canada’s Horn River basin. Hypocenters and densities of the microseismic signals were used to identify the hydraulic fracture distributions and the stimulated reservoir volume. In addition, the fracture width decreased because of fluid pressure reduction during production, which was integrated with the dynamic permeability change of the hydraulic fractures. We also incorporated the geometric change of hydraulic fractures to the 3D reservoir simulation model and established a new shale gas modeling procedure. Results demonstrate that the accuracy of the predictions for shale gas flow improved. We believe that this technique will enrich the community’s understanding of fluid flows in shale gas reservoirs.


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.


2018 ◽  
Vol 2018 ◽  
pp. 1-14 ◽  
Author(s):  
Yongliang Wang ◽  
Yang Ju ◽  
Yongming Yang

Hydrofracturing technology of perforated horizontal well has been widely used to stimulate the tight hydrocarbon reservoirs for gas production. To predict the hydraulic fracture propagation, the microseismicity can be used to infer hydraulic fractures state; by the effective numerical methods, microseismic events can be addressed from changes of the computed stresses. In numerical models, due to the challenges in accurately representing the complex structure of naturally fractured reservoir, the interaction between hydraulic and pre-existing fractures has not yet been considered and handled satisfactorily. To overcome these challenges, the adaptive finite element-discrete element method is used to refine mesh, effectively identify the fractures propagation, and investigate microseismic modelling. Numerical models are composed of hydraulic fractures, pre-existing fractures, and microscale pores, and the seepage analysis based on the Darcy’s law is used to determine fluid flow; then moment tensors in microseismicity are computed based on the computed stresses. Unfractured and naturally fractured models are compared to assess the influences of pre-existing fractures on hydrofracturing. The damaged and contact slip events were detected by the magnitudes, B-values, Hudson source type plots, and focal spheres.


Sign in / Sign up

Export Citation Format

Share Document