Evaluation of Transportation Network Reliability under Emergency Based on Reserve Capacity

There are differences between the requirements for traffic network for traffic demand in daily and emergency situations. In order to evaluate how the network designed for daily needs can meet the surging demand for emergency evacuation, the concept of emergency reliability and corresponding evaluation method is proposed. This paper constructs a bilevel programming model to describe the proposed problem. The upper level problem takes the maximum reserve capacity multiplier as the optimization objective and considers the influence of reversible lane measures taken under emergency conditions. The lower level model adopts the combined traffic distribution/assignment model with capacity limits, to describe evacuees’ path and shelter choice behavior under emergency conditions and take into account the traits of crowded traffic. An iterative optimization method is proposed to solve the upper level model, and the lower level model is transformed into a UE assignment problem with capacity limits over a network of multiple origins and single destination, by adding a dummy node and several dummy links in the network. Then a dynamic penalty function algorithm is used to solve the problem. In the end, numerical studies and results are provided to demonstrate the rationality of the proposed model and feasibility of the proposed solution algorithms.


Introduction
In recent decades, the frequent occurrence of emergencies around the world has caused certain casualties and property losses, posing a great challenge to public security.Especially when these emergencies occur in densely populated urban areas, they are undoubtedly becoming a burden to the already crowded urban traffic.This premise puts forward higher requirements for traffic network under emergency.
Apparently, different from the daily demand of the traffic network, the traffic demand under emergency conditions has an evident feature of unidirection and will surge within a short time.As a result, the existing capacity of road infrastructures is usually not able to meet the demand of evacuation.
To investigate the maximum emergency evacuation demand that can be satisfied by the network designed for daily needs, this paper analyzes the traffic demand under emergency conditions and takes the emergency reliability of traffic network as the research target.
Meanwhile, so as to take full advantage of the finite road resources, ease the traffic pressure which exceeds the load of the road network, increase the evacuation capacity of the road network, and improve the evacuation efficiency, traffic managers will take certain emergency traffic management measures, such as lane reversal.Hence, when studying the traffic network emergency reliability, we should not only consider the particularity of emergency traffic demand, but also consider the impact of these measures on the network service capacity at the same time.As a result, the network under this situation can be deemed as a variable one with capacity constraints, which is more in line with the traits of traffic supply and demand in an emergency.
The remainder of this paper is organized as follows: Section 2 is the literature review.Section 3 is the problem statement.Section 4 establishes the mathematical model and Section 5 proposes the solution algorithm.Section 6 is the case study.Section 7 summarizes the conclusions.

Literature Review
Traffic reliability research began in the 1980s.Over the past decade, transportation reliability has been a new avenue for advancing transportation research in terms of both theories and practical applications, and has attracted tremendous effort over the past few years [1].Chen first introduced capacity reliability as a new network performance index [2].In general, it is defined as one of the important indexes to measure traffic network quality and evaluate network operation state.The higher the reliability of the road network, the more stable the overall operation of the road network, indicating that the structure of the road network is reasonable, with strong anti-interference ability.Existing reliability studies of road networks at home and abroad mainly include the following three categories: connectivity reliability, travel time reliability, and capacity reliability.
Chang [3] analyzed with different paths failure strategy quantificationally by measuring the betweenness centrality of the paths so as to study the reliability of highway and transportation network.Li [4] proposed a method to calculate connectivity reliability of shortest path in the network model of complex system.Based on the concept of network connectivity reliability, Harun [5] discussed the assessment of road maintenance and network capacity expansion planning, revealed the relationship between network capacity expansion and connectivity reliability.
Wael [6] investigated the impact of roadworks undertaken on a given road link over wider parts of the network and assessed network travel time reliability by setting up a network assignment model.Sabyasachee [7] proposed a method based on empirically observed travel time data, to measure and determine the value of travel time reliability, estimate changes in it, and incorporate reliability in the transportation planning process.Vincenza [8] proposed a methodology for estimating travel time reliability of an extended road network, through the calibration of empirical relations and the finding's representation on GIS maps with a dynamic simulation model.Taking into account the influence of travel time variability on the congestion profile, Roxana [9] addressed the valuation of travel time reliability in the presence of endogenous congestion and the role of scheduling preferences.
Researches in capacity reliability mainly concern about the max-flow of the network in early stage, especially in transportation network.Capacity reliability is defined as the largest multiplier applied to a given existing (or basic) OD demand matrix that can be allocated to a network without violating the link capacities [10].Chen [2] extended the capacity reliability analysis by providing a comprehensive methodology, which combines reliability and uncertainty analysis, network equilibrium models, and sensitivity analysis of equilibrium network flow, to assess the performance of a degradable road network.To evaluate capacity reliability, Kuang [11] built up a bilevel programming model based on travel time reliability, concerning OD traffic demand multiplier, and investigated the interaction between two kinds of reliability.Lee [12] proposed a capacity reliability algorithm with Monte Carlo simulation in communication network based on the shortest delay which can indicate the delay's transient variation caused by flow fluctuation.Qin [13] considered quantity of emergency resources, reserve capacity, and number of locations and built an optimization model for emergency resource layout considering service reliability, in which the objective function is to minimize the system cost and various constraints are considered.
The reliability of road network under extreme conditions is very important.In the case of not considering the road damage, how much margin a network has to accommodate emergency traffic demand is an important question that needs to be considered in traffic network design.Based on capacity reliability, this paper evaluates traffic network reliability under emergency conditions from planning perspective.

Problem Statements
This section first analyzes the characteristics of evacuation traffic flow, then summarizes emergency traffic management measures under emergency conditions, and finally leads to the establishment of emergency reliability concept.

Traffic Flow Characteristics under Emergency Conditions.
Different from the traffic state under daily conditions, the traffic state under emergency evacuation has the characteristics of suddenness, high risk, contingency, and so on.In general, emergency evacuation traffic has the following characteristics.
(1) Large and Centralized Demand.From the perspective of space, the origin node of emergency evacuation traffic demand is concentrated in the influenced area.And from the perspective of time, traffic demand for emergency evacuation is basically concentrated in a short period of time before and after the disaster, and people need to be evacuated to a safe area in a short time.This leads to the uneven spatial and temporal distribution of traffic flow directly [14].In the case of emergencies, the peak of travel demand is several times higher than the peak of daily travel demand, which brings great difficulties to the traffic management of emergency evacuation.
(2) Characteristics of Crowded Traffic Flow.Due to the large and centralized traffic demand, a large number of crowded points and congested sections will inevitably appear in the traffic network.The traffic flow on the road is basically saturated or oversaturated, and the evacuation traffic flow is basically in the following state, so it is difficult to find an acceptable gap for overtaking [15].Therefore, emergency evacuation traffic network is a network with congestion influence.
(3) High Reliability Requirement.Sudden disasters threaten people's life and property safety, so emergency evacuees have higher requirements on the reliability of whether they can be quickly evacuated to a safe area.The reliability requirements are as follows: network connectivity reliability, road capacity reliability, travel time reliability, etc.
(4) Uneven Flow of Both Directions.In the process of evacuation, the traffic flow from the affected area to other refuge areas is usually large, resulting in oversaturation and serious traffic congestion.However, the traffic flow into the evacuation zone is small, and most lanes are idle.

Emergency Traffic Management strategies.
Depending on traffic flow characteristics under emergency conditions, there are several emergency traffic management measures that are commonly used in practice.The purpose is to improve evacuation efficiency and reduce casualties and property losses.
(1) Staged Evacuation Strategy.Staged evacuation strategy is widely used in emergency situations when different parts of the target network may suffer different levels of severity over different time windows.By evacuating those populations in the network via an optimized sequence, the staged evacuation strategy can best utilize available roadway capacity, optimally distribute the total demand over the evacuation time horizon, and thus minimize the network congestion level [16,17].
(2) Lane Reversal.In order to alleviate the one-way traffic demand in emergency evacuation which is far beyond the capacity of the road network, the direction of one or more lanes of the road can be reversed to increase the capacity of evacuation direction [15].
(3) Crossing Elimination Strategy.The delay of evacuation vehicles mainly occurs at intersections, which is caused by the conflict and interweaving of evacuation flows from different directions.The interference of evacuation traffic flow can be reduced by eliminating or reducing the confluence point and conflict point of intersections or changing the conflict point into confluence point, which can improve the maximum traffic flow of evacuation network significantly.
The most widely used emergency traffic control strategy is lane reversal.

Establishment of Emergency Reliability
Concept.Given a traffic network which is designed for daily traffic demand, in order to evaluate the extent how the network can meet the evacuation requirements under emergency conditions, the concept of emergency reliability is proposed here.In this study, the concept of network reserve capacity is used to evaluate emergency reliability of traffic network.Network reserve capacity refers to the maximum traffic demand multiplier that can be applied to the given existing (or basic) traffic OD matrix in the process of trip assignment when the road capacity constraints can be satisfied [9].
In this paper, traffic demand multiplier  is used to express emergency reliability of the traffic network.First, the given emergency evacuation demand is expanded  times to obtain a new evacuation demand.Then it is assigned to the traffic network according to evacuees' travel behavior.The impact of lane reversal implementation on road capacity is considered during trip assignment.Since the capacity of each road section has a certain upper limit, there exists a maximum multiplier .
With the constant adjustment of traffic demand multiplier , the trip assignment will approach the network equilibrium state continuously, until the traffic flows allocated to certain road sections reach the upper limit of their capacity.At that time, the maximum value of traffic demand multiplier  is obtained.The maximum value of  reflects the extent of network reliability.If multiplier  ≥ 1, it indicates that the traffic network is reliable.And the larger the multiplier is, the better the traffic network works, and the more reliable the traffic network is.At the same time, it also means that the traffic network can accommodate greater extra traffic demand.
In the process of trip assignment, it is assumed that the evacuation demand from each origin node is known, and the evacuation demand to each destination depends on the choice of evacuees.Each evacuation destination corresponds to a shelter.And the evacuation demand that each shelter can accommodate is limited.It is assumed that, in case of emergency evacuation, evacuees tend to choose the fastest routes and the nearest shelter in order to reach the safe area as soon as possible, considering the congestion on the routes and in the shelters.

Mathematical Formulation
In this section, the mathematical description for the problem of interest is formally proposed with a bilevel programming model to analyze the aforementioned multiplier  quantitatively.

Upper Level Model.
Given an existing traffic network (, ), where  denotes the set of nodes and  denotes the set of existing links.,   ∈ , represent for the two opposite links of the same road segment respectively.  and    denote the number of lanes on the corresponding links.   is the total number of lanes of the segment.
The upper level model aims to maximum the variable  while searching for a feasible lane distribution scheme under unbalanced distribution of two-way traffic flow.The formulation is as follows: In this model, the total number of lanes of a certain segment,    , should be a constant on the basis of the existing traffic network.And the number of lanes of all the segments in each direction in the lane distribution scheme could not exceed the corresponding upper bounds, according to the constraint conditions.
The upper level model influences the lower level model through two sets of variables.The variable  changes the actual evacuating demand calculated in the lower level model, and the number of lanes on each link,   and    , influences the capacity of each link and ultimately affects the capacity of the entire road network.

Lower Level Model.
The lower level model is a traveler behavior model describing route and shelter choices of evacuees in emergency situations.The total number of trips generated at origin nodes is given, and the capacity of a single lane and each shelter is also known.Evacuees will choose the shelter that they can arrive as soon as possible and choose the route by which they can reach the safe area in the shortest time, with the consideration of congestion influence on the routes and in the shelters under emergency circumstance.The number of evacuees reaching each shelter is finally decided by the lower level model.Therefore the lower level model is formulated as a joint UE distribution/assignment model with capacity constraints: ≥ 0 (2f) The definitional constraints are as follows: where In this model, (2a) is the objective function, which is the sum of the integrals of the link performance functions.This function does not have any intuitive economic or behavioral interpretation.It should be viewed strictly as a mathematical construct that is utilized to solve equilibrium problems.Equations (2b) and (2c) represent a set of flow conservation constraints.Equation (2b) states that the flow on all paths connecting each OD pair has to equal the OD trip rate, and (2c) indicates that the total number of trips from each origin node to all the destination nodes has to equal the amplified traffic demand with variable  at each origin node.Equation (2d) indicates that the total number of trips end with a certain shelter is constrained by the capacity of the shelter.Equation (2e) states that the capacity of a link constrains the flow on it.Equations (2f) and (2g) are nonnegative conditions to ensure that the solution of the program will be physically meaningful.
The objective function of program (2a)-(2h), (, , ), is formulated in terms of link flows, traffic demand multiplier and number of lanes on each link, whereas the flow conservation constraints are formulated in terms of path flows.Equation (2h) expresses the incident relationships between link flow and path flow, which is determined by the network structure itself.

Equivalency Conditions of Lower Level Model.
To demonstrate the equivalence conditions of lower level model established above, it has to be shown that any flow pattern that solves the mathematical program (2a)-(2h) also satisfies the equilibrium conditions.This equivalency is demonstrated in this section by proving that the first-order conditions for the minimization program are identical to the equilibrium conditions.These conditions can be derived by forming and analyzing the Lagrangian [15,18,19].
The Lagrangian of the equivalent minimization problem with respect to all the equality and inequality constraints can be formulated as where   ,   ,   , and   are Lagrangian multipliers.
According to Karush-Kuhn-Tucker (KKT) conditions, Lagrangian function (3) must satisfy the following conditions at its extreme points: One can solve the partial derivatives in expression (3) and substitute them into expressions (4a) to (4f), which yields (  −   +   )   = 0 (5c) where Here,   and    both denote the generalized travel time, while   denotes the penalty coefficient related to shelter capacity, and   denotes queuing delay due to limited link capacity respectively.
Equations (5a) and (5b) spell out the familiar userequilibrium condition over the network.With (5k), it can be inferred that the flow on path connecting any origin and destination satisfies That is to say, the paths connecting any O-D pair can be divided into two categories: those carrying flow, on which the generalized travel time equals the minimum O-D generalized travel time; and those not carrying flow, on which the generalized travel time is greater than (or equal to) the minimum O-D generalized travel time.This is satisfied with the first equilibrium principle of Wardrop.
Similarly, with (5c), (5d), and (5l), the following conditions can be speculated: That is to say, when shelter  is chosen by evacuees from origin , the generalized impedance from origin  to shelter  is the smallest one among the generalized impedances from origin  to all the shelters.Otherwise, if none of evacuees choose shelter  as a destination, it means the generalized impedance from origin  to shelter  is greater than (or equal to) the minimum generalized travel time from origin  to all the shelters.This is consistent with the assumption that all evacuees tend to choose the nearest shelter.
From (5e) and (5f), it can be inferred that, for shelter , its capacity penalty coefficient satisfies The penalty coefficient of shelter  depends on the amount of demand reaching this destination.It only exists when shelter  has reached its capacity.Penalty coefficient of shelter  is caused by congestion in the shelter, but the penalty coefficient will not increase indefinitely, as evacuees will change their destinations to avoid excessive delays.
Similarly, from (5g) and (5h), it can be inferred that, for link , its queuing delay satisfies the following: That is, the queuing delay on link  depends on its traffic flow.It only exists when the link traffic volume achieves its capacity.Queue delay is caused by congestion on the link and it causes the traffic flow to be redistributed spatially so that the traffic flow of all links does not exceed their capacity.
By solving the lower programming model, the link flows and path flows (the result of path flows is not unique) that meet the above conditions can be obtained.But if the demand is too large and exceeds the capacity of the whole traffic network, the model will have no solution.
Step 2: Generate initial lane distribution schemes.
Generate  lane distribution schemes randomly, including the original lane distribution scheme.All the schemes are represented by u = { 1 ,  2 , ⋅ ⋅ ⋅ ,   }, which is the change in the number of lanes in each pair of sections.
Step 3: Find the feasible lane distribution schemes.
For each lane distribution scheme: Let  =   − .The algorithm ends.

Solution Algorithm
This section aims to design an algorithm to effectively solve the proposed bilevel model.Since bilevel programming is a NP-hard problem, which is difficult to solve by traditional optimization algorithms or analytical methods, the usual solving methods are based on heuristic algorithms.Therefore, in this paper, iterative optimization algorithm is adopted in the upper model.Besides, the addition of constraint conditions in the lower model changes the original Descartes form of UE model, which results in the fact that the lower model cannot be solved by traditional F-W method.Commonly used methods include penalty function method and Lagrange multiplier method.In this paper, the lower model is transformed and then solved by dynamic penalty function algorithm proposed by Zhang [19].
Step-by-step procedures of two algorithms are also given below, respectively.

Solution Algorithm for Upper Level Model.
The specific steps of the iterative optimization algorithm are shown in Algorithm 1.
The process of local searching algorithm in Algorithm 1 is as follows.Lane schemes are generally generated in a monotonous way, such as random generating method, which can easily get trapped in local optimality.In order to obtain better newly generated schemes, this paper starts from historically feasible lane distribution schemes in set P and searches their neighborhoods to get new schemes based on the characteristics of specific problems.
For each road section, judge the flow distribution results of a certain historical feasible scheme; if   /  ≥    /   ,   increases by 1 with probability  while ensuring that it does not exceed its upper limit; otherwise,   decreases by 1 with probability  while ensuring it does not exceed its lower limit.This is the principle of generating a new lane distribution scheme.

Solution Algorithm for Lower Level Model.
As for the lower level model, the topology of the origin network is firstly transferred.Augment the original network with a dummy node, denoted   , and connect all the destination nodes (i.e., shelters) to node   , deriving several dummy links.The total number of dummy links is the product of the number of shelters.Assume that free flow travel time on each dummy link is 0, and capacity of each dummy link equals the capacity of the corresponding shelter.
This simple modification of the network topology can transform the combined distribution/assignment problem with constraint of traffic capacity into an equivalent UE problem with constraint of traffic capacity.Based on the modification, a dynamic penalty function algorithm is used to solve the problem.
The basic idea of the dynamic penalty function algorithm is to transform equivalent UE problem with constraint into traditional UE problem by adding penalty function.In UE problem with constraint, the generalized travel time    at equilibrium is composed of travel time calculated by impedance function and queuing delay denote by Lagrangian multipliers   .However, the Lagrangian multipliers cannot be calculated directly, thus a penalty function (  ,   ) =   (  /  ) is used to replace it.Here,   is penalty function coefficient related to link a.Let  denotes the ratio of   to Step 1: Initialization.
During the problem solving process, both   and (  /  ) will be constantly updated according to the new link traffic flow obtained from each iteration, so that the penalty function (  ,   ) will finally approach Lagrange multiplier   .Update rule for   is as follows: Due to the poor convergence performance of the traditional Frank-Wolfe method for solving UE model, the linear decomposition algorithm proposed by Aashitiani and Magnanti [20,21] is used here to solve UE problem.Based on the path flow, the linear decomposition algorithm starts from the nonlinear complementarity problem (NCP) of UE problem and transforms it into a series of linear complementarity problems.
The specific steps of the dynamic penalty function algorithm for the lower level model are shown in Algorithm 2.

Numerical Experiments
In this section, the above proposed model and algorithm will be tested and evaluated in an experimental network.First the topology of network and relevant information is given, then the analyzation of the results.

Topology of Network.
Take the test network in the Figure 1 as the example, which has 9 nodes and 24 directed links.Here node 1 and node 4 are the origin points of evacuation, and node 3, 6, and 9 are shelters.Information of link properties in the network, including link free flow travel time and link capacities, is listed in Table 1.Information of evacuation demand of each origin point is listed in Table 2, and Table 3 gives the information of the capacity of shelters.Modify the original network by adding a dummy node 10 and connect it with three shelters by three dummy links.The topology of modified network is shown in Figure 2. As mentioned above, free flow travel time of each dummy link is always 0, and capacity of each dummy link equals to the capacity of the corresponding shelter.On the basis of the modification above, the problem of interest is transformed into UE assignment problem of a multiorigin and singledestination network with link capacity constraints.

Result Analysis
(1) Result Analysis of the Lower Model.The trip assignment results are analyzed when  = 1 without changing the number of lanes in each link.In the experiment, values of each parameter are as follows:  = 0.01,  1 =  2 = 0.001, the initial penalty function coefficient  0  = 0.1 0 , where  0 denotes the average of free flow travel time on all links.The BPR formula is employed as impedance function and two parameters in the function are set as 0.15 and 4, respectively.
(a) From the Perspective of Route Choice.Tables 4 and 5 show the path details of OD pair 1-10 and 4-10 at equilibrium.From the perspective of route choice, for OD pair 1-10, three paths are chosen by evacuees, and the generalized travel time for these three paths is all equal to 4.8443, which is the minimum generalized travel time of this OD pair.The generalized travel time of unselected paths including path 4, 5, and 6 is 6.1449, 6.0435, and 6.9712, respectively, all larger than 4.8443.Similarly, for OD pair 4-10, three paths are chosen by evacuees, and the generalized travel time for these three paths is all equal to 3.6940.The generalized travel time of unselected paths including path 4, 5, and 6 is 4.5875, 5.6443, and 6.0504, respectively.
Although the path flow results of traffic assignment in crowded traffic networks with queuing delays are not necessarily unique, the result satisfies that generalized travel time of all the selected paths is equal, which meets the conditions of UE equilibrium.
(b) From the Perspective of Shelter Choice.According to the path flow results, the demand from origins to all the shelter can be referred, shown in Table 6.For origin node 1, the number of evacuees choosing shelter 3, 6, and 9 is 900, 42.3956, and 557.6044, respectively.And the generalized travel time from node 1 to all the shelters is equal to 4.8443.For origin node 4, the number of evacuees choosing shelter 6 and 9 is 757.6044 and 442.3956, respectively.No one chooses shelter 3 because its generalized travel time is 3.8375, a little larger than other shelters.
That means for a certain origin, generalized travel time of all the selected shelters is equal, and generalized travel time of shelters which are not chosen is larger than the chosen shelters.
Table 7 shows the final shelter demands.The number of evacuees choosing the three shelters is 900, 800, and 1000, respectively.And among them, only shelter 6 reaches its capacity, and the queuing delay is 1.1064 .The other two shelters are not saturated and the queuing delays are both 0. This is consistent with the previous conclusion that the queuing delay associated with shelters only exists when the shelter has reached its capacity.All the experimental results are consistent with the above analysis of first-order conditions of the Lagrangian function.
(2) Results Analysis of Bilevel Model.In order to illustrate the rationality of the model and the effectiveness of the algorithm proposed in this paper, three different demand scenarios and two different shelter capacity scenarios (shown in Table 8) are constructed, and four cases of combination of different scenarios are generated (shown in Table 9).
(a) Results Analysis of Different Cases.In case 1, the maximum multiplier  that the solution algorithm can obtain is 1.11.The corresponding feasible lane distribution scheme and link flow are shown in Table 10.Only links 5, 15, and 19 are saturated, but there is still plenty of potential room for the capacity of these links.Table 11 shows the corresponding shelter demands.Three shelters are chosen by 997, 800, and 1200 evacuees, respectively.All the shelters are almost saturated.It turns out that the emergency reliability of the entire road network in the case here is limited to the shelter capacity.In case 2, by increasing the shelter capacity, the maximum multiplier  coms up to 2.00.The corresponding results at optimal are shown in Tables 12 and 13.Although shelter 6 is  In conclusion, the emergency reliability of the entire road network can be limited by both shelter capacity and road capacity.Although improving road capacity is a good way to improve network emergency reliability, sometimes, increasing the shelter capacity, there would also be a significant increase in the network emergency reliability.
(b) Performance Evaluation of the Proposed Algorithm.Case 3 and 4 are obtained by adjusting the pattern of evacuation demand.Both cases can be solved effectively by the proposed algorithm.
By analyzing the solution results of all the above cases and combining the enumeration method, the true optimal value of multiplier  in different cases can be obtained.The true optimal can also be checked by the network topology under feasible lane distribution schemes.The results are shown in Table 14.It can be seen that, in different cases, the optimal multiplier calculated by the proposed algorithm is almost equal to its true optimal value, and the slight difference in some cases is caused by the step size.Therefore, the algorithm proposed in this paper is proved to be of high accuracy.
The computation time of the proposed algorithm is mainly consumed in solving the lower level joint UE distribution/assignment problem.Therefore, in this paper, the times of solving the lower assignment problem are used to evaluate the solving efficiency of the proposed algorithm.Each case runs 10 times, and the average times to solve the lower assignment problem are listed in Table 14.It can be seen that the calculation times are acceptable in those numerical experiments.In addition, the research results of this paper are applicable to evaluate the existing road network at the level of management planning, and the solving efficiency of the model can meet the actual needs.

Conclusions
Reliability is an important feature of transportation network.In order to evaluate the reliability of road network under emergency evacuation, a bilevel programming model is proposed in this paper.The concept of reserve capacity is adopted, and the influence of reversible lane measures taken under emergency conditions, as well as the capacity limits of links and shelters, are considered.An iterative optimization method is proposed to solve the upper level model, and the lower level model is transformed and then solved by a dynamic penalty function algorithm.Finally, the numerical example demonstrates the rationality of the proposed model and feasibility of the proposed solution algorithms.
The proposed lower level model can effectively reflect the route choice and shelter choice behavior of evacuees under emergency condition.Evacuees may choose to evacuate to another shelter if the shelter with smaller travel time is full.Besides, when the shortest path has reached its maximum capacity, evacuees are forced to choose other subprime paths.The assignment result accords with UE equilibrium.
Although the reserve capacity of the network can be improved by reversible lane measures, but both link capacity and shelter capacity could be the restriction of emergency reliability of the entire road network.In addition, the efficiency of the proposed algorithms is also verified in the case study.
The model proposed here is a little complicated and many practical problems have been simplified.In further study, the proposed model can be not difficult to contain congestion and delays at intersections.Additionally, uncertainty factors should be considered, including demand uncertainty, disaster location uncertainty, link capacity uncertainty, and so on.
3.1) Calculate the number of lanes on each link ,   , and its capacity   ; (3.2) Run the lower programming model according to demand     ; (3.3)If any solution exists, store the corresponding lane distribution scheme into the set of historical feasible schemes P. If all schemes can't find a feasible solution in the lower level model, they should be renewed according to following procedure: (4.1) Renew the update count of lane schemes, , and if  > , turn to Step 6; (4.2) If there are more than  schemes in set P, randomly pick  schemes and do mutation processing by using local search algorithm, then obtain  new schemes.(4.3)If there are less than  schemes in set P, suppose it is , do mutation processing by using local searching algorithm to all of these schemes.Besides, generate  −  schemes at random.(4.4) Turn to Step 3. Step 5: Let  +1 =   + ,  =  + 1, and turn to Step 3.
2) Set valid paths of each O-D pair    = ⌀.Step 2: Calculate initial solution (2.1) Calculate the generalized travel time   according to the initial link flow {  = 0, ∀}, and find the shortest path between each OD pair; (2.2) Apply all-or-nothing assignment to get the initial path flow .Step 3: For each OD pair : (3.1) Update the set of valid paths    and calculate the generalized travel time   , obtain the shortest path   and its travel time   between each OD pair .If the convergence criterion Eq. ( * ) is not satisfied, then    =    ∪ {  }; Update the link flow   = () and penalty function parameter   .

Table 1 :
Information of link properties.

Table 2 :
Information of evacuation demand.

Table 3 :
Information of shelter capacity.

Table 4 :
Path details of OD pair 1-10 at equilibrium.

Table 5 :
Path details of OD pair 4-10 at equilibrium.

Table 6 :
Detailed information for OD pairs.

Table 7 :
Detail information for shelter choice.

Table 8 :
Information of different demand and shelter capacity scenarios.

Table 9 :
Information of combination of different scenarios.

Table 10 :
Feasible lane distribution scheme and link flow in Case 1 when  = 1.11.

Table 11 :
Detailed information for shelter choice in Case 1 when  = 1.11.

Table 12 :
Feasible lane distribution scheme and link flow in Case 2 when  = 2.00.

Table 13 :
Detailed information for shelter choice in Case 2 when  = 2.00.
Meanwhile, the opposite links of link 1, 3, and 5 have no remaining capacity, which means the capacity of these three links cannot be enlarged by applying reversal lane strategy.In other words, the emergency reliability of the entire road network is limited.

Table 14 :
Comparison between calculated optimal solution and true optimal solution.