The Analysis of a SEIRS Epidemic Model with Time Delay on Complex Networks

Abstract

In this paper, a new epidemic SEIRS model with time delay on complex networks is proposed. Based on the mean field theory, the basic reproductive number and equilibriums of the model are derived. Moreover, the impact of the network topology and time delay on the basic reproductive number is analyzed. Theoretical analyses indicate that the basic reproductive number is dependent on the topology of the underlying networks. The time delay cannot change the basic reproductive number, but it can reduce the endemic level and weaken the epidemic spreading. The global asymptotically stability of the disease-free equilibrium and the permanence of epidemic are proved in detail. Numerical simulations confirm the analytical results.

Share and Cite:

Dong, J. , Li, T. , Wan, C. and Liu, X. (2017) The Analysis of a SEIRS Epidemic Model with Time Delay on Complex Networks. Open Access Library Journal, 4, 1-10. doi: 10.4236/oalib.1103901.

1. Introduction

Mathematical models describing the population dynamics of infectious diseases have played an important role in better understanding epidemiological patterns and disease controls for a long time. The basis of modeling dynamics and evolution of infectious diseases are summarized by Anderson and May [1] . It is well known that the spread of a communicable disease involves disease-related factors such as infectious agent, mode of transmission, incubation periods, infectious periods, susceptibility, and resistance.

Continuous time deterministic epidemic models are traditionally elaborated as systems of ordinary differential equations. More realistic models should include some states of these systems, and ideally, a real system should be modeled by delay differential equation. Time delay plays an important part in propagation process of the epidemic, and we can simulate the exposed period of infectious diseases, the infections period of patients and the immunity period of recovery of the disease with time delay. However, compared with studies of the dynamical behaviors of the epidemic models with time delays on complex networks, only a few attentions about early studies have been paid to them on complex networks. In recent years, the dynamics of the SIS, SIR, SEI and SEIR epidemic models have received considerable attention [2] [3] [4] [5] [6] . In order to study the effects of disease latency or immunity in real life, the delay is incorporated in such models [7] - [12] . Liu and Deng et al. discussed epidemic SIS model with discrete time delay which represents the infectious period [13] , and they obtained the basic reproduction number and discussed the persistence of the disease, but they failed to give the proof in detail. Wang and Wang et al. discussed an epidemic SIR model with discrete time delay which represented the exposed period [14] . Based on the above, in this paper, we will present a suitable epidemic SEIRS model with time delay which represents the infectious period on complex networks using functional differential equations to investigate the epidemic spreading.

The rest of the paper is organized as follows: Section 2 derives a SEIRS model with time delay mechanism on scale-free networks. Then, Section 3 obtains two equilibriums and basic reproductive number. In Section 4, numerical simulations are performed. Finally, we conclude the paper in Section 5.

2. Model Formulation

One of the most effective interventions to contain the spread of epidemic diseases is the delay mechanism as discussed above. In order to investigate the efficiency of delay mechanism, we consider the new SEIRS model with delay mechanism on complex networks. On the complex networks, each individual is represented by a node of the network and the edges are the connections between individuals along which the infection may spread. Taking the connectivity among different nodes into consideration, let S k ( t ) , E k ( t ) , I k ( t ) and R k ( t ) be the relative densities of susceptible, exposed, infected and recovered nodes of degree k at time t respectively. In the course of disease transmission, a susceptible individual will be infected with probability β 1 ( β 2 ) if it connects to a exposed (infected) one. r is the state transition rate from exposed to infectious individual. The rate constant of recovery for infected individuals is denoted by α . τ , for recovered nodes, is the effectively constant of immunity, and there are still recovered nodes, when t moment from t τ moment, which will lose the immunity ability to become the susceptible nodes again. Here, we assume that the birth rate equals the death rate, and the rate constant is l . Thus, the dynamic mean-field reaction rate equations can be written as

{ d S k ( t ) d t = l l S k ( t ) k β 1 Θ 1 ( t ) S k ( t ) k β 2 Θ 2 ( t ) S k ( t ) + α I k ( t τ ) e l τ d E k ( t ) d t = k β 1 Θ 1 ( t ) S k ( t ) + k β 2 Θ 2 ( t ) S k ( t ) ( l + r ) E k ( t ) d I k ( t ) d t = r E k ( t ) ( l + α ) I k ( t ) d R k ( t ) d t = α I k ( t ) α I k ( t τ ) e l τ l R k ( t ) (1)

where the probability Θ 1 ( t ) denotes a link pointing to an exposed individual, and satisfies

Θ 1 ( t ) = k k P ( k ) E k ( t ) s s P ( s ) = 1 k k k P ( k ) E k ( t ) .

The probability Θ 2 ( t ) denotes a link pointing to an infected individual, and satisfies

Θ 2 ( t ) = k k P ( k ) I k ( t ) s s P ( s ) = 1 k k k P ( k ) I k ( t ) .

Where k describes the average degree and P ( k ) describes the connectivity distribution. E ( t ) = k P ( k ) E k ( t ) is the total density of exposed individuals, and the I ( t ) = k P ( k ) I k ( t ) is the total density of infected individuals in the whole network.

For simplicity, through this paper, let ρ ( t ) = β 1 Θ 1 ( t ) + β 2 Θ 2 ( t ) and P k = k P ( k ) / k . Thus, the system (1) is equivalent to the following model

{ d S k ( t ) d t = l l S k ( t ) k ρ ( t ) S k ( t ) + α I k ( t τ ) e l τ d E k ( t ) d t = k ρ ( t ) S k ( t ) ( l + r ) E k ( t ) d I k ( t ) d t = r E k ( t ) ( l + α ) I k ( t ) d R k ( t ) d t = α I k ( t ) α I k ( t τ ) e l τ l R k ( t ) (2)

With the normalization condition, there is

S k ( t ) + E k ( t ) + I k ( t ) + R k ( t ) = 1.

Obviously, the initial conditions for system satisfy

0 S k ( 0 ) , E k ( 0 ) , I k ( 0 ) , R k ( 0 ) 1.

3. Analysis of the Novel SEIRS Model

In this part, we put forward the analytic solution of system (2) for describing the dynamic behavior.

Theorem 1. Denote

R 0 = [ β 1 ( l + α ) + r β 2 ] k 2 ( l + α ) ( l + r ) k .

The disease-free equilibrium E 0 ( 1 , 0 , 0 , 0 ) of system (2) always exists and there exists a positive endemic equilibrium E + ( S k , E k , I k , R k ) of system (2) when R 0 > 1 .

Proof. To get the equilibrium solution E + ( S k , E k , I k , R k ) , we need to make the right side of system (2) equal to zero. Then, the equilibrium E + ( S k , E k , I k , R k ) should satisfy

{ l + α I k e l τ k ρ S k l S k = 0 k ρ S k ( l + r ) E k = 0 r E k ( l + α ) I k = 0 α I k α I k e l τ l R k = 0. (3)

Thus, we have

{ S k = l + α I k e l τ l + k ρ E k = l + α r I k R k = α α e l τ l I k (4)

According to the normalization condition S k + E k + I k + R k = 1 , we can get

{ S k = ( r l 2 + k ρ r l ) α e l τ + l ( l + k ρ ) [ l ( l + α ) + α r ( 1 e l τ ) + r l ] B k E k = ( l + α ) k ρ l B k I k = k ρ r l B k R k = ( 1 e l τ ) α k ρ r B k (5)

where B k = r l α e l τ + ( l + k ρ ) [ l ( l + α ) + α r ( 1 e l τ ) + r l ] .

For ρ = β 1 Θ 1 ( t ) + β 2 Θ 2 ( t ) = k 1 k k P ( k ) ( E k β 1 + I k β 2 ) β 1 + β 2 , we can find that ρ = 0 satisfies Equation (5). Therefore, S k = 1 and E k = I k = R k = 0 is a disease-free equilibrium of system (2). Substituting I k and E k of Equation (5) into ρ , we can obtain

ρ = k k 2 P ( k ) k l ρ [ r β 2 + ( l + α ) β 1 ] B k . (6)

We suppose

f ( ρ ) = 1 k k 2 P ( k ) B k l [ r β 2 + ( l + α ) β 1 ] k ,

so that, Equation (6) can be written

ρ f ( ρ ) = 0.

Due to d f ( ρ ) d ρ > 0 and f ( β 1 + β 2 ) > 0 , the equation f ( ρ ) = 0 has a non-trivial solution if and only if f ( 0 ) < 0 .

So, we have

β 1 ( l + α ) + r β 2 ( l + α ) ( l + r ) k k 2 P ( k ) k > 1 .

Let’s define k 2 = k k 2 P ( k ) , thus, we can compute the base reproduction number as follows

R 0 = [ β 1 ( l + α ) + r β 2 ] k 2 ( l + α ) ( l + r ) k . (7)

Namely, there is the unique nontrivial solution if and only if R 0 > 1 . From Equations ((4) and (5)), we are able to obtain

0 < S k < 1 , 0 < E k < 1 , 0 < I k < 1 , 0 < R k < 1 .

So, we define the equilibrium E + ( S k , E k , I k , R k ) . Hence, when R 0 > 1 , system (2) has a unique positive endemic equilibrium E + ( S k , E k , I k , R k ) . The proof is completed.

Remark. The basic reproductive number is obtained by (7), which depends on different model parameters and fluctuations of the degree distribution. The delay parameter cannot change the basic reproductive number. Clearly, in the eco-limits of infinite size network the number of nodes grows to infinity, i.e.,

N , then k 2 k grows to infinity, so the increase of the basic reproductive number, i.e., R 0 , is obvious.

Theorem 2. If R 0 < 1 , the disease-free equilibrium E 0 is globally asymptotically stable. When R 0 > 1 , the epidemic disease is permanent, which means there exists ξ > 0 , such that

lim inf t { S k ( t ) , E k ( t ) , I k ( t ) , R k ( t ) } k = 1 n ξ ,

where ( S k ( t ) , E k ( t ) , I k ( t ) ) is any solution of system (2), satisfying E k ( 0 ) > 0 or I k ( 0 ) > 0 .

Proof. For S k ( t ) + E k ( t ) + I k ( t ) + R k ( t ) = 1 , i.e., the values of S k ( t ) , E k ( t ) and I k ( t ) are fixed, a unique corresponding R k ( t ) exists. We can discuss the first three equations of system (2) with R k ( t ) = 1 S k ( t ) E k ( t ) I k ( t ) at all t . System (2) satisfies

{ d E k ( t ) d t = k ρ ( 1 E k ( t ) I k ( t ) R k ( t ) ) ( l + r ) E k ( t ) d I k ( t ) d t = r E k ( t ) ( l + α ) I k ( t ) d R k ( t ) d t = α I k ( t ) l R k ( t ) α I k ( t τ ) e l τ ,

Thus, the jacobian matrix of the virus-free equilibrium for the system (2) can be obtained, the jacobian matrix is a 3 n × 3 n matrix, such that

J = ( A 11 A 1 n A n 1 A n n ) ,

A 11 = ( ( l + r ) + k β 1 P 1 k β 2 P 1 0 r ( l + α ) 0 0 α ( 1 e l τ ) l ) ,

A 1 n = ( k β 1 P n k β 2 P n 0 0 0 0 0 0 0 ) ,

A n 1 = ( n k β 1 P 1 n k β 2 P 1 0 0 0 0 0 0 0 ) ,

A n n = ( ( l + r ) + k β 1 P n n β 2 P n 0 r ( l + α ) 0 0 α ( 1 e l τ ) l ) .

We can put forward the characteristic equation of the disease-free equilibrium from jacobian matrix, such that

( λ + l + α ) n 1 ( λ + l + r ) n 1 ( λ 2 + p λ + q ) = 0 ,

where p = ( l + r ) + ( l + α ) n β 1 P n and q = ( l + α ) ( l + r ) [ ( l + α ) β 1 + r β 2 ] n P n .

Note that R 0 < 1 is equivalent to q > 0 which

( l + r ) ( l + α ) [ ( l + α ) β 1 + r β 2 ] i = 1 n i P i > 0 , so that

( l + r ) ( l + α ) > [ ( l + α ) β 1 + r β 2 ] i = 1 n i P i , which means that p > 0 . Therefore, the

eigenvalue λ of J are all negative if R 0 < 1 ; otherwise, if and only if R 0 < 1 , here exists a unique positive eigenvalue λ of J . Using the Perron-Feobenius theorem, this suggests that the maximal real part of all eigenvalues of J is positive only if R 0 > 1 . Thus, we can obtain the results of this theorem by using a theorem presented of Lajmanovich and York [15] . The proof is completed.

4. Numerical Simulations and Sensitivity Analysis

To support and explain our theoretical analysis results, we present several numerical simulations in this section. The system (2) is considered on complex networks, whose the degree distribution is P ( k ) = c k 3 , where the parameter satisfies k = 1 n c k 3 = 1 , n = 1000 . We suppose S k ( 0 ) = 0.9 , E k ( 0 ) = 0.1 , I k ( 0 ) = 0 and R k ( 0 ) = 0 .

In Figure 1, the parameters are chosen as β 1 = 0.01 , β 2 = 0.1 , l = 0.2 , u = 0.4 , r = 0.5 , thus the threshold value R 0 = 0.4230 < 1 . According to R 0 < 1 , we can see that when R 0 < 1 , ( E k + I k ) grows to zero, i.e., the infectious individuals will ultimately disappear. If R 0 < 1 , the equilibrium E 0 is globally attractive and the disease eventually disappear.

In Figure 2, the parameters are chosen as β 1 = 0.3 , β 2 = 0.4 , l = 0.2 , u = 0.4 , r = 0.5 , thus the threshold value R 0 = 2.8705 > 1 . We can see that when R 0 > 1 , ( E k + I k ) grows to a constant, i.e., the epidemic disease is permanent and the number of infected individuals will converge to a positive constant.

In Figure 3, the parameters are chosen as β 1 = 0.01 , β 2 = 0.1 , l = 0.2 , u = 0.4 , r = 0.5 , R 0 = 0.4230 . We can see that the corresponding ( E k + I k )

Figure 1. With R 0 = 0.4230 < 1 , the prevalence of E k ( t ) and I k ( t ) versus t corresponding for k = 50 , 100 , 200 , 500.

Figure 2. With R 0 = 2.8705 > 1 , the prevalence of E k ( t ) and I k ( t ) versus t corresponding for k = 50 , 100 , 200 , 500.

Figure 3. With R 0 < 1 , the prevalence of E k ( t ) and I k ( t ) versus t corresponding for k = 100 and different parameter τ .

Figure 4. With R 0 > 1 , the prevalence of E k ( t ) and I k ( t ) versus t corresponding for k = 100 and different parameter τ .

decreases significantly as the delay parameter τ increases, i.e., a larger delay parameter can weaken the spreading of disease. The larger delay is, the lower corresponding ( E k + I k ) level is.

In Figure 4, the parameters are chosen as β 1 = 0.3 , β 2 = 0.4 , l = 0.2 , u = 0.4 , r = 0.5 , R 0 = 0.4230 . We can see that the corresponding ( E k + I k ) decreases significantly as the delay parameter τ increases, i.e., a larger delay parameter can reduce the endemic level.

5. Conclusion

In this paper, a novel SEIRS epidemic model of delay mechanism on complex networks has been presented. Using the field theory, we calculated two epidemiologically relevant quantities, i.e., the basic reproductive number for invasion and the endemic prevalence of epidemic diseases. We have proved that if the basic reproductive number R 0 < 1 , the disease-free equilibrium is globally asymptotically stable; while R 0 > 1 , the endemic equilibrium is permanent. Moreover, increasing delay parameters can result in the weakness of the diseases spreading and the decrease of population infected. Numerical simulations show that the endemic equilibrium E + is globally asymptotically stable when R 0 > 1 (as shown in Figure 2). It is interesting but challenging to discuss the stability of equilibrium E + . We leave it for our future work. This study has valuable guiding significance in effectively predicting epidemic spreading.

Acknowledgements

This work is supported in part by the National Natural Science Foundation of China under Grants 61672112, the Teaching Research Project of Hubei Province under Grant 2014260, and the Teaching Research Project of Yangtze University under Grant JY2015010.

Conflicts of Interest

The authors declare no conflicts of interest.

References

[1] Anderson, R.M. and May, R.M. (1991) Infectious Diseases of Humans: Dynamics and Control. Oxford University Press Inc., New York.
[2] Wei, J. and Zou, X. (2006) Bifurcation Analysis of a Population Model and the Resulting SIS Epidemic Model with Delay. Journal of Computational and Applied Mathematics, 197, 169-187.
https://doi.org/10.1016/j.cam.2005.10.037
[3] Huang, G., Takeuchi, Y., Ma, W., et al. (2010) Global Stability for Delay SIR and SEIR Epidemic Models with Nonlinear Incidence Rate. Bulletin of Mathematical Biology, 72, 1192-1207.
https://doi.org/10.1007/s11538-009-9487-6
[4] Xiao, Y., Chen, L. and ven den Bosch, F. (2002) Dynamical Behavior for a Stage-Structured SIR Infectious Disease Model. Nonlinear Analysis: Real World Applications, 3, 175-190.
https://doi.org/10.1016/S1468-1218(01)00021-9
[5] Gao, S., Chen, L. and Teng, Z. (2007) Impulsive Vaccination of an SEIRS Model with Time Delay and Varying Total Population Size. Bulletin of Mathematical Biology, 69, 731-745.
https://doi.org/10.1007/s11538-006-9149-x
[6] Li, T., Liu, X.D., Wu, J., et al. (2016) An Epidemic Spreading Model on Adaptive Scale-Free Networks with Feedback Mechanism. Physica A: Statistical Mechanics and Its Applications, 450, 649-656.
https://doi.org/10.1016/j.physa.2016.01.045
[7] Li, J. and Ma, Z. (2006) Analysis of Two SEIS Epidemic Models with Fixed Period of Latency. Journal of Systems Science and Mathematical Sciences, 26, 228-236.
[8] Enatsu, Y., Nakata, Y. and Muroya, Y. (2012) Global Stability of SIRS Epidemic Models with a Class of Nonlinear Incidence Rates and Distributed Delays. Acta Mathematica Scientia, 32, 851-865.
https://doi.org/10.1016/S0252-9602(12)60066-6
[9] Nakata, Y. (2011) On the Global Stability of a Delayed Epidemic Model with Transport-Related Infection. Nonlinear Analysis: Real World Applications, 12, 3028-3034.
https://doi.org/10.1016/j.nonrwa.2011.05.004
[10] Zhang, T. and Teng, Z. (2008) Global Behavior and Permanence of SIRS Epidemic mOdel with Time Delay. Nonlinear Analysis: Real World Applications, 9, 1409-1424. https://doi.org/10.1016/j.nonrwa.2007.03.010
[11] Liu, Q., Jiang, D., Shi, N., et al. (2016) Asymptotic Behavior of a Stochastic Delayed SEIR Epidemic Model with Nonlinear Incidence. Physica A: Statistical Mechanics and Its Applications, 462, 870-882.
[12] Sharma, N. and Gupta, A.K. (2017) Impact of Time Delay on the Dynamics of SEIR Epidemic Model using Cellular Automata. Physica A: Statistical Mechanics and Its Applications, 471, 114-125.
[13] Liu, Q.M., Deng, C.S. and Sun, M.C. (2014) The Analysis of an Epidemic Model with Time Delay on Scale-Free Networks. Physica A: Statistical Mechanics and Its Applications, 410, 79-87.
[14] Wang, J., Wang, J., Liu, M., et al. (2014) Global Stability Analysis of an SIR Epidemic Model with Demographics and Time Delay on Networks. Physica A: Statistical Mechanics and Its Applications, 410, 268-275.
[15] Lajmanovich, A. and Yorke, J.A. (1976) A Deterministic Model for Gonorrhea in a Nonhomogeneous Population. Mathematical Biosciences, 28, 221-236.

Copyright © 2024 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.