Modeling_传染病的数学建模与分析报告

1、传染病的数学建模与分析时间 ：2010 年 9 月 7 日 地点 ：2 楼阶梯教室一、 传染病建模的意义传染病历来就是威胁人类健康的大敌 .人类征服传染病的道路依然曲折漫长 。近 20 年来像 AIDS 病、SARS、 禽流感等重大传染病相继爆发 .在全球蔓延 。2008 手足口病的爆发曾给婴幼儿的健康带来了极大的危害 。2009 年的 H1N1 又来侵害年轻的我们 。结核 、白喉 、鼠疫 、 登革热等一些老的传染病也重新抬头 .给人们工作 、生活和国民经济的发展带来了极大的影响 。2003 年突发的 SARS 传染病给我们的公共卫生体系应对突发性传染病提出了新的要求 .也给数学在研究传染病动

2、力学性态和预测等方面提出了一系列新问题 。因此.研究和分析传染病传播的数量规律 .建立有效的防控机制既是摆在我们面前的一个困难问题 .也是一项紧迫任务 。建立数学模型的目的是 ： 描述传染病的传播过程 ； 分析受感染人数的变化规律 ；预报传染病高潮到来的时刻 ； 预防传染病蔓延的手段 。二、 基本的传染病动力学模型在传染病动力学中 .长期以来主要使用的数学模型是所谓的 “仓室 ”（compartment ）模型.它的基本思想由 Kermack 与 McKendrick 创立于 1927 年.但一直到现在仍然被广泛的使用和不断地发展着 。下面我们以他们提出的一个经典的基本模型为例 .来阐述建立仓

3、室模型的基本思想和有关基本概念 .并显示由模型能得到的主要结论 。Kermack-McKendrick 的 S I R仓室模型所谓 S I R 仓室模型就是针对某类传染病将该地区的人群分成以下三类 （即三个仓室）：易感者 （Susceptibles ）类 其数量记为 S(t) .表示 t 时刻未染病但有可能被该类疾病传染的人数 。染病者 （Infectives ）类 其数量记为 I (t ) .表示 t 时刻已被感染成病人而且具有传染力的人数 。. .移出者 （Removed ）类 其数量记为 R(t ) .表示 t 时刻已从传染病者类移出的人数 。设总人口为 N(t) .则有 N (t) S

4、(t) I (t) R(t ) 。K-M 的 SIR模型是一个十分简单粗糙的模型 。它的建立基于以下三个基本假设 ：（1）不考虑人口的出生 、 死亡 、流动等种群动力因素 。 这意味着考虑一个封闭环境而且假定疾病随时间的变化要比出生 、 死亡随时间变化显著得多 .从而后者可以忽略不计 。这样.此环境的总人口始终保持为一个常数 .即 N(t) K .或S(t) I (t) R(t ) K 。（2） 一个病人一旦与易感者接触就必然具有一定的传染力 。 这里假设 t 时刻单位时间内.一个病人能传染的易感者数目与此环境内易感者总数 S(t ) 成正比 .比例系数为 （称为传 染系 数 ） .从 而在

5、t 时 刻 单位 时 间内 被 所有 病 人传 染的人 数（ 即新 病 人数 ） 为S(t) I (t) 。（3） t 时刻 .单位时间内从染病者类移出的人数与病人数量成正比 .比例系数为 .从而单位时间内移出者的数量为 I (t )。显然 . 是单位时间内移出者在病人中所占的比例 .称为移出率系数 .当不致混淆时也简称为移出率 。当移出者中仅包括康复者时 .移出率系数又称为恢复率系数 或简称为 恢复率 。在以上三个基本假设下 .易感者从患病到移出的过程可用下述框图描述 ：对每一个仓室的人口变化率建立平衡方程式 .便得到以下模型 ：dSdtSIdIdtSI IdRdtI（1-1 ）下面.我们通

6、过对模型 （1-1 ） 的分析和解的渐近性态研究来初步显示动力学模型对认识传染病流行规律所起的作用 。. .将（1-1 ） 中三个方程两端分别相加 .得d(S I R)dt0.从而S(t ) I (t) R(t) K（常数）由于 （1-1 ） 中前两个方程中不含 R.故实际上我们只需先讨论前两个方程 ：dSdtSIdIdtI ( S )（1-2 ）dS由于 0dt. S(t ) 单调递减且有下界 （为 0）.故极限lim S(t) S t存在 。由（1-2 ）有dI1dS S. （1-3 ）可见.当 S 时. I 达到极大值 。可见.当初始时刻易感者数量S(0) S 时.随时间增长 .0染病者

7、数量 I (t) 将先增加达到最大值 I ( ) .然后再逐渐减少而最终消亡 。 这一现象表明 .只要S .即01 dIS 1.就有 I ( S ) 00dt.疾病就会流行 。令1 S0R S （1-4 ）0 0则当R0 1 时.疾病流行 ；当 R0 1时.疾病不会流行 .染病者数量 I (t) 将单调下降而趋向于零。R0 1是区分疾病流行与否的 阈值。R 表示在发病初期 .一个病人在传染期内所传染的人数 .称为 基本再生数 （ 具有很强的0生物学意义 ）。. .1 表示 平均移出时间 .也就是 平均患病期 。事实上 .由移出率系数应当指出 .（1-4 ）中的1的定义可见 .若病人数量为 n

8、.则单位时间内移出者的数目为 n.故经过时间.病人全部移出。要防止疾病流行 .必须减少R 使它小于 1.由表达式 （1-4）可知.这可以通过加强治疗以01缩短染病期 或采取杀菌等措施以减少疾病的传染力 .或通过隔离措施以减少与患病者可能接触的人数即这里的易感者S 来实现 。更为有效的方法是通过疫苗接种以使易感者成为0免疫者而直接进入移出者类 R .从而减少初始时刻易感者的数量 S0 。根据疾病的病理特点 .我们可以建立各种模型 ：1） SIS模型 。患病后可以治愈 .然后直接变为易感者 。2） SIRS模型 。 病人康复后只有暂时免疫力 .单位时间内将有 R 的康复者丧失免疫而可能再次被感染

9、。3） SEIR模型 。 在被感染后成为患病者 I (t) 之前有一段病菌潜伏期 .并且假定在潜伏期内的感染者没有传染力 。记t 时刻潜伏期的人数为 E(t) .疾病的平均潜伏期为1。. .4） SEIRS模型.病人康复后仅有暂时免疫力5） SIR无垂直传染模型 。 即母体的疾病不会先天传染给新生儿 .故新生儿均为易感者。6） SIR（有垂直传染且康复者的新生儿不具有免疫力 ）模型三、 SARS的建模与预测SARS 是由一种冠状病毒引起的传染性很强的呼吸道传染病 .它主要通过近距离空气飞沫以及接触病人呼吸道分泌物和密切接触进行传播 .也可能通过病人飞沫污染物 、 如通过手、 衣物 、食物 、

10、水或环境等途径传播 。SARS 潜伏期一般为 2-11 天.在潜伏期无感染 。SARS 患者的主要症状有 ：发热 （体温 38以上 ）为首发症状 .多为高热 .并可持续 1-2 周以上.可伴有寒战或其他症状 .包括头痛 、 全身酸痛和不适 、乏力.部分病人在早期也会有轻度的呼吸道症状 (如咳嗽 、 咽痛等 ）。 SARS 患者治愈后不会再被感染 。. .在前面的模型中 .感染率 等都是常量 .也没有考虑因病死亡率 .这样的模型对感染情况变化不大 、因病死亡极少的情况比较适用 （如 HFMD ）。 但对 SARS而言 .从 2003 年4 月下旬开始 .我国各地普遍采取了十分严格的控制措施 .使

11、得一个患者能感染他人的情况大大降低 .传染率有了明显的变化 。另外 .SARS 传染者的死亡率相对很高 .是正常人口死亡率的十多倍.SARS 患者的治疗和康复情况也随着经验的增加而大大改进 。所以.我们在模型中需引入随时间变化的感染率 、 恢复率和因病死亡率 。由于 SARS 感染者在潜伏期内不易被发现 .且潜伏期中不感染或感染很少 .所以在这个简易模型中我们不考虑潜伏期 。 在建立模型时我们讲总人口分为易感者 、 感染者和恢复者 3 类.并做以下假设 ：（1） 单位时间内感染 人数与现有感染者与现有感染者与易感者的乘积成比例 ；（2） 单位时间内治愈恢复的感染者与现有的感染者成比例 ；（3）

12、 单位时间内死亡的感染者与现有的感染者成比例 ；（4） SARS 感染者治愈后不会再被感染 ；（5） 所有的新生儿都是易感者 。根据 SARS 传播的规律和这些假设 .我们可以得到下面的框图 （图 1-3 ） 和模型 。图 1-3 SARS 传播的 SIR框图dS(t )dtbN(t ) S(t ) (t )S(t )I (t)dI (t)dt(t)S(t )I (t ) ( (t) (t) I (t )dR(t )dt(t) I (t) R(t )（1-5 ）其中 N(t) S(t) I (t) R(t ) 为 t时刻的人口总数 . b是总人口的出生率 . (t ) 是 SARS 感染. .

13、者的因病死亡率 .其它记号和参数的含义同前 。 模型中随时间变化的感染率 (t) 、 治愈率(t) 和因病死亡率 (t) 可以使我们更加方便和实际地描述 SARS 的传播规律 。模型 （1-5 ） 应用的最大困难是易感者的处理 。理论上讲 .每个未被 SARS感染的人都是易感者 .我国 13 亿人口中除过当时已感染的几千人外都是易感者 。 而实际情况使每个人的活动能力和区域有限 .SARS是通过接触传播的 .每个病人根本不可能与全国 13 亿人等可能地接触 。另一方面 .如果将 13 亿人作为易感者带入模型计算时 .由于 13 亿易感者和几千人的感染者的差距如此之大 .很难得到合理的计算结果

14、。克服这一困难的途径是选取有可能和易感者接触的人员作为易感者 .这些人包括医护人员 、 感染者的亲属 、朋友 、同事 、 旅伴和同去过一些公共场所的人员 。 由于这些人员也很难统计或估计 .所以我们采用另外一条途径 。在一个疾病的传播过程中 .人们最感兴趣的问题是当前有多少感染者 .以后将会有多少感染者 .所以我们将注意力集中在描述感染者数量变化的模型 （1-5 ） 中的第二个方程上 。根据国家每天公布疫情数据这一特点 .我们以天为时间单位将 （1-5 ）的第二个方程离散化得I (t 1) I (t) (t )S(t )I (t) ( (t) (t) I (t)（1-6 ）其中 I (t) 是

15、第 t 天 SARS 感染者数量 . (t)S(t)I (t) 是第 t 天新增的感染人数 . (t )S(t ) 可以看作每个 SARS 感染者在第 t 天所感染的人数 .( (t) (t) I (t) 为第 t 天 SARS 感染者治愈和死亡的人数 。记 f (t) (t) S(t) ( (t) (t) .我们可以将 （1-6 ）简化为I (t 1) I (t) f (t )I (t)（1-7 ）如果我们知道了某一天 SARS 感染者的数量 .并且知道了函数 f (t) 的具体表达式 .我们就可以用 （1-7 ） 来递推预测 SARS 以后随时间的变化规律 。 因为国家每天公布疫情数据 .

16、（1-7 ） 中所用的初始值很容易得到 .困难在于如何确定 f (t) 。由（1-7 ） 可以得到 f (t) 具有下面的表达式. .f (t)I (t 1) I(t)I(t)（1-8 ）如果我们知道了若干天 SARS 感染者的数据 .就可以从 （1-8 ） 中得到 f (t) 。即当 SARS感染的数据公布后 . f (t) 就是已知的函数了 。 我们需要的是要对未来 SARS 的感染情况进行预测 .不能等实际统计结果出来以后进行回顾 。 所以我们根据已有的 SARS 感染数据用（1-8 ） 计算出 f (t) 在这一段时间内的值 .再根据 f (t) 在这一段时间内的值进行回归 .得到f

17、(t) 的函数形式后再代到 （1-7 ） 就可以预测未来 。 通过搜索得到我国在 2003 年 4 月 21日的 SARS 感染人数为 2158 人.再根据附件中的数据可以计算出各天的 SARS感染人数 .利用（1-8 ） 中的公式计算得到三周内 f (t) 的值 。对这三周内 f (t) 的值利用指数曲线回归可以得到 f (t )的函数表达式 .其回归的结果见图 1-4 。图 1-4 每天每个 SARS病人平均感染的人数 .光滑曲线为拟合曲线 .折线为实际计算结果将这样回归得到的 f (t) 代入到 （1-7 ）中.就可以得到以后各天 SARS 感染者的人数和新增感染者人数 。 其预测结果见

18、图 1-5 和图 1-6 。. .图 1-5 每天在医病人数随时间变化的关系 .绿线 ：实际统计值 .蓝线 ：预测值图 1-6 每天新增病人随时间变化的关系 .绿线 ：实际统计值 .蓝线 ：预测值这一预测结果我们于 2003 年 5 月 21 日向媒体公布 .到 SARS 在我国的流行结束后 .我们收集了 2 个月 SARS 感染者的数据 .对比我们的模型和 3 周的数据得到的 f (t) 所给出的预测是相当成功的 。我们的预测人数和实际统计人数的对比见图 1-7 。. .图 1-7 每天在医病人随时间变化的关系 .紫线 ：实际统计值 .绿线 ：预测值对建立的 SARS 在我国的传播模型 .进行了理论研究并研制了预测与控制分析软件 .于2003 年5 月21 日向新闻媒体发布了我们的研究结果 .预测按世界卫生组织 （WHO ） 的标准 .我国将于 6 月下旬解除旅游警告 。 届时全国累计感染者人数为 6000 人左右 .这与后来的实际情况是吻合的 .而世界卫生组织宣布解除北京市旅行警告的日期也正是 2003 年 6 月 24 日。. .欢 迎 您的光临 ，Word 文 档下载 后可修 改编辑 双击可删 除页眉 页脚 谢谢 ！希 望您提 出您宝 贵的意