Abstract
In order to study the impact of continuous vaccination and voluntary isolation for the COVID-19, a susceptible-exposed-infected-recovered-quarantine-vaccines (SEIR-QV) model is proposed. A basic regeneration number is defined to determine the extinction or persistence of the disease. We numerically analyze the impact of key parameters based on actual parameters of COVID-19, such as the vaccination rate, population importation rate, and natural (or causal) mortality transmission rate on the dynamics of disease transmission. Then we obtain sensitivity indices of some parameters on by sensitivity analysis. Finally, the stability of the system and the effectiveness of the optimal control strategy are verified by numerical simulation.
1. Introduction
Since the outbreak of the coronavirus disease 2019 (COVID-19), it has been significantly impacting the world economy and our lives. Globally, as of 4:00 pm CET, 28 February 2022, there have been 434,154,739 confirmed cases of COVID-19, including 5,944,342 deaths, reported to the WHO. As of 26 February 2022, a total of 10,585,766,316 vaccine doses have been administered.
With the emergence of multiple mutant strains of novel coronavirus, including Alpha, Beta, Gamma, Delta, and Omicron, many scholars further explore the transmission and epidemic pattern of COVID-19 from different view. These research results are helpful to learn about the transmission and mode of infectious diseases and provide reliable information for the prediction and control of infectious diseases.
In the study of the dynamics of infectious diseases, the susceptibility-infection-recovery (SIR) model was mainly established by Kermack and McKendrick in 1927 with the kinetic method. In order to better study the characteristics of infectious diseases and obtain the best prevention strategies, scholars propose many improved models based on the traditional SIR model for different situations. By introducing time variation factors, Chen et al. [1] present a new SIR model, which has stronger adaptability and robustness in predicting the number of confirmed cases and inflection points of infectious diseases. Wang et al. [2] propose a SIR model with time-varying isolation protocol. In this model, the dynamic evolution of infectious diseases can be comprehensively analyzed and evaluated. In addition, Hota et al. [3] introduce a closed-loop framework combined with the SIR model and expound the significance of early detection through two feasible optimization questions.
In view of the latent characteristics of some infectious diseases, SEIR models are further proposed based on the SIR model. Li et al. [4] fit the basic regeneration number curve based on the SEIR model. Trends in COVID-19 outbreaks in China, the USA, India, and Iran are predicted and analyzed. The results show that the rapid and efficient isolation measures adopted in China are significant in suppressing COVID-19. In addition, some new models take into account the influence of different factors on infectious diseases, such as the presence of temperature and humidity [5, 6]. It is concluded that COVID-19 is more infectious and lethal at low temperatures and humidity. Besides, many research studies show that the different epidemic prevention measures, such as mandatory isolation, wearing masks, vaccination, and government control policies, have significant effects on preventing the rapid spread of COVID-19 [7–9].
Due to the problem of asymptomatic infected persons in infectious diseases, Yu et al. [10] propose an SEIR-AQ model by considering factors including prevention and control efforts, isolation strategies, and asymptomatic infected persons. The SEIR-AQ model can well anticipate the spread trend of COVID-19, providing technical support for the scientific assessment of the infectious disease situation. Analysis from the SEIR-AQ model shows that prevention and control isolation, medical follow-up isolation, and other measures have a significant inhibitory impact on the spread of COVID-19. In addition, some scholars use deep learning methods to study the impact of the COVID-19 outbreak on human society [11]. Some other scholars have analyzed and studied the transmission trends of the COVID-19 outbreak. The spatiotemporal evolution and transmission trends of the COVID-19 epidemic are analyzed using an ontological modeling approach by Liu et al. [12].
Although control measures, vaccination, latent infectivity, and other factors are considered in these studies, the stability and existence of the model are not discussed from the mathematical view point. Therefore, we establish a new susceptible-exposed-infected-recovered-quarantine-vaccines (SEIR-QV) model based on the latent period and the influence of prevention and control measures of COVID-19. In mathematical theory, the extinction or persistence of a disease can be determined based on the basic reproduction number and the Lyapunov function. We numerically simulate and analyze the impact of key real-world parameters of COVID-19. What’s more, our sensitivity analysis indicates that the related parameters have significant effects on the stability and existence of the SEIR-QV model.
2. The SEIR-QV Model
In this paper, the susceptible-exposed-infectious-recovered-asymptomatic-quarantine-vaccines (SEIR-AQV) model is developed based on the traditional SEIR model. The SEIR-AQV model takes into account the effect of isolation based on the traditional model and divides the population into unisolated susceptible people (), unisolated exposed people (), unisolated infected people (), isolated susceptible people (), isolated exposed people (), asymptomatic infected people (), hospitalized people (), and recovered people ().
The SEIR-AQV model has the following assumptions: first assumption: medical resources are adequate, i.e., isolated exposed people can be directly converted to hospitalized people after diagnosis; second assumption: once infected with COVID-19, one must undergo inpatient treatment to recover, i.e., infected people are unlikely to recover on their own; third assumption: both isolated exposed people and hospitalized people are isolated from the outside world, i.e., they are not infectious; fourth assumption: the infected people will be immune after recovery, i.e., he or she will not become susceptible again after recovery.
The bin transformation relationship is shown in Figure 1.
The equations of the SEIR-AQV model are as follows:where means the rate of exposure; means the natural mortality rate; represents the proportion of isolated; represents the vaccine coverage rate; represents the probability of transmission; stands for the ratio of relative to the transmission capacity of ; stands for the ratio of relative to the transmission capacity of ; stand for the conversion of , , and, respectively to ; represent the rate of conversion of into ; stands for the rate of conversion of into ; stands for the ratio of to conversion and ; means the ratio of to conversion; mean the rate of conversion of , , and the rate of cause-specific death of ; means the rate of recovery of ; and means the rate of population importation.
However, given the reality of the COVID-19 epidemic, the status of isolated susceptible people is unlikely to persist over time, and isolated susceptible people have no translational relationship to hamlets other than unisolated susceptible people. Therefore, the population type of isolated susceptible people and the parameters involved in it are removed. Considering also that the proportion of asymptomatic infected people is very small and has minimal effect on the overall transmission trend of the COVID-19 outbreak, the population type of asymptomatic infected people and its involved parameters are also removed. Therefore, we improved the SEIR-AQV model in the context of the COVID-19 outbreak. We construct a new seasonal susceptible-exposed-infected-removed-quarantine-vaccines (SEIR-QV) model with the population transformation relationship bin view shown in Figure 2.
The equations of the SEIR-QV model are as follows:where represent the rate of cause-specific death of unisolated infected and hospitalized , respectively, and represent the rate of conversion of unisolated infected and isolated exposed to hospitalized , respectively.
2.1. Proof of Stability of the SEIR-QV Model
To study the stability of equation (2), it is sufficient to study the stability of the following three formulas in (2) [13, 14]:
Considering the biological significance of the system, the dynamical properties of equation (3) are discussed mainly in the closed set , where denotes the first trigonometric limit of and contains the boundary.
2.1.1. Existence of Equilibrium Point
The equilibrium point of (3) satisfies
When , (3) has the point of diseased equilibrium .
When , from , we get ; from , we get , from , we get , so .
Bring and into to obtain the endemic equilibrium point , where
2.1.2. Basic Regeneration Number
The basic regeneration number () indicates the number of people infected by a patient during the average disease period when all are susceptible at the beginning of the disease. can be used as a threshold to decide whether the disease is extinguished or not. When , the disease will become extinct. When , the disease will persist. The basic regeneration number is closely related to the stability of the endemic equilibrium point.
Next, we study the basic regeneration number of equation (3). Let , then equation (3) can be rewritten aswhere
The Jacobi matrices of and at the disease-free equilibrium are, respectively,where
The basic regeneration number, denoted by , is thus given by the following equation:
2.1.3. Proof of Stability of the SEIR-QV Model
Theorem 1. For equation (3), the disease-free equilibrium is locally asymptotically stable if .
Proof. Linearizing equation (3) at the disease-free equilibrium point , we obtain the linearization matrix at point as the following equation:The characteristic equation of this matrix is , where is a unit matrix. Expanding it givesObviously, this characteristic equation has a negative characteristic root . The other characteristic roots satisfy the following equation:It is known that .
Let , thenSo we can obtainIt is easy to verify for and that the roots of the quadratic equation are as follows:All have negative real parts, i.e., the disease-free equilibrium point of equation (3) is proved to be locally asymptotically stable.
Theorem 2. For equation (3), the disease-free equilibrium is globally asymptotically stable if .
Proof. By considering the Lyapunov function,Clearly, the solution of along equation (3) has and when and only when and . The full derivative of the solution of the function along equation (3) is as follows:Take , .
Then,where the inequality sign is obtained based on . When , there is . Thus, when and only when . Therefore, if , the maximum tight invariant set is the single point set . Therefore, if , the disease-free equilibrium point of equation (3) is globally asymptotically stable in .
Theorem 3. are positive when .
Proof. It is known that the parameters are positive and , so obviously holds the following equation:It is known that , and the parameters are positive, so .It is known that , and the parameters are positive, so .
Theorem 4. For equation (3), the endemic equilibrium is locally asymptotically stable if .
Proof. Linearize equation (3) at the endemic equilibrium point and obtain the linearization matrix at point as follows:The characteristic equation of this matrix is , where is a unit matrix. Expanding it gives , whereThen,where and .
Due to and , .
Therefore, the endemic equilibrium point of equation (3) is locally asymptotically stable on .
Theorem 5. For equation (3), the endemic equilibrium is globally asymptotically stable if .
Proof. Considering the Lyapunov function, we can obtainThe full derivative of the solution of the function along equation (3) is as follows:Bringing in , yieldsTake and .
Then,Since and when , .
At this time,and when and , , and holds, so . Therefore, if , the endemic equilibrium point of equation (3) is globally asymptotically stable on .
3. Numerical Experiments
3.1. Numerical Simulation
The stability of disease-free and endemic equilibrium points is discussed in this section, and this section uses the software Matlab to numerically simulate equation (3) to verify the above conclusions.
The initial values of the system are (0.7,0.2,0.1), and the parameters are assigned as , , , , , , , and , where and are two variables. In this paper, we will take some random values to simulate the change of , , and ratio when and , respectively. When , takes a random value between 0 and 1 so that holds. When , takes a random value between 0 and 1 so that holds. Its value of R0 is less than 1 () when q and delta take the values in Table 1. Its value of R0 is greater than 1 () when q and delta take the values in Table 2.
The global asymptotic stability of both disease-free and endemic equilibrium points can be seen from Theorems 1–5. Figures 3–6 illustrates the correctness of the results obtained from the above theorem.
(a)
(b)
(c)
(a)
(b)
(c)
(a)
(b)
(c)
(a)
(b)
(c)
If , by numerically simulating the stability of the disease-free equilibrium point, we can analytically conclude the following points. Figure 3 shows the variation curves of , , and as the parameter varies when is constant. The variation of increases, and the rate of decrease and the rate of increase both increase with the increase of . The rate of decrease of and decreases with the increase of . Finally, , , and all stabilized (see Figures 3(a)–3(c)). The change of decreases, and the rate of decrease of decreases with the increase of . But the change of has no obvious pattern with the increase of . Finally, , , and all tend to be stable (see Figures 4(a)–4(c)). Figures 3 and 4 illustrate that when , the changes of and have significant effects on the trends of , , and . And when , the disease-free equilibrium point is tending to be stable.
In Figures 5 and 6, the curve of , , and variation with and is simulated. By comparing the trend of , , and with in Figure 5 and the trend of , , and with in Figure 6, we can draw the following conclusions. Increases in both and lead to larger fluctuations in and delayed stabilization times. It will lead to an increase in the value of as it stabilizes. What is more, as increases, the rate of increase and decrease of decreases. The monotonically increasing interval of decreases and the rate of decrease of the monotonically decreasing interval decreases with the increase of . The maximum peak value of both and decreases with the increase of . And the time of reaching the maximum peak value of is delayed and the time of reaching the maximum peak value of is advanced. Furthermore, as increases, the rate of increase and decrease of decreases and reaches its maximum peak earlier and at a lower value. As changes, there is no obvious pattern of change in , but it can be concluded that has a significant effect on the change in the value of . Figures 5 and 6 illustrate that when , the changes of and have significant effects on the trends of , , and . And when , the endemic equilibrium point is tending to be stable.
3.2. Model Simulation
The real data of China in this paper are from the official website of the WHO. We preprocess the data, remove some bad points, and make the time continuous so that the data can be better adapted to the SEIR-QV model. The Chinese COVID-19 outbreak is fitted based on the SEIR-QV model. The parameters used in the model are listed as follows: , , , , , , , , , , , , and . The simulation results and the actual data are shown in Figure 7.
The Chinese government has decisively taken strict preventive and control measures following the emergence of the COVID-19 outbreak. These prevention and control measures, along with vigorous advocacy by the Chinese government, have reduced exposure rates and increased isolation rates and minimized the development of aggregated outbreaks. The COVID-19 outbreak has become manageable and stable. As can be seen in Figure 7, the simulated data fit better with the real data from January 20, 2020 to May 20, 2021. However, there is a small outbreak of the epidemic in China in June 2021. The main reason is that the epidemic is out of control in many places abroad, and people at home are taking it lightly and weakening the efforts of prevention and control, thus giving the opportunity for the epidemic to spread.
As can be seen in Figure 8, both the isolation ratio and vaccine coverage play an important role in the control of the outbreak. The highest number of confirmed cases per day in China decreases as both and increase.
(a)
(b)
The correlation coefficient is used to measure the accuracy of the fit, and the value of R is found to be 0.9394. Therefore, the results show that the SEIR-QV model can be used for COVID-19 epidemic development status assessment and has important implications.
4. Conclusion
An SEIR-QV model based on vaccination, isolation strategies, and the impact of different parameters on the development of infectious diseases is developed. The equilibrium point and stability of the new model are proved by using the basic regeneration number and the Lyapunov function theory. Simulation experiments show that the new method has certain theoretical value for analyzing and predicting the development of the COVID-19 epidemic. We can get the conclusion that there is a significant impact on the development of infectious diseases by different vaccination and isolation strategies.
Data Availability
The dataset used to support the findings of this study is available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (no. 62171203).