Abstract

The traditional rock burst estimation method is usually based on the (rock strength) in practice, while the actual occurrence of rock burst depends more on the structure and strength of the rock mass. In this work, the actual rock bursts occurred in a railway tunnel project in Western China, and the (rock mass strength) was calculated by the generalized Hoek–Brown criterion. According to the actual situation of rock bursts, a modified rock burst estimation criterion using the ratio of σcm to (maximum geostress) was proposed. The influence of randomness on the reliability of rock burst estimation criterion was considered. The estimation results based on the traditional and modified method were furthermore compared with those of the actual rock bursts. The results show that calculated by the generalized Hoek–Brown criterion may be considered well in the rock type and strength, construction condition, and structure features of the rock mass; the estimation results of rock burst using the ratio of to are quite different from the actual situation, while those using the ratio of the to coincided relatively with the actual rock bursts; the ratios of to , which are greater than 0.167, 0.066 to 0.167, 0.012 to 0.066, and less than 0.012, are corresponded to the slight, medium, strong, and violent grades of the rock bursts, respectively; the randomness of data selection has certain influence on the rock burst estimation criterion, but the variation range is small; the modified estimation criterion of rock burst proposed in this work has a good reliability. The results presented herein are important for tunnel construction and the prevention of rock burst in the high geostress areas.

1. Introduction

When excavating in ground prone to rock burst, it is essential to estimate the intensity of rock burst in order to implement suitable ground control measures. A rock burst is a visible rock mass damage to excavations and represents a major concern in tunneling projects since it is associated with casualties and accidents. This deadly phenomenon initiates in high geostress and brittle rock conditions with extremely complex mechanical features. This may lead to a violent expulsion of rock from the surrounding rock mass. Usually, the rock burst occurs with a sudden release of elastic strain energy in different ways including slabbing, spalling, ejecting, or throwing [1, 2].

At present, rock burst is usually estimated based on the (rock strength) in practice, and there are many achievements that have been made in this field. For example, Tao [3] proposed a rock burst criterion based on the ratio of to (maximum geostress), while Gu et al. [4] proposed a comprehensive criterion for the occurrence of rock burst based on σc and rock integrity. Later on, Zhang and Fu [5] modified the Tao–Gu criterion. Liu et al. [6] introduced a membership function to investigate the influence of on rock burst occurrence. Miao et al. [7] used the Tao discriminant criterion and combined with other rock parameters to comprehensively estimate the rock burst grades. Zhou et al. [8], Li et al. [9], and Afraei et al. [10] regarded the as one of the main factors to estimate rock burst intensity. Xue et al. [11] adopted the ratio of to geostress as an indicator to estimate the rock burst.

However, the actual occurrence of rock burst depends not only on the , but also on the structure of the rock mass and the . Rock mass structure, which may change the evolution mode of rock burst activity, is the important factor causing rock burst [12, 13] There are also some reports about the influence of rock mass structure on rock burst. For example, Zhou et al. [14] found that the rock mass structure plane plays an important role in controlling rock burst and analyzed the mechanism of rock burst caused by the structure plane. Mohamad et al. [15] found that when rock burst occurred, a large number of flying rocks are produced in the rock mass plane with small joint spacing and large aperture. Du et al. [16] recognized that the weakening feature of rock mass structure surface is the main factor to estimate whether rock burst is instantaneous or lagged.

At present, the Hoek–Brown (H-B) strength criterion for calculating based on rock mass structure is widely recognized. Hoek and Brown [17] gave a nonlinear empirical relationship between the rock mass type, rock mass quality criterion, and . Hoek and Brown [18] obtained the relationship between the ultimate principal stress of rock mass and the rock mass failure through statistical analysis of rock triaxial test data and the rock mass tests. Sharifzadeh et al. [19] calculated the shear strength of rock mass using Hoek–Brown failure criterion and geological strength criterion (GSI). Wu et al. [20] proposed a method to analyze the GSI and the disturbance factor of the rock mass and calculated the RMS based on this method. Ma et al. [21] thereafter estimated the rock burst by a ratio of the H-B strength-based to the horizontal stress perpendicular to the tunnel axis.

The traditional rock burst estimation method usually estimates the rock burst grades based on the ratio of to , and the estimated results are often different from the actual situation. Nevertheless, the actual rock burst depends not only on the geostress state and the strength of rock, but also on the rock mass structure and the construction disturbance. In this work, the geostress state was determined by the hollow inclusion stress relief method. The σc was obtained by the laboratory uniaxial compressive test. The rock mass structure was characterized by GSI according to the actual situation. The σcm was calculated by the generalized H-B strength criterion. The 18 measurement points were randomly selected from the 22 field measurement points to compare the actual rock burst situation. According to the actual construction situation and rock burst situation, a modified rock burst estimation criterion based on σmax and σcm was thereafter investigated. The reliability of the modified rock burst estimation criterion was verified by the remaining 4 measurement points, and the influence of randomness of data selection on the stability of the modified rock burst estimation criterion was considered. The research results have important reference value for tunnel construction and determination of rock burst prevention measures.

2. Measurement of Geostress

2.1. Engineering Background

The geostress measurements located at the construction site of a railway tunnel in Western China. The measurement area has extremely harsh climate where the mountain ranges extend longitudinally in the north-south direction and the terrain fluctuates greatly. There are many regional great and active faults with frequent earthquakes, high tectonic stresses, and well-developed joints. Due to the fact that the surrounding rock mass is hard and brittle, it is expected that rock burst may occur. The 22 geostress measurements are all taken from the railway tunnel.

2.2. Process of Geostress Measurement

The geostress data were obtained by using the hollow inclusion stress relief method. As shown in Figure 1, a great hole inclined by 3° to 5° with a diameter of 130 mm was drilled during the test. A concentric small hole with a diameter of 36 mm was then drilled from the bottom of the great hole.

As shown in Figure 2, after completing the small hole, a dry towel was used to wipe the hole and to arrange the adhesion. After mixing the binder, the butter was smeared evenly on the orientation instrument, the surface of the hollow inclusion stress gauge was grinded with sandpaper, the binder was poured into the inner cavity of the stress meter, and the stress gauge was propelled into the borehole with the mounting rod. The entire installation process was controlled within 20 minutes to maintain the fluidity of the colloid.

After the colloid was completely solidified, the inclination and orientation of the borehole, the position of the orifice (earth coordinates), and the mounting angle were measured. As shown in Figure 3, when the stress gauge was released, the data were collected by the YJZ-16+ intelligent digital strain gauge, and the core with the strain gauge is removed. The data are received by the KJ327-F type of the mine pressure monitoring system, and the core with the stress gauge was placed into the confining calibration instrument.

As shown in Figure 4, after gradually applying the confining pressure on the core, the stress-strain curve was obtained and the elastic modulus and Poisson’s ratio were thereafter computed.

In this work, the geostress data from the 22 in situ measurement points were used in the later analysis, where 18 of them are randomly selected to compare the actual rock burst situation and estimation results by using the traditional method. The improvement of the traditional estimation method was conducted by analyzing the actual rock burst grades. The correctness of the modified method will be verified by the remaining 4 points. The field results of the σmax of these 18 measurement points and the corresponding laboratory results from the uniaxial compressive tests are shown in Table 1.

3. Estimation of Rock Burst Tendency

3.1. Actual Situation of Rock Burst

During the tunnel construction, the appearances of the 18 measurement points are shown in Figure 5.

Five parameters, i.e., motion, sound, aging, impact on construction, and influence depth, were used to comprehensively classify the rock burst grades by the code for hydropower engineering geological investigation (CHEGI), suggested by the National Standards Compilation Group of People’s Republic of China [22]. The measurement point 1 was taken as an example to concisely describe the details in estimating rock burst grade. At the measurement point 1, many large rocks flew out rapidly, accompanied by the rock powder ejection; a strong burst sound was heard; the rock burst lasted for a long time; there is a great influence on the construction of the tunnel; the rock burst pit is distributed continuously with the influencing distances of more than 2 m. Therefore, the rock burst at this point was determined at the strong level. Based on the actual situations of rock bursts, the rock burst grades at 18 measurement points were determined by the above method and are shown in Table 2.

3.2. Rock Burst Estimation Results Using Traditional Criterion

The CHEGI criterion based on the ratio of σc to σmax is often used to estimate the rock burst in practice. The corresponding criterion is shown in Table 3.

In this work, the σcm, σc, and rock burst grades were obtained by using the hollow inclusion stress relief method, laboratory uniaxial compressive tests, and the ratio of σc to σmax, respectively. The results are shown in Table 4. For comparison purposes, the observation results of field rock bursts are also listed in Table 4. From Table 4, it can be seen that results by using the CHEGI rock burst estimation criterion are quite different from those at the field situations.

3.3. Modified Criterion of Rock Burst Estimation

In the traditional criteria for rock burst classification, the strength in the strength-stress ratio method generally refers to the uniaxial compressive strength. However, the actual occurrence of rock burst in practice relates much with the rock mass structure and the strength of the corresponding rock mass. Therefore, it is necessary to improve the traditional method to pay attention both to the rock mass structure features and to the σcm in the rock burst estimation.

The rock mass structure features are generally characterized by the geological strength criterion, or GSI. Based on GSI and other relevant parameters (such as the disturbance coefficient and uniaxial compressive strength), the generalized H-B strength criterion is often used to calculate the σcm, and the ratio of σcm to σc may be used as estimation criterion of the rock burst. Therefore, the rock burst estimation based on σcm will be proposed in the following.

3.3.1. Determination of GSI

The geological strength criterion, or GSI, introduced by Marinos and Hoek [23] and considering the structural features, weathering condition, and the surface features of the rock mass, could better reflect the geological situation of rock mass. As shown in Figure 6, based on the discontinuity, structure, and surface condition of rock mass, the average value of GSI may be estimated. In this figure, “N/A” means that it is not applicable within this range.

As shown in Figure 6, the surface quality of rock mass structure may be divided into five categories using the weathering condition of rock mass and the surface features of joints, which are Very Good, Good, Fair, Poor, and Very Poor. Among them, Very Good means very rough, fresh, and unweathered surface of rock mass; Good means rough, mild weathered, iron surface of rock mass; Fair means medium weathered and altered surface of rock mass; Poor means smooth, highly weathered rock mass surface with a dense overburden or filler or angular fragments; Very Poor means smooth, severely weathered rock mass surface with a soft clay coating or filler. The corresponding range of values is from 100 to 0 in order of high to low, and the higher the value, the better the quality grade of rock mass surface. The structural features of the rock mass are divided into six categories using the order of the integrity of the rock mass structure surface, which are Intact or Massive, Blocky, Very Blocky, Blocky/Disturbed/Seamy, Disintegrated, and Laminated or Sheared. Wherein the Intact or Massive means a complete rock mass or a large rock mass structure with few large spacing and discontinuity; the Blocky means a good and original rock mass structure composed of cubic blocks formed by three mutually orthogonal joint faces; Very Blocky means a partially disturbed rock mass structure, which composed of multifaceted angular blocks formed by at least 4 sets of joints; Blocky/Disturbed/Seamy means a rock mass structure in which a plurality of sets of discontinuous surfaces are mutually cut to form an angular rock mass, and undergoes a fold activity, a layer, or a flank plane continuous; Disintegrated means a severely fractured rock mass, which contains a mixture of angular and circular rock blocks; Laminated/Sheared means the lack of a massive rock mass structure due to the weak schistosity or the close spacing of the shear planes. The corresponding range of values is from 100 to 0 in order of high to low, and the higher the value, the better the integrity of the rock mass.

The Marinos–Hoek method in estimating GSI was used in this work. The measurement point 1 was taken as an example to concisely describe the process in estimating GSI. In the measurement point 1, the structural fissures and unloading fissures are well developed; the structural fissures are often distributed in an “X” shape; the joints generally do not extend; the fissures are 1 to 10 mm wide with a maximum width of 15 mm; the fissures are mostly half-open and filled visibly; the joint spacing is more than 2 m. Accordingly, the surface condition of the rock mass at this point was determined as “Good,” and the value was estimated as 62; the rock mass structure was determined as a blocky structure with a value of 71.

As shown in Figure 7, when the intersection of the vertical line of the rock surface condition and the horizontal line of the rock structure characteristic was between the two GSI values, linear interpolation was often used for calculation. Accordingly, the GSI of the measurement point 1 was estimated as 58. According to the above method, the GSI values of the remaining measurement points were estimated, and the estimated GSIs at 18 measurement points are listed in Table 5.

3.3.2. Determination of Relevant Parameters

Other relevant parameters mainly relate to the rock mass disturbance and σc. The related parameter of the rock mass disturbance degree is coefficient D, which represents the disturbance of the rock mass. Considering the actual situation of the tunnel excavation, D’s of all the measurement points are assumed to be 0.5. Rock uniaxial compressive strength was determined: the σc at each measurement point is determined by the laboratory uniaxial compression test (see early-mentioned Table 1).

3.3.3. Calculation of Rock Mass Strength Using Hoek–Brown Criterion

Based on Griffith’s theory, Hoek et al. [24] investigated the relations between the ultimate principal stress in a rock mass and the rock mass failure from the statistical results of the rock triaxial tests and the rock mass tests and obtained the generalized H-B strength guidelines.

In the generalized H-B strength criterion, the σcm is calculated using the following formula:where σc is the uniaxial compressive strength of intact rock; s is the empirical parameter reflecting the fracture degree of the rock mass, ranging from 0 to 1, respectively, taking 0 for the fully fractured rock mass and 1 for the intact rock mass; a is the empirical parameter that reflects the features of the rock mass; GSI is the geological strength criterion of the rock mass; D is a parameter reflecting the disturbance degree of the field rock mass influencing by external factors and ranges from 0 to 1, respectively, taking 0 for the undisturbed rock mass and 1 for the completely disturbed rock mass; σcm is the uniaxial compressive strength of the rock mass.

At measurement point 1, substituting GSI = 58 and D = 0.5 into equations (1) and (2), s and a will be 0.003698 and 0.5033, respectively. Substituting s, a, and σc into equation (3), the σcm is 2.42 MPa. The calculation results of the rock mass strength σcm of 18 measurement points are shown in Table 6.

3.3.4. Rock Burst Estimation Based on Rock Mass Strength

In computing the ratio of the rock mass strength σcm to the geostress, Ma et al. [21] represented the geostress as the horizontal stress perpendicular to the tunnel axis. In the current study, the maximum stress σmax was still used to represent the geostress, considering the availability in practice and the comparability with the existing specifications (i.e., the code for hydropower engineering geological investigation, 2016).

The ratios of σcm to σmax at different measurement points are listed in Table 7 and shown in Figure 5. The on-site field rock burst grades from the actual observations are also shown in Table 7.

As shown in Table 7, among the 18 randomly chosen points, the number of measurement points for slight rock burst, medium rock burst, strong rock burst, and violent rock burst is 1, 11, 6, and 0, respectively. Because the point number of violent rock burst is zero, it is difficult to accurately determine the extent of the ratio for the violent rock burst. In order to facilitate the division of the rock burst, the ratio of σcm to σmax is taken as 0 to be temporarily used to estimate the violent rock burst.

From Figure 8, it can be seen that the rock burst is lighter when the ratio of σcm to σmax is greater, and vice versa. Therefore, three boundary lines may be used to divide the ratio of these measurement points into four parts, corresponding to the grades of the slight, medium, strong, and violent rock bursts. The determination processes of these three boundary lines are as follows:(a)The determination of the boundary line between the slight and medium rock bursts. In this case, one and 12 measurement points belong, respectively, to the slight and medium rock bursts. Because  = 0.1850 and  = 0.1497 represent the minimum and maximum ratios of the slight and medium rock bursts, respectively, at the measurement points 3 and 16, the midvalue of and , or y = 0.16735, is used as the boundary ratio between these two grades of rock bursts.(b)The determination of the boundary line between the medium and strong rock bursts. Because  = 0.0738 and  = 0.0573 represent the minimum and maximum ratios of the medium and strong rock bursts, respectively, at the measurement points 8 and 12, the midvalue of and , or y = 0.06555, is used as the boundary ratio between the grades of the rock bursts.(c)The determination of the boundary line between the strong and violent rock bursts. In this case,  = 0.0236 and y = 0 are taken as the minimum and maximum ratios for the strong and violent rock bursts, respectively, at the measurement point 6 and others. The midvalue of and 0, or y = 0.01180, is regarded as the boundary ratio between the strong and violent rock bursts.

To consider the facility in applications, it is better to transfer these threshold values into ones with three digits after decimal points. Accordingly, the ratio intervals are more than 0.167, , , and no more than 0.012, respectively, representing the slight, medium, strong, and violent rock bursts. A rock burst estimation method (see Table 8) based on the RMS-to-MG ratio is thereafter obtained.

3.3.5. Verification of Estimation Criterion for Rock Burst

In order to address the reliability of the abovementioned estimation criterion for rock bursts, the remaining 4 field measurement points are used for verification. The appearances after the rock bursts at these points are shown in Figure 9.

Using the traditional method based on the ratio of σc to σmax, the estimation results at these four verification points are obtained (see Table 9).

Using the modified criterion, the ratios of σcm to σmax at these four verification points were obtained (see Table 10).

Comparing the index in Table 9 and the index in Table 10 with the rock burst grade estimation criteria in Table 3 and Table 8, respectively, the rock burst estimation results of the CHEGI’s criterion and modified criterion at four verification points were obtained (see Table 11). Table 11 also lists the field observations.

As shown in Table 11, the estimation results using the CHEGI criterion based on the ratio of σc to σmax are quite different from those in the actual situations; if the criterion based on the ratio of σcm to σmax is used, the estimation results for the rock bursts will be consistent with the field observations. Accordingly, the rock burst estimation only considering the rock strength regardless of the rock mass structure is not suitable. In the modified criterion, much attention is paid on the rock mass structure in computing the ratio of σcm to σmax. The estimation results based on the modified criterion are consistent with the actual situation and may be used for rock burst estimation.

3.4. The Influence of Randomness of Data Selection on Rock Burst Estimation Criterion

Using the results of 22 field measurement points, 18 of them were randomly selected to match the actual situation of rock bursts, and a new rock burst criterion was proposed. In order to investigate the influence of the randomness of data selection on the estimation criterion of rock burst, 12 sets of data, each of which includes 18 measurement points, were randomly selected from 22 field measurement points. Using the above method, the corresponding rock burst estimation criteria were obtained and compared. The results are shown in Table 12.

It can be seen from Table 12 that among the estimation criteria of slight rock burst grade in 12 sets, 11 of them are more than 0.167, and the other is greater than 0.156, with little change in range; the estimation criterion of medium rock burst grade is 0.066 to 0.167 in 8 sets, 0.062 to 0.167 in 3 sets, and 0.066 to 0.156 in 1 set; the change of estimation criterion of strong rock burst grade is similar to that of medium rock burst grade, with 8 sets ranging from 0.012 to 0.066, 3 sets ranging from 0.012 to 0.062, and 1 set ranging from 0.019 to 0.066; the estimation criterion of violent rock burst grade is no more than 0.012 in 9 sets and less than 0.019 in 3 sets; the range of change is a little big, but the probability of occurrence of a situation less than 0.019 is small. Accordingly, the criteria of slight, medium, strong, and violent rock bursts were determined to be greater than 0.167, 0.066 to 0.165, 0.012 to 0.066, and less than 0.012, respectively, which have a good stability.

Comparing the data in Tables 12 and 8, it is found that the rock burst estimation criterion obtained in this work has a relatively high stability. It can be seen from Tables 7, 10, and 11 that there is only one slight rock burst measurement point in the 22 field measurement points. If the slight rock burst measurement point was not selected, the maximum ratio of RMS to MG in medium rock burst measurement points will have to be taken as the boundary value between slight and medium rock bursts, which will have a certain impact on these two rock burst estimation criteria. However, the probability of not selecting this slight measurement point from 22 points is 18.18%, so the probability value is relatively small. Moreover, measurement points of the medium rock burst were relatively large, which will weaken this impact to some extent. In the later studies, the accuracy of the boundary value between the slight and medium rock bursts may be modified by increasing the number of slight rock burst measurement points.

4. Discussions

(1)Considering the rock mass structure could change the evolution mode of rock burst activity and the rock mass (GSI) degradation could be used for avoiding the risk of rock burst [25, 26], the effect of GSI on the estimation of rock burst grades was explored. The verification point 3 was taken as an example (see Figure 10), where various GSI values (50, 51, 52, 53, 54, 55, 56, 57, 58, 59, and 60) were selected.From Figure 10, it can be seen that the estimation index of rock burst grade increases approximately linearly with the increase in GSI, and (goodness of fit) is 0.9916; the effect of GSI on the estimation of rock burst grades is relatively great, and increases by 7.1% when GSI increases by 1; the rock burst grade of verification point 3 is strong if GSI is between 50 and 53, while the rock burst grade is medium if GSI is between 54 and 60. Therefore, the effect of GSI on the rock burst grade needs to be paid more attention.(2)Considering the Hoek–Brown criterion may overestimate the strength of rock mass [27, 28], various methods for estimating rock mass strength were conducted. Taking the measurement point 1 as an example, the values of (the weight of rock), Q (rock mass quality rating), and RMR (rock mass rating) are 26 kN/m3, 0.12, and 25, respectively. The estimation results of rock mass strength are shown in Table 13.From Table 13, it can be seen that the rock mass strength is mainly distributed in the range of 2.4 to 9.0 MPa and 5 to 20 times less than the rock mass strength; the estimations using Hoek’s and Kalamaras’s methods are relatively close, while those using Barton’s and Singh’s methods seem to be too high. Hence, the rock mass strength estimated by Hoek’s method is relatively reasonable.(3)In order to explore the reliability of the modified method in this paper, the rock burst grade estimation results were compared using various criteria (see Table 14). As can be seen from Table 14, the classification of Barton’s criterion, which simply classified the rock burst into the mild and heavy rock burst grades, was not accurate enough; even though Rehman et al. [32] considered the rock is jointed and proposed an index of SRFQ, the classification of rock burst was also inaccurate; based on Barton’s criterion, Tao [3] divided the mild rock burst into the medium and slight ones, while there is no subdivisions in heavy rock bursts; the rock burst grades were divided into four grades (slight, medium, strong, and violent) by the CHEGI’s criterion, but this criterion ignored the effect of the rock mass structure on the rock burst. In the current study, the rock burst grades are divided into four grades using the modified criterion where the rock mass structure was well considered.The estimation results at four verification points are shown in Table 15. As can be seen from Table 15, the rock burst grades of verification points 1 and 4 estimated by various criteria are relatively consistent, while the rock burst grades at verification points 2 and 3 estimated by Barton’s, Rehman’s and Tao’s criteria are all greater than those estimated by the modified criterion; for these four verification points, the rock burst grades estimated by the CHEGI’s criterion are greater (in an one grade) than those estimated by the modified criterion.Estimated results using various criteria with the actual rock burst were compared and are shown in Table 16.From Table 16, it can be seen that the estimation results using Barton’s, Rehman’s and Tao’s criteria are overestimated and have a lower accuracy; compared with the actual rock burst results, the estimated ones using the CHEGI’s criterion are also overestimated; the estimated rock burst grades using the modified criterion in the current study are in good agreement with the actual ones.(4)When determining the boundary value of adjacent rock burst grades, it is not enough to consider only the maximum or minimum value of the rock burst grades as the boundary value, and the reliability of the rock burst estimation index obtained by this method is greatly influenced by the factors of artificial selection of data. In this work, the method of taking the midvalue may well consider the data of two adjacent rock bursts, which makes the estimation criterion of rock burst closer to the real value and has better generalization ability. In the future, the boundary value of rock burst grades may be further optimized by increasing the number and diversity of rock burst data.(5)Taking the Neelum-Jhelum Hydroelectric Project in Pakistan as an example ([12] and [33]), the estimation results using various criteria were compared to further analyze the reliability of the modified criterion. The results are shown in Table 17.

It can be seen from Table 17 that the rock burst estimation grades using traditional criteria are quite different from the actual ones; the estimation results using the modified criterion in this study are much closer to the actual rock burst; the estimation results using the modified criterion have a good reliability for the measurement points with GSI greater than 60.

5. Conclusions

In this work, after using the generalized H-B strength criterion to calculate the σcm, a modified rock burst estimation criterion based on the ratio of σcm to σmax was proposed. The estimation results of this method were compared with those of the traditional method based on the ratio of σc to σmax. It shows that(1)The generalized H-B criterion in calculating σcm may reflect both the structural features of the rock masses and actual situations of the excavation disturbances(2)The estimation results by using the ratio of σc to σmax in the traditional method are quite different from the actual situation, while the estimation results by using the ratio of σcm to σmax in this work are consistent with the actual situation(3)If a ratio of σcm to σc is used for the rock burst estimation, the ratio intervals in the slight, medium, strong, and violent rock burst grades are more than 0.167, , , and no more than 0.012, respectively.(4)The randomness of the selection of rock burst measurement points has a certain influence on the determination of the estimation criterion of rock burst grades, but the overall change of the range is not large. The rock burst estimation criterion proposed in this work has a good reliability.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

Financial supports for the study were provided by the Natural Sciences Foundation Committee of China under Grant no. 41472254 and the Science and Technology Research and Development Program of China Railway Construction Corporation Limited under Grant no. 17-C13. These supports are gratefully acknowledged.