Next Article in Journal
Development and Validation of an Interdisciplinary Worker’s Health Approach Instrument (IWHAI)
Previous Article in Journal
Geospatial Analysis of Mass-Wasting Susceptibility of Four Small Catchments in Mountainous Area of Miyun County, Beijing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Key Factors Identification and Risk Assessment for the Stability of Deep Surrounding Rock in Coal Roadway

1
State Key Laboratory of Mining Disaster Prevention and Control Co-Founded by Shandong Province and the Ministry of Science and Technology, Shandong University of Science and Technology, Qingdao 266590, China
2
College of Mining and Safety Engineering, Shandong University of Science and Technology, Qingdao 266590, China
*
Authors to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2019, 16(15), 2802; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16152802
Submission received: 8 July 2019 / Revised: 24 July 2019 / Accepted: 30 July 2019 / Published: 6 August 2019

Abstract

:
In order to evaluate the stability of deep surrounding rock, all of the affecting factors should be theoretically identified. However, some factors have slight impacts on the stability of deep surrounding rock compared with others. To conduct an effective risk assessment, key factors should be first extracted. The analytic hierarchy process (AHP) and grey relation analysis (GRA) methods are integrated to determine the key factors. First, the AHP method is applied to sort the factors by calculating the weights of them. Seven out of fifteen factors are extracted as the key factors, which account for 80% of the weights. Further, the GCA method is used to validate the effects of these key factors by analyzing the correlation between the performance of each factor and that of the reference. Considering the influence of these key factors and experts’ judgements, the multilevel fuzzy comprehensive evaluation method is adopted to obtain the risk level of the deep surrounding rock stability. Finally, the risk assessment of the deep surrounding rock in the E-Zhuang coal mine of Chinese Xinwen Mining Area illustrates the operability of the proposed method.

1. Introduction

The geological environment of coal mine is complex and its ecological environment is fragile. In recent years, the conditions of coal mines in China have gradually developed into deep and complex ones. The depth and breadth of mining have been increasing and rock engineering has been more difficult. Meanwhile, the gradual enlargement of the mining scale and the continuous improvement of the mechanization level have brought great difficulties to the stability control of surrounding rock in mines. In addition, factors such as stress and mining influence have led to the damage of weak rock surfaces in the process of excavation and artificial disturbance. Once the whole rock mass is destroyed, the instability of the surrounding rock could lead to a series of environmental disasters, such as collapse and roof separation, which result in great harm to the geological environment. From 2013 to 2017, nearly 39% of coal mine accidents in China were caused by roof accidents. The particular environment of the mine and these disasters indicate that the stability of the deep surrounding rock of the coal roadway is an urgent problem for safe mining. Therefore, it is of great significance to evaluate the stability of deep surrounding rock, which is not only related to the safety of people’s lives and property, but also to an urgent problem of safety in mining and environmental protection.
Deep surrounding rock plays an important role in the coal mining process. The instability of the deep surrounding rock will result in accidents such as roof caving and wall caving. During the years 2013–2017, nearly 39% of coal mining fatalities in China were caused by roof accidents. Substantial disasters have indicated that the stability of the deep surrounding rock in the coal roadway is an urgent problem for safe mining. Therefore, it is significant to conduct the evaluation of the deep surrounding rock stability.
Many works on surrounding rock stability evaluation are now available. For example, Kentil et al [1] studied the engineering geological properties of exposed rocks on the Ankara-Pozanti Expressway and evaluated the stability of excavated slopes. Andonov [2] evaluated the geomechanical conditions for establishing a free channel irrigation system and proposed some suggestions for improving the stability of rock mass. Saffari et al [3] applied a rock engineering systems approach to evaluate and classify the coal spontaneous combustion potential in the Eastern Alborz coal mines. Wang et al [4] used the ideal point method to categorize the affecting indices of surrounding rock stability and combined the objective and subjective weighting method to obtain the grade of the surrounding rock. Mishra et al [5] proposed the geotechnical risk management concept and applied it to the intelligent deep mines; Galvin [6] emphasized the importance of risk management in ground engineering and provided many management strategies. Babaeian et al [7] proposed a new framework for rock evaluation in open pit mines, integrating the multivariate regression method and the decision-making trial. Santos et al [8] proposed a method to predict the stability of the mine rock slope based on the principal component analysis and discriminant analysis. It can be seen that the stability of surrounding rock is determined by multiple indices and there is no consensus on stability affecting factors.
In fact, the stability of the deep surrounding rock is affected by not only geological features, but environmental factors and engineering factors. Theoretically, all of these affecting factors should be considered for a comprehensive risk assessment of the deep surrounding rock stability. However, it could be time-consuming to include all of the affecting factors in the assessment model. In fact, some of these factors only have a small impact on the stability of the deep surrounding rock. Therefore, a scientific method is required to extract those with a high important degree. One efficient method is to rank these factors and choose the top-ranked ones [9,10,11]. The analytic hierarchy process (AHP) has been proven to be a useful method for selecting indicators [12,13,14,15,16,17,18], but the relative importance between every two factors in AHP is determined by subjective judgment. There lacks the verification of these affecting factors. On the other hand, with the complexity of affecting factors, it is impracticable to obtain a definite value to express the stability of the deep surrounding rock, especially in situations where many of these factors cannot be measured precisely and subjective judgment from experts is inevitable. Therefore, a method to synthesize the key affecting factor and the experts’ knowledge is required.
The novelty of this paper is to identify and verify the key factors of the deep surrounding rock stability and to establish a model for a comprehensive risk assessment. The rest of this paper is organized as follows: Section 2 introduces the key factors identification method and a fuzzy risk assessment method. Section 3 is the application of these methods with the E-Zhuang coal mine of Chinese Xinwen Mining Area as a case study. Section 4 provides conclusions.

2. Methodology

2.1. Key Factors Identification Method

Theoretically, all of the affecting factors should be identified to evaluate the stability of the deep surrounding rock, but with the limits of time and resource, only those contributing significantly to the stability of the deep surrounding rock can be considered. We named these factors as key factors. On the other hand, different affecting factors of the stability of the deep surrounding rock have different dimensions. Besides, some factors such as the seepage quantity of groundwater can be quantified, while others such as the purposes and age limit of the roadway cannot. All of these features of the affecting factors of the deep surrounding rock stability result in the inadaptability of any mathematical function. Therefore, experts’ knowledge and experience should be considered in a good way.
In order to determine the key factors contributing significantly to the stability of the deep surrounding rock in the coal roadway, an integrated method combining the AHP and GCA is proposed. The principle of the integrated method is shown as follows:
(1) Identify all the affecting factors (Fi, i = 1, 2, …, n) and construct a hierarchy structure. The analytical object is named the target (top) layer. Then the categories of factors are classified in the criterion (middle) layer and all of the affecting factors compose the index (bottom) layer.
(2) Divide units for the hierarchy structure. Each factor in the upper layer and its affecting factors are included in one unit. For example, there is a three-layered structure, and the unit division is shown in Figure 1.
(3) Represent the relative importance of different factors in each unit with the numbers 1–9 or the multiplicative inverse. Then the judgment matrices are determined. An example is shown in Equation (1). Assume B is in the criterion layer. C1, C2, and C3 are the affecting factors of B, which are in the index layer. According to expert judgement, C2 is relatively more important than C1 but less important than C3, and C3 is much more important than C1. Then the value assigned to C3 is bigger than the value of C2 and the value assigned to C2 is bigger than the value of C1. Note that the relative value determines the relative importance of different factors. The consistency of each judgment matrix should be checked, and the details can be seen in reference [19].
R = B     C 1   C 2   C 3 C 1 C 2 C 3 [ 1 2 3 1 / 2 1 1 / 2 1 / 3 2 1 ]
(4) Calculate the maximum characteristic root and the feature vector. In MATLAB, we can obtain the characteristic roots and the feature vector with the function eig, i.e., [x, y] = eig(A). The feature vector to which the maximum characteristic root corresponds is the weight vector. Further, the weight vector should be normalized.
(5) Calculate the combination weights. With the judgment matrices, we can only get the weight of each factor in each individual unit. The objective is to obtain the weight of each factor in the index layer for the target layer. Therefore, the weights should be combined. Assume that the weight vectors in each unit shown in Figure 1 are represented as WCB1, WCB2, WCB3, WBA, then the combination weights can be calculated as:
W C B 1 A = W B A × W C B 1
W C B 1 A = W B A × W C B 1
W C B 1 A = W B A × W C B 1
(6) Sort the affecting factors according to their weights and extract the key factors, which account for 80% of the weights.
(7) Validate the results of AHP with GRA. The correlation between the measured value and reference value can be reflected with GRA (Grey Relation Analysis). Accordingly, the measured value of the key factors should be highly correlated with the reference value. To verify, the indicators reflecting the measured value and the reference value should be first defined.
① Define the correlation indicators and sequences.
For the affecting factors of the deep surrounding rock stability, the indicators include the possibility of accidents, the time of exposure to hazardous conditions, consequence, and current control measures. The grading standards of these indicators are defined as shown in Table 1. The score is determined by considering the comprehensive correlation factors in the MLS method. MLS is one of the risk assessment methods for the working conditions. The control correlation strength indices of the surrounding rock stability of the deep roadway factors are chosen, which separately are: damage possibility (accident possibility), cycle strength (exposure time to dangerous environment), disaster degree (possible consequence), and control measures (existing prevention and control measures). The value range of the four indices is also determined by this method, and then the correlation coefficient and correlation degree are calculated using MATLAB. Here the relevance degree is a comprehensive representation of the importance level of each factor. The smaller the degree, the greater the importance.
For each affecting factor, there is a parameter sequence consisting of these four indicators, which is Xi = {value 1 (possibility of accidents), value 2 (time of exposure to hazardous conditions), value 3 (consequence), and value 4 (current control measures)}. In addition, a reference sequence is defined as X0 = {0.1, 0.5, 1, 1} based on the optimum value of each indicator.
② Calculate the correlation coefficient.
Assume that X 0 = { x 0 ( t k ) , k = 1 , 2 , , n } represents the reference sequence and X i = { x i ( t k ) , k = 1 , 2 , , n } represents the parameter sequence. Calculate the absolute difference > i ( n ) between the value of each indicator x i ( n ) and the reference value x 0 ( n ) [20,21]. The equation is as follows:
> i ( n ) = | x 0 ( n ) x i ( n ) |
Then the correlation coefficient ξ i ( n ) can be calculated based on Equation (6):
ξ i ( n ) = > min + ρ > max > i ( n ) + ρ > max = min i | x 0 x i | + ρ min i | x 0 x i | | x 0 x i | + ρ max i | x 0 x i | , n = 1 , 2 , 3 , 4 ;   i = 1 , 2 , , n
The smaller that ρ is, the bigger the correlation coefficient difference. In general, ρ = 0.5 .
③ Calculate the correlation degree.
Since there are four indicators for each factor, the correlation degree should reflect all of these four aspects. Therefore, the mean value of correlation coefficients is defined as the correlation degree. The smaller the correlation degree, the more dangerous is the affecting factor. The equation is shown as follows:
γ i = 1 n k = 1 n ξ i ( k )

2.2. Fuzzy Risk Assessment Method

The risk of the deep surrounding rock stability can be divided into five levels. The risk evaluation set is shown as follows:
V = (V1, V2, V3, V4, V5) = {very low, relatively low, general, relatively high, very high}
(1)
Establish the risk assessment matrix ( M i j ) k × m . K experts are invited to evaluate the key factors of the deep surrounding rock stability. The matrix ( M i j ) k × m is shown as follows:
( M i j ) k × m = ( M 11 M 12 M 1 m M 21 M 22 M 2 m M k 1 M k 2 M k m )
(2)
Determine the membership matrix.
Segmented grey whitening weight functions are used to transform the evaluation value into membership degrees. Therefore, the evaluation value is the independent variable, and the memberships are the dependent variables. In response to five risk evaluation levels, five-segmented functions are defined.
(1) The grey whitening weight function is shown in Equation (10) when the risk evaluation level is “very high”:
f 1 ( d i ) = { 0 d i / 5 1    d [ 0 , 10 ] d [ 0 , 5 ] d [ 5 , 10 ]
(2) The grey whitening weight function is shown in Equation (11) when the risk evaluation level is “relatively high”:
f 2 ( d i ) = { 0 d i / 4 ( 8 d i ) / 4    d [ 0 , 8 ] d [ 0 , 4 ] d [ 4 , 8 ]
(3) The grey whitening weight function is shown in Equation (12) when the risk evaluation level is “general”:
f 3 ( d i ) = { 0 d i / 3 ( 6 d i ) / 3    d [ 0 , 6 ] d [ 0 , 3 ] d [ 3 , 6 ]
(4) The grey whitening weight function is shown in Equation (13) when the risk evaluation level is “relatively low”:
f 4 ( d i ) = { 0 d i / 2 ( 4 d i ) / 2    d [ 0 , 4 ] d [ 0 , 2 ] d [ 2 , 4 ]
(5) The grey whitening weight function is shown in Equation (14) when the risk evaluation level is “very low”:
f 5 ( d i ) = { 0 1 2 d i    d [ 0 , 2 ] d [ 0 , 1 ] d [ 1 , 2 ]
Then after normalization, each key factor’s membership degree to each risk level can be obtained with the grey whitening weight functions. For the factor j, its membership degree to risk level e can be represented with Equation (15):
s e = i = 1 k f e ( M i j ) e = 1 5 i = 1 k f e ( M i j )
where i refers to the experts and i = 1, 2, …, k; j represents the factor and j = 1, 2, …, m; e means the risk level; and e = 1, 2, 3, 4, 5. For factor j, the membership degree vector is shown below:
S 11 = ( s 111 , s 112 , s 113 , s 114 , s 115 )
Similarly, the membership degree vectors of other factors can be obtained. Then the membership matrix is shown in Equation (17):
S = ( S 11 S 12 S i j ) ( s 111 s 112 s 115 s 121 s 122 s 125 s i j 1 s i j 2 s i j 5 )
(3)
Multi-level fuzzy comprehensive evaluation.
To evaluate a system with affecting factors interrelated in a hierarchy structure, the multi-level fuzzy comprehensive evaluation has proven to be an effective method [22,23]. The membership degree of each factor in the lower level to its corresponding upper level factor can be obtained with a fuzzy comprehensive evaluation [24,25,26]. The weight set consists of weights of key factors, represented as {a,a,a}. By integrating the weight set with membership matrix through a fuzzy operator, the fuzzy comprehensive evaluation results are available. For a three-layered structure, the fuzzy comprehensive evaluation should be conducted twice. First, calculate the membership degrees of one factor in the intermediate layer to each risk level. The equation is as follows:
B i = A i · S i = ( a i 1 , a i 2 , a i j ) ( s i 11 s i 1 k s i j 1 s i j k ) = ( b i 1 , b i 2 , b i k )
Similarly, the membership degrees of other factors in the intermediate layer to each risk level can be obtained, and all these membership degree sets form a membership degree matrix:
S = ( B 1 B 2 B m ) ( A 1 · S 1 A 2 · S 2 A m · S m )
Then the risk level of the object can be calculated according to Equation (20):
B = A · S = ( a 1 , a 2 , , a m ) ( b 11 b 1 k b m 1 b m k ) = ( b 1 , b 2 , , b k )

3. Applications

3.1. Affecting Factors of Deep Surrounding Rock Stability

The stability of the deep surrounding rock is affected by geological features such as rock strength, environmental factors such as ground stress, and engineering factors such as the thickness of strengthened rock [27,28,29]. To evaluate the risk of the deep surrounding rock, all these aspects should be considered. Based on field experience and experts’ knowledge, 15 affecting factors are identified for the risk evaluation of an E-Zhuang coal mine in Chinese Xinwen Mining Area. The hierarchy structure model for the deep surrounding rock stability affecting factors is shown in Figure 2. However, some of them have little impact on the stability of the deep surrounding rock. The methods introduced above are applied to extract the key ones.

3.2. Key Factors Identification Based on AHP and GCA

3.2.1. Determine the Weights of Each Factor with AHP

According to AHP, pairwise comparisons are conducted for the 15 affecting factors. The scoring results are shown in Table 2.
According to the calculation rules introduced in Section 2.1, the weights of each factor can be calculated. The analysis and ranking results are shown in Table 3. The bigger the weight, the more important the factor. It can be seen that the top 7 factors’ weights account for 0.8151, which means that these 7 factors are key factors for the deep surrounding rock stability.
The results are in consistent with the actual conditions of the deep surrounding rock in the coal roadway. In fact, the rock strength, rock mass quality indicator RQD, and coal seam dip angle are the direct causes of the deep surrounding rock failure according to accident statistics. The mining influence coefficient and ground stress are the precipitating factors because they determine the environment of the surrounding rock. The sectional size of the roadway and the thickness of the strengthened rock are also important contributors since great uncertainties exist during engineering excavation.

3.2.2. Validate the Results of AHP with GRA

GRA (Grey Relation Analysis) extracts the key factors by calculating the correlation between the measured value and reference value. Four correlation indicators have been defined in Section 2.1, which are the possibility of accidents, the time of exposure to hazardous conditions, consequence, and current control measures. According to the grading standards of these indicators presented in Table 1, experts in the E-Zhuang coal mine of Chinese Xinwen Mining Area assigned evaluation values for each factor. Accordingly, there are 15 parameter sequences and one reference sequence, shown in Table 4, and corresponding normalized data are shown in Table 5.
According to Equations (5)–(8) mentioned in Section 2.1, the correlation coefficients and correlation degrees for each factor can be calculated, shown in Table 6. Then these 15 factors are sorted by their correlation degrees.
From Table 5 we can conclude that the rock strength, rock mass quality indicator RQD, coal seam dip angle, mining influence coefficient, ground stress, sectional size of the roadway, and thickness of strengthened rock are dominant factors of the deep surrounding rock stability in the coal roadway. The results validate the results of AHP.

3.3. Fuzzy Risk Evaluation

(1) Establish the evaluation index system.
According to the key factors identification results, a structural evaluation index system can be established, shown in Figure 3. Based on the index system, two layers of fuzzy comprehensive evaluation are defined. The first layer is composed of the bottom factors and the intermediate indices, and the second layer consists of the intermediate indices and the top object. The evaluation factor set for the second layer is represented as T = { T 1 , T 2 , T 3 } , and the evaluation factor set for the first layer is represented as M 1 = { M 11 , M 12 , M 13 } , M 2 = { M 21 , M 22 , M 23 } and M 3 = { M 31 , M 32 , M 33 } .
(2) Establish the weight set.
According to the AHP results, the weights of these seven key factors are available. Further, the normalized weights of them are calculated, as shown in Table 7. Take rock strength as an example; the calculation process is as follows:
Normalized   weight   ( rock   strength ) = 0.2948 0.2948 + 0.1150 + 0.0761 = 0.6067
Accordingly, the weight sets are represented as follows:
C 1 = ( 0.6067 0.2367 0.1566 )
C 2 = ( 0.5068 0.4932 )
C 3 = ( 0.6188 0.3812 )
(3) Determine the membership matrix.
Further, 10 experts were employed to evaluate these seven key factors. For each key factor, every expert assigns a number between 0 and 10 according to the interval definition for each risk evaluation level.
The risk of the deep surrounding rock stability can be divided into five levels. The risk evaluation set is V = (V1, V2, V3, V4, V5) = {very low, relatively low, general, relatively high, very high}. Then the evaluation matrix ( M i j ) 10 × 7 of these seven factors is obtained, shown in Equation (25), where M i j represents the risk evaluation value of factor j made by expert i, and where i refers to the experts and i = 1, 2, …, k; j represents the factor and j = 1, 2, …, m:
( M i j ) 10 × 7 = ( M 11 M 12 M 1 , 7 M 21 M 22 M 2 , 7 M 10 , 1 M 10 , 2 M 10 , 7 ) = [ 4 2 2 6 7 6 4 6 3 3 7 8 5 5 5 2 2 6 7 4 4 4 3 3 6 6 5 5 5 2 4 7 7 4 4 4 4 4 6 8 4 6 6 3 2 8 6 6 5 5 3 3 6 7 5 4 4 4 4 7 7 4 5 5 2 2 6 8 6 6 ]
For example, for the factor j, its membership to risk evaluation level “very low” can be calculated as:
s 1 = i = 1 k f 1 ( M i 1 ) e = 1 5 i = 1 k f e ( M i 1 ) = i = 1 10 f 1 ( M i 1 ) e = 1 5 i = 1 10 f e ( M i 1 ) = f 1 ( 4 ) + f 1 ( 6 ) + f 1 ( 5 ) + f 1 ( 4 ) + f 1 ( 5 ) + f 1 ( 4 ) + f 1 ( 6 ) + f 1 ( 5 ) + f 1 ( 4 ) + f 1 ( 5 ) i = 1 10 f 1 ( M i 1 ) + i = 1 10 f 2 ( M i 1 ) + i = 1 10 f 3 ( M i 1 ) + i = 1 10 f 4 ( M i 1 ) + i = 1 10 f 4 ( M i 1 ) = 9.2 9.2 + 8.0 + 4.0 + 0 + 0 = 0.4340
Similarly, the memberships of factor j to other risk evaluation levels can be calculated. Consequently, the membership degree set of factor j is shown below:
S 11 = [ 0.4340 0.3774 0.1887 0.0000 0.0000 ]
The other factors’ membership sets can be established and represented as:
S 1 = [ 0.4340 0.3774 0.1887 0.0000 0.0000 0.2105 0.2632 0.3008 0.2256 0.0000 0.2212 0.2765 0.2925 0.2098 0.0000 ]
S 2 = [ 0.7273 0.2727 0.0000 0.0000 0.0000 0.8163 0.1837 0.0000 0.0000 0.0000 ]
S 3 = [ 0.4462 0.3758 0.1780 0.0000 0.0000 0.4340 0.3774 0.1887 0.0000 0.0000 ]
S1 represents the risk evaluation memberships of three key factors that reflect geological features to risk levels. S2 represents the risk evaluation memberships of two key factors that reflect environmental features to risk levels. S3 represents the risk evaluation memberships of three key factors that reflect engineering features to risk levels.
(4) First layer fuzzy comprehensive evaluation.
According to the comprehensive evaluation principles, the first level risk evaluation of geological features (B1), environmental features (B2), and engineering features (B3) can be available. The calculation process is shown below:
B 1 = C 1 S 1 = ( 0.6067 0.2367 0.1566 ) [ 0.4340 0.3774 0.1887 0.0000 0.0000 0.2105 0.2632 0.3008 0.2256 0.0000 0.2212 0.2765 0.2925 0.2098 0.0000 ] = ( 0.3131 0.3012 0.2084 0.0777 0 )
B 2 = C 2 S 2 = ( 0.5068 0.4932 ) [ 0.7273 0.2727 0.0000 0.0000 0.0000 0.8163 0.1837 0.0000 0.0000 0.0000 ] = ( 0.5868 0.1741 0 0 0 )
B 3 = C 3 S 3 = ( 0.6188 0.3812 ) [ 0.4462 0.3758 0.1780 0.0000 0.0000 0.4340 0.3774 0.1887 0.0000 0.0000 ] = ( 0.3048 0.2598 0.1257 0 0 )
(5) Second layer fuzzy comprehensive evaluation.
The first level evaluation results (B1, B2, B3) are used as input data for the second layer fuzzy comprehensive evaluation. The calculating process is as follows:
S = ( B 1 B 2 B 3 ) ( C 1 · S 1 C 2 · S 2 C 3 · S 3 ) = [ 0.3131 0.3012 0.2084 0.0777 0 0.5868 0.1741 0 0 0 0.3048 0.2598 0.1257 0 0 ]
B = A S = ( 0.5396 0.1634 0.2970 ) [ 0.3131 0.3012 0.2084 0.0777 0 0.5868 0.1741 0 0 0 0.3048 0.2598 0.1257 0 0 ] = ( 0.3554 0.2681 0.1498 0.0419 0 )
(6) Result analysis.
The calculation results indicate the risk of the deep surrounding rock stability. According to the maximum membership principle (Li et al., 2015), the risk level of the deep surrounding rock in the E-Zhuang coal mine of Chinese Xinwen Mining Area is very low and the membership degree is 0.3554. This case illustrates that the combined approach can help identify key factors of the deep surrounding rock stability and obtain the risk level.

4. Conclusions

(1) The first and most important step of risk evaluation is to identify risk affecting factors. The stability of the deep surrounding rock is affected not only by geological features, but also by environmental factors and engineering factors. Since not all of these factors can be measured or quantified, especially man-made ones, a semi-quantitative approach is necessary. In order to extract the key affecting factors and evaluate the risk level of the deep surrounding rock stability, a method combing AHP, GCA, and FCE is presented.
(2) The traditional AHP has been a useful approach to obtain the weight of each factor. With GCA, the result of AHP can be verified. Therefore, the key affecting factors are identified, including rock strength, rock mass quality indicator RQD, coal seam dip angle, mining influence coefficient, ground stress, sectional size of roadway, and thickness of strengthened rock, which are dominant factors of the deep surrounding rock stability in the coal roadway. In addition, the weights of these main factors are calculated.
(3) Further, a multilevel FCE method is used to include both the influence of these key affecting factors and the experts’ judgements. The comprehensive method was applied to the risk assessment of the deep surrounding rock in the E-Zhuang coal mine of Chinese Xinwen Mining Area to obtain the risk level. The results show that the risk level of the deep surrounding rock is extremely low and support measures are necessary. The results can also provide a reference for other similar deep surrounding rock.
(4) Although the stability of surrounding rock in the deep mine has been well controlled with the gradual improvement of underground rock engineering construction, the traditional geotechnical engineering still faces a series of new environmental protection problems. With the gradual development of coal mines in China to deeper and more complex conditions, the problems of environmental effects and ecological environment protection have not been paid enough attention. The rock mass stability evaluation can provide theoretical support for ecological environment protection.

Author Contributions

D.H. contributed to the conception of the study; W.L. contributed to translation and visualization; X.C. contributed to analysis; and Y.T. helped perform the analysis with constructive discussions.

Funding

This research was funded by National Natural Science Foundation of China (51704179), Natural Science Foundation of Shandong Province (ZR2019BEE018, ZR2016EEB23), Scientific Research Foundation of Shandong University of Science and Technology for Recruited Talents (2017RCJJ002), and the First-class discipline construction projects (01AQ03703, 01CK05902). And The APC was funded by National Natural Science Foundation of China (51704179), Natural Science Foundation of Shandong Province (ZR2016EEB23).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kentli, B.; Topal, T. Evaluation of rock excavatability and slope stability along a segment of motorway, pozanti, turkey. Environ. Geol. 2004, 46, 83–95. [Google Scholar] [CrossRef]
  2. Andonov, E. Evaluation of Rock Mass Stability around Tunnels of the Irrigation System “Afrin”. J. Min. Sci. 2005, 41, 320–331. [Google Scholar] [CrossRef]
  3. Saffari, A.; Sereshki, F.; Ataei, M.; Ghanbari, K. Applying rock engineering systems (res) approach to evaluate and classify the coal spontaneous combustion potential in eastern alborz coal mines. Int. J. Min. Geo Eng. 2013, 47, 115–127. [Google Scholar]
  4. Wang, Y.; Zhang, J.; Guo, E.; Sun, Z. Fuzzy Comprehensive Evaluation-Based Disaster Risk Assessment of Desertification in Horqin Sand Land, China. Int. J. Environ. Res. Public Health 2015, 12, 1703–1725. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Mishra, R.; Janiszewski, M.; Uotinen, L.; Szydlowska, M.; Siren, T.; Rinne, M. Geotechnical Risk Management Concept for Intelligent Deep Mines. Procedia Eng. 2017, 191, 361–368. [Google Scholar] [CrossRef]
  6. Galvin, J. Critical role of risk management in ground engineering and opportunities for improvement. Int. J. Min. Sci. Technol. 2017, 27, 725–731. [Google Scholar] [CrossRef]
  7. Babaeian, M.; Ataei, M.; Sereshki, F.; Sotoudeh, F.; Mohammadi, S. A new framework for evaluation of rock fragmentation in open pit mines. J. Rock Mech. Geotech. Eng. 2019, 11, 325–336. [Google Scholar] [CrossRef]
  8. Santos, A.E.M.; Lana, M.S.; Cabral, I.E.; Pereira, T.M.; Naghadehi, M.Z.; da Silva, D.F.S.; dos Santos, T.B. Evaluation of rock slope stability conditions through discriminant analysis. Geotech. Geol. Eng. 2018, 37, 775–802. [Google Scholar] [CrossRef]
  9. Li, W.; He, M.; Sun, Y.; Cao, Q. A proactive operational risk identification and analysis framework based on the integration of ACAT and FRAM. Reliab. Eng. Syst. Saf. 2019, 186, 101–109. [Google Scholar] [CrossRef]
  10. Racicot, M.; Zanabria, R.; Paradis, M.; Gaucher, M.L.; Arsenault, J.; Leroux, A.; Letellier, A.; Tiwari, A.; Griffiths, M.; Charlebois, S.; et al. Identification of risk factors to be considered for food establishments’ risk assessment models. Microb. Risk Anal. 2019, 11, 1–10. [Google Scholar] [CrossRef]
  11. Yan, L.; Zhang, L.; Liang, W.; Li, W.; Du, M. Key factors identification and dynamic fuzzy assessment of health, safety and environment performance in petroleum enterprises. Saf. Sci. 2017, 94, 77–84. [Google Scholar] [CrossRef]
  12. Acharya, V.; Sharma, S.K.; Gupta, S.K. Analyzing the factors in industrial automation using analytic hierarchy process. Comput. Electr. Eng. 2018, 71, 877–886. [Google Scholar] [CrossRef]
  13. Hadidi, L.A.; Khater, M.A. Loss prevention in turnaround maintenance projects by selecting contractors based on safety criteria using the analytic hierarchy process (AHP). J. Loss Prev. Process. Ind. 2015, 34, 115–126. [Google Scholar] [CrossRef]
  14. Ho, W.; Ma, X. The state-of-the-art integrations and applications of the analytic hierarchy process. Eur. J. Oper. Res. 2018, 267, 399–414. [Google Scholar] [CrossRef]
  15. Lee, M. Strategies for Promoting the Medical Device Industry in Korea: An Analytical Hierarchy Process Analysis. Int. J. Environ. Res. Public Health 2018, 15, 2659. [Google Scholar] [CrossRef] [PubMed]
  16. Saaty, T.L. How to make a decision: The analytic hierarchy process. Eur. J. Oper. Res. 1990, 48, 9–26. [Google Scholar] [CrossRef]
  17. Saaty, T.L. Analytic hierarchy process. In Encyclopedia of Operations Research and Management Science; Springer: Boston, MA, USA, 2001. [Google Scholar]
  18. Tian, J.; Huo, Z.; Ma, F.; Gao, X.; Wu, Y. Application and Selection of Remediation Technology for OCPs-Contaminated Sites by Decision-Making Methods. Int. J. Environ. Res. Public Health 2019, 16, 1888. [Google Scholar] [CrossRef]
  19. Zhang, X.; Yang, J.; Zhao, X. Optimal Study of the Rural House Space Heating Systems Employing the AHP and FCE Methods. Energy 2018, 150, 631–641. [Google Scholar] [CrossRef]
  20. Li, S.Y.; Tao, G.; Zhang, L.J. Fire Risk Assessment of High-rise Buildings Based on Gray-FAHP Mathematical Model. Procedia Eng. 2018, 211, 395–402. [Google Scholar] [CrossRef]
  21. Liu, S.C. Network Security Risk Assessment Method Based on HMM and Attack Graph Model. In Proceedings of the 17th IEEE/ACIS International Conference on Software Engineering, Artificial Intelligence, Networking and Parallel/Distributed Computing (SNPD), Shanghai, China, 30 May–1 June 2016; pp. 517–522. [Google Scholar]
  22. Eboli, L.; Fu, Y.; Mazzulla, G. Multilevel Comprehensive Evaluation of the Railway Service Quality. Procedia Eng. 2016, 137, 21–30. [Google Scholar] [CrossRef]
  23. Lv, C.; Wu, Z.; Liu, Z.; Shi, L. The multi-level comprehensive safety evaluation for chemical production instalment based on the method that combines grey-clustering and EAHP. Int. J. Disaster Risk Reduct. 2017, 21, 243–250. [Google Scholar] [CrossRef]
  24. Li, W.; Liang, W.; Zhang, L.; Tang, Q. Performance assessment system of health, safety and environment based on experts’ weights and fuzzy comprehensive evaluation. J. Loss Prev. Process. Ind. 2015, 35, 95–103. [Google Scholar] [CrossRef]
  25. Wang, Y.; Zhao, N.; Jing, H.; Meng, B.; Yin, X. A Novel Model of the Ideal Point Method Coupled with Objective and Subjective Weighting Method for Evaluation of Surrounding Rock Stability. Math. Probl. Eng. 2016, 2016. [Google Scholar] [CrossRef]
  26. Yang, D.; Mak, C.M. An assessment model of classroom acoustical environment based on fuzzy comprehensive evaluation method. Appl. Acoust. 2017, 127, 292–296. [Google Scholar] [CrossRef]
  27. Jiang, L.; Wu, Q.; Wu, Q.; Wang, P.; Xue, Y.; Kong, P.; Gong, B. Fracture failure analysis of hard and thick key layer and its dynamic response characteristics. Eng. Fail. Anal. 2019, 98, 118–130. [Google Scholar] [CrossRef]
  28. Zhang, Z.; Wang, W.; Li, S.; Bai, J.; Hao, S.; Wu, H.; Yu, X. An Innovative Approach for Gob-Side Entry Retaining With Thick and Hard Roof: A Case Study. Teh. Vjesn. Tech. Gaz. 2018, 25, 1028–1036. [Google Scholar]
  29. Zhang, Z.; Yu, X.; Wu, H.; Deng, M. Stability Control for Gob-Side Entry Retaining with Supercritical Retained Entry Width in Thick Coal Seam Longwall Mining. Energies 2019, 12, 1375. [Google Scholar] [CrossRef]
Figure 1. The unit division in a three-layered structure.
Figure 1. The unit division in a three-layered structure.
Ijerph 16 02802 g001
Figure 2. Hierarchy structure model for the deep surrounding rock stability affecting factors.
Figure 2. Hierarchy structure model for the deep surrounding rock stability affecting factors.
Ijerph 16 02802 g002
Figure 3. Deep surrounding rock stability evaluation index system.
Figure 3. Deep surrounding rock stability evaluation index system.
Ijerph 16 02802 g003
Table 1. The grading standards of correlation indicators.
Table 1. The grading standards of correlation indicators.
Accident PossibilityScoreTime of ExposureScoreConsequenceScoreControl MeasuresScore
Absolutely impossible0.1Very rare0.5Minor consequences and no casualties1Preventive measures1
Extremely unlikely0.2Several times per year1Relatively serious consequences and injuries3Emergency measures3
Imaginable but highly unlikely0.5One time per month2Serious consequences and injuries7No control measures5
Far less likely1One time per week or accidentally3Very serious and one death15
Not often but likely3Exposure to hazardous conditions every day6At least two deaths40
Very likely6Continual exposure to hazardous conditions10Multiple fatalities100
Predictable10
Table 2. Scoring table with AHP.
Table 2. Scoring table with AHP.
AB1B2B3
B1132
B21/311/2
B31/221
B1C1C2C3C4C5
C114578
C21/41246
C31/51/2135
C41/71/41/312
C51/81/61/51/21
B2C6C7C8C9C10
C614517
C71/4121/42
C81/51/211/52
C914516
C101/71/21/21/61
B3C11C12C13C14C15
C1111/71/51/41/3
C1271234
C1351/2123
C1441/31/212
C1531/41/31/21
Table 3. Factors’ weights and ranks.
Table 3. Factors’ weights and ranks.
FactorsB1B2B3ARank
C10.5463--0.29481
C20.2131--0.11503
C30.1410--0.07615
C40.0609--0.03299
C50.0387--0.020911
C6-0.3856-0.06306
C7-0.1103-0.018012
C8-0.0767-0.012514
C9-0.3753-0.06137
C10-0.0520-0.008515
C11--0.04610.013713
C12--0.42710.12682
C13--0.26310.07814
C14--0.16220.04828
C15--0.10150.030110
Table 4. Parameter sequence of each factor and reference sequence.
Table 4. Parameter sequence of each factor and reference sequence.
FactorsAccident PossibilityTime of ExposureConsequenceControl Measures
Rock strength10101003
Rock mass quality indicator RQD66405
Coal seam dip angle36155
Structural surface strength coefficient1313
Volumetric joint count3315
Mining influence coefficient36155
Seepage quantity of groundwater1171
Geological structural influencing coefficient3233
Ground stress66401
Buried depth0.50.515
Cross-section shape of roadway3133
Sectional size of roadway63405
Thickness of strengthened rock36155
Support form and strength3673
Purposes and age limit of roadway3231
Reference value0.10.511
Table 5. Parameter sequence and reference sequence after normalization.
Table 5. Parameter sequence and reference sequence after normalization.
FactorsAccident PossibilityTime of ExposureConsequenceControl Measures
Rock strength2.75232.43905.15460.8491
Rock mass quality indicator RQD1.65141.46342.06191.4151
Coal seam dip angle0.82571.46340.77321.4151
Structural surface strength coefficient0.27520.73170.05160.8491
Volumetric joint count0.82570.73170.05161.4151
Mining influence coefficient0.82571.46340.77321.4151
Seepage quantity of groundwater0.27520.24390.36080.2830
Geological structural influencing coefficient0.82570.48780.15460.8491
Ground stress1.65141.46342.06190.2830
Buried depth0.13760.12200.05161.4151
Cross-section shape of roadway0.82570.24390.15460.8491
Sectional size of roadway1.65140.73172.06191.4151
Thickness of strengthened rock0.82571.46340.77321.4151
Support form and strength0.82571.46340.36080.8491
Purposes and age limit of roadway0.82570.48780.15460.2830
Reference value0.02750.12200.05160.2830
Table 6. Correlation coefficient and sorting of correlation degree.
Table 6. Correlation coefficient and sorting of correlation degree.
FactorsAccident PossibilityTime of ExposureConsequenceControl MeasuresCorrelation DegreeRank
Rock strength0.48360.52410.33330.81840.53991
Rock mass quality indicator RQD0.61110.65540.55930.69270.62962
Coal seam dip angle0.76170.65540.77950.69270.72237
Structural surface strength coefficient0.91150.87071.00000.81840.884312
Volumetric joint count0.76170.80711.00000.69270.81549
Mining influence coefficient0.76170.65540.77950.69270.72237
Seepage quantity of groundwater0.91150.95440.89191.00000.939415
Geological structural influencing coefficient0.76170.87460.96120.81840.854010
Ground stress0.61110.65540.55931.00000.70654
Buried depth0.95861.00001.00000.69270.912814
Cross-section shape of roadway0.76170.95440.96120.81840.873911
Sectional size of roadway0.61110.80710.55930.69270.66763
Thickness of strengthened rock0.76170.65540.77950.69270.72237
Support form and strength0.76170.65540.89190.81840.78198
Purposes and age limit of roadway0.76170.87460.96121.00000.899413
Table 7. The normalized weights of key factors.
Table 7. The normalized weights of key factors.
Index LayerKey FactorsWeightsNormalized Weights
Geological featuresRock strength0.29480.6067
Rock mass quality indicator RQD0.11500.2367
Coal seam dip angle0.07610.1566
Sum0.48591
Environmental features Mining influence coefficient0.06300.5068
Ground stress0.06130.4932
Sum0.12431
Engineering featuresSectional size of roadway0.12680.6188
Thickness of strengthened rock0.07810.3812
Sum0.20491

Share and Cite

MDPI and ACS Style

Huang, D.; Li, W.; Chang, X.; Tan, Y. Key Factors Identification and Risk Assessment for the Stability of Deep Surrounding Rock in Coal Roadway. Int. J. Environ. Res. Public Health 2019, 16, 2802. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16152802

AMA Style

Huang D, Li W, Chang X, Tan Y. Key Factors Identification and Risk Assessment for the Stability of Deep Surrounding Rock in Coal Roadway. International Journal of Environmental Research and Public Health. 2019; 16(15):2802. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16152802

Chicago/Turabian Style

Huang, Dongmei, Weijun Li, Xikun Chang, and Yunliang Tan. 2019. "Key Factors Identification and Risk Assessment for the Stability of Deep Surrounding Rock in Coal Roadway" International Journal of Environmental Research and Public Health 16, no. 15: 2802. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16152802

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop