SARS流行病模型及其对未来走势的预测

1/1SARS流行病模型及其对未来走势的预测SARS流行病模型及其对未来走势的预测摘要SARS是21世纪第一个在世界范围内传播的疾病，它的爆发和蔓延给我国的经济发展和人民生活带来了巨大影响，因此定量地研究传染病的传播规律具有十分重要的意义。

本文在一系列合理假设的基础上，通过对问题的详细分析，首先对于附件（1）中的模型给予客观公正的评价。

该模型简单明了，具有一定的合理性和实用性，但是又存在很多的缺陷，具有一定的片面性，对影响SARS疫情传播的因素考虑的不够全面。

在此基础上，我们提出了SIR模型，该模型将某一地区的人们看作一个动力系统，再将其划分为易感染者、患病者和移出者，通过建立微分方程组来加以解决。

在此模型中，我们考虑了个体的免疫率、死亡率以及被感染者的康复率等因素。

为了验证模型、算法的正确性和有效性，通过计算机模拟，对北京、内蒙古、广东、河北、山西各地的SARS疫情进行拟和，并给出了我们的预测数据，结果表明我们的模型与实际情况相当的吻合。

另外，我们还分别针对得病后入院时间以及隔离强度的不同，对SARS疫情传播所造成的影响做出估计，其结果表明：

SARS病人1.5天后入院与2天后入院相比，SARS发病总人数可能会减少1500人，SARS疫情得到控制的时间可能会提前1个月；而得病1.5天后入院与2.5天后入院相比，SARS发病总人数可能会减少2400人，SARS疫情得到控制的时间可能会提前1个半月；隔离措施强度60%与50%相比，SARS发病总人数可能会减少700人，SARS疫情得到控制的时间可能会提前半个月；隔离措施强度60%与隔离措施强度40%相比，SARS发病总人数可能会减少1100人，SARS疫情得到控制的时间可能会提前1个月。

一、问题的重述一、问题的重述SARS（SevereAcuteRespiratorySyndrome，严重急性呼吸道综合症,俗称：

非典型性肺炎）是21世纪第一个在世界范围内传播的传染病。

SARS的爆发和蔓延给我国的经济发展和人民生活带来了很大影响，我们从中得到了许多重要的经验和教训，认识到定量地研究传染病的传播规律，为预测和控制传染病蔓延创造条件的重要性。

请对SARS的传播建立数学模型，具体要求如下：

（1）对附件1所提供的一个早期的模型，评价其合理性和实用性。

（2）建立你们自己的模型，说明为什么优于附件1中的模型；特别要说明怎样才能建立一个真正能够预测以及能为预防和控制提供可靠、足够的信息的模型，这样做的困难在哪里？对于卫生部门所采取的措施做出评论，如：

提前或延后5天采取严格的隔离措施，对疫情传播所造成的影响做出估计。

（3）给当地报刊写一篇通俗短文，说明建立传染病数学模型的重要性。

二、已有模型的分析对于附件1中的模型，我们可以很直观地看出，此模型简单明了，没有涉及到太多的变量和参数，读者很容易理解。

模型不但包括了每个病人可以传染他人的期限L，而且还考虑了不同阶段、不同环境下传染概率的变化，并对香港和广东的疫情进行计算和分析，其结果论证了该模型具有一定的合理性。

为了便于分析，其假定参数L为20天，此时该模型能较准确地反映出香港疫情的变化趋势。

用此模型对广东疫情进行测算时也较准确地反映了疫情上升期间的趋势，但在随后的下降过程中反映的数据却与实际情况有一定出入，这可能与后期广东对SARS疫情的预防或控制状况有一定的关系，需要更进一步的分析。

经过对广东及香港两地的疫情的计算和分析，他们又利用该模型来计算和预测北京地区的病例情况，发现其预测结果与实际情况大致相同，且基本上能够反映疫情的走势。

从文中的实际数据和模拟结果的对比情况来看，该模型比较合理地反映了上述三处的SARS疫情，并对后期的实际操作及预防起到了一定的作用。

但是我们看到，该模型还是存在着某些不足。

比如期限L被固定在20天，这就有了一定的局限性，因为随着社会的重视程度及医疗水平的不断提高，还有其它诸多因素的影响，期限L将会缩短，而不是一个恒定不变的值；另外该模型考虑的影响SARS流行病的因素较少，从而导致其未能反映出一些必要的数据，比如病人的康复数及康复率等等。

因此在此基础上，经过合理的假设和必要的分析，我们提出了如下的SIR模型：

三、SIR模型的建立及求解大多数流行病如天花、肝炎、麻疹等治愈后有很强的免疫力，它们已退出传染系统。

严重的急性呼吸道综合症SARS就是21世纪出现的第一个既严重且又易于传播的疾病。

但是它又不同于其他传染病，它引起了社会的极大恐慌，并造成社会运作方式发生相应的变化，因此其本质已经超出了一般传染病的发展规律。

但在政府、各单位以及每个公民的共同努力下，如传染者与其他人接触率的减少、改进医院的控制措施、患者接触的检疫以及公众自动减少的接触等有效措施的实行，必然使非典疫情的发展满足更一般意义的规律。

因此，我们结合SARS流行病的特点，建立SARS流行病动力学的SIR模型[2]。

我们将人群分为易感染者（Susceptible）、患病者（Infectious）、移出者（Removed）。

所谓易感染者是那些健康且与病人接触过就容易感染疾病的人，他们中的一部分可能转化为患病者；患病者是那些能够将疾病传染给易感染者的人，在一定期限内（感染阶段），他们仍属于患病者；移出者则包括已治愈的患病者和死亡者[3]，我们可以通过图1来表示他们之间的转化关系：

图1.SIR模型中各群体间的相互转换图1.SIR模型中各群体间的相互转换各类t时刻的人数分别用tS、tI、tR表示，且t时刻因病死亡人数也属于移出者一类。

若用t表示t时刻人口总规模，则有：

tNtRtEtS假设每个个体都具有因感染SARS而死去的相同的死亡率；单位时间内每个患者传染的人数为kIS（非线性传染率），其中k为大于0的常数；一个易受感染者在与病人进行有效接触后立即被感染且被感染后的康复率；患者病愈后具有一定的免疫能力，其免疫率为h。

我们再假设只考虑因SARS死亡、康复的人，而不考虑人的自然出生与自然死亡人数（因数据不全，且影响小），则NtRtEtS（常数）由此我们得出带有免疫率的SARS流行病的数学模型为：

tRtItStNRhIdtdRIhIkISdtdINBrRSIkISdtdS)()()()()(（1）其中NB为出生率，它是总人口数N的函数，初值0S=0S＞0，00II＞0，00R，Tt,0，T为流行病的流行时间。

若假定出生率等于死亡率，即RISN保持不变，记为0N，则0NNB，从（1）式中消去S可得：

RIdtdRIhRINkIdtdI)()()(02（2）相对于附件（1）中的模型，我们的模型较全面地考虑了影响SARS的因素，包括复发率、感染率、免疫率等。

为了进一步验证模型的合理性和有效性，我们根据世界卫生组织、中国卫生部以及香港卫生署等公布的疫情数据（包括确诊数、死亡数和出院人数），通过计算机，利用MatLab语言[5]（程序详见附录一），对北京、内蒙古、广东、河北、山西各地的SARS疫情进行拟和，并给出了我们的预测数据（详见附录二、三）。

拟和结果表明我们的模型与实际情况相当的吻合，以下是各地区的疫情变化曲线图：

4月20日至7月1日北京SARS累计感染人数曲线图4月20日至7月1日内蒙古SARS累计出院人数曲线图4月25日至7月14日广东SARS累计出院人数曲线图4月25日至7月14日广东SARS累计感染人数曲线图4月21日至7月1日河北SARS累计确诊人数曲线图4月21日至7月1日山西SARS累计死亡人数曲线图为了建立一个真正能够预测以及能为预防和控制提供可靠、足够信息的模型，我们收集了SARS发病较严重的部分地区的资料，其难点是：

(1)数据收集相当困难，而且我们根本无法保证其可靠性;(2)为了使模型的参数尽量符合实际，我们从多方面加以考虑，但其困难在于程序的调试比较麻烦,而且我们无法保证模拟时数据的选择是否合理。

四、评价与分析政府的措施我们选择得病后入院时间变化的模拟分析和隔离措施强度两个参数疫情传播所造成的影响做出估计。

1、得病后入院时间变化的模拟分析在SIE模型其他参数不变的情况下，我们针对得病后入院时间的不同，在得病后入院时间为1.5天、2.0天与2.5天情况下，对模型进行模拟，得到三种不同的曲线，图形如下：

得病后入院时间不同的SARS发病总人数曲线图对上述的三条曲线进行分析，可获得：

SARS病人1.5天后入院与2天后入院相比，SARS发病总人数可能会减少1500人，SARS疫情得到控制的时间可能会提前1个月；而SARS病人1.5天后入院与2.5天后入院相比，SARS发病总人数可能会减少2400人，SARS疫情得到控制的时间可能会提前1个半月。

2、隔离强度不同的模拟分析在SIE模型其他参数不变的情况下，我们针对隔离强度的不同，对隔离强度分别为40%、50%、和60%的情况进行模拟，得到三种不同的曲线，如下图所示：

隔离措施强度不同SARS发病总人数的变化曲线对以上三条曲线进行分析，可获得：

隔离措施强度60%与隔离措施强度50%相比，SARS发病总人数可能会减少700人，SARS疫情得到控制的时间可能会提前半个月；隔离措施强度60%与隔离措施强度40%相比，SARS发病总人数可能会减少1100人，SARS疫情得到控制的时间可能会提前1个月。

五、模型分析与评价我们在建立模型时没有将病例加以细分，因为根据实际情况，病例可分为显性病例和隐性病例（即潜伏患者），我们只是将其视为患病者。

表面上看我们似乎考虑的不够全面，而实际上我们是根据易感染者的感染比率来确定患病人数的，这就说明对于任何的显性患者或是隐性患者，其感染率是通过易感染者的自身免疫力来确定的，跟患者数目没有直接关系。

根据全国非典科技攻关组公布的七大科研进展可知，已出院的非典患者都不具有传染性，密切接触者无隐性感染可能，潜伏期患者传染可能性很小，而SARS则是一种潜伏期较短的传染病。

当然，我们不否认所建立的模型并不是十分完美的，因为还有很多因素都可能对预测结果造成影响，诸如SARS潜伏期天数，患者被发现前的时间长短，住院后治愈时间的长短，隔离措施的强度，还有个人每天接触的人数等等。

另外，我们并不能保证避免地区外的病例输入，而这一因素对预测结果也有很大的影响。

若有可能，可以通过建立系统动力学模型对其进行预测。

六、建立传染病数学模型，势在必行山雨欲来风满楼，中国古代传说中就有能够根据某些迹象预测灾难发生的圣贤之辈，但现实中我们不能再用迷信的方法来欺骗自己，尤其是在划时代的今天，我们要用科学的理论来预防任何灾难的发生。

如今，SARS是21世纪第一个在世界范围内传播的传染病，它的爆发和蔓延给我国的经济发展和人民生活带来了很大影响，我们从中得到了许多重要的经验和教训，如果SARS能够被预测的话，中国就不会陷入动荡不安的社会状况，就不会有数百个生灵离开人间，就不会有那么多不幸的事情发生所以，痛定思痛之后，我们有必要建立数学模型来解决这些难题。

建立传染病数学模型的优点在于它具有科学的理论依据，根据已有的相关信息定量地研究传染病的传播规律，用抽象的数学方法能够计算出不同条件下的各种可能性结果，进而采取有效措施来预防灾害的发生。

然而，传染病数学模型的建立并不是一件轻松容易的事情，它需要对疫情传播所造成的影响做出估计，需要研究人员具备严谨的数学理念，通过计算机编程来预测和控制传染病蔓延，并对各种结果进行误差分析，使传染病产生的危害降到最低限度。

我们有信心，有能力建立一个最佳的传染病数学模型来解决我们身边的任何一类疾病传播问题。

虽然我们的认识是有限的，但我们的不断创新能力是无限的，我们渴望幸福生活的拼搏精神是永远不会磨灭的。

参考文献：

[1]R.M.AndersonandR.M.May,Populationbiologyofinfectious:PartI,Nature280,361-367(1797)[2]吴子牛、杨清红等自控制起效后非典滞留患者数目变化规律的预测2002年5月23[3]/issue4/fertures/disease[4]黄运新一种非线形传染病模型的定性分析数理医药学杂志1999年第12卷第一期第7页1999年出版[5]薛定宇、陈阳泉基于MATLAB/Simulink的系统仿真技术与应用北京清华大学出版社2004、4[6]阻止SARS疫病在农村传播的建议，中国科学院研究生院管理学院SARS管理与控制研究小组《研究简报》第3期附录一：

Matlab程序fd.dydq.dydq2为自编的matlab识别的m文件，依照这三个m文件内的%说明文件输入数据，然后在matlab工作区间运行dydq2即可。

下面为这三个m文件：

fd:function[output1,output2,output3]=fd(t,x,flag,theta)ifnargin4|isempty(theta),%theta=[110.1111],theta=[110.111],endifstrcmp(flag,’’)%output1=[-theta(1)*x(1)*x(3);...%theta(1)*x(1)*x(3)-theta(2)*x(2);...%theta(2)*x(2)-theta(3)*x(3)-theta(4)*x(3);...%theta(3)*x(3);...%theta(4)*x(3)];output1=[theta(1)*(x(1)).*(N0-x(1)-x(2))-(theta(3)+theta(4))*x(1);...theta(2)*x(1)-(theta(3)+theta(2))*x(2);0;0];%output1=[-theta(1)*x(1)*x(3);...%theta(2)*(theta(1)*x(1)*x(3)-x(2)*(theta(3)-theta(4)*x(2)))/(theta(3)-2*theta(4)*x(2));...%x(2)*(theta(3)-theta(4)*x(2))-theta(5)*x(3);...%theta(5)*x(3)];elseifstrcmp(flag,’init’),output1=[0,59];output2=[69999007000000-23-6999900230];output3=odeset;end%/(1+theta(2)*7000000+sqrt(1+2*theta(2)*7000000))dydq:functionE=dydq(theta)formatlong;x=[];%输入数据长度B=[]；%对应地区死亡人数C=[]；%对应地区出院人数A=[]；%对应地区累计个案y=A-B-C;s0=theta(4);n0=theta(5);[t,s]=ode45(‘fd’,[0:93],[s0n0230],[],theta);%fun=inline(‘s(:,3)’,’theta’,[]);%quad8(fun,0,73)%fun=inline(‘‘,’theta’,’x’);%theta=lsqcurvefit(fun,[111],x,y)E=sum((y-s(:,3)).+(B+C-s(:,4)).)%+(C-s(:,5)).)%(y-s(:,3)).%(B-s(:,4)).+(C-s(:,5)).+(A-s(:,3)%plot(t,s(:,3),t,y,’o’)%figure%plot(t,s(:,1),t,s(:,2),t,s(:,3),t,s(:,4),t,s(:,5))%[X,FVAL,EDITFLAG]=fminsearch(‘HK’,[0.00C=[]；%对应地区出院人数A=[]；%对应地区累计个案y=A-B-C;%mm=max(s(:,3))[t,s]=ode45(‘fd’,[0:100],[theta(4)theta(5)230],options,theta);a=theta(1)*(theta(4)+theta(5)+23)a1=1/theta(2)a2=1/theta(3)a3=a*a2%fori=1:80%ss(i)=s(i,3)/mm;%end%ss;%T=0:79;%plot(T,ss)%figureplot(t,s(:,3))holdonT=0:93;plot(T,y,’o’)xlabel(‘时间’);ylabel(‘染病人数’);title(‘对应地区：

对应时间’);%输入对应数据figureplot(t,s(:,4))holdonT=0:93;plot(T,B+C,’o’)xlabel(‘时间’);ylabel(‘移除人数’);title(‘对应地区：

对应时间’);%输入对应数据figure%plot(t,s(:,5))%holdon%T=0:93;%plot(T,C,’o’)%xlabel(‘时间’);%ylabel(‘出院人数’);%title(‘对应地区：

对应时间’);%输入对应数据%figureplot(t,(s(:,3)+s(:,4)))%+s(:,5))holdonT=0:93;plot(T,(y+B+C),’o’)xlabel(‘时间’);ylabel(‘累计人数’);title(‘对应地区：

对应时间’);%输入对应数据figureplot(t,s(:,1))xlabel(‘时间’);ylabel(‘易感人数’);title(‘对应地区：

对应时间’);%输入对应数据figureplot(t,s(:,2))xlabel(‘时间’);ylabel(‘潜伏人数’);title(‘对应地区：

对应时间’);%输入对应数据%figure%plot(t/93,s(:,3)/(theta(4)+theta(5)+23))%holdon%plot(t/93,s(:,4)/(theta(4)+theta(5)+23))%holdon%plot(t/93,s(:,1)/(theta(4)+theta(5)+23))%holdon%plot(t/93,s(:,2)/(theta(4)+theta(5)+23))%legend(‘实际染病人数’);%text(65,s(65,3),’\leftarrow拟合人数为s(65,3)’);%text(65,s(65,3),’\leftarrow实际人数309’);附录二：

北京地区SARS疫情表日期已确诊病例累计现有疑似病例死亡累计治愈出院累计预测结果（累计病例）4月20日33940218332874月21日48261025434114月22日58866628465334月23日69378235556454月24日77486339647504月25日87795442738524月26日988109348769614月27日11141255567810894月28日11991275597811774月29日13471358668313224月30日14401408759014215月01日155314158210015265月02日163614689110916115月03日174114939611517225月04日1803153710011817765月05日1897151010312118765月06日1960152310713419345月07日2049151411014120225月08日2136148611215221115月09日2177142511416821545月10日2227139711617522055月11日2265141112018622515月12日2304137812920822905月13日2347133813424423335月14日2370130813925223595月15日2388131714025723785月16日2405126514127323965月17日2420125014530724115月18日2434125014733224235月19日2437124915034924305月20日2444122515439524385月21日2444122115644724395月22日2456120515852824515月23日2465117916058224635月24日2490113416366724805月25日2499110516770424925月26日2504106916874724995月27日2512100517282825075月28日251494117586625105月29日251780317692825155月30日2520760177100625175月31日2521747181108725196月01日2522739181112425216月02日2522734181115725236月03日2522724181118925236月04日2522718181126325256月05日2522716181132125256月06日2522713183140325286月07日2523668183144625286月08日2522550184154325286月09日2522451184165325306月10日2522351186174725306月11日2523257186182125306月12日2523155187187625306月13日252271187194425316月14日25224189199425316月15日25223189201525326月16日25213190205325326月17日25215190212025326月18日25214191215425326月19日25213191217125326月20日25213191218925326月21日25212191223125326月22日25212191225725326月23日2521219122772532附录三广东省非典疫情累计表附录三广东省非典疫情累计表====================日期确诊死亡出院（预测）7.1151258145114956.30151258145114956.29151258145114956.28151258145014956.27151258145014956.26151258145014956.25151258145014956.24151158145014956.23151158145014956.22151158145014946.21151158145014946.20151158145014946.19151158145014946.18151158145014946.17151158145014946.16151158145014946.15151158145014946.14151158145014946.13151158145014946.12151158145014946.11151158145014946.10151158144914946.9151158144914946.8151157144914946.7151157144714946.6151157144714946.5151157144714946.4151157144714946.3151157144514946.2151157144514946.1151157144214945.31151157144114945.30151157143914945.29151157143814945.28151157143614945.27151257143414945.26151357143314945.25151357143214945.24151357142914955.2315135714281495日期确诊死亡出院（预测）5.22151357142314965.21151356141914965.20151456141514945.19151456138814955.18151456136314955.17151456135514945.16151556135314935.15151356134214935.14151056132714925.13150956132214895.12150956131914875.11150656131414865.10150456130514835.9150256128814745.8148556127314655.7147855125614585.6146555125114445.5145354124514325.4144151123714265.3143251123314175.2142151120814085.1141251120613994.30140551120113914.29139951120113824.28139551119113744.27138251118513634.26137951117813554.2513745011731351内蒙古非典疫情累计表====================日期确诊死亡出院7.1282282546.30282282546.29282282546.28282282546.27282282546.26282282546.25282282546.24282282526.23282282526.22282282526.21282282526.20282282526.1928228252617282282526.162822825261428228245612282282456.11283292446.10283292316.9283292246.8283292166.7284292076.6284291966.5284291906.4284291806.3284291716.2284291666.1284281655.31284281565.30284281415.29284281315.28285281275.27285281105.2628628975.2528627805.2428626755.232862567日期确诊死亡出院5.2228525615.2128525585.2028725565.1928625525.1828925415.1729025405.1628824365.1528824235.1428823235.1328823235.1228820215.1128818205.1028917185.928417165.826616135.725816135.625114105.523014105.42251465.31901465.21801355.11541124.30127924.29120804.28114704.2776604.2674604.257060河北省非典疫情累计表====================日期确诊死亡出院（预测）7.1215122032646.30215122032626.29215122032626.28215122032616.27215122032616.26215122032606.25215122032606.24215122032596.23215122032586.22215122032586.21215122032576.20215122032566.19215122032556.18215122032546.17215122032546.16215122022546.15215122022536.14215122022536.13215122022536.12215122022516.11215122002496.10215121992496.9215121982486.8215121972486.7215121912476.6215121862476.5216121822466.4216121792446.3216121692436.2216121592416.1216121552395.31215121522395.30216121462375.29216121402365.28216111372365.27216111342355.26216111332315.25222111202295.24223111152265.2322811108224日期确诊死亡出院（预测）5.2222711972205.2122111882165.2021711792155.1921411632115.1821010572075.1720210432025.1620210352015.152009322005.141948271985.131918221945.121778211925.111668131875.101578111815.91568111765.8147691635.713465