Abstract
This paper studies an SEIR-type epidemic model with time delay and vaccination control. The vaccination control is applied when the basic reproduction number \(R_0>1\). The vaccination strategy is expressed as a state delayed feedback which is related to the current and previous state of the epidemic model, and makes the model become a linear system in new coordinates. For the presence and absence of vaccination control, we investigate the nonnegativity and boundedness of the model, respectively. We obtain some sufficient conditions for the eigenvalues of the linear system such that the nonnegativity of the epidemic model can be guaranteed when the vaccination strategy is applied. In addition, we study the stability of disease-free equilibrium when \(R_0<1\) and the persistent of disease when \(R_0>1\). Finally, we use the obtained theoretical results to simulate the vaccination strategy to control the spread of COVID-19.
Similar content being viewed by others
1 Introduction
The outbreak of COVID-19 has caused a large number of deaths and economic losses to the world [1,2,3,4]. At present, many countries have not effectively controlled the epidemic, and the worst is that there are many asymptomatic infections [5,6,7,8]. Asymptomatic infections can transmit the virus to others, but they do not have any symptoms. In fact, the presence of asymptomatic infections is not unique to COVID-19, and many diseases have this property, such as Epstein–Barr virus, HIV, typhoid fever, etc. However, unlike other diseases having asymptomatic infections, the high infectiousness and a large number of asymptomatic infections make it difficult to control the spread of COVID-19 [8]. In order to effectively control the epidemic, various effective vaccines are the best choice. For COVID-19, many vaccines have entered phase III clinical trials [9,10,11]. Appropriate mathematical models and vaccination strategies are currently needed.
This paper introduces an SEIR type epidemic model with time delay, and a vaccination control is provided to make the removed-by-immunity population asymptotically track the whole population. The model assumes that there are many asymptomatic infections and exists an incubation period during which the strain of disease develops in the asymptomatic infections and it is only after that time that the asymptomatic infections become contagious. These two assumptions are suitable for COVID-19. Time delay is a common phenomenon in the spread of various diseases [12,13,14,15]. Since Cooke studied an SIR model with time delay which is caused by a vector [12], many epidemic models began to consider the factor of time delay in order to better simulate the spread of different diseases, such as SIS [16,17,18], SEIR [19,20,21] and SEIRS models [22,23,24]. The above epidemic models cannot describe the case that there are many asymptomatic infections. The two main factors causing time delay are (1) the latency of infected vector becoming infectious [12] and (2) the latency of the infected host becoming infectious. In addition, the time delay may occur in other scenarios, such as vaccination period and immune period [24].
Vaccination is the most effective way to completely control the spread of various diseases, such as smallpox [25] and measles [26]. When there exist various effective vaccines, appropriate vaccination strategy is needed to improve the efficacy of vaccination [27,28,29]. Since the epidemic models are nonlinear, many nonlinear control methods are used to design appropriate vaccination strategy [30,31,32,33,34,35,36]. For example, in [30], the authors considered an SEIR epidemic model and proposed a vaccination strategy based on feedback linearization techniques. A continuous linear vaccination control was proposed for SEIR model [31]. In [32], the impulsive vaccination was designed for cholera epidemics. The impulsive vaccination strategy was also studied for other cases, such as SIR epidemic model [33] and SEIR epidemic model with time delay [34]. The sliding mode control has been used to design vaccination strategy for SEIR epidemic model [36]. Although there are many vaccination strategies for various epidemic models, they cannot be directly used for proposed epidemic model with time delay. This paper will use feedback linearization technique to design a switching vaccination strategy for the proposed model.
The main contributions of this paper include the following: (1) We propose an SEIR type epidemic model with time delay based on the properties of COVID-19 (high infectiousness and a large number of asymptomatic infections). When there do not exist vaccination control, we analyze the stability of the disease-free equilibrium and persistence of the disease. (2) Based on feedback linearization technique, we design a switching vaccination strategy for the proposed model, and obtain sufficient conditions for the nonnegativity of the model when the vaccination strategy is applied. (3) We simulate the vaccination strategy to control the spread of COVID-19 and show the effectiveness of the obtained results.
The paper is organized as follows: Section 2 presents the problem statement and some basic properties of the model. Section 3 shows the stability of disease-free equilibrium when \(R_0<1\), and the persistence of disease when \(R_0>1\). Section 4 provides a vaccination strategy for the proposed model. Section 5 simulates the vaccination control of COVID-19. Section 6 collects our conclusions.
2 Problem statement and basic properties
Motivated by Refs. [15, 31], suppose that the total population N(t) can be divided into four compartments, that is, S(t), E(t), I(t), and R(t), where S(t), E(t), I(t), R(t) denote the susceptible, exposed, infectious, recovered individuals at time t, respectively. Consider the SEIR-type epidemic model with time delay and vaccination control which is presented as following functional differential equation:
where all parameters are nonnegative (see Table 1 for details), u(t) is a vaccination control needed to design. u(t) can be seen as a vaccination rate [13, 28]. The flow diagram for model (1) is presented in Fig. 1. This epidemic model involves the following assumptions:
-
(1)
There exists a vector between susceptible and exposed individuals. The vector will become infectious after a time \(\tau >0\). The susceptible population in S enter into E when they come into contact with infected vectors.
-
(2)
The population in E move to I after the average latent time \(\beta ^{-1}\).
-
(3)
Part of the recovered people become susceptible again to the disease.
-
(4)
The death rate in I is large than the natural death rate.
Note that the vector may be an animal (mosquito) [12] or a human (asymptomatic infection) [5]. Based on Assumption 1), in model (1), the total number of infections per unit of time at time t is given by \(\frac{\alpha S(t)I(t-\tau )}{N_0}\), which is called “true-mass action model” [31]. The model (1) can be used to describe the case that there exist many asymptomatic infections, such as COVID-19. In Sect. 5, we will use our model to simulate the second outbreak of COVID-19 in Italy and show the effectiveness of the proposed vaccination control.
Let \({\mathbb {R}}^n\) denote the set of all real n-dimensional column vectors, and \({\mathbb {R}}^{n}_{+}:=\{x\in {\mathbb {R}}^n: x_i\ge 0,\, i=1,\ldots ,n\}\). Let \({\mathscr {C}}:={\mathscr {C}}([-\tau ,0],{\mathbb {R}}^4)\) be the Banach space of continuous functions mapping the interval \([-\tau ,0]\) into \({\mathbb {R}}^4\) with norm \(||\phi ||=\sum _{i=1}^4\left( \max _{\theta \in [-\tau ,0]}|\phi _i(\theta )|\right) \), where \(\phi =[\phi _1,\ldots ,\phi _4]\in {\mathscr {C}}\). Let \({\mathscr {C}}^+:={\mathscr {C}}([-\tau ,0],{\mathbb {R}}_+^4)\) which is a cone of \({\mathscr {C}}\). We consider the solutions of system (1) with initial conditions
where \(\phi (s)\in {\mathscr {C}}^+\).
One can know from (1) that
If the death rate \(\mu =\mu ^I\), then \({\dot{N}}=\mu (N_0-N)\) and the total population \(N(t)=N(0)=N_0, \forall t\ge 0\). Generally speaking, \(\mu ^I\) is slightly larger than the death rate \(\mu \). Hence, \(N\le \frac{\mu -e^{-\psi t}(\mu -\psi )}{\psi }N_0\), where \(\psi =\min \left\{ \mu ,\mu ^I\right\} \). Since \(\mu -\psi \approx 0\), the total population N can be regarded as constant in a short time.
In this paper, we want to design a vaccination strategy to make the infective population approach zero, and the removed-by-immunity population asymptotically track the whole population.
When \(u(t)\equiv 0\), that is, there does not exist vaccination control, the following theorem shows that the system (1) has positively invariant set
Theorem 1
Consider the system (1) with initial condition (2). Then, the system (1) has positively invariant set \(\varGamma \).
Proof
Let \(g(\phi )=[g_1(\phi ),g_2(\phi ),g_3(\phi ),g_4(\phi )]^{\top }\), where \(\forall \phi \in {\mathscr {C}}^+\), and
Then, for each compact set \(\varOmega \in {\mathbb {R}}\times {\mathscr {C}}^+\), \(g(\phi )\) is continuous and Lipschitz in \(\varOmega \). Hence, there exists a unique solution of system (1) through \((0,\phi )\) [37, Theorem 2.2.3]. Since \(\phi =[0,0,0,0]^{\top }\) is a solution of system (1), all solutions of system (1) are nonnegative at maximal interval. One can see from the above discussion that the total population \(N(t)\le N_0, \forall t\ge 0\). Hence, all solutions of system (1) are bounded, and exist on the interval \([0,\infty )\). Thus, the set \(\varGamma \) is positively invariant for system (1). \(\square \)
3 Stability analysis
This section will introduce the basic reproduction number of system (1), present the stability analysis of the disease-free equilibrium, and show that the disease is persistent when the basic reproduction number is greater than one.
When \(u(t)\equiv 0\), the disease-free equilibrium of model (1) is \(P_0=[N_0,0,0,0]^{\top }\). Based on Ref. [14], let \({\mathscr {F}}\) and \({\mathscr {Y}}\) denote the rate of secondary infections increase and the rate of disease progression, death, recovery decrease, respectively. Then,
The next generation matrix of system (1) can be presented as \(FY^{-1}\), where \(F=(\frac{\partial {\mathscr {F}}}{\partial x}(P_0)), Y=(\frac{\partial {\mathscr {Y}}}{\partial x}(P_0)),x=[E,I]^{\top }\).
Hence,
and the next-generation matrix is \(K=FY^{-1}\),
Then, the basic reproduction number equals the spectral radius of matrix K, that is,
Theorem 2
Consider the system (1). If \(R_0<1\), then the disease-free equilibrium \(P_0\) is globally asymptotically stable. If \(R_0>1\), then the disease-free equilibrium \(P_0\) is unstable.
Proof
We will use the method of fluctuations [38, Lemma 4.2] to show that
For simplicity, let \({\bar{S}}=\lim _{t\rightarrow \infty }\sup S(t)\), \({\bar{E}}=\lim _{t\rightarrow \infty }\sup E(t)\), \({\bar{I}}=\lim _{t\rightarrow \infty }\sup I(t)\), \({\bar{R}}=\lim _{t\rightarrow \infty }\sup R(t)\). By Lemma 4.2 [38], there exist four sequences \(\pi _n^i\rightarrow \infty , n\rightarrow \infty , i=1,2,3,4\) such that
Then, one has the following inequalities from (1)
The third inequality of (4) implies \({\bar{E}}\le \frac{(\gamma +\mu ^{I}){\bar{I}}}{\beta }\). Since \({\bar{S}}\le N_0\), the second inequality of (4) implies \(0\le \alpha {\bar{I}}-(\beta +\mu ){\bar{E}}\). Hence, \(0\le (R_0-1){\bar{I}}\). Since \(R_0<1\), \({\bar{I}}={\bar{E}}=0\). The fourth inequality of (4) implies \({\bar{R}}=0\). The first inequality of (4) implies \({\bar{S}}=N_0\). Hence, the disease-free equilibrium \(P_0\) is globally attractive.
One can know from Theorem 2.1 in [14] that the disease-free equilibrium \(P_0\) is locally stable when \(R_0<1\), and the disease-free equilibrium \(P_0\) is unstable when \(R_0>1\). Hence, if \(R_0<1\), then the disease-free equilibrium \(P_0\) is globally asymptotically stable. \(\square \)
When \(R_0>1\), the disease will persist in the population. The following theorem will provide the detail proof.
Theorem 3
Consider the system (1). If \(R_0>1\), then the disease is uniformly persistent, that is, there exists a positive number b such that \(\lim _{t\rightarrow \infty }\inf I(t)\ge b\).
Proof
For the system (1), let \(x(t,\phi )=[S(t),E(t),I(t),\)\(R(t)]^{\top }\) be the solution of system (1) with initial condition \(\phi \in {\mathscr {C}}^+\). Because \(S(t)\rightarrow N_0\) as \(x(t,\phi )\rightarrow P_0\), for sufficiently small \(\nu >0\) there exists a neighborhood of \(P_0\) such that \(N_0-\nu<S(t)<N_0+\nu \).
First, we will show that there exists \(\varepsilon (\nu )>0\) such that \(\lim _{t\rightarrow \infty }\sup ||x(t,\phi )-P_0||\ge \varepsilon (\nu )\) for any \(\phi \in {\mathscr {C}}^+\). We will prove it by contradiction. If the above condition does not hold, then there exists \(\phi \in {\mathscr {C}}^+\) such that \(||x(t,\phi )-P_0||<\varepsilon (\nu )\). Hence, there exists \(T>0\) such that \(|S(t)-N_0|<\varepsilon (\nu )\), \(|E(t)|<\varepsilon (\nu )\), \(|I(t)|<\varepsilon (\nu )\), and \(|R(t)|<\varepsilon (\nu )\) for \(t>T\). One can obtain from the second and third equation of (1) that
By Lemma 4.2 [38], there exist a sequence \(\pi _n\rightarrow \infty , n\rightarrow \infty \) such that \(\lim _{n\rightarrow \infty }E(\pi _n)=0,\;{\dot{E}}(\pi _n)=0\). Hence, the following inequality holds
Consider the following equation
The characteristic equation of (7) is
When \(\nu =0\) and \(R_0>1\), \(\varTheta (0)=(\beta +\mu )(\gamma +\mu ^{I})(1-R_0)<0\). It is easy to see that \(\varTheta (\lambda )\) is a continuous function for \(\lambda \ge 0\) and increases to \(\infty \) as \(\lambda \rightarrow \infty \). Hence, for Eq. (8) there exists a positive root. Suppose that \({\bar{\lambda }}(\nu )\) is the principle eigenvalue. Thus, \({\bar{\lambda }}(0)>0\). Since \({\bar{\lambda }}(\nu )\) is a continuous function for \(\nu \), one can get \({\bar{\lambda }}(\nu )>0\) for sufficiently small \(\nu >0\). Hence, \(ce^{{\bar{\lambda }}(\nu )t}\) (\(c>0\)) is a solution of Eq. (7). By the comparison theorem, one can obtain that there exists \(a>0\) such that \(I(\pi _n)\ge ace^{{\bar{\lambda }}(\nu )\pi _n}\). Hence, \(I(\pi _n)\rightarrow \infty \) as \(n\rightarrow \infty \), which is a contraction for the fact that I(t) is bounded.
Let \(\omega (\phi )\) be the omega limit set of \(x(t,\phi )\). Define \(M=\{\phi \in {\mathscr {C}}^+:x(t,\phi )\in \partial {\mathscr {C}}^+,t\ge 0\}\). One can obtain from the third equation of (1) that
Hence, if there exists \(T\ge 0\) such that \(I(T)>0\), then \(I(t)>0\) for all \(t\ge T\). Thus, \(I(t)=0\) when \(\phi \in M\). For any \(\phi \in M\) \(\omega (\phi )=P_0\), and there does not exist a cycle in M. With the similar proof of Theorem 3.1 [14], there exists \(b>0\) such that \(\lim _{t\rightarrow \infty }\inf I(t)\ge b\). \(\square \)
4 Vaccination control strategy
This section will design a vaccination control u(t) which makes the removed-by-immunity population asymptotically track the whole population. u(t) will present as a state delayed feedback.
Since \(N=S+E+I+R\), we can obtain that
Hence, \(N=e^{-\mu t}N_0+\int _{0}^{t}e^{-\mu (t-s)}[\mu N_0+(\mu -\mu ^I)I]ds\). Because \(R=N-S-E-I\) and N can be presented as a function of I, we just need to consider the first three equations about S, E, I in model (1).
The system (1) can be rewritten as the following nonlinear delay system:
where \(x(t)=[S(t),\;E(t),\;I(t)]^\top \in {\mathbb {R}}_{+}^3\),
In the sequel, let \(x_{i\tau }(t)=x(t-i\tau )\), \(u_{i\tau }(t)=u(t-i\tau )\), \(i=0,1,\ldots \). \(x_{i\tau }\) is defined for \(t\ge (i-1)\tau \), and \(u_{i\tau }\) is defined for \(t\ge i\tau \).
Definition 1
Assume that for system (9), there exists an integer r such that, for every X in \(\varOmega _r\subseteq {\mathbb {R}}^{3(r+1)}\), the following conditions are satisfied
where
then the delay relative degree of system (9) is equal to r in \(\varOmega _r\).
By a simple calculation, \(L_{G}L_{F}^{k}H(X)=0, k=0,1\) and \(L_{G}L_{F}^{2}H(X)[1,0,0]^{\top }=-\mu \beta \alpha I_{\tau }\). Hence, the delay relative degree of system (9) is 3 in \(\varOmega _3:=\{X:x=[S,\;E,\;I]^\top ,\;I_{\tau }\ne 0\}\). Then, we can obtain that
where
Let \(z(t)=[y(t),\;{\dot{y}}(t),\;\ddot{y}(t)]^\top \). Then, we have
If we choose
then the equation (12) becomes
where the parameters \(k_1, k_2\) and \(k_3\) are needed to design.
The above analysis process can be presented in the following proposition:
Proposition 1
Consider the system (1) with time delay. If the state delayed feedback control law u(t) defined as in (13), then one can obtain the closed-loop system (14).
Remark 1
Note that the control law u(t) is well defined all \(x\in \varOmega _3:=\{X:x=[S,\;E,\;I]^\top ,\;I_{\tau }\ne 0\}\). However, the disease may be considered eradicated from the population once the infectious population is smaller than one individual. Hence, the vaccination strategy can be switched off when \(0<I(t)\le \varepsilon <1\). In this sense, the efficient control law is given by
where \(t_{f}\) is defined as follows:
Remark 2
Note that although the vector X is defined for \(t\ge 2\tau \), the control u(t) is a function of x and \(x(t-\tau )\). Hence, the control \(u_{e}(t)\) can be applied at \(\tau \).
Remark 3
For linear system (14), its characteristic polynomial \(P(s)=s^3-k_3s^2-k_2s-k_1=0\). Making \(P(s)=(s+r_1)(s+r_2)(s+r_3)\). In order to make the linear system (14) asymptotically stable, the parameters \(k_1, k_2, k_3\) should be chosen to make \(r_i>0 ,i=1,2,3\). It can be get that \(k_1=-r_1r_2r_3<0, k_2=-(r_{1}r_{2}+r_{1}r_{3}+r_{2}r_{3})<0, k_3=-(r_{1}+r_{2}+r_3)<0\).
When the control \(u_{e}(t)\) is applied at \(\tau \), we must guarantee that S, E, I, R of system (1) are nonnegative. Hence, the parameters \(k_1, k_2, k_3\) should be chosen to make \(r_i>0 ,i=1,2,3\) satisfy the following theorem.
Theorem 4
Consider the SEIR type epidemic model (1) with an initial condition (2), and \(I(0)\le I(\tau )\). Suppose that the positive real numbers \(r_i, i=1,2,3\) satisfy the following conditions:
-
(1)
\(r_1<\min \{\mu +\beta ,\mu ^I+\gamma \}\), \(r_2=\mu ^I+\gamma \), and \(r_3>\max \{\mu +\beta ,\mu ^I+\gamma \}\).
-
(2)
\(r_1\) and \(r_3\) satisfy the inequalities:
$$\begin{aligned} {\left\{ \begin{array}{ll} (r_{3}-r_{1})(r_{3}-\mu ^I-\gamma )>\beta \alpha ,\\ r_1+r_3-\alpha -\mu -(\beta +\mu +\mu ^I+\gamma )\ge 0,\\ (\beta +\mu )^2+r_1r_3-(r_1+r_3)(\beta +\mu )\ge 0. \end{array}\right. } \end{aligned}$$(17)Then, the application of control law (15) guarantees the nonnegativity of the SEIR-type epidemic model (1), and \(u(t)\ge 1, \forall t\in [\tau ,t_f]\).
Proof
Because \(z(t)=[y(t),\;{\dot{y}}(t),\;\ddot{y}(t)]^\top \) and \(y(t)=I(t)\), \(I(t)=c_{1}e^{-r_{1}(t-\tau )}+c_{2}e^{-r_{2}(t-\tau )}+c_{3}e^{-r_{3}(t-\tau )}, \forall t\ge \tau \). The following equations can be obtained
The above equations can be expressed in matrix form \(R\cdot C=M\), where
The matrix R is the Vandermonde matrix and its determinant is equal to \((r_{2}-r_{1})(r_{3}-r_1)(r_{3}-r_{2})\). Hence, R is an invertible matrix when \(r_{1}\ne r_{2}\ne r_{3}\). By direct calculation,
If the condition (1) holds, then \(c_{1}=\frac{\beta (r_{3}-\mu -\beta )E(0)+\beta \frac{\alpha S(0)I(-\tau )}{N}}{(r_{3}-r_{1})(\mu +\gamma -r_{1})}\) and \(r_{3}>r_{2}>r_{1}\). Hence, \(c_{1}>0\). In the sequel, we show that \(I(t)\ge 0\).
Since the sign of \(c_{2}\) and \(c_{3}\) is uncertain, there are four cases: \((i)\ c_{2}\ge 0\) and \(c_{3}\ge 0\), \((ii)\ c_{2}\ge 0\) and \(c_{3}<0\), \((iii)\ c_{2}<0\) and \(c_{3}\ge 0\), and \((iv)\ c_{2}<0\) and \(c_{3}<0\).
For the case (i) and (ii), since \(c_{2}\ge 0\), one can obtain
For the case (iii), that is \(c_{2}<0\) and \(c_{3}\ge 0\), the following equation holds:
Since \(I(0)\le I(\tau )\), \(S(\tau )\le N_0\), one has
When the condition (2) holds, then \(I(\tau )-c_{3}\ge 0\). Hence, for the case (iii), \(I(t)\ge 0, \forall t\ge 0\).
For the case (iv), that is \(c_{2}<0\) and \(c_{3}<0\), we have
Hence, for the cases (i)–(iv), \(I(t)\ge 0, \forall t\ge 0\).
We will show that \(E(t)\ge 0\) and \(S(t)\ge 0\). We can obtain from the third equation of (1) that
Taking \(r_{2}=\mu ^I+\gamma \), then
Hence, \(c_{3}(\mu ^I+\gamma -r_{3})=\beta E(\tau )-c_{1}(\mu ^I+\gamma -r_{1})\). Then, one can obtain from (20) that
Similarly, one can obtain from the second equation of (1) that
Combining (20) and (21), one has
From the above equation, it follows that
Hence,
The nonnegativity of R(t) is proven by using continuity arguments. If R(t) has a negative value, then there is a point \(t_0\) such that \(R(t_0)=0\) when \(R(0)\ge 0\). Then, it follows from the fourth equation of (1) that
Because \(k_1=-r_1r_2r_3<0, k_2=-(r_{1}r_{2}+r_{1}r_{3}+r_{2}r_{3})<0, k_3=-(r_{1}+r_{2}+r_3)<0\) and \(r_2=\mu ^I+\gamma \), then one has
If the condition (17) holds, then \({\dot{R}}(t_0)\ge 0\) and there does not exist a point \(t_1>t_0\) such that \(R(t_1)<0\). Hence, \(R(t)\ge 0, \forall t\ge \tau \).
It follows from (13) that
When \(k_1=-r_1r_2r_3<0, k_2=-(r_{1}r_{2}+r_{1}r_{3}+r_{2}r_{3})<0, k_3=-(r_{1}+r_{2}+r_3)<0\) and \(r_2=\mu ^I+\gamma \), one can get
If the condition (17) holds, then \(u(t)\ge 1\). \(\square \)
Theorem 4 presents some conditions for the nonnegativity of the SEIR-type epidemic model (1). The following theorem shows that the control law (15) can make the symptomatic infections disappear after a finite time \(t_f\).
Theorem 5
Consider the SEIR-type epidemic model (1) with an initial condition (2), and \(I(0)\le I(\tau )\). Suppose that the conditions (1) and (2) of Theorem 4 hold. Then, the application of control law (15) guarantees that disease is asymptotically eradicated from the symptomatic infections after a finite time \(t_f\).
Proof
For the linear system (14), \(z(t)=[y(t),\;{\dot{y}}(t),\;\ddot{y}(t)]^\top \) and \(y(t)=I(t)\). Hence,
When the control law (15) is applied, \(I(t-\tau )>\varepsilon \). Thus,
If the conditions (1) and (2) of Theorem 4 hold, then \(I(t_f)=z_1(t_f)\le \varepsilon <1\). If we choose \(\varepsilon <\min \{\frac{\beta }{1+\mu ^I+\gamma },1\}\), then \(E(t_f)<1\) which is obtained from the second equation of (23). Hence, disease is asymptotically eradicated from the host population (eventually becoming symptomatic infections) after a finite time \(t_f\). \(\square \)
Remark 4
Since we suppose that there exist many asymptomatic infections, and we see the asymptomatic infections as vectors, the control law (15) just makes symptomatic infections and the host population (eventually becoming symptomatic infections) disappear after a finite time \(t_f\). Hence, when there are some effective vaccines, we do not need to consider a large number of asymptomatic infections. As long as the population is vaccinated according to strategy, the virus will not cause a large number of deaths.
5 Simulating the vaccination control of COVID-19
This section will provide a simulation of the second outbreak of COVID-19 in Italy, and shows that the symptomatic infections can be eradicated by using the obtained vaccination control. According to reports [39], the incubation period of COVID-19 was 2–14 days. We take 5.8 days as time units which is the mean incubation period of COVID-19 [40]. The rate of exposed class who become infective \(\beta \) is equal to \(1/5.2=0.1923\). Suppose that the time delay \(\tau =3\).As pointed out in [41], the natural death rate in Italy is equal to \(12.04\%\). The case fatality rate of COVID-19 \(\mu ^I=7.7\%\) [42]. There are lack of real data about the proportion of recovery population who become susceptible again, transmission rate and treatment rate. Hence, we assume that \(\delta =1\%\), \(\alpha =0.5\), and \(\gamma =0.151\).
We set August 1, 2020, as time zero. Around August 1, the second outbreak of COVID-19 occurred in Italy, and a large number of asymptomatic infections appeared. There are about 60.39 million permanent residents in Italy. The number of susceptible people is approximately equal to that of its permanent residents. Since the epidemic in Lombardy is very serious, we assume that the number of exposed people is equal to the permanent residents of Lombardy. The initial value of the infective and the recovered can obtain from published data [43]. Hence, \(S(\omega )=60.09\) million, \(E(\omega )=33000\), \(I(\omega )=12646\), \(R(\omega )=200229\), \(-\tau \le \omega \le 0\). Figure 2 depicts the actual number of active infections in Italy according to the data of coronavirus resource center [43].
Based on the above data, it is easy to obtain that \(R_{0}=2.0642>1\). As pointed out in [44], the obtained basic reproduction number \(R_{0}\) is reasonable for COVID-19 in Italy. Figure 3a depicts the time evolution of system (1) when there does not exist control. Figure 3b depicts the number of active infections with no control. Comparing with the actual number of active infections in Italy (see Fig. 2), the number of active infections in Italy is much higher than the actual data. The reason for this phenomenon is that we assume the medical conditions are very poor and no other methods of epidemic control have been implemented, such as wearing a mask, not going outside and keeping social distance.
If there is an effective vaccine, then the government can apply vaccination strategy (15). If we choose \(\varepsilon =0.9\), \(r_1=\beta +0.99\mu \), \(r_2=\mu ^I+\gamma \) and \(r_3=3.86(\mu ^I+\gamma )\), then the conditions of Theorem 4 hold, and the nonnegativity of system (1) is guaranteed. The vaccination strategy (15) is applied when \(t=3\) days. By direct calculation, \(t_{f}\approx 66\) days. Figure 4a depicts the time evolution of system (1) with vaccination control (15). Figure 4b depicts the number of active infections with control \(u_e(t)\). One can see that the active infections are eradicated from the population. Moreover, the removed-by-immunity population asymptotically tracks the whole population as time tends to infinity. Figure 5 depicts the time evolution of vaccination control (15).
When the vaccination strategy (15) is applied, the proportion of the cured population to the total population will soon decline. Figure 6a illustrates the percentage of the total number of people who have recovered through treatment under different situations. As pointed out in [44], if there is no effective control strategy and the disease has spread to 350 days, then percentage of the cured population to the total population will approximate \(0.6\%\) (red line in Fig. 6a). The blue line is obtained from system (1) with no vaccination control. The green line is obtained from system (1) when the vaccination strategy (15) is applied. Because vaccination has been implemented, the number of people cured by treatment will become very small (Fig. 6b). Most susceptible gain immunity directly through vaccination and do not need treatment. This shows that vaccination can greatly reduce the burden on medical facilities.
In addition, if we set \(\varepsilon =0.9\), \(r_1=\beta -\mu \), \(r_2=\mu ^I+\gamma \) and \(r_3=3(\mu ^I+\gamma )\), then the conditions of Theorem 4 do not hold. However, one can see from Fig. 7 that the nonnegativity of system (1) is guaranteed when one applies vaccination strategy (15). This is because the condition of Theorem 4 is a sufficient condition. The vaccination strategy (15) is applied when \(t=3\) days. By direct calculation, \(t_{f}\approx 71\) days. Figure 7a depicts the time evolution of system (1) with vaccination control (15). Figure 7b depicts the number of active infections with control \(u_e(t)\). Figure 8 depicts the time evolution of vaccination control (15).
6 Concluding remarks
This paper introduced an SEIR-type epidemic model with time delay and vaccination control which assumed that the disease is high infectiousness and there exists a large number of asymptomatic infections. We have shown that the disease will persist in population when there does not exist vaccination control and the basic reproduction number \(R_0>1\). The vaccination control was presented as a switching control based on feedback linearization technique. The proposed vaccination control can make the active infections disappear in finite time. Since COVID-19 satisfies the properties of the proposed model (high infectiousness and a large number of asymptomatic infections), we apply the proposed control to simulate the vaccination of COVID-19, and illustrate the effectiveness of the obtained results.
References
Chan, J.F.W., Yuan, S., Kok, K.H., To, K.K.W., Chu, H., Yang, J., Xing, F., Liu, J., Yip, C.C.Y., Poon, R.W.S., Tsoi, H.W., Lo, S.K.F., Chan, K.H., Poon, V.K.M., Chan, W.M., Ip, J.D., Cai, J.P., Cheng, V.C.C., Chen, H., Hui, C.K.M., Yuen, K.Y.: A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating person-to-person transmission: a study of a family cluster. The Lancet 395(10223), 514–523 (2020)
Badr, H.S., Du, H., Marshall, M., Dong, E., Squire, M.M., Gardner, L.M.: Association between mobility patterns and COVID-19 transmission in the USA: a mathematical modelling study. Lancet. Infect. Dis 20(11), 1247–1254 (2020)
Sjodin, H., Wilder-Smith, A., Osman, S., Farooq, Z., Rocklov, J.: Only strict quarantine measures can curb the coronavirus disease (COVID-19) outbreak in Italy, 2020. Eurosurveillance 25(13), 7–12 (2020)
Zhai, S., Gao, H., Luo, G., Tao, J.: Control of a multigroup COVID-19 model with immunity: treatment and test elimination. Nonlinear Dyn. (2020). https://doi.org/10.1007/s11071-020-05961-4
Zheng, X., Luo, S., Sun, Y., Han, M., Liu, J., Sun, L., Zhang, L., Ling, P., Ding, Y., Jin, T., Liu, Z., Weng, J.: Asymptomatic patients and asymptomatic phases of coronavirus disease 2019 (COVID-19): a population-based surveillance study. Natl. Sci. Rev. 7(10), 1527–1539 (2020)
Gao, M., Yang, L., Chen, X., Deng, Y., Yang, S., Xu, H., Chen, Z., Gao, X.: A study on infectivity of asymptomatic SARS-CoV-2 carriers. Respir. Med. 169, 106026 (2020)
Arons, M.M., Hatfield, K.M., Reddy, S.C., Kimball, A., James, A., Jacobs, J.R., Taylor, J., Spicer, K., Bardossy, A.C., Oakley, L.P., Tanwar, S., Dyal, J.W., Harney, J., Chisty, Z., Bell, J.M., Methner, M., Paul, P., Carlson, C.M., McLaughlin, H.P., Thornburg, N., Tong, S., Tamin, A., Tao, Y., Uehara, A., Harcourt, J., Clark, S., Brostrom-Smith, C., Page, L.C., Kay, M., Lewis, J., Montgomery, P., Stone, N.D., Clark, T.A., Honein, M.A., Duchin, J.S., Jernigan, J.A., County, P.H.S.K., Team, C.C.I.: Presymptomatic SARS-CoV-2 infections and transmission in a skilled nursing facility. N. Engl. J. Med. 382(22), 2081–2090 (2020)
Qiu, J.: Covert coronavirus infections could be seeding new outbreaks. Nature (2020). https://doi.org/10.1038/d41586-020-00822-x
Zhu, F.C., Guan, X.H., Li, Y.H., Huang, J.Y., Jiang, T., Hou, L.H., Li, J.X., Yang, B.F., Wang, L., Wang, W.J., Wu, S.P., Wang, Z., Wu, X.H., Xu, J.J., Zhang, Z., Jia, S.Y., Wang, B.S., Hu, Y., Liu, J.J., Zhang, J., Qian, X.A., Li, Q., Pan, H.X., Jiang, H.D., Deng, P., Gou, J.B., Wang, X.W., Wang, X.H., Chen, W.: Immunogenicity and safety of a recombinant adenovirus type-5-vectored COVID-19 vaccine in healthy adults aged 18 years or older: a randomised, double-blind, placebo-controlled, phase 2 trial. The Lancet 396(10249), 479–488 (2020)
Cook, T.M., Farrar, J.J.: COVID-19 vaccines: one step towards the beginning of the end of the global impact of the pandemic. Anaesthesia 76(4), 435–443 (2021)
Tanne, J.H.: COVID-19: Pfizer-BioNTech vaccine is rolled out in US. BMJ 371, 4836 (2020)
Cooke, K.: Stability analysis for a vector disease model. Rocky Mount. J. Math. 9(1), 31–42 (1979)
Al-Darabsah, I.: Threshold dynamics of a time-delayed epidemic model for continuous imperfect-vaccine with a generalized nonmonotone incidence rate. Nonlinear Dyn. 101, 1281–1300 (2020)
Zhao, X.Q.: Basic reproduction ratios for periodic compartmental models with time delay. J. Dyn. Diff. Equ. 29, 67–82 (2017)
Kumar, A., Goel, K.: Nilam: A deterministic time-delayed SIR epidemic model: mathematical modeling and analysis. Theory Biosci. 139(1), 67–76 (2020)
Agaba, G., Kyrychko, Y., Blyuss, K.: Time-delayed SIS epidemic model with population awareness. Ecol. Complex. 31, 50–56 (2017)
Liu, J., Bai, Z., Zhang, T.: A periodic two-patch SIS model with time delay and transport-related infection. J. Theor. Biol. 437, 36–44 (2018)
Zhu, L., Guan, G., Li, Y.: Nonlinear dynamical analysis and control strategies of a network-based SIS epidemic model with time delay. Appl. Math. Model. 70, 512–531 (2019)
Gao, S., Teng, Z., Xie, D.: The effects of pulse vaccination on SEIR model with two time delays. Appl. Math. Comput. 201(1), 282–292 (2008)
De la Sen, M., Alonso-Quesada, S., Ibeas, A.: On the stability of an SEIR epidemic model with distributed time-delay and a general class of feedback vaccination rules. Appl. Math. Comput. 270, 953–976 (2015)
Gao, S., Chen, L., Teng, Z.: Pulse vaccination of an SEIR epidemic model with time delay. Nonlinear Anal. Real World Appl. 9(2), 599–607 (2008)
Wang, W.: Global behavior of an SEIRS epidemic model with time delays. Appl. Math. Lett. 15(4), 423–428 (2002)
Yuan, Y., Belair, J.: Threshold dynamics in an SEIRS model with latency and temporary immunity. J. Math. Biol. 69, 875–904 (2014)
Yang, P., Wang, Y.: Dynamics for an SEIRS epidemic model with time delay on a scale-free network. Phys. A 527, 121290 (2019)
Delacruz, W.P., Savona, M.R., Thornton, J.A., Danaher, P.J.: Evidence of vaccinia dissemination despite lack of major reaction following smallpox vaccination. Vaccine 38(7), 1589–1592 (2020)
Agaba, G., Kyrychko, Y., Blyuss, K.: Dynamics of vaccination in a time-delayed epidemic model with awareness. Math. Biosci. 294, 92–99 (2017)
Farrington, C.: On vaccine efficacy and reproduction numbers. Math. Biosci. 185(1), 89–109 (2003)
Lv, W., Ke, Q., Li, K.: Dynamical analysis and control strategies of an SIVS epidemic model with imperfect vaccination on scale-free networks. Nonlinear Dyn. 99, 1507–1523 (2020)
Zhao, Z., Pang, L., Chen, Y.: Nonsynchronous bifurcation of SIRS epidemic model with birth pulse and pulse vaccination. Nonlinear Dyn. 79, 2371–2383 (2015)
Alonso-Quesada, S., Sen, M.D.L., Agarwal, R.P., Ibeas, A.: An observer-based vaccination control law for an SEIR epidemic model based on feedback linearization techniques for nonlinear systems. Adv. Differ. Equ. 2012(1), 1–32 (2012)
De la Sen, M., Ibeas, A., Alonso-Quesada, S.: On vaccination controls for the SEIR epidemic model. Commun. Nonlinear Sci. Numer. Simul. 17(6), 2637–2658 (2012)
Sisodiya, O.S., Misra, O., Dhar, J.: Dynamics of cholera epidemics with impulsive vaccination and disinfection. Math. Biosci. 298, 46–57 (2018)
Jiao, J., Cai, S., Li, L.: Impulsive vaccination and dispersal on dynamics of an SIR epidemic model with restricting infected individuals boarding transports. Phys. A 449, 145–159 (2016)
Zhao, Z., Chen, L., Song, X.: Impulsive vaccination of SEIR epidemic model with time delay and nonlinear incidence rate. Math. Comput. Simul. 79(3), 500–510 (2008)
Liu, X., Takeuchi, Y., Iwami, S.: SVIR epidemic models with vaccination strategies. J. Theor. Biol. 253(1), 1–11 (2008)
Jiao, H., Shen, Q.: Dynamics analysis and vaccination-based sliding mode control of a more generalized SEIR epidemic model. IEEE Access 8, 174507–174515 (2020)
Hale, J., Lunel, S.: Introduction to Functional Differential Equations, vol. 99. Springer, New York (1993)
Hirsch, W.M., Hanisch, H., Gabriel, J.P.: Differential equation models of some parasitic infections: methods for the study of asymptotic behavior. Commun. Pure Appl. Math. 38(6), 733–753 (1985)
Yang, Z., Zeng, Z., Wang, K., Wong, S.S., Liang, W., Zanin, M., Liu, P., Cao, X., Gao, Z., Mai, Z., Liang, J., Liu, X., Li, S., Li, Y., Ye, F., Guan, W., Yang, Y., Li, F., Luo, S., Xie, Y., Liu, B., Wang, Z., Zhang, S., Wang, Y., Zhong, N., He, J.: Modified SEIR and AI prediction of the epidemics trend of COVID-19 in China under public health interventions. J. Thorac. Dis. 12(3) (2020)
Li, Q., Guan, X., Wu, P., Wang, X., Zhou, L., Tong, Y., Ren, R., Leung, K.S., Lau, E.H., Wong, J.Y., Xing, X., Xiang, N., Wu, Y., Li, C., Chen, Q., Li, D., Liu, T., Zhao, J., Liu, M., Tu, W., Chen, C., Jin, L., Yang, R., Wang, Q., Zhou, S., Wang, R., Liu, H., Luo, Y., Liu, Y., Shao, G., Li, H., Tao, Z., Yang, Y., Deng, Z., Liu, B., Ma, Z., Zhang, Y., Shi, G., Lam, T.T., Wu, J.T., Gao, G.F., Cowling, B.J., Yang, B., Leung, G.M., Feng, Z.: Early transmission dynamics in Wuhan, China, of novel coronavirus-infected pneumonia. N. Engl. J. Med. 382(13), 1199–1207 (2020)
Carcione, J.M., Santos, J.E., Bagaini, C., Ba, J.: A simulation of a COVID-19 epidemic based on a deterministic SEIR model. Front. Public Health 8, 230 (2020)
Lazzerini, M., Putoto, G.: COVID-19 in Italy: momentous decisions and many uncertainties. The Lancet Global Health 8(5), e641–e642 (2020)
CENTER, C.R.: COVID-19 dashboard by the center for systems science and engineering (CSSE) at Johns Hopkins University (JHU). https://coronavirus.jhu.edu/map.html (2020)
Giordano, G., Blanchini, F., Bruno, R., Colaneri, P., Colaneri, M.: Modelling the COVID-19 epidemic and implementation of population-wide interventions in Italy. Nat. Med. 26, 855–860 (2020)
Acknowledgements
This work was supported in part by the Natural Science Foundation of Chongqing of China under Grant cstc2019jcyj-msxmX0109, the Scientific and Technological Research Program of Chongqing Municipal Education Commission under Grant KJQN202000608, the Science and Technology Innovation Project of “Construction of Chengdu-Chongqing Double City Economic Circle” under Grant KJCX2020029.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The authors declare that they have no conflict of interest.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
About this article
Cite this article
Zhai, S., Luo, G., Huang, T. et al. Vaccination control of an epidemic model with time delay and its application to COVID-19. Nonlinear Dyn 106, 1279–1292 (2021). https://doi.org/10.1007/s11071-021-06533-w
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11071-021-06533-w