埃博拉传播动力学模型及其应用摘要:埃博拉出血热,是由埃博拉病毒引起的急性疾病,致死率为25%到90%,在西非国家已经引起了人群恐慌,直接影响着人类身心健康。论文根据埃博拉病毒在三种不同情形下的传播原理,分别建立了不同的DSEIR模型,并且最后给出了埃博拉出血热的流行趋势以及防控策略。关键字:尸体的传播率;基本再生数;有效的防控措施;药物和疫苗模型1.引言埃博拉出血热,是一种致死率很高的急性传染病,若未能在短时间内及时采取相关的治疗措施可能会直接导致死亡。埃博拉病毒通过野生动物,比如果蝠(可能为该病毒的自然宿主),黑猩猩,豪猪等传播至人,并且通过人际间的传播蔓延开来。目前埃博拉出血热在西非传播广泛,尤其在几内亚,利比里亚和塞拉利昂,已经成为公共健康问题。自从西非在2014年3月报告出现首批病例,截至到目前为止,各项指标均显示此次疫情是1976年第一次发现埃博拉病毒以来发生的规模最大且病理最复杂的埃博拉疫情,本次疫情的病例和死亡数字已经超过了所有其他疫情的总和,并且疫情还在各个国家之间继续传播。因此,寻求埃博拉出血热的预防和控制措施显得尤为重要,然而当前的口服补液或者静脉注射等早期的支持性治疗护理方法以及对于特定症状而采取的症状疗法仅可以改善当前的生存率,目前仍未有获得官方认证并且实验证明可中和病毒的治疗措施,但正在开发各种血液,免疫,和药物疗法,当前也没有获得认可的针对埃博拉出血热的疫苗,但在2014年年底经过初步临床实验得到的有效数据的基础之上,世界卫生组织将联合各个有关部门在病例数最高的几内亚地区启动埃博拉疫苗效力试验。数学模型是历年来人们研究传染病的主要工具,它可以清晰的描述疾病的传播规律,结合已经获取的有效数据,能够在理论上进一步探讨不同防治措施下的效果,同时建立传染病动力学模型可以对传染病进行理论性的定性分析和定量研究,它可以描述疾病的传播,流行规律,在此之后根据实验数据预测疾病流行的最终趋势,并给出预警,为已经得出的防控措施提供理论基础和数量依据。所以,建立以及分析埃博拉出血热的动力学模型,有利于揭示该疾病的流行机制。因此,为了优化根除埃博拉,或者至少能够减轻目前的压力,在本文中,根据埃博拉出血热的传播特点,构建一个可行的数学模型;然后根据是否使用药物或疫苗,对模型进行二次优化,使得到的模型更加符合实际生活中疾病的传播情况;最后,基于对模型的分析与研究,能够预测埃博拉出血热在疫区的流行趋势并且给出预防和控制的最优策略。2.主要结果2.1基本假设i)易感个体从染病到死亡的时间较短,因此在这段时间内不考虑出生和自然死亡。ii)不考虑垂直传播。iii)埃博拉病毒对紫外线和伽马射线很敏感,因此不考虑病毒排泄物对疾病传播的影响。iv)潜伏期的人群不具有传染性。v)已免疫的个体不会被二次传染。vi)相关参数如下所示:几内亚:1因病致死率(1)0.38972康复率(2)0.19413尸体存放时间(T)0.1429(周)4治疗时间(RT)2(周)5总人口(N)112000006死亡率(Q)0.61197平均死亡时间(P)1.57(周)8潜伏期(W)1.64(周)9接触患者的传染率(1)5.5x81010接触尸体的传染率(2)1.3x71011及症状治疗等早期支持性医护办法可改善生存率。1111111潜伏期到染病期的转换率()0.6098利比里亚:`1因病致死率(1)0.27262康复率(2)0.28583尸体存放时间(T)0.1429(周)4治疗时间(RT)2(周)5总人口(N)54859956死亡率(Q)0.42847平均死亡时间(P)1.57(周)8潜伏期(W)1.64(周)9接触患者的传染率(1)5.5x81010接触尸体的传染率(2)1.3x7101111潜伏期到染病期的转换率()0.6098塞拉利昂:1因病致死率(1)0.23132康复率(2)0.31853尸体存放时间(T)0.1429(周)4治疗时间(RT)2(周)5总人口(N)37867646死亡率(Q)0.36317平均死亡时间(P)1.57(周)8潜伏期(W)1.64(周)9接触患者的传染率(1)5.5x81010接触尸体的传染率(2)1.3x710111111潜伏期到染病期的转换率()0.60982.2模型建立2.2.1不使用新药和疫苗条件下,传染模型的建立及分析在埃博拉疫区,埃博拉出血热主要在人际间传播,由于其传染性很强,因此易感人群可以看作疫区全部人口。一般情况下,埃博拉出血热的传播方式为直接接触感染者或者死者尸体,且当易感个体感染埃博拉病毒后进入潜伏期直到出现临床症状后患者才具有传染性。康复后的个体不会被二次传染,即不会再次成为易感个体,此外死者尸体的传染率极高,因此有关部门会在短时间内对死者尸体采取干预措施,即对尸体进行强制性埋葬或焚烧。基于目前的治疗和干预措施,在不考虑使用新药和疫苗的情况下,根据埃博拉病毒的传播特点以及模型假设,在疫区埃博拉出血热在人际间的传播关系见图1。SI1EIRIDSD2图1.仅采用当前治疗措施下的埃博拉出血热传染流程图在该模型中,分别用)()()()()(tDtRtItEtS,,,,表示t时刻易感,潜伏,染病,康复,死亡人群数量。对于易感仓室S,其输出为感染SDSI21,易感个体可以通过与染病个体接触而致病,且传染系数为1;也可以通过与尸体接触而致病,且传染系数为2。因此,单位时间内易感仓室S的变化如下所示:SDSIdtdS21,易感个体被感染将会进入潜伏仓室E,为SDSI21,潜伏仓室E的输出为临床爆发E。单位时间内潜伏仓室E的变化为:ESDSIdtdE21,对于染病仓室I,临床爆发E为输入,死亡1I和康复2I为输出。所以单位时间内染病仓室I的变化如下所示:IIEdtdI21,感染个体康复后将进入康复仓室R。因此单位时间内康复仓室的变化为:IdtdR2,SEIRD因病死亡1I将进入死亡仓室,其输出为尸体的焚烧或掩埋TD,所以单位时间内的变化如下所示:TDIdtdD1,由以上描述可得疫区仅在当前治疗措施下埃博拉出血热动力学模型为:,,,,,12212121TDIdtdDIdtdRIIEdtdIESDSIdtdESDSIdtdS2.2.2使用新药后的模型建立与分析在2.2.1中建立的模型基础上,如果考虑使用新药,那么不妨假设在单位时间(周)内的给药量为k(I)IM,其中M表示人均所需药量,)(Ik是比例系数,药效为l,且k)(I,l[0,1]。则根据当前的治疗措施并且结合新药治疗,此时康复率为R,致死率为D。结合2.2.1和模型假设,在采取当前治疗措施并且结合使用新药的情况下,疫区埃博拉出血热在人际间的传播关系见图2。SI1EIRIDSD2图2.当前治疗措施结合新药治疗下的埃博拉病毒病传染流程是否使用新药对易感仓室和潜伏仓室的输入和输出情况无影响,因此对于易感仓室S和潜伏仓室E与2.2.1中建立的模型完全一致,即有,,21SDSIdtdSSEIRDESDSIdtdE21,临床爆发E为染病仓室I的输入,其输出包括使用新药康复RTIkl,未使用新药康复(1-kl)2I和因病死亡Ikl1)1(。此时在使用新药后与2.2.1中的模型相比,单位时间内康复个体增加(既包括未使用新药康复又包括使用新药康复),死亡个体减少。所以单位时间内染病仓室I的变化如下所示:,)1()1(12IklIklTklIEdtdIR对于康复仓室R,其输入为使用新药康复RTIkl以及未使用新药康复Ikl2)1(。因此单位时间内康复仓室的变化为:,)1(2IklTklIdtdRR因病死亡Ikl1)1(将进入死亡仓室,其输出为尸体的焚烧或掩埋TD,所以单位时间内的变化如下所示:,)1(1TDIkldtdD由以上描述可得疫区在当前治疗措施结合新药的前提下埃博拉出血热动力学模型为:,)1(,)1(,)1()1(,,12122121TDIkldtdDIklTklIdtdRIklIklTklIEdtdIESDSIdtdESDSIdtdSRR2.2.3使用疫苗和新药条件下,传染模型的建立及分析在2.2.2的基础上,如果考虑使用疫苗,那么不妨假设疫苗的有效率为,并且假设疫苗长期有效,即易感个体注射疫苗康复后不会再次被传染。结合2.2.2和模型假设,在采取当前治疗措施并结合新药以及疫苗的情况下,疫区埃博拉出血热在人际间的传播关系见图3。SSI1EIRRDSD2图3采取当前治疗措施并结合新药以及疫苗情况下的埃博拉病毒病传染流程图对于易感仓室S,其输出包括感染SDSI21以及易感个体的免疫S。此时在使用疫苗后与2.2.2中的模型相比,单位时间内易感个体的输出增加,既包括被传染,也包括使用疫苗后直接进入康复仓室。因此,单位时间内易感仓室S的变化如下所示:,21SSDSIdtdS是否使用疫苗对潜伏仓室和染病仓室的输入和输出情况无影响,因此对于潜伏仓室S和染病仓室I与2.2.2中建立的模型完全一致,即有,ESDSIdtdE21,,)1()1(12IklIklTklIEdtdIR对于康复仓室R,其输入包括使用新药康复,RTIkl未使用新药康复Ikl2)1(以及易感个体的免疫康复S,因此单位时间内康复仓室的变化为:,)1(2SIklITklIdtdRR因病死亡Ikl1)1(将进入死亡仓室,其输出为尸体的焚烧或掩埋TD,所以单位时间内的变化如下所示:SEIRD,)1(1TDIkldtdD由以上描述可得疫区在当前治疗措施结合新药以及疫苗的前提下埃博拉出血热动力学模型为:,)1(,)1(,)1()1(,,12122121TDIkldtdDSIklTklIdtdRIklIklTklIEdtdIESDSIdtdESSDSIdtdSRR3模型预测及防控措施3.1疾病预测与防控措施在此次疫情中,几内亚,利比里亚和塞拉利昂的疫情最为严重,因此下面仅就这三个西非国家的实际情况给出疾病预测以及防控措施。在2.2建立的三个传染模型中,接触患者的传染率1和接触尸体的传染率2是未知的,其余参数均可由已知数据求出。对于1与2,则可以利用MATLAB软件对不同时间段累计病例进行拟合,从而得到具体的1和2的数值。通过对这三个国家不同时间段累计病例的数据收集,数据拟合的结果分别如下图所示:图4几内亚不同时间段累计病例数据拟合图(data1表示实际累计病例,data2表示拟合之后的累计病例)(下同)图5利比里亚不同时间段累计病例数据拟合图图6塞拉利昂不同时间段累计病例数据拟合图21和拟合数值分别为:国家12几内亚5.4x8101.3x710利比里亚8.5x8107.5x710塞拉利昂7.8x7105.7x710利用分别拟合出的三个国家接触患者的传染率1和接触尸体的传染率2,再结合2.2.1,2.2.2和2.2.3中建立的传染模型,能够预测出2015年三个国家在仅采取当前治疗措施,当前治疗措施结合新药以及当前治疗措施结合新药和疫苗三种情况下累计病例发展的最终趋势。图7,图8,图9分别给出了三个国家2015年埃博拉出血热在三种情况下累计病例数的最终趋势,通过对累计病例的预测能够反映出三个国际埃博拉出血热传播的最终趋势。图7几内亚在三种情况下的累计病例数的最终趋势(蓝色曲线:在仅采取当前治疗措施的情况下对疾病的预测。中间红色曲线:在采取当前治疗措施结合新药的情况下对疾病的预测。底部红色曲线:在采取当前治疗措施结合新药以及疫苗的情况下对