# A Multiobjective Optimization Model for Prevention and Control of Coronavirus Disease 2019 in China

Authors: Xiaocheng Li; Zhaoli Liu; Fangzhen Ge
Journal: Discrete Dynamics in Nature and Society
 (2022)

---

## Abstract

It is a global issue to set up a practical, sensitive, and useful model to eradicate or mitigate the coronavirus disease 2019 (COVID-19). Taking Central China’s Hubei Province for example, three models were established. Firstly, a susceptible-probable-infectious-recovered (SPIR) model was proposed to predict the monthly number of confirmed and susceptible cases in each city. Next, an epidemic prefecture clustering model was set up to find proper vaccine delivery sites, according to the distance of each city. Finally, a dynamic material delivery optimization model was established for multiple epidemic prefectures, aiming to speed up vaccine production and storage in each delivery site.


---

## Body

## 1. Introduction

Coronavirus disease 2019 (COVID-19) is a serious and fatal disease, which threatens the public health of the entire world. Some susceptible cases were found in Wuhan, China, as early as December 2019 [1]. On January 12, 2020, the World Health Organization (WHO) named this highly contagious disease 2019-novel coronavirus (2019-nCoV) [2] and announced that the disease has no specific medicine, but routine antiviral treatment. Therefore, it is a global issue to set up a practical, sensitive, and useful model to eradicate or mitigate the disease. Otherwise, COVID-19 will upgrade and spread like Ebola, put the fundamental survival of the society at stake, and might even paralyze countries. China’s experience in upgrading personal protective equipment, building mobile labs, mobilizing medical staff, and real-time tracking of epidemic spread has played an important role in COVID-19 prevention and control [3].Li and Bai [4] established an infectious disease model constrained by nonlinear incidence and maximum treatment capacity. Comparing recent cumulative and predicted cases, Kumar and Bhagavathula [5] forecasted that the confirmed cases, deaths, and recoveries will double in all the observed countries except China, Switzerland, and Germany, and the death and recovery rates will increase faster than confirmed cases in the next two months. Agusto [6] modeled the transmission dynamics of Ebola with relapse and reinfection and discovered that high efficacy control measures can effectively curtail Ebola in the community. Yi et al. [7] proved that the joint prevention between the government, the masses, and hospitals is a scientific way to prevent and control the epidemic, using the five-early model (early detection, early screening, early diagnosis, early quarantine, and early recovery). Kennedy et al. [8] suggested that the magnitude and duration of social distancing required to control the pandemic can be reduced, if people exercise caution in public by protecting themselves. Using the S(usceptible) U(nidentified infectious) I(dentified infectious) H(ospitalized) R(ecovered) model (SUIHR), Zhen et al. [9] justified the extremely strict border control in entities of China. Rui and Tian [10] divided the population into susceptible, infective, recovered, and dead populations and then reveal the law of epidemic transmission by the SIRD model. Wu et al. [11] established a generalized-growth model to present the evolution of the number of the total confirmed cases changing with time.In the light of the above, this article aims at modeling COVID-19 trend based on factors such as spread rate, vaccine demand, delivery site, and vaccine production rate, to name but a few, and optimize the model by bringing in other key factors. Specifically, three models were established to provide the optimal scheme to allay the severe COVID-19 epidemic. Firstly, the susceptible-probable-infectious-recovered (SPIR) model was employed to predict the situation at the outbreak of the disease. Then, the k-means clustering (KMC) model [12–14] was introduced to optimize the delivery sites of vaccines and medical staff. Finally, a dynamic material delivery optimization model was constructed to deliver sufficient vaccines and medical staff in the shortest time and at the minimum cost. Figure 1 is the roadmap of this research.Figure 1
Research roadmap.


## 2. Hypotheses and Symbols

### 2.1. Hypotheses

H1: the objects form an ideal, uniformly mixed group. The total number of objects remains at N, without death caused by immigration, emigration, or other reasons.H2: those infected with COVID-19 will eventually recover and acquire long-term immunity through treatment.H3: the infected are immediately contagious, because COVID-19 is mildly infectious, despite its long incubation period.H4: the objects are divided into four classes: the susceptible, the probable, the infected, and the removed.H5: once the epidemic breaks out in a prefecture, the prefecture will be immediately quarantined, cutting off the flow of population.H6: the storagesf emergency materials in each prefecture have enough vehicles to deliver the materials in time.


### 2.2. Symbols

The main variables involved in this article are shown in Table1.Table 1
Symbols of key variates.
SymbolUnitDescriptionNPersonTotal populations(t)%Percentage of the infected at timet in the total populationi(t)%Percentage of patients at timet in the total populatione(t)%Percentage of the susceptible at timet in the total populationr(t)%Percentage of the recovered at timet in the total populationθ1%Percentage of the susceptible diagnosed as normally infected per dayθ2%Percentage of the susceptible diagnosed as infected per dayh%Remove ratio (the sum of death rate and recover rate per day)λPersonDaily contracting ratio (the mean number of people infected by a patient through contact)Ci$Fixed cost of setting upi storagesCs$Subsidy of medical staffCd$Transport cost of each dose per unit of distancetijdayTravel time fromi to jTkdayTime threshold of vaccine delivery (violating the threshold will bring a penalty cost)MijDoseQuantity of vaccines delivered fromi to jXijkmDistance betweeni and jDjDoseVaccine demand ofjNijGroupNumber of medical groups delivered fromi to jβ—Ratio of the infected to medical staffIPersonNumber of the infectedV0km/hVaccine delivery speed from storages to the epidemic prefectureVpDose/dayVaccine production speedYijPersonNumber of medical staff delivered fromi to jSiDoseThe drug storage capacity of thei-th inventory point


## 2.1. Hypotheses

H1: the objects form an ideal, uniformly mixed group. The total number of objects remains at N, without death caused by immigration, emigration, or other reasons.H2: those infected with COVID-19 will eventually recover and acquire long-term immunity through treatment.H3: the infected are immediately contagious, because COVID-19 is mildly infectious, despite its long incubation period.H4: the objects are divided into four classes: the susceptible, the probable, the infected, and the removed.H5: once the epidemic breaks out in a prefecture, the prefecture will be immediately quarantined, cutting off the flow of population.H6: the storagesf emergency materials in each prefecture have enough vehicles to deliver the materials in time.


## 2.2. Symbols

The main variables involved in this article are shown in Table1.Table 1
Symbols of key variates.
SymbolUnitDescriptionNPersonTotal populations(t)%Percentage of the infected at timet in the total populationi(t)%Percentage of patients at timet in the total populatione(t)%Percentage of the susceptible at timet in the total populationr(t)%Percentage of the recovered at timet in the total populationθ1%Percentage of the susceptible diagnosed as normally infected per dayθ2%Percentage of the susceptible diagnosed as infected per dayh%Remove ratio (the sum of death rate and recover rate per day)λPersonDaily contracting ratio (the mean number of people infected by a patient through contact)Ci$Fixed cost of setting upi storagesCs$Subsidy of medical staffCd$Transport cost of each dose per unit of distancetijdayTravel time fromi to jTkdayTime threshold of vaccine delivery (violating the threshold will bring a penalty cost)MijDoseQuantity of vaccines delivered fromi to jXijkmDistance betweeni and jDjDoseVaccine demand ofjNijGroupNumber of medical groups delivered fromi to jβ—Ratio of the infected to medical staffIPersonNumber of the infectedV0km/hVaccine delivery speed from storages to the epidemic prefectureVpDose/dayVaccine production speedYijPersonNumber of medical staff delivered fromi to jSiDoseThe drug storage capacity of thei-th inventory point


## 3. COVID-19 Modeling

According to China’s National Health Commission (NHC), the nationwide distribution of COVID-19 epidemic on February 1, 2020, was obtained. The number of COVID-19 cases in Chinese provinces is shown in Figure2.Figure 2
Top 10 provinces by the number of COVID-19 cases.As shown in Figure2, Hubei, Zhejiang, and Guangdong led the country in terms of the number of COVID-19 cases. The number of infected people in Hubei province accounted for 63.1 percent of the total, higher than any other province in China. Therefore, Hubei was taken as an example in our research.


### 3.1. Models

As mentioned before, the objects form an ideal, uniformly mixed group. The total number of objects remains at N, without death caused by immigration, emigration, or other reasons. The objects were divided into four groups: the susceptible S, the probable P, the infected I, and the removed R. Figure3 shows the relationship between the four groups.Figure 3
Relationship between different groups.The traditional susceptible, infected, recovered (SIR) model [15–18] can be expressed as(1)didt=λsi−μi,i0=i0,dsdt=−λsi,s0=s0,drdt=μi,r0=r0,where i is the number of the infected person, s is the number of the susceptible person, r is the number of the recovered person, and μ is the recovered ratio.


### 3.2. Models Optimization

The transmission mechanism of COVID-19 differs from that of general infectious diseases, which involve the normally infected, patients, and the removed. COVID-19 has an additional group of objects: the susceptible. This group has similar symptoms with COVID-19, but cannot be diagnosed as infected with COVID-19 or general flus. Some of these objects might be transformed into patients or normally infected. Hence, the susceptible was added to the traditional SIR model:(2)s+p+i+r=1,dsdt=θ1p−λsi,s0=s0,didt=λsi−hi+θ2p,i0=i0,drdt=hi,r0=r0,dedt=−θ1p−θ2p,e0=e0.


### 3.3. Parameter Evaluation

The above model contains the following key parameters: percentage of the susceptible diagnosed as normally infectedθ1 per day; percentage of the susceptible diagnosed as infected θ2 per day; and the sum of death rate and recover rate per day (remove ratio) h. The values of these parameters can be approximated according to the official data released each day:(3)θ1=newly added removal of suspected per daytotal suspected on the current day,θ2=newly added transformation of susceptible to patients per daytotal suspectible on the current day,h=newly added recovered and deaths per daytotal suspectible on the current day.By our research hypotheses, the contracting probability of two persons is closely related to the population size. Therefore, it is assumed that the effective contracting ratio is the reciprocal of the population,λ=1/N.


### 3.4. Model Solution

The epidemic severity in Hubei varies from prefecture to prefecture. It is necessary to optimize the material delivery sites for the province [19]. Before that, the number of COVID-19 patients should be predicted for each prefecture.As shown in Figure4, from 31 December 2019 to 1 February 2020, Wuhan had more infected than any other prefecture in Hubei. Several other prefectures also faced a serious epidemic: Huanggang, Xiaogan, Jingzhou, Ezhou, Suizhou, and Huangshi. In the said period, Wuhan and Jingzhou were the epicenters of the epidemic. Hence, this subsection attempts to predict the COVID-19 situation in each prefecture of Hubei with the SPIR model. However, it is impossible to directly drive S, P, I, and R values with the SPIR model. To solve this problem, the Runge–Kutta method of MATLAB was called to obtain these solutions [21, 22]. And in the Runge–Kutta method, the initial value of parameters, θ1, θ2, h, and λ, are set as 0.001, 0.001, 0.0001, and 0.001, respectively.Figure 4
Epidemic distribution in Hubei [20].The key parameters for Hubei were solved by differential equations (2). The solutions are given in Table 2. Then, the COVID-19 cases in the next month were predicted.Table 2
Values of key parameters.
Parametersθ1θ2hλValues0.00350.00370.00030.00017Figure5 shows the number of diagnosed patients in Hubei. It can be seen that the epidemic reached its peak in the 10th week after the outbreak, as the patients took up 0.035% of the total population. In the 28th week, the epidemic was under control, when the percentage of patients dropped to 0.Figure 5
Number of patients over time.


### 3.5. Results and Analysis

The above model was applied to predict the COVID-19 situation of other preferences in four weeks. The results are given in Tables3 and 4.Table 3
Predicted new susceptible in four weeks.
PrefectureWeek1234TotalWuhan1008.24035.78658.812510.426213.1Huanggang66.8212.4604.4878.61762.1Xiaogan78.0244.5576.31045.71944.5Xiangyang26.491.0257.3335.2709.9Suizhou31.688.4248.6336.8705.5Yichang23.268.5185.5239.1516.3Jingzhou34.6114.3304.6492.0945.5Jingmen24.868.8209.6249.3552.5Ezhou29.688.1280.5362.0760.1Huangshi25.867.8219.5252.2565.2Xianning17.257.0175.3229.6479.0Shiyan15.845.5123.5156.5341.3Xiantao14.045.7116.3168.4344.3Enshi5.018.144.772.9140.6Tianmen9.333.288.0121.2251.6Qianjiang4.616.939.459.0120.0Shennongjia0.30.92.13.46.6Table 4
Predicted new patients in four weeks.
PrefectureWeek134TotalWuhan789.96506.69392.119341Huanggang45.1371.8536.71105.1Xiaogan55.1453.5654.61347.9Xiangyang18.8155.0223.7460.7Suizhou20.0164.7237.7489.5Yichang14.7121.0174.7359.8Jingzhou25.5210.2303.4624.9Jingmen16.0132.1190.7392.8Ezhou21.0172.7249.2513.2Huangshi16.0131.9190.4392.0Xianning13.1108.1156.1321.4Shiyan10.384.8122.4252.0Xiantao8.973.2105.6217.5Enshi3.931.745.894.4Tianmen7.259.085.1175.3Qianjiang3.125.636.975.9Shennongjia0.21.42.04.2


## 3.1. Models

As mentioned before, the objects form an ideal, uniformly mixed group. The total number of objects remains at N, without death caused by immigration, emigration, or other reasons. The objects were divided into four groups: the susceptible S, the probable P, the infected I, and the removed R. Figure3 shows the relationship between the four groups.Figure 3
Relationship between different groups.The traditional susceptible, infected, recovered (SIR) model [15–18] can be expressed as(1)didt=λsi−μi,i0=i0,dsdt=−λsi,s0=s0,drdt=μi,r0=r0,where i is the number of the infected person, s is the number of the susceptible person, r is the number of the recovered person, and μ is the recovered ratio.


## 3.2. Models Optimization

The transmission mechanism of COVID-19 differs from that of general infectious diseases, which involve the normally infected, patients, and the removed. COVID-19 has an additional group of objects: the susceptible. This group has similar symptoms with COVID-19, but cannot be diagnosed as infected with COVID-19 or general flus. Some of these objects might be transformed into patients or normally infected. Hence, the susceptible was added to the traditional SIR model:(2)s+p+i+r=1,dsdt=θ1p−λsi,s0=s0,didt=λsi−hi+θ2p,i0=i0,drdt=hi,r0=r0,dedt=−θ1p−θ2p,e0=e0.


## 3.3. Parameter Evaluation

The above model contains the following key parameters: percentage of the susceptible diagnosed as normally infectedθ1 per day; percentage of the susceptible diagnosed as infected θ2 per day; and the sum of death rate and recover rate per day (remove ratio) h. The values of these parameters can be approximated according to the official data released each day:(3)θ1=newly added removal of suspected per daytotal suspected on the current day,θ2=newly added transformation of susceptible to patients per daytotal suspectible on the current day,h=newly added recovered and deaths per daytotal suspectible on the current day.By our research hypotheses, the contracting probability of two persons is closely related to the population size. Therefore, it is assumed that the effective contracting ratio is the reciprocal of the population,λ=1/N.


## 3.4. Model Solution

The epidemic severity in Hubei varies from prefecture to prefecture. It is necessary to optimize the material delivery sites for the province [19]. Before that, the number of COVID-19 patients should be predicted for each prefecture.As shown in Figure4, from 31 December 2019 to 1 February 2020, Wuhan had more infected than any other prefecture in Hubei. Several other prefectures also faced a serious epidemic: Huanggang, Xiaogan, Jingzhou, Ezhou, Suizhou, and Huangshi. In the said period, Wuhan and Jingzhou were the epicenters of the epidemic. Hence, this subsection attempts to predict the COVID-19 situation in each prefecture of Hubei with the SPIR model. However, it is impossible to directly drive S, P, I, and R values with the SPIR model. To solve this problem, the Runge–Kutta method of MATLAB was called to obtain these solutions [21, 22]. And in the Runge–Kutta method, the initial value of parameters, θ1, θ2, h, and λ, are set as 0.001, 0.001, 0.0001, and 0.001, respectively.Figure 4
Epidemic distribution in Hubei [20].The key parameters for Hubei were solved by differential equations (2). The solutions are given in Table 2. Then, the COVID-19 cases in the next month were predicted.Table 2
Values of key parameters.
Parametersθ1θ2hλValues0.00350.00370.00030.00017Figure5 shows the number of diagnosed patients in Hubei. It can be seen that the epidemic reached its peak in the 10th week after the outbreak, as the patients took up 0.035% of the total population. In the 28th week, the epidemic was under control, when the percentage of patients dropped to 0.Figure 5
Number of patients over time.


## 3.5. Results and Analysis

The above model was applied to predict the COVID-19 situation of other preferences in four weeks. The results are given in Tables3 and 4.Table 3
Predicted new susceptible in four weeks.
PrefectureWeek1234TotalWuhan1008.24035.78658.812510.426213.1Huanggang66.8212.4604.4878.61762.1Xiaogan78.0244.5576.31045.71944.5Xiangyang26.491.0257.3335.2709.9Suizhou31.688.4248.6336.8705.5Yichang23.268.5185.5239.1516.3Jingzhou34.6114.3304.6492.0945.5Jingmen24.868.8209.6249.3552.5Ezhou29.688.1280.5362.0760.1Huangshi25.867.8219.5252.2565.2Xianning17.257.0175.3229.6479.0Shiyan15.845.5123.5156.5341.3Xiantao14.045.7116.3168.4344.3Enshi5.018.144.772.9140.6Tianmen9.333.288.0121.2251.6Qianjiang4.616.939.459.0120.0Shennongjia0.30.92.13.46.6Table 4
Predicted new patients in four weeks.
PrefectureWeek134TotalWuhan789.96506.69392.119341Huanggang45.1371.8536.71105.1Xiaogan55.1453.5654.61347.9Xiangyang18.8155.0223.7460.7Suizhou20.0164.7237.7489.5Yichang14.7121.0174.7359.8Jingzhou25.5210.2303.4624.9Jingmen16.0132.1190.7392.8Ezhou21.0172.7249.2513.2Huangshi16.0131.9190.4392.0Xianning13.1108.1156.1321.4Shiyan10.384.8122.4252.0Xiantao8.973.2105.6217.5Enshi3.931.745.894.4Tianmen7.259.085.1175.3Qianjiang3.125.636.975.9Shennongjia0.21.42.04.2


## 4. KMC Model

Under the condition of ensuring the safety of the epidemic, in order to quickly transport medical workers, medical supplies, and related materials to the affected areas, Hubei was divided into three regions according to the geographical coordinates of the cities and airports. For each region, a material delivery site was established [23, 24]. Once the epidemic breaks out, the local site will transfer vaccines and medical staff to the affected prefectures. If the local force is insufficient to control the epidemic, the other prefectures will transfer vaccines and medical staff to this prefecture.


### 4.1. KMC

The distance from pointi to point j can be defined as(4)dijq=∑k=1pxik−xjkq1/q.Ifq=2, then the Euclidean distance dij2=∑k=1pxik−xjk21/2.The general steps of KMC [25, 26] are given in Figure 6:Step 1: select a number of sample points as cluster heads by experience or basic knowledge of dataStep 2: compute the Euclidean distance between each sample point and each cluster head, and allocate the sample point to the cluster of the nearest cluster headStep 3: classify one sample point at a time, and compute the mean of each cluster to update the cluster headStep 4: repeat Step 2 with the new cluster head, and classify the sample points againStep 5: repeat Steps 2–4 cyclically and adjust the classes until no cluster head changes to obtain the final classesFigure 6
Flow of KMC.


### 4.2. Solution of the Model

Table5 sums up the coordinates of the central cities of the prefectures in Hubei.Table 5
Coordinates of central cities of the prefectures in Hubei.
PrefecturexyWuhan21.832.04Huanggang24.871.75Xiaogan20.003.94Xiangyang11.549.84Suizhou17.498.20Yichang8.312.28Jingzhou12.510.29Jingmen11.914.02Ezhou24.261.16Huangshi25.610Xianning22.20−2.12Shiyan5.6612.96Xiantao18.100.45Enshi00Tianmen16.912.25Qianjiang15.450.79Shennongjia4.768.07According to the clustering results, Hubei was divided into three regions. The first region includes Yichang, Jingzhou, and Enshi; the second contains Wuhan, Huanggang, Xiaogan, Suizhou, Ezhou, Huangshi, Xianning, Xiantao, Tianmen, and Qianjiang; and the third covers Xiangyang, Jingmen, Shiyan, and Shennongjia.The coordinates of the head of each cluster were solved and taken as the location of the delivery site in the corresponding region, because the head is the closest point to any sample point in the cluster. The coordinates of the three heads were Yichang (8.31, 2.28), Wuhan (21.83, 2.04), and Xiangyang (11.54, 9.84).


## 4.1. KMC

The distance from pointi to point j can be defined as(4)dijq=∑k=1pxik−xjkq1/q.Ifq=2, then the Euclidean distance dij2=∑k=1pxik−xjk21/2.The general steps of KMC [25, 26] are given in Figure 6:Step 1: select a number of sample points as cluster heads by experience or basic knowledge of dataStep 2: compute the Euclidean distance between each sample point and each cluster head, and allocate the sample point to the cluster of the nearest cluster headStep 3: classify one sample point at a time, and compute the mean of each cluster to update the cluster headStep 4: repeat Step 2 with the new cluster head, and classify the sample points againStep 5: repeat Steps 2–4 cyclically and adjust the classes until no cluster head changes to obtain the final classesFigure 6
Flow of KMC.


## 4.2. Solution of the Model

Table5 sums up the coordinates of the central cities of the prefectures in Hubei.Table 5
Coordinates of central cities of the prefectures in Hubei.
PrefecturexyWuhan21.832.04Huanggang24.871.75Xiaogan20.003.94Xiangyang11.549.84Suizhou17.498.20Yichang8.312.28Jingzhou12.510.29Jingmen11.914.02Ezhou24.261.16Huangshi25.610Xianning22.20−2.12Shiyan5.6612.96Xiantao18.100.45Enshi00Tianmen16.912.25Qianjiang15.450.79Shennongjia4.768.07According to the clustering results, Hubei was divided into three regions. The first region includes Yichang, Jingzhou, and Enshi; the second contains Wuhan, Huanggang, Xiaogan, Suizhou, Ezhou, Huangshi, Xianning, Xiantao, Tianmen, and Qianjiang; and the third covers Xiangyang, Jingmen, Shiyan, and Shennongjia.The coordinates of the head of each cluster were solved and taken as the location of the delivery site in the corresponding region, because the head is the closest point to any sample point in the cluster. The coordinates of the three heads were Yichang (8.31, 2.28), Wuhan (21.83, 2.04), and Xiangyang (11.54, 9.84).


## 5. Dynamic Material Delivery Optimization Model

By Model 1, it is possible to predict the number of patients in each prefecture in the near future and determine the newly added patients according to the distribution of the infected. From the perspective of emergency rescue, the number of patients and the growth of the infected must be suppressed by sufficient supply of vaccines. This is the only way to effectively control the epidemic [27, 28].Based on SPIR model, this article establishes a dynamic emergency delivery optimization model for multiple epidemic prefectures [29, 30]. The model has two objective functions: the minimum delivery time F1t and the minimum transport cost F2t.


### 5.1. Objective Functions

#### 5.1.1. Minimum Delivery Time

Considering the high death rate of COVID-19, vaccines and medical staff should be delivered as soon as possible, that is, within the threshold of the arrival time. Otherwise, the severity and scope of the epidemic will go out of control. The penalty cost incurred by the violation of the threshold of the arrival time can be described as an exponential function:(5)minF1t=emax0,tij−Tk.


#### 5.1.2. Minimum Transport Cost

Besides delivery time, the transport cost must also be minimized to deliver vaccines and medical staff to the epidemic prefecture. The transport cost covers delivery cost, staff salary, and storage cost:(6)minCt=∑i=1m∑j=1nCdXijMij+CsYij+C0,where C0 is the total fixed cost of setting up i storages; Cs is the subsidy of each medical staff member; Cd is the transport cost of a dose per unit distance; Yij is number of medical staff delivered from i to j; Mij is the quantity of vaccines delivered from i to j; and Xij is the distance between i and j.


#### 5.1.3. Composite Objective Function

The dynamic material delivery optimization model has two objective functions: the minimum delivery timeF1t and the minimum transport cost F2t. The two functions need to be normalized to unify their dimensionality:(7)minF1t¯=F1t−F1tminF1tmax−F1tmin,minF2t¯=F2t−F2tminF2tmax−F2tmin,where F1tmax and F1tmin are the maximum and minimum expectations of F1t, respectively, and F2tmax and F2tmin are the maximum and minimum expectations of F2t, respectively.Since two objectives have different priorities, weight coefficients were introduced to combine them into one composite objective function, turning the multiobjective optimization problem into a single-objective optimization problem:(8)minFt=αF1t¯−1−αF2t¯,α∈0,1.


### 5.2. Constraints

#### 5.2.1. Vaccine Demand

To control epidemic dissemination, the vaccine demand of each epidemic prefecture must be satisfied by the sum of vaccines and medical staff delivered to that prefecture:(9)∑i=1mMij=Dj,i=1,2,3,⋯m;j=1,2,3⋯n.


#### 5.2.2. Vaccine Supply

Practically, it is impossible to deliver an infinite number of materials to epidemic prefectures. The total vaccines transferred from storagei to each epidemic prefecture should not surpass the vaccine supply capability:(10)∑j=1nMij≤Si,i=1,2,3,⋯m;j=1,2,3⋯n.Every medical staff member involved in the rescue mission must be trained professionally, as medical staff have the highest infection rate of COVID-19 among all groups of people. According to the Health Commission of Hubei Province, there are a total of 54 million medical staff that could support the epidemic control in Hubei. This number should not be exceeded by the total medical staff delivered to epidemic prefectures. LetNij be a group of eight medical staff. In addition, there is a correspondence between medical staff groups and the infected: one group is needed to treat every 15 infected. The two numbers should grow in proportionality. If not, the rounding-plus-one method was adopted to confirm the number of groups:(11)∑i=1m∑j=1nYij≤540000,Yij=8Nij,Nij=βI,βI∈N+,βI+1,βI∉N+.


#### 5.2.3. Arrival Time

The period from the epidemic outbreak to the arrival of vaccines and medical staff at the epidemic prefecture covers the following stages: the period from vaccine production to delivery at each storage; the period for vaccines and medical staff to move from the storage to the epidemic prefecture; and the response period of the emergency system to the epidemic:(12)tij=XijV0+∑i=1m∑j=1nMijVp+t0,where t0 is the response period of the system; Vp is the vaccine production speed; and V0 is the delivery speed from storage to epidemic prefecture.


### 5.3. Solution of the Model

Referring to epidemic situations in previous years [31], the values of key parameters can be obtained as shown in Table 6.Table 6
Value of each key parameter.
ParametersC0CdCsTkValue32004.7251002.5Withα=0.8, the dynamic material delivery optimization model can be expressed as(13)minF1t=0.8F1t¯+0.2F2t¯,s.t.∑i=1mMij=Dj,i=1,2,3,⋯m;j=1,2,3⋯n,∑j=1mMij≤Si,i=1,2,3,⋯m;j=1,2,3⋯n,∑i=1m∑j=1nYij≤4116,Mij≥0,tij=XijV0+∑i=1m∑j=1nMijVp+t0.The distance between prefectures and delivery sites is shown in Table7.Table 7
Distance between prefectures and delivery sites.
PrefectureSiteYichangWuhanXiangyangWuhan288.910257.99Huanggang344.5757.05311.79Xiaogan253.2351.59206.9Xiangyang168.31257.990Suizhou228.69152.67119.86Yichang0288.91168.31Jingzhou99.81200.52187.95Jingmen95.11206.03110.11Ezhou346.159.81315.77Huangshi363.4982.89340.53Xianning306.6483.36317.67Shiyan221.79403.47146.45Xiantao210.985.31221.1Enshi178.91463.85318.53Tianmen180.64108.32178.53Qianjiang156.23137.41190.85Shennongjia130.28368142.83Through LINGO modeling, the optimal vaccine distribution scheme was obtained (Table8). The vaccine quantities delivered from the three sites were S1 = 2560, S2 = 27160, and S3 = 4520.Table 8
The drug distribution scheme.
Vaccine quantityS1-Jingzhou1840S1-Enshi720S2-Huanggang11040S2-Xiaogan4840S2-Suizhou3200S2-Ezhou2040S2-Huangshi1720S2-Xianning1600S2-Xiantao1720S2-Tianmen680S2-Qianjiang320S3-Jingmen3120S3-Shiyan1400S3-Shennongjia0The vaccine production speed wasVp=∑i=1m∑j=1mMij/t0=34240/28≈1223 doses/day.Then, the number of medical staff delivered from the three sites to each epidemic prefecture was calculated (Table9). The coordinates of the three delivery sites were S1 (8.31, 2.28), S2 (21.83, 2.04), and S3 (11.54, 9.84).Table 9
Medical staff delivery scheme.
Number of medical staffS1-Jingzhou30S1-Enshi12S2-Huanggang180S2-Xiaogan79S2-Suizhou52S2-Ezhou33S2-Huangshi28S2-Xianning26S2-Xiantao28S2-Tianmen11S2-Qianjiang5S3-Jingmen51S3-Shiyan23S3-Shennongjia0


## 5.1. Objective Functions

### 5.1.1. Minimum Delivery Time

Considering the high death rate of COVID-19, vaccines and medical staff should be delivered as soon as possible, that is, within the threshold of the arrival time. Otherwise, the severity and scope of the epidemic will go out of control. The penalty cost incurred by the violation of the threshold of the arrival time can be described as an exponential function:(5)minF1t=emax0,tij−Tk.


### 5.1.2. Minimum Transport Cost

Besides delivery time, the transport cost must also be minimized to deliver vaccines and medical staff to the epidemic prefecture. The transport cost covers delivery cost, staff salary, and storage cost:(6)minCt=∑i=1m∑j=1nCdXijMij+CsYij+C0,where C0 is the total fixed cost of setting up i storages; Cs is the subsidy of each medical staff member; Cd is the transport cost of a dose per unit distance; Yij is number of medical staff delivered from i to j; Mij is the quantity of vaccines delivered from i to j; and Xij is the distance between i and j.


### 5.1.3. Composite Objective Function

The dynamic material delivery optimization model has two objective functions: the minimum delivery timeF1t and the minimum transport cost F2t. The two functions need to be normalized to unify their dimensionality:(7)minF1t¯=F1t−F1tminF1tmax−F1tmin,minF2t¯=F2t−F2tminF2tmax−F2tmin,where F1tmax and F1tmin are the maximum and minimum expectations of F1t, respectively, and F2tmax and F2tmin are the maximum and minimum expectations of F2t, respectively.Since two objectives have different priorities, weight coefficients were introduced to combine them into one composite objective function, turning the multiobjective optimization problem into a single-objective optimization problem:(8)minFt=αF1t¯−1−αF2t¯,α∈0,1.


## 5.1.1. Minimum Delivery Time

Considering the high death rate of COVID-19, vaccines and medical staff should be delivered as soon as possible, that is, within the threshold of the arrival time. Otherwise, the severity and scope of the epidemic will go out of control. The penalty cost incurred by the violation of the threshold of the arrival time can be described as an exponential function:(5)minF1t=emax0,tij−Tk.


## 5.1.2. Minimum Transport Cost

Besides delivery time, the transport cost must also be minimized to deliver vaccines and medical staff to the epidemic prefecture. The transport cost covers delivery cost, staff salary, and storage cost:(6)minCt=∑i=1m∑j=1nCdXijMij+CsYij+C0,where C0 is the total fixed cost of setting up i storages; Cs is the subsidy of each medical staff member; Cd is the transport cost of a dose per unit distance; Yij is number of medical staff delivered from i to j; Mij is the quantity of vaccines delivered from i to j; and Xij is the distance between i and j.


## 5.1.3. Composite Objective Function

The dynamic material delivery optimization model has two objective functions: the minimum delivery timeF1t and the minimum transport cost F2t. The two functions need to be normalized to unify their dimensionality:(7)minF1t¯=F1t−F1tminF1tmax−F1tmin,minF2t¯=F2t−F2tminF2tmax−F2tmin,where F1tmax and F1tmin are the maximum and minimum expectations of F1t, respectively, and F2tmax and F2tmin are the maximum and minimum expectations of F2t, respectively.Since two objectives have different priorities, weight coefficients were introduced to combine them into one composite objective function, turning the multiobjective optimization problem into a single-objective optimization problem:(8)minFt=αF1t¯−1−αF2t¯,α∈0,1.


## 5.2. Constraints

### 5.2.1. Vaccine Demand

To control epidemic dissemination, the vaccine demand of each epidemic prefecture must be satisfied by the sum of vaccines and medical staff delivered to that prefecture:(9)∑i=1mMij=Dj,i=1,2,3,⋯m;j=1,2,3⋯n.


### 5.2.2. Vaccine Supply

Practically, it is impossible to deliver an infinite number of materials to epidemic prefectures. The total vaccines transferred from storagei to each epidemic prefecture should not surpass the vaccine supply capability:(10)∑j=1nMij≤Si,i=1,2,3,⋯m;j=1,2,3⋯n.Every medical staff member involved in the rescue mission must be trained professionally, as medical staff have the highest infection rate of COVID-19 among all groups of people. According to the Health Commission of Hubei Province, there are a total of 54 million medical staff that could support the epidemic control in Hubei. This number should not be exceeded by the total medical staff delivered to epidemic prefectures. LetNij be a group of eight medical staff. In addition, there is a correspondence between medical staff groups and the infected: one group is needed to treat every 15 infected. The two numbers should grow in proportionality. If not, the rounding-plus-one method was adopted to confirm the number of groups:(11)∑i=1m∑j=1nYij≤540000,Yij=8Nij,Nij=βI,βI∈N+,βI+1,βI∉N+.


### 5.2.3. Arrival Time

The period from the epidemic outbreak to the arrival of vaccines and medical staff at the epidemic prefecture covers the following stages: the period from vaccine production to delivery at each storage; the period for vaccines and medical staff to move from the storage to the epidemic prefecture; and the response period of the emergency system to the epidemic:(12)tij=XijV0+∑i=1m∑j=1nMijVp+t0,where t0 is the response period of the system; Vp is the vaccine production speed; and V0 is the delivery speed from storage to epidemic prefecture.


## 5.2.1. Vaccine Demand

To control epidemic dissemination, the vaccine demand of each epidemic prefecture must be satisfied by the sum of vaccines and medical staff delivered to that prefecture:(9)∑i=1mMij=Dj,i=1,2,3,⋯m;j=1,2,3⋯n.


## 5.2.2. Vaccine Supply

Practically, it is impossible to deliver an infinite number of materials to epidemic prefectures. The total vaccines transferred from storagei to each epidemic prefecture should not surpass the vaccine supply capability:(10)∑j=1nMij≤Si,i=1,2,3,⋯m;j=1,2,3⋯n.Every medical staff member involved in the rescue mission must be trained professionally, as medical staff have the highest infection rate of COVID-19 among all groups of people. According to the Health Commission of Hubei Province, there are a total of 54 million medical staff that could support the epidemic control in Hubei. This number should not be exceeded by the total medical staff delivered to epidemic prefectures. LetNij be a group of eight medical staff. In addition, there is a correspondence between medical staff groups and the infected: one group is needed to treat every 15 infected. The two numbers should grow in proportionality. If not, the rounding-plus-one method was adopted to confirm the number of groups:(11)∑i=1m∑j=1nYij≤540000,Yij=8Nij,Nij=βI,βI∈N+,βI+1,βI∉N+.


## 5.2.3. Arrival Time

The period from the epidemic outbreak to the arrival of vaccines and medical staff at the epidemic prefecture covers the following stages: the period from vaccine production to delivery at each storage; the period for vaccines and medical staff to move from the storage to the epidemic prefecture; and the response period of the emergency system to the epidemic:(12)tij=XijV0+∑i=1m∑j=1nMijVp+t0,where t0 is the response period of the system; Vp is the vaccine production speed; and V0 is the delivery speed from storage to epidemic prefecture.


## 5.3. Solution of the Model

Referring to epidemic situations in previous years [31], the values of key parameters can be obtained as shown in Table 6.Table 6
Value of each key parameter.
ParametersC0CdCsTkValue32004.7251002.5Withα=0.8, the dynamic material delivery optimization model can be expressed as(13)minF1t=0.8F1t¯+0.2F2t¯,s.t.∑i=1mMij=Dj,i=1,2,3,⋯m;j=1,2,3⋯n,∑j=1mMij≤Si,i=1,2,3,⋯m;j=1,2,3⋯n,∑i=1m∑j=1nYij≤4116,Mij≥0,tij=XijV0+∑i=1m∑j=1nMijVp+t0.The distance between prefectures and delivery sites is shown in Table7.Table 7
Distance between prefectures and delivery sites.
PrefectureSiteYichangWuhanXiangyangWuhan288.910257.99Huanggang344.5757.05311.79Xiaogan253.2351.59206.9Xiangyang168.31257.990Suizhou228.69152.67119.86Yichang0288.91168.31Jingzhou99.81200.52187.95Jingmen95.11206.03110.11Ezhou346.159.81315.77Huangshi363.4982.89340.53Xianning306.6483.36317.67Shiyan221.79403.47146.45Xiantao210.985.31221.1Enshi178.91463.85318.53Tianmen180.64108.32178.53Qianjiang156.23137.41190.85Shennongjia130.28368142.83Through LINGO modeling, the optimal vaccine distribution scheme was obtained (Table8). The vaccine quantities delivered from the three sites were S1 = 2560, S2 = 27160, and S3 = 4520.Table 8
The drug distribution scheme.
Vaccine quantityS1-Jingzhou1840S1-Enshi720S2-Huanggang11040S2-Xiaogan4840S2-Suizhou3200S2-Ezhou2040S2-Huangshi1720S2-Xianning1600S2-Xiantao1720S2-Tianmen680S2-Qianjiang320S3-Jingmen3120S3-Shiyan1400S3-Shennongjia0The vaccine production speed wasVp=∑i=1m∑j=1mMij/t0=34240/28≈1223 doses/day.Then, the number of medical staff delivered from the three sites to each epidemic prefecture was calculated (Table9). The coordinates of the three delivery sites were S1 (8.31, 2.28), S2 (21.83, 2.04), and S3 (11.54, 9.84).Table 9
Medical staff delivery scheme.
Number of medical staffS1-Jingzhou30S1-Enshi12S2-Huanggang180S2-Xiaogan79S2-Suizhou52S2-Ezhou33S2-Huangshi28S2-Xianning26S2-Xiantao28S2-Tianmen11S2-Qianjiang5S3-Jingmen51S3-Shiyan23S3-Shennongjia0


## 6. Conclusions

Based on the characteristics of the spread of COVID-19 and treatment procedures, three models of organic linkage were established to address the epidemic prevention and control issues systematically. The results reflect the stability and reliability of our models.Model 1 considers the driving effect of COVID-19 diffusion on the distribution of materials and the inverse control of material flow on disease spread. The modeling results are reasonable and practical, because various factors were coupled with the SPIR infectious diseases model.Model 2 relies on KMC to cluster the epidemic prefectures by the Euclidean distance. The modeling results lay the basis for optimizing the material delivery of Model 3 and help to choose the best delivery sites.Model 3 takes account of the competition between epidemic prefectures for medical resources, and the distance between epidemic prefecture and delivery site, as well as the severity of the epidemic. In this way, the authors derived a practical distribution scheme for emergency materials.



---

*Source: 6329601-2022-01-05.xml*