# Equilibrium Decision Method for Earthquake First-Aid Medicine Allocation Based on Demand Information Updating

Authors: Yong Ye; Shalei Zhan; Shih-Yung Wei
Journal: Mathematical Problems in Engineering
 (2017)

---

## Abstract

The allocation of rescue resources after an earthquake has become a popular research topic in the field of emergency management. The allocation of first-aid medicine for earthquake rescue has stronger time sensitivity than that of general rescue materials. This study focuses on the problem of first-aid medicine allocation in earthquake response. First, we consider the incompleteness and renewal of decision information in an emergency environment, as well as the balance between the risk of decision error and delay. Second, we propose an equilibrium decision method for the allocation of first-aid medicine in earthquake rescue based on information update. This method attempts to realize a fair allocation to all disaster places and minimize total transport time loss. Third, a simulation analysis is performed in which the proposed method is applied to the first-aid medicine allocation problem in the Wenchuan earthquake response. Results show that the method can be used to create a good allocation plan in an earthquake rescue situation.


---

## Body

## 1. Introduction

Earth has been subjected to earthquakes during its entire geological period. Historically, the frequent occurrence of catastrophic or major earthquakes has resulted in a considerable loss of life and property worldwide. The magnitude 8.3 earthquake that hit San Francisco on April 18, 1906, took more than 60,000 lives. The great Kanto earthquake that occurred in Japan on September 1, 1923, destroyed 600,000 buildings, took 143,000 lives, and left 200,000 injured and 500,000 homeless. The terrible fire that occurred in the business district of Tokyo at the same time took 240,000 lives in Yokohama Park. Earthquakes also trigger landslides, debris flow, and other secondary disasters, all of which cause large numbers of casualties. A magnitude 8.9 earthquake with an epicentral intensity of more than 11 hit Chile on May 22, 1980, and hundreds of stronger earthquakes followed continuously. Among the earthquakes that followed, 3 were more than 8.0 in magnitude and 10 were more than 7.0. The great Kanto earthquake also triggered a violent tsunami and a severe volcanic eruption, taking 5700 lives. A magnitude 8.1 earthquake occurred in the Pacific Ocean, off the coast of the Mexican state, on September 19, 1985. A magnitude 7.5 aftershock occurred the next day, taking more than 7000 lives, causing injury to 11,000 people, and leaving more than 0.3 million homeless; the economic loss totaled USD 1.1 billion. The magnitude 7.2 Great Hanshin earthquake occurred on January 17, 1995, taking 5466 lives and leaving more than 30,000 injured and hundreds of thousands of others homeless; the victims reached more than 1.4 million. A magnitude 7.4 earthquake occurred in the western part of Turkey on August 17, 1999, taking more than 13,000 lives. A magnitude 8.7 earthquake occurred on December 26, 2004, and a tsunami was triggered off the northwest coast of Sumatra, Indonesia, killing 188,149 people in the countries around the Indian Ocean and with 43,861 still missing. A magnitude 7.8 earthquake occurred in Pakistan on October 8, 2005, affecting Pakistan, India, Afghanistan, and Indonesia. More than 79,000 Pakistanis lost their lives, over 60,000 were injured, and 3.3 million were left homeless. A magnitude 7.4 earthquake occurred in Haiti on January 12, 2010, taking 222,500 lives and injuring 196,000 people. A magnitude 9.0 earthquake occurred in Japan on March 11, 2011, and triggered a tsunami, causing unit numbers 1–4 of the Fukushima I Nuclear Power Plant to develop a nuclear leak, taking 13,498 lives, and leaving 14,734 people still missing. Clearly, global earthquake relief work is arduous.Unlike other disasters, earthquakes are characterized by strong suddenness, huge destructiveness, and high defense difficulty. After an earthquake, the demands for drugs, food, tents, and other supplies in disaster areas increase sharply within a short period. The construction of emergency resources in China has many weak links, and the demands of disaster emergency cannot be fully met. For example, a shortage of emergency rescue materials occurred during the Wenchuan (WC), Yushu, and Ludian earthquakes. The shortage or delay of emergency medical supplies, particularly for sterilization, anti-infection, anesthesia, and stanch bleeding, prevents the timely treatment of victims. For example, numerous victims died in Indonesia earthquake because of drug shortage.In medical rescue, treatment during the first three days after an earthquake is critical; the effect of medical rescue will gradually weaken after this period [1]. Drug allocation for earthquake first aid has higher time sensitivity than that for general rescue supplies. On the one hand, the acquired decision-making information is insufficient if a decision is made instantly, and thus decision-making efficiency may be low. This phenomenon can reduce the urgency of decision-making time in the late stage, resulting in invalid execution and the complexity of problems. Thus, the risk of decision error is high and considerable manpower and material resources are consumed. On the other hand, the acquired decision-making information is more robust if a decision is made later, and thus decision-making efficiency may be better. This situation can reduce the risk of decision error but may delay decision making because of the further deterioration of the disaster, thereby increasing the risk of delayed implementation. An ideal allocation plan can be obtained once a simple program is designed if complete information on the demand for rescue drugs, resource supply, and transportation is obtained. However, the occurrence and development of earthquake disaster are accompanied by numerous uncertainties, and acquiring complete disaster information is difficult. Therefore, a decision maker needs to utilize the available information, which includes sample information obtained at the current stage and historical information of related disasters, to make a scientific and reasonable allocation plan. The sample information at the current stage can be continuously updated through observation, that is, information renewal of drug allocation for earthquake rescue. Multiple vectors are used to describe the related information in the decision-making process of emergency resource allocation for earthquake disaster, namely, group information. The group information for the allocation decision of rescue drugs is continuously acquired and updated with the constant disaster evolution, and the related allocation decision is dynamically optimized.Therefore, a decision maker should consider several problems in the decision-making process, including when to allocate, where to allocate, how much to allocate, and how to allocate. He or she needs to focus on the balance between the decision time and risk of rescue medicine allocation under an uncertain environment and the dynamic renewal of demand information. Based on the summary of results of previous studies and using Bayesian analysis theory, this study comprehensively considers the balance of risks of decision error and delay and constructs a decision-making method for rescue medicine allocation under the renewal of demand information. This method is applied to an example of rescue medicine allocation for the WC earthquake to perform a simulation analysis. The result shows that an emergency decision maker can update data by inputting the corresponding demand information, acquiring the optimal decision time under the balance between decision time and risk of rescue medicine allocation, and creating an allocation plan with the minimum transport time under a fair allocation principle.


## 2. Research Status

The following paragraphs summarize the related studies on the supply and demand for earthquake first-aid medicine and decision-making method for the allocation of disaster rescue supplies.


### 2.1. Study on the Supply and Demand for Earthquake First-Aid Medicine

Studies on emergency medicine have achieved fruitful results, mainly by focusing on the development of emergency medicine, emergency medical service, and analysis of all types of common diseases after a disaster [2]. At present, studies on the supply and demand for earthquake first-aid medicine utilize the practical data of earthquake medical rescue to analyze the main injury types caused by earthquake disasters. Moreover, these studies use the supply, demand, and application of related drugs to raise suggestions on enhancing the reservation, management, and application of earthquake first-aid drugs.Studies on the supply of earthquake first-aid drugs mainly include the supply characteristics, supply security, donation management of rescue drugs, production purchasing, dispatching, transportation, and existing problems. Donation is the major source of earthquake first-aid drugs. However, some donation drugs are expired, of poor quality, or useless; thus, the use efficiency of drugs is low. A total of 5000 tons of drugs and medical tools were received in the rescue process of the San Francisco earthquake in 1989, but only approximately 30% could be used immediately [3]. Xiao et al. [4] indicated that the dispatch, distribution, donation rules, and use management of rescue drugs in China need to be standardized further. Mori et al. [5] also analyzed the significant effect of drug shortage on the rescue during the large earthquake in Japan in 2011.Studies on the demand for earthquake first-aid drugs according to historical earthquake first-aid practices have analyzed the demand for earthquake first-aid drugs on the basis of the studies on earthquake injury and medication. Fracture, soft tissue injury, and crush syndrome are the most significant factors in earthquake-disaster-caused casualties [6]. Kazancioglu et al. [7] analyzed the change in the demand for blood of crush syndrome patients by using first-aid data from the Marmara earthquake. Their result showed that the demand for blood was highest during the first week after the earthquake. Kudo et al. [8] investigated the reservation and demand for drugs, medical tools, and other medical supplies of 14 main hospitals after the East Japan earthquake. Their result showed that the demands of all hospitals for medicine were accommodated on the fourth day after the earthquake.


### 2.2. Study on the Decision-Making Method for the Allocation of Disaster Rescue Supplies

The allocation of disaster rescue supplies is the main content of an emergency logistics system operation. Existing related studies are mainly focused on solving emergency logistics, transshipment center location, and distribution route selection in an uncertain environment through the model construction and design of solving algorithms.Studies on the decision-making model of the allocation of disaster rescue supplies have constructed an allocation and dispatch model to satisfy the requirements of a disaster area maximally. Meanwhile, efficiency and fairness have been considered under the conditions of sharply increasing demand, uncertain supply, and road damage by stochastic programming, robust optimization, Bayesian analysis, space-time network, mixed integer programming, and other theories. Sheu [9] constructed an allocation model of emergency supplies with the maximization of survivor fitness as the objective by starting from the demand and using stochastic programming theory. Haghani et al. [10] established a time-space network model for selecting the distribution route of disaster rescue supplies. Özdamar and Yi [11] established a mixed integer programming model based on emergency logistics combined with crowd evacuation, resource allocation, and route selection.Studies on the solving algorithm of the allocation model of disaster rescue supplies have considered that the allocation of disaster rescue supplies is a nondeterministic polynomial time-hard problem that involves various factors, such as multiple objectives, starting points, demand points, transportation modes, materials, and agents. The solving algorithms of related decision-making models include Lagrangian relaxation technique, genetic algorithm, bionic algorithm, ant colony optimization, and particle swarm optimization. Chang et al. [12] introduced a multiobjective genetic algorithm based on greedy search to solve a large-scale emergency resource allocation model. Zhang et al. [13] designed a bionic algorithm against emergency logistic network route programming problems. Yuan and Wang [14] solved a multiobjective emergency logistics route selection model by using ant colony optimization.Studies on the supply and demand for earthquake first-aid medicines and the allocation of disaster rescue supplies are relatively abundant. By contrast, studies on the supply and demand for an adaptation strategy of the allocation of rescue medicines are limited. In particular, no study has evaluated the equilibrium decision method between decision time for rescue medicine allocation and risk of decision error from the perspective of information updating in an uncertain environment. In previous studies, a particular algorithm was required to design the solutions to a model because of the model complexity. The solving efficiency was not high, and only a satisfactory solution could be obtained. In consideration of the urgency in an emergency environment, the related model can not only increase the solving efficiency, but also obtain the optimum solution if it is solved using a mature mathematical software. On the basis of the preceding analysis, this study integrates an information updating method with the balanced allocation of earthquake first-aid medicine using Bayesian analysis theory to achieve the balance between the risk of decision error and risk of decision delay while staying true to the principle of fairness. The study likewise aims to use a decision method for earthquake first-aid medicine allocation, which is solved using mature mathematical software.


## 2.1. Study on the Supply and Demand for Earthquake First-Aid Medicine

Studies on emergency medicine have achieved fruitful results, mainly by focusing on the development of emergency medicine, emergency medical service, and analysis of all types of common diseases after a disaster [2]. At present, studies on the supply and demand for earthquake first-aid medicine utilize the practical data of earthquake medical rescue to analyze the main injury types caused by earthquake disasters. Moreover, these studies use the supply, demand, and application of related drugs to raise suggestions on enhancing the reservation, management, and application of earthquake first-aid drugs.Studies on the supply of earthquake first-aid drugs mainly include the supply characteristics, supply security, donation management of rescue drugs, production purchasing, dispatching, transportation, and existing problems. Donation is the major source of earthquake first-aid drugs. However, some donation drugs are expired, of poor quality, or useless; thus, the use efficiency of drugs is low. A total of 5000 tons of drugs and medical tools were received in the rescue process of the San Francisco earthquake in 1989, but only approximately 30% could be used immediately [3]. Xiao et al. [4] indicated that the dispatch, distribution, donation rules, and use management of rescue drugs in China need to be standardized further. Mori et al. [5] also analyzed the significant effect of drug shortage on the rescue during the large earthquake in Japan in 2011.Studies on the demand for earthquake first-aid drugs according to historical earthquake first-aid practices have analyzed the demand for earthquake first-aid drugs on the basis of the studies on earthquake injury and medication. Fracture, soft tissue injury, and crush syndrome are the most significant factors in earthquake-disaster-caused casualties [6]. Kazancioglu et al. [7] analyzed the change in the demand for blood of crush syndrome patients by using first-aid data from the Marmara earthquake. Their result showed that the demand for blood was highest during the first week after the earthquake. Kudo et al. [8] investigated the reservation and demand for drugs, medical tools, and other medical supplies of 14 main hospitals after the East Japan earthquake. Their result showed that the demands of all hospitals for medicine were accommodated on the fourth day after the earthquake.


## 2.2. Study on the Decision-Making Method for the Allocation of Disaster Rescue Supplies

The allocation of disaster rescue supplies is the main content of an emergency logistics system operation. Existing related studies are mainly focused on solving emergency logistics, transshipment center location, and distribution route selection in an uncertain environment through the model construction and design of solving algorithms.Studies on the decision-making model of the allocation of disaster rescue supplies have constructed an allocation and dispatch model to satisfy the requirements of a disaster area maximally. Meanwhile, efficiency and fairness have been considered under the conditions of sharply increasing demand, uncertain supply, and road damage by stochastic programming, robust optimization, Bayesian analysis, space-time network, mixed integer programming, and other theories. Sheu [9] constructed an allocation model of emergency supplies with the maximization of survivor fitness as the objective by starting from the demand and using stochastic programming theory. Haghani et al. [10] established a time-space network model for selecting the distribution route of disaster rescue supplies. Özdamar and Yi [11] established a mixed integer programming model based on emergency logistics combined with crowd evacuation, resource allocation, and route selection.Studies on the solving algorithm of the allocation model of disaster rescue supplies have considered that the allocation of disaster rescue supplies is a nondeterministic polynomial time-hard problem that involves various factors, such as multiple objectives, starting points, demand points, transportation modes, materials, and agents. The solving algorithms of related decision-making models include Lagrangian relaxation technique, genetic algorithm, bionic algorithm, ant colony optimization, and particle swarm optimization. Chang et al. [12] introduced a multiobjective genetic algorithm based on greedy search to solve a large-scale emergency resource allocation model. Zhang et al. [13] designed a bionic algorithm against emergency logistic network route programming problems. Yuan and Wang [14] solved a multiobjective emergency logistics route selection model by using ant colony optimization.Studies on the supply and demand for earthquake first-aid medicines and the allocation of disaster rescue supplies are relatively abundant. By contrast, studies on the supply and demand for an adaptation strategy of the allocation of rescue medicines are limited. In particular, no study has evaluated the equilibrium decision method between decision time for rescue medicine allocation and risk of decision error from the perspective of information updating in an uncertain environment. In previous studies, a particular algorithm was required to design the solutions to a model because of the model complexity. The solving efficiency was not high, and only a satisfactory solution could be obtained. In consideration of the urgency in an emergency environment, the related model can not only increase the solving efficiency, but also obtain the optimum solution if it is solved using a mature mathematical software. On the basis of the preceding analysis, this study integrates an information updating method with the balanced allocation of earthquake first-aid medicine using Bayesian analysis theory to achieve the balance between the risk of decision error and risk of decision delay while staying true to the principle of fairness. The study likewise aims to use a decision method for earthquake first-aid medicine allocation, which is solved using mature mathematical software.


## 3. Equilibrium Decision-Making Method without Information Updating

### 3.1. Problems and Symbols

Many disaster points have required emergency rescue medicine supplies after the 2011 East Japan earthquake and other major earthquake disasters. However, local emergency reserve drugs are difficult to maintain during such a sudden increase in demands. Therefore, we need to supply related drugs from other supply points outside the disaster areas. Supposing thatm supply points and n disaster points exist, the other parameters are shown in Notations: (i) Symbols.


### 3.2. Hypotheses and Definitions

Hypothesis 1.
Compared with life-type emergency rescue supplies, rescue medicine is characterized by a relatively small demand despite its importance. The distribution of rescue medicine can be prioritized if necessary. Therefore, this study assumes that the supply point for rescue medicine has an adequate distribution capability. Hence, transport capacity is not considered in the model.Hypothesis 2.
The requirement of pointj in t can be expressed by the population transfer rate. The requirement is denoted by the stochastic variable θjt and submits to the detailed distribution. The probability density function is π(θjt).Hypothesis 3.
The supply in the supply point is confirmed.Definition 1.
The requirement of pointj in t is the function of population transfer rate θjt, which is defined as follows:(1)djt=dθjt.Definition 2.
The loss of decision errorRtError in t is a function related to the prediction accuracy pt, which is defined as follows:(2)RtError=1-pt,where p(t) is a monotone increasing function in decision time t.Definition 3.
The loss of decision delayRtDelay in t is a monotone increasing function related to decision time, which is defined as follows:(3)RtDelay=Rt.Definition 4.
The loss of transport timeLtT in t is the total time of medicine distribution from each supply point to each disaster point, which is defined as follows:(4)LtT=∑i=1m∑j=1ntij0·sijt.Definition 5.
The fair allocation principle ensures that the expected demand satisfaction rates of all disaster points are equal. Thus, the expected demand satisfaction rate of pointj in t is defined as follows:(5)ujt=min⁡∑j=1nEdjt,∑i=1mSit∑j=1nEdjt.
In (5), when ∑j=1nE(djt)>∑i=1mSit, the expected demand satisfaction rate ujt is less than 1 because the supply is not enough. However, when ∑j=1nE(djt)≤∑i=1mSit, the expected demand satisfaction rate ujt equals 1, which means that the supply is enough to satisfy all the demands.


### 3.3. Decision-Making Model

According to the above hypotheses and definitions, an equilibrium configuration model of earthquake first-aid medicine can be established in an uncertain environment. The specific objectives are as follows: (I) balanced risks of decision error and delay, (II) similar expected demand satisfaction rates of all disaster points, and (III) smallest total loss of transportation time. The model can be established as follows:(6)min⁡z=∑i=1m∑j=1ntij0·sijt.Thus,(7)uj′t=uj′′t,∀j′,j′′∈J,(8)∑j=1nsijt=Sit,∀i∈I,(9)RtError-RtDelay2≤α,(10)0≤t≤T0,sijt≥0.Equation (6) refers to the objective function; constraint condition (7) indicates that the expected demand satisfaction rates of all disaster points are the same; constraint condition (8) meets the total supply limit of all supply points; constraint condition (9) refers to the equilibrium condition between the risks of decision error and delay; that is, the square of the difference between both is not more than α; and constraint condition (10) refers to the range of the variable. The model can be solved using Lingo software.


## 3.1. Problems and Symbols

Many disaster points have required emergency rescue medicine supplies after the 2011 East Japan earthquake and other major earthquake disasters. However, local emergency reserve drugs are difficult to maintain during such a sudden increase in demands. Therefore, we need to supply related drugs from other supply points outside the disaster areas. Supposing thatm supply points and n disaster points exist, the other parameters are shown in Notations: (i) Symbols.


## 3.2. Hypotheses and Definitions

Hypothesis 1.
Compared with life-type emergency rescue supplies, rescue medicine is characterized by a relatively small demand despite its importance. The distribution of rescue medicine can be prioritized if necessary. Therefore, this study assumes that the supply point for rescue medicine has an adequate distribution capability. Hence, transport capacity is not considered in the model.Hypothesis 2.
The requirement of pointj in t can be expressed by the population transfer rate. The requirement is denoted by the stochastic variable θjt and submits to the detailed distribution. The probability density function is π(θjt).Hypothesis 3.
The supply in the supply point is confirmed.Definition 1.
The requirement of pointj in t is the function of population transfer rate θjt, which is defined as follows:(1)djt=dθjt.Definition 2.
The loss of decision errorRtError in t is a function related to the prediction accuracy pt, which is defined as follows:(2)RtError=1-pt,where p(t) is a monotone increasing function in decision time t.Definition 3.
The loss of decision delayRtDelay in t is a monotone increasing function related to decision time, which is defined as follows:(3)RtDelay=Rt.Definition 4.
The loss of transport timeLtT in t is the total time of medicine distribution from each supply point to each disaster point, which is defined as follows:(4)LtT=∑i=1m∑j=1ntij0·sijt.Definition 5.
The fair allocation principle ensures that the expected demand satisfaction rates of all disaster points are equal. Thus, the expected demand satisfaction rate of pointj in t is defined as follows:(5)ujt=min⁡∑j=1nEdjt,∑i=1mSit∑j=1nEdjt.
In (5), when ∑j=1nE(djt)>∑i=1mSit, the expected demand satisfaction rate ujt is less than 1 because the supply is not enough. However, when ∑j=1nE(djt)≤∑i=1mSit, the expected demand satisfaction rate ujt equals 1, which means that the supply is enough to satisfy all the demands.


## 3.3. Decision-Making Model

According to the above hypotheses and definitions, an equilibrium configuration model of earthquake first-aid medicine can be established in an uncertain environment. The specific objectives are as follows: (I) balanced risks of decision error and delay, (II) similar expected demand satisfaction rates of all disaster points, and (III) smallest total loss of transportation time. The model can be established as follows:(6)min⁡z=∑i=1m∑j=1ntij0·sijt.Thus,(7)uj′t=uj′′t,∀j′,j′′∈J,(8)∑j=1nsijt=Sit,∀i∈I,(9)RtError-RtDelay2≤α,(10)0≤t≤T0,sijt≥0.Equation (6) refers to the objective function; constraint condition (7) indicates that the expected demand satisfaction rates of all disaster points are the same; constraint condition (8) meets the total supply limit of all supply points; constraint condition (9) refers to the equilibrium condition between the risks of decision error and delay; that is, the square of the difference between both is not more than α; and constraint condition (10) refers to the range of the variable. The model can be solved using Lingo software.


## 4. Equilibrium Decision-Making Method under Demand Information Updating

### 4.1. Symbols and Hypotheses

The decision of earthquake first-aid medicine allocation cannot be in place in just one step but is continuously updated with the development of an earthquake disaster. Such a decision is a dynamic allocation based on the related information updates. The balanced allocation of earthquake first-aid medicine based on the demand information update can be modeled using Bayesian analysis [15], which is used to establish related models. The symbols are presented in Notations: (ii) New Symbols.In an earthquake first-aid medicine allocation model under demand information updating, a decision maker can adopt various rescue medicine allocation plan setsA={a} (the plan mainly refers to the quantity of medicine distributed from all supply points to all disaster points) to describe various possible status sets of population transfer rate Θ={θ}. The decision-making rule δ(x) shows the action taken [δ(x) refers to the function from X to A] and X=x is the observation value of the sample information. The loss and risk functions are based on the following hypotheses.Hypothesis 4.
The density function of population transfer rate in prior distributionπ(θ) is known, and the sample information related to θ can be continuously observed and obtained in cycle T. The observation can be conducted numerous times in cycle T. The comprehensive observation value is obtained.The Bayes decision method comprehensively uses the prior information of parameters and sample information to make a decision. The density function of disaster information in prior distributionπ(θ) can be derived, where θ is a multidimensional random variable, based on the disaster databases from government departments (e.g., the National Climatic Data Center of the National Oceanic and Atmospheric Administration, the National Geophysical Data Center, and the National Oceanographic Data Center).The population transfer rate information observation is expressed byX1,X2,…. Xk=(X1k,X2k,…) is the sample observed in stage k, and xk refers to the sample information. The condition density function is fk(xk∣θk), which is defined in Xk=X1×X2×⋯×Xk. The posterior probability density while θk is in stage k is πk∗(θk∣xk). δ(xk)=δ(x1k,x2k,…) is the allocation plan for all disaster points.Hypothesis 5.
A difference exists between the emergency supply allocation and actual demand, causing a decrease in the utility of the limited material, which is defined as the first-aid medicine mismatch error losses (MMELs).A difference always exists between first-aid medicine allocation and actual demand because of the incompleteness of information. Excess or deficiency will reduce the effectiveness of drug allocation. The effectiveness of drug allocation will be significant if the allocation error is high. If the demand for resources is set tod(θ) and the supply is s, the first-aid MMEL function can be expressed as follows:(11)LFDMELθ,δ=dθ-sβ,β∈R.


### 4.2. Fair Allocation Principle Based on Demand Information Updating

Equation (11) indicates that the first-aid MMEL function in demand point j in t is(12)LjtFMMELθjt,δjt=dθjt-sijtβ,β∈R.Int, the first-aid MMEL function of all demand points is(13)LFMMELθ,δ=∑j=1ndθjt-∑i=1msitβ,β∈R.Equations (12) and (13) do not consider the observation sample information Xk. Consequently, the effectiveness of decision is affected and the updated data cannot be used. According to the formula of Bayes risk defined in literature (Berger, 1980), the Bayes risk function of first-aid MMEL in demand point j in t is defined as follows:(14)rjtFMMELπ∗,δjt=Eπ∗EθjtxjtLjtFMMELθjt,δjt.Int, the Bayes risk function of first-aid MMEL in all demand points is(15)rFMMELπ∗,δ=Eπ∗EθjtxjtLFMMELθ,δ.Definition 6.
The fair allocation principle under demand information updating ensures that the Bayes risks of first-aid MMEL in all disaster points are the same; that is,(16)rjtFMMELπ∗,δjt=1n·rFMMELπ∗,δ.


### 4.3. Decision Model under Demand Information Updating

The optimal decision scheme can be obtained by minimizing the Bayes risk. Based on literature (Berger, 1980), when the posterior distribution isπ(θ∣x), the posterior Bayes expected loss of action a is(17)ρπθ∣x,a=∫ΘLθ,adFπθ∣xθ.Theorem 7 (Berger, 1980).
The posterior Bayes action obtained after minimizing (17) is equal to the Bayes action obtained after minimizing (18)∫ΘLθ,afx∣θdFπθ.Theorem 8 (Berger, 1980).
Whenδ is a nonrandomized estimator,(19)rπ,θ=∫x:mx>oρπθ∣x,δxdFmx.The integral sign in Theorems7 and 8 shows that when random variables are continuous, they are calculated by an integer; when random variables are discrete, they are calculated by a summation formula. Theorems 7 and 8 are applied to (16) to simplify its calculation.On the basis of the preceding discussion, the equilibrium decision model without information updating can be improved. The model discussed in this paragraph is the equilibrium decision model under demand information updating. Thus, the target function is (6), with unchanged constraints (8), (9), and (10). The related improvement is mainly for constraint (7), and the specific improvement is as follows:(20)∫Δ∫ΘLjtFMMELθjt,δjtftxjt∣θjtdFπθjtdFuxjtΔ=xjt:uxjt>o,=∫Δ∫ΘLFMMELθ,δftx∣θdFπθdFuxΔ=x:ux>0.Meanwhile, the Bayes-risk-related calculation formula is added as follows:(21)uxjt=∫Θfxjt∣θjtdFπθjt,∀j∈J,(22)hxjt,θjt=πθjtfxjt∣θjt,∀j∈J,(23)π∗θjt∣xjt=hxjt,θjtuxjt,∀j∈J.Equation (20) refers to the mathematical expression of the fair allocation principle under demand information updating, (21) is the marginal probability density formula, (22) is a joint probability density formula, and (23) is the posterior probability density formula.


## 4.1. Symbols and Hypotheses

The decision of earthquake first-aid medicine allocation cannot be in place in just one step but is continuously updated with the development of an earthquake disaster. Such a decision is a dynamic allocation based on the related information updates. The balanced allocation of earthquake first-aid medicine based on the demand information update can be modeled using Bayesian analysis [15], which is used to establish related models. The symbols are presented in Notations: (ii) New Symbols.In an earthquake first-aid medicine allocation model under demand information updating, a decision maker can adopt various rescue medicine allocation plan setsA={a} (the plan mainly refers to the quantity of medicine distributed from all supply points to all disaster points) to describe various possible status sets of population transfer rate Θ={θ}. The decision-making rule δ(x) shows the action taken [δ(x) refers to the function from X to A] and X=x is the observation value of the sample information. The loss and risk functions are based on the following hypotheses.Hypothesis 4.
The density function of population transfer rate in prior distributionπ(θ) is known, and the sample information related to θ can be continuously observed and obtained in cycle T. The observation can be conducted numerous times in cycle T. The comprehensive observation value is obtained.The Bayes decision method comprehensively uses the prior information of parameters and sample information to make a decision. The density function of disaster information in prior distributionπ(θ) can be derived, where θ is a multidimensional random variable, based on the disaster databases from government departments (e.g., the National Climatic Data Center of the National Oceanic and Atmospheric Administration, the National Geophysical Data Center, and the National Oceanographic Data Center).The population transfer rate information observation is expressed byX1,X2,…. Xk=(X1k,X2k,…) is the sample observed in stage k, and xk refers to the sample information. The condition density function is fk(xk∣θk), which is defined in Xk=X1×X2×⋯×Xk. The posterior probability density while θk is in stage k is πk∗(θk∣xk). δ(xk)=δ(x1k,x2k,…) is the allocation plan for all disaster points.Hypothesis 5.
A difference exists between the emergency supply allocation and actual demand, causing a decrease in the utility of the limited material, which is defined as the first-aid medicine mismatch error losses (MMELs).A difference always exists between first-aid medicine allocation and actual demand because of the incompleteness of information. Excess or deficiency will reduce the effectiveness of drug allocation. The effectiveness of drug allocation will be significant if the allocation error is high. If the demand for resources is set tod(θ) and the supply is s, the first-aid MMEL function can be expressed as follows:(11)LFDMELθ,δ=dθ-sβ,β∈R.


## 4.2. Fair Allocation Principle Based on Demand Information Updating

Equation (11) indicates that the first-aid MMEL function in demand point j in t is(12)LjtFMMELθjt,δjt=dθjt-sijtβ,β∈R.Int, the first-aid MMEL function of all demand points is(13)LFMMELθ,δ=∑j=1ndθjt-∑i=1msitβ,β∈R.Equations (12) and (13) do not consider the observation sample information Xk. Consequently, the effectiveness of decision is affected and the updated data cannot be used. According to the formula of Bayes risk defined in literature (Berger, 1980), the Bayes risk function of first-aid MMEL in demand point j in t is defined as follows:(14)rjtFMMELπ∗,δjt=Eπ∗EθjtxjtLjtFMMELθjt,δjt.Int, the Bayes risk function of first-aid MMEL in all demand points is(15)rFMMELπ∗,δ=Eπ∗EθjtxjtLFMMELθ,δ.Definition 6.
The fair allocation principle under demand information updating ensures that the Bayes risks of first-aid MMEL in all disaster points are the same; that is,(16)rjtFMMELπ∗,δjt=1n·rFMMELπ∗,δ.


## 4.3. Decision Model under Demand Information Updating

The optimal decision scheme can be obtained by minimizing the Bayes risk. Based on literature (Berger, 1980), when the posterior distribution isπ(θ∣x), the posterior Bayes expected loss of action a is(17)ρπθ∣x,a=∫ΘLθ,adFπθ∣xθ.Theorem 7 (Berger, 1980).
The posterior Bayes action obtained after minimizing (17) is equal to the Bayes action obtained after minimizing (18)∫ΘLθ,afx∣θdFπθ.Theorem 8 (Berger, 1980).
Whenδ is a nonrandomized estimator,(19)rπ,θ=∫x:mx>oρπθ∣x,δxdFmx.The integral sign in Theorems7 and 8 shows that when random variables are continuous, they are calculated by an integer; when random variables are discrete, they are calculated by a summation formula. Theorems 7 and 8 are applied to (16) to simplify its calculation.On the basis of the preceding discussion, the equilibrium decision model without information updating can be improved. The model discussed in this paragraph is the equilibrium decision model under demand information updating. Thus, the target function is (6), with unchanged constraints (8), (9), and (10). The related improvement is mainly for constraint (7), and the specific improvement is as follows:(20)∫Δ∫ΘLjtFMMELθjt,δjtftxjt∣θjtdFπθjtdFuxjtΔ=xjt:uxjt>o,=∫Δ∫ΘLFMMELθ,δftx∣θdFπθdFuxΔ=x:ux>0.Meanwhile, the Bayes-risk-related calculation formula is added as follows:(21)uxjt=∫Θfxjt∣θjtdFπθjt,∀j∈J,(22)hxjt,θjt=πθjtfxjt∣θjt,∀j∈J,(23)π∗θjt∣xjt=hxjt,θjtuxjt,∀j∈J.Equation (20) refers to the mathematical expression of the fair allocation principle under demand information updating, (21) is the marginal probability density formula, (22) is a joint probability density formula, and (23) is the posterior probability density formula.


## 5. Numerical Simulation Analysis

### 5.1. Background and Data

A shallow-source great earthquake at magnitude 8.0 suddenly occurred at 2:28:04 PM, China Standard Time, on May 12, 2008, in WC and Beichuan (BC), Sichuan Province. The maximum intensity was 11. This earthquake was the strongest earthquake in the area with the greatest destruction since the founding of new China. A total of 417 counties, 4,656 towns, and 47,789 villages in Sichuan, Gansu, Shanxi, Chongqing, and other provinces (districts and cities) were affected by the earthquake. Among these areas, 10 disaster areas in the Sichuan Province were seriously affected. The WC earthquake disaster situation is shown in Table1. The total casualties in the 10 severe disaster areas are shown in Table 2. In consideration of the significant impact of the WC earthquake and data availability, a simulation analysis of the first-aid medicine allocation in the 10 severe disaster areas was conducted by using the earthquake first-aid medicine equilibrium decision model under an uncertain environment. Figure 1 shows that the rescue supplies across the country were first distributed to these four points and then sent to WC, BC, Mianzhu (MZ), Shifang (SF), Qingchuan (QC), Maoxian (MX), Anxian (AX), Dujiangyan (DJY), Pingwu (PW), Pengzhou (PZ), and other disaster areas with Chengdu (CD), Deyang (DY), Mianyang (MY), and Guangyuan (GY) as the rescue points.Table 1
Victims and casualties in the WC earthquake (unit: people, square kilometer).
Scope Death Missing Injured Transferred population Affected population Affected area Country 69,181 18,498 374,171 1,510,000 above 46,260,000 440,000 Sichuan 68,669 18,498 360,352 1,460,000 above 29,830,000 280,000 Data resources: report on the 2008 Sichuan big earthquake.Table 2
Total casualties and disaster index in 10 severe disaster areas.
Disaster points WC BC MZ SF QC MX AX DJY PW PZ Population (10 thousand people) 11.2 16.1 51.6 43.3 25.0 10.4 48.4 62.2 18.6 77.0 Death/missing (people) 23,871 20,047 11,380 6132 4819 4088 3295 3388 6565 1131 Injured (people) 34,583 9693 36,468 31,990 15,453 8183 13,476 4388 32,145 5770 Average intensity 8.89 9.16 9.14 8.68 8.74 7.91 8.89 9.13 8.15 8.53 Comprehensive disaster index 0.8675 0.7050 0.6612 0.5953 0.5146 0.5107 0.4993 0.4910 0.4424 0.4333 Data resources: Ministry of Civil Affairs, Ministry of Finance, Ministry of Science and Technology, China Earthquake Administration, and China National Commission for Disaster Reduction.Figure 1
Distribution of severe disaster areas and main rescue points in the WC earthquake. Data resources: China Seismological Bureau, National Wenchuan Earthquake Reconstruction Planning Group, Google Earth.In the most severely affected areas in the WC earthquake, the main highway transportation tools participating in material distribution were various trucks, and the main air tools were helicopters (Black Hawk, Zhi 8, Mi 17 series, Mi 26, etc.). In consideration of the importance of first-aid medicine, a helicopter is used for transport in the simulation analysis. The driving distance among all nodes in the distribution network and the driving time are shown in Table3. The medicine supplies in all supply points are shown in Table 4.Table 3
Driving distance among all nodes in the distribution network and the driving time.
WC BC MZ SF QC MX AX DJY PW PZ CD 100/35 135/45 76/28 53/21 240/77 116/40 108/37 57/22 200/65 40/17 DY 85/30 78/28 30/14 22/12 180/59 80/29 48/20 76/28 143/48 45/19 MY 103/36 45/19 48/20 62/24 135/46 82/30 13/10 115/40 105/37 86/31 GY 238/76 146/49 198/65 215/70 60/23 205/67 156/56 265/85 124/42 240/77 Note: the unit is km/minute; the flying speed is 200 km/h.Table 4
Driving distance among all nodes in the distribution network and the driving time.
Supply points CD DY MY GY Quantity of drug A 35000 14000 13999 6999


### 5.2. Equilibrium Time Point between the Risks of Decision Error and Delay

As time passes, the forecast accuracy increases, the risk of decision error decreases, and the risk of decision delay increases. Therefore, the calculation formula of forecast accuracyp(t) adopted in this study is as follows:(24)pt=1T0-t+1.The calculation formula of the risk of decision error is(25)RtError=1-pt=1-1T0-t+1.The calculation formula of the risk of decision delay is(26)RtDelay=Rt=tT0,where 0<t<T0 and T0 refers to the length of decision cycle, which can be set to 24 h. The calculation formulas of forecast accuracy (24), the risk function of decision error (25), and the risk of decision delay (26) can be modified based on the specific conditions. According to the above formulas, when decision time t∈[14.1903,22.3909] can be obtained as α=0.1, the risk of decision error balances the risk of decision delay. When decision time t=19.5756, the difference between the risks of decision error and delay is 0. The function graphs of the risk of decision error, risk of decision delay, and the balance between them are shown in Figure 2.Figure 2
Function graphs of risk of decision error, risk of decision delay, and the balance between them.


### 5.3. Distribution of Random Variables for Demand and Setting of Other Calculation Formulas

The calculation formula of the demand in all disaster areas is(27)djt=dθjt=Nj·γj·λ·θjt,where Nj refers to the population in disaster point j, γj refers to the casualty ratio in disaster point j, and λ refers to the demand per casualty for a certain drug.The historical seismic data published by the China National Committee for Disaster Reduction and the Ministry of Science and Technology can be used as reference to check if the various status sets of population transfer rateΘ={θ} obey the distribution by population transfer rate and road destructive rate in the 10 severe disaster areas and 36 severe disaster counties in the WC earthquake. The data are shown as follows.Population Transfer Rate. 1.0223, 0.8625, 0.9029, 0.8529, 0.9811, 1.2824, 0.9720, 0.7174, 0.7767, 1.0008, 0.5150, 0.8476, 0.5229, 0.4261, 0.7970, 0.5351, 0.6398, 0.8000, 0.1151, 0.2666, 0.5037, 0.4656, 0.5241, 0.5811, 0.3169, 0.3687, 0.1085, 0.1375, 0.2570, 0.663, 0.3151, 0.2138, 0.9789, 0.5754, 1.2792, 0.5511, 0.4385, 0.8119, 0.3401, 0.1708, 0.3565, 0.9075, 0.5313, 0.5220, 0.3713, 0.3350, —, —.On the basis of the above data, SPSS is used to make a statistical analysis, as shown in Figure3 and Table 5. The examination result is shown in Table 6. The population transfer rate θjt obeys the following distribution in [0,1]:(28)θjt·Nμ,σ2,∀j∈J,t∈T.Table 5
Statistical table of population transfer rate.
Description Population transfer rate N Valid 46 Missing 0 Mean .596971739 Std. deviation .3013077317 Skewness .362 Std. error of skewness .350 Kurtosis −.546 Std. error of kurtosis .688Table 6
One-sample Kolmogorov-Smirnov test result.
Description Population transfer rate N 46 Normal parameters (a, b) Mean .596971739 Std. deviation .3013077317 Most extreme differences Absolute .108 Positive .108 Negative −.073 Kolmogorov-SmirnovZ .732 Asymp. sig. (2-tailed) .657 a: the test distribution is normal. b: calculated from data.Figure 3
Frequency histogram of population transfer rate.Table7 shows the statistics of sample data under the assumption that the population transfer rate of this earthquake in all severe disaster areas is sampled.Table 7
Statistical value of population transfer rate sample data in severe disaster areas.
Demand points WC BC MZ SF QC MX AX DJY PW PZ Total Mean value 0.61 0.53 0.45 0.47 0.51 0.48 0.46 0.57 0.53 0.41 0.502 Variance 0.09 0.10 0.12 0.09 0.11 0.09 0.12 0.10 0.11 0.09 0.10


### 5.4. Calculation Result

On the basis of the above settings, letγj=0.02 (j=1,2,…,10), λ=2,β=2; then the equilibrium decision method of the first-aid medicine allocation without information updating can be used to obtain the decision-making time t∈[14.1903,22.3909]. The detailed allocation scheme is shown in Table 8.Table 8
Allocation scheme of drug A.
WC BC MZ SF QC MX AX DJY PW PZ CD 2145 0 0 6055 0 0 0 11939  14861 DY 0 0 9704 2426  1870 0 0 0 0 MY 0 3368 0 0 0 169 9120 0 1342 0 GY 0 0 0 0 4786 0 0 0 2213 0  Total 2145 3368 9704 8481 4786 2039 9120 11939 3555 14861The equilibrium decision method of the first-aid medicine allocation under demand information updating is used to obtain the decision-making timet∈[14.1903,22.3909]. The detailed allocation scheme is shown in Table 9.Table 9
Allocation scheme of drug A.
WC BC MZ SF QC MX AX DJY PW PZ CD 3311 0 0 7884 0 0 0 11435  12370 DY 0 0 12384 187  1429 0 0 0 0 MY 0 4054 0 0 0 545 8653 0 747 0 GY 0 0 0 0 4578 0 0 0 2421 0  Total 3311 4054 12384 8071 4578 1974 8653 11435 3168 12370Meanwhile, we can compare the population, allocation amount without information updating, allocation amount, and population multiplied by the comprehensive disaster index in each disaster area of the total. The comparison results are shown in Figure4. The allocation amount without information updating fits the population in each disaster area well. By contrast, the allocation amount with information updating fits the population multiplied by the comprehensive disaster index in each disaster area better. Therefore, the decision method with demand information updating better fits the real disaster situation based on historical and sample information.Figure 4
Comparison of the population and allocation amount.


## 5.1. Background and Data

A shallow-source great earthquake at magnitude 8.0 suddenly occurred at 2:28:04 PM, China Standard Time, on May 12, 2008, in WC and Beichuan (BC), Sichuan Province. The maximum intensity was 11. This earthquake was the strongest earthquake in the area with the greatest destruction since the founding of new China. A total of 417 counties, 4,656 towns, and 47,789 villages in Sichuan, Gansu, Shanxi, Chongqing, and other provinces (districts and cities) were affected by the earthquake. Among these areas, 10 disaster areas in the Sichuan Province were seriously affected. The WC earthquake disaster situation is shown in Table1. The total casualties in the 10 severe disaster areas are shown in Table 2. In consideration of the significant impact of the WC earthquake and data availability, a simulation analysis of the first-aid medicine allocation in the 10 severe disaster areas was conducted by using the earthquake first-aid medicine equilibrium decision model under an uncertain environment. Figure 1 shows that the rescue supplies across the country were first distributed to these four points and then sent to WC, BC, Mianzhu (MZ), Shifang (SF), Qingchuan (QC), Maoxian (MX), Anxian (AX), Dujiangyan (DJY), Pingwu (PW), Pengzhou (PZ), and other disaster areas with Chengdu (CD), Deyang (DY), Mianyang (MY), and Guangyuan (GY) as the rescue points.Table 1
Victims and casualties in the WC earthquake (unit: people, square kilometer).
Scope Death Missing Injured Transferred population Affected population Affected area Country 69,181 18,498 374,171 1,510,000 above 46,260,000 440,000 Sichuan 68,669 18,498 360,352 1,460,000 above 29,830,000 280,000 Data resources: report on the 2008 Sichuan big earthquake.Table 2
Total casualties and disaster index in 10 severe disaster areas.
Disaster points WC BC MZ SF QC MX AX DJY PW PZ Population (10 thousand people) 11.2 16.1 51.6 43.3 25.0 10.4 48.4 62.2 18.6 77.0 Death/missing (people) 23,871 20,047 11,380 6132 4819 4088 3295 3388 6565 1131 Injured (people) 34,583 9693 36,468 31,990 15,453 8183 13,476 4388 32,145 5770 Average intensity 8.89 9.16 9.14 8.68 8.74 7.91 8.89 9.13 8.15 8.53 Comprehensive disaster index 0.8675 0.7050 0.6612 0.5953 0.5146 0.5107 0.4993 0.4910 0.4424 0.4333 Data resources: Ministry of Civil Affairs, Ministry of Finance, Ministry of Science and Technology, China Earthquake Administration, and China National Commission for Disaster Reduction.Figure 1
Distribution of severe disaster areas and main rescue points in the WC earthquake. Data resources: China Seismological Bureau, National Wenchuan Earthquake Reconstruction Planning Group, Google Earth.In the most severely affected areas in the WC earthquake, the main highway transportation tools participating in material distribution were various trucks, and the main air tools were helicopters (Black Hawk, Zhi 8, Mi 17 series, Mi 26, etc.). In consideration of the importance of first-aid medicine, a helicopter is used for transport in the simulation analysis. The driving distance among all nodes in the distribution network and the driving time are shown in Table3. The medicine supplies in all supply points are shown in Table 4.Table 3
Driving distance among all nodes in the distribution network and the driving time.
WC BC MZ SF QC MX AX DJY PW PZ CD 100/35 135/45 76/28 53/21 240/77 116/40 108/37 57/22 200/65 40/17 DY 85/30 78/28 30/14 22/12 180/59 80/29 48/20 76/28 143/48 45/19 MY 103/36 45/19 48/20 62/24 135/46 82/30 13/10 115/40 105/37 86/31 GY 238/76 146/49 198/65 215/70 60/23 205/67 156/56 265/85 124/42 240/77 Note: the unit is km/minute; the flying speed is 200 km/h.Table 4
Driving distance among all nodes in the distribution network and the driving time.
Supply points CD DY MY GY Quantity of drug A 35000 14000 13999 6999


## 5.2. Equilibrium Time Point between the Risks of Decision Error and Delay

As time passes, the forecast accuracy increases, the risk of decision error decreases, and the risk of decision delay increases. Therefore, the calculation formula of forecast accuracyp(t) adopted in this study is as follows:(24)pt=1T0-t+1.The calculation formula of the risk of decision error is(25)RtError=1-pt=1-1T0-t+1.The calculation formula of the risk of decision delay is(26)RtDelay=Rt=tT0,where 0<t<T0 and T0 refers to the length of decision cycle, which can be set to 24 h. The calculation formulas of forecast accuracy (24), the risk function of decision error (25), and the risk of decision delay (26) can be modified based on the specific conditions. According to the above formulas, when decision time t∈[14.1903,22.3909] can be obtained as α=0.1, the risk of decision error balances the risk of decision delay. When decision time t=19.5756, the difference between the risks of decision error and delay is 0. The function graphs of the risk of decision error, risk of decision delay, and the balance between them are shown in Figure 2.Figure 2
Function graphs of risk of decision error, risk of decision delay, and the balance between them.


## 5.3. Distribution of Random Variables for Demand and Setting of Other Calculation Formulas

The calculation formula of the demand in all disaster areas is(27)djt=dθjt=Nj·γj·λ·θjt,where Nj refers to the population in disaster point j, γj refers to the casualty ratio in disaster point j, and λ refers to the demand per casualty for a certain drug.The historical seismic data published by the China National Committee for Disaster Reduction and the Ministry of Science and Technology can be used as reference to check if the various status sets of population transfer rateΘ={θ} obey the distribution by population transfer rate and road destructive rate in the 10 severe disaster areas and 36 severe disaster counties in the WC earthquake. The data are shown as follows.Population Transfer Rate. 1.0223, 0.8625, 0.9029, 0.8529, 0.9811, 1.2824, 0.9720, 0.7174, 0.7767, 1.0008, 0.5150, 0.8476, 0.5229, 0.4261, 0.7970, 0.5351, 0.6398, 0.8000, 0.1151, 0.2666, 0.5037, 0.4656, 0.5241, 0.5811, 0.3169, 0.3687, 0.1085, 0.1375, 0.2570, 0.663, 0.3151, 0.2138, 0.9789, 0.5754, 1.2792, 0.5511, 0.4385, 0.8119, 0.3401, 0.1708, 0.3565, 0.9075, 0.5313, 0.5220, 0.3713, 0.3350, —, —.On the basis of the above data, SPSS is used to make a statistical analysis, as shown in Figure3 and Table 5. The examination result is shown in Table 6. The population transfer rate θjt obeys the following distribution in [0,1]:(28)θjt·Nμ,σ2,∀j∈J,t∈T.Table 5
Statistical table of population transfer rate.
Description Population transfer rate N Valid 46 Missing 0 Mean .596971739 Std. deviation .3013077317 Skewness .362 Std. error of skewness .350 Kurtosis −.546 Std. error of kurtosis .688Table 6
One-sample Kolmogorov-Smirnov test result.
Description Population transfer rate N 46 Normal parameters (a, b) Mean .596971739 Std. deviation .3013077317 Most extreme differences Absolute .108 Positive .108 Negative −.073 Kolmogorov-SmirnovZ .732 Asymp. sig. (2-tailed) .657 a: the test distribution is normal. b: calculated from data.Figure 3
Frequency histogram of population transfer rate.Table7 shows the statistics of sample data under the assumption that the population transfer rate of this earthquake in all severe disaster areas is sampled.Table 7
Statistical value of population transfer rate sample data in severe disaster areas.
Demand points WC BC MZ SF QC MX AX DJY PW PZ Total Mean value 0.61 0.53 0.45 0.47 0.51 0.48 0.46 0.57 0.53 0.41 0.502 Variance 0.09 0.10 0.12 0.09 0.11 0.09 0.12 0.10 0.11 0.09 0.10


## 5.4. Calculation Result

On the basis of the above settings, letγj=0.02 (j=1,2,…,10), λ=2,β=2; then the equilibrium decision method of the first-aid medicine allocation without information updating can be used to obtain the decision-making time t∈[14.1903,22.3909]. The detailed allocation scheme is shown in Table 8.Table 8
Allocation scheme of drug A.
WC BC MZ SF QC MX AX DJY PW PZ CD 2145 0 0 6055 0 0 0 11939  14861 DY 0 0 9704 2426  1870 0 0 0 0 MY 0 3368 0 0 0 169 9120 0 1342 0 GY 0 0 0 0 4786 0 0 0 2213 0  Total 2145 3368 9704 8481 4786 2039 9120 11939 3555 14861The equilibrium decision method of the first-aid medicine allocation under demand information updating is used to obtain the decision-making timet∈[14.1903,22.3909]. The detailed allocation scheme is shown in Table 9.Table 9
Allocation scheme of drug A.
WC BC MZ SF QC MX AX DJY PW PZ CD 3311 0 0 7884 0 0 0 11435  12370 DY 0 0 12384 187  1429 0 0 0 0 MY 0 4054 0 0 0 545 8653 0 747 0 GY 0 0 0 0 4578 0 0 0 2421 0  Total 3311 4054 12384 8071 4578 1974 8653 11435 3168 12370Meanwhile, we can compare the population, allocation amount without information updating, allocation amount, and population multiplied by the comprehensive disaster index in each disaster area of the total. The comparison results are shown in Figure4. The allocation amount without information updating fits the population in each disaster area well. By contrast, the allocation amount with information updating fits the population multiplied by the comprehensive disaster index in each disaster area better. Therefore, the decision method with demand information updating better fits the real disaster situation based on historical and sample information.Figure 4
Comparison of the population and allocation amount.


## 6. Conclusion

Unlike the allocation of general first-aid supplies, earthquake first-aid medicine allocation is specific. The transport constraint is relatively small and the scientific nature of the allocation scheme must be high. In this study, according to the change in the supply and demand for first-aid drugs, the balance between the risks of decision error and delay reflects that the square of their difference is not beyond the constraint on a certain expected value. The expected demand satisfaction rates of all disaster areas should be the same, and the Bayes risk of first-aid medicine allocation loss in all demand points should be equal to the total Bayes risk to achieve fair allocation to all disaster areas. This study also introduces a set of equilibrium decision methods for earthquake first-aid medicine allocation based on information updating to minimize the total transport time loss. Finally, the simulation analysis is applied to a certain drug allocation in the WC earthquake to verify the effectiveness of the constructed method. This study can provide enlightenment on the following management decisions for emergency decision makers:(1)
In earthquake first-aid medicine allocation decision, the relation between decision-making time and risk is important. A decision maker must realize the balance between the risks of decision error and delay to realize the maximization of allocation efficiency.(2)
The decision-making information in an earthquake disaster environment is incomplete. Bayes analysis is an effective method for realizing the integration of the related history and sample information.However, the method proposed in this study has numerous limitations. Future studies on earthquake first-aid medicine allocation can comprehensively consider the changes in supply and demand, the balance between the risks of decision error and delay under information updating, time window constraint, and comprehensive allocation in a time-space network.



---

*Source: 6326938-2017-01-31.xml*