scholarly journals Numerical Simulation of Artificial Fracture Propagation in Shale Gas Reservoirs Based on FPS-Cohesive Finite Element Method

Geofluids ◽  
2019 ◽  
Vol 2019 ◽  
pp. 1-16 ◽  
Author(s):  
Xiaoqiang Liu ◽  
Zhanqing Qu ◽  
Tiankui Guo ◽  
Ying Sun ◽  
Zhifeng Shi ◽  
...  

The simulation of hydraulic fracturing by the conventional ABAQUS cohesive finite element method requires a preset fracture propagation path, which restricts its application to the hydraulic fracturing simulation of a naturally fractured reservoir under full coupling. Based on the further development of a cohesive finite element, a new dual-attribute element of pore fluid/stress element and cohesive element (PFS-Cohesive) method for a rock matrix is put forward to realize the simulation of an artificial fracture propagating along the arbitrary path. The effect of a single spontaneous fracture, two intersected natural fractures, and multiple intersected spontaneous fractures on the expansion of an artificial fracture is analyzed by this method. Numerical simulation results show that the in situ stress, approaching angle between the artificial fracture and natural fracture, and natural fracture cementation strength have a significant influence on the propagation morphology of the fracture. When two intersected natural fractures exist, the second one will inhibit the propagation of artificial fractures along the small angle of the first natural fractures. Under different in situ stress differences, the length as well as aperture of the hydraulic fracture in a rock matrix increases with the development of cementation superiority of natural fractures. And with the increasing of in situ horizontal stress differences, the length of the artificial fracture in a rock matrix decreases, while the aperture increases. The numerical simulation result of the influence of a single natural fracture on the propagation of an artificial fracture is in agreement with that of the experiment, which proves the accuracy of the PFS-Cohesive FEM for simulating hydraulic fracturing in shale gas reservoirs.

2020 ◽  
Vol 10 (8) ◽  
pp. 3319-3331 ◽  
Author(s):  
Belladonna Maulianda ◽  
Cindy Dhevayani Savitri ◽  
Aruvin Prakasan ◽  
Eziz Atdayev ◽  
Twon Wai Yan ◽  
...  

Abstract Hydraulic fracturing has been around for several decades since 1860s. It is one of the methods used to recover unconventional gas reservoirs. Hydraulic fracturing design is a challenging task due to the reservoir heterogeneity, complicated geological setting and in situ stress field. Hence, there are plenty of fracture modelling available to simulate the fracture initiation and propagation. The purpose of this paper is to provide a review on hydraulic fracturing modelling based on current hydraulic fracturing literature. Fundamental theory of hydraulic fracturing modelling is elaborated. Effort is made to cover the analytical and numerical modelling, while focusing on eXtended Finite Element Modelling (XFEM).


2020 ◽  
Vol 10 (8) ◽  
pp. 3333-3345
Author(s):  
Ali Al-Rubaie ◽  
Hisham Khaled Ben Mahmud

Abstract All reservoirs are fractured to some degree. Depending on the density, dimension, orientation and the cementation of natural fractures and the location where the hydraulic fracturing is done, preexisting natural fractures can impact hydraulic fracture propagation and the associated flow capacity. Understanding the interactions between hydraulic fracture and natural fractures is crucial in estimating fracture complexity, stimulated reservoir volume, drained reservoir volume and completion efficiency. However, because of the presence of natural fractures with diffuse penetration and different orientations, the operation is complicated in naturally fractured gas reservoirs. For this purpose, two numerical methods are proposed for simulating the hydraulic fracture in a naturally fractured gas reservoir. However, what hydraulic fracture looks like in the subsurface, especially in unconventional reservoirs, remain elusive, and many times, field observations contradict our common beliefs. In this study, the hydraulic fracture model is considered in terms of the state of tensions, on the interaction between the hydraulic fracture and the natural fracture (45°), and the effect of length and height of hydraulic fracture developed and how to distribute induced stress around the well. In order to determine the direction in which the hydraulic fracture is formed strikethrough, the finite difference method and the individual element for numerical solution are used and simulated. The results indicate that the optimum hydraulic fracture time was when the hydraulic fracture is able to connect natural fractures with large streams and connected to the well, and there is a fundamental difference between the tensile and shear opening. The analysis indicates that the growing hydraulic fracture, the tensile and shear stresses applied to the natural fracture.


1982 ◽  
Vol 22 (02) ◽  
pp. 209-218 ◽  
Author(s):  
Sunder H. Advani ◽  
J.K. Lee

Abstract Recently emphasis has been placed on the development and testing of innovative well stimulation techniques for the recovery of unconventional gas resources. The design of optimal hydraulic fracturing treatments for specified reservoir conditions requires sophisticated models for predicting the induced fracture geometry and interpreting governing mechanisms. This paper presents methodology and results pertinent to hydraulic fracture modeling for the U.S. DOE's Eastern Gas Shales Program (EGSP). The presented finite-element model simulations extend available modeling efforts and provide a unified framework for evaluation of fracture dimensions and associated responses. Examples illustrating the role of multilayering, in-situ stress, joint interaction, and branched cracks are given. Selected comparisons and applications also are discussed. Introduction Selection and design of stimulation treatments for Devonian shale wells has received considerable attention in recent years1-3. The production of natural gas from such tight eastern petroliferous basins is dependent on the vertical thickness of the organically rich shale matrix, its inherent fracture system density, anisotropy, and extent, and the communication-link characteristics of the induced fracture system(s). The investigation of stimulation techniques based on resource characterization, reservoir property evaluation, theoretical and laboratory model simulations, and field testing is a logical step toward the development of commercial technology for optimizing gas production and related costs. This paper reports formulations, methodology, and results associated with analytical simulations of hydraulic fracturing for EGSP. The presented model extends work reported by Perkins and Kern,4 Nordgren,5 Geertsma and DeKlerk,6 and Geertsma and Haafkens.7 The simulations provide a finite-element model framework for studying vertically induced fracture responses with the effects of multilayering and in-situ stress considered. In this context, Brechtel et al.,8 Daneshy,9 Cleary,10 and Anderson et al.11 have done recent studies addressing specific aspects of this problem. The use of finite-element model techniques for studying mixed-mode fracture problems encountered in dendritic fracturing and vertical fracture/joint interaction also is illustrated along with application of suitable failure criteria. Vertical Hydraulic Fracture Model Formulations Coupled structural fracture mechanics and fracture fluid response models for predicting hydraulically induced fracture responses have been reported previously.12,13 These simulations incorporate specified reservoir properties, in-situ stress conditions, and stimulation treatment parameters. One shortcoming of this modeling effort is that finite-element techniques are used for the structural and stress intensity simulations, while a finite-difference approach is used to evaluate the leakoff and fracture-fluid response in the vertical crack. A consistent framework for conducting all simulations using finite-element modeling is formulated here.


2002 ◽  
Vol 5 (03) ◽  
pp. 249-254 ◽  
Author(s):  
Colleen A. Barton ◽  
Mark D. Zoback

Summary Natural fractures and drilling-induced wellbore failures provide critical constraints on the state of in-situ stress and the direct applicability to problems of reservoir production, hydrocarbon migration, and wellbore stability. Acoustic, electrical, and optical wellbore images provide the means to detect and characterize natural fracture systems and to distinguish them from induced wellbore failures. We present new techniques and criteria to measure and characterize attributes of natural and induced fractures in borehole image data. These techniques are applied to the characterization of fracture permeability in two case studies. Introduction Wellbore image logs are extremely useful for identifying and studying a variety of modes of stress-induced wellbore failures. We present examples of how these wellbore failures appear in different types of image data and how they can be discriminated from natural fractures that intersect the wellbore. We then present brief overviews of two studies, which illustrate how the techniques have been applied to address specific issues of fracture permeability. Drilling-induced failures are ubiquitous in oil and gas and geothermal wells because the process of drilling a well causes a concentration of the far-field tectonic stress close to the wellbore, which often can exceed rock strength. Through the use of wellbore imaging and other logging techniques, stress-induced failures can be detected and categorized (compressive, tensile, or shear) and then used to estimate the unknown components of the stress field. We demonstrate how these modes of wellbore failures appear in different types of image data and the pitfalls in their interpretations. The most valuable use of drilling-induced features is to constrain the orientations and magnitudes of the current stress field. The use of drilling-induced features as stress indicators has become routine in the oil and gas industry.1–8 The detection of these features at the wellbore wall has become a primary target for Logging While Drilling/Measurement While Drilling (LWD/ MWD) real-time operations.9 A strong correlation between critically stressed fractures (fractures optimally oriented to the stress field for frictional failure) and hydraulic conductivity has been documented in a variety of reservoirs worldwide.10–12 When faults are critically stressed, permeabilities are increased, and the movement of fluid along faults is possible. We present examples of how knowledge of the stress state and natural fracture population may be used to access reservoir permeability. Drilling-Induced Tensile Wall Fractures Compressive and tensile failure of a wellbore is a direct result of the stress concentration around the wellbore, which results from drilling a well into an already stressed rock mass.13 Compressive wellbore failures (wellbore breakouts), first identified with caliper data, are useful for determining stress orientation in vertical wells.14–16 The study of such features with acoustic and electrical imaging devices makes it possible to clearly identify such features and to use them to determine stress magnitude and stress orientation.15,17–19 It is well known that if a wellbore is pressurized, a hydraulic fracture will form at the azimuth of the maximum horizontal stress.20 The formation of drilling-induced tensile wall fractures is the result of the natural stress state, perhaps aided by drilling-related perturbations, that causes the wellbore wall to fail in tension. The general case of tensile and compressive failure of arbitrarily inclined wellbores in different stress fields is described by Peska and Zoback,1 who demonstrate that there is a wide range of stress conditions under which drilling-induced tensile fractures occur in wellbores, even without a significant wellbore-fluid overpressure. We call these fractures tensile wall fractures because they occur only in the wellbore wall as a result of the stress concentration. These failures form in an orientation of the maximum principal horizontal stress in a vertical borehole (Fig. 1a) and as en echelon features in deviated wells (Fig. 1b). Because drilling-induced tensile wall fractures are very sensitive to the in-situ stress, they can be used to constrain the present state of stress.1,2,21–23 Pitfalls in Interpretation of Tensile Wall Fractures in Wellbore Image Data In cases in which drilling-induced tensile fractures form at an angle to the wellbore axis, it can be difficult to distinguish them from natural fractures, especially in electrical image logs that do not sample the entire wellbore circumference. Because misinterpretation of such features could lead to serious errors in the characterization of a fractured (or possibly not fractured!) reservoir, as well as the assessment of in-situ stress orientation and magnitude, we present criteria that are useful for discriminating natural from induced tensile fractures when observed in wellbore image logs. This is especially important because the wellbore stress concentration can have a significant effect on the appearance of natural fractures that intersect the wellbore. It is well known that fractures are mechanically weakened at their intersection with the borehole. This erosion causes the upper and lower peak and trough of the fracture sinusoid to be enlarged and subsequently enhanced in the standard 2D unwrapped view of wellbore image data (Fig. 2). Where the borehole hoop stress is tensile, the intersection of a natural fracture or foliation plane with the tensile region of the borehole may be preferentially opened in tension (Fig. 3a). These drilling-enhanced natural fractures can be mistaken easily for inclined tensile wellbore failures (Fig. 1b), thus resulting in serious errors in geomechanical modeling. Incipient wellbore breakouts are the early stages of wellbore breakout development, in which the borehole compressive stress concentration has exceeded the rock strength and initiated breakout development. The failed material within the breakout, however, has not yet spalled into the borehole (Fig. 3b). In a vertical borehole, these failures may appear as thin "fractures" that propagate vertically in the borehole and may be confused with drilling induced tensile wall cracks.


2015 ◽  
Vol 52 (7) ◽  
pp. 926-946 ◽  
Author(s):  
N. Zangeneh ◽  
E. Eberhardt ◽  
R.M. Bustin

Hydraulic fracturing is the primary means for enhancing rock mass permeability and improving well productivity in tight reservoir rocks. Significant advances have been made in hydraulic fracturing theory and the development of design simulators; however, these generally rely on continuum treatments of the rock mass. In situ, the geological conditions are much more complex, complicated by the presence of natural fractures and planes of weakness such as bedding planes, joints, and faults. Further complexity arises from the influence of the in situ stress field, which has its own heterogeneity. Together, these factors may either enhance or diminish the effectiveness of the hydraulic fracturing treatment and subsequent hydrocarbon production. Results are presented here from a series of two-dimensional (2-D) numerical experiments investigating the influence of natural fractures on the modeling of hydraulic fracture propagation. Distinct-element techniques applying a transient, coupled hydromechanical solution are evaluated with respect to their ability to account for both tensile rupture of intact rock in response to fluid injection and shear and dilation along existing joints. A Voronoi tessellation scheme is used to add the necessary degrees of freedom to model the propagation path of a hydraulically driven fracture. The analysis is carried out for several geometrical variants related to hypothetical geological scenarios simulating a naturally fractured shale gas reservoir. The results show that key interactions develop with the natural fractures that influence the size, orientation, and path of the hydraulic fracture as well as the stimulated volume. These interactions may also decrease the size and effectiveness of the stimulation by diverting the injected fluid and proppant and by limiting the extent of the hydraulic fracture.


2021 ◽  
Author(s):  
Mostafa Gorjian ◽  
Sepidehalsadat Hendi ◽  
Christopher D. Hawkes

Abstract. This paper presents selected results of a broader research project pertaining to the hydraulic fracturing of oil reservoirs hosted in the siltstones and fine grained sandstones of the Bakken Formation in southeast Saskatchewan, Canada. The Bakken Formation contains significant volumes of hydrocarbon, but large-scale hydraulic fracturing is required to achieve economic production rates. The performance of hydraulic fractures is strongly dependent on fracture attributes such as length and width, which in turn are dependent on in-situ stresses. This paper reviews methods for estimating changes to the in-situ stress field (stress shadow) resulting from mechanical effects (fracture opening), poro-elastic effects, and thermo-elastic effects associated with fluid injection for hydraulic fracturing. The application of this method is illustrated for a multi-stage hydraulic fracturing operation, to predict principal horizontal stress magnitudes and orientations at each stage. A methodology is also presented for using stress shadow models to assess the potential for inducing shear failure on natural fractures. The results obtained in this work suggest that thermo and poro-elastic stresses are negligible for hydraulic fracturing in the Bakken Formation of southeast Saskatchewan, hence a mechanical stress shadow formulation is used for analyzing multistage hydraulic fracture treatments. This formulation (and a simplified version of the formulation) predicts an increase in instantaneous shut-in pressure (ISIP) that is consistent with field observations (i.e., ISIP increasing from roughly 21.6 MPa to values slightly greater than 26 MPa) for a 30-stage fracture treatment. The size of predicted zones of shear failure on natural fractures are comparable with the event clouds observed in microseismic monitoring when assumed values of 115°/65° are used for natural fracture strike/dip; however, more data on natural fracture attributes and more microseismic monitoring data for the area are required before rigorous assessment of the model is possible.


2020 ◽  
Vol 2020 ◽  
pp. 1-9
Author(s):  
Peilun Li ◽  
Yan Dong ◽  
Sheng Wang ◽  
Peichao Li

Natural fractures usually develop in shale reservoirs. Thereby, in the process of hydraulic fracturing, it is inevitable that hydraulic fractures will intersect with natural fractures. In order to reveal the interaction mechanism between hydraulic-induced fractures and natural fractures, a two-dimensional fracture intersection model based on the extended finite element method (XFEM) is proposed, and the different types of intersecting criteria reported in the literature are compared. Then, the effects of natural fracture azimuth, fluid pressure in hydraulic fracture, and in situ principal stress difference on hydraulic fracturing are studied in detail. The results show that the fracture morphology is different under different criteria and working conditions. And the stress concentration phenomenon mainly concentrates on the tip in the obtuse angle side of natural fracture. Meanwhile, different fluid pressures in hydraulic fracture can also induce different intersection patterns. The obtained results in this work are of great benefit to understand the intersection mechanism between hydraulic fractures and natural fractures.


Sign in / Sign up

Export Citation Format

Share Document