传染病模型数学建模论文

甲型H1N1流感传播模型研究

摘要

本文采用了SIR模型对的甲型h1n1流感病毒的传播规律进行了研究和预测，文章收集了美国地区的甲流实验室确认病例数量的数据，对模型进行了验证，并提出了如何降低流感在人群中发病率的俩种可靠方法。

一、问题重述

1

邹健：甲型H1N1流感传播模型研究

近年来由墨西哥发端的甲型h1n1型流感（又称猪流感）正成为人们关注的焦点，通过相关网站获得数据，建立一个模型对甲型h1n1流感的走势进行预测。

二、问题分析

甲型h1n1流感的传播是一道传染病问题。在数学建模领域已经有很多关于这方面的研究，其中SIR模型是比较完整的模型。SIR模型通过建立微分方程组，按照一般的传播机理建立集中模型。本文选取美国地区的甲流实验室确认病例数量，建立SIR模型，对甲型h1n1流感的传播规律进行预测。

美国甲型H1N1流感实验室确认病例数量：

时间 确诊（包括死亡病例） 死亡（累计） 4月23日 5 0 4月24日 8 0 4月25日 11 0 4月26日 20 0 4月27日 40 0 4月28日 0 4月29日 91 0 4月30日 109 1 5月1日 141 1 5月2日 160 1 5月3日 226 1 5月4日 279 1 5月5日 403 1 5月6日 2 2 5月7日 6 2 5月8日 1639 2 5月9日 22 2 5月10日 2532 3 5月11日 2600 3 5月12日 3009 3 5月13日 3352 4 5月14日 4298 4 5月15日 4714 4 三、建立模型

（一）、不考虑潜伏期的数学模型

1、模型假设

（1）、在甲型H1N1流感传播期内，美国境内的总人数为N亿不变，既不考虑生死，也不

2

数学建模

考虑迁移，人群分为易感染者S，发病人群I和退出人群R(括死亡者和治愈者)四类，时刻t内这三类人在总人数中所占比例分别为s(t)、i(t)、r(t)。 （2）、i(t)关于时间的增长率与s(t)成正比，比例常数为。

病人的数量减少速度与当时的病人总人数成正比，比例常数为。治愈

的病人具有了免疫力，即治愈后不再会成为二次患者。 （3）、s(t)、r(t)、i(t)之和是一个常数1。

2、模型构成

易感者和发病者有效接触后成为发病者者。设每个发病者平均每天有效接触的易感者数为S(t)，NI(t)个发病者平均每天能使S(t)NI(t)个易感者成为病毒潜伏者。所以有：

dS(t)S(t)I(t) （1） dt单位时间内退出者的变化等于发病人群的减少，即

dR(t)I(t) （2） dt发病人群的变化等于易感人群转入的数量，即

dI(t)S(t)I(t)I(t) （3） dt 记初始时刻的健康者和病人的比例分别为S0、R0（不妨设R0=0）。

3、模型求解

方程组（1）、（2）、（3）无法求出解析解，我们定义一个新的变量 /，于是可以求出方程的解为：

i(s0i0)s1lns （4） s0 下面分析s(t)、i(t)、r(t)的变化情况：

a、不论初始条件S0、R0如何，病人最终将消失，即i0。 b、最终未被感染者的健康者的比例是s，是方程

(s0i0)s1lns0在(0,1/)内的根。 s0C、若s01/，则开始有：i(t)先增加。当s01/时，i(t)达到最大值，然后i(t)减小且趋于零，s(t)则单调减小至s。

d、若s01/，则i(t)单调减小至5，s(t)则单调减小至s。

我们发现人们的卫生水平越高，日接触率越小；医疗水平越高，日治愈率越高，于是越

3

邹健：甲型H1N1流感传播模型研究

小，所以提高卫生水平和医疗水平有利于传染病的蔓延。

结合美国的具体情况和假设条件进行分析：

根据所得的数据画出美国患病人数变化曲线和治愈人数变化曲线：

4

数学建模

根据图形来看，甲型h1n1流感在美国呈现出蔓延的形式，即现在属于s01/ 的情况，即

/1/s0。由假设条件可知的取值范围在1.4~1.6之间。现在我们取

=1.6，则表示/(1/s0)1.6，即美国每天平均治愈的人数最多为1.6人，这与美国

疾病预防与控制中心所发布的数据不同。如果美国平均每天治愈1.6个人的话，那么从4月23日期，治愈的总人数为1.6*2336.8人，这与实际的情况相差甚远。产生这个问题的原因有以下几个方面：

第一：对每个病人每天有效接触的平均人数估计值偏小。不是简单的成正比关系，应该是成多次方关系，甚至是指数关系。

第二：美国疾病预防与控制中心所得到的数据具有滞后性。

第三：在美国s00不一定成立。可以把那些身体强壮的、注意自己个人卫生的人排除在外。

（二）、考虑潜伏期的数学模型

1、模型假设

（1）、在甲型H1N1流感传播期内，美国境内的总人数为N亿不变，既不考虑生死，也不考虑迁移，人群分为易感染者S，病毒潜伏人群E,发病人群I和退出人群R(括死亡者和治

5

邹健：甲型H1N1流感传播模型研究

愈者)四类，时刻t内这三类人在总人数中所占比例分别为st、et、it、rt。 （2）、每个病人每天有效接触的平均人数为，称为日接触率，当已感染者与易感染者有效接触时，使易感染者变为病毒潜伏人群，病毒潜伏人群过一段时间再转换成发病人群，发病人群被治愈。

2、模型构成

易感者和发病者有效接触后成为病毒潜伏者。设每个发病者平均每天有效接触的易感者数为(t)S(t)，NI(t)个发病者平均每天能使(t)S(t)NI(t)个易感者成为病毒潜伏者。

dS(t)化简得： (t)S(t)NI(t)dt dS(t)(t)S(t)I(t) dt病毒潜伏人群的变化等于易感人群转入数量减去转化为发病人群的数量，即

dE(t)(t)S(t)I(t)(t)E(t) dt所以有N其中(t)表示潜伏期日发病率，即每个潜伏者平均有效发病的人数。 单位时间内退出者的变化等于发病人群的减少，即

dR(t)(t)I(t) dt其中(t)表示日退出率，即每个病人平均有效病情结束的人数。 发病人群的变化等于潜伏人群转入的数量，即

dI(t)(t)E(t)(t)I(t) dts(t)i(t)r(t)1

初始时刻易感染者，已感染者与病愈免疫者的比例分别是

s0(s00),i0(i00),r00

3、模型求解

由于潜伏期的人群数量不能确定，所以可视为是易感人群的一部分，因此求解过程跟忽略潜伏期的一样。

四、模型的改进

就如何确定日接触率的值。就如何确定日接触率可以进行改进，根据以前的流感疫

6

数学建模

情治愈率，加权平均得到值，而不是简单的是一个正比关系。病毒在人群中的传播刚开始阶段一个有一个爆发阶段，该阶段的日接触率很大，可设为是一个冲激变量。

参考文献：

[1]姜启源 谢金星 叶俊 数学建模（第四版）高等教育出版社 [2]数据来源：美国疾病预防控制中心

7