Advances in Civil Engineering

Advances in Civil Engineering / 2020 / Article
Special Issue

Advances in Transportation Geotechnics

View this Special Issue

Research Article | Open Access

Volume 2020 |Article ID 9631804 |

Jinlong Liu, Hanxi Jia, Junqi Lin, Heng Hu, "Seismic Damage Rapid Assessment of Road Networks considering Individual Road Damage State and Reliability of Road Networks in Emergency Conditions", Advances in Civil Engineering, vol. 2020, Article ID 9631804, 16 pages, 2020.

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

Academic Editor: Xue Zhang
Received23 Feb 2020
Revised28 Jun 2020
Accepted21 Jul 2020
Published06 Aug 2020


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.

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 [911]. 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 [1618]. 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 (Table 1). 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].

NumberHeight difference (m)Road gradesSupporting and retaining typeSeismic fortification intensityPractical intensity (degree)Slope height (m)Soil type of the site

G213: K1020 + 96010.43Gravity-type groutedVIIXI13.6II
G213: K1018 + 6001.23Gravity-type groutedVIIXI40.0I
G213: K1022 + 90015.63Reinforced retaining wallVIIXI0I
G213: K1029 + 8503.13Gravity-type groutedVIIX33.0I
G213: K1008 + 9003.03Gravity-type groutedVIIXI20.0I
G213: K1014 + 1756.03Gravity-type groutedVIIXI0I
G213: K1008 + 4002.03Cutting slopeVIIXI50.0III
G213: K1008 + 5809.93Gravity-type concreteVIIXI50.0III
G213: K1012 + 40013.63Gravity-type concreteVIIXI15.5III
G213: K1029 + 7006.23Gravity-type concreteVIIX50.0III
G213: K1023 + 7002.23Prestressed cableVIIXI19.0III
G213: K1009 + 0807.13Cutting slopeVIIXI0III
G213: K1008 + 9804.03Gravity-type concreteVIIXI0III
G213: K916 + 5083.63Cutting slopeVIIIX7.0III
G213: K37 + 3501.82Gravity-type groutedVIIX294.0II
G213: K38 + 5001.02Gravity-type groutedVIIX326.3II
G213: K40 + 1003.62Gravity-type groutedVIIX121.2I
G213: K42 + 4004.42Facing wallVIIX175.5II
G213: K26 + 8001.52Cutting slopeVIIXI24.3I
G213: K73 + 0002.02Cutting slopeVIIIX175.9I
G213: K29 + 9501.02Facing wallVIIXI19.7I
G213: K34 + 1401.52Cutting slopeVIIXI122.7III
G213: K35 + 6003.92Cutting slopeVIIXI100.0III
G213: K50 + 2003.02Cutting slopeVIIX480.1III
G213: K58 + 1603.42Gravity-type groutedVIIX15.0III
G2131.53Gravity-type groutedVIIIX7.1III
G2131.73Gravity-type groutedVIIIX92.1II
G2134.03Gravity-type groutedVIIIX27.0II
G2133.43Cutting slopeVIIIX93.4III
G2132.33Facing wallVIIIX71.0III
G2131.53Gravity-type groutedVIIIX0II
G21313Gravity-type groutedVIIIX10.0II
G2131.23Cutting slopeVIIIX0II
S303: K10 + 3411.82Gravity-type groutedVIIX5.0III
S303: K28 + 9410.52Gravity-type groutedVIIX50.0II
S303: K40 + 5000.52Gravity-type groutedVIIIX27.4III
S303: K16 + 4412.62Cutting slopeVIIX237.5II
S303: K31 + 8411.02Cutting slopeVIIIX28.7I
S303: K36 + 6001.02Cutting slopeVIIIX100.0III
S303: K43 + 2001.02Cutting slopeVIIIX30.0II
S303: K26 + 9003.12Gravity-type groutedVIIIX25.0III
S3022.02Cutting slopeVIIIX27.0III
S302: K732 + 5601.52Gravity-type groutedVIIX7.4III
S3022.52Facing wallVIIX47.5I
S1052.23Gravity-type groutedVIIXI2.8I
S1053.13Gravity-type groutedVIIXI25.3III
S1053.53Facing wallVIIXI36.7I
S105: K167 + 4605.93Cutting slopeVIIXI8.4III
S105: K168 + 9007.43Cutting slopeVIIXI9.0III
XU092.34Gravity-type concreteVIX44.3II
S2101.33Gravity-type groutedVIIVII324.4I
S2103.73Cutting slopeVIIVII63.6II
S2102.33Cutting slopeVIIVII24.4III
G2124.23Gravity-type groutedVIIVIII34.3II
G2131.02Gravity-type groutedVIIVIII20.3I
S3033.84Cutting slopeVIIVII194.3III
S3034.24Gravity-type groutedVIIVIII35.0III
G3171.53Cutting slopeVIIVII173.3II
G3172.52Gravity-type groutedVIIVI2.5II
G3172.42Gravity-type groutedVIIVI5.0I
G3176.02Gravity-type groutedVIIVI0III

There 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. Figure 1 demonstrates some typical damages of the road in the Wenchuan earthquake, which were divided into some road segments and had the individual number [24].

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 (Figure 2) (Tables 2 and 3):

Mean shear wave velocity (m/s)Site category

ClassificationSlope protection type

Level INo protection, wall protection, plant protection, spraying (no net)
Level IINet hanging, shotcrete, active and passive protection net
Level 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 Table 3; 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 (Table 4).


Highway classificationExpresswayBetter seismic performance
First-classBetter seismic performance
Third-classFailure occurs more frequently
Fourth-classFailure occurs more frequently
Subgrade typeNo digging and no fillingNo damage
CuttingDestruction accounts for 26% of the total damage
EmbankmentDestruction accounts for 28% of the total damage
Half-filled and half-dugFailure occurs most frequently (48%)
Supporting and retaining typeReinforced retaining wall, prestressed cable, facing wallBeneficial to mitigating damage
Gravity-type concreteFailure occurs most frequently
Gravity-type wet masonryMore serious damage
NoNot conducive to mitigating damage
Soil type of the siteI
IIIIntensifying the damage slightly
IVIntensifying the damage
Slope height (h)0 mThe least number of damages occurred
0 m < h ≤ 15Less destruction occurs
15 m < h ≤ 40 m
40 m < h ≤ 120 mMore susceptible to damage
120 m < hMore susceptible to damage
Slope angle (θ)The least number of damages occurred
0° < θ ≤ 35°
35° < θ ≤ 65°Failure occurs most frequently
65° < θMore susceptible to damage
Protective measureNo, steening retaining wallNot conducive to resisting seismic, small quantity
Frame beam, precast blockBeneficial to seismic
Hanging net and guniting, safety netting systemCan resist part of seismic
Shotcrete, mortarNot conducive to resisting seismic
Height difference of subgradeh ≤ 1 m
1 m < h ≤ 3 m
3 m < h ≤ 5 m
5 m < h
Seismic fortification intensityVI
Practical intensityVI

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

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 Section 3.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: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:where y is the value of the cost function, n presents the amount of the roads, is the damage state of the ANN model, and t is the actual damage state.where is the ReLU function, and x is the input of the ReLU function.

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 ( 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 Section 3.3 gave only the damage state and could not directly obtain the pass probability, , 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 .

Damage stateDamage index (I)Pass statePass probability ()


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, , can be computed by [32]where is the average traffic interval between consecutive traffic. When the city does not have the value of , the road capacity can use the value in Table 6.

Calculated driving speed (km/h)50403020

NP (pcu/h)1690164015501380

According 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, , can be computed bywhere is the road classification coefficient of the vehicle lane (Table 7).

Road classificationExpresswayArterial roadMinor arterial roadBranch road


The initial RNC, , before the earthquake is defined as follows:where 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 is usually calculated based on the betweenness centrality. However, the betweenness centrality involves calculation in a complex network [33]. Hence, the index can be expressed as (7) for the convenience of calculation and speed of the assessment:where is the number of road links connected to the i-th road link and is the total number of road links. If the road link contains bridges, tunnels, etc., the pass probability is obtained based on the ANN model, and the RNC after the earthquake can be obtained by (8). The reliability of the road network, , is defined as the ratio of to (9). Table 8 presents the reliability state of the road network, including reliable, medium reliable, and unreliable.

Reliability stateReliableMedium reliableUnreliable

Reliability of the road network

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 Figure 5. The road network contains 64 road links and 39 nodes, totaling 60 km (Table 9).

Road linksVertex iVertex jNameRoad classificationLength (m)Year built

1139Xihuan roadArterial road11122000
21739Kunlun roadArterial road6822003
33839Xihuan roadArterial road6972000
41638Xingfu roadArterial road7012002
53738Xihuan roadArterial road16252000
61537Nanxin roadArterial road9032002
73637Xihuan roadArterial road8642000
81436Shiji roadArterial road9422003
91415Youjian roadMinor arterial road8462003
101314Shiji roadArterial road6852003
11415Nanxin roadArterial road6752002
121516Youjian roadMinor arterial road11902003
13516Xingfu roadArterial road6512002
141617Youjian roadMinor arterial road6702003
15617Kunlun roadArterial road6812003
161718Youjian roadMinor arterial road11082003
17118Zhunger roadArterial road6802003
181819Zhunger roadArterial road7002003
19619Youyi roadArterial road11062003
2026Kunlun roadArterial road9142003
2156Youyi roadArterial road6402003
2235Xingfu roadArterial road9642002
2345Yingbin roadArterial road8842000
24411Nanxin roadArterial road8022002
25413Yingbin roadArterial road8332000
261213Shiji roadArterial road9292003
271330Yingbin roadArterial road7172000
283035Yingbin roadArterial road9552000
293536Xihuan roadArterial road23222000
302930Baoshi roadMinor arterial road11242012
311920Zhunger roadArterial road9162003
32720Shengli roadArterial road5922003
3378Ashan roadMinor arterial road7022004
3427Shengli roadArterial road5122003
3529Kunlun roadArterial road8022003
3623Shengli roadArterial road5152003
37211Shengli roadArterial road8432003
381011Nanxin roadArterial road7662002
391112Shengli roadArterial road6962003
401229Shengli roadArterial road7362003
412934Shengli roadArterial road13092003
423435Nanhuan roadArterial road16332012
432729Baoshi roadMinor arterial road8892012
44910Dongjiao roadMinor arterial road11712004
451026Nanxin roadArterial road7072002
4689Dongjiao roadMinor arterial road5122004
472021Zhunger roadArterial road10132003
482122Changzheng roadArterial road5992004
492224Jinyuan roadArterial road7542000
50923Kunlun roadArterial road3702003
512326Baihua roadMinor arterial road13382004
522526Nanxin roadArterial road8222002
532627Baihua roadMinor arterial road8202004
542733Ruyi roadMinor arterial road11602004
553334Nanhuan roadArterial road7242012
562728Baoshi roadMinor arterial road12442012
573233Nanhuan roadArterial road15032012
582425Jinyuan roadArterial road10092000
592528Jinyuan roadArterial road10092000
602832Jinyuan roadArterial road6082000
613132Donghuan roadArterial road30912000
622131Zhunger roadArterial road23032003
632324Kunlun roadArterial road4112003
64822Ashan roadMinor arterial road3162004

According to the artificial neural network assessment model obtained in Section 3, 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.

Road linksVertex iVertex jNameVIVIIVIIIIX

1139Xihuan road0.050.110.320.63
21739Kunlun road0.
33839Xihuan road0.
41638Xingfu road0.
53738Xihuan road0.
61537Nanxin road0.110.340.560.87
73637Xihuan road0.
81436Shiji road0.
91415Youjian road0.
101314Shiji road0.
11415Nanxin road0.140.360.630.95
121516Youjian road0.
13516Xingfu road0.110.330.610.91
141617Youjian road0.
15617Kunlun road0.
161718Youjian road0.
17118Zhunger road0.060.130.320.57
181819Zhunger road0.120.380.590.93
19619Youyi road0.130.390.620.95
2026Kunlun road0.060.150.330.61
2156Youyi road0.060.140.320.60
2235Xingfu road0.
2345Yingbin road0.
24411Nanxin road0.
25413Yingbin road0.
261213Shiji road0.
271330Yingbin road0.
283035Yingbin road0.
293536Xihuan road0.
302930Baoshi road0.
311920Zhunger road0.120.380.590.93
32720Shengli road0.
3378Ashan road0.
3427Shengli road0.
3529Kunlun road0.060.150.330.61
3623Shengli road0.
37211Shengli road0.
381011Nanxin road0.060.150.330.61
391112Shengli road0.
401229Shengli road0.
412934Shengli road0.
423435Nanhuan road0.
432729Baoshi road0.080.250.480.77
44910Dongjiao road0.
451026Nanxin road0.
4689Dongjiao road0.
472021Zhunger road0.190.400.710.98
482122Changzheng road0.
492224Jinyuan road0.
50923Kunlun road0.
512326Baihua road0.
522526Nanxin road0.
532627Baihua road0.
542733Ruyi road0.
553334Nanhuan road0.
562728Baoshi road0.080.250.480.77
573233Nanhuan road0.
582425Jinyuan road0.
592528Jinyuan road0.
602832Jinyuan road0.
613132Donghuan road0.
622131Zhunger road0.
632324Kunlun road0.
64822Ashan road0.

The reliability of each road link under different seismic intensities is obtained based on Section 4.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.

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 Table 11.

Seismic intensityVIVIIVIIIIX

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

Data Availability

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

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

The work in this manuscript was done with the participation of all the authors. Dr. Jinlong Liu was the main planner of the manuscript. The artificial neural network assessment model for an individual road and the analysis model for the reliability of road networks were mainly built by Dr. Liu. However, Ph.D. student Hanxi Jia has also done important work. She was mainly responsible for collecting data of road damage in the Wenchuan earthquake and accomplishing the assessment of the factor importance using the AdaBoost Algorithm. Professor Lin gave them a lot of advice on road damage assessment, and he reviewed the case study results in detail. Postgraduate Heng Hu provided them with detailed data of the road network in Section 5 and helped them to complete the road network reliability analysis.


The authors are very grateful to Professor Junqi Lin for his help. This work was supported by the National Key Technology R&D Program of China (grant no. 2018YFC1504401) and the Heilongjiang Provincial Key Laboratory of Underground Engineering Technology Open Project of China (grant no. 2017-HXYKF-06).


  1. C. A. Davis, A. Mostafavi, and H. Wang, “Establishing characteristics to operationalize resilience for lifeline systems,” Natural Hazards Review, vol. 19, no. 4, Article ID 04018014, 2018. View at: Publisher Site | Google Scholar
  2. P. Bocchini and D. M. Frangopol, “A stochastic computational framework for the joint transportation network fragility analysis and traffic flow distribution under extreme events,” Probabilistic Engineering Mechanics, vol. 26, no. 2, pp. 182–193, 2011. View at: Publisher Site | Google Scholar
  3. A. Kiremidjian, J. Moore, Y. Y. Fan, O. Yazlali, N. Basoz, and M. Williams, “Seismic risk assessment of transportation network systems,” Journal of Earthquake Engineering, vol. 11, no. 3, pp. 371–382, 2007. View at: Publisher Site | Google Scholar
  4. A. Guo, Z. Liu, S. Li, and H. Li, “Seismic performance assessment of highway bridge networks considering post-disaster traffic demand of a transportation system in emergency conditions,” Structure and Infrastructure Engineering, vol. 13, no. 12, pp. 1523–1537, 2017. View at: Publisher Site | Google Scholar
  5. L. Cheng, M. Du, X. Jiang, and H. Rakha, “Modeling and estimating the capacity of urban transportation network with rapid transit,” Transport, vol. 29, no. 2, pp. 165–174, 2014. View at: Publisher Site | Google Scholar
  6. Y. Ji, R. G. Mishalani, and M. R. McCord, “Estimating transit route OD flow matrices from APC data on multiple bus trips using the IPF method with an iteratively improved base: method and empirical evaluation,” Journal of Transportation Engineering, vol. 140, no. 5, Article ID 04014008, 2014. View at: Publisher Site | Google Scholar
  7. X. Zhou, M. Venigalla, and S. Zhu, “Bounding box approach to network pruning for efficient path search through large networks,” Journal of Computing in Civil Engineering, vol. 31, no. 5, Article ID 04017033, 2017. View at: Publisher Site | Google Scholar
  8. H. Heng, “Research on assessment method of road seismic damage based on the Wenchuan earthquake,” Institute of Engineering Mechanics, China Earthquake Administration, Beijing, China, 2018. View at: Google Scholar
  9. Q. Yu and R. Tian, “Research on reversal lane application method of urban road network based on the Bi-level programming,” Advances in Intelligent Systems and Computing, vol. 279, Springer, Berlin, Germany, 2014. View at: Publisher Site | Google Scholar
  10. K. Pitilakis, S. Argyroudis, K. Kakderi, and J. Selva, “Systemic vulnerability and risk assessment of transportation systems under natural hazards towards more resilient and robust infrastructures,” in Proceedings of the 6th Transportation Research Arena, vol. 14, pp. 1335–1344, Warsaw, Poland, April 2016. View at: Publisher Site | Google Scholar
  11. P. Gauthier, A. Furno, and N.-E. El Faouzi, “Road network resilience: how to identify critical links subject to day-to-day disruptions,” Transportation Research Record: Journal of the Transportation Research Board, vol. 2672, no. 1, pp. 54–65, 2018. View at: Publisher Site | Google Scholar
  12. H. Hu, J. Lin, and J. Liu, “Improved methodology for seismic estimation of highway based on Wenchuan earthquake,” in Proceedings of the 2018 International Conference on Modeling, Simulation and Optimization, pp. 291–296, Xiamen, China, November 2018. View at: Publisher Site | Google Scholar
  13. Y. Lecun, Y. Bengio, and G. Hinton, “Deep learning,” Nature, vol. 521, no. 7553, pp. 436–444, 2015. View at: Publisher Site | Google Scholar
  14. P.-j. Chun, H. Yamashita, and S. Furukawa, “Bridge damage severity quantification using multipoint acceleration measurement and artificial neural networks,” Shock and Vibration, vol. 2015, Article ID 789384, 11 pages, 2015. View at: Publisher Site | Google Scholar
  15. K. Morfidis and K. Kostinakis, “Approaches to the rapid seismic damage prediction of r/c buildings using artificial neural networks,” Engineering Structures, vol. 165, pp. 120–141, 2018. View at: Publisher Site | Google Scholar
  16. N. Wang, W. Cheng, M. Zhao, Q. Liu, and J. Wang, “Identification of the debris flow process types within catchments of Beijing mountainous area,” Water, vol. 11, no. 4, p. 638, 2019. View at: Publisher Site | Google Scholar
  17. S. Park and J. Kim, “Landslide susceptibility mapping based on random forest and boosted regression tree models, and a comparison of their performance,” Applied Sciences, vol. 9, no. 5, p. 942, 2019. View at: Publisher Site | Google Scholar
  18. H. Jia, J. Lin, and J. Liu, “An earthquake fatalities assessment method based on feature importance with deep learning and random forest models,” Sustainability, vol. 11, no. 10, p. 2727, 2019. View at: Publisher Site | Google Scholar
  19. B. H. Menze, B. M. Kelm, R. Masuch et al., “A comparison of random forest and its Gini importance with standard chemometric methods for the feature selection and classification of spectral data,” BMC Bioinformatics, vol. 10, pp. 1–16, 2009. View at: Publisher Site | Google Scholar
  20. D. P. Solomatine and D. L. Shrestha, “Adaboost.RT: a boosting algorithm for regression problems,” in Proceedings 2004 IEEE International Joint Conference, Budapest, Hungary, July 2004. View at: Google Scholar
  21. Code of China, “Chinese specification of seismic design for highway engineering (JTG B02-2013),” Code of China, Beijing, China, 2013. View at: Google Scholar
  22. Code of China, “Chinese seismic intensity scale (GBT 17742-2008),” Code of China, Beijing, China, 2008. View at: Google Scholar
  23. L. Chen, Report on the Highways’ Damage in the Wenchuan Earthquake, China Communications Press, Beijing, China, 2012.
  24. Code of China, “Code for design of urban road engineering (CJJ 37-2012),” Code of China, Beijing, China, 2012. View at: Google Scholar
  25. Code of China, “Chinese code for design of highway subgrades (JTG D30-2015),” Code of China, Beijing, China, 2015. View at: Google Scholar
  26. D.-s. Ling, Y. Zhao, B. Huang, F. Zhang, and Y. Zhou, “Analysis of dynamic stress path in inhomogenous subgrade under moving aircraft load,” Soil Dynamics and Earthquake Engineering, vol. 111, pp. 65–76, 2018. View at: Publisher Site | Google Scholar
  27. I. Chowdhury, R. Tarafdar, A. Ghosh, and S. P. Dasgupta, “Seismic response of rectangular liquid retaining structures resting on ground considering coupled soil-structure interaction,” Bulletin of Earthquake Engineering, vol. 15, no. 9, pp. 3695–3726, 2017. View at: Publisher Site | Google Scholar
  28. A. Imtiaz, V. Perron, F. Hollender et al., “Wavefield characteristics and spatial incoherency: a comparative study from argostoli rock-and soil- site dense seismic arrays,” Bulletin of the Seismological Society of America, vol. 108, no. 5A, pp. 2839–2853, 2018. View at: Publisher Site | Google Scholar
  29. D. L. Shrestha and D. P. Solomatine, “Experiments with AdaBoost.RT, an improved boosting scheme for regression,” Neural Computation, vol. 18, no. 7, pp. 1678–1710, 2006. View at: Publisher Site | Google Scholar
  30. China Earthquake Administration, “Classification of earthquake damage to lifeline engineering GB/T 24336-2009,” China Standard Press, Beijing, China, 2009. View at: Google Scholar
  31. Standardization Administration and State Administration for Market Regulation, “The Chinese seismic intensity scale GB/T 17742-2008,” China Standard Press, Beijing, China, 2008. View at: Google Scholar
  32. Ministry of Housing and Urban-Rural Development of the People’s Repulic of China, “Code for design of urban road engineering CJJ 37-2012,” China Architecture and Building Press, Beijing, China, 2016. View at: Google Scholar
  33. Z. Safra and U. Segal, “Risk aversion in the small and in the large: calibration results for betweenness functionals,” Journal of Risk and Uncertainty, vol. 38, no. 1, pp. 27–37, 2009. View at: Publisher Site | Google Scholar
  34. M. S. Asadi, R. P. Orense, M. B. Asadi, and M. J. Pender, “Post-liquefaction behavior of natural pumiceous sands,” Soil Dynamics and Earthquake Engineering, vol. 118, pp. 65–74, 2019. View at: Publisher Site | Google Scholar
  35. S. Hatfield-Dodds, H. Schandl, P. D. Adams et al., “Australia is “free to choose” economic growth and falling environmental pressures,” Nature, vol. 527, no. 7576, pp. 49–53, 2015. View at: Publisher Site | Google Scholar

Copyright © 2020 Jinlong Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.