Next Article in Journal
Compressed Stabilized Earthen Blocks and Their Use in Low-Cost Social Housing
Previous Article in Journal
Common Space and Behavior at the Border between Slum and Metropolitan Area: The Case of “Catambor” and “Alvalade”
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Study on Seismic Damage Risk Assessment of Mountain Tunnel Based on the Extension Theory

1
Key Laboratory of Transportation Tunnel Engineering, Ministry of Education, Southwest Jiaotong University, Chengdu 610031, China
2
Shudao Investment Group Co., Ltd., Chengdu 610041, China
3
Sichaun Lu Shi Expressway Group Co., Ltd., Chengdu 610041, China
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(6), 5294; https://doi.org/10.3390/su15065294
Submission received: 20 February 2023 / Revised: 10 March 2023 / Accepted: 13 March 2023 / Published: 16 March 2023

Abstract

:
In this paper, the extension method is applied to the assessment of seismic damage risk in mountain tunnels. Based on various considerations such as the geological conditions of the tunnel site, the state of the tunnel, the structural situation and the earthquake, the seismic damage risk assessment index for mountainous tunnels is proposed. The range of values of different assessment indexes under each grade is quantified. The seismic damage risk of the tunnel is divided into four classes, and assessed based on the extension theory. The specific assessment process is as follows: Firstly, five tunnels affected by the Wenchuan earthquake were selected to construct the matter element to be assessed, and construct the classical domain and nodal domain; then, the entropy weight method and analytical hierarchy process (AHP) were used to determine the weights of each index; finally, the correlation function values were calculated to obtain the seismic damage risk level of the tunnels. Among them, the seismic damage risk level of the Longchi tunnel and Longxi tunnel is “high risk”, and the seismic damage risk level of the Longdongzi tunnel, Youyi tunnel and Maanshi tunnel is “moderate risk”. The five tunnels are sorted according to the eigenvalues of the seismic damage risk level as follows: Longchi tunnel > Longxi tunnel > Longdongzi tunnel > Youyi tunnel > Maanshi tunnel. Combined with the seismic damage investigation of the Wenchuan earthquake, the actual seismic damage of tunnel is consistent with the assessment results, which verifies the accuracy and effectiveness of this method. The seismic damage risk assessment model for mountain tunnels proposed in this paper has certain reference value for the future seismic damage risk assessment of tunnels.

1. Introduction

The investigation of seismic damage caused by the Tangshan earthquake, Chi-Chi earthquake, Hanshin earthquake, Wenchuan earthquake and other earthquakes shows that the tunnels have suffered varying degrees of damage [1,2,3,4,5]. The most common forms of damage include tunnel deformation and cracking, lining leakage, concrete block falling, and collapse in the tunnel portal section and the cavern section. Once the tunnel collapses, the losses incurred are immeasurable. With the increasing use of underground space, the number of tunnels located in high seismic intensity areas is also increasing. The seismic intensity represents the strength of the impact of an earthquake event on buildings. The higher the seismic intensity, the more severe the earthquake damage will be when it occurs. How to quantify the seismic damage risk of tunnel is a complex work.
Tunnel safety assessment includes many aspects, such as the risk of collapse during the construction period, the risk of water and mud inrush in karst tunnels, etc. Ou et al. selected 11 risk indexes such as excavation span, grade of surrounding rock to evaluate the collapse risk of deep buried tunnel during construction, and successfully applied it to Yuxi tunnel [6]. Kim et al. collected 24 tunnel collapse cases, established a database, and quantified the influencing factors of tunnel collapse using an analytic hierarchy process (AHP) [7]; Chen et al. proposed a mountain tunnel collapse risk assessment model integrating case-based reasoning, rough set theory and um-spa theory, and evaluated the Xiucun tunnel crossing the shattered fault zone [8]. Lin et al. developed a simple and practical software package for assessing the risk of water gushing in tunnels using cloud model and variable weight theory [9]; Luo et al. established a shield tunneling risk assessment model based on the characteristics of shield tunnels in soft soil, and successfully applied it to Fuzhou Metro Line 6 [10]. Zhen and Gong et al. carried out the damage and safety assessment of the tunnel during the operation period, which provided a reference for the reinforcement and maintenance of the tunnel [11,12]. In summary, a significant number of studies have been carried out for collapse risk, water gushing risk and excavation risk in tunnels; however, almost all studies have not considered the seismic damage risk factor when assessing the safety of tunnels. Scholars have also carried out a large number of researches on the dynamic response of tunnels under earthquake action through an investigation method [13,14], theoretical analysis [15,16], numerical simulation [17,18,19] and shaking table tests [20,21], respectively. Their studies have focused on the investigation and classification of seismic damage of tunnels, mechanical behavior of tunnels, and anti-seismic measures, as they are more concerned with the tunnel itself. However, the seismic damage risk of tunnels is also related to geological conditions, tunnel status, tectonic conditions and many other factors. With the increasing number and scale of tunnels in high intensity areas, the safety of tunnels under earthquakes poses a great challenge; a reasonable seismic damage risk assessment of the tunnel is helpful to understand the safety of the tunnel, and can strengthen the tunnel with high seismic damage risk in a targeted manner so as to reduce the loss of life and property caused by the tunnel damage when the earthquake occurs.
In recent years, a large number of assessment methods have been used in engineering safety assessment, including fault tree analysis [22], risk matrix method [23,24], analytic hierarchy process [25], Bayesian network [26,27], artificial neural network method [28], and fuzzy comprehensive assessment method [29], among others; however, these methods have some shortcomings when in use. For example, the fault tree analysis tends to omit potential risk factors. The risk matrix method and analytic hierarchy process need to rank the risk factors, and the selection of weight depends on the experience of experts. The artificial neural network method requires a large amount of data as a learning sample. Bayesian models can make poor predictions due to the assumed prior model. When using the fuzzy comprehensive assessment method, a larger number of indexes in the same index level will lead to smaller weight coefficients and affect the assessment accuracy.
The shortcomings of the above-mentioned assessment methods make them somewhat limited in engineering applications. The risk of seismic damage in tunnels is a multi-objective decision problem that requires the consideration of numerous factors, and the magnitude of the impact of these factors on the tunnel is uncertain. The extension method can transform incompatible conflicts into compatible relationships, thus achieving the best global decision. It can quantitatively describe the seismic damage risk of tunnel and make up for the deficiency of single subjective or objective weighting. It can also consider the fuzziness and uncertainty of risk factors. In addition, to avoid the occurrence of invalid correlation functions, which are caused by the magnitudes of the matter-element out of the range of the nodal domain, the conventional extension method is improved. In the present work, the seismic damage risk factors of mountain tunnels are summarized and a seismic damage risk assessment index system is established based on a large amount of literature research and seismic damage research. Then, the analytic hierarchy process and the improved entropy weight method are used to determine the weight of each risk assessment index, five damaged tunnels in the Wenchuan earthquake are selected to form the matter-elements to be evaluated. The correlation function value is calculated, and the seismic damage risk level of the tunnel is finally obtained. The assessment results are compared with the actual earthquake damage to verify the correctness of the model proposed in this paper.

2. Risk Assessment Model Based on Extension Theory

The assessment model based on extension theory can transform multi-objective decision problems into single-objective decision problems, and quantitatively express the assessment results, as it has strong practicability and has been widely used in the field of engineering safety assessment [30].
The construction process of matter-element extension assessment model generally includes five steps.
  • Establishing the classical domain and nodal domain matter-element matrix of assessment indexes;
  • Constructing the correlation function;
  • Determining the index weights and calculating the comprehensive correlation degree;
  • Determining the assessment level.

2.1. Determination of Classic Domains and Nodal Domains

For a particular risk assessment index, the classical domain is a range of values describing the characteristics of the assessment object for different assessment levels, and the nodal domain is the maximum range that the risk assessment index can take. According to the assessment standard, the tunnel seismic damage risk assessment index is divided into several grades; the classical domain is the range of values of the seismic damage risk assessment index for the tunnel at the specified level, which is expressed as shown in Equation (1):
R j = N j , C i , U i j = N j c 1 U 1 j c 2 U 2 j c n U n j = N j c 1 a 1 j , b 1 j c 2 a 2 j , b 2 j c n a n j , b n j
where R j ( j = 1,2 , , m ) is the j-th matter-element; N j is the j-th assessment level of the risk index; C i ( i = 1,2 , , n ) is the i-th assessment index; U i j is N j on the range of quantitative values specified by index C i , U i j = a i j , b i j is the range taken by the assessment index C i at the j-th level, that is, the classical domain of N j .
Referring to the way of determining the classical domain, the nodal domain can be obtained as shown in Equation (2):
R p = P , C i , U i p = P c 1 U 1 p c 2 U 2 p c n U n p = P c 1 a 1 p , b 1 p c 2 a 2 p , b 2 p c n a n p , b n p
where P is all assessment levels of tunnel risk indexes; U i p = a i p , b i p , U i p is all assessment levels which are in proximity the range of values taken by the assessment index C i , that is, the nodal domain of P .
The value of the classic domain is generally determined according to the existing specifications in the industry norms, objective laws, industry expert experience judgment, data statistical analysis, etc. The classical domain cannot extend beyond the nodal domain and the value range of the nodal domain is also determined after the classical domain is determined.

2.2. Matter-Element to Be Assessed

For a tunnel to be assessed, the basic information or data analysis results of the tunnel are represented by matter-element R; then, R is called matter-element to be assessed, as shown in Equation (3).
R = P c 1 u 1 c 2 u 2 c n u n
where P is all assessment levels for an index of a tunnel to be evaluated; u i is the value of the tunnel about the assessment index c i , that is, the specific value of the tunnel to be assessed.

2.3. Correlation Functions

The correlation function expresses the degree of affiliation of each index of the object to be evaluated with respect to each assessment level, and the association degree of the index is calculated according to the definition of distance in extension theory; its calculation formula is as shown in Equation (4).
K j ( u i ) = ρ ( u i , U i j ) D ( u i , U i j , U i p )
where ρ u i , U i j represents the distance between u i and the classical domain U i j = a i j , b i j ; its calculation formula is as shown in Equation (5).
ρ u i , U i j = u i a i j + b i j 2 b i j a i j 2
where D ( u i , U i j , U i p ) represents the place value of the point with respect to the interval U i j and the interval U i p , U i j U i p , place value and distance can quantitatively express the positional relationship between points and intervals; its calculation formula is as shown in Equation (6).
D u i , U i j , U i p = U i j = b i j a i j u i U i j ρ u i , U i p ρ u i , U i j u i U i j

2.4. Determination of the Assessment Level

According to the actual situation, choose the appropriate method to determine the weight coefficient of each index. Commonly used methods include Delphi method, frequency statistical analysis method, entropy method, AHP, etc. The comprehensive correlation degree K j ( p ) under the risk level J is obtained through the correlation degree and the corresponding weight coefficient, as shown in Equation (7).
K j p = i = 1 n ω i K j ( u i )
where ω i is the weight coefficient of the i-th assessment index; K j ( u i ) is the degree of association of the i-th assessment index at level J.
The larger the K j p value of the matter-element p to be evaluated, the smaller the deviation of the index from the assessment level; if K j 0 p = max [ K j p ( j = 1,2 , , m ) ] , then the seismic damage risk of the tunnel to be evaluated is level J0.
The risk level eigenvalues j * are calculated as follows.
K j p = K j p m i n [ K j p ] m a x [ K j p ] m i n [ K j p ]
j * = j = 1 m j × K j p j = 1 m K j p
j * corresponds to the extent to which p is biased towards the adjacent level; for example, j * = 2.6 , which indicates that p belongs to the second level and is biased towards the third level.

3. Improvement of Matter-Element Extension Model for Mountain Tunnel

When the conventional matter-element extension method is used for quantitative analysis, once the value of the matter-element u i is out of the range of the nodal domain, it leads to the invalid of the correlation function to be calculated, and to inaccurate assessment results. On the other hand, in determining the weights of each index, both the subjective and objective weighting methods have certain drawbacks. The subjective weighting method relies too much on subjective judgment, and the weight calculation results will have obvious differences due to different subjective experiences of researchers. The calculation results of the objective weighting method are based on actual data and lack the correction of subjective experience. Therefore, it is necessary to integrate the advantages of the subjective and objective weighting method, combining subjective judgment based on practical problems and theoretical calculation based on objective data to determine the comprehensive weight of indexes, try to make the weight of each index match the actual importance, improve the credibility of the index weight, and improve the accuracy of the assessment results.

3.1. Standardized Processing of Classical Domains and Matter Elements

In response to the problem that the correlation function fails due to the value of the matter-element being outside the range of the nodal domain, this paper normalizes the matter-element matrix, classical domain and nodal domain matrix to be evaluated; the specific improvement method is as follows:
R = 1 W h e n   t h e   v a l u e   o f   t h e   r i s k   i n d e x   i s   l a r g e r ,   t h e   r i s k   i s   l o w e r . 0   W h e n   t h e   v a l u e   o f   t h e   r i s k   i n d e x   i s   s m a l l e r ,   t h e   r i s k   i s   l o w e r .
Then the matter-element value to be evaluated is normalized as shown in Equation (11).
u i = R u i b i p ( b i p u i ) 1 R ( u i a i p ) R b i p a i p u i ( a i p , b i p ) 1 R u i a i p
where u i is normalized matter-element value, u i is the actual value of the matter-element to be assessed; the above processing also solves the problem of non-uniform dimensions.
Normalize the classical domain and nodal domain functions in the same way. The normalized classical domain and nodal domain are as shown in Equations (12) and (13).
U i j = a i j a i p b i p a i p , b i j a i p b i p a i p
U i p = [ 0,1 ]

3.2. Comprehensive Weight of Improved Extension Method Index

The method of determining the weight of indexes can be divided into subjective weighting method and objective weighting method; subjective weighting method is mainly employs on subjective experience and judgment to assign weights to assessment indexes, such as AHP, Delphi method, etc. The objective weighting method calculates the weight value of each index based on the information reflected by the sample data, such as the entropy weighting method. When using the subjective weighting method to determine the weights, it involves more a reflection of human intention and the decision result may be subjective and arbitrary. The entropy weighting method follows only mathematical theory. Both the subjective and objective weighting methods have certain drawbacks. This paper adopts a comprehensive subjective–objective weighting method which can effectively avoid the excessive influence of subjective factors in the assignment process and fully exploit the information in the sample data.

3.2.1. Analytic hierarchy Process to Determine the Weight of Assessment Index

First, a judgment matrix Q = x i i n × n ( i = 1,2 , 3 , , n ) is constructed by inviting a number of experienced experts to score. When constructing the judgment matrix, two related factors are compared according to the scale h a b ; specific standards are as shown in Table 1.
Consistency check based on the constructed judgment matrix; check formula is as shown in Equation (14).
C R = C I / R I
where C I = ( λ m a x n ) / ( n 1 ) C I is consistency check standard, λ m a x is maximal characteristic root of the judgment matrix.
R I is an assessment of the average random consistency index, related to the order n of the judgment matrix (see Table 2 for details).
When the ratio C R = C I / R I < 0.10 of random consistency is less than 0.10, it can be considered that the single sorting results of the hierarchy have good consistency; otherwise, the values of the elements of the judgment matrix must be readjusted.
The weight vector is obtained using the geometric average approach, and the elements of each row of the judgment matrix are multiplied separately, as shown in Equation (15).
W i = t = 1 n x t n
The weights of each index W s are obtained by normalizing W i , as shown in Equation (16).
W i = W i t = 1 n W i

3.2.2. Determination of the Weight of Assessment Index using Improved Entropy Weight Method

The entropy weighting method is an objective weighting method [31,32], which is a method to determine the weight of an index from the objective data obtained. The entropy value of each index is calculated according to the degree of variation of each index; the index information entropy is inversely proportional to the degree of variation and weight. There are x assessment indexes, y tunnels to be assessed, and the normalized original assessment matrix X = r i j x × y ,   i = 1,2 , 3 , x ,   j = 1,2 , 3 , y is obtained, as shown in Equations (17) and (18). First, the entropy value of each index needs to be found.
E i = 1 l n x j = 1 y P i j × l n P i j
P i j = r i j j = 1 y r i j
where P i j is the proportion of the j-th tunnel under the i-th assessment index. The entropy weight of each assessment index is calculated according to Equation (19).
w i = 1 E i i = 1 x ( 1 E i )
When calculating using the above formula, if the entropy value E i is close to 1, the slight difference of different indexes may cause a large change in the entropy weight. Some improvements to the above equation can reduce the errors arising from similar indexes [33], as shown in Equation (20).
w i = e x p   k = 1 x E k + 1 E i e x p   E i l = 1 x [ e x p   k = 1 x E k + 1 E l e x p   E l ]

3.2.3. Determination of Comprehensive Weights

In order to make the difference degree of the index weight calculated by the AHP and the improved entropy weight method consistent with the difference of the corresponding distribution coefficient, in this paper, the distance function is used to calculate the comprehensive weight of the assessment index [34]. Specific steps are as follows:
d ( w i , W i ) = i = 1 x ( w i W i ) 2
d ( w i , W i ) = α β α + β = 1
where the comprehensive weight of the assessment index i is the linear combination of the results calculated by the above AHP and the improved entropy weight method.
W i * = α w i + β W i

4. Seismic Damage Risk Assessment Index System of Mountain Tunnel

4.1. The Assessment Index System

The seismic damage in the tunnel is the result of the combined action of many factors, insufficient surrounding rock strength or integrity caused by hazardous environments such as geological structure and unfavorable geology, which can easily destabilize the surrounding rock during earthquakes and cause disasters such as tunnel collapse. The stability of the surrounding rock is greatly affected by its own strength and integrity. Generally speaking, the lower the strength and integrity of the surrounding rock, the more likely it is to become unstable and produce a seismic hazard when an earthquake occurs [35]. Weathering will reduce the mechanical properties of the surrounding rock, and the degree of weathering has a great relationship with the stability of the surrounding rock [36]. In the case of loose stratum or large slope angle at the entrance area of the tunnel, the collapse of the slope and rockfall caused by the earthquake may lead to the cracking and damage of the tunnel portal and the blockage of the tunnel entrance [37].
The condition of tunnel lining itself, including tunnel defect, tunnel burial depth and size, etc., all determine whether the lining can resist the earthquake and maintain its integrity when an earthquake occurs. Tunnel defects, including cavities behind the lining, cracking, water seepage, etc. [38], will reduce the bearing capacity of the structure and thus reduce the seismic performance of the tunnel. It can be seen from the statistics of seismic damage of tunnels around the world that the magnitude of the seismic damage to the tunnel is related to the burial depth of the tunnel [39]. Tunnel earthquake damage survey shows that when the burial depth is greater than 100 m and the tunnel does not cross an active fault, the tunnel will suffer minor or no damage. Most of the tunnels with burial depths of less than 50 m have suffered some degree of damage. The span of the tunnel will also affect its seismic characteristics [16]. The earthquake intensity in the tunnel site area reflects the average level of the ground vibration intensity in the area when the earthquake occurs. The higher the earthquake intensity, the greater the potential for damage in the event of an earthquake, and unfavorable geology can exacerbate the destructive power of earthquakes [40].
In summary, the seismic damage risk assessment factors of mountain tunnels are divided into three categories: geological conditions of the tunnel site, the condition of the tunnel lining, geological structure and earthquakes. The engineering geological conditions include the surrounding rock classification, the weathering of surrounding rock and the slope angle of tunnel entrance. The condition of the tunnel includes the maximum burial depth of tunnel, tunnel defect, tunnel width, structural form and strength of the lining. Unfavorable geology and earthquakes include the length of shattered fault zone, the earthquake intensity of the tunnel site and the unfavorable geology, as shown in Figure 1.

4.2. The Seismic Damage Risk Assessment Level

Based on a large number of tunnel seismic damage investigations [3,4,5,39], combined with relevant specifications for tunnel seismic design and engineering practice, this paper defines the seismic damage risk assessment grades as four categories: low risk (I), moderate risk (II), high risk (III), and extremely high risk (IV), as shown in the Table 3.

4.3. Normalization of Assessment Index

The assessment indexes are normalized. In Table 3, there are many qualitative descriptions in the assessment indexes, such as c 1 , c 2 , c 5 , c 7 , c 9 , which quantify the qualitative indexes; the specific quantitative criteria are shown in Table 4. When quantifying seismic damage risk level according to surrounding rock classification, c 1 is the sum of the surrounding rock proportion multiplied by the grade.
In order to make different assessment indexes comparable with each other and facilitate scientific generalization, according to Equations (11)–(13), the indexes were normalized and dimensionless, and the processed ones are shown in Table 5.

5. Application Example of Seismic Damage Risk Assessment Model for Mountain Tunnel Based on Extension Theory

On 21 May 2008, am earthquake with a magnitude of 8.0 occurred in Yingxiu Town, Wenchuan County, Sichuan Province (31.0° N, 103.4° E), which caused severe damage to the surrounding buildings. In this section, the seismic damage risk assessment method proposed in this paper is used to assess five tunnels (Longchi tunnel, Longxi tunnel, Longdongzi tunnel, Youyi tunnel, Maanshi tunnel) affected by the Wenchuan earthquake.

5.1. Matter-Element to Be Assessed, Classical Domain, Nodal Domain

Taking the five mountain tunnels affected by the Wenchuan earthquake as the research objects, the seismic damage risk assessment of the mountain tunnels proposed in this paper is carried out, the specific indexes of the five tunnels to be assessed are shown in the Table 6 below.
According to Table 3, Table 4 and Table 5, and Equations (1)–(3), the value range corresponding to the seismic damage risk index of grades I to IV is taken as the classical domain Rj; then, the nodal domain R p and the matter-element to be assessed, Longchi Tunnel R 1 , Longxi Tunnel R 2 , Longdongzi Tunnel R 3 , Youyi Tunnel R 4 , and Maanshi Tunnel R 5 , respectively, are as follows:
R p = P c 1 0,1 c 2 0,1 c 3 0,1 c 4 0,1 c 5 0,1 c 6 0,1 c 7 0,1 c 8 0,1 c 9 0,1 R 1 = p 1 c 1 0.600 c 2 0.580 c 3 0.556 c 4 1.000 c 5 0.640 c 6 0.480 c 7 0.840 c 8 0.550 c 9 0.625 R 2 = p 2 c 1 0.440 c 2 0.350 c 3 0.500 c 4 1.000 c 5 0.000 c 6 0.416 c 7 0.620 c 8 0.650 c 9 0.625 R 3 = p 3 c 1 0.560 c 2 0.600 c 3 0.611 c 4 0.432 c 5 0.320 c 6 0.416 c 7 0.380 c 8 0.600 c 9 0.375 R 4 = p 4 c 1 0.700 c 2 0.660 c 3 0.500 c 4 0.283 c 5 0.320 c 6 0.340 c 7 0.350 c 8 0.015 c 9 0.375 R 5 = p 5 c 1 0.520 c 2 0.380 c 3 0.444 c 4 0.567 c 5 0.280 c 6 0.420 c 7 0.140 c 8 0.050 c 9 0.375

5.2. Correlation Function Value and Comprehensive Weight Calculation

According to Equations (4)–(6), the correlation function value of the seismic damage risk level of the tunnel to be assessed is obtained as shown in Table 7 below.
Using AHP and improved entropy weight method, the subjective weight and objective weight of nine assessment indexes of seismic damage risk level of mountain tunnel are obtained, where the weights of the AHP were calculated from the results of a questionnaire survey carried out by 20 tunneling and geological experts. There is a certain difference in the assessment index weights determined by the improved entropy weight method and the analytic hierarchy process. It is too unilateral to rely solely on pure mathematical theory or pure engineering experience for the seismic damage risk assessment of such complex mountain tunnels. Therefore, combining the subjective weight and the objective weight, the comprehensive weight is calculated according to Equations (20)–(22), which has positive significance for improving the rationality of the assessment results. From the comprehensive weight, it can be seen from Table 8 that the weight of tunnel defect, unfavorable geology, length of shattered fault zone and earthquake intensity is relatively high, and they are more important in assessing seismic damage risk.

5.3. Correlation Function Value and Comprehensive Weight Calculation

According to Equation (7), the correlation degree of the seismic damage risk of each tunnel i with respect to the level j is obtained:
K i j p = 0.558 0.232 0.019 0.184 0.427 0.180 0.169 0.378 0.328 0.166 0.173 0.344 0.151 0.101 0.353 0.549 0.096 0.145 0.412 0.605
The larger the K j p value of the matter-element p to be evaluated, the smaller the deviation of the index from the assessment level; if K j 0 p = max [ K j p ( j = 1,2 , , m ) ] , then the seismic damage risk of the tunnel to be evaluated is level J0. According to the calculation results in Table 9, the seismic damage risk levels of the five tunnels to be assessed can be obtained. Among them, the seismic damage risk level of Longxi tunnel and Longchi tunnel is Level III (high risk), and the seismic damage risk level of Longdongzi tunnel, Youyi tunnel and Maanshi tunnel is Level II (moderate risk), According to Equations (8) and (9), the eigenvalues of the seismic damage risk level of the tunnels can be obtained: j 1 * = 3.04 ,   j 2 * = 2.64 ,   j 3 * = 2.22 ,   j 4 * = 1.84 ,   j 5 * = 1.78 ; this indicates that the seismic damage risk level of Longchi tunnel is “high risk” biased towards “extremely high risk”, the seismic damage risk level of the Longxi tunnel is “high risk” biased towards “moderate risk”, the seismic damage risk level of the Longdongzi tunnel is “moderate risk” biased towards “high risk”, and the seismic damage risk level of the Youyi tunnel and Maanshi tunnel is “moderate risk” biased towards “low risk”. The five tunnels are sorted according to the eigenvalues of the seismic damage risk level as follows: Longchi tunnel > Longxi tunnel > Longdongzi tunnel > Youyi tunnel > Maanshi tunnel.
Longchi tunnel crosses four large-scale shattered fault zones, and large cracks appeared at the retaining wall. Longitudinal cracks, circumferential cracks and diagonal cracks appeared in the lining of the portal entrance, with localized falling blocks and exposed steel reinforcement (Figure 2a). Longitudinal cracks appeared in the vault, arch foot and side walls of the tunnel body section, and diagonal cracks appeared near the tunnel crossing the shattered fault zone (Figure 2b,c). Pavement protrusions and dislocations were observed in some sections (Figure 2d).
Longxi tunnel is located between the Yingxiu Fault (F3) and the Longxi Fault (F2), crossing the F8 Fault, as well as a number of minor secondary shattered faults. The tunnel was damaged for a length of 4500 m. The lining of the portal section was severely cracked (Figure 3a), the reinforcement was exposed (Figure 3b), and obvious dislocations up to 20 cm were produced at the construction joints (Figure 3c). Circular and longitudinal cracks appeared in the tunnel body section (Figure 3d), accompanied by water leakage (Figure 3e). The second lining and waterproof plate collapsed (Figure 3g), and the surrounding rock was destabilized in the severely damaged areas, producing an overall collapse (Figure 3h). The seismic damage in Longchi tunnel and Longxi tunnel is in accordance with the evaluation result of “high risk”.
The slope at the portal of Longdongzi tunnel collapsed; thus, the damage to the entrance of the tunnel was relatively large, and cracks and water leakage appeared in some parts of the lining (Figure 4a–c), which essentially conformed to the assessment results of “moderate risk” biased to “high risk”. The Youyi tunnel (Figure 4d–f) and Maanshi tunnel (Figure 4g–i) were slightly affected by the earthquake, which was mainly manifested in local cracks and water leakage, cracking and spalling of the secondary lining. The seismic damage is in accordance with the assessment result of “moderate risk”.
In general, the seismic damage investigation is essentially consistent with the assessment results. Therefore, the conclusions obtained using the seismic damage risk assessment method proposed in this paper are consistent with the actual situation and have great practical value, which can also be applied to the seismic damage risk assessment of similar projects.

6. Conclusions

In this paper, the extension theory is introduced to evaluate the seismic damage risk of mountain tunnels, taking into account the geological conditions of the tunnel site, the tunnels’ state, structural conditions and seismic damage risk factors. The assessment process includes the establishment of the risk index system, the calculation of the weight of the risk index, the division of the risk level, the quantification of the tunnel to be assessed, and the classification of the assessment results. The following conclusions were obtained:
  • A comprehensive approach to weighting using a combination of improved entropy weighting method and AHP method can effectively avoid the excessive influence of subjective factors in the assignment process and fully exploit the information in the sample data. Numerically speaking, tunnel defect, unfavorable geology, length of shattered fault zone and earthquake intensity have higher weights and are more important in assessing seismic damage risk.
  • Five tunnels affected by the Wenchuan earthquake were assessed for seismic damage risk. Among them, the Longchi tunnel is locally defective and large deformation of the surrounding rock may occur, leaving it at a significantly higher risk of seismic damage. The seismic damage risk level of Longchi tunnel is “high risk” biased towards “extremely high risk”. The Longxi tunnel crosses the longest shattered fault zone and has high ground stresses; thus, it has a “high risk” biased towards “moderate risk” for seismic damage. The seismic damage risk level of the Longdongzi tunnel is “moderate risk” biased towards “high risk”. The seismic damage risk level of the Youyi tunnel and Maanshi tunnel is “moderate risk” biased towards “low risk”. The five tunnels are sorted according to the eigenvalues of the seismic damage risk level as follows: Longchi tunnel > Longxi tunnel > Longdongzi tunnel > Youyi tunnel > Maanshi tunnel.
  • Combined with the seismic damage investigation of the Wenchuan earthquake, the Longchi tunnel and the Longxi tunnel suffered the most serious seismic damage during the earthquake, which conforms to the assessment result of “high risk”. Longdongzi tunnel, Youyi tunnel and Maanshi tunnel suffered less seismic damage, which conforms to the assessment result of “moderate risk”. The assessment results of this paper are essentially consistent with the actual tunnel damage. The accuracy and effectiveness of the method are verified.

Author Contributions

Conceptualization, C.C. and P.G.; methodology, C.C. and W.G.; validation, W.Z. and P.G.; investigation, P.G. and C.H.; data curation, F.Y.; writing—original draft preparation, C.C. and F.Y.; writing—review and editing, P.G. and W.G.; supervision, W.Z.; project administration, C.H.; funding acquisition, P.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China (Key Program, Grant No. 52130808); and National Key R&D Programme of China (Key Projects for International Science and Technology Innovation Cooperation between Governments, Grant No. 2022YFE0104300).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fu, X.; Ma, Y.; Cui, Z. Probabilistic Seismic Hazard Analysis for Fault Dislocation Magnitude Induced by Strong Earthquakes: A Case Study of the Sichuan-Yunnan Region. Sustainability 2021, 13, 2383. [Google Scholar] [CrossRef]
  2. Hu, X.; Zhou, Z.; Chen, H.; Ren, Y. Seismic Fragility Analysis of Tunnels with Different Buried Depths in a Soft Soil. Sustainability 2020, 12, 892. [Google Scholar] [CrossRef] [Green Version]
  3. Azadiab, M. The Seismic Behavior of Urban Tunnels in Soft Saturated Soils. Procedia Eng. 2011, 14, 3069–3075. [Google Scholar] [CrossRef] [Green Version]
  4. Tao, S.J.; Gao, B.; Wen, Y.M.; Zhou, X. Investigation and Analysis on Seismic Damages of Mountain Tunnels Subjected to Wenchuan Earthquake. Appl. Mech. Mater. 2011, 99–100, 273–281. [Google Scholar] [CrossRef]
  5. Wang, W.L.; Wang, T.T.; Su, J.J.; Lin, C.H.; Seng, C.R.; Huang, T.H. Assessment of damage in mountain tunnels due to the Taiwan Chi-Chi Earthquake. Tunn. Undergr. Space Technol. 2001, 16, 133–150. [Google Scholar] [CrossRef]
  6. Ou, G.; Jiao, Y.; Zhang, G.; Zou, J.; Tan, F.; Zhang, W. Collapse risk assessment of deep-buried tunnel during construction and its application. Tunn. Undergr. Space Technol. 2021, 115, 104019. [Google Scholar] [CrossRef]
  7. Kim, J.; Kim, C.; Kim, G.; Kim, I.; Abbas, Q.; Lee, J. Probabilistic tunnel collapse risk evaluation model using analytical hierarchy process (AHP) and Delphi survey technique. Tunn. Undergr. Space Technol. 2022, 120, 104262. [Google Scholar] [CrossRef]
  8. Chen, W.; Zhang, G.; Jiao, Y.; Wang, H. Unascertained Measure-Set Pair Analysis Model of Collapse Risk Evaluation in Mountain Tunnels and Its Engineering Application. KSCE J. Civ. Eng. 2021, 25, 451–467. [Google Scholar] [CrossRef]
  9. Lin, C.; Zhang, M.; Zhou, Z.; Li, L.; Shi, S.; Chen, Y.; Dai, W. A new quantitative method for risk assessment of water inrush in karst tunnels based on variable weight function and improved cloud model. Tunn. Undergr. Space Technol. 2020, 95, 103136. [Google Scholar] [CrossRef]
  10. Luo, Q.; Li, W.; Su, H.; Chen, X. Evaluating Construction Risks of Modified Shield Machine Applicable to Soft Soils Based on Fuzzy Comprehensive Evaluation Method. Math. Probl. Eng. 2020, 2020, 8861801. [Google Scholar] [CrossRef]
  11. Zhen, H.; Fu, H.L.; Zhang, J.; Wei, C.; Yue, S. Structural Damage Evaluation Method for Metro Shield Tunnel. J. Perform. Constr. Fac. 2018, 33, 04018097. [Google Scholar] [CrossRef]
  12. Gong, Z.X.; Yang, J.G. Fuzzy Matter-Element Evaluation Method for Reliability Analysis of an Existing Highway Tunnel. Adv. Mater. Res. 2010, 163–167, 3110–3113. [Google Scholar] [CrossRef]
  13. Lai, J.; He, S.; Qiu, J.; Chen, J.; Wang, L.; Wang, K.; Wang, J. Characteristics of seismic disasters and aseismic measures of tunnels in Wenchuan earthquake. Environ. Earth Sci. 2017, 76, 94. [Google Scholar] [CrossRef]
  14. Li, T. Damage to mountain tunnels related to the Wenchuan earthquake and some suggestions for aseismic tunnel construction. Bull. Eng. Geol. Environ. 2012, 71, 297–308. [Google Scholar] [CrossRef]
  15. Wang, Z.Z.; Jiang, Y.J.; Zhu, C.A. Seismic energy response and damage evolution of tunnel lining structures. Eur. J. Environ. Civ. Eng. 2019, 23, 758–770. [Google Scholar] [CrossRef]
  16. Xie, H.T. Numerical Analysis on Influence of Cross Section Shape on Earthquake Resistant Capability of Shallow-Buried Tunnel. Appl. Mech. Mater. 2013, 405–408, 1292–1296. [Google Scholar] [CrossRef]
  17. Li, P.; Song, E. Three-dimensional numerical analysis for the longitudinal seismic response of tunnels under an asynchronous wave input. Comput. Geotech. 2015, 63, 229–243. [Google Scholar] [CrossRef]
  18. Sazid, M.; Ahmed, H.A. Stability Analysis of Shallow Depth Tunnel in Weak Rock Mass: 3D Numerical Modeling Approach. J. City Dev. 2019, 1, 18–22. [Google Scholar]
  19. Sazid, M.; Singh, T.N.; Saharan, M. Risk Analysis of Mine Dump Slope Stability—A Case Study. Min. Eng. J. 2012, 12, 11–15. [Google Scholar]
  20. Zheng, H.B.; Li, P.F.; Ma, G.W.; Zhang, Q.B. Experimental investigation of mechanical characteristics for linings of twins tunnels with asymmetric cross-section. Tunn. Undergr. Space Technol. 2022, 119, 104209. [Google Scholar] [CrossRef]
  21. Cheng, X.; Zhou, X.; Liu, H.; Zhou, Y.; Shi, W. Numerical Analysis and Shaking Table Test of Seismic Response of Tunnel in a Loess Soil Considering Rainfall and Traffic Load. Rock Mech. Rock Eng. 2021, 54, 1005–1025. [Google Scholar] [CrossRef]
  22. MacLeod, J.; Tan, S.; Moinuddin, K. Reliability of fire (point) detection system in office buildings in Australia—A fault tree analysis. Fire Saf. J. 2020, 115, 103150. [Google Scholar] [CrossRef]
  23. Moradi, M.R.; Farsangi, M. Application of the Risk Matrix Method for Geotechnical Risk Analysis and Prediction of the Advance Rate in Rock TBM Tunneling. Rock Mech. Rock Eng. 2014, 47, 1951–1960. [Google Scholar] [CrossRef]
  24. Xie, G.Y.; Yue, Z. Seismic Safety Evaluation for Mountainous Highway Bridge Based on Risk Matrix. Adv. Mater. Res. 2011, 255–260, 4212–4216. [Google Scholar] [CrossRef]
  25. Darban, S.; Ghasemzadeh Tehrani, H.; Karballaeezadeh, N.; Mosavi, A. Application of Analytical Hierarchy Process for Structural Health Monitoring and Prioritizing Concrete Bridges in Iran. Appl. Sci. 2021, 11, 8060. [Google Scholar] [CrossRef]
  26. Li, T.; Zhou, Y.; Zhao, Y.; Zhang, C.; Zhang, X. A hierarchical object oriented Bayesian network-based fault diagnosis method for building energy systems. Appl. Energy 2022, 306, 118088. [Google Scholar] [CrossRef]
  27. Choi, H. A Causal Network-Based Risk Matrix Model Applicable to Shield TBM Tunneling Projects. Sustainability 2021, 13, 4846. [Google Scholar]
  28. Ma, L.Y.; Gu, B.X.; Liu, G.L. A Study on High-Rise Building Structure Selections Using Artificial Intelligence Methods. Appl. Mech. Mater. 2013, 351–352, 1198–1201. [Google Scholar] [CrossRef]
  29. Zhang, G.; Wang, C.; Jiao, Y.; Wang, H.; Qin, W.; Chen, W.; Zhong, G. Collapse Risk Analysis of Deep Foundation Pits in Metro Stations Using a Fuzzy Bayesian Network and a Fuzzy AHP. Math. Probl. Eng. 2020, 2020, 4214379. [Google Scholar] [CrossRef]
  30. Ju, Y.; Yu, Y.; Ju, G.; Cai, W. Extension set and restricting qualifications of matter-elements’ extension. In Proceedings of the Third International Conference on Information Technology and Applications (ICITA’05), Sydney, NSW, Australia, 4–7 July 2005. [Google Scholar]
  31. Gao, C.; Li, S.; Wang, J.; Li, L.; Lin, P. The Risk Assessment of Tunnels Based on Grey Correlation and Entropy Weight Method. Geotech. Geol. Eng. 2018, 36, 1621–1631. [Google Scholar] [CrossRef]
  32. Li, Y.; Zhao, L.; Suo, J. Comprehensive Assessment on Sustainable Development of Highway Transportation Capacity Based on Entropy Weight and TOPSIS. Sustainability 2014, 6, 4685–4693. [Google Scholar] [CrossRef] [Green Version]
  33. Zhang, C.; Qing, W.; Ping, C.J.; Guang, G.F.; Wen, Z. Evaluation of debris flow risk in Jinsha River based on combined weight process. Rock Soil Mech. 2011, 3, 831–836. [Google Scholar]
  34. Lang, F.Y.; Li, X.G. Multi-Sensors Information Fusion Based on Momentis Method and Euclid Distance. Adv. Mater. Res. 2011, 383–390, 5447–5452. [Google Scholar] [CrossRef]
  35. Liu, K.; Liu, B.; Fang, Y. An intelligent model based on statistical learning theory for engineering rock mass classification. Bull. Eng. Geol. Environ. 2019, 78, 4533–4548. [Google Scholar] [CrossRef]
  36. Gurocak, Z.; Yalcin, E. Excavatability and the effect of weathering degree on the excavatability of rock masses: An example from Eastern Turkey. J. Afr. Earth Sci. 2016, 118, 1–11. [Google Scholar] [CrossRef]
  37. Niu, J.; Jiang, X.; Wang, F. Stability Analysis of Rock Slope with Small Spacing Tunnel under Earthquakes and Influence of Ground Motion Parameters. Geotech. Geol. Eng. 2018, 36, 2437–2453. [Google Scholar] [CrossRef]
  38. Xu, S.; Ma, E.; Lai, J.; Yang, Y.; Liu, H.; Yang, C.; Hu, Q. Diseases failures characteristics and countermeasures of expressway tunnel of water-rich strata: A case study. Eng. Fail. Anal. 2022, 134, 106056. [Google Scholar] [CrossRef]
  39. Sharma, S.; Judd, W.R. Underground opening damage from earthquakes. Eng. Geol. 1991, 30, 263–276. [Google Scholar] [CrossRef]
  40. Niu, F.; Cai, Y.; Liao, H.; Li, J.; Tang, K.; Wang, Q.; Wang, Z.; Liu, D.; Liu, T.; Liu, C.; et al. Unfavorable Geology and Mitigation Measures for Water Inrush Hazard during Subsea Tunnel Construction: A Global Review. Water 2022, 14, 1592. [Google Scholar] [CrossRef]
  41. Specifications for Design of Highway Tunnels Section 1 Civil Engineering (JTG 3370.1-2018); Ministry of Transport of the People’s Republic of China: Beijing, China, 2018.
Figure 1. Classification of seismic damage risk assessment factors for mountain tunnels.
Figure 1. Classification of seismic damage risk assessment factors for mountain tunnels.
Sustainability 15 05294 g001
Figure 2. Seismic damage map of Longchi tunnel ((a). Broken side wall at portal (b). Longitudinal cracking and water seepage in the vault (c). Falling concrete blocks (d). Pavement protrusions).
Figure 2. Seismic damage map of Longchi tunnel ((a). Broken side wall at portal (b). Longitudinal cracking and water seepage in the vault (c). Falling concrete blocks (d). Pavement protrusions).
Sustainability 15 05294 g002
Figure 3. Seismic damage map of Longxi Tunnel ((a). Crack in the tunnel portal (b). Exposed reinforcement (c). Dislocation (d). Cracks appeared in the tunnel body section (e). Water leakage (f). Pavement protrusions (g). Second lining and waterproof plate collapsed (h). Overall collapse of tunnel).
Figure 3. Seismic damage map of Longxi Tunnel ((a). Crack in the tunnel portal (b). Exposed reinforcement (c). Dislocation (d). Cracks appeared in the tunnel body section (e). Water leakage (f). Pavement protrusions (g). Second lining and waterproof plate collapsed (h). Overall collapse of tunnel).
Sustainability 15 05294 g003
Figure 4. Seismic damage map of Longdongzi Tunnel, Youyi Tunnel and Maanshi Tunnel ((a). Slope collapse (b). Portal door cracking (c). Rockfall at the entrance of the tunnel (d). Falling concrete blocks Cracks in different directions (e,f). Water leakage (g). Longitudinal cracks (h). Cracks in side walls (i). Cracks in vault.).
Figure 4. Seismic damage map of Longdongzi Tunnel, Youyi Tunnel and Maanshi Tunnel ((a). Slope collapse (b). Portal door cracking (c). Rockfall at the entrance of the tunnel (d). Falling concrete blocks Cracks in different directions (e,f). Water leakage (g). Longitudinal cracks (h). Cracks in side walls (i). Cracks in vault.).
Sustainability 15 05294 g004
Table 1. Meaning of scale h a b .
Table 1. Meaning of scale h a b .
Scale habMeaning
1Indexes a and b are equally important
3Index a is slightly more important than b
5Index a is obviously more important than b
7Index a is strongly more important than b
9Index a is extremely more important than b
2, 4, 6, 8The degree of importance takes the middle value of two adjacent judgments
1, 1/2, …, 1/9The ratio of the effects of a to b is opposite to the above
Table 2. The value of the random consistency index RI.
Table 2. The value of the random consistency index RI.
n123456789
RI000.580.901.121.241.321.411.45
Table 3. Classification of the seismic damage risk level.
Table 3. Classification of the seismic damage risk level.
NO.Assessment IndexSeismic Damage Risk Level
Grade I
(Low Risk)
Grade II
(Moderate Risk)
Grade III
(High Risk)
Grade IV
(Extremely High Risk)
c 1 Surrounding rock classification (JTG 3370.1-2018) [41]I, IIIIIIVV
c 2 Degree of weatheringNo weathering, no impact on structureWeak weathering, little impact on structureStrong weathering, great impact on the structureFull weathering, extreme impact on the structure
c 3 Slope angle of tunnel entrance/°0~2020~4040~6060~90
c 4 Maximum burial depth of tunnel/m300~200200~100100~5050~0
c 5 Tunnel defectStructure in good conditionStructure is damaged in some places, and the degree of damage is smallStructure is damaged in some places, and the degree of damage is moderateStructure is damaged in many places, the degree of damage is serious
c 6 Tunnel width/m0~88~1212~1818~25
c 7 Unfavorable geological conditionNo unfavorable geological conditionThere is unfavorable geological condition, but the scale is small, and the impact on the safety of the tunnel is smallThe scale of unfavorable geology is large, which has a great impact on the safety of the tunnelThe scale of unfavorable geology is large, which has an extreme impact on the safety of the tunnel
c 8 Length of shattered fault zone/m0~55~2020~5050~100
c 9 Earthquake intensity≤VIVIIVIII≥IX
Table 4. Quantitative criteria for qualitative indexes.
Table 4. Quantitative criteria for qualitative indexes.
Number of Risk Assessment IndexRisk Level
Grade IGrade IIGrade IIIGrade IV
c 1 1~2.52.5~3.53.5~4.54.5~6
c 2 , c 5 , c 7 , c 9 0~2.52.5~5.05.0~7.57.5~10
Table 5. Normalization of seismic damage risk assessment index parameters.
Table 5. Normalization of seismic damage risk assessment index parameters.
NO.Assessment IndexSeismic Damage Risk Level
Grade I
(Low Risk)
Grade II
(Moderate Risk)
Grade III
(High Risk)
Grade IV
(Extremely High Risk)
c 1 Surrounding rock classification (GB/T 50218-2014)0~0.300.30~0.500.50~0.700.70~1.00
c 2 Degree of weathering0~0.250.25~0.500.50~0.750.75~1.00
c 3 Slope angle of tunnel entrance/°0~0.330.33~0.670.67~0.830.83~1.00
c 4 Maximum burial depth of tunnel/m0~0.250.25~0.500.50~0.750.75~1.00
c 5 Tunnel defect0~0.250.25~0.500.50~0.750.75~1.00
c 6 Tunnel width/m0~0.320.32~0.480.48~0.720.72~1.00
c 7 Unfavorable geological condition0~0.250.25~0.500.50~0.750.75~1.00
c 8 Length of shattered fault zone/m0~0.050.05~0.200.20~0.500.50~1.00
c 9 Earthquake intensity0~0.250.25~0.500.50~0.750.75~1.00
Table 6. Indexes of the tunnel to be assessed.
Table 6. Indexes of the tunnel to be assessed.
Assessment IndexLongchi TunnelLongxi TunnelLongdongzi TunnelYouyi TunnelMaanshi Tunnel
Surrounding rock classificationIII, IV, VIII, IV, VIII, IV, VIV, VIII, IV, V
Degree of weatheringStrong weathering, great impact on the structureWeak weathering, little impact on the structureStrong weathering, great impact on the structureStrong weathering, great impact on the structureWeak weathering, little impact on the structure
Slope angle of tunnel entrance40~6035~5540~7030~6025~55
Maximum burial depth of tunnel470830173215130
Tunnel defectStructural is damaged in some places, and the degree of damage is moderate, local leakage water in tunnelThe tunnel is still in the construction stage. No damageStructural is damaged in some places, and the degree of damage is smallStructural equipment is damaged in some places, and the degree of damage is smallStructural equipment is damaged in some places, and the degree of damage is small
Tunnel width1210.410.48.510.5
Unfavorable geological conditionThe tunnel traverses a 254 m long coal seam, where moderate-to-large deformation of the surrounding rock may occurHigh ground stress, and moderate rockburst may occurDangerous rock collapse at the entrance of the tunnelLarge deformation of surrounding rock may occur locally in the shattered fault zoneSmall collapse area at the entrance of the tunnel
Length of shattered fault zone5565601.58
Earthquake intensity88777
Table 7. The value of the correlation function.
Table 7. The value of the correlation function.
Indicators c 1 c 2 c 3 c 4 c 5 c 6 c 7 c 8 c 9
Tunnel
LongchiKI−0.429−0.440−0.337−1.000−0.520−0.250−0.787−0.526−0.500
KII−0.200−0.1600.3350.240−0.2800.000−0.680−0.438−0.250
KIII0.5000.320−0.204−1.0000.4400.000−0.360−0.1000.500
KIV−0.200−0.288−0.382−1.000−0.234−0.3330.3600.100−0.250
LongxiKI−0.241−0.222−0.254−1.0000.000−0.188−0.493−0.632−0.500
KII0.3000.4000.5000.240−1.0000.400−0.240−0.563−0.250
KIII−0.120−0.300−0.254−1.000−1.000−0.1330.480−0.3000.500
KIV−0.371−0.533−0.398−1.000−1.000−0.422−0.2550.300−0.250
LongdongziKI−0.371−0.467−0.419−0.290−0.179−0.188−0.255−0.579−0.250
KII−0.120−0.2000.1740.2400.2800.4000.480−0.5000.500
KIII0.3000.400−0.132−0.310−0.360−0.133−0.240−0.200−0.250
KIV−0.241−0.273−0.360−0.467−0.573−0.422−0.4930.200−0.500
YouyiKI−0.571−0.547−0.254−0.104−0.179−0.056−0.2220.300−0.250
KII−0.400−0.3200.5000.2400.2800.1250.400−0.7000.500
KIII0.0000.360−0.254−0.402−0.360−0.292−0.300−0.925−0.250
KIV0.000−0.209−0.398−0.567−0.573−0.528−0.533−0.970−0.500
MaanshiKI−0.314−0.255−0.204−0.423−0.097−0.1920.4400.000−0.250
KII−0.0400.4800.3350.2400.1200.375−0.4400.0000.500
KIII0.100−0.240−0.337−0.305−0.440−0.125−0.720−0.750−0.250
KIV−0.273−0.493−0.465−0.461−0.627−0.417−0.813−0.900−0.500
Table 8. Determination of comprehensive weights.
Table 8. Determination of comprehensive weights.
Assessment IndexesWeights
AHPImproved Entropy Weight MethodComprehensive Weights
Surrounding rock classification0.0260.1040.057
Degree of weathering0.0280.1050.059
Slope angle of tunnel entrance0.0560.1050.075
Maximum burial depth of tunnel0.0840.1100.094
Tunnel defect0.1430.1230.135
Tunnel width0.0600.1040.077
Unfavorable geological condition0.1740.1130.150
Length of shattered fault zone0.1920.1310.168
Earthquake intensity0.2370.1050.185
Table 9. Correlation degree of tunnel seismic damage risk.
Table 9. Correlation degree of tunnel seismic damage risk.
No.TunnelKI (p)KII (p)KIII (p)KIV (p)
1Longchi tunnel−0.558−0.2320.019−0.184
2Longxi tunnel−0.427−0.180−0.169−0.378
3Longdongzi tunnel−0.3280.166−0.173−0.344
4Youyi tunnel−0.1510.101−0.353−0.549
5Maanshi tunnel−0.0960.145−0.412−0.605
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chen, C.; Zou, W.; Geng, P.; Gu, W.; Yuan, F.; He, C. Study on Seismic Damage Risk Assessment of Mountain Tunnel Based on the Extension Theory. Sustainability 2023, 15, 5294. https://doi.org/10.3390/su15065294

AMA Style

Chen C, Zou W, Geng P, Gu W, Yuan F, He C. Study on Seismic Damage Risk Assessment of Mountain Tunnel Based on the Extension Theory. Sustainability. 2023; 15(6):5294. https://doi.org/10.3390/su15065294

Chicago/Turabian Style

Chen, Changjian, Wei Zou, Ping Geng, Wenqi Gu, Feiyun Yuan, and Chuan He. 2023. "Study on Seismic Damage Risk Assessment of Mountain Tunnel Based on the Extension Theory" Sustainability 15, no. 6: 5294. https://doi.org/10.3390/su15065294

APA Style

Chen, C., Zou, W., Geng, P., Gu, W., Yuan, F., & He, C. (2023). Study on Seismic Damage Risk Assessment of Mountain Tunnel Based on the Extension Theory. Sustainability, 15(6), 5294. https://doi.org/10.3390/su15065294

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