Next Article in Journal
Effects of Urbanization on Watershed Evapotranspiration and Its Components in Southern China
Previous Article in Journal
Application of Artificial Neural Networks in Forecasting a Standardized Precipitation Evapotranspiration Index for the Upper Blue Nile Basin
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
College of Hydraulic and Civil Engineering, Xinjiang Agricultural University, Urumqi 830052, China
2
Sinohydro Engineering Bureau 4 Co., Ltd., Xining 810007, China
3
Xinjiang Design Institute of Water Conservancy and Hydroelectric Power, Urumqi 830001, China
*
Author to whom correspondence should be addressed.
Water 2020, 12(3), 644; https://doi.org/10.3390/w12030644
Submission received: 29 January 2020 / Revised: 20 February 2020 / Accepted: 24 February 2020 / Published: 27 February 2020
(This article belongs to the Section Hydraulics and Hydrodynamics)

Abstract

:
Tunnel water inrush is complex, fuzzy, and random, and it is affected 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 final 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.

1. 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.

2. Study Area

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.

3. 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.

3.1. 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 (Q1), geological factors (Q2), construction factors (Q3), and dynamic monitoring (Q4). In this paper, ten factors were selected to establish the assessment index system for water inrush, including landform and physiognomy (C1), unfavorable geological conditions (C2), strata inclination (C3), strength of the surrounding rock (C4), groundwater level (C5), water supply (C6), excavation disturbance (C7), supporting measures (C8), monitoring measurement (C9), and geological prediction (C10). C1~C4 are geological factors, C5 and C6 are hydrological factors, C7 and C8 are construction factors, and C9 and C10 correspond to dynamic monitoring. In the established evaluation system, C3, C4, and C5 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.

3.1.1. Landform and Physiognomy (C1)

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.

3.1.2. Unfavorable Geological Conditions (C2)

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.

3.1.3. Strata Inclination (C3)

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.

3.1.4. Strength of Surrounding Rock (C4)

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.

3.1.5. Groundwater Level (C5)

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].

3.1.6. Water Supply (C6)

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.

3.1.7. Excavation Disturbance (C7)

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.

3.1.8. Supporting Measure (C8)

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.

3.1.9. Monitoring Measurement (C9)

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.

3.1.10. Geological Prediction (C10)

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.

3.2. 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].

3.2.1. 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]:
μ : U [ 0 , 1 ] x U μ ( x )
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.

3.2.2. 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):
E x = C m a x + C m i n 2
E n = C m a x C m i n 6
H e = k
where C m a x   and   C m i n 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 m i n ] or [ c m a x , + ) , the boundary parameters can be determined by the maximum and minimum values of the variable.

3.2.3. 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].
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 , E n ´ 2 ) , and E n ´ obeys the Gaussian distribution E n ´ ~ N ( E n , H e 2 ) . The certainty degree of x to A , μ ( x ) , satisfies [27]:
μ ( x ) = e ( x E x ) 2 2 E n ´ 2
where E 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:
  • Calculate expectation E x , entropy E n , and hyper entropy H e ;
  • Generate a normally distributed random number E n i ´ with expectation E n i and variance H e i : E n i ´ ~ N ( E n i , H e i 2 ) ;
  • Generate a normally distributed random number x i with expectation E x i and variance E n i ´ : x i ~ N ( E x i , E n i ´ 2 ) ;
  • Calculate the certainty degree of x i , μ ( x i ) = e ( x i E x i ) 2 2 E n i ´ 2 ;
  • Generate a cloud drop ( x i ,   μ ( x i ) ), and repeat steps 1 to 5 until N cloud drops are generated.

3.3. Evaluation Index Weight Calculation

3.3.1. 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).
M = [ M 1 M 1 M 1 M 2 M 1 M n M 2 M 1 M 2 M 2 M 2 M n M n M 1 M n M 2 M n M n ] = [ M 11 M 12 M 1 n M 21 M 22 M 2 n M n 1 M n 2 M n n ]
Secondly, the weight vector and the maximum eigenvalues of M can be calculated as:
w i = j = 1 n M i j n / i = 1 n j = 1 n M i j n ( i = 1 , 2 , n )
λ m a x = i = 1 n ( M w ) i n w i
where w i is the weight of the i th assessment index, M is the judgment matrix, n is the dimension of judgment matrix M , and λ m a x 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).
C R = ( λ m a x n ) / ( n 1 ) R I
where R I is the random consistency index value. The values of R I are shown in Table 3.
If C R 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.

3.3.2. 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 i j ) m × n can be constructed as:
X = [ x 11 x 12 x 1 n x 21 x 22 x 2 n x m 1 x m 2 x m n ]
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:
b i j = x i j m i n { x i j } m a x { x i j } m i n { x i j }
The smaller an index value, the lower the risk, and the standardization can be calculated as:
b i j = m a x { x i j } x i j m a x { x i j } m i n { x i j }
With Formulas (11) and (12), the normalized matrix B = ( b i j ) m × n can be obtained:
B = [ b 11 b 12 b 1 n b 21 b 22 b 2 n b m 1 b m 2 b m n ]
According to the definition of entropy, the entropy value of the j th evaluation index is calculated by Equations (14) and (15):
r i j = ( 1 + b i j ) / i = 1 m ( 1 + b i j )
E j = ( l n   m ) 1 i = 1 m r i j   l n   r i j
The weight of the j th evaluation index can be determined by the following formula:
w j = 1 E j j = 1 n ( 1 E j )

3.3.3. 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:
p = β w i + ( 1 β ) w j
where p is the combination weight of the evaluation indices, and β is the weight coefficient. β can be determined by expert scoring. If the decisions prefer expert experience, β [ 0.5 ,   1 ] . If the decisions prefer measured data, β [ 0 ,   0.5 ] [25,44]. For example, if experts give [ β 1 , β 2 ,   β 3 ] = [ 0.6 ,   0.7 ,   0.6 ] , then the weight coefficient is β = β 1 + β 2 + β 3 3 = 0.63 .

3.4. 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).
U = i n μ ( x ) p i
where p i is the combination weight of the i th evaluation index, and μ ( x ) is the certainty degree of the i th assessment index.

4. Results and Discussion

In this study, the SS tunnel from SD52+160 to SD50+617 was selected as the research area. According to the actual geological survey and the hydrological conditions of the SS tunnel, the research area was divided into ten sections. Moreover, based on the actual hydrogeological investigation, construction, management, advance geological forecasting, and the expert scoring method, the risk factor parameter values of water inrush in the ten samples could be obtained, as shown in Table 4.

4.1. Weight Calculation

4.1.1. Subjective Weight Calculation Based on AHP

Based on the established assessment index system of water inrush in this study, the judgment matrix M can be obtained by AHP.
M = [ 1 1 / 4 5 4 1 / 3 1 / 3 5 7 6 1 4 1 6 5 2 3 6 8 7 4 1 / 5 1 / 6 1 1 / 3 1 / 5 1 / 5 2 2 2 1 / 6 1 / 4 1 / 5 3 1 1 / 5 1 / 4 3 4 5 1 / 4 3 1 / 2 5 5 1 2 6 7 7 3 3 1 / 3 5 4 1 / 2 1 5 6 6 3 1 / 5 1 / 6 1 / 2 1 / 3 1 / 6 1 / 5 1 2 3 1 / 5 1 / 7 1 / 8 1 / 2 1 / 4 1 / 7 1 / 6 1 / 2 1 2 1 / 6 1 / 6 1 / 7 1 / 2 1 / 5 1 / 7 1 / 6 1 / 3 1 / 2 1 1 / 5 1 1 / 4 6 4 1 / 3 1 / 3 5 6 5 1 ]
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 λ m a x = 10.8765. In addition, R I = 1.49 and C R = 0.0654 < 0.1.

4.1.2. Objective Weight Calculation Based on Entropy

For the evaluation indices of C1~C3 and C5~C7, the smaller the value, the lower the risk of water inrush. For the evaluation indices of C4 and C8~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).
B = [ 0.0000 0.0000 0.0000 0.2857 0.8889 0.0000 0.2000 0.0000 0.0000 0.2000 0.1429 0.0714 0.3333 0.2857 1.0000 0.5000 0.3000 0.1429 0.2500 0.2000 0.7143 0.5714 0.5000 0.8571 0.2778 1.0000 0.8000 0.8571 0.7500 0.9000 1.0000 0.7143 0.2500 0.3571 0.4444 1.0000 0.7000 0.7143 0.7500 0.8000 1.0000 0.9286 0.3333 0.7143 0.2778 1.0000 1.0000 1.0000 1.0000 1.0000 0.8571 0.0000 0.5000 0.1429 0.1111 1.0000 0.0000 0.2857 0.5000 0.0000 0.7143 0.5000 0.9167 0.0000 0.3889 1.0000 0.4000 0.4286 0.5000 0.5000 1.0000 0.3571 1.0000 0.4286 0.6667 1.0000 0.3000 0.5714 0.5000 0.6000 1.0000 1.0000 0.0833 0.7143 0.1667 1.0000 0.9000 1.0000 0.7500 1.0000 1.0000 0.7143 0.3333 1.0000 0.0000 1.0000 0.8000 0.7143 0.7500 0.8000 ]
With Formulas (14)–(16), the entropy-based weight of each assessment index can be calculated, as presented in Table 5.

4.1.3. 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.

4.2. 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.

4.3. 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.

5. 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 m3/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 m3/h and 950 m3/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 m3/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.

6. 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.

Author Contributions

X.W. conceived and wrote the paper; K.S. and Q.S. supervised the findings of this work; and H.D. and M.C. analyzed the data. All authors discussed the results and contributed to the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the Scientific Research and Innovation Project of Xinjiang Agricultural University (XJAUGRI2019002) and the Science and Technology Support Plan Project of Xinjiang (201233132).

Acknowledgments

The authors appreciate the aid provided by the reviewers and editors to improve the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liu, G.; Ma, F.; Liu, G.; Zhao, H.; Guo, J.; Cao, J. Application of Multivariate Statistical Analysis to Identify Water Sources in A Coastal Gold Mine, Shandong, China. Sustainability 2019, 11, 3345. [Google Scholar] [CrossRef] [Green Version]
  2. Xue, Y.; Liu, Y.; Dang, F.N.; Liu, J. Assessment of the Nonlinear Flow Characteristic of Water Inrush Based on the Brinkman and Forchheimer Seepage Model. Water 2019, 11, 855. [Google Scholar] [CrossRef] [Green Version]
  3. Wu, H.N.; Shen, S.L.; Liao, S.M.; Yin, Z.Y. Longitudinal structural modelling of shield tunnels considering shearing dislocation between segmental rings. Tunn. Undergr. Sp. Tech. 2015, 50, 317–323. [Google Scholar] [CrossRef]
  4. Yang, W.M.; Fang, Z.D.; Yang, X.; Shi, S.S. Experimental Study of Influence of Karst Aquifer on the Law of Water Inrush in Tunnels. Water 2018, 10, 1211. [Google Scholar] [CrossRef] [Green Version]
  5. Shen, S.L.; Wu, H.N.; Cui, Y.J.; Yin, Z.Y. Long-term settlement behavior of metro tunnels in the soft deposits of Shanghai. Tunn. Undergr. Sp. Tech. 2014, 40, 309–323. [Google Scholar] [CrossRef]
  6. Zhang, Q.S.; Li, S.C.; Ge, Y.H.; Xu, Z.H.; Liu, R.T. Study on Risk Evaluation Method of Water Inrush and Integrated Geological Prediction Technology in High-Risk Karst Tunnel. In Proceedings of the Geohunan International Conference 2011, Hunan, China, 9–11 June 2011. [Google Scholar]
  7. Lyu, H.M.; Shen, J.S.; Arulrajah, A. Assessment of Geohazards and Preventative Countermeasures Using AHP Incorporated with GIS in Lanzhou, China. Sustainability 2018, 10, 304. [Google Scholar] [CrossRef] [Green Version]
  8. Zhang, J.; Yang, T. Study of a roof water inrush prediction model in shallow seam mining based on an analytic hierarchy process using a grey relational analysis method. Arab. J. Geosci. 2018, 11, 153. [Google Scholar] [CrossRef]
  9. Chen, Y.F.; Hong, J.M.; Zheng, H.K.; Li, Y. Evaluation of Groundwater Leakage into a Drainage Tunnel in Jinping-I Arch Dam Foundation in Southwestern China: A Case Study. Rock Mech. Rock Eng. 2016, 49, 961–979. [Google Scholar] [CrossRef]
  10. Hao, Y.Q.; Rong, X.L.; Ma, L.J.; Fan, P.X.; Lu, H. Uncertainty Analysis on Risk Assessment of Water Inrush in Karst Tunnels. Math Probl. Eng. 2016, 2, 1–11. [Google Scholar] [CrossRef] [Green Version]
  11. Lyu, H.M.; Sun, W.J.; Shen, S.L. Flood risk assessment in metro systems of megacities using a GIS-based modeling approach. Sci. Total Environ. 2018, 626, 1012–1025. [Google Scholar] [CrossRef]
  12. Li, X.P.; Li, Y.N. Research on risk assessment system for water inrush in the karst tunnel construction based on GIS: Case study on the diversion tunnel groups of the Jinping II hydropower station. Tunn. Undergr. Sp. Tech. 2014, 40, 182–191. [Google Scholar] [CrossRef]
  13. Li, S.C.; Zhou, Z.Q.; Li, L.P.; Xu, Z.H.; Zhang, Q.Q.; Shi, S.S. Risk assessment of water inrush in karst tunnels based on attribute synthetic evaluation system. Tunn. Undergr. Sp. Tech. 2013, 38, 50–58. [Google Scholar] [CrossRef]
  14. Chu, H.D.; Xu, G.L. Risk assessment of water inrush in karst tunnels based on two-class fuzzy comprehensive evaluation method. Arab. J. Geosci. 2017, 10, 179. [Google Scholar] [CrossRef]
  15. Yang, X.L.; Zhang, S. Risk assessment model of tunnel water inrush based on improved attribute mathematical theory. J. Cent. South Univ. 2018, 25, 379–391. [Google Scholar] [CrossRef]
  16. Li, T.Z.; Yang, X.L. Risk assessment model for water and mud inrush in deep and long tunnels based on normal grey cloud clustering method. KSCE J. Civ. Eng. 2018, 22, 1991–2001. [Google Scholar] [CrossRef]
  17. Li, L.P.; Lei, T.; Li, S.C.; Zhang, Q.Q.; Xu, Z.H. Risk assessment of water inrush in karst tunnels and software development. Arab. J. Geosci. 2015, 8, 1843–1854. [Google Scholar] [CrossRef]
  18. Li, S.C.; Wu, J. A multi-factor comprehensive risk assessment method of karst tunnels and its engineering application. Bull. Eng. Geol. Environ. 2019, 78, 1761–1776. [Google Scholar] [CrossRef]
  19. Li, S.C.; Wu, J.; Xu, Z.H. Unascertained Measure Model of Water and Mud Inrush Risk Evaluation in Karst Tunnels and Its Engineering Application. KSCE J. Civ. Eng. 2017, 4, 1170–1182. [Google Scholar] [CrossRef]
  20. Wang, J.; Li, S.C.; Li, L.P.; Lin, P.; Xu, Z.H. Attribute recognition model for risk assessment of water inrush. Bull. Eng. Geol. Environ. 2019, 78, 1057–1071. [Google Scholar] [CrossRef]
  21. Moon, J.; Fernandez, G. Effect of excavation-induced groundwater level drawdown on tunnel inflow in a jointed rock mass. Eng. Geol. 2010, 110, 33–42. [Google Scholar] [CrossRef]
  22. Shin, J.H. Analytical and combined numerical methods evaluating pore water pressure on tunnels. Geotechnique 2010, 60, 141–145. [Google Scholar] [CrossRef]
  23. Wang, Y.C.; Jing, H.W.; Yu, L.Y.; Su, H.J. Set pair analysis for risk assessment of water inrush in karst tunnels. Bull. Eng. Geol. Environ. 2017, 76, 1199–1207. [Google Scholar] [CrossRef]
  24. Wang, X.T.; Li, S.C.; Xu, Z.H.; Lin, P. Analysis of Factors Influencing Floor Water Inrush in Coal Mines: A Nonlinear Fuzzy Interval Assessment Method. Mine Water Environ. 2019, 38, 81–92. [Google Scholar] [CrossRef]
  25. Xu, Z.H.; Li, S.C.; Li, L.P. Risk assessment of water or mud inrush of karst tunnels based on analytic hierarchy process. Rock Soil Mech. 2011, 32, 1757–1766. [Google Scholar]
  26. Li, D.Y.; Meng, H.J.; Shi, X.M. Membership clouds and membership cloud generators. J. Comput. Res. Dev. 1995, 6, 15–20. [Google Scholar]
  27. Li, D.Y.; Liu, C.G.; Gan, W.Y. A new cognitive model: Cloud model. Int. J. Intell. Syst. 2009, 24, 357–375. [Google Scholar] [CrossRef]
  28. Wang, J.C.; Guo, J. Research on Rock Mass Quality Classification Based on An Improved Rough Set-Cloud Model. IEEE Access 2019, 7, 123710–123724. [Google Scholar] [CrossRef]
  29. Guan, X.J.; Qin, H.D.; Meng, Y.; Wu, Z.N. Comprehensive evaluation of water-use efficiency in China’s Huai river basin using a cloud-compound fuzzy matter element-entropy combined model. J. Earth Syst. Sci. 2019, 128, 179. [Google Scholar] [CrossRef] [Green Version]
  30. Yang, S.M.; Han, X.Q.; Cao, B.; Li, B.; Yan, F. Cloud-Model-Based Method for Risk Assessment of Mountain Torrent Disasters. Water 2018, 10, 830. [Google Scholar] [CrossRef] [Green Version]
  31. Wang, G.Y.; Xu, C.L.; Li, D.Y. Generic normal cloud model. Inf. Sci. 2014, 280, 1–15. [Google Scholar] [CrossRef]
  32. Liu, Z.B.; Shao, J.F.; Xu, W.Y.; Xu, F. Comprehensive Stability Evaluation of Rock Slope Using the Cloud Model-Based Approach. Rock Mech. Rock Eng. 2014, 47, 2239–2252. [Google Scholar] [CrossRef]
  33. Wang, X.T.; Li, S.C.; Ma, X.Y.; Xue, Y.G. Risk Assessment of Rockfall Hazards in a Tunnel Portal Section Based on Normal Cloud Model. Pol. J. Environ. Stud. 2017, 26, 2295–2306. [Google Scholar] [CrossRef]
  34. Li, X.C.; Zhong, D.H.; Ren, B.Y. Study on evaluation of rock mass groutability of dam foundation based on Fuzzy RES-cloud model. J. Hydraul. Eng. 2017, 48, 1311–1323. [Google Scholar] [CrossRef]
  35. Deng, W.H.; Wang, G.Y.; Zhang, X.R. A novel hybrid water quality time series prediction method based on cloud model and fuzzy forecasting. Chemom. Intell. Lab. Syst. 2015, 149, 39–49. [Google Scholar] [CrossRef]
  36. Yang, X.J.; Yan, L.L.; Peng, H.; Gao, X.D. Encoding words into Cloud models from interval-valued data via fuzzy statistics and membership function fitting. Knowl. Based Syst. 2014, 55, 114–124. [Google Scholar] [CrossRef]
  37. Wang, Y.C.; Yin, X.; Jing, H.W.; Liu, R.C.; Su, H.J. A novel cloud model for risk analysis of water inrush in karst tunnels. Environ. Earth Sci. 2016, 75, 1450. [Google Scholar] [CrossRef]
  38. Lin, Y.; Zhou, K.; Li, J. Application of cloud model in rock burst prediction and performance comparison with three machine learning algorithms. IEEE Access 2018, 6, 30958–30968. [Google Scholar] [CrossRef]
  39. Zhang, T.; Yan, L.S.; Yang, Y. Trust evaluation method for clustered wireless sensor networks based on cloud model. Wirel. Netw. 2018, 24, 777–797. [Google Scholar] [CrossRef]
  40. Saaty, T.L. How to Make a Decision: The Analytic Hierarchy Process. Inf. J. Appl. Anal. 1994, 24, 19–43. [Google Scholar] [CrossRef]
  41. Saaty, T.L. Decision-making with the AHP: Why is the principal eigenvector necessary. Eur. J. Oper. Res. 2003, 145, 85–91. [Google Scholar] [CrossRef]
  42. Long, Y.; Yang, Y.; Lei, X.; Tian, Y.; Li, Y. Integrated Assessment Method of Emergency Plan for Sudden Water Pollution Accidents Based on Improved TOPSIS, Shannon Entropy and a Coordinated Development Degree Model. Sustainability 2019, 11, 510. [Google Scholar] [CrossRef] [Green Version]
  43. Won, K.; Chung, E.-S.; Choi, S.-U. Parametric Assessment of Water Use Vulnerability Variations Using SWAT and Fuzzy TOPSIS Coupled with Entropy. Sustainability 2015, 7, 12052–12070. [Google Scholar] [CrossRef] [Green Version]
  44. Li, Z.; Liu, W.B.; Dong, F. Research on Capability Evaluation of Equipment Support Based on Comprehensive Weighting Method and Cloud Model. Command Control Simul. 2018, 40, 69–74. [Google Scholar] [CrossRef]
  45. Xu, C.; Liu, X.L.; Wang, E.Z.; Zheng, Y.L.; Wang, S.J. Rockburst prediction and classification based on the ideal-point method of information theory. Tunn. Undergr. Sp. Tech. 2018, 81, 382–390. [Google Scholar] [CrossRef]
  46. Tahmasebi, S.; Moradi, P.; Ghodsi, S. An ideal point based many-objective optimization for community detection of complex networks. Inf. Sci. 2019, 502, 125–145. [Google Scholar] [CrossRef]
  47. Chen, L.Y.; Deng, Y. A new failure mode and effects analysis model using Dempster-Shafer evidence theory and grey relational projection method. Eng. Appl. Artif. Intel. 2018, 76, 13–20. [Google Scholar] [CrossRef]
Figure 1. The geographical location of the study area.
Figure 1. The geographical location of the study area.
Water 12 00644 g001
Figure 2. The geological profile of the SS (Shuang-san) tunnel from SD52+500 to SD50+300. TBM, tunnel boring machine.
Figure 2. The geological profile of the SS (Shuang-san) tunnel from SD52+500 to SD50+300. TBM, tunnel boring machine.
Water 12 00644 g002
Figure 3. The measured data of the TBM monthly tunneling length and monthly pumping volume of the SS (Shuang-san) tunnel during the 19 months from August 2017 to February 2019.
Figure 3. The measured data of the TBM monthly tunneling length and monthly pumping volume of the SS (Shuang-san) tunnel during the 19 months from August 2017 to February 2019.
Water 12 00644 g003
Figure 4. The risk assessment process of tunnel water inrush based on a normal cloud model.
Figure 4. The risk assessment process of tunnel water inrush based on a normal cloud model.
Water 12 00644 g004
Figure 5. A normal cloud model with parameters E x = 0 , E n = 2 , and H e = 0.2 .
Figure 5. A normal cloud model with parameters E x = 0 , E n = 2 , and H e = 0.2 .
Water 12 00644 g005
Figure 6. Two different types of cloud generators. (a) forward cloud generator; (b) backward cloud generator.
Figure 6. Two different types of cloud generators. (a) forward cloud generator; (b) backward cloud generator.
Water 12 00644 g006
Figure 7. Cloud of each assessment index generated by the forward cloud generator. (a) landform and physiognomy (C1), unfavorable geological conditions (C2), water supply (C6), and excavation disturbance (C7); (b) strata inclination (C3); (c) strength of the surrounding rock (C4); (d) groundwater level (C5); (e) supporting measures (C8), monitoring measurement (C9), and geological prediction (C10).
Figure 7. Cloud of each assessment index generated by the forward cloud generator. (a) landform and physiognomy (C1), unfavorable geological conditions (C2), water supply (C6), and excavation disturbance (C7); (b) strata inclination (C3); (c) strength of the surrounding rock (C4); (d) groundwater level (C5); (e) supporting measures (C8), monitoring measurement (C9), and geological prediction (C10).
Water 12 00644 g007
Figure 8. The measured maximum water inrush volume of the ten samples in the SS (Shuang-san) tunnel.
Figure 8. The measured maximum water inrush volume of the ten samples in the SS (Shuang-san) tunnel.
Water 12 00644 g008
Figure 9. The advanced geological prediction result of the complex frequency conductivity method at SD52+160.
Figure 9. The advanced geological prediction result of the complex frequency conductivity method at SD52+160.
Water 12 00644 g009
Figure 10. Practical situation of water inrush in the SS (Shuang-san) tunnel. (a) water inrush at SD52+160; (b) water inrush at SD52+090; (c) water inrush at SD51+265; (d) water inrush at SD50+660.
Figure 10. Practical situation of water inrush in the SS (Shuang-san) tunnel. (a) water inrush at SD52+160; (b) water inrush at SD52+090; (c) water inrush at SD51+265; (d) water inrush at SD50+660.
Water 12 00644 g010
Table 1. Assessment indices and classification criteria for water inrush.
Table 1. Assessment indices and classification criteria for water inrush.
Assessment IndicesI (Low Risk)II (Medium Risk)III (High Risk)IV (Higher Risk)
Landform and physiognomy (C1)<2525~5050~7575~100
Unfavorable geological conditions (C2)<2525~5050~7575~100
Strata inclination (C3)<1010~3535~7575~90
Strength of surrounding rock (C4)(BQ) > 450350 < (BQ) ≤ 450250 < (BQ) ≤ 350(BQ) ≤ 250
Groundwater level (C5)<1010~3030~60>60
Water supply (C6)<2525~5050~7575~100
Excavation disturbance (C7)<2525~5050~7575~100
Supporting measures (C8)85~10070~8560~70<60
Monitoring measurement (C9)85~10070~8560~70<60
Geological prediction (C10)85~10070~8560~70<60
C5 refers to the height difference between the tunnel floor and the groundwater level, and C3 and C5 are measured in meters and degrees, respectively. BQ refers to the basic quality of surrounding rock.
Table 2. The 1~9 scale method of analytic hierarchy process (AHP).
Table 2. The 1~9 scale method of analytic hierarchy process (AHP).
Value   of   M i j Two-by-Two Comparison
1 M i and M j are equally important
3 M i is slightly more important than M j
5 M i is more important than M j
7 M i is substantially more important than M j
9 M i is absolutely more important than M j
2, 4, 6, 8the middle value of two adjacent judgements
ExplanationWhen M j and M i are compared, the value is
the reciprocal of M i and M j scalar.
Table 3. The random consistency index values.
Table 3. The random consistency index values.
n12345678910
R I 000.580.901.121.241.321.411.451.49
Table 4. The risk factor parameter values of water inrush in the SS (Shuang-san) tunnel.
Table 4. The risk factor parameter values of water inrush in the SS (Shuang-san) tunnel.
SampleC1C2C3C4C5C6C7C8C9C10
SD52+160–SD52+1354585663201043070556550
SD52+135–SD52+0604080623201022065607050
SD52+060–SD51+9802045604001151040858085
SD51+980–SD51+9171035633301121045808080
SD51+917–SD51+2801020623801151030908590
SD51+280–SD51+2641585603001181080657540
SD51+264–SD51+2122050552801131060707565
SD51+212–SD51+1701060543401081065757570
SD51+170–SD50+6601015653801171035908090
SD50+660–SD50+6171035624201201040808080
Table 5. Weights of the assessment indices of water inrush.
Table 5. Weights of the assessment indices of water inrush.
IndexC1C2C3C4C5
AHP-based Weight0.10840.26720.03200.05410.2043
Entropy-based Weight0.10290.12870.10050.09850.1074
Combination Weight0.10620.21180.05940.07190.1655
IndexC6C7C8C9C10
AHP-based Weight0.15910.02900.02100.01800.1067
Entropy-based Weight0.07710.10100.10370.07200.1082
Combination Weight0.12630.05780.05410.03960.1073
Table 6. Cloud numerical characteristics ( E x , E n , H e ) of each assessment index.
Table 6. Cloud numerical characteristics ( E x , E n , H e ) of each assessment index.
IndexIIIIIIIV
C1(12.5, 4.1667, 0.01)(37.5, 4.1667, 0.01)(62.5, 4.1667, 0.01)(87.5, 4.1667, 0.01)
C2(12.5, 4.1667, 0.01)(37.5, 4.1667, 0.01)(62.5, 4.1667, 0.01)(87.5, 4.1667, 0.01)
C3(5, 1.6667, 0.01)(22.5, 4.1667, 0.01)(55, 6.6667, 0.01)(82.5, 2.5, 0.01)
C4(550, 33.3333, 0.01)(400, 16.6667, 0.01)(300, 16.6667, 0.01)(125, 41.6667, 0.01)
C5(5, 1.6667, 0.01)(20, 3.3333, 0.01)(45, 5, 0.01)(90, 10, 0.01)
C6(12.5, 4.1667, 0.01)(37.5, 4.1667, 0.01)(62.5, 4.1667, 0.01)(87.5, 4.1667, 0.01)
C7(12.5, 4.1667, 0.01)(37.5, 4.1667, 0.01)(62.5, 4.1667, 0.01)(87.5, 4.1667, 0.01)
C8(92.5, 2.5, 0.01)(77.5, 2.5, 0.01)(65, 1.6667, 0.01)(30, 10, 0.01)
C9(92.5, 2.5, 0.01)(77.5, 2.5, 0.01)(65, 1.6667, 0.01)(30, 10, 0.01)
C10(92.5, 2.5, 0.01)(77.5, 2.5, 0.01)(65, 1.6667, 0.01)(30, 10, 0.01)
Table 7. Certainty degrees of assessment indices C1~C10.
Table 7. Certainty degrees of assessment indices C1~C10.
SampleC1C2C3C4C5
SD52+160–SD52+1350.19870.83500.25790.48640.3761
SD52+135–SD52+0600.83660.19590.57760.48640.4871
SD52+060–SD51+9800.19440.19870.75571.00000.0441
SD51+980–SD51+9170.83590.83550.48750.19850.0888
SD51+917–SD51+2800.83590.19440.57760.48690.0441
SD51+280–SD51+2640.83550.83500.75571.00000.0199
SD51+264–SD51+2120.19440.01121.00000.48730.0707
SD51+212–SD51+1700.83590.83520.98880.05620.1985
SD51+170–SD50+6600.83590.83550.32640.48690.0259
SD50+660–SD50+6170.83590.83550.57760.48600.0110
SampleC6C7C8C9C10
SD52+160–SD52+1350.20220.19750.04421.00000.1355
SD52+135–SD52+0600.19440.83520.01130.01110.1355
SD52+060–SD51+9800.83590.83550.01110.60360.0111
SD51+980–SD51+9170.83590.19840.60360.60360.6036
SD51+917–SD51+2800.83590.19660.60380.01110.6038
SD51+280–SD51+2640.83590.19791.00000.60470.6064
SD51+264–SD51+2120.83590.83520.01110.60471.0000
SD51+212–SD51+1700.83590.83570.60470.60470.0118
SD51+170–SD50+6600.83590.83610.60380.60360.6038
SD50+660–SD50+6170.83590.83550.60360.60360.6036
Table 8. Calculation results of the different evaluation methods applied to the SS (Shuang-san) tunnel.
Table 8. Calculation results of the different evaluation methods applied to the SS (Shuang-san) tunnel.
SampleSynthetic Certainty DegreeThe Cloud Model MethodIdeal Point MethodGray Relation Projection Method
U(I)U(II)U(III)U(IV)
SD52+160–SD52+13500.04660.10130.2561IVIVIV
SD52+135–SD52+0600.02450.08890.11800.1373IVIIIIV
SD52+060–SD51+9800.12600.18800.04490.0073IIIIIII
SD51+980–SD51+9170.19430.30980.04320.0147IIIIII
SD51+917–SD51+2800.33340.04640.03430.0073IIIII
SD51+280–SD51+2640.19430.02400.17090.2567IVIVIV
SD51+264–SD51+2120.12620.02400.25300.0117IIIIIIIII
SD51+212–SD51+1700.19430.05670.28930.0329IIIIIIIII
SD51+170–SD50+6600.46880.10720.01940.0043IIIII
SD50+660–SD50+6170.19430.38160.03430.0018IIIIII

Share and Cite

MDPI and ACS Style

Wang, X.; Shi, K.; Shi, Q.; Dong, H.; Chen, M. 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. Water 2020, 12, 644. https://doi.org/10.3390/w12030644

AMA Style

Wang X, Shi K, Shi Q, Dong H, Chen M. 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. Water. 2020; 12(3):644. https://doi.org/10.3390/w12030644

Chicago/Turabian Style

Wang, Xin, Kebin Shi, Quan Shi, Hanwei Dong, and Ming Chen. 2020. "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" Water 12, no. 3: 644. https://doi.org/10.3390/w12030644

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