Abstract
Background In December 2019, a novel coronavirus-induced pneumonia (COVID-19) broke out in Wuhan, China. On 19 January and on 8 June 2020, there were two wave COVID-19 epidemics happened in Beijing. Modelling, simulations and analysis for the two wave epidemics are important issues.
Methods This study introduces a symptomatic-asymptomatic-recoverer-death differential equation model (SARDDE). It presents the conditions of the asymptotical stability on the disease-free equilibrium of the SARDDE. It proposes the necessary conditions of disease spreading for the SARDDE. Based on the reported data of the first and the second COVID-19 epidemics in Beijing and numerical simulations, it determines the parameters of the SARDDE, respectively.
Results Numerical simulations of the SARDDE describe well the outcomes of the current symptomatic and asymptomatic individuals, the recovered symptomatic and asymptomatic individuals, and the died infected individuals, respectively. The numerical simulations obtain the following results.
The transmission rate of the symptomatic infections caused by the symptomatic individuals in the second Beijing epidemic is about two times higher than the one in the first Beijing epidemic.
Both the symptomatic and the asymptomatic individuals cause lesser asymptomatic spread than symptomatic spread.
The blocking rates of 89.32% and 97.48% (reaching the infection turning points) to the symptomatic infections cannot prevent the spreads of the first and the second COVID-19 epidemics in Beijing, respectively.
That on the day 28, the symptomatic infection blocking rates reached to 100% has made the second Beijing epidemics epidemic end on day 56.
That on the day 98, the symptomatic infection blocking rates reached to 100% has made the the first Beijing epidemics epidemic end on day 140.
Keeping the blocking rates, the recovery rates and the death rates reaching the infection turning points would make the numbers of current hospitalized infected individuals reach, on day 140, 208958 individuals and 25017 individuals for the two Beijing epidemics, respectively.
Conclusions Virtual simulations suggest that the strict prevention and control strategies implemented by Beijing government are effective and necessary; using the data from the beginning to the days after about two weeks after the turning points can estimate well and approximately the following outcomes of the two COVID-19 epidemics, respectively. To avoid multiple epidemic outbreaks, a recommendation is that the authorities need to have maintained the prevention and control measures implemented, at least, 7 days after reaching the turning point until new current infection cases disappear. It is expected that the research can provide better understanding, explaining, and dominating for epidemic spreads, prevention and control measures.
1 Introduction
In December 2019, a novel coronavirus-induced pneumonia (COVID-19) broke out in Wuhan, China. As of 15 July 2022, over 557 million confirmed cases of COVID-19 and over 6.35 million deaths have been reported globally [1]. Now COVID-19 affects more than 220 countries and regions including Antarctica.
One of the reasons of such a tragedy is that people in some countries do not pay attentions to theoretical analysis and estimations for COVID-19 epidemics. In fact mathematical models for epidemic infectious diseases have played important roles in the formulation, evaluation, and prevention of control strategies. Modelling the dynamics of spread of disease can help people to understand the mechanism of epidemic diseases, formulate and evaluate prevention and control strategies, and predict tools for the spread or disappearance of an epidemic [2].
Since the outbreak of COVID-19 in Wuhan, many scholars have published a large numbers of articles on the modeling and prediction of COVID-19 epidemic (for examples see [3–15]). It is difficult to describe well the dynamics of COVID-19 epidemics. In a Lloyd-Smith et al’s paper, it described nine challenges in modelling the emergence of novel pathogens, emphasizing the interface between models and data [16].
On 19 January, two Beijingers returning from Wuhan were diagnosed with COVID-19. That triggered the first wave of COVID-19 in Beijing. During the first wave of COVID-19, a total of 420 locally diagnosed cases were reported. After 140 days, on 8 June, 411 COVID-19 infected individuals were cured and 9 died [17]. However three days later, there was one new confirmed individual infected a kind of different COVID-19 in an agricultural products wholesale market in Tongzhou district, Beijing. This triggered a second wave of COVID-19 epidemic in Beijing. After 56 days, totals of 335 and 50 locally symptomatic and asymptomatic COVID-19 individuals were reported, and all individuals were cured. Medical staff have achieved zero infection [17].
This paper introduces a SARDDE. It gives the conditions of the asymptotical stability on the disease-free equilibrium of the SARDDE. Using simulations determines the parameters of the SARDDE based on the reported data of the two COVID-19 epidemics in Beijing [17]. Numerical simulations of the SARDDE describe well the practical outcomes of current infected symptomatic and asymptomatic individuals, recovered infected symptomatic and asymptomatic individuals, and died infected individuals. Virtual simulations are given to estimate the effectiveness of the prevention and control strategies. This paper improves the simulation results implemented in the previous version of this paper [18], and increases a comparing subsection.
The rest of this article is organized as follows. Subsections 2.2.1 and 2.3.1 establish the SARDDEs with 5 and 4 variables, respectively. Subsections 2.2.2 and 2.3.2 provide the criterions of the asymptotical stability of the disease-free equilibriums of the two SARDDEs. Subsections 2.2.3 and 2.3.3 determine the necessary conditions of disease spreading. Subsections 3.1.1 and 3.2.1 implement the dynamic simulations of the SARDDEs to describe the reported data of the first and second COVID-19 epidemic in Beijing, respectively. Subsections 3.1.2 and 3.2.2 discuss and analyze the simulation results. Subsections 3.1.3 and 3.2.3 implement virtual simulation examples to predict the outcomes of the two Beijing epidemics on different prevention and control measures. Subsection 3.3 compares the two COVID-19 epidemics in Beijing. Section 4 concludes the paper. A recommendation to avoid multiple epidemic outbreaks is proposed in the same Section.
2 SARDDEs and Dynamic Properties
2.1 Data collection
The original COVID-19 data sets of the first and second epidemics in Beijing were downloaded from the Beijing authority’s Website [17]. The author edited the original COVID-19 data sets for using in this study.
2.2 SARDDE with death variable
2.2.1 Model
For SARDDE model with death variable, there are five states. I(t), Ia(t), Ir(t), Ira(t) and D(t) represent the fraction of current symptomatic infected individuals, and current asymptomatic infected individuals, cumulative recovered symptomatic infected individuals, cumulative recovered asymptomatic infected individuals and cumulative died individuals, respectively. The transition among these states is governed by the following rules (Flowchart of the rules is shown in Fig.1, where S represents susceptible population.).
Flowchart of disease transmission among susceptible population S, current symptomatic infected individuals I, current infected individuals Ia, recovered symptomatic infected individuals Ir, recovered asymptomatic but infected individuals Ira, and died individuals D.
First, the symptomatic infected individuals (I) and the asymptomatic infected individuals (Ia) infect the susceptible population (S) with the transmission rates of β11 and β21, respectively, making S become symptomatic infected individuals, and with the transmission rates of β12 and β22, respectively, making S become asymptomatic individuals. Then, a symptomatic individual is cured at a rate κ, an asymptomatic individual returns to normal at a rate κa. An infected individual dies at a rate α. Here all parameters are positive numbers.
Assume that the dynamics of an epidemic can be described by m-time intervals, which correspond different prevention and control measures, and medical effects. At lth time interval, the SARDDE has the form:
where Θ1(l) = (1 − θ1(l)) and Θ2(l) = (1 − θ2(l)) (l = 1, …, m) represent the blocking rates to symptomatic and asymptomatic infections, respectively. Then SARDDE (1) has a disease-free equilibrium:
2.2.2 Stability of disease-free equilibrium
The stability of SARDDE (1) is determined by the first two equations (1a) and (1b). Denote in (1a) and (1b)
Then at the disease-free equilibrium of SARDDE (1), the Jacobian matrix of (1a) and (1b) is
Solving the corresponding eigenequation obtains two eigenvalues:
Therefore we obtain the following:
Suppose that a11, a12, a21 and a22 are defined by (3)-(6). Then the disease-free equilibrium E (2) of SARDDE (1) is globally asymptotically stable if, and only if, the following inequalities hold:
2.2.3 The necessary condition of disease spreading
If an epidemic can occur, then
This implies that
Solving the above inequalities gives the following
If SARDDE (1) satisfies the following inequalities
then a disease transmission will occur.
2.3 SARDDE Model without death variable
2.3.1 Model
Similar to Section 2.2.1, the transition among these states is governed by the following rules (Flowchart of the rules is shown in Fig. 2, where S represents susceptible population.)
Flowchart of disease transmission among susceptible population S, current symptomatic infected individuals I, current asymptomatic infected individuals Ia, recovered symptomatic infected individuals Ir, recovered asymptomatic infected individuals Ira.
Assume that the dynamics of an epidemic can be described by m time intervals. At lth time interval, SARDDE (1) becomes the form:
Then SARDDE (11) has a disease-free equilibrium:
2.3.2 Stability of disease-free equilibrium
The stability of SARDDE (11) is determined by the first two equations (11a) and (11b). Denote in (11a) and (11b):
Then at the disease-free equilibrium of SARDDE(11), the Jacobian matrix of (11a) and (11b) has the form
Solving the corresponding eigenequation obtains 2 eigenvalues:
Therefore we obtain the following:
Suppose that a11, a12, a21 and a22 are defined by (13)-(16). Then the disease-free equilibrium E of SARDDE (11) is globally asymptotically stable if, and only if, the following inequalities hold:
2.3.3 The necessary condition of disease spreading
If an epidemic can occur, then
This implies that
Solving the above inequalities gives the following
If SARDDE (11) satisfies the following inequalities
Then a disease transmission will occur.
3 Applications
3.1 Simulation and prediction of the first COVID-19 epidemic in Beijing
This Section will discuss the applications of Theorem 1, Theorem 2 and SARDDE (1) to simulate the real-word COVID-19 epidemic data from 19 January to 8 June 2020 in Beijing [17]. Numerical simulations and drawings are performed by using MATLAB software programs.
3.1.1 Modeling and Simulation
Figures 3(a) and 3(b) show the first 55 days’ reported clinical data on the current confirmed infection cases, and the reported clinical data on recovered cases of the COVID-19 epidemic in Beijing [17]1. The number of the current symptomatic infected individuals is shown in Fig4(a) by circles. The numbers of the cumulative recovered symptomatic infected individuals, and the cumulative died infected individuals are shown in Fig4(b) by circles and stars respectively.
(a) Outcome of the number of the current infected individuals. (b) Outcome of the number of the cumulative recovered individuals.
The number of the current infected individuals was risen rapidly in the first 4 days (see Fig.3(a) and Fig. 4(a)). The number of the current infected individuals reached the highest 295 on day 24, February 12. Then after day 31, February 19, they declined rapidly (see Fig.3(a) and 4(a)).
(a) Outcome of the number of the current symptomatic individuals (CSI), representing by circles. Solid lines and dash lines are outcomes of the corresponding simulations of SARDDE (1). (b) Outcomes of the numbers of the cumulative recovered symptomatic individuals (CRSI) and the cumulative died individuals (CDI), representing by circles and stars, respectively. Solid lines and dash lines are the corresponding simulations of SARDDE (1).
On day 86 (15 April), there were 3 Chaoyang district infected people coming back Beijing form foreign country which made calculated blocking rate rise. After day 97, the number of the current infected individuals were reducing continually (in piecewise linear way) until day 140 (8 June) reaching to zero.
Observe from Figs. 4(a) and 4(b) that the overall changes in the number of the current confirmed infections are not subject to the law of exponential changes, but the data can be approximated in good agreement with 9 straight lines in log scale (see Fig. 4). This phenomenon can be explained as different medical measures prevention and control strategies have been adopted over the different 9-time intervals. Therefore the l in SARDDE model (1) should be chosen as l = 1, 2, …, 9.
Denote
[tl−1, tl] is the lth time interval. Denote Ic(tl) to be the number of the reported current symptomatic individuals, and Denote Icr(tl) to be the number of the reported current cumulative recovered symptomatic individuals. Dc(tl) be the number of the reported current cumulative died individuals. Using the minimization error square criterion determines the βij(l)′s, κ(l)′s, κa(l)′s), θ1(l)′s, θ2(l)′s, and α(l)′s.
(1) One can assume that S = 1 because the effects of S can be deleted by calculated 
. This makes the calculated 
have general sense.
(2) In the previous version of this paper [18], κ(l), α(l) were determined via using the number of the recovered patients and died patients to divide the days of patients stayed in the hospital during the lth time interval, respectively. However It finds that with the helps of the formulas κ(l), and α(l) [18], we determine the search ranges of κ(l) and α(l) and can obtain better simulation results as those described later.
(3) The website [17] did not reported the data of asymptomatic infected individuals and announced only that there were no asymptomatic infected individuals after day 73 (April 1). Hence we cannot estimate the errors between the simulated Ia(l) and Ira(l) with the reported ones.
(4) The search ranges of the infection rate 
s are difficult to be determined. However if we know the number of current symptomatic infected individuals at day 0 and day T when no prevention and control measures had been implemented. Additional the number of current asymptomatic infected individuals was very small. Then we can estimate 
from equation (1a) as follows
Using the real world data of the first COVID-19 epidemic in Beijing [17] selects following initial condition
The calculated parameters are shown in Table 1. The corresponding simulation results of equation (1) are shown in Fig. 1 and Fig. 2. Observe that the simulation results of equation (1) were in good agreement with the data of the COVID-19 epidemics. At the end points (see solid dots in Figs. 1 and 2) of the 9 investigated time-interval [tl−1, tl]′s, the simulated numbers and the actual reported numbers were approximate the same∼ errors were less than one, except on day 140, the numbers of the simulated and the reported cumulative recovered infected individuals was one difference. (See Figs. 4(a) and 4(b) and legends in them).
The data of the first COVID-19 epidemic in Beijing on different days and corresponding calculated SARDDE parameters. Where NCSII and NCDI represent the numbers [17] of the current symptomatic infected individuals and the cumulative died individuals, respectively; NRSII the number [17] of the recovered symptomatic infected individuals over the lth time interval.β11 = 0.75454, β21 = 0.1119, β12 = 0.48757, β22= 0.07993.
3.1.2 Discussions
On the end points of the 9 investigated time-interval [tl−1, tl]′s, that is on day 0, day 3, day 12, day 24, day 30, day 243, day 73, day 97, day 116, and day 140, we obtain the following results (also see Figs. 4).
On the end point days, the reported and simulated current symptomatic individuals are are approximate the same.
On day 140, the numbers of the reported and simulated cumulative recovered individuals were one difference. On the other end point days, the numbers were the same.
On the end point days, the numbers of the reported and simulated cumulative died individuals are approximate the same.
There is no information on the current symptomatic infected and recovered symptomatic infected individuals. But it has reported that after the 73 day, April 1, there is no symptomatic infected individuals until day 143, June 11 [17]. Our simulation results shows that on day 73, the number of the simulated current symptomatic infected individuals was less than one (≈ 0.43612), the cumulative asymptomatic infected individuals discharged from observations were about 104 individuals, which seem to explain the actual missing report asymptomatic infected individual data.
Computed results (see Table 1) of the transmission rates
show that the ratio of the transmission rates of asymptomatic and symptomatic individuals infecting susceptible population to become symptomatic individuals is about 0.159 (β21:β11). It suggests that the asymptomatic individuals cause lesser symptomatic spread than the symptomatic individuals do.The computed results (see Table 1) also show that the ratios of the transmission rates of the asymptomatic and the symptomatic individuals infecting susceptible population to become the asymptomatic and the symptomatic individuals are about 0.646 (β12 : β11) and 0.667 (β22 : β21), respectively. It suggests that both the symptomatic and the asymptomatic individuals cause lesser asymptomatic spreads than symptomatic spreads.
The criteria (7) and (8) of the asymptotical stability of the disease-free equilibrium of SARDDE model (1) over the 8-time intervals are shown in the 5th ∼ 8th columns of Table 2. It is shown that until the blocking rates (Θ1, Θ2) reach about (98.44%, 99%), the disease-free equilibrium becomes globally asymptotical stability. The conditions (9) and (10) of disease spreading are listed in the last two columns in Table 2. It shows also that if the blocking rates (θ1, θ2) reach about (98.44%, 99%), the spreading of COVID-19 epidemic can be blocked. Table 2 suggests also that 90% and 96% blocking rates to symptomatic infections and asymptomatic infections cannot prevent spreads to transmissions of symptomatic infections and asymptomatic infections.
Criteria of the asymptotical stability and disease spreading of disease-free equilibrium of SARDDE (1) over 9-time intervals.
3.1.3 Virtual Simulations
Now assume that after day 24, 12 February, it still keeps the blocking rates (Θ1(3), Θ2(3)) ≈ (89.93%, 96%), the recoverr rates (κ(3), κa(3)), and the died rate α(3) until day 140, 8 June. The simulation results of SARDDE (1) are shown in Figs 5(a) and 5(b). Observe that the numbers of current symptomatic and the asymptomatic infected individuals reach about 208958 and 22847, respectively. The numbers of cumulative recovered symptomatic and died individuals reach about 90436 and 3028, respectively.
Virtual simulations: (a) Outcome of the number of the current symptomatic individuals (CSI), representing by circles. Solid lines and dash lines are outcomes of the corresponding simulations of SARDDE (1). (b) Outcomes of the numbers of the cumulative recovered symptomatic individuals (CRSI) and the cumulative died individuals (CDI), representing by circles and stars, respectively. Solid lines and dash lines are the corresponding simulations of SARDDE (1).
Furthermore assume that after day 52, 11 March, it still keeps the blocking rates (Θ1(5), Θ2(4)) ≈ (98.93%, 99%), the recovery rates (κ(5), κa(5)), and the died rate α(5) until day 140, 8 June. The simulation results of SARDDE (1) are shown in Figs 6(a) and 6(b). Observe that the numbers of current symptomatic and asymptomatic infected individuals are about one and zero, respectively; the numbers of cumulative recovered symptomatic and died individuals are about 415 and 11, respectively. The results suggest that using the data before day 52 (about 17 days after the turning point) can estimate well the following outcome of the first COVID-19 epidemic in Beijing.
Virtual simulations: (a) Outcome of the number of the current symptomatic individuals (CSI), representing by circles. Solid lines and dash lines are outcomes of the corresponding simulations of SARDDE (1). (b) Outcomes of the numbers of the cumulative recovered symptomatic individuals (CRSI) and the cumulative died individuals (CDI), representing by circles and stars, respectively. Solid lines and dash lines are the corresponding simulations of SARDDE (1).
In summary, SARDDE (1) can simulate the outcomes of the first COVID-19 epidemic in Beijing. The calculated equation parameters can help us to understand and explain the mechanism of epidemic diseases and control strategies for the event of the practical epidemic.
3.2 Simulation and prediction of the second COVID-19 epidemic in Beijing
This Section will discuss the applications of Theorem 3, Theorem 4 and SARDDE (11) to simulate the real-word COVID-19 epidemic data from 11 June to 6 August 2020 in Beijing [17]. This event of the second epidemic in Beijing provides a valuable example of accurate preventing and controlling strategies and excellent clinical treatments. The data of the symptomatic and asymptomatic infected and recovered individuals provide a footstone to modeling the event.
3.2.1 Modeling and Simulation
Figure 7(a) show that the reported data on the current confirmed symptomatic infection cases. Figure 7(b) show that the reported data on the cumulative recovered symptomatic infection cases. Figure 8(a) show that the reported data on the current confirmed asymptomatic infection cases. Figure 8(b) show that the reported data on the cumulative recovered asymptomatic infection cases. All data set were edited from [17].
(a) Outcome of the number of the current symptomatic infected individuals. (b) Outcome of the number of the cumulative recovered symptomatic infected individuals.
(a) Outcome of the number of the current asymptomatic infected individuals. (b) Outcome of the number of the cumulative recovered asymptotic infected individuals.
The evolution of the current symptomatic infected individuals, and the current asymptomatic infected individuals are shown in Fig. 9(a) by circles and diamonds, respectively. The evolution of the cumulative recovered symptomatic infected individuals, and the cumulative recovered asymptomatic infected individuals are shown in Fig. 9(b) by circles and diamonds, respectively.
Outcomes of the numbers of: (a) the current symptomatic individuals (CSI) and the current asymptomatic individuals (CAI), representing by circles and diamonds. Solid lines and dash lines are the corresponding simulations of SARDDE (11). Outcomes of the numbers of: (b) the cumulative recovered symptomatic individuals (CRSI) and the cumulative recovered asymptomatic individuals (CRAI), representing by circles and diamonds, respectively. Solid lines and dash lines are the corresponding simulations of SARDDE (11).
Observe from Fig.9 that the overall changes in the number of the current confirmed infections are not subject to the law of exponential changes, but the data can be approximated in good agreement with 8 straight lines in log scale (see Fig.9). This phenomenon can be explained as different medical measures and prevention and control strategies have been adopted over different 8-time intervals. Therefore the l in SARDDE (11) satisfies i = 1, 2, …, 8.
Denote
Denote Ic(tl) to be the number of the reported current symptomatic individuals, and Ia(tl) be the number of the reported current asymptomatic individuals charged in medical medical observations. Denote Icr(tl) to be the number of the reported current cumulative recovered symptomatic individuals. and Icra(tl) be the number of the reported current cumulative asymptomatic individuals discharged from medical medical observations.
Using the minimization error square criterion (23) determines the βij(l)′s, κ(l)′s, κa(l)′s), θ1(l)′s and θ2(l)′s.
The calculated results are shown in Table 3. The corresponding simulation results of SARDDE (11) are shown in Fig. 9(a) and 9(b). Observe that the simulation results of model (11) describe well the dynamics of the second COVID-19 epidemic in Beijing.
The data of the second wave COVID-19 epidemics on different days and the corresponding calculated parameters of SARDDE model (11). Where NCSII and NCAII represent the numbers of the current symptomatic infected individuals and the current asymptomatic infected individuals, respectively; NRSII and NRAII represent the numbers of the cumulative recovered symptomatic infected individuals and the asymptomatic infected individuals over the ith interval. β11 = 1.4521, β21 = 0.072824, β12 = 0.072824, β22 = 0.72824.
(1) Similar the description in Section 3.1.1, we have
Using the practical data of the second COVUD-19 epidemic (also see the second line in Table 3) selects following initial condition
The simulations of SARDDE (11) with the above model parameters are shown in Figs. 9(a) and 9(b). Observe that the simulation results are in good agreement with the reported clinical data On the nine time interval endpoint days, the errors between the reported data and the simulated data were less than one (see the solid and dash lines and legends in Figs. 9(a) and 9(b)).
3.2.2 Discussions
On the end points of the 8 investigated time-interval [tl−1, tl]′s, that is on day 0, day 3, day 10, day 19, day 27, day 34, day 50, and day 55, we obtain the following results (also see Figs. 9 and legends).
On the end point days, the numbers of the reported and simulated current symptomatic individuals were approximate the same.
On the end point days, the numbers of the reported and simulated current asymptomatic individuals were approximate the same.
On the end point days, the numbers of the reported and simulated cumulative recovered symptomatic individuals were approximate the same.
On the end point days, the numbers of the reported and simulated cumulative recovered asymptomatic individuals were approximate the same.
On the end point days, the numbers of the reported and simulated cumulative died individuals are approximate the same.
Computed results (see Table 3) show that the ratio of the transmission rates of the asymptomatic and the symptomatic individuals infecting susceptible population to become the symptomatic individuals is about 5% (β21 : β11). It suggests that the asymptomatic individuals cause lesser symptomatic spread than the symptomatic individuals do.
Computed results (see Table 3) also show that the ratios of the transmission rates of the asymptomatic and symptomatic individuals infecting susceptible population to become the asymptomatic and symptomatic individuals are about 5% (β12:β11). It suggests that the symptomatic individuals cause lesser asymptomatic spread than symptomatic spread.
The criteria of the stability of the disease-free equilibrium of SARDDE (11) on 8-time intervals are shown in Table 4. It shows that the blocking rates reach about (97.47%, 96.6%) cannot prevent the spread of the second COVID-19 epidemic in Beijing; the blocking rates (Θ1, Θ2) reach about (99.80%, 98.19%), the disease-free equilibrium becomes globally asymptotical stable.
Criteria of the asymptotical stability and disease spreading of disease-free equilibrium of SARDDE (11) over 8-time intervals.
3.2.3 Virtual Simulation
Now assume that it keeps still the blocking rates (Θ1(3), Θ2(3)) ≈ (97.47%, 96.6%) and the cure rates (κ(3), κa(3)) until day 56, 6 August. The simulation results of SARDDE (11) are shown in Fig.10. Observe that on day 56, the numbers of the current symptomatic and asymptomatic infected individuals reach about 1230 and 71, respectively; The numbers of the cumulative recovered symptomatic and asymptomatic infected individuals reach 23 and 67, respectively. If kept the above equation parameters until day 140, the numbers of the current symptomatic and asymptomatic infected individuals reach about 25017 and 1312, respectively; the numbers of the cumulative recovered symptomatic and asymptomatic infected individuals reach 574 and 1190, respectively.
Virtual simulations. Outcomes of the numbers of: (a) the current symptomatic individuals (CSI) and the current asymptomatic individuals (CAI), representing by circles and diamonds. Solid lines and dash lines are the corresponding simulations of SARDDE (11). Outcomes of the numbers of: (b) the cumulative recovered symptomatic individuals (CRSI) and the cumulative recovered asymptomatic individuals (CRAI), representing by circles and diamonds, respectively. Solid and dash lines are the corresponding simulations of SARDDE (11).
Furthermore assume that after the day 34th, 14 July, it still keeps the blocking rates (Θ1(5), Θ2(5)), the cure rates (κ(5), κa(5)) until day 56, 6 August. The simulation results of SARDDE (11) are shown in Figs. 11(a) and 11(b). Observe that on day 56, the numbers of the current symptomatic and asymptomatic infected individuals reach about 36 and 5, respectively. The numbers of the cumulative recovered symptomatic and the asymptomatic individuals are about 300 and 46, respectively. The results suggest that using the data before the day 34 (about two weeks after the turning point) can approximately estimate the following outcome of the second COVID-19 epidemic in Beijing.
Virtual simulations. Outcomes of the numbers of: (a) the current symptomatic individuals (CSI) and the current asymptomatic individuals (CAI), representing by circles and diamonds. Solid lines and dash lines are the corresponding simulations of SARDDE (11). Outcomes of the numbers of: (b) the cumulative recovered symptomatic individuals (CRSI) and the cumulative recovered asymptomatic individuals (CRAI), representing by circles and diamonds, respectively. Solid and dash lines are the corresponding simulations of SARDDE (11).
3.3 Comparing
The model main parameters of the first and second COVID-19 in Beijing and related data are listed in Table 5. It follows:
Model parameters and related data of the first and second COVID-19 epidemics in Beijing.
The transmission rate β11 of the symptomatic infections caused by the symptomatic individuals in the second Beijing epidemic is much higher than the one in the first Beijing epidemic (1.4521:0.75454).
The numbers of the highest hospitalized symptomatic individuals (NCSII) in Beijing second and first epidemics are 326 (on day 24) and 295 (on day 19), respectively. The corresponding blocking rates (Θ1) to symptomatic infections reach about 97.48% and 89.32%, respectively.
The days (denoted by Day100%) that the blocking rates reached to 100% in Beijing first and second epidemics were on day 98 and on day 28, respectively.
The days (denoted by Dayclear) that all infected individuals were cleared up in Beijing first and second epidemics were on day 140 and on day 56, respectively.
The numbers of cumulative symptomatic infected individuals in Beijing first and second epidemics are 420 individuals and 335 individuals, respectively,
Although the baseline transmission rate β11 in the second Beijing epidemic is much higher than the one in the first Beijing epidemic, the higher blocking rate and the recovery rate (after day 27) to the symptomatic infections make that the the second Beijing epidemic ended 84 days earlier than the first epidemic in Beijing and no died infected patients were reported.
4 Conclusions
The main contributions of this paper are summarized as follows:
It proposes SARDDE models ((1) and (11)) which describe the evolutions of infectious diseases, estimates transmission rates, recovery rates, and blocking rates rate to symptomatic and asymptomatic infections, and symptomatic infected individuals’ death rates.
It provides the criterion inequalities for the asymptotical stability of the disease free equilibrium point of SARDDE (see Theorem 1 and Theorem 3).
It presents the criterion inequalities for epidemic transmission (see Theorem 2 and Theorem 4) of the symptomatic and asymptomatic infections.
It uses two models to simulate the first and second COVID-19 epidemics in Beijing. The simulation results on the end points of transmission intervals were in good agreement with the real word data [17].
The simulation results can provide possible interpretations and estimations of the prevention and control measures, and the effectiveness of the treatments to the two epidemics in Beijing.
Numerical simulation results suggest that
The transmission rate β11 of the symptomatic infections caused by the symptomatic individuals in the second Beijing epidemic is much higher than the one in the first Beijing epidemic.
The blocking rates of 89.32% and 97.48% to the symptomatic infections cannot prevent the spreads of first and second COVID-19 epidemics in Beijing.
The days that the symptomatic infection blocking rates reached to 100% in Beijing first and second epidemics were on day 98 and on day 28, respectively.
Using the data form the beginning to the days after about two weeks from the turning points, we can estimate well or approximately the following outcomes of the first or second COVID-19 epidemics in Beijing.
Numerically interpreted that both symptomatic and asymptomatic individuals cause lesser asymptomatic spread than symptomatic spread.
The lack of the data of the asymptomatic infected individuals in the first COVID-19 in Beijing may prevent the accurate determinations of the number of asymptomatic infected individuals.
The higher blocking rate and the higher recovery rate (after day 27) to the symptomatic infections make that the the second Beijing epidemic ended 84 days earlier than the first epidemic in Beijing.
Keeping the blocking rates, recovery rates and death rates reaching the infection turning points would make the numbers of current hospitalized infected individuals reach, on day 140, 208958 individuals and 25017 individuals for the two Beijing epidemics, respectively.
Different combinations of the parameters of Equations (1) and (11) may generate similar simulation results. Therefore it needs further study to obtain better parameter combinations to interpret COVID-19 epidemics.
The proposed two SARDDEs are simpler than those given in Ref. [11]. However they can better describe and explain the real world data [17].
Because not all infected people can go to the hospital for treatment and be confirmed at the first time. In some cases: adequate resources, no shortage of beds and medical treatment advantages, patients may be left behind when they are discharged from the hospital. Therefore, it does not have very important that the simulation results of the model are required accurately describe every datum reported on the epidemic. Long-term accumulated data, such as the total number of patients and the number of deaths may eliminate short-term deviations. Therefore, the accuracy of predicting long-term epidemics should be the standard for evaluating the rationality of the selected model and unknown model parameters.
When a society cannot withstand high transmission blocking rate load (for example 99%), how can administrative authorities take measures to avoid multiple outbreaks? A recommendation is that the authorities need to at least maintain the prevention and control measures implemented 7 days after reaching the turning point. Administrative authorities need to persuade their society to implement a prevention and control strategy to be withstood until all infection cases are cleared up. It is not a wise strategy to withdraw all prevention and control measures before the number of the all infected people have been cleared up. 100% blocking rate to transmission of COVID-19 infection is the keystone strategy to clear up or reduce spread of an epidemic as early as possible.
The strict prevention and control strategies implemented by Beijing government is not only effective but also necessary. It is expected that the research can provide better understanding, explanation, and dominating the spreads and control measures of epidemics.
Data Availability
All data cite in this article selected from the web ref. [18] (http://wjw.beijing.gov.cn/
Funding
The author has not declared a specific grant for this research from any funding agency in the public, commercial or not for profit sectors.
Conflict of Interest
The author declares no potential conflict of interest.
Data availability statement
Data are available on reasonable request. Please email the author (13501029489@163.com).
Ethical Statement
Not applicable/No human participants included.
Note The previous versions of the manuscript have been published online [11, 18].
Footnotes
E-mail: 13501029489{at}163.com
Dear The medRxiv team, My previous submitted manuscript is another amended medRxiv paper (ID#: MEDRXIV/2022/2732250). Now the changes in the amended version are given as follows: (1) Three new references are added. (2) A few printed/grammar errors have been corrected.. (3) The styles of sevral references are improved. Lequan
↵1 In the cases that some reported data crossed one day, we assign approximately the numbers according to the ratios of time intervals.
