Dynamic Route Network Planning Problem for Emergency Evacuation in Restricted-Space Scenarios

We investigate a dynamic route planning problem in restricted-space evacuation, namely, theMultiobjectiveDynamic RouteNetwork Planning (MODRNP) problem. It models the multisource to multidestination evacuation in restricted-space scenarios, with the objectives of minimizing the whole evacuation delay and maximizing the evacuation efficiency. We study the problem in 3D scenarios, which can provide intuition vision for the geographic space and contribute to the evacuation plan and implementation. Based on the auxiliary graph transformation, we propose a heuristic algorithm referred to the classical problem, Minimum Weighted Set Cover. We finally conduct extensive experiments to evaluate the performance of the proposed algorithm and give an application instance on a typical kind of restricted-space scenarios. The results indicate that the proposed algorithm outperforms the existing alternatives in terms of the utilization as well as timeliness.


Introduction
As a significant research problem, disaster prevention and relief became one of the new requirements for social development in recent years, and emergency evacuation planning is an important part of rescue and relief work.Emergency evacuation generally contains the evacuee evacuation [1] and the relief resource deployment [2], and the evacuee evacuation is the most urgent work of disaster relief.An effective evacuation path planning strategy can not only effectively reduce the evacuation time consumption, but also prevent unnecessary congestion in the evacuation process, even secondary damage.Thus, the optimization design and theoretical analysis for evacuation path planning are of great practical significance [3].
According to the disaster scenarios, evacuation path planning can be classified into free-space (e.g., urban streets and squares) planning and restricted-space (e.g., underground pipeline, tunnel, and mine laneway) planning [4,5].The former one has most applications in city-security and antiterrorist, which has been well studied based on one-source to one-destination model or one-source to multidestination model.
The later one involves in a large range of industrial applications, i.e., the evacuation planning in fire disaster, gas explosion, or water inrush in underground mine.The restricted-space has several complex characteristics which are different from the free space: (a) the space and time consuming limitations on entrance and evacuation.The density of personnel is sparser; (b) the unsuitability of long operating time.The reliability of temporary shelters is weaker; (c) the trapped-vulnerability of internal structure in disasters.When a disaster occurs, the number of available evacuation paths is much less.Thus, the restricted-space has higher accident rate and brings more difficulties and challenges to rescue and relief work, which has been urgent to be studied.To improve the quality of emergency response in restrictedspace planning, we focus on dynamic path planning based on the multisource to multidestination model.
In practical, restricted-space emergency evacuations have different evacuation requirements according to the accident type, nature, and distribution.On the one hand, considering 2 Journal of Advanced Transportation the limitations and unsuitability of restricted spaces, it may cause local or global congestions when evacuation traffic load is too large, e.g., the evacuees intensively rushing to a certain exit.Thus, evacuation requirements do not only include the shortest evacuation time consumption, but also the unobstructed connectivity of route network.The goals of the restricted-space emergency evacuation are minimizing the global evacuation time consumption and maximizing the utilization rate of each exit (balancing among the exits' traffic loads).To minimize the global evacuation time consumption, each evacuee tends to choose the shortest path and nearest exit, which may cause the imbalance of exit utilization, i.e., most evacuees intensively rushing to a certain exit via their shortest paths.To maximize the utilization rate of each exit, the evacuees should dispersedly choose the exits, which may cause that some evacuee chooses the suboptimal path rather than the shortest one and its evacuation time consumption will be lengthened.Thus, the two objectives are conflicting and they need a tradeoff to balance their realization, i.e., each evacuee should choose the global optimum rather than local optimum.On the other hand, due to the trapped-vulnerability of restricted spaces, the static planning strategies cannot adapt to the dynamic spreading of disasters or accidents and other environment changes.Thus, it is necessary to design dynamic evacuation planning based on environment changes.Based on the above considerations, we propose new evacuation routes planning problem, the Multiobjective Dynamic Route Network Planning (MODRNP) problem, and aim to design dynamic route network planning strategies and give theoretical analysis.
In this paper, firstly we introduce the restricted-space model and formulate the Multiobjective Dynamic Route Network Planning (MODRNP) problem in restricted-space evacuation.Secondly, we propose dynamic and globaloptimized strategies to construct the evacuation route network, whose performance will be evaluated by the simulations and instance experiments at last.The rest of the paper is organized as follows.We summarize the related work in Section 2. Section 3 introduces the restricted-space model and problem definitions.Section 4 introduces the strategy framework for MODRNP Problem and describes the algorithms.Simulation results and corresponding discussions are given in Section 5. Section 6 concludes this paper.

Related Works
As a crucial kind of combinational optimization problem, the path planning problem has been widely studied such like recent research works on traveling salesman problem (TSP) [6], traffic assignment/vehicle routing problem [7][8][9], and evolutionary optimization problems, e.g., general transportation planning problems, facility location problem, and roadway repair problem [10].Based on the existing theoretical achievements, the research on evacuation path planning can be classified as static planning and dynamic planning.
Research on the static path planning problem provides theoretical guideline for the emergency analysis and escaping exercise preparation.Most algorithms for it are based on Dijkstra algorithm, improved Dijkstra algorithm, Floyd algorithm, the first  shortest path algorithm,  * algorithm, and dynamic programming.Under the static problem model, the solutions generally simplified the disaster type, evacuee psychological pattern, and other factors by introducing the path equivalent length, with the assumption that the disaster would be spreading steadily.
The studies on the dynamic path planning problem are the basis of real-time decision making in emergency response and management in disaster/accident.The solutions can be classified into two kinds, the approximate algorithms with high theoretical precision and the heuristic ones with high efficiency.
To design approximate algorithms for the dynamic problem model, most researchers firstly transformed the problem into the classical network flow problem model [11][12][13][14] and then utilized the classical polynomial time algorithms to design the path selection strategy.By applying the maximum flow algorithm, Dunn and Newton [15] constructed a path network with the goal of guiding the most survivors in accordance with the constructed path network under the condition of permission capacity.By applying the minimum cost flow algorithm, Yamada [16] assigned evacuation traffic aimed at minimizing the overall length of the chosen evacuation paths, and Cova and Johnson [17] focused on complex path network with the multiple objectives of minimizing the intersections between the evacuation paths and minimizing the overall evacuation length.By applying the fast flow algorithm, the authors in [18,19] aimed at the application scenarios modeled as directed graph.But the above algorithms could be only applied to the case of grid structure path network and unified path capacity.By applying the contraflow algorithm, the authors in [20] considered the setup time of the contraflow operation and proposed a two-layer algorithm to minimize the evacuation time consumption and the setup time for flow patterns.Since the algorithms for the flow problems have high time or space complexity, most of these strategies lack high efficiency in practice.
To design heuristic algorithms for the dynamic problem model, researchers in [21,22] designed simulated evolutionary algorithms.They took full advantages of the positive feedback information mechanism of ant colony algorithm and the fast convergence of the genetic algorithm [23].For the airport emergency evacuation scenarios, the authors in [24] performed a series of simulations based on an agentbased model (ABM) to determine the collective behaviors and the overall evacuation time, which are affected by the environment and complex structures.And [25] was motivated by the improvement of the marine evacuation system plan in a practical case, Dalian offshore airport of China, and developed a modeling framework to evaluate and optimize the emergency evacuation capability.For the crowded building emergency evacuation scenarios, the authors in [26] studied the human behavior of crowd evacuation to guide the design of buildings and the reinforcement of security management.The proposed algorithms can solve the path planning problem in densely populated scenarios.Due to the sparsity of personnel in restricted spaces, these strategies cannot be applied directly in the evacuation planning.
The above related work had limitations at three aspects.Firstly, most research works considered application scenarios for two-dimensional (2D) planes, in which the proposed algorithms cannot adapt to three-dimensional (3D) application scenarios.Secondly, most existing solutions for the path network construction problem were with high time or space complexity, which may bring low efficiency in restrict space evacuation.Thirdly, the existing researches on dynamic planning of emergency evacuation paths were still limited; e.g., the authors in [27] conducted research on the vehicle routing problem for postdisaster scenarios.And their application effect in practical is not ideal in terms of dynamic adjustment and timely feedback.Therefore, we intend to introduce a new dynamic planning problem model and utilized new technology in solution for restricted-space scenarios.

Problem Formulation
3.1.Restricted-Space Model.We focus on the emergency evacuation in confined space works and model such a densely populated restricted-space, as a 3D connected graph  = (, ), where  = {V 1 , ⋅ ⋅ ⋅ , V  } is the set of predetermined observation vertices in the space,  is the set of bidirectional edges, and for each V  ∈ , V  = (  ,   ,   ,   ) and for each   ∈ ,   =<   ,   > (1 ≤  <  ≤ ).For any vertex subset   ⊆ , [  ] is the subgraph of  induced by   .Similarly, [  ] is the subgraph of  induced by an edge subset   ⊆ .The source set  and the destination set  are important subsets of the vertex set in , which are composed of the evacuees' initial positions and the feasible exits, respectively.Here we assume  and  are predefined.
When a disaster or an accident happens, the evacuees need to escape via the edge/segment in the restricted-space in the shortest possible time.And the traffic capacity of each edge is affected by several influence factors, which need to be considered in the decision of the most reliable evacuation paths.We list the main influence factors of the restrictedspace as follows.
(a) Edge types:  1 (  ).The properties like the shape and slope of the edge/section decide the different types of the edge.Different edge types have different level effects on evacuee evacuation.Here we consider the scope of the edge, and for any pair of edges <   ,   > and (b) Passibility:  2 (  ).The passibility of the edge is greatly influenced by the obstacles and the vehicle of the evacuees.We assume that the evacuees adopt a uniform activity manner and have the same traveling speed.
For any pair of edges <   ,   > and (e) Travel length: ℎ(  ).The geometric length of the edge.For any pair of edges <   ,   > and Based on these influence factors, the traffic capacity on an edge can be expressed in terms of the equivalent length/edge weight: Thus, the constrained space can be further modeled as an edge-weighted graph  = (, , ).Here we adopt travel length instead of traveling time to be the basis of traffic capacity for the reason that the traveling speed of each evacuee is regarded as a uniform value.

Disaster Spreading Model.
When the disaster or accident happens, the spreading speed and mode have great influence on the decision of evacuation route.Generally, the disaster spreading pattern is determined by several influence elements, e.g., the fire spreading pattern is effected by carbon monoxide in the volume flow of smoke and fire gases , oxygen in the fire place  2 , heat flow, and so on.And the spreading pattern is directly reflected in the speed of disaster on each edge.Generally the spreading pattern can be modeled as a statistical model depending on the disaster itself.
To this end, we value the maximum spreading time duration of the disaster as  s ; i.e., it takes T from the disasters happening to its whole coverage.Note that we assume that all the evacuees can find the feasible evacuation path and their evacuation time consumption is not exceeding   .We also assume that   can be splitted as several timeslots according to the spreading rule and the evacuation time can be measured as the same timeslot.Thus, the disaster spreading pattern can be modeled as a function of F(  , ), where   = (  ,   ,   ) is the temporal space subgraph of  determined by disaster spreading speed in the -th timeslot.And the influence factor of the disaster is as follows.

Problem Definitions.
We consider a Multiobjective Dynamic Route Network Planning problem in a restrictedspace evacuation, which is to design a staged and adjusted construction scheme of evacuation path network from all the evacuees in  to the emergency exit set  in a 3D space graph  = (, , ).
Firstly, the problem is motivated from two points of view: the one is for each evacuee and the goal is to decide its globally ideal escape exit based on evacuation priority and find the escaping path to the ideal exit; the other one is for each escape exit.In most of practical applications, the feasible escape exits have different evacuation priorities; e.g., in underground coal mine evacuation, the auxiliary and main shafts have higher priority than chambers or mobile capsules.Thus, to enhance the utilization rate of each exit, it is needed to assign relatively balanced count of evacuees to the exit such that the case of congesting in some an overdrawn exit can be avoided.
Secondly, the problem has two objectives; one is in time dimension, i.e., minimizing the whole evacuation travel length.Here the whole evacuation travel length is defined as the travel length consumed by the latest evacuee, which can be referred to as the length of the longest escaping path, max ∀  ∈,∀  ∈ ℎ(  ,   ).And the other one is in space dimension, i.e., maximizing the priority-oriented evacuation efficiency.The priority-oriented evacuation efficiency is a global measure representing the balance performance among each exit's utilization efficiency; i.e., a smaller efficiency stands for a higher balance performance.And the priority-oriented evacuation efficiency is defined as the maximum value among the utilization efficiency differences between any pairs of exits, max   ,  ∈ |(  )/|  | − (  )/|  ||, where   represents two meanings, the exit   itself and the set of the sources in  escaping from   finally.
Thirdly, the problem has dynamic constrains, i.e., disaster spreading status, which is the most influential one among the factors for traffic capacity and is varying along with the evacuee evacuation.In each timeslot, the disaster spreading pattern is regarded as an input item of the path planning process and it can be modeled as a function of F(  , ), where   is the temporal version of .
Our problem in restricted-space evacuation is defined as follows, which is a staged problem model.
Find: Find the route network in each timeslot  from all the sources to their optimal destinations, and the whole route network is constructed with two goals, minimizing the whole evacuation travel length max ∀  ∈,∀  ∈ ℎ(  ,   ) and maximizing the priority-oriented evacuation efficiency max   ,  ∈ |(  )/|  | − (  )/|  ||.

Theorem 2. MODRNP Problem in restricted-space evacuation is NP-hard.
Proof.In order to proof the hardness of MODRNP Problem, we firstly decompose the problem in each timeslot as a subproblem and rewrite the subproblem into a mathematical formulation as below: In a 3D connected and edge-weighted graph   = (  ,   ,   ), given two vertex subsets of   , the Then we consider a special case of it: Case (a).∀  ∈ , ∀  ∈   , ℎ(  ,   )s are a uniform value; i.e., the maximum value in (ii) is unified as a constant.Thus, the minimization of the maximum value can be regarded a default implement.Based on Case (a), we transform these isometric paths into the corresponding one-hop edges; i.e.,   is transformed into a bipartite graph with the two vertex sets   and ,   [,   ] as shown in Figure 1.
Case (b).∃ 0 ∈ , y( 0 ) = 0, i.e., the objective measure in (i) minmax Based on the above problem variant, the subproblem in Cases (a) and (b) is to find a subcollection C ⊆ D (here D = { 1 ,  2 , . . .,  || }) such that ⋃   ∈C   =  and the total weight of   in C is minimized.Here we review an important and classical NP-hard problem, Minimum Weighted Set Cover (MWSC) Problem: Given a set X composed of  elements, a collection D of  subsets of X (D = {X 1 , X 2 , . . ., X  }), and each X  ⊆ X (1 ≤  ≤ ) has a weight (X  ), the problem is to find a subcollection C ⊆ D such that ⋃ X  ∈C X  = X and ∑ X  ∈C (X  ) is minimized.It is easy to discover that this special version of each subproblem is equivalent to MWSC problem, which was proved to be NP-hard [28].Therefore, MODRNP Problem is NP-hard in general.

Multiobjective Dynamic Route Network Planning Algorithm
To solve the NP-hard problem MODRNP, we propose a partition strategy (the framework is shown in Figure 2), which has a uniform algorithm and progressive inputs in each timeslot.The division of timeslot is determined by the disaster spreading function F(  , ) and it generates the temporal graph   in timeslot .In each timeslot, we construct the route network according to the following iterations.
In each timeslot , (i) Phase 1: update the space graph according to the spreading pattern,   ← F(  , ), where   can be regarded as a progressive input for the algorithm.(ii) Phase 2: search for the current route network P  in   , which connects the current source set   to the set .
Since the evacuees from the current sources may not reach the destinations  in this timeslot limited by the evacuee speed, it is also needed to extract the path set P *  from   to their destinations in the current timeslot, collected by   .(iii) Phase 3: update   to   as the new source set in the next timeslot and enter timeslot  + 1.
In Phase 2, the searching strategy consists of three steps.Firstly, we find all the shortest paths from each node of   to each node of  in   based on Dijkstra Algorithm and record all the obtained paths as a subgraph    .Secondly, we sort all the paths of    in the nondecreasing order of their lengths to obtain an ordered path set Thirdly, we construct a route network from the current source set   to , whose details are as follows: (1) construct a subgraph    composed of all the shortest paths from each node of   to the whole destinations of  via Dijkstra Algorithm; (2) design a greedy strategy to select the optimal whole destination for each node of   in    based on Minimum Weighted Set Cover and record the selected shortest paths as ; (3) for each node in   , decide its current destination on its path in    and store the corresponding paths into .The construction is under the goals of minimizing the longest path and maximizing the priority-oriented evacuation efficiency.
The former objective is realized by a binary search as shown by the while loop in steps ( 5)- (17).In each iteration of the binary search, for each ℎ  ∈    , starting from  = ⌊(1 + )/2⌋, we temporarily remove all the ℎ ∈    whose length satisfies ℎ(ℎ) ≥ ℎ(ℎ  ) and check if the subgraph of    can still construct a route network without those edges, which is verified by a searching algorithm.If a route network can be constructed, we obtain a route network  on the subgraph of    and proceed with the increased .Otherwise, we keep those edges in    and proceed with the decreased .The latter objective, balancing the evacuation efficiency among all the nodes of , is implemented by the route network construction in each iteration of the binary search: we propose a path planning algorithm, RN-Greedy-MWSC(  , ) (Algorithm 2).To compute evacuation paths from multiple sources to their most suitable destinations, we referred to the classical approximation algorithm with the ratio of ln( + 1) for MWSC problem: select a subset T ∈ F that minimizes the whole redundant weight ((T)/|T ⋂ U|) in each loop, where U stands for the uncovered elements in X; repeat the greedy selection until the set of Ts cover X.Thus, the criteria of greedy selection here are the redundant utilization, which is defined in step (8) of Algorithm 2; i.e., we repeat selecting a best   with the minimum |(  )/|  ⋂ S|| value in each iteration until each node   can be reached by the selected   s.

Theorem 3. The time complexity of Route Network Planning Algorithm for MODRNP Problem is O(n 2 ).
Proof.For the four phases for Algorithm 1, the first one (Step (2)) is constructing all the candidate paths between   and  which has O(n 2 ) time complexity.In the second phase (Step (3)), sorting the candidate paths based on their lengths takes O((|S| * |D|) 2 ) time.Thus, the time cost in this phase is ( log ).For Phase 3 (Steps (5)-( 17)), the binary search for the current route network construction of the auxiliary graph has the time complexity of O(|D| ⋅ log(|S| ⋅ |D|)).In the final phase (Steps (18)-( 21)), the worst-case time complexity is O(|S|).To conclude, the whole algorithm has the time complexity of O(n 2 ).coordinate values.And the edge set  is randomly generated based on the deployment of  to guarantee that the 3D graph is connected: firstly allocate the connecting relationship to each pair of nodes in ; secondly determine whether the graph is connected.If being connected, the edge generation is successful; otherwise, repeat the edge generation process.We consider the following parameters and their variation modes are described as follows:

Performance Evaluation
(i)  = ||: the number of observation nodes in the constrained space; (ii)  = ||: the number of source nodes; (iii)  = ||: the number of destination nodes; (iv) : an edge's priority is decided by the relevance between the edge and the destinations with different For the evaluation of algorithm performances, we compare our algorithm with the shortest path (SP) algorithm in dynamic and static patterns; i.e., we apply the algorithms in each timeslot under F(  , ) for dynamic pattern; meanwhile, we apply these algorithms in static pattern only under the latest status of disaster spreading F(  , ) and without the consideration of the disaster changes.We discuss about the average behavioral characteristics of these algorithms (denoted as Dynamic RN-Greedy-MWSC, Dynamic SP, Static RN-Greedy-MWSC, and static SP).
We will investigate the influence of all the relevant parameters on our algorithms and the parameter setting is based on the scale proportions between the pairs of the space size, escaping evacuee (source) and exits (destination) in practical: (a)  varies from 200 to 550 by the step of 50 and  = 30,  = 6.
For each parameter setting, we run 100 instances and compute their average for evaluation.Based on the parameters and their variation modes, we study the average characteristic behavior of the proposed algorithm and evaluate their average performance in terms of three items in Table 1.

Performance on Priority
Factor.We first observe the impact of , , and  on the priority factor.Overall in the three subfigures of Figure 3, the priority factors of the dynamic and static SP with the variation of parameters are all higher than those of the dynamic and static RN-Greedy-MWSC.And dynamic RN-Greedy-MWSC is obviously lower than the static one during the second half of the variation.
(a) In Figure 3(a), comparing our algorithms with SP algorithms, with the growth of the space scale, the factors in the three compared algorithms are shown in parabolic shapes and the gap between RN-Greedy-MWSCs and SPs becomes bigger.It can be explained as that enlarging the space scale supplies more feasible evacuation paths.Thus, RN-Greedy-MWSCs can perform their advantages on balancing the utilization efficiency among these paths.
Comparing our algorithm in dynamic and static patterns, we can observe that, with the increasing of , the factor's fluctuation in static RN-Greedy-MWSC is more obvious than that in dynamic RN-Greedy-MWSC and the dynamic one is trending downward on the whole.When  > 250, the advantage of dynamic one becomes obvious.It is because when the density of observation vertices deployment becomes larger, the connectivity of the candidate route network is improved, which makes more influence of disaster spreading on each edge.Dynamic RN-Greedy-MWSC considers the influence changes and thus can select paths with better priorityoriented utilization.
(b) As shown in Figure 3(b), comparing our algorithms with SP algorithms, SP algorithms perform higher than ours on the priority factor and present state of rises and falls.It is because that SP algorithms can guarantee minimizing the route lengths but cannot provide the utilization balance of the routes.Thus, with the changes of source number, the maintaining of the utilization balance may have instability.
Comparing our algorithm in dynamic and static patterns, we also find that when  is larger than 15, the variation trends for the two algorithms become similar and dynamic one maintains the advantage.The number of evacuee  raises the fluctuation in the two algorithms with the peak values at  = 30 in dynamic RN-Greedy-MWSC and  = 35 in static RN-Greedy-MWSC.
(c) From Figure 3(c), comparing our algorithms with SP algorithms, under a small number of destinations, SP algorithms have undesirable priority factor values but show improvement when  > 10.It can be explained that the number of the available and shorter paths can be increased with the growth of destination number, which can raise the efficiency of path utilization.The increasing of the destination number is beneficial to improve the utilization balance of the path planning algorithms.
Comparing our algorithm in dynamic and static patterns, we can learn that the priority factor gained from static RN-Greedy-MWSC can be guaranteed into range [1,2] and the dynamic one shows a trend of decline and becomes much closer to 0 along with the rise of the exit number.The reason is that with the adjustment of exit utilization based on disaster spreading, the priority-oriented utilization for each exit can be enhanced.

Performance on Average Length
Factor.Secondly, we investigate the impact of , , and  on the average length factor.Note that, for the SP algorithm, the average length factors of dynamic and static SP are both 0 (= ℎ(  ,   ) − (  ,   )) which are the lower bounds of all the pathfinding algorithms.These three subfigures show that dynamic and static RN-Greedy-MWSCs have gaps with SP on the average length performance.In Figure 4(a), with the growth of due to the consideration of the balanced utilization of each destination, the path assignment in RN-Greedy-MWSC can be adjusted under the change of the path network topological structure.Based on Figures 4(b) and 4(c), it is worth noting that, with the increasing of the number of sources or destinations, the gap between the average length of the escaping network generated by RN-Greedy-MWSC and that by SP tends to be gradually narrowing.
Comparing dynamic RN-Greedy-MWSC with static RN-Greedy-MWSC, the average length factors of the dynamic one with the variation of parameters are all higher than the static one.Based on Figure 4(a), the average length factor of dynamic RN-Greedy-MWSC is in volatility and that of static RN-Greedy-MWSC is gradually raising along with the rise of .It is because, with the enlarging of the scale of restrictedspace, the global lengths of paths between each pair of source and destination will be increased.From Figures 4(b) and 4(c), we can observe that the number of source nodes  has more influence on this factor of our solutions than the number of destination nodes .With the increasing of , the factor of two algorithms presents parabolic shape in the range [10,25] and the gap between the two algorithms reaches its maximum when  = 20.And with the growth of , the factors of two algorithm both show a parabolic shape in a reasonable range [0, 0.35].

Performance on Global Length
Factor.Similar to the average length factor, the global length factors of dynamic and static SPs are both 0 (= max ℎ(  ,   ) − max (  ,   )) which are the lower bounds of our algorithms.Seen from Figure 5(a), the difference between the global length obtained from RN-Greedy-MWSCs and that from SPs is getting smaller with the rising of the space scale.And when  is larger than 450, the difference is getting 0 which indicates our strategy can enhance the utilization efficiency of the escaping exits with the minimization guarantee of the global length.On the contrary, as shown in Figures 5(b) and 5(c), the difference becomes larger with the increasing of the source number and could be controlled in a steady state when  ≥ 35 or  ≥ 10.This reason is that, with the rise of the source number or destination number, the utilization efficiency of shorter paths can be increased, which can narrow the gap of global path length between SP algorithms and RN-Greedy-MWSCs.
Comparing dynamic RN-Greedy-MWSC with static RN-Greedy-MWSC, seen from Figures 5(a) and 5(c), the difference between the global escaping time consumption obtained from dynamic solution and that from static scheme is getting smaller with the enlarging of space scale or the rising of the evacuee number.And when  is larger than 450 and  ≥ 14, the difference is approaching 0 which indicates our strategies can enhance the priority utilization efficiency of multiple escaping exits with the complication of the space structure.On the contrary, as shown in Figure 5(b), the difference between the two algorithms is getting larger since the source count is larger than 25, which is also caused by the dynamic influence of the disaster.
Via the performance evaluations on the three measurement factors, we can draw the conclusions.
(A) The path planning algorithms: firstly, RN-Greedy-MWSC algorithms have more advantages on path utilization than SP algorithms.And the dynamic RN-Greedy-MWSC outperforms the static one in terms of priority factor when the scale of a restricted-space, the scale of sources, and destinations are relatively large; secondly, when evaluating on the average and global path length, the advantage of dynamic RN-Greedy-MWSC maintains signification.
(B) The influence factors: the influence of the scale of sources has more influence on the algorithm performance than that of destinations and space scale.

Application Instance Experiment and Discussion
. We further perform an instance experiment on our algorithm in a representative kind of restricted-space application scenarios, underground coal mine water inrush.Here we adopt the 3D model data of a mine area in Zhangjiakou, Hebei Province of China from Kailuan Group Company.We preprocess and transform the original spatial data of the model; i.e., the vertex set in  is composed of all the traverse points, and the edge set of  maintains consistent direction with the original laneway.Water spreading status is obtained by the water spreading algorithm in [29]; i.e., in each timeslot, the prediction data of the water level on all observation points are the input of our algorithms.According to practical requirements, the priorities of three feasible exits, auxiliary shaft (d1), main shaft (d2), and air shaft(d3), are valued as follows: (1) = 5, (2) = 2, and (3) = 1.5.
In the visualization results in Figure 6, we present the paths from the 10 evacuation source positions to their temporal destinations from timeslots 4 to 11, and finally all the evacuees successfully escape from the 3 red exits in the last subfigure.Note that the points in blue stand for the submerged part by water inrush, which are all avoided from by the escapers in our algorithms.In Table 2, the evaluation criteria, priority factor, average length factor, and global length factor are represented by PF, ALF, and GLF in short.The number of paths reaching these exits is denoted as d1, d2, and d3, respectively.The parameter results show that the dynamic algorithm outperforms the static one in terms of priority factor in a practical application scene with large scale.
To conclude, our algorithm has advantages on the performance of the escaping exit priority utilization and the evacuation time consumption.And when the scale of space, evacuee, and exits are relatively high, the advantage will be significant.Furthermore, the parameter, the number of evacuees, has more influence on our algorithm than the number of exits on the length factor.Therefore, our strategy can enhance the priority utilization efficiency of the escaping  exits and guarantee a tolerable range of the global escaping time consumption.

Conclusions
In this paper, we propose an optimization problem for emergency evacuation planning in restricted-space scenarios, the Multiobjective Dynamic Route Network Planning (MOD-RNP) problem.To minimize the whole evacuation delay and maximize the priority-oriented evacuation efficiency, we design a dynamic and global-optimized strategy to construct the evacuation route network, based on Minimum Weighted Set Cover (MWSC).Based on the performance evaluation, our algorithm outperforms the shortest path algorithm in dynamic pattern on the priority utilization efficiency of the escaping exits, and it can guarantee that the deviation of the global travel length from the shortest one is in a reasonable range.Furthermore, we believe there are lots of room for further investigation and we plan to design distributed dynamic route network planning strategies to contribute or reorganize route network in real time.

Figure 2 :
Figure 2: The framework of Route Network Planning Algorithm.

Figure 3 :
Figure 3: Performance on the priority factor.

Figure 4 :
Figure 4: Performance on the average length factor.

Figure 5 :
Figure 5: Performance on the global length factor.

Figure 6 :
Figure 6: The visualization results of dynamic algorithm.
The factor is generated from the incident relationship with the destinations in , i.e., if edge   is an incident edge with more than one   in , (  ) = max   ∈ and   ∈  r(  ).For any pair of edges <   ,   > and <   ,   >, (<   , c) Exit escaping priority: (  ).Due to the exits' security levels, structural characteristics, or functional characteristics, different exits have different escaping priorities {(  ) | ∀  ∈ } which have influence on the edge's traffic.(d)Evacuation priority: (  ).
set   and the destination set , and 's priority set of (), the problem is to find a path network from   to  such that (i) the maximum difference among 's utilization efficiencies, max   ,  ∈ |(  )/|  |−(  )/ | source  ||, can be minimized, where   is not only a destination in , but also a set of nodes in   and connected to   by the constructed paths, (ii) the maximum path length among all sources and destinations, max ∀  ∈,∀  ∈  ℎ(  ,   ), can be minimized.