Quantitative Risk Assessment for Deep Tunnel Failure Based on Normal Cloud Model: A Case Study at the ASHELE Copper Mine, China

The spatial and temporal distribution of tunnel failure is very complex due to geologic heterogeneity and variability in both mining processes and tunnel arrangement in deep metal mines. In this paper, the quantitative risk assessment for deep tunnel failure was performed using a normal cloud model at the Ashele copper mine, China. This was completed by considering the evaluation indexes of geological condition, mining process, and microseismic data. A weighted distribution of evaluation indexes was determined by implementation of an entropy weight method to reveal the primary parameters controlling tunnel failure. Additionally, the damage levels of the tunnel were quantitatively assigned by computing the degree of membership that different damage levels had, based on the expectation normalization method. The methods of maximum membership principle, comprehensive evaluation value, and fuzzy entropy were considered to determine the tunnel damage levels and risk of occurrence. The application of this method at the Ashele copper mine demonstrates that it meets the requirement of risk assessment for deep tunnel failure and can provide a basis for large-scale regional tunnel failure control in deep metal mines.


Introduction
At the current rate of mine expansion, more than 30% of the metal mines in China will enter into the deep mining category (i.e., a mining depth greater than 1000 m) by the end of 2025. This indicates that the exploitation of deep mineral resources will become commonplace in China [1][2][3]. Under high ground stress conditions in deep metal mines, frequent and strong blast disturbances can lead to the spalling and collapse of the rock surrounding tunnels. For example, more than 20 m of the roadway in the Hongtoushan copper mine, one of the deepest nonferrous metal mines in China, was damaged by blasting disturbance at a mining depth of 1000 m [4]. In the Dongguashan copper mine, large-scale collapses in deep tunnels were induced by high ground stress and frequent blasting. In the Jinchuan No. 2 mine, the combined effects of high geo-stress and horizontal tectonic stress resulted in the large-scale deformation and failure of the rock mass, and more than 56% of the deep tunnels have serious damage from collapse, where the maximum convergent deformation can reach 500 mm [5]. In the Ashele copper mine, tunnel collapse and the burial of drilling equipment lead to a long interruption of production. As such, a quantitative risk assessment for deep tunnel failure should be carried out to ensure the safety of workers and equipment (e.g., targeted support measures and systematic mining process optimization).
The deep metal mines have many unique characteristics. First, the mining system is very complex, with the tunnel and shaft arranged in three-dimensional space that can

Entropy Weight Method
In the application of the normal cloud model, the different evaluation indexes have contributions to the evaluation object. The weight, expressed in a quantitative value that reflects the relative importance of each evaluation index in the set of factors, directly affects the resulting comprehensive evaluation. As such, scientific and reasonable assessment of the evaluation index weights is very important for building an objective evaluation model.
The entropy weight method can be used to weigh multiple evaluation indicators and does not depend on the subjective judgment or experience of decision-makers, which is more objective and accurate for weight calculation. Moreover, the calculation process of this method is simple, and it easily solves the weighting problem in multi-index target evaluation scenarios. The entropy weight method assigns weights according to degree of difference of each evaluation index. If there is a large difference in a specific index, it means that the index contains a lot of information and has a great influence on the comprehensive evaluation results. Thus, it is given a larger weight value. There are three main steps of the entropy weight method, as follows: (1) Dimensionless processing of the data. According to the property of each index, the profitability index and the cost index can be divided as: maxx ij (i = 1, 2, . . . , n) the profitability index minx ij x ij (i = 1, 2, . . . , n) the cos t index (1) where j = 1, 2, . . . , m; n is the number of data, and m is the number of evaluation indicators. The larger the profitability index and the smaller the cost index, the better they are. Then, the data can be normalized as: (2) Calculate the information entropy of each indicator. The information entropy of the j-th evaluation index is expressed as: u ij ln u ij , (j = 1, 2, . . . , m; i = 1, 2, . . . , n) where , l = 1 ln n , and u ij = 0, u ij lnu ij = 0.
(3) Calculate the weight. The weight of each evaluation index is calculated as:

Cloud Model 2.2.1. Definition and Numerical Characteristics of the Cloud Model
The cloud model considers the uncertainty of converting between qualitative concepts and quantitative numerical representations [47]. Let U be a quantitative set, then U = {x}. C is a qualitative concept in U, and the definite parameter x∈U is a random occurrence in C. For any element x of U, the membership degree of x in C is u(x)∈[0,1], and u(x) is a stable random number. The distribution of x in U is called the cloud, and every x is called a Appl. Sci. 2021, 11, 5208 4 of 19 cloud droplet. A cloud that is composed of a large number of cloud droplets can represent a qualitative concept in the domain space.
Three numerical characteristics (Ex, En, and He) are introduced to express the qualitative concept of cloud theory. Expectation (Ex) is the expectation of the spatial distribution of cloud droplets in the domain space and the mean value of the set. Ex is also the point that best represents the qualitative concept. Entropy (En) is determined by the randomness and fuzziness of the qualitative concept. Specifically, En is a measure of the randomness and fuzziness of a qualitative concept. Hyper entropy (He) is a measure of the uncertainty of the entropy and reflects the cohesion of the uncertainty at all points in the domain space. The value of hyper entropy He indirectly reflects the thickness of cloud droplets.

Cloud Generator
Cloud models are executed by cloud generators. There are four commonly used cloud generators: the forward cloud generator, the backward cloud generator, the twodimensional cloud generator, and the x-condition cloud generator ( Figure 1). The forward cloud generator (CG) is used to generate the cloud drops with a membership degree that conforms to the three numerical characteristics. On the contrary, in the backward generator (CG −1 ), the transferring process is used to derive the qualitative concept represented by the three descriptors from cloud drops. The two-dimensional cloud generators can establish a three-dimensional cloud droplet (x 1 , x 2 , y) that satisfies the two-dimensional cloud distribution on the basis of the one-dimensional cloud generator. Lastly, the x-condition cloud generator (CG x ) uses the three numerical characteristics of cloud and the specified x = x 0 to generate the cloud drops.

Definition and Numerical Characteristics of the Cloud Model
The cloud model considers the uncertainty of converting between qualitative concepts and quantitative numerical representations [47]. Let U be a quantitative set, then U = {x}. C is a qualitative concept in U, and the definite parameter x∈U is a random occurrence in C. For any element x of U, the membership degree of x in C is u(x)∈[0,1], and u(x) is a stable random number. The distribution of x in U is called the cloud, and every x is called a cloud droplet. A cloud that is composed of a large number of cloud droplets can represent a qualitative concept in the domain space.
[ ] ( ) Three numerical characteristics (Ex, En, and He) are introduced to express the qualitative concept of cloud theory. Expectation (Ex) is the expectation of the spatial distribution of cloud droplets in the domain space and the mean value of the set. Ex is also the point that best represents the qualitative concept. Entropy (En) is determined by the randomness and fuzziness of the qualitative concept. Specifically, En is a measure of the randomness and fuzziness of a qualitative concept. Hyper entropy (He) is a measure of the uncertainty of the entropy and reflects the cohesion of the uncertainty at all points in the domain space. The value of hyper entropy He indirectly reflects the thickness of cloud droplets.

Cloud Generator
Cloud models are executed by cloud generators. There are four commonly used cloud generators: the forward cloud generator, the backward cloud generator, the twodimensional cloud generator, and the x-condition cloud generator ( Figure 1). The forward cloud generator (CG) is used to generate the cloud drops with a membership degree that conforms to the three numerical characteristics. On the contrary, in the backward generator (CG −1 ), the transferring process is used to derive the qualitative concept represented by the three descriptors from cloud drops. The two-dimensional cloud generators can establish a three-dimensional cloud droplet (x1, x2, y) that satisfies the two-dimensional cloud distribution on the basis of the one-dimensional cloud generator. Lastly, the x-condition cloud generator (CGx) uses the three numerical characteristics of cloud and the specified x = x0 to generate the cloud drops. By analyzing the determination of membership function in both natural and social sciences, Li et al. [48] determined that the normal distribution is the most suitable distri- By analyzing the determination of membership function in both natural and social sciences, Li et al. [48] determined that the normal distribution is the most suitable distribution for the membership function. According to the central-limit theorem, when a random variable is affected by many small and independent random factors, and each factor does not play an overwhelming role in the controlling distribution, this random variable has a normal distribution. In the case of tunnel collapse, there are many factors that can lead to tunnel damage, and each factor does not play a leading role. Therefore, the normal distribution cloud generator is used for the quantitative risk assessment of deep tunnel failure.
where w j is the weight of the j-th evaluation index, j = 1, 2, . . . , m.

Calculation Process of Cloud Model
In this paper, the comprehensive membership degree and comprehensive evaluation value of sample data are calculated based on the cloud model (shown in Figure 2). The calculation process includes three steps.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 6 of 21 where Pi is the comprehensive evaluation value of the i-th sample data and vk is the eigenvalue of the k-th damage level.

The Index System of Risk Assessment for Tunnel Failure
The evaluation index system is the foundation of risk assessment for tunnel failure. In this paper, the classification standards of tunnel failure are determined based on the investigation of deep tunnel failure in the Ashele copper mine. The evaluation indexes used here were derived by analyzing the index's reliability, weight value, and characteristic value.

Damage Level of the Tunnel
The tunnel damage level was determined by considering the existing standards [50][51][52] and the results of site surveys. Four damage levels were created: no damage, slight damage, moderate damage, and serious damage (Table 1). (1) The numerical characteristics of the cloud associated with each evaluation index are calculated by using the backward cloud generator.
(2) Based on x-conditional cloud model, the membership degree u(x ij ) of sample data belonging to a certain damage level is calculated. After the weight of each evaluation index, w j is determined by the entropy weight method, and the comprehensive membership degree b i of sample data at certain damage levels is expressed as: The damage level of the sample data is then obtained according to the maximum membership rule. However, when the comprehensive membership degree of different damage levels is similar, using the maximum membership rule may provide incorrect values. The fuzzy entropy, H i , can provide additional information in these cases and reflect the similarity of the membership degree [49]. The fuzzy entropy term is defined as: where b ik is the comprehensive membership degree of the k-th damage level, s is the number of damage levels, and f is the normalization coefficient.
(3) The integrated cloud model is established according to Equations (6)- (8). After the eigenvalues of damage levels are determined from the expectation, Ex, the comprehensive evaluation value p of the sample data can be calculated based on the weighted mean method as: where P i is the comprehensive evaluation value of the i-th sample data and v k is the eigenvalue of the k-th damage level.

The Index System of Risk Assessment for Tunnel Failure
The evaluation index system is the foundation of risk assessment for tunnel failure. In this paper, the classification standards of tunnel failure are determined based on the investigation of deep tunnel failure in the Ashele copper mine. The evaluation indexes used here were derived by analyzing the index's reliability, weight value, and characteristic value.

Damage Level of the Tunnel
The tunnel damage level was determined by considering the existing standards [50][51][52] and the results of site surveys. Four damage levels were created: no damage, slight damage, moderate damage, and serious damage (Table 1). Table 1. Description of the different damage levels of the tunnel.

Damage Level Description
No damage There is no damage or only slight deformation in the tunnel. The current support is able to keep the rock mass stable.

Slight damage
The deformation of rock mass in the tunnel can be observed by the naked eye. A small range of rock mass has fallen or is bulking. Parts of the support appear damaged but can keep being used after small repairs.

Moderate damage
The support and rock mass in the tunnel appear to be damaged and fall in a range greater than a meter in length. Secondary support should be put in place.

Serious damage
The integrity of the support and rock mass in tunnel is severely damaged along several meters and has had serious impact on the normal production work. It is difficult to clean up the collapsed rock. Secondary reinforcement or support should be put in place.
By systematic investigation of tunnel damage in the Ashele copper mine, it was found that there was a variety of damage in the tunnel in the 350 m level to the 0 m level (i.e., the mining depth was approximately from 565 m to 915 m), including spalling and bulking of the sidewall or the roof caving (Figures 3 and 4). The majority of the damage belonged to the slight damage and moderate damage index; however, slight damage accounted for the highest proportion. There were four typical characteristics of the tunnel damage in the Ashele copper mine. First, damage degree increased with buried depth. For example, there was a large range of tunnel damage in strike-drift at the 50 m level. Second, tunnel damage is closely related to the maximum principal stress direction. The majority of the damage occurred in the strike-drift that was oriented perpendicular to the maximum principal stress, rather than in the transverse drift, which was oriented parallel to the maximum principal stress. Third, tunnel damage was greatly affected by mining disturbance. For example, some transverse drifts were located close to the stope and are prone to damage. Fourth, the lithology and the development of joints and fractures had an important influence on the damage of tunnel. For example, tunnels located in rock with significant pyrite developed joints and fractures in the hanging wall and are prone to damage.
By systematic investigation of tunnel damage in the Ashele copper mine, it was found that there was a variety of damage in the tunnel in the 350 m level to the 0 m level (i.e., the mining depth was approximately from 565 m to 915 m), including spalling and bulking of the sidewall or the roof caving (Figures 3 and 4). The majority of the damage belonged to the slight damage and moderate damage index; however, slight damage accounted for the highest proportion. There were four typical characteristics of the tunnel damage in the Ashele copper mine. First, damage degree increased with buried depth. For example, there was a large range of tunnel damage in strike-drift at the 50 m level. Second, tunnel damage is closely related to the maximum principal stress direction. The majority of the damage occurred in the strike-drift that was oriented perpendicular to the maximum principal stress, rather than in the transverse drift, which was oriented parallel to the maximum principal stress. Third, tunnel damage was greatly affected by mining disturbance. For example, some transverse drifts were located close to the stope and are prone to damage. Fourth, the lithology and the development of joints and fractures had an important influence on the damage of tunnel. For example, tunnels located in rock with significant pyrite developed joints and fractures in the hanging wall and are prone to damage.

Evaluation Index
Three types of evaluation indexes (geological parameters, mining disturbances, and MS data) ( Figure 5) were selected for the risk assessment of tunnel failure in the Ashele copper mine by considering the mining engineering, monitoring data, and tunnel damage characteristics. The evaluation indexes of geological parameters were angled between the tunnel and the maximum principal stress, depth, and lithology. The evaluation index of mining disturbances was the distance between the tunnel and its nearest stope. The evaluation indexes of MS data included the accumulative number of MS events, the accumu-

Evaluation Index
Three types of evaluation indexes (geological parameters, mining disturbances, and MS data) ( copper mine by considering the mining engineering, monitoring data, and tunnel damage characteristics. The evaluation indexes of geological parameters were angled between the tunnel and the maximum principal stress, depth, and lithology. The evaluation index of mining disturbances was the distance between the tunnel and its nearest stope. The evaluation indexes of MS data included the accumulative number of MS events, the accumulative released energy, the density of MS events, apparent stress, and displacement. The definitions of these nine evaluation indexes are introduced in Table 2.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 9 of 21 Accumulated energy Total energy released by MS events within a certain time period and region. The accumulated energy is used to reflect the degree of energy released by cracks generated inside the rock mass.
Event density Number of MS events per unit volume of rock mass. The event density is used to describe the cluster extent of MS events.

Apparent stress σa
The apparent stress is the ratio of the total radiated seismic energy to the seismic moment. It assesses the amount of energy released per unit of deformation and is defined as: σa = μ/M0, where u is the modulus of rigidity of the source and M0 is the seismic moment.

Displacement u
It is the average displacement of the source, where r0 is the source radius and is defined as: In total, 100 observations between April 2017 to June 2018 were used to populate a risk assessment model for tunnel failure ( Table 3). The numbers of no damaged samples, slightly damaged samples, and moderately damaged samples were 24, 44, and 22, respectively. In total, 90 of the observations were randomly selected for modeling, while 10 observations were used for the generalization ability test.   Table 2. Description of the evaluation indexes for tunnel failure in the Ashele copper mine.

Evaluation Index Description
Angle Defined as the angle between the tunnel and the maximum principal stress. The majority of the damage occurred in the strike-drift that, which is oriented perpendicular to the maximum principal stress.

Depth
The deeper the tunnels are buried, the greater the vertical stress and the more serious damage that occurs.

Lithology
The lithology and the development of joints and fractures have an important influence on the damage of tunnel.

Distance
Rock mass is inevitably subjected to changes in stress field after the ore is mined and during blasting disturbances. The closer the tunnel is to the stope, the more the mining disturbances and unloading effects affect the tunnel stability.
Accumulated events Total number of MS events generated within a certain time period and region. The number of accumulated events is used to evaluate the change in regional fracture activities.
Accumulated energy Total energy released by MS events within a certain time period and region. The accumulated energy is used to reflect the degree of energy released by cracks generated inside the rock mass.
Event density Number of MS events per unit volume of rock mass. The event density is used to describe the cluster extent of MS events.

Apparent stress σ a
The apparent stress is the ratio of the total radiated seismic energy to the seismic moment. It assesses the amount of energy released per unit of deformation and is defined as: σ a = µ/M 0 , where u is the modulus of rigidity of the source and M 0 is the seismic moment. In total, 100 observations between April 2017 to June 2018 were used to populate a risk assessment model for tunnel failure ( Table 3). The numbers of no damaged samples, slightly damaged samples, and moderately damaged samples were 24, 44, and 22, respectively. In total, 90 of the observations were randomly selected for modeling, while 10 observations were used for the generalization ability test.

Correlation Between Evaluation Indicators and Damage Levels
The correlation between the damage levels and the nine evaluation indexes was analyzed to verify the reliability of the evaluation indexes. Table 4 displays the number of angles and lithology occurrences at different damage levels. It was observed that more slight damage observations occurred in tuff and the moderate damage observations were prone to occur in pyrite. In addition, 77 damage observations occured where the tunnel was oriented perpendicular to the maximum principal stress (angle = 90 • ), and only 13 damage observations occurred when the tunnel was oriented parallel to the maximum principal stress (angle = 0 • ). This indicates that tunnel orientation has a strong correlation to tunnel damage.  Figure 6 displays the distribution of the other seven evaluation indexes, including the accumulated MS events, the accumulated released energy from the MS events, the MS events density, the displacement, the apparent stress, distance between the tunnel and its nearest stope, and the buried depth. There are many outliers within the upper part of the boxplot for accumulated MS events and the accumulated energy. This may indicate that some sample data of the two indexes deviate from the data distribution center. If these data were used in the tunnel damage risk assessment, the accuracy of the results would be significantly affected. Due to this dispersion, the uppermost and lowermost 5% of the data were removed from each of the seven indexes prior to implementation in tunnel damage risk assessment. The data distribution of these seven evaluation indexes and their mean values after truncation are shown in Table 5. The boundaries of the statistical parameters of each evaluation index overlap between the three damage levels. However, from the scatter diagram of the distribution of evaluation indexes with different damage levels (Figure 7), it can be seen that, although the distribution of the MS events, MS events density, displacement, the apparent stress overlap in each damage category to some extent, each damage level has a distinct cluster of observations. The boundaries of distributions of the other three indexes (the angle, distance, and buried depth) are indistinct at the three damage levels. Because of this, it is not possible to conduct the tunnel damage risk assessment by using a single evaluation index, and a comprehensive analysis by using multi-indexes is necessary to improve the accuracy of the evaluation results.

Weight of Evaluation Indexes
The weight of the nine evaluation indexes was calculated by using the entropy weight method, according to Equations (1)-(4). Any single index does not play a dominant role in risk assessment compared to other indexes ( Figure 8). Thus, multi-indexes models should be applied during the risk assessment of tunnel damage. Among the nine indexes, the accumulated energy and accumulated events have the largest weights (0.211 and 0.21, respectively) and have the greatest impact on the results of tunnel failure risk assessment. The weight of the buried depth index is smallest, which is consistent with the scattered distribution and indistinct boundary of this index at different damage levels (Figure 7). It is worth noting that the evaluation indexes with the five largest weights are MS parameters and account for a cumulative 83.8% of the total weight. This indicates that the MS parameters play a leading role in the risk assessment of tunnel failure.

Weight of Evaluation Indexes
The weight of the nine evaluation indexes was calculated by using the entropy weight method, according to Equations (1)-(4). Any single index does not play a dominant role in risk assessment compared to other indexes ( Figure 8). Thus, multi-indexes models should be applied during the risk assessment of tunnel damage. Among the nine indexes, the accumulated energy and accumulated events have the largest weights (0.211 and 0.21, respectively) and have the greatest impact on the results of tunnel failure risk assessment. The weight of the buried depth index is smallest, which is consistent with the scattered distribution and indistinct boundary of this index at different damage levels ( Figure 7). It is worth noting that the evaluation indexes with the five largest weights are MS parameters and account for a cumulative 83.8% of the total weight. This indicates that the MS parameters play a leading role in the risk assessment of tunnel failure.

Weight of Evaluation Indexes
The weight of the nine evaluation indexes was calculated by using the entropy weight method, according to Equations (1)-(4). Any single index does not play a dominant role in risk assessment compared to other indexes ( Figure 8). Thus, multi-indexes models should be applied during the risk assessment of tunnel damage. Among the nine indexes, the accumulated energy and accumulated events have the largest weights (0.211 and 0.21, respectively) and have the greatest impact on the results of tunnel failure risk assessment. The weight of the buried depth index is smallest, which is consistent with the scattered distribution and indistinct boundary of this index at different damage levels (Figure 7). It is worth noting that the evaluation indexes with the five largest weights are MS parameters and account for a cumulative 83.8% of the total weight. This indicates that the MS parameters play a leading role in the risk assessment of tunnel failure.

Numerical Characteristics of Evaluation Indexes
The numerical characteristics of the evaluation indexes were obtained by using a backward cloud generator (Table 6). Based on the numerical characteristics of the evaluation indexes, the membership relationships between each evaluation index and different damage levels were obtained by the two-dimensional cloud generator, in which each cloud contains 10,000 cloud drops. The results were visualized through clustering diagrams (Figure 9). The boundaries of the cloud drops that represent the buried depth and angle indexes at different damage levels are not clear, while boundaries of the cloud drops of the accumulated events and accumulated energy indexes are distinct. This result is in agreement with the weight value of the evaluation indexes ( Figure 8). cloud contains 10,000 cloud drops. The results were visualized through clustering diagrams ( Figure 9). The boundaries of the cloud drops that represent the buried depth and angle indexes at different damage levels are not clear, while boundaries of the cloud drops of the accumulated events and accumulated energy indexes are distinct. This result is in agreement with the weight value of the evaluation indexes ( Figure 8).

Eigenvalues of Damage Levels
During tunnel damage risk assessment, using a weighted mean method, it is necessary to determine the eigenvalues of the different damage levels. The expectation (Ex) is the best parameter to quantitatively describe the damage levels (i.e., it can reflect the difference between the damage levels). Ex was used to determine the eigenvalues of damage 0 Figure 9. Clustering of the evaluation indexes of different damage levels generated by the cloud models.

Eigenvalues of Damage Levels
During tunnel damage risk assessment, using a weighted mean method, it is necessary to determine the eigenvalues of the different damage levels. The expectation (Ex) is the best parameter to quantitatively describe the damage levels (i.e., it can reflect the difference between the damage levels). Ex was used to determine the eigenvalues of damage levels to avoid subjectivity in using the weighted mean method. When the weight of the evaluation indexes were obtained, the numerical characteristics can be calculated by applying the integrated cloud model ( Table 7). The membership degrees of the three damage levels are displayed in Figure 10. In this paper, the Ex of the three damage levels were proceeded by normalization. The Ex of the three levels is 15.46, 22.17, and 38.03 for no damage, slight damage, and moderate damage, respectively. After normalization, the eigenvalues of the three damage levels are 1.00, 1.43, and 2.46, respectively. levels to avoid subjectivity in using the weighted mean method. When the weight of the evaluation indexes were obtained, the numerical characteristics can be calculated by applying the integrated cloud model ( Table 7). The membership degrees of the three damage levels are displayed in Figure 10. In this paper, the Ex of the three damage levels were proceeded by normalization. The Ex of the three levels is 15.46, 22.17, and 38.03 for no damage, slight damage, and moderate damage, respectively. After normalization, the eigenvalues of the three damage levels are 1.00, 1.43, and 2.46, respectively.

Determination of the Damage Levels
Based on the weights of the evaluation indexes and the degree of membership obtained by the x-condition cloud generator, the comprehensive membership degree B, fuzzy entropy H, and comprehensive evaluation value p of the 100 data samples were calculated according to Equations (9)-(12) ( Table 8). When the maximum membership rule was used to determine damage levels, portions of the evaluation information might be lost. Thus, the comprehensive evaluation value p was used to determine the damage level by using the weighted mean method. Moreover, the correctness of the eigenvalues for each damage level calculated by normalized expectations was verified ( Figure 11). The comprehensive evaluation results were interrogated by using the two assignment methods for the damage levels: the first method defined the eigenvalues of different damage levels as 1, 2, and 3, while the second method defined them as 1.00, 1.43, and 2.46. The results indicate that the distribution of the later method is smaller, which verifies that using the normalized expectations to define the eigenvalues of the three damage levels was more accurate. The average value of the comprehensive evaluation values for each damage level, P

Determination of the Damage Levels
Based on the weights of the evaluation indexes and the degree of membership obtained by the x-condition cloud generator, the comprehensive membership degree B, fuzzy entropy H, and comprehensive evaluation value p of the 100 data samples were calculated according to Equations (9)-(12) ( Table 8). When the maximum membership rule was used to determine damage levels, portions of the evaluation information might be lost. Thus, the comprehensive evaluation value p was used to determine the damage level by using the weighted mean method. Moreover, the correctness of the eigenvalues for each damage level calculated by normalized expectations was verified ( Figure 11). The comprehensive evaluation results were interrogated by using the two assignment methods for the damage levels: the first method defined the eigenvalues of different damage levels as 1, 2, and 3, while the second method defined them as 1.00, 1.43, and 2.46. The results indicate that the distribution of the later method is smaller, which verifies that using the normalized expectations to define the eigenvalues of the three damage levels was more accurate. The average value of the comprehensive evaluation values for each damage level, P No , P Slight , and P Moderate , were 1.50, 1.74, and 1.87, respectively.
Comparing the known damage levels of 100 data samples to the predicted damage levels, it was found that the accuracy rate of damage level determined based on the principle of maximum affiliation was 81.1% (Table 9). The accuracy rate of slight damage was the lowest, with 13 of 44 sample data misclassified as moderate damage, which lead to a lower accuracy rate of the model. It is worth noting that, for the majority of samples that were misclassified, the comprehensive membership degrees of the correct damage level and the chosen damage level were very close. For example, sample data No. 7 was identified as slight damage level in the mine but was classified as moderate damage. In the comprehensive membership, degrees of the slight damage and moderate damage were 0.39 and 0.40. According to the maximum membership rule, the evaluation result of this sample data is moderate damage. However, the comprehensive evaluation value p of this sample is 1.75, which is very close to the mean value of slight damage samples, P Slight = 1.74. In this case, classification of the damage level based on the comprehensive evaluation value p was more correct. By analyzing samples that were misclassified by the maximum membership rule, it was found that the fuzzy entropy of each was larger than 0.89. Therefore, the comprehensive evaluation value p was used to recalculate the damage level for misclassified samples when the fuzzy entropy was larger than 0.89. This caused four sample data that were originally misclassified to be reclassified as the correct damage level. This additional classification improved the accuracy rate of total sample data improved to 86%. Thus, the implementation of the maximum membership rule and the comprehensive evaluation value, where the maximum membership rule was used, when fuzzy entropy H < 0.89, and the comprehensive evaluation value is used when fuzzy entropy H > 0.89, can improve the evaluation accuracy of tunnel damage level.   Comparing the known damage levels of 100 data samples to the predicted damage levels, it was found that the accuracy rate of damage level determined based on the principle of maximum affiliation was 81.1% (Table 9). The accuracy rate of slight damage was the lowest, with 13 of 44 sample data misclassified as moderate damage, which lead to a

Tunnel Failure Risk Assessment Method
To ensure the safety of mining process, it is necessary to establish a risk assessment method of tunnel damage. Using the inputs and methods of the normal cloud model above, a systematic assessment process was formulated that fit the mining conditions in the Ashele copper mine ( Figure 12). First, the damage level was calculated according to the maximum membership rule and the comprehensive evaluation value. When the two results were consistent, the damage levels were assigned to risk levels (no damage = no risk, slight damage = low risk, moderate damage = medium risk, and serious damage = high risk). Second, if the two results were inconsistent, the evaluation threshold T of the fuzzy entropy was determined by comparing the actual damage level and the result obtained from the maximum membership rule. The tunnel failure risk was then retrieved from an evaluation matrix ( Figure 13). Third, in the process of the application of the evaluation matrix, two cases should be considered. When the fuzzy entropy H was larger than threshold T, that is the comprehensive membership degrees of different damage levels were very close, the higher damage level was output as the result according to the conservative principle. Otherwise, when the fuzzy entropy H was smaller than threshold T, the result depended on the comprehensive evaluation value.

Application
Using the assessment method, regular risk assessment for tunnel failure was carried out every 6 months in the Ashele copper mine. Using the model created above, the average values of the comprehensive evaluation values of the three damage levels P No , P Slight , and P Moderate were 1.50, 1.74, and 1.87, respectively, and the threshold of fuzzy entropy was 0.89. The risk assessment results are visually displayed on the 3D geological model, which was used for the optimization of the mining process and support.
Results of the risk assessment performed in June 2019 show that the medium risk areas were dominantly concentrated at the footwall at the 150 m level and 100 m level (Figure 14), accounting for 2% of the total length of the tunnels. Low risk areas were mainly located in four regions: the footwall of the 200 m, 150 m, 100 m, and 50 m levels and the hang wall of the 150 m and 100 m levels. The length of these tunnels accounted for 9% of the total tunnel length. The modeled results were verified by site investigations, and secondary reinforcements were put in place through areas of medium risk. This integrated model and site investigation workflow effectively reduced the risk of ground pressure hazards during deep mining in the Ashele copper mine.

Application
Using the assessment method, regular risk assessment for tunnel failure was carried out every 6 months in the Ashele copper mine. Using the model created above, the average values of the comprehensive evaluation values of the three damage levels P -No, P -Slight, and P -Moderate were 1.50, 1.74, and 1.87, respectively, and the threshold of fuzzy entropy was 0.89. The risk assessment results are visually displayed on the 3D geological model, which was used for the optimization of the mining process and support.
Results of the risk assessment performed in June 2019 show that the medium risk areas were dominantly concentrated at the footwall at the 150 m level and 100 m level (Figure 14), accounting for 2% of the total length of the tunnels. Low risk areas were mainly located in four regions: the footwall of the 200 m, 150 m, 100 m, and 50 m levels and the hang wall of the 150 m and 100 m levels. The length of these tunnels accounted for 9% of the total tunnel length. The modeled results were verified by site investigations,

Conclusions
Rock mass failure in deep metal mines is complex and uncertain in terms of time, space, frequency, and intensity. Considering the complex characteristics of tunnel failure in the Ashele copper mine, a quantitative method to assess risk of tunnel failure was carried out by considering three aspects of tunnel failure (geological conditions, mining dis-