《658.甲型H1N1流感预测模型研究》由会员分享，可在线阅读，更多相关《658.甲型H1N1流感预测模型研究（17页珍藏版）》请在金锄头文库上搜索。
1、甲型H1N1流感预测模型研究 【摘要】传染病（甲型H1N1流感等）的控制预防历来都是人们所关注的焦点。每当传染性很强的流感爆发后，人们首先想到的是注射疫苗。但是据美国疾病预防控制中心工作人员透露，此次流感疫苗的大规模生产面市最快的时间是9月份，在此之前人类处于免疫的真空状态。在这期间，流感会以某种速度在某种程度范围内传播是社会各界比较关注的问题，也是关系到患者生命和应疫苗数量的多少的问题。传染病问题数学建模的研究是按照一般的传播机理建立集中模型，主要是建立微分方程模型，一般有SI模型、SIS模型、SIR模型。本文首先拟采用SIR模型1分析、预测美国甲型H1N1流感的传播趋势。我们发现SIR模型
2、在实际应用中有很多的局限，特别是整个传染期间内每个病人有效接触的平均人数称为接触数的（）的确定极其困难，由此得出的阈值在确定传染病是否会蔓延以及蔓延的程度问题方面也带来了很多技术性的困难。本文第二个数学模型采用时间序列计量经济学模型来分析甲型H1N1的传播趋势。流感传播的方式由两个：一是自发，二是由患者传染给健康的人，并且以由患者传染给健康的人为主要方式。每天新增流感患者的数量与前期患流感患者的人数有直接的关系，两者以一定得数学模型联系在一起的。通过时间序列分析，在没有疫苗出现的情况和其他条件不变的假设下，我们建立了模型：。这样就可以预测近期的甲型H1N1经实验室确认患者的数量。我们预测5月1
3、6号美国经实验室确认的流感患者将达到人。问题二和问题三都是基于SIR模型进行修正，并且得出结论。1、我们认为随着的增大，减小得越来越快。在现实中的意义是，随着投放疫苗数量的增加，患病的人数减少的速度变得越来越快。2、如果病毒变异导致潜伏期变长，一个患者将流感传染到别人身上的可能性大大增加，也就是说平均感染人数将增加。由此可以看出流感蔓延的程度将增加，流感疫情将持续更长的时间，流感疫情也变得更加严重。1 问题重述由墨西哥发端的甲型H1N1流感疫情目前正成为人们关注的焦点。据5月13日的最新统计，全球现有33国家（或地区）出现确诊病例，共有确诊病例5887人。你们可以通过有关网站获得最新的数据，如
4、新浪网新闻中心（网址：http:/ 问题分析甲型H1N1流感的传播是一道传染病问题。在数学建模领域已经有很多有关于这方面的研究，其中SIR模型是比较完整的模型。SIR模型通过建立微分方程模型，按照一般的传播机理建立集中模型。本文先采用SIR模型来分析，传染病的一般传播过程，即从蔓延到消失的过程。相关资料1也指出：我们看到在SIR模型中是一个重要参数。实际上m、k很难估计，而当一次传染病结束后，可以获得和，根据相关的公式得。当同样的传染病到来时，如果估计量没有多大的变化，那么就可以用上面得到的分析这次传染病的蔓延过程。问题就出现在这里：一、世界上没有完全同样的传染病，即与患者接触后，染病的概率肯
5、定不同。二、m会随着交通方式的发展、经济一体化的发展而增加，k却会随着医疗技术的发展而增加，这样=k/m的变化方向就不确定了。本文第二个数学模型采用时间序列计量经济学模型来分析甲型H1N1的传播趋势。问题二和问题三都是基于SIR模型进行修正，并且得出结论。1、我们认为随着的增大，减小得越来越快。在现实中的意义是，随着投放疫苗数量的增加，患病的人数减少的速度变得越来越快。2、如果病毒变异导致潜伏期变长，一个患者将流感传染到别人身上的可能性大大增加，也就是说平均感染人数将增加。由此可以看出流感蔓延的程度将增加，流感疫情将持续更长的时间，流感疫情也变得更加严重。3 建模准备3.1 数据收集与分析3.
6、1.1 美国甲型H1N1流感实验室确认病例数量时间确诊（包括死亡病例） 死亡(累计)4月23日504月24日804月25日1104月26日2004月27日4004月28日6404月29日9114月30日10915月1日14115月2日16015月3日22615月4日27915月5日40315月6日64225月7日89625月8日163925月9日225425月10日253235月11日260035月12日300935月13日335245月14日429845月15日47144感染甲型H1N1流感的患者平均可感染1.4至1.6人。H1N1, with an estimated human mort
7、ality rate of only 2.5 to 10%, but with much higher infectivity , produced an estimated 50 million deaths in the 1918 pandemic. 在美国每年大约3.6万人死于普通流感，20万人因此住院治疗。墨西哥调确诊病例死亡率1.2%。世界卫生组织流感专家进藤奈邦子12日说，大多数病情较轻的甲型H1N1流感患者无需服用抗病毒药即可痊愈，因此没必要给所有感染者服用抗病毒药，这些药最好留给那些容易发病严重甚至死亡的高风险人群使用，高风险人群指的是那些本来就患有心血管疾病、糖尿病以及其他慢
8、性病的人或者孕妇等。这些人一旦感染甲型H1N1流感病毒容易发病严重。在墨西哥和美国的甲型H1N1流感确诊病例中，有约10%需要住院治疗，这个比例远高于季节性流感。 4 数学模型 SIR模型用以估计长期传播趋势4.1 模型假设(1)病人数占总人数关于时间的增长率与当时的未受感染人数成正比,比例常数为。k：每个病人每天有效接触的平均人数称为日接触率。(2病人数占总人数治愈率是常数m( 即病人数减少速度与当时的病人数成正比,比例常数为)。m：每天被治愈的病人数占病人总数的比例为常数，称为日治愈率。(3)治愈后的病人(其人数占总人口数 记为)有免疫力, 即治愈后的病人不再被视为未受感染之人。(4) 、
9、与之和是一个常数1。4.2 模型建立 : .4.3 模型求解由第一式得；第三式两边求导：；再将第二式代入上式消去，并利用第一式，最后可得：。利用进行求解：这个微分方程没有解析解，即不能用初等函数表示出、。 利用相平面法，即转而考虑的函数关系。由上面三个方程联立可得: 通过数学软件可解得： 4.4模型分析 4.4.1微分方程分析（一） 当 S ( t ) m / k 时, t 增大时, I ( t ) , 传染病为 蔓延期 ；当 S ( t ) 1/,则开始有：先增加,当s=1/时, 达到最大值：即 ， 然后减小且趋于零,s(t)则单调减小至。4.若 1/,则恒有 ，i(t)单调减小至零,s(t
10、)单调减小至。,可以看出,如果仅当病人比例有一段增长的时期才认为传染病在蔓延,那么是一个阈值，当 (即)时传染病就会蔓延。如果能够减小传染期接触数,即提高阈值使得(即),传染病就不会蔓延(健康者比例的初始值 是一定的,通常可认为 接近1)。我们注意到在中,人们的卫生水平越高,日接触率k越小；医疗水平越高,日治愈率没m越大，于是越小,所以提高卫生水平和医疗水平有助于控制传染病的蔓延。4.4.3 结合模型的假设条件和美国具体情况进行分析从相关报道来看，甲型H1N1流感在美国呈现出蔓延的形势，并且有加快的迹象。根据上面的相平面分析。即现在是属于的情况，即。由假设条件可以知道k的取值在1.4-1.6范
11、围内。现在我们取，则表示，即美国每天平均治愈的人数最多有1.6人，这与美国疾病预防与控制中心所发布的数据不同。美国疾病预防与控制中心表示，患有此次流感的患者中只有大约10%的人需要到医院就诊。这就表明：大约有90%的人不必到医院就诊，只需去医院以外的医疗机构寻求医疗帮助，即可痊愈。如果美国平均每天治愈1.6人的话，那么从发现第一例流感病例（4-23至5-15）以来总共治愈的人数只有人，这与人（截至15日有4714经实验室确认患者，4人死亡）相距甚远。产生这个问题的原因有以下几个：第一：对每个病人每天有效接触的平均人数称为日接触率k的估计值偏小。随着现代交通方式的发展和经济一体化的发展。如果有一
12、人患病，流动到其他的地方。在这个流动的过程中只要感染了k个人，那么这k个人又会感染其他的k*k个人。那么计算病人这一天接触率就应当包括那些间接接触感染的病人，即。第二，美国疾病预防与控制中心所得到的经实验室确认的患者数量具有滞后性。如果一个患者今天住院，但是经过实验室正式确认则往往要经过一段时间。那么这个数据就具有滞后性。第三：在美国并不一定成立。可以把那些身体强壮的、注意个人卫生人群排出在之外，由于，所以。相关资料也指出：我们看到在SIR模型中=k/m是一个重要参数，实际上m、k很难估计，而当一次传染病结束后，可以获得和，根据相关的公式得 。当同样的传染病到来时，如果估计量m、k没有多大的变化，那么就可以用上面得到的分析这次传染病的蔓延过程。尽管有以上技术方面数据获取的缺陷，SIR在理论上具有无可比拟的优势。它能够适用于解释已经发生的传染病的解释，预测流感传染的趋势，也就是在抽象意义上具有绝对的优势。如果要在微观上预测传染病的传播趋势则要采用其他的模型时间序列计量模型。5.模型二：时间序列计量模型用以估计短期患者人数符号说明：表示第期确诊的患者数量。假设:：只与有关系（）对时间序列，进行平稳性分析：152831142054066479181099141101601122612279134031464215896161639172254
