A Numerical Investigation on the Hydraulic Fracturing Effect of Water Inrush during Tunnel Excavation

Zhang, Quan; Huang, Bingxiang; He, Manchao; Guo, Shan

doi:https://doi.org/10.1155/2020/6196327

Geofluids

On this page

Abstract Introduction Methods Results Discussion Conclusions Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2020 | Article ID 6196327 | https://doi.org/10.1155/2020/6196327

A Numerical Investigation on the Hydraulic Fracturing Effect of Water Inrush during Tunnel Excavation

Quan Zhang,^1,2,3Bingxiang Huang ,¹Manchao He,³and Shan Guo³

Academic Editor: Umberta Tinivella

Received27 Jan 2020

Revised25 Oct 2020

Accepted09 Nov 2020

Published09 Dec 2020

Abstract

When a high-pressure water source is located near a tunnel under excavation, water inrush is commonly associated with a hydraulic fracturing effect. To study the hydraulic fracturing effect of water inrush (HFEWI), flow-rock failure process analysis (F-RFPA^2D) was adopted to simulate the water inrush process. The simulated results indicated that a stress disturbance area formed in front of the excavation face and that a hydraulic fracture zone formed in front of the karst cavity. Similarly, stress concentrations formed in front of the excavation face and the karst cavity. The hydraulic fracturing effect was characterized by stress concentration, and the local hydraulic crack propagation was the result of stress concentration. In addition, a pore pressure gradient formed in the crack-free area of the surrounding rock, and the occurrence of hydraulic cracking was the root cause of the significant change in water flow. When the hydraulic cracks initially formed and expanded, the zone of crack activity was large. As the cracks continued to expand, the range of activity decreased and finally concentrated directly in front of the excavation face. Additionally, the shapes of the water inrush channel obtained by the experimentation and numerical simulation were basically the same: semielliptical. During the evolution of hydraulic crack initiation, expansion, and penetration, the bottom of the excavated borehole was initially dry and then experienced seepage and water inrush. Finally, the minimum safe thickness of the rock wall was calculated to provide a safety guideline for this type of water inrush.

1. Introduction

Water inrush disasters occur frequently during the mining of mineral resources and the excavation of traffic tunnels and hydropower stations [1–4]. With the reduction in shallow mineral resources [5], mining depths have been increasing year over year [6]; thus, the water pressures expected at these mining depths have gradually increased as well. When the mining depth is more than 1000 m, karst water pressure can exceed 7 MPa [7]. For example, the water pressure in the Xinyi Mine, Yi’an Mine, and Mengjin Mine of the Yimei Group is 7.5 MPa, and the water pressure of the Ordovician limestone confined water in the Shandong Xingcun Minefield can reach 13 MPa [8]. The water pressure in deep tunnels is also high. For example, the maximum burial depth of the diversion tunnel of the Jinping II Hydropower Station is approximately 2525 m, and the maximum water pressure exceeds 10 MPa [9]. Therefore, when a rock mass excavation project becomes deep, the pressure of the pressurized water is high and hydraulic fracturing is prone to occur during the water inrush process [10].

Hydraulic fracturing is one of the core mechanisms of the water inrush process under high water pressure [11]. Water can continuously induce hydraulic fracturing at crack tips [12], so the rock failure mechanism [13] of the water inrush channel caused by hydraulic fracturing is different from that caused by solid stress [14]. The mechanical action of groundwater is induced by its hydrostatic pressure and hydrodynamic pressure, both of which can cause the hydraulic fracturing effect of water inrush (HFEWI) [15]. There are two HFEWI types: pure splitting and stress splitting. The pure splitting type reflects hydraulic fracturing of a fractured rock mass, while the stress splitting reflects the hydraulic fracturing of an intact rock mass [16]. The HFEWI is associated with a high risk of fracturing, rapid flow rate, and high water flow pressure [17].

At present, the HFEWI has mainly been analyzed theoretically. From the perspective of fracture mechanics, the critical water pressure of the HFEWI from an unidentified fault in a coal floor was calculated [18]. The HFEWI occurred in the sidewall of a collapse column when the collapse column was far from the mining face or roadway, and the intrinsic cracks in the coal and rock mass were opened by the groundwater via hydraulic fracturing, which caused water inrush [19]. Additionally, the critical burst pressure of water inrush was derived for this case [19, 20]. For the water inrush from a fractured rock mass, the formulas for calculating the critical water pressure of the HFEWI from a deep tunnel under the two failure modes of tensile-shear composite fracturing and compression-shear composite fracturing were deduced [21]. In addition, the critical water pressure for compression-shear failure was much less than that for tension-shear failure [19, 22], so the rock mass was more likely to undergo compression-shear failure. Moreover, crack propagation under the HFEWI presented intermittent characteristics [22]. To prevent water inrush in the case in which a high-pressure water source is in front of the excavation face, the minimum safe thickness of the rock wall was calculated [23, 24]. However, most of these theoretical analyses focused on fractured rock masses. In addition, few numerical simulations have been carried out for the HFEWI. A functional relationship between the permeability coefficient and the damage variable was established and substituted into a numerical model to analyze the HFEWI from a mine floor [14]. The finite difference continuum code FLAC3D was used to analyze the HFEWI of a bottom collapse column. When the water pressure in a collapse column was greater than the minimum principal stress of the key layer of the bottom plate, the hydraulic fracturing phenomenon occurred in the water-blocking rock mass and formed a pipe surge, which led to mine water inrush [25]. By introducing a seepage-stress coupling model and employing ABAQUS to simulate the water inrush in deeply buried tunnels during excavation, it was considered that the high hydraulic gradient triggered the HFEWI [26]. Regarding the numerical simulation of the HFEWI, the floor water inrush during mining and the water inrush in a tunnel with high-pressure water have been studied, but research on water inrush from karst cavities with high water pressures during tunnel or roadway excavation has rarely been carried out. The water pressure of the current water inrush experiments is still relatively low and has not reached the level of hydraulic fracturing. Therefore, there has been no experimental study of the HFEWI.

This is the first study to investigate the HFEWI when a karst cavity with a high water pressure is encountered in front of the face of a tunnel (Figure 1) with a true triaxial water inrush experiment. However, the experimental setup could not support observations of the whole process of crack initiation, expansion, and penetration through the test block, and only the final shape of the water inrush channel was inspected. The numerical simulation considered can visually reproduce the dynamic process of crack evolution under the disturbance action of the excavation and hydraulic splitting action of the high-pressure water. In addition, flow-rock failure process analysis (F-RFPA2D), which is based on the flow-stress-damage (FSD) coupling model, can describe the evolution of the stress field (maximum principal stress, minimum principal stress, and shear stress) and the seepage field (pore water pressure and volume of flow) during the process of water inrush channel formation. These models compensate for the deficiencies of the physical experiment. Therefore, F-RFPA2D was used to simulate the HFEWI of an intact rock mass. This work is conducive to the prevention and control of water inrush disasters, thereby increasing the safety of tunnel excavation.

2. Methods

2.1. Introduction of the Numerical Model

Considering the nonuniformity of a rock mass, F-RFPA2D transforms a complex macroscopic nonlinear problem into a simple mesoscopic linear problem that can accurately simulate the FSD coupling mechanism of rock during crack initiation and expansion. This software has been widely used in fluid-solid coupling studies [28], and its established numerical models have been verified many times [29–33].

Seepage-stress coupling follows classical Biot consolidation theory [34]. The main governing equations are as follows.

Equilibrium equation is as follows:

Strain-displacement equation is as follows:

Constitutive equation is as follows:

Seepage equation is as follows:where is the sum of the normal stress in the - plane, is the body force in the th direction, is the rock density, is the body strain ( is a positive strain), is the Kronecker delta, is the shear modulus, is a Lamé constant, is the Biot constant, and is the Laplacian operator.

The mesostructural mechanical parameters of the material obey a Weibull distribution [35], which represents nonuniformity, that is,where is a certain parameter of the mesoscopic unit, is the statistical average of the parameter of the mesoscopic unit, and is the heterogeneous coefficient. The strength distribution changes with five different heterogeneous coefficients, as shown in Figure 2. The elastic modulus follows the same distribution.

2.2. Numerical Simulation Scheme

The numerical simulation considered a plane strain model. As shown in Figure 3, the model size was , which was divided into . An karst cavity with a water pressure of 3 MPa was represented in the model. The sample parameters [26] are shown in Table 1. The uniaxial compressive strength and uniaxial tensile strength of the rock were 20 MPa and 2 MPa, respectively. The residual strength factor was 0.1. The elastic modulus was 10000 MPa. The parameters such as elastic modulus and rock mass strength were randomly distributed according to a Weibull distribution and expressed by the homogeneity coefficient . The greater was, the more uniform the parameter. The homogeneity coefficient was 4. The permeability coefficient of the rock was 0.01 m/d, and the permeability coefficient jump rate was 20. In addition, the pore water pressure coefficient was 0.1. Poisson’s ratio and internal friction angle were 0.25 and 30°, respectively. The vertical stress was set to 5 MPa, and the horizontal stress was 1.5 MPa. The water inrush process was triggered by an excavation disturbance. The diameter of the excavation was 20 mm, and the number of excavation steps was set to 3. In addition, the excavation distances were 30 mm, 25 mm, and 20 mm.

3. Results

3.1. The Evolution Pattern of the Stress Field

Excavation caused stress redistribution in the surrounding rock, and a stress disturbance zone was formed in front of the excavation face. A hydraulic fracture zone was created near the karst cavity. Stress concentrations formed in front of the excavation face and near the karst cavity. The drilling centerline (DCL) was selected to analyze the stress change process (Figure 4(a)). When the excavation face was 80 mm from the karst cavity, the maximum shear stress reached a peak of 1.84 MPa (Figure 4(a)) at 4 mm in front of the excavation face, the maximum value of the maximum principal stress was 4.48 MPa at 7 mm in front of the excavation (Figure 5(a)), and the maximum value of the minimum principal stress was 3.10 MPa at 23 mm in front of the excavation (Figure 6(a)). Furthermore, the shear stress and the minimum principal stress near the karst cavity also exhibited the phenomenon of stress concentration. When the distance between the excavation face and the karst cavity was 55 mm, the stress disturbance zone in front of the excavation face became larger, and the peak values of the shear stress, maximum principal stress, and minimum principal stress along the DCL were 4.19 MPa (Figure 4(b)), 9.08 MPa (Figure 5(b)), and 4.43 MPa (Figure 6(b)), respectively, which were greater than the corresponding results when the excavation face was 80 mm from the karst cavity. Moreover, as the distance between the excavation face and karst cavity decreased, the stress gradient in the surrounding rock increased. The peak shear stress, maximum principal stress, and minimum principal stress near the karst cavity were 2.52 MPa, 4.86 MPa, and -1.28 MPa, respectively, which were lower than the peak values of the stresses near the excavation face. When the excavation face was 35 mm from the karst cavity, the peak shear stress, maximum principal stress, and minimum principal stress along the DCL were 6.84 MPa (Figure 4(e)), 17.62 MPa (Figure 5(e)), and 3.95 MPa (Figure 6(e)), respectively. Additionally, the peak shear stress and maximum principal stress were higher than those and the peak minimum principal stress was slightly lower than that when the distance between the excavation face and the karst cavity was 55 mm. However, the distance between the excavation face and the karst cavity was further decreased, so the stress gradient in the surrounding rock still increased.

(a) Crack initiation stage

(b) Crack propagation stage (Step 2-2)

(c) Crack propagation stage (Step 2-11)

(d) Crack propagation stage (Step 2-21)

(e) Crack penetration stage

(a) Crack initiation stage

(b) Crack propagation stage (Step 2-2)

(c) Crack propagation stage (Step 2-11)

(d) Crack propagation stage (Step 2-21)

(e) Crack penetration stage

(a) Crack initiation stage

(b) Crack propagation stage (Step 2-2)

(c) Crack propagation stage (Step 2-11)

(d) Crack propagation stage (Step 2-21)

(e) Crack penetration stage

There was no hydraulic crack initiation in the surrounding rock during the initial excavation, but energy was concentrated at the stress concentration. As the tunnel or roadway continued to be excavated, the rock at the stress concentration fractured, and hydraulic cracks were generated there. As the hydraulic cracks propagated forward, the position of the stress peak near the hydraulic fracture zone continuously moved forward (Figures 4(b)–4(d), Figures 5(b)–5(d), and Figures 6(b)–6(d)). This indicated that the stress peak near the karst cavity was located at the tip of the hydraulic crack. After the hydraulic cracks were generated, the water flowed into these cracks and created obvious stress concentrations and accumulated energy in the rock around the crack tips. When the energy at a crack tip exceeded the threshold required for hydraulic crack propagation, the crack expanded, releasing the energy and reducing the stress. Then, the water flowed into the extended crack tip and formed another stress concentration zone. Therefore, the hydraulic fracturing effect during the crack propagation stage is characterized by stress concentration, and local crack propagation is the result of this stress concentration. Under the splitting action of water pressure, the entire rock mass does not reach the threshold of damage simultaneously; instead, the phenomenon of stress concentration appears locally. When the local concentrated energy exceeds the energy required for rock failure, hydraulic cracks propagate.

The similarity of the three types of stress evolution was that the three types of stress all formed stress concentrations in front of the excavation face and karst cavity, causing the surrounding rock to rupture and the water inrush channel to form. There were three differences between the three types of stress. First, the peak stress values of the shear stress and the minimum principal stress were smaller than the peak stress value of the maximum principal stress at the same excavation step. Second, the peak stress values of the shear stress and the maximum principal stress in front of the excavation face were all larger than those in front of the karst cavity during the whole excavation process. However, the peak stress value of the minimum principal stress in front of the excavation surface was greater than that in front of the karst cavity when the tunnel excavated to Step 1-1, Step 2-2, and Step 2-11. As the tunnel continued to be excavated to Step 2-21, the peak stress value of the minimum principal stress in front of the excavation face was less than that in front of the karst cavity. Third, the sign of stress represented the direction of stress in Figures 4–6. The shear stress and the maximum principal stress at different positions were always positive throughout the water inrush process, indicating that the direction of the stress did not change. However, the minimum principal stress changed from a positive value to a negative value at some locations, and from a negative value to a positive value at other locations, so the direction of the minimum principal stress direction changed.

3.2. The Evolution Pattern of the Seepage Field

The change in pore water pressure (Figure 7) showed that the pore water pressure was relatively high at the crack but low in the crack-free portion along the DCL. This was because the pressurized water flowed into the crack. When the distance from the excavation face to the karst cavity was 80 mm (Figure 7(a)), no crack was generated along the DCL, and the pore pressure reduced almost linearly from the edge of the karst cavity (0.3 MPa) to the face of the excavation (0 MPa), so a pore water pressure gradient formed in the surrounding rock. When the excavation face was 55 mm (Figures 7(c) and 7(d)) and 35 mm (Figure 7(e)) from the karst cavity, the pore water pressure at the crack-free area gradually decreased from the karst cavity to the excavation face but exhibited a peak pressure at the crack.

(a) Crack initiation stage

(b) Crack propagation stage (Step 2-2)

(c) Crack propagation stage (Step 2-11)

(d) Crack propagation stage (Step 2-21)

(e) Crack penetration stage

The flow rate (Figure 8) reflects two patterns. One is that the horizontal flow was much greater than the vertical flow, indicating that the flow was mainly parallel to the direction of the DCL. Figure 8(a) shows that the flow in the horizontal direction of the surrounding rock along the DCL ranged from 322 to 449 m³/(d·m), while that in the vertical direction ranged from 0.01 to 5.48 m³/(d·m). Figures 8(b)–8(e) show the same trends. The other rule is that the flow along the DCL was low before hydraulic cracks were generated, while the flow sharply increased after the cracks occurred. When there was no crack in the surrounding rock, the horizontal and vertical flow rates were low (Figure 8(a)). With the initiation (Figure 8(b)), expansion (Figures 8(c) and 8(d)), and penetration (Figure 8(e)) of the hydraulic cracks, the horizontal and vertical flow rates dramatically increased. After formation of the water inrush channel, the horizontal flow along the DCL reached 53050 m³/(d·m) (Figure 8(e)), which was 118 times the maximum flow rate (449 m³/(d·m)) in the horizontal direction when no crack occurred (Figure 8(a)). Additionally, the maximum vertical flow rate reached 17908 m³/(d·m), which was much higher than the maximum vertical flow rate (5.48 m³/(d·m)) when no crack was generated (Figure 8(a)).

(a) Crack initiation stage

(b) Crack propagation stage (Step 2-2)

(c) Crack propagation stage (Step 2-11)

(d) Crack propagation stage (Step 2-21)

(e) Crack penetration stage

3.3. The Evolution Pattern of the Fracture Field

Under the joint effect of high-pressure water and excavation unloading, crack initiation, propagation, and penetration occurred in the surrounding rock, thus gradually forming a crack field. Acoustic emission analysis can record the process of damage and energy release inside a rock mass [36]. Therefore, it was used to describe the evolution process of the fracture field. The pattern of acoustic emissions under different excavation steps is shown in Figure 9. The sound waves in the acoustic emission diagram are represented by points and circles, and the positions of the points and circles illustrate the locations where the rock mass was destroyed.

(a) Crack initiation stage

(b) Crack expansion stage

(c) Crack penetration stage

In the crack initiation stage (Figure 9(a)), microcracks gradually developed near the karst cavity. When the microcracks were generated [37], the rupture energy was released in the form of elastic energy, and an acoustic emission signal appeared near the karst cavity. The crack originated mainly near the karst cavity directly in front of the excavation face. However, there was no crack initiation near the karst cavity far ahead of the excavation face.

During the crack propagation phase (Figure 9(b)), the surrounding rock experienced multiple progressive failures in a total of 26 small steps. When a crack initiated, the water pressure at the crack tip was 0. After the water reached the crack tip, the water pressure increased continuously. When the energy at the tip of the crack reached that required for crack propagation, the energy accumulation ended, and the crack began to expand. Then, the crack evolution entered the period of crack propagation. This process was called the energy accumulation period. Before water inrush occurred, the crack propagation required passing through multiple energy accumulation periods, and the crack propagation range needed to exceed the minimum safe thickness of the rock wall, so the HFEWI exhibited hysteresis. Therefore, this characteristic can be used to ensure the safety of equipment and workers.

In the crack penetration stage (Figure 9(c)), the cracks near the karst cavity directly in front of the excavation face extended to the excavation face. The cracks near the karst cavity far ahead of the excavation face generally did not extend far from the karst cavity; their extents were short and did not reach the excavation face. Since the water pressure is the same everywhere in the karst cavity, the stress disturbance directly in front of the excavation face between the excavation face and the karst cavity is greatest. The acoustic emission signal significantly increased during the crack penetration phase due to the sudden release of energy that had accumulated in the rock during an earlier stage. In addition, the acoustic emission signal was mainly concentrated in the water inrush channel and was more sparsely distributed in other locations.

4. Discussion

4.1. Comparison of Simulation and Experimental Results

A water inrush experiment under true triaxial stress was carried out. The size of the cement mortar sample was (Figure 10). The sample contained a karst cavity with high-pressure water. The karst cavity was cylindrical, with a disc radius of 120 mm and a thickness of 18 mm. The upper and lower circular plates were 4 mm thick. The edges between the upper and lower circular plates were supported by hard plastic and sealed with glue. Therefore, the aquifer was 10 mm high. A water injection pipe was connected behind the aquifer. The stress in the excavation direction of the test block was 6 MPa, and the stresses in the other two directions were 4 MPa and 2 MPa. The test block was placed in the true triaxial water inrush experimental system and excavated from the top down (Figure 11). To observe the shape of the water inrush channel, red dye with a certain particle size was added to the water. There were five excavations in total. The first three excavation depths were 10 cm, 4.5 cm, and 5 cm, and the water pressure during these three excavations was 6.5 MPa. Subsequently, the water pressure was increased to 7 MPa, and the fourth excavation was 2 cm. Next, the water pressure was increased to 8 MPa, but no excavation occurred. Finally, the water pressure was increased to 8.5 MPa, and the fifth excavation was 1 cm.

(a) Schematic diagram [38]

(b) Physical diagram

When the water pressure was 6.5 MPa, the bottom of the hole was initially dry (Figure 12(a)) and then seeped (Figure 12(b)). However, the water did not turn red, since water could penetrate to the excavation surface, but the red pigment particles were large and could not penetrate there. This result indicated that the surrounding rock did not generate hydraulic cracks, only seepage. When the water pressure was 7 MPa, the water turned red (Figure 12(c)), indicating that a hydraulic fracturing effect had occurred and that a water inrush channel had formed. To analyze the influence of the increase of water pressure on the hydraulic fracturing effect in the water inrush process, the water pressure was increased to 8 MPa and 8.5 MPa. After the water pressure rose, the flow of water inrush became larger, and the amount of water inrush increased significantly.

(a)

(b)

(c)

After the experiment, the test block was opened. The red area in Figure 13 shows the shape of the water inrush channel. Numerous cracks were observed near the karst cavity, while fewer cracks were observed near the excavation face; those near the excavation face were concentrated directly in front of the excavation face. The simulation results (Figure 9) also show that when hydraulic cracks initiated and expanded, the hydraulic crack zone near the karst cavity was large. However, when the hydraulic cracks extended to the excavation face, the crack zone decreased, mainly located directly in front of the excavation borehole. Therefore, the morphologies of the water inrush channel obtained by the experiment and that obtained by the simulation are basically the same, and the shapes are all semielliptical.

Under FSD coupling, the hydraulic crack initiation and expansion of water inrush is the same as that of hydraulic fracturing. Therefore, under the action of high-pressure water, the water inrush during excavation has a hydraulic fracturing effect.

4.2. Minimum Safe Thickness of the Rock Wall

According to the simulation and test results, the minimum safe thickness of the rock wall is the sum of two distances (Figure 14): the maximum length () of the stress disturbance zone in front of the excavation face and the maximum length () of the hydraulic fracture zone in front of the karst cavity, that is,

According to the classical hydraulic fracturing crack extension model, the Khristianovich-Geertsma-de Klerk (KGD) model is as follows [39]:where is the water pressure, is the minimum in situ stress, is the shear modulus, is the displacement of water, is the fluid viscosity, is Poisson’s ratio of the rock, and is the crack length. Therefore, the maximum length of the hydraulic fracture zone is as follows:

The maximum length of the stress disturbance zone in front of the excavation face is the thickness of the loose ring [40], that is,where is the radius of the roadway or tunnel, is the internal friction angle of the surrounding rock, is the peak of the supporting pressure, is the slope of the stress-strain curve of the plastic zone, is the adhesion force of the surrounding rock, is the residual cohesive force of the surrounding rock, is the reaction force of the roadway support, is the distance from the center of the circular roadway, is the radius of the plastic zone, , is the elastic modulus of the rock, and .

Therefore, when water inrush exhibits a hydraulic fracturing effect, the minimum safe thickness of the rock wall is as follows:

5. Conclusions

It was previously believed that the mechanism of water inrush was solid stress and water seepage. A new mechanism of water inrush, which is the hydraulic fracturing effect, was proposed. A simulation method based on the FSD coupling model was used to study this mechanism. In addition, the true triaxial water inrush experiment method was adopted for the first time to investigate the HFEWI. Therefore, the following conclusions can be obtained.

The stress concentration caused by excavation and high water pressure makes the cracks in the surrounding rock form, expand, and penetrate, thus creating a water inrush channel. Moreover, this type of water inrush first undergoes three stages: seepage, hydraulic fracturing, and large-scale water inrush. Interestingly, the HFEWI exhibits hysteresis since the crack propagation process requires multiple energy accumulation periods. In addition, the shape of the water inrush channel is semielliptical. Finally, the minimum safety of the rock wall for fracturing water inrush is given quantitatively, which provides a theoretical basis for tunnel safety.

During the experiment, the pattern of the stress field and seepage field could not be studied. Therefore, pore water pressure sensors and skeleton stress sensors will be buried in the test blocks to investigate the patterns of solid stress and water pressure during water inrush. In addition, the HFEWI under high water pressure will need to be verified on site. This study focuses on the HFEWI in intact rock masses. Thus, experimental and simulation studies will be carried out on the HFEWI in fractured rock masses.

Data Availability

Data can be obtained by the email zhangquan@cumt.edu.cn.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors gratefully acknowledge the financial support from the National Key Research and Development Program of China (No. 2016YFC0600901), and the National Natural Science Foundation of China (No. 41941018).

References

D. Ma, X. Cai, Z. Zhou, and X. Li, “Experimental investigation on hydraulic properties of granular sandstone and mudstone mixtures,” Geofluids, vol. 2018, Article ID 9216578, 13 pages, 2018.
View at: Publisher Site | Google Scholar
J. L. Shao, F. Zhou, and W. Sun, “Evolution model of seepage characteristics in the process of water inrush in faults,” Geofluids, vol. 2019, Article ID 4926768, 14 pages, 2019.
View at: Publisher Site | Google Scholar
J. Wang, Y. Zhang, Z. Qin, S. G. Song, and P. Lin, “Analysis method of water inrush for tunnels with damaged water-resisting rock mass based on finite element method-smooth particle hydrodynamics coupling,” Computers and Geotechnics, vol. 126, article 103725, 2020.
View at: Publisher Site | Google Scholar
D. Ma, M. Rezania, H. S. Yu, and H. B. Bai, “Variations of hydraulic properties of granular sandstones during water inrush: effect of small particle migration,” Engineering Geology, vol. 217, pp. 61–70, 2017.
View at: Publisher Site | Google Scholar
K. Peng, Z. P. Liu, Q. L. Zou, Q. H. Wu, and J. Q. Zhou, “Mechanical property of granite from different buried depths under uniaxial compression and dynamic impact: an energy-based investigation,” Powder Technology, vol. 362, pp. 729–744, 2020.
View at: Publisher Site | Google Scholar
Q. Wang, Z. H. Jiang, B. Jiang, H. K. Gao, Y. B. Huang, and P. Zhang, “Research on an automatic roadway formation method in deep mining areas by roof cutting with high-strength bolt-grouting,” International Journal of Rock Mechanics and Mining Sciences, vol. 128, article 104264, 2020.
View at: Publisher Site | Google Scholar
M. C. He, H. P. Xie, S. P. Peng, and Y. D. Jiang, “Study on rock mechanics in deep mining engineering,” Chinese Journal of Rock Mechanics and Engineering, vol. 24, pp. 2803–2812, 2005.
View at: Google Scholar
D. Xu, Mechanism of fault reactivation-induced disasters during coal mining above high water-pressure confined aquifers, [Ph.D. thesis], China University of Mining and Technology, Geological Resources and Geological Engineering, 2012.
L. P. Li, W. Lu, S. C. Li, Q. S. Zhang, Z. H. Xu, and S. S. Shi, “Research status and developing trend analysis of the water inrush mechanism for underground engineering construction,” Journal of Shandong University (Engineering Science), vol. 40, pp. 104–112, 2010.
View at: Google Scholar
V. N. Odintsev and N. A. Miletenko, “Water inrush in mines as a consequence of spontaneous hydrofracture,” Journal of Mining Science, vol. 51, no. 3, pp. 423–434, 2015.
View at: Publisher Site | Google Scholar
J. Q. Guo and C. S. Qiao, “Study on water-inrush mechanism and safe thickness of rock wall of karst tunnel face,” Journal of the China Railway Society, vol. 34, pp. 105–111, 2012.
View at: Google Scholar
Y. L. Zhao, L. Y. Zhang, J. Liao, W. J. Wang, Q. Liu, and L. Tang, “Experimental study of fracture toughness and subcritical crack growth of three rocks under different environments,” International Journal of Geomechanics, vol. 20, no. 8, article 04020128, 2020.
View at: Publisher Site | Google Scholar
Q. Zhang, M. C. He, J. Wang et al., “Instantaneous expansion with a single fracture: A new directional rock- breaking technology for roof cutting,” International Journal of Rock Mechanics and Mining Sciences, vol. 132, article 104399, 2020.
View at: Publisher Site | Google Scholar
X. H. Xie, B. Y. Su, Y. F. Gao, and X. B. Duan, “Numerical study on water inrush above a confined aquifer in coal mining using hydro-fracturing,” Chinese Journal of Rock Mechanics and Engineering, vol. 24, pp. 987–993, 2005.
View at: Google Scholar
R. Q. Huang, X. N. Wang, and L. S. Chen, “Hydro splitting off analysis on underground water in deep lying tunnels and its effect on water gushing out,” Chinese Journal of Rock Mechanics and Engineering, vol. 19, pp. 573–576, 2000.
View at: Google Scholar
B. C. Liu, Z. K. Xiang, Y. S. Lin, C. M. Lu, and B. H. Zhang, “Karst tunnel water gushing mechanism in Jinping stage 2 hydropower station,” Journal of Guilin University of Technology, vol. 28, pp. 484–488, 2008.
View at: Google Scholar
X. Y. Li, D. L. Zhang, Q. Fang, and H. R. Song, “On water burst patterns in underwater tunnels,” Modern Tunnelling Technology, vol. 52, pp. 24–31, 2015.
View at: Google Scholar
Z. H. Chen, Z. P. Hu, H. Li, and Q. F. Chen, “Fracture mechanical model and criteria of insidious fault water inrush in coal mines,” Journal of China University of Mining & Technology, vol. 40, pp. 673–677, 2011.
View at: Google Scholar
S. X. Yin and Q. Wu, “Generalized modes and academic criterions of water inrush from paleo-sinkholes,” Journal of University of Science and Technology Beijing, vol. 28, pp. 812–817, 2006.
View at: Google Scholar
Z.-H. Li, G. R. Feng, and C. Z. Zhai, “Study on “triangle” water-inrush mode of strong water-guide collapse column,” Journal of Central South University, vol. 23, no. 9, pp. 2402–2409, 2016.
View at: Publisher Site | Google Scholar
Z. L. Li, H. C. Zhang, Q. W. Ren, and Y. H. Wang, “Analysis of hydraulic fracturing and calculation of critical internal water pressure of rock fracture,” Rock and Soil Mechanics, vol. 26, pp. 1216–1220, 2005.
View at: Google Scholar
X. Dong, H. Lu, H. X. Huang, Y. Hao, and Y. Xia, “Failure mode of the water-filled fractures under hydraulic pressure in karst tunnels,” Open Geosciences, vol. 9, no. 1, pp. 186–196, 2017.
View at: Publisher Site | Google Scholar
S. C. Li, Y. C. Yuan, L. P. Li, Z. H. Ye, Q. Q. Zhang, and T. Lei, “Water inrush mechanism and minimum safe thickness of rock wall of karst tunnel face under blast excavation,” Chinese Journal of Geotechnical Engineering, vol. 37, pp. 313–320, 2015.
View at: Google Scholar
Y. L. Zhao, S. L. Luo, Y. X. Wang, W. J. Wang, L. Y. Zhang, and W. Wan, “Numerical analysis of karst water inrush and a criterion for establishing the width of water-resistant rock pillars,” Mine Water and the Environment, vol. 36, no. 4, pp. 508–519, 2017.
View at: Publisher Site | Google Scholar
S. X. Yin and Q. Wu, “Simulation and mechanism analysis of water inrush from karstic collapse columns in coal floor,” Chinese Journal of Rock Mechanics and Engineering, vol. 23, pp. 2551–2556, 2004.
View at: Google Scholar
J. X. Wang, Q. Li, W. X. Wu et al., “On water burst in deep-buried tunnelling based on hydrofracturing method,” Modern Tunnelling Technology, vol. 55, pp. 164–169, 2018.
View at: Google Scholar
Z. Huang, W. Zeng, Y. Wu, S. J. Li, and J. Zhao, “Experimental investigation of fracture propagation and inrush characteristics in tunnel construction,” Natural Hazards, vol. 97, no. 1, pp. 193–210, 2019.
View at: Publisher Site | Google Scholar
T. H. Yang, Study on infiltrate character and coupling analysis of seepage and stress in rock failure process, [Ph.D. thesis], Northeastern University, Engineering Mechanics, 2001.
L. C. Li, Q. M. Meng, S. Y. Wang, G. Li, and C. A. Tang, “A numerical investigation of the hydraulic fracturing behaviour of conglomerate in glutenite formation,” Acta Geotechnica, vol. 8, no. 6, pp. 597–618, 2013.
View at: Publisher Site | Google Scholar
S. Y. Wang, S. W. Sloan, S. G. Fityus, D. V. Griffiths, and C. A. Tang, “Numerical modeling of pore pressure influence on fracture evolution in brittle heterogeneous rocks,” Rock Mechanics and Rock Engineering, vol. 46, no. 5, pp. 1165–1182, 2013.
View at: Publisher Site | Google Scholar
Y. Wang, X. Li, R. Q. Zhou, B. Zheng, B. Zhang, and Y. F. Wu, “Numerical evaluation of the effect of fracture network connectivity in naturally fractured shale based on FSD model,” Science China-Earth Sciences, vol. 59, no. 3, pp. 626–639, 2016.
View at: Publisher Site | Google Scholar
T. Xu, P. G. Ranjith, A. S. K. Au et al., “Numerical and experimental investigation of hydraulic fracturing in kaolin clay,” Journal of Petroleum Science and Engineering, vol. 134, pp. 223–236, 2015.
View at: Publisher Site | Google Scholar
H. Q. Zhang, Y. N. He, C. A. Tang, B. Ahmad, and L. J. Han, “Application of an improved flow-stress-damage model to the criticality assessment of water inrush in a mine: a case study,” Rock Mechanics and Rock Engineering, vol. 42, no. 6, pp. 911–930, 2009.
View at: Publisher Site | Google Scholar
M. A. Biot, “General theory of three-dimensional consolidation,” Journal of Applied Physics, vol. 12, no. 2, pp. 155–164, 1941.
View at: Publisher Site | Google Scholar
W. Weibull, “A statistical distribution function of wide applicability,” Journal of Applied Mechanics, vol. 18, pp. 293–297, 1951.
View at: Google Scholar
X. Y. Shang and H. Tkalčić, “Point‐source inversion of small and moderate earthquakes from P‐wave polarities and P/S amplitude ratios within a hierarchical bayesian framework: implications for the geysers earthquakes,” Journal of Geophysical Research: Solid Earth, vol. 125, no. 2, 2020.
View at: Publisher Site | Google Scholar
K. Peng, J. Q. Zhou, Q. L. Zou, Y. J. Zhang, and G. W. Tan, “Deformation characteristics and failure modes of sandstones under discontinuous multi-level cyclic loads,” Powder Technology, vol. 373, pp. 599–613, 2020.
View at: Publisher Site | Google Scholar
Q. Zhang, Hydraulic fracturing effect of water inrush in excavation of roadways or tunnels, [M.S. thesis], China University of Mining and Technology, Mining Engineering, 2018.
J. Geertsma and F. de Klerk, “A rapid method of predicting width and extent of hydraulically induced fractures,” Journal of Petroleum Technology, vol. 21, pp. 1571–1581, 1969.
View at: Publisher Site | Google Scholar
M. X. Chu and J. G. Yan, “Definition and control of thickness of relaxed zone around mining gateway,” Coal Science and Technology, vol. 31, pp. 53–55, 2003.
View at: Google Scholar

Copyright

Copyright © 2020 Quan Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies