留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

两个异质城市间具有路途感染的SEIR传染病模型

上一篇

下一篇

康萍, 刘贤宁. 两个异质城市间具有路途感染的SEIR传染病模型[J]. 西南师范大学学报(自然科学版), 2020, 45(5): 7-12. doi: 10.13718/j.cnki.xsxb.2020.05.002
引用本文: 康萍, 刘贤宁. 两个异质城市间具有路途感染的SEIR传染病模型[J]. 西南师范大学学报(自然科学版), 2020, 45(5): 7-12. doi: 10.13718/j.cnki.xsxb.2020.05.002
Ping KANG, Xian-ning LIU. On Dynamics of an SEIR Epidemic Model with Transport-Related Infection between Two Heterogeneous Cities[J]. Journal of Southwest China Normal University(Natural Science Edition), 2020, 45(5): 7-12. doi: 10.13718/j.cnki.xsxb.2020.05.002
Citation: Ping KANG, Xian-ning LIU. On Dynamics of an SEIR Epidemic Model with Transport-Related Infection between Two Heterogeneous Cities[J]. Journal of Southwest China Normal University(Natural Science Edition), 2020, 45(5): 7-12. doi: 10.13718/j.cnki.xsxb.2020.05.002

两个异质城市间具有路途感染的SEIR传染病模型

  • 基金项目: 国家自然科学基金项目(11671327)
详细信息
    作者简介:

    康萍(1993-), 女, 硕士研究生, 主要从事生物数学及动力系统理论及其应用研究 .

    通讯作者: 刘贤宁, 教授, 博士研究生导师
  • 中图分类号: O175

On Dynamics of an SEIR Epidemic Model with Transport-Related Infection between Two Heterogeneous Cities

计量
  • 文章访问数:  2073
  • HTML全文浏览数:  2073
  • PDF下载数:  352
  • 施引文献:  0
出版历程
  • 收稿日期:  2019-05-20
  • 刊出日期:  2020-05-20

两个异质城市间具有路途感染的SEIR传染病模型

    通讯作者: 刘贤宁, 教授, 博士研究生导师
    作者简介: 康萍(1993-), 女, 硕士研究生, 主要从事生物数学及动力系统理论及其应用研究
  • 西南大学数学与统计学院, 重庆 400715
基金项目:  国家自然科学基金项目(11671327)

摘要: 建立并分析了两个异质城市间具有路途感染的SEIR传染病模型.计算得出了基本再生数R0,证得当R0 < 1时,存在一个全局渐近稳定的无病平衡点;当R0>1时,存在一个地方病平衡点且系统是一致持久的.

English Abstract

  • 随着跨地区交流的日益密切,大量传染病通过种群迁移迅速扩散到世界各地,如流感、非典、甲型H1N1等[1].为了研究人群流动对传染病传播的影响,文献[2]建立了两斑块间具有路途感染的SIS传染病模型,文献[3-4]研究了具有路途感染和出入境检查的SIQS传染病模型.鉴于流感、水痘、麻疹等传染病的宿主从感染者变成患者期间会有一个潜伏期,并且康复后对疾病具有免疫力,文献[5-6]分析了具有路途感染的SEIS,SEIRS传染病模型的动力学行为.这些模型均采用相同的流行病学参数来描述不同斑块间传染病的传播.然而,各城市的人群密度、结构组成、资源分配等都是不同的.本文在文献[6]的基础上,进一步研究空间异质性及路途感染对传染病传播的影响.建立如下具有路途感染的SEIR模型:

    1) 将每个城市的总人口分为4类:易感者S、潜伏者E、患者I和免疫者R,则城市it时刻的总人口为Ni(t)=Si(t)+Ei(t)+Ii(t)+Ri(t).城市i单位时间内新增加的人口为ai,每类个体的自然死亡率为bi,且单位时间内每个潜伏者的发病率为ci,患者的因病死亡率为ei,患者的恢复率为dii=1,2.

    2) 城市i中疾病发生率为标准发病率$\frac{\beta_{i} S_{i} I_{i}}{N_{i}}$,其中βi为城市i中的接触率.两城市之间有直接的迁移,城市i到城市j的人均迁移率为αi(jii=1,2),迁移路途中的疾病发生率为

    其中γi为路途感染率,并且个体在迁移途中无出生和死亡.

    3) 易感者在迁移途中可能变成患者,但患者在迁移途中不会恢复,且个体在潜伏期没有传染能力,不能直接恢复成易感者.

    初始条件为

    注1  迁移过程中部分易感者会变成感染者,但剩余易感者数量必须非负,即对所有的SiEiIiRi≥0,αiSi-$\frac{\gamma_{i} \alpha_{i} S_{i} I_{i}}{N_{i}}$≥0,所以始终要求0≤γ≤1,i=1,2.

    注2  迁移途中的新感染不能大于每个城市内部的感染,即对所有的SiEiIiRi≥0,当0≤γi≤1时,βi>γiαii=1,2.

  • 若没有个体的迁移,即α1=α2=0,则模型(1)降维成简单的SEIR模型:

    由问题的实际含义,仅在区域

    中讨论系统(3)解的性态,容易验证D是系统(3)的可行域.易得系统(3)存在唯一无病平衡点

    利用下一代矩阵理论[7],可以得到基本再生数

    R0i>1时,系统(3)存在唯一的地方病平衡点P*i=(S*iE*iI*iR*i),其中,

    其中Λ=bici(R0i-1)(bi+ei)+R0ibi[bi(bi+di+ei)+cidi]+bi2cii=1,2.

    根据文献[6]中定理3.1和定理3.2的证明过程可得如下结论:

    定理1  当R0i < 1时,系统(3)的无病平衡点P0i局部渐近稳定;当R0i>1时,系统(3)的无病平衡点P0i不稳定,但是系统(3)的地方病平衡点P*i局部渐近稳定.

  • 定理2  在初始条件(2)下,系统(1)的解始终非负,并且是一致最终有界的.

      易证系统(1)的满足初值条件的解(S1E1I1R1S2E2I2R2)始终非负.根据系统(1)可得

    y(t)=N1(t)+N2(t),则

    从而$\mathop {\lim \sup }\limits_{t \to \infty } y(t)\frac{{{a_1} + {a_2}}}{{\min \left\{ {{b_1}, {b_2}} \right\}}}$$y(t) \leqslant \max \left\{\frac{a_{1}+a_{2}}{\min \left\{b_{1}, b_{2}\right\}}, y(0)\right\}$,因此系统(1)的满足初值条件的解(S1E1I1R1S2E2I2R2)是一致最终有界的.证毕.

    下面讨论平衡点的存在性及稳定性.

    系统(1)始终存在一个无病平衡点P0=(S10,0,0,0,S20,0,0,0),其中

    根据下一代矩阵理论可得系统(1)的基本再生数:

    其中Δ=(m3m4n1R01-m1m2n2R02+α1α2ϕ12)2+4α1α2(m4n1R01+ϕ13)(m2n2R02+ϕ14)>0,定义

    如果两城市间没有个体迁移,则

    对城市1,R0=R01;对城市2,R0=R02.

    如果在迁移途中没有路途感染,即γ1=γ2=0,则$\hat{\phi}_{11}$=β1c2+β2c1$\hat{\phi}_{12}$=β1c2-β2c1$\hat{\phi}_{13}$=β1c2m1$\hat{\phi}_{14}$=β2c1m3,显然

    即路途感染使得基本再生数增大,从而加剧了传染病的传播.

    定理3  当R0 < 1时,系统(1)的无病平衡点P0全局渐近稳定;当R0>1时,系统(1)的无病平衡点P0不稳定.

      根据文献[7]中定理2,很容易验证系统(1)满足假设(A1)-(A5),从而当R0 < 1时,无病平衡点P0是局部渐近稳定的;当R0> 1时,无病平衡点P0是不稳定的.下证P0的全局吸引性.由于

    考虑辅助的线性系统

    系统(6)有唯一的平衡点E*=(u1*u2*),其中

    平衡点E*处对应的特征方程的所有根均具有负实部,因此E*是局部渐近稳定的.由于系统(6)是一个合作的不可约系统,则$\lim\limits_{t \rightarrow \infty}\left(u_{1}(t), u_{2}(t)\right)=\left(S_{1}^{0}, S_{2}^{0}\right)$.

    u1(0)≥S1(0),u2(0)≥S2(0),由比较定理可得,对任意足够小的ε1> 0,存在t0> 0,当t>t0时,

    所以对t>t0,可以得到E1I1R1E2I2R2的方程

    显然,系统(7)的辅助线性系统有一个平衡点E0=(0,0,0,0,0,0).且辅助线性系统右边的系数矩阵为J=F-V,根据文献[7]中定理2证明过程,如果R0 < 1,则s(J)=s(F-V) < 0,所以矩阵J的所有特征值具有负实部,因此系统(7)的辅助线性系统的每个解均趋于零.由比较定理可得:Ei(t)0,Ii(t)0,Ri(t)0,t→∞,所以,系统(1)的极限系统为

    系统(8)有一个正平衡点E1=(S10S20),由系统(6)的分析可知E1是全局渐近稳定的,因此根据渐近自治系统理论[8],当R0 < 1时,系统(1)的无病平衡点P0是全局渐近稳定的.

    定理4  当R0>1时,系统(1)至少存在一个地方病平衡点P*=(S1*E1*I1*R1*S2*E2*I2*R2*),且存在一个常数ε> 0,使得系统(1)的每个满足初值条件Si(0)≥0,Ei(0)≥0,Ii(0)≥0,Ri(0)≥0且满足I1(0)+I2(0)> 0的解(S1(t),E1(t),I1(t),R1(t),S2(t),E2(t),I2(t),R2(t))满足

      由系统(1)可知,

    $\liminf\limits_{t \rightarrow \infty} S_{1}(t) \geqslant \frac{a_{1}}{2\left(\beta_{1}+b_{1}+\alpha_{1}\right)}$,同理可得$\liminf\limits_{t \rightarrow \infty}$S2(t)≥$\frac{a_{2}}{2\left(\beta_{2}+b_{2}+\alpha_{2}\right)}$.根据系统(1)可知,如果Ii(t)是持久的,则Ei(t)和Ri(t)也是持久的,从而只需证明$\liminf\limits_{t \rightarrow \infty}$Ii(t)≥εi=1,2.令

    则可以证明传染病关于(X0∂X0)是一致持久的.

    显然XX0都是系统(1)的正不变集,且∂X0X上是相对闭的,由定理2可知,系统(1)是点耗散的.定义

    根据文献[9],容易证明M={(S1E1,0,R1S2E2,0,R2)|Si≥0,Ei≥0,Ri≥0,i=1,2}.

    定义系统(1)从初值ϕ(0)=(S1(0),E1(0),I1(0),R1(0),S2(0),E2(0),I2(0)R2(0))∈X开始的解的ω极限集为ω(S1(0),E1(0),I1(0),R1(0),S2(0),E2(0),I2(0),R2(0)),则有

    对任意ϕ(0)∈MIi(tϕ)=0,∀t≥0.则对系统(1),$\lim\limits_{t \rightarrow \infty} E_{i}(t, \phi)=0$$\lim\limits_{t \rightarrow \infty} R_{i}(t, \phi)=0$.根据渐近自治半流理论[8],可得Si(t)满足系统(8),由系统(8)的分析可知$\lim\limits_{t \rightarrow \infty} S_{i}(t, \phi)=S_{i}^{0}$,所以Ω={P0}.显然P0Ω的一个孤立覆盖,且由于在M上不存在连接P0到它自身的解,所以是非周期的,下证

    这里Ws(P0)表示P0的稳定流形.假设(9)式不成立,则系统(1)存在一个满足初值条件的解{(S1(t),E1(t),I1(t),R1(t),S2(t),E2(t),I2(t),R2(t))∈X0,∀t≥0},使得

    从而对任意小的δ>0,且δ < Si0,存在t1> 0,使得对所有t>t1

    因此对t>t1,根据系统(1)可得

    其辅助的线性系统右边的系数矩阵为Jδ,且Jδ=J-δJ1,其中

    由于R0>1,则s(J)=s(F-V)> 0,取充分小的δ>0,使得当R0>1时,s(J-δJ1)>0.因此Jδ有一个具有正特征向量的正特征值s(Jδ),由比较定理可得

    与(10)式矛盾,故Ws(P0)∩X0=∅.所以系统(1)关于(X0∂X0)是一致持久的[10].根据文献[11],系统(1)至少存在一个正平衡点P*=(S1*E1*I1*R1*S2*E2*I2*R2*).定理得证.

  • 本文基于文献[6]研究了两异质城市间具有路途感染的SEIR传染病模型.证明了当R0 < 1时,两城市间的传染病不能入侵人群,当R0>1时,存在地方病平衡点且系统是一致持久的.

参考文献 (11)

目录

/

返回文章
返回