Optimization and Planning of Emergency Evacuation Routes Considering Traffic Control

Emergencies, especially major ones, happen fast, randomly, as well as unpredictably, and generally will bring great harm to people's life and the economy. Therefore, governments and lots of professionals devote themselves to taking effective measures and providing optimal evacuation plans. This paper establishes two different emergency evacuation models on the basis of the maximum flow model (MFM) and the minimum-cost maximum flow model (MC-MFM), and proposes corresponding algorithms for the evacuation from one source node to one designated destination (one-to-one evacuation). Ulteriorly, we extend our evaluation model from one source node to many designated destinations (one-to-many evacuation). At last, we make case analysis of evacuation optimization and planning in Beijing, and obtain the desired evacuation routes and effective traffic control measures from the perspective of sufficiency and practicability. Both analytical and numerical results support that our models are feasible and practical.


Introduction
Major emergencies, concerning accidents such as natural disasters, public health events, and social abrupt affairs, occur in a very short period of time and always cause serious damage to the society. They are all characterized by suddenness, uncertainty, and serious harmfulness, and have been a big challenge for sustainable development of the human society.
In recent years, various disasters occur frequently. Therefore, countries all over the world attach great importance to emergency management. China has entered into the highincidence season of emergencies and will face the baptism brought by emergencies in a long period of time [1]. With the further acceleration of urbanization, the population is more concentrated, so how to respond effectively and timely to emergencies is especially important. In most emergency management, how to evacuate people to safety zone becomes a crucial step. For example, cyclohexane spill that occurred in 1976 in Seveso, Italy, rained a cloud of dioxin on surrounding communities, and the local authority organized a evacuation of 220,000 people. Similarly, when liquid chlorine cylinder explosion accident occurred in southern China in 1979, emergency evacuation was also implemented to move 60,000 people away. Likewise, in 1984, after Mexico City gas storage tank exploded, 350,000 people were evacuated [2].
The evacuation issues can be divided into two different types: small scale evacuation and long-distance regional evacuation [3]. Small scale evacuation generally refers to a kind of evacuation when emergencies affect small area swiftly and violently, such as explosion within a finite range, house collapsing, fire breaking out in shopping malls, and so on. Generally, evacuation of this kind mainly deals with evacuation on foot, rather than using vehicles. In contrast, longdistance regional evacuation refers to a form of evacuation implemented along with the appearance of a wide range of emergencies, such as leakage and diffusion of harmful gases and earthquakes. Since the evacuation route is long-distance, this kind of evacuation generally needs vehicles for transport. In addition, according to the difference of subject and actions taken after emergencies happen, the long-distance regional evacuation can be divided into autonomous evacuation, recommended evacuation, and mandatory evacuation [4]. Autonomous evacuation and recommended evacuation are 2 The Scientific World Journal carried out spontaneously and preparedly by people when they get messages or notifications of disasters, in which case the alert time is comparatively longer, while mandatory evacuation responds to emergencies having the need of urgent evacuation and generally commanded by governments or related departments. All of the evacuations considered in this paper belong to the latter situation.
The rest of our paper is organized as follows. Section 2 presents a review of related literature. Section 3 deals with one-to-one evacuation optimization problem based on MFM. Section 4 develops one-to-one evacuation model based on MC-MFM. Section 5 explores the one-to-many evacuation problem. Case analysis is reported in Section 6, and finally some concluding remarks are given in Section 7.

Literature Review
Generally, scholars mainly tend to consider long-distance regional evacuation or mandatory evacuation first, figuring evacuation routes as network diagrams to achieve shortest routes, maximum flow, or minimum-cost maximum flow, such as adopting contraflow transportation network reconfiguration [5], establishing the cell-based evaluation network model [6], analyzing multifactors with the construction of a system which involves the evacuation network [7], building the network flow model for lane-based evacuation routing [8], and solving evacuation problems by using dynamic network flows [9,10]. Lots of researchers have not only optimized evacuation route, but also explored many methods. The representatives refer to Choi et al. [11], Han [12], Hoppe and Tardos [13], Klüpfel et al. [14], Theodoulou [15], and Adewmi and Garba [16]. With regard to multiobjective evacuation routing, Stepanov and Smith [17] present a methodology for designing optimal routing policies for emergency evacuation planning (EEP) through building an integer program (IP) model. On the problem of roadway deployment, Pal and Bose [18] propose an integer model to find the best location and assign response vehicles to those depots under reliability constraints. Saeed and Ram [19] examine the evacuation routing problems by proposing a two-step approach that consists of an incapacitated integer multi-commodity network model and a computational algorithm. Considering nonlinear relations, Zhang et al. [20] explore a min-costpursued swapping dynamic (NMSD) system to model the evolution of selfish routing games on the traffic network.
Besides, many topics related to emergency evacuation also have been discussed. Regarding the evacuation problem in surrounding areas of nuclear power plant, Dunn and Newton [21] investigate how to evacuate people as far as possible within the prescribed time, looking upon the evacuation problem as a maximum flow problem to find the optimal evacuation route. As to the evacuation in an earthquake, Yamada [22] describes route optimization of neighborhoods as the shortest route problem and the minimum cost flow problem, respectively. Campos and da Silva [23] try to reduce the conflicts between finding the evacuation flows and increasing the traffic capacity of evacuation route. They treat the ratio of traffic capacity and travel time as a measurable index of traffic performance. On the issue of hurricane evacuation, Dixit and Radwan [24] provide a new and innovative technique with a "network breathing strategy" at destinations after dictating when to schedule evacuation orders and capacities required on different routes. Based on the data of earthquakes in 1994, Northridge and 1995, Kobe, Koike et al. [25] provide a probabilistic approach to simulate the evacuation scenario along the streets crowded with evacuation people. Specifically, the ratio of fire-proof structures, ignition rates from fire gutted houses, and population density are all assumed in a probabilistic manner. In view of blocking effects on crowd movement, Luh et al. [26] establish a new macroscopic network-flow model assuming that fire, smoke, and psychological factors can evoke crowd's desire to escape at the expected flow rate and then develop a divide-andconquer approach to reduce computational complexity and reflect psychological changes.
As for effective solutions for route optimization, Hamacher and Tjandra [27] all systematically summarize the dynamic network flow models, which are widely used in the optimization of human organization evacuation planning, including the dynamic maximum flow model, the fastest dynamic flow model, and the global maximum dynamic flow model. Besides, they analyze the utilization potentiality of these models in regional and extensive evacuation. Mamada et al. (2003) [28] treat the optimization of evacuation routes and the assignment of departure times as the fastest flow problem. With the assumption that the traffic starting from the source node could only depart for a designated destination, they provide a algorithm to find the fastest flow in the tree network diagram. Similarly, Lu et al. [29] regard the evacuation planning that includes departure time, evacuation routes, and ending point selection as the fastest transshipment problem, and propose a heuristic algorithm on the basis of network figural representation to find the optimal organization planning directly. Comprehensively, Lin [30] constructs a framework integrated with optimizations, evaluations, and adjustments, aiming to achieve the optimization of evacuation organization planning. Fuellerer et al. [31] consider that, on two-dimensional level, successive routing of the vehicles in the road network by a colony algorithm can satisfy customers' demands with the consideration of the freight loading factor. Likewise, Leung et al. [32] develop an extended guided tabu search (EGTS) method and a new heuristic packing algorithm for the two-dimensional loading vehicle routing problem, which can help tabu search to escape from local optimum effectively. From the perspective of information quality and evolution, Pillac et al. [33] present a general description of dynamic routing and introduce the notion of degree of dynamism, and then bring out a comprehensive review of applications and solution methods for dynamic vehicle routing problems. When capacity constraints exist, Xu et al. [34] propose a stochastic user equilibrium assignment model for a schedule-based transit network, through which they can simultaneously predict how passengers choose their transit vehicles to minimize their travel costs and estimate the associated costs. And they also find that when a connection segment reaches its capacity level, the Lagrange multipliers of the mathematical problem The Scientific World Journal 3 give the equilibrium passenger overload delays in this transit network.
In addition, some scholars formulate models to find the concrete flow and path from source nodes to designated destinations. To minimize the end-to-end delays for the specific routing mechanism, Grimmell and Rao [35] discuss the problem that deals with the transmission of a message from a source node to a designated destination over a network with propagation delays and dynamic bandwidth constraints on the links. They make available bandwidth for each link specified as a piecewise constant function and present for each message forwarding. In order to identify the network paths accurately, Murtaza et al. [36] present a new mechanism for detecting shared bottlenecks between end-to-end paths in a network, which is based on the wellknown linear algebraic approach-singular value decomposition (SVD). Specific to the evacuation in the feeder-bus network, Deng et al. [37] extend the demand pattern M-to-1 (i.e., multiple bus stops and a single station) to M-to-M, considering the passenger travel cost. Moreover, they present a new genetic algorithm to determine the optimal feederbus operating frequencies under strict constraints, and finally find that demand distributions also should be considered when designing a feeder-bus network.
The literature mentioned above mostly converts evacuation problems to dynamic or static network model by finding shortest routes, maximum flow, and minimum-cost maximum flow. However, intuitive calculation makes the entire emergency management on evacuation too rigid, and there are many actual-world factors that cannot be fully reflected in the model.
First of all, intersections are everywhere in reality, so the impacts of them should be considered during evacuation process in avoidance of possible time delays. Although it may make computing process more complex, it is better than ignoring intersections' impact. By the way, the complexity of intersections does not only reflect on the limitation of certain straight road sections, but also on that of different turning routes. Being different from classic emergency evacuation management, this paper fully examines the influence of intersections.
Secondly, previous researches on impacts of traffic management and control measures are very few. The actual evacuation may be affected by many aspects of constraints, which are all objective conditions that constrain the evacuation carried on. Therefore, control measures should be considered in the model, such as duplicated row, single row, and forbidding for passing as well as some uncertain real-time management and control measures, like signal control and real-time directing.
Finally, regarding the problem of establishing evacuation models, two different kinds of evacuation, namely, evacuations from one source node to one designated destination(one-to-one evacuation) and one source node to multiple designated destinations (one-to-many evacuation), also are taken into account according to the actual situation.

Evacuation Optimization Model Based on MFM
When dangerous situations are unknown to people, the maximum number of evacuations can only be implemented within the capacity limitations, and each batch of evacuation must be the maximum flow. In addition, in order to facilitate the description and calculation, the maximum flow below generally refers to the flow.

Assumptions and Prerequisites.
In the network diagram model with the given time and the goal of maximum flow, its basic assumptions and prerequisites are as follows.
(1) There is only one source node and one designated destination, and vehicle number in the source node are known.
(2) Allowed evacuation time can be obtained by forecasts, and thus evacuation time is given.
(3) Traffic capacity of each road section and each intersection's different turnings are known.
(4) Limitations of traffic capacity are only valid for one batch, and next batch has new limitations on its traffic capacity.
(5) The goal is to maximize evacuation vehicle number under limitations of evacuation time and traffic capacity, namely, to find the maximum flow in a network diagram.
In the above prerequisites, assumption (4) also reflects a characteristic of this modeling. Maximum-flow model generally refers to the model for achieving one-time maximum flow under traffic capacity constraints. Due to the limitations of traffic capacity, maximum flow in multiple times is taken into account in this model. Under this circumstance, batches will be limited, and each batch has the same limitation of traffic capacity. Moreover, the maximum flow includes vehicle number of the th batch and the ( − 1)th batch, where = ⌈ / ⌉.

Notations.
Nodes include source nodes, intersections, and evacuation destinations. Arcs represent road sections between the intersections, source nodes, destinations, and adjacent intersections, which convert the entire road network to directive network diagrams = ( , , , ), where represents a vertex in the figure, expressed by characters , , . The node 1 and node represent the source node and the designated destination, respectively.
represents the arc, while ( , ) represents straight line arc from to , and ( , , ) represents steering arc from ( , ) to ( , ) through . stands for the node weight set, and represents the traffic capacity of from road section ( , ) to road section ( , ) through .
is the arc weight set, and stands for the traffic capacity of road section ( , ). represents the total number of evacuation flow.
represents the flow on arc ( , ).
represents the flow from the arc ( , ) to arc ( , ) through .

Formulations Based on MFM.
Based on the above assumptions and definitions, this evacuation can be expressed as the following models: While formula (1) is the objective function, formula (2) is its decomposition type; formulas (3) and (4) are constraints of traffic capacity of road sections and intersections, respectively; formula (5) means that the flow of each batch reaching the destination equals the flow departing from the source node; formulas (6) and (7) are the flow conservative constraints of each intersection; formula (8) represents that the flow of the final batch is not more than the maximum evacuation flow of front batches; formula (9) is nonnegative constraints.
3.4. Solutions. The model established above is not linear, so linear programming cannot be used directly to find the solution. Certainly, this model can be converted into one or multiple linear programming models, but obviously the solution process will be more complex. Comparatively speaking, it is much more convenient to establish the network diagram model.
From the view of formula 1, the key point of the solution is to obtain , 1 , and .
In fact, 1 can be obtained by the MFM in the network diagram, using labeling method. Traditional solution to find maximum flow in network diagram is to label points in sequence to obtain the augmented chain. However, the solution of this model is to label arcs. Besides, arcs can be labeled including straight line arcs and steering arcs at the intersection on road sections. In this model, let represents the set of the augmented routes; let + and − represent the set of forward arcs and backward arcs in the augmented route, respectively. Let [ ] or [ ] represent the prior label of arc or arc on augmented chain. Steps of the calculation are as follows.
Step 1. Give an initial feasible flow in the network. Zero flow can be treated as the initial feasible flow.
Step 2. Label the arc and find an augmented chain.
(1) Label a random arc that regards the source node as the starting point with {∞, ℎ }, while ℎ is the remaining possible maximum flow of the arc. (3) If a labeled arc's endpoint is the destination, it represents the augmented chain has been found. When continuous labeling cannot make the labeled arc point to the destination, it indicates that the augmented chain does not exist, so the calculation ends.
Step 3. Adjust the flow.
(1) Find the minimum value of arc ℎ or arc ℎ on the augmented chain, and denote it with ℎ.
(2) Adjust the flow according to the following formulas: Then obtain the new possible flow 1 .
Step 4. Cancel all labels and repeat Steps 2 and 3 constantly until new augmented chain cannot be found. At this moment, sets consisting of of each arc are the maximum flow sets.
The Scientific World Journal 5 In fact, when maximum flows of each batch are implemented, we can know the time spent in each evacuation, as well as the evacuation time of each possible route. In the final evacuation, it is obvious that not all routes were able to evacuate in place, which means that evacuation routes consuming shorter time will be chosen preferentially. For example, before th batch evacuation begins, it remains about time after prior -1 batches, and then the last evacuation batch only aim at the route with the evacuation time less than , and is the flow of final batch evacuation.

Evacuation Model Based on MC-MFM
Our evacuation model investigates a comparatively less clamant emergency situation, such as slight toxic gas leak.
Everyone is required to be evacuated out of the danger zone. Here, people can rationally implement evacuation route planning and establish models with the goal of minimizing total evacuation time and maximizing the evacuation vehicle number. So evacuation problems in practical can be abstracted as a kind of minimum-cost maximum flow model (MC-MFM). Evacuations in multiple batches are still considered in this model, and only subsequent evacuation routes are implicit in the diagram. Specific model is as follows.

Premise and Hypothesis.
In the network model with the object of obtaining the maximum flow within the given evacuation time, its basic premise and hypothesis are as follows.
(1) Only one source node and one designated destination exist, and the number of evacuation vehicles are known.
(2) Evacuation intervals of each batch are given and known.
(3) Traffic capacity of each road section and each intersection are known.
(4) The average traffic time on each road section and the average delay of different turnings of each intersection are known.
(5) The goal is to minimize the total evacuation time and maximize evacuation vehicle number.

Symbol Definition.
Nodes include source node, intersections, and evacuation destination. Similarly, arcs include road sections between the intersections, evacuation source node, and destination. Therefore, the network can be converted to network diagram = ( , , , , , ).
is the vertex, expressed by characters , , .
is the arc, while ( , ) is the arc from to .
is a point weight set, while represents the delay of road section ( , ) to the road section ( , ) through in th batch.
is also a point weight set, while represents the traffic capacity of road section ( , ) through to the road section ( , ).
represents the arc weight set, while represents average travel time of road section ( , ) in th batch.
represents the arc weight set, while represents the traffic capacity of road section ( , ). 0 is the total number of vehicles at source node.
stands for the flow of arc ( , ) in th batch.
represents the flow of arc ( , ) to arc ( , ) through in th batch.
is the evacuation interval in each batch, also known as the valid time of traffic.
Besides, node 1 and node stand for the source node and the designated destination, respectively.
While formula (12) is the objective function; formulas (13) and (14) are the traffic capacity constraints of road sections and intersections, respectively, which means that the flow of road sections and intersections does not exceed the traffic capacity; formula (15) is the constraints of aggregate demand of evacuation, which illustrates that the total flow to the destination is equal to the total flow departing from source node; formula (16) refers to the flow of evacuation road section ( , ) in th batch evacuation is equal to the total flow turning to other road sections at node ; formula (17) means that the flow of road sections ( , ) in the th batch evacuation is equal to the total flow turning from other road sections to this road section through node , and formula (18) is nonnegative constraints.

Solutions.
Likewise, this model can also be converted to linear programming and solved by simplex method, but obviously it is easier to be solved by our MC-MFM. 6 The Scientific World Journal If we find the first minimum cost flow with the cost of , then we assume that there are numerous other cost flows existing with costs of + , + 2 , and so on. Each time a feasible flow with minimum cost is found, we assume that there are numerous virtual routes. And hereby we declare that traffic capacities of these virtual routes are recovered and have the same value as the first batch route. In the process of looking for smaller cost flows, if a smaller cost is not less than + , then we take flow in the second batch as the minimum cost possible flow, and so on. Let represents the augmented route or the minimum cost route in the process of solution. + and − are the set of forward arcs and backward arcs in the augmented route.
[ ] is the prior node of node in the augmented route or the minimum cost route. The steps of calculation are as follows.
Step 1. Give an initial feasible flow 0 , which can be zero.
Step 2. Find a minimum cost road from node 1 to node in this network diagram.
Step 3. Augment the unit flow along route and add flow ℎ in augmented route + or intersection turnings. Otherwise, reduce flow ℎ in augmented route − or intersection turnings. Then a new possible flow 1 is formed.
Step 4. After one round of augmentation, wipe off the arcs which have been saturated and then repeat Steps 2 and 3 and continue to augment until 0 has been fully evacuated out. Then calculate the total time.
Specific calculation methods of this model will be displayed in the case analysis.

Extension to One-to-Many Evacuation
The models mentioned above are the evacuation from one source node to one designated destination. To make the study go further and deeper, the evacuation from single source node to many designated destinations will be considered in this part.
As for the evacuation from single source node to many designated destinations, the destination is no longer a single parameter as but expressed by many parameters such as 1 , 2 , . . . , . We only need to change formula (15) to ∑ 1 = ∑ ∑ , which represents the total flow arriving at multiple destinations is equal to the total flow of the source node.
We have two ways as follows to solve this problem. The first one is to set a virtual evacuation destination and then assume the delay time or travel time of all turning arcs and line arcs toward to arc ( , ) are zero, and the traffic capacity is ∞. In this way, evacuation from one source node to many designated destinations can be turned back to one source node to one designated destination.
The second one is to solve the model of one-to-many evacuation directly. In the process of calculation, the domain of the minimum cost route becomes broader. For example, if destinations 1 , 2 , and 3 exist, then we must find minimum cost routes from the source node to 1 , 2 , and 3 , respectively and make comparisons to find the minimum cost route which will be regarded as the augmented route at this time. Though this method has differences with the first one, the essence is the same.
Certainly, emergency evacuation problem in one-tomany evacuation has a little difference from one-to-one evacuation. In other words, multiple evacuation destinations may have multiple capacity limitations, which require us to consider whether the evacuation flow toward certain destination will reach the maximum capacity.
According to Figure 1, we abstract its main routes to form a network route map, see Figure 2.

Case Analysis
As we all know, Beijing is a big city with high concentration of population, especially in the downtown. Once some hazardous emergencies occur, it is likely to cause very serious losses. Under such circumstances, this paper will consider the specific route map within the second ring of Beijing, and the second ring road will be abstracted into network diagrams. Since the routes of the second ring road are comparatively complex, the abstracted route network diagrams are simplified according to importance, and more actual values for each parameter are given to ensure the accuracy of case analysis.

Selection of Network Diagrams.
This case analysis will calculate separately according to two different kinds of evacuation model, which is on the basis of actual second ring road map of Beijing, as shown in Figure 1.
Assume that the black spot in Figure 2 is the source node. It is located between the imperial palace and the Beihai Park, which is also the area with the most intensive stream of people. Small white circles represent the assumed possible destinations. We assume that the big circle is the range of evacuation, and it is required that people should be evacuated to a destination or multiple destinations out of this circle. In order to calculate more conveniently, this paper ignores some minor roads and only considers main roads, so this road map will be simplified to Figure 3.
According to the situation of actual evacuation, some road sections and intersection turnings should be regulated, and some possible traffic control measures should be implemented. Possible routes of each road section and accessible intersection turning will be described, and the traffic capacity and the amount of time consumed, including the travel time at road sections and the delay time at intersection turning, will be illustrated, see Tables 1 and 2. Possible road sections and intersection turnings described in the above two tables are designed based on the traffic control measures. For example, there is no turning (3, 2, 11) at intersection 2 and turning (8,9,13) at intersection 9. In addition, we also assume that destinations will not be taken as the midpoint of other evacuation routes. For instance, there are no similar turnings as (19,3,17) and (15,12,13), and so forth. Of course, limitations of these control measures make calculations more convenient but do not affect the accuracy of the entire case. Figure 3 is marked with multiple destinations, hereby we consider one-to-one evacuation at first and select node 13 as the destination. In order to calculate conveniently, we still assume that evacuation routes cannot pass other alternative destinations beside destination node 13. Thus, evacuation route network diagram can be simplified in Figure 4.

Problem Description.
As the intersections 1, 6, 11, and 12 are no longer T-roads and intersections 2, 7, and 9 are also no longer crossroads, the delay time of each intersection turning will change inevitably. The traffic capacity and the amount of time consumed, including travel time and delay time of each road section and intersection turnings, are showed in Table 3.
In this paper, the total evacuation time is limited to 60 mins and the allowed evacuation time in each batch is 16 mins.

Solutions.
The first thing is to obtain the value of and = ⌈ / ⌉ = ⌈60/16⌉ = 4. The evacuation time in the fourth batch is only 12 mins.
The second thing is to get the static maximum flow 1 .
Step 1. Let initial flow 0 be zero.
Label the arc with the next three rounds with adjustment.
Step 3. Since the limit time in the fourth batch is 12 mins. Therefore some routes, in which the travel time and the delay time are significantly more than 12 mins, will not be considered. Only route (0-1-9-13) is possible, and its flow is 8.
Due to the evacuation only towards a fixed destination, the delay time of some intersections is shorter than one-tomany evacuations. In our case, the traffic control measures have been considered. For example, road section (9, 10) is a two-way traffic, and most road sections are only one-way traffic.

Emergency Evacuation Based on MC-MFM.
Then we continue to consider one-to-one evacuation with node 13 as the destination. For the simplicity of calculation, we still assume that the evacuation route cannot pass other alternate evacuation destinations except node 13. Thus, network diagram of evacuation routes is also illustrated in Figure 4. In addition, the total number of vehicles in the source node 0 is 60 hundred vehicles. We need to find the minimum cost route and the augment route on this way until evacuation flow reaches the target value. The approach used here is almost like the Dijkstra algorithm, but the review for the routes with minimum cost has been changed into the review for arcs.
The first thing is to obtain the minimum cost route, see Table 5.
It can be seen from Table 5 that the minimum cost of route 0-1-9-13 is 10. As the bottleneck traffic flow on this route is arc (1,9,13), which is 8, this route is the augmented chain. When we increase its traffic flow about 8, arc (1,9,13) of the first batch is saturated at this moment, and the total route cost is 80, namely 8000 mins. In addition, subsequent batches of this route should be considered in the course of finding the minimum cost route.
As illustrated in Table 5, since the first batch route 0-2-11-10-12-13 does not include arc (1,9,13), the total cost of this route is 16, which is more than the virtual route's minimum value in the preamble, and thus the minimum cost route is route 0-1-9-13 in the second batch. As the bottleneck of traffic flow is arc (1,9,13), which is 8, the traffic capacity on this augmented chain is increased with 8, and arc (1,9,13) in the  second batch is saturated at this moment with the total route cost of 112, namely 11200 mins. We repeat the above steps according to Table 5 and obtain the following results in Table 6.

Case Analysis of One-to-Many Evacuation.
Simple case calculations are implemented in the above section for the one-to-one evacuation problem on the basis of MC-MFM. In this section, the form of evacuation is extended to oneto-many evacuation. According to the second ring road map of Beijing, the evacuation destinations include eight nodes, which are nodes 13, 14, 15, 16, 17, 18, 19, and 20. In specific evacuation, the problems of route optimization are ensured based on MFM.   To keep the comparability with the case in the previous calculation, the total number of vehicles in the source node is still 60. The numerical values of specific capacity and the amount of time consumed are shown in Table 1. Source node is node 0, and evacuation destinations are 13, 14, 15, 16, 17, 18, 19, and 20. In spite of being a different form, we still use label censorship method for calculation. Specific calculation process is shown in Table 7.
As illustrated in Table 7, we can find the minimum cost route 0-6-18 with the cost of 8. As the bottleneck of traffic capacity on arc (0, 6, 18) is 8, this route is augmented chain and traffic flow should be increased with 8. After that, arc (0, 6,18) in the first batch is saturated, and the total route cost is 64. In addition, in the process of finding the minimum cost routes in the next step, subsequent batches should be considered.
Similarly, route 0-1-4-16 in the first batch does not include the arc (0, 6, 18), and its cost is 10, which is the minimum cost in the remaining viable routes, so this route is augmented chain. Since the bottleneck of traffic flow is arc (4,16), which is 10, traffic flow on this augmented chain can be increased with 10. At this point, arc (4,16) in the second batch is saturated, and the total route cost is 100.
We repeat the above steps according to Table 7 and obtain the following results, see Table 8.
Compared with one-to-one evacuation, one-to-many evaluation will take just 13 minutes, which is much faster and more efficient. Therefore, we draw conclusions as follows: under normal circumstances, when a certain emergency occurs and the emergency evacuation is required, multiple destinations for the evacuation are quicker and more efficient than only one destination. Therefore, in order to prevent heavy casualties caused by major emergencies and improve evacuation efficiency, more evacuation roads and shelters should be built in the places with high population density and high occurrence of emergencies.

Conclusions
Appropriate traffic control measures must be taken into account in evacuation route optimization. Overall, this paper explores evacuation optimization and planning of evacuation routes, considering some traffic control measures. We abstract the road network as directive network, taking the vehicles waiting for evacuation in the source node as flows and evacuation road sections along the road as arcs which are under traffic control, such as one-way, two-way, delay, and so forth. Besides, the traffic capacity of intersections is expressed as the weights of network diagram. After model formulation, we use graph theory method to solve our model and then verify its feasibility.
Based on the case analysis to the second ring of Beijing, we present evacuation optimization and planning after emergencies happen. Although we make a lot of prerequisites and assumptions for modeling and calculation, there are still some issues remained to be explored as follows.