Research Article  Open Access
Jian Wang, Danqing Shen, Mingzhu Yu, "Multiobjective Optimization on Hierarchical Refugee Evacuation and Resource Allocation for Disaster Management", Mathematical Problems in Engineering, vol. 2020, Article ID 8395714, 18 pages, 2020. https://doi.org/10.1155/2020/8395714
Multiobjective Optimization on Hierarchical Refugee Evacuation and Resource Allocation for Disaster Management
Abstract
This paper studies a locationallocation problem to determine the selection of emergency shelters, medical centers, and distribution centers after the disaster. The evacuation of refugees and allocation of relief resources are also considered. A mixedinteger nonlinear multiobjective programming model is proposed to characterize the problem. The 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.
1. 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. The 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 lifesaving resources for refugees [3–5].
In recent years, three different categories of evacuation/refuge ESs are widely used after a natural disaster, including immediate shelters, shortterm shelters, and longterm shelters. Immediate shelter is used for less than 1 day temporary evacuations, shortterm shelter for less than 1 week duration of stay, and longterm shelter for more than 1 week and less than 1 month of stay. This 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]. Thus, 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 72h 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 (SAPSO) is applied to solve the complex model. The group search, crossover, and mutation operators of GA are combined into SA to optimize and expand the searching space of our proposed algorithm. Through the realcase computational studies, we find that (1) compared to SAPSO, the improved SA algorithm and the PSO algorithm (ISAPSO) 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 timevarying demand estimation, emergency evacuation planning under stable demand would cause more evacuation distance and more unmet resources.
The 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. These analyses provide useful suggestions for the government on shelters and resource planning management.
This 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. Then, 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.
2. Literature Review
Emergency disaster management including site selection for emergency facilities and the corresponding resource allocation has been extensively studied in the literature.
The 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 Pcenter problem and the Pmedian problem and the covering model like the location set covering problem (LSCP) and the maximal coverage location problem (MCLP). Solution of the Pcenter problem and the Pmedian problem seeks to set the total number of selected ESs in advance. The objective of the Pcenter problem is to minimize the longest distance, while the objective of the Pmedian problem is to minimize the total travel distance from the AAs to the ESs [9, 10]. The LSCP and the MCLP searches for the optimal location of ESs to cover all demands within a certain service distance. The 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]. These singleobjective 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 Pcenter problem and developed a doubletarget (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 threelevel 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. Thus, the location of a longterm shelter is not a complex problem. But, the location of a shortterm 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 shortterm 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 lifesaving project [7]. In recent years, some studies focused on simultaneous site selection as well as relief resource allocation during the response phase. A threestage 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. The 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íguezEspí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 tradeoffs between the biobjective of evacuation time and system cost, considering uncertainty in refugee demand. They 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]. BozorgiAmiri 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 tradeoff 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. The first is the formulation of a decisionmaking 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. The other aspect is the allocating of resources considering insufficient relief resources and hierarchical requirement of each demand point.
In this paper, we consider the 72hours 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 twolevel 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, twolevel 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 tradeoffs among three targets and parametric settings.
3. 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. The 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 shortterm 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, AESs (Apriority emergency shelters) and BESs (Bpriority emergency shelters), are the destination nodes. EMCs provide services for injuries, and AESs provide higher priority services than BESs for refugees. A final set of nodes of EDCs represent the source nodes of resources which provide resources to demand points.
3.1. Assumptions
The proposed model is based on the following assumptions: Assumption 1. The timevarying 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, Apriority evacuated people, Bpriority evacuated people, and injured people, and their corresponding places are AAs, AESs, BESs, and EMCs. We set the order of priority from high to low as injured, Apriority evacuated, Bpriority evacuated, and nonevacuated. Assumption 3. The resources are classified into two categories: commodities and medical supplies. Injuries and different hierarchy of refugees require different kinds of resources. The refugees with higher priority can receive more and closer resources. Assumption 4. The EMCs are selected from the existing hospitals. Therefore, regardless of the construction cost of the EMCs, all EMCs can be the selected points. The cost involves the construction cost of AESs, BESs, and EDCs and the fixed and transportation cost of resources. Assumption 5. The number of evacuation and transportation vehicles is unlimited. Therefore, all refugees, injured, and resources can be transported to the corresponding nodes.
3.2. Model Formulation
3.2.1. Sets and Parameters
Sets and parameters of this model, along with their brief description, are given as follows: Sets and indices: AA, : set and index of the affected areas, AES, : set and index of Apriority emergency shelter locations, BES, : set and index of Bpriority emergency shelter locations, EMC, : set and index of candidate existing emergency medical centers, EDC, : set and index of emergency distribution centers, , : set and index of the time period, , : set and index of the type of relief resources, as commonly indexed by ; refers to commodities; refers to medical supplies, Parameters: : the total number of people in the affected area , : the number of refugees in the affected area a at time , , : the number of Apriority refugees in the affected area a at time , , : the number of Bpriority refugees in the affected area a at time , , : the number of injured people in the affected area a at time , , : the quantity of relief resources r in the emergency distribution center , : the maximum capacity (unit: person) of node , BES EMC : the construction cost (unit: RMB) of node , BES EDC : the fixed cost (unit: RMB) of per unit resources , : the transportation cost (unit: RMB/km) of per unit resources , : the distance (unit: km) from the affected area a to node , BES EMC : the distance (unit: km) from the emergency distribution center c to node , AES BES EMC : the priority weight of node , AES BES EMC : the number of resources r consumed by per person per unit time in node , AES BES EMC, : the ratio of total resources in the emergency distribution center to the total required resources (considering the insufficient feature)
3.2.2. Variables
Dependent variables and decision variables in this model are given as follows: Dependent variables: : the number of priority refugees in Apriority emergency shelter at time , , : the number of priority refugees in Bpriority emergency shelter at time , , : the number of injured people in emergency medical center at time , , : the quantity of resources required by node at time , , AES BES EMC, : the total quantity of resources required by all resource demand points throughout the rescue period of 72h, : the quantity of unmet resources in node at time , , AES BES EMC, : the loss for the unmet resources of node at time , Decision variables: : 1 if node is selected, 0 otherwise, BES EDC. : the number of priority refugees transported from affected area to Apriority emergency shelter at time , , , : the number of priority refugees transported from affected area to Apriority emergency shelter at time , , , : the number of injured people transported from affected area to emergency medical center at time , , , : the quantity of resources transported from emergency distribution center to node at time , , , AES BES EMC,
3.2.3. Formulation
The proposed hierarchical multiobjective nonlinear mathematical formulation is presented as follows: Objective function: is subject to
The programming model involves three objective functions (minimization of the total evacuation distance, total cost, and total unmet resources), specified by equations (2)–(4). Constraints (5)–(7) make sure that all refugees requiring hierarchical evacuation in AAs can be evacuated into the corresponding nodes. The maximum capacities of AESs, BESs, 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 ( 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.
4. Solution Methodology
4.1. 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].
The 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]. Therefore, the SA algorithm can be a good solution to solve the site selection of AESs, BESs, EMCs, and EDCs.
The 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]. The 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. The 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.
4.2. Fundamental SAPSO Algorithm
4.2.1. SA Algorithm
The main steps of the SA algorithm are as follows: Step 1. Initialize the parameters: start temperature , end temperature , decay rate of temperature , receive tolerance for inferior solutions , initial solution , optimal solution unchanged times key, and maximum optimal solution unchanged times keymax. Step 2. Iterate times in each temperature, generate a new neighbor from the previous solution s based on certain rules in the current temperature , according to the Metropolis guidelines [27], and decide if is accepted as the current solution with the following probability: where is the current temperature of the th iteration, , , and is the target value of solution . Step 3. Outer loop: if is lower than or , end the algorithm; else, if , 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. Thus, the objective function of the SA algorithm is based on all objective functions , , and . By combining the multiobjective into a singleobjective function, we use a weighted sum approach that sets the weights of different targets. Therefore, can be represented as the following equation, and :
Minimum and maximum objective function values are shown in Table 1.

The random initialization site selection of the AES, BES, 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). The initial site selection solution of the SA algorithm is shown in Figure 3, where , , and mean the candidates AES, BES, and EMC are selected and , , and mean the candidates AES, BES, and EMC are not selected.
The SA algorithm randomly selects a new solution from the neighborhood feasible solution space of the current solution s. In our site selection problem, the method of generating a new solution , , 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.
Taking the location selection of the four candidate points as an example, assume that the initial solution is as shown in Figure 4, where 1 means the candidate point is selected, 0 otherwise.
Then, based on the previous method of generating a new solution, the neighborhood feasible solution space is shown in Figure 5.
4.2.2. PSO Algorithm
The main steps of the PSO algorithm are as follows: Step 1. Initialize the parameters: set as the particle number, set as the initial position of the th particle (), set as the initial velocity of the ith particle (), and calculate the fitness function of the initial position of the ith particle (). 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 of the ith particle in the kth iteration (), personal best and global best are obtained, and if , return and end the algorithm, else return to step 2.
The inner PSO algorithm is to find the best solution to allocate insufficient resources, and the fitness function only optimizes the resource allocation part. Thus, based on the objective function of and , can be represented as follows:
The random initialization resource allocation solution from EDCs to each resource demand points at time t is based on constraints (11) and (25). The initial PSO resource allocation solution of the ith particle () can be represented by a fourdimensional array of all , , , and in the model in Section 3.2. The updating method of PSO is based on equations (31) and (32). If does not satisfy constraints (11) and (25), (s, AESs, BESs, and AAs) will be set to 0.
4.3. Improved SAPSO Algorithm
The 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 realworld 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. Thus, the outer site selection problem is solved by the ISA algorithm.
4.3.1. GA Algorithm
The main steps of the GA algorithm are as follows: Step 1. Initialize the parameters: population size , maximal number of generations , crossover rate , and mutation rate , set as initial solution of the th individual (), calculate the fitness function of the th individual, and find . 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. Step 3. Calculate the fitness function of the th individual in the jth iteration and find (, ), and if , return and end the algorithm, else return to step 2.
4.3.2. ISA Algorithm
The main steps of the GA algorithm are as follows: Step 1. Initialize the parameters of SA: start temperature , end temperature , 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 , crossover rate , and mutation rate . Set as the initial solution of the th individual (), calculate the fitness function of the th individual, and find . Step 2. Iterate times in each temperature, use roulette wheel selection to select the new population from prior population Pop, update the new population Pop' by randomly roulette wheel selection, crossover, and mutation, calculate the fitness function of the individual in the iteration , and find , . According to the Metropolis guidelines, decide if is accepted as the best solution and if is accepted as the new population with the following probability: Step 3. Outer loop: if is lower than or , end the algorithm; else if , 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 . Combining the characteristics of GA evolution, the update solution of ISA includes selection, crossover, and mutation operation.
5. 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 tradeoff 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 to and a longitude of to , mainly in the North China earthquake zone Yinmiddle Yanshan seismic zone. The 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. There 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.
5.1. Data Processing
We divide the first 3day gold rescue time after a disaster into 6 hours of time intervals. Thus, the emergency evacuation and resource allocation network are designed in 12 (=72/6) intervals.
5.1.1. 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. Thus, based on the Beijing Central City Earthquake and Emergency Refuge Site Planning Outline and the Chaoyang Government Network, another sixteen candidate shelters were selected. The EMCs need professional medical equipment and nursing staffs. Thus, 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 assume that candidate EDCs are located at street/district offices. So, there are 42 candidate EDCs.
5.1.2. Population
Chaoyang district governs 24 streets and 19 districts, and the population size of each street/district is extracted from the Chaoyang District Beijing Statistical Yearbook 2017. The number of refugees in each AA of each time interval is roughly estimated [17], and Figure 8 presents the ratio of to . We assume that 10% of the population is injured , 20% of the population is assumed to be Apriority refugees , and 30% is assumed to be Bpriority refugees .
5.1.3. 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.

5.1.4. Capacity
We assume that each ES is divided into two parts: AES and BES, where the capacity of AESs is smaller than BES. The capacities of AES, BES, and EMCs are presented in Tables 3 and 4.


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–33]. The capacity of each EDC can be derived from the total number of EDCs, total quantity of required resources, and the ratio .
5.1.5. 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 AESs, BESs, and EMCs and from EDCs to AAs, AESs, BESs, and EMCs using version 10.1 of ArcGIS.
5.1.6. Cost
The total cost involves the construction cost of candidate AESs, BESs (existing shelters have no construction cost), and EDCs and the fixed and transportation cost of per unit resource . The construction cost of each site is determined by its capacity, which is estimated by the government report and previous research studies [5, 18, 19]. The 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.

5.2. 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 AESs BESs AAs and (s, AESs, BESs, and AAs), we set the priorities of different places , , , and 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. 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 . 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, , , and 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. is the start temperature, 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, and are the learning factors, indicates the possibility of the particle updating according to its local optimal position, indicates the possibility of the particle updating according to the global optimal position, is the inertia factor, which reflects the particle’s ability to inherit the previous speed, and and are the weight setting of the objective functions and . For parameters of GA, G represents the number of groups, is the crossover rate, and 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 of SA, we calculate the minimum/maximum value of each objective of our two test cases (as shown in Table 1).
We firstly test the case with these existing shelters as case 1, in which the construction cost of AESs and BESs set to be zero. Then, these candidate shelters are added in the second test case as case 2, using the parameters shown in Table 6.
Through 50 experiments, the average value of objective functions of the two test cases is shown in Table 7.

The dividend is calculated by the equation , 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. Thus, 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 AESs and BESs at time of the two cases are shown in Figure 9.
(a)
(b)
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 AESs and six BESs were selected at time t = 1 (one more AES and BES were selected at time ) and optimized the longdistance allocation result. Through 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.
5.3. 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. Therefore, all the scenarios in this section are based on case 2.
Scenario 1. The first scenario is to compare the efficiency of SAPSO and ISAPSO. Parameters of the GA are shown in Table 8. We analyze these two algorithms by comparing the results of case 2 under SAPSO and ISAPSO. The 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 ISAPSO is always better than that of SAPSO. Table 9 presents the objective values calculated by ISAPSO. The dividend is calculated by the equation , and it can be seen that compared with SAPSO, the objective function values by ISAPSO can achieve Pareto optimization. But, because of the characteristics of population search of ISAPSO, the time consumption of ISAPSO is correspondingly longer.
The algorithm used in the following calculations will use ISAPSO. And, the following comparison of objective values will be based on the results of ISAPSO.


Scenario 2. The second scenario is to analyze the impact of the objective weight. Thus, we set two of the objective weights to 50% of the original and the other target weight to twice the original: S21: ; S22: ; S13: . Through 50 experiments, the average value of three objective functions is presented in Table 9.
It can be seen from Table 9 that is greatly affected by the objective weight. When the weight increases, the objective value is optimized, while the objective value of and is affected by the weight to a small extent. In addition, it can be seen that and are negatively correlated with , and increases as and decrease. The corresponding dividend comparison of three objective values of Scenario 2 is shown in Figure 11.
The dividend is calculated by the equation , and it can be inferred that is the most sensitive to the objective weight. Therefore, provided that the budget is sufficient, increasing the objective weight of appropriately can better optimize the result.
Scenario 3. The third scenario is to analyze the effect of parameter on three objective functions. The 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 S31: , S32: , and S33: , respectively. Through 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 and . It implies that resources and their transportation costs account for only a small part of the total cost (). And, obviously, as increases, decreases. Figure 12 illustrates the corresponding dividend comparison of Scenario 2. The dividend is calculated by the equation . It is observed that as increases, decreases rapidly and there were almost no change in and . This scenario illustrates that adjusting the parameter significantly changes objective values of in a linear relationship, but does not affect the other objects. This implies that increasing the quantity of resources is an effective method that improves the target value of . It can be seen that preparing sufficient resources in advance can play a vital role in the disaster relief process at a small cost.

Scenario 4. This scenario is to analyze the impact of population in each priority. Here, we increase the number of Apriority and Bpriority refugees to be 30% and 40% of the population in S41 and 40% and 50% of the population in S42.
In this scenario, the total capacity of the existing shelters cannot meet the refugee demand. Thus, the selecting of candidate shelters is essential when the refugee demand increases.
Through 50 experiments, the average value of three objective functions of Scenario 4 is presented in Table 11. The dividend is calculated by the equation . It is obvious that when the number of A and Bpriority 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 AESs and BESs 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 . 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 A and Bpriority refugees increases, the number of refugees in AESs and BESs increases as well. This result demonstrates that refugees in AESs and BESs are more sensitive to unmet resources than those in AAs.

(a)
(b)
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. Through 50 calculations, the average value of three objective functions of Scenario 5 is presented in Table 12. The dividend is calculated by the equation . The 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 timevarying demand. The reason of the evacuation distance dividend is that the remaining capacity of AESs and BESs becomes smaller as time goes by, but the number of A and Bpriority refugees has not become smaller. Therefore, refugees in the later period will have to evacuate to farther AESs and BESs.
Similarly, most of the resources in EDCs were allocated to the AAs in the early period. The remaining resources in EDCs are getting less and less as time goes by. However, the demand for resources of AESs and BESs is increasing much faster than the timevarying demand model. Therefore, the resource requirements of AESs and BESs 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. Thus, it is obvious that the timevarying 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.

6. Conclusions
This paper studies decisionmaking on the dynamic emergency shelter and resource allocation for disaster management, considering the hierarchical requirement of refugees and allocation of insufficient relief resources. The 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. Thus, a multiobjective locationallocation 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. The 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 SAPSO 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. The computational experiment results indicate that selecting candidate shelters for refugees is obviously more worthwhile and more effective than evacuating refugees to existing shelters. The 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. The 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 realtime road network like congestion and destruction should be considered in the refugee evacuation planning. Thus, 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
The numerical study data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The work described in this article was supported by National Natural Science Foundation of China (NSFC) under project no. 71671071, no. 71771154, and no. 91846301.
References
 N. Laframboise and B. Loko, “Natural disasters: mitigating impact, managing risks,” IMF Working Papers, vol. 12, no. 245, p. 32, 2012. View at: Publisher Site  Google Scholar
 S. Roh, S. Pettit, I. Harris, and A. Beresford, “The prepositioning of warehouses at regional and local levels for a humanitarian relief organisation,” International Journal of Production Economics, vol. 170, pp. 616–628, 2015. View at: Publisher Site  Google Scholar
 X. Zhao, W. Xu, Y. Ma, and F. Hu, “Scenariobased multiobjective optimum allocation model for earthquake emergency shelters using a modified particle swarm optimization algorithm: a case study in Chaoyang district, Beijing, China,” PLoS One, vol. 10, no. 12, Article ID e0144455, 2015. View at: Publisher Site  Google Scholar
 C. L. Hu, X. Liu, and Y. K. Hua, “A biobjective robust model for emergency resource allocation under uncertainty,” International Journal of Production Research, vol. 54, no. 24, pp. 7421–7438, 2016. View at: Publisher Site  Google Scholar
 H. Li, L. Zhao, R. Huang, and Q. Hu, “Hierarchical earthquake shelter planning in urban areas: a case for Shanghai in China,” International Journal of Disaster Risk Reduction, vol. 22, pp. 431–446, 2017. View at: Publisher Site  Google Scholar
 Z. Chen, X. Chen, Q. Li, and J. Chen, “The temporal hierarchy of shelters: a hierarchical location model for earthquakeshelter planning,” International Journal of Geographical Information Science, vol. 27, no. 8, pp. 1612–1630, 2013. View at: Publisher Site  Google Scholar
 H. S. Na and A. Banerjee, “A disaster evacuation network model for transporting multiple priority evacuees,” IIE Transactions, vol. 47, no. 11, pp. 1287–1299, 2015. View at: Publisher Site  Google Scholar
 A. BozorgiAmiri and M. Khorsi, “A dynamic multiobjective location–routing model for relief logistic planning under uncertainty on demand, travel time, and cost parameters,” The International Journal of Advanced Manufacturing Technology, vol. 85, no. 5–8, pp. 1633–1648, 2016. View at: Publisher Site  Google Scholar
 S. Elloumi, “A tighter formulation of the pmedian problem,” Journal of Combinatorial Optimization, vol. 19, no. 1, pp. 69–83, 2010. View at: Publisher Site  Google Scholar
 C.C. Lu and J.B. Sheu, “Robust vertex pcenter model for locating urgent relief distribution centers,” Computers & Operations Research, vol. 40, no. 8, pp. 2128–2137, 2013. View at: Publisher Site  Google Scholar
 Y. Ma, W. Xu, L. Qin, and X. Zhao, “Site selection models in natural disaster shelters: a review,” Sustainability, vol. 11, no. 2, p. 399, 2019. View at: Publisher Site  Google Scholar
 A. Trivedi and A. Singh, “A hybrid multiobjective decision model for emergency shelter locationrelocation projects using fuzzy analytic hierarchy process and goal programming approach,” International Journal of Project Management, vol. 35, no. 5, pp. 827–840, 2017. View at: Publisher Site  Google Scholar
 W. Xu, Y. Ma, X. Zhao, Y. Li, L. Qin, and J. Du, “A comparison of scenariobased hybrid bilevel and multiobjective locationallocation models for earthquake emergency shelters: a case study in the central area of Beijing, China,” International Journal of Geographical Information Science, vol. 32, no. 2, pp. 236–256, 2018. View at: Publisher Site  Google Scholar
 C.H. Tsai and Y.L. Yeh, “The study of integrating geographic information with multiobjective decision making on allocating the appropriate refuge shelters: using kengting national park as an example,” Natural Hazards, vol. 82, no. 3, pp. 2133–2147, 2016. View at: Publisher Site  Google Scholar
 F. Hu, S. Yang, and W. Xu, “A nondominated sorting genetic algorithm for the location and districting planning of earthquake shelters,” International Journal of Geographical Information Science, vol. 28, no. 7, pp. 1482–1501, 2014. View at: Publisher Site  Google Scholar
 J. Xu, X. Yin, D. Chen, J. An, and G. Nie, “Multicriteria location model of earthquake evacuation shelters to aid in urban planning,” International Journal of Disaster Risk Reduction, vol. 20, pp. 51–62, 2016. View at: Publisher Site  Google Scholar
 J.B. Sheu and C. Pan, “A method for designing centralized emergency supply network to respond to largescale natural disasters,” Transportation Research Part B: Methodological, vol. 67, pp. 284–305, 2014. View at: Publisher Site  Google Scholar
 O. RodríguezEspíndola and J. Gaytán, “Scenariobased preparedness plan for floods,” Natural Hazards, vol. 76, no. 2, pp. 1241–1262, 2015. View at: Publisher Site  Google Scholar
 J. Dalal and H. Üster, “Combining worst case and average case considerations in an integrated emergency response network design problem,” Transportation Science, vol. 52, no. 1, pp. 171–188, 2017. View at: Publisher Site  Google Scholar
 H. Su, “Siting and sizing of distributed generators based on improved simulated annealing particle swarm optimization,” Environmental Science and Pollution Research, vol. 26, no. 18, pp. 17927–17938, 2019. View at: Publisher Site  Google Scholar
 W.S. Li, W.P. Hsiung, O. Po, K. S. Candan, and D. Agrawal, “Evaluations of architectural designs and implementation for databasedriven web sites,” Data & Knowledge Engineering, vol. 43, no. 2, pp. 151–177, 2002. View at: Publisher Site  Google Scholar
 T.H. Wu, C.C. Chang, and S.H. Chung, “A simulated annealing algorithm for manufacturing cell formation problems,” Expert Systems with Applications, vol. 34, no. 3, pp. 1609–1617, 2008. View at: Publisher Site  Google Scholar
 S. Bandyopadhyay, S. Saha, U. Maulik, and K. Deb, “A simulated annealingbased multiobjective optimization algorithm: Amosa,” IEEE Transactions on Evolutionary Computation, vol. 12, no. 3, pp. 269–283, 2008. View at: Publisher Site  Google Scholar
 H. Zheng, X. Hu, M. Zheng, and R.j. Liu, “An improved hybrid algorithm based on particle swarm optimization and simulated annealing and its application,” Computer Technology for Developing Areas, vol. 23, no. 7, pp. 26–30, 2013. View at: Google Scholar
 Q. Cai, M. Gong, B. Shen, L. Ma, and L. Jiao, “Discrete particle swarm optimization for identifying community structures in signed social networks,” Neural Networks, vol. 58, pp. 4–13, 2014. View at: Publisher Site  Google Scholar
 P.Y. Yin, S.S. Yu, P.P. Wang, and Y.T. Wang, “Multiobjective task allocation in distributed computing systems by hybrid particle swarm optimization,” Applied Mathematics and Computation, vol. 184, no. 2, pp. 407–420, 2007. View at: Publisher Site  Google Scholar
 D. P. Culhane and S. Metraux, “Rearranging the deck chairs or reallocating the lifeboats? homelessness assistance and its alternatives,” Journal of the American Planning Association, vol. 74, no. 1, pp. 111–121, 2008. View at: Publisher Site  Google Scholar
 J. H. Zhang, “Research of improved simulated annealing optimization algorithm based on the global harmony search mechanism,” Advanced Materials Research, vol. 482–484, pp. 2500–2503, 2012. View at: Publisher Site  Google Scholar
 A. K. Mollah, S. Sadhukhan, P. Das, and M. Z. Anis, “A cost optimization model and solutions for shelter allocation and relief distribution in flood scenario,” International Journal of Disaster Risk Reduction, vol. 31, pp. 1187–1198, 2018. View at: Publisher Site  Google Scholar
 G. Nastasi, V. Colla, S. Cateni, and S. Campigli, “Implementation and comparison of algorithms for multiobjective optimization based on genetic algorithms applied to the management of an automated warehouse,” Journal of Intelligent Manufacturing, vol. 29, no. 7, pp. 1545–1557, 2018. View at: Publisher Site  Google Scholar
 A. BozorgiAmiri, M. S. Jabalameli, and S. M. J. Mirzapour AleHashem, “A multiobjective robust stochastic programming model for disaster relief logistics under uncertainty,” OR Spectrum, vol. 35, no. 4, pp. 905–933, 2013. View at: Publisher Site  Google Scholar
 D. Khayal, R. Pradhananga, S. Pokharel, and F. Mutlu, “A model for planning locations of temporary distribution facilities for emergency response,” SocioEconomic Planning Sciences, vol. 52, pp. 22–30, 2015. View at: Publisher Site  Google Scholar
 Z. Yang, H. Zhou, X. Gao, and S. Liu, “Multiobjective model for emergency resources allocation,” Mathematical Problems in Engineering, vol. 2013, Article ID 538695, 6 pages, 2013. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Jian Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.