Research Article  Open Access
Yonghong Liu, Yucheng Li, De Huang, "A Multiobjective Optimization Model for Continuous Allocation of Emergency Rescue Materials", Mathematical Problems in Engineering, vol. 2020, Article ID 5693182, 15 pages, 2020. https://doi.org/10.1155/2020/5693182
A Multiobjective Optimization Model for Continuous Allocation of Emergency Rescue Materials
Abstract
Emergency rescue operations play a vital role in alleviating human suffering, reducing casualties, and cutting down economic losses. One key aspect in the management of these operations is the rational allocation of emergency relief materials, where the allocation is continuous, dynamic, and concurrent. This allocation should be made not only to minimize the emergency rescue losses, but also to reduce the cost of emergency rescue work. A reasonable and effective allocation scheme for emergency relief materials can be established to adapt to the continuity, dynamics, and concurrency of material distribution. In this work, we propose a multiobjective optimization model of emergency material allocation with continuous timevarying supply and demand constraints, where the objective is to minimize the losses and the economic cost incurred by the emergency rescue operations. The constrained optimization problem is handled through sequential unconstrained minimization techniques, and the multiobjective optimization is carried out by the fast nondominated sorting genetic algorithm (NSGAII) with an elite strategy to obtain a Pareto solution set with fairness and balance of loss and cost. The loss and cost associated with the Pareto frontier are employed to find an appropriate noninferior solution and its corresponding material allocation scheme. We verify through several simulations the model feasibility and the effectiveness of the proposed method, which can provide decision support for continuous material allocation in emergency rescue operations.
1. Introduction
Human populations are frequently impacted by natural disasters such as earthquakes, tsunamis, floods, typhoons, volcanic eruptions, and debris flow, as well as major accidents such as those of mines, traffic, and fires. Such disasters and accidents not only harm the lives and properties of people, but also seriously hinder economic growth and society development [1, 2]. The emergency rescue operations, in response to such natural disasters and accidents, should be highly active, effective, and quick. These requirements promoted extensive research and investigations on the issues of coordination among emergency rescue units, rescue route selection, and material distribution. In particular, effective allocation of emergency materials is one of the most essential and key problems in rescue operations [3, 4]. A reasonable and effective material allocation scheme can complete the emergency rescue work at the lowest cost and loss. More importantly, such a scheme can also reduce, to a certain extent, the personal and property losses caused by disasters and accidents [5, 6].
Fiedrich et al. [7] introduced a dynamic optimization model to formulate the problem of minimizing the total number of deaths after an earthquake. This problem was solved using simulated annealing and Tabu search methods to determine the best resource allocation plan for earthquake sites. Balcik and Beamon [8] used mathematical models to investigate the facility location problem in humanitarian relief chains. Hence, they determined the numbers and locations of distribution centres in relief networks, as well as the quantities of needed relief supplies for disasteraffected areas. Li et al. [9] reviewed coverage models and optimization techniques for the location and planning of emergency response facilities. The authors particularly reviewed stateoftheart models including the hypercube queuing model, the dynamic allocation model, the progressive coverage model, and the cooperative coverage model. Wex et al. [10] developed a decisionsupport model for rescue unit allocation and arrangement. The model seeks to minimize the sum of rescue completion times weighted by incident severity levels. A metaheuristic algorithm was used to solve this optimization problem and significantly reduce rescue times. Coordination among rescue units and independent agencies for fast material delivery was also discussed as a potential research direction. Yang et al. [11] proposed a decisionsupport system for public safety and emergency management in patrol service centres. This system shows a high potential for the allocation, deployment, and scheduling of resources in order to improve the efficiency of patrolrelated resources and services and boost emergency response capabilities. However, this system is not suitable for realtime or dynamic scenarios. Zhou et al. [12] proposed a multiobjective optimization model for the multiperiod dynamic emergency resource scheduling problem. Schryen et al. [13] proposed a method to support decision makers in sudden disaster situations with high risk, time pressure, uncertainty, conflict, or lack of information. The experimental results showed small differences between heuristicbased and optimal solutions over short periods of time. Chen and Yu [14] used network partitioning and integer programming to find temporary locations for the transportation infrastructure of emergency medical services (EMS). Taking into consideration worstcase EMS demands that arise after the onset of a disaster, the authors expect that their analysis results can help decision makers with site selection issues and can serve as a benchmark for planning postdisaster emergency services. Duhamel et al. [15] proposed a mathematical model and a heuristic algorithm for a multiperiod locationallocation problem in postdisaster relief operations. They also related the proposed quantitative model to another model for disasterresistant systems. Pradhananga et al. [16] proposed a threelevel network model for integrated emergency preparedness, response planning, and emergency material distribution. This model seeks to minimize social costs through determining a set of potential resource supply points at the highest level. The results showed that prearrangement of some items and postdisaster procurement can reduce the shortage of emergency supplies. Alem et al. [17] introduced a new twostage stochastic network traffic model for providing rapid humanitarian assistance to victims after a disaster. This model takes into account typically overlooked elements in rescue operations, such as budget allocation, fleet sizing, and procurement and lead times of dynamic multicycle changes. Sung and Lee [18] modelled the ambulance routing problem as a resourceconstrained classification problem. Then, they determined the source and destination hospitals for evacuating patients. A column generation algorithm was used to solve the formulated fixedpriority resource allocation problem. Though the best obtained solution may require complicated operations, it may still be easy to follow and implement in a postdisaster environment. Bai [19] proposed a new twostage optimization method for emergency material distribution with multiple vendors, disaster areas, relief resources, and vehicles. Criteria of fairness, timeliness, and economic efficiency are considered in the formulation of this multiobjective optimization model for facility location, vehicle routing, and resource allocation decisions. The goals of this model are to minimize the proportion of unsatisfied demands, minimize the response time of emergency relief, and minimize the total cost of the entire process. Deng et al. [20] investigated the problem of twostage randomcapacitylimit locationallocation in emergency logistics. To solve this problem, the authors established a random expectedvalue model and its deterministic counterpart. Then, they developed an improved particle swarm optimization algorithm with restart of the Gaussian cloud operator and an adaptive parameter selection strategy.
Boonmee et al. [21] studied facility location problems in emergency humanitarian logistics, where they analyzed facility location types, disaster characteristics, and decisionmaking approaches. They also investigated data models, typical formulations of objectives and constraints, solution methods, practical application procedures, and case studies. Liu et al. [22] addressed the allocation of emergency relief materials after chemical leakage accidents in rivers. The authors proposed a material allocation framework that minimizes the response time for dynamic emergency material demands. They established a multiobjective emergency resource allocation model that seeks to minimize both the total resource allocation cost and the total loss caused by resource insufficiency. Particle swarm optimization was used to obtain the optimal solutions. Nevertheless, the emergency resource allocation process did not consider the response time, resource urgency, and disaster severity. Guo et al. [23] addressed the problem of distributing marine emergency rescue materials. The authors formulated this multipleresource allocation problem within a multiobjective optimization framework that considers the longrange nature of maritime rescue operations and sustainable recycling of emergency resources. The fast nondominated sorting genetic algorithm (NSGAII) was used for solving this optimization problem, and the candidate solution set approaching the Pareto frontier was exploited to find a suitable noninferior resource allocation scheme. The authors adopted a discrete resource allocation problem formulation with decreasing proportional constraints and proposed expertsystem mathematical models for solving this problem. This formulation has the advantage of a decreasing ratio which can prevent unfair material allocation, but it did not account for continuous material distribution.
Therefore, material allocation is generally formulated as a multiobjective and multiconstrained optimization problem, for which reasonable and effective solutions are crucially needed to ensure the success of emergency rescue work. At present, existing approaches for emergency relief material allocation focus on selecting the locations of the relief materials. An allocation scheme is typically determined based on known locations of material storehouses and material storage capacities. However, a welldesigned allocation scheme for emergency relief materials needs to consider the emergency rescue time, economic costs, the loss of materials that are undeliverable on time, and the balance of the material distribution plan. Therefore, this paper comprehensively considers the continuity, dynamics, and concurrency of emergency relief material allocation. Continuity of material allocation is defined for materials that are genuinely necessary and continuously consumed in emergency rescue sites, such as water, medicines, and food. Towards continuous material allocation, we construct a multiobjective optimization model, based on a timevarying order of material transportation, for the minimization of the conflicting factors of the cost and loss of the emergency rescue operations. Based on existing methods and algorithms for emergency material allocation with timevarying order of material transportation, we choose to minimize two conflicting objectives, namely, the economic costs and the losses. Also, we adopt the NSGAII heuristic algorithm to solve this optimization problem and obtain Pareto frontier solutions. From this frontier, a noninferior solution is selected to get an economical, reasonable, and effective emergency relief material allocation plan that achieves fairness and balance between cost and loss. Such a plan can provide decision support for decision makers on the best continuous allocation schemes for emergency rescue materials. The experimental results show that the solution of the emergency material distribution problem is feasible based on criteria of fairness and balance of material distribution. Also, the solution shows characteristics of timeliness, timevarying dynamics, and continuity.
The remainder of the paper is structured as follows. The proposed multiobjective material allocation optimization model is discussed in Section 2. Section 3 presents the algorithmic solution steps of that optimization problem. Section 4 provides case studies of the emergency rescue material allocation. In Section 5, the experimental results are presented and discussed. Finally, conclusions are drawn in Section 6.
2. Multiobjective Optimization for Emergency Material Allocation
2.1. Problem Description
After a natural disaster or a major accident, emergency rescue decision makers should direct rescue teams for people rescue missions and material allocation. A reasonable and effective allocation plan of relief materials should have the shortest time and the lowest cost, while minimizing the losses caused by the natural disaster or accident. Indeed, creating a multiobjective optimization model with theses optimality criteria can lead to sought a reasonable and effective emergency rescue plan.
The types, total quantities, and the total number of storage locations of emergency relief materials are typically known. Let n be the number of the candidate storage units for the emergency rescue materials. This number is set based on the estimated emergency rescue time, the time when each storage unit arrives at the rescue site, and the number of units that fail to deliver materials within the expected rescue time. The rescue time , the total number of rescue material types , and the quantity of each material type () are estimated based on the emergency damage level. Disregarding interruptions in material allocation, let the shortest time for transporting the material storage unit to the rescue location be . We seek to determine candidate storage units for emergency rescue materials and complete the allocation of the materials in time . Based on timevarying sorting, the times at which the candidate storage units (i = 1, 2, …, n) arrive at the rescue locations are sorted such that . The capacity for the material type provided by the storage unit is denoted by , where . Let the decision variable indicate whether a quantity of the material type is assigned from the material storage unit to the accident location . The problem of the emergency relief material allocation should be solved for in order to achieve the shortest rescue time, the lowest cost of material allocation, and the minimization of losses caused by the emergency. The notation used in the model is presented in Table 1.

2.2. Proposed Mathematical Model
Based on the above problem description, the rescue time is limited by the distance of the material rescue storage unit to the accident location. As a practical allocation scheme for emergency rescue materials at the minimum cost while also minimizing the emergency losses, we propose the following multiobjective mathematical optimization model for continuous allocation of emergency rescue materials. The objective functions arewhere
The constraints are
The objective function (1) uses the coefficient C_{ij} as a variable objectivedependent cost. In particular, the choice of this coefficient depends on the relative magnitude of the relief material cost and the minimum material deployment time. On the one hand, if the objective is to minimize the deployment time, the transportation cost will be higher. On the other hand, if the objective is to minimize the transportation cost, the deployment time may be longer, and hence the loss can be larger. Practically, the appropriate transportation mode should be selected according to the actual situation, in order to appropriately determine the coefficient C_{ij}. The objective function (2) represents the losses inflected by the emergency on the people, environment, or economy. These losses are determined based on the nature and degree of the emergency. For the emergency rescue operations, the main influencing factors include the amount of available materials, the loss caused by the lack of materials, and the consumption of materials. In the model mentioned above, the loss includes not only the economic and property losses, but also the damage caused by the lack of lifesaving essentials. The loss caused by the lack of materials is based on the magnitude of the material influence on people, animals, or environments. This loss is quantified as a value in the unit interval [0, 1].
For the objective (2), represents the amount of missing material of type j after the storage unit B_{i} arrives at the emergency site. The period represents the time difference between the two latest (current and previous) arrival times of the materials from the storage unit B_{i}. In this period, the amount of consumed material of type j is , the remaining amount of this material at the last moment is , and the amount of delivered material is . From these three amounts, the amount of lost material can be obtained. The amount of the remaining material is given by (4). When the sum of the amount of the delivered material and the amount of the surplus material is greater than the amount of the consumed material , then the amount of missing materials is zero. Otherwise, the amount of missing materials is the difference between the consumed materials and the sum of the delivered and remaining materials. The constraint (5) indicates that the demand for a material of type j must be more than or equal to the consumption within the rescue time T. Also, under conditions of material sufficiency, the material demand needs to be less than or equal to the sum of the materials in the storage units. The constraint (6) indicates that the amounts of various rescue materials deployed from different material storage units cannot exceed the amounts in the storage units. The constraint (7) shows that the total sum of the relief supplies should be more than the demand of the relief materials. The constraint (8) indicates that the longest time for a material storage unit to arrive at a rescue site cannot exceed the expected rescue time, where the material storage units are sorted by the material transportation times. Side constraints on the variables are given in (9).
3. Solution of the Multiobjective Material Allocation Problem
According to the established mathematical model (1)–(9), the two objectives (1) and (2) conflict with each other. Hence, it is difficult to achieve the minimum transportation cost and loss at the same time. For such multiobjective optimization problems, Paretooptimal solutions can be obtained through swarm intelligence heuristicsearch stochastic optimization algorithms. We analyzed and compared the improved harmony search algorithm (IHS) [24], the multiobjective evolutionary algorithm based on decomposition (MOEA/D) [25], and the fast nondominated sorting genetic algorithm (NSGAII) with an elite strategy [26]. The NSGAII method appears more appropriate and effective in solving the established material allocation problem. Indeed, the NSGAII method typically shows the best performance and has been widely used in resource allocation and inventory management problems [27–29]. As well, this method has been effective in solving multiobjective optimization problems, especially for 23 objective functions [30, 31]. This method essentially uses fast nondominated sorting and a diversitypreserving elite strategy to sort the parent and progeny populations at different levels of the nondominated solution frontier. A density estimation method based on the crowding distance is adopted to maintain population diversity.
3.1. NSGAII Algorithm
Fast nondominated sorting of NSGAII reduces the computational complexity of this algorithm. Moreover, this lowcomplexity algorithm uses a congestion degree and a congestion degree comparison operator to replace the fitness sharing strategy which needs to specify a shared radius. As well, NSGAII sets the winning criterion in the same level of comparison after sorting, so that individual solutions in the Pareto domain can be extended to the whole domain in an evenly distributed diversified manner. In addition, NSGAII introduces an elite strategy, which enlarges the sampling space, combines the parent population with its offspring populations, and creates competition for producing a bestfit nextgeneration population. By layering all individuals in the population, the best individuals will survive, and the population quality will be raised rapidly [32].
Based on the NSGAII algorithm, the steps for solving the multiobjective optimization problem of emergency material allocation are given as follows: Step 1. Model Parameter Setting and Initialization. Let T be the predicted emergency rescue time, m be the number of rescue material types, n be the number of candidate material storage units, ν_{j} be the consumption of the jth rescue material per unit time, and L_{j} be the loss caused by the lack of the jth rescue material. Based on the maximum time t_{i} required for each material repository to reach the rescue site, the material repositories B_{i} are ranked. Thus, for each repository B_{i} and material type j, the storage capacity S_{ij} and the storage cost C_{ij} are obtained. The material demand q_{j} is initialized according to the material storage, material consumption, and estimated rescue completion time. Step 2. Initialization of the NSGAII Parameters. The following NSGAII parameters are initialized: the number of individuals , the number of generations , the crossover probability , the crossover distribution index , the mutation probability , and the mutation distribution index . Step 3. Coding Strategy and Chromosome Design. The sought design variables are all of the decision variables x_{ij} for allocating a quantity of the jth material from the ith rescue material storage unit. According to the parameter setting and initialization [30, 31, 33], the chromosome associated with each individual solution of the overall population is encoded as . Each chromosome is initialized by real numbers according to the range set by the constraints in (6). Step 4. Nonlinear Constraint Processing. The sequential unconstrained minimization technique is used to deal with the nonlinear constraints in (7), and hence obtain a multiobjective unconstrained minimization problem [34]. The unconstrained objective is obtained from the constrained objective and the nonlinear constraints as follows [34]: where is the jth constraint, is the number of constraints, is the constrained objective function, is the objective function value corresponding to the worst feasible solution in the population, <> is a penalty function of , and the operator <·> denotes the absolute value of the operand. The definition of the penalty function here is different from those of the traditional methods and earlier genetic algorithms. The following rules of solution preference are adopted: (1) any feasible solution is superior to any infeasible solution; (2) among two feasible solutions, the solution with the better objective function value is preferred; (3) among two infeasible solutions, the solution with less conflict of constraints is preferred. Step 5. Fitness Calculation. Before calculating the fitness of a solution, the chromosome coding in Step 3 needs to be decoded into an matrix, as shown in Figure 1. The matrix rows represent the amounts of demand for each emergency supply repository, and the matrix columns indicate the amounts of each material that needs to be allocated from each repository. According to the decoded decision , the fitness values are obtained by calculating the objective functions (1) and (2). Step 6. Fast Nondominated Sorting. According to the order of dominance (or preference) of the Pareto solutions, the population is divided into subpopulation levels. The similarity between subpopulations is evaluated on the Pareto frontier, and the resulting subpopulation and similarity measures are used to promote the diversity of nondominant solutions. The calculated fitness values are stratified by the fast nondominated sorting method. Thus, the population density near a specific solution can be judged. The crowding degree of each individual in the nondominated frontier is calculated according to the difference between the fitness values. So, the population density near a specific solution in a population is determined, and the density can be obtained through pairwise competition and selection. Step 7. Crossmutation Genetic Operations. According to the initial crossover and mutation rates, the crossover and mutation genetic operations are carried out on the new population with the smallest population density. New individuals are generated, subjected to a cross mutation operation, and then merged with the parent population to obtain a new population. Step 8. Termination and Output. The fitness value of Step 5 is calculated for the individual population in Step 7, and the ranking method by Step 6 is used to reorder and stratify the population, in order to generate a new population in Step 7. This process is iteratively repeated until either the preset termination condition or the maximum number of iterations is reached. Hence, the optimal Pareto frontier is output.
3.2. IHS Algorithm
The IHS multiobjective optimization method [24] is used to solve nonlinear nonconvex problems. This method generates new solution vectors to improve the accuracy and convergence speed of the harmony search algorithm (HS). The HISbased solution steps for the optimization problem of emergency rescue multitarget material allocation are as follows: Step 1. The same as that of NSGAII. Step 2. Initialization of the IHS Parameters. The following IHS parameters are initialized: the number of individuals , the number of generations , the probability of choosing from memory , the minimum pitch adjustment rate , and the maximum pitch adjustment rate . Step 3 to Step 5. The same as those of NSGAII. Step 6. Improvisation of a New Harmony. Based on memory consideration, pitch adjustment, and random selection, the new harmony is generated as follows:(1)Memory consideration:(2)Pitch adjustment decision:(3)Random selection: If the pitch adjustment decision for , then where is an arbitrary distance bandwidth and r_{d} is a random floating number between 0 and 1. Step 7. Solution Update. Iterate over x_{i}, where i = 1, 2, …, n, and exclude the worst existing harmony from the solution X. Step 8. Termination and Output. If the preset termination condition or the maximum number of iterations is reached, the algorithm is terminated. Otherwise, Steps 6 to 8 are repeated.
3.3. MOEA/D Algorithm
In this algorithm, a multiobjective optimization problem is decomposed into several scalar optimization subproblems, which are then simultaneously optimized. The solution of each optimization subproblem uses only information from its adjacent subproblems. This significantly reduces the computational complexity of each MOEA/D generation. The MOEA/Dbased solution steps for the optimization problem of emergency rescue multitarget material allocation are as follows: Step 1. The same as that of NSGAII. Step 2. Initialization of the MOEA/D Parameters. The following MOEA/D parameters are initialized: the number of individuals , the number of generations , the size of the weight neighborhood , and set the external population . Step 3. Distance Computation. For each i = 1, 2, …, N_{I}, set , where are the n_{d} closest weight vectors to . Step 4. Initial Population Generation. Generate randomly an initial population , set , and initialize using equations (1) and (2). Step 5. Reproduction. Randomly select two serial numbers b and l from . Use genetic operators to generate a new solution y from and . Based on equations (1) and (2), apply the improvement heuristic on y to produce . Step 6. Update of the Objective Values z and Neighbouring Solutions. For each , then set . If , , then set and . Step 7. External Population Update. For j = 1 and 2, remove from the external population (EP) all of the vectors dominated by . Also, add to EP if no vectors in EP dominate . For each i = 1, 2, …, N_{I}, loop through the execution of Step 4 to Step 7. Step 8. Termination and Output. If any of the stopping criteria are satisfied, then stop and output EP. Otherwise, go to Step 3.
4. Experimental Setup
In this section, we verify the validity of the multiobjective mathematical optimization model for the allocation of rescue materials. We check as well the feasibility of the solutions obtained by the NSGAII algorithm. Given n rescue material repositories, m types of rescue materials, an estimated rescue time T, and the constraints (5)–(9), we generated randomly the material demands, the storage capacities, the material transportation costs, and the loss due to material shortage. From this simulation, an example allocation scheme for emergency rescue materials is generated. Optimal solutions of the established optimization problem are found using numerical simulations of the NSGAII algorithm.
4.1. Experimental Design
According to actual situations, the number of the relief material storage units n that reach the rescue site within the estimated rescue time is chosen in the interval [3, 30]. Also, the range for the number of types of rescue materials m is set to [15, 120]. The minimum rescue time T is expected to be 3 hours, while the maximum rescue time should be 72 hours. Given the total number of rescue material repositories n and the number of types of the rescue materials m, the rescue time is randomly generated, and 45 test cases are obtained. As shown in Table 2, for the number of the material storage units n, the maximum and minimum values are 30 and 5, respectively. The 25% and 75% percentile values are, respectively, 10 and 28, covering the upper and lower limits of the set. Also, the 45 test cases are more evenly distributed within the set range. For the number of recue material types m, the maximum and minimum values are 115 and 43, respectively. The 25% and 75% percentile values are, respectively, 61 and 86. Though the values are mostly distributed in [61, 86], they still cover most of the preset range. For the estimated rescue time T, the maximum and minimum values are 70 and 4, respectively. The 25% and 75% percentile values are 14 and 49, respectively. For the time variables, the maximum and minimum values are close to the preset upper and lower limits. In addition, the time in the test cases is concentrated between 14 and 49 hours.

The minimum time for a material storage unit to reach a rescue site is assumed to be 0.2 hours, while the longest time is the corresponding predicted rescue time T under each test condition. The demand q_{j} of each type of material is in the range of [1, 2000]. For each material type j and material repository B_{i}, the storage capacity S_{ij} is in the range [100, 1000], and the unit transportation cost C_{i j} is in [1, 50]. As well, the loss L_{j} caused by the lack of the jth material per unit time is in [0, 1]. According to the data and variable ranges given in Table 2, the required variable values corresponding to each test are randomly generated under the constraints (5)–(9). The simulation data are generated as follows:(1)For each test case, the time t_{i} required to send rescue materials from each material storage B_{i} to the rescue site A is generated randomly, while satisfying the side constraint and sorting t_{i} so that , i = 1, 2, …, n.(2)The demand q_{j} of the jth material type is generated randomly. The demand is related to the emergency impact and the rescue time. Especially, the demand is positively correlated with the rescue time T through , .(3)The loss L_{j} caused by the lack of the jth material type in a unit time is generated randomly within the range [0, 1], .(4)The quantity ν_{j} of the jth material type consumed per unit time is randomly generated. Here, the material consumption is considered as a fixed value that is related to the rescue time and the material demand by , .(5)For each material type j and material repository B_{i}, the unit transportation cost C_{ij} is generated randomly within the range [1, 50], i = 1, 2, …, n, .(6)For each material type j and material repository B_{i}, the storage capacity S_{ij} is generated randomly, such that the total amount of all materials in the material storage units meets the rescue needs, i.e., , i = 1, 2, …, n, .
4.2. Evaluation Indices
4.2.1. Hypervolume Index
The hypervolume index is a set measure used to evaluate the convergence and diversity performance of multiobjective optimization evolutionary algorithms [35]. In our model, there are two objective functions. So, when the number of the Pareto solutions is N, the kth solution is represented by , k = 1, 2, …, N. If we set the reference point as , then the hypervolume index value is defined as
The larger this value is, the better the algorithm convergence and diversity are. Inversely, the worse the algorithm convergence and diversity are, the lower the value of the hypervolume index is.
4.2.2. Diversity Metric
To evaluate the diversity of the Paretooptimal region in multiobjective optimization, the diversity metric is typically used [32]. This metric measures the extensibility and nonuniformity of the distribution of the Paretooptimal region. The diversity metric is defined aswhere N is the number of solutions on the Pareto frontier; is the Euclidean distance of the extreme solution; is the Euclidean distance of the obtained boundary solution; is the Euclidean distance between the ith point and the (i + 1)th point in the obtained nondominated set of solutions (i = 1, 2, …, N − 1); is the average of all distances d_{i}. A perfect solution distribution will make all the Euclidean distances equal to and satisfy . Therefore, the better the distribution is, the smaller the diversity metric is. Inversely, a larger diversity metric will lead to a worse solution distribution.
5. Results and Discussion
In order to verify the feasibility of the proposed mathematical model of continuous material allocation as well as the effectiveness of the algorithmic solution of the multiobjective optimization problem, the improved harmony search algorithm (IHS) [36], the MOEA/D algorithm [25], and NSGAII were used to compare and analyze the solutions for 45 test cases, where the iteration count is set to and the population size . For the three aforementioned algorithms, the average results of the hypervolume index, the diversity metric, and the fitness values are shown in Figures 2–4, respectively. The parameters of each group of test solutions are consistent. The parameters of the IHS algorithm are set as follows: the probability of choosing from the memory , the minimum pitch adjustment rate , and the maximum pitch adjustment rate . The parameters of the MOEA/D algorithm are set as follows: the size of the weight neighborhood , the crossover parameter in the differential evolution operator , the parameter for the differential evolution operator , and the distribution index used by the polynomial mutation . The NSGAII parameters are set as follows: the crossover probability , the crossover distribution index , the mutation probability , and the mutation distribution index .
(a)
(b)
Using the IHS algorithm, we found the mean and variance of the hypervolume index to be 1.93 and 0.099, respectively. Based on the MOEA/D algorithm, the mean and the variance of the hypervolume index are 2.18 and 0.137, respectively. For the NSGAII method, the corresponding mean and variance of the hypervolume index are 2.42 and 0.003, respectively. The evaluation of the hypervolume index shows that the best algorithmic solution of the continuous material allocation problem is obtained by NSGAII, followed by the MOEA/D algorithm and then the IHS one (see Figure 2).
Using the IHS method, the mean and variance of the diversity metric values are 1.7 and 0.185, respectively. Based on the MOEA/D method, the corresponding mean and variance of the diversity metric are 1.51 and 0.112, respectively. For the NSGAII method, we obtained the mean and the variance of the diversity metric as 1.23 and 0.109, respectively. The evaluation of the diversity metric shows the same relative performance of the three algorithms, i.e., the NSGAII method outperforms the MOEA/D one, which is in turn better than the IHS algorithm (see Figure 3).
Using the IHS method, we found the fitness mean and variance of the objective f_{1} to be and , respectively. The corresponding values for the objective f_{2} are and , respectively. Using the MOEA/D method, the fitness average and variance of the objective f_{1} were found to be and , respectively. The corresponding values for f_{2} are and , respectively. Using the NSGAII method, we computed the fitness mean and variance of f_{1} to be and , respectively. The corresponding values for f_{2} are and , respectively. The results in Figure 4(b) show fitness values below 15 using the NSGAII method, while the IHS method got three larger fitness values in the cases , , and , respectively. Also, the MOEA/D method resulted in two larger fitness values in the cases , , and , respectively. In the case , , and , the methods of NSGAII, IHS, and MOEA/D got the largest fitness values of 12.7, 166.4, and 166.4, respectively. Due to the limitations on the number of materials and the quantity of each material, the disasterinduced loss is relatively large no matter what algorithm is used. Nevertheless, IHS and MOEA/D are more likely to fall into local optima than NSGAII during the optimization process. Also, it is difficult to minimize the loss in the shortest time. The results show a consistent ranking of the three algorithms where the NSGAII method gives the best objective value, followed by the MOEA/D method and then the IHS one (see Figure 4).
In general, for the given parameters and among the methods considered for solving the multiobjective continuous material allocation problem, the NSGAII method has the best performance in terms of convergence, diversity, and extensiveness, while the IHS method has the worst performance.
By using the NSGAII method, the convergence and distribution hypervolume measure of each group was calculated. For the 45 test cases, the average and variance of the hypervolume index were 2.42 and 0.003, respectively. The average and variance of the hypervolume index for each case are shown in Table 3.

The average, minimum, maximum variance, and minimum variance of the hypervolume index are 2.494, 2.271, 2.12E − 03, and 1.33E − 08, respectively.
Using the NSGAII method, we calculated the diversity metric of each test case. For the 45 test cases, the average and variance of the diversity metric were 1.23 and 0.109, respectively. The average results for ,, and for multiple runs are shown in Table 4.

The maximum values of and are 0.281 and 0.155, respectively. The average difference between and is 0.013, while the maximum absolute difference is 0.155. Test cases with a difference of less than 0.0009 accounted for approximately 57% of the total test cases, while test cases with a difference of less than 0.01 accounted for approximately 80% of the total number of cases.
Analysis of the distribution of the Pareto dominant set shows that this set has sufficient extensibility and uniform distribution for all test cases. The measurement magnitude has nothing to do with the number of rescue material storage units n, the types of rescue materials m, and the estimated time T. That demonstrates the rationality of the material allocation model and the feasibility of the solution method.
Four test cases were selected to analyze the distribution of the Pareto solution set in a single run. The Pareto frontier results are shown in Figure 5. The abscissa is the fitness value of the objective function f_{1} in (1), and the ordinate represents the fitness value of the objective function f_{2} in (2). In order to exhibit particular details legibly and clearly, the target values are reduced by 10^{2}–10^{6}. As shown in Figure 5, the Pareto noninferior solution is represented by a dotted curve, and the Pareto frontier is represented by a broken line.
(a)
(b)
(c)
(d)
According to the objective functions (1) and (2), the emergency rescue material allocation is a multiobjective nonlinear optimization problem. The Pareto frontier distribution demonstrates the validity and reliability of the constraint processing method. The mathematical model described in this paper is a multiobjective optimization model with sufficient materials. Under this condition, the Pareto frontier results show that when the number of rescue storage units n and the number of rescue material types m are fixed, the longer the estimated rescue time T is and the greater the emergency rescue cost and the resulting losses are. When the number of the emergency rescue storage units n and the estimated rescue time T are ascertained, the more the number of emergency rescue material types m is and the more the emergency rescue costs and caused losses are. When the number of the emergency rescue material types m and the estimated rescue time T are constant, the more the number of the emergency rescue storage units n is and the smaller the emergency rescue cost and the resulting losses are. Therefore, in order to minimize the economic costs and rescue losses, we should establish a sufficient number of emergency rescue material repositories and ensure that there are enough material types and quantities. Any point on the Pareto frontier distribution is a noninferior solution of the multiobjective optimization problem, and there is no distinction for such solutions. There is a contradiction between the objective values f_{1} and f_{2}, which stems from the fact that the smaller the emergency rescue cost is, the greater the loss is and the higher the emergency rescue cost is. Therefore, in the emergency rescue operations, decision makers should consider the actual situation, especially when the lack of materials has a little impact on the damage caused by the rescue operations. When decision makers seek to reduce the economic cost of the emergency rescue operations, the satisfactory emergency rescue solution should be selected from the lower right corner in Figure 5. When the economic cost of the emergency rescue has a little impact, then the lack of materials will cause little rescue damage, and the satisfactory solution should be obtained from the upper left corner in Figure 5. When the economic cost of the emergency rescue has the same impact as the loss of materials, a compromise could be made, and the satisfactory solution should be obtained from the middle part in Figure 5.
A noninferior solution of the bottom right corner of the Pareto frontier is taken as an example to analyze the proposed scheme. In this example, there are 30 material storage warehouses that can provide materials of 86 types to the rescue site within the estimated rescue time of 22 hours (). The time distribution after sorting is shown in Figure 6. The abscissa shows the number of each emergency rescue material repository after sorting the times to reach the rescue site. The ordinate is the time for each repository to arrive at the rescue site. The shortest time is 0.4 hours, and the longest time is 20 hours.
The number of material types needed is 86. The distribution of the loss caused by the lack of each material type per unit time is shown in Figure 7, where the abscissa represents the material types, the left ordinate represents the demand of each material type (the black solid line), and the right ordinate represents the loss caused by the lack of each material type (the red dotted line). The least demand is that of the 80th material m_{80}, for which 88 units of supplies are needed. The highest demand is that of the 49th material m_{49}, for which 968 units of materials are required. The loss caused by the lack of materials hits a value of 1 for the 25th material m_{25} and the 56th material m_{56}. The minimum loss of 0.003 is attained for the 77th material.
The NSGAII method is used to solve the established multiobjective optimization problem. Under the above conditions, 234 noninferior solutions are obtained by the NSGAII method. Different test cases may lead to equal objective function values. So, 255 alternative material allocation schemes are obtained. The material allocation scheme corresponding to the selected noninferior solution is shown in Figure 8. Table 5 is a partial result corresponding to the noninferior material allocation scheme.

In Figure 8, the abscissa indicates the number of 30 emergency relief material repositories sorted by time, while the ordinate indicates 86 rescue material types. The quantity of the emergency rescue material is expressed by an indicative color, which has a few variations from dark to light settings. The darker the color is, the more the allocated materials are. Inversely, the lighter the color is, the less the allocated materials are. Because the value of the objective function f_{2} in (2) is related to the time when each material storage unit reaches the rescue site and the material storage units are sorted by that time, the allocation of various material types corresponding to B_{1} in Figure 7 is the most appropriate. According to the demand size (the black solid line in Figure 6) and the loss caused by the lack of unit materials per unit time (the red dotted line in Figure 6), not all materials are dispatched in the first material storage unit B_{1} to ensure that the balance between the two objective values is achieved. From the material allocation scheme in Figure 7, it can be confirmed that if the demand for the 40th material type m_{40} is q_{40} = 660, the loss caused by the lack of this material per unit time is L_{40} = 0.06, and its allocation amounts in B_{1} and B_{2} are 21 and 258, respectively. Therefore, these results demonstrate that the emergency relief material allocation scheme meets the requirements of the established mathematical optimization model.
6. Conclusions
Observing the constraints of continuous supply and demand of timevarying sequences, we establish a multiobjective mathematical optimization model for the continuous allocation of emergency relief materials. We consider the dynamic characteristics of disasters and accidents and the demand continuity for emergency relief materials, in order to reduce the losses and the economic cost of the emergency relief operations. By comparing the convergence, distribution, and diversity of the solutions based on the IHS, MOEA/D, and NSGAII methods of multiobjective optimization, we conclude that the Pareto solution set can be obtained by the NSGAII method with fairness and balance of loss and cost. The solution set has the advantages of higher scalability, uniformity of distribution, and convergence. According to the results of the loss caused by the lack of rescue materials and the impact of the economic costs on the rescue work, the decision makers can select the appropriate noninferior material allocation scheme from the frontier of the Pareto set. Forty five sets of emergency relief material allocation examples were simulated and tested. It is found that the NSGAII method is more effective than its competing algorithms for solving the multiobjective optimization problem of continuous rescue material allocation. For the NSGAII method, the mean and variance of the hypervolume index are 2.42 and 0.003, respectively. The corresponding values for the diversity metric are 1.23 and 0.109, respectively. These results verify the validity of the mathematical model. Through the analysis of the noninferior solutions, the results show that the emergency rescue material allocation scheme agrees with the established mathematical optimization model. The number of storage units of the emergency rescue materials, the types of emergency rescue materials, and the estimated rescue time are independent of the solution results, an observation that demonstrates the feasibility of using the NSGAII method. Thus, the multiobjective mathematical optimization problem of continuous material allocation for emergency relief can be solved by the NSGAII method, and the results can provide decision support for the continuous allocation and planning of emergency relief materials.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This work was supported by the National Key Research and Development Program of China (2017YFC0804301); the National Natural Science Foundation of China (51774168); the Foundation of Liaoning Educational Committee (LJYL002).
References
 L. Zhang, Y. Lin, G. Yang, and H. Chang, “Emergency resources scheduling based on adaptively mutate genetic algorithm,” Computers in Human Behavior, vol. 27, no. 5, pp. 1493–1498, 2011. View at: Publisher Site  Google Scholar
 J. Raikes and G. McBean, “Responsibility and liability in emergency management to natural disasters: a Canadian example,” International Journal of Disaster Risk Reduction, vol. 16, pp. 12–18, 2016. View at: Publisher Site  Google Scholar
 Z. Su, G. Zhang, Y. Liu, F. Yue, and J. Jiang, “Multiple emergency resource allocation for concurrent incidents in natural disasters,” International Journal of Disaster Risk Reduction, vol. 17, pp. 199–212, 2016. View at: Publisher Site  Google Scholar
 Z. Tang and J. Sun, “Multi objective optimization of railway emergency rescue resource allocation and decision,” International Journal of System Assurance Engineering and Management, vol. 9, no. 3, pp. 696–702, 2017. View at: Publisher Site  Google Scholar
 B. S. Manoj and A. H. Baker, “Communication challenges in emergency response,” Communications of the ACM, vol. 50, no. 3, pp. 51–53, 2007. View at: Publisher Site  Google Scholar
 P. Ghasemi, K. KhaliliDamghani, A. Hafezalkotob, and S. Raissi, “Uncertain multiobjective multicommodity multiperiod multivehicle locationallocation model for earthquake evacuation planning,” Applied Mathematics and Computation, vol. 350, pp. 105–132, 2019. View at: Publisher Site  Google Scholar
 F. Fiedrich, F. Gehbauer, and U. Rickers, “Optimized resource allocation for emergency response after earthquake disasters,” Safety Science, vol. 35, no. 1–3, pp. 41–57, 2000. View at: Publisher Site  Google Scholar
 B. Balcik and B. M. Beamon, “Facility location in humanitarian relief,” International Journal of Logistics Research and Applications, vol. 11, no. 2, pp. 101–121, 2008. View at: Publisher Site  Google Scholar
 X. Li, Z. Zhao, X. Zhu, and T. Wyatt, “Covering models and optimization techniques for emergency response facility location and planning: a review,” Mathematical Methods of Operations Research, vol. 74, no. 3, pp. 281–310, 2011. View at: Publisher Site  Google Scholar
 F. Wex, G. Schryen, S. Feuerriegel, and D. Neumann, “Emergency response in natural disaster management: allocation and scheduling of rescue units,” European Journal of Operational Research, vol. 235, no. 3, pp. 697–708, 2014. View at: Publisher Site  Google Scholar
 B. Yang, Z.H. Hu, and J.X. Zhou, “Configuration, deployment, and scheduling models for management and optimization of patrol services,” Mathematical Problems in Engineering, vol. 2015, Article ID 738578, 13 pages, 2015. View at: Publisher Site  Google Scholar
 Y. Zhou, J. Liu, Y. Zhang, and X. Gan, “A multiobjective evolutionary algorithm for multiperiod dynamic emergency resource scheduling problems,” Transportation Research Part E: Logistics and Transportation Review, vol. 99, pp. 77–95, 2017. View at: Publisher Site  Google Scholar
 G. Schryen, G. Rauchecker, and T. Comes, “Resource planning in disaster response,” Business & Information Systems Engineering, vol. 57, no. 4, pp. 243–259, 2015. View at: Publisher Site  Google Scholar
 A. Y. Chen and T.Y. Yu, “Network based temporary facility location for the Emergency Medical Services considering the disaster induced demand and the transportation infrastructure in disaster response,” Transportation Research Part B: Methodological, vol. 91, pp. 408–423, 2016. View at: Publisher Site  Google Scholar
 C. Duhamel, A. C. Santos, D. Brasil, E. Châtelet, and B. Birregah, “Connecting a population dynamic model with a multiperiod locationallocation problem for postdisaster relief operations,” Annals of Operations Research, vol. 247, no. 2, pp. 693–713, 2016. View at: Publisher Site  Google Scholar
 R. Pradhananga, F. Mutlu, S. Pokharel, J. HolguínVeras, and D. Seth, “An integrated resource allocation and distribution model for predisaster planning,” Computers & Industrial Engineering, vol. 91, pp. 229–238, 2016. View at: Publisher Site  Google Scholar
 D. Alem, A. Clark, and A. Moreno, “Stochastic network models for logistics planning in disaster relief,” European Journal of Operational Research, vol. 255, no. 1, pp. 187–206, 2016. View at: Publisher Site  Google Scholar
 I. Sung and T. Lee, “Optimal allocation of emergency medical resources in a mass casualty incident: patient prioritization by column generation,” European Journal of Operational Research, vol. 252, no. 2, pp. 623–634, 2016. View at: Publisher Site  Google Scholar
 X. Bai, “Twostage multiobjective optimization for emergency supplies allocation problem under integrated uncertainty,” Mathematical Problems in Engineering, vol. 2016, Article ID 2823835, 13 pages, 2016. View at: Publisher Site  Google Scholar
 Y. Deng, W. Zhu, J. Tang, and J. Qin, “Solving a twostage stochastic capacitated locationallocation problem with an improved PSO in emergency logistics,” Mathematical Problems in Engineering, vol. 2017, Article ID 6710929, 15 pages, 2017. View at: Publisher Site  Google Scholar
 C. Boonmee, M. Arimura, and T. Asada, “Facility location optimization model for emergency humanitarian logistics,” International Journal of Disaster Risk Reduction, vol. 24, pp. 485–498, 2017. View at: Publisher Site  Google Scholar
 J. Liu, L. Guo, J. Jiang, D. Jiang, and P. Wang, “Emergency material allocation with timevarying supplydemand based on dynamic optimization method for river chemical spills,” Environmental Science and Pollution Research, vol. 25, no. 18, pp. 17343–17353, 2018. View at: Publisher Site  Google Scholar
 Y. Guo, “A multiobjective INLP model of sustainable resource allocation for longrange maritime search and rescue,” Sustainability, vol. 11, no. 3, 2019. View at: Publisher Site  Google Scholar
 M. Mahdavi, M. Fesanghary, and E. Damangir, “An improved harmony search algorithm for solving optimization problems,” Applied Mathematics and Computation, vol. 188, no. 2, pp. 1567–1579, 2007. View at: Publisher Site  Google Scholar
 Z. Qingfu and L. Hui, “MOEA/D: a multiobjective evolutionary algorithm based on decomposition,” IEEE Transactions on Evolutionary Computation, vol. 11, no. 6, pp. 712–731, 2007. View at: Publisher Site  Google Scholar
 K. Mitra, “Genetic algorithms in polymeric material production, design, processing and other applications: a review,” International Materials Reviews, vol. 53, no. 5, pp. 275–297, 2008. View at: Publisher Site  Google Scholar
 S. M. Mousavi, J. Sadeghi, S. T. A. Niaki, and M. Tavana, “A biobjective inventory optimization model under inflation and discount using tuned Paretobased algorithms: NSGAII, NRGA, and MOPSO,” Applied Soft Computing, vol. 43, pp. 57–72, 2016. View at: Publisher Site  Google Scholar
 G. W. Mann and S. Eckels, “Multiobjective heat transfer optimization of 2D helical microfins using NSGAII,” International Journal of Heat and Mass Transfer, vol. 132, pp. 1250–1261, 2019. View at: Publisher Site  Google Scholar
 M. Laszczyk and P. B. Myszkowski, “Improved selection in evolutionary multiobjective optimization of multiskill resourceconstrained project scheduling problem,” Information Sciences, vol. 481, pp. 412–431, 2019. View at: Publisher Site  Google Scholar
 S. H. R. Pasandideh, S. T. A. Niaki, and K. Asadi, “Biobjective optimization of a multiproduct multiperiod threeechelon supply chain problem under uncertain environments: NSGAII and NRGA,” Information Sciences, vol. 292, pp. 57–74, 2015. View at: Publisher Site  Google Scholar
 Z. Sazvar, S. M. J. Mirzapour Alehashem, A. Baboli, and M. R. Akbari Jokar, “A biobjective stochastic programming model for a centralized green supply chain with deteriorating products,” International Journal of Production Economics, vol. 150, pp. 140–154, 2014. View at: Publisher Site  Google Scholar
 K. Deb, A. Pratap, S. Agarwal, and T. Meyarivan, “A fast and elitist multiobjective genetic algorithm: NSGAII,” IEEE Transactions on Evolutionary Computation, vol. 6, no. 2, pp. 182–197, 2002. View at: Publisher Site  Google Scholar
 P. Mittal and K. Mitra, “Determining layout of a wind farm with optimal number of turbines: a decomposition based approach,” Journal of Cleaner Production, vol. 202, pp. 342–359, 2018. View at: Publisher Site  Google Scholar
 K. Deb, “An efficient constraint handling method for genetic algorithms,” Computer Methods in Applied Mechanics and Engineering, vol. 186, no. 24, pp. 311–338, 2000. View at: Publisher Site  Google Scholar
 A. Auger, “Theory of the hypervolume indicator: optimal μdistributions and the choice of the reference point,” in Proceedings of the Tenth ACM SIGEVO Workshop on Foundations of Genetic Algorithms, ACM, Orlando, FL, USA, January 2009. View at: Publisher Site  Google Scholar
 M. NazariHeris, A. F. Babaei, B. MohammadiIvatloo, and S. Asadi, “Improved harmony search algorithm for the solution of nonlinear nonconvex shortterm hydrothermal scheduling,” Energy, vol. 151, pp. 226–237, 2018. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Yonghong Liu 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.