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

/传染病的数学建模与分析时间：2010年9月7日地点：2楼阶梯教室传染病建模的意义传染病历来就是威胁人类健康的大敌.人类征服传染病的道路依然曲折漫长。近20年来像AIDS病、SARS、禽流感等重大传染病相继爆发.在全球蔓延。2008手足口病的爆发曾给婴幼儿的健康带来了极大的危害。20XX的H1N1又来侵害年轻的我们。结核、白喉、鼠疫、登革热等一些老的传染病也重新抬头.给人们工作、生活和国民经济的发展带来了极大的影响。20XX突发的SARS传染病给我们的公共卫生体系应对突发性传染病提出了新的要求.也给数学在研究传染病动力学性态和预测等方面提出了一系列新问题。因此.研究和分析传染病传播的数量规律.建立有效的防控机制既是摆在我们面前的一个困难问题.也是一项紧迫任务。建立数学模型的目的是：描述传染病的传播过程；分析受感染人数的变化规律；预报传染病高潮到来的时刻；预防传染病蔓延的手段。基本的传染病动力学模型在传染病动力学中.长期以来主要使用的数学模型是所谓的"仓室"〔compartment模型.它的基本思想由Kermack与McKendrick创立于1927年.但一直到现在仍然被广泛的使用和不断地发展着。下面我们以他们提出的一个经典的基本模型为例.来阐述建立仓室模型的基本思想和有关基本概念.并显示由模型能得到的主要结论。Kermack-McKendrick的仓室模型所谓仓室模型就是针对某类传染病将该地区的人群分成以下三类〔即三个仓室：易感者〔Susceptibles类其数量记为.表示时刻未染病但有可能被该类疾病传染的人数。染病者〔Infectives类其数量记为.表示时刻已被感染成病人而且具有传染力的人数。移出者〔Removed类其数量记为.表示时刻已从传染病者类移出的人数。设总人口为.则有。K-M的模型是一个十分简单粗糙的模型。它的建立基于以下三个基本假设：〔1不考虑人口的出生、死亡、流动等种群动力因素。这意味着考虑一个封闭环境而且假定疾病随时间的变化要比出生、死亡随时间变化显著得多.从而后者可以忽略不计。这样.此环境的总人口始终保持为一个常数.即.或。〔2一个病人一旦与易感者接触就必然具有一定的传染力。这里假设时刻单位时间内.一个病人能传染的易感者数目与此环境内易感者总数成正比.比例系数为〔称为传染系数.从而在时刻单位时间内被所有病人传染的人数〔即新病人数为。〔3时刻.单位时间内从染病者类移出的人数与病人数量成正比.比例系数为.从而单位时间内移出者的数量为。显然.是单位时间内移出者在病人中所占的比例.称为移出率系数.当不致混淆时也简称为移出率。当移出者中仅包括康复者时.移出率系数又称为恢复率系数或简称为恢复率。在以上三个基本假设下.易感者从患病到移出的过程可用下述框图描述：对每一个仓室的人口变化率建立平衡方程式.便得到以下模型：〔1-1下面.我们通过对模型〔1-1的分析和解的渐近性态研究来初步显示动力学模型对认识传染病流行规律所起的作用。将〔1-1中三个方程两端分别相加.得.从而由于〔1-1中前两个方程中不含.故实际上我们只需先讨论前两个方程：〔1-2由于.单调递减且有下界〔为0.故极限存在。由〔1-2有.〔1-3可见.当时.达到极大值。可见.当初始时刻易感者数量时.随时间增长.染病者数量将先增加达到最大值.然后再逐渐减少而最终消亡。这一现象表明.只要.即.就有.疾病就会流行。令〔1-4则当时.疾病流行；当时.疾病不会流行.染病者数量将单调下降而趋向于零。是区分疾病流行与否的阈值。表示在发病初期.一个病人在传染期内所传染的人数.称为基本再生数〔具有很强的生物学意义。应当指出.〔1-4中的表示平均移出时间.也就是平均患病期。事实上.由移出率系数的定义可见.若病人数量为.则单位时间内移出者的数目为.故经过时间.病人全部移出。要防止疾病流行.必须减少使它小于1.由表达式〔1-4可知.这可以通过加强治疗以缩短染病期或采取杀菌等措施以减少疾病的传染力.或通过隔离措施以减少与患病者可能接触的人数即这里的易感者来实现。更为有效的方法是通过疫苗接种以使易感者成为免疫者而直接进入移出者类.从而减少初始时刻易感者的数量。根据疾病的病理特点.我们可以建立各种模型：1模型。患病后可以治愈.然后直接变为易感者。2模型。病人康复后只有暂时免疫力.单位时间内将有的康复者丧失免疫而可能再次被感染。3模型。在被感染后成为患病者之前有一段病菌潜伏期.并且假定在潜伏期内的感染者没有传染力。记时刻潜伏期的人数为.疾病的平均潜伏期为。4模型.病人康复后仅有暂时免疫力5无垂直传染模型。即母体的疾病不会先天传染给新生儿.故新生儿均为易感者。6〔有垂直传染且康复者的新生儿不具有免疫力模型SARS的建模与预测SARS是由一种冠状病毒引起的传染性很强的呼吸道传染病.它主要通过近距离空气飞沫以及接触病人呼吸道分泌物和密切接触进行传播.也可能通过病人飞沫污染物、如通过手、衣物、食物、水或环境等途径传播。SARS潜伏期一般为2-11天.在潜伏期无感染。SARS患者的主要症状有：发热〔体温38℃以上为首发症状.多为高热.并可持续1-2周以上.可伴有寒战或其他症状.包括头痛、全身酸痛和不适、乏力.部分病人在早期也会有轻度的呼吸道症状<如咳嗽、咽痛等。SARS在前面的模型中.感染率等都是常量.也没有考虑因病死亡率.这样的模型对感染情况变化不大、因病死亡极少的情况比较适用〔如HFMD。但对SARS而言.从20XX4月下旬开始.我国各地普遍采取了十分严格的控制措施.使得一个患者能感染他人的情况大大降低.传染率有了明显的变化。另外.SARS传染者的死亡率相对很高.是正常人口死亡率的十多倍.SARS患者的治疗和康复情况也随着经验的增加而大大改进。所以.我们在模型中需引入随时间变化的感染率、恢复率和因病死亡率。由于SARS感染者在潜伏期内不易被发现.且潜伏期中不感染或感染很少.所以在这个简易模型中我们不考虑潜伏期。在建立模型时我们讲总人口分为易感者、感染者和恢复者3类.并做以下假设：单位时间内感染人数与现有感染者与现有感染者与易感者的乘积成比例；单位时间内治愈恢复的感染者与现有的感染者成比例；单位时间内死亡的感染者与现有的感染者成比例；SARS感染者治愈后不会再被感染；所有的新生儿都是易感者。根据SARS传播的规律和这些假设.我们可以得到下面的框图〔图1-3和模型。图1-3SARS传播的框图〔1-5其中为时刻的人口总数.是总人口的出生率.是SARS感染者的因病死亡率.其它记号和参数的含义同前。模型中随时间变化的感染率、治愈率和因病死亡率可以使我们更加方便和实际地描述SARS的传播规律。模型〔1-5应用的最大困难是易感者的处理。理论上讲.每个未被SARS感染的人都是易感者.我国13亿人口中除过当时已感染的几千人外都是易感者。而实际情况使每个人的活动能力和区域有限.SARS是通过接触传播的.每个病人根本不可能与全国13亿人等可能地接触。另一方面.如果将13亿人作为易感者带入模型计算时.由于13亿易感者和几千人的感染者的差距如此之大.很难得到合理的计算结果。克服这一困难的途径是选取有可能和易感者接触的人员作为易感者.这些人包括医护人员、感染者的亲属、朋友、同事、旅伴和同去过一些公共场所的人员。由于这些人员也很难统计或估计.所以我们采用另外一条途径。在一个疾病的传播过程中.人们最感兴趣的问题是当前有多少感染者.以后将会有多少感染者.所以我们将注意力集中在描述感染者数量变化的模型〔1-5中的第二个方程上。根据国家每天公布疫情数据这一特点.我们以天为时间单位将〔1-5的第二个方程离散化得〔1-6其中是第天SARS感染者数量.是第天新增的感染人数.可以看作每个SARS感染者在第天所感染的人数.为第天SARS感染者治愈和死亡的人数。记.我们可以将〔1-6简化为〔1-7如果我们知道了某一天SARS感染者的数量.并且知道了函数的具体表达式.我们就可以用〔1-7来递推预测SARS以后随时间的变化规律。因为国家每天公布疫情数据.〔1-7中所用的初始值很容易得到.困难在于如何确定。由〔1-7可以得到具有下面的表达式〔1-8如果我们知道了若干天SARS感染者的数据.就可以从〔1-8中得到。即当SARS感染的数据公布后.就是已知的函数了。我们需要的是要对未来SARS的感染情况进行预测.不能等实际统计结果出来以后进行回顾。所以我们根据已有的SARS感染数据用〔1-8计算出在这一段时间内的值.再根据在这一段时间内的值进行回归.得到的函数形式后再代到〔1-7就可以预测未来。通过搜索得到我国在20XX4月21日的SARS感染人数为2158人.再根据附件中的数据可以计算出各天的SARS感染人数.利用〔1-8中的公式计算得到三周内的值。对这三周内的值利用指数曲线回归可以得到的函数表达式.其回归的结果见图1-4。图1-4每天每个SARS病人平均感染的人数.光滑曲线为拟合曲线.折线为实际计算结果将这样回归得到的代入到〔1-7中.就可以得到以后各天SARS感染者的人数和新增感染者人数。其预测结果见图1-5和图1-6。图1-5每天在医病人数随时间变化的关系.绿线：实际统计值.蓝线：预测值图1-6每天新增病人随时间变化的关系.绿线：实际统计值.蓝线：预测值这一预测结果我们于20XX5月21日向媒体公布.到SARS在我国的流行结束后.我们收集了2个月SARS感染者的数据.对比我们的模型和3周的数据得到的所给出的预测是相当成功的。我们的预测人