一类带接种疫苗的SEIR传染病模型的定性分析

1、一类带接种疫苗的SEIR传染病模型的定性分析摘要众所周知，传染病严重影响人类的健康，并在一定程度上会阻碍社会经济的发展。因此,寻求疾病的最优预防与控制策略是当今社会重要的研究课题。而传染病动力学则是对传染病进行定量研究的一门重要学科。早在 1927年，Kermack与McKendrick利用动力学的方法建立 了 SIR仓室模型。当考虑易感者（S）在感染病菌后不会立即发病，例如： HIV病毒，这样即可建立SEIR模型。本文主要研究了一类带有接种疫苗的SEIR传染病模型的定性分析。首先，通过研究模型本身的等价系统，给出了基本再生数的显式表达式；其次，应用常 微分方程的平衡点局部稳定性及全局稳定性的

2、判断依据，讨论了无病平衡点以及地方病平衡点的存在性和稳定性。关键词：SEIR传染病模型，基本再生数.无病平衡点,地方病平衡点引言传染病是由细菌、?菌或病毒等病原体或蠕虫等寄生虫感染人或者其他生物后所产生的能在种群中相互传播的疾病。历史上一次又一次的大规模传染病盛行给人类生存和国民生计带来了巨大的灾难。长期以来，人类都在与传染病进行不屈不挠的斗争。20世纪以来，人 类在征服传染病的路上取得了辉煌的成就。 肆虐千年的天花病毒被消火了： 麻风、脊髓灰质 炎消失的日子 也不远了；百日咳、白喉等疾病已经在许多国家得到了遏制； 众多抗生素的问 世，也使得一度令 人闻风丧胆的瘟疫不能再危害社会。然而世界卫生

3、组织发表的报告中表明，传染病依旧是危害人类健康的第一杀手。以95年的数据为例，全世界5200万死亡人口中，丧命于传染病的有 1700万，将近三分之一。近 20年 来，AIDS、霍乱、疯牛病、SARS和甲型Hll流感 等恶性突发疾病给人类社会带来了巨大危 害，具至一些老的传染病如鼠疫等也死灰复燃，这种悄况已经引起了全球的高度重视，如何对这类传染病发展做出科学 预测，并实施有效手段的问题已经得到业界工作者的普遍关注。数学作为一门基础学科，到如今已经渗透到了科学研究领域的各个方面。近些年来，数学1模型已经被广泛的运用到了对生态系统的研究中来。生物学中有众多有趣的模型，这些模型在数学中大多可以被归纳为

4、微分方程组或方程，其中传染病动力学模型与我们的生活息息相关。对传染病传播特性的的研究，特别具有里程碑式意义的工作是1927年 Kermack和Mckendrick提出的仓室建模方法，该法一直沿用至今，并发展成为一种专门的理论：传染病动力学。传染病动力学是大众研究传染病的非常重要的方法，传染病动力学模型依据种群生长的特点、疾病产生及在种群内传播的规律，并且考虑周圉环境的变化等 因素，建立起能反映传染病动力学特征的数学模型，而通过对模型动力学性态的 定性、定量和数值模拟，来显示疾病的发展过 程，揭示疾病的流行规律同时还可 以预测其发展趋势，以便找到对其预防和控制的最优方案，为人类制定防治决策 提供

5、有效的理论和数量方面的依据。近 30年来，国际上传染病动力学的研究 进 展已经十分迅速了。从已有的研究结果看，主要集中在常微分方程形式的SI, SIR, SEIR,MSEIR等模型。SIR （易感者，感染者，恢复类）是常见的传染病模型。可是，许许多多的 传染病在发病 前都会在宿主体内潜伏一段时间。比如，麻疹病毒的潜伏期是8-13天，狂犬病毒的潜伏期少则儿天多则达数年。因此我们顾忌到传染病毒的潜伏期,就需要考虑一个更适合贴切的传染病模型，这个模型引进了一个新的传染病群 体，即潜伏者E （t）,即该传染病模型为SEIR （易感者，潜 伏者，感染者，恢复 类）模型，这个模型具有十分重要的实际意义。在

6、传染病的实际传播过程中，不同类型的疾病通常有不同的传播特点，部分疾病的传播和年龄有着莫大的关系，而有些疾病对不同年龄的人的传染程度不同；有些疾病的传播与预防疾病的描施，比如接种疫苗有着很大关系。在我国，通过接种疫苗，许多传染病如麻疹、破伤风、口日咳等的预防和治疗已经取得了不错的成绩。当然，事实说明部分接种疫苗的人群在接种疫苗之后有可能会再次感染该疾病。文中建立和研究了带有接种疫苗的 SEIR传染病模型，并且考虑到了接种疫 苗之后有可能 被再次感染的情况，运用常微分方程中稳定性理论的知识和思想， 在一定的条件假设下，对这一 类传染病模型的动力学性质进行了粗浅的分析。 求出了基本再生数的表达式，

7、证明了时，系统存 在唯一的无病平衡点并讨论了其稳 定性；及时，无病平衡态不稳定，此时存在地方病平衡点， 并 研究分析了其稳定 性。将其投入到实际中，希望能对疾病的预防消除起到一定作用。本文安排如下，第一节在一定的假设条件下， 建立了带有接种疫苗的 SEIR传染病模型，并 且简要的介绍了数学建模的理论知识；第二节分析并证明了模型的初值问题和一些基本性质；第三节求出了再生数的表达式，得到了无病平衡点和地方病平衡点的存在条件，并且在一定条件下分析了其稳定性。第一章模型的建立假设是在时刻种群内易感者数量，是在时刻种群内已经接种疫苗的人数，是在时刻的时候群体内的疾病潜伏者人数，是在时刻种群内感染者的人数

8、，康复者为错误！未找到引用源。，则在时刻种群的总量可表示为进一步，本文所研究模型的建立基丁 ？以下假设：(1)易感染者口然增加量包括新生儿和移入人口为，11 一旦易感染者与潜伏者或感染者存在有效接触，即有可能染上疾病，儿率是错误！未找到引用源“错误！未找到引用 源。是易感 染者与感染者接触患病的儿率，错误！未找到引用源。是易感染者与潜伏者 接触后可能染病的儿 率与错误！未找到引用源。的比值：易感染者接种疫苗的速度是错 误！未找到引用源。；种群的口 然死亡率为，因此，易感者满足假设疫苗的成功率弁不是百分之百，接种疫苗的人群还有可能患病，几率是未找到引 用潦.？错误味找到引用源。(错误！未找到引用

9、源。)是接种人群成功的概率，故此时在与感染者与潜伏者接触时，人群中的潜伏者数量增多，而潜伏者数量还受自然死亡率及传染病影响，假设错误！未找到引用源。为潜伏者成为感染者的速度，故此时(4)感染者来源于病变的潜伏者，同时受限于疾病和自然死亡以及治愈率，假设该病死亡 率为d,治愈率为错误！未找到引用源。，则有上文各量之间转换关系用图表呈现为图1(pS(5)设错误！未找到引用源。为该病治愈率，则康复者满足dihv从而，有如下传染病模型=0a )s_( i在（1? 1）中，前四个等式并未涉及康复者人数，因此，只需讨论如下传染病模dS di dV孑=声_ (1 _ 7)0(也尹 _ 0 P一 =0 Q +

10、 sE)S - (1 - b )0Q + sEE bE_k di Q=&E - (d+b?Ldt(1.2)第二章初值问题2.1解的正性及有界性从实际情况出发，可假设错误！未找到引用源。，错误！未找到引用源。.错误！未找到引用源。，错误味找到引用源。？错误！未找到引用源。令错诫未找到引用以定理1对于一切错误！未找到引用源。，模型 （1.2）的解错误！未找到引用源。 也是非负的。同时，lim$upN (r) ocU证明：假设错误！未找到引用源。，显然 错误！未找到引用源。由一个方程可知dS(A(t) + b + O,V O,EOJ 0在此定义下，有如下结论：定理2对于（1.2）,如果初始值在 错误

11、！未找到引用源。中时，则可行域错误 ！未找到引用 源。是不变的。第三章基本再生数及无病平衡点易知系统（1.2）存在无病平衡点错误！未找到引用源。）其 中6,错误！未找到引用源。，错误！未找到引用源令错误！未找到引用源。，则(1.2)可写为一=八(x)-X(x)其中PG + 畔 4-(1- O) p(7 + 定)V000定义错误！未找到引用源v=(灯+b o ) -ki b+ d + d )v-=l(幻+b)+d)则(1.2)的基本再生数为_ L?0+Ci-ffX f ,+- a)VoJ/?v丁 (灯十 b)(b + d 十 6)错误！未找到引用源。(3.1)对于传染病模型的基本再生数是刻画疾病

12、是否传播的重要域值。山于错误味找到引用源。的每一部分都是非负的，易知特别地，可证 明如下结论: 定理3当错误！未找到引用源。时，系统（1.2）的无病平衡点错误！未 找到引用源。是局部渐近稳定的，“后昔误！未找到引用源。，E是不稳定的。证明：已知错误味找到引用源。，且错误味找到引用源。，则由（1.2）的最后两个方程可 得为了得到&0,对所有的心，有错误味找到引用源。则由（1.2）的笫一个方程立即得 到于是有另外，由（1? 2）的第二个方程可得从而山错误！未找到引用源。的任意性可得到于是，当错误！未找到引用源。时，可得到因此、勺错误味找到引用源。时，错误 ！未找到引用源。在可行域错误！未找到引 用

13、源。是局部渐近稳定的，反之则不稳定，得证第四章地方病平衡点的存在性及稳定性注意到，当时，系统 （1.2）存在地方病平衡点令错误！未找到引用源。，则地方病平衡点满足错误！未找到引用源。（2其中则方程（2.5）的正实根即为（1.2）的地方病平衡点，容易看出1 ?当错误！未找到引用源。，模型 （1.2）存在唯一地方病平衡点；2? X错误味找到引用源。，且错误味找到引用源。或者错误味找到引用源。时，则模型存在一个地方病平衡点；3?半错误味找到引用源。，错误！未找到引用源。，则模型存在两个地方病平衡点；4?当错误！未找到引用源。且错误！未找到引用源。时，模型没有地方病平衡点但是，事实上，当错误！未找到引用源。时，模型没有地方病平衡点。定理4当错误！未找到引用源。是，系统 （1.2）不存在地方病平衡点。证明：假设错误！未找到引用源。，则有如下关系式两式相比可以得到10cr(b + 初、b + r(b + 如 (b + ocp)o化简可得到b2 + bop + J:( p(b + 0) 0该不等式与参数为正相矛盾，因此，错误 ！未找到引用源。时模型没有地方病平衡点。当错误第五未找到引用源。时，可以大致确定二次函数图像，没有正根，故这种情况也不会有地方病平衡 点，得证章总结本文主要研究了一类