关于病毒扩散与传播的控制模型
摘 要
本文主要以某种传染病疫情为例,利用微分方程来研究和讨论一般传染性病
毒扩散与传播的的控制模型。
模型一:针对问题一,本文在考虑人群分为五类:确诊患者、疑似患者、治
愈者、死亡和正常人前提下,建立传染性病毒扩散与传播的控制模型,即控制后
模型。
针对问题二,在模型一的基础上且满足问题中的四个条件要求对问题二进行
模拟求解,患者人数随时间变化的曲线见图 1 明确标识图中的一些特殊点的具体
数据,并且分析结果的合理性。
针对问题三,在问题二的基础上,对问题二的条件 4 作微调之后进行模拟求
解具体求解见图 2。
针对问题四,在问题二的基础上,对问题二的条件 3 作修改之后进行模拟求
解,模拟结果见图 3。
针对问题五,在问题二的基础上,仅对问题二的条件 1 改动之后进行模拟求
解,模拟结果见图 4。
针对问题六,在模型一的大前提下,对问题二的四个条件作单变量调整并得
出模拟结果图,通过对各个结果的分析,本文得到参数对计算结果的敏感性。
针对问题七,依据如上数据,模型的求解和参数的敏感性分析,本文结合实
际情况给政府部门一个建议报告。
关键词:退出率 过渡期 稳定期 控制末期 微分方程 自由带菌者
一、 问题重述
已知某种不完全确知的具有传染性病毒的潜伏期为 d1~d2 天,病患者的治愈
时间为 d3 天。该病毒可通过直接接触、口腔飞沫进行传播、扩散,该人群的人
均每天接触人数为 r。为了控制病毒的扩散与传播将该人群分为五类:确诊患者、
疑似患者、治愈者、死亡和正常人,可控制参数是隔离措施强度 p(潜伏期内的
患者被隔离的百分数)。 要求:
1.在合理的假设下试建立该病毒扩散与传播的控制模型;
2. 利用你所建立的模型针对如下数据进行模拟
条件 1:d1=1, d2=11, d3=30, r=10,
条件 2:已经知道的初始发病人数为 890、疑似患者为 2000
条件 3:隔离措施强度 p=60%
条件 4:患者 2 天后入院治疗,疑似患者 2 天后被隔离,试给出患者人数随时间
变化的曲线图,并明确标识图中的一些特殊点的具体数据,分析结果的合理性。
3.若将 2 中的条件 4 改为条件:患者 1.5 天后入院治疗,疑似患者 1.5 天后被
隔离,模拟结果有何变化?
4.若仅将 2 中的条件 3 改为条件:隔离措施强度 p=40%,模拟结果有何变化?
5.若仅将 2 中的条件 1 改为条件:d1=1, d2=11, d3=30, r=250,模拟结果有
何变化?
6.分析问题中的参数对计算结果的敏感性。
7.针对如上数据给政府部门写一个不超过 400 字的建议报告。
二、问题分析
关于问题一,为了建立该病毒扩散与传播的控制模型,我们将人群分为易感
人群 ( )S t 、疑似患者 ( )E t 、确诊患者 ( )
I t 、治愈者 ( )R t 和死亡者。结合问题的要
求,我们采用微分方程的思想建立含有两个时滞的 SEIRS 模型,其中,易感人群
一天内人口的出入包括人口的输入、自然死亡的人数、转入疑似患者的人数、转
到确诊患者的人数、疑似患者转入的人数和治愈者过了免疫期的人数;疑似患者
包括由易感人群转入的人数、转入易感人群的人数(经过 2d 天后确定没有患病
离开隔离室的人数)、转入确诊患者的人数和自然死亡的人数;确诊患者包括由
易感人群转入的人数、疑似患者转入的人数、治愈者、因病死亡和自然死亡的人
数;治愈者包括由确诊患者转入的人数、转入易感人群的人数。问题二、问题三、
问题四和问题五都是在问题一的基础上进行模型的变更求解的。
三、模型假设
1)将人群分为五类:健康者(易受感染者),用 ( )S t 表示;确诊病人(己被隔
离的病人),用 ( )
l t 表示;疑似病人(当日没有发病的隔离者),用 ( )D t 表示;
退出者(包括“死亡者”和“治愈者”),用 ( )R t 表示;自由带菌者(没有收治
的传染病患者),用 ( )F t 表示。
2)潜伏期为一确定常数,且由于潜伏期时间较隔离、发病到治愈或死亡时间,
因此,这里不再考虑潜伏期。
3)为了计算方便,假设自由带菌者一经发现就被确定为确诊病人,而不必隔离
观察;而被自由带菌者感染的人被发现必须隔离观察,即归为疑似病人;未被发
现则归为自由带菌者。
4)隔离人群完全断绝与外界接触,不再具有传染性。
5)根据目前的医学调查资料,传染病康复者尚未有复发情况,因为对于一个传
染病康复者而言,一定时期内身边的人会主动远离他。因此,可以假设一个传染
病康复者二度感染传染病的概率为 0,这些人既不是健康者(易受感染者),也
不是病人(已感染者),他们已经退出传染体系。另外,还有一部分确诊病人由
于传染病或其他并发症、自然情况等因素而死亡,他们同样也退出了传染体系,
因此,模型中不予考虑。
6)人群中自由带菌者和健康者均匀混合。
四、符号约定
在这一阶段,每日自由带菌者都会感染一定数量的健康人群。他们中的一些
人被医院收治,成为新增疑似患者;另一些成为新增自由带菌者。而新增确诊人
数的组成则是从疑似转化而来的患者以及新增被医院收治的自由带菌者。基于上
述分析,相关参数设定如下:
1l :疑似病人中每日被确诊人数占疑似病人总数比例
2l :疑似病人中每日被排除的人数占疑似病人总数比例
:每日新增被感染人群中被隔离的概率
:每日自由带菌者被确诊为病人的概率
M :平均每日每隔离一个疑似患者同时会隔离健康者的人数
( )S t :健康者(易受感染者)
( )
l t :确诊病人(己被隔离的病人)
( )D t :疑似病人(当日没有发病的隔离者)
( )R t :退出者(包括“死亡者”和“治愈者”)
( )F t :自由带菌者(没有收治的传染病患者)
五、模型的建立与求解
5.1 模型Ⅰ
开始控制后每日现有确诊病例数目和疑似病例数目以及自由带菌者数目随
时间 t(单位:天)的关系是:
)
(1
q
F
l D
1
1
1
t
t
(1
)
(1
l
l
F k M
1
1
2
1
t
)
(1
(1
)
F k
1
t
I
t
1
t
D D
t
t
F F
t
t
I
1
)
其中, 0
I D F 为初值,即为刚开始控制时的各类人群的数目。
,
,
0
0
1)
q 每日死亡人数 每日治愈人数
前一天的确诊人数
+
2) 1l 每日疑似病人转为确诊病人数
疑似病人总数
3) 2l 每日被排除的疑似病人数
疑似病人总数
4) ,
, M 都与政府的控制力度和公众的防范意识有关,在一定的社会环境下,
可以认为是常量,可以通过相关数据统计而得。
根据官方提供的数据,对过渡期、平稳期和控制末期分别给出各个参数的估
计取值:
过渡期:
q
0.01,
l
1
0.025,
l
2
0.03,
0.4,
0.3,
k
0.338,
M
5
平稳期:
q
0.03,
l
1
0.005,
l
2
0.03,
0.8,
0.2,
k
0.338,
M
4
控制末期:
q
0.01,
l
1
0.001,
l
2
0.2,
08,
0.2,
k
0.338,
M
4
5.2 问题二的求解
在模型一的基础上,我们代入数据
1,
30,
10,
11,
p
d
r
d
3
2
0.6,
( ) 890,
I t
0
E t
( ) 2000,
0
m
经
2
d
1
查 相 关 文 献 可 知 :
N t
6
( ) 2 10 ,
75.5%,
9.3%,
90.7%
运 用
MATLAB 软件求解(程序见附件 1),得到患者人数随时间变化的曲线图为:
图 1 患者人数随时间变化
由图可知,首先,可能由于该传染性病毒是不完全确知的,所以在该病毒发
生的初期,即前 20 天左右,病毒的传播速度更快。其次,随着对该传染病认识
的不断深入,采取的措施及医疗水平都有所提高,因而该病毒又得到短暂的控制,
另外 30 天后部分病人治愈出院,使的患者有所减少。然后,随着染病人数的递
增,每天被传播的人数也在递增,并在第 70 天达到高峰,另外考虑到治愈期比
较长,所以患病者又在一段时间内急剧增长(35~70 天)。最后,在政府的相关
隔离措施的基础上,病情终于得到控制,所以在 92 天该传染病最终消失。
5.3 问题三的求解
在问题二的基础上,我们把 2m 改为 1.5
m
,得到患者人数随时间变化
的曲线图为:
图 2 患者人数随时间的变化
与问题二进行对比可知,患者越早入院治疗,疑似患者越早被隔离,病情越
早得到控制,在第 86 天的时候病情就已经得到了控制。
5.4 问题四的求解
在问题二的基础上,我们把隔离措施强度 60%
p
改为 40%
p
,得到患者
人数随时间变化的曲线图为:
图 3 患者人数随时间的变化
由上图我们发现,在措施强度 40%的时候,该传染病不能被控制,隔离的
措施会对该传染病造成一定的影响,使它周期变化,但是却不能使这个传染病消
除。
5.5 问题五的求解:
在问题二的基础上,把 10
r 改为 250
r ,得到患者人数随时间变化的曲线图
为:
图 4 患者人数随时间的变化
由上图我们发现,若仅将 2 中的条件 1 改为条件:d1=1, d2=11, d3=30,
r=250,,该传染病不能被控制,隔离的措施会对该传染病造成一定的影响,使它
产生持续上升的趋势,不能使这个传染病消除。
六、传染病建议报告
随着社会的进步、卫生设施和医疗水平的改善、人类文明的不断进步,虽然
诸如霍乱、天花等曾经肆虐全球的传染病已经得到了有效的控制,但是在世界的
某些地区,特别是在较贫困的发展中国家,很可能出现传染病流行的情况。因此,
找出一个能够预测并为预防和控制提供足够可靠信息的方法是人们所期望的,建
立传染病数学模型可以预测、预防以及控制传染病的流行,阻止其蔓延。
传染病模型是根据传染病传播的一般机理,利用以观测到的数据,对传染病
的传播规律进行模拟,从而找出他的传播方式和途径,预测在以后较短一段时间
内传染病的传播程度,为有关部门提供预防这种传染病的方案,现将有关建议提
供如下:
开展我国新发传染病现状调查;建立国家病原微生物和新发传染病实验室;
建立以实验室为基础的传染病监测系统;加强新发传染病的诊断技术研究;开展
新发传染病的发生机理和预警技术研究;积极引进在我国还没有出现、但在临近
国家已经出现的新发传染病的病原学诊断技术和方法;开展有关新发传染病的群
众性宣传教育工作,提高公众对新发传染病的认识及如何预防的措施和方法;建
立新发传染病的专业疾病控制队伍;开展新发传染病的培训工作,使基层卫生防