Modeling of Hydraulic Fracture Propagation in Shale Gas Reservoirs: A Three-Dimensional, Two-Phase Model

2016 ◽  
Vol 139 (1) ◽  
Author(s):  
Chong Hyun Ahn ◽  
Robert Dilmore ◽  
John Yilin Wang
Keyword(s):  
Hydraulic Fracture ◽  
Three Dimensional ◽  
Fracture Network ◽  
Gas Reservoirs ◽  
Two Phase ◽  
Coupled Equations ◽  
Detailed Method ◽  
Network Propagation ◽  
Hydraulic Fracture Propagation ◽  
Fracture Height

A three-dimensional, two-phase, dual-continuum hydraulic fracture (HF) propagation simulator was developed and implemented. This paper presents a detailed method for efficient and effective modeling of the fluid flow within fracture and matrix as well as fluid leakoff, fracture height growth, and the fracture network propagation. Both a method for solving the system of coupled equations, and a verification of the developed model are presented herein.

Three-Dimensional Hydraulic Fracture Propagation in the Presence of Stress Variations

1983 ◽  
Vol 23 (06) ◽  
pp. 870-878 ◽  
Author(s):  
Ian D. Palmer ◽  
H.B. Carroll
Keyword(s):  
Boundary Condition ◽  
Hydraulic Fracture ◽  
Three Dimensional ◽  
Fracture Propagation ◽  
Rock Properties ◽  
Vertical Fracture ◽  
Hydraulic Fracture Propagation ◽  
Net Pressure ◽  
Fracture Height ◽  
Stress Contrast

Abstract Models of three-dimensional (3D) fracture propagation are being developed to study the effect of variations of stress and rock properties on fracture height and bottomhole pressure (BHP). Initially a blanket sand bounded by zones of higher minimum in-situ stress is considered, with stresses symmetrical about both the pay-zone axis and the wellbore. An elliptical fracture perimeter is assumed. Fluid flows are one-dimensional (1D) Newtonian in the direction of the pay zone. Two models, FL1 and FL2, are developed. In FL1, a discontinuous stress variation is approximated by a y2 variation in the vertical coordinate, and the fracture criterion, Ki = Kc, is satisfied at both major and minor axes. The net pressure at the tip, Lf, of the long axis required by the boundary condition Ki = Kc does not seem crucial in determining fracture height or BHP (compare with one group of published models that assumes p = 0 at Lf). Model FL2 properly represents the discontinuous stresses, and satisfies Ki = Kc at the wellbore but not at the tip of the long axis. A parametric study is made, with both models, of the comparative effects of stress contrast, Kc, pay-zone height, h, and Young's modulus, E, on fracture height and BHP. Results indicate that Kc does not have as much effect as either E or, at least for large stress contrasts. Model FL2 suggests the possibility of a rapid growth in fracture height as is reduced. Such modeling may be able to give an upper or "safe" limit on the pumping parameters ( and ) to ensure good containment. When the stress contrast is high, 700 psi [4826 kPa], an analytic derivation of BHP appears to be a good approximation for the parameters we use, if everywhere the fracture height is assumed equal to the pay zone height. Although leakoff is neglected here, subsequent modeling results show that, for leak off coefficients 0.001 ft- min [3.9 × 10 -5 m.s ], the results herein are a good approximation to the case when leak off is included. Introduction In their essence, models of hydraulic fracture propagation involve elasticity theory and fluid mechanics. The first is concerned with the fracture opening or width, w(p), as a function of net pressure on the fracture faces, while the second is concerned with the pressure drop, p(w), caused by the flow of viscous fluids in the fracture. Simultaneous solution of these equations includes a boundary condition that often takes the form Ki = Kc, where Ki is the stress-intensity factor at a point on the fracture tip, and Kc is the fracture toughness. The final solution is very complex in 3D, when a vertical fracture can expand vertically as well as horizontally along the pay zone. Thus, the first solutions were essentially two-dimensional (2D), and they assumed that the fracture height, hf, was fixed at the pay zone height, h. The 2D solutions were clustered in two groups as summarized by Nordgren, Perkins, and Geertsma and Haafkens. The first grouping, based on a model by Christianovich and Zheltov, assumed that the sides of an elongated, vertical fracture were parallel (i.e., free slippage between the pay and bounding zones, or no vertical stiffness). Other papers in this grouping included Geertsma and de Klerk, Daneshy and Settari. SPEJ P. 870^

No Need for Radio Active Tracer When Mother Nature Helps Detecting the Effective Hydraulic Fracture Height

2022 ◽  
Author(s):  
Rinat Lukmanov ◽  
Said Jabri ◽  
Ehab Ibrahim
Keyword(s):  
Hydraulic Fracture ◽  
Natural Radioactivity ◽  
Gamma Ray ◽  
Fracture Propagation ◽  
Radioactive Isotopes ◽  
Main Element ◽  
Tight Gas ◽  
Well Test ◽  
Gas Reservoirs ◽  
Fracture Height

Abstract The tight gas reservoirs of Haima Supergroup provide the majority of gas production in the Sultanate of Oman. The paper discusses a possibility of using the anomalies from natural radioactivity to evaluate the fracture height for complex tight gas in mature fields of Oman. The standard industry practice is adding radioactive isotopes to the proppant. Spectral Gamma Ray log is used to determine near wellbore traced proppant placement. Spectral Noise log in combination with Production logs helps to identify the active fractures contributing to production. These methods complement each other, but they are obviously associated with costs. Hence, majority of wells are fracced without tracers or any other fracture height diagnostics. However, in several brown fields, an alternative approach to identify fracture height has been developed which provides fit-for-purpose results. It is based on the analysis of naturally occurring radioactive minerals (NORM) precipitation. The anomalies were observed in the many gas reservoirs even in cases when tracers were not used. At certain conditions, these anomalies can be used to characterize fracture propagation and optimize future wells hydraulic Fracture design. A high number of PLTs and well test information were analyzed. Since tight formations normally don't produce without fracturing, radioactive anomalies flag the contributing intervals and hence fracture propagation. The main element of analysis procedure is related to that fact that if no tracers applied, the discrepancy between normalized Open Hole Gamma Ray and Gamma Ray taken during PLT after 6-12 months of production can be used instead to establish fracture height. This method cannot be applied for immediate interpretation of fracture propagation because time is required to precipitate NORM and using the anomalies concept. The advantage of this method is that it can be used in some fields to estimate the frac effectiveness of wells without artificial tracers. It is normally assumed that the Natural radioactivity anomalies appear mainly due to co-production of the formation water. However, in the fields of interest the anomalies appear in wells producing only gas and condensate. This observation provides an opportunity for active fracture height determination at minimum cost.

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 ◽  
...  
Keyword(s):  
Fluid Flow ◽  
Hydraulic Fracture ◽  
Fracture Propagation ◽  
Propagation Model ◽  
Rock Deformation ◽  
Hydraulic Fractures ◽  
Fracture Geometry ◽  
Fracture Width ◽  
Hydraulic Fracture Propagation ◽  
Fracture Height

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.

Sign in / Sign up

Export Citation Format

Share Document