Multiobjective Optimization on Hierarchical Refugee Evacuation and Resource Allocation for Disaster Management

+is paper studies a location-allocation problem to determine the selection of emergency shelters, medical centers, and distribution centers after the disaster. +e evacuation of refugees and allocation of relief resources are also considered. A mixedinteger nonlinear multiobjective programming model is proposed to characterize the problem. +e hierarchical demand of different refugees and the limitations of relief resources are considered in the model. We employ a combination of the simulated annealing (SA) algorithm and the particle swarm optimization (PSO) algorithm method to solve the complex model. To optimize the result of our proposed algorithm, we absorb the group search, crossover, and mutation operator of GA into SA. We conduct a case study in a district of Beijing in China to validate the proposed methodology. Some computational experiments are conducted to analyze the impact of different factors, such as the target weight setting, selection of candidate shelters, and quantity of relief resources.


Introduction
A natural disaster is defined as a natural occurrence or activity that may cause injuries and deaths, loss of property and normal life, public disorder, economic recession, or environmental damage. e growing of population has escalated both the severity and frequency of natural disasters [1]. Research shows that natural disasters impact people with increasing frequency [2]. Site selection of emergency shelters (ESs) and emergency resource allocation have been studied in many disaster risk reduction projects, as they can provide safe havens, medical care, and life-saving resources for refugees [3][4][5].
In recent years, three different categories of evacuation/ refuge ESs are widely used after a natural disaster, including immediate shelters, short-term shelters, and long-term shelters. Immediate shelter is used for less than 1 day temporary evacuations, short-term shelter for less than 1 week duration of stay, and long-term shelter for more than 1 week and less than 1 month of stay. is hierarchical design of ESs has been widely used by many research studies [5,6]. However, the urgency for evacuation and resource demand of refugees in the same period are often different, and refugees should be divided into different priorities [7]. us, ESs should be divided into different categories in correspondence with different priorities of refugees in the same period. Furthermore, relief resources are typically not sufficient to meet all demands in the early days of disaster [8]. Hence, the goals of the emergency relief resource allocation planning policy should not only focus on efficiency (e.g., minimizing the evacuation distance and total cost), but also consider the fairness with the limitations of relief resources.
In this paper, we consider the facility site selection and relief resource allocation in the response phase of a disaster, especially during the critical 72-h time window after the disaster happens. We propose a hierarchical multiobjective model to determine site selection and the corresponding relief resource allocation. We consider not only traditional objectives, e.g., minimizing the hierarchical evacuation distance, minimizing the total costs (including construction costs of ESs and emergency distribution centers (EDCs), and fixed and transportation costs of emergency resources), and minimizing the total hierarchical unmet resources, but also maximizing the fairness of refugees. A combination of the SA algorithm and the PSO algorithm (SA-PSO) is applied to solve the complex model. e group search, crossover, and mutation operators of GA are combined into SA to optimize and expand the searching space of our proposed algorithm.
rough the real-case computational studies, we find that (1) compared to SA-PSO, the improved SA algorithm and the PSO algorithm (ISA-PSO) can optimize the results by expanding its searching space, but population search also increases time consumption; (2) when the objective weight increases, the objective value of the total evacuation distance improves the most; (3) increasing the resources in EDCs can greatly optimize the target of total hierarchical unmet resources with a small increase in the cost; (4) when the evacuation demand of refugees is increased, it is necessary to select more candidate shelters; (5) compared to the results of time-varying demand estimation, emergency evacuation planning under stable demand would cause more evacuation distance and more unmet resources. e main contributions of our research are as follows. First, we formulate a hierarchical multiobjective locationallocation model to determine the selection of hierarchical ESs, emergency medical centers (EMCs), and EDCs, along with the corresponding hierarchical evacuation and relief resource allocation. Second, we proposed a heuristic algorithm to solve the complex nonlinear model. Furthermore, we conduct a case study to validate our model and analyze the impacts of different algorithms and factors. ese analyses provide useful suggestions for the government on shelters and resource planning management. is paper is organized as the following. Section 2 provides review of the related literature, while Section 3 describes the problem and presents a multiobjective nonlinear programming model. en, Section 4 proposes two heuristic algorithms to solve this complex problem. We present computational results of the Chaoyang district, Beijing, China, as our case study in Section 5 and conclude in Section 6.

Literature Review
Emergency disaster management including site selection for emergency facilities and the corresponding resource allocation has been extensively studied in the literature. e site selection of ESs and EMCs has long been a critical and difficult issue. Nowadays, many models have been developed, including the minimal distance model like the P-center problem and the P-median problem and the covering model like the location set covering problem (LSCP) and the maximal coverage location problem (MCLP). Solution of the P-center problem and the P-median problem seeks to set the total number of selected ESs in advance. e objective of the P-center problem is to minimize the longest distance, while the objective of the P-median problem is to minimize the total travel distance from the AAs to the ESs [9,10]. e LSCP and the MCLP searches for the optimal location of ESs to cover all demands within a certain service distance. e LSCP is to minimize the total number of selected facilities as well as to minimize the construction cost, and the MCLP is to minimize the number of refugees that will not be served within the service distance [11]. However, in an actual response period, the project manager must deal with multiple conflicting objectives in the varying environment [12]. ese single-objective models are too simplistic to solve the ES and EMC site selection problems because of the ignorance of important objectives [13]. To solve the complex problem, multiobjective models based on the minimal distance model and the covering model have been developed by many studies. Tsai and Yeh [14] integrated the MCLP with the P-center problem and developed a double-target (minimizing the farthest distance and maximizing the total number of refugees in these candidate ESs) model to decide the site selection problem. Hu et al. [15] considered multiple criteria like the safety of victims, minimization of total evacuation distance, cost reduction, and capacity constraint satisfaction to solve the ES selection problem. Herein, site selection of ESs should consider multiple criteria such as minimizing the total distance/time, minimizing the total cost, and maximizing the coverage.
In addition, varying environment corresponds to different levels of service need [5]. Chen et al. [6] developed a three-level hierarchical ES location model based on the varying need of different periods after a disaster, with the objective of minimizing the total hierarchical distance. Previous classification of ESs is based only on the timevarying demand of different periods; however, no study has researched on the varying need of different levels of refugees to different levels of ESs in the same period. And, among these classified earthquake evacuation ESs, according to Xu et al. [16], the immediate shelter is not necessary to determine a specific model for its site selection because it is only used for an immediate evacuation. And, because the longterm shelter is used for a longer period, there are limited candidate locations of these shelters and are normally away from the city center. us, the location of a long-term shelter is not a complex problem. But, the location of a short-term shelter must meet multiple criteria and constraints as previously introduced. Since the service demand of refugees varies in the same period, the most critical site selection of the short-term shelter should consider hierarchical division based on different demand of refugees.
In the early stage of emergency disaster response, site selection of ESs and allocation of refugees are of great importance to emergency management. However, effective distribution of relief resources is also an indispensable part of this life-saving project [7]. In recent years, some studies focused on simultaneous site selection as well as relief resource allocation during the response phase. A three-stage programming model is designed by Sheu and Pan [17] for disaster management, attempting to identify the integration of the ES location network and the EMC location network. e objective functions of each stage involve not only the most fundamental criteria such as minimizing the total distance and cost, but also consider the psychological cost perceived by refugees in AAs, ESs, and EMCs. Rodríguez-Espíndola and Gaytán [18] developed a multicriteria approach to solve the simultaneous site selection of ESs and EDCs, along with the allocation of human resources and goods by minimizing the total costs and total distance. Dalal andÜster [19] proposed an overall disaster response model including supply (relief ) and demand (evacuation) with the trade-offs between the biobjective of evacuation time and system cost, considering uncertainty in refugee demand. ey presented a set of disaster scenarios to optimize the average and worst case costs. Furthermore, the demand for emergency resources is uncertain in emergency resource allocation problems [4]. Bozorgi-Amiri and Khorsi [8] pointed out that the relief resources are often not enough to meet the demand and presented a relief logistic planning model making the trade-off between three objectives. Hence, considering the insufficient feature of relief resources and the hierarchical demand of different refugees, the objective of this project should contain not only multiple criteria like cost, distance, and resource satisfaction, but also maximize the fairness of refugees. Fairness can be measured by the hierarchical demand for evacuation and relief resources.
Our model differs from most existing studies [4,7,16,17] in two aspects. e first is the formulation of a decision-making project that contemplates a novel model on the assignment of different priority refugees to different levels of ESs, using hierarchy to evaluate the urgent demand of refugees and improve efficiency of management. e other aspect is the allocating of resources considering insufficient relief resources and hierarchical requirement of each demand point.
In this paper, we consider the 72-hours gold rescue time after the disaster. Based on these previous research studies, we propose a hierarchical multiobjective nonlinear model, which considers different demands of different levels of refugees, simultaneous site selection of two-level ESs, EMCs, and EDCs, and the allocation of different level refugees and insufficient relief resources. Our study is distinguished from others by bringing a hierarchical design of refugees and ESs to optimize evacuation efficiency, using different hierarchy of refugees and ESs to measure their urgency for evacuation and resource demand. Furthermore, we use hierarchical management to allocate insufficient resources to different priority refugees in AAs, two-level ESs, and EMCs to achieve fairness. In this research, we use GIS to generate distance data to test the effectiveness and efficiency of this model in a case study. Finally, we use ArcGIS and visualize the comparison of target values to examine the trade-offs among three targets and parametric settings.

Problem Definition and Model Formulation
In this section, we provide the problem definition and model formulation on the emergency evacuation and the resource allocation network. e model considers the specific evacuation and the distribution network problem with hierarchical ESs, EMCs, and AAs ( Figure 1) in an advanced or a developed city because such a city has basic buildings like public parks that can provide short-term emergency evacuation and adequate hospitals for injuries.
In the evacuation and distribution network, we define the AAs as the source nodes whose populations dictate the required outflow. Another set of nodes of EMCs, A-ESs (Apriority emergency shelters) and B-ESs (B-priority emergency shelters), are the destination nodes. EMCs provide services for injuries, and A-ESs provide higher priority services than B-ESs for refugees. A final set of nodes of EDCs represent the source nodes of resources which provide resources to demand points.

Assumptions.
e proposed model is based on the following assumptions: Assumption 1. e time-varying demand of refugees and resources is given. Based on the basic population data of AAs, the number of different priority refugees and injuries is roughly estimated. Besides, the quantity of required resources can also be estimated from the number of people in each demand point. Assumption 2. Refugees are classified as nonevacuated, A-priority evacuated people, B-priority evacuated people, and injured people, and their corresponding places are AAs, A-ESs, B-ESs, and EMCs. We set the order of priority from high to low as injured, A-priority evacuated, B-priority evacuated, and nonevacuated. Assumption 3.
e resources are classified into two categories: commodities and medical supplies. Injuries and different hierarchy of refugees require different kinds of resources. e refugees with higher priority can receive more and closer resources.   4. e EMCs are selected from the existing hospitals. erefore, regardless of the construction cost of the EMCs, all EMCs can be the selected points. e cost involves the construction cost of A-ESs, B-ESs, and EDCs and the fixed and transportation cost of resources. Assumption 5. e number of evacuation and transportation vehicles is unlimited. erefore, all refugees, injured, and resources can be transported to the corresponding nodes. E a : the total number of people in the affected area a, a ∈ AA E t a : the number of refugees in the affected area a at time t, a ∈ AA, t ∈ T E t Aa : the number of A-priority refugees in the affected area a at time t, a ∈ AA, t ∈ T E t Ba : the number of B-priority refugees in the affected area a at time t, a ∈ AA, t ∈ T E t Ia : the number of injured people in the affected area a at time t, a ∈ AA, t ∈ T R rc : the quantity of relief resources r in the emergency distribution center c, r ∈ R MA ψ : the maximum capacity (unit: person) of node ψ,

Model
CR r : the fixed cost (unit: RMB) of per unit resources r, r ∈ R CT r : the transportation cost (unit: RMB/km) of per unit resources r, r ∈ R D aψ : the distance (unit: km) from the affected area a to node ψ, ψ ∈ A − ES ∪ B-ES ∪ EMC D cΦ : the distance (unit: km) from the emergency distribution center c to node Φ, the ratio of total resources in the emergency distribution center to the total required resources (considering the insufficient feature)

Variables. Dependent variables and decision variables
in this model are given as follows: RT r : the total quantity of resources r required by all resource demand points throughout the rescue period of 72h, r ∈ R R t urΦ : the quantity of unmet resources r in node Φ at time t, r ∈ R, Φ ∈ AA ∪ A-ES ∪ B-ES ∪ EMC, t ∈ T FR t urΦ : the loss for the unmet resources of node Φ at time t, Decision variables: aj : the number of A-priority refugees transported from affected area a to A-priority emergency shelter j at time t, a ∈ AA, j ∈ A − ES, t ∈ T EB t ak : the number of B-priority refugees transported from affected area a to A-priority emergency shelter k at time t, a ∈ AA, k ∈ B − ES, t ∈ T EH t ah : the number of injured people transported from affected area a to emergency medical center h at time t, a ∈ AA, h ∈ EMC, t ∈ T R t rcΦ : the quantity of resources r transported from emergency distribution center c to node Φ at time t,

Formulation.
e proposed hierarchical multiobjective nonlinear mathematical formulation is presented as follows: Objective function: is subject to e programming model involves three objective functions (minimization of the total evacuation distance, total cost, and total unmet resources), specified by equations (2)-(4).

Mathematical Problems in Engineering 5
Constraints (5)-(7) make sure that all refugees requiring hierarchical evacuation in AAs can be evacuated into the corresponding nodes. e maximum capacities of A-ESs, B-ESs, EMCs, and EDCs are limited in constraints (8)- (11). Constraints (12)- (19) reflect the number of people in each node at time t. Constraints (20)-(23) define the quantity of required resources in each node at time t; thus, constraint (24) defines the total quantity of relief resources r required from all resource demand points throughout the rescue period of 72 h. To avoid the wastage of insufficient resources, we define that the quantity of received resources of each node is no more than the quantity of required resources of each node by equation (25). Considering the hierarchical resource demand of each node and insufficient resources of EDCs, to reduce the secondary damage caused by insufficient resources, constraint (26) sets the total quantity of resources in all EDCs (α ≥ 1 means the total quantity of resources is sufficient). Constraint (27) indicates the quantity of unmet resources in each node. Constraint (28) defines the loss for the unmet resources of each node at time t, where the parameter μ is used to show the importance placed on large losses.

Algorithm Frameworks.
In the previous section, we proposed the hierarchical multiobjective nonlinear mathematical model. Considering our research on site selection and resource allocation during the critical 72 h after disaster, the general commercial solver like CPLEX and LINGO cannot solve this complex nonlinear problem [5]. e SA algorithm is the expansion of the local search algorithm. Different from the local search algorithm, SA accepts a worse solution with a certain probability in iterations to prevent local optimum [20]. Based on this optimization feature, SA has been used in the site selection research [21], manufacturing problem [22], and multiobjective problem [23]. SA can optimize to the global optimum with very slow annealing [24]. erefore, the SA algorithm can be a good solution to solve the site selection of A-ESs, B-ESs, EMCs, and EDCs. e PSO algorithm is a stochastic evolutionary based on birds' group behavior. Due to its fast convergence, high efficiency, and robustness, PSO has become a popular tool to solve complex problems in the social network [25] and the computer science [26]. e algorithm has also been used in solving largescale emergency resource allocation problems and has proved to be time efficient [4]. In our research, PSO is an effective algorithm to solve the resource allocation problem.
In this section, we use the combination of two heuristic algorithms to solve our programming model, and outer SA and inner PSO methods are applied to site selection and resource allocation. e flow of our combination algorithms is shown in Figure 2. Next, we would explain how the SA and PSO methods are applied to our multiobjective problem with constraints.

Fundamental SA-PSO Algorithm
e main steps of the SA algorithm are as follows: Step 1. Initialize the parameters: start temperature t s , end temperature t e , decay rate of temperature ρ, receive tolerance for inferior solutions ε, initial solution s 0 , optimal solution unchanged times key, and maximum optimal solution unchanged times keymax.
Step 2. Iterate L times in each temperature, generate a new neighbor s ′ from the previous solution s based on certain rules in the current temperature t k , according to the Metropolis guidelines [27], and decide if s ′ is accepted as the current solution with the following probability: where t k is the current temperature of the kth iteration, t k > 0, ∀k, and f(s) is the target value of solution s.
Step 3. Outer loop: if t k is lower than t e or key ≥ keymax, end the algorithm; else, if t k � ρ * t k , return to step 2.
In our combination algorithms of solving site selection and resource allocation problems, we get the optimal objective value through the outer SA algorithm.
us, the objective function f(s) of the SA algorithm is based on all objective functions f 1 , f 2 , and f 3 . By combining the multiobjective into a single-objective function, we use a weighted sum approach that sets the weights of different targets. erefore, f(s) can be represented as the following equation, and 3 ∇�1 ω ∇ � 1: Minimum and maximum objective function values are shown in Table 1.
e random initialization site selection of the A-ES, B-ES, and EMC solution is based on constraints (5)- (10), which ensure that all refugees are evacuated into corresponding nodes with limited capacity. And, the random initialization site selection of the EDC solution is based on constraint (26). e initial site selection solution s 0 of the SA algorithm is shown in Figure 3, where x j � 1, x k � 1, and x c � 1 mean the candidates A-ES, B-ES, and EMC are selected and x j � 0, x k � 0, and x c � 0 mean the candidates A-ES, B-ES, and EMC are not selected. e SA algorithm randomly selects a new solution from the neighborhood feasible solution space N(s) of the current solution s. In our site selection problem, the method of generating a new solution s ′ , s ′ ∈ N(s), can be divided into three types.
Randomly adding an unselected point; randomly removing a selected point on the premise of meeting the above constraints; randomly changing the selected point (randomly adding an unselected point while randomly removing a selected point) on the premise of meeting the above constraints. 6 Mathematical Problems in Engineering Taking the location selection of the four candidate points as an example, assume that the initial solution s is as shown in Figure 4, where 1 means the candidate point is selected, 0 otherwise. en, based on the previous method of generating a new solution, the neighborhood feasible solution space N(s) is shown in Figure 5.

PSO Algorithm.
e main steps of the PSO algorithm are as follows: Step 1. Initialize the parameters: set N as the particle number, set ps 0 i as the initial position of the ith particle (i ∈ [1, N]), set v 0 i as the initial velocity of the ith particle (i ∈ [1, N]), and calculate the fitness function  f(ps 0 i ) of the initial position of the ith particle (i ∈ [1, N]).
Step 2. Updating the velocity and position of the ith particle in the kth iteration with the following equations: Step 3. Calculating the fitness function f(ps k i ) of the ith particle in the kth iteration (k ∈ [1, K]), personal best pbest k i and global best gbest k are obtained, and if k > K, return gbest k and end the algorithm, else return to step 2. e inner PSO algorithm is to find the best solution to allocate insufficient resources, and the fitness function f(ps k i ) only optimizes the resource allocation part. us, based on the objective function of f 2 and f 3 , f(ps k i ) can be represented as follows: e random initialization resource allocation solution from EDCs to each resource demand points at time t is based on constraints (11) and (25). e initial PSO resource allocation solution of the ith particle (i ∈ [1, N]) ps 0 i can be represented by a four-dimensional array of all R t rch , R t rcj , R t rck , and R t rca in the model in Section 3.2. e updating method of PSO is based on equations (31) and (32). If ps k i does not satisfy constraints (11) and (25), R t rc * ( * ∈ EMCs, A-ESs, B-ESs, and AAs) will be set to 0.

Improved SA-PSO Algorithm.
e SA algorithm has the advantages of jumping out of local optimal and fast convergence speed. But, its single individual evolution method cannot grasp the entire search space, which makes its search efficiency not high [20,28]. GA is a global search algorithm based on the theory of evolution and the natural rule of "survival of the fittest" [29]. GA has been widely applied in real-world problems as it is very robust and has strong global search capabilities [30]. Based on the characteristics of SA and GA, we absorb the group search, crossover, and mutation operator of GA into SA to optimize and expand its population diversity. us, the outer site selection problem is solved by the ISA algorithm.

GA Algorithm.
e main steps of the GA algorithm are as follows: Step 1. Initialize the parameters: population size G, maximal number of generations M, crossover rate p c , and mutation rate p m , set s 0 g as initial solution of the gth individual (g ∈ [1, G]), calculate the fitness function f(s 0 g ) of the gth individual, and find Fitness Best .
Step 2. Use roulette wheel selection to select the new population, and use the randomly crossover and mutation to optimize and widen the solution space.

ISA Algorithm.
e main steps of the GA algorithm are as follows: Step 1. Initialize the parameters of SA: start temperature t s , end temperature t e , decay rate of temperature ρ, receive tolerance for inferior solutions ε, optimal solution unchanged times key, and maximum optimal solution unchanged times keymax; initialize the parameters of GA: population size G, crossover rate p c , and mutation rate p m . Set Step 3. Outer loop: if t k is lower than t e or key ≥ keymax, end the algorithm; else if t k � ρ * t k , return to step 2.
Based on the initial site selection and the update solution of the SA algorithm in Section 4.2.1, the chromosome, which represents the initial site selection solution of the ISA algorithm, is the same as s 0 . Combining the characteristics of GA evolution, the update solution of ISA includes selection, crossover, and mutation operation.

Case Study
In this section, we intend to (1) present a case study to illustrate the efficiency of our mathematical model to realize quick location selection and resource allocation during disasters; (2) consider some scenario analyses for further investigation on emergency refugee evacuation and resource allocation decisions; (3) analyze the trade-off between different factors through the result of the computational experiments.
Beijing is located in the northern part of the North China Plain, with a latitude of 39°26 ′ N to 41°03 ′ N and a longitude of 115°25 ′ E to 117°30 ′ E, mainly in the North China earthquake zone Yin-middle Yanshan seismic zone. e earthquakes with a magnitude greater than 4 in the Beijing area have occurred nearly 200 times, and more than 10 earthquakes with magnitude 5 have occurred.
ere is a possibility of a moderately devastating earthquake in the Beijing area, which makes the situation of earthquake prevention and mitigation faced by Beijing very serious. Chaoyang district, affiliated to Beijing, located in the eastern part of Beijing, covers an area of 470.8 km 2 and is the largest district in the Beijing central city. In recent years, some emergency shelters have been built in Beijing, in which the Chaoyang district has the largest emergency shelters. We choose the Chaoyang district ( Figure 6) as our study area to demonstrate the effectiveness and efficiency of our proposed model and algorithms.

Data Processing.
We divide the first 3-day gold rescue time after a disaster into 6 hours of time intervals. us, the emergency evacuation and resource allocation network are designed in 12 (�72/6) intervals.

Facility Data.
According to the database of the Beijing Earthquake Agency, the Chaoyang district has built fifteen earthquake disaster shelters. However, these shelters may not accommodate so many refugees or cause more evacuation distances.
us, based on the Beijing Central City Earthquake and Emergency Refuge Site Planning Outline and the Chaoyang Government Network, another sixteen candidate shelters were selected.
e EMCs need professional medical equipment and nursing staffs. us, in our research, we select 26 candidate EMCs from the hospitals located in the Chaoyang district, as shown in Figure 7. Since there are no EDCs in the Chaoyang District, Beijing, we  Mathematical Problems in Engineering assume that candidate EDCs are located at street/district offices. So, there are 42 candidate EDCs.

Population.
Chaoyang district governs 24 streets and 19 districts, and the population size of each street/district E a is extracted from the Chaoyang District Beijing Statistical Yearbook 2017. e number of refugees in each AA of each time interval E t a is roughly estimated [17], and Figure 8 presents the ratio of E t a to E a . We assume that 10% of the population is injured E t Ia , 20% of the population is assumed to be A-priority refugees E t Aa , and 30% is assumed to be B-priority refugees E t Ba .

Demand.
Quantity of required resources of different refugees is estimated by previous research [17]. Each refugee requires two types of relief resources: commodities and medical supplies. Commodities per unit are herein 1 kg, and the medical supplies per unit are 0.5 kg. Different refugees have different quantity of demand for different resources, and the specific resource requirements are presented in Table 2. Tables 3 and 4.

Capacity. We assume that each ES is divided into two parts: A-ES and B-ES, where the capacity of A-ESs is smaller than B-ES. e capacities of A-ES, B-ES, and EMCs are presented in
Because of the lack of EDC data in Chaoyang District, Beijing, combining the total demand of resources of all refugees in each node, we assume here that 9 EDCs were selected from 42 candidate EDCs based on previous EDC selection research studies [4,8,[31][32][33]. e capacity of each EDC can be derived from the total number of EDCs, total quantity of required resources, and the ratio α.

Distance.
We obtain the road networks from the map of Beijing published by the China Cartographic Publishing House in 2016. We calculate the actual distances from AAs to A-ESs, B-ESs, and EMCs and from EDCs to AAs, A-ESs, B-ESs, and EMCs using version 10.1 of ArcGIS.

Cost.
e total cost involves the construction cost of candidate A-ESs, B-ESs (existing shelters have no construction cost), and EDCs and the fixed and transportation cost of per unit resource r. e construction cost of each site is determined by its capacity, which is estimated by the government report and previous research studies [5,18,19]. e fixed cost of resources is also estimated by previous research studies, and the transportation cost of resources is estimated by the Beijing logistic price. All these parameters are shown in Table 5.

Model Testing and Results.
In this section, the proposed model and algorithms are tested with the abovementioned parameter settings. As the priority of refugees in different places is EMCs > A-ESs > B-ESs > AAs and * ω * � 1 ( * ∈ EMCs, A-ESs, B-ESs, and AAs), we set the priorities of different places ω h , ω j , ω k , and ω a as following. We consider that the quantity of resources in all EDCs is insufficient and higher priority places have a higher minimum resource satisfaction degree. Different satisfaction degree settings are based on the ratio of total resources in DCs to total required resources. α � 0.9 means that the total resources of the emergency distribution center account for 90% of the total resources needed for the entire rescue cycle. μ is a measure of the fairness loss function, and μ < 0. When the amount of unsatisfied resources increases, the fairness loss function grows like an exponential function as μ decreases and grows like a linear function as μ increases. For parameters of SA, ω 1 , ω 2 , and ω 3 are the weights of the three objective functions in the model, and also the weights of the three subobjective functions in the fitness function of the outer simulated annealing algorithm. t s is the start temperature, t e is the end temperature, ρ is the decay rate of temperature, ε is the received tolerance for inferior solutions, key is the optimal solution unchanged times, keymax is the maximum optimal solution unchanged times, and L is the iterate times in each temperature. For parameters of PSO, K is the maximum iteration times, N is the particle number, c 1 and c 2 are the learning factors, c 1 indicates the possibility of the particle updating according to its local optimal position, c 2 indicates the possibility of the particle updating according to the global optimal position, w is the inertia factor, which reflects the particle's ability to inherit the previous speed, and ω 2 ′ and ω 3 ′ are the weight setting of the objective functions f 2 ′ and f 3 . For parameters of GA, G represents the number of groups, p c is the crossover rate, and p m is the mutation rate. According to relevant research [20-22, 25, 26, 29, 30], parameters of the model and algorithms are also shown in Table 6.
By setting the weight of each objective function equal to 1, with the objective of minimizing/maximizing f(s) of SA, we calculate the minimum/maximum value of each objective of our two test cases (as shown in Table 1).

Existing B-ESs Candidate A-ESs
Candidate B-ESs Urban road We firstly test the case with these existing shelters as case 1, in which the construction cost of A-ESs and B-ESs set to be zero. en, these candidate shelters are added in the second test case as case 2, using the parameters shown in Table 6.
rough 50 experiments, the average value of objective functions of the two test cases is shown in Table 7. e dividend is calculated by the equation (� (f(case 1)/f(case 2)) − 1), as shown in Table 7, and the results reveal that the total evacuation distance in case 2 decreases by 40.92%, while the total cost increases by 3.85% with selected candidate shelters compared to case 1. And, we see that the total unmet resources in case 2 also decreases.   us, it is obvious that selecting candidate shelters can reduce the total evacuation distance to a large extent and also reduce the total unmet resources while only adding a small amount of cost.
Taking the corresponding refugee evacuation in which the objective function value is the closest to the average objective function value as an example, the allocation results of refugees from AAs to A-ESs and B-ESs at time t � 1 of the two cases are shown in Figure 9.
One can see from Figure 9(a) that many AAs have no nearby shelter, leading to the relatively long evacuation distance of refugees in these AAs. However, as is shown in Figure 9(b), six A-ESs and six B-ESs were selected at time t � 1 (one more A-ES and B-ES were selected at time t � 2) and optimized the long-distance allocation result.
rough the comparison of Figures 9(a) and 9(b), it is obvious that selecting candidate shelters can greatly reduce the evacuation distance of some AAs.

Scenario Analyses.
In this section, five scenarios with added candidate shelters are designed to compare the characteristics of our proposed model and change in disaster situation of refugees. As can be seen from the results of Section 5.2, selecting candidate shelters is a more efficient planning method. erefore, all the scenarios in this section are based on case 2.

Scenario 1.
e first scenario is to compare the efficiency of SA-PSO and ISA-PSO. Parameters of the GA are shown in Table 8. We analyze these two algorithms by comparing the results of case 2 under SA-PSO and ISA-PSO. e operational iteration graph of the two algorithms is shown in Figure 10, and through 50 experiments, the average value of three objective functions is presented in Table 8.
It can be seen from Figure 10 that the iterative result of ISA-PSO is always better than that of SA-PSO. Table 9 presents the objective values calculated by ISA-PSO. e dividend is calculated by the equation (� (f(scenario 1)/f(case 2)) − 1), and it can be seen that compared with SA-PSO, the objective function values by ISA-PSO can achieve Pareto optimization. But, because of the characteristics of population search of ISA-PSO, the time consumption of ISA-PSO is correspondingly longer. e algorithm used in the following calculations will use ISA-PSO. And, the following comparison of objective values will be based on the results of ISA-PSO.

Scenario 2.
e second scenario is to analyze the impact of the objective weight.
It can be seen from Table 9 that f 1 is greatly affected by the objective weight. When the weight increases, the objective value is optimized, while the objective value of f 2 and f 3 is affected by the weight to a small extent. In addition, it can be seen that f 1 and f 3 are negatively correlated with f 2 , and f 2 increases as f 1 and f 3 decrease. e corresponding dividend comparison of three objective values of Scenario 2 is shown in Figure 11. e dividend is calculated by the equation (� (f(S2 − * )/f(scenario 1)) − 1)( * ∈ 1, 2, 3), and it can be inferred that f 1 is the most sensitive to the objective weight. erefore, provided that the budget is sufficient, increasing the objective weight of f 1 appropriately can better optimize the result.

Scenario 3.
e third scenario is to analyze the effect of parameter α on three objective functions. e purpose of considering this scenario is to enable managers to improve the efficiency of resource allocation under uncertain or  insufficient conditions. In our research, we assume the quantity of relief resources is insufficient, and here we change the ratio α as S3-1: α � 0.75, S3-2: α � 0.85, and S3-3: α � 0.95, respectively. rough 50 experiments, the average value of three objective functions of Scenario 3 is presented in Table 10. Table 10 presents that the change in α has little effect on f 1 and f 2 . It implies that resources and their transportation costs account for only a small part of the total cost (f 2 ). And, obviously, as α increases, f 3 decreases. Figure 12 illustrates the corresponding dividend comparison of Scenario 2. e dividend is calculated by the equation (� (f(S3 − * )/f (scenario 1)) − 1)( * ∈ 1, 2, 3). It is observed that as α increases, f 3 decreases rapidly and there were almost no change in f 1 and f 2 . is scenario illustrates that adjusting the parameter α significantly changes objective values of f 3 in a linear relationship, but does not affect the other objects. is implies that increasing the quantity of resources is an effective method that improves the target value of f 3 . It can be seen that preparing sufficient          A-priority and B-priority refugees to be 30% and 40% of the population in S4-1 and 40% and 50% of the population in S4-2.
In this scenario, the total capacity of the existing shelters cannot meet the refugee demand. us, the selecting of candidate shelters is essential when the refugee demand increases.
rough 50 experiments, the average value of three objective functions of Scenario 4 is presented in Table 11. e dividend is calculated by the equation (� (f(s4 − * )/f(scenario 1)) − 1)( * ∈ 1, 2). It is obvious that when the number of A-and B-priority refugees increases, the total evacuation distance increases evidently.
Taking the corresponding refugee evacuation in which the objective function value is the closest to the average objective function value as an example, the allocation results of refugees from AAs to A-ESs and B-ESs at t � 1 of Scenario 4 are shown in Figure 13. Combining the results in Table 9, all these imply that the construction cost has only very small impact on the total cost f 2 . Besides, when the number of refugees who need to evacuate increases, the total evacuation distance will also increase even if more candidate shelters are selected. Furthermore, when the evacuation demand of Aand B-priority refugees increases, the number of refugees in A-ESs and B-ESs increases as well. is result demonstrates that refugees in A-ESs and B-ESs are more sensitive to unmet resources than those in AAs.

Scenario 5.
In this scenario, we set the population of each period to be stable and compare the results with timevarying demand and stable demand. Figure 14 displays the affected population in each time interval with a stable refuge demand. rough 50 calculations, the average value of three objective functions of Scenario 5 is presented in Table 12. e dividend is calculated by the equation (� (f(scenario 5)/f(scenario 1)) − 1). e results in Table 12 show that the total evacuation distance increment is 14.8% higher, while the total cost is only 0.54% less than the result under time-varying demand. e reason of the evacuation distance dividend is that the remaining capacity of A-ESs and B-ESs becomes smaller as time goes by, but the number of A-and B-priority refugees has not become smaller. erefore, refugees in the later period will have to evacuate to farther A-ESs and B-ESs.
Similarly, most of the resources in EDCs were allocated to the AAs in the early period. e remaining resources in EDCs are getting less and less as time goes by. However, the demand for resources of A-ESs and B-ESs is increasing much faster than the time-varying demand model. erefore, the resource requirements of A-ESs and B-ESs may not be met in the later period. Recall the implication of Scenario 4, and this might be the reason that the objective value of unmet resources is also 41.88% higher than the results of the timevarying demand. us, it is obvious that the time-varying demand estimation is a more efficient evacuation and allocation planning model because it reduces evacuation distance and unmet resources. Furthermore, it fits the actual refuge demand better.

Conclusions
is paper studies decision-making on the dynamic emergency shelter and resource allocation for disaster management, considering the hierarchical requirement of refugees and allocation of insufficient relief resources. e planning of evacuating hierarchical refugees to hierarchical shelters and EMCs along with the corresponding resource allocation to these refugees is a difficult problem with multicriteria restrictions, especially when the quantity of relief resources is insufficient.
us, a multiobjective location-allocation model of minimizing the total hierarchical evacuation distance, total cost, and total hierarchical unmet resources is proposed to characterize the problem. It considers both traditional objectives and fairness perceived by refugees. e SA algorithm was applied to site selection in the outer layer, and the PSO algorithm in the inner layer is to solve the   corresponding resource allocation. And, an improved SA-PSO is proposed in order to expand the local exploration scope of SA and enhance its population diversity. A case study of Chaoyang District, Beijing, is conducted to verify the capability of the proposed model and algorithms. e computational experiment results indicate that selecting candidate shelters for refugees is obviously more worthwhile and more effective than evacuating refugees to existing shelters. e results of this study can be considered as a shelter selection proposal for Chaoyang District, Beijing, emergency facility planning. Further, some scenario analyses are considered for further investigation on emergency refugee evacuation and resource allocation decisions. e scenario analyses not only provide useful suggestions on site selection of candidate emergency facilities and corresponding resource allocation, but also perceive a more costeffective and realistic estimation on the refugee demand.
For future research, the real-time road network like congestion and destruction should be considered in the refugee evacuation planning. us, the evacuation planning and resource allocation planning would be different over time, with the change in the road distance. In addition, integration with physicians and nurse allocation and multitype vehicle allocation could provide a more comprehensive disaster reduction system.

Data Availability
e numerical study data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that they have no conflicts of interest.