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.

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

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):

(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).

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

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.

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

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.

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

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.

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.

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.

Acknowledgments

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