Abstract

Zonal disintegration phenomenon with alternative distribution of fracture zones and nonfracture zones is a characteristic failure in deep rock masses, especially for deep tunnel excavated by drill and blast method. To investigate the mechanism of zonal disintegration under coupled high axial geostress and blasting load, elastic stress field distribution for a circular tunnel in deep isotropic rock masses is obtained. Furthermore, Hoek-Brown criterion is amended by considering blast-induced damage effect. Both radial blasting load and blast-induced damage effect are assumed to decay in a negative exponential function. Taking the deep tunnel of Dingji coal mine in China as engineering background, the number and width of fracture zones are determined by deduced elastic stress field and modified Hoek-Brown criterion. Then numerical computation is conducted. Numerical computation results indicate that both the number and width of fracture zones mainly depend on high axial geostress and mechanical parameters of deep rock masses, and peak radial blasting load plays an important role in determining the width of fracture zone near the excavation.

1. Introduction

With the ever-increasing demand of energy resources, deep mineral resources exploitation is inevitable and imperative due to the continuing consumption of shallow mineral resources [1], and deep mining has been identified as an important topic for research under China’s State Key Research and Development Program [2]. During deep mineral resources exploitation, some special rock failure phenomena, such as zonal disintegration and rockburst, are frequently reported and quite different from those in shallow rock Engineering [35]. Among these special rock failure phenomena, zonal disintegration is such a unique phenomenon with alternative distribution of fracture zones and nonfracture zones around an excavation in deep rock mass [6, 7].

Zonal disintegration phenomenon is considered as one of the characteristics of deep rock masses. In 1970s, zonal disintegration phenomenon was firstly observed ahead of stope faces by periscope in Doornfontein gold mine of South Africa at about -2300 m level [8]. Then zonal disintegration phenomenon was also monitored around underground workings by electrometric techniques in Oktyabrsky Mine and Taymyrsky Mine of Russia at about -1000 m level [9]. However, an obvious zonal disintegration phenomenon in China was detected by bore TV video in Huainan mine at -910 m level [10].

Formation condition and mechanism of zonal disintegration under in situ stresses have been studied subsequently by theoretical analysis, model experiment, numerical simulation, and in situ monitoring. Shemyakin et al. [11] considered that the fracture zone is the splitting failure in rock mass. Odintsev [12] found that zonal disintegration can be predicted by spatial periodicity of splitting failure around an excavation. Zhang et al. [13] established an elastic damage-softening model for zonal disintegration according to strain gradient theory and continuum damage mechanics. Chen et al. [5] revealed that the circulation of fracture zone is mainly caused by the peak stress intensity factor of a “new” plastic boundary. Jia et al. [14] stated that failure patterns of surrounding rock mass are influenced by multiaxial stress state and a high axial geostress is essential for zonal disintegration. Yuan and Gu [15] indicated that the major radial tension strain following tension breakage is the key factor of zonal disintegration. While traditional static mechanics theory cannot illustrate the formation of zonal disintegration well, therefore the unloading effect during excavation process is considered. Zhu et al. [16] indicated that a slowly unloading P-wave reflecting from a free surface may induce zonal disintegration. By assuming deep rock masses containing a set of parallel cracks, zonal disintegration mechanism was analyzed and fracture zone occurred when cracks satisfy the crack grown criterion or the strain energy density factor theory [17, 18].

Based on non-Euclidean model, zonal disintegration mechanism around a deep circular tunnel was investigated and elastic stress field of surrounding rock was derived under various conditions, such as nonhydrostatic pressure state and dynamic unloading. In non-Euclidean model, radial stress and tangential stress around deep tunnel were obviously fluctuant, which was quite different with classic elastic solution [19, 20]. And zonal disintegration around an excavation is the result of the alternate appearance of stress wave crests and troughs [21, 22]. By establishing a new non-Euclidean model in a nonhydrostatic pressure state, it was found that the magnitude and position of fracture zones depend on the actual geostress state and non-Euclidean parameters [23]. Considering the characteristics of rock mass, both crack-weakened rock mass and cross-anisotropic rock mass were analyzed for zonal disintegration with a non-Euclidean model [24, 25]. Moreover, zonal disintegration around a deep circular tunnel under dynamic unloading condition was also studied with non-Euclidean model [26, 27]. When considering the excavation unloading effect, it was found that the magnitude and position of fracture zone mainly depend on unloading rate, and the quantity of fracture zones increases with decreasing excavation time [18, 22].

In consideration of complex geological environment and blasting excavation, rock failure in deep rock mass can be induced by blasting load [28, 29]. General Particle Dynamics (GPD) method was put forward to study the zonal disintegration mechanism of deep rock mass in blasting excavation, and it was found that dynamic loads and high geostress are two dominant factors for occurrence of zonal disintegration [30]. Blasting excavation not only applies a dynamic load on deep rock masses but also induces blasting damage around an excavation [31, 32]. By considering both high axial geostress and blasting load in deep rock masses, failure mechanism of zonal disintegration is investigated by elastic stress field analyses and modified Hoek-Brown criterion. Then the effects of high axial geostress, mechanical parameters of deep rock masses, and radial blasting load on zonal disintegration are analyzed.

2. Elastic Stress Field under Coupled High Axial Geostress and Blasting Load

2.1. Decomposition of Elastic Stress Field

It is found that the axial stress of tunnel is vital for zonal disintegration [33]. Therefore, high axial geostress should be considered during elastic stress field analyses.

It is assumed that a deep circular tunnel, whose radius is r0, is subjected to a vertical geostress p0, a horizontal geostress p0 at infinity, and a high axial geostress along the tunnel axis. Under coupled high axial geostress and blasting load, elastic stress field around a circular excavation in deep rock mass can be divided into two parts, elastic stress field caused by high axial geostress and transient elastic stress field induced by blasting load, which can be seen from Figure 1.

2.2. Elastic Stress Field Caused by High Axial Geostress

To simplify the analyses, three-dimensional problem is converted into a quasi-plane strain problem with an initial axial strain [34]. According to simplified quasi-plane strain state, axial stress component can be expressed aswhere E is Young’s modulus, μ is Poisson’s ratio, and and are radial stress component and tangential stress component, respectively.

As , , , and , the initial axial strain can be obtained as

According to the research of Qian et al. [34], the first invariant of stress tensor I1 in surrounding rock with a distance r from tunnel center can be obtained aswhere q in non-Euclidean model is the fitting parameter determined by experimental data, J0(ρ1/2r) and J1(ρ1/2r) are zero-order and first-order Bessel function, N0(ρ1/2r) and N1(ρ1/2r) are zero-order and first-order Neumann function, K0(ρ1/2r) and K1(ρ1/2r) are zero-order and first-order modified Bessel function of second kind, and C coincides with the Wronskian of linearly independent solutions J0(ρ1/2r) and N0(ρ1/2r).

According to elasticity mechanics, the first invariant of stress tensor I1 in surrounding rock with a distance r from tunnel center can also be written as

Then, the tangential stress component σθ1 can be given as follows:

Under high axial geostress, the equation of equilibrium can be expressed as

Integrating (9), the radial stress component can be expressed as

Substituting (10) into (8), the tangential stress component σθ1 can be expressed as

2.3. Transient Elastic Stress Field Induced by Radial Blasting Load

Elastic stress field induced by blasting load can be treated as an axial symmetric problem in plain strain state for long circular tunnel. As an impact load, blasting load usually increases firstly and then decreases. But the main cause of zonal disintegration phenomenon in deep rock mass is the tensile failure of surrounding rock in radial direction [15]. And the tensile failure in radial direction usually occurs in unloading stage rather than in loading stage. Therefore, a monotone decreasing radial blasting load p(t) is applied on the boundary of circular tunnel in consideration of blasting load in excavation [35]. As shown in Figure 2, radial blasting load decays in a negative exponential function with action time. And the expression for time-varying radial blasting load p(t) is shown as follows.where is the peak pressure of radial blasting load, a is exponential decay coefficient, and t0 is the action time of radial blasting load. ρ0 and D are density and detonation velocity of explosive.

Under blasting load, the equations of equilibrium can be written aswhere is longitudinal wave velocity: .

With the help of Laplace transform, the Laplace transform of u(r,t) can be expressed aswhere p is parameter of Laplace transform.

Then the equations of equilibrium can be written as

Assuming that , (16) can be obtained as

Equation (17) is modified Bessel equation, and its general solution iswhere I1(η) and K1(η) are first-order modified Bessel function of first kind and second kind and M and N are the parameters determined by boundary conditions.

As , . While U(η, p) is bounded, the parameter N is equal to 0. According to the asymptotic property of Bessel function, the first-order modified Bessel function of first kind I1(η) can be expressed as

Therefore, (18) can be written as

According to generalized Hooke’s law, the strain for plain strain problem in axisymmetric condition can be calculated as

By adding (23) multiplied by (1-μ)/μ with (24), the following can be obtained:

Then, the radial stress component can be rewritten as

By adding (23) multiplied by μ with (24) multiplied by (1-μ), the following can be obtained:

Then, the tangential stress component σθ2 can be rewritten as

Assuming that and , the stress-strain relation can be rewritten as

On the boundary of circular tunnel, the stress boundary condition can be expressed as

With the help of Laplace transform, the stress boundary condition can be written aswhere .

Substituting (22) into (32), the parameter M can be obtained as

According to and (34), (22) can be expressed as

According to convolution theorem of Laplace transform,

With the help of inverse Laplace transform, the radial displacement around circular tunnel is obtained as

Therefore, according to (29) and (30), the radial stress component and tangential stress component induced by radial blasting load can be obtained as

2.4. Elastic Stress Field under Coupled High Axial Geostress and Blasting Load

Consequently, the radial stress component and tangential stress component under coupled high axial geostress and blasting load can be expressed as

Substituting (2), (42), and (43) into (1), the axial stress component can be expressed as

3. Hoek-Brown Criterion Considering Blast-Induced Damage Effect

Assume that when deep rock mass meets the nonlinear Hoek-Brown criterion, deep rock mass fracture occurs. And the size and location of fracture zone can be obtained according to principal stress components. The Hoek-Brown criterion can be expressed as [36, 37]where σ1 and σ3 are the maximum and minimum principal stress, is uniaxial compressive strength of intact rock, and and s are strength parameters.

Strength parameters and s can be calculated by GSI (geological strength index). During blasting excavation, blast-induced damage plays down the integrity of rock mass and weakens the mechanical properties of rock mass. Hence, blast-induced damage factor D is introduced for cumulative damage effects of blasting excavation [38]. Then, the value of strength parameters and s can be obtained aswhere is the value of for intact rock.

With the increase of distance, blast-induced damage factor D decays in a negative exponential function [38], and it can be assumed as

4. Validation of Present Results

In China, an obvious zonal disintegration phenomenon was observed in Dingji coal mine at -910 m level. The distribution of facture zones is shown in Figure 3 [10]. According to research of Zhang et al. [39], an essential condition for zonal disintegration is maximum principal stress parallel to tunnel axis, exceeding 1.5 times uniaxial compressive strength of surrounding rock mass.

According to in situ condition, density of sandstone is 2620 kg/m3, Young’s modulus is 12.97 GPa, Poisson’s ratio is 0.268, and uniaxial compressive strength is 88.55 MPa. The height and width of straight wall arch tunnel are 3.88 m and 5.00 m, respectively. According to equivalent circular method [40], the equivalent circle radius is 2.54 m. In consideration of overburden rock mass, vertical geostress and horizontal geostress p0 at infinity are 25.10 MPa. High axial geostress is assumed to be uniaxial compressive strength of surrounding rock mass, 88.55 MPa. Non-Euclidean model q = 3.28×107 N·m2; C = 50 m−2. Considering rock structure, joints, and initial cracks, Hoek-Brown parameter is 10, and GSI is 75. Density and detonation velocity of explosive are 1100 kg/m3 and 3262 m/s, respectively. Peak pressure, decay coefficient, and action time of radial blasting load are 73.15 MPa, 0.87, and 8 ms, respectively. Blast-induced damage factor D decays in a negative exponential function.

Figure 4 presents the number and location of fracture zones from present numerical computation results. As seen from Figure 4, zonal disintegration arises in both conditions, while the width of fracture zones is different. On one hand, radial blasting load leads to stress redistribution around the deep tunnel. Then the elastic stress field is changed. On the other hand, radial blasting load produces some microcracks in surrounding rock with deterioration of its mechanical properties and damage to its integrity. It can be found that blasting load plays an important role in determining the width of fracture zones near the excavation. Comparing Figure 3 with Figure 4, the number and location of fractured zones under coupled high axial geostress and blasting load are in good accordance with the results observed in Dingji coal mine.

5. Numerical Computation Results

5.1. Effect of Geological Strength Index (GSI) on Distribution of Fracture Zones

Distribution of fractures zones under different geological strength index (GSI) conditions is shown in Figure 5. The other parameters are the same as mentioned in Section 4.

As seen from Figure 5, both the number and the width of fracture zones increase with the reduction of geological strength index (GSI).

5.2. Effect of Hoek-Brown Parameter on Distribution of Fracture Zones

Figure 6 presents the distribution of fractures zones under various Hoek-Brown parameter conditions when geological strength index (GSI) is 80. The other parameters are the same as mentioned in Section 4.

It can be seen from Figure 6 that when Hoek-Brown parameter is greater than 10, only the width of fracture zones increases with the reduction of Hoek-Brown parameter . While when Hoek-Brown parameter is less than 10, both the number and the width of fracture zones increase with the reduction of Hoek-Brown parameter .

5.3. Effect of High Axial Geostress on Distribution of Fracture Zones

As geological strength index (GSI) is 90 and Hoek-Brown parameter is equal to 15, the distribution of fractures zones under various high axial geostress conditions is shown in Figure 7. The other parameters are the same as mentioned in Section 4. High axial geostress varies from 1 times to 1.6 times of uniaxial compressive strength of surrounding rock mass.

As seen from Figure 7, both the number and the width of fracture zones increase with the increase of high axial geostress .

5.4. Effect of Peak Radial Blasting Load on Distribution of Fracture Zones

When geological strength index (GSI) is equal to 80 and Hoek-Brown parameter is equal to 10, the distribution of fractures zones under various peak radial blasting load conditions is shown in Figure 8. The other parameters are the same as mentioned in Section 4.

As seen from Figure 8, peak radial blasting load , which decays in a negative exponential function, shows a small influence on the number of fracture zones but a relatively big effect on the width of fracture zone near the excavation.

6. Conclusions

As a characteristic failure in deep rock masses, zonal disintegration phenomenon is such a unique failure phenomenon with alternative distribution of fracture zones and nonfracture zones. To study the formation mechanism of zonal disintegration under coupled high axial geostress and blasting load, the elastic stress field distribution for a circular tunnel in deep isotropic rock masses is deduced by assuming that radial blasting load decays in a negative exponential function. And Hoek-Brown criterion is modified by assuming that blast-induced damage effect decays in a negative exponential function. With the deep tunnel of Dingji coal mine in China as engineering background, the number and width of fracture zones are determined and they are in good accordance with the results observed in Dingji coal mine. Some numerical computations are carried out to study the effect of mechanical parameters of deep rock masses, high axial geo-stress, and radial blasting load on zonal disintegration. It is found that both the number and width of fracture zones mainly depend on high axial geostress and mechanical parameters of deep rock masses, and peak radial blasting load plays an important role in determining the width of fracture zone near the excavation. Both the number and the width of fracture zones increase with the increase of high axial geostress or with the deterioration of mechanical parameters of deep rock masses, while only the width of fracture zone near the excavation is extended with the increase of peak radial blasting load.

Data Availability

The datasets generated and analyzed during the current study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This research is funded by the National Natural Science Foundation of China (no. 51174004 and no. 51774011), Anhui Natural Science Foundation (no. 1808085QE148), Natural Science Research Project of Colleges and Universities in Anhui Province (no. KJ2017A097), Young Teacher Scientific Research Project of Anhui University of Science and Technology (no. QN201607), and Science and Technology Project of Department of Housing and Urban-Rural Development of Anhui Province (no. 2017YF-08).