Abstract
Apart from private traffic, the evacuation of transitdependent population is also an essential component of emergency preparedness, especially under nonotice evacuation scenarios with limit evacuation horizon. In literature, most busbased evacuation models for nonotice evacuation are established under implicit assumptions of uniform evacuation horizon among different pickup locations or fixed bus fleet in the evacuation area. These constraints will distance their models from realworld situations, where evacuation horizon is various due to spatial distribution of pickup locations and fleet size of bus available for allocation will increase over time in nonotice evacuation. This research presents a riskbased bus schedule model which is differentiated from the vehicle routing problem (VRP) and bus evacuation problem (BEP) in literature, including the objective and the timedependent parameters. A quantified definition of evacuation risk for pickup location with concerns of disaster dynamics and timevarying supplydemand conditions is proposed in this paper as a criterion for bus allocation, also acting as a reflection of social equity to some extent. A notion of timeevolving disadvantageous evacuation units (DEU) is introduced to represent the pickup locations selected for bus allocation with limited resource. The binary integer linear programming (BILP) named riskbased bus schedule model incorporated with DEU notion can provide a reference for resource allocation in stage of both evacuation planning and operation for transitdependent population. The proposed model structure can effectively capture the changes of evacuation risk among pickup locations over time to realize realtime bus schedule. Numerical experiments are conducted using the transportation network of the city of Xi’an, China, to test the performance of the model. The applicability and comparison of different bus evacuation models are also discussed in this paper. This research provides insights into dealing with disaster dynamics and timevarying supply conditions in realistic busbased nonotice evacuation operations.
1. Introduction
Apart from private traffic, the evacuation of transitdependent population is also an essential component of emergency preparedness, especially under nonotice evacuation scenarios with dynamic disaster impacts. The researches on busbased evacuation have drawn worldwide attention after Hurricane Katrina [1–10]. Compared with short notice evacuation, the planning of busbased evacuation under nonotice scenario is more difficult because of limited evacuation horizon and insufficient preparedness. As noted by the TRB special report 294 [11], although the probability of occurrence of nonotice disasters is low, the spatiotemporal randomness and catastrophic impact justify the rationality of research on such disasters. Modeling for busbased evacuation planning under nonotice evacuation with limited resources has drawn more attention in recent years [1, 7, 12–20]. The literature on busbased evacuation mainly solves problems of pickup location selection (Kulshrestha et al., 2014) [4], route optimization (Tunc et al., 2011) [7], and resource allocation (Kulshrestha et al., 2014) [1]. The busbased evacuation problem in this paper is more like a dynamic resource allocation (DRA) problem, which has become an active area in recent years. He et al. (2015) modeled the problem of DRA of movable resources for largescale evacuation as a mixed integer linear program and developed a solution algorithm based on the Benders decomposition technique to solve the proposed model.
Existing busbased evacuation models solve the bus evacuation problem (BEP) by minimizing evacuation duration [7], Goerigk and Grün 2013 [14]. While under nonotice evacuation, the primary concern is not to reduce the network clearance time, but instead to maximize the safe evacuation numbers or minimize the total amount of casualties or exposure to the threat area [13, 21]. Meanwhile, since the available number of buses is limited, the social equity should also be considered in the bus schedule model. Recently, Aalami and Kattan [1] proposed a DRA model in busbased evacuation with consideration of “proportional fairness” and developed an algorithm based on Lagrange dual approach to find proportional fair allocation of resources. However, both the bus fleet size and bus schedule are unchangeable and a complex objective function is utilized in their model; therefore only the demand aspect is considered in resource allocation. Further, the bus schedule scheme should be adjustable to realtime information. It is noteworthy that, in busbased nonotice evacuation, the supply conditions and disaster impacts may change dynamically during evacuation procession based on evacuation response and disaster analysis [18, 20]. However, most existing busbased evacuation models are established under implicit assumptions of uniform evacuation horizon among different pickup locations and fixed bus fleet in the evacuation network. These constraints will distance their models from realworld situations, where evacuation horizon is various due to spatial distribution of pickup locations and fleet size of bus available for allocation will increase over time in nonotice evacuation. Therefore, the spread of disaster and bus supply condition combined with evacuation demand should be incorporated in the busbased evacuation modeling to better adapt to realistic situation.
To bridge the aforementioned gaps in incorporation of dynamic disaster characteristics and supplydemand conditions in busbased evacuation, the busbased evacuation problem under nonotice evacuation is addressed as a riskbased nonotice bus schedule problem (RNBSP) proposed in this paper. The quantified evacuation risk of pickup locations derived from the interactions of aforementioned factors acts as criterion of bus allocation with consideration of social equity. To the best of the authors’ knowledge, there is no research on busbased evacuation modeling under nonotice evacuation incorporated with dynamics of disaster impacts and supplydemand conditions.
The rest of the paper is organized as follows. The preliminaries for modeling the problem are provided in the next section. Section 3 describes the evacuation risk assessment of busbased evacuation in this paper. The method of determination of initial disadvantageous evacuation units (DEU) is proposed in Section 4; Section 5 introduces the modeling and solution method of RNBSP. Numerical results and discussions are presented in Section 6. The study is summarized in Section 7.
2. Preliminaries
2.1. Evacuation Scenarios and Scope of Research
Consider the scenario when a disaster such as a hazmat explosion strikes an urban area: citizens in the affected area (predefined by the disaster characteristics and wind directions) are required to evacuate the area within a very short time horizon, e.g., less than an hour. Under this situation, evacuees with access to automobile will evacuate by personal vehicles, while, for transitdependent citizens, public transportation systems coordinated by governmental agencies might be their only available means of evacuation. Although there might be other mode choices available for the transitdependent evacuees (e.g., walking and waterway in 9/11 and carpools in Hurricane Katrina), this study focuses on busbased evacuation under nonotice evacuation scenario where evacuees are guided to designated pickup locations (e.g., bus termini, landmarks like park, and schools) and depend on public transit service. It is understandable that deadlines of different locations within the area are different through the endangered area because of spatial distribution and disaster propagation characteristics, interpreted as spatial difference of evacuation risk associated with characteristics of disaster among pickup locations within the endangered area.
Under this circumstance, buses within and around evacuation area will dismiss their scheduled path and follow different paths to the temporary scheduled pickup locations as required. It is noteworthy that not only the buses within the evacuation area is available for allocation to associated pickup locations, the buses operated outside the evacuation area would also be assigned to endangered area with a time delay. That is, at the starting point of nonotice evacuation, the available bus fleet only contains the buses currently in the evacuation area, but the size of fleet will increase over time as more buses out of the evacuation area arriving in the evacuation area. Therefore, a dynamic supply condition is considered based on available resources for allocation over time.
2.2. Analysis of Compliance Rate and DecisionMaking Process of Evacuees under Emergencies with Different Type
Evacuation demand for each pickup location should be estimated first before resource allocation. A possible direction for demand estimation for each predefined pickup location from demand side is discussed in this subsection. Although the evacuation demand at each pickup location is assumed as a model input in this study, it is noteworthy that the uncertainty of the evacuees’ compliances in the evacuation plan is a main concern in real operation. Currently, the main flow of studies on compliance of evacuees is focused on evacuation of private vehicles, i.e., compliance behavior of car owners [22–25]. Unless it is a mandatory evacuation performed by the law enforcement units, it is very likely to have a certain proportion of the population that would not comply with the order or the recommendations. Compliance rate is an aggregate behavioral characteristic that indicates the percentage of evacuees that follow a recommendation or order by EMAs. It is a vital influential factor for the effectiveness of control strategies during evacuation operation. According to Abdelgawad and Abdulhai [8], noncompliance percentage acts as a sensitivity parameter in their traffic flow optimization model for mass evacuation.
Compliance of evacuees in evacuation includes compliance with evacuation warnings, departure time, destination, and evacuation route, etc. Evacuees’ compliance with evacuation warning would impact the determination of total amount of traffic in the network. Risk perception is one of the key factors in understanding the evacuation decisionmaking process. However, the heterogeneity among population will turn out different level of perceived risk with given warning massage and circumstance. Furthermore, current practices in conducting evacuations are very general; for example, same set of instructions are given to evacuees in the endangered area containing available evacuation time and evacuation route (for carowners). However, evacuees need much more detailed instructions, which are thus more difficult to disseminate [26]. Current advances in communications offer several feasible options in this arena, where providing customized instructions is a feasible way to improve the level of compliance in the case of evacuations. Normally, under evacuations with relative longer notice, evacuees will select the departure time and pickup location reasonably by service rate, i.e., the number of evacuees transported out from a pickup location per unit time, queue size i.e., the number of evacuees waiting at the pickup location currently, and arrival rate, i.e. the number of evacuees arriving at the pickup location per unit time; it is regarded as a tradeoff between perceived cost and risk exposure. Under this scenario, the behavior of evacuees might be not in accordance with the recommendation of emergency management agencies (EMAs) because of heterogeneity among population. Evacuees will compare the recommendation with the information that they received from other sources and make their own decision after deep consideration. The accuracy of demand estimation in terms of evacuees’ compliance of recommendations can directly influence the performance of evacuation plan. Therefore, the evacuees’ compliance is probably one of main factors impacting the evacuation risk. However, under scenario of nonotice evacuation, evacuees intend to get out of the endangered area as soon as possible out of the fear of risk; thus the departure time can be perceived as the time of disaster occurrence. Meanwhile, the decisionmaking process of evacuees may depend on simple principle with approximate judgment on a few factors and alternatives, due to the time pressure associated with high level of emergency caused by the disaster [25]. For evacuees with access to automobiles, the decisionmaking process is mainly about route choice; here, evacuees’ compliance to route recommendation will impact the traffic assignment of evacuation network and further impact the evacuation operation, while, for transitdependent evacuees, the decision is mainly about selection of pickup locations. Transitdependent population is vulnerable under nonotice evacuation, especially for situations with low information penetration, where disasterrelated information provided by media associated with personal vehicles is not accessible to these groups of people. Compliance of transitdependent population will impact the initial demand distribution among pickup locations, which acts as one of the determinants of bus resource allocation. Therefore, given the limitations of mobility of transitdependent population, EMAs should assign evacuees with most convenient pickup locations to obtain higher compliance rate. The initial distribution of transitdependent population can be obtained from historical statistics and daily trip data. It is assumed that evacuees have no idea of the current situation of each pickup location (e.g., queue size, service rate, and arrival rate) near their current location under nonotice evacuation. Whether a pickup location selected by an evacuee or not is only relied on the intuitionistic perceived cost comparing to other pickup locations, which can be simply interpreted as seeking nearest pickup locations of no information received. Since the decisionmaking process of evacuees may depend on simple principle with approximate judgment on a few factors and alternative, we circumvent the estimation of compliance rate of evacuees by directly assigning them corresponding pickup locations based on the minimum distance principle. That is, it is assumed that the evacuees would choose the nearest pickup location as no other valuable information is available.
The decisionmaking mechanism highlights the possibility of using simple principles for decisionmaking during nonotice evacuation. However, under shortnotice evacuation scenario, transitdependent evacuees’ compliance would be a vital factor that impacts demand level of each pickup location. Therefore, compliance remains a practical issue to be addressed in the future. In evacuation operation, the buses usually have fixed routes compared to private vehicles. In future, the compliance of bus drivers on evacuation route should be studied, as drivers might change the assigned route based on their judgement on traffic conditions.
2.3. Parameters and Problem Setting
The evacuation area can be abstracted to topological structure of a directed network G (N, A), where N and A, respectively, denote the set of nodes and arcs. Each node , representing a location, is assumed to cover a predefined geographical area in its vicinity. is composed of two subsets of nodes; P, a set of predefined pickup locations serving associated geographical area in its vicinity; and S, a set of shelters . The network is considered at a set T of discrete time intervals , where , is a constant increment of time, and let represent the starting point of evacuation. The end of evacuation is defined by . It is reasonable to set equivalent to the longest evacuation time horizon of the pickup locations in the endangered network. Thus, there are periods in total evacuation horizon. Here, to capture the dynamics of supply condition and guarantee relative accuracy of travel time, a time increment of 1 minute is taken in this paper.
Since the evacuation demand increases explosively under nonotice evacuation, the arrival rate of evacuees is higher than service rate because of limit resources available for allocation, and thus the total demand is given as a constant instead of a monotone increasing function over time without violation of reality. We denote the total number of evacuees with , where is demands of pickup location . The value of is acquired by assigning evacuees to certain pickup locations based on the minimum distance. The approximate distribution of population is obtained by historical statistics of the area. Since each pickup location is expected to receive a relatively large number of evacuees, it is assumed that each vehicle will always be fully loaded at its scheduled pickup location, and thus it will always travel directly between its pickup facility and the shelter (instead of traveling to other locations for additional pickups). We also assume the number of evacuees at every pickup location is known, in terms of integer multiples of bus capacities denoted by c. The capacity of each shelter is uncapacitated as the aim of evacuation under this circumstance is to transport the evacuees out of endangered area timely; thus the destination is not necessary a place that can provide accommodation to evacuees. is a given parameter denoting the total number of buses available for assignment at time , determined by the initial number of buses in endangered area and arrival rate of buses outside the area. Although the evacuation time horizon for nonotice evacuation is limited, multiple trips between shelters and pickup locations should be concerned instead of a onetime transportation; represents the estimated round trip time between pickup location i and shelter j, in terms of integer multiples of time increment. In the case of various over time (multimodal evacuation with influence of private vehicles), the notation is used to refer to at time t. is introduced to denote whether bus operates between and at time . If the answer is yes, equals one; otherwise, equals zero. can denote the number of evacuees who can leave in time unit t, i.e., the service rate in time unit t.
3. Evacuation Risk Assessment for TransitDependent Population
3.1. Related Work
In this paper, evacuation risk is utilized to rank the priority of pickup locations for bus allocation with limited resources. The principle for bus assignment in this paper is to allocate limit resources where most in need is based on current situation (disaster dynamics, supply condition, level of demand). The definition of evacuation risk is still obscure in literature, especially for busbased evacuation. Generally, the evacuation risk is notional according to the conceptual danger as research needs. The evacuation risk in this paper is defined in terms of whether the transitdependent population of a pickup location can be safely evacuated before the impact of disaster on it under current supply condition. It is an integration of disaster dynamics, supply, and demand conditions. MurrayTuite and Wolshon [27] provide a comprehensive overview of evacuation transportation modeling. The majority of researches are focused on traffic management strategies from supply side, while the disaster characteristics are treated as a notional concept, except for Hsu and Peeta [28], they incorporate disaster characteristics into risk measurement for evacuation of private vehicles in a dynamic manner. The risk of a location is mainly determined by the potential disaster impact on it associated with geographical features in the natural hazards management literature. Method of “Hazard Scores” has been proposed to rank level of risk for planning purposes [29]. Church and Cova [30] quantify evacuation risk in terms of network and demographic characteristics to produce maps of evacuation risk or vulnerability. However, the measurement method of risk with concerns of various factors combination is questionable. Social vulnerability has also been introduced into risk calculation to represent the level of needs for evacuation [31]. Most evacuation risk proposed in literature is addressed in a static planning context mainly based on disaster analysis and has key gaps to adapt its application on dynamic evacuation operations from transportation side. Further, as the measures are developed as research needs, most existing measures for risk assessment are specific to particular disaster types and lack general applicability. As mentioned in Hsu and Peeta [28], it is a challenge to measure risk robustly to the affected region using a unified model considering different disaster types and dynamics. The evacuation risk measurement for busbased evacuation in this paper comprehensively considering dynamics of disaster, demand, and supply condition is still missing in literature. The riskbased bus schedule principle is of significance as it put social equity in priority instead of convenience and efficiency, i.e., serving the population in need first rather than population with shortest roundtrip cost.
3.2. Evacuation Risk Formulation
The main challenge of nonotice busbased evacuation is formed by the evolving characteristics of three factors: disaster dynamics, demand level, and supply condition. The disaster dynamics have tempospatial influence on evacuation risk among the pickup locations within the evacuation area as the spread of disaster impact. The demand level at different pickup locations will determine the time duration for evacuation with given supply condition, and vice versa. With aim of developing a formulation of evacuation risk measurement that is consistent with the busbased evacuation scenario, a formulation of evacuation risk that can capture the dynamics of disaster characteristics, demand, and supply condition is proposed in this subsection. Therefore, to obtain the evacuation risk of transitdependent population in a comprehensive manner, the evacuation risk should be formulated by proxy parameters of the aforementioned three factors.
The proposed measure for evacuation risk in the RNBSP context accounts for disaster characteristics using the notion of a residual evacuation time. Residual evacuation time of a location refers to the estimated length of time available for evacuation before the arrival of disaster impact on that location. In this paper, means the available evacuation time left for location at time . Residual evacuation time is assumed to be known based on analysis of disaster characteristics, which is obtained by meteorological department or other related agencies with assistance of relevant prediction models of disaster impact and spread, such as meteorological, hydrologic, or plume dispersion models, etc. It characterizes disasterrelated risk, as the potential for significantly increased casualties if evacuation of transitdependent population is not finished by the end of its residual evacuation time, i.e. . It is understandable that the risk increases as the end of residual evacuation time approaching. At the starting point of evacuation, , the residual evacuation time of can be expressed as follows:If the propagation speed is constant, the residual evacuation time of at time unit , is formulated as follows:orIf the propagation speed is various over time, the residual evacuation time of at time unit , can be written as follows:where is distance between and location of disaster occurrence. , , and are the propagation speed of disaster impact at different time unit. denotes tolerable exposure time, which is predefined based on disaster type.
In the RNBSP context, utilization of residual evacuation time as proxy of disaster impact not only reflects the time available to evacuate of a disaster relative to the location of interest, but also manifests the spatial differentiation of disasterinduced risk across various pickup locations in the endangered area. Residual evacuation time is negative related to evacuation risk of a location.
The evacuation risk formulation considers the effects of demand and supply conditions by a notion of evacuation completion time. Evacuation completion time for a location is defined as the time required to evacuate all the evacuees left in that location at current time. The estimated evacuation completion time required to evacuate citizens left in at time , , can be written as follows:
where is the number of evacuees left at at time .
Both evacuation completion time and residual evacuation time are timedependent itself measured by same unit, thereby enabling the incorporation of the evolving traffic dynamics resulting from the demand–supply–performance interactions. As mentioned before, the evacuation risk in this paper is defined in terms of whether the transitdependent population of a pickup location can be safely evacuated before the impact of disaster on it with current supply condition. Therefore, evacuation risk of at time without adding new arrived buses, , can be expressed as follows:
It is noticeable that is a nonlinear component in risk formulation, which will dramatically increase the computation complexity of the problem. Thus, instead of representing risk in terms of time, we transform (7) into (8) without changing the generality. where denotes the demands left in at time , and specifies the number of evacuees that can escape from the disaster before based on current supply condition. The proxy of evacuation risk in (8) is in terms of the number of evacuees instead of evacuation time. If the value of is nonpositive, then all the evacuees in are able to leave the area in time under current supply condition. Otherwise, the value of is equal to the number of casualties. After assignment of available bus fleet at , the evacuation risk will change into :
The linear structure of the risk formulation can capture the effect of bus assignment on evacuation risk.
4. Disadvantageous Evacuation Units (DEU) Determination
In this paper, a notion of disadvantageous evacuation units (DEU) is introduced to represent the pickup locations selected for bus assignment based on ranking of evacuation risk; in other words, consists of a set of with highest evacuation risk at time unit t. Busbased evacuation in this paper is more like a dynamic resource allocation (DRA) problem; it is a supply side operation instead of action taken by demand side (evacuation of private vehicles), although a pickup location represents associated geographical area in its vicinity, it is not a graph partition problem, thus contiguity constraint is not a necessity for DEU determination [32]. The confusion of evacuees in automobile evacuation scenarios mentioned in Hsu and Peeta [28] is also perceived not happening in busbased evacuation scenarios because (a) the distance between two pickup locations is relatively far for walking distance and the available information is limited under nonotice evacuation; (b) evacuation of these citizens completely depends on transit system instead of themselves. Thus, the areas covered by pickup locations are allowed to be discrete as long as they can better capture the evacuation risk of different pickup locations. Unlike automobile evacuation operation [28], the busbased evacuation is totally conducted by supply side (EMA); thus it can perfectly adapt the bus schedule scheme to minimize evacuation risk of the total endangered area for whole evacuation horizon.
4.1. Initial DEU Determination
In this section, DEU formulation at the starting point of evacuation is developed, acting as the initial resource allocation under RNBSP context. The dynamic expanding DEU over time is equivalent to the dynamic bus schedule scheme in evacuation operation. is set of available buses that can be assigned to pickup locations to assist evacuation operation at initial time stage, is the total number of bus available for assignment and , and is number of pickup locations. Superscript is omitted for convenience.
Decision Variables. is binary variable that equals one if pickup location is selected into DEU at ; otherwise, it equals zero. .
is binary variable that equals 1 if bus is assigned to operate between and at ; otherwise, it equals zero. .
The determination of DEU can be formulated in the following manner:where
Equation (11) requires to be greater than or equal to the maximum risk of any , which is then minimized by the objective function (10), i.e., minimizing the highest risk. Because of this, we will refer to this objective function as the “min–max” objective. The above model can be utilized to identify pickup locations with highest evacuation risk based on limited supply condition of the endangered area. The aim of the model is to select a set of with higher risk compared to others and decrease the evacuation risk by assigning buses to them. It can be treated as an iteration procedure that continuously decrease the current highest risk to second highest risk with current available buses . Therefore, the objective function is developed by minimizing the maximization of evacuation risk of pickup locations within the evacuation area after bus assignment. For better understanding the objective function, the set of pickup locations can be settled in xy rectangular coordinates, and let x coordinates be pickup location number and y coordinates be corresponding evacuation risk. Linking the nodes by number, we can get a line chart. The objective is to shave the upper peak value of the chart and obtain a chart with lowest upper peak by assignment of bus. Equation (12) allows a bus to make at most one trip at a time. Equation (13) specifies that the total number of buses allocated to pickup locations does not exceed . Equation (14) is a social equity constraint, which specifies that at least one bus should be assigned to each pickup location regardless of risk. Equation (15) ensures that supply should not surpass demand. Equations (16) and (17) together explain that is selected into DEU if any bus is assigned to it. is a big number bigger than . Equations (18) and (19) are the logical binary restrictions on variables and .
4.2. Numerical Experiments
A simple network is given in Figure 1 to test the validity of the model. Only one shelter is determined for simplicity. The travel time for each bold link in Figure 1 is assumed to be identical as 5 min. The evacuation demand, residual evacuation time, number of buses, and roundtrip time of pickup locations at initial stage are given in Table 1. The schedule scheme of different supply condition is also shown in Table 1, calculated by the model above; the bold and italic numbers denote pickup locations for each added bus to the evacuation network.
The value of is nonnegative for all ; i.e., all evacuees can be safely evacuated, after assignment of 32 buses into the current network. The initial DEU determined for different initial supply condition (given number of ) is shown as the bold and italic numbers in Table 1. For example, for =17, the initial DEU consists of ,,,, and DEU=, when =10. The bold numbers interpret the pickup locations for added buses under different supply condition. The DEU formulation in this section is a static resource allocation deployment to understand the principle and mechanism of DEU in busbased evacuation. However, in real operation, the bus schedule scheme will change dynamically by the dynamics of disaster, supply, and demand conditions reflected on risk. Actually, the real number of buses required is less than the static model as a bus can be assigned to other pickup locations after the evacuation of current served pickup location is finished under dynamic situation. The principle for DEUbased bus schedule is to assign available buses to pickup locations with highest evacuation risk in any time unit considering dynamics of disaster, supply, and demand conditions. The dynamic model formulation of RNBSP and the solution approach will be introduced in next section.
5. RNBSP Formulation and Solution Method
5.1. RNBSP Formulation
This section develops the RNBSP formulation of a dynamic bus schedule scheme for busbased evacuation operation. The model aims to optimize the bus schedule scheme according to the dynamic changing evacuation risk among pickup locations.
5.1.1. Objective Function
With consideration of time attribute in dynamic resource assignment, the evacuation risk denotes the evacuation risk of pickup location i at t after assignment of current available buses in the evacuation network. Since the evacuation operation lasts for the entire evacuation horizon, represented by , the RNBSP is equivalent to determination of dynamic expanding DEU during time units by assignment of the buses available to pickup locations with highest evacuation risk over t. After considering time attribute, the objective function is developed as follows:where ; the expression of is given in (9). The objective function is to minimize the risk of with highest evacuation risk for from the start of evacuation operation.
5.1.2. Dynamics of Supply Condition
In the model proposed in Section 4, the supply of bus B is treated as a given input to identify the initial DEU and determine the number of buses required to evacuate the transitcitizen in time. In realistic operation, the bus will be assigned to the evacuation network continuously over time. At the initial stage, only the buses within the endangered area can serve the pickup locations immediately, with process of the evacuation operation, buses relative far from the evacuation area can also arrive in the pickup locations and provide service. The dynamic set of supply is represented by . contains the set of new arrival buses and the buses already in the evacuation network . It is a timedependent parameter based on the travel time determined by distance between the assigned pickup location and buses outside the endangered area and evolving traffic conditions over time that obtained by exogenous models. , , the relationship can be expressed as follows:This is the core equation of supply dynamics in proposed bus schedule model.
5.1.3. Summarized RNBSP Formulation
In summary, the RNBSP for an endangered area is performed by solving the following formulation.
Objectivewhere can refer to (9), and refers to (6).
Equation (23) is constraint of dynamic supply condition. Equation (24) requires to be greater than or equal to the maximum . Equation (25) ensures that a bus can make at most one trip at a time. Equation (26) specifies that at least one bus should be assigned to each pickup location when evacuation starts. Equation (27) implies that all evacuees might not depend on bus to escape from the endangered area because of high level of urgency under nonotice evacuation scenario. Equation (28) specifies that a bus can be assigned again only after it finishes the current tour. Equations (29) and (30) explain the relationship between and . Equation (31) specifies that no bus will be assigned to after ; we can decrease the search range effectively by adding this constraint. The last two constraints are logical binary restrictions on variables and .
5.2. Solution Method
For , if , then all the evacuees are safely evacuated. Actually, for , as long as when , the safe evacuation is guaranteed. The mechanism of the objective function is divided into two steps: (a) identifying with highest evacuation risk before its deadline at current stage; (b) consecutively, assigning buses to it until it is not the riskiest one, then assigning the remainder to the pickup location with currently highest risk (the former secondhighest one), and iterating this procedure until no bus is available at current time stage. This operation can ensure that the highest casualties of different are always the lowest. In scenario with ideal bus fleet size, the evacuation risk is almost identical among pickup locations. Therefore, each pickup location is treated equally with consideration of social equity.
The aim of the RNBSP in this paper is to provide recommendations of realtime bus schedule scheme for EMA. The problem changes into a BILP after the transformation of risk formulation in Section 3. Under nonotice evacuation scenario, bus schedule scheme is required to be realtime adjustable. The solution steps of the RNBSP formulation can be summarized as in Algorithm 1.

Based on the calculated evacuation risk over the endangered area, the solution procedure for DEU determination is implemented in a timedependent operational framework as illustrated in Figure 2. The evacuation operation begins in the first stage (t = 1) with the initial conditions of the evacuation network. The EMAs determine DEU continuously by solving the corresponding RNBSP. The shadowed part in Figure 2 is the mechanism of DEU determination for each t.
6. Numerical Experiments and Discussion
6.1. Objectives of Numerical Experiments
Numerical experiments are conducted in this section to test the performance of the dynamic riskbased bus schedule scheme under nonotice evacuation scenarios. These experiments focus on illustrating the practical applications of the proposed approach for DEU determination and bus schedule scheme in the operational context, where evacuation risk over the study region needs to be captured in a dynamic manner. In the experiments, different disaster impact scenarios are created to test the robustness of the proposed approach to the related variability.
6.2. Numerical Experiments Scenario Descriptions and Operations
The transportation network of Xi’an city, China, is used in numerical experiment (see in Figure 3(a)); the research area is located in the southern part of the Xi’an city. In total, 238 bus stops are located on the daily route of transit vehicles in the research area. As mentioned in Bish [7], the runtime of minmax objective grows quickly with problem size. Thus, the number of pickup locations for selection should be controlled; i.e., the evacuation demands in each pickup location are large enough for multiple times of bus capacity. This precondition is in accordance with the singlevisit singletrip assumption. Three scenarios are created to test the robustness of the proposed riskbased bus schedule model under disasters with different dissemination pattern and locations of occurrence, described as follows: Scenario A: A vehicle carrying hazardous materials explodes and creates a hazmat disaster in the southeast part of the network with diffusive impact uniformly spreading throughout the network in a radial pattern.10 among 51 bus stops are selected as pickup locations in Figure 3(b). Scenario B: The same incident happened in the central part of the network. 9 among 42 bus stops are selected as pickup locations in Figure 3(c). Scenario C: The same incident happened in the middle part of the network. The spread of hazmat is affected by the wind direction (from left to right, 0.7km/h); 9 among 44 bus stops are selected as pickup locations in Figure 3(d).
(a)
(b)
(c)
(d)
The allowable exposure time duration of hazmat is assumed to be 15 min, after which the risk of poisoning and casualties will increase significantly. The diffusion speed of hazmat is assumed to be 3km/h and the effect radius is 5 km in this numerical experiment. Residual evacuation horizon for different pickup locations in the affected area can be calculated by distance between the pickup locations and the location of disaster occurrence divided by the spread speed of the disaster impact, then plus the value with allowable exposure time duration. Evacuation demand of transitdependent population can be assumed as model input. The density of traffic in the network is collected from an output file that is generated by DynasmartP when simulating the abovementioned disaster scenarios. The density output file contains the traffic density on all road segments in every time unit. The travel time required to traverse a path is calculated by speed of traveling on the path and the length of the path.
The time interval is set as 5 minutes for each t, and extra buses will be assigned to the evacuation network within ; i.e., no more buses can be assigned to assist evacuation after 30min. The value of is flexible based on the user’s need or the realistic network supply condition. is generated by the software DYNASMARTP with given input of personal vehicles in the network. The population of transitdependent evacuees of each pickup location is generated by six demand levels of 1000, 2000, 3000, 4000, 5000, and 6000. The initial number of buses available for assignment is 20, and for each the bus assigned to evacuation network is 3, 5, 7, 9 for different level of resource. The capacity of bus is 50.
6.3. Criterion of Performance and Results
Since the RNBSP is focused on evacuation with high level of emergency, i.e., the residual evacuation time is so short that all of the population might not be evacuated using the current available bus fleet, the number of evacuees that can be saved by bus fleet during the evacuation time horizon is a valid criterion of model performance, similar to the percentage of total evacuees saved by available bus fleet:
It is noteworthy that the social equity is an important principle in evacuation resource allocation, especially in evacuation with high level of urgency. Social equity is reflected in terms of evacuation risk in this paper. We introduce a notion of social injustice index as another criterion for model performance, expressed as follows:where and , respectively, represent the number of evacuees saved by bus fleet and total demands at . The formulation of can reflect the degree of unevenness in terms of variance in percentage of evacuees saved by bus fleet among various pickup locations. It is easy to obtain the maximum value of is 0.5, under an extreme situation of two pickup locations with, respectively, 100% and 0% of safely evacuated citizens. Total percentage of evacuees saved for different level of supply condition is shown in Figure 4. Illustration of social injustice index is shown in Figure 5. From Figure 4, it is evident that is positive with the increment of bus fleet and negative with increase of demand level. The social injustice index is stable and near to 0 even under situations where all evacuees are not able to escape the location in time according to Figure 5. Actually, the same pattern is also observed in Scenarios B and C.
6.4. Discussion
The applicability of the general used objective function for busbased evacuation that minimizes travel cost of total evacuation network (MINEC), minimizes the duration of the evacuation (MINED), maximizes the total number of rescued citizens (MAXEN), and maximizes evacuation rate (MAXER), i.e., maximizes the number of evacuees who reach safety by any given deadline after the evacuation, is discussed in this part. The comparison of RNBSP, MAXEN, and MAXER is also discussed since similar evacuation scenarios are faced by EMA under these contexts.
6.4.1. Applicability Analysis of Different Objective Functions in Bus Evacuation Models
The objective function of the aforementioned bus evacuation models of MINEC, MINED, MAXEN, and MAXER is shown in following formulation, respectively:
According to output by the model proposed in Section 5, value of objective function can be obtained accordingly with (36), (37), (38), and (39) under RNBSP context.
Though seldom considered in evacuation modeling, the disaster characteristics affect the level of evacuation risk spatially and temporally as the disaster impact spread through the affected area in terms of different value of in this paper. With consideration of tempospatial difference of pickup locations caused by disaster characteristics, the applicability of the aforementioned four objective functions for busbased evacuation model is various due to different value of . For deterministic within the endangered area, if is long enough compared with the network clearance time, i.e., for , , it can be solved as a vehicle routing problem (VRP), represented by model of MINEC. The cost of evacuation operation is the main consideration for bus schedule as no casualties needs to worry about. If is smaller or slightly larger than network clearance time, it is a bus evacuation problem (BEP) represented by model of MINED. Whether the transitdependent citizens can be safely evacuated before the impacts of disaster can be determined by comparing the value of and . This is also the criterion of whether extra bus resources are requirement. If the bus fleet cannot complete evacuation before the impacts of disaster, in other words, the deadline is so short that not all of the population can be evacuated using the current available bus fleet; generally, the objective is either to maximize the evacuation rate represented by MAXER or to maximize the total number of citizens saved by available bus fleet represented by MAXEN. Normally, MAXEN is a special situation of MAXER when , where is the deadline of entire evacuation region. After introducing dynamic supply condition into the problem, the information provision of supply condition will distinguish these two objective functions by bus assignment. Both MAXER and MAXEN maximize the number of evacuees who leave pickup locations in different time horizon without identifying the level of risk endured by evacuees at different pickup locations.
6.4.2. RNBSP vs MAXER & MAXEN
Based on the aforementioned discussion, the evacuation scenario faced by EMA in RNBSP is similar to that in MAXER or MAXEN as the supply might not be sufficient to evacuate all the demands within the endangered area. Let be total number of saved citizens. Let denote the bus operating from to for m times. The comparison of RNBSP, MAXER, and MAXEN is introduced by the following three examples.
Example 1. An evacuation scenario is given in Figure 6(a). Only one bus is available in the area with capacity of 50. In RNBSP context, the trip order of the bus is , which is the same as that in MAXEN context. , =87.5% and =0.125 can obtained in both contexts. Under MAXER context, the trip order is , and is only 250, =62.5% and =0.375.
(a)
(b)
Example 2. Same evacuation scenario with Example 1 is given in Example 2, only the capacity of bus is reduced to 20. The trip order in RNBSP context is still , while, in MAXEN and MAXER context, the trip exists only between and for entire evacuation time horizon. i.e., . We can obtain =140, =35%, and =0.05 in RNBSP context and =200, =50%, and =0.5 in MAXEN and MAXER context. It is noteworthy that achieves theoretical maximum under this scenario. Although the number of saved citizens in MAXER context is higher than that in RNBSP, the evacuees in are ignored by models of MAXEN and MAXER. The problem of social equity is exposed under this bus schedule scheme. The evacuees in different pickup locations should have equal opportunity to be rescued in terms of human nature. Therefore, the limited resource should be sent to locations that are most in need.
Example 3. For evacuation scenarios in Figure 6(b) with bus capacity of 40, in RNBSP context, the trip order of the bus is , the trip order in MAXEN context is . It is interesting to find that the bus route in MAXER is various as assignment to or can provide the same amount evacuation rate. Here we find the route that can also provide maximum of total rescued population on premise of maximizing evacuation rate; it turns out that the trip order is the same as that in MAXEN. =320, =53.3% and =0.305 in RNBSP context, and =400, =66.7% and =0.471 in MAXEN and MAXER context. Actually, for the most adverse condition in MAXER context, the trip order is , we can obtain =200, =33.3% and =0.471. is ignored in context of MAXEN and MAXER; further, is also ignored in context of MAXER for the most adverse condition.
Consider a scenario in which the whole population of a city must be evacuated to shelters outside of the city because of a disaster expanding from the downtown area. MAXER suggests to evacuate people who live in the suburb since they are closer to the shelters without considering the fact that people in downtown may be in a greater danger based on Aalami and Kattan [1], while, under MAXEN context, the buses available will also be allocated to pickup locations near the city center with higher danger as long as the time left is enough to evacuate the citizens in the suburban area as discussed in Example 1. While if both pickup locations are suffering from risk of insufficient evacuation time, as discussed in Example 2, the limited resource will be allocated to pickup locations with highest evacuation rate in MAXEN context. The problem of social equity is exposed under bus schedule scheme both in MAXER and in MAXEN context. The evacuees in different pickup locations should have equal opportunity to be rescued. Therefore, the limited resource should be sent to locations with higher evacuation risk. We can conclude from the three examples that, under situation of pickup locations with different deadlines, the MAXER is not an appropriate objective function for busbased evacuation under evacuation scenarios with short and various deadlines. Further comparison between MAXEN and RNBSP is given in next part.
6.4.3. Influence of Information Provision on Bus Schedule in MAXEN Context
In this part, the bus schedule scheme of MAXEN is discussed under different level of information provision. It is noteworthy that, in any case, can only influence the bus schedule in and after under MAXER or RNBSP context, which means the computation complexity of MAXER and RNBSP is not sensitive to the length of evacuation time. Three scenarios of information provision are listed as follows: (a) complete information scenario, i.e., for , is known at the starting point ; (b) bounded information scenario, i.e., is known a bounded time space before t; (c) none information scenario, i.e., is informed at t. Let be the time point that bus b is available for assignment. Again, take the scenario in Figure 6(a) as an example; let , and , . The trip order and for the aforementioned three information scenarios in MAXER context combined with RNBSP context are, respectively, shown in Figure 7. In MAXEN context, under complete information scenario, =65%, =0.35, under bounded information scenario with 4 time unit advance, =65%, =0.35, under none information scenario, =60%, =0.4. In RNBSP context, =60%, =0.2. Comparing MAXER with RNBSP under different information scenarios in terms of and , we can conclude that the complete information scenario has a largest , followed by bounded information scenario. Although ignoring the importance of social equity, the advantage of MAXEN under none information scenario is no better than RNBSP even in terms of . We can conclude that information provision of supply condition has a critical influence on bus schedule in MAXEN context because it aims to maximize the number of safely evacuated citizens for entire time horizon instead of a certain time stage. The percentage of rescued population MAXERN is negatively related to degree of delay. However, the information has no influence on RNBSP as it is determined by current bus fleet size. Therefore, the robustness of proposed model is better compared with MAXEN under uncertain supply condition.
(a)
(b)
(c)
(d)
6.4.4. Solution Method of MAXEN
A solution method of MAXEN with consideration of different residual evacuation time among pickup locations is developed in this part, which can fast identify the bus schedule and estimate total number of evacuees that can be rescued.
The MAXEN model can be solved by dynamic programing approach of back reasoning. First, the pickup locations are sorted by its residual evacuation time at the starting point of evacuation. . The total evacuation time horizon of is divided into a series of stages with time intervals equal to the difference of two adjacent elements in . Similar to finding the shortest path, here, the “shortest path” is equivalent to the “path” that possesses highest value of rescued population, the “cost” of link is expressed by , denoting the number of rescued evacuees during . The approach searches backward from the last time period to the start of evacuation. If we seek the maximum value of , which equals as the total number of rescued evacuees at time . By back reasoning method used in shortest path identification, the following equations are obtained:where is evacuation demand of , and denotes the demand left in pickup locations with residual evacuation time larger than . The illustration of this approach is shown in Figure 8 with a given example, where denotes the evacuation efficiency of , interpreted as number of evacuees being rescued per unit time, equivalent to .
6.4.5. Summaries
This study quantifies evacuation risk in terms of number of casualties caused by impacts of disaster based on current bus schedule scheme. Modeled as a timedependent variable, it is treated as a proxy of social equity and characterized by disaster dynamics, demand, and supply conditions. Since the objective of RNBSP is to shave the peak of evacuation risk among different pickup locations in the network, the difference of casualties might be caused by disaster is lowest among pickup locations. Not surprisingly, the runtime for model with either objective grows quickly with problem size, especially with the minmax objective. This is the reason that the problem size should be controlled in order to obtain feasible answer for realistic operation in dynamic manner. Out of consideration of social equity, bus schedule of RNBSP is recommended in evacuation operation. Actually, the disaster characteristics are not always predictable or just shortterm predictable in realistic situation; thus the resource allocation based on current situation is more reliable than global allocation which covers entire evacuation time horizon. Based on the discussion above, MAXER is not suitable for busbased evacuation with consideration of pickup locations with different deadlines, and the efficiency of MAXEN is highly dependent on degree of advance and accuracy of predictable information of supply conditions.
The benefits of the riskbased dynamic bus schedule scheme are evident when considering social equity. The formulation of evacuation risk enables the DEU concept to be seamlessly applied to bus schedule scheme considering pickup locations with different residual evacuation time; further, the estimation of evacuation risk is applicable for different types of disasters, providing a generalized framework for busbased evacuation with dynamic supply condition and disaster impacts. It is interesting to find possibility that the performance of MAXER and MAXEN is no better than RNBSP even in aspect of number of rescued citizens. The results also highlight the importance of considering dynamics of disaster, demand, and supply conditions.
7. Conclusions
This paper presented a methodology of ability to model a riskbased dynamic bus schedule scheme for evacuation of transitdependent citizens under nonotice scenarios. The proposed concept of evacuation risk can reflect the dynamic situation of a pickup location within the endangered area considering disaster characteristics, demand patterns, and supply conditions. The concept of DEU introduced in this paper is a novel way for EMAs to design bus schedule scheme in a dynamic manner. Incorporation of disaster dynamics into evacuation modeling can better adapt the model to realistic situation.
The results of this research are of value to emergency management agencies. EMA can use the models proposed in this paper both in process of planning for a disaster and during disaster response. In the planning phase for deterministic disaster, the initial DEU determination model can be applied to identify the number of buses needed to evacuate all transitdependent citizens in a static manner. In the operations phase of nonotice evacuation, the RNBSP model has the ability to identify optimal schedule scheme to minimize the evacuation risk of transitdependent evacuees and maximize the vehicle utilization with limit resource in a dynamic manner. A future study direction is to take psychological experience into evacuation planning for transitdependent populations. Under nonotice evacuation, normally, evacuees will care less about their physical experience during evacuation operation, but the psychological experience will have longterm influence on them. On premise of ensuring evacuation safety, seeking methods to improve psychological experience of transitdependent evacuees during nonotice evacuation is an efficient way to take care of special needs and disadvantageous evacuees.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The opportunity to explore this topic was made possible by funding provided by the National Key Research and Development Program of China (Grant no. 2016YFC0802208), the National Natural Science Foundation of China (Grant no. 71101185), and the Natural Science Foundation of Shaanxi Province (Grant no. 2017JQ5122).