Next Article in Journal
Simulation Analysis of Artificial Intelligence Technology Diffusion under Market Competition and Policy Incentives Based on Complex Network Evolutionary Game Models
Next Article in Special Issue
The Impact of Digital Technology Innovation on the Supply Chain Position: Micro Evidence from the Chinese New Energy Vehicle Companies
Previous Article in Journal
Research on Collaborative Risk Management Mechanism of Mega Projects: A Tripartite Evolutionary Game Model Considering the Participation of Insurance Institution
Previous Article in Special Issue
Comparative Analysis of Key Success Factors in S&OP: Focusing on Manufacturing and Retail Industry in South Korea
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Facilities Sites Selection Optimization for Food Emergency Logistics to Meet Urgent Demands

School of Economics and Management, Changsha University of Science and Technology, Changsha 410075, China
*
Authors to whom correspondence should be addressed.
Systems 2024, 12(7), 241; https://doi.org/10.3390/systems12070241
Submission received: 6 June 2024 / Revised: 3 July 2024 / Accepted: 4 July 2024 / Published: 5 July 2024
(This article belongs to the Special Issue Multi-criteria Decision Making in Supply Chain Management)

Abstract

:
Effective emergency logistics facility site selection is vital for ensuring prompt and fair food supply during crises. This study tackles the intricate task of choosing optimal sites for emergency food logistics facilities by considering varying urgency levels of needs, uncertain demands, and potential facility interruptions. A novel weighted Mahalanobis distance–gray relational analysis–TOPSIS method is devised to evaluate demand urgency and guide site selection decisions. The proposed location model aims to minimize total cost and unmet demand while integrating discrete scenario strategies to address interruption events. Leveraging the Social Network Search (SNS) algorithm, the model is solved, and its effectiveness is validated through a case study analysis. The results highlight the accuracy of the urgency level determination method in capturing demand characteristics and the model’s provision of an objective and practical framework for formulating rational facility location strategies. This approach holds significant promise for enhancing the promptness and fairness of food supply assurance during emergencies.

1. Introduction

Frequent sudden events worldwide have a profound impact on residents, economies, and social development. Post disasters, there is a surge in food demand in affected regions, requiring swift supply under resource constraints [1,2,3]. Optimizing the spatial layout of emergency logistics facilities is a widely employed and effective strategy to reduce material distribution time, manage costs, and enhance material supply efficiency [4,5,6]. It serves as a crucial practical foundation for emergency logistics system operations and has the potential to significantly enhance the effectiveness of material supply efforts. The optimization process should adhere to the principle of prioritizing urgent demands first. Facilities should be strategically located near areas with high demand urgency, and their capacity should align with the urgency of material requirements. By accounting for demand variations and fully acknowledging the urgency’s impact, the development of scientifically driven site selection strategies for emergency logistics facilities can facilitate the establishment of a robust and efficient emergency logistics system. This, in turn, can enhance allocation efficiency, bolster material supply operations during crises, and ensure reliable support [7,8,9,10,11,12,13,14].
Various research endeavors have been pursued to tackle this challenge. Ye Feng et al. introduced an Emergency Warehouse Location Problem (EWLP) model, prioritizing effectiveness in emergency management over efficiency, to investigate the site selection problem for national emergency warehouses in China [15]. Liu Jin et al. advocated emergency material reserve warehouse site selection based on coverage satisfaction and cost-effectiveness [16]. Mansoor et al. presented a comprehensive bilateral location model considering center and customer utilities to enhance service quality and center density [17]. Naranjo et al. proposed a hierarchical facility location model for multi-category facilities, emphasizing storage capacity and coverage rate to optimize overall coverage satisfaction [18]. Regis-Hernández et al. assert that combining location and allocation factors can strengthen the scientific foundation of location decisions [19]. To capture the intricacies and uncertainties of emergencies, Michael et al. developed a multi-capability P-median facility location model that accounts for the diverse needs of individuals following disasters [20]. Guo et al. addressed the intermittent unavailability of established emergency logistics facilities due to capacity constraints or disaster impact in their research on emergency logistics facility location issues [21]. Several groups of researchers have extensively investigated the effects of uncertain demand on location problems [22,23,24,25]. Several groups of scholars emphasize the significance of transportation and transit costs in location problems [26,27,28,29]. In contrast, other groups of scholars advocate incorporating humanitarian considerations into site selection decisions [30,31,32,33].
These studies have given less emphasis to variations in demand urgency. To tackle this gap, Pamucar et al. introduced the fuzzy FUCOM-D’Bonferroni model to prioritize transportation demand management (TDM) issues [34]. Song et al. introduced a novel material dispatching model that accounts for distinct disaster severity levels, ensuring prioritization and concentration on emergency material dispatching in severely affected areas while also considering equity in emergency material distribution across all levels [35]. Furthermore, Song et al. recommended that scheduling decisions should factor in rescue time window constraints across different urgency levels [36]. Zhang et al. proposed a dynamic optimization model that considers demand urgency, arguing that conventional emergency material supply models focus on meeting quantity requirements but overlook the subjective urgent needs of diverse emergency materials and disaster victims, hindering the fair and rational maximization of limited resources’ utility [37]. For complex location models, Mišković et al. proposed a variable neighborhood search method to tackle practical problems that the CPLEX solver cannot handle [38]. Filippi et al. introduced a kernel search heuristic to improve the effectiveness of model solutions [39]. Presently, research on emergency material dispatch primarily focuses on medical supplies and specialized rescue equipment, often overlooking the importance of food provision. Emergency food plays a crucial role in not only sustaining life but also maintaining social order, facilitating swift responses, preventing diseases, supporting long-term rescue efforts, and catering to vulnerable populations. Unlike medical supplies that have specific usage guidelines, various types of food can address basic survival requirements. Additionally, food presents distinct challenges like perishability and specific transportation needs. Therefore, in situations with limited resources, it is essential to prioritize various food requirements and make strategic decisions regarding facility locations.
This paper explores the problem of food emergency logistics facility location by fully considering the impact of demand urgency and determining the optimal layout and configuration of these facilities. Firstly, it designs a method for assessing demand urgency and establishes a location model aimed at minimizing total cost and unmet demand, with considerations for potential facility disruptions. Secondly, it develops a solution method for the model and verifies it through case analysis. Finally, this paper summarizes the findings, aiming to improve China’s crisis management capabilities and enhance emergency response effectiveness.

2. Model Establishment

2.1. Problem Description and Assumptions

Within an emergency logistics system, emergency logistics facilities primarily comprise emergency logistics centers and distribution centers. The material demand points, representing the disaster-affected areas, can be fulfilled directly by the emergency logistics centers or indirectly through the distribution centers. Emergency logistics centers are divided into permanent and temporary types, with permanent ones established during the prevention and preparation phase. In the aftermath of a sudden incident, the emergency response phase may witness substantial casualties and losses, rendering existing logistics centers inadequate to meet the demand. Hence, it becomes imperative at this juncture to address the two-tier food emergency logistics facility location problem encompassing “logistics centers–distribution centers”. The model integrates uncertain demand and demand urgency to formulate the location model for food emergency logistics facilities. To streamline the model, it focuses solely on domestic logistics and considers a single mode of transportation, be it road, waterway, or air transport. Furthermore, the model operates under the following assumptions:
(1)
The positions of potential emergency logistics facilities and demand points are predetermined.
(2)
Each emergency logistics facility can cater to multiple demand points, and conversely, demand points can be serviced by multiple facilities.
(3)
Factors such as demand, damage rate, and damage cost of various material types are not taken into account.
(4)
The capacity of emergency logistics facilities is restricted, and can be adjusted as needed, in proportion to construction and operational expenses.
Following a sudden disaster, uncertainty parameters may not exhibit randomness due to information scarcity, and the outcomes of diverse disasters can differ. It is plausible that uncertainty parameters cannot be deduced from historical data or conform to a probability distribution. In such instances, the prerequisites for employing stochastic programming are not fulfilled, prompting the creation of a robust optimization model. To tackle uncertainty parameters in the location selection issue, uncertainty sets are employed to represent demand uncertainty. We devise a food emergency logistics facility location model centered on demand urgency. The symbols utilized in the model are elucidated below:
I : a set of demand points, I = 1 , 2 , , | I | .
M : a set of established and candidate logistics centers, M = { 1 , 2 , , | M | } .
N : a set of candidate distribution centers, N = { 1 , 2 , , | N | } .
k : the total categories of emergency food.
q i , q i k ¯ : these, respectively, represent the expected demand of the demand point ( I i ) and the actual demand of the demand point ( I i ) for materials ( k ), k K .
λ i , λ i k : these denote the urgency of emergency demand and the urgency of emergency food ( k ) demand at the demand point ( I i ).
F C M , F C N : the unit construction cost of the logistics center and distribution center, F C M > F C N .
V C M , V C N : the unit operating cost of the logistics center and distribution center, V C M > V C N .
S C : the unit transportation cost refers to the cost of transporting goods per unit distance.
T C : the unit transit cost is the charge for moving each unit of goods from the logistics center to the demand point through the distribution center.
v N max , ν N max : the maximum and minimum construction capacity (i.e., the volume of materials they can handle) of logistics centers.
s m i , s n i , s m n : the distance from the demand point to the logistics center, the distance from the demand point to the distribution center, and the distance from the logistics center to the distribution center.
w M , w N : the decision variable, i.e., the number of emergency logistics centers and distribution centers opened.
x m M { 0 , 1 } , x n N { 0 , 1 } : The decision variable signifies the opening status of a logistics center ( M m ) or distribution center ( N n ). A value of 1 indicates that the center is open, while any other value indicates closure. m M , n N .
v m M , v n N : The decision variable. The capacity of a logistics center and distribution center denotes the volume of goods they can hold.
z m n k , z m i k , z n i k : the decision variable representing the allocation quantity of food supply ( k ) among logistics centers, distribution centers, and demand points.

2.2. Determination of Demand Urgency

Following a sudden event, the urgency of demand for emergency food varies across demand points due to factors like the severity of the disaster, economic conditions, and population demographics. The urgency of demand is influenced by factors such as food importance, production capabilities, transportation challenges, and emergency response efficiency. As a result, demand points exhibit varying levels of urgency for different food types. Therefore, the analysis of demand urgency should consider two aspects: the urgency at demand points and the urgency across food categories. Priority should be given to meeting the needs of high-urgency demand points while also addressing those with lower urgency. This approach allows for centralized coordination by emergency material management centers, rational selection of emergency logistics sites, and efficient resource utilization. It enhances decision-making precision, ensures timely food supply during emergencies with limited resources, reduces management expenses, and boosts rescue effectiveness.
To assess the urgency of food demand at various points, we consider the following indicators: degree of disaster (X1), extent of building or road damage (X2), energy demand density (X3), primary industry GDP (X4), population density (X5), and ratio of elderly and young population (X6). For evaluating the urgency of food demand itself, we choose the following indicators: demand level (Y1), nutritional demand (Y2), demand frequency (Y3), transportation temperature and humidity index (Y4), lead time (Y5), demand price elasticity (Y6), self-sufficiency capacity (Y7), edible days (Y9), food energy (Y9), shelf life (Y10), and convenience of consumption (Y11). To ensure objectivity, we first construct an initial decision matrix and preprocess it. Then, we use the entropy weight method to assign weights. Next, we apply the weighted Mahalanobis distance–gray relational–TOPSIS method to obtain urgency coefficients. Finally, we use hierarchical clustering to determine the urgency levels of demand. The steps are as follows:
Step 1: Interval number type conversion. Liu et al. introduced a method for converting attribute values represented as interval numbers [40].
x i e = x + i e ( x + x )
where i e are expected indicative values, and the general value is 0.5; x between [ x , x + ] x and x + belongs to ; x is less than x + .
Step 2: Normalization of indicator values. Define set n as the evaluation objects, associated with set m as the evaluation indicators, and create an initial evaluation indicator matrix A = ( a i j ) n × m , where a i j represents the j -th indicator of the i -th object, 1 i n , 1 j m . Matrix A is normalized using the min–max method, resulting in matrix B = ( b i j ) n × m .
Step 3: Calculate the information entropy H j :
H j = i = 1 n c i j ln c i j / ln n
where c i j = b i j / i = 1 n b i j represents the specific gravity of the indicator. To prevent calculation errors in cases where it equals zero, the lower limit of the normalization interval is set to 0.002.
Step 4: Use w j = 1 H j / n j = 1 n H j to calculate index entropy weight w j .
Step 5: Calculate the positive ideal solution A + and negative ideal solution A :
A + = { a 1 + , a 2 + , , a m + } = ( max x i j 1 i n | j J + , min x i j 1 i n | j J )
A = { a 1 , a 2 , , a m } = ( min x i j 1 i n | j J + , max x i j 1 i n | j J )
where J + and J are benefit-type and cost-type index sets, respectively.
Step 6: Calculate the dimensionless Markov distance D j + between each evaluation scheme and the positive ideal solution, d i + = d ( a i , A + ) = ( a i A + ) Ω Σ 1 Ω T ( a i A + ) T , D i + = d i + max 1 i n d i + . Similarly, the dimensionless Markov distance from the negative ideal solution D i = d i max 1 i n d i , Ω = diag ( w 1 , w 2 , , w m ) ; Σ is the covariance matrix of matrix A ; and d i + is the Mahalanobis distance with dimension.
Step 7: Calculate the dimensionless gray correlation degree R i + and R i of each scheme with positive and negative ideal solutions.
Firstly, calculate the weighted normalized matrix C = ( ( b i j w j ) n × m ) Τ . Then, calculate the gray correlation coefficient matrix r + = ( r i j + ) n × m of each scheme and positive ideal solution, where r i j + = min i min j | c j + c i j | + ρ max i max j | c j + c i j | | c j + c i j | + ρ max i max j | c j + c i j | . Finally, the gray correlation degree R i + = r i + max 1 i n r i + with the positive ideal solution is calculated. Here, r i + = 1 m j = 1 m r i j + and ρ are the resolution coefficients, and ρ = 0.5 is generally selected; c j + = max ( c j ) ; r i j + represents the gray coefficient of the i scheme and the positive ideal solution with respect to the j index; and r i + is the gray correlation degree with dimension.
Similarly, the dimensionless gray correlation degree between each scheme and the negative ideal solution is R i = r i max 1 i n r i .
Step 8: Combine the Mahalanobis distance and the gray correlation degree. The larger D i and R i + are, the closer they are to the positive ideal scheme; conversely, the larger D i + and R i are, the closer they are to the ideal solution.
H i + = α D i + β R i +
H i = α D i + + β R i
where 1 i n ; α + β = 1 , generally 0.5; H i + and H i are the weighted synthetic values of D i and R i + and the weighted synthetic values of D i + and R i , respectively.
Step 9: Utilize H i * = H i + H i + + H i to compute the urgency value H i * . A higher proximity value indicates a more urgent demand, correlating to a heightened level of demand urgency.
Step 10: Use the hierarchical clustering method H i * to cluster urgent values.

2.3. Construction of Site Selection Model Considering the Urgency of Demand

2.3.1. Location Model

Utilizing the maximum coverage model, a food emergency logistics facility location model, denoted as S A M , is developed with the aim of minimizing total cost and unmet demand while considering demand urgency.
Total construction cost: T F C = m M x m M v m M F C M + n N x n N v n N F C N .
Total operating cost: T V C = m M x m M v m M V C M + n N x n N v n N V C n .
Total transportation cost:
T S C = m M ( i I k K z m i k s m i S C + n N k K z m n k s m n S C ) + n N i I k K z n i k s n i S C .
Total transit cost: T T C = n N i I k K x n N z n i k T C .
( S A M ) :   min f 1 = T F C + T V C + T S C + T T C = m M ( x m M v m M ( F C M + V C M ) + i I k K z m i k s m i S C + n N k K z m n k s m n S C ) + n N ( x n N v n N ( F C N + V C N ) + i I k K ( z n i k s n i S C + x n N z n i k T C ) )
min f 2 = i I k K λ i k ( q i k ¯ m M z m i k x m M n M z n i k x m M ) + i I λ i ( k K q i k ¯ m M k K z m i k x m M n N k K z n i k x n n )
s . t . n N x n N = w N
m M x m M = w M
i I k K z m i k + n N k K z m n k v m M , m M
m M k K z m n k v n N , n N
m M z m n k = i I z n i k , n N
z m n k , z m i k = { 0 ,   x m M = 1 0       , e l s e z n i k = { 0 ,   x n N = 1 0       , e l s e
x m M , x n N { 0 , 1 }
v m M v M min , 0 v n N v N max
i I , m M , n N , k K
Objective function (7) aims to minimize the total cost, encompassing construction, operating, operating, transportation, and transit costs. Objective function (8) focuses on minimizing the material shortage at demand points based on their urgency, considering both total demand shortfall and specific material deficiencies. Constraint (9) restricts the number of allowable distribution centers. Constraint (10) limits the number of new logistics centers that can be established. Constraints (11) and (12) ensure that material supply at logistics and distribution centers, respectively, does not exceed their capacities. Constraint (13) enforces flow conservation at distribution centers. Constraint (14) mandates that supplies are only provided if the facility is operational. Constraints (15) to (17) define variable type limitations.
In the model S A M , the objective function (8) incorporates an uncertain parameter related to uncertain demand q i k ¯ for food supply k . This paper utilizes the robust counterpart optimization method, employing a “box” uncertainty set to characterize the uncertainty level in demand q i k ¯ . It integrates the robust counterpart model introduced by Bertsimas and Sim to convert the objective function (8) into a robust equivalent model.
Assume q i k ¯ [ q i k a i k u i k , q i k + a i k u i k ] , i I , where q i k is the nominal demand of demand point I , a i k = ε i k q i k is the demand perturbation, ε i k is the perturbation ratio, and u i k is the uncertainty factor. In uncertainty set U k = { u k : i I u i k Γ u k , i I , the parameter Γ u k signifies the level of uncertainty within the demand set. This parameter quantifies the conservatism of the constraints, mirroring the decision-maker’s risk preference. A higher value of Γ u k indicates a more conservative approach. Here, Γ u k is a subset of [0, I ], while I denotes the count of demand points, suggesting that not all demand points will encounter fluctuations in material ( k ) requirements. This implies that the demands of up to Γ i k demand points can fluctuate within their intervals, with each demand perturbation being ( Γ i k Γ i k ) q i k ¯ , while the demands for the remaining demand points remain at their nominal level q i k . Under this condition, the resulting robust solution remains viable [41]. Part i I k K λ i k q i k ¯ + i I k K λ i q i k ¯ of objective function (8) containing uncertain demand is sorted out to obtain i I k K ( λ i + λ i k ) q i k ¯ . After applying the method described above, i I k K ( λ i + λ i k ) ( q i k + a i k u i k ) is derived. According to the Bertsimas and Sim equivalent model, the inner maximization problem is denoted as max i I k K ( λ i + λ i k ) a i k u i k , and its constraints are represented as i I u i k Γ u k , i I , k K , 0 u i k 1 [42].
Applying duality theory involves introducing dual variables φ k and γ i k to address this maximization problem, resulting in the following outcome:
min i I k K ( λ i + λ i k ) q i k + k K Γ u k φ k + i I k K γ i k
φ k + γ i k ( λ i + λ i k ) ε i k q i k
φ k , γ i k 0 , i I , k K
Ultimately, Equation (8) is converted into Equation (21).
min f 2 = i I k K λ i k ( q i k m M z m i k x m M n M z n i k x m M ) + i I λ i ( k K q i k m M k K z m i k x m M n N k K z n i k x n n ) + k K Γ u k φ k + i I k K γ i k

2.3.2. Site Selection Model Considering Facility Interruption

In light of potential sudden events, particularly natural disasters, that could damage infrastructure, established emergency facilities may face interruptions. When a facility is disrupted, it may lose partial or complete functionality, leading to unmet demands at various points. This necessitates the construction of new facilities or the reassignment of existing ones to address the situation. However, the temporary facilities may be distant and lack sufficient inventory capacity, causing delays in the supply of emergency food materials and impacting rescue operations. A new assumption is introduced, considering the risk of facility interruptions in both emergency logistics centers and distribution centers, with interruptions only considered in cases of complete functional loss. Discrete random scenarios are employed to depict facility interruption events [43]. In addition to the current site selection model, new variables are introduced based on the characteristics of reliability site selection issues. The symbols’ meanings are as follows:
Ξ : the scenario set for facility disruptions, Ξ = { 1 , 2 , , | Ξ | } , ξ Ξ .
p ξ : the probability of occurrence of scenario ξ , p ξ [ 0 , 1 ] .
y m M ξ , y n N ξ { 0 , 1 } : In scenario ξ , whether logistics center ( M m ) and distribution center ( N n ) is interrupted. If it is 1, it is interrupted.
z m n k ξ , z m i k ξ , z n i k ξ : the decision variable, denoted as, signifies the allocation of food material ( k ) among logistics centers, distribution centers, and demand points within scenario ξ .
Building upon model S A M , a location selection model based on discrete scenarios is developed and denoted as model I S A M .
Total construction cost: T F C = m M x m M v m M F C M + n N x n N v n N F C N
Total operating cost: T V C = ξ Ξ m M p ξ x m M y m M ξ v m M V C M + ξ Ξ n N p ξ x n N y n N ξ v n N V C N
Total transportation cost:
T S C = ξ Ξ m M p ξ ( i I k K z m i k ξ s m i S C + n N k K z m n k ξ s m n S C ) + ξ Ξ n N i I k K p ξ z n i k ξ s n i S C
Total transit cost: T S C = ξ Ξ n N i I k K p ξ x n N y n N ξ z n i k ξ T C
( I S A M ) :   min f 1 = T F C + T V C + T S C + T T C = m M x m M v m M F C M + n N x n N v n N F C N + ξ Ξ n N p ξ x n N y n N ξ v n N V C N + ξ Ξ m M p ξ ( x m M y m M ξ v m M V C M + i I k K z m i k ξ s m i S C + n N k K z m n k ξ s m n S C ) + ξ Ξ n N i I k K p ξ z n i k ξ ( s n i S C + x n N y n N ξ T C )
min f 2 = ξ Ξ i I k K P ξ λ i k ( q i k ¯ m M z m i k ξ x m M y m M ξ n N z n i k ξ x m M y m M ξ ) + ξ Ξ i I P ξ λ i ( k K q i k ¯ m M k K z m i k ξ x m M y m M ξ n N k K z n i k ξ x m M y m M ξ )
s . t . n N x n N = w N
m M x m M = w N
i I k K z m i k ξ + n N k K z n i k ξ v m M , m M , ξ Ξ
m M k K z m n k ξ v n N , n N , ξ Ξ
m M z m n k ξ = i I z n i k ξ , n N , ξ Ξ
m M x m M y m M + n N x n N y n N 1
z m n k ξ = { 0 ,   x m M y m M ξ x n N y n N ξ = 1 0       , e l s e z m i k ξ = { 0 ,   x m M y m M ξ = 1 0       , e l s e z n i k ξ = { 0 ,   x n N y n N ξ = 1 0       , e l s e
x m M , x n N { 0 , 1 }
v m M v M min , 0 v n N v N max
i I , m M , n N , ξ Ξ
Constraint (29) guarantees the availability of at least one operational emergency facility to deliver emergency supplies to the point of need in all scenarios. The remaining objective functions and constraints convey the same meaning as those in the model S A M .
For the part ξ Ξ i I k K p ξ ( λ i + λ i k ) q i k ¯ of the objective function (23) that contains uncertain parameter q i k ¯ , the uncertain set description of “box” is adopted, which is robust and transformed by introducing dual variables. The dual variables θ k and ρ i k are introduced and converted into Equations (34)–(36), which are put into the objective function, and finally the final form (37) of the objective function (23) is obtained.
min ξ Ξ i I k K p ξ ( λ i + λ i k ) q i k + k K Γ u k θ k + i I k K ρ i k
θ k + ρ i k ( λ i + λ i k ) ε i k q i k ξ Ξ p ξ
θ k , ρ i k 0 , i I , k K
min f 2 = ξ Ξ i I k K P ξ λ i k ( q i k m M z m i k ξ x m M y m M ξ n N z n i k ξ x m M y m M ξ ) + ξ Ξ i I P ξ λ i ( k K q i k m M k K z m i k ξ x m M y m M ξ n N k K z n i k ξ x m M y m M ξ ) + k K Γ u k θ k + i I k K ρ i k

3. Location Model Solving

3.1. Multiobjective Processing

Balancing cost and unmet value as dual objectives poses challenges in real rescue operations due to differing units of measurement, making the traditional linear weighting method impractical. This paper adopts a method to address multi-objective decision-making, where f 1 * and f 2 * represent the optimal values for the respective sub-problems, and α denotes the cost preference weight determined by the decision-maker [44,45].
min f = α f 1 f 1 * + ( 1 α ) f 2 f 2 *

3.2. Algorithm Design

This paper employs a novel meta-heuristic algorithm, the Social Network Search (SNS), to tackle the model. The SNS algorithm, introduced by SIAMAK TALATAHARI et al. in 2021, mimics user behaviors in social networks to gain recognition, replicating actions like imitation, dialog, argument, and innovation [46]. These behaviors serve as optimization operators, mirroring how users are influenced and driven to share fresh perspectives. The algorithm accounts for interpersonal relationships, incorporates additional random elements and population diversity, and boasts low computational overhead and rapid convergence speed.
(1)
Encoding
The facilities (logistics center, distribution center) for the supply of materials accepted by the demand point are taken as the code, that is, the emergency logistics center and the corresponding rescue provided by the distribution center are taken as the code. Assuming that there are m demand points, each of which can be supplied by n emergency facilities, a vector m × n dimensional vector (called user individual) is constructed using symbol coding.
X = ( x 1 , x 2 , , x k , , x 2 m )
where the value x k ( 1 k 2 m ) of gene location k represents the number of the logistics center or distribution center that provides materials to the corresponding demand point, and the location k = n ( i 1 ) + j corresponds to the number of the Class j facility that supplies materials to the i demand point.
(2)
SNS algorithm steps
Step 1: Initialize scale ( n U s e r ), maximum number of evaluations ( M a x E v a l ), maximum number of iterations ( M a x E v a l ), and other parameters.
Step 2: Initialize individual users and evaluate them.
Step 3: Users are randomly selected to enter four behavioral modes of imitation, dialog, argument, and innovation to obtain a new solution. Since the encoded solution is discrete, the SNS algorithm needs to reprocess the solution after it is updated, and the processing method is Formula (40), where U B and L B represent the maximum and minimum values of the variable value range, and r o u n d ( ) represents the nearest integer value.
x i   n e w = { L B , x i   n e w < L B r o u n d ( x i   n e w ) , L B x i   n e w U B U B , U B < x i   n e w
Step 4: If the new solution satisfies the condition, it replaces the old solution and updates the function value. Otherwise, return to Step 3.
Step 5: If the termination conditions are not met, which typically involve reaching the maximum number of iterations or attaining the global optimum, the process returns to Step 3.
The algorithm flowchart is depicted in Figure 1.

4. Example Analysis

4.1. Data Description

This case study delves into a region in China impacted by the novel coronavirus pneumonia event in October 2022. It involves 20 demand points, with their food demand urgency data and demand urgency analysis data detailed in Table 1 and Table 2, respectively. Table 2 specifically outlines the food demand data for demand point 14. Assuming logistics centers 1 and 2 are already operational, two additional candidate logistics centers (numbered 3 and 4) and six candidate distribution centers can be added. The food demand at demand points follows a specific distribution, with probabilities of 0.1, 0.15, and 0.19 for three scenarios. The emergency food categories required consist of Material 1, Material 2, and Material 3. The coordinates and relevant demand data for each demand point are provided in Table 3, Table 4 and Table 5, respectively. The unit construction cost for logistics centers and distribution centers is 2 and 1, respectively. The unit operating costs are 0.05 and 0.03, respectively. The unit transportation cost for goods is 0.012, with a unit transfer cost of 0.001. The minimum construction capacity for logistics centers is 800, while for distribution centers it is 200, with a maximum construction capacity of 1000. The facility interruption scenarios are detailed in Table 6.

4.2. Results and Analysis

4.2.1. Model Solution Results

(1)
Solution for demand urgency
Using the above method and data for determining demand urgency, we obtained the numerical values for demand urgency and identified the emergency demand urgency level for each demand point through cluster analysis. Given that China categorizes emergency response to sudden events into four levels, we have similarly classified demand urgency into four levels, as shown in Table 7. For demand point 14, which has the highest demand urgency, detailed evaluations of demand urgency for different food categories are provided in Table 8. Additionally, the demand urgency for each demand point across different food categories is presented in Table 9.
According to Table 7, demand point number 14 has an emergency demand urgency level of Level I. This area experiences severe disasters, high per capita energy demand, weak self-supply capabilities, and high population density. There are eight demand points with Level II urgency levels, characterized by relatively lower disaster severity but high energy demand density and lower primary industry GDP. Demand point number 3, despite having a smaller population density, has a higher proportion of elderly and young populations, resulting in a relatively higher urgency level. This demonstrates that the model fully considers various scenarios. There are six demand points with Level III urgency levels, which have higher disaster severity and higher primary industry GDP, but lower demand density and population density. Lastly, there are five demand points with Level IV urgency levels. These points generally have lower disaster severity, low energy demand density, high primary industry GDP, and low population density. They can somewhat meet their own emergency requirements, justifying their classification as Level IV.
Referring to Table 8, Level I comprises 1 class, Level II includes 5 classes, Level III consists of 12 classes, and Level IV encompasses 5 classes. Through the analysis, it is observed that emergency food items sharing the same demand urgency level exhibit consistent attributes. Class I emergency food stands out as a vital resource crucial for human survival and development, representing a life-sustaining element. Class II emergency food embodies emergency characteristics fully. Class III emergency food exhibits high demand, nutritional requirements, and frequency, albeit with subpar performance, making it suitable for scenarios with specific cooking conditions. Grade IV emergency food primarily consists of less essential items with high interchangeability. Thus, it can be inferred that emergency food items with identical demand urgency levels exhibit uniform attributes.
(2)
Solve the model without considering the facility interruption
According to the above model ( S A M ), the SNS algorithm is used to solve the location problem in the prevention and preparation stage. The user group size ( n U s e r ) is set to 10; the maximum evaluation time ( M a x E v a l ) is 2000; the maximum iteration time ( M a x I t e r ) is 200; the probability of constraint violation ( M a x I t e r ) is 0.1; the stocks of three kinds of materials in the built logistics center 2 and logistics center 3 are, respectively, (1000,700, 400) and (450,400,800); the uncertainty levels of the three different food materials ( Γ u ) were 7, 10, and 8, respectively; the disturbance proportion ( ε ) was 0.1; and the cost target weight ( α ) was 0.3. The solution results are shown in Table 10 and Table 11.
Table 10 shows that the optimal solution involves establishing one new logistics center and three distribution centers, resulting in a minimum comprehensive objective function value of 1.0588. The total cost amounts to 4238.869, with a relatively low unmet demand of 2645.678, indicating economic viability during the rescue phase. In Table 11, it is evident that the optimal solution with the minimum comprehensive objective value designates the new logistics center as number 4, with a capacity of 152. The newly opened distribution centers are numbered 6, 7, and 9, with capacities of 554.2330, 838.6398, and 412.5656, respectively. Figure 2 illustrates the layout of supply points and the supply dynamics between supply points and demand points in this scenario.
(3)
Consider facility interruption model solution
According to the ISAM model, the results obtained are shown in Table 12 and Table 13.
Table 12 illustrates that with two new emergency logistics centers and two distribution centers, the minimum comprehensive target value of 1.2356 is achieved, representing the optimal outcome. At this configuration, the total cost amounts to 1934.6315, with unmet demands of 1146.8449 and dimensionless values of 1.6809 and 1.0447, respectively, where the latter is minimized. These solution outcomes effectively mirror the requirements of emergency material support operations. In Table 13, logistics centers 3 and 4 are selected with capacities set at 137.1711 and 61.1074, while distribution centers 7 and 8 have capacities of 442.7607 and 532.5178, respectively. Figure 3 depicts the location layout and the demand–supply dynamics of various materials under this scenario.

4.2.2. Impact Analysis of Cost Weight

Figure 4 depicts the correlation between cost weight (CW) and the comprehensive target value (CTV) alongside the dimensionless target value (DTV). As CW approaches 0, CTV converges to 0.9764; at CW 0.32, CTV peaks at 1.3469; conversely, as CW nears 1, CTV declines to its lowest point at 0.6829.
Distinct intervals exhibit diverse curve behaviors. Within the range of CW from 0 to 0.32, CTV rises with increasing CW; the dimensionless total cost (C-DTV) fluctuates notably but generally decreases, while the dimensionless unmet demand (D-DTV) stabilizes, narrowing the gap between them.
As CW spans from 0.32 to 0.4, CTV inversely correlates with CW; D-DTV sees a marked upsurge, while C-DTV maintains a downward trajectory at a slower rate. The disparity between them diminishes with rising CW, aligning at a certain threshold before gradually diverging.
In the CW range of 0.4 to 1, CTV sustains its decline with minor oscillations. C-DTV and D-DTV uphold their respective trends, with C-DTV’s descent decelerating as D-DTV consistently surpasses C-DTV. The gap between them widens as CW values increase.
Figure 5 illustrates the impact of facility interruption. For CW values between 0 and 0.2, CTV rises with CW, indicating a gradual decline in solution quality. Within this range, C-DTV decreases notably while D-DTV increases slowly, narrowing the gap between them until they equalize around CW 0.18 before diverging again.
In the CW range of 0.2 to 1, CTV initially rises and then falls, peaking at 0.36 for CW 0.36, representing the least favorable solution, and dropping to its minimum at CW 1. Changes in C-DTV and D-DTV follow similar patterns but with varying rates. C-DTV decreases rapidly with CW, while D-DTV increases. The rate of change accelerates beyond CW 0.72.
Comparing Figure 4 and Figure 5 reveals overall consistency with notable differences at key points. CTV’s trend is similar but peaks earlier in Graph 5. D-DTV exhibits a gradual increase followed by a sharp rise, hitting its threshold earlier than C-DTV and after their equalization, underscoring the importance of considering facility interruption in site selection models for decision-makers to grasp the impact of diverse disasters on target values.
Decision-makers can leverage these trends to select CW, factoring in the emergency response stage’s characteristics. During the response phase emphasizing demand satisfaction, this aspect should be proritized with CW between 0 and 0.72. As the situation improves, focus should be shifted to cost, opting for CW between 0.72 and 1.

4.2.3. Impact Analysis of Uncertainty Level

Figure 6 illustrates the relationship between the uncertainty level (UL) and the Critical Threshold Value (CTV) of three different emergency food supplies. The ranges between the maximum and minimum CTV values under varying ULs for these supplies are 0.1859, 0.1386, and 0.1727, respectively. Category 1 exhibits the most significant fluctuation in CTV, while Category 3 shows the least, with a gradual and steady upward trend. Category 2, however, experiences the most frequent fluctuations. Consequently, decision-makers should enhance the accuracy of their assessments for Category 1 and 2 emergency supplies prior to making decisions.
Figure 7 presents the maximum and minimum comprehensive objective values under different UL variations for categories 1, 2, and 3 as 1.4692, 1.5110, and 1.5368, and 1.2508, 1.2694, and 1.2998, respectively. The differences between the maximum and minimum values are 0.2184, 0.2416, and 0.2370. The analysis of curve changes reveals that Category 1 has the smallest difference between maximum and minimum values but shows relatively high frequency in curve fluctuations. For Category 2, a notable observation is that the CTV is generally lower in the earlier part than in the latter when dividing the graph at UL 10, suggesting that priority should be given to whether UL exceeds half of its upper limit when locating facilities for this type of supply. Category 3, with the poorest minimum comprehensive value, exhibits significant fluctuations between the maximum and minimum values, and the uncertainty in CTV increases as UL rises. Therefore, decision-makers should exercise increased caution when setting the UL for this category.
The sensitivity analysis comparing scenarios with and without facility interruptions highlights significant differences in how different types of supplies are affected, underscoring the importance of developing separate models for facility interruption occurrences.

4.2.4. Influence Analysis of Disturbance Proportion

Figure 8a shows the relationship between the disturbance proportion (DP) and the comprehensive target value (CTV). The disturbance proportion varies from 0 to 2. The change in the target value depends on the range of the disturbance proportion.
The relationship changes with DP values. When DP is between 0 and 1, the CTV rises with the DP. This rise is first rapid and then more gradual. When DP is between 1 and 2, the CTV becomes more sensitive to changes in it. CTV oscillates wildly between highs and lows in this range. So, when DP reaches this range, the decision-maker should carefully consider the parameters.
Figure 8b considers the disruption of emergency facilities. In this context, it shows the trend in the CTV with DP. When DP changes in the interval (0, 1), CTV first increases and then decreases as DP increases. The maximum value in this interval is 1.5121, which occurs when DP is 0.32.
When DP is between 1 and 2, CTV varies erratically as DP rises. The fluctuation frequency is high. CTV is more sensitive to DP in this range. The maximum value is 1.6309 when the weight is 1.58, and the minimum value is 1.3627 when the weight is 0.83. As a result, the minimum and maximum values differ by 0.2689, which is a large value.
Compare the results of the sensitivity analysis. First, look at the results in the 0–1 interval without considering facility outages. Without facility outages, CTV always maintains a upward trend. Second, look at the results with facility outages. CTV first rises and then declines. But the interval between 1 and 2 gives different results. The relationship shows large fluctuations and high frequency with or without outages. This confirms that when DP reaches this level, the parameter should be calculated precisely. Decision-makers should not rely on experience alone.

5. Conclusions

This paper presents several innovative aspects in exploring the location problem of emergency logistics facilities. Firstly, food is essential for human survival, making its demand both necessary and important. This paper treats food as emergency supplies, enriching previous studies that primarily focused on medical supplies and specialized rescue materials, thereby expanding the research scope in this field. Secondly, by recognizing the urgency of food demands, the model incorporates this critical factor and devises an effective method for its determination, enhancing the model’s realism. Thirdly, it thoroughly examines the location problem within the context of potential facility disruptions, establishing location models for various scenarios and enriching the model’s application scope. The main conclusions of this paper are as follows:
(1)
Fully considering demand urgency in food emergency logistics facility location is crucial. The method proposed for determining demand urgency can identify both the numerical value and the urgency level of demand. The results are straightforward and practical. By observing that demand points or foods with the same urgency level share common attributes, decision-making in the location and configuration of emergency logistics facilities can be made more rational.
(2)
The model extensively addresses the decision-making environment involving facility disruptions. Case studies highlight the importance of creating distinct location models and factoring in facility disruptions. The model’s solution algorithm can generate comprehensive objective values for different scenarios and determine the optimal decision plan, encompassing the type, number, location, and capacity of emergency logistics facilities, along with material supply. This capability aids decision-makers in promptly devising effective strategies.
(3)
Analyzing variations in comprehensive objective values under different parameter settings and comparing results across various decision-making scenarios reveal distinct trends and frequencies of change for specific parameters within different intervals. Some parameters exhibit clear thresholds, and their changes vary across decision-making contexts. This insight can guide decision-makers in identifying relevant parameters effectively before making informed decisions.
The methods and models presented in this paper effectively tackle the location problem of emergency logistics facilities for food. However, there are areas that can be improved. The method for assessing demand urgency relied solely on objective weighting, overlooking subjective factors that could enhance future research. Furthermore, the model was constructed on specific assumptions without considering real-world factors such as cargo damage costs and model dynamics. Subsequent research could explore these dimensions to further improve the analysis.

Author Contributions

Conceptualization, X.Z. and Y.C.; methodology, Y.C. and X.Z.; software, Y.C.; validation, Y.C.; formal analysis, Y.C. and X.Z.; investigation, X.Z. and Y.C.; resources, X.Z. and Y.C.; data curation, Y.C.; writing—original draft preparation, Y.C.; writing—review and editing, Y.C. and X.Z.; visualization, Y.C. and L.L.; supervision, X.Z.; project administration, X.Z.; funding acquisition, X.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Social Science Foundation of China (Grant No. 21BGL176).

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding authors.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Ni, H.; Yao, X.; Chu, Y.; Shi, B. Integrated “Capability-Timeliness-Decision” Emergency Response Effectiveness Pre-Assessment Model. China Saf. Sci. J. 2020, 30, 148–156. [Google Scholar]
  2. Zhang, H. Comprehensive Balance in Emergency Management: A New Topic. Chin. Public Adm. 2020, 35, 123–130. [Google Scholar]
  3. Wang, B.; Li, D. Research on Emergency Management Mechanism for Sudden Disasters from the Perspective of Intelligence. CSTPCD 2022, 37, 192–197. [Google Scholar]
  4. Zhou, L.; Chen, X.; Chen, H.; Peng, M. Study on the Characteristics of Disaster Relief Supplies under Emergency Situations—Taking Food Supply during the Wenchuan Earthquake as an Example. Acad. Manag. Rev. 2008, 20, 25–29. [Google Scholar]
  5. Wang, D.; Chen, S. Research on the Emergency Management System for Food Safety Incidents in China and Analysis of Environmental Pollution Cases. J. Food Sci. 2016, 37, 283–289. [Google Scholar]
  6. Zhao, Z. Discussion on the Framework for Evaluating the Effectiveness of Government Emergency Response—An Analysis Based on Investigation Reports of Major Accidents and Disasters since 2010. Leadersh. Sci. 2023, 38, 118–121. [Google Scholar]
  7. Li, Y.; Xie, X.; Xie, Y. Research on Grading of Emergency Ordnance Supplies Based on PPSVM during Emergency Tasks. CSTPCD 2014, 34, 123–125. [Google Scholar]
  8. Zhou, Y.; Chen, N.; Li, Z.; Gong, Y. Optimization Design of Emergency Logistics Network in the Early Stage of Post-Earthquake Rescue Considering Facility Disruption Scenarios. Oper. Res. 2020, 29, 107–112. [Google Scholar]
  9. Zhao, J.; Han, W.; Zheng, W.; Zhao, Y. Urban Emergency Medical Supplies Distribution during Major Public Health Emergencies. J. Transp. Eng. 2020, 20, 168–177. [Google Scholar]
  10. Yan, X.; Hou, H.; Yang, J.; Fang, J. Site Selection and Layout of Material Reserve Based on Emergency Demand Graduation under Large-Scale Earthquake. Sustainability 2021, 13, 1236. [Google Scholar] [CrossRef]
  11. Wang, N.; Zuo, T.; Feng, P. Research on Demand Grading of Military Supplies Emergency Distribution Based on Principal Component Clustering Analysis. J. Acad. Mil. Transp. 2021, 23, 45–50. [Google Scholar]
  12. Wang, W.; Du, Y.; Song, X.; Fan, J.; Wang, Y. Research on Priority of Deep-Sea Emergency Supplies Based on Fuzzy Support Vector Machine. J. Health Saf. Environ. 2022, 22, 3321–3325. [Google Scholar]
  13. Saldanha-da-Gama, F. Facility Location in Logistics and Transportation: An Enduring Relationship. Transp. Res. Part E Logist. Transp. Rev. 2022, 166, 102903. [Google Scholar] [CrossRef]
  14. Che, A.; Li, J.; Chu, F.; Chu, C. Optimizing Emergency Supply Pre-Positioning for Disaster Relief: A Two-Stage Distributionally Robust Approach. Comput. Oper. Res. 2024, 166, 106607. [Google Scholar] [CrossRef]
  15. Ye, F.; Zhao, Q.; Xi, M.; Dessouky, M. Chinese National Emergency Warehouse Location Research based on VNS Algorithm. Electron. Notes Discret. Math. 2015, 47, 61–68. [Google Scholar] [CrossRef]
  16. Liu, J.; Zou, R.; Han, Q.; Wang, W.; Qi, D. Research on the Optimization of Emergency Material Storage Location and Material Allocation based on Adaptive Genetic Algorithm. J. Saf. Environ. 2021, 21, 295–302. [Google Scholar]
  17. Mansoor, D.; Jafar, R. Bi-Sided Facility Location Problems: An Efficient Algorithm for k-Centre, k-Median, and Travelling Salesman Problems. Int. J. Syst. Sci. Oper. Logist. 2023, 10, 2235814. [Google Scholar]
  18. Regis-Hernández, F.; Lanzarone, E.; Bélanger, V.; Ruiz, A. Solving jointly districting and resource location and allocation problems: An application to the design of Emergency Medical Services. Comput. Ind. Eng. 2023, 179, 109232. [Google Scholar] [CrossRef]
  19. Naranjo, M.B.; Merino, L.I.M.; Chía, A.M.R. Multi-Product Maximal Covering Second-Level Facility Location Problem. Comput. Ind. Eng. 2024, 189, 109961. [Google Scholar] [CrossRef]
  20. Widener, M.J.; Horner, M.W. A Hierarchical Approach to Modeling Hurricane Disaster Relief Goods Distribution. J. Transp. Geogr. 2010, 19, 821–828. [Google Scholar] [CrossRef]
  21. Guo, Y.; Hu, D.; Zhu, L.; Duan, C. Modeling Research on Location Allocation of Emergency Logistics Facilities Considering Reliability Factors. J. Saf. Sci. Technol. 2017, 13, 85–89. [Google Scholar]
  22. Sarma, D.; Das, A.; Bera, U.K. Uncertain demand estimation with optimization of time and cost using Facebook disaster map in emergency relief operation. Appl. Soft Comput. 2020, 87, 105992. [Google Scholar] [CrossRef]
  23. Gourtani, A.; Nguyen, T.D.; Xu, H. A distributionally robust optimization approach for two-stage facility location problems. EURO J Comput. Optim. 2020, 8, 141–172. [Google Scholar] [CrossRef]
  24. Cheng, C.; Adulyasak, Y.; Rousseau, L.M. Robust facility location under demand uncertainty and facility disruptions. Omega 2021, 103, 102429. [Google Scholar] [CrossRef]
  25. Byrne, T.; Kalcsics, J. Conditional facility location problems with continuous demand and a polygonal barrier. Eur. J. Oper. Res. 2022, 296, 22–43. [Google Scholar] [CrossRef]
  26. Karatas, M.; Yakıcı, E. A multi-objective location analytics model for temporary emergency service center location decisions in disasters. Decis. Anal. 2021, 1, 100004. [Google Scholar] [CrossRef]
  27. Manupati, V.K.; Schoenherr, T.; Wagner, S.M.; Soni, B.; Panigrahi, S.; Ramkumar, M. Convalescent plasma bank facility location-allocation problem for COVID-19. Transp. Res. Part E Logist. Transp. Rev. 2021, 156, 102517. [Google Scholar] [CrossRef]
  28. Shehadeh, K.S. Reducing disparities in transportation distance in a stochastic facility location problem. Transp. Res. Part C Emerg. Technol. 2023, 153, 104199. [Google Scholar] [CrossRef]
  29. Emami, A.; Hazrati, R.; Delshad, M.M.; Pouri, K.; Khasraghi, A.S.; Chobar, A.P. A novel mathematical model for emergency transfer point and facility location. J. Eng. Res. 2024, 12, 182–191. [Google Scholar] [CrossRef]
  30. An, S.; Cui, N.; Bai, Y.; Xie, W.; Chen, M.; Ouyang, Y. Reliable emergency service facility location under facility disruption, en-route congestion and in-facility queuing. Transp. Res Part E Logist. Transp. Rev. 2015, 82, 199–216. [Google Scholar] [CrossRef]
  31. Koca, E. Two-stage stochastic facility location problem with disruptions and restricted shortages. Comput. Ind. Eng. 2023, 183, 109484. [Google Scholar] [CrossRef]
  32. Dönmez, Z.; Kara, B.Y.; Karsu, Ö.; Saldanha-Da-Gama, F. Humanitarian facility location under uncertainty: Critical review and future prospects. Omega 2021, 102, 102393. [Google Scholar] [CrossRef]
  33. Boonmee, C.; Arimura, M.; Asada, T. Facility location optimization model for emergency humanitarian logistics. Int. J. Disaster Risk Reduct. 2017, 24, 485–498. [Google Scholar] [CrossRef]
  34. Pamučar, D.; Deveci, M.; Canıtez, F.; Bozanic, D. A fuzzy Full Consistency Method-Dombi-Bonferroni model for prioritizing transportation demand management measures. Appl. Soft Comput. 2020, 87, 105952. [Google Scholar] [CrossRef]
  35. Song, Y.; Bai, M.; Ma, Y.; Lyu, W.; Huo, F. An Optimal Model for Fair Dispatch of Emergency Materials Considering Regional Disaster Classification. China Saf. Sci. J. 2022, 32, 172–179. [Google Scholar]
  36. Song, Y.; Sun, P. Considering the Urgency of the Disaster Area, the Dispatching Strategy of Emergency Supplies Air Drop Force. Sci. Technol. Eng. 2023, 23, 4011–4018. [Google Scholar]
  37. Zhang, J.; Huang, J.; Wang, T.; Zhao, J. Dynamic Optimization of Emergency Logistics for Major Epidemic Considering Demand Urgency. Systems 2023, 11, 303. [Google Scholar] [CrossRef]
  38. Mišković, S.; Stanimirović, Z.; Grujičić, I. An efficient variable neighborhood search for solving a robust dynamic facility location problem in emergency service network. Electron. Notes Discret. Math. 2015, 47, 261–268. [Google Scholar] [CrossRef]
  39. Filippi, C.; Guastaroba, G.; Huerta-Muñoz, D.; Speranza, M. A kernel search heuristic for a fair facility location problem. Comput. Oper. Res. 2021, 132, 105292. [Google Scholar] [CrossRef]
  40. Liu, Z.; Li, G. Interval-Based Benefit-Cost Multi-Attribute Decision Making Algorithm. Math. Pract. Theory 2021, 51, 93–100. [Google Scholar]
  41. Yu, D.; Gao, L.; Zhao, S. Maximum Coverage Location Model for Emergency Facilities Considering Shared Uncertainties. Oper. Res. Manag. Sci. 2020, 29, 43–50. [Google Scholar]
  42. Bertsimas, D.; Sim, M. The price of robustness. Oper. Res. 2004, 52, 35–53. [Google Scholar] [CrossRef]
  43. Sun, J.; Chen, Z.; Chen, Z.; Li, X. Research on Robust Optimization of Closed-Loop Supply Chain Network based on Prim-DMGA Algorithm. Appl. Res. Comput. 2023, 40, 2984–2992. [Google Scholar]
  44. Karatas, M.; Yakıcı, E. An iterative solution approach to a multi-objective facility location problem. Appl. Soft Comput. 2018, 62, 272–287. [Google Scholar] [CrossRef]
  45. Yan, S.; Qi, J. Study on Location Selection of Multi-Level Emergency Logistics Facilities Considering Demand Uncertainty. Oper. Res. Manag. Sci. 2022, 31, 7–13. [Google Scholar]
  46. Talatahari, S.; Bayzidi, H.; Saraee, M. Social Network Search for Global Optimization. IEEE Access 2021, 9, 92815–92863. [Google Scholar] [CrossRef]
Figure 1. Algorithm flowchart.
Figure 1. Algorithm flowchart.
Systems 12 00241 g001
Figure 2. Location layout and material supply relationship.
Figure 2. Location layout and material supply relationship.
Systems 12 00241 g002
Figure 3. Location layout and material supply relationship.
Figure 3. Location layout and material supply relationship.
Systems 12 00241 g003
Figure 4. Changes in target values for different cost weights.
Figure 4. Changes in target values for different cost weights.
Systems 12 00241 g004
Figure 5. Changes in target values considering different cost target weights for facility outages.
Figure 5. Changes in target values considering different cost target weights for facility outages.
Systems 12 00241 g005
Figure 6. Changes in comprehensive target values at different levels of uncertainty.
Figure 6. Changes in comprehensive target values at different levels of uncertainty.
Systems 12 00241 g006
Figure 7. Changes in comprehensive target values at different levels of uncertainty.
Figure 7. Changes in comprehensive target values at different levels of uncertainty.
Systems 12 00241 g007
Figure 8. Changes in comprehensive target values of different disturbance proportions.
Figure 8. Changes in comprehensive target values of different disturbance proportions.
Systems 12 00241 g008
Table 1. Urgency of food demand at demand points.
Table 1. Urgency of food demand at demand points.
Demand Point NumberCoordinateX1X2X3X4X5X6
1(20,85)28.57%02172.8611.271460.0431.95%
2(5,45)14.29%02182.095.422439.8232.70%
3(42,15)45.45%02172.779.111416.1133.92%
4(38,5)7.14%02196.796.302449.4433.55%
5(95,35)98.00%02213.940.20106.5725.35%
6(85,25)89.00%02220.31084.7625.40%
7(62,80)82.00%02238.932.7054.0322.79%
8(58,75)83.00%02235.850.5093.5822.61%
9(55,85)76.00%02241.315.4046.9327.96%
10(18,80)92.00%02236.3410.246.7726.84%
11(25,30)4.00%02158.2080.87410.8540.48%
12(15,10)5.00%02161.9848.58375.8744.10%
13(45,65)10.27%02201.0338.65248.6229.68%
14(65,20)47.37%02178.8621.742529.3132.19%
15(31,52)32.00%02136.6071.20315.4238.78%
16(2,60)38.46%02158.1623.73336.0744.75%
17(5,5)60.00%02116.6412.47985.8733.70%
18(57,29)8.33%02107.9218.29159.7440.81%
19(4,18)11.76%02158.0639.04202.6441.31%
20(26,35)57.14%02150.4725.96152.1741.93%
Table 2. Analysis data of food demand urgency.
Table 2. Analysis data of food demand urgency.
Broad CategoryMiddle ClassY1Y2Y3Y4X5Y6Y7Y8Y9Y10Y11
grainflour225030.360.50.3270.634.00364.002702
rice225030.692.00.390.474.00120.002702
corn225030.640.50.3057.254.00112.00162
soya bean53010.341.00.430.5933.33390.009002
other225030.300.50.706.354.00325.00452
vegetablefresh340030.561.50.42194.62.5043.4262
dehydrated34030.130.55.344.0425.00273.162254
other340030.580.50.8510.772.5040.0092
fruitfresh327530.461.00.6541.413.6485.97255
dried327.510.330.57.373.5436.36338.322105
nutnut53010.362.05.359.7233.331076.43005
poultry eggegg44530.521.00.721.2522.22147.0020.04
other54530.521.00.8911.8822.22204.8325.04
edible saltedible salt62.530.110.50.540.00400.00.0014601
edible oilvegetable627.530.040.50.6310.6536.36829.383601
animal627.530.042.00.740.0136.36867.573601
sugarsugar65.000.50.131.00.560.00200400.006303
meatlivestock meat4552.50.531.00.797.6918.18236.091.51
poultry4552.50.530.50.851.9018.18284.902.51
processed foodconvenience food225030.430.27.31166.74.00473.003304
dairy product54001.50.562.02.0632.822.5313.2514.05
beverage116080.040.20.6878,5050.6337.154055
processed aquatic products457.50.50.882.54.487.2217.39310.451381
Table 3. Information of candidate logistics centers.
Table 3. Information of candidate logistics centers.
Number of the Built Logistics Centerx-Coordinatey-CoordinateAdded Candidate Logistics Center Numberx-Coordinatey-Coordinate
1558032050
220504438
Table 4. Information of candidate distribution centers.
Table 4. Information of candidate distribution centers.
Candidate Numberx-Coordinatey-CoordinateCandidate Numberx-Coordinatey-Coordinate
5171588021
6257091025
745461023.352
Table 5. Food demand at each demand point.
Table 5. Food demand at each demand point.
Demand Point NumberThe Demand for Different SituationsNominal Demand for Different MaterialsDemand Point NumberThe Demand for Different SituationsNominal Demand for Different Materials
123123123123
12001802205013020111502101701103010
2180240300706050122201901601503040
321020017090309013160230210602080
4220190210506011014210140190806070
5230220190705011015170190210704060
6230240210308012016210210210707070
719018016060706017180240230806040
821026026090507018190220170409060
9210150230509070192301702101005080
10170170150707030202401902206060120
Table 6. Sets of disruption scenarios of logistics facilities.
Table 6. Sets of disruption scenarios of logistics facilities.
SituationOpportunityLogistics Facility Code
12345678910
situation 10.21111111111
situation 20.0510 *11111111
situation 30.11111011111
situation 40.11111111110
* 0 indicates facility interruption. Source of data.
Table 7. Comprehensive parameters for food demand urgency evaluation at demand points.
Table 7. Comprehensive parameters for food demand urgency evaluation at demand points.
Demand Point NumberEvaluate the Comprehensive Parameter ValuesH*SortUrgency Level
D+DR+RH+H
141.8072 2.4434 0.8303 0.6966 1.6369 1.2519 0.5666 1 I
51.7104 1.7563 0.9471 0.7033 1.3517 1.2069 0.5283 2 II
61.7247 1.9923 0.8152 0.8123 1.4038 1.2685 0.5253 3 II
22.1442 2.1346 0.7912 0.7322 1.4629 1.4382 0.5043 4 II
42.1624 2.1447 0.8028 0.7360 1.4738 1.4492 0.5042 5 II
101.4862 1.6830 0.7350 0.9473 1.2090 1.2167 0.4984 6 II
81.8329 1.8761 0.7487 0.8877 1.3124 1.3603 0.4910 7 II
71.8452 1.8531 0.7302 0.8877 1.2917 1.3665 0.4859 8 II
162.3692 2.0854 0.7810 0.7159 1.4332 1.5426 0.4816 9 II
91.9514 1.7347 0.7363 0.9988 1.2355 1.4751 0.4558 10 III
132.5738 1.8462 0.9829 0.8365 1.4146 1.7052 0.4534 11 III
202.7843 1.8132 0.7801 0.7368 1.2966 1.7606 0.4241 12 III
122.6218 1.7583 0.7158 0.8020 1.2371 1.7119 0.4195 13 III
32.7077 1.5804 0.7591 0.6953 1.1697 1.7015 0.4074 14 III
172.7785 1.5384 0.7513 0.7449 1.1449 1.7617 0.3939 15 III
183.0372 1.3996 0.6711 0.8513 1.0353 1.9442 0.3475 16 IV
193.1233 1.3375 0.6887 0.8046 1.0131 1.9639 0.3403 17 IV
13.0537 1.2186 0.7226 0.7451 0.9706 1.8994 0.3382 18 IV
112.9690 1.2881 0.6530 0.8503 0.9706 1.9097 0.3370 19 IV
152.9692 1.2551 0.6589 0.8044 0.9570 1.8868 0.3365 20 IV
Table 8. Comprehensive parameters of urgency evaluation of emergency food demand.
Table 8. Comprehensive parameters of urgency evaluation of emergency food demand.
CategoryEvaluate the Comprehensive Parameter ValuesH*SortUrgency Level
D+DR+RH+H
beverage7.28216.44590.79260.80753.61924.04480.47221I
nut8.73833.87610.74760.82622.31184.78220.32592II
convenience food8.64393.73710.71020.8492.22374.74640.31903II
sugar8.66693.70710.70210.83742.20464.75210.31694II
edible salt8.57583.14730.75760.85121.95254.21350.31675II
dairy product8.76213.56230.73560.83822.1494.80020.30926II
rice9.02113.49560.73130.83022.11344.92570.30027III
flour8.98223.39370.70480.85342.04934.91780.29418III
other vegetables8.85493.3250.69420.87522.00964.86510.29239III
fresh fruit9.08023.34960.72920.84622.03944.96320.291210III
fresh vegetable8.93593.28660.70720.85841.99694.89720.289711III
corn9.00383.22490.70420.86991.96454.93680.284712III
animal fat9.16663.2730.70370.86981.98845.01820.283813III
soya bean9.65823.40510.70080.85162.0535.25490.280914III
vegetable oil9.10863.21610.68160.89311.94895.00080.280415III
other grain9.13213.13910.6970.87111.91815.00160.277216III
other eggs9.33583.1470.70440.86521.92575.10050.274117III
egg9.40113.12570.70860.86061.91725.13080.27218III
dried fruit9.52422.9120.70450.86621.80825.19520.258219IV
dehydrated vegetable9.59482.90740.69230.87051.79985.23270.255920IV
processed aquatic products9.45642.83580.71310.87291.77455.16470.255721IV
poultry9.42322.78130.67910.90231.73025.16270.25122IV
livestock meat9.50732.69170.68370.8951.68775.20110.24523IV
Table 9. Urgency of demand for different food materials.
Table 9. Urgency of demand for different food materials.
Demand PointDifferent Urgency of Food DemandDemand PointDifferent Urgency of Food Demand
Category 1Category 2Category 3Category 1Category 2Category 3
10.65670.67520.5786110.56880.46370.4312
20.65990.64350.5416120.64230.53810.5498
30.5760.51730.579130.57600.57180.5981
40.53850.53740.5971140.51430.53810.5237
50.46930.47200.5122150.51970.44200.5124
60.43910.44380.4685160.46370.46670.4328
70.45990.46210.4610170.68940.59730.1357
80.56750.51730.5578180.56870.56370.5173
90.43680.48390.4764190.6420.61200.6389
100.47120.57130.3751200.53640.53290.6210
Table 10. The model solution results under the number of different logistics facilities.
Table 10. The model solution results under the number of different logistics facilities.
w m w n f f 1 f 2 f 1 * f 2 * f 1 / f 1 * f 2 / f 2 *
011.66852383.1993269.633919.2032569.7402.5927 1.2724
21.78522291.4923620.558926.5952429.3202.4730 1.4904
31.74572661.9053171.0831097.2742180.7022.4259 1.4542
41.29524057.2082586.0002465.9552258.3231.6453 1.1451
51.65903533.5953047.0851476.862266.2432.3926 1.3446
61.76822786.6003183.4051200.0332079.4922.3221 1.5309
111.24473632.9542866.1382585.5182437.3171.4051 1.1759
21.07693993.4072716.5313107.2022750.4831.2852 0.9877
31.05584238.8692645.6783514.8282668.2401.2060 0.9915
41.18324361.4492939.3363368.3382589.0611.2948 1.1353
51.15364583.6032547.3323372.0582390.7431.3593 1.0655
61.23443358.1123091.2012294.0032720.8941.4639 1.1361
211.19864057.6992884.9462950.9362569.1421.3751 1.1229
21.20144560.5392470.4003396.6102165.4191.3427 1.1408
31.09174602.1882596.6783365.6422667.1171.3674 0.9736
41.07565248.0562412.3753629.8122631.0001.4458 0.9169
51.06265285.1082427.9103860.4932607.2001.3690 0.9312
61.31173153.1513094.1612293.7502408.3681.3747 1.2848
Table 11. Solution results of opening facilities and capacity under different number of logistics facilities.
Table 11. Solution results of opening facilities and capacity under different number of logistics facilities.
w m w n Emergency Logistics Facility Number
345678910
01121.5065
21.0163110.6026
3122.8886220.7189298.8664
4313.8970405.2732189.1734313.4505
5574.6851297.219989.475154.50375.3834
6312.6386224.1756185.2262121.0599109.1487121.7259
1123.4350277.6299
2154.0637488.9404234.6332
3152554.2330838.6398412.5656
4117.4460497.70701258.24340.2762129
5100.7196691.0052408.6187760.94480177.2688292.2753
669.6243150.7925493.5004210.1093359.934027.209275.5644
2112158588.3328
212575804.0603844.2429
3263.265713.32539982.3554286.5072159.9360
4200.696536.786731005.940704.9396613.777494.0270
525551159.382778.1607841.4426103.23623.9527
653.168121.2980237.9712222.5209250.5111224.936786.991198.2520
Table 12. Solution results of target values under different number of emergency facilities in the case of facility interruption.
Table 12. Solution results of target values under different number of emergency facilities in the case of facility interruption.
w m w n f f 1 f 2 f 1 * f 2 * f 1 / f 1 * f 2 / f 2 *
011.37411250.90441222.5636635.52311092.09051.96831.1195
21.40311193.61641359.3207632.96881136.26761.88571.1963
31.28361016.30251400.7349738.07471126.35951.37701.2436
41.47381234.83631425.5614633.79221122.11731.94831.2704
51.30401272.70251373.4383819.51631147.16801.55301.1972
61.3275939.70731440.9201702.65691088.86111.33741.3233
111.39962069.31081324.89511425.5944961.91711.45151.3773
21.45562809.16571172.04701195.71951092.81682.34941.0725
31.5308 3009.3252 1246.6189 1251.9059 1077.7028 2.4038 1.1567
41.24092461.42201213.13921689.28331056.46091.45711.1483
51.38962756.37151125.57121298.82361046.42382.12221.0756
61.50192228.57911334.53661179.5523999.03181.88931.3358
212.06172085.08471425.5944546.44251088.25353.81571.3100
21.23561934.63151146.84491150.93691097.75881.68091.0447
31.48203027.56591150.93691372.3786982.33052.20611.1716
43.89821870.87501689.2833201.06791068.37059.30471.5812
52.21472151.82981298.8236485.34371027.82304.43361.2637
61.43342442.53971179.55231310.0228944.63231.86451.2487
Table 13. Solution results of facility establishment and capacity configuration under different number of emergency facilities.
Table 13. Solution results of facility establishment and capacity configuration under different number of emergency facilities.
w m w n Emergency Facility Number
345678910
01327.50960
2250.896686.48520
3187.2122399.767225.6657
4358.6288301.449454.79758107.7734
5243.1211273.8265121.4983136.4630358.3099
694.8141304.828250.9340122.0936130.320363.0702
114.91139.1049
249.530094.5167631.6002
3103.4061942.6161162.2912165.3176
4161.4599306.2215180.289255.6302161.5663
5417.4045138.8492224.1519559.4366197.902886.3463
6183.9976156.3549109.3003201.502487.5131112.0054107.1257
2139.93442.784188.3463
2137.171161.1074442.7607532.5178
3206.611424.2078155.2727875.1702232.9581
47.54151.06956.798517.317724.764012.7616
5171.905769.8868119.522789.8118164.595049.5461208.8138
653.386511.0256118.9475275.9588150.777222.582627.9950123.1602
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zeng, X.; Chen, Y.; Liu, L. Facilities Sites Selection Optimization for Food Emergency Logistics to Meet Urgent Demands. Systems 2024, 12, 241. https://doi.org/10.3390/systems12070241

AMA Style

Zeng X, Chen Y, Liu L. Facilities Sites Selection Optimization for Food Emergency Logistics to Meet Urgent Demands. Systems. 2024; 12(7):241. https://doi.org/10.3390/systems12070241

Chicago/Turabian Style

Zeng, Xiaoqing, Yanping Chen, and Liming Liu. 2024. "Facilities Sites Selection Optimization for Food Emergency Logistics to Meet Urgent Demands" Systems 12, no. 7: 241. https://doi.org/10.3390/systems12070241

APA Style

Zeng, X., Chen, Y., & Liu, L. (2024). Facilities Sites Selection Optimization for Food Emergency Logistics to Meet Urgent Demands. Systems, 12(7), 241. https://doi.org/10.3390/systems12070241

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