A Normal Cloud Model-Based Method for Risk Assessment of Water Inrush and Its Application in a Super-Long Tunnel Constructed by a Tunnel Boring Machine in the Arid Area of Northwest China

: Tunnel water inrush is complex, fuzzy, and random, and it is a ﬀ ected by many factors, such as hydrology, geology, and construction. However, few papers have considered the impact of dynamic monitoring on water inrush in previous research. In this study, considering geological, hydrological, and construction factors, as well as dynamic monitoring, a new multi-index evaluation method is proposed to analyze the risk of tunnel water inrush based on the normal cloud model. A new weight algorithm combining analytic hierarchy process and entropy method is used to calculate the index weight. The certainty degree of each evaluation index belonging to the corresponding cloud can be obtained by the cloud model theory. The ﬁnal level of tunnel water inrush is determined via the synthetic certainty degree. The proposed method is applied to analyze the risk of water inrush in the SS (Shuang-san) tunnel constructed by a tunnel boring machine in the arid area of Northwest China. The evaluation results are not only basically identical to the results calculated by the ideal point and gray relation projection methods, but also agree well with the actual excavation results. This demonstrates that this new risk assessment method of water inrush has high accuracy and feasibility. Simultaneously, it also provides a new research idea to analyze the probability of tunnel water inrush and can provide a reference for related projects.


Introduction
Water inrush has great perniciousness in tunnel excavation [1,2]; particularly for a tunnel constructed by a tunnel boring machine (TBM), it may not only submerge the mechanical equipment of the TBM, which is worth tens of millions of dollars, but also seriously endanger the safety and life of the builders in the tunnel [3,4]. For the construction of tunnels in arid areas, water inrush also leads to a decline in groundwater level, surface vegetation degradation, and land desertification [5]. Therefore, predicting the risk of water inrush and reducing its harm are of great significance for regional sustainable development.
Water inrush in tunnels is affected by many factors, such as hydrology, geology, and construction, and it is complex, fuzzy, and random. In recent years, many researchers and engineers have made great efforts to assess the risk of tunnel water inrush, and have developed various methods of analysis, including the geological analysis method [6], the analytic hierarchy process (AHP) method [7,8], the neural network method [9,10], the geographic information system (GIS) method [11,12], and the fuzzy synthetic evaluation method [13,14]. Although these methods provide a theoretical basis to analyze the risk of water inrush, there still have some shortcomings. For example, due to the limitation of previous exploration technology, the geological analysis method cannot accurately obtain the hydrogeological conditions along the tunnel. The neural network method needs a large amount of historical data to construct learning samples in the evaluation process. The GIS method is mainly based on topography, geomorphology, and other factors to determine the risk of water inrush, which ignores the complexity of water inrush. Although the fuzziness and randomness of water inrush are taken into account in the above methods, few scholars have considered the impact of dynamic monitoring on water inrush in the evaluation process. In this study, based on detailed geological, hydrological, construction, and dynamic monitoring data of the SS (Shuang-san) tunnel, which is a super-long tunnel constructed by a TBM in the arid area of Northwest China, a new multi-index evaluation method is proposed to analyze the risk of tunnel water inrush based on the normal cloud model. The aim of this method is to further improve the accuracy of the risk prediction of water inrush.
The structure of this paper is designed as follows: Section 2 introduces the engineering background of the study area. Section 3 proposes a new risk assessment method of water inrush based on the cloud model and the comprehensive weight algorithm. Section 4 discusses the evaluation results of the proposed cloud model method by comparing them with other methods. Section 5 verifies that the actual excavation situations are consistent with the evaluation results. Finally, Section 6 summarizes the research of this paper.

Engineering Background
The SS tunnel is located at the junction of the northern region of the Tianshan Mountains and the Jungger Basin in Xinjiang Province, China, as shown in Figure 1. It is a typical long tunnel with a length of 31.5 km and a longitudinal slope of 1/5000, 22.7 km of which was constructed by a TBM. The excavation diameter is 5.5 m, the ground elevation along the line is 775-1005 m, the topography is not undulating, the relative height difference is 5-35 m and the overburden thickness is 80-295 m. The geological structure of the SS tunnel mainly consists of Tuffaceous sandstone and lithic tuff from SD52+500 to SD50+300, small-scale faults and joint fissures are well developed, and the width of the fault zones is 0.3-5 m. The geological section of the SS tunnel is presented as Figure 2. The surface water above the tunnel is not abundant, and there is no river passing nearby. The groundwater is mainly fissure water stored in the rock mass, but the water is abundant. It is very likely that large water inrush would occur during the TBM construction, which would seriously threaten construction safety. Actually, water inrush occurred many times in the SS tunnel during the 19 months from August 2017 to February 2019, which seriously affected the tunneling efficiency of the TBM, as shown in Figure 3. Accordingly, there is a need to assess the risk level of water inrush in this section of the SS tunnel and provide risk warning for the subsequent tunneling of the TBM.

Risk Assessment of Water Inrush
The occurrence of tunnel water inrush not only endangers engineering construction in arid areas, but also has an important irreversible impact on the ecological environment at the top of the tunnel. Consequently, analyzing the risk of tunnel water inrush has important engineering practical value in arid areas.
A new risk analysis method, based on the normal cloud theory, is proposed to analyze the risk of water inrush, as presented in Figure 4. First, a complete evaluation index system and evaluation criteria of water inrush risk assessment are determined. Second, the numerical characteristics of the cloud are calculated by the normal cloud model theory. Third, the weight of each evaluation index is calculated by a comprehensive weight calculating algorithm. Finally, the comprehensive certainty degree is determined by the comprehensive weight and certainty degree belonging to the corresponding level. Then, the risk level of water inrush can be determined.

Risk Assessment of Water Inrush
The occurrence of tunnel water inrush not only endangers engineering construction in arid areas, but also has an important irreversible impact on the ecological environment at the top of the tunnel. Consequently, analyzing the risk of tunnel water inrush has important engineering practical value in arid areas.
A new risk analysis method, based on the normal cloud theory, is proposed to analyze the risk of water inrush, as presented in Figure 4. First, a complete evaluation index system and evaluation criteria of water inrush risk assessment are determined. Second, the numerical characteristics of the cloud are calculated by the normal cloud model theory. Third, the weight of each evaluation index is calculated by a comprehensive weight calculating algorithm. Finally, the comprehensive certainty degree is determined by the comprehensive weight and certainty degree belonging to the corresponding level. Then, the risk level of water inrush can be determined.

Evaluation Index System
In actual engineering, there are many factors that cause water inrush, and the reasons are very complicated. According to previous research results [15][16][17], these factors can be categorized into hydrological factors (Q 1 ), geological factors (Q 2 ), construction factors (Q 3 ), and dynamic monitoring (Q 4 ). In this paper, ten factors were selected to establish the assessment index system for water inrush, including landform and physiognomy (C 1 ), unfavorable geological conditions (C 2 ), strata inclination (C 3 ), strength of the surrounding rock (C 4 ), groundwater level (C 5 ), water supply (C 6 ), excavation disturbance (C 7 ), supporting measures (C 8 ), monitoring measurement (C 9 ), and geological prediction Water 2020, 12, 644 5 of 20 (C 10 ). C 1~C4 are geological factors, C 5 and C 6 are hydrological factors, C 7 and C 8 are construction factors, and C 9 and C 10 correspond to dynamic monitoring. In the established evaluation system, C 3 , C 4 , and C 5 are quantitative indices, and the index values were obtained by the measured data of the project. The other indices are qualitative indices, and the index values are determined by the expert scoring method. Combined with water inrush classification of the tunnel, the ten risk assessment indices were divided into four levels: low risk (I), medium risk (II), high risk (III), and higher risk (IV), as shown in Table 1. Table 1. Assessment indices and classification criteria for water inrush.

Assessment Indices I (Low Risk) II (Medium Risk) III (High Risk) IV (Higher Risk)
Landform and physiognomy (C 1 ) C 5 refers to the height difference between the tunnel floor and the groundwater level, and C 3 and C 5 are measured in meters and degrees, respectively. BQ refers to the basic quality of surrounding rock.

Landform and Physiognomy (C 1 )
Landform and physiognomy reflect the fluctuation characteristics of the terrain above the tunnel. Under different geographic and geomorphic conditions, the possibility of tunnel water inrush is different. If the terrain above the tunnel fluctuates greatly and there is a large area of depression, the possibility of tunnel water inrush is large. Conversely, if the terrain above the tunnel fluctuates little and the landform is relatively flat, the possibility of tunnel water inrush is small. In the present study, the proportion of negative landform area above the tunnel was quantified [18], and it was divided into four levels: Without negative landform (0~25), small-scale negative landform (25~50), medium-scale negative landform (50~75), and large-scale negative landform (75~100), as shown in Table 1.

Unfavorable Geological Conditions (C 2 )
Faults and large fissures not only bring difficulties to construction, but also are important unfavorable geological conditions that cause water inrush. Because of the existence of faults and fissures, the spatial distribution balance of water originally stored in the strata is destroyed in the process of tunnel construction, which easily leads to groundwater inrush in the tunnel along these paths. According to the degree of fault and fissure development in the strata [19], the unfavorable geological conditions was divided into four levels in this paper: non-developed strata (0~25), weakly developed strata (25~50), medially developed strata (50~75), and strongly developed strata (75~100), as shown in Table 1.

Strata Inclination (C 3 )
Strata development is complex and changeable in a tunnel, and the trend of strata has an effect on the flow of groundwater. The permeability of groundwater is anisotropic in rock strata, which is smaller in the vertical direction than in the horizontal direction. In this study, the angle (0~90 • ) between the trend of rock strata and the axis of the tunnel, α, was selected as the strata inclination, and the strata inclination was divided into four grades: 0~10 • , 10~35 • , 35~75 • , and 75~90 • [15], as shown in Table 1. In the process of tunnel construction, the strength of the surrounding rock is an important reference index for the risk of water inrush and construction safety. If the integrity of the surrounding rock is poor, the surrounding rock is prone to collapse and result in deformation in tunnel construction, and the risk of water inrush is high. According to the size of the surrounding rock classification standard (BQ), the strength of surrounding rock was divided into four levels in this paper: (BQ) > 450, 350 < (BQ) ≤ 450, 250 < (BQ) ≤ 350 and (BQ) ≤ 250 [20], as shown in Table 1.

Groundwater Level (C 5 )
For underground engineering excavation, the distribution characteristics of groundwater are not only closely related to water inrush, but are also one of the important driving factors. Meanwhile, when water inrush occurs, groundwater not only erodes the surrounding rock, but is also the carrier of material movement [21,22]. When the groundwater level is high, the water pressure is higher, and the risk of water inrush is greater. In this paper, the height difference h between the tunnel floor and the groundwater level are defined as the effective groundwater level, and h is classified into four grades: 0 < h ≤ 10 m, 10 < h ≤ 30 m, 30 < h ≤ 60 m and h > 60 m [23].

Water Supply (C 6 )
Generally, rivers, surface precipitation, lakes, and groundwater in rock mass are the main supply of the water sources of water inrush, and adequate water supply may lead to large scale and long duration of water inrush. Therefore, when analyzing the cause of water inrush, we should not only analyze the flow characteristics of groundwater in the rock mass, but also consider the potential threat of water supply to tunnel water inrush. If the water supply capacity is strong, the possibility of tunnel water inrush is high [16,24]. Based on the supply capacity of water resources, it was divided into four levels in this paper: weaker supply capacity (0~25), weak supply capacity (25~50), strong supply capacity (50~75), and stronger supply capacity (75~100), as shown in Table 1.

Excavation Disturbance (C 7 )
Construction factors have great uncertainty in the process of tunnel construction. The construction parameters, technical level of builders and engineering management may affect the construction quality of the tunnel. If the construction scheme is not reasonable, the excavation disturbance to the tunnel will be larger, and the possibility of water inrush will be greater. Therefore, it is necessary to formulate a reasonable excavation scheme to minimize the disturbance to surrounding rock. According to the disturbance degree of the surrounding rock during tunnel excavation, it was divided into four grades in this paper: smaller disturbance (0~25), small disturbance (25~50), large disturbance (50~75), and larger disturbance (75~100), as shown in Table 1.
The original balance of the strata around the tunnel is destroyed during tunnel construction, causing the deformation or collapse of the surrounding rock. The timeliness and rigidity of the supporting are important indicators to measure whether the supporting measures are reasonable. Reasonable supporting measures are critical to ensuring the safety of tunnel construction, in addition, it can prevent the weathering of the surrounding rock and reduce the impact of groundwater on the tunnel construction [25]. Based on the rationality of the supporting measure, it was divided into four grades: unreasonable (0~60), basically reasonable (60~70), reasonable (70~85), and more reasonable (85~100), as shown in Table 1. The hydrogeological conditions and the deformation of the surrounding rock dynamically change in time and space in the process of tunnel construction, and monitoring measurement is an important means to capture these dynamic changes. According to continuous field observation and measurement, we can obtain a lot of feedback information, such as seepage and displacement of the surrounding rock, which is of great significance for us to analyze the potential risk of water inrush. According to the rationality of monitoring measurement scheme, it was divided into four grades in this paper: unreasonable (0~60), basically reasonable (60~70), reasonable (70~85), and more reasonable (85~100), as shown in Table 1.

Geological Prediction (C 10 )
Advanced geological forecast is one of the important technical methods to detect the risk of water inrush when the tunnel crosses the water-rich section [18]. A reasonable geological prediction method can not only effectively detect the water storage situation of the rock mass in front of the tunnel face, but can also provide an early risk warning for the builders. This gives technicians enough time to formulate a detailed pretreatment scheme to minimize the probability of water inrush, and then reduce the adverse impact of water inrush on tunnel construction. According to the rationality of the geological prediction method, it was divided into four grades in this paper: unreasonable (0~60), basically reasonable (60~70), reasonable (70~85), and more reasonable (85~100), as shown in Table 1.

The Normal Cloud Model
The cloud model proposed by Li is a mathematical theory that is used to analyze the uncertainty problem, which is based on the probability statistics method and the fuzzy set theory [26,27]. It can achieve a good conversion between a qualitative concept and its quantitative instantiations. So far, it has been widely applied in the stability classification of surrounding rocks, evaluation of water-use efficiency, prevention and control of disaster, and other areas [28][29][30].

Cloud and Cloud Droplets
Let U be a quantitative universe of discourse, which can be described by precise numerical values, and A is a qualitative concept of U. Assume that there is a quantitative number x ∈ U, and x is a random realization of concept A. For arbitrary element x, there is a random variable with stabilization tendency, called the certainty degree µ(x) [31,32]. The parameter µ(x) can be estimated as [33,34]: where the distribution of x on U is called a cloud, and each (x, µ(x)) is a cloud drop. According to the definition of cloud, cloud is composed of numerous cloud drops, and it can transform the qualitative concept into a quantitative value. Each cloud drop is a quantitative representation of a qualitative concept in the U space. The certainty degree of a cloud drop reflects the fuzziness and randomness of a concept.

Numerical Characteristics of Cloud
Theoretically, the distribution of x on U can be described by the cloud numerical characteristics (E x , E n , H e ). The expected value (E x ) represents the central value of the concept in the U space. Entropy (E n ) is the uncertainty measurement of the qualitative concept, which reflects the value range of the distribution of cloud drops. The dispersion degree of cloud drops can be indicated by Hyper entropy (H e ), which can reflect the uncertainty degree of E n and control the thickness of cloud drops [35,36]. For instance, Figure 5 shows a normal cloud diagram, in which the cloud numerical characteristics are E x = 0, E n = 2, and H e = 0.2. The three cloud numerical characteristics are calculated by formulas (2)-(4): where and are the maximum and minimum boundary values of the corresponding grade standard, respectively. Based on the practical situation, can be set as an appropriate constant (k < 0.5) to ensure that the cloud thickness is appropriate. In this paper, was 0.01 [37]. If a variable has only an upper or lower limit, such as (−∞, ] or [ , ∞), the boundary parameters can be determined by the maximum and minimum values of the variable.

Cloud Generator
Cloud generator is a tool to implement the conversion between the qualitative index and the quantitative index. The cloud generator is divided into forward cloud generator (FCG) and backward cloud generator (BCG), as shown in Figure 6. The FCG maps from qualitative to quantitative, it can generate a large number of cloud drops to form a cloud diagram by inputting the cloud numerical characteristics ( , , ). In contrast, the BCG can convert quantitative values into the qualitative a concept [38,39]. In this study, the forward cloud generator was used to analyze the risk of water inrush. If a variable ∈ is a random realization of the concept and obeys the Gaussian The three cloud numerical characteristics are calculated by Formulas (2)- (4): where C max and C min are the maximum and minimum boundary values of the corresponding grade standard, respectively. Based on the practical situation, k can be set as an appropriate constant (k < 0.5) to ensure that the cloud thickness is appropriate. In this paper, k was 0.01 [37]. If a variable has only an upper or lower limit, such as (−∞, c min ] or [c max , +∞) , the boundary parameters can be determined by the maximum and minimum values of the variable.

Cloud Generator
Cloud generator is a tool to implement the conversion between the qualitative index and the quantitative index. The cloud generator is divided into forward cloud generator (FCG) and backward cloud generator (BCG), as shown in Figure 6. The FCG maps from qualitative to quantitative, it can generate a large number of cloud drops to form a cloud diagram by inputting the cloud numerical characteristics (E x , E n , H e ). In contrast, the BCG can convert quantitative values into the qualitative a concept [38,39]. The three cloud numerical characteristics are calculated by formulas (2)-(4): where and are the maximum and minimum boundary values of the corresponding grade standard, respectively. Based on the practical situation, can be set as an appropriate constant (k < 0.5) to ensure that the cloud thickness is appropriate. In this paper, was 0.01 [37]. If a variable has only an upper or lower limit, such as (−∞, ] or [ , ∞), the boundary parameters can be determined by the maximum and minimum values of the variable.

Cloud Generator
Cloud generator is a tool to implement the conversion between the qualitative index and the quantitative index. The cloud generator is divided into forward cloud generator (FCG) and backward cloud generator (BCG), as shown in Figure 6. The FCG maps from qualitative to quantitative, it can generate a large number of cloud drops to form a cloud diagram by inputting the cloud numerical characteristics ( , , ). In contrast, the BCG can convert quantitative values into the qualitative a concept [38,39]. In this study, the forward cloud generator was used to analyze the risk of water inrush. If a variable ∈ is a random realization of the concept and obeys the Gaussian In this study, the forward cloud generator was used to analyze the risk of water inrush. If a variable x ∈ A is a random realization of the concept A and obeys the Gaussian distribution x ∼ N(E x ,É n 2 ), and E n obeys the Gaussian distributionÉ n ∼ N(E n , H 2 e ). The certainty degree of x to A, µ(x), satisfies [27]: whereÉ n is a normally distributed random number, the distribution of x in U is called a normal cloud. The specific algorithm of the forward normal cloud model is as follows: 1.
Calculate expectation E x , entropy E n , and hyper entropy H e ; 2.
Generate a normally distributed random numberÉ ni with expectation E ni and variance H ei : Generate a normally distributed random number x i with expectation E xi and varianceÉ ni : Generate a cloud drop (x i , µ(x i )), and repeat steps 1 to 5 until N cloud drops are generated.

Subjective Weight Calculation
The analytic hierarchy process (AHP), proposed by Saaty, is an effective method to analyze complicated decisions [40]. It can analyze the underlying logical relationship of the study objects in quantitative terms. Moreover, AHP has been widely applied in many fields. In this study, the specific algorithm of AHP is shown below: Firstly, according to the established evaluation index system, the relative importance of any two indices at each level is analyzed by the 1~9 scale method (see Table 2) [41]. Then, the judgment matrix M can be constructed by the Formula (6).
Secondly, the weight vector and the maximum eigenvalues of M can be calculated as: where w i is the weight of the ith assessment index, M is the judgment matrix, n is the dimension of judgment matrix M, and λ max is the maximum eigenvalue of M. Finally, to ensure the accuracy of the results, the consistency of judgment matrix M should be verified by Formula (9). CR = (λ max −n)/(n−1) RI where RI is the random consistency index value. The values of RI are shown in Table 3.  If CR is less than 0.1, it indicates that the judgment matrix satisfies the consistency requirement, and therefore the normalized eigenvector can be used as the weight vector.

Objective Weight Calculation
Water inrush is a dynamic process affected by many factors, and the process is highly uncertain. In the information theory, entropy is an index to measure the disorder and degree of confusion of a system. It can not only objectively reflect the amount of effective information provided by the data, but can also measure the uncertainty of the system. The smaller the entropy value of an evaluation index, the greater the amount of effective information provided by the index, and the greater the weight of the index. This method can avoid human interference in the calculation process as much as possible so that the weight assignment is more in line with engineering practice [42,43].
Assuming that there are m cases of tunnel water inrush and n evaluation indices, the original data matrix X = (x ij ) m×n can be constructed as: In order to eliminate the influence of dimension in the different types of data, the original data should be standardized. The larger an index value, the lower the risk, and the standardization is calculated by the following formula: The smaller an index value, the lower the risk, and the standardization can be calculated as: With Formulas (11) and (12), the normalized matrix B = (b ij ) m×n can be obtained: According to the definition of entropy, the entropy value of the jth evaluation index is calculated by Equations (14) and (15): The weight of the jth evaluation index can be determined by the following formula:

Combination Weight Calculation
The calculation of the combination weight is used to integrate the above two weights so that the weight reflects both the subjective experience and the objective facts. The calculation method is as follows: where p is the combination weight of the evaluation indices, and β is the weight coefficient.

Calculation of the Synthetic Certainty Degree
According to the forward normal cloud generator, the normal cloud model for the risk assessment of water inrush can be generated, and the certainty degree of x belonging to the corresponding cloud can also be calculated. Then, combined with the combination weights calculated by Formula (17), the synthetic certainty degree, U, can be calculated by Formula (18).
where p i is the combination weight of the ith evaluation index, and µ(x) is the certainty degree of the ith assessment index.
Then, the weight vector and the maximum eigenvalues of M can be calculated by Formulas (7) and (8). The AHP-based weight of each assessment index is presented in Table 5, where λ max = 10.8765. In addition, RI = 1.49 and CR = 0.0654 < 0.1. For the evaluation indices of C 1~C3 and C 5~C7 , the smaller the value, the lower the risk of water inrush. For the evaluation indices of C 4 and C 8~C10 , the larger the value, the lower the risk of water inrush. Based on the evaluation index values of water inrush in Table 4, the normalized matrix B can be obtained by Formulas (11) and (12).
With Formulas (14)-(16), the entropy-based weight of each assessment index can be calculated, as presented in Table 5.

Combination Weight Calculation
According to the calculation results of the AHP-based weight and the entropy-based weight, the combination weight of each assessment index can be calculated by Formula (17). In this study, β was taken as 0.6. The calculation results of the index weight are presented in Table 5.

Cloud Model of Water Inrush
According to the assessment index values and classification criteria of water inrush, the cloud numerical characteristics (E x , E n , H e ) of different risk levels can be calculated by Formulas (2)-(4), as shown in Table 6. Then, the certainty degree of each assessment index belonging to different risk levels can be calculated by the forward cloud generator, as shown in Table 7. Combined with the combination weight of each evaluation index, the synthetic certainty degree can be calculated by Formula (18), as shown in Table 8. The risk level of water inrush was determined by the maximum synthetic certainty degree, and the distributions of the certainty degrees of each evaluation index in the four risk levels of water inrush are illustrated in Figure 7. In Figure 7, the X-axis is the value of each assessment index, and the Y-axis is the corresponding value of the certainty degree.

Risk Level Assessment
In order to assess the reliability of the cloud model, the ideal point and gray relation projection methods were compared in this paper [45][46][47]. The evaluation results calculated by different methods are shown in Table 8.
Based on the comparative analysis of the evaluation results in Table 8, the results calculated by  the three methods are exactly same in the sections of SD52+160-SD52+135, SD51+980-SD51+917,  SD51+280-SD51+264, SD51+264-SD51+212, SD51+212-SD51+170 and SD50+660-SD50+617. The final  levels of SD52+160-SD52+135 and SD51+280-SD51+264 are IV, the risk levels of SD51+264-SD51+212  and SD51+212-SD51+170 are III, and the final levels of SD51+980-SD51+917 and SD50+660-SD50+617  are II. Moreover, the synthetic certainty degree of the SD52+135-SD52+060 section belonging to the different risk levels is: U(IV) = 0.1373, U(III) = 0.1180, U(II) = 0.0889, and U(I) = 0.0245. Hence, the probability of water inrush in this section is level IV (very high), which is the same as the evaluation result calculated by the gray relation projection method. The evaluation result calculated by the ideal point method is level III, which is lower than the above two methods, and the final level of SD52+135-SD52+060 is IV. The maximum synthetic certainty degree of SD52+060-SD51+980 is: U(II) = 0.1880. This shows that the result calculated by the normal cloud model is of level II, which is the same as the result of the ideal point method and is lower than the result of the gray relation projection method. Therefore, the final level of SD52+060-SD51+980 is II. The synthetic certainty degree of SD51+917-SD51+280 and SD51+170-SD50+660 belonging to the different risk levels is U(I) > U(II) > U(III) > U(IV). Compared with the results calculated by the ideal point and gray relation projection methods, the result calculated by the normal cloud model is lower, and the final levels of SD51+917-SD51+280 and SD51+170-SD50+660 are I.
In summary, the 70% (7/10) final levels are the same as the ideal point and gray relation projection method. It can be seen that the final results calculated by the three methods have high consistency.

Excavation Verification
According to the statistical analysis of the location and the volume of water inrush in the SS tunnel, within the study area from SD52+160 to SD50+617, the measured maximum water inrush volumes of the ten samples are presented as Figure 8. In view of the advanced geological prediction result of the complex frequency conductivity method, the risk of water inrush is very high for SD52+160, as seen in Figure 9. Obvious water inrush appears on the tunnel face when the TBM reaches SD52+160, as seen in Figure 10a, and the measured maximum water inrush volume is 758 m 3 /h. Moreover, large-scale water inrush also appears on the right wall of SD52+090 and SD51+265, as seen in Figure 10b,c, and the measured maximum water inrush volume is 587 m 3 /h and 950 m 3 /h, respectively. Several outflow points of water inrush occur at the TBM shield of SD50+660, as seen in Figure 10d, and the measured maximum water inrush volume is 220 m 3 /h. Therefore, the evaluation results calculated by the proposed cloud model method In view of the advanced geological prediction result of the complex frequency conductivity method, the risk of water inrush is very high for SD52+160, as seen in Figure 9. Obvious water inrush appears on the tunnel face when the TBM reaches SD52+160, as seen in Figure 10a, and the measured maximum water inrush volume is 758 m 3 /h. Moreover, large-scale water inrush also appears on the right wall of SD52+090 and SD51+265, as seen in Figure 10b,c, and the measured maximum water inrush volume is 587 m 3 /h and 950 m 3 /h, respectively. Several outflow points of water inrush occur at the TBM shield of SD50+660, as seen in Figure 10d, and the measured maximum water inrush volume is 220 m 3 /h. Therefore, the evaluation results calculated by the proposed cloud model method show high consistency with the actual excavation situations, which demonstrates that this risk assessment method is of high accuracy and reliability for practical engineering. In view of the advanced geological prediction result of the complex frequency conductivity method, the risk of water inrush is very high for SD52+160, as seen in Figure 9. Obvious water inrush appears on the tunnel face when the TBM reaches SD52+160, as seen in Figure 10a, and the measured maximum water inrush volume is 758 m 3 /h. Moreover, large-scale water inrush also appears on the right wall of SD52+090 and SD51+265, as seen in Figure 10b,c, and the measured maximum water inrush volume is 587 m 3 /h and 950 m 3 /h, respectively. Several outflow points of water inrush occur at the TBM shield of SD50+660, as seen in Figure 10d, and the measured maximum water inrush volume is 220 m 3 /h. Therefore, the evaluation results calculated by the proposed cloud model method show high consistency with the actual excavation situations, which demonstrates that this risk assessment method is of high accuracy and reliability for practical engineering.

Conclusions
Considering geological, hydrological, and construction factors, as well as dynamic monitoring, a new multi-index evaluation method was proposed in this paper to analyze the risk of tunnel water inrush based on the normal cloud model. The weight coefficients of the assessment indicators were calculated by the objective and subjective weights, and the certainty degree of each assessment indicator belonging to the corresponding cloud was obtained by the normal cloud generator. The risk level of water inrush was determined by the comprehensive determination degree.
The proposed evaluation method was applied to assess the risk of water inrush in the SS tunnel

Conclusions
Considering geological, hydrological, and construction factors, as well as dynamic monitoring, a new multi-index evaluation method was proposed in this paper to analyze the risk of tunnel water inrush based on the normal cloud model. The weight coefficients of the assessment indicators were calculated by the objective and subjective weights, and the certainty degree of each assessment indicator belonging to the corresponding cloud was obtained by the normal cloud generator. The risk level of water inrush was determined by the comprehensive determination degree.
The proposed evaluation method was applied to assess the risk of water inrush in the SS tunnel constructed by a TBM in the arid area of Northwest China. The evaluation results are not only consistent with the results of the ideal point and gray relation projection methods, but also agree well with the actual excavation situations. This demonstrates that this new risk evaluation method of water inrush has high accuracy and reliability and can provide a reference for related projects.
Since water inrush is characterized as complex, fuzzy, and random, the application of cloud theory to the risk prediction of water inrush is only a preliminary attempt, and there are still some issues to be further studied, such as the flow mechanism and distribution characteristics of groundwater in the fractures of the surrounding rock. In addition, it is worth exploring the introduction of the multi-dimensional cloud model into the risk analysis of water inrush.