Impact of the stress state and the natural network of fractures/faults on the efficiency of hydraulic fracturing operations in the Goldwyer Shale Formation

2020 ◽  
Vol 60 (1) ◽  
pp. 163
Author(s):  
Partha Pratim Mandal ◽  
Reza Rezaee ◽  
Joel Sarout

Cost-effective hydrocarbon production from low-permeability unconventional reservoirs requires multi-stage hydraulic fracturing (HF) operations. Each HF stage aims to generate the most spatially extended fracture network, giving access to the largest volume of reservoir possible (stimulated volume) and allowing hydrocarbons to flow towards the wellbore. The size of the stimulated volume, and therefore, the efficiency of any given HF stage, is governed by the rock’s deformational behaviour and presence of pre-existing natural fractures/faults. Naturally elevated pore pressures at depth not only help to reduce the injection energy required to generate hydraulic fractures but can also induce slip along pre-existing fractures/faults, and therefore, enhance production rates. Here we analyse borehole image, density, resistivity and sonic logs available from a vertical exploration well in the Goldwyer Shale Formation (Canning Basin) to (i) characterise the pre-existing network of natural fractures; and (ii) estimate the in-situ pore pressure and stress state at depth. The aim of such an analysis is to evaluate the possibility of fracture/fault reactivation (slip) during and following HF operations. Based on this analysis, we found that an increase in the formation's pore pressure by only a few MPa (typically ~5–10 MPa) could lead to slip along pre-existing fractures/faults, provided they are favourably oriented with respect to the prevalent stress field for future production. We also found that slip along the horizontal or sub-horizontal bedding of the Goldwyer Formation is unlikely in view of the prevalent strike-slip faulting regime, unless an extremely large overpressure exists within the reservoir.

Geofluids ◽  
2020 ◽  
Vol 2020 ◽  
pp. 1-13 ◽  
Author(s):  
Yintong Guo ◽  
Lei Wang ◽  
Xin Chang ◽  
Jun Zhou ◽  
Xiaoyu Zhang

Refracturing technology has become an important means for the regeneration of old wells reconstruction. It is of great significance to understand the formation mechanism of hydraulic fracturing fracture for the design of hydraulic fracturing. In order to accurately evaluate and improve fracturing volume after refracturing, it is necessary to understand the mechanism of refracturing fracture in shale formation. In this paper, a true triaxial refracturing test method was established. A series of large-scale true triaxial fracturing experiments were carried out to characterize the refracturing fracture initiation and propagation. The results show that for shale reservoirs with weak bedding planes and natural fractures, hydraulic fracturing can not only form the main fracturing fracture, which is perpendicular to horizontal minimum principal stress, but it can also open weak bedding plane or natural fractures. The characteristics of fracturing pump curve indicated that the evolution of fracturing fractures, including initiation and propagation and communication of multiple fractures. The violent fluctuation of fracturing pump pressure curve indicates that the sample has undergone multiple fracturing fractures. The result of refracturing shows that initial fracturing fracture channels can be effectively closed by temporary plugging. The refracturing breakdown pressure is generally slightly higher than that of initial fracturing. After temporary plugging, under the influence of stress induced by the initial fracturing fracture, the propagation path of the refracturing fracture is deviated. When the new fracturing fracture communicates with the initial fracturing fracture, the original fracturing fracture can continue to expand and extend, increasing the range of the fracturing modifications. The refracturing test results was shown that for shale reservoir with simple initial fracturing fractures, the complexity fracturing fracture can be increased by refracturing after temporary plugging initial fractures. The effect of refracturing is not obvious for the reservoir with complex initial fracturing fractures. This research results can provide a reliable basis for optimizing refracturing design in shale gas reservoir.


2020 ◽  
pp. 014459872096083
Author(s):  
Yulong Liu ◽  
Dazhen Tang ◽  
Hao Xu ◽  
Wei Hou ◽  
Xia Yan

Macrolithotypes control the pore-fracture distribution heterogeneity in coal, which impacts stimulation via hydrofracturing and coalbed methane (CBM) production in the reservoir. Here, the hydraulic fracture was evaluated using the microseismic signal behavior for each macrolithotype with microfracture imaging technology, and the impact of the macrolithotype on hydraulic fracture initiation and propagation was investigated systematically. The result showed that the propagation types of hydraulic fractures are controlled by the macrolithotype. Due to the well-developed natural fracture network, the fracture in the bright coal is more likely to form the “complex fracture network”, and the “simple” case often happens in the dull coal. The hydraulic fracture differences are likely to impact the permeability pathways and the well productivity appears to vary when developing different coal macrolithtypes. Thus, considering the difference of hydraulic fracture and permeability, the CBM productivity characteristics controlled by coal petrology were simulated by numerical simulation software, and the rationality of well pattern optimization factors for each coal macrolithotype was demonstrated. The results showed the square well pattern is more suitable for dull coal and semi-dull coal with undeveloped natural fractures, while diamond and rectangular well pattern is more suitable for semi-bright coal and bright coal with more developed natural fractures and more complex fracturing fracture network; the optimum wells spacing of bright coal and semi-bright coal is 300 m and 250 m, while that of semi-dull coal and dull coal is just 200 m.


SPE Journal ◽  
2013 ◽  
Vol 19 (03) ◽  
pp. 443-462 ◽  
Author(s):  
Sahar Ghannadi ◽  
Mazda Irani ◽  
Rick Chalaturnyk

Summary Inductive methods, such as electromagnetic steam-assisted gravity drainage (EM-SAGD), have been identified as technically and economically feasible recovery methods for shallow oil-sands reservoirs with overburdens of more than 30 m (Koolman et al. 2008). However, in EM-SAGD projects, the caprock overlying oil-sands reservoirs is also electromagnetically heated along with the bitumen reservoir. Because permeability is low in Alberta thermal-project caprock formations (i.e., the Clearwater shale formation in the Athabasca deposit and the Colorado shale formation in the Cold Lake deposit), the pore pressure resulting from the thermal expansion of pore fluids may not be balanced with the fluid loss caused by flow and the fluid-volume changes resulting from pore dilation. In extreme cases, the water boils, and the pore pressure increases dramatically as a result of the phase change in the water, which causes profound effective-stress reduction. After this condition is established, pore pressure increases can lead to shear failure of the caprock, the creation of microcracks and hydraulic fractures, and subsequent caprock integrity failure. It is typically believed that low-permeability caprocks impede the transmission of pore pressure from the reservoir, making them more resistant to shear failure (Collins 2005, 2007). In cases of induced thermal pressurization, low-permeability caprocks are not always more resistant. In this study, analytical solutions are obtained for temperature and pore-pressure rises caused by the constant EM heating rate of the caprock. These analytical solutions show that pore-pressure increases from EM heating depend on the permeability and compressibility of the caprock formation. For stiff or low-compressibility media, thermal pressurization can cause fluid pressures to approach hydrostatic pressure, and shear strength to approach zero for low-cohesive-strength units of the caprock (units of the caprock with high silt and sand percentage) and sections of the caprock with pre-existing fractures with no cohesion.


Geophysics ◽  
2020 ◽  
pp. 1-49
Author(s):  
Aleksei Titov ◽  
Gary Binder ◽  
Youfang Liu ◽  
Ge Jin ◽  
James Simmons ◽  
...  

Optimization of well spacings and completions are key topics in research related to the development of unconventional reservoirs. In 2017, a vertical seismic profiling (VSP) survey using fiber-optic-based distributed acoustic sensing (DAS) technology was acquired. The data include a series of VSP surveys taken before and immediately following the hydraulic fracturing of each of 78 stages. Scattered seismic waves associated with hydraulic fractures are observed in the seismic waveforms. Kinematic traveltime analysis and full-wavefield modeling results indicate these scattered events are converted PS-waves. We tested three different models of fracture-induced velocity inhomogeneities that can cause scattering of seismic waves: single hydraulic fracture, low-velocity zone, and tip diffractors. We compare the results with the field observations and conclude that the low-velocity zone model has the best fit for the data. In this model, the low-velocity zone represents a stimulated rock volume (SRV). We propose a new approach that uses PS-waves converted by SRV to estimate the half-height of the SRV and the closure time of hydraulic fractures. This active seismic source approach has the potential for cost-effective real-time monitoring of hydraulic fracturing operations and can provide critical constraints on the optimization of unconventional field development.


Author(s):  
Shike Zhang ◽  
Jincheng Lv ◽  
Xinsheng Yuan ◽  
Shunde Yin

Although lots of ways can be used to estimate geo-stress state, estimation of geo-stress state without knowing geomechanical parameters such as pore pressure, tensile strength and Poisson’s ratio, etc., still remains one of the most challenging tasks in geotechnical engineering. The main contribution of this paper is to present a back-propagation neural network (BPNN) with genetic algorithm (GA) optimization to predict the geo-stresses based on wellbore pressures of hydraulic fracturing tests during drilling. In the suggested hybrid model, the BPNN is used establish a mapping between the recording pressures and the geo-stress state. Also the GA is used to carry out the optimization of the weights and thresholds of BPNN model for improving accuracy of prediction. Finally, based on the record pressures in hydraulic fracturing (HF) tests, the BPNN model with genetic algorithm optimization successfully predicts the geo-stresses at the corresponding formation in the event that these parameters such as pore pressure, tensile strength and Poisson’s ratio are unavailable. In the meantime, the geo-stress state has been calculated using the theoretical formula by assuming pore pressure and tensile strength of rock mass are known. Then results from theoretical equation, BPNN and BPNN with GA optimization are compared, which shows that the degree accuracy of geo-stresses predicted by using GA-BPNN model is more obviously improvement than the predicted results by the basic BPNN model.


2020 ◽  
Vol 2020 ◽  
pp. 1-10
Author(s):  
Qingying Cheng ◽  
Bingxiang Huang ◽  
Xinglong Zhao

Rock directional fracturing is one of the difficult problems in deep mines. Directional fracturing controlled by hydraulic fracturing in dense linear multiboreholes is a novel directional fracturing technology of rock mass, which has been applied to the ground control in mines. In this paper, a physical model experiment was performed to study the fracture propagation process between multiboreholes. The results show that the intersecting of fractures between boreholes caused the sharp fluctuation of injecting water pressure. A directional fracturing plane was formed along with the direction of boreholes layout, and the surface of the fracturing plane is relatively flat. The dynamic initiation and propagation process of cracks between boreholes during directional hydraulic fracturing were simulated. The evolution of poroelastic stress and pore pressure between multiboreholes was analyzed. The numerical results indicated that a poroelastic stress concentration zone and pore pressure increase zone appeared between boreholes in the direction of boreholes layout. The pore pressure distribution is generally an elliptical seepage water pressure zone with the long axis along the direction of the boreholes layout. After the hydraulic fractures are initiated along the direction of the boreholes layout, the poroelastic stress on both sides of fractures decreases.


2015 ◽  
Author(s):  
Mark W. McClure ◽  
Mohsen Babazadeh ◽  
Sogo Shiozawa ◽  
Jian Huang

Abstract We developed a hydraulic fracturing simulator that implicitly couples fluid flow with the stresses induced by fracture deformation in large, complex, three-dimensional discrete fracture networks. 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 relations are used to relate normal stress, fracture opening, and fracture sliding to fracture aperture and transmissivity. Fluid leakoff is treated with a semianalytical one-dimensional leakoff model that accounts for changing pressure in the fracture over time. Fracture propagation is treated with linear elastic fracture mechanics. Non-Darcy pressure drop in the fractures due to high flow rate is simulated using Forchheimer's equation. 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. The code is efficient enough to perform field-scale simulations of hydraulic fracturing with a discrete fracture network containing thousands of fractures, using only a single compute node. 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 which 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.


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


Geofluids ◽  
2018 ◽  
Vol 2018 ◽  
pp. 1-23 ◽  
Author(s):  
Zhaohui Chong ◽  
Qiangling Yao ◽  
Xuehua Li

The presence of a significant amount of discontinuous joints results in the inhomogeneous nature of the shale reservoirs. The geometrical parameters of these joints exert effects on the propagation of a hydraulic fracture network in the hydraulic fracturing process. Therefore, mechanisms of fluid injection-induced fracture initiation and propagation in jointed reservoirs should be well understood to unleash the full potential of hydraulic fracturing. In this paper, a coupled hydromechanical model based on the discrete element method is developed to explore the effect of the geometrical parameters of the joints on the breakdown pressure, the number and proportion of hydraulic fractures, and the hydraulic fracture network pattern generated in shale reservoirs. The microparameters of the matrix and joint used in the shale reservoir model are calibrated through the physical experiment. The hydraulic parameters used in the model are validated through comparing the breakdown pressure derived from numerical modeling against that calculated from the theoretical equation. Sensitivity analysis is performed on the geometrical parameters of the joints. Results demonstrate that the HFN pattern resulting from hydraulic fracturing can be roughly divided into four types, i.e., crossing mode, tip-to-tip mode, step path mode, and opening mode. As β (joint orientation with respect to horizontal principal stress in plane) increases from 0° to 15° or 30°, the hydraulic fracture network pattern changes from tip-to-tip mode to crossing mode, followed by a gradual decrease in the breakdown pressure and the number of cracks. In this case, the hydraulic fracture network pattern is controlled by both γ (joint step angle) and β. When β is 45° or 60°, the crossing mode gains dominance, and the breakdown pressure and the number of cracks reach the lowest level. In this case, the HFN pattern is essentially dependent on β and d (joint spacing). As β reaches 75° or 90°, the step path mode is ubiquitous in all shale reservoirs, and the breakdown pressure and the number of the cracks both increase. In this case, β has a direct effect on the HFN pattern. In shale reservoirs with the same β, either decrease in k (joint persistency) and e (joint aperture) or increase in d leads to the increase in the breakdown pressure and the number of cracks. It is also found that changes in d and e result in the variation in the proportion of different types of hydraulic fractures. The opening mode of the hydraulic fracture network pattern is observed when e increases to 1.2 × 10−2 m.


Sign in / Sign up

Export Citation Format

Share Document