scholarly journals Slickwater hydraulic fracture propagation: near-tip and radial geometry solutions

2019 ◽  
Vol 880 ◽  
pp. 514-550 ◽  
Author(s):  
Brice Lecampion ◽  
Haseeb Zia

We quantify the importance of turbulent flow on the propagation of hydraulic fractures (HF) accounting for the addition of friction reducing agents to the fracturing fluid (slickwater fluid). The addition in small quantities of a high molecular weight polymer to water is sufficient to drastically reduce friction of turbulent flow. The maximum drag reduction (MDR) asymptote is always reached during industrial-like injections. The energy required for pumping is thus drastically reduced, allowing for high volume high rate hydraulic fracturing operations at a reasonable cost. We investigate the propagation of a hydraulic fracture propagating in an elastic impermeable homogeneous solid under a constant (and possibly very high) injection rate accounting for laminar and turbulent flow conditions with or without the addition of friction reducers. We solve the near-tip HF problem and estimate the extent of the laminar boundary layer near the fracture tip as a function of a tip Reynolds number for slickwater. We obtain different propagation scalings and transition time scales. This allows us to easily quantify the growth of a radial HF from the early-time turbulent regime(s) to the late-time laminar regimes. Depending on the material and injection parameters, some propagation regimes may actually be bypassed. We derive both accurate and approximate solutions for the growth of radial HF in the different limiting flow regimes (turbulent smooth, rough, MDR) for the zero fracture toughness limit (corresponding to the early stage of propagation of a radial HF). We also investigate numerically the transition(s) between the early-time MDR regime to the late-time laminar regimes (viscosity and toughness) for slickwater fluid. Our results indicate that the effect of turbulent flow on high rate slickwater HF propagation is limited and matters only at early times (at most during the first minutes for industrial hydraulic fracturing operations).

2021 ◽  
pp. 014459872198899
Author(s):  
Weiyong Lu ◽  
Changchun He

Directional rupture is one of the most important and most common problems related to rock breaking. The goal of directional rock breaking can be effectively achieved via multi-hole linear co-directional hydraulic fracturing. In this paper, the XSite software was utilized to verify the experimental results of multi-hole linear co-directional hydraulic fracturing., and its basic law is studied. The results indicate that the process of multi-hole linear co-directional hydraulic fracturing can be divided into four stages: water injection boost, hydraulic fracture initiation, and the unstable and stable propagation of hydraulic fracture. The stable expansion stage lasts longer and produces more microcracks than the unstable expansion stage. Due to the existence of the borehole-sealing device, the three-dimensional hydraulic fracture first initiates and expands along the axial direction in the bare borehole section, then extends along the axial direction in the non-bare hole section and finally expands along the axial direction in the rock mass without the borehole. The network formed by hydraulic fracture in rock is not a pure plane, but rather a curved spatial surface. The curved spatial surface passes through both the centre of the borehole and the axial direction relative to the borehole. Due to the boundary effect, the curved spatial surface goes toward the plane in which the maximum principal stress occurs. The local ground stress field is changed due to the initiation and propagation of hydraulic fractures. The propagation direction of the fractures between the fracturing boreholes will be deflected. A fracture propagation pressure that is greater than the minimum principle stress and a tension field that is induced in the leading edge of the fracture end, will aid to fracture intersection; as a result, the possibility of connecting the boreholes will increase.


2015 ◽  
Author(s):  
Manhal Sirat ◽  
Mujahed Ahmed ◽  
Xing Zhang

Abstract In-situ stress state plays an important role in controlling fracture growth and containment in hydraulic fracturing managements. It is evident that the mechanical properties, existing stress regime and the natural fracture network of its reservoir rocks and the surrounding formations mainly control the geometry, size and containments of produced hydraulic fractures. Furthermore, the three principal in situ stresses' axes swap directions and magnitudes at different depths giving rise to identifying different mechanical bedrocks with corresponding stress regimes at different depths. Hence predicting the hydro-fractures can be theoretically achieved once all the above data are available. This is particularly difficult in unconventional and tight carbonate reservoirs, where heterogeneity and highly stress variation, in terms of magnitude and orientation, are expected. To optimize the field development plan (FDP) of a tight carbonate gas reservoir in Abu Dhabi, 1D Mechanical Earth Models (MEMs), involving generating the three principal in-situ stresses' profiles and mechanical property characterization with depth, have been constructed for four vertical wells. The results reveal the swap of stress magnitudes at different mechanical layers, which controls the dimension and orientation of the produced hydro-fractures. Predicted containment of the Hydro-fractures within the specific zones is likely with inevitable high uncertainty when the stress contrast between Sv, SHmax with Shmin respectively as well as Young's modulus and Poisson's Ratio variations cannot be estimated accurately. The uncertainty associated with this analysis is mainly related to the lacking of the calibration of the stress profiles of the 1D MEMs with minifrac and/or XLOT data, and both mechanical and elastic properties with rock mechanic testing results. This study investigates the uncertainty in predicting hydraulic fracture containment due to lacking such calibration, which highlights that a complete suite of data, including calibration of 1D MEMs, is crucial in hydraulic fracture treatment.


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.


2021 ◽  
Author(s):  
Debotyam Maity ◽  
Jordan Ciezobka

Abstract In this case study, we apply a novel fracture imaging and interpretation workflow to take a systematic look at hydraulic fractures captured during thorugh fracture coring at the Hydraulic Fracturing Test Site (HFTS) in Midland Basin. Digital fracture maps rendered using high resolution 3D laser scans are analyzed for fracture morphology and roughness. Analysis of hydraulic fracture faces show that the roughness varies systematically in clusters with average cluster separation of approximately 20' along the core. While isolated smooth hydraulic fractures are observed in the dataset, very rough fractures are found to be accompanied by proximal smoother fractures. Roughness distribution also helps understand the effect of stresses on fracture distribution. Locally, fracture roughness seems to vary with fracture orientations indicating possible inter-fracture stress effects. At the scale of stage lengths however, we see evidence of inter-stage stress effects. We also observe fracture morphology being strongly driven by rock properties and changes in lithology. Identified proppant distribution along the cored interval is also correlated with roughness variations and we observe strong positive correlation between proppant concentrations and fracture roughness at the local scale. Finally, based on the observed distribution of hydraulic fracture properties, we propose a conceptual spatio-temporal model of fracture propagation which can help explain the hydraulic fracture roughness distribution and ties in other observations as well.


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.


2020 ◽  
Vol 38 (6) ◽  
pp. 2466-2484
Author(s):  
Jianguang Wei ◽  
Saipeng Huang ◽  
Guangwei Hao ◽  
Jiangtao Li ◽  
Xiaofeng Zhou ◽  
...  

Hydraulic fracture initiation and propagation are extremely important on deciding the production capacity and are crucial for oil and gas exploration and development. Based on a self-designed system, multi-perforation cluster-staged fracturing in thick tight sandstone reservoir was simulated in the laboratory. Moreover, the technology of staged fracturing during casing completion was achieved by using a preformed perforated wellbore. Three hydraulic fracturing methods, including single-perforation cluster fracturing, multi-perforation cluster conventional fracturing and multi-perforation cluster staged fracturing, were applied and studied, respectively. The results clearly indicate that the hydraulic fractures resulting from single-perforation cluster fracturing are relatively simple, which is difficult to form fracture network. In contrast, multi-perforation cluster-staged fracturing has more probability to produce complex fractures including major fracture and its branched fractures, especially in heterogeneous samples. Furthermore, the propagation direction of hydraulic fractures tends to change in heterogeneous samples, which is more likely to form a multi-directional hydraulic fracture network. The fracture area is greatly increased when the perforation cluster density increases in multi-perforation cluster conventional fracturing and multi-perforation cluster-staged fracturing. Moreover, higher perforation cluster densities and larger stage numbers are beneficial to hydraulic fracture initiation. The breakdown pressure in homogeneous samples is much higher than that in heterogeneous samples during hydraulic fracturing. In addition, the time of first fracture initiation has the trend that the shorter the initiation time is, the higher the breakdown pressure is. The results of this study provide meaningful suggestions for enhancing the production mechanism of multi-perforation cluster staged fracturing.


Geofluids ◽  
2021 ◽  
Vol 2021 ◽  
pp. 1-24
Author(s):  
Xin Zhang ◽  
Yuqi Zhang ◽  
Bingxiang Huang

Hydraulic fracturing applications have shown a stress disturbance effect during hydraulic fracture propagation, which is often ignored. Using laboratory and discrete element numerical simulation tests, hydraulic fracture propagation under this stress disturbance is systematically studied. The results show that during hydraulic fracturing, the bedding plane is damaged by the stress disturbance, forming a bedding fracture zone (BFZ). The nonlinear fracture characteristics of the formation process of the disturbed fracture zone are revealed, and two indexes (the number of fractures in the disturbed fracture zone and the size of the disturbed fracture zone) are proposed to evaluate the fracturing effect of the stress disturbance. Based on these indexes, multifactor sensitivity tests are conducted under different geological conditions and operational factors. When the principal stress ( σ 1 ) difference is large, the number of hydraulic fractures gradually decreases from many to one, and the direction of the hydraulic fractures gradually approaches the vertical direction of σ 3 , but the change in the in situ stress condition has no obvious effect on the stress disturbance effect. The weaker the bonding strength of the bedding plane, the more significant the stress disturbance effect is, and the easier it is for the fractures to expand along the bedding plane. With increasing injection rate, the stress disturbance effect first increases and then decreases, and the hydraulic fracture propagates from along the bedding plane to cross the bedding plane. With increasing relative distance between the injection hole and bedding plane, the stress disturbance effect presents a linearly increasing trend, and the hydraulic fractures along the bedding planes extend. Based on the experimental results, the relationship between the fracturing effect of the stress disturbance and the extension mode of the hydraulic fracture is determined, and an optimization method for hydraulic fracturing in composite rock reservoirs is given. The research results can provide a theoretical basis for controlling the formation of complex fracture networks in composite rock reservoirs.


Fractals ◽  
2018 ◽  
Vol 26 (02) ◽  
pp. 1840009 ◽  
Author(s):  
KAI ZHANG ◽  
XIAOPENG MA ◽  
YANLAI LI ◽  
HAIYANG WU ◽  
CHENYU CUI ◽  
...  

Hydraulic fracturing is an important measure for the development of tight reservoirs. In order to describe the distribution of hydraulic fractures, micro-seismic diagnostic was introduced into petroleum fields. Micro-seismic events may reveal important information about static characteristics of hydraulic fracturing. However, this method is limited to reflect the distribution area of the hydraulic fractures and fails to provide specific parameters. Therefore, micro-seismic technology is integrated with history matching to predict the hydraulic fracture parameters in this paper. Micro-seismic source location is used to describe the basic shape of hydraulic fractures. After that, secondary modeling is considered to calibrate the parameters information of hydraulic fractures by using DFM (discrete fracture model) and history matching method. In consideration of fractal feature of hydraulic fracture, fractal fracture network model is established to evaluate this method in numerical experiment. The results clearly show the effectiveness of the proposed approach to estimate the parameters of hydraulic fractures.


Sign in / Sign up

Export Citation Format

Share Document