Evaluation of Excavation-Damaged Zone around Underground Tunnels by Theoretical Calculation and Field Test Methods

: Excavation-damaged zones (EDZs) induced in underground mining and civil engineering potentially threaten tunnel safety and stability, and increase construction and support costs. In this paper, an investigation of the excavation damaged zone (EDZ) around roadways in Fankou lead-zinc mine in Guangzhou, China is performed by applying a seismic velocity method accompanied by SET-PLT-01 nonmetallic ultrasonic detector. Meanwhile, the in situ stress in the mining area was measured based on the stress relief method with the Swedish high-precision LUT system. The results indicate that the stress ﬁeld is dominated by the maximum horizontal tectonic stress, and the extents of the EDZ on the roof-ﬂoor region are greater than that on the sidewall. In addition, both of the in situ stresses and EDZs show an increasing trend with an increase of depth. Analytical solutions of EDZ around circular openings in the brittle rock mass subjected to non-hydrostatic stress ﬁelds are presented in terms of the Mohr–Coulomb and generalized Hoek–Brown criteria, and validated by several cases mentioned above. The extents of EDZ solved by closed-form solutions were found to be in a great agreement with those obtained in the ﬁeld. Finally, a series of parametric studies are conducted to investigate the e ﬀ ects of cohesion ( c ), friction angle ( ϕ ), geological strength index ( GSI ), m i , uniaxial compressive strength ( σ c ), and disturbance factor ( D ) on EDZ. It is shown that the e ﬀ ects of c , ϕ , GSI , and σ c are signiﬁcant; however, more attention should be paid to consider the dynamic disturbances induced by mechanical drilling, blasting, and seismic waves in tunnel excavations or operations. elastoplastic rock mass is subjected to a stress ﬁeld of vertical stress, σ 0 , and horizontal stress, λσ 0 . The support resistance and ﬂuid containment distributed on the tunnel wall are assumed to be zero and can be simpliﬁed as a plane strain problem in the cylindrical coordinate system.


Introduction
The supply of energy is related to the lifeblood of global economic development. As one of the important components of energy, mineral resources are rapidly decreasing in the shallow crust. Therefore, in order to ensure sustainable human production activities, a number of mines are gradually shifting to deep mining around the world. It is well known that in situ stress redistribution will be generated during the excavation progress of a man-made opening, then the tangential stress gradually increases while radial stress on the surrounding rock mass decreases. Both of them reach an extreme state on the surface of the openings and therefore stress concentrations occur in the surrounding of the openings. Once the concentrated stress exceeds the failure strength of the surrounding rocks, excavation-induced cracks are initiated and gradually expand deep into the rock mass until the stress reaches a new equilibrium condition. Eventually, an undesirable excavation-damaged zone (EDZ) where the basic physical, thermal-mechanical, and hydrogeological properties of surrounding rocks are greatly changed is induced in the rock mass surrounding the opening during the construction process [1][2][3][4]. Previous studies have indicated that geohazards, such as roof collapse, rheology, and floor sinking are related to EDZs [5][6][7][8], which threaten the safety and efficiency of deep underground mining. Therefore, it is important to develop a solution and measurement method to determine the mechanism of EDZ development and accurately evaluate the thickness of an EDZ.
Extensive studies have been conducted on evaluating the EDZs based on theoretical analysis, numerical simulation, and field test [9][10][11][12]. Wu et al. [13] introduced a simplified formula for assessing an EDZ according to the strain energy density released due to excavation. However, closed-form solutions of stress and displacement distributions around circular openings are mainly based on different constitutive models and failure criteria. The Mohr-Coulomb and Hoek-Brown failure criteria are simultaneously employed to analyze circular tunnel problems. For instance, the theoretically consistent closed-form solution of displacements around a circular opening under an elastic-brittle-plastic rock mass satisfying the Hoek-Brown failure criterion is proposed [14,15]. Carranza-Torres [16] developed a perfectly elastic-plastic constitutive model and effectively combined it with the Hoek-Brown criterion. Park et al. [17] proposed analytical solutions for predicting displacements around a circular opening in terms of both the Mohr-Coulomb and Hoek-Brown failure criteria. However, these criteria do not consider the effect of intermediate principal stress. Because of this deficiency, researchers have started using the Drucker-Prager criterion and unified strength theory to solve the elastic-plastic problem of surrounding rocks in circular openings. Based on the unified strength theory and the hypothesis of plane strain, Roatesi et al. [18] proposed a new closed-form solution and a specific expression for the circular tunnel problem in finite and infinite rock masses, respectively. Moreover, they analyzed the effect of the intermediate principal stress, radius-dependent Young's modulus, and dilatation on the stress and displacement fields of EDZs. The results indicate that stress and radial displacement decrease with the increase in the intermediate principal stress parameter, b. More recently, Singh et al. [19] proposed analytical solutions for the circular opening problem in an elastoplastic rock mass under the Mogi-Coulomb, inscribed Drucker-Prager, middle circumscribed Drucker-Prager, circumscribed Drucker-Prager, and Mohr-Coulomb failure criteria. In addition, considering the effect of the intermediate principal stress and selection of the appropriate failure criteria have a significant effect on optimizing the support system. Several studies have investigated the effect of out-of-plane stresses on tunnel problems in elastoplastic rock masses [20][21][22][23]. Based on the nature of a non-plane strain problem, an analytical solution of the cylindrical opening subjected to various in situ axial stresses can be derived by employing the Mohr-Coulomb failure criterion [24].
In conjunction with the development of high-speed computer algorithms, numerical simulation, as an effective method in geotechnical engineering [25][26][27][28], has been widely applied to investigate EDZs. Golshani et al. [29] simulated the excavation damage of a circular opening in a brittle rock mass in terms of an extended micromechanics-based continuum damage model. The results clearly show that the development of an EDZ is a function of time and that moisture has a significant effect on its development. Moreover, an EDZ formed in the sidewalls is larger when the vertical stress is greater than the horizontal stress. Yang et al. [30] studied the relationship between EDZ size and the excavation process, lining, rock mass structure, initial stress, and tunnel radius. Perras et al. [31] evaluated the validity of the damage initiation and spalling limit approach for predicting the depth of the damaged area, as well as the predictive equations for the EDZ range of circular tunnels in the brittle rock mass. Based on the Flac3d platform, Ren et al. [32] investigated the surrounding rock displacement field, stress field, and plastic damage zone by developing an elastic-plastic damage numerical method for fractured rock mass, which can be effectively applied to optimize the excavation scheme of URL. Other numerical studies show that the effect of lateral stress ratio and high in-situ stress are closed related to the excavation damaged zone [33,34]. As discussed above, several analytical and numerical solutions have been introduced to assess and characterize EDZs around circular roadways under hydrostatic pressure. However, knowledge available in the literature concerning closed-form solutions of elastoplastic problems around circular openings subjected to asymmetrical stress conditions is still limited. Therefore, the present study attempted to evaluate the EDZ in Fankou lead-zinc mine in terms of theoretical calculation and field test methods. The investigation of the distribution law of in situ stress field and EDZ in the mining area is conducted. Closed-form solutions are derived for predicting the radius of EDZ around a circular opening in brittle rock mass under non-hydrostatic stress field based on the Mohr-Coulomb and generalized Hoek-Brown failure criteria. Then the accuracy of presented analytical solutions are validated by comparing the thickness of EDZ measured in the field with the single-hole seismic velocity method. In addition, a parametric study of c, ϕ, GSI, m i , D, and σ c on the damaged zone is carried out. The results of this study show that both the in situ stresses and the thickness of EDZ increase with the depth. Great agreement between the measured and calculated thickness of EDZ is obtained. The findings can provide measures to control damage zones and can also provide theoretical guidance for underground tunnel excavation and support.

Closed-Form Solutions of EDZ around Circular Tunnels
In general, closed-form solutions for circular tunnel problems are derived on the assumption that a tunnel is subjected to a hydrostatic stress field. Nevertheless, the initial stress in the surrounding rock mass has to be considered as inhomogeneous in most practical situations; that is, the lateral pressure coefficient is not equal to (λ 1) one. Two expressions are presented in this paper for predicting the extent of an EDZ around a circular tunnel in an infinite elastoplastic rock mass, which is governed by a non-hydrostatic pressure and obeys both the Mohr-Coulomb and Hoek-Brown failure criteria.
As shown in Figure 1, where r is the radial distance from the center of the opening surface, and R p is the plastic zone radius, a circular tunnel of radius R 0 in an infinite elastoplastic rock mass is subjected to a stress field of vertical stress, σ 0 , and horizontal stress, λσ 0 . The support resistance and fluid containment distributed on the tunnel wall are assumed to be zero and can be simplified as a plane strain problem in the cylindrical coordinate system. ethods. The investigation of the distribution law of in situ stress field and EDZ in the m is conducted. Closed-form solutions are derived for predicting the radius of EDZ arou lar opening in brittle rock mass under non-hydrostatic stress field based on the Mohr-Cou eneralized Hoek-Brown failure criteria. Then the accuracy of presented analytical solution ated by comparing the thickness of EDZ measured in the field with the single-hole se ity method. In addition, a parametric study of c, φ, GSI, mi, D, and σc on the damaged zo d out. The results of this study show that both the in situ stresses and the thickness of ase with the depth. Great agreement between the measured and calculated thickness of ED ned. The findings can provide measures to control damage zones and can also pro etical guidance for underground tunnel excavation and support.

sed-Form Solutions of EDZ around Circular Tunnels
n general, closed-form solutions for circular tunnel problems are derived on the assum a tunnel is subjected to a hydrostatic stress field. Nevertheless, the initial stress in unding rock mass has to be considered as inhomogeneous in most practical situations; th teral pressure coefficient is not equal to (λ ≠ 1) one. Two expressions are presented in this p redicting the extent of an EDZ around a circular tunnel in an infinite elastoplastic rock m is governed by a non-hydrostatic pressure and obeys both the Mohr-Coulomb and H n failure criteria. s shown in Figure 1, where r is the radial distance from the center of the opening surface the plastic zone radius, a circular tunnel of radius R0 in an infinite elastoplastic rock ma cted to a stress field of vertical stress, σ0, and horizontal stress, λσ0. The support resistance containment distributed on the tunnel wall are assumed to be zero and can be simplified strain problem in the cylindrical coordinate system. n fact, when r > Rp, the radial and tangential stresses are independent of the rock mass fa ion and are only related to the initial stress state of the rock mass. Therefore, the stresses i c zone can be obtained directly based on the theory of elasticity: In fact, when r > R p , the radial and tangential stresses are independent of the rock mass failure criterion and are only related to the initial stress state of the rock mass. Therefore, the stresses in the elastic zone can be obtained directly based on the theory of elasticity: where σ e r and σ e θ are the radial and tangential stresses in the elastic zone, respectively, σ i represents the radial stress on the elastic-plastic interface, and S 1 = (1 + λ) σ 0 /2, S 2 = (1 − λ) σ 0 /2, and µ = R p /r.
The stresses in the plastic zone can be solved in terms of the failure criteria and equilibrium Equation [2]:

Extent of the Plastic Zone in Terms of the Mohr-Coulomb Failure Criterion
In 1900, Mohr and Coulomb proposed the Mohr-Coulomb strength theory, which can characterize changes in rock strength with the surrounding rock pressure. It has been employed extensively in expressing the shear failure of rock materials because of its remarkable advantages, specifically, its simplicity and wide variety of classical plasticity solutions. For an active loading condition, the Mohr-Coulomb failure criterion can be written as [35] where ; and c and ϕ represent the cohesion and friction angle of the rock mass, respectively. By substituting Equation (3) into (2), then combining them with the boundary condition σ r = 0 when r = R 0 , the expressions for radial and tangential stresses in the plastic zone are derived as On the elastic-plastic interface (r = R p ), the stress states σ i = σ r and σ θ = 2S 1 − σr + 4S 2 cos2θ, satisfy the Mohr-Coulomb failure criterion. Thus, by substituting σ θ into Equation (3), σ i is given by With the continuous conditions of stress on the elastic-plastic zone interface (σ p r = σ e r , σ p θ = σ e θ ), the extent of the plastic zone around the circular tunnel in the rock mass governed by the Mohr-Coulomb failure can be obtained by combining Equations (1) and (4) with (5): where A = 1−sinϕ and B = 2c·cotϕ. concept of disturbed and undisturbed rock masses in 1992. The generalized Hoek-Brown criterion is given by [36][37][38].

Extent of Plastic
where σ 1 and σ 3 represent the major principal stress and minor principal stress at failure, respectively; σ c is the uniaxial compressive strength of the intact rock; m g , δ, and α are the Hoek-Brown constants, which depend on the properties of the rock mass. In cylindrical coordinates, for an active condition, Equation (7) can be written as The coefficients m g , δ, and α are calculated based on the empirical index called the geological strength index (GSI). This index can be quantified from charts based on the quality of the rock structure and condition of rock surfaces. The rock mass parameters are expressed according to the Carranza-Torres model [14] In Equations (9) and (10), m i is the coefficient, D is the factor that reflects the degree of disturbance to which the rock mass has been subjected because of blast damage and stress relaxation. It varies between 0 and 1. Actually, in the analytical solution of a circular tunnel problem, deducing when α > 0.5 is not straightforward. Therefore, in a relatively good quality rock mass, it is assumed that α = 0.5.
The stresses in the radial and tangential directions in the plastic region can be described according to Equations (2) and (8)- (11). Then combining them with the boundary condition σ r = 0 when r = R 0 , the following is obtained: where C = m g ln(r/R 0 ) and E = σ c /m g . At the elastic-plastic interface (r = R p ), σ i = σ r , and σ θ = 2S 1 − σ r + 4S 2 cos2θ. The radial stress is obtained in terms of the Hoek-Brown failure criterion through the substitution of σ θ into Equations (8)-(11): where K = S 1 + 2S 2 cos2θ and J = (m 2 g + 16δ+16m g K/σ c ) 0.5 . With the continuous conditions of stress at the elastic-plastic interface (σ p r = σ e r , σ p θ = σ e θ ), a closed-form solution of the extent of damage in the plastic zone within the surrounding rock mass, which satisfies the Hoek-Brown failure criterion, can be obtained by combining Equations (1) and (12) with (13): where H = δ+m 2 g /8 and F = 2/m g .

Fankou Lead-Zinc Mine in Situ Stress Distribution Laws
An investigation of the distribution law of in situ stress field was performed in Fankou lead-zinc mine in Guangzhou, China where the average elevation is approximately 100 to 200 m. The mine is currently operating with the vertical crater retreat (VCR) mining method at depths from 80 m to 880 m below the ground surface. The ore contains a large amount of yellow iron lead-zinc ores incorporating massive pyrite, the main components of the surrounding rock are massive and porphyritic limestone, dolomitic marl, and nodular limestone. The basic physical and mechanical properties of the surrounding rocks and ores were determined by laboratory tests of specimens extracted from exploration drilling cores and rock blocks taken directly from Fankou lead-zinc mine, experimental results showed that the uniaxial compressive strength, Poisson's ratio, and Platts' coefficient of rocks in this mining area were approximately 130 MPa, 0.25, and 8 to 13, respectively. Evidently, their stability was generally good, and typically, there was no need to design a support system for the upper ore body.
Numerous measurements of in situ stresses in the mining areas were carried out by applying the Swedish high-precision LUT rock triaxial stress measurement system which is with a basic principle of overcoring. Based on the assumption of linear and elastic deformations in the surrounding rock mass around boreholes, as well as the coincidence of loading and unloading curves, the main steps of field measurement are briefly described below. Boreholes deepened into the rock mass is firstly drilled, then twelve strain gauges are stuck firmly on the surface of drilling borehole. Meanwhile, the initial strain of each strain gauge is set to zero. Secondly, recovering the elastic deformation of the cores by releasing the surrounding rock stress with the overcoring method and measuring the elastic recovery strain of the cores. Finally, according to the measured elastic recovery strain and elastic constant of the rock mass, the stress conditions of each measurement site, including the stress component and principal stress, are analyzed by employing a specific LUT-str rock mass stress calculation program, the magnitude and direction of in situ stresses can be solved out with high accuracy.
The maximum horizontal principal stress, σ h-max ; minimum horizontal principal stress, σ h-min ; and vertical principal stress, σ v , were calculated by inputting the drilling parameters, elastic parameters, and recovered strain into the LUT-str calculation systems. Moreover, the ratio of σ h-max to σ h-min and the ratio of σ h-max to σ v were also calculated, as summarized in Table 1. Table 1. In situ stresses of each measurement sites. It can be seen from Table 1 that the maximum horizontal, minimum horizontal, and vertical principal stresses present an increasing trend in the mining area, and the in situ stress field in the Fankou lead-zinc mine was dominated by horizontal tectonic stresses. At the same overburden depth, because there was no evident difference among the magnitudes of the same type of stresses at different measurement points, the in situ stresses in the mining area were relatively uniform. In addition, the ratio of σ h-max to σ v in the mining area was approximately 1.80, and it gradually decreased as the excavation depth increased. The results indicate that as the depth increased, the effect of the gravitational stress on in situ stresses became more significant. On the other hand, it is noteworthy that the ratio of σ h-max to σ h-min was approximately 1.81 to 2.44. This revealed that the distribution law of horizontal tectonic stress in the rock mass is related to the depth. In order to visualize the relationship between in situ stresses and depth, the least-square method was employed to conduct a linear regression analysis of the stress data measured in the field. The regression equations of the relationship between σ h.max , σ h,min , and σ v , and the depth of the measurement stations are expressed as

No. Depth/m σ h-max /MPa
where H is the overburden depth of measurement sites.

Test Method and Apparatus
So far, many in situ test methods of EDZs have been invented to adapt to complex geological conditions [13,[39][40][41]. Whereas, because of better adaptability and operability, nonmetallic ultrasonic detection techniques are still widely used in mining and civil engineering [42]. In this study, the single-hole acoustic wave monitoring method incorporating a SET-PLT-01 nonmetallic ultrasonic detector was applied. This testing system is mainly comprised of a host computer, a water container, an ultrasonic source probe, and two receiving probes. The nonmetallic ultrasonic detector is relatively easy to operate and highly adaptable. In particular, it improves the testing environment and overcomes the shortcoming of the fully water-coupled acoustic wave test method, which cannot be used in a severely fractured rock mass. A simple illustration of the principle of the single-hole ultrasonic detection technique is shown in Figure 2. as the excavation depth increased. The results indicate that as the depth increased, the effect of the gravitational stress on in situ stresses became more significant. On the other hand, it is noteworthy that the ratio of σh-max to σh-min was approximately 1.81 to 2.44. This revealed that the distribution law of horizontal tectonic stress in the rock mass is related to the depth. In order to visualize the relationship between in situ stresses and depth, the least-square method was employed to conduct a linear regression analysis of the stress data measured in the field. The regression equations of the relationship between σh.max, σh,min, and σv, and the depth of the measurement stations are expressed as where H is the overburden depth of measurement sites.

Test Method and Apparatus
So far, many in situ test methods of EDZs have been invented to adapt to complex geological conditions [13,[39][40][41]. Whereas, because of better adaptability and operability, nonmetallic ultrasonic detection techniques are still widely used in mining and civil engineering [42]. In this study, the single-hole acoustic wave monitoring method incorporating a SET-PLT-01 nonmetallic ultrasonic detector was applied. This testing system is mainly comprised of a host computer, a water container, an ultrasonic source probe, and two receiving probes. The nonmetallic ultrasonic detector is relatively easy to operate and highly adaptable. In particular, it improves the testing environment and overcomes the shortcoming of the fully water-coupled acoustic wave test method, which cannot be used in a severely fractured rock mass. A simple illustration of the principle of the single-hole ultrasonic detection technique is shown in Figure 2. The source and receiving probes were inserted into the boreholes. This triggered the source probes to generate ultrasonic waves, which were refracted several times and received by probes S1 and S2. The waves were finally converted into electrical signals and stored in the host computer. Then data from all monitored points along each borehole were processed to obtain the depth and wave velocity curves. The formulation of the wave velocity at each monitored point can be expressed as The source and receiving probes were inserted into the boreholes. This triggered the source probes to generate ultrasonic waves, which were refracted several times and received by probes S1 and S2. The waves were finally converted into electrical signals and stored in the host computer. Then data from all monitored points along each borehole were processed to obtain the depth and wave velocity curves. The formulation of the wave velocity at each monitored point can be expressed as where t 1 and t 2 represent the initial times when receiving probes S 1 and S 2 start receiving signals, respectively, and L is the distance between the two receiving probes. The velocity of the ultrasonic waves decreased with the development of fractures. That is, the wave velocity was high if the surrounding rocks were intact; otherwise, the velocity was low. Hence, the thicknesses of EDZs around roadways can be determined by analyzing wave velocity variations at various locations.

Field Measurements
A total of five monitoring roadways with 20 test sites in the mining area were selected for EDZ assessment. These roadways were arch-shaped and had neither support nor serious disturbance. The locations and borehole parameters of these detection stations are summarized in Table 2.
where t1 and t2 represent the initial times when receiving probes S1 and S2 start receiving signals, respectively, and L is the distance between the two receiving probes.
The velocity of the ultrasonic waves decreased with the development of fractures. That is, the wave velocity was high if the surrounding rocks were intact; otherwise, the velocity was low. Hence, the thicknesses of EDZs around roadways can be determined by analyzing wave velocity variations at various locations.

Field Measurements
A total of five monitoring roadways with 20 test sites in the mining area were selected for EDZ assessment. These roadways were arch-shaped and had neither support nor serious disturbance. The locations and borehole parameters of these detection stations are summarized in Table 2. Four test boreholes were symmetrically drilled on the left and right sides of the roadways. The borehole depths were from 2.0 to 3.2 m, with a diameter of 51 mm and an inclination angle of 30°. The measurements were performed along the detecting borehole at an interval of 0.1 m from top to bottom. The arrangement of the boreholes monitored with the SET-PLT-01 nonmetallic ultrasonic detector around the roadways is shown in Figure 3. The main steps involved in field-testing with the application of the dry single-hole ultrasonic detection technology are briefly introduced below. The drilling boreholes are firstly flushed with high-pressure water on account of the rough surface and a considerable amount of debris, then water is injected into the charging pipe by compressing the container until water overflowed continuously from the vent hole. Simultaneously, the vent is sealed with screws, and gas is prevented from entering  The main steps involved in field-testing with the application of the dry single-hole ultrasonic detection technology are briefly introduced below. The drilling boreholes are firstly flushed with high-pressure water on account of the rough surface and a considerable amount of debris, then water is injected into the charging pipe by compressing the container until water overflowed continuously from the vent hole. Simultaneously, the vent is sealed with screws, and gas is prevented from entering the pipe. Secondly, the valve is opened to allow water to return to the container. Thereafter, the monitoring detector is thrust into the borehole and water is supplied into it until the rubber capsule is fully in contact with the surface of the borehole. Thirdly, to ensure accurate data collection, parameters such as the pulse width and frequency were adjusted. Finally, the data should be collected and saved when the signals become stable. It should be noted that the first measurement point is set at a distance of 0.2 to 0.4 m away from the sidewall because of the extreme fragmentation of rocks near the roadway surface.

Measurement Results and EDZ Analysis
Variation curves of the wave velocity with depth, shown in Figure 4, were obtained by loading the saved data into the special analysis software. According to these variation curves, the ultrasonic wave velocity increased as the borehole depth increased. In fact, the wave velocity was low because of the large number of cracks induced near the orifice; in contrast, it was high and stable when the monitoring points were located considerably far from the orifice until the joints closed. the pipe. Secondly, the valve is opened to allow water to return to the container. Thereafter, the monitoring detector is thrust into the borehole and water is supplied into it until the rubber capsule is fully in contact with the surface of the borehole. Thirdly, to ensure accurate data collection, parameters such as the pulse width and frequency were adjusted. Finally, the data should be collected and saved when the signals become stable. It should be noted that the first measurement point is set at a distance of 0.2 to 0.4 m away from the sidewall because of the extreme fragmentation of rocks near the roadway surface.

Measurement Results and EDZ Analysis
Variation curves of the wave velocity with depth, shown in Figure 4, were obtained by loading the saved data into the special analysis software. According to these variation curves, the ultrasonic wave velocity increased as the borehole depth increased. In fact, the wave velocity was low because of the large number of cracks induced near the orifice; in contrast, it was high and stable when the monitoring points were located considerably far from the orifice until the joints closed.    As seen in Figure 4, there was an evident low wave velocity region near the roadway surface. Figure 4a shows that the velocity at depths of 0.2 to 0.8 m varied in the low range from 2800 to 4000 m/s with a relatively slow increase rate, indicating that the rocks were severely damaged, and fractures were well developed in the area. However, there was a slight decline in the acoustic velocity, and subsequently, the wave velocity increased sharply within a short period. When the detection As seen in Figure 4, there was an evident low wave velocity region near the roadway surface. Figure 4a shows that the velocity at depths of 0.2 to 0.8 m varied in the low range from 2800 to 4000 m/s with a relatively slow increase rate, indicating that the rocks were severely damaged, and fractures were well developed in the area. However, there was a slight decline in the acoustic velocity, and subsequently, the wave velocity increased sharply within a short period. When the detection devices were set further down the borehole, the velocity slightly fluctuated and stabilized between 5200 and 5400 m/s. It can be concluded that a damaged zone where the wave velocity had an overall rising stage existed. In addition, a few micro-cracks may have been initiated in rock masses located considerably far from the excavation, but these remained intact. The other variation curves show that changes in the regulations of wave velocity were similar to those shown in Figure 4a. These results indicate that the extent of excavation damage gradually increased as the embedded depth of the roadway increased; that is, the thickness of the EDZ is positively related to the buried depth. On the other hand, it can be shown that the EDZ monitored in boreholes 1# and 4# was smaller than that in boreholes 2# and 3#. The variation of the damage thickness with buried depth was fitted by using the least square method, and the result is shown in Equation (19): where T EDZ is the EDZ thickness, and H is the embedded depth of the roadway.

Theoretical Calculations of EDZ Extent in Terms of the Mohr-Coulomb and Hoek-Brown Criteria
In Section 2, a closed-form solution of EDZ by the Mohr-Coulomb and Hoek-Brown criteria was proposed. Based on the practical conditions of the Fankou lead-zinc mine, the theoretical EDZ thickness was calculated for comparison with field measurement results. The surrounding rock masses were blocky, porphyritic limestone, and dolomitic marl. The parameters of the rock mass and material constants for the Mohr-Coulomb and Hoek-Brown failure criteria, determined by site investigations and laboratory experiments, are listed in Tables 3 and 4, respectively.

Results Based on the Mohr-Coulomb Failure Criterion
Based on Equation (6), the EDZ thickness is equal to the deviator radius (R p − R 0 ). The calculated EDZ ranges in the surrounding rock mass of underground roadways are presented in Figure 5. It is shown from Figure 5 that the EDZ size was related to inclination; the EDZ in the vertical direction was always larger than that in the horizontal direction because of the non-hydrostatic stresses field. Corresponding to the field roadways in the Fankou lead-zinc mine, EDZs on the sidewalls are smaller than those on the roof. In addition, the EDZ extent gradually increased as the embedded depth increased. When the depth was 780 m, the extents of damage in the vertical and horizontal directions were 1.79 and 1.19 m, respectively, and their ratio was 1.50. On the other hand, the thicknesses in the vertical direction were 0.835, 0.836, 0.886, 0.816, and 0.783 times the average values measured in the field at mining depths of 370, 450, 530, 680, and 780 m, respectively. It was found that the EDZ thicknesses calculated based on the Mohr-Coulomb criterion were smaller than the actual values measured in rock masses under various geological and geometrical conditions. In fact, once a tunnel is excavated in a deep underground mine, a damaged zone will occur in the surrounding rock masses as a result of redistribution of in situ stress and dynamic disturbance caused by excavation unloading. Our results demonstrated that predicting the extent of an EDZ in terms of the Mohr-Coulomb failure criterion has a high reliability.

Results Based on the Hoek-Brown Failure Criterion
By substituting the Hoek-Brown constants and parameters of roadways listed in Tables 3 and 4 into Equation (14), the ranges of EDZs in the surrounding rock mass around underground roadways were obtained, as shown in Figure 6. The results also showed that the EDZ contours obtained in terms of the Hoek-Brown criterion were elliptical. In other words, the EDZ in the vertical direction was always larger than that in the It is shown from Figure 5 that the EDZ size was related to inclination; the EDZ in the vertical direction was always larger than that in the horizontal direction because of the non-hydrostatic stresses field. Corresponding to the field roadways in the Fankou lead-zinc mine, EDZs on the sidewalls are smaller than those on the roof. In addition, the EDZ extent gradually increased as the embedded depth increased. When the depth was 780 m, the extents of damage in the vertical and horizontal directions were 1.79 and 1.19 m, respectively, and their ratio was 1.50. On the other hand, the thicknesses in the vertical direction were 0.835, 0.836, 0.886, 0.816, and 0.783 times the average values measured in the field at mining depths of 370, 450, 530, 680, and 780 m, respectively. It was found that the EDZ thicknesses calculated based on the Mohr-Coulomb criterion were smaller than the actual values measured in rock masses under various geological and geometrical conditions. In fact, once a tunnel is excavated in a deep underground mine, a damaged zone will occur in the surrounding rock masses as a result of redistribution of in situ stress and dynamic disturbance caused by excavation unloading. Our results demonstrated that predicting the extent of an EDZ in terms of the Mohr-Coulomb failure criterion has a high reliability.

Results Based on the Hoek-Brown Failure Criterion
By substituting the Hoek-Brown constants and parameters of roadways listed in Tables 3 and 4 into Equation (14), the ranges of EDZs in the surrounding rock mass around underground roadways were obtained, as shown in Figure 6. It is shown from Figure 5 that the EDZ size was related to inclination; the EDZ in the vertical direction was always larger than that in the horizontal direction because of the non-hydrostatic stresses field. Corresponding to the field roadways in the Fankou lead-zinc mine, EDZs on the sidewalls are smaller than those on the roof. In addition, the EDZ extent gradually increased as the embedded depth increased. When the depth was 780 m, the extents of damage in the vertical and horizontal directions were 1.79 and 1.19 m, respectively, and their ratio was 1.50. On the other hand, the thicknesses in the vertical direction were 0.835, 0.836, 0.886, 0.816, and 0.783 times the average values measured in the field at mining depths of 370, 450, 530, 680, and 780 m, respectively. It was found that the EDZ thicknesses calculated based on the Mohr-Coulomb criterion were smaller than the actual values measured in rock masses under various geological and geometrical conditions. In fact, once a tunnel is excavated in a deep underground mine, a damaged zone will occur in the surrounding rock masses as a result of redistribution of in situ stress and dynamic disturbance caused by excavation unloading. Our results demonstrated that predicting the extent of an EDZ in terms of the Mohr-Coulomb failure criterion has a high reliability.

Results Based on the Hoek-Brown Failure Criterion
By substituting the Hoek-Brown constants and parameters of roadways listed in Tables 3 and 4 into Equation (14), the ranges of EDZs in the surrounding rock mass around underground roadways were obtained, as shown in Figure 6. The results also showed that the EDZ contours obtained in terms of the Hoek-Brown criterion were elliptical. In other words, the EDZ in the vertical direction was always larger than that in the The results also showed that the EDZ contours obtained in terms of the Hoek-Brown criterion were elliptical. In other words, the EDZ in the vertical direction was always larger than that in the horizontal direction. It was concluded that the maximum calculated EDZ thicknesses were 1.064, 1.367, 1.519, 1.763, and 1.917 m at each depth, from shallow to deep. On the other hand, the average EDZ thicknesses measured in the field were 1.118, 1.086, 1.033, 1.013, and 1.129 times the average values calculated based on the Hoek-Brown criterion in the vertical direction. The deviations between the theoretical calculated EDZ thickness and field measured EDZ thickness did not exceed 12%. Thus, this failure criterion has high accuracy in predicting EDZ thickness. Because the Hoek-Brown criterion is a semi-empirical strength criterion, the EDZ radius was determined based on rock mass parameters. Consequently, changes in the rock mass parameters were related to significant differences among EDZs.
The calculated thicknesses based on the Hoek-Brown criterion were 6.6 to 11.7% larger than those calculated based on the Mohr-Coulomb criterion, which demonstrates that the Hoek-Brown criterion achieved better performance than the Mohr-Coulomb criterion in evaluating EDZs surrounding the Fankou lead-zinc mine. In addition, the EDZs assessed with the single-hole acoustic wave monitoring method in the field were larger than those calculated with both the Mohr-Coulomb and Hoek-Brown criteria. These differences may be attributed to the deficient working conditions in the field and the fact that the rock masses in the mining area were not perfectly elastoplastic materials. Additionally, the rock masses were inevitably subjected to dynamic disturbances, such as mechanical drilling, blasting vibration, and stress relief, during excavation. Simultaneously, with the redistribution of in situ stress, the deformations and fractures in the surrounding rock masses continuously deteriorated. Nevertheless, it is important to note that both theoretical solutions in terms of the Hoek-Brown and Mohr-Coulomb failure criteria did not consider the effect of the creep behavior of the rock. Moreover, only a few exact analytical solutions for non-circular tunnel problems have been presented in the literature. Therefore, the difference was inevitable when employing closed-form solutions, derived on the basis of circular openings to evaluate EDZs around these roadways, which have arch and rectangular shapes.

Affecting of EDZ Thickness
A parametrical study was carried out for the sake of thoroughly investigating the influence of the parameters of Mohr-Coulomb and Hoek-Brown criteria on the stability of excavations in a rock mass with different physical and mechanical properties. Therefore, based on those findings which can provide guidance for practical engineering and provide insight into the characteristics of the rock mass, more rational excavation methods and more stable support systems can be developed. By analyzing Equations (6) and (14), but without considering project conditions, such as λ and the tunnel radius, it was found that two key physical coefficients, namely, cohesion and friction angle, affected the calculated EDZ thickness, based on the Mohr-Coulomb criterion. Moreover, four variables, namely, GSI, m i , D, and σ c affect the EDZ thickness according to the Hoek-Brown criterion. In the following, the effect of six variables, specifically, c, ϕ, GSI, m i , D, and σ c , on the calculated EDZ thickness are discussed in relation to tunnel construction and support. The other characteristics and stress states of rock masses were assumed to be the same as those of the five roadways in the Fankou lead-zinc mine. Moreover, the following radius of the plastic zone was calculated at θ = 90 • . The calculated thicknesses based on the Hoek-Brown criterion were 6.6 to 11.7% larger than those calculated based on the Mohr-Coulomb criterion, which demonstrates that the Hoek-Brown criterion achieved better performance than the Mohr-Coulomb criterion in evaluating EDZs surrounding the Fankou lead-zinc mine. In addition, the EDZs assessed with the single-hole acoustic wave monitoring method in the field were larger than those calculated with both the Mohr-Coulomb and Hoek-Brown criteria. These differences may be attributed to the deficient working conditions in the field and the fact that the rock masses in the mining area were not perfectly elastoplastic materials. Additionally, the rock masses were inevitably subjected to dynamic disturbances, such as mechanical drilling, blasting vibration, and stress relief, during excavation. Simultaneously, with the redistribution of in situ stress, the deformations and fractures in the surrounding rock masses continuously deteriorated. Nevertheless, it is important to note that both theoretical solutions in terms of the Hoek-Brown and Mohr-Coulomb failure criteria did not consider the effect of the creep behavior of the rock. Moreover, only a few exact analytical solutions for non-circular tunnel problems have been presented in the literature. Therefore, the difference was inevitable when employing closed-form solutions, derived on the basis of circular openings to evaluate EDZs around these roadways, which have arch and rectangular shapes.

Affecting of EDZ Thickness
A parametrical study was carried out for the sake of thoroughly investigating the influence of the parameters of Mohr-Coulomb and Hoek-Brown criteria on the stability of excavations in a rock mass with different physical and mechanical properties. Therefore, based on those findings which can provide guidance for practical engineering and provide insight into the characteristics of the rock mass, more rational excavation methods and more stable support systems can be developed. By analyzing Equations (6) and (14), but without considering project conditions, such as λ and the tunnel radius, it was found that two key physical coefficients, namely, cohesion and friction angle, affected the calculated EDZ thickness, based on the Mohr-Coulomb criterion. Moreover, four variables, namely, GSI, mi, D, and σc affect the EDZ thickness according to the Hoek-Brown criterion. In the following, the effect of six variables, specifically, c, φ, GSI, mi, D, and σc, on the calculated EDZ thickness are discussed in relation to tunnel construction and support. The other characteristics and stress states of rock masses were assumed to be the same as those of the five roadways in the Fankou lead-zinc mine. Moreover, the following radius of the plastic zone was calculated at θ = 90°.
(a) (b)  As shown in Figure 7a, when the cohesion c is 2 MPa, the range of the excavation damage zone is approximately 2.3 to 4.3 times the radius of excavation. In this instance, an increase of c results in a rapid decrease to Rp/R0. It is indicated that the stability of the roadway can be significantly improved with the increase of cohesion. Nevertheless, once c exceeded 6.0 MPa, the increase of c will no longer have an important effect on Rp/R0. Thus, it is concluded that a minimal damaged zone in surrounding rocks masses with high cohesion c will be induced.
It can be seen in Figure 7b that Rp/R0, which is from 1.9 to 2.78 when the friction angle was 10°, decreases exponentially with the increasing friction angle φ of the rocks. Moreover, the extent of EDZ increases with the mining depth. Occasionally, when φ > 30°, the effect of the friction angle on Rp/R0 diminished, although a continuous increase in φ further reduced the damaged zone in the surrounding rock. A comparison of the effects of c and φ on the integrity of the surrounding rock showed that the effect of c was evidently greater than that of φ. In underground mining, the instability of rock masses around excavations caused by disturbances of the cohesion and internal friction angle should be minimized. Figure 7c shows that Rp/R0 is greatly reduced with an increasing geological strength index GSI value, which is similar to each of the other four field conditions. The ratio of Rp to R0 approaches 5.15 when GSI = 10. In this case, an increase in GSI rapidly decreased Rp/R0, which means that the stability of the surrounding rocks improved. However, once GSI increases to 40, the effect on Rp/R0 was no longer remarkable. Therefore, several techniques, such as deep grouting, spray anchor supporting, and bolt grouting provide remarkable effects on the rock mass integrity and tunnel stability, especially in inadequately jointed rock masses.
The influence of mi on the extent of the damage zone is illustrated in Figure 7d. Here, there is a negative correlation between the extent of EDZ and mi of the rock mass. The effect of mi, which is compared with the influence of GSI on the extent of damage, is quite small. On the other hand, when mi increases to 25, the Rp/R0 ratios of the five types of cases did not approach 1, and there is a certain 8   As shown in Figure 7a, when the cohesion c is 2 MPa, the range of the excavation damage zone is approximately 2.3 to 4.3 times the radius of excavation. In this instance, an increase of c results in a rapid decrease to R p /R 0 . It is indicated that the stability of the roadway can be significantly improved with the increase of cohesion. Nevertheless, once c exceeded 6.0 MPa, the increase of c will no longer have an important effect on R p /R 0 . Thus, it is concluded that a minimal damaged zone in surrounding rocks masses with high cohesion c will be induced.
It can be seen in Figure 7b that R p /R 0 , which is from 1.9 to 2.78 when the friction angle was 10 • , decreases exponentially with the increasing friction angle ϕ of the rocks. Moreover, the extent of EDZ increases with the mining depth. Occasionally, when ϕ > 30 • , the effect of the friction angle on R p /R 0 diminished, although a continuous increase in ϕ further reduced the damaged zone in the surrounding rock. A comparison of the effects of c and ϕ on the integrity of the surrounding rock showed that the effect of c was evidently greater than that of ϕ. In underground mining, the instability of rock masses around excavations caused by disturbances of the cohesion and internal friction angle should be minimized. Figure 7c shows that R p /R 0 is greatly reduced with an increasing geological strength index GSI value, which is similar to each of the other four field conditions. The ratio of R p to R 0 approaches 5.15 when GSI = 10. In this case, an increase in GSI rapidly decreased R p /R 0 , which means that the stability of the surrounding rocks improved. However, once GSI increases to 40, the effect on R p /R 0 was no longer remarkable. Therefore, several techniques, such as deep grouting, spray anchor supporting, and bolt grouting provide remarkable effects on the rock mass integrity and tunnel stability, especially in inadequately jointed rock masses.
The influence of m i on the extent of the damage zone is illustrated in Figure 7d. Here, there is a negative correlation between the extent of EDZ and m i of the rock mass. The effect of m i , which is compared with the influence of GSI on the extent of damage, is quite small. On the other hand, when m i increases to 25, the R p /R 0 ratios of the five types of cases did not approach 1, and there is a certain amount of deviation. This result may be attributed to the fact that as the depth of embedment increased, and fractures are inevitably caused in the surrounding rock masses around tunnels under high in situ stress.
The effect of the disturbance factor D on R p /R 0 is depicted in Figure 7e where a curvilinear trend is seen. The slope increases with the intensity of the dynamic disturbance of rock masses. When D > 0.6, the influence on the EDZ becomes even more evident, especially for deeply embedded rock masses. These results demonstrate that the dynamic disturbance had an evident weakening effect on the stability and integrity of rock masses under high in situ stress. Hence, for practical mining engineering, the disturbance factor, D, should not be larger than 0.6. Otherwise, the thickness of the EDZ will be considerably affected. In the construction of underground projects, it is necessary to include measures to decrease D, such as smooth and presplit blasting, decoupled charging, and delayed detonation. Figure 7f shows that the range of EDZ increases with the increase of σ c . When the uniaxial compressive strength of the rock is 70 MPa, the ratio of the damaged radius to tunnel radius is approximately 3.14, and the variation trend of R p /R 0 at 780 m is similar to those of the other four types of rock masses. When the uniaxial compressive strength of the rock increase and reaches to 115 MPa, the effect of σ c on R p /R 0 diminishes. Thus, in the construction of openings in rock masses with σ c less than 115 MPa, reinforcing supports for the excavation should be considered.
In conclusion, the stability and integrity of a rock mass are positively correlated with c, ϕ, GSI, m i , and σ c . On the contrary, the parameter D varied from 0 to 1 has a weakening effect on the stability of tunnels. In other words, thanks to the action of dynamic disturbance, an extensive EDZ will be induced in rock mass subjected to high in situ stresses. It is noteworthy that ϕ, GSI, and D had more significant effects on the extent of EDZ than the other three characteristics. Thus, in tunnel construction, measures to reduce the extent of EDZs in rock masses with varying physical and mechanical characteristics should be carefully considered. In addition, some measures should be adopted to control the expansion of damaged zones, such as pre-grouting, bolt installation coupled with cables, decoupled charging, and so forth.

Conclusions
During the excavation of underground tunnels in mining, hydropower and radioactive waste disposal operations, an unexpected excavation damaged zone is induced around the tunnels. In this paper, closed-form solutions based on the Mohr-Coulomb and Hoek-Brown failure criteria are derived for predicting the extent of EDZ in the surrounding rock mass around a circular opening subjected to a non-hydrostatic stress field. In addition, the distribution law of in situ stress field in the Fankou lead-zinc mine was analyzed by the overcoring method (stress relief method) accompanying the Swedish high-precision LUT system. The results indicate that σ h-max , σ h-min , and σ v of the mining area increase approximately linearly with the depth, and the in situ stress field is dominated by maximum horizontal tectonic stress, followed by vertical principal stress. An evaluation of EDZ is carried out by using a SET-PLT-01 nonmetallic ultrasonic detector. It is shown that the thickness of EDZ in the surrounding rock mass shows a linearly increasing trend with depth, and they are greater on the roof-floor region than those on the sidewall. A case study is presented to verify the accuracy of the proposed closed-form solutions in terms of field measurements. The results show that the thicknesses calculated based on Hoek-Brown criterion were 6.6 to 11.7% larger than those based on the Mohr-Coulomb criterion; thus, it reveals that Hoek-Brown criterion achieved a better performance than the Mohr-Coulomb criterion in EDZ evaluation. The damage thickness assessed with the single-hole acoustic wave monitoring method was approximately 12 to 21% larger than that calculated under the two derived analytical solutions. The results of the parametric studies illustrate that parameters c, ϕ, GSI, and σ c have a significant effect on EDZ. Inevitable dynamic disturbance D shows a more critical influence on the stability of underground rock mass. Therefore, the study on evaluating the damage zone around underground openings by theoretical calculation and field test methods will be a great help in understanding the evolution law of in situ stresses and EDZ with the depth. The work could also find its application in the field to predict the extent of the excavation-damaged zone under the premise that the physical and mechanical properties of rock mass have been obtained.