A fractional-order SEIHDR model for COVID-19 with inter-city networked coupling effects.

In the end of 2019, a new type of coronavirus first appeared in Wuhan. Through the real-data of COVID-19 from January 23 to March 18, 2020, this paper proposes a fractional SEIHDR model based on the coupling effect of inter-city networks. At the same time, the proposed model considers the mortality rates (exposure, infection and hospitalization) and the infectivity of individuals during the incubation period. By applying the least squares method and prediction-correction method, the proposed system is fitted and predicted based on the real-data from January 23 to March 18 - m where m represents predict days. Compared with the integer system, the non-network fractional model has been verified and can better fit the data of Beijing, Shanghai, Wuhan and Huanggang. Compared with the no-network case, results show that the proposed system with inter-city network may not be able to better describe the spread of disease in China due to the lock and isolation measures, but this may have a significant impact on countries that has no closure measures. Meanwhile, the proposed model is more suitable for the data of Japan, the USA from January 22 and February 1 to April 16 and Italy from February 24 to March 31. Then, the proposed fractional model can also predict the peak of diagnosis. Furthermore, the existence, uniqueness and boundedness of a nonnegative solution are considered in the proposed system. Afterward, the disease-free equilibrium point is locally asymptotically stable when the basic reproduction number R 0 ≤ 1 , which provide a theoretical basis for the future control of COVID-19. © Springer Nature B.V. 2020.

Introduction

A infectious disease caused by a new coronavirus was discovered firstly in Wuhan, China, in the end of 2019. Contrary to the initial report [1], it is indeed spread from person to person through frequent interpersonal communication [2]. Then, it quickly spread all over the word. By 5 February, 2020, more than 24,550 had been confirmed. Since January 23, 2020, the Wuhan government has gradually implemented city-wide quarantine, which has greatly suppressed the spread of the virus.

As the main transportation hub in central China, Wuhan has a large population movement during the Spring Festival. During the transmission of the new coronavirus, population migration caused by geographical factors has added fundamental difficulties for the prevention and intervention of this epidemic, such as the domestic passenger flow from Wuhan is estimated at 5 million during the Spring Festival, which makes it difficult for following the infected people. Prasse et al. studied a discrete epidemic model with network to describe this infectious disease and their work pointed out that the epidemic model with network helped to accurately predict the outbreak of epidemics [3]. Peng et al. found that the first appearance of COVID-19 could be dated back to the beginning of December 2019 in China by a generalized SEIR model [4]. Furthermore, the basic reproduction number represents number of cases in which an infected individual is expected to have a second infection in a fully susceptible population. If , an infected person, on average, infects fewer than one infection in its infectious period, and the epidemic will not exist. Conversely, if , each infected individual infects more than one infected individual, and the disease will affect the entire population [5]. At present, there are a large number of articles that estimate the basic reproduction number to further determine whether COVID-19 will widespread in the population [6-8]. Meanwhile, during the transmission of infectious diseases, susceptible people come into contact with infected people and are infected with a probability. A great deal of evidence shows that incidence rates are an important method for describing infectious diseases [9-11]. Bilinear incidence rate is often used in the modeling process of COVID-19 because of its high infectivity. The incubation period of the infected person in the traditional SEIR model is not infectious, but there is already evidence that for COVID-19, the incubation period is very high infectious [12]. Not only that, due to various conditions, all patients infected with COVID-19 could not be sent to the hospital, and infected individuals have a high mortality rate before hospitalization.

It is worth noting that in the classic integer-order epidemic model, the state of the model does not depend on the history. However, not only does the infectious diseases depends on its current state, but also on its past state in real life. At the same time, the fractional order system can describe the change of the system whose instantaneous rate of change depended on the past state. This way of description is called the memory effect [13]. And it is confirmed that the fractional model perfectly fits the test data of memory phenomena in different disciplines by using numerical least square method. That is, a physical meaning of the fractional-order is an index of memory [14]. A large number of articles have discussed the stability of fractional order systems in different ways [15-17], such as Jan et al. concerned with some stability properties of linear autonomous fractional differential and difference systems involving derivative operators of the Riemann–Liouville type and they showed discretizations based on backward differences can retain the key qualitative properties of underlying fractional differential systems [15], and Deng et al. derived the sufficient conditions for the local asymptotical stability of nonlinear fractional differential equations, including Riemann–Liouville fractional order and Caputo fractional order derivative [16]. Furthermore, Smethurst et al. found that the waiting time for patients follows the power law model [18]. And the power law distribution generates Caputo fractional-order derivative of the same order [19]. When dealing with practical problems, Caputo fractional order and integer order derivatives have the same initial conditions, which has a specific physical meaning. At the same time, the Caputo derivative of its constant is zero. More importantly, the accuracy of Caputo derivative can supersede the integer derivative caused by its change and non-local behavior. Angstmann et al. derived a fractional-order infectivity SIR model from a stochastic process and they found the fractional derivative appears in the generalized master equations of a continuous time random walk through SIR compartments, with a power-law function in the infectivity [20]. Khan et al. described the interaction between bats and unknown hosts, people to people, and the source of infection (seafood market), and their work shown fractional model played an important role in limiting the number of infected people [21]. Chen et al. established a fractional delay dynamic system (FTDD) to describe COVID-19, and they used reconstruction coefficients to predict the spread of new coronaviruses COVID-19 [22]. Amjad et al. constructed the fractional order COVID-19 model. In their work, the effects of preventive measures, various mitigation measures are estimated, and future outbreaks and potential control strategies are predicted [23]. A fractional SEIQRD model is analyzed by Xu et al., which is of guiding significance for predicting the possible outbreaks of some epidemic diseases [24].

To incorporate the time fractional order and the coupling effect between cities, a SEIHDR epidemic model is established to study the dynamic behavior of COVID-19. Many studies on COVID-19 shown that not only infected individuals are infectious, but people in the incubation period have same infectivity with infected individuals. Therefore, the infectivity of the incubation period is considered here. Further, the mortality of different individuals (exposed, infected and hospitalization) is considered in this model. Furthermore, during the outbreak of COVID-19, due to the limitation of various resources, the patients with infection cannot be diagnosed immediately and the patients with diagnosis cannot be hospitalized immediately, all need waiting time. Thus, from the analysis made above, a fractional inter-city network SEIHDR model for COVID-19 is established. Then, the fractional inter-city network model is verified to be rational by several numerical examples with the official data [25]. Compared with the integer system, the non-network fractional model has been verified to fit the data of Beijing, Shanghai, Wuhan and Huanggang better. The analysis with inter-city network shows that the proposed system may not have a positive effect on the spread of viruses in China, but may be important to other countries. Finally, the stability of the proposed system is analyzed through , which has theoretical significance for further intervention and prevention of this infectious disease.

The rest of the article is organized as follows. COVID-19’s SEIHDR fractional model is established in Sect. 2. Then, some dynamic behaviors of the proposed system are analyzed in Sect. 3. In Sect. 4, numerical simulations are provided to illustrate the theoretical results. Finally, the conclusion is held in Sect. 5.

Model development

It has been found that the fractional-order derivatives have a wide range of applications in the modeling of many asynchronous dynamic processes, such as engineering, biology, medicine and many other fields [26-30]. Before presenting the fractional epidemic model, some necessary preliminaries are introduced.

Preliminaries

This section begins with some definitions and results.

Definition 2.1

[31] The Gamma function satisfies the following equation:

Definition 2.2

[31] For , the Caputo fractional-order derivative of order ( for a function is defined by

Remark 2.1

If ,

Definition 2.3

[32] Consider the Caputo fractional dynamical system:a constant is an equilibrium point of the above system if and only if .

Lemma 2.1

[33] The fractional-order system is considered as followed:with the initial condition , where and . If g(t, x) satisfies the Lipschitz condition on x, the above system has a unique solution.

Lemma 2.2

[34] Consider the following fractional-order system:where , . The equilibrium points are locally asymptotically stable if and only if all eigenvalues of the Jacobian matrix satisfy the following equation:

System description

In December 2019, a new coronavirus was first discovered in Wuhan, China. There are a large number of epidemic models that describe the spread of this infectious disease [35]. And the spread of COVID-19 started during the Chinese New Year, and the massive population movement had made the connection between cities closer. Meanwhile, Wuhan, the most severely affected area, as a major transportation hub in central China, COVID-19 is mainly spread by means of transportation. Therefore, it is necessary to establish an inter-city network to observe the transmission of this infectious disease. Due to the limitation of medicine, latent individuals and infection individuals cannot be hospitalized immediately, which will increase the mortality rate. In addition, Tang et al. considered that the incubator of COVID-19 infection is highly infectious [12]. From the analysis made above, a fractional inter-city network SEIHDR epidemic model is constructed as follows:with the initial conditionThe city k’s total population is presented by that is classified into , , , , and denoted the class of susceptible individuals, exposed individuals, infective individuals (infected but not hospitalized), hospitalization individuals, death individuals and recovered individuals, respectively. The susceptible individual contact with and then infected by , where is the transmission coefficient. The dead individual includes death during susceptible , exposure , infection , and hospitalization , where denotes the natural death rates, and imply the disease-related mortality. The parameter denotes the inflow number of susceptible individuals; be the recovery rate; imply the transit rate of the exposed class ; denote hospitalization rate. Furthermore, , , and are positive constants; bounded function and satisfy and for , where and are positive constants. Without loss of generality, , , , , , and are still denoted as , , , , , and , respectively.

However, during the outbreak of COVID-19, individual migration has been tightly controlled, especially for the city of Wuhan to implement the closure measures, and the natural death rate of individuals is negligible compared to the death caused by COVID-19. By the above assumptions, the improved SEIHDR model is given as followed:Then, the transmission diagram of a fractional SEIHDR model without inter-city network for COVID-19 is shown in Fig. 1.

Fig. 1

Transmission diagram for system (2)

Basic properties about the model

The dynamical analysis of the proposed system (1) is studied in this section. Here, it can be seen that the susceptible individual , the exposed individual , the infected individual and the hospitalized individual of system (1) are not effected by the recovered individual and the death class , so the following system is considered:with the initial condition

Nonnegativity and boundedness

Before the numerical process, the existence, uniqueness and boundedness of the nonnegative solution of system (4) must be proved. Therefore, this subsection discusses the above properties of system (4).

Theorem 3.1

Consider the following initial condition:system (4) has a unique bounded and nonnegative solution , , .

Proof

Let . Adding all equations of system (4) yieldswhich , then one hasSo one has , , and .

Let , and whereObviously, one haswhere , , , , . So satisfies the Lipschitz condition on . Then, system (4) has a unique bounded solution .

Further, the auxiliary system is considered as followed:Through the comparison theorem, it is found thatThus, it can be concluded that , , and for .

Stability analysis

Now, exploring the stability of the system (4) through the basic reproduction number . Obviously, the disease-free equilibrium for system (4) is 0).

Theorem 3.2

The basic reproduction number of system (4) iswhere , , , and the spectral radius .

Let , , and . Then, taking the derivative of , , and with respect to , and at , it can be concluded thatwhere , , , , , , , and . Then, according to [5], it can be yielded that

It is obvious to get the following theorem:

Theorem 3.3

If the basic reproduction number , system (4) is locally stability at the disease-free equilibrium point 0).

Consider the following Jacobian matrix of system (4) at where . One can calculate that the eigenvalues are and andwhere , and are all eigenvalues of the matrix , and , respectively. Then, if , it can be yielded that . Thus, system (4) is locally stability at the disease-free equilibrium point .

Remark 3.1

When city k and city j have no communication (ie, city k is closed), the basic reproduction number of city k can be expressed as followed:So, is local stability when the basic reproduction number .

Remark 3.2

It is evident that are dependent on and . The sensitivity of to the other parameters , , , and is calculated as follows:where , , and represent the normalized sensitivity on , , and , respectively. It is worth noting that the n times’ increase on leads to the n times’ increase on , but the n times’ increase on , , and leads to the n times’ decrease on .

Numerical analysis

It can be seen from the previous discussion that system (1) has a unique bounded solution, and the disease-free equilibrium point is locally asymptotically stable, which can provide theoretical support for the further control and prediction of COVID-19. Then, the numerical solutions of system (1) are analyzed by least squares method [36] and predictor-correctors scheme. In order to predictions, it is possible to move the data by a fixed date m before March 18. Determine predictions are made based on data of from January 23 to March 18, 2020. Then, it can be predict the number of people diagnosed on March 18 and m represents the predicted number of days. Like in [12], the disease-related mortality and the cure rates are time-varying functions aswhere and are initial cure rate and mortality. Considering without the network, it can be seen that the fitting effect of the fractional system (3) on the peak value and the peak time is a little bit better than that of the integer system (ie ) from Figs. 2, 3, 4, 5 and 6. And the relative error between the real-date and the numerical solutions of and can be seen from Table 1. Particularly, corresponding the basic reproduction number and the fractional-order are shown in Table 2. Then, consistent with Remark 3.2, Fig. 7 indicates the sensitivity of . Furthermore, according to Fig. 8, it is observed that the peak time and the peak value are given as the coefficient and the hospitalization rate change. The numerical results show that strict control of individual exposure and increased detection of nucleic acid reagents play a vital role in reducing confirmed cases and delaying the peak period. However, from Fig. 8, it can also be found that the peak time reaches earlier and the peak value becomes larger as the disease coefficient increases. So, the non-closed policy adopted by many countries are also reasonable

Fig. 2

The number of recovered and confirmed for Beijing (m = 6)

Fig. 3

The number of recovered and confirmed for Shanghai ()

Fig. 4

The number of recovered and confirmed for Hubei ()

Fig. 5

The number of recovered and confirmed for Wuhan ()

Fig. 6

The number of recovered and confirmed for Huanggang ()

Table 1

Error with the real-data and the numerical solutions

Error (%)	Different country
	Beijing	Shanghai	Hubei	Wuhan	Huanggang
α=1	23.6	34.3	21.8	19.2	78.3
α≠1	5.9	14.1	16.2	19.1	49.7

Table 2

Performance comparison with different city

Different country	Different index
	α	R_0
Beijing	1.1913	0.8855
Shanghai	0.9844	0.8833
Hubei	1.4343	0.9848
Wuhan	0.9995	1.2675
Huanggang	0.0332	0.9568
America	0.7945	1.1443
Japan	0.3279	1.6346
Italy	0.7056	1.3368×10^3

Fig. 7

The sensitivity of the basic reproduction number

Fig. 8

Prediction of peak in Hubei province

When considered Hubei and other cities interacting with each other, the communication between cities is not clear and can be identified from system (1). Then, Figs. 9 and 10 describe the network coupling between Beijing, Shanghai, Wuhan and Huanggang. It can be seen that since China closed the city on January 23, the epidemic model with network is no longer suitable for China. However, regarding the population movements in other countries, a network-based model can be established to describe the spread of the virus, which will show better performance, and this will also be the focus of future research

Fig. 9

With and without inter-city network in Beijing and Shanghai

Fig. 10

With and without considering inter-city network in Wuhan and Huanggang

Meanwhile, considering system (1) without the natural death rate , the short-term prediction for America and Japan is shown in Figs. 11, 12 when considered the real-date dating from 22 January and 1 February to 16 April, 2020 [37], respectively. And it can be shown from Figs. 11 and 12 that the number of confirmed and recovery individuals will keep increasing for America and Japan. Finally, the prediction is considered in Fig. 13, according to the Italian real-data of hospitalization, death and recovery from February 24 to March 31 [38]. And it is obviously that the peak time of symptomatic individuals and hospitalizations reach on 30 June and 20 April, respectively. Then, the basic reproduction number and the fractional order can be seen from Table 2, which suggests that the epidemic model of Italy, America and Japan will not disappear in a short time.

Fig. 11

American report and forecast from 22 January to 16 April

Fig. 12

Japan report and forecast from 1 February to 16 April

Fig. 13

Italy report and forecast from 24 February to 31 March

Discussion

Based on the coupling effects of inter-city network in Beijing, Shanghai, Wuhan and Huanggang, this paper proposes a fractional SEIHDR model. By applying the method of least squares and prediction-correction, the numerical solution of the system (3) is compared with the actual value of Beijing, Shanghai, Hubei, Wuhan and Huanggang from January 23 to March 18, 2020. It can be found from this study that fractional systems have a better fitting effect than integer systems. Besides, based on the inter-city network, the analysis shows that the system with inter-city network could not be a better description of infectious transmission in China. At the same time, system (1) fits the data of Japan, the USA from January 22 to April 16 and Italy from February 24 to March 31, and predicts the peak value and the peak time, respectively. Moreover, the existence, uniqueness and boundedness of a nonnegative solution are established. Further, the local stability at the disease-free equilibrium point is studied by . And the sensitivity of without network are researched to provide a theoretical basis for infectious disease control.

Although the above study may not be able to show the impact of the inter-city network-based epidemic model on disease transmission in China because of strict locking measures, it may on the whole indicates that the fractional epidemic model with network may predict more accurate than independent. Although the situation of COVID-19 can be effectively predicted in China, this study raises some questions that require further research: whether the network model can be used to predict epidemics in other regions; how medical and other factors affect the virus transmission trend. Further study how to predict and intervene in the development of epidemic situations in other countries when the dynamic model includes network effects.