A Fast Calculation Method for Natural Fractures Activation Considering Stress Shadow and Study on the Law of Natural Fractures Activation State Changes

During hydraulic fracturing process of the Permian Basin in North America, the cluster spacing has been shortened to 3m, and stress shadow can no longer be ignored. Many scholars have studied the influence of stress shadows to optimize cluster spacing. For reservoirs with natural fractures, how to activate more natural fractures through hydraulic fracturing has become the purpose. However, few scholars have studied changes in the activation law of natural fractures under stress shadow conditions. This paper establishes stress change value around single fracture according to Sneddon formula, and calculates the maximum and minimum principal stress according to plane principal stress calculation formula. Considering attenuation of net pressure, stress field of multiple fractures is established, and influence of various factors on stress re-orientation is studied. Finally, considering attenuation of net pressure with distance, according to discriminant formulas of tension & shear activation, the proportion of natural fractures that are easily activated is calculated. By designing orthogonal experiments, the influence of different factors on the proportion of activated natural fractures was studied. The stress increase in the direction of the minimum principal stress is much greater than the increase in the direction of the maximum principal stress. The stress increases in the direction of the maximum principal stress at the tip of the hydraulic fracture. The tip position between hydraulic fractures is "neutralized" due to the superposition of shear stress. Stress-fracture angle and the in-situ stress difference are the common main influencing factors for both tensile and shear activation, but the net pressure has little effect on the tensile activation of natural fracture. The fracture spacing has little effect on the activation of natural fractures. When formulating the fracturing scheme, we should pay more attention to the net pressure rather than the fracture spacing. This article provides a fast calculation method for the activation state of natural fractures considering the stress shadow, which provides a reference index for activating more natural fractures and increasing the production of a single well.

Impact of geomechanical heterogeneity on multiple hydraulic fracture propagation

To extract more gas from shale gas reservoirs, the spacing among hydraulic fractures should be made smaller, resulting in a significant stress shadow effect. Most studies regarding the stress shadow effect are based on the assumption of homogeneity in rock properties. However, strong heterogeneity has been observed in shale reservoirs, and the results obtained with homogeneous models can be different from practical situations. A series of case studies have been conducted in this work to understand the effects of mechanical heterogeneity on multiple fracture propagation. Fracture propagation was simulated using the extended finite element method. A sequential Gaussian simulation was performed to generate a heterogeneous distribution of geomechanical properties. According to the simulation results, the difficulty of fracture propagation is negatively correlated with the Young's modulus and Poisson's ratio, and positively correlated with tensile strength. When each of the multiple fractures propagates in a homogeneous area with different mechanical properties, the final geometry of the fracture is similar to homogeneous conditions. When the rock parameter is a random field or heterogeneity perpendicular to the propagation direction of fracture, the fracture will no longer take the wellbore as the center of symmetry. Based on the analysis of fracture propagation in random fields, a small variance of elastic parameters can result in asymmetrical propagation of multiple fractures. Moreover, the asymmetrical propagation of hydraulic fractures is more sensitive to the heterogeneity of Poisson's ratio than Young's modulus. This study emphasises the importance of considering geomechanical heterogeneity and provides some meaningful suggestions regarding hydraulic fracturing designs.

Effect of Stress Shadow Caused by Multistage Fracturing from Multiple Well Pads on Fracture Initiation and Near-Wellbore Propagation from Infill Wells

Summary Multistage fracturing with multiwell pads (MSFMP) is an essential technology for the efficient development of unconventional oil and gas reservoirs, but the reservoir area between two well pads is often not stimulated. Fracture initiation and near-wellbore propagation from infill horizontal wells drilled with different azimuth from the optimal azimuth in the unstimulated area is poorly understood, largely because of the stress shadow (or induced stress) caused by MSFMP. In this study, we propose an integrated method for calculating the stress shadow caused by MSFMP and then determine optimal completion parameters for infill horizontal wells in the unstimulated connecting area between two well pads. First, we develop a theoretical stress shadow model caused by MSFMP on the basis of the dislocation theory. Considering two extreme cases, fully open and completely closed propped fractures, the range of stress shadow in the unstimulated area after MSFMP of 20 horizontal wells in Platform H of tight reservoirs in the Changqing Oilfield, China, is considered as an example. Second, we import the calculated stress shadow into a 3D perforated fracturing model that is built based on the discrete lattice method. Then, we investigate the influence of perforation technology, horizontal wellbore azimuth, phase angle, and injection rate on fracture initiation and near-wellbore propagation. Our results show that this model is capable of calculating stress shadow at any position and then can be used to optimize the fracturing interval for the middle unstimulated area. We find that appropriate perforation and fracturing parameters significantly decrease the complexity of near-wellbore fractures. The models and results presented in this paper provide a new method and new insight for quantifying and optimizing fracture initiation and propagation for infill horizontal wells to maximize reservoir stimulation efficiency.

Study on the Rock Fracture during Fluid Injection Using a Coupled Flow-Stress-Damage (FSD) Model: Insight into the Stress Shadow Effect

In order to investigate the influence of pore pressure on hydraulic fracturing behavior in the local and whole model, the coupled flow-stress-damage (FSD) analysis system RFPA-Flow was used to study the influence of rock heterogeneity, natural stress ratio, double-hole spacing, and water pressure gradient on the stress shadow effect. The numerical results show that the tensile crack induced by pore water pressure is significantly affected by the pore water pressure and water pressure gradient. The larger the pore pressure gradient is, the more asymmetrical the crack development pattern and the smaller the instability pressure of the model. In addition, the shape of hydraulic fracture becomes much more irregular with the increase in rock heterogeneity. The number and shape of tip microcracks under the influence of local water pressure are closely related to the homogeneity of rock. Moreover, when the natural stress difference is large, the hydraulic fracture propagates parallel to the maximum principal stress; when the stress field is close and the spacing between two holes is less than 5 times the diameter, the propagation direction of hydraulic fractures between holes is perpendicular to the maximum principal stress. It is found that no hydraulic fractures occur between the two holes when the distance between holes is greater than 5 times the diameter.

Numerical Simulation of Multi-Cluster Fracture Propagation In Horizontal Wells With Limited-Entry Designs

The effective propagation of multi-cluster fractures in horizontal wells is the key to the development of unconventional reservoirs. Due to the influence of pressure drops at perforating holes and the stress shadow effect, it is difficult to predict the fracturing fluid distribution and fracture dimensions in a fracturing stage. In this paper, a two-dimensional fluid-solid coupling model for simultaneous propagation of multiple fractures is established, and fluid distributions and dimensions of multiple fractures are studied with respect to different perforation designs. The model combines the User Amplitude Curve Subroutine (UAMP) in ABAQUS and the cohesive zone model (CZM) to calculate the perforating friction, fluid distribution and fracture propagation behaviors. After the accuracy of this model is verified by the analytical solution, a group of simulation is conducted to compare fracture propagations when the conventional limited-entry method (CLE) and extreme limited-entry method (less than 5 perforations per cluster, XLE) are used. Simulation results show that the edge and sub-central fractures in CLE cases almost get all the fluid and effectively propagate; central fractures receive little fluid and hardly propagate. In XLE cases, the fluid distribution of each fracture is relatively uniform, but the fracture lengths within one fracturing stage is still uneven; however, only reducing numbers or radii of perforation holes cannot achieve the uniform fracture propagation, where diverters might be further needed. Findings of this study provide a reference for the perforation optimization of multi-cluster horizontal wells in the field.

A Time-Continuation Solver for Hydraulic Fracture Propagation

We present an efficient time-continuation scheme for fluid-driven fracture propagation problems in the frame-work of the extended finite element method (XFEM). The fully coupled, fully implicit hydro-mechanical system is solved in conjunction with the linear elastic fracture propagation criterion by the Newton-Raphson method. Therefore, at the end of each time-step solve, the model ensures the energy release rate of weakest fracture tips within the equilibrium propagation regime. Besides, an initialization procedure for newly created fracture space as well as a priori estimate of stress intensity factor (SIF) growth rates are also developed to further improve the solver performance. We validate the model by the analytical solution and extend the problem to the multiple fracture propagation where stress shadow phenomenon occur.

Interaction Between Historical Earthquakes and the 2021 Mw7.4 Maduo Event and Their Impacts on the Seismic Gap Areas Along the East Kunlun Fault

On May 21, 2021 (UTC time), a Mw7.4 earthquake struck Maduo County, Qinghai Province, China. The rupture of this typical strike-slip event and its aftershocks along the Kunlun-Jiangcuo fault (JCF) propagated approximately 170 km from the epicenter. In this study, we calculated the coseismic and postseismic Coulomb stress changes induced by 14 historical earthquakes and investigated their impacts on the 2021 Maduo source area. We found that the JCF is in the stress shadow of these historical events with a combined ΔCFS range of approximately -0.4 to -0.2 MPa. Since the seismogenic fault of the 1937 event is nearly parallel and close to the JCF, the rupture of the 1937 event had the greatest inhibitory effect on Maduo source area. We hypothesize that the actual loading rate at the depth of the seismogenic layer in the Maduo source area is much higher than the simulated value (0.3 kPa/a). Consequently, the Maduo earthquake still occurred despite the considerable delaying effect of these historical earthquakes (especially the 1937 event). Our findings also indicate that the tectonic stress in the eastern Bayanhar block is still rapidly accumulating and adjusting. Our investigation further reveals the enhanced stress induced by the historical and Maduo events with ΔCFS values of approximately 30~300 kPa and 50~300 kPa on the XDS and the eastern end of the EKF, respectively, not only on the MMS but also at the eastern end of each branch segment of the EKF. Hence, considering the accumulation of tectonic stress, we suggest that the seismic hazard in these two regions has been promoted.

3D Extended Finite Element Modeling of Hydraulic Fracturing Design and Operation in Shale Gas Reservoir and Validation with Micro-Seismicity Data

Hydraulic fracturing is vital in unconventional shale gas development in order to produce economically from the reservoir. An optimum hydraulic fracturing design and operation can be the key difference between good and poor producing well and economics of the well. One of the most common hydraulic fracturing designs is ball drop system. Using ABAQUS software with XFEM method, a three layers model is used to represent overburden formation, shale gas formation and underburden formation. Rock properties, pore pressure and stress data are used as inputs for the generated model. A horizontal well is created in the middle shale gas formation with three fracture stages and 100m perforation spacing between them. Each hydraulic fracture stage is pressurized sequentially based on the treatment plan of ball drop sliding sleeve completion. The simulated hydraulic fractures are evaluated and compared with the measured field data. The comparison of the average wellbore pressure is good as they all showed within 10% of the measured data. The comparison of the hydraulic fracture geometry with the micro-seismicity data is reasonable overall in view of the data evaluation showing considerable uncertainties in the data. The hydraulic fracturing results also show that at 100m perforation spacing and using sequential hydraulic fracturing method (such as ball drop system), the effect of stress shadow is minimal and does not inhibit the fractures growth. However, the stress shadow effect is found to be pronounced for closer spacing between hydraulic fractures. For future application of the developed XFEM hydraulic fracturing model, it can be utilized to design new hydraulic fracturing completion in order to recommend the optimum completion, including perforation spacing, of development wells in unconventional shale gas field.

On the simultaneous growth of multiple hydraulic fractures emanating from an inclined well

The primary focus of this paper is to investigate the interaction between simultaneously propagating multiple fractures, initiated from an inclined well. In particular, the aim is to better understand the influence of the well inclination angle on the stress shadow between the fractures and on the overall resulting geometry of individual cracks. To simplify the analysis, the paper assumes the limit of large perforation friction, which leads to a uniform flux distribution between the fractures. The mathematical model for multiple hydraulic fractures is constructed by coupling together the respective models for individual fractures, each representing a single planar fracture model. In this approach, the fracture induced stress or stress shadow from a previous time step is used as an input for a given single hydraulic fracture to propagate independently. Further, to reduce computational burden, the effects associated with tangential stresses and displacements are neglected, whereby the stress interaction between the fractures is solely described by the normal opening and the normal stress component. Numerical results are presented for the storage viscosity dominated regime, whereby the effects of toughness and leak-off are negligible. An interesting behaviour is observed, demonstrating that the well inclination angle plays a significant role on the overall fracture symmetry. For zero inclination, all the fractures are nearly symmetrical and identical. However, once well inclination is introduced, this breaks the symmetry, making a profound effect on the final result.