Abstract
Modeling and simulation are important tools that can be used to control, prevent and understand an epidemic spread. This paper 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 simulations, it determines the parameters of the SARDDE, respectively. Numerical simulations of the SARDDE describe well the outcomes of current symptomatic and asymptomatic individuals, recovered symptomatic and asymptomatic individuals, and died individuals, respectively. The numerical simulations suggest that both symptomatic and asymptomatic individuals cause lesser asymptomatic spread than symptomatic spread; the blocking rates of about 90% and 97% to the symptomatic individuals cannot prevent the spread of the first and second COVID19 epidemics in Beijing, respectively. Virtual simulations suggest that the strict prevention and control strategies implemented by Beijing government are effective and necessary. The numerical simulations suggest also that using the data from the beginning to the days after about two weeks from the turning points can estimate approximately the following outcomes of the two COVID-19 academics, respectively. A recommendation to avoid multiple epidemic outbreaks is proposed. 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. Globally by 8 September 2021 there have been about 22.1 million confirmed cases of COVID-19, including 4.5 million deaths, reported to WHO (see https://covid19.who.int/). 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 [1].
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 [2–11]). 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 [12].
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 [13]. 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 [13].
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 [13]. 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.
The rest of this article is organized as follows. Sections 2.1 and 2.4 establish the SARDDEs with 5 and 4 variables. Sections 2.2 and 2.5 provide the criterions of the asymptotical stability of the disease-free equilibriums of the two SARDDEs. Sections 2.3 and 2.6 determine the necessary conditions of disease spreading. Sections 3.1 and 3.4 implement the dynamic simulations of the SARDDEs to describe the reported data of the first and second COVID-19 epidemic in Beijing, respectively. Sections 3.2 and 3.5 discuss and analyze the simulation results. Sections 3.3 and 3.6 implement two virtual simulation examples to predict the outcomes of the two Beijing epidemics on different prevention and control measures. 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 Web [13]. 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 ith time interval, the SARDDE has the form:
where (1 − θ1(i)) and (1 − θ1(i)) (i = 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 ith 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 [13]. 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 [13]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.
(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).
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)).
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 8 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 8-time intervals. On day 86, 15 April, there were 3 Chaoyang district infected people coming back Beijing form foreign country which made calculated blocking rate to rise. Therefore the i in SARDDE model (1) should be chosen as i = 1, 2, …, 8.
First it needs to determine the parameters κ(i), κa(i) and α(i). There are different methods for calculating the recovery rate κ(i) in a specific time interval. It seems to be reasonable that we represent κ(i) via using the number of the recovered patients to divide the days of patients stayed in the hospital during the ith time interval.
Denote s1(i) and s2(i) to be the days that the old patients and the new patients stayed in the hospital during ith time interval. Denote R(i) and d(i) to be the numbers of the recovered patients and died patients during ith time interval. Then we have
Denote ti and ti+1 to be the first date and the last date in the ith time interval. Denote
I1(ti) = the number of the current symptomatic infected individuals on day ti,
I2(ti) = the number of the daily recovered symptomatic infected individuals,
I3(ti) = the number of the daily died infected individuals,
I4(ti) = the number of the daily increased symptomatic infected individuals. Then we obtain
Since there is no information on recovered asymptomatic infected individuals, we take
That is, an asymptomatic infected individual will recover average seven days. The calculated κ(i)′s and α(i)′s are shown in the 7 ∼ 8 columns in Table 1.
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 [13] of the current symptomatic infected individuals and the current died individuals, respectively; NRSII the number [13] of the recovered symptomatic infected individuals.
Second it needs to determine the parameters 
in SARDDE (1). One can assume that S = 1 because the effects of S can be deleted by calculated 
. This makes the calculated 
have general sense. Using the real world data of the first COVID-19 epidemic in Beijing [13] (also see the second row in Table 1) selects following initial condition
Note The infection rate 
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
Substitute parameters κ(1), α(1), θ1(1) and θ2(1) given in Table 1 into SARDDE (1). Using a minimization error square criterion
determines 
.
A group (β11, β12, β21, β22) that makes δ be “smallest” (considering continued simulations) are
The first 4 days’ simulations of SARDDE (1) with the above equation parameters are shown in Figs. 4(a) and 4(b). The simulation results are in good agreement with the reported clinical data (see the first solid line and the first dash line, and the legends in Figs. 4(a) and 4(b)).
Third it needs to determine: θ1(i), θ2(i), i = 2, 3, …, 8. Denote
Let Ic(ti) to be the number of the reported current symptomatic infected individuals at ti in the first CONVID-19 epidemic in Beijing. Let Icr(ti) and Dc(ti) to be the numbers of the reported cumulative recovered infected and the died individuals at ti, respectively.
Using the minimization error square criterion:
determines the θ1(i)′s and θ2(i)′s. The calculated results are shown in Table 1. The corresponding simulation results of SARDDE (1) are shown in Figs. 4(a) and 4(b). Observe that the simulation results of SARDDE (1) describe well the dynamics of the first COVID-19 epidemic in Beijing (see the solid lines and the dash lines and the legends in Figs. 4(a) and 4(b)).
3.1.2 Discussions
On day 0, day 3, and day 140, the numbers of the reported and simulated current symptomatic individuals are are approximate the same. On day 12, day 36, day 73 and day 97, there are only one or two differences. On day 24 and day 52, there are 3 (relative error 1.02%) and 5 (relative error 6.85%) differences.
On day 0, day 3, day 12 and day 140, the numbers of the reported and the simulated current recovered symptomatic individuals are approximate the same. On day 52 and day 97, there are 2 and 3 differences, respectively. On day 24, day 36 and day 73, there are 7, 5, and 6 differences (relative error -10.3%, -2.33%, +1.51%), respectively.
On day 0, day 3, day 12, day 24, day 36, day 52, day 73, day 97 and day 140, 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 [13]. Our simulation results shows that on day 73, the number of the simulated current symptomatic infected individuals was less than one (≈ 0.7), which seems to explain the actual report data.
Computed results (see (29)) 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 (29)) 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 − θ1, 1 − θ2) reach about (98.09%, 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.09%, 99%), the spreading of COVID-19 epidemic can be blocked.
3.1.3 Virtual Simulations
Now assume that after day 24, 12 February, it still keeps the blocking rates (1−θ1(3), 1−θ2(3)) ≈ (89.65%, 99%), the cure 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 1.899e5 and 4679, respectively. The numbers of cumulative recovered symptomatic and died individuals reach about 74278 and 2359, 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 − θ1(5), 1 − θ2(5)) ≈ (98.3%, 99%), the cure 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 both less than one, respectively; the numbers of cumulative recovered symptomatic and died individuals are about 411 and 11, respectively. The results suggest that using the data before day 52 (about 17 days after the turning point) can approximately estimate 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 realword COVID-19 epidemic data from 11 June to 6 August 2020 in Beijing [13]. 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 [13].
(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 7 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 7-time intervals. Therefore the i in SARDDE (11) satisfies i = 1, 2, …, 7.
First it needs to determine the parameters κ(i), κa(i). Denote s1(i) and s2(i) to be the days that the old patients and the new patients stayed in the hospital during ith time interval. sa1(i) and sa2(i) to be the days that the old asymptomatic individuals and the new asymptomatic individuals stayed in the hospital during ith time interval. Denote R(i) and Ra(i) to be the numbers of the recovered symptomatic patients and asymptomatic patients over the ith time interval, respectively. Then κ(i) and κa(i) can be defined by
Denote ti and ti+1 to be the first date and the last date in the ith time interval. Denote
I1(ti) = the number of the current symptomatic infected individuals on day ti,
I2(ti) = the number of the daily recovered symptomatic infected individuals,
I3(ti) = the number of the daily increased symptomatic infected individuals.
Then we obtain
Similarly, we can define Ra(i), sa1(i), sa2(i). The calculated κ(i)′s and κa(i)′s are shown in Table 3.
The data of the second wave COVID-19 epidemics on 7 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.
Second it needs to determine the parameters 
in SARDDE (11). One can assume that S = 1 because the effects of S can be deleted by calculated 
. This makes the calculated 
have general sense. Similar to (36), we have
Using the practical data of the second COVUD-19 epidemic (also see the second line in Table 3) selects following initial condition
Substitute parameters κ(1), κa(1), θ1(1) and θ2(1) showed in Table 3 into SARDDE (11). Using a minimization error square criterion
determines 
. A group 
that makes δ be “smallest” are
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 first 4 days’ clinical data (see the solid and dash lines in Figs. 9(a) and 9(b)).
Third it needs to determine: θ1(i), θ2(i), i = 2, 3, …, 7. Denote
Let Ic(ti) and Ica(ti) to be the numbers of the current symptomatic and asymptomatic infected individuals at time ti, respectively; Icr(ti) and Icra(ti) to be the numbers of the cumulative recovered symptomatic and asymptomatic infected individuals at time ti, respectively.
Using the minimization error square criterion:
determines the θ1(i) and θ2(i). 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.
3.2.2 Discussions
On day 0, day 3, day 10, day 19, and day 44, the numbers of the reported and simulated current symptomatic individuals are approximate the same. On day 27 and day 56, they have only one difference. On day 34, they have 9 differences (relative error -5.3%).
On day 0, day 3, day 10, day 19, day 27, day 44 and day 56, the numbers of the practical and the simulated current asymptotic individuals are approximate the same. On day 34, it has only two differences.
On 0 day, 3, 10 and 19, the numbers of the practical and the simulated cumulative recovered symptomatic individuals are approximate the same. On day 56 and day 44, they has two and three differences, respectively. On day 27, it has 4 differences (relative error +6.7%). On day 34, it has 12 differences (relative error +7.2%).
On day 0, day 3, day 19, and day 27, the numbers of the practical and the simulated cumulative recovered asymptomatic individuals are approximate the same. On day 10 and day 34, they have only one difference. On day 44 and day 56, they has two and three differences, respectively.
Computed results (see (39)) 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 (39)) 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) at 7-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 − θ1, 1 − θ2) reach about (99.71%, 98.19%), the disease-free equilibrium becomes globally asymptotical stable.
3.2.3 Virtual Simulation
Now assume that it keeps still the blocking rates (1 − θ1(3), 1 − θ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 1236 and 72, respectively; The numbers of the cumulative recovered symptomatic and asymptomatic infected individuals reach 23 and 72, 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 − θ1(5), 1 − θ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 29 and 3, respectively. The numbers of the cumulative recovered symptomatic and the asymptomatic individuals are about 309 and 46, respectively. The results suggest that using the data before day 34 (about two weeks after the turning point) can approximately to estimate the following outcome of the second COVID-19 academic 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).
4 Conclusions
The main contributions of this paper are summarized as follows:
Proposed the SARDDE models ((1) and (11)) which describe epidemic transmissions and controls.
Provided the criterion inequalities for the asymptotical stability of the disease free equilibrium point of SARDDE (see Theorem 1 and Theorem 3).
Presented the criterion inequalities for epidemic transmission (see Theorem 2 and Theorem 4) of the symptomatic and asymptomatic infections.
Determined model parameters interpreted quantitatively the outcomes of the first and second epidemics in Beijing.
Numerical simulation results suggest that
100% block rates were implemented after infection appeared 73 days and 34 days for the first and second epidemics, respectively.
Block rates of 90% and 97.5% to symptomatic infections cannot prevent the spreads of the first and second epidemics, respectively (see Fig.5 and Fig.10).
The transmission rate of symptomatic individuals of the first epidemic was about two time larger than that of the second epidemic (1.4492:0.75454 (see (29) and (39)).
Using the data form the beginning to the day after about two weeks from the turning points, we can estimate well or approximately the following outcomes of the first or second COVID-19 academics in Beijing (see Fig.6 and Fig.11).
Numerically interpreted that both symptomatic and asymptomatic individuals cause lesser asymptomatic spread than symptomatic spread (see (29) and (39)).
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.
On day 24–the turning point day of the first epidemics in Beijing, the calculated blocking rate of the symptomatic infections was about 89.65%. On day 19–the turning point day of the second epidemics in Beijing, the calculated blocking rate of the symptomatic infections was about 97.47%. The calculated blocking rates provide a quantitative interpretation that the clearing up days of the hospitalized individuals of the first and second epidemics in Beijing were about 87 days and 25 days, respectively.
The selections of the transmission rates
are difficult because different combinations of can produce very closed simulation errors. Agreements of followed simulations are used to judge the reasonableness of the selected .The proposed two SARDDEs are simpler than those given in Ref. [10]. However they can better describe and explain the real world data [13].
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 practicality, that the simulation results of the model are required accurately describe every datum reported on the epidemic. Longterm 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 100% of the transmission blocking load, 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.
Funding
The author has not declared a specific grant for this research.
Conflict of Interest
The author declares no potential conflict of interest.
Ethical Statement
Not applicable/No human participants included.
Footnotes
E-mail: 13501029489{at}163.com
Dear Editor, The main changes of the amended version are listed as follows: Add a recommendation paragraphs and a comparassiom paragraphs in Conclusion Section. Please see item (vii) and the last paragraph in page 20. Best wishes! Lequan
↵1 In the cases that some reported data crossed one day, we assign approximately the numbers according to the ratios of time intervals.
References
Subject Area
Reviews and Context
0
Comment
0
TRIP Peer Reviews
0
Community Reviews
0
Automated Services
0
Blogs/Media
Author Videos