# Seismic Damage Rapid Assessment of Road Networks considering Individual Road Damage State and Reliability of Road Networks in Emergency Conditions

Authors: Jinlong Liu; Hanxi Jia; Junqi Lin; Heng Hu
Journal: Advances in Civil Engineering
 (2020)

---

## Abstract

Road networks are one of the vital components of a transportation system that influence the traffic capacity and disaster losses after the earthquakes. The road network reliability is crucial for the postearthquake emergency decision-making. In this study, a method is proposed to assess the seismic damage of road networks considering individual road damage state and reliability of road networks in emergency conditions. First, we assess the importance of the factors that affect seismic road damage using the AdaBoost algorithm. In addition, artificial neural networks are used to evaluate the damage state of an individual road based on the factors that are selected with higher importance. Then, the improved estimation for the reliability of road networks is adopted to assess the damage of road networks. Last, the method is demonstrated using the road networks in Karamay, China.


---

## Body

## 1. Introduction

Lifeline systems refer to the infrastructure networks, including transportation, communication, water supply, drainage, power supply, gas supply, and oil transportation, which have a major impact on social life and production [1]. The road networks are one of the critical components of the transportation system that transmit information in the lifeline systems. They are also the life passage for seismic relief after the earthquake. Rapidly and accurately assessing the damage state of the road networks after the earthquake can provide reference information for the deployment of emergency rescue works and rationally distribute the manpower and material resources for disaster relief. A timely rescue can reduce the casualties and property losses in the disaster areas.In recent years, the seismic damage evaluations of the transportation networks have been widely investigated by some studies. In these past studies, most of them regard bridge failure as the main factor affecting the capacity of the road network, while ignoring the impact of road failure on the overall performance of the road network [2]. Kiremidjian et al. proposed an assessment model of the transportation networks based on the loss evaluation of the bridges [3]. Guo et al. considered the roads as the main components of the transportation systems in their method, which were developed according to postdisaster traffic demand [4]. However, the focus of the research proposed by Guo et al. remained the bridges. Some researchers have focused studies on the road networks. The method proposed by Cheng et al. was based on the road network capacity (RNC) model and combined travel demand model, which was mainly calculated by the original-destination (OD) matrix [5]. It is only required to solve the shortest path traffic reliability between all OD sets and OD sets of the road networks [6]. The expected value of the traffic probability of all the shortest paths [7] is the connectivity reliability of the road network, which can be used to evaluate the overall operation of the road networks after the earthquake. The method can reflect the structure and evaluate the traffic capacity of the road network. However, as road network vertices increase, the calculation of the number of OD sets will increase exponentially [8], which will greatly increase the workload of assessment. Hence, a faster calculation method is needed to assess the road capacity after the earthquake.Furthermore, the damage assessments of the individual road were ignored in most previous studies [9–11]. Pitilakis et al. investigated the dynamic transportation performance directly based on the whole road networks [10]. Some studies are used to evaluate the road seismic damage based on the linear regression model [12]. The causes of the road damages are very complicated, and they are not only influenced by the earthquakes. However, the linear models are still flawed in solving the problem of very complicated internal mechanism, and they often need some assumptions. In recent years, the rises of machine learning algorithms have greatly advanced the development of the prediction problems [13]. For example, the ANN model was used to assess the bridge damage severity quantification [14] and buildings’ damage seismic state [15]. The ANN model also has almost no application in the road network seismic damage assessment. Before using the ANN model to predict, the features should be screened first rather than by empirically selecting the influencing factors [12]. One of the best ways to choose the factors that influence road damage is to assess the importance of them [16]. The method of ensemble learning algorithms has been proved to have a good performance in the analysis of feature importance [16–18]. There are many algorithms in the ensemble learning models, such as the most commonly used random forest algorithm [19] and adaptive boosting algorithm (AdaBoost) [20]. However, ensemble learning models have not been used to assess the importance of different influencing factors of seismic road damage.This paper shows a method for the seismic damage assessment of road networks considering the damage state of an individual road and reliability of road networks. The first part of this method consists of two models: (1) the assessment of factor importance applying the AdaBoost algorithms [20] and (2) the estimation of the damage state of an individual road using the ANN model based on the results of the AdaBoost model. The second part is to evaluate the reliability of road networks considering the capacity of road networks (RNC) and pass probability based on the results of the ANN model. Finally, we use a case study to verify the applicability of the method proposed in this study, and it is proved that the method has a good performance on the problem of the seismic damage assessment of road networks.


## 2. Data

Because the data of the road are less than the bridge in an earthquake, some studies have to make some assumptions so as to verify the efficiency of their method. The data used in this article are derived from actual seismic damage and do not contain any assumptions. The road damage data were collected by the government and some institutions (Sichuan Highway Planning, Survey and Design Institute, China Highway Planning and Design Institute, and Shanxi Highway Planning Survey and Design Institute). They were all derived from the 2008 Ms 8.0 Wenchuan earthquake in China (Table1). In Table 1, the seismic fortification intensity and soil type refer to the “Chinese Specification of Seismic Design for Highway Engineering” (JTG B02-2013) [21], and the practical intensity refers to the Chinese seismic intensity scale (GBT 17742-2008) [22].Table 1
The damage roads of the Wenchuan earthquake with some characteristics; the number consists of the road name and the mileage number.
NumberHeight difference (m)Road gradesSupporting and retaining typeSeismic fortification intensityPractical intensity (degree)Slope height (m)Soil type of the siteG213: K1020 + 96010.43Gravity-type groutedVIIXI13.6IIG213: K1018 + 6001.23Gravity-type groutedVIIXI40.0IG213: K1022 + 90015.63Reinforced retaining wallVIIXI0IG213: K1029 + 8503.13Gravity-type groutedVIIX33.0IG213: K1008 + 9003.03Gravity-type groutedVIIXI20.0IG213: K1014 + 1756.03Gravity-type groutedVIIXI0IG213: K1008 + 4002.03Cutting slopeVIIXI50.0IIIG213: K1008 + 5809.93Gravity-type concreteVIIXI50.0IIIG213: K1012 + 40013.63Gravity-type concreteVIIXI15.5IIIG213: K1029 + 7006.23Gravity-type concreteVIIX50.0IIIG213: K1023 + 7002.23Prestressed cableVIIXI19.0IIIG213: K1009 + 0807.13Cutting slopeVIIXI0IIIG213: K1008 + 9804.03Gravity-type concreteVIIXI0IIIG213: K916 + 5083.63Cutting slopeVIIIX7.0IIIG213: K37 + 3501.82Gravity-type groutedVIIX294.0IIG213: K38 + 5001.02Gravity-type groutedVIIX326.3IIG213: K40 + 1003.62Gravity-type groutedVIIX121.2IG213: K42 + 4004.42Facing wallVIIX175.5IIG213: K26 + 8001.52Cutting slopeVIIXI24.3IG213: K73 + 0002.02Cutting slopeVIIIX175.9IG213: K29 + 9501.02Facing wallVIIXI19.7IG213: K34 + 1401.52Cutting slopeVIIXI122.7IIIG213: K35 + 6003.92Cutting slopeVIIXI100.0IIIG213: K50 + 2003.02Cutting slopeVIIX480.1IIIG213: K58 + 1603.42Gravity-type groutedVIIX15.0IIIG2131.53Gravity-type groutedVIIIX7.1IIIG2131.73Gravity-type groutedVIIIX92.1IIG2134.03Gravity-type groutedVIIIX27.0IIG2133.43Cutting slopeVIIIX93.4IIIG2132.33Facing wallVIIIX71.0IIIG2131.53Gravity-type groutedVIIIX0IIG21313Gravity-type groutedVIIIX10.0IIG2131.23Cutting slopeVIIIX0IIS303: K10 + 3411.82Gravity-type groutedVIIX5.0IIIS303: K28 + 9410.52Gravity-type groutedVIIX50.0IIS303: K40 + 5000.52Gravity-type groutedVIIIX27.4IIIS303: K16 + 4412.62Cutting slopeVIIX237.5IIS303: K31 + 8411.02Cutting slopeVIIIX28.7IS303: K36 + 6001.02Cutting slopeVIIIX100.0IIIS303: K43 + 2001.02Cutting slopeVIIIX30.0IIS303: K26 + 9003.12Gravity-type groutedVIIIX25.0IIIS3022.02Cutting slopeVIIIX27.0IIIS302: K732 + 5601.52Gravity-type groutedVIIX7.4IIIS3022.52Facing wallVIIX47.5IS1052.23Gravity-type groutedVIIXI2.8IS1053.13Gravity-type groutedVIIXI25.3IIIS1053.53Facing wallVIIXI36.7IS105: K167 + 4605.93Cutting slopeVIIXI8.4IIIS105: K168 + 9007.43Cutting slopeVIIXI9.0IIIXU092.34Gravity-type concreteVIX44.3IIS2101.33Gravity-type groutedVIIVII324.4IS2103.73Cutting slopeVIIVII63.6IIS2102.33Cutting slopeVIIVII24.4IIIG2124.23Gravity-type groutedVIIVIII34.3IIG2131.02Gravity-type groutedVIIVIII20.3IS3033.84Cutting slopeVIIVII194.3IIIS3034.24Gravity-type groutedVIIVIII35.0IIIG3171.53Cutting slopeVIIVII173.3IIG3172.52Gravity-type groutedVIIVI2.5IIG3172.42Gravity-type groutedVIIVI5.0IG3170.52NoVIIVI0IIIG3176.02Gravity-type groutedVIIVI0IIIThere were five national highways (G108, G212, G213, G317, and G318) and eleven provincial highways (S101, S105, S106, S202, S205, S210, S211, S301, S302, S303, and S306) in the disaster-stricken areas of Sichuan province damaged to varying degrees. National road G213, provincial road S302, and provincial road S303 crossing the area of earthquake intensity VIII were seriously damaged, and the total damage of rural roads was about 24,103 kilometers. There were two national highways (G212 and G316) and nine provincial highways (S205, S206, S208, S219, S306, S307, S313, X482, and X484) in Gansu province which were damaged to varying degrees, and the total damage was about 5518 kilometers including rural roads. There were only one national highway (G108) and two provincial highways (S210 and S309) in Shanxi province which were damaged to varying degrees, and the total damage was about 1,791 kilometers including rural roads. The above roads were divided into 63 links based on the damage and mileage. Figure1 demonstrates some typical damages of the road in the Wenchuan earthquake, which were divided into some road segments and had the individual number [24].Figure 1
The typical damage of the road in the Wenchuan earthquake, including the retaining wall being cut, the wall body being inclined outward, the subgrade being buried by the rock, and the subgrade partially collapsed [23].


## 3. Assessment of Seismic Road Damage

### 3.1. Select Factors

There are many factors affecting road seismic damage, but not every factor is critical. In the rapid assessment of road damage, the most important factors should be selected for evaluation, which can save time in obtaining data and running procedures.The road is mainly composed of subgrade and pavement, retaining structure, and slope. The road damage is divided into three types based on the failure parts: subgrade damage, supporting structure damage, and slope damage. The seismic damage data show that the three seismic types are not independent. For example, the retaining structure or subgrade will be destroyed when the slope collapses. The following is a detailed description of the selection factors (Figure2) (Tables 2 and 3):Figure 2
Diagrammatic drawing of the subgrade type. (a) Embankment. (b) Cutting. (c) Half-filled and half-dug. (d) No digging and no filling.
(a)(b)(c)(d)Table 2
Classification of engineering sites.
Mean shear wave velocity (m/s)Site categoryIIIIIIIVvse>5000———500≥vse>250<5≥5——250≥vse>140<3≥3,≤50>50—vse≤140<3≥3,≤15>15,≤80>80Table 3
Classification of slope protection types.
ClassificationSlope protection typeLevel INo protection, wall protection, plant protection, spraying (no net)Level IINet hanging, shotcrete, active and passive protection netLevel IIIFrame beam, prestressed anchor (cable) protection, composite measures(1)
Different road grades have different importance in seismic design. According to the Chinese Design Code of Highway Subgrades (JTG D30-2015) [25], the road can be divided into 4 grades. For each grade, the road has different materials and construction methods, and so, the earthquake resistance ability of each grade road is different.(2)
Subgrade-related factors included the subgrade type and height difference of the subgrade [26]. The type of subgrade was classified according to the form of excavation and filling. The subgrade height difference is the difference between the design elevation of the subgrade centerline and the original ground elevation, which was used to describe the stability of the subgrade in the earthquake. The diagrammatic drawing of the subgrade type is shown in Figure 2.(3)
The retaining structures are used to support and stabilize the subgrade filling or slope to prevent soil slipping and thus to maintain soil stability [27]. The retaining structure can be divided into gravity retaining wall and flexible retaining wall. Gravity retaining wall can be divided into three types according to the construction method: dry masonry, wet masonry, and concrete. The seismic damage data show that the dry masonry retaining wall is more prone to serious damage, and its actual seismic capacity is obviously weaker than that of wet and concrete, so the dry masonry retaining wall and the nonretaining wall are combined into one category. The reinforced retaining wall, antislide pile, and anchor cable structure with good seismic performance are classified as advantageous retaining, so the retaining types are divided into four categories: nonretaining, wet masonry gravity, concrete gravity, and advantageous retaining.(4)
The site soil has a great influence on the bearing capacity of the foundation [28]. In the Wenchuan earthquake, a large number of shoulder walls were inclined and deformed due to insufficient bearing capacity of the foundation. According to the Specification of Seismic Design for Highway Engineering (JTG B02-2013) [21], the engineering site is divided into four categories, as shown in Table 2.(5)
Slope-related factors select slope angle, height, and protective measures. As the slope height and slope increase, the number of earthquake damage increases simultaneously. The higher the slope is, the more obvious the acceleration amplification effect is under the action of the earthquake, which is easy to produce earthquake damage such as rock projectile. Similarly, the steeper the slope is, the worse the stability will be. In the Wenchuan earthquake, the slope damage mainly concentrates on the slope between 35° and 65°. In addition, the characteristics of seismic damage show that the seismic damage of the slope mainly occurs on the soil slope without protection, while the failure of the slope with the combination of the frame, anchor bolt, and multimeasures is less. This indicates that appropriate protective measures can effectively reduce the degree and quantity of damage to the slope. According to the slope protection method, slope protection can be divided into three categories, as shown in Table3; from type I to type III, the seismic protection capability of the slope is gradually enhanced.(6)
The seismic intensity describes the intensity of the earthquake. The damage data of the Wenchuan earthquake also show that the amount and extent of damage increase with the increase of intensity, which is consistent with our understanding. The seismic fortification level is an important factor of road’s seismic resistance ability. In the engineering design stage, the seismic resistance of the road will be guaranteed according to the corresponding specifications [21]. Therefore, the fortification intensity can be taken as one of the evaluation indexes.Above all, we estimated the importance of ten factors as follows: highway classification, subgrade type, supporting and retaining type, soil type of the site, slope height, slope angle, protective measure, height difference of subgrade, seismic fortification intensity, and practical intensity (Table4).Table 4
Factors of seismic road damage.
FactorsClassificationFeaturesHighway classificationExpresswayBetter seismic performanceFirst-classBetter seismic performanceSecond-classThird-classFailure occurs more frequentlyFourth-classFailure occurs more frequentlySubgrade typeNo digging and no fillingNo damageCuttingDestruction accounts for 26% of the total damageEmbankmentDestruction accounts for 28% of the total damageHalf-filled and half-dugFailure occurs most frequently (48%)Supporting and retaining typeReinforced retaining wall, prestressed cable, facing wallBeneficial to mitigating damageGravity-type concreteFailure occurs most frequentlyGravity-type wet masonryMore serious damageNoNot conducive to mitigating damageSoil type of the siteIIIIIIIntensifying the damage slightlyIVIntensifying the damageSlope height (h)0 mThe least number of damages occurred0 m <h ≤ 15Less destruction occurs15 m <h ≤ 40 m40 m <h ≤ 120 mMore susceptible to damage120 m <hMore susceptible to damageSlope angle (θ)0°The least number of damages occurred0° <θ ≤ 35°35° <θ ≤ 65°Failure occurs most frequently65° <θMore susceptible to damageProtective measureNo, steening retaining wallNot conducive to resisting seismic, small quantityFrame beam, precast blockBeneficial to seismicHanging net and guniting, safety netting systemCan resist part of seismicShotcrete, mortarNot conducive to resisting seismicHeight difference of subgradeh ≤ 1 m1 m <h ≤ 3 m3 m <h ≤ 5 m5 m <hSeismic fortification intensityVIVIIVIIIIXPractical intensityVIVIIVIIIIXXXI


### 3.2. Assessment of the Factor Importance Applying the AdaBoost Algorithm

There are many algorithms for assessing the importance of features in the machine learning domain, such as support vector machine, AdaBoost algorithm, and decision tree model. This study selected the AdaBoost algorithm in ensemble learning. The essence of the ensemble learning model is to combine many weak processors to get a better predictive processor. Ensemble learning model works well for large data set and insufficient data set. The data in this study are not sufficient. Common ensemble learning models are of two types: bagging algorithms and boosting algorithms. The boosting algorithm predicts through a series of aggregated estimated model weighted averages. AdaBoost is a representative algorithm in the boosting algorithm [29]. In this study, the input parameters were ten factors and five damage states. The damage states of individual roads are divided into five levels: none, slight, moderate, extensive, and complete [30].The procedures of the AdaBoost algorithm can be summarized as follows:(i)
Step 1: initialize the sample weights and perform equal weight processing.(ii)
Step 2: train the base classifier. The weights are updated according to the results of each decision tree and then trained until the conditions are met.(iii)
Step 3: the weighted average method is used to combine the base classifiers into strong classifiers, and the decision trees with smaller errors are more weighted.Base classifier: classification and regression tree classifier (maximum depth = 10, minimum sample leaf = 5, algorithm = SAMME.R, number of estimators = 200, and learning rate = 0.5).Due to the small data set, we chose the cross validation function in the sklearn function library to validate the model rather than dividing the data set into training set and testing set. The mean average accuracy of the output classification was 0.8204, which was higher than the accuracy of the AdaBoost algorithm in other models [16]. Figure 3 presents the importance of the ten factors, and there are five factors that have the importance value more than 10%.Figure 3
The importance values of factors of the seismic road damage based on the AdaBoost algorithm.


### 3.3. Artificial Neural Network Assessment Model

In the rapid assessment after the earthquake, the assessment speed is one of the most important points. Reducing the parameters can reduce the amount of data and speed up the evaluation. Therefore, we used the AdaBoost algorithm to evaluate the importance of different influencing factors. According to the results of Section3.2, we selected the five most important factors as input parameters. We chose practical intensity, supporting and retaining type, height difference of subgrade, slope height, and soil type of the site as the input parameters. The sum of the importance value of these five parameters reached 80%. Hence, we ignore other parameters when establishing the evaluation model. ANN has strong adaptive, self-learning, and nonlinear mapping capabilities, which can solve the problem of less data and uncertainty, and is not limited by nonlinear models. Earthquake is an accidental event, and destructive earthquakes have a lower probability of occurrence, such as the Wenchuan earthquake. Therefore, there are small data set on road damages and great uncertainties. Moreover, the relationship between various influencing factors of road damage and the results of damage is not clear. Above all, artificial neural networks are suitable for road damage assessment. A typical neural network includes an input layer, an implicit layer, and an output layer. The layers are fully connected, and there is no connection between the layers. The hidden layer can have one or more layers. For a typical network, a single layer of hidden layers is sufficient [13]. The input and output layers have been determined. The hidden layer selects a layer based on the actual amount of data. The number of hidden layer neurons is determined according to the following empirical equation:(1)L=n+m+a,where L is the number of hidden layers, n and m are the number of input and output layers, respectively, and a is a constant between 1 and 10. n = 10 and m = 1. We determined the number of hidden layer neurons as ten according to the convergence speed and accuracy during the training process. Therefore, this study selected the most typical three-layer backpropagation neural network, with 5 input layer neurons, 10 hidden layer neurons, and 1 output layer neuron (Figure 4). The maximum epochs were 10,000, and the learning rate was 0.1 based on the running processing and the data set. The mean square error loss function and sigmoid function were selected as the cost function (2) and activation function (3) of the hidden layer, respectively:(2)y=1n∑y^−t2,where y is the value of the cost function, n presents the amount of the roads, y^ is the damage state of the ANN model, and t is the actual damage state.(3)fx=11+e−z,where fx is the ReLU function, and x is the input of the ReLU function.Figure 4
Layout for artificial neural networks for the assessment model.We selected 50 roads as the training set and the remaining 12 roads as the testing set. The model was run in the Jupyter Notebook of the Anaconda Navigator environment with python 3 languages (https://www.anaconda.com/). The accuracy of the linear model was 87% [8], while the accuracy of the ANN model was 92%. It can be seen that the ANN model had a good performance in assessing road damage problems.


## 3.1. Select Factors

There are many factors affecting road seismic damage, but not every factor is critical. In the rapid assessment of road damage, the most important factors should be selected for evaluation, which can save time in obtaining data and running procedures.The road is mainly composed of subgrade and pavement, retaining structure, and slope. The road damage is divided into three types based on the failure parts: subgrade damage, supporting structure damage, and slope damage. The seismic damage data show that the three seismic types are not independent. For example, the retaining structure or subgrade will be destroyed when the slope collapses. The following is a detailed description of the selection factors (Figure2) (Tables 2 and 3):Figure 2
Diagrammatic drawing of the subgrade type. (a) Embankment. (b) Cutting. (c) Half-filled and half-dug. (d) No digging and no filling.
(a)(b)(c)(d)Table 2
Classification of engineering sites.
Mean shear wave velocity (m/s)Site categoryIIIIIIIVvse>5000———500≥vse>250<5≥5——250≥vse>140<3≥3,≤50>50—vse≤140<3≥3,≤15>15,≤80>80Table 3
Classification of slope protection types.
ClassificationSlope protection typeLevel INo protection, wall protection, plant protection, spraying (no net)Level IINet hanging, shotcrete, active and passive protection netLevel IIIFrame beam, prestressed anchor (cable) protection, composite measures(1)
Different road grades have different importance in seismic design. According to the Chinese Design Code of Highway Subgrades (JTG D30-2015) [25], the road can be divided into 4 grades. For each grade, the road has different materials and construction methods, and so, the earthquake resistance ability of each grade road is different.(2)
Subgrade-related factors included the subgrade type and height difference of the subgrade [26]. The type of subgrade was classified according to the form of excavation and filling. The subgrade height difference is the difference between the design elevation of the subgrade centerline and the original ground elevation, which was used to describe the stability of the subgrade in the earthquake. The diagrammatic drawing of the subgrade type is shown in Figure 2.(3)
The retaining structures are used to support and stabilize the subgrade filling or slope to prevent soil slipping and thus to maintain soil stability [27]. The retaining structure can be divided into gravity retaining wall and flexible retaining wall. Gravity retaining wall can be divided into three types according to the construction method: dry masonry, wet masonry, and concrete. The seismic damage data show that the dry masonry retaining wall is more prone to serious damage, and its actual seismic capacity is obviously weaker than that of wet and concrete, so the dry masonry retaining wall and the nonretaining wall are combined into one category. The reinforced retaining wall, antislide pile, and anchor cable structure with good seismic performance are classified as advantageous retaining, so the retaining types are divided into four categories: nonretaining, wet masonry gravity, concrete gravity, and advantageous retaining.(4)
The site soil has a great influence on the bearing capacity of the foundation [28]. In the Wenchuan earthquake, a large number of shoulder walls were inclined and deformed due to insufficient bearing capacity of the foundation. According to the Specification of Seismic Design for Highway Engineering (JTG B02-2013) [21], the engineering site is divided into four categories, as shown in Table 2.(5)
Slope-related factors select slope angle, height, and protective measures. As the slope height and slope increase, the number of earthquake damage increases simultaneously. The higher the slope is, the more obvious the acceleration amplification effect is under the action of the earthquake, which is easy to produce earthquake damage such as rock projectile. Similarly, the steeper the slope is, the worse the stability will be. In the Wenchuan earthquake, the slope damage mainly concentrates on the slope between 35° and 65°. In addition, the characteristics of seismic damage show that the seismic damage of the slope mainly occurs on the soil slope without protection, while the failure of the slope with the combination of the frame, anchor bolt, and multimeasures is less. This indicates that appropriate protective measures can effectively reduce the degree and quantity of damage to the slope. According to the slope protection method, slope protection can be divided into three categories, as shown in Table3; from type I to type III, the seismic protection capability of the slope is gradually enhanced.(6)
The seismic intensity describes the intensity of the earthquake. The damage data of the Wenchuan earthquake also show that the amount and extent of damage increase with the increase of intensity, which is consistent with our understanding. The seismic fortification level is an important factor of road’s seismic resistance ability. In the engineering design stage, the seismic resistance of the road will be guaranteed according to the corresponding specifications [21]. Therefore, the fortification intensity can be taken as one of the evaluation indexes.Above all, we estimated the importance of ten factors as follows: highway classification, subgrade type, supporting and retaining type, soil type of the site, slope height, slope angle, protective measure, height difference of subgrade, seismic fortification intensity, and practical intensity (Table4).Table 4
Factors of seismic road damage.
FactorsClassificationFeaturesHighway classificationExpresswayBetter seismic performanceFirst-classBetter seismic performanceSecond-classThird-classFailure occurs more frequentlyFourth-classFailure occurs more frequentlySubgrade typeNo digging and no fillingNo damageCuttingDestruction accounts for 26% of the total damageEmbankmentDestruction accounts for 28% of the total damageHalf-filled and half-dugFailure occurs most frequently (48%)Supporting and retaining typeReinforced retaining wall, prestressed cable, facing wallBeneficial to mitigating damageGravity-type concreteFailure occurs most frequentlyGravity-type wet masonryMore serious damageNoNot conducive to mitigating damageSoil type of the siteIIIIIIIntensifying the damage slightlyIVIntensifying the damageSlope height (h)0 mThe least number of damages occurred0 m <h ≤ 15Less destruction occurs15 m <h ≤ 40 m40 m <h ≤ 120 mMore susceptible to damage120 m <hMore susceptible to damageSlope angle (θ)0°The least number of damages occurred0° <θ ≤ 35°35° <θ ≤ 65°Failure occurs most frequently65° <θMore susceptible to damageProtective measureNo, steening retaining wallNot conducive to resisting seismic, small quantityFrame beam, precast blockBeneficial to seismicHanging net and guniting, safety netting systemCan resist part of seismicShotcrete, mortarNot conducive to resisting seismicHeight difference of subgradeh ≤ 1 m1 m <h ≤ 3 m3 m <h ≤ 5 m5 m <hSeismic fortification intensityVIVIIVIIIIXPractical intensityVIVIIVIIIIXXXI


## 3.2. Assessment of the Factor Importance Applying the AdaBoost Algorithm

There are many algorithms for assessing the importance of features in the machine learning domain, such as support vector machine, AdaBoost algorithm, and decision tree model. This study selected the AdaBoost algorithm in ensemble learning. The essence of the ensemble learning model is to combine many weak processors to get a better predictive processor. Ensemble learning model works well for large data set and insufficient data set. The data in this study are not sufficient. Common ensemble learning models are of two types: bagging algorithms and boosting algorithms. The boosting algorithm predicts through a series of aggregated estimated model weighted averages. AdaBoost is a representative algorithm in the boosting algorithm [29]. In this study, the input parameters were ten factors and five damage states. The damage states of individual roads are divided into five levels: none, slight, moderate, extensive, and complete [30].The procedures of the AdaBoost algorithm can be summarized as follows:(i)
Step 1: initialize the sample weights and perform equal weight processing.(ii)
Step 2: train the base classifier. The weights are updated according to the results of each decision tree and then trained until the conditions are met.(iii)
Step 3: the weighted average method is used to combine the base classifiers into strong classifiers, and the decision trees with smaller errors are more weighted.Base classifier: classification and regression tree classifier (maximum depth = 10, minimum sample leaf = 5, algorithm = SAMME.R, number of estimators = 200, and learning rate = 0.5).Due to the small data set, we chose the cross validation function in the sklearn function library to validate the model rather than dividing the data set into training set and testing set. The mean average accuracy of the output classification was 0.8204, which was higher than the accuracy of the AdaBoost algorithm in other models [16]. Figure 3 presents the importance of the ten factors, and there are five factors that have the importance value more than 10%.Figure 3
The importance values of factors of the seismic road damage based on the AdaBoost algorithm.


## 3.3. Artificial Neural Network Assessment Model

In the rapid assessment after the earthquake, the assessment speed is one of the most important points. Reducing the parameters can reduce the amount of data and speed up the evaluation. Therefore, we used the AdaBoost algorithm to evaluate the importance of different influencing factors. According to the results of Section3.2, we selected the five most important factors as input parameters. We chose practical intensity, supporting and retaining type, height difference of subgrade, slope height, and soil type of the site as the input parameters. The sum of the importance value of these five parameters reached 80%. Hence, we ignore other parameters when establishing the evaluation model. ANN has strong adaptive, self-learning, and nonlinear mapping capabilities, which can solve the problem of less data and uncertainty, and is not limited by nonlinear models. Earthquake is an accidental event, and destructive earthquakes have a lower probability of occurrence, such as the Wenchuan earthquake. Therefore, there are small data set on road damages and great uncertainties. Moreover, the relationship between various influencing factors of road damage and the results of damage is not clear. Above all, artificial neural networks are suitable for road damage assessment. A typical neural network includes an input layer, an implicit layer, and an output layer. The layers are fully connected, and there is no connection between the layers. The hidden layer can have one or more layers. For a typical network, a single layer of hidden layers is sufficient [13]. The input and output layers have been determined. The hidden layer selects a layer based on the actual amount of data. The number of hidden layer neurons is determined according to the following empirical equation:(1)L=n+m+a,where L is the number of hidden layers, n and m are the number of input and output layers, respectively, and a is a constant between 1 and 10. n = 10 and m = 1. We determined the number of hidden layer neurons as ten according to the convergence speed and accuracy during the training process. Therefore, this study selected the most typical three-layer backpropagation neural network, with 5 input layer neurons, 10 hidden layer neurons, and 1 output layer neuron (Figure 4). The maximum epochs were 10,000, and the learning rate was 0.1 based on the running processing and the data set. The mean square error loss function and sigmoid function were selected as the cost function (2) and activation function (3) of the hidden layer, respectively:(2)y=1n∑y^−t2,where y is the value of the cost function, n presents the amount of the roads, y^ is the damage state of the ANN model, and t is the actual damage state.(3)fx=11+e−z,where fx is the ReLU function, and x is the input of the ReLU function.Figure 4
Layout for artificial neural networks for the assessment model.We selected 50 roads as the training set and the remaining 12 roads as the testing set. The model was run in the Jupyter Notebook of the Anaconda Navigator environment with python 3 languages (https://www.anaconda.com/). The accuracy of the linear model was 87% [8], while the accuracy of the ANN model was 92%. It can be seen that the ANN model had a good performance in assessing road damage problems.


## 4. Analysis on the Reliability of Road Networks

According to the results of the ANN model, we could obtain the damage state of each individual road. However, it is also necessary to consider the damage of the road network so that the road could not affect the emergency rescue after the earthquake. The combination of the damage state on the individual road and reliability of road networks can ensure that the lifeline system of emergency rescue after the earthquake is unimpeded. Therefore, we proposed the postearthquake evaluation method of the reliability of road networks.


### 4.1. The Pass Possibility of an Individual Road

A strong earthquake caused damage to the roads. Earthquake damage index is usually used to describe the damage degree of road, but the damage degree of earthquake is only the description of the physical damage state of road, and it cannot reflect the status of road traffic function. Therefore, it is necessary to establish the relationship between seismic damage index and pass probability. The results of Section3.3 gave only the damage state and could not directly obtain the pass probability, Pr,i, of an individual road. Therefore, we used the damage index [31] as the conversion indicator to translate the damage state into pass possibility. Table 5 shows the relationship between the damage index and the damage state. The damage index and pass probability had a certain range, and we decided to use the median value of them as the value of the damage index and Pr,i.Table 5
The relationship between the damage index and the damage state.
Damage stateDamage index (I)Pass statePass probability (p)None0.00≤I<0.10Allowable0.90<p≤1.00Slight0.10≤I<0.30Limitation0.70<p≤0.90Moderate0.30≤I<0.55Limitation0.45<p≤0.70Extensive0.55≤I<0.85Unallowedp≤0.001Complete0.85≤I≤1.00Unallowedp≤0.001


### 4.2. Reliability of the Road Network Model

According to the results of the ANN model, we could obtain the damage state of each individual road. However, it is also necessary to consider the damage of the road network so that the road could not affect the emergency rescue after the earthquake. Therefore, we proposed the postearthquake evaluation method of the reliability of road networks. Urban RNC is an important decision indicator in urban road network planning and construction [4]. The RNC represents the maximum traffic volume and provides important information for effective traffic flow control and demand management. The RNC is the maximum flow rate determined according to the geometric characteristics of the road, traffic conditions, and specified operational characteristics, and its value has relative stability and prescriptiveness. The simplest consideration for RNC is to add up the maximum capacity of all road links to get the physical capacity of an urban transportation network. Under the general road and traffic conditions of the city and when not affected by the intersection of planes, the possible RNC, NP, can be computed by [32](4)NP=3600δi,where δi is the average traffic interval between consecutive traffic. When the city does not have the value of δi, the road capacity can use the value in Table 6.Table 6
Possible capacity of the road [32].
Calculated driving speed (km/h)50403020NP (pcu/h)1690164015501380According to the code for design of urban road engineering, in road capacity design, the actual capacity of the road is not allowed to reach the maximum capacity value. Instead, it tends to reduce the maximum capacity as the design value for safety. The higher the level of the road, the greater the reduction coefficient. When the RNC is not affected by a plane intersection, the design RNC,Ni0, can be computed by(5)Ni0=ac×NP,where ac is the road classification coefficient of the vehicle lane (Table 7).Table 7
Coefficient of the road classification [24].
Road classificationExpresswayArterial roadMinor arterial roadBranch roadac0.750.800.850.90The initial RNC,N0, before the earthquake is defined as follows:(6)N0=∑i=1nZi×Ni0,where Zi is the importance index of the i-th road, which mainly reflects the importance or influence of the road in the network. The urban road network consists of several road sections, but each road section has different importance in the whole road network. Some road sections have high utilization rate, and some are very low. The impact of these unit sections on the connectivity of the entire road network is not the same. Therefore, the weight of the unit road link needs to be calibrated. The index Zi is usually calculated based on the betweenness centrality. However, the betweenness centrality involves calculation in a complex network [33]. Hence, the index Zi can be expressed as (7) for the convenience of calculation and speed of the assessment:(7)Zi=BiBtd,where Bi is the number of road links connected to the i-th road link and Btd is the total number of road links. If the road link contains bridges, tunnels, etc., the pass probability Pr,i is obtained based on the ANN model, and the RNC after the earthquake can be obtained by (8). The reliability of the road network, R, is defined as the ratio of N to N0 (9). Table 8 presents the reliability state of the road network, including reliable, medium reliable, and unreliable.(8)N=∑i=1nPr,i×Zi×Ni0,(9)R=NN0.Table 8
The reliability state of the road network.
Reliability stateReliableMedium reliableUnreliableReliability of the road networkR≥0.80.5≤R<0.8R<0.5


## 4.1. The Pass Possibility of an Individual Road

A strong earthquake caused damage to the roads. Earthquake damage index is usually used to describe the damage degree of road, but the damage degree of earthquake is only the description of the physical damage state of road, and it cannot reflect the status of road traffic function. Therefore, it is necessary to establish the relationship between seismic damage index and pass probability. The results of Section3.3 gave only the damage state and could not directly obtain the pass probability, Pr,i, of an individual road. Therefore, we used the damage index [31] as the conversion indicator to translate the damage state into pass possibility. Table 5 shows the relationship between the damage index and the damage state. The damage index and pass probability had a certain range, and we decided to use the median value of them as the value of the damage index and Pr,i.Table 5
The relationship between the damage index and the damage state.
Damage stateDamage index (I)Pass statePass probability (p)None0.00≤I<0.10Allowable0.90<p≤1.00Slight0.10≤I<0.30Limitation0.70<p≤0.90Moderate0.30≤I<0.55Limitation0.45<p≤0.70Extensive0.55≤I<0.85Unallowedp≤0.001Complete0.85≤I≤1.00Unallowedp≤0.001


## 4.2. Reliability of the Road Network Model

According to the results of the ANN model, we could obtain the damage state of each individual road. However, it is also necessary to consider the damage of the road network so that the road could not affect the emergency rescue after the earthquake. Therefore, we proposed the postearthquake evaluation method of the reliability of road networks. Urban RNC is an important decision indicator in urban road network planning and construction [4]. The RNC represents the maximum traffic volume and provides important information for effective traffic flow control and demand management. The RNC is the maximum flow rate determined according to the geometric characteristics of the road, traffic conditions, and specified operational characteristics, and its value has relative stability and prescriptiveness. The simplest consideration for RNC is to add up the maximum capacity of all road links to get the physical capacity of an urban transportation network. Under the general road and traffic conditions of the city and when not affected by the intersection of planes, the possible RNC, NP, can be computed by [32](4)NP=3600δi,where δi is the average traffic interval between consecutive traffic. When the city does not have the value of δi, the road capacity can use the value in Table 6.Table 6
Possible capacity of the road [32].
Calculated driving speed (km/h)50403020NP (pcu/h)1690164015501380According to the code for design of urban road engineering, in road capacity design, the actual capacity of the road is not allowed to reach the maximum capacity value. Instead, it tends to reduce the maximum capacity as the design value for safety. The higher the level of the road, the greater the reduction coefficient. When the RNC is not affected by a plane intersection, the design RNC,Ni0, can be computed by(5)Ni0=ac×NP,where ac is the road classification coefficient of the vehicle lane (Table 7).Table 7
Coefficient of the road classification [24].
Road classificationExpresswayArterial roadMinor arterial roadBranch roadac0.750.800.850.90The initial RNC,N0, before the earthquake is defined as follows:(6)N0=∑i=1nZi×Ni0,where Zi is the importance index of the i-th road, which mainly reflects the importance or influence of the road in the network. The urban road network consists of several road sections, but each road section has different importance in the whole road network. Some road sections have high utilization rate, and some are very low. The impact of these unit sections on the connectivity of the entire road network is not the same. Therefore, the weight of the unit road link needs to be calibrated. The index Zi is usually calculated based on the betweenness centrality. However, the betweenness centrality involves calculation in a complex network [33]. Hence, the index Zi can be expressed as (7) for the convenience of calculation and speed of the assessment:(7)Zi=BiBtd,where Bi is the number of road links connected to the i-th road link and Btd is the total number of road links. If the road link contains bridges, tunnels, etc., the pass probability Pr,i is obtained based on the ANN model, and the RNC after the earthquake can be obtained by (8). The reliability of the road network, R, is defined as the ratio of N to N0 (9). Table 8 presents the reliability state of the road network, including reliable, medium reliable, and unreliable.(8)N=∑i=1nPr,i×Zi×Ni0,(9)R=NN0.Table 8
The reliability state of the road network.
Reliability stateReliableMedium reliableUnreliableReliability of the road networkR≥0.80.5≤R<0.8R<0.5


## 5. Case Study Results

The methodology proposed in this paper is applied to the Karamay district in Karamay city, China. Due to the lack of information on the roads of all cities and for the propose of demonstrating the method proposed in this study, the reliability of road networks is applied to the area enclosed by the Zhunger road, Xihuan road, Nanhuan road, and Donghuan road, as shown in Figure5. The road network contains 64 road links and 39 nodes, totaling 60 km (Table 9).Figure 5
Road network plan.Table 9
Properties of the road network.
Road linksVertexiVertexjNameRoad classificationLength (m)Year built1139Xihuan roadArterial road1112200021739Kunlun roadArterial road682200333839Xihuan roadArterial road697200041638Xingfu roadArterial road701200253738Xihuan roadArterial road1625200061537Nanxin roadArterial road903200273637Xihuan roadArterial road864200081436Shiji roadArterial road942200391415Youjian roadMinor arterial road8462003101314Shiji roadArterial road685200311415Nanxin roadArterial road6752002121516Youjian roadMinor arterial road1190200313516Xingfu roadArterial road6512002141617Youjian roadMinor arterial road670200315617Kunlun roadArterial road6812003161718Youjian roadMinor arterial road1108200317118Zhunger roadArterial road6802003181819Zhunger roadArterial road700200319619Youyi roadArterial road110620032026Kunlun roadArterial road91420032156Youyi roadArterial road64020032235Xingfu roadArterial road96420022345Yingbin roadArterial road884200024411Nanxin roadArterial road802200225413Yingbin roadArterial road8332000261213Shiji roadArterial road9292003271330Yingbin roadArterial road7172000283035Yingbin roadArterial road9552000293536Xihuan roadArterial road23222000302930Baoshi roadMinor arterial road11242012311920Zhunger roadArterial road916200332720Shengli roadArterial road59220033378Ashan roadMinor arterial road70220043427Shengli roadArterial road51220033529Kunlun roadArterial road80220033623Shengli roadArterial road515200337211Shengli roadArterial road8432003381011Nanxin roadArterial road7662002391112Shengli roadArterial road6962003401229Shengli roadArterial road7362003412934Shengli roadArterial road13092003423435Nanhuan roadArterial road16332012432729Baoshi roadMinor arterial road889201244910Dongjiao roadMinor arterial road11712004451026Nanxin roadArterial road70720024689Dongjiao roadMinor arterial road5122004472021Zhunger roadArterial road10132003482122Changzheng roadArterial road5992004492224Jinyuan roadArterial road754200050923Kunlun roadArterial road3702003512326Baihua roadMinor arterial road13382004522526Nanxin roadArterial road8222002532627Baihua roadMinor arterial road8202004542733Ruyi roadMinor arterial road11602004553334Nanhuan roadArterial road7242012562728Baoshi roadMinor arterial road12442012573233Nanhuan roadArterial road15032012582425Jinyuan roadArterial road10092000592528Jinyuan roadArterial road10092000602832Jinyuan roadArterial road6082000613132Donghuan roadArterial road30912000622131Zhunger roadArterial road23032003632324Kunlun roadArterial road411200364822Ashan roadMinor arterial road3162004According to the artificial neural network assessment model obtained in Section3, the failure states of all units in the road network under different earthquake intensities can be obtained, and the seismic damage index for each road is shown in Table 10.Table 10
Damage state of elements in the network.
Road linksVertexiVertexjNameVIVIIVIIIIX1139Xihuan road0.050.110.320.6321739Kunlun road0.010.030.040.0933839Xihuan road0.020.040.080.1541638Xingfu road0.010.050.120.3053738Xihuan road0.040.090.230.4361537Nanxin road0.110.340.560.8773637Xihuan road0.040.070.210.3681436Shiji road0.030.060.090.1891415Youjian road0.040.070.220.37101314Shiji road0.030.060.090.1811415Nanxin road0.140.360.630.95121516Youjian road0.040.070.240.4213516Xingfu road0.110.330.610.91141617Youjian road0.010.040.090.2715617Kunlun road0.040.070.200.43161718Youjian road0.040.080.270.4617118Zhunger road0.060.130.320.57181819Zhunger road0.120.380.590.9319619Youyi road0.130.390.620.952026Kunlun road0.060.150.330.612156Youyi road0.060.140.320.602235Xingfu road0.020.080.130.292345Yingbin road0.010.030.080.1824411Nanxin road0.010.020.060.1325413Yingbin road0.010.020.060.13261213Shiji road0.010.020.070.15271330Yingbin road0.010.020.060.14283035Yingbin road0.020.040.080.18293536Xihuan road0.020.040.080.18302930Baoshi road0.020.050.130.30311920Zhunger road0.120.380.590.9332720Shengli road0.040.080.220.413378Ashan road0.040.080.220.413427Shengli road0.010.020.040.103529Kunlun road0.060.150.330.613623Shengli road0.040.070.220.3737211Shengli road0.010.030.060.15381011Nanxin road0.060.150.330.61391112Shengli road0.010.020.050.14401229Shengli road0.010.020.050.14412934Shengli road0.010.020.050.14423435Nanhuan road0.020.040.080.18432729Baoshi road0.080.250.480.7744910Dongjiao road0.010.040.090.27451026Nanxin road0.010.020.040.104689Dongjiao road0.010.040.090.27472021Zhunger road0.190.400.710.98482122Changzheng road0.010.020.060.12492224Jinyuan road0.010.020.060.1250923Kunlun road0.010.020.060.12512326Baihua road0.010.040.090.27522526Nanxin road0.010.020.040.10532627Baihua road0.010.040.090.27542733Ruyi road0.010.040.090.27553334Nanhuan road0.020.040.080.18562728Baoshi road0.080.250.480.77573233Nanhuan road0.020.040.080.18582425Jinyuan road0.010.020.060.12592528Jinyuan road0.010.020.060.12602832Jinyuan road0.010.020.060.12613132Donghuan road0.010.020.050.14622131Zhunger road0.010.020.050.14632324Kunlun road0.010.020.060.1264822Ashan road0.010.020.060.12The reliability of each road link under different seismic intensities is obtained based on Section4.1. Figure 6 presents the road network reliability with the intensity of VI, VII, VIII, and IX, respectively. From the figure, it can be observed that the reliability of the road network in intensity VI is very large. With the increase of the seismic intensity, the probabilities of the road network at medium reliable and unreliable become more obvious.(1)
When the seismic intensity is VI, only one road link is in medium reliable. The RNC is almost unaffected, and the traffic order is normal, which can ensure that people can quickly resume normal work and life.(2)
When the seismic intensity is VII, there are 14% of the road links which are in medium reliable, and the rest are in the reliable state. There are no unreliable road links, and the RNC is affected to some extent. The whole transportation function is intact, and relief supplies and ambulance personnel can enter the disaster area at the first time. The normal traffic can be restored after a short period of time.(3)
When the seismic intensity is VIII, there are 59%, 30%, and 11% of the road links which are in the reliable, medium reliable, and unreliable state, respectively. Some road links have lost the RNC, but they can still meet emergency rescue needs.(4)
When the seismic intensity is IX, there are 44%, 33%, and 23% of the road links which are in the reliable, medium reliable, and unreliable state, respectively. Some road links are blocked and difficult to pass, which is mainly due to the safety distance between the buildings on both sides, and the road is not enough. The debris piling after the collapse of the buildings has a great impact on the RNC. Fortunately, most of the buildings and the roads have sufficient safety distance, and the road conditions are good due to the reasonable planning of the Karamay district. Hence, there are about half of the road links which are in the reliable state under intensity IX.Figure 6
Road network reliability diagram under different seismic intensities: (a) VI; (b) VII; (c) VIII; (d) IX. Green lines represent the road links being in the reliable state, yellow represents medium reliable, and red represents unreliable.
(a)(b)(c)(d)We calculated the number of roads connected to each road and the importance index and obtained the RNC before and after the earthquake. Then, the road network reliability under different seismic intensities can be obtained as shown in Table11.Table 11
The reliability value corresponding to the seismic intensity of the road network in the Karamay district.
Seismic intensityVIVIIVIIIIXR0.92990.89410.75430.6631Reliable stateReliableReliableMedium reliableMedium reliable


## 6. Discussion and Future Works

The rescue time after the earthquakes is crucial. The shorter the time to predict the roads, the faster the rescuer can reach the disaster area, which can greatly reduce casualties and economic losses. Therefore, it is necessary to select suitable factors to assess the roads, which can save the time of collecting data and obtaining suitable accuracy. In this study, we first select ten factors that affect road damage and analyzed their importance using the AdaBoost algorithm. Among ten factors of the road damage, the practical intensity, supporting and retaining type, height difference of subgrade, slope height, and soil type of the site are more important than the others. We chose the five factors as the input variables of the ANN assessment model. The accuracy is more than the regression model [8]. However, the model only considers the data of the Wenchuan earthquake and ignores some factors, such as soil liquefaction, secondary disasters, and economic conditions. Whether the foundation soil in the earthquake is liquefied, it is related to the gravel content, the buried depth of gravel soil, the depth of groundwater, and the seismic intensity [34]. The liquefaction phenomenon does not occur until the earthquake occurs, which is a complex phenomenon, and there is no suitable method to evaluate in a few minutes [3]. Most of the current methods are field surveys. Similarly, secondary disasters also have such problems, and it is unreasonable to use only numerical simulations to predict them [18]. In areas with more developed economies, the more developed the construction technology, the higher the level of safety of the corresponding roads. The problems in the evaluation of economic conditions are the lack of information, and some data come from many years ago. Annual inflation affects economic conditions [35], so previous data may not apply to the present. Moreover, the small data sets influence the accuracy and applicability of the ANN model. The damage data of the road networks should be further collected in the future earthquakes.The methodology proposed is more convenient for calculating the reliability of the networks, the speed is faster, and the parameters considered are less than the previous studies [2, 3, 5]. In the evaluation of the road network, we use the median value of the range given by the standard [27] for both the seismic damage index and the pass probability. This part should be optimized in future research to make the predictions more accurate. According to the data, the actual number of earthquake damages is positively correlated with the increase of seismic intensity. The intensity is determined by the macroscopic phenomenon of earthquake damage and the statistical results. The assessment is subjective and ambiguous, and the model using seismic parameters as seismic inputs should be developed, such as the occurrence time, magnitude, and earthquake peak acceleration [10]. Moreover, the method proposed in this study can only evaluate static damage states and road reliability. The dynamic reliability of the road networks ought to be further investigated to reflect the real-time traffic behavior after earthquakes.


## 7. Conclusions

This paper presents a method to assess the reliability of the road networks considering the factors’ importance, damage index, pass possibility, and traffic capacity after the earthquakes. This method is developed based on the machine learning algorithms and linear model. The conclusions of the methodology are as follows: (1) the changes in the practical intensity, supporting and retaining type, height difference of subgrade, slope height, and soil type of the site have a greater impact on the damage of the road. These factors should be considered for the assessment model of the road damage; (2) the seismic intensity has a strong correlation with the damage of the road networks. The greater the seismic intensity, the lower the reliability of the road network and the more serious the damage to the road; and (3) the proposed methodology has a good performance on the evaluation of the road networks. It can be seen that the assessment on the network damages should first consider the individual structure and then calculate the networks.



---

*Source: 9631804-2020-08-06.xml*