SARS传播地数学模型和对经济地影响文档格式.docx

1、两者之间有实质的区别，文中的说法显然不妥。从预测思想来看，该模型是借鉴先发地区、的有关参数对的疫情进行预测的。由于、的疫情和控制都在之前，已经过了高峰期，到5月8日为止每日新增病例已降至10来例，基本处于后期控制阶段。而当时的疫情刚过了高峰期，正处于社会剧烈调整时期，数据较为凌乱，略有下降趋势，但不明显。可见在当时，采取这种借鉴是无奈之举。但是由于城市之间的政策，风俗习惯等不同，城市之间的可比性不强，借鉴存在很大的局限性。如在，由于对传播机制认识不足，中途又出现高度感染的特殊情况。另外使用借鉴法无法对首发城市进行预测。三、 传播模型（一）问题的分析在SARS爆发的初期，由于潜伏期的存在，人们对

2、病毒传播速度和危害程度的认识不够，未能及时识别这一传染病的存在。但当病患数不断增加，政府开始采取应对措施对其进行控制，同时社会舆论加大宣传力度，人们的警觉性提高，病毒的传播速度下降。因此，我们通常把传染病的传播模式近似分为两个阶段：第一、自由传播阶段（即控前阶段）：在采取切实有效的控制措施之前的一段时间。第二、控后阶段：介入人为因素之后的一段时间。由于SARS的传播涉及的因素很多，如潜伏期、人群的迁入迁出，感病者的数量、易感者的数量、传染率和治愈率的大小等。而且在以上因素中，潜伏期的大小、传染率和治愈率的大小因人而易，具有一定的随机性。不可能一开始就把所有的因素全部考虑在建立模型。对此，我们将

3、作出相应的假设进行简化。分析附件2所给出的数据，发现6月1日至15日，已确症病例累计数为2522人，其中夹杂3天累计数为2523人。但6月16日后累计数降至2521人，认为累计数的减少可能是误诊引起的。由于误诊的可能性很低，故在这里忽略不计。（二）基本假设1、 国家卫生部提供的疫情统计真实可信。（误诊数仅为1，可忽略不计）。2、 由于非典的主要传播途径是近距离接触，通过受感染者咳嗽或打喷嚏时产生的飞沫传播，这里将所有传播途径都视为与病源的直接接触。3、 不考虑出生与自然死亡的过程和人群的迁入迁出（或认为迁入和迁出基本平衡），认为疾病传播期间所考察地区的总人数为常数。4、 根据国家卫生部资料可知

4、处于潜伏期的SARS病人不具有传染性。5、 目前尚不清楚康复患者是否具有免疫力，但据国家卫生部资料可知康复后的病患无一例复发。故假设康复患者退出传染系统。6、 根据资料显示，SARS病毒的潜伏期一般为27天，平均约为5天。（这一条件将在后期的模型中有所改动）（三）常用基本模型目前常用的传染病模型，通常将传染病流行围的人群分为三类： S类：易感者，指未得病者，但与感病者接触后容易受到感染。I类：感病者，指染上传染病的人。R类：移出者，指因患病而被隔离，或因病愈而具有免疫力的人，他们即非感病者，也非易感者，实际上他们已经退出了传染病系统。并通过三类之间的互相转化关系建立微分方程组进行求解： （1）

5、变量和符号说明： 传染率：每个病人平均每天有效接触（足以使被解除者感染）的人数。 退出率：单位时间治愈和死亡人数占感病者人数的百分数。S（t）易感人群的总数。I（t）感病者总数。R（t）退出者总数。N一个城市总人口数。观察附件二中给出的数据，我们发现截至6月23日，感病者累计为2521人，远远小于城市的总人口数150万人，故认为感病者和退出者对易感人群的总数影响不大，易感者总人数I为一常数。原方程变形为： （2）注意到退出者不是我们研究的围，故方程组（2）实际上是一个常微分方程 （3）其中，不难用分离变量法解出： （4）其中I0为初始值。根据以上分析我们可以看出，常微分方程的传染病模型只适用于

6、病例数与总人口数具有可比性的情况。当病例数远小于总人口数时，常微分方程模型的实质与附件1的模型相同，感病人数将随时间以指数增长。考虑这一特点，我们用计算机跟踪病毒的个体传播情况，建立了模拟模型。（四）计算机模拟模型：在该模型中，我们将传染系统中的人分为五类：自由携带者（）身上携带病毒并均匀散布在人群中的患者，根据基本假设自由携带者在潜伏期不具有传染力， 日增患者（）每天被医疗部门发现并加以隔离的感病者 被隔离者（）因曾与自由携带者接触而被怀疑携带SARS病毒的人 有效接触者（）每日与自由携带者接触并感染上病毒的人 无效接触者（）每日与自由携带者接触但未染上病毒的人 并作出如下假设：1、由于传染

7、性SARS最初（12天）的症状通常为发热（），发热通常为高热1。症状明显，易于辨认，故可认为自由携带者发病后当天或第二天就立即入院治疗，入院后不会再参与疾病的传播。2、根据实际情况，假设SARS病人被发现的三天，有关部门将采取措施，将部分与病源有效接触者隔离，这部分人即使发病后也不会参与疾病的传播。3、与病源有效接触者必然发病。根据基本假设，潜伏期一般为2至7天，这里取为5天。（这一假设在改进模型中有进一步的讨论。另外，对模拟模型中出现的符号变量说明如下： 有效接触率，表示一个自由携带者平均每天有效接触的人数。 无效接触率，表示一个自由携带者平均每天无效接触的人数。与自由携带者接触后（包括有效

8、接触和无效接触）的人群中可以控制的人数所占的百分比。模拟模型中个体传播情况如图1所示：自由携带者15天处于潜伏期，不具有传染能力；5天后发病，发病后每天有效接触人，2天后（第7天）被隔离，在隔离前每天无效接触人。与病源接触后的可控人群（占接触者总人数的）在3天后被视为疑似病人。疑似病人中的有效接触者在接触病源的第7天被发现确认为日增患者，而有效接触者中其他人作为自由携带者留在人群中，继续这之前的个体传播。图1：个体传播示意图用数学模型描述各个变量之间的关系如下： （5）由于在4月20日才开始建立每日疫情报道制度，故认为政府采取严格的隔离措施开始于4月20日，以这一天为分界线，之前属于自由传播阶

9、段。根据这一模型，用计算机模拟5月10日之前SARS的传播情况，并对5月10日以后的传播情况进行预测。图2：5月10日以前数据拟合图图3：5月10日以后的预测曲线通过图2两条曲线的拟合，得到控前的有效接触率（表征病毒的传染力）k=1.351,可控率（表征政府的控制力度）；控后k1=0.8， =0.7。根据这两个参数作出5月10日的预测曲线见图3。，根据预测，将在第97天（6月下旬）实现零增长，累计病例数2448人。经过分析，以上确定型的模拟模型存在以下两点问题：第一、SARS病毒的潜伏期一般是27天，模型将潜伏期确定为5天，从感染到被发现的时间确定为7天，可以明显地看到以7天为传播周期的曲线变

10、动，而在实际曲线中，虽然数据有上下波动的趋势，但周期性并不明显。第二、在采用隔离政策时，模型假定与病源接触的人群以固定的比例受到控制，这种假设加大了人为主观因素的影响。基于以上两点，在拟合5月10日之前的数据时，得到的有效接触率与实际统计数据有所偏差。这种偏差降低了模型的可信度。基于此，我们查找有关统计数据，为参数的确定寻求医学上的支持，并以随机模拟取代完全确定性的模拟，对原模型进行改进，建立了随机模拟模型。（五）改进后的随机模拟模型1、根据上文的分析，在原有模型的基础上作出如下改进假设：假设一：假设潜伏期的长短服从以5天为期望值，2天为方差的正态分布，即认为自由携带者将以68%的概率在37天

11、发病，以95%的概率在19天发病。假设二：假设与病源接触的人群每天均以一定的比例被隔离。该比例服从以为期望值，0.05为方差的正态分布。2、参数的确定：控前k15月23日，美国科学杂志及其上发表了两篇有关SARS的流行病学论文及一篇长篇评论2，研究确定SARS病患的基本传播率（即文中的有效接触率）为3。控前表示有效接触者中在发病前未能受到有效隔离的人所占的比例。由于在SARS传播前期，对疾病未能作良好的监控，而在潜伏期病患无法得知自己的真实情况。所以在这一阶段不妨设，即几乎没有病患在潜伏期被隔离。这与控前的实际情况是吻合的。控后k1控后的有效接触率实际上受多种因素综合影响。比如人们不愿外出，减

12、少了与病源的接触机会，更注意卫生等，这些因素均使有效接触率明显降低。英国和的一个研究小组的研究表明，传染率（即本文的有效接触率）控后仍有较长一段时间在1左右徘徊。为简化模型，我们取k1=1。控后控后的隔离措施被细化为几个等级，措施的级别越高，控制力度越大，社会付出代价越高，难度越大3。经研究发现，采取隔离措施后短期使未隔离的病患在总自由携带者中占40%90%的比例3，即在0.4在0.9之间。由于各地的控制力度不同，对应的不同，对4月20日以前的数据进行拟合后得到当时，拟合较好。从中可以看出的措施是比较严厉的。3、预测结果由于缺乏4月20日之前的数据，不妨设初始6天的病例数分别为1，2，3，4，

13、5，6（之后将对初始值的合理性进行讨论）。为了较好地拟合数据，采用3月15日作为病毒传播的开始时间，认为在病毒传染前期，主要以输入型传播为主，传播机理不明，输入病例数量无法确定。采用以上数据和相关参数，通过计算机编程，产生正态分布的随机数，并对传染情况进行500次模拟，取其平均数得到预测结果见表1。表1：实际值与预测值的比较日期5.115.125.135.145.155.165.175.18实际值3839432318171514预测值43.40373026205.195.205.215.225.235.245.255.2637129255161311105.275.285.295.305.31

14、6.16.26.3821646.46.56.66.76.86.96.106.11图4：改进后的拟合曲线图5：改进后的预测曲线图6：提前5天采取措施的预测曲线若卫生部门提前5天采取严格的隔离措施，预测将会出现图4的情况：累计病例总数为1958人，到第100天就可实现零增长。若滞后5天才采取严格的隔离措施，预测曲线如图5所示，疫情将延迟至第121天才能基本消除，届时累计病例总数将达到4487人。若卫生部门的隔离措施滞后15天才得到全面实施，将出现图6的情况，到第139天疫情才能基本消除，累计病例数将增至10208人,比原来增长了约8000多名患者。 图8：滞后15天的预测曲线若社会及政府等相关部门

15、加大宣传力度，人们警觉性提高，使有效接触率降低，设k=0.8，得到预测曲线图7，预计到第92天疫情解除，累计病例2395人。若卫生部门实施严格隔离措施，加大对与自由携带者相接触的人群的控制力，设可控率提高到。得到预测曲线图8，预计到第80天疫情解除，累计病例2545人。表2：改变有关参数后的预测值改变相关参数累计病例数实现零增长的时间与预测病例数的差值与预测实现零增长的日期相比最初预测2816109提前5天1958100-858-9滞后5天44871211671+12滞后15天102081397392+30加大控制力度254580-271-29加大宣传力度239592-421-18以上预测结果

16、表明，加大宣传力度，可以明显减少感染总数；而加大控制力度，可以显著减少疫情延续的时间。根据以上预测结果，对SARS疫情的控制提出以下建议：1、应尽早采取有效的控制措施。提早采取严格的隔离措施，将会使疫情提早得到控制或解除。2、隔离和宣传并重。一方面加大控制力度，提高接触病源的人群的可控率，提前结束疫情；另一方面加大宣传力度，提高全民健康意识，切断传染途径，从而减少自由携带者的有效接触人数。这样即使防措施晚了一些，出现了高峰期，仍然能够在一定的时间解除疫情。四、 SARS影响经济模型从题中数据可以看出SARS对旅游行业造成了严重的影响，影响最大的是五月，海外旅游人数下降到了1.78万人，比去年下

17、降了27.22万人，占99.34%。SARS过后海外旅游人数有了较快的回升，在近期有望达到疫前水平。从历年各月旅游人数的数据可以看出，近年来旅游行业保持着良好势头，中国旅游业的进步是有目共睹的。如果没有SARS，2003年旅游业仍然会保持持续增长。SARS只是一个突发事件，不会对旅游业造成长期的影响，不过在近几个月还有一个恢复的过程。因此认为对2003年今后几月旅游人数的影响因素有：1、旅游业在的发展动力。如不断增加旅游业投资，为2008年奥运会做的建设等。2、SARS流行期间对旅游业的影响。两种因素的影响特点有：第一种是旅游业进步的在本质因素，具有根本性、长期性。第二种是突发事件，是偶然因素

18、，具有暂时性，作用效果看有突出重要性。分别分析两种力量对2003年旅游业的作用，两种力量的效果是独立影响的，最后把两种效果叠加起来。（二）模型的建立和求解1、不考虑SARS对2003年各月旅游人数的预测从19972002年的数据可以看出，各月人数的分布有明显的规律性，见图9规律一：历年各月都在不断增长，从1999年到2002年4年的平均增幅为27.2万人。增长过程见图10。我们做了每月与历年该月的相关性检验，发现各月与前两年的该月份显著相关。规律二：每年中各月的大小关系基本不变，表现为图9中的同时上升，同时下降，且各月的增幅比例大致相同。我们也作了12个月中每个月与相邻几个月的相关性检验，发现

19、有一定相关性。如果有足够的数据，我们希望得出这样的相关关系：F（i,j）=F（i,j）表示第i年，j月的人数；反映历年该月的影响；反映相邻月份的影响。但是目前得到的数据太少，不能统计出相关系数。下面从数据的特点预测没有SARS影响的各月旅游人数。近几年来，旅游业的不断进步，必然吸引越来越多的国外游客。不过各年的增幅有一些波动，其实有点波动是很正常的，因为事物都是在曲折中前进，螺旋式上升的。从图10中可见，2001到2002年（4点到5点）的增长轨迹与1998-1999-2000年（1-2-3点）的增长轨迹很相似。从1998年开始经过两年快速增长期后，进入一个低速增长期，两个时期的总长约为3年。

20、2001年开始的快速增长，很可能又要进入下一个快速增长期，可推测2003年在2002年的基础上会增长40万人。这40万总增量分配到12个月中，分配比例按表3中的统计比例。得到2003年各月的预计值并加入图9中比较。图9图10：近几年旅游人数的年增长数表3：各月份增幅比较月份月2月3月4月5月6月平均增幅1.0252.171.8252.252.3752.4增幅比例0.0412140.0872540.0733820.090470.0954970.0965027月8月9月10月11月12月2.052.2252.52.0252.2751.750.0824290.0894650.1005230.0814

21、230.0914760.070366表4：2003年预测值1月,2月,3月,4月,5月,15.431,33.36466,26.18203,32.69976,33.01086,31.453087月,8月,9月,10月,11月,29.462,32.2,35.62195,36.01978,33.04198,25.855372、考虑SARS对旅游人数的影响用2003年的预测值减去2003年8月以前的实际值，可看出SARS减少的旅游人数，见图11。其中对2003年2月作了处理，2月时没有SARS，而且1月、3月的人数比2002年都有增加，符合增长规律。则2月份旅游人数的降低（近一半），由其他突发因素引起

22、。用2002年2月的数据替换。考虑到SARS的蔓延是每天变化的，把每月的旅游人数离散到每天。把每月人数作平均，作三次样条插值，得到光滑的曲线，曲线上的点反映每天访问人数。差值曲线反映了SARS没天减少的访问人数。我们发现每天减少的访问人数有个峰值，这个峰值出现在5月，6月之间，到8月时已经过。说明SARS对5、6月影响最大，到8月已逐渐减弱。关于SARS的指标有多个，如果有个指标能反映对旅游者心理的影响，则SARS对旅游业影响减弱时这个指标也下降。经过对众多指标（死亡率，死亡人数，治愈人数，治愈率，疑似患者数，累计患者数）的筛选，选出能较好地反映对旅游人数的影响的指标每日现有患者数，并和旅游人

23、数变化曲线画在一起观察，如图12。图11：曲线比较图图12分别选取两条曲线中的两个序列作相关分析，得到每日现有患者数对旅游人数的影响有个时延，这个值为8天，也正反映了两条曲线的峰值相差10天左右。到2003年6月23日，患者还有53个，2个疑似患者，可以确信SARS即将消失，因为已经没有疑似患者，一段时期多个危险指标没有增长。那么SARS对旅游的影响也即将过去。目前患者数和日访问人数都一直保持恒定的下降速度。可从数据得出日访问人数以0.019的斜率降低，10月25日以后SARS已没有影响了。累加SARS对每天访问人数的减少量，得出9月比预计值减少11.5944万，10月比预计值减少2.95万。

24、最后我们综合两种因素后得到预测值见表5，图13。表5：1月15.429.723.511.61.782.618.816.224.0275533.0688833.04198图13以上只分析了旅游业的在增长因素和SARS的打击作用。在SARS时期没有旅游的人，可能累积到SARS以后旅游，这会使实际旅游人数比预测值大。可见以上只是保守的预测，不过都是很乐观的预测。五、 短文传染病模型的研究任重道远传染病是危及人类身体健康的重要因素之一，长期以来一直受到世界各国的关注。随着社会的发展、卫生设施和医疗水平的改善和人类文明的不断发展，诸如霍乱、天花等曾经肆虐全球的传染病已经得到了有效的控制。但是在世界的某些地区，由于战争对抗、人为失误或病毒变异等因素，还不时出现传染病流行的情况，如82年英国的全国性沙门氏菌病暴发流行，以及最近迅速蔓延全球的SARS等。目前很多学