Numerical Modeling of Multistranded-Hydraulic-Fracture Propagation: Accounting for the Interaction Between Induced and Natural Fractures

SPE Journal ◽  
2011 ◽  
Vol 16 (03) ◽  
pp. 575-581 ◽  
Author(s):  
Arash Dahi-Taleghani ◽  
Jon E. Olson

Summary Recent examples of hydraulic-fracture diagnostic data suggest that complex, multistranded hydraulic-fracture geometry is a common occurrence. This reality is in stark contrast to the industry-standard design models based on the assumption of symmetric, planar, biwing geometry. The interaction between pre-existing natural fractures and the advancing hydraulic fracture is a key condition leading to complex fracture patterns. Performing hydraulic-fracture-design calculations under these less-than-ideal conditions requires modeling fracture intersections and tracking fluid fronts in the network of reactivated fissures. Whether a hydraulic fracture crosses or is arrested by a pre-existing natural fracture is controlled by shear strength and potential slippage at the fracture intersections, as well as potential debonding of sealed cracks in the near-tip region of a propagating hydraulic fracture. We present a complex hydraulic-fracture pattern propagation model based on the extended finite-element method (XFEM) as a design tool that can be used to optimize treatment parameters under complex propagation conditions. Results demonstrate that fracture-pattern complexity is strongly controlled by the magnitude of anisotropy of in-situ stresses, rock toughness, and natural-fracture cement strength, as well as the orientation of the natural fractures relative to the hydraulic fracture. Analysis shows that the growing hydraulic fracture may exert enough tensile and shear stresses on cemented natural fractures that the latter may be debonded, opened, and/or sheared in advance of hydraulic-fracture-tip arrival, while under other conditions, natural fractures will be unaffected by the hydraulic fracture. Detailed aperture distributions at the intersection between fracture segments show the potential for difficulty in proppant transport under complex fracture-propagation conditions.

2015 ◽  
Author(s):  
Hisanao Ouchi ◽  
Amit Katiyar ◽  
John T. Foster ◽  
Mukul M. Sharma

Abstract A novel fully coupled hydraulic fracturing model based on a nonlocal continuum theory of peridynamics is presented and applied to the fracture propagation problem. It is shown that this modeling approach provides an alternative to finite element and finite volume methods for solving poroelastic and fracture propagation problems and offers some clear advantages. In this paper we specifically investigate the interaction between a hydraulic fracture and natural fractures. Current hydraulic fracturing models remain limited in their ability to simulate the formation of non-planar, complex fracture networks. The peridynamics model presented here overcomes most of the limitations of existing models and provides a novel approach to simulate and understand the interaction between hydraulic fractures and natural fractures. The model predictions in two-dimensions have been validated by reproducing published experimental results where the interaction between a hydraulic fracture and a natural fracture is controlled by the principal stress contrast and the approach angle. A detailed parametric study involving poroelasticity and mechanical properties of the rock is performed to understand why a hydraulic fracture gets arrested or crosses a natural fracture. This analysis reveals that the poroelasticity, resulting from high fracture fluid leak-off, has a dominant influence on the interaction between a hydraulic fracture and a natural fracture. In addition, the fracture toughness of the rock, the toughness of the natural fracture, and the shear strength of the natural fracture also affect the interaction between a hydraulic fracture and a natural fracture. Finally, we investigate the interaction of multiple completing fractures with natural fractures in two-dimensions and demonstrate the applicability of the approach to simulate complex fracture networks on a field scale.


SPE Journal ◽  
2013 ◽  
Vol 19 (01) ◽  
pp. 161-171 ◽  
Author(s):  
Arash Dahi Taleghani ◽  
Jon E. Olson

Summary Hydraulic fracturing is recognized as the main stimulating technique to enhance recovery in tight fissured reservoirs. These fracturing treatments are often mapped by use of hypocenters of induced microseismic events. In some cases, the microseismic mapping shows asymmetry of the induced-fracture geometry with respect to the injection well. In addition, the conventional theories predict fracture propagation along a path normal to the least compressive in-situ stresses, whereas in some cases the microseismic data suggest fracture propagation parallel to the minimum compressive stress. In this paper, we present an extended-finite-element-method (XFEM) model that can simulate asymmetric fracture-wing development as well as diversion of the fracture path along natural fractures. Simulation results demonstrate the sensitivity of the fracture-pattern geometry to differential stress and natural-fracture orientation with respect to the in-situ maximum compressive stress. We examine the properties of sealed natural fractures that are common in formations such as the Barnett shale and show that they may still serve as weak paths for hydraulic-fracture beginning and/or diversion. The presented model predicts faster fracture propagation in formations where natural fractures are favorably aligned with the tectonic stresses.


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.


Processes ◽  
2018 ◽  
Vol 6 (8) ◽  
pp. 113 ◽  
Author(s):  
Shen Wang ◽  
Huamin Li ◽  
Dongyin Li

To investigate the mechanism of hydraulic fracture propagation in coal seams with discontinuous natural fractures, an innovative finite element meshing scheme for modeling hydraulic fracturing was proposed. Hydraulic fracture propagation and interaction with discontinuous natural fracture networks in coal seams were modeled based on the cohesive element method. The hydraulic fracture network characteristics, the growth process of the secondary hydraulic fractures, the pore pressure distribution and the variation of bottomhole pressure were analyzed. The improved cohesive element method, which considers the leak-off and seepage behaviors of fracturing liquid, is capable of modeling hydraulic fracturing in naturally fractured formations. The results indicate that under high stress difference conditions, the hydraulic fracture network is spindle-shaped, and shows a multi-level branch structure. The ratio of secondary fracture total length to main fracture total length was 2.11~3.62, suggesting that the secondary fractures are an important part of the hydraulic fracture network in coal seams. In deep coal seams, the break pressure of discontinuous natural fractures mainly depends on the in-situ stress field and the direction of natural fractures. The mechanism of hydraulic fracture propagation in deep coal seams is significantly different from that in hard and tight rock layers.


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.


Author(s):  
Yunsuk Hwang ◽  
Jiajing Lin ◽  
David Schechter ◽  
Ding Zhu

Multiple hydraulic fracture treatments in reservoirs with natural fractures create complex fracture networks. Predicting well performance in such a complex fracture network system is an extreme challenge. The statistical nature of natural fracture networks changes the flow characteristics from that of a single linear fracture. Simply using single linear fracture models for individual fractures, and then summing the flow from each fracture as the total flow rate for the network could introduce significant error. In this paper we present a semi-analytical model by a source method to estimate well performance in a complex fracture network system. The method simulates complex fracture systems in a more reasonable approach. The natural fracture system we used is fractal discrete fracture network model. We then added multiple dominating hydraulic fractures to the natural fracture system. Each of the hydraulic fractures is connected to the horizontal wellbore, and some of the natural fractures are connected to the hydraulic fractures through the network description. Each fracture, natural or hydraulically induced, is treated as a series of slab sources. The analytical solution of superposed slab sources provides the base of the approach, and the overall flow from each fracture and the effect between the fractures are modeled by applying the superposition principle to all of the fractures. The fluid inside the natural fractures flows into the hydraulic fractures, and the fluid of the hydraulic fracture from both the reservoir and the natural fractures flows to the wellbore. This paper also shows that non-Darcy flow effects have an impact on the performance of fractured horizontal wells. In hydraulic fracture calculation, non-Darcy flow can be treated as the reduction of permeability in the fracture to a considerably smaller effective permeability. The reduction is about 2% to 20%, due to non-Darcy flow that can result in a low rate. The semi-analytical solution presented can be used to efficiently calculate the flow rate of multistage-fractured wells. Examples are used to illustrate the application of the model to evaluate well performance in reservoirs that contain complex fracture networks.


Author(s):  
Chong Hyun Ahn ◽  
Robert Dilmore ◽  
John Yilin Wang

The most effective method for stimulating shale gas reservoirs is horizontal drilling with successful multi-stage hydraulic fracture treatments. Recent fracture diagnostic technologies have shown that complex fracture networks are commonly created in the field. The interaction between preexisting natural fractures and the propagating hydraulic fracture is a critical factor affecting the complex fracture network. However, many existing numerical models simulate only planar hydraulic fractures without considering the pre-existing fractures in the formation. The shale formations already contain a large number of natural fractures, so an accurate fracture propagation model needs to be developed to optimize the fracturing process. In this paper, we first understood the interaction between hydraulic and natural fractures. We then developed a new, coupled numerical model that integrates dynamic fracture propagation, reservoir flow simulation, and the interactions between hydraulic fractures and pre-existing natural fractures. By using the developed model, we conducted parametric studies to quantify the effects of rock toughness, stress anisotropy, and natural fracture spacing on the geometry and conductivities of the hydraulic fracture network. Lastly, we introduced new parmeters Fracture Network Index (FNI) and Width Anistropy (Wani) which may describe the characteristics of the fracture network in shale gas reservoirs. This new knowledge helps one understand and optimize the stimulation of shale gas reservoirs.


Energies ◽  
2019 ◽  
Vol 12 (23) ◽  
pp. 4477 ◽  
Author(s):  
Heng Zheng ◽  
Chunsheng Pu ◽  
CHOE TONG IL

Hydraulic fracturing is an essential technique for the development of shale gas, due to the low permeability in formation. Abundant natural fractures contained in a formation are indispensable for the development of a fracture network. In this paper, a damage-stress-seepage coupled hydraulic fracture expansion model, based on the extended finite element method, is established. The simulation results show that shear failure occurs when the hydraulic fracture interacts with a frictional natural fracture, while tensile failure occurs when it interacts with a cement natural fracture. Low interaction angles and high tensile strength of the rock are beneficial for the generation of a complex fracture network. Furthermore, under the same geological conditions and injection parameters, frictional natural fractures are more beneficial for the generation of a complex fracture network, when compared with cement natural fractures. This can not only effectively increase the propagation length of the natural fracture, but also effectively reduce its reactive resistance. This research is of great significance for the efficient exploitation of unconventional oil and gas resources.


2021 ◽  
Author(s):  
Ghazal Izadi ◽  
Colleen Barton ◽  
Pierre-Francois Roux ◽  
Tebis Llobet ◽  
Thiago Pessoa ◽  
...  

Abstract For tight reservoirs where hydraulic fracturing is required to enable sufficient fluid mobility for economic production, it is critical to understand the placement of induced fractures, their connectivity, extent, and interaction with natural fractures within the system. Hydraulic fracture initiation and propagation mechanisms are greatly influenced by the effect of the stress state, rock fabric and pre-existing features (e.g. natural fractures, faults, weak bedding/laminations). A pre-existing natural fracture system can dictate the mode, orientation and size of the hydraulic fracture network. A better understanding of the fracture growth phenomena will enhance productivity and also reduce the environmental footprint as less fractures can be created in a much more efficient way. Assessing the role of natural fractures and their interaction with hydraulic fractures in order to account for them in the hydraulic fracture model is achieved by leveraging microseismicity. In this study, we have used a combination of borehole and surface microseismic monitoring to get high vertical resolution locations and source mechanisms. 3D numerical modelling of hydraulic fracturing in complex geological conditions to predict fracture propagation is essential. 3D hydraulic fracturing simulation includes modelling capabilities of stimulation parameters, true 3D fracture propagation with near wellbore 3D complexity including a coupled DFN and the associated microseismic event generation capability. A 3D hydraulic fracture model was developed and validated by matching model predictions to microseismic observations. Microseismic source mechanisms are leveraged to determine the location and geometry of pre-existing features. In this study, we simulate a DFN based on the recorded seismicity of multi stage hydraulic fractures in a horizontal well. The advanced 3D hydraulic fracture modelling software can integrate effectively and efficiently data from a variety of multi-disciplinary sources and scales to create a subsurface characterization of the unconventional reservoir. By incorporating data from 3D seismic, LWD/wireline, core, completion/stimulation monitoring, and production, the software generates a holistic reservoir model embedded in a modular, multi-physics software platform of coupled numerical solvers that capture the fundamental physics of the processes being modelled. This study illustrates the importance of a powerful software tool that captures the necessary physics of stimulation to predict the effects of various completion designs and thereby ensure the most accurate representation of an unconventional reservoir response to a stimulation treatment.


Sign in / Sign up

Export Citation Format

Share Document