#### 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.

#### 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.

#### 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 that supply points and 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 point in can be expressed by the population transfer rate. The requirement is denoted by the stochastic variable and submits to the detailed distribution. The probability density function is .

*Hypothesis 3. *The supply in the supply point is confirmed.

*Definition 1. *The requirement of point in is the function of population transfer rate , which is defined as follows:

*Definition 2. *The loss of decision error in is a function related to the prediction accuracy , which is defined as follows:where is a monotone increasing function in decision time .

*Definition 3. *The loss of decision delay in is a monotone increasing function related to decision time, which is defined as follows:

*Definition 4. *The loss of transport time in is the total time of medicine distribution from each supply point to each disaster point, which is defined as follows:

*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 point in is defined as follows:In (5), when , the expected demand satisfaction rate is less than 1 because the supply is not enough. However, when , the expected demand satisfaction rate 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:Thus,

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 sets (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 shows the action taken [ refers to the function from to ] and 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 . The observation can be conducted numerous times in cycle . 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 by . is the sample observed in stage , and refers to the sample information. The condition density function is , which is defined in . The posterior probability density while is in stage is . 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 to and the supply is , the first-aid MMEL function can be expressed as follows:

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

Equation (11) indicates that the first-aid MMEL function in demand point in is

In , the first-aid MMEL function of all demand points is

Equations (12) and (13) do not consider the observation sample information . 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 in is defined as follows:

In , the Bayes risk function of first-aid MMEL in all demand points is

*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,

##### 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 , the posterior Bayes expected loss of action is

Theorem 7 (Berger, 1980). *The posterior Bayes action obtained after minimizing (17) is equal to the Bayes action obtained after minimizing *

Theorem 8 (Berger, 1980). *When is a nonrandomized estimator,*

The integral sign in Theorems 7 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:Meanwhile, the Bayes-risk-related calculation formula is added as follows:

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 Table 1. 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.

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 Table 3. The medicine supplies in all supply points are shown in Table 4.

##### 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 accuracy adopted in this study is as follows:

The calculation formula of the risk of decision error is

The calculation formula of the risk of decision delay iswhere and 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 can be obtained as , the risk of decision error balances the risk of decision delay. When decision time , 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.

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

The calculation formula of the demand in all disaster areas iswhere refers to the population in disaster point , refers to the casualty ratio in disaster point , 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 Figure 3 and Table 5. The examination result is shown in Table 6. The population transfer rate obeys the following distribution in :

Table 7 shows the statistics of sample data under the assumption that the population transfer rate of this earthquake in all severe disaster areas is sampled.

##### 5.4. Calculation Result

On the basis of the above settings, let (), ; then the equilibrium decision method of the first-aid medicine allocation without information updating can be used to obtain the decision-making time . The detailed allocation scheme is shown in Table 8.

The equilibrium decision method of the first-aid medicine allocation under demand information updating is used to obtain the decision-making time . The detailed allocation scheme is shown in Table 9.

Meanwhile, 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 Figure 4. 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.

#### 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.

#### Notations

*(i) Symbols*

*Set*

: | Supply point set, |

: | Disaster point set, |

*Stochastic Variables*

: | Population transfer rate |

: | Population transfer rate in disaster point in |

*Variables*

: | Requirement of point in |

: | Demand prediction accuracy in |

: | Loss of error decision in |

: | Loss of delay decision in |

: | Loss of transportation time of the decision-making program in |

: | Expected demand satisfaction rate of point in |

*Decision Variables*

: | Decision time |

: | Requirement of point in from supply point |

*Parameters*

: | Transport time in point from supply point |

: | Requirement in point in . |

*(ii) New Symbols*

*Stochastic Variables*

: | Independent observation value of |

: | Specific realization value of |

: | Specific allocation plan of rescue drugs. |

*Set*

: | All possible set |

: | All possible set |

: | All possible set |

: | All positive integer set. |

*Functions and Equations*

: | Prior probability density function of |

: | Edge probability density function of |

: | Joint probability density function of and |

: | Condition probability density function of |

: | Posterior probability density function of |

: | Cumulative distribution function of |

: | Bayesian risk of loss of rescue drug allocation error in |

: | Bayesian risk of loss of rescue drug allocation error in disaster point in |

: | Marginal cumulative distribution function of |

: | Decision in space is a plan under the condition of . |

#### Competing Interests

The authors declare that they have no competing interests.

#### Acknowledgments

This work is supported by the National Natural Science Foundation of China (Grant no. 71601146), Humanities and Social Sciences Foundation of Ministry of Education of China (Grant no. 15YJCZH211), Zhejiang Provincial Philosophy and Social Science Program of China (Grant no. 16ZJQN025YB), Zhejiang Provincial Natural Science Foundation of China (Grant no. LQ16G010005), and Zhejiang Provincial College Major Program of Humanities and Social Science of China (Grant no. 2014QN006).