A prediction model on rockburst intensity grade based on variable weight and matter-element extension

Abstract

Rockburst is a common dynamic disaster in deep underground engineering. To accurately predict rockburst intensity grade, this study proposes a novel rockburst prediction model based on variable weight and matter-element extension theory. In the proposed model, variable weight theory is used to optimize the weights of prediction indexes. Matter-element extension theory and grade variable method are used to calculate the grade variable interval corresponding to the classification standard of rockburst intensity grade. The rockburst intensity grade of Engineering Rock Mass is predicted by rock burst intensity level variable and the interval. Finally, the model is tested by predicting the rockburst intensity grades of worldwide several projects. The prediction results are compared with the actual rockburst intensity grades and the prediction results of other models. The results indicate that, after using variable weight theory and grade variable method, the correct rate of prediction results of matter-element extension model is improved, and the safety of the prediction results is also enhanced. This study provides a new way to predict rock burst in underground engineering.

About the Authors:

Jianhong Chen

Contributed equally to this work with: Jianhong Chen, Yi Chen

Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Supervision

Affiliation: School of Resource and Safety Engineering, Central South University, Changsha, Hunan, China

Yi Chen

Contributed equally to this work with: Jianhong Chen, Yi Chen

Roles Investigation, Methodology, Validation, Visualization, Writing – original draft

Affiliation: School of Resource and Safety Engineering, Central South University, Changsha, Hunan, China

Shan Yang

Roles Funding acquisition, Supervision, Writing – review & editing

* E-mail: [email protected]

Affiliation: School of Resource and Safety Engineering, Central South University, Changsha, Hunan, China

ORCID logo http://orcid.org/0000-0002-8249-0183

Xudong Zhong

Roles Investigation, Methodology, Software, Validation

¶‡ These authors also contributed equally to this work.

Affiliation: School of Resource and Safety Engineering, Central South University, Changsha, Hunan, China

Xu Han

Roles Investigation, Validation, Visualization

¶‡ These authors also contributed equally to this work.

Affiliation: School of Resource and Safety Engineering, Central South University, Changsha, Hunan, China

Introduction

Rockburst is the common energy release of elastic deformation with sudden, uncontrollable and destructive damage in deep mining and underground space development, which is manifested as the fracture and instantaneous projection of surrounding rock mass [1–3]. Strong rockbursts may cause significant losses to underground engineering, such as the failure of the supporting mechanism and equipment, casualty and construction delay [4,5]. Therefore, it is significant to predict and alarm the rockburst in underground engineering [6,7]. The research on the mechanism of rockburst has been studied from different perspectives with scientific rockburst criteria and prediction method [8]. Xie H P [9] and Cai M F [10,11] obtained the rockburst criteria based on the energy principle by studying the internal relationship between energy accumulation, dissipation, release, rock strength and overall failure in the process of rock deformation and failure. Arteiro A [12] summarized the critical energy value of rockburst through triaxial test and used it as the criteria of rockburst prediction. Jiang Q [13]et al. proposed that local energy release rate (LERR) is a new energy index to simulate rockburst conditions, contributing to understanding the rockburst mechanism from the perspective of energy. He M C [14,15], Kuksenko V S [16] et al. obtained the rockburst prediction criteria by analyzing the characteristics of acoustic emission signal in the rock failure. Feng X T [17], Ma T H [18], Liu X H [19], Chen Z H [20] et al. used microseismic monitoring technology to explore the relationship between the spatiotemporal evolution of microseismic signals and rockburst, and obtained the criteria for rockburst prediction and early warning based on the features of microseismic events. Sharan S K [21] proposed a finite element model for predicting the potential occurrence of rockburst in underground caverns. Mole-Coulomb failure criteria[22] and Hoek–Brown failure criteria[23] was used as rockburst criteria. These studies on rockburst criteria enrich the understanding of rockburst mechanism. However, due to the complex rock mass structure and geological environment in the underground engineering, rock burst prediction by a single rock burst criterion is usually less accurate.

To improve the correct rate of prediction, the information of rockburst reflected by various rockburst criteria have been synthesized by the relevant mathematical theory, and many rockburst prediction and early warning models have been established. For example, Zhou J [24,25] used Fisher discriminant theory to construct a Fisher analysis discriminant model for rockburst classification prediction, which can accurately predict rockburst in deep tunnels and some coal mines. Guo L [26] used the theory of neural network to construct the BP neural network prediction model of rockburst prediction. A large amount of rockburst data were required for the model training, otherwise the correct rate of prediction and the application scope were limited. Liu Z J [27] used the fuzzy probability theory to establish a new fuzzy probability model for predicting the occurrence and intensity of rockburst. The fuzzy weight was introduced in this model, and the limitation of traditional fuzzy probability model was overcame to a certain extent in practical engineering. Gao W [28] constructed a rockburst prediction model by ant colony clustering algorithm. Engineering analogy was performed to realize the automatic classification of rockburst in this model, while prediction accuracy still needed to be further improved. Pei Q T [29] considered the complex relationship between rockburst and its influencing factors as a grey system, optimized the grey whitening weight function, and established an improved grey evaluation model for rockburst prediction. The model greatly overcame the problems of multiple intersections of grey types and abnormality, and improved the prediction accuracy. Based on Projection Pursuit (PP), Particle Swarm Optimization (PSO) and Logistic Curve Function (LCF), a particle swarm projection pursuit model for rockburst prediction was constructed by Xu F [30], ZHOU Xuanchi [31] et al. Particle swarm optimization was used to optimize the parameters of projection index function and LCF, guaranteeing the accuracy of model parameters and prediction accuracy. However, with the increasing number of prediction indexes, the optimization of model parameters becomes more difficult. Zhou J [32–34] used a large number of rockburst data to compare the learning ability of 11 supervised learning algorithms, including linear discriminant analysis (LDA), partial least-squares discriminant analysis (PLSDA), quadratic discriminant analysis (QDA), multilayer perceptron neural network (MLPNN), support vector machine (SVM), random forest (RF), gradient-boosting machine (GBM) and so on. Limitations of these algorithms were also analyzed by Zhou J. In addition, the commonly used models for rockburst prediction include evidence theory model [35], set pair analysis theory model [36], efficacy coefficient method model [37], normal cloud model [38]. Although these models have been applied to the engineering, there are still some defects in the models, such as the difficult fusion of conflicting indicators in evidence chain model; the randomness of forecasting cannot be reflected by the failure of set pair analysis model; the difficult obtaining of satisfactory and unsatisfactory values in efficacy coefficient method; and the gap between forecasting results and actual ones caused by the over ideal distribution of indicators in normal cloud model. Therefore, it is necessary to improve the existing rockburst prediction theories and models.

Compared with the extension prediction model in the previous study, this paper proposes a new rockburst intensity measurement model for underground engineering by variable weight theory and matter-element extension theory. In the proposed model, variable weight theory is introduced to optimize the constant weight of predictors so as to improve the rationality of weighting deamination. Grade variable method and matter-element extension method are combined to establish the grade variable interval of rockburst intensity prediction, corresponding to the classification standard of rockburst intensity grade. The rockburst intensity grade is determined in line with grade variable interval of rockburst intensity. Because the grade variable comprehensively reflects the comprehensive extension correlation information of each rockburst intensity grade, and the correct rate of rockburst intensity prediction can be improved.

Theoretical framework

2.1 The matter-element extension theory

In the matter-element extension theory, the change of objects indicates the exploitation, and the change possibility is called extensibility, and the extensibility of objects is described by the matter-element extension. If an object S has the feature y, the eigenvalue of feature y is expressed by v, then the ordered triple R = (S, y, v) is used to describe the basic-element of the object, abbreviated as matter-element. If a matter S has n features, the corresponding eigenvalues of features y1, y2, ⋯, yn are v1, v2, ⋯, vn respectively. The matter-element describing the Subject S is denoted as R [39].

(1)

The matter-element classical domain of object S is Rj.

(2)

The matter-element joint domain of object S is R0.

(3)

Where S is the object, vi is the eigenvalue, Sj is the object corresponding to feature greades of j, j∈{1,⋯,m}, m is the number of feature greades, yi is the features of object i, i∈{1,⋯,n}, n is the number of features, Vji is the eigenvalues range of Sj about yi, S0 is the object corresponding to overall feature grades, V0i is the eigenvalues range of S0 about yi.

According to the extension set and extension distance [40], extending "distance" to "extension distance", the real variable function in classical mathematics is extended to correlation function then the extension distance formula of feature yi in the object S in regard to feature grade j can be expressed as:(4)where vi is the object feature i, Vji is the value range of Sj in regard to feature yi, Vij = [aij,bij], V0i is the value range of S0 with respect to feature yi, V0i = [a0i,b0i], ρ(vi,Vji) is the extension distance of feature yi with respect to the classical domain of j feature grade, and ρ(vi,V0i) is the extension distance of feature yi with respect to joint domain.

If vi∊V0i, the elementary function is obtained for calculating the feature correlation [41].(5)where Sj(vi) is the feature correlation of the feature yi of the object S with respect to the feature grade j. Combining with the feature weight vector W, the comprehensive correlation degrees of the object S with respect to the feature grade j is calculated as follows [42].

(6)

Where wi is the weight coefficient of the feature yi, Sj(Vj) is the comprehensive correlation degrees between the object S and the feature grade j, Sj(vi) is the single index correlation of the feature yi of the object S with respect to the feature grade j.

2.2 Variable weight theory

Since the weight value is always fixed in the process of constant weighting, it can only reflect the relative importance of each feature of objects, ignoring the influence caused by eigenvalue change to the feature weight [43]. In this regard, Zhu Y Z and Li H X improved the variable weight theory [44]. According to the variable weight theory, constructing the variable weight vector can optimize the constant weight of objects, and obtain the feature variable weight.

According to the axiomatic system, the definition of variable weight vectors [45] as follows. There are m mappings Wj, j∈{1,⋯,m}, Wj: [0,1]m→[0,1] (x1, ⋯, xm)→Wj (x1, ⋯, xm). if Wj satisfies the following three conditions: Condition a. W1+ ⋯ + Wm = 1. Condition b. Wj for variable xj continuous. Condition c. Wj is monotone decreasing (or increment) for variable xj. W(X) = {W1(X), ⋯, Wm(X)} is called the penalty (or excitation) variable weight vector. The definition of feature variable weight vector [45] as follows. There is a mapping P:[0,1]m→[0,1]m, X→P(X) = {P1(X), ⋯ Pm(X)}. P is called an m dimensional feature variable weight vector, if P satisfies the following three conditions: Condition 1. xi≥xj→Pi(X)≤Pj(X) or xi≥xj→Pi(X)≥Pj(X). Condition 2. Pj(X) is continuous for each variable xj. Condition 3. for any constant weight vector W0 = {w01,⋯,w0m}, Eq (7) satisfies the condition a, condition b and condition c, P is called a penalty (or excitation) feature variable weight vector.

(7)

Essentially, feature variable weight vector Pj(X) is a gradient vector of m-dimensional real function B(x) (also known as balanced function) [46], whose calculation equation is as follows.

(8)

According to Eqs (7) and (8), the variable weight vector W(X) of the object X can be obtained.

Prediction model for rockburst intensity grade

3.1 Modeling solutions

Due the complex environment around the rock masses in underground engineering, many factors can result in the rockburst. Firstly, based on the analysis of the existing achievement, the scientific predictive index system of rockburst intensity grading was constructed, and the rockburst intensity predictive index grading standard was established. Secondly, the extreme value method was used to standardize the grading standards and the indicators of Engineering Rock Mass. The dimensionality of the forecasting indicators was achieved, and the normalized grading standards for rockburst prediction indicators were obtained. Thirdly, the each grades in the grading standard was regarded as an Ideal Rock Mass; the constant weight of the predictive indexes were calculated by using entropy method to synthesize the index information of the Ideal Rock Masses and the Engineering Rock Masses, and the difference of the predictive indexes was fully considered by using variable weight theory to calculate index variable weight. Finally, based on the matter-element extension theory, the comprehensive correlation degrees between the normalized grading and the Ideal Rock Masses, Engineering Rock Masses, were calculated, and whose grading variables of rockburst were calculated. The prediction grading interval of rockburst intensity was constructed by using the rockburst grading variables of Ideal Rock Masses. The rock burst intensity grade of Engineering Rock Masses can be obtained by judging the grading interval of the variable value of Engineering Rock Masses.

3.2 Predictive index system and grading standard

According to a large number of test results and practical experience, the factors causing rockburst can be divided into internal factors and external factors. Internal factors mainly includes compressive strength, shear strength, brittleness, hardness and Poisson's ratio. External factors mainly includes blasting disturbance, tunneling depth, size and shape of underground space. These external factors can damage the integrity and spatial structure of the rock mass, thus changing the stress distribution of the surrounding rock. Rock burst is caused by both internal and external factors. Based on the existing research [47–50], the rockburst intensity predictive index system was established based on the ratio σθ/σc of shear stress to uniaxial compressive strength, the ratio σc/σt of rock uniaxial compressive strength to tensile strength, rock brittleness index Is and rock elastic energy index Wet. The Classification standard of rock burst intensity grading features is shown in Table 1, and the corresponding prediction index grading standard is shown in Table 2.

[Figure omitted. See PDF.]

Table 1. Classification standard of rockburst intensity grading features.

https://doi.org/10.1371/journal.pone.0218525.t001

[Figure omitted. See PDF.]

Table 2. The corresponding prediction index grading standard.

https://doi.org/10.1371/journal.pone.0218525.t002

To concisely calculate the variable interval corresponding to the rockburst intensity grade, grading criteria of the rockburst intensity is classified into five kinds of Ideal Rock Masses. The predictive indexes of rockburst intensity and their normalized values of these five Ideal Rock Masses are shown in Table 3, and the normalized method of rockburst intensity prediction index is shown in Section 3.3.

[Figure omitted. See PDF.]

Table 3. The rockburst intensity predictive index values and corresponding normalized values of Ideal Rock Masses.

https://doi.org/10.1371/journal.pone.0218525.t003

3.3 Weighting of prediction indicators

Firstly, the Ideal Rock Masses and the Engineering Rock Masses were subjected to dimensionless treatment by the extremum method, so that the dimension of each predictive index was consistent. Then the entropy method was used to determine the constant weight of each predictive index. Finally, the variable weight vector was constructed based on the variable weight theory to optimize the constant weight, and the variable weight of predictive index was obtained. Weighting steps of the rockburst intensity predictive index are as follows:

1. The extreme value method is used to standardize the predictive index xij of rock mass Xi and transform it into dimensionless value vij. The forward predictive index xij of rock mass Xi is disposed by Eq (9), and the reverse predictive index xij of rock mass Xi is disposed by Eq (10) [51]. n is the number of Engineering Rock Masses, m is the number of predictive indexes, i∈{1, ⋯, n}, j∈{1, ⋯, m}.(9)(10)

Where xij is the predictive index of rock mass Xi, mj is the maximum value of predictive index xij, mj represents the minimum value of predictive index xij, vij is the dimensionless value of xij.

2. The normalized predictive index value vij of rock mass Xi is processed and the predictive index information entropy ej of rock mass Xi is calculated.(11)(12)

Where vij is the normalized predictive index of the rock mass Xi, ej is the information entropy of the predictive index of the rock mass Xi, rij is the normalized predictive index value of the rock mass Xi.

3. the constant weight of predictive index of rock mass Xi is calculated.(13)

Where wj is the weight of the predictive index of the rock mass Xi, ej is the information entropy of the predictive index of the rock mass Xi, and m is the number of the predictive indexes of the rock mass Xi.

The constant weight vector W0 = {w01,⋯,w0m} can be obtained.

4. According to the variable weight theory, the construction and selection of state variable weight vectors are the core of variable weight implementation. The common feature variable weight vectors include sum type, product type, exponential type and mixed type. Since the exponential feature variable weight vector can determine the appropriate parameters by given different equilibrium forces, it is more flexible [52]. The feature variable weight vectors Pj(X) = (Pj(x1),⋯,Pj(xm)) of exponential type are used. The calculation equation of Pj (xi) is as follows.(14)

Where xij is the predictive index of rock mass Xi, β indicates the adjustment degree of variable weight vector to the vector equilibrium. β is selected according to the actual situation and β > 0.

3.4 Prediction of the rockburst intensity grade

According to matter-element extension theory and rockburst intensity predictive index standard, the classical domain of rockburst intensity predictive index can be expressed as Rj, and the nodal domain as R0. The classical domain Rj expresses the change range of normalized index value in each feature grade of slope rock mass evaluation index, and the nodal R0 expresses the whole range of normalized value of rock mass classification index.

(15)(16)

According to the grading standard of rockburst intensity, if a predictive index of Engineering Rock Mass is greater than the maximum index, the maximum index is used as the index of Engineering Rock Mass. if it is less than the minimum index, the minimum index is used as the index of Engineering Rock Mass. Then the rockburst intensity predictive indexes of each Engineering Rock Massare normalized according to the Eqs (9) and (10), and normalized predictive index values are obtained. Finally, by Eqs (4), (5) and (6), the comprehensive correlation degrees between each Ideal Rock Mass, each Engineering Rock Mass and the rockburst intensity prediction grading standard are calculated. The grading variable values of rockburst intensity grade of Ideal Rock Masses and Engineering Rock Masses are calculated based on the grading variable method [53]. The equation for calculating grading variables k is as follows.(17)(18)where Sj(Vj) is the comprehensive correlation degree between S and feature grading j, Pj(Vj) is the normalized value of comprehensive correlation, k is the value of feature grading of S,Bs = min{Sj(Vj)}, As = max{Sj(Vj)}.

Finally, the grading interval of grading variable is constructed by using the grading variables of Ideal Rock Masses. The rockburst intensity grade of a Engineering Rock Mass is determined by distinguishing the grading interval corresponding to the grading variable of the Engineering Rock Mass.

Model verification and results analysis

4.1 Cases selection

To test the accuracy and safety of the proposed model, several classical rockburst cases of underground engineering at home and abroad [54] are selected for the prediction and analysis.

1. Case 1: Tianshengqiao II Hydropower Station Headrace Tunnels;

2. Case 2: Ertan hydropower station 2;

3. Case 3: Underground Tunnels of Lubuge Hydropower Station;

4. Case 4: Yuzixi Hydropower Station Headrace Tunnels;

5. Case 5: Taipingyi Hydropower Station Headrace Tunnels;

6. Case 6: Pingjin II Hydropower Station Headrace Tunnels;

7. Case 7: Parallel adit K261+9398 of Erlangshan Tunnels;

8. Case 8: Zhongnanshan Extra Highway Tunnels of Qinling Mountains;

9. Case 9: Jiuhuashan Tunnels of Fu’ning Expressway in Fujian Province

10. Case 10: Heggura Highway Tunnels of NorwayVietas Hydropower Station Headrace Tunnels of Sweden

11. Case 11: Vietas Hydropower Station Headrace Tunnels of Sweden

12. Case 12: Japanese Kan-Etsu Tunnels

13. Case 13: Taipingyi Hydropower Station Headrace Tunnels.

The measured values and normalized values of rockburst prediction in each case are shown in Table 4.

[Figure omitted. See PDF.]

Table 4. Predictive index values of rockburst cases in underground engineering at home and abroad.

https://doi.org/10.1371/journal.pone.0218525.t004

Rockburst cases in Table 4 are occurred in different countries and regions, so there is a large difference in the engineering geological conditions. This difference can be reflected in the distribution of the forecasting indicators (Fig 1). At the same time, As shown in Fig 1, there is no regularity between the predictive indexes of different cases in the model verification, and there is a great difference between the predictive indexes of the same case. It is impossible to directly determine the rockburst intensity grade of each case based on a single index. Therefore, these cases can be used to test the correctness of the prediction results of the model.

[Figure omitted. See PDF.]

Fig 1. The distribution of predictive indexes for each case.

Index 1 is σθ/σc, the ratio of rock shear stress and uniaxial compressive strength. Index 2 is σc/σt, the ratio of rock uniaxial compressive strength and tensile strength. Index 3 is Is, the brittleness index of rock. Index 4 is Wet, the elastic energy index of rock.

https://doi.org/10.1371/journal.pone.0218525.g001

4.2 Calculation and analysis of the weights of predictive indexes

In the exponential variable weight vectors, ?? is the penalty grade, which reflects the adjustment degree of the variable vectors weighting to vector equilibrium of the predictive index. The bigger the value, the more the optimal result is biased toward equilibrium among indexes. ?? = 1 is selected in this paper.

(19)

According to the theory of extreme value entropy weighting and variable weight, the constant weight of each rockburst prediction index and the variable weight value of each engineering case are calculated by using Eqs (7), (9)–(13) and (19). The results are shown in Table 5. The constant weight of prediction index synthesizes the prediction indexes in Table 5, which reflects the relative importance of each prediction index. In Fig 2, the relative importance of each prediction index is shown as: that index 1 is more important than Indicator 2, Index 2 is more important than Index 3, and Index 3 is more important than Index 4. Considering the different values of rocks in the prediction index, the variable weight of prediction index is obtained by adjusting the constant weight. It reflects the influence of different values of prediction index on the weight of prediction index. As shown in Fig 2, the relative importance of prediction indexes of different Engineering Rock Masses changes after introducing variable weights. There is no regularity in these changes, indicating the extensive representativeness of the selected rockburst cases in this paper.

[Figure omitted. See PDF.]

Fig 2. The distribution of variable weights for each case.

Index 1 is σθ/σc, the ratio of rock shear stress and uniaxial compressive strength. Index 2 is σc/σt, the ratio of rock uniaxial compressive strength and tensile strength. Index 3 is Is, the brittleness index of rock. Index 4 is Wet, the elastic energy index of rock.

https://doi.org/10.1371/journal.pone.0218525.g002

[Figure omitted. See PDF.]

Table 5. Weight values of predictive indexes for ideal rock mass and engineering cases.

https://doi.org/10.1371/journal.pone.0218525.t005

4.3 Prediction results and analysis

According to the Eqs (4)–(6) and (15)-(18), the comprehensive correlation degrees and the rockburst intensity grading variables between each Ideal Rock Mass and each prediction grading standard are calculated as shown in Table 6. The grading variable interval for constructing the rockburst intensity grading is shown in Table 7.As shown in Table 6, prediction index values of rockburst intensity of Ideal Rock Mass come from classification criteria. In this case, the comprehensive correlation degree of a Ideal Rock Mass is a one-to-one correspondence with its rockburst intensity grade. Therefore, according to the comprehensive correlation degree of each ideal rock mass in Table 6, the value of the characteristic grade variable of the calculated rockburst intensity is also a one-to-one correspondence relationship with the rockburst intensity grade, namely, the corresponding relationship in Table 7 is correct.

[Figure omitted. See PDF.]

Table 6. The comprehensive correlation degrees and rockburst intensity grading variables of each ideal rock mass.

https://doi.org/10.1371/journal.pone.0218525.t006

[Figure omitted. See PDF.]

Table 7. Rockburst intensity feature grades and grading variable interval.

https://doi.org/10.1371/journal.pone.0218525.t007

According to the Eqs (4)—(6), (15)—(18), the comprehensive correlation between the Engineering Rock Mass and the standard grade of rockburst intensity and the variable value of the characteristic grade of rockburst intensity are calculated. Table 8 shows the final prediction results. Model 1 in Table 8 refers to the cloud model based on index distance and uncertainty measurement [50]. Model 2 refers to the normal cloud model [38], Model 3 refers to the set pair theory model [36], and Model 4 refers to the efficiency coefficient method model [37]. Method 1 refers to the matter-element extension model based on variable weight theory with the discrimination criteria of the characteristic grade variables of rockburst intensity. Method 2 refers to the matter-element extension model that uses the variable weight with the discrimination criteria of the comprehensive relevance. Method 3 refers to the matter-element extension model that uses the constant weight with the discrimination criteria of the comprehensive relevance.

[Figure omitted. See PDF.]

Table 8. Prediction results 1.

https://doi.org/10.1371/journal.pone.0218525.t008

In Table 8, prediction results show that the prediction grades of Method 1 on Projects 2, 11 and 13 are lower than the actual grade, and the misjudgment rate is about 23%. In other projects, the prediction results are completely consistent with the actual situation. Method 2 has a lower prediction grade than the actual grade on Projects 4, 10 and 13, a lower prediction grade than the actual grade on Project 2; and a risk misjudgment (higher prediction grade than the actual grade) on Project 1, with a misjudgment rate of about 38%. Method 3 has a lower prediction grade than the actual grade on Project 13; and there are dangerous misjudgments on Projects 1, 9, 10 and 12, the false positive rate is about 38%. Thus, after applying variable weight theory and grading variable method in matter-element extension model, the misjudgment rate of the model is reduced by about 15%, and the occurrence of dangerous misjudgment is reduced.

In Table 9, predicted results show that Model 1 has dangerous misjudgments in Projects 1 and 8. All predictive grades of Model 2 are consistent with the actual grade. The prediction grade of model 3 on Projects 10 and 13 is one grade lower than the actual grade. The predictive grade of Model 4 in Project 2 is lower than the actual grade. Therefore, Method 1 is more secure than Model 1 in predicting results. Compared with Models 2, 3 and 4, the correct rate of Method 1 needs to be improved. However, Method 1 has a relatively high correct rate of prediction results, and the predicted grade is not greater than the actual grade. Hence, Method 1 is of great significance for guiding the safe construction of underground engineering.

[Figure omitted. See PDF.]

Table 9. Prediction results 2.

https://doi.org/10.1371/journal.pone.0218525.t009

Conclusions

1. Taking the ratio of tangential stress to uniaxial compressive strength of rock, uniaxial compressive strength to tensile strength of rock, brittleness index and elastic energy index of rock as prediction index system, the dimension of prediction index is unified by extreme value method. The objective constant weight of prediction index is calculated by entropy weight method, which reflects the relative importance of prediction index. On this basis, the variable weight theory is introduced to fully consider the influence of the difference of rockburst prediction indexes on index weighting, and the variable weight function of penalty characteristic variable weight function is selected to calculate the variable weight of prediction index, so as to improve the rationality of the weighting of prediction index.

2. Based on the grading standards of rockburst intensity index, five ideal rock masses are constructed. The matter-element extension theory and the grade variable method are used to calculate the characteristic grade variables of rockburst intensity for each ideal rock mass. Then value range of the rockburst intensity is constructed, corresponding to the grading standard. Based on this interval, the rockburst grade can be predicted by the variables of characteristic grade of rockburst intensity. This method overcomes the deficiency of matter-element extension model in predicting the rockburst grade with the maximum comprehensive correlation degree.

3. This model is also used to predict the rockburst intensity of several typical underground engineering cases at home and abroad. Results show that in a few projects, the prediction grade of this model is lower than the actual grade, accounting for about 23%. In most projects, the prediction grade is consistent with the actual grade, accounting for about 77%, and there is no large risk misjudgment or the difference between the prediction grade and the actual grade. Compared with the traditional matter-element extension model, the method of distinguishing rockburst intensity grade by the characteristic grade variable of rockburst intensity is proposed in this paper. The correct rate of prediction results is improved by about 15%, and the safety of prediction results is also higher.

4. Compared with normal cloud model, set pair theory model and efficacy coefficient method, the correct rate of the prediction results of this model needs to be further improved. Meanwhile, although the engineering cases collected in this paper are representative, the number is limited, and the reliability of the model needs to be further tested in engineering practice. In addition, in the process of calculating the variable weight of prediction index, the selection of variable weight function is significantly important to the weight calculation. A more reasonable variable weight is required to construct the function according to the characteristics of rockburst index.

Supporting information

[Figure omitted. See PDF.]

S1 File.

https://doi.org/10.1371/journal.pone.0218525.s001

(DOCX)

Citation: Chen J, Chen Y, Yang S, Zhong X, Han X (2019) A prediction model on rockburst intensity grade based on variable weight and matter-element extension. PLoS ONE 14(6): e0218525. https://doi.org/10.1371/journal.pone.0218525

References

1. Mansurov V A. Laboratory Experiments: Their Role in the Problem of Rock Burst Prediction. Rock Testing & Site Characterization. 1993: 745–771.

2. Adoko A C, Gokceoglu C, Wu L, Zuo Q J. Knowledge-based and data-driven fuzzy modeling for rockburst prediction. International Journal of Rock Mechanics & Mining Sciences. 2013; 61(4): 86–95.

3. Diederichs M S. Early assessment of dynamic rupture hazard for rockburst risk management in deep tunnel projects. Journal of the Southern African Institute of Mining & Metallurgy. 2018; 118(3): 193–204.

4. Kabwe E, Wang Y. Review on Rockburst Theory and Types of Rock Support in Rockburst Prone Mines. Open Journal of Safety Science & Technology. 2015; 05(4):104–121. Available from: http://file.scirp.org/Html/3-1480140_61983.htm.

5. Mazaira A, Konicek P. Intense rockburst impacts in deep underground construction and their prevention 1. Canadian Geotechnical Journal. 2015; 52(10): 1426–1439.

6. Li N, Jimenez R. A logistic regression classifier for long-term probabilistic prediction of rock burst hazard. Natural Hazards. 2018; 90(1):197–215.

7. Roohollah S F, Abbas T. Long-term prediction of rockburst hazard in deep underground openings using three robust data mining techniques. Engineering with Computers. 2018: 1–17.

8. Zhang J J, Fu B J. Rock burst and its criteria and control. Chinese Journal of Rock Mechanics and Engineering. 2008; 27 (10): 2035–2042.

9. Xie H P, Ju Y, Li L Y. Criteria for strength and structural failure of rocks based on energy dissipation and energy release principles. Chinese Journal of Rock Mechanics and Engineering. 2005; 24(17): 3003–3010.

10. Cai M F, Ji D, Guo Q F. Study of rockburst prediction based on in-situ stress measurement and theory of energy accumulation caused by mining disturbance. Chinese Journal of Rock Mechanics and Engineering. 2013; 32(10): 1973–1980.

11. Cai M F. Prediction and prevention of rockburst in metal mines—A case study of Sanshandao gold mine. Journal of Rock Mechanics & Geotechnical Engineering. 2016; 8(2): 204–211.

12. Heib M A, Lafont G, Didier C, Gunzburger Y, Piguet J P, Pytel W, et al. Critical energy and mine rockburst potential. ISRM International Symposium (EUROCK). Wroclaw, POLAND. 2013 Sep 23–26. Rock Mechanics for Resources, Energy and Environment. 2013: 551–556.

13. Jiang Q, Feng X T, Xiang T B, Su G S. Rockburst characteristics and numerical simulation based on a new energy index: a case study of a tunnel at 2,500 m depth. Bulletin of Engineering Geology & the Environment. 2010; 69(3): 381–388.

14. He M C, Miao J L, Feng J L. Rock burst process of limestone and its acoustic emission characteristics under true-triaxial unloading conditions. International Journal of Rock Mechanics & Mining Sciences. 2010; 47(2): 286–298.

15. He M C, Zhao F, Zhang Y, Du S, Guan L. Feature evolution of dominant frequency components in acoustic emissions of instantaneous strain-type granitic rockburst simulation tests. Rock and Soil Mechanics. 2015; 36(1): 1–8.

16. Kuksenko V S, Inzhevatkin I E, Manzhikov B T, Stanchits S A, Tomilin N G, Frolov D I. Physical and methodological principles of rock burst prediction. Soviet Mining. 1987; 23(1): 6–17.

17. Feng X T, Liu J B, Chen B R, Xiao Y X, Feng G L, Zhang F P. Monitoring warning and control of rockburst in deep metal mines. Engineering. 2017; 3(4): 538–545.

18. Ma T H,Tang C A,Tang L X, Zhang W D, Wang L.Mechanism of rock burst forcasting based on micro-seismic monitoring technology. Chinese Journal of Rock Mechanics and Engineering. 2016; 35(3): 470–483.

19. Liu X H, Wu A X, Wang C L, Wang H J,Wang Y M. Study on rock burst forecasting prediction in a deep mine. Journal of Mining & Safety Engineering. 2012; 29(1): 78–83.

20. Chen Z H, Tang C A, Huang R Q. A double rock mass sample model for rockbursts. International Journal of Rock Mechanics & Mining Sciences. 1997; 34(6): 991–1000.

21. Sharan S K. A finite element perturbation method for the prediction of rockburst. Computers & Structures. 2007; 85(17): 1304–1309.

22. Labuz J F, Zang A. Mohr–Coulomb Failure Criterion. Rock Mechanics and Rock Engineering. 2012, 45(6): 975–979.

23. Hoek E, Brown E T. Underground excavations in rock. The Institute of Mining and Metallurgy, 1980.

24. Zhou J, Shi X Z, Dong L, Hu H Y, Wang H Y. Fisher discriminant analysis model and its application for prediction of classification of rockburst in deep-buried long tunnel. International Journal of Coal Science & Technology. 2010; 16(2): 144–149.

25. Jian Z, Shi X Z. Fisher discriminant analysis method for prediction of classification of rock burst risk. Journal of China Coal Society. 2010; 35(S1): 22–27.

26. Guo L, Li X B, Yan X M, Xiong L H. Rock Burst Prediction Methods Based on BP Network Theory. Industrial Safety & Dust Control. 2005(10):32–35.

27. Liu Z J,Yuan Q P,Li J L. Application of fuzzy probability model to prediction of classification of rockburst intensity.Chinese Journal of Rock Mechanics and Engineering. 2008; 27(S1): 3095–3103.

28. Gao W. Prediction of rock burst based on ant colony clustering algorithm. Chinese Journal of Geotechnical Engineering. 2010; 32(6): 874–880.

29. Pei Q T, Li H B, Liu Y Q, Liu J T. Rockburst prediction based on a modified grey evaluation model. Chinese Journal of Rock Mechanics and Engineering. 2013; 32 (10): 2088–2093. 2013.10.017.

30. Xu F, Xu W Y. Projection pursuit model based on particle swarm optimization for rock burst prediction. Chinese Journal of Geotechnical Engineering. 2010; 32(5): 718–723.

31. Zhou X C, Bai C H, Wang Z Q, Ling D C.Application of projection pursuit model and particle swarm optimization in rockburst prediction. Journal of Shanghai Jiaotong University. 2012; 46 (12): 1956–1961.

32. Zhou J, Shi X Z, Huang R D, Qiu X Y, Chen C. Feasibility of stochastic gradient boosting approach for predicting rockburst damage in burst-prone mines. Transactions of Nonferrous Metals Society of China. 2016; 26(7): 1938–1945. 0.1016/s1003-6326(16)64312-1.

33. Zhou J, Li X, Shi X. Long-term prediction model of rockburst in underground openings using heuristic algorithms and support vector machines. Safety Science. 2012; 50(4): 629–644.

34. Zhou J, Li X B, Mitri H S. Classification of rockburst in underground projects: comparison of ten supervised learning methods. Journal of Computing in Civil Engineering. 2016; 30(5): 04016003.

35. Jia Y P, Lu Q, Shamg Y Q, Zhi M M, Du L L. Rockburst prediction based on evidence theory. Chinese Journal of Geotechnical Engineering. 2014; 6(6): 1079–1086.

36. Chen H H, Li X B, Zhang Y. Study on application of set pair analysis method to prediction of rockburst. Journal of University of South China. 2008; 22(4): 10–14.

37. Wang Y C, Shang Y Q, Sun H Y, Yan X S. Study of prediction of rockburst intensity based on efficacy coefficient method. Rock and Soil Mechanics. 2010; 31 (2): 529–534.

38. Wang Y C, Jing H W, Zhang Q, Yu L Y, Xu Z M. A normal cloud model-based study of grading prediction of rockburst intensity in deep underground engineering. Rock and Soil Mechanics. 2015; 36(4): 1189–1194.

39. Cai W. Introduction of Extenics. System Engineering Theory and Practice. 1998(1): 76–84.

40. Cai W, Yang C Y. Basic theory and methodology on extenics. Chinese Science Bulletin. 2013; 58(13): 1190–1199.

41. Yang C Y, Cai W. Recent research progress in dependent functions in extension sets. Journal of Guangdong University of Technology. 2012; 29(2): 7–14.

42. Shi Z P, Shan T H, Liu W F, Zhang X P. Comprehensive evaluation of power network operation risk based on matter-element extensible model. Power System Technology. 2015; 39 (11): 3233–3239.

43. Tang X W, Zhou Z F, Shi Y. The variable weighted functions of combined forecasting. Computers and Mathematics with Applications. 2003; 45(2): 723–730.

44. Zhu Y Z, Li H X. Axiomatic system ofstate variable weights and construction of balance functions. Systems Engineering-Theory and Practice. 1999; 19(7): 116–131.

45. Li D Q, Li H X. The properties and construction of state variable weight vectors. Journal of Beijing Normal University:Natural Science. 2002; 38(4): 455–461.

46. Li H X. Factor spaces and mathematical frame of knowledge represention (Ⅷ). Fuzzy Systems and Mathematical. 1995; 9(3): 1–9.

47. Gong F Q, Li X B. A distance discriminant analysis method for prediction of possibility and classification of rockburst and its application. Chinese Journal of Rock Mechanics and Engineering. 2007; 26(5): 1012–1018.

48. Jia Z Q, Zhang R, Zhang Y F, Gao M Z, Liu Q Y, Li A Q et al. Integrated prediction for rockburst of underground powerhouse of Houziyan hydropower station. Chinese Journal of Geotechnical Engineering. 2016; 38(S2): 110–116.

49. Wang X B, Tao S, Pan Y S, Wang W. Numerical simulation of the distribution of rockburst regions with different magnitudes based on loop stress criteria. Journal of Sichuan University (Engineering Science Edition). 2011; 43(5): 7–12.

50. Zhang B, Dai X G. A cloud model for predicting rockburst intensity grade based on index distance and uncertainty measure. Rock and Soil Mechanics. 2017; 38(S2): 257–265.

51. Xu H J, Zhao B F, Zhou Y, Liu S X. Evaluation on water disaster from roof strata based on the entropy-weight and matter-element extension model. Journal of Mining&Safety Engineering. 2018; 35(01): 112–117.

52. Li D Q, Hao F L. Weights transferring effect of state variable weights vector. Systems Engineering: Theory & Practice. 2009; 29(6): 127–131.

53. Hu B Q, Zhang X. The improvement and application of extension evaluation method. Journal of Wuhan University. 2003; 36(5): 79–84.

54. Xiong X B, Gui G Q, Zheng M Q, Xu J C, Wang B J, Zhao Q H, et al. Study on application of extension method to prediction of rockburst in deep and over-length highway tunnel. Symposium of the 10th Learning Congress About Rock Mechanics and Engineering. Beijing: China Electric Power Press, 2008:403–411.