Next Article in Journal
Lane Optimization of Highway Reconstruction and Expansion Work Zone Considering Carbon Dioxide Emission Factors
Previous Article in Journal
Business Models and Sustainability Plans in the FinTech, InsurTech, and PropTech Industry: Evidence from Spain
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Two-Stage Stochastic Programming Model for Emergency Supplies Pre-Position under the Background of Civil-Military Integration

School of Economics and Management, Yanshan University, Qinhuangdao 066004, China
*
Author to whom correspondence should be addressed.
Sustainability 2022, 14(19), 12080; https://doi.org/10.3390/su141912080
Submission received: 6 August 2022 / Revised: 5 September 2022 / Accepted: 20 September 2022 / Published: 24 September 2022

Abstract

:
The pre-positioning of emergency supplies is of great significance for a timely and effective rescue after a disaster. Using the background of civil-military integration in China, this paper puts the military storage facilities into the layout scheme of emergency supplies reserve and builds a two-stage stochastic programming model of emergency supplies location–allocation optimization, the aim of which is to effectively use the reserve resources of both military and civilian sides to reduce the reserve cost. Then, an improved whale optimization algorithm (IWOA) with more strategies is designed to solve the model. The applicability of the model is proved via a real-world case study in Tangshan, China. The case study shows that when the unit storage cost of military storage facilities is less than 1.5 times of that of civilian emergency storage, the military and civilian joint reserve mode can reduce the reserve cost-effectively; decision-makers can set different maximum rescue times according to different preferences to adapt to different emergency decision-making needs.

1. Introduction

In recent years, natural disasters such as earthquakes and floods have occurred frequently across the whole world, seriously threatening the lives and safety of people [1]. These natural disasters are sudden, urgent and highly uncertain, so a common strategy for dealing with natural disasters is to pre-position emergency supplies so that they can be delivered to the affected area in the first place [2]. Pre-positioning of emergency supplies involves a preparation stage and response stage. The core of the preparation stage is to determine the location of emergency reserve facilities and the type and quantity of supplies to be stored. The core of the response stage is to deliver the materials present in emergency facilities to the disaster site in time and at the lowest cost when a disaster occurs. Therefore, it is important to get global optimization of facility location and resource allocation problems under limited resources and costs.
In times of disaster, the military reserve system plays a key role in emergency response. However, when the military and the government participate in the rescue, they often face problems such as resource waste and low efficiency caused by multiple leaderships and repetitive construction of facilities. This is because the core idea of military–civilian integration is the mutual transformation of military and civilian technology and sharing between military resources and civilian resources. Emergency logistics under the background of civil–military integration emphasizes the sharing of emergency logistics resources, such as the joint reserve of emergency supplies and joint construction of emergency facilities, the aim of which is to increase emergency efficiency and reduce logistics costs. Inevitably, this integration makes pre-positioning more difficult.
In this paper, we put the military storage facilities into the layout scheme of emergency in order to build a civil-military emergency pre-position model, which is to determine the locations of military storage facilities and national civilian emergency facilities and the allocations of emergency supplies, such that the pre-positioning cost is decreased. Then, an improved whale optimization algorithm (IWOA) is designed to solve the model. In order to give full play to the good search performance of the whale optimization algorithm, logistic chaos mapping, nonlinear convergence factor and the Cauchy-Gaussian variation are introduced to adjust the global search ability and local search ability of the algorithm. The main contributions of the paper can be summarized as follows:
  • Put the military storage facilities into the layout scheme of the national emergency reserve; the surplus military warehouse is used to store emergency supplies in order to save reserve resources and improve emergency efficiency.
  • Build a two-stage stochastic programming model of emergency supplies pre-position; both the facility location of pre-disaster and supply allocation post-disaster are simultaneously considered to solve the suboptimal problem caused by separate optimization.
  • Design the IWOA to cope with the proposed model; compared with other classical algorithms, experiments show that IWOA has better optimization performance.
The remainder of this paper is organized as follows. Section 2 reviews the closely relevant literature on emergency facilities’ location and resource allocation operations management. Section 3 presents the problem description and formulates a two-stage stochastic programming model to optimize the emergency operations. Section 4 designs the IWOA with more strategies to solve the model. Section 5 implements the model to carry out a case study in Tangshan, China with the illustrations of collected data and computational results. Finally, we conclude the paper in Section 6.

2. Literature Review

At present, the emergency supply pre-positioning literature is mainly divided into two stages: the preparation stage (pre-disaster) and the response stage (post-disaster). The research related to our problem can be categorized into three streams: pre-disaster emergency decision-making [3,4,5,6,7,8,9], post-disaster emergency decision-making [10,11,12,13,14] and the combination of pre-and-post disaster emergency decision-making [15,16,17,18,19].
The pre-disaster emergency decision-making has focused mainly on the location of emergency facilities, the determination of reserves and the inventory problem. It was first used by the U.S. military to study logistics in a potential war using a pre-position model of emergency supplies [3]. Grass et al. [4] proposed a mixed integer linear programming to determine the facilities’ locations considering the uncertainty of disaster occurrence. Paul and MacDonald [5] considered the location and reserve capacity of the emergency distribution centers to prevent facility damage and casualties. Yang et al. [6] formulated a multi-period dynamic distributionally robust pre-positioning model of emergency supplies under uncertainty of demand. For rapid-onset predictable disasters, Jon and Subodha [7] considered the influence of incorporating returns into pre-disaster deployments. Shu et al. [8] proposed a design model for the humanitarian relief supply network under large-scale natural disasters such as floods and hurricanes. Turkes et al. [9] assumed uncertainties in demand, transportation network availability and inventory depletion, and utilized iterative local search techniques to find good locations and inventory allocation.
The post-disaster emergency decision-making includes the location of emergency shelters, the allocation of emergency supplies and the optimization of vehicle routes. Sheu [10] developed a hybrid fuzzy clustering-optimization method for the operation of emergency logistics co-distribution responding to urgent relief demands during the crucial rescue period. Najafi et al. [11] studied the transport of supplies and casualties in emergency response after an earthquake when resource shortages occur. Zhang et al. [12] proposed a collaborative truck-and-drone system to perform the post-disaster assessment task in the humanitarian relief networks. Lu et al. [13] developed a rolling horizon-based prediction and optimization framework for real-time relief distribution after disasters. The objective of the proposed model is to minimize the total time needed to deliver relief supplies to beneficiaries while considering the risk-averse attitude of the decision-makers. Hu et al. [14] established a bi-objective robust model for emergency resource allocation under uncertainty and developed a heuristic particle swarm optimization algorithm to solve the Pareto frontier.
Facility location of pre-disaster and supply allocation post-disaster are two key factors affecting rescue efficiency, which are interrelated and mutually restricted. Therefore, considering operations of pre-disaster and post-disaster, simultaneously, is conducive to realizing the overall coordination and rational layout of emergency supplies, and meeting the maximization of emergency supplies fill rate in the early stage of the disaster. Rawls and Turnquist [15] presented a two-stage stochastic mixed integer programming model under natural disaster scenarios to solve the problems of location and quantity of pre-disaster storage and the distribution of post-disaster materials. Wang et al. [16] explored how mobile phone location data were applied to capture accurate disaster information rapidly and proposed a preparedness-response two-stage scenario-based stochastic programming model with mobile phone location data for integrated pre-positioning and real-time response operation optimization. Aslan and Celik [17] considered the influence of the demand uncertainty for disaster relief items and the vulnerability of roads and facilities on the pre-disaster decision-making of warehouse location and item reservation; moreover, a two-stage stochastic programming model is established. Chen et al. [18] presented a two-stage delivery process model to represent the delivery of post-disaster relief materials through the combination of social donation and on-site emergency procurement and proposed a risk-sharing scheme in which the state and suppliers jointly reserve relief materials. Davis [19] proposed a stochastic programming model to acquire the amount of supply and the cooperative reallocation scheme of reserve network inventory by using the forecast information. An overview of studies on emergency decision-making is shown in Table 1. Whether it is a pre-disaster or post-disaster emergency, there are few quantitative studies on the introduction of the idea of “civil–military integration” into emergency decision-making.
The last important stream related to our work is the methodologies of dealing with large-scale optimization models. For the recent large-scale optimization models, the solving speed and accuracy of the traditional optimal method cannot perform well. Nature provides great inspiration to solve complex optimization problems. Scholars design some optimization algorithms to deal with large-scale optimization problems by simulating the foraging or evolution behavior of organisms in nature. Some widely used heuristic algorithms are particle swarm optimization algorithm, sparrow algorithm, ant colony optimization algorithm and so on [20,21,22]. Australian researchers Mirjalili and Lewis designed a novel swarm intelligence optimization algorithm, namely the whale optimization algorithm (WOA), which is inspired by the hunting mechanism of humpback whales in nature in 2016 [23]. This algorithm simulates the shrinking encircling, spiral updating position and random hunting mechanisms of humpback whale pods. Nowadays, WOA has been widely applied in engineering, clustering, classification, robot path, image processing, network, task scheduling and other fields. Liu and Zhang [24] proposed a differential evolution chaotic whale optimization algorithm and applied the algorithm to the distribution network fault location of IEEE-33 nodes, which improved the effectiveness and accuracy of the distribution network fault location. Yan et al. [25] improved the whale optimization algorithm to solve the 3D path planning of autonomous underwater vehicles based on forward-looking sonar. Bo et al. [26] proposed a multi-stage mine reuse scheduling method based on an improved whale optimization algorithm. Qian et al. [27] testified the efficiency of the improved multi-objective whale optimization algorithm in the field of vehicle crashworthiness. Tair et al. [28] applied a chaotic oppositional WOA with a firefly search to medical diagnostics. Zhang et al. [29] proposed an energy minimization whale optimization algorithm for workflow scheduling in clouds. In terms of facility location and supply allocation, Yan et al. [30] used the improved WOA to solve the multi-objective optimization problem of water resource allocation in Handan, China. Javid et al. [31] used WOA to solve the robust fuzzy mathematical programming model of a closed-loop supply chain network. Jiang et al. [32] combined WOA and particle-swarm optimization algorithms to cope with the site selection problem of municipal solid waste incineration plants based on certain rules. An overview of studies on optimization algorithms is shown in Table 2.
At present, most of the research on civil-military integration is qualitative. Scholars pay more attention to the construction of a national emergency reserve strategy from the perspective of national policy and industrial planning. There are few quantitative studies on civil-military integration and, as far as we know, there is no report on the location–allocation optimization of emergency rescue in the civilian and military joint reserve. Our work is also the first to propose the IWOA for this kind of combinatorial optimization problem. The results indicate that the IWOA yields higher convergence precision and rates.

3. Model Formulation

3.1. Problem Description

The basic idea of the military and civilian joint reserve is that the local government should make overall planning for the layout of the military and civilian joint reserve, and optimize the reserve scale and structure according to the characteristics of public emergencies in different regions. The military uses its idle storage capacity and specialized storage function to reserve, manage and maintain the stored national emergency and disaster relief materials, which can also be applied to war supplies.
There are local civilian emergency storage facility candidates set I, military storage facility candidates set D and affected areas set J. Scenario set S is used to describe the potential disaster scenario, and the demand for supplies is generated according to the scenario. Considering distributive fairness, if the minimum demand in the affected area is not met, there will be a higher penalty cost due to the lack of stock or the transfer of goods from outside. Maximum rescue time is set according to the urgency of disaster relief. This emergency pre-position decision has obvious two-stage characteristics. When determining the location and reserve capacity of civilian emergency facilities and military storage facilities, future disasters are unknown, so the demand of each affected area is random. Therefore, a two-stage stochastic programming model is established. The objective of the first stage is to minimize the sum of costs for facility location and inventory pre-positioning as well as the expected value of the objective function of the second stage in regard to disaster scenarios, so as to determine the location and capacity of local emergency storage facilities and military storage facilities. After the occurrence of random variables, the objective of the second stage is to minimize the sum of total travel costs and penalty cost of unsatisfied demands for deciding the optimal allocation scheme of emergency supplies in the corresponding scenario after the occurrence of an emergency. The decision of the first stage needs to consider the expected loss of the second stage, and the decision of the second stage depends on the facility’s location and capacity decision of the first stage.
To simplify the problem, we consider that the following assumptions are satisfied:
  • The locations and numbers of civilian emergency facility candidates and military storage facility candidates are known.
  • Use the same mode of transportation and the same type of vehicle to transport emergency supplies from the storage facilities to the affected area, and the unit material transport cost depends on the transport time between the two places.
  • Considering the emergency rescue of the early stage of disaster, the supply needs of the affected areas can only be met by emergency storage and military storage facilities, and social donations are not considered.
  • The objects of the reserve are general emergency supplies with low confidentiality, such as supplies, bedding, medicinal materials and general parts.

3.2. Two-Stage Stochastic Programming Model

A two-stage stochastic programming model for the pre-position of emergency supplies considering military storage facilities is established in this section to minimize the total cost of the system. The first stage: when the demand is unknown, determine the facility location and capacity of the civilian emergency facility and the military storage facility. The second stage: when the location and capacity of each storage depot are known, the disaster situation and the demand of the affected area are also known, and the emergency supplies allocation plan of emergency supplies under the corresponding scenario is determined.
min f = m M i I c m x i m + i I k K b i k y i k + d D k K b d k y d k + E S Q x , y , s
The objective function of the first stage (1) includes the fixed cost of locating and operating facilities, storage cost of emergency supplies pre-positioned at civilian emergency facilities and military storage facilities and the expected value of the second stage solution in regard to disaster scenarios E S Q x , y , s .
  k K y i k v k m M x i m V m ¯ ,   k K y d k v k V ¯ d ,       i I , d D , m M
Constraint (2) indicates that the sum of the volume of supplies stored in the civilian emergency facilities and the military storage facilities does not exceed its maximum reserve capacity.
m M x i m 1         i I ,   m M
Constraint (3) ensures each potentially affected area can only be equipped with a maximum of one emergency storage.
m M i I x i m h ,   d D x d g ,         i I ,   d D ,   m M
Constraint (4) indicates that the number of the selected civilian emergency facilities is at most h, and the number of the selected military storage facilities is at most g.
x i 0 , 1 , x d 0 , 1 ,         i I , d D
Constraints (5) indicates that the x i and x d are binary variable.
y i k 0 , y d k 0 ,         i I , d D , k K
Constraint (6) ensures y i k and y d k are non-negative.
Q x , y , s = i I j J k K μ k t i j z i j k s + d D j J k K μ k t d j z d j k s + j J k K ω k l j k s i I z i j k s + d D z d j k s
The objective function of the second stage problem (7) incorporates the total travel time cost (includes the time of distributing emergency supplies from civilian emergency facilities to affected areas and from military storage facilities to affected areas) and the penalty cost of unsatisfied demand.
j J z i j k s y i k ,     j J z d j k s y d k ,     i I , d D , k K
Constraint (8) limits the total amount of emergency supplies sent from depots to the affected area by the number of emergency supplies pre-positioned at the storage facilities.
s i g n ( z i j k ( s ) ) t i j T ,     s i g n ( z d j k ( s ) )   t d j T ,     i I , d D , j J
Constraint (9) indicates that the travel time should meet the requirements of emergency distribution.
i I z i j k s + d D z d j k s l j k s ,       j J , k K
Constraint (10) shows that the unmet quantity of material demand is nonnegative
z i j k s 0 , z d j k s 0 ,     i I , d D , j J , k K
Constraint (11) ensures z i j k ( s ) and z d j k ( s ) are nonnegative.

3.3. Assessment of Demand in Affected Areas

In different disaster scenarios, there is usually a disaster center with the highest disaster level, and the disaster level of other affected areas decreases according to the distance from the disaster center. The disaster level includes four levels: extremely major (Level I), major (Level II), general (Level III) and no disaster (Level V). It is assumed that other disaster-affected points decrease by one level every 50 km from the disaster center. If the disaster center is located at point j 1 , the disaster level of j 2 can be approximated by the following formula:
ω j 2 ( s ) = min 4 , r o u n d ( ω j 1 ( s ) + o j 1 j 2 50 )
In Formula (12), o j 1 j 2 is the straight-line distance between the two affected areas, and r o u n d (   ) is the rounded whole function.
Assume that half of the population at each Level I affected area needs emergency supplies. Based on the demand quantity of Level I disasters, the demand quantity decreases with the decrease in disaster grade. With each decrease in disaster grade, the affected population decreases by half. The demand for the first emergency supplies in the disaster scenario is as follows:
l j k ( s ) = n j 2 ω ( s ) l k max ,     ω ( s ) 3 0 ,     ω ( s ) > 3  
In Formula (13), n j represents the population in the affected areas, and l k max represents the demand of each person for the kth emergency supplies when a Level I disaster occurs.

4. Materials and Methods

Due to the large scale and many other variables of the two-stage stochastic programming problem, traditional precise solutions such as the branch and bound method and the cut plane method are difficult to obtain better results quickly. In order to solve the optimization model quickly, the IWOA is considered as the algorithm in this paper. WOA is a new swarm intelligence bionic optimization algorithm proposed by Australian scholars Mirjalili and Lewis in 2016 [23]. As a novel heuristic optimization algorithm, WOA has a simple updating mechanism, few adjustment parameters and certain randomness. Compared with other heuristic optimization algorithms, WOA has higher search ability and search speed, and shows good performance in facility location and allocation problems [30,31,33,34]. In order to prevent the algorithm from falling into local optimum, logistic chaos mapping, nonlinear convergence factor and the Cauchy-Gaussian variation are introduced to adjust the global and local searching ability of the algorithm, and give full play to the good searching performance of WOA.

4.1. Whale Optimization Algorithm

This algorithm is inspired by the hunting mechanism of humpback whales in nature and simulates the shrinking encircling, spiral updating position and random hunting mechanisms of humpback whale pods. This model consists of the following three stages: encircling prey, bubbling-net attacking and search for prey.
  • Encircling Prey
WOA simulates that humpback whales hunt by surrounding the contracted prey. Assuming that the current candidate solution is the target prey, other whales update their positions towards the target location, and their behaviors can be expressed as follows:
D ˜ = C X t ¯ X t ¯
X t ¯ + 1 = X t ¯ A D ˜
A = 2 a × r a
C = 2 r
where t ¯ represents the current iteration, X t ¯ is the position vector of the individual whale at iteration t ¯ , X t ¯ is the position vector of the optimal solution at iteration t ¯ , D ˜ represents the distance between the above two, A and C are coefficient vectors, r is a random vector in [0, 1] and a is a control parameter that is linearly decreased from 2 to 0 over the course of iterations.
2.
Bubbling-Net Attacking
As whales catch their prey, they shrink, encircle and spiral updating position. The shrinking encircling mechanism is mainly achieved by decreasing the value of control parameter a. For the spiral updating position, a logarithmic spiral equation is created between the position of the whales and the current best solution to simulate the spiral movement of humpback whales. Assume that the shrinking encircling mechanism and the spiral position mechanism have the same renewal probability during optimization ( P = 0.5 ) . The model can be expressed as
X t ¯ + 1 = X t ¯ A D ˜                                         P < 0.5 D e b ˙ l cos ( 2 π l ) + X t ¯           P 0.5
where D = X t ¯ X t ¯ , p is a random probability in [0, 1], l is a random number in [−1, 1] and b ˙ is a constant for defining the shape of the logarithmic spiral.
When A < 1 , the whale’s next position can be anywhere between its restricted position and the position of its prey.
3.
Search for Prey
In the search for prey, when A 1 is set, the WOA selects random whale individuals instead of the current optimal individuals as target prey, so that they can search in the global scope. The mathematical model is as follows:
D ˜ = C X r a n d X t ¯
X t ¯ + 1 = X r a n d A D ˜
where X r a n d is a position vector randomly selected from the current whale population.

4.2. Improved Whale Optimization Algorithm

In order to give full play to the good search performance of WOA, an IWOA is proposed in this section. Logistic chaos mapping, nonlinear convergence factor and the Cauchy-Gaussian variation are introduced to adjust the global search ability and local search ability of the algorithm. The detailed steps are as follows:
  • Logistic Chaos Mapping
The initial position of the individual population is very important to the optimization performance of the swarm intelligence algorithm, so the logistic chaotic mapping is introduced to initialize the population, so that the population is evenly distributed, and the convergence speed and optimization accuracy of the algorithm are improved [35]. The formula is as follows:
X t ¯ + 1 = α X t ¯ 1 X t ¯
where α is a control parameter, α = 4 .
2.
Nonlinear Convergence Factor
The exploration and development capability of WOA largely depends on the change of convergence factor a . Since the original convergence factor decreases linearly from 2 to 0, the convergence accuracy of the algorithm will be seriously affected in the late iteration and the algorithm will easily fall into local optimal. In order to balance the relationship between global search and local search of the algorithm, the nonlinear improved convergence factor method is introduced to IWOA [36]. The improved convergence factor is:
a = a i n i t i a l a f i n a l + 1 t ¯ / T max 1 μ t ¯ / T max
where a i n i t i a l and a f i n a l are the start value and end value of a , respectively, μ is the linear weight factor, μ = 25.
3.
The Cauchy-Gaussian Variation
At the late iterate stage of WOA, local optimal stagnation tends to occur. In order to solve this problem, the Cauchy-Gauss mutation strategy is adopted to select the individual with the best fitness at present for mutation, and then compare the position before and after mutation to select the better position for the next iteration [37]. The formula is as follows:
X n e w b e s t t = X best t 1 + λ 1 C a u c h y 0 , σ 2 + λ 2 G a u s s 0 , σ 2
σ = 1 , f X b e s t < f X i exp f X b e s t f X i f X b e s t , o t h e r w i s e
where X newbest t represents the position of the optimal individual after mutation, σ 2 represents the standard deviation of the Cauchy-Gaussian variation strategy, C a u c h y 0 , σ 2 is a random variable that satisfies a Cauchy distribution, G a u s s 0 , σ 2 is a random variable that satisfies a Gaussian distribution, λ 1 = 1 t 2 / T max 2 , λ 2 = t 2 / T max 2 , in the optimization process, λ 1 gradually decrease, λ 2 gradually increase, so that the algorithm can jump out of the current stagnation, and coordinate its local development and global exploration ability.

4.3. Solving Steps

Pseudocode of IWOA solution step is as Algorithm 1.
Algorithm 1. Pseudocode of IWOA solution.
Begin
 Set the maximum iteration number and population size, then initialize the population based on logistic chaos mapping according to Equation (21);
 Calculate the fitness value of individual whale and obtain the current optimal solution;
w h i l e   ( t t max )
f o r   i = 1 : N
  Calculate the value of the nonlinear convergence factor a according to Equation (22), update parameter A , C , P , l ;
   i f ( p < 0.5 )
    i f A < 1
    Get a new individual by encircling the prey according to Equation (15);
    e l s e   i f   A 1
    Get a new individual by searching for prey according to Equation (20);
   end if
   e l s e   i f   p 0.5
    Get a new individual by bubbling-net attacking according to Equation (18);
  end if
 end if
  The Cauchy-Gaussian variation is performed according to Equation (23);
  Update the current optimal individual and fitness location;
t = t + 1 ;
end while
end
As shown in Figure 1, the solving process of IWOA is as follows.

5. Case Study

5.1. Background of Study Area

An earthquake is one of the major disasters in China, which is also regarded as one of the most recurrent and harmful events that could cause massive damages and severe economic losses. The eastern Hebei region is located in the seismic belt around the Bohai Sea, which is not only a place where earthquakes and other natural disasters happen frequently but also an important military strategic location. Tangshan city in the east of Hebei Province is selected to conduct the case study. Earthquakes caused by the compression of tectonic plates in these areas may occur unpredictably, which makes the case study more meaningful as it may serve as a decision-supporting tool for such events.

5.2. Model Parameter Determination

This subsection includes two parts: the data and parameters related to (i) potential civilian emergency facilities and military storage facilities, (ii) seismic scenario generation and (iii) affected area needs. These are designed closely to follow the realistic situation, yet should be used for illustration only.

5.2.1. Potential Civilian Emergency Facilities and Military Storage Facilities

We set the central point of each district as the civil emergency facility candidates, which are also the potential affected areas. Therefore, the time distance between the civilian emergency facilities and the locally affected areas is 0 for each region. Four military storage facility candidates were assumed because of confidentiality reasons. Figure 2 shows the locations of potential affected areas and storage facility candidates.
In response to a possible earthquake, a number of civilian emergency facilities will be set up among 14 potential affected areas, and several military storage facilities will be selected from the four local options for military storage facilities of common material. Drinking water and food are considered as two kinds of emergency supplies. It is assumed that the reserve cost of military storage facilities is 10% higher than that of civilian emergency storage, and the related costs are shown in Table 3. The construction costs and capacity of different types of facilities are shown in Table 4 below.

5.2.2. Scene Collection Generation

According to the Hebei Seismological Bureau, from the beginning of historical records to 2021, 44 earthquakes with a magnitude of 5 or above occurred within 100 km of Tangshan. Among them, 37 with a magnitude of 5.0–5.9, 5 with a magnitude of 6.0–6.9, and two with a magnitude of 7.0–7.9. Disasters are classified into three levels according to earthquake magnitude: Level III (7.0–7.9), Level II (6.0–6.9) and Level I (5.0–5.9) (Table 5).
According to the probability of earthquake occurrence of different levels in historical data, several typical disaster scenarios of each level are selected to construct a scenario set of 8 scenarios, as shown in Table 6.

5.2.3. Affected Areas Needs

According to The Statistical Yearbook of Hebei Province 2020, the total population of 14 potential disaster sites at the end of 2020 is shown in Table 7. The actual road network data from Baidu Map was used to calculate the travel distance between 14 potential affected areas. The average transport speed of military and civilian relief vehicles is 80 km/h and 50 km/h respectively, and the maximum rescue time is 2 h because the initial stage of rescue is considered only.
The disaster level and demand of affected areas is determined according to Formulas (12) and (13). The per capita demand in the affected areas is calculated according to the standard that 1.5 L of drinking water and 0.5 kg of food are provided per person per day. The demand in the affected areas under various scenarios is shown in Table 8 below.

5.3. Results Analysis

5.3.1. Model Results Analysis

With MATLAB2018a as the operating platform, the algorithm is executed on Intel(R), Core(TM)i7-6500ucpu, 2.50 GHz, 16.00 GB memory, Windows10 operating system PC. We set the population size as 100 and the maximum number of iterations as 100. By solving the pre-position model of emergency supplies considering the military storage facilities, the optimal site selection and storage sites are obtained, respectively, Qian’an (medium), Zunhua (small), Lunan (medium), Guye (small), 2nd military storage facility and 3rd military storage facility. As shown in Table 9, these six facilities reserve 16.1523 million liters of drinking water and 5.2689 million kilograms of food, with the average utilization rate of the warehouse reaching 84.53%.
When the disaster reaches Levels I and II, drinking water and food can basically meet all the needs for seven days after the disaster. When the disaster reaches Level III, the actual total demand for drinking water and food in Tangshan is 28.24 million liters and 9.41 million kilograms. According to the method in this paper, the optimal reserves of drinking water and food can meet 61% and 64% of the initial demand of disaster respectively. Due to the magnitude of the earthquake, the optimal food storage plan could not meet all the needs of the first seven days of the disaster, but it could still meet all the needs of the first four days of the earthquake. Within four days, emergency supplies can be sent to the affected area through inter-provincial allocation to meet follow-up needs.
Due to space limitations, only the allocation scheme of scenario 5 is shown, as shown in Figure 3 below.
In order to verify the influence of military storage facilities on the results, IWOA with the same parameters is used to acquire the facility location scheme, reserve scheme and distribution scheme when local emergency supplies are stored independently. As shown in Table 10, compared with the local independent reserve mode, the military and civilian joint reserve mode reduce the total cost of the system by 5.02% and the weighted transportation time by 37.45%, and the material satisfaction rate under major disaster scenarios (level III) has little change on the whole. The results show that the plan of military and civilian joint reserve can greatly improve the economy of supplies reserve and the timeliness of supplies distribution, and improve the emergency ability to deal with major emergencies in this region.

5.3.2. Parameter Sensitivity Analysis

  • Maximum rescue time
In the emergency reserve decision-making, if the emergency rescue time is too long, the restriction on the transportation distance will not play a role of constraint; if the emergency rescue time is too small, emergency facilities can only cover a small number of affected areas, which may lead to high reserve cost. In order to analyze the impact of the maximum rescue time limit on the total system cost and weighted transport time, the maximum rescue time was increased from 0.5 to 3. The experimental results are shown in Figure 4. The results show that there is a clear inflection point in the curve when the maximum rescue time is 1 h or 2 h and the total cost of the system changes obviously.
2.
Unit reserve cost of military storage facilities
The military storage facilities can greatly save the construction cost of the national reserve by using the idle storage capacity of the army to reserve emergency supplies. However, due to the high degree of specialization and strong management ability of the military reserve, the unit reserve cost of the military agent storage point is higher than that of the ordinary civilian emergency reserve. The change in unit reserve cost of military storage facilities will directly affect the total cost and total reserve of the system. Analyzing the change of unit reserve costs of the military storage facilities will affect the total costs and total reserve capacity of the system, The reserve cost of the civilian emergency facilities remains unchanged, and the ratio of unit storage cost of military storage facilities to civilian emergency facilities gradually changes from 1 to 2 for testing, The experimental results are shown in Figure 5. The results show that when the unit reserve cost of military storage facilities is less than 1.5 times that of civilian emergency facilities, the optimal result of the model will contain the maximum number of military storage facilities; when the unit reserve cost of military storage facilities is 1.65 times higher than that of civilian emergency facilities, the optimal results of the model will not include military storage facilities.

5.4. Comparison with WOA and Other Algorithms

In order to verify the performance of the IWOA, the genetic algorithm (GA), particle swarm optimization (PSO), simulated annealing (SA) and WOA are employed to solve the model by using the same parameters. The convergence curve in Figure 6 shows that GA, PSO and SA tend to fall into local optimum prematurely for the pre-position model proposed in this paper. IWOA integrates logistic chaotic population initialization, nonlinear convergence factor and improved points of the Cauchy-Gaussian variation, which improves the global search ability of the algorithm and shows better convergence speed and accuracy.
Compared with the results of the five algorithms, the minimum cost is 59.3349 million yuan, 62.4123 million yuan, 620.17 million yuan, 61.7649 million yuan and 61.6366 million yuan, respectively. As shown in Figure 7, compared with other algorithms, IWOA can not only calculate the minimum total cost of the system, but also obtain a better value in the emergency rescue time of the obtained allocation scheme and the emergency supplies satisfaction rate in the case of level III disasters.

6. Managerial Implications

In this section, we will discuss and summarize some interesting managerial insights that emerged from our model formulation and case study analytical results.
(I) We studied the pre-positioning of emergency supplies problem from the perspective of civil-military integration, the aim of which is to effectively make use of the logistics resources of both military and civilian sides to reduce emergency reserve costs. The model of using idle resources of the military to provide storage support for national emergency reserve activities provides a thinking point for managers. In our case study of Tangshan, we first compared the joint reserves model with the independent reserves model. The experimental results show that the joint reserves model obtains a relatively optimal solution under the same parameters. Through the model, it can be verified that the joint reserves model proposed in this paper has practical significance in management practice.
(II) We studied the influence of parameters T on the model results. In order to analyze the impact of the maximum rescue time limit on the total system cost and weighted transport time, the maximum rescue time was increased from 0.5 to 3. The results show that the total cost of the system changes obviously when the maximum rescue time is 1 h and 2 h. If decision-makers are inclined to consider the total cost of pre-disaster storage and response, the maximum rescue time can be 2 h. If decision-makers prefer to consider the timeliness of post-disaster distribution, they can choose the maximum rescue time of 1 h.
(III) Finally, we adjusted the cost of the military storage facilities. The results indicate that when the unit reserve cost of military storage facilities is less than 1.5 times that of civilian emergency facilities, selecting the largest number of military storage facilities can effectively reduce the reserve cost. When the unit reserve cost of military storage facilities is 1.65 times higher than that of civilian emergency facilities, it is more economical to reserve all civilian emergency storage facilities for managers.

7. Discussion

The military has always been the backbone of emergency relief, and the implementation of the military and civilian joint reserve is an important way to effectively solve the problem of emergency relief supplies reserve and improve the ability of the country and the military to deal with natural disasters and emergencies. For the military and civilian general emergency supplies, the military and civilian joint reserve can effectively share the idle storage capacity resources of the army, improve the utilization rate of resource elements and avoid repeated construction. Based on this, this paper proposes a pre-position problem of emergency supplies considering the military storage facilities and builds a two-stage stochastic programming model based on the scenario: the model in the first stage determines the optimal location and the optimal amount of storage, the model in the second stage determines the optimal supply distribution scheme. Aimed at the problem that the initial WOA solution is random and the late iteration is easy to fall into local optimum, the logistic chaotic mapping is introduced to initialize the population and the non-linear convergence factor method and the Cauchy-Gauss mutation strategy are introduced to balance the relationship between global search and local search of the algorithm. To sum up, we developed the IWOA with multi-strategies when searching for solutions to the problem. Based on the case analysis of the Tangshan earthquake disaster scenario, the optimal decision of location selection, supplies storage and post-disaster distribution of military and civilian storage was obtained. The results show that: (1) the total cost of the system can be slightly reduced and the post-disaster transportation time can be effectively optimized by using the military and civilian joint reserve mode for the general emergency supplies. (2) Decision-makers can choose the maximum rescue time according to different emergency decision-making needs. (3) When the unit storage cost of military storage facilities is less than 1.5 times that of civilian emergency storage, the military and civilian joint reserve mode can reduce the reserve cost. (4) For the model in this paper, compared with WOA, the IWOA based on chaotic initialization, nonlinear convergence factor and the Cauchy-Gaussian variation shows better convergence accuracy and convergence speed. In terms of calculation results of the model, IWOA can not only obtain the minimum total cost of the system, but also get a better value in the emergency rescue time of the obtained allocation scheme and the emergency supplies satisfaction rate in the case of level III disasters.
However, the location–allocation model in this paper only considers the civilian emergency relief scenarios and does not consider the military response scenarios. In subsequent studies, both civilian emergency and military response efficiencies can be considered as the objectives of the model to conduct joint optimization. In the optimization of resource deployment under the background of civil–military integration, there are still many new problems that have not been deeply explored, such as the unified scheduling and route planning of military–civilian vehicles in emergency rescue, and the design of military–civilian logistics integration network. These problems can be studied systematically in the future to provide decision support for military and civilian managers.

Author Contributions

Q.L.: conceptualization, software, data curation, methodology, original draft preparation, visualization. J.W.: supervision, writing—reviewing editing, validation. Y.W.: re-viewing, software. J.L.: reviewing, editing English, literature survey and review. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Hebei Province Military and Civil Integration Development Strategy Research and Consulting Project (HBJMRH202104), The Key Research Base of Humanities and Social Sciences of Colleges and Universities of Hebei Province (639000264), Youth Project of Natural Science Foundation of Hebei Province (G2021203003).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data will be available upon request.

Acknowledgments

The authors would like to thank Yanshan University for their support in developing this work.

Conflicts of Interest

The authors declare no conflict of interest.

Glossary

Sets:
I Set of local civilian emergency facility candidates i
D Set of military storage facility candidates d
J Set of affected areas j
S Set of scenarios s
K Set of the types of emergency supplies k
M Set of the level of civilian emergency facility m
Parameters:
p s The probability of the scenario s occurring
h The maximum number of civilian emergency facility
g The maximum number of military storage facility
t i j Unit travel time from civilian emergency facilities i to affected area j
t d j Unit travel time from military storage facilities i to affected area j
c m Fixed cost of locating and operating a level m civilian emergency facility
b i k Unit storage cost of emergency supplies k in civilian emergency facility i
b d k Unit storage cost of emergency supplies k in military storage facilities d
μ k Unit time cost of distribution vehicles of emergency supplies k
ω k Unit penalty cost of unsatisfied demands for emergency supply k
l j k s Demand for emergency supplies k for each affected areas j in scenario s
v k Unit volume of emergency supplies k
V m ¯ Spatial capacity of level m civilian emergency facility i
V d ¯ Spatial capacity of military storage facilities d
T Maximum emergency rescue time
Decision variables:
x i m 1, if a level m civilian emergency facility is selected at location i; 0, otherwise
x d 1, if a military storage facility is selected at location d; 0, otherwise
y i k The quantity of the emergency supplies k pre-positioned at civilian emergency facilities i
y d k The quantity of the emergency supplies k pre-positioned at military storage facilities d
z i j k s The quantity of emergency supplies k delivered from civilian emergency facilities i to affected area j in scenario s
z d j k s The quantity of emergency supplies k delivered from military storage facilities i to affected area j in scenario s

References

  1. Caunhye, A.M.; Nie, X.; Pokharel, S. Optimization models in emergency logistics: A literature review. Socio Econ. Plan. Sci. 2012, 46, 4–13. [Google Scholar] [CrossRef]
  2. Sahebjamnia, N.; Torabi, S.A.; Mansouri, S.A. A hybrid decision support system for managing humanitarian relief chains. Decis. Support Syst. 2017, 95, 12–26. [Google Scholar] [CrossRef]
  3. Johnstone, D.P.; Hill, R.R.; Moore, J.T. Mathematically modeling munitions prepositioning and movements. Math. Comput. Model. 2004, 39, 759–772. [Google Scholar] [CrossRef]
  4. Grass, E.; Fischer, K.; Rams, A. An accelerated L-shaped method for solving two-stage stochastic programs in disaster management. Ann. Oper. Res. 2018, 284, 557–582. [Google Scholar] [CrossRef]
  5. Paul, J.A.; MacDonald, L. Location and capacity allocations decisions to mitigate the impacts of unexpected disasters. Eur. J. Oper. Res. 2016, 251, 252–263. [Google Scholar] [CrossRef]
  6. Yang, M.; Liu, Y.; Yang, G. Multi-period dynamic distributionally robust pre-positioning of emergency supplies under demand uncertainty. Appl. Math. Model. 2021, 89, 1433–1458. [Google Scholar] [CrossRef]
  7. Stauffer, J.M.; Kumar, S. Impact of Incorporating Returns into Pre-Disaster Deployments for Rapid-Onset Predictable Disasters. Prod. Oper. Manag. 2020, 30, 451–474. [Google Scholar] [CrossRef]
  8. Shu, J.; Lv, W.; Na, Q. Humanitarian relief supply network design: Expander graph based approach and a case study of 2013 Flood in Northeast China. Transp. Res. Part E-Logist. Transp. Rev. 2021, 146, 102178. [Google Scholar] [CrossRef]
  9. Turkes, R.; Sorensen, K.; Cuervo, D.P. A matheuristic for the stochastic facility location problem. J. Heuristics 2021, 27, 649–694. [Google Scholar] [CrossRef]
  10. Sheu, J.-B. An emergency logistics distribution approach for quick response to urgent relief demand in disasters. Transp. Res. Part E-Logist. Transp. Rev. 2007, 43, 687–709. [Google Scholar] [CrossRef]
  11. Najafi, M.; Eshghi, K.; Dullaert, W. A multi-objective robust optimization model for logistics planning in the earthquake response phase. Transp. Res. Part E-Logist. Transp. Rev. 2013, 49, 217–249. [Google Scholar] [CrossRef]
  12. Zhang, G.; Zhu, N.; Ma, S.; Xia, J. Humanitarian relief network assessment using collaborative truck-and-drone system. Transp. Res. Part E-Logist. Transp. Rev. 2021, 152, 102417. [Google Scholar] [CrossRef]
  13. Lu, C.-C.; Ying, K.-C.; Chen, H.-J. Real-time relief distribution in the aftermath of disasters—A rolling horizon approach. Transp. Res. Part E-Logist. Transp. Rev. 2016, 93, 1–20. [Google Scholar] [CrossRef]
  14. Hu, C.L.; Liu, X.; Hua, Y.K. A bi-objective robust model for emergency resource allocation under uncertainty. Int. J. Prod. Res. 2016, 54, 7421–7438. [Google Scholar] [CrossRef]
  15. Rawls, C.G.; Turnquist, M.A. Pre-positioning of emergency supplies for disaster response. Transport. Res. B-Meth. 2010, 44, 521–534. [Google Scholar] [CrossRef]
  16. Wang, J.; Cai, J.; Yue, X.; Suresh, N.C. Pre-positioning and real-time disaster response operations: Optimization with mobile phone location data. Transp. Res. Part E-Logist. Transp. Rev. 2021, 150, 102344. [Google Scholar] [CrossRef]
  17. Aslan, E.; Celik, M. Pre-positioning of relief items under road/facility vulnerability with concurrent restoration and relief transportation. IISE Trans. 2019, 51, 847–868. [Google Scholar] [CrossRef]
  18. Chen, J.; Liang, L.; Yao, D.-Q. Pre-positioning of relief inventories for non-profit organizations: A newsvendor approach. Ann. Oper. Res. 2017, 259, 35–63. [Google Scholar] [CrossRef]
  19. Davis, L.B.; Samanlioglu, F.; Qu, X.; Root, S. Inventory planning and coordination in disaster relief efforts. Int. J. Prod. Econ. 2013, 141, 561–573. [Google Scholar] [CrossRef]
  20. Katoch, S.; Chauhan, S.S.; Kumar, V. A review on genetic algorithm: Past, present, and future. Multimed. Tools Appl. 2021, 80, 8091–8126. [Google Scholar] [CrossRef]
  21. Mirjalili, S.; Mirjalili, S.M.; Lewis, A. Grey Wolf Optimizer. Adv. Eng. Softw. 2014, 69, 46–61. [Google Scholar] [CrossRef]
  22. Yang, X.S.; Hossein Gandomi, A. Bat algorithm: A novel approach for global engineering optimization. Eng. Comput. 2012, 29, 464–483. [Google Scholar] [CrossRef]
  23. Mirjalili, S.; Lewis, A. The Whale Optimization Algorithm. Adv. Eng. Softw. 2016, 95, 51–67. [Google Scholar] [CrossRef]
  24. Liu, L.; Zhang, R. Multistrategy Improved Whale Optimization Algorithm and Its Application. Comput. Intell. Neurosci. 2022, 2022, 3418269. [Google Scholar] [CrossRef]
  25. Yan, Z.; Zhang, J.; Zeng, J.; Tang, J. Three-dimensional path planning for autonomous underwater vehicles based on a whale optimization algorithm. Ocean Eng. 2022, 250, 111070. [Google Scholar] [CrossRef]
  26. Bo, L.; Li, Z.; Liu, Y.; Yue, Y.; Zhang, Z.; Wang, Y. Research on Multi-Level Scheduling of Mine Water Reuse Based on Improved Whale Optimization Algorithm. Sensors 2022, 22, 5164. [Google Scholar] [CrossRef] [PubMed]
  27. Qian, L.J.; Yu, L.X.; Huang, Y.Z.; Jiang, P.; Gu, X.G. Improved whale optimization algorithm and its application in vehicle structural crashworthiness. Int. J. Crashworthiness 2022, 27, 1–15. [Google Scholar] [CrossRef]
  28. Tair, M.; Bacanin, N.; Zivkovic, M.; Venkatachalam, K. A Chaotic Oppositional Whale Optimisation Algorithm with Firefly Search for Medical Diagnostics. CMC Comput. Mater. Con. 2022, 72, 959–982. [Google Scholar] [CrossRef]
  29. Zhang, L.X.; Wang, L.; Xiao, M.S.; Wen, Z.C.; Peng, C. EM_WOA: A budget-constrained energy consumption optimization approach for workflow scheduling in clouds. Peer-to-Peer Netw. Appl. 2022, 15, 973–987. [Google Scholar] [CrossRef]
  30. Yan, Z.; Sha, J.; Liu, B.; Tian, W.; Lu, J. An Ameliorative Whale Optimization Algorithm for Multi-Objective Optimal Allocation of Water Resources in Handan, China. Water 2018, 10, 87. [Google Scholar] [CrossRef] [Green Version]
  31. Ghahremani-Nahr, J.; Kian, R.; Sabet, E. A robust fuzzy mathematical programming model for the closed-loop supply chain network design and a whale optimization solution algorithm. Expert Syst. Appl. 2019, 116, 454–471. [Google Scholar] [CrossRef]
  32. Jiang, S.; Li, Z.; Gao, C. Study on site selection of municipal solid waste incineration plant based on swarm optimization algorithm. Waste Manag. Res. 2022, 40, 205–217. [Google Scholar] [CrossRef] [PubMed]
  33. Chen, H.; Xu, Y.; Wang, M.; Zhao, X. A balanced whale optimization algorithm for constrained engineering design problems. Appl. Math. Model. 2019, 71, 45–59. [Google Scholar] [CrossRef]
  34. Pham, Q.-V.; Mirjalili, S.; Kumar, N.; Alazab, M.; Hwang, W.-J. Whale Optimization Algorithm With Applications to Resource Allocation in Wireless Networks. IEEE Trans. Veh. Technol. 2020, 69, 4285–4297. [Google Scholar] [CrossRef]
  35. Ma, Z. Chaotic populations in genetic algorithms. Appl. Soft. Comput. 2012, 12, 2409–2424. [Google Scholar] [CrossRef]
  36. Ding, H.; Wu, Z.; Zhao, L. Whale optimization algorithm based on nonlinear convergence factor and chaotic inertial weight. Concurr. Comput. 2020, 32, e5949. [Google Scholar] [CrossRef]
  37. Wang, W.-C.; Xu, L.; Chau, K.-W.; Xu, D.-M. Yin-Yang firefly algorithm based on dimensionally Cauchy mutation. Expert Syst. Appl. 2020, 150, 113216. [Google Scholar] [CrossRef]
Figure 1. Solving Steps of IWOA.
Figure 1. Solving Steps of IWOA.
Sustainability 14 12080 g001
Figure 2. Location map of potentially affected areas and reserve facility options in Tangshan.
Figure 2. Location map of potentially affected areas and reserve facility options in Tangshan.
Sustainability 14 12080 g002
Figure 3. Supplies distribution options in scenario 5. (The numbers in parentheses indicate the allocation amount of drinking water and food, respectively, and the arrows indicate the distribution relationship rather than the transportation route).
Figure 3. Supplies distribution options in scenario 5. (The numbers in parentheses indicate the allocation amount of drinking water and food, respectively, and the arrows indicate the distribution relationship rather than the transportation route).
Sustainability 14 12080 g003
Figure 4. Maximum rescue time sensitivity.
Figure 4. Maximum rescue time sensitivity.
Sustainability 14 12080 g004
Figure 5. Unit reserve cost sensitivity of military storage facilities.
Figure 5. Unit reserve cost sensitivity of military storage facilities.
Sustainability 14 12080 g005
Figure 6. Convergence curve.
Figure 6. Convergence curve.
Sustainability 14 12080 g006
Figure 7. Comparison results of different algorithms.
Figure 7. Comparison results of different algorithms.
Sustainability 14 12080 g007
Table 1. Overview of studies on emergency decision-making.
Table 1. Overview of studies on emergency decision-making.
TypeAuthorYearMethodTopic
Pre-disaster emergency decision-making[3]2004Integer programmingPosition and configure pre-positioned assets
[4]2018Two-stage stochastic programsDisaster management
[5]2016Stochastic programsLocation and capacity allocations
[6]2021Multi-period dynamic distributionally robust optimizationPre-positioning of emergency supplies
[7]2020Stochastic programsPre-Disaster Deployments
[8]2021Nonlinear integer programmingEmergency facility location and relief supply pre-positioning
[9]2021Linear programmingFacility location
[10]2007Hybrid fuzzy clustering-optimization approachThe operation of emergency logistics co-distribution
Post-disaster emergency decision-making[11]2013Multi-objective robust optimizationLogistics planning in the earthquake response phase
[12]2021Mixed-integer linear programmingHumanitarian relief network assessment
[13]2016Rolling horizon-based frameworkRelief distribution in the aftermath of disasters
[14]2016Bi-objective robust optimizationEmergency resource allocation under uncertainty
Combination of pre-and-post disaster emergency decision-making[15]2010Two-stage stochastic mixed integer programPre-positioning of emergency supplies
[16]2021Two-stage scenario-based stochastic programming modelIntegrated pre-positioning and real-time response operation optimization
[17]2019Two-stage stochastic programmingPre-disaster decisions of warehouse location and item pre-positioning
[18]2017Newsvendor approachPre-positioning of relief inventories
[19]2013Stochastic programmingInventory planning and coordination in disaster relief efforts
Table 2. Overview of studies on optimization algorithm.
Table 2. Overview of studies on optimization algorithm.
AlgorithmAuthorYearTopicImprovement Strategy
Genetic algorithm[20]2021Review of Genetic algorithm-
Grey wolf optimizer[21]2014The first study of GWO-
Bat algorithm[22]2012Global engineering optimization-
Whale optimization algorithm[23]2016The first study of WOA-
[24]2022Layer recognitionChaotic logistic map,
Exploitation and exploration,
Lévy flight mechanism,
Evolutionary population dynamics
[25]2022Three-dimensional path planning of autonomous underwater vehicles.-
[26]2022Mine water reuseOpposition-based learning strategy,
Lévy flight strategy,
Nonlinear convergence factor
[27]2022Deterministic optimization of vehicle structural crashworthinessEvolution operators
[28]2022Medical diagnostics-
[29]2022Workflow scheduling in clouds-
[30]2018Allocation of water resourcesLogistic mapping,
Inertia weighting
[31]2019The closed-loop supply chain management-
[32]2019Constrained engineering design problemsLevy flight,
Chaotic local search
Table 3. Unit volume and associated costs of emergency supplies.
Table 3. Unit volume and associated costs of emergency supplies.
Emergency SuppliesUnit Volume (m3)Unit Storage Cost in Civilian Emergency Facility
(CNY)
Unit Storage Cost in Military Storage Facility
(CNY)
Unit Travel Time (CNY/h)Unit Penalty Cost (CNY)
Drinking water (L)0.0151.121.231.525.4
Food (KG)0.0191.781.961.135.6
Table 4. The construction costs and capacity of storage facilities.
Table 4. The construction costs and capacity of storage facilities.
Storage FacilitiesQuantitativeConstruction Cost
(104 CNY)
Reserve Capacity (104 m3)
Large civilian emergency facilities 650015
Medium civilian emergency facilities30010
Small civilian emergency facilities2005
Military storage facilities 205
Table 5. Frequency of earthquakes occurring within 100 km of Tangshan city.
Table 5. Frequency of earthquakes occurring within 100 km of Tangshan city.
MagnitudeLevelFrequencyProbability
7.0–7.9III20.045
6.0–6.9II50.114
5.0–5.9I370.841
Total 441
Table 6. Definition of disaster scenario.
Table 6. Definition of disaster scenario.
ScenarioDisaster CenterLevelProbability
1LuanzhouI0.022
2FengnanI0.023
3GuyeII0.061
4LuannanII0.053
5KaipingIII0.232
6YutianIII0.361
7LuanzhouIII0.248
Table 7. Population at each potentially affected area.
Table 7. Population at each potentially affected area.
Potential Affected AreaFengrunLubeiQian’anZunhuaYutianFengnanLuanzhou
Population (104)80.0778.4677.6870.766.4955.2552.01
Potential Affected AreaLuannanLaotingQianxiCaofeidianLunanGuyeKaiping
Population (104)50.8538.9236.5635.2133.4231.7927.94
Table 8. Drinking water and food requirements in different situations (104 L, 104 KG).
Table 8. Drinking water and food requirements in different situations (104 L, 104 KG).
Scene1234567
SuppliesWaterFoodWaterFoodWaterFoodWaterFoodWaterFoodWaterFoodWaterFood
Fengrun420.37140.12210.1870.06210.1870.06105.0935.03105.0935.03105.0935.03105.0935.03
Lubei205.9668.65411.92137.31205.9668.65102.9834.33102.9834.33102.9834.3300
Qian’an407.82135.94203.9167.97203.9167.97101.9633.99101.9633.9900101.9633.99
Zunhua185.5961.86185.5961.8692.7930.9392.7930.930092.7930.9300
Yutian174.5458.18174.5458.1887.2729.0987.2729.09000000
Fengnan145.0348.34290.0696.6972.5224.1772.5224.1772.5224.170000
Luanzhou136.5345.51136.5345.51136.5345.51136.5345.5168.2622.750068.2622.75
Luannan266.9688.99133.4844.49133.4844.49133.4844.4966.7422.250066.7422.25
Laoting204.3368.11102.1734.0651.0817.03102.1734.06000051.0817.03
Qianxi95.9731.9995.9731.9995.9731.9947.991647.99160000
Caofeidian92.4330.81184.8561.6246.2115.492.4330.81000000
Lunan175.4658.49175.4658.4987.7329.2487.7329.2443.8614.620043.8614.62
Guye166.955.6383.4527.8283.4527.8283.4527.8241.7213.910041.7213.91
Kaiping146.6948.9146.6948.973.3424.4573.3424.4536.6712.220036.6712.22
Table 9. Facilities location and optimal amount of reserve.
Table 9. Facilities location and optimal amount of reserve.
SuppliesDrinking Water (104 L)Food (104 KG)Warehouse Utilization
Location
Qian’an Medium civilian emergency facilities371.47115.4677.66%
ZunhuaSmall civilian emergency facilities194.341.6474.11%
LunanMedium civilian emergency facilities447.47172.6399.92%
GuyeSmall civilian emergency facilities253.353.1677.21%
Military storage facilities 2208.1185.3294.85%
Military storage facilities 3140.51108.6883.45%
Table 10. Schemes comparison of local independent reserve and joint reserve.
Table 10. Schemes comparison of local independent reserve and joint reserve.
The Reserve TypeThe Number of Location FacilitiesSupplies Satisfaction Rate in Level IIIWeighted Transit Time (h)The Total Cost (104 CNY)
Drinking WaterFood
Independent reserves464.65%59.17%358.806247.01
Joint reserves661.08%64.31%224.425933.49
Rate of change −5.52%+8.69%−37.45%−5.02%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, Q.; Wang, J.; Wang, Y.; Lv, J. A Two-Stage Stochastic Programming Model for Emergency Supplies Pre-Position under the Background of Civil-Military Integration. Sustainability 2022, 14, 12080. https://doi.org/10.3390/su141912080

AMA Style

Li Q, Wang J, Wang Y, Lv J. A Two-Stage Stochastic Programming Model for Emergency Supplies Pre-Position under the Background of Civil-Military Integration. Sustainability. 2022; 14(19):12080. https://doi.org/10.3390/su141912080

Chicago/Turabian Style

Li, Qingwen, Jiuhe Wang, Yinggang Wang, and Jian Lv. 2022. "A Two-Stage Stochastic Programming Model for Emergency Supplies Pre-Position under the Background of Civil-Military Integration" Sustainability 14, no. 19: 12080. https://doi.org/10.3390/su141912080

APA Style

Li, Q., Wang, J., Wang, Y., & Lv, J. (2022). A Two-Stage Stochastic Programming Model for Emergency Supplies Pre-Position under the Background of Civil-Military Integration. Sustainability, 14(19), 12080. https://doi.org/10.3390/su141912080

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop