An Optimization Method for Distributing Emergency Materials Which Balances Multiple Decision Criteria

: In emergency rescue operations, it is very important for emergency management to arrange emergency materials scientiﬁcally and reasonably according to emergency needs. However, the decision criteria for emergency material distribution are often diverse, coupled with the lack of information acquisition due to the suddenness and development complexity of disasters, which makes the scientiﬁc decisions about emergency material distribution face great challenges. This paper proposes a multi-criteria decision optimization model of emergency material distribution under uncertainty that combines the 4E criteria of performance evaluation (economic criterion, efﬁciency criterion, effectiveness criterion and equity criterion). Then, the solution method of the multiple decision criteria model is designed. Finally, a case study of emergency medical material distribution during COVID-19 in China was conducted for model veriﬁcation. The results show that different decision criteria have an important impact on the distribution of emergency materials, and the proposed model has signiﬁcant advantages in the allocation of large-scale disaster relief materials, which can provide a useful reference for emergency material distribution under uncertain information conditions.


Introduction
In recent years, various large-scale emergencies, such as earthquakes, floods and infectious diseases, have occurred frequently, causing huge casualties and social and economic losses. Therefore, how to make and choose scientific emergency rescue decisions has attracted great attention [1][2][3][4][5]. The distribution of emergency materials is an important guarantee for emergency response and post-disaster recovery, which has an important impact on the success of the overall emergency relief work [6][7][8][9]. However, there are many decision criteria such as efficiency, equity, economy and effectiveness in the distribution of relief materials in large-scale emergencies. The selection of different decision criteria has an important and different impact on the actual distribution of emergency materials [10,11]. In addition, the information related to the distribution of relief materials in emergency situations is usually uncertain, incomplete or inaccurate [12,13]. During major outbreaks of infectious diseases, for example, in the process of emergency rescue, a large number of emergency materials such as masks, protective clothing, disinfectants and medicines are needed, but during the rescue process, problems such as insufficient initial supply, lagging demand information acquisition, dynamic changes of demand with the epidemic trend, and limited distribution capacity are often faced, which seriously affects the rescue work. Therefore, how to maximize the balance of multiple decision indicators to obtain a scientific and reasonable emergency material allocation plan under the condition of uncertain disaster information is an important practical problem to be solved in the decision field of large-scale emergency rescue. Based on the above background, this study proposes a multi-period distribution optimization model of emergency materials based on multiple decision criteria under uncertain information. The purpose of this study is to provide new ideas and

Literature Review
There is a growing literature that explores emergency material distribution in humanitarian logistics in recent years due to the increased frequency of disasters and the severity of the consequences [14][15][16][17][18][19]. Meanwhile, the issue of different decision criteria/objectives for the distribution of emergency materials has also received more attention, mainly focusing on the efficiency criterion (the shortest total time) [20][21][22], economic criterion (the lowest total cost) [23][24][25][26][27][28], and equity criterion (the lowest total loss) [29,30]. For example, Yan and Shih [20] proposed an optimization model for disaster emergency material allocation with the goal of minimizing the total allocation time, including path transportation time and road repair time. Tang and Ye [22] developed an emergency medical material allocation model facing multiple distribution centers and multiple demand points with the goal of minimizing the total time of loading, unloading and transportation. Zhu et al. [26] took the social cost of emergency rescue scheduling as the optimization objective, and presented an optimization model of post-disaster emergency materials dynamic scheduling considering heterogeneous behavior. Liu et al. [27] established a real-time dispatching and distribution model of emergency materials for epidemic relief based on the data-driven idea with the goal of minimizing emergency costs. Ge et al. [29] proposed an equity distribution model of emergency materials with the optimization objective of minimizing the loss of victims on the basis of considering the severity of disasters, demand points and attributes of emergency materials. Pang et al. [30] constructed an allocation optimization decision model based on a three-level emergency material transportation network with the goal of minimizing system loss.
There are also studies that combine different decision criteria/objectives to put forward corresponding emergency material allocation models [31][32][33][34][35][36]. For example, Tzeng et al. [31] developed a multi-stage emergency material distribution model with the goal of minimizing the total distribution time, minimizing the total cost and maximizing the satisfaction of victims. Liu et al. [33] took minimizing the overall loss comparison effect and distribution time comparison effect as the goal and constructed a double-layer cooperative optimization emergency material distribution model considering the irrational comparison psychology of victims. Zahedi et al. [35] proposed a multi-objective optimization decision model for emergency material allocation with the goal of minimizing unmet demand and minimizing the total cost, including vehicle costs and goods costs. Xue et al. [36] presented a multiobjective decision model for emergency material allocation based on capacity limitation, taking the minimum total network cost and the shortest average waiting time as the optimization objectives. However, it should be noted that the above studies are based on deterministic disaster information in the model construction.
With the continuous deepening of research, some studies on emergency material distribution began to gradually consider relevant uncertain factors [37][38][39][40]. For example, Gao [37] proposed a multi-commodity allocation rebalancing stochastic optimization model for earthquake relief considering uncertain transportation network and demand, aiming at minimizing the overall unmet level and expected transportation time. Zhang et al. [39] developed a two-stage emergency material distribution model combining disaster preparedness and disaster relief with the goal of minimizing the total cost on the basis of considering demand uncertainty.
In summary, a brief review of related research shows that it is crucial to consider multiple decision criteria and uncertain disaster information in emergency material distribution. However, there are still gaps in the existing research.

•
Firstly, previous research on emergency material distribution mostly focused on a single decision criterion; although some studies began to consider the combination of the different decision criteria, most of them focused on deterministic disaster information, ignoring the large amount of uncertain information in rescue and its impact on material distribution. • Secondly, in the relevant studies of emergency material distribution considering uncertain information, the consideration and trade-off of multiple decision criteria are usually ignored.
Therefore, this study attempts to develop an optimization model of emergency material distribution that balances multiple decision criteria under uncertain conditions. Compared with most existing works, four different decision criteria in the case of emergency material distribution-efficiency, economy, effectiveness and equity-are simultaneously considered and balanced in this paper. We also discussed the overall and periodization effects of different decision criteria on the distribution strategy of emergency materials to help obtain high-quality or optimal material distribution plan and provide support and reference for the decision of emergency material distribution in actual large-scale disaster relief responses.

Notation
The selection of sets, parameters and variables in this study is related to the reality of emergency material distribution.
First of all, regarding the selection of the sets, emergency material distribution generally refers to a logistics activity that distributes various types of emergency materials (such as medicines, daily necessities and rescue equipment) from supply centers to affected sites within the emergency response time period [41]. Therefore, this study chooses supply centers M, affected sites N, emergency materials I and emergency periods T as the sets for model building.
Secondly, the selection of parameters, especially the selection of uncertain parameters, can be input into the model in the form of fuzzy numbers to make them obey certain distribution rules. The reason is that the sudden and destructive nature of the disaster makes it impossible to obtain completely accurate information during the emergency rescue process. The demand and supply of materials, as well as transportation costs, are highly uncertain. This uncertain information can be input into the model as fuzzy numbers and this has been successfully done by previous researchers [42]. Therefore, this paper chooses to use the uncertain parameters with disturbance coefficient to represent the incompleteness of material demand information acquisition and the fuzzy uncertainty of transportation cost in emergency situations and uses the normal distribution to represent the uncertainty of material supply caused by inaccurate or delayed information in emergency contexts. Thirdly, for the selection of variables, the ultimate purpose of model construction and solution is to obtain the optimal emergency material distribution scheme, including whether to conduct distribution, the quantity of material distribution, the quantity of shortage, the quantity of inventory, the satisfaction of distribution, etc. Therefore, these variables are selected for model building in this paper.
The sets and indices, decision variables and parameters used in the model formulation are shown in Table 1. Binary variable indicating whether materials are distributed to affected site n ∈ N from supply center m ∈ M during time period t ∈ T. If so, the value is 1; otherwise, it is 0 S t ni Shortfalls of material i ∈ I at affected site n ∈ N at the end of each time period t ∈ T x t mni Amount of material i ∈ I allocated to affected site n ∈ N from supply center m ∈ M during time period t ∈ T J t mi Inventory of material i ∈ I at supply center m ∈ M at the end of each time period t ∈ T η t ni Satisfaction rate of emergency material i ∈ I at affected site n ∈ N during time period t ∈ T

Parameters D mn
Shortest distance from supply center m ∈ M to affected site n ∈ N V t mn Vehicle traveling speed from supply center m ∈ M to affected site n ∈ N during time period t ∈ T L t mni Loading and uploading time of unit material i ∈ I from supply center m ∈ M to affected site n ∈ N during time period t ∈ T ∂ t ni Urgency of affected site n ∈ N for emergency material i ∈ I during time period t ∈ T. The higher the value, the more urgent the affected site's demand for such materials β t n Vulnerability of affected site n ∈ N during time period t ∈ T, and it can be comprehensively determined according to disaster intensity, disaster degree, characteristics of victims and disaster carrying capability. The greater the vulnerability value is, the greater the possible loss caused by disasters at the affected site ρ Confidence level that makes constraints hold a i Space required to store per unit of material i ∈ I U t m Available inventory capacity at supply center m ∈ M during period t ∈ T Pr{·} Probability that constraint conditions in {·} hold C t i Unit purchase cost of material i ∈ I during time period t ∈ T H t mni Unit cost of loading and uploading material i ∈ I from supply center m ∈ M to affected site n ∈ N during time period t ∈ T G t mn Fixed cost of transportation from supply center m ∈ M to affected site n ∈ N during time period t ∈ T R t ni New forecast demand of emergency material i ∈ I at affected site n ∈ N during time period t ∈ T expressed by uncertain parameters with disturbance coefficient µ. R t ni = r t ni ,r t ni + µ ·r t ni , wherer t ni represents the minimum nominal amount of new forecast demand, µ ·r t ni represents the maximum disturbance value of the new forecast demand, and the disturbance coefficient µ can be determined according to the actual emergency situation Q t mi Fuzzy quantity of the new raised of emergency material i ∈ I at the supply center m ∈ M during time period t ∈ T. It obeys a normal distribution with mean E Q t mi and variance Unit mileage transportation cost of per unit of material i ∈ I from supply center m ∈ M to affected site wherek t mni represents the lowest nominal unit mileage transportation cost, µ ·r t ni represents the maximum disturbance value of the unit mileage transportation cost, and the disturbance co-efficient µ can be determined according to the actual emergency objective conditions Q t mi Actual supply of material i ∈ I at supply center m ∈ M during time period t ∈ T R t ni Actual demand for material i ∈ I at affected site n ∈ N during time period t ∈ T Y t mn Maximum total amount of material i ∈ I that can go from supply center m ∈ M to affected site n ∈ N during time period t ∈ T Processes 2022, 10, 2317 5 of 16

Mathematical Model
This paper constructs an optimization model of emergency material distribution that balances efficiency, economy, effectiveness and equity. The objective Function (1) minimizes the total time of material distribution in all emergency periods to pursue the efficiency criterion. The objective Function (2) minimizes the total costs for distributing materials in all emergency periods to pursue the economic criterion. The objective Function (3) maximizes the overall satisfaction rate at each affected site in all emergency periods to pursue the effectiveness criterion. The objective Function (4) minimizes the total losses caused by unmet demand at each affected site in all emergency periods to pursue the equity criterion. Constraint (5) represents the inventory capacity constraints of the supply center. Constraints (6) and (7) show demand and supply constraints, respectively. Constraint (8) indicates transportation capacity constraints. Constraint (9) is the constraint to satisfying as much demand as possible. Constraint (10) is a binary variable indicating whether the connection from the supply center to the affected site is used during time period t. Constraints (11) and (12) are, respectively, the expressions of actual demand and actual supply. Constraints (13) and (14) are, respectively, the expressions of shortfalls and satisfaction rate of materials. Constraint (15) is the value constraint of binary variables, and constraint (16) is a nonnegative constraint of decision variables. The objective functions and constraints of the proposed model are as follows.

•
Objective Function

Solution Method
The solution method of the proposed model is divided into two steps: the first step is to transform the uncertain parameters in the model with certainty, and the second step is to convert the proposed multi-objective model into a single-objective form for solution.
According to relevant research [43], it is assumed that there is a disturbance parameter µ that changes around the nominal value E, W is the fuzzy uncertain parameter value, then introducing the optimization constraint level σ (σ ∈ [0, 1], which can be set according to the disaster situation, disaster area and rescue reality), and let Taking the uncertain parameter K t mni with disturbance coefficient in the constructed model as an example, K t mni ≤ σ · k t mni +k t mni . The objective Function (2) can be transformed as follows: (17) Similarly, other uncertain parameters with disturbance coefficients can be transformed into deterministic forms.
According to the chance constrained programming problem proposed by Charnes and Cooper [44], if the function g(x, ψ) = h(x) − ψ, then for any given confidence level Assume that the random distribution vector is ψ = (k 1 , k 2 , · · · , k n , b), and the form of the random constraint function is g(x, ψ) = k 1 x 1 + k 2 x 2 + · · · + k n x n − b. If k ω and b follow a normal distribution and are independent of each other, then for any given confidence level ρ(0 ≤ ρ ≤ 1), if and only if: where E[k ω ] is the expectation, and Var[k ω ] is the variance. Based on the above analysis, the Constraint (7) of the proposed model can be transformed into the following deterministic form: • Step 2. Multi-objective transformation. We first normalize the objective functions of benefit type and cost type, where F min j (x t mni ) and F max j (x t mni ) are the minimum and maximum value of the objective function, respectively, and j is the number of objective functions.
Then we determine the weight coefficients of each objective function. In the case of emergency rescue, the importance of each decision objective in different rescue stages is different. The weight coefficient can be determined by relevant decision makers and emergency experts according to the degree of disaster, the vulnerability of victims, the urgency of demand, and should be satisfied Finally, the following single-objective model is constructed.

Case Description
This part takes the distribution of emergency medical materials during the COVID-19 epidemic as the background to design a research case to verify the effectiveness and feasibility of the proposed model. Wuhan (WH), Xiaogan (XG), Huanggang (HG) and Jingzhou (JZ), which were seriously affected by the epidemic, are selected as the affected sites for emergency medical materials, and Changsha City (CS) and Hefei City (HF) were selected as supply centers. Disposable surgical masks (KZ) and medicine (Lianhua Qingwen capsules (YP)) were selected as the urgently needed materials in the epidemic area, and the units are the thousand pieces and the thousand boxes, respectively. As part of the data could not be accurately obtained in the emergency rescue phase of the epidemic, the relevant parameters were set by combining the actual data with some simulation data. The demand for emergency materials in each period can be estimated according to the daily usage per capita of various emergency materials, as shown in Table 2. The latest amount of emergency materials that can be raised in each period under the urgency of emergency rescue is estimated according to the epidemic situation, material production and mobilization capacity of emergency relief centers, as shown in Table 3. The average speed of highway transportation is set to be 100 km/h, the shortest distance from the supply center to the affected site is obtained by Baidu Maps, and the corresponding fixed transportation cost in the case of disasters is estimated, as shown in Table 4. The procurement cost, transportation cost, loading and unloading cost, loading and unloading time and required storage space per unit of emergency materials are shown in Table 5. The available inventory capacity at the supply center at the end of each time period is shown in Table 6. The maximum transportable volume of emergency materials from the supply center to the affected site during each time period is shown in Table 7, and the vulnerability of each affected site and the urgency of the demand for emergency material are shown in Table 8. This computational case was solved by using Lingo 11.0 on a computer with an Intel(R) Core(TM) 1.90 GHz processor with 16.0 GB of RAM, and the confidence level and constraint optimization level are, respectively, ρ = 0.95 and σ = 0.9.  a, b), where a is the nominal value and b is the disturbance coefficient. Table 3. Supply of emergency materials at the supply centers at the beginning of each time period.     The influence of different decision criteria on the overall strategy of emergency material distribution was explored by calculating the value of objective functions under the five decision criteria, as shown in Figure 1. The five decision criteria are as follows: the efficiency criterion W 1 = (w 1 , w 2 , w 3 , w 4 ) = (1, 0, 0, 0), that is, the decision weights of the four objective functions under the efficiency criterion are (1, 0, 0, 0) respectively, denoted as W 1 = (1, 0, 0, 0). By analogy, the economic criterion is denoted as W 2 = (0, 1, 0, 0), the effectiveness criterion is denoted as W 3 = (0, 0, 1, 0), the equity criterion is denoted as W 4 = (0, 0, 0, 1), the balance-criterion (considering four objectives simultaneously) is denoted as W * = (0.25, 0.25, 0.25, 0.25).

Influence of Different Decision Criteria on Multi-Period Emergency Material Distribution
The time, distribution cost, satisfaction rate and loss of emergency material distribu tion, as well as the evolution trend in each time period under the condition of adhering to the efficiency criterion (W1), economic criterion (W2), effectiveness criterion (W3), equity criterion (W4) and balance-criterion (W*), are shown in Figures 2-5, respectively.
It can be seen from Figures 2-5 that the five decision criteria all reflect the following trends: with the passage of time period, the time, cost and satisfaction rate of emergency material distribution in each time period gradually increase, and the system loss of mate rial distribution gradually decreases until it reaches zero. The main reasons are as follows At the beginning of emergency rescue, the supply centers are able to provide limited emer gency materials, and the time and cost required to distribute the limited materials to af fected sites are relatively low. However, at this time, all the affected sites tend to need large number of rescue materials. Due to only a small amount of materials being obtained the satisfaction rate of material distribution in the early period of rescue is low, which As can be seen from Figure 1, the total time, total allocation cost, total satisfaction rate and total loss of emergency material distribution are significantly different under different decision criteria, indicating that different decision criteria have an important impact on the overall strategy selection of emergency material distribution. We can find that balance-criterion has advantages in the process of large-scale multi-period emergency material distribution. As shown in Figure 1, compared with other decision criteria, under balance-criterion (W*), the distributing time, cost, satisfaction rate and loss of emergency materials allocation are above the medium level, indicating that the decision scheme formed based on the balance-criteria can properly balance and take into account various indicators in combination with the specific disaster situation and rescue information of different emergency periods in the material distribution process, thus avoiding more extreme distribution situations as much as possible.

Influence of Different Decision Criteria on Multi-Period Emergency Material Distribution
The time, distribution cost, satisfaction rate and loss of emergency material distribution, as well as the evolution trend in each time period under the condition of adhering to the efficiency criterion (W 1 ), economic criterion (W 2 ), effectiveness criterion (W 3 ), equity criterion (W 4 ) and balance-criterion (W*), are shown in Figures 2-5, respectively.           Further analysis of Figures 2-5 shows that the time of material distribution in each period is the shortest under adherence to the efficiency criterion, but the system loss is the largest and the coverage satisfaction rate is almost the lowest. The cost of distribution per period is the lowest based on the economic criterion, but the satisfaction rate is very low Under the effectiveness criterion, the satisfaction rate of material distribution in each pe riod is the highest, but the distribution time and cost are higher. Under the criterion o equity, the loss of material distribution in each period is the least, but it has the longes time and the highest cost. However, compared with other decision criteria, under adher ence to the balance-criterion, the time, cost, satisfaction rate and loss of emergency mate rial distribution in each time period have obvious advantages and meet the requirements of multi-period sustainable material distribution.

Multi-period Distribution Scheme of Emergency Materials under the Balance-Criterion of "Efficiency-Economy-Effectiveness-Equity"
Through calculation, the multi-period distribution scheme of emergency materials under balance-criterion is obtained, as shown in Table 9. The multi-objective evolution trend of multi-period distribution of emergency materials is shown in Figure 6. The time and cost required by each affected site to obtain emergency materials in each time period the satisfaction rate obtained and the system loss generated are shown in Figures 7-10 respectively. Table 9. Multi-period distribution scheme of emergency materials based on balance-criterion. It can be seen from Figures 2-5 that the five decision criteria all reflect the following trends: with the passage of time period, the time, cost and satisfaction rate of emergency material distribution in each time period gradually increase, and the system loss of material distribution gradually decreases until it reaches zero. The main reasons are as follows: At the beginning of emergency rescue, the supply centers are able to provide limited emergency materials, and the time and cost required to distribute the limited materials to affected sites are relatively low. However, at this time, all the affected sites tend to need a large number of rescue materials. Due to only a small amount of materials being obtained, the satisfaction rate of material distribution in the early period of rescue is low, which leads to large system losses. In the middle period of emergency rescue, with the continuous supply of all types of emergency materials, the quantity of materials needed to be distributed in each time period gradually increases, the time, cost and the satisfaction rate gradually increase, and the system loss gradually decreases. In the later period of emergency rescue, when the supply of materials is greater than or equal to the demand, the demand at all affected sites can be met, and the satisfaction rate reaches the maximum value, while the system loss decreases to the minimum; that is, when the demand is fully met, there is no loss caused by material shortage.
Further analysis of Figures 2-5 shows that the time of material distribution in each period is the shortest under adherence to the efficiency criterion, but the system loss is the largest and the coverage satisfaction rate is almost the lowest. The cost of distribution per period is the lowest based on the economic criterion, but the satisfaction rate is very low. Under the effectiveness criterion, the satisfaction rate of material distribution in each period is the highest, but the distribution time and cost are higher. Under the criterion of equity, the loss of material distribution in each period is the least, but it has the longest time and the highest cost. However, compared with other decision criteria, under adherence to the balance-criterion, the time, cost, satisfaction rate and loss of emergency material distribution in each time period have obvious advantages and meet the requirements of multi-period sustainable material distribution.

Multi-period Distribution Scheme of Emergency Materials under the Balance-Criterion of "Efficiency-Economy-Effectiveness-Equity"
Through calculation, the multi-period distribution scheme of emergency materials under balance-criterion is obtained, as shown in Table 9. The multi-objective evolution trend of multi-period distribution of emergency materials is shown in Figure 6. The time and cost required by each affected site to obtain emergency materials in each time period, the satisfaction rate obtained and the system loss generated are shown in Figures 7-10, respectively.  Figure 6. The multi-objective evolution trend of multi-period distribution of emergency materials based on the balance-criterion.
As can be seen from Figure 6, in general, the multi-period material distribution scheme formed based on the balance-criterion shows a trend of gradual increase in the total time, cost and satisfaction rate, and a trend of gradual decrease in the loss, which is in line with the reality of multi-period emergency materials distribution.  As can be seen from Figure 6, in general, the multi-period material distribution scheme formed based on the balance-criterion shows a trend of gradual increase in th total time, cost and satisfaction rate, and a trend of gradual decrease in the loss, which i in line with the reality of multi-period emergency materials distribution.       Figures 7-10 show that the distribution at each affected site is generally in line wit the above trend. However, it should be noted that when the affected site for a type o emergency material decreases, the required distribution time is shortened accordingly and the growth range of the required distribution cost is also slowed down. Meanwhil   Figures 7-10 show that the distribution at each affected site is generally in line with the above trend. However, it should be noted that when the affected site for a type o emergency material decreases, the required distribution time is shortened accordingly and the growth range of the required distribution cost is also slowed down. Meanwhile As can be seen from Figure 6, in general, the multi-period material distribution scheme formed based on the balance-criterion shows a trend of gradual increase in the total time, cost and satisfaction rate, and a trend of gradual decrease in the loss, which is in line with the reality of multi-period emergency materials distribution. Figures 7-10 show that the distribution at each affected site is generally in line with the above trend. However, it should be noted that when the affected site for a type of emergency material decreases, the required distribution time is shortened accordingly, and the growth range of the required distribution cost is also slowed down. Meanwhile, the satisfaction rate gradually rises, and the system loss gradually decreases (e.g., WH (KZ) after the third period in Figures 7-10).
In addition, as can be seen from Figures 7-10, compared with the material distribution of other affected sites, the distribution of material KZ for affected site WH requires a longer time and higher cost, but its satisfaction is lower and its loss is greater in the early period of emergency. The main reason is that in the early period of emergency rescue, compared with the demand at other affected sites, the affected site WH has the largest demand for material KZ, so the supply centers distribute more materials to this affected site on the basis of weighing all decision indicators, resulting in longer distribution time and higher distribution cost. However, the amount of materials obtained could not completely cover all the demand of the affected site WH, and its vulnerability and demand urgency coefficient were the highest among all affected sites at the initial period of emergency relief, which resulted in relatively high system loss. It should be noted that this is already the best global distribution strategy under the condition of balancing all affected sites. Therefore, it is particularly necessary to scientifically weigh the decision criteria based on the multi-period global optimization perspective and consider the vulnerability and demand urgency of different affected sites for emergency material distribution under the condition that supply falls short of demand.

Conclusions and Directions for Future Research
This paper proposed an optimization model for emergency material distribution that balances multiple decision criteria-efficiency, economy, effectiveness and equity, which are simultaneously considered and balanced. A case study was conducted to validate the proposed model. Accordingly, the following conclusions were reached.
(1) In the initial period of emergency relief, when multiple affected sites are faced with emergency peak demand, paying attention to the effectiveness criterion and equity criterion can avoid the low coverage satisfaction rate and high system loss caused by material shortage as much as possible. (2) When the initial rescue task is urgent and the time requirement is extremely high, the efficiency criterion can be concerned. In the case of sufficient supply of materials in the middle and later emergency periods, paying attention to economic criteria can improve the satisfaction rate and reduce system losses while completing the task of material distribution at a lower cost. However, from the perspective of long-term multi-period sustainable emergency rescue, the balance-criteria should be highly concerned and valued, and the weight selection of various decision criteria should be scientifically grasped. (3) A decision scheme based on balance-criteria can trade-off the efficiency, economic, effectiveness and equity of emergency material distribution, achieving in a shorter time and lower cost for higher satisfaction rate and smaller system loss and thus avoiding the extreme or non-global optimization of distribution schemes caused by only paying attention to a certain decision criterion.
In the future, how to obtain real-time disaster information to distribute emergency materials more scientifically is a direction worthy of thinking and studying. Furthermore, as multi-period distribution of emergency supplies is a complex and systematic process, with the increase of the scale and complexity of the research problem, it is necessary to design more complex models and more effective model-solving methods with the increase of the scale and complexity of the research problem. Therefore, how to use big data technology to obtain real-time information, such as supply, demand and transportation capacity, and to design scientific and efficient solutions in the multi-time distribution process of large-scale disaster emergency relief operations still needs further discussion and consideration. In short, it is hoped that the proposed model can have certain reference and promotion values for large-scale disaster emergency material allocation under uncertain conditions. Funding: This research was funded by the National Social Science Foundation of China, grant number 21CGL044.