广东省SARS传播模型实证研究

1、广东省SARS传播模型实证研究广东省sars传播模型实证研究【摘要】目的讨论严重急性呼吸系统综合征severeauterespiratrysyndre，sars在人群中的传播规律，研究预防措施对疾病的影响，为预防控制sars提供科学根据。方法用确定性的sir模型和seidr模型等传播动力学数学模型对sars流行过程进展估计与分析。结果通过对卫生部2022年4月22日公布的广东省的数据，用确定性sir模型估计消除比。广州市sars的根本繁殖率为3.5左右。结论确定性的sir模型、确定性seidr模型能比较好模拟并估计出sars传播过程，并可以对将来的sars流行趋势和流行规模进展预测。【关键词】

2、流行病学;疾病传播sirdel，allparaetersithepideilgialeaninginludingtransis-sinrate，revalratiandthreshldvaluehadbeenestiated.itasnsideredthatthebasirepr-dutinratebrrsarstransissin.【keyrds】epideilgy;diseasetransissinhinjdisntrlprev2022，106:560-563传染性非典型肺炎是一种sars冠状病毒引起的急性呼吸道传染病，是21世纪人类遭遇的第一种能引起世界范围内流行的新发传染玻该疾病埋伏期

3、短，发病急，病死率高，抗生素根本治疗无效。sars在我国以及全世界多个国家流行，给社会带来很大的危害和冲击。广东省对sars病例流行病学调查，根本上认为sars传播方式为接触传播或短间隔 的飞沫传播。对于未知的并且没有很好治疗措施的一种急性传染病，最有效的控制措施是疫情公示和隔离，通过媒体将疫情向公众公布，同时将感染人群和疫区与外界隔离，防止疾病的扩散更显得尤其重要。但是隔离措施的施行也需要根据疾病传播的速度和强度而定，为此本文欲讨论sars的传播强度和应用动力学模型预测sars的流行趋势。1模型1.1确定性的sir模型假定研究对象为理想人群并且均匀混合，总人数保持在固定程度n，无迁入、迁出及

4、其他原因引起的死亡现象。并将总人群分为3类:易感者ssuseptible、患病者iinfe-tive;移除者rrevedlass。其中移除者包括疾病恢复的、感染后获得免疫的以及被隔离的。用st、it和rt表示每一类在t时刻的人数，同时假设:产生感染者速度是与易感者以及感染的人数成比例，即rsir0。易感者以同样的速度移除。感染者向移除类转移的速度与感染的人数成比例，即aia0。埋伏期足够短并可以忽略。dsdt=-rsi;didt=rsi-ai;drdt=ai以上就是经典的kerak-kendrik模型1。其初始状态为，s0=s00，i0=i00，r0=0，st+it+rt=ndidt=i0rs

5、0-，因此根本繁殖数r0=rs0，由sir模型微分方程组的前两个方程可以解出:is=i0+s0-s+lnss0。令=/r，表示阈值didt=i=n-r-s=n-r-s0exp-r/，r0=0假如疾病流行不是很大，即r/很小，取泰勒级数e-r/=1-r+12r2+的前3项近似得:drdt=n-s0+s0-1r-s0r222rt=2s0s0-1+tanht2-，s0-12+2s0n-s021/2，=tanh-1s0-1这样就可以近似把kerak和kendrik的移除率做如下近似估计2:drdt222s0seh2t2-;drdtaseh2t-;sehx=2ex+e-x;a=222s0;tanhx=e

6、x-e-xex+e-x;=2对于移除类人数，用以下计算公式来估计:rt=累计确诊人数+累计疑似人数-累计确诊的医务人员数，drdtrt+1-rt移除类人数的函数表达式，通过对该函数做最小二乘法拟合，估计模型中的参数。1.2确定性的seidr模型根据sars的疾病生态学和埋伏期研究，同时考虑到广东省sars较高的病死率，对传播动力学过程有较大的影响。故在一般的sir三分室模型的根底上增加埋伏期分室e、病死者分室d，设定五分室的模型如下3:考虑的模型的实用性和简单性以及广东省广州市sars埋伏期较短5d左右和死亡率较高的特点，对模型进展如下变形:s+e+d+i+r=p;dedt=hs-ie;did

7、t=ie-ri-ki;drdt=ri;dddt=ki;h=abi，p为总人口数在应用以上模型时，做了如下几个假设:疾病的传播是由人到人为唯一的传播方式。人群是均匀混合，人群的每一个个体接触率是一样的。由于sars没有有效的治疗药物，在研究期间恢复率和病死率不变。假设同一个群体中每一个个体每天与病人接触率和一次亲密接触感染成功概率相等。假设模拟时间很短，不考虑迁入和迁出率。模型中根本参数的含义如下:接种率h:人群中平均每人每天受到有效感染的概率，是由易感者s向埋伏期者e转移的速率。接触率a:人群中平均每人每天与一个病人发生近间隔 接触的概率。该指标反映隔离措施的估计。感染效率b:近间隔 接触成为

8、有效感染的概率，影响这个概率的因素有病毒的致病力。该指标也反映不同程度防护措施。转阳率i:埋伏期天数的倒数，是埋伏期e向发病者i转移的速率。恢复率r:病人每天得到恢复且获得免疫的概率，是有传染性期限天数的倒数。病人日死亡率k:是sars的病死率与恢复率的乘积。根本繁殖率r0=abpr，维持最低传播的根本繁殖率的临界值为1。假如1，即1个感染者产生1个二代病例，也就意味着传播即将终止。估计在整个流行过程中根本繁殖率是如何变化的，以及通过对的变化估计干预措施和防护措施的效果，并且预测疾病的开展趋势。2模型的实证分析本次研究的sars患者是广东省疾病预防控制中心搜集的并经过回忆性调查的从2002年1

9、1月26日2022年6月3日符合sars临床诊断标准确实诊的sars病例。2022年4月14日前执行临床诊断标准。2022年4月14日后执行卫生部?传染性非典型肺炎临床诊断标准?4。2.1确定性的sir模型由于卫生部发布的广东省累计报告病例、累计疑似病例以及医务人员数据是从2022年4月22日开场，所以每日治愈或消除人数只能从2022年4月22日开场计算。根据前述假设，笔者将把确诊的普通患者以及疑似患者作为排除人员对待。使用atlab的最优化工具箱中提供的lsqurvefit函数进展最小二乘法拟合图1，确定了dr/dt的参数。其中状态初始值设为40，0.1，1。得出:a=66.2158，=0.

10、0887，=0.8030。图3发病者的拟合曲线与原始曲线比照图2。可以通过移除类人数r计算出易感人群s，通过易感人群s推导出发病者的拟合曲线图3。2.2seidr模型根据模型的假设，考虑到人群中任何两个个体之间发生接触的概率与人群的大小即人口数亲密相关，因此设接触率a等于总人口数的倒数，即=1/在医院这个特定的群体，接触率可能会更高。感染效率b:一些流行病学专家认为sars的感染效率为0.5左右。埋伏期:sars的埋伏期，广东sars病例的埋伏期为5d，所以转阳率i为1/5。恢复率r:有关sars病人有传染性的期限，h专家估计为10d左右，本文取sars的传染性期限为7d，即r=1/75。日死

11、亡率k:根据广东省报道的sars患者的病死率在5%左右，所以病人的日死亡率为r5%，这同美国hell等6对sars研究的日死亡率估计值一致。对广州市从2022年1月2日5月31日的发病数的拟合结果如图4，拟合的广州市初始人口数为1000万。考虑到sars在医务人员中的聚集现象，笔者将广州市sars每日发病数据中减除医务人员的每日发病数据，并对数据进展拟合图5。3讨论从以上分析来看，sir模型描绘出sars在人群中传播的特征，从理论上说明了采取预防措施的重要性和必要性，为sars防制提供了理论根据7。然而sars具有一定的埋伏期，大约在5d左右，因此该模型对发病者的拟合曲线与原始曲线比较不是很理

12、想。seidr确定性模型微分方程组考虑到sars较高的病死率和大约5d左右的埋伏期的影响。在sir模型的根底上加以改进8，能更好的符合疾病的实际发病过程，同时seidr模型的计算和相关概念与疾病的控制策略和预测有亲密关系。从整个流行时间段内，根本繁殖数在整体上呈现波动下降趋势恰好也提示政府采取的控制措施以及加强人群的防护意识的有效性。广东省建立报告监测系统，发出诊断标准及控制措施指引时间为2022年2月5日左右即35d;发布医院收治患者指引在63d左右;发布学校、公共场所指引在76d左右;而在92109d即4月4日21日这个时间段内根本繁殖数反而上升，这种反常变化恰好发生在sars诊断标准发生

13、变化的时间前后，所以笔者认为这种改变可能由于不同诊断标准造成的。而从整个发病过程来看sars病毒毒力是逐渐下降的。在对排除医务人员的sars数据以及广州市7区数据分析中，根本繁殖数也都是呈现波动下降，并且其波动时间与全市发病拟合的波动时间比较接近。综上所述，本课题通过对广东省和广州市sars发病数据采用不同的数学模型分析，进而理解sars在人群中的传播过程，发现本文所用的几种模型都能较好的模拟疾病的发生过程，并且能估计sars所可能造成的危害，并且能为卫生行政部门和疾控部门对于疾病的传播以及预防措施的制定具有很重要的参考价值。【参考文献】2黄德生，关鹏，周宝森.sir模型对北京市sars疫情流

14、行规律的拟合研究j.疾病控制杂志，2022，85:399-400.3吴开琛，吴开录，陈文江，等.sars传播数学模型与流行趋势预测研究j.中国热带医学，2022，34:421-423.4彭国文，何剑峰，林锦炎，等.广东省传染性非典型肺炎流行病学特征初步调查j.中华流行病学杂志，2022，245:353-354.5lipsith，hent，perb，etal.transissindynaisandntrlfsevereauterespiratrysyndrej.siene，2022，3005267:1966-1970.6hellg，fenirep，astill-garsa，etal.sarsutbreakinntari，hngkngandsingapre:therlefdiagn-sisandislatinasantrl