Abstract

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.

1. Introduction

Directional fracturing is a common and significant problem in rock engineering [13]. Directional hydraulic fracturing is an effective directional fracturing technology, which has been widely applied in oil and gas engineering and mining engineering, such as hydraulic fracturing with oriented perforation and preslotting directional hydraulic fracturing [47]. The oriented perforation is a key technology in reservoir stimulation of oil and gas wells, which can control the initiation direction of fractures in near-wellbore [813].

At present, preslotting directional hydraulic fracturing technology includes two types according to the slotting method. One is the mechanical slotting method, and another one is the hydraulic slotting method. The mechanical preslotting directional hydraulic fracturing is achieved by using a mechanical device to produce a wedge ring slot at the bottom of the borehole, which induced the borehole to initially fractured along the wedge ring slot direction due to the stress concentration [1418]. The hydraulic method for preslotting directional hydraulic fracturing is achieved by pre-water jet slotting to produce a fracturing plane along the radial or axial of the borehole. A stress concentration will appear at the tip of the fracturing plane, which caused the borehole to initiate along the desired direction during the subsequent hydraulic fracturing [1921].

However, these current directional hydraulic fracturing technologies only focus on controlling the directional initiation of hydraulic fractures by taking some directional measure to ensure the fracture initiation along the desired direction. The subsequent directional propagation is indirectly controlled through directional initiation. The related research indicated that when the perforation or preslotting is not perpendicular to the direction of the minimum principal in situ stress, the initiation direction of hydraulic fracture will be first along the direction of the perforation or preslotting direction, but the final propagation direction of hydraulic fracture will gradually turn to be perpendicular to the direction of the minimum principal in situ stress [22]. The hydraulic fracture reorientation took place within 2.2–5.8 wellbore diameters. Therefore, the effective range of orientation is very limited [10].

This study focuses on a rock directional fracturing method controlled by hydraulic fracturing in dense linear multiholes (Figure 1) [23, 24]. This method considered the superimposing and coupling of poroelastic stress and pore pressure between the dense linear boreholes to produce a stress concentration area. Both the initiation and propagation of hydraulic fractures were controlled to be along the direction of boreholes layout, which can form a better directional fracturing plane. In the process of directional fracturing of the hard roof in underground coal mines, rapid drilling construction can be carried out with anchor rigs. The spacing of dense boreholes can be 0.5–1 m, and the number of boreholes in a group can be 5–10 to perform the controlling directional hydraulic fracturing. Directional fracturing controlled by hydraulic fracturing in dense linear multi-boreholes is a novel directional fracturing technology of rock mass, which has been applied to the ground control in mines. It also can be applied to the directional fracturing and excavation of rock in underground geotechnical engineering, tunnel engineering, and directional excavation of slope engineering. As long as the rock is hard and has good integrity, this technology can be applied.

In this paper, a physical model experiment was performed to study the fracture propagation process between multiboreholes. The dynamic initiation and propagation process of cracks between multiboreholes during directional hydraulic fracturing was simulated. The evolution of poroelastic stress and pore pressure between multiboreholes was analyzed.

2. Physical Model Experiment

2.1. Experimental Design

A true triaxial hydraulic fracturing experimental system was used in this study (Figure 2), and a cubic cement mortar sample with a size of 500 × 500 × 500 mm3 was used to simulate the coal rock mass (Figure 3). The physical and mechanical properties of the sample are shown in Table 1.

A three-dimensional stress field with  = 6 MPa,  = 5 MPa, and  = 4 MPa was loaded to the sample to simulate the in situ stress. During the experiment, one channel was used to inject high-pressure water into the three boreholes with a pumping rate of 600 mL/min. The pumping pressure is equal in each borehole.

2.2. Experiment Results and Analysis

The water pressure curve during hydraulic fracturing is shown in Figure 4, which can be divided into five parts: the rise of water-injection pressure (AB), cracks initiation (BC), cracks independent propagation (CD), cracks intersecting and overall propagation (DE), and pressure release due to the hydraulic fracture propagating to the sample surface (EF), respectively.

In the first 100 s, the water output from the pump is used to fill the pipeline and boreholes, after which the water pressure rises rapidly (AB). At 173 s, the water pressure rose to the maximum value of 9.62 MPa and crack initiation occurs. After crack initiation, the effective volume used for filling the initial hydraulic cracks decreases, which leads to the water pressure decreased slightly by 21.04% (BC). After that, the hydraulic fracture propagates independently, and the water pressure is stable at around 8.3 MPa (CD). At 370 s, the hydraulic fractures of the three boreholes intersect and connect with each other, and the water pressure curve decreased sharply, with a drop of 38.43%. Subsequently, the hydraulic fracture propagates overall. The water pressure rises again to the second peak and then falls again. It continuously fluctuates within the range of 4–6 MPa (DE). The pumping stopped when the hydraulic fracture propagates through the surface of the sample. The water pressure drops rapidly and then gradually decreases ‵(EF).

After the experiment, the sample was taken out from the experimental frame. The morphology of the hydraulic fracture was shown in Figure 5. A directional fracturing plane was formed along the direction of the boreholes layout, and the surface of the fracturing plane is relatively flat. All three boreholes were fractured, the cracks propagated towards each other between holes finally intersected and connected.

The poroelastic stress concentration zone and pore pressure increase zone occurred between the holes after injecting water simultaneously into the multiholes. It created a maximal effective tensile stress at the direction of boreholes layout, which priority to exceed the fracturing strength of the sample. Therefore, the fractures initiated along the direction of boreholes layout due to the superimposing of poroelastic stress and pore pressure, and hydraulic-mechanics coupling effect between boreholes. Because the fractures propagation area was still located at the superimposing area of poroelastic stress and pore pressure, the fractures did not alter the extending direction. They continue to propagate along the initiation direction and extend toward each other between boreholes. A wide range of ruptures occurred when the cracks tip intersects and connect each other, which caused a sharp water pressure relief. After that, a connected water pressure area was formed. The hydraulic fractures continued to extend as an overall fracturing plane until it extended to the surface of the sample. The initiation and propagation direction of the hydraulic fractures in each hole is consistent and all along the direction of boreholes layout, which is the guarantee for forming a single flat directional fracturing plane.

3. Numerical Simulation Model Using RFPA2D-Flow

3.1. Background of RFPA2D-Flow

The Coupled Analysis of Flow and Solid Mechanics in Rock Failure Process Analysis software RFPA2D-Flow was developed by Dalian Mechanics Software Co., Ltd. (China) [2527]. This is a numerical test software that can simulate the progressive failure of a quasi-brittle material. The calculation method is based on finite element theory and statistical damage theory. This method takes into account the non-uniformity of material properties and the random distribution of defects [25].

The formulation of seepage-stress coupling is based on the following basic assumptions, which are described in detail elsewhere [26, 28, 29]:(1)The seepage process of rock mass satisfies Biot’s consolidation theory and the modified Terzaghi effective stress principle.(2)The rock mass is fully saturated and is a brittle elastic material with residual strength. Besides, its loading and unloading behaviours are following the elastic damage mechanics.(3)Regard the maximum tensile strength criterion and the Mohr–Coulomb failure envelope as the damage threshold condition of element.(4)Variation in the permeability of rock as a function of the stress states in elastic deformation. The variation increases dramatically when the element fails.(5)The local heterogeneity in the properties of rock masses is defined by the Weibull function, which is given as follows:where is the macroscopic magnitude of the property, is the scale parameter, and is the homogeneity.

The seepage-stress coupling process in a saturated rock mass can be explained by Biot’s consolidation theory [29], RFPA2D-Flow extended Biot’s consolidation theory by considering the influence of stress on permeability, through which its basic governing equations were formed:

3.1.1. Balance Equation

where is the stress and is the body force in the direction.

3.1.2. Geometrical Equation

where is the strain and is the displacement in the direction.

3.1.3. Constitutive Equation

where is the coefficient of pore-fluid pressure, is the pore pressure, is the Kronecker constant, is the Lame coefficient, and is the shear modulus.

3.1.4. Seepage Equation

where is the coefficient of permeability and is Biot’s constant.

3.1.5. Coupling Equation

where is the initial coefficient of permeability and is material constant.

Equations (2)∼(5) are based on Biot’s theory of consolidation [30]. Equation (6) represents the influence of stress on permeability, which is achieved by assuming that both permeability and stress follow a negative exponential function.

3.2. The Design of Numerical Model

The numerical model is a cubic rock sample with a cross section of 1000 mm × 1000 mm. A two-dimensional model was established for the sample section, and the model was divided into 500 × 500 units by using the plane strain method. The program of numerical simulation is shown in Table 2 and Figure 6. The modified Mohr–Coulomb strength criterion is taken as the failure criterion of rock, where the tension fractures of rocks are taken into consideration. The parameter selection is important for the numerical model. Many studies have been carried out by using a neural network model to predict the mechanical parameters of rock, such as unconfined compressive strength, Young’s modulus, and shear strength [3137]. The spacing and diameter of the boreholes are designed according to the field construction parameters. Five boreholes with a spacing of 141 mm and a diameter of 18 mm are drilled to perform the hydraulic fracturing synchronously, and each fractured borehole is equal in water pressure. The angle between the boreholes layout direction and the maximum principal stress direction is 15°. The water pressure of boreholes is loaded step by step from 0 MPa at an increment of 0.3 MPa, which simulates the process of hydraulic fracturing.

4. Numerical Results and Analysis

4.1. Evolution Process of Pore Pressure

Pore pressure plays an important role during hydraulic fracturing [38]. The contour variation of pore pressure in the process of directional hydraulic fracturing controlled by dense linear multiboreholes is shown in Figure 7. As water pressure grows, the pressurized water infiltrates into the surrounding rocks of boreholes, which produce pore pressure and its gradient. The distribution of pore pressure builds an oval region of infiltration whose macro axis is the ligature of these boreholes. The axial ratio becomes smaller as the distance from the axes of the boreholes increases. It means that the infiltration coefficients along the direction of boreholes layout are relatively small, which leads to the heavy attenuation of water pressure and large gradient of pore pressure. On the contrary, the infiltration coefficients whose direction is perpendicular to the direction of boreholes layout experience the opposite variation compared with the case above, which leads to slight attenuation of water pressure and a small gradient of pore pressure (Figure 6, Step 2-1). The existence of pore pressure contributes to the prior crack initiation of the borehole walls along the direction of the boreholes layout. In the central area of the boreholes layout, there is an oval region with high pore pressure around every hole. Due to the superimposing effect of pore pressure between dense linear multiboreholes, the long axis of the region is along the direction of boreholes layout. Meanwhile, every region stands alone and does not connect (Figure 6, Step 2-1). When the water pressure goes up to 9.9 MPa (Figure 7, Step 34-1), cracks initiation begins. In this case, the infiltration coefficient along the boreholes layout grows. The high-pressure regions of three middle-positioned boreholes are connected (Figure 8, Step 34-1). When the water pressure continues to rise at 10.8 MPa, hydraulic cracks further propagate, and the highest pressure regions are connected (Figure 7, Step 37-2), which results in a connected region of peak pore pressure. As the water pressure reaches 11.2 MPa, the hydraulic cracks within boreholes gradually propagate closer and tend to propagate in the direction of the boreholes layout (Figure 7, Step 38-3). Finally, a connected fracturing plane formed along the direction of the boreholes layout (Figure 7, Step 38-6).

4.2. Evolution Process of Poroelastic Stress

Due to the fluid-structure coupling, the pore pressure changes will also induce poroelastic stress change [39]. The contour of horizontal poroelastic stress in the process of hydraulic fracturing is shown in Figure 9. As the pressurized water injects the multiboreholes simultaneously, poroelastic stress concentration regions occur around every borehole. Some regions parallel to the direction of boreholes layout, while the others perpendicular to it. The poroelastic stress concentration regions are connected between boreholes. The poroelastic stress around the middle boreholes is larger than that of the outside two boreholes. Besides, the poroelastic stress grows with the decrease of its distance from boreholes (Figure 10). After the crack initiation along with the direction of boreholes layout, the stress drops from both sides of fractures, which presents an elliptical contour of poroelastic stress. The drop speed of poroelastic stress goes up as the distance from the boreholes layout decreases. With the cracks propagating towards each other, the distance between cracks decreases. The augment of poroelastic stress superposition leads to the poroelastic stress at the tips of cracks increase. The superposing of poroelastic stress fields at the tips of cracks contribute to the connection of fractures during the propagation towards each other, which finally results in a directional fracturing plane along the direction of boreholes layout.

4.3. The Role of Controlling Directional Fracturing

Directional hydraulic fracturing controlled by dense linear multi-boreholes is achieved by controlling the spacing and size of the boreholes to make stresses superimposed between boreholes. The stress superposing causes the stress concentration zone to appear in the direction of the boreholes layout, where the broken occurs preferentially. The role of controlling directional fracturing includes the superimposing of poroelastic stress and pore pressure, and the hydraulic-mechanics coupling effect.

4.3.1. Superimposing of Poroelastic Stress

The rock mass is in equilibrium with stress before hydraulic fracturing, borehole drilling can change the stress distribution inside the rock mass, and an area of stress concentration will occur around the borehole. Besides, the borehole will be squeezed by high-pressure water when hydraulic fracturing, which increases the area of stress concentration around the borehole.

The boreholes spacing of dense linear multiboreholes hydraulic fracturing must be smaller than the stress concentration range around boreholes, which keep the adjacent borehole within the stress concentration zone of the borehole. The poroelastic stress concentration zone of adjacent holes will affect each other and superimposed. Finally, a connected poroelastic stress concentration zone between the boreholes is formed. Poroelastic stress gradient will occur between the poroelastic stress concentration zone and the in situ stress zone, which will control the fracture initiation and propagation. The poroelastic stress gradient direction (tangential tensile stress direction) is perpendicular to the direction of the boreholes layout. Therefore, the initiation and propagation direction of hydraulic fractures will be along the direction of the boreholes layout.

4.3.2. Superimposing of Pore Pressure

Rock is a porous medium. Pressurized water can penetrate the rock around the borehole to form pore water pressure during hydraulic fracturing. An area of pore pressure increase will also form around the borehole.

The layout of dense linear multiboreholes keeps each borehole within the pore pressure increase zone of its adjacent borehole. The pore pressure increase zone of adjacent boreholes will affect each other and superimposed. Finally, a connected pore pressure increase zone between the boreholes is formed. Pore pressure gradient will occur between the pore pressure concentration zone and the initial pore pressure zone. The pore pressure gradient direction is also perpendicular to the direction of boreholes layout, which will also induce the initiation and propagation direction of hydraulic fractures along the direction of boreholes layout.

4.3.3. Coupling Effect of Hydraulic-Mechanics

The coupling effect of hydraulic-mechanics between multiboreholes also plays an important role in directional hydraulic fracturing controlled by dense linear multiboreholes. The superposing of poroelastic stress causes the pore structure changes and affects the distribution of pore pressure. Meanwhile, the superposition of pore pressure also causes pore structure changes and affects the distribution of poroelastic stress.

5. Conclusions

This study presents the dynamic directional initiation and propagation process of hydraulic fractures between dense linear multi-boreholes through numerical simulation and physical experiment. The dynamic evolution of pore pressure and poroelastic stress during fracturing was also analyzed. The primary conclusions are as follows:(1)The water pressure curve during the physical model experiment reflected the fractures propagation, intersection, and connection process between boreholes. A directional fracturing plane was formed along the direction of the boreholes layout, and the surface of the fracturing plane is relatively flat.(2)The initiation and propagation direction of the hydraulic fractures in each borehole are consistent and all along the direction of boreholes layout, which is the guarantee for forming a single flat directional fracturing plane.(3)Directional fracturing of dense linear multiholes hydraulic fracturing is achieved by controlling the spacing and size of the boreholes to make stresses between boreholes superimposed. The initiation and propagation direction of hydraulic fractures between boreholes can be controlled along the direction of boreholes layout by the superimposing of poroelastic stress and pore pressure, and hydraulic-mechanics coupling effect.(4)The poroelastic stress concentration zone and pore pressure increase zone appears between holes in the direction of the boreholes layout during directional hydraulic fracturing. 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 boreholes layout, the poroelastic stress on both sides of the direction of boreholes layout decreases.(5)The future work will focus on the physics model experiment to analyze the influence of principal stress difference, injection rate, borehole aperture, and spacing using real rock samples.

Data Availability

The data used in the study will be available on request.

Conflicts of Interest

The authors declare no conflicts of interest.

Acknowledgments

The authors express their gratitude to the National Natural Science Foundation of China. This research was funded by the National Natural Science Foundation of China (nos. 51774272 and 52004269).