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:
: a set of demand points, .
: a set of established and candidate logistics centers, .
: a set of candidate distribution centers, .
: the total categories of emergency food.
, : these, respectively, represent the expected demand of the demand point () and the actual demand of the demand point () for materials (), .
, : these denote the urgency of emergency demand and the urgency of emergency food () demand at the demand point ().
, : the unit construction cost of the logistics center and distribution center, .
, : the unit operating cost of the logistics center and distribution center, .
: the unit transportation cost refers to the cost of transporting goods per unit distance.
: 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.
, : the maximum and minimum construction capacity (i.e., the volume of materials they can handle) of logistics centers.
, , : 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.
, : the decision variable, i.e., the number of emergency logistics centers and distribution centers opened.
, : The decision variable signifies the opening status of a logistics center () or distribution center (). A value of 1 indicates that the center is open, while any other value indicates closure. , .
, : The decision variable. The capacity of a logistics center and distribution center denotes the volume of goods they can hold.
, , : the decision variable representing the allocation quantity of food supply () 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].
where
are expected indicative values, and the general value is 0.5;
between [
,
]
and
belongs to
;
is less than
.
Step 2: Normalization of indicator values. Define set as the evaluation objects, associated with set as the evaluation indicators, and create an initial evaluation indicator matrix , where represents the -th indicator of the -th object, , . Matrix is normalized using the min–max method, resulting in matrix .
Step 3: Calculate the information entropy
:
where
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 to calculate index entropy weight .
Step 5: Calculate the positive ideal solution
and negative ideal solution
:
where
and
are benefit-type and cost-type index sets, respectively.
Step 6: Calculate the dimensionless Markov distance between each evaluation scheme and the positive ideal solution, , . Similarly, the dimensionless Markov distance from the negative ideal solution , ; is the covariance matrix of matrix ; and is the Mahalanobis distance with dimension.
Step 7: Calculate the dimensionless gray correlation degree and of each scheme with positive and negative ideal solutions.
Firstly, calculate the weighted normalized matrix . Then, calculate the gray correlation coefficient matrix of each scheme and positive ideal solution, where . Finally, the gray correlation degree with the positive ideal solution is calculated. Here, and are the resolution coefficients, and is generally selected; ; represents the gray coefficient of the scheme and the positive ideal solution with respect to the index; and is the gray correlation degree with dimension.
Similarly, the dimensionless gray correlation degree between each scheme and the negative ideal solution is .
Step 8: Combine the Mahalanobis distance and the gray correlation degree. The larger
and
are, the closer they are to the positive ideal scheme; conversely, the larger
and
are, the closer they are to the ideal solution.
where
;
, generally 0.5;
and
are the weighted synthetic values of
and
and the weighted synthetic values of
and
, respectively.
Step 9: Utilize to compute the urgency value . A higher proximity value indicates a more urgent demand, correlating to a heightened level of demand urgency.
Step 10: Use the hierarchical clustering method 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 , is developed with the aim of minimizing total cost and unmet demand while considering demand urgency.
Total construction cost: .
Total operating cost: .
Total transportation cost:
Total transit cost:
.
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 , the objective function (8) incorporates an uncertain parameter related to uncertain demand for food supply . This paper utilizes the robust counterpart optimization method, employing a “box” uncertainty set to characterize the uncertainty level in demand . It integrates the robust counterpart model introduced by Bertsimas and Sim to convert the objective function (8) into a robust equivalent model.
Assume
,
, where
is the nominal demand of demand point
,
is the demand perturbation,
is the perturbation ratio, and
is the uncertainty factor. In uncertainty set
, the parameter
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
indicates a more conservative approach. Here,
is a subset of [0,
], while
denotes the count of demand points, suggesting that not all demand points will encounter fluctuations in material (
) requirements. This implies that the demands of up to
demand points can fluctuate within their intervals, with each demand perturbation being
, while the demands for the remaining demand points remain at their nominal level
. Under this condition, the resulting robust solution remains viable [
41]. Part
of objective function (8) containing uncertain demand is sorted out to obtain
. After applying the method described above,
is derived. According to the Bertsimas and Sim equivalent model, the inner maximization problem is denoted as
, and its constraints are represented as
[
42].
Applying duality theory involves introducing dual variables
and
to address this maximization problem, resulting in the following outcome:
Ultimately, Equation (8) is converted into Equation (21).
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, , .
: the probability of occurrence of scenario , .
, : In scenario , whether logistics center () and distribution center () is interrupted. If it is 1, it is interrupted.
, , : the decision variable, denoted as, signifies the allocation of food material () among logistics centers, distribution centers, and demand points within scenario .
Building upon model , a location selection model based on discrete scenarios is developed and denoted as model .
Total construction cost:
Total operating cost:
Total transportation cost:
Total transit cost:
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 .
For the part
of the objective function (23) that contains uncertain parameter
, the uncertain set description of “box” is adopted, which is robust and transformed by introducing dual variables. The dual variables
and
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.