摘要:目的:探讨并分析在新型冠状病毒肺炎疫情中基于SEIR模型干预措施的影响。方法:收集2019-nCoV疫情数据,分析2019-nCoV的传播特点,利用python建立SEIR模型,加入潜伏期传染率、感染人群变化率等新参数,尝试分析防控手段的有效性。结果:给出多种相关干预措施下疫情发展趋势的预测。模型显示,基于严格限制出行的隔离措施能够减2019-nCoV发病曲线发展的趋势,使潜伏和感染人群的峰值降低45.71%和29.90%,潜伏人群数量较快到达群峰;模拟发现,本次疫情的拐点将出现在3月初,但对比2019-nCoV疫情实际数据,疫情拐点或将提前出现。结论:在2019-nCoV肺炎疫情中基于SEIR模型干预措施是有效的,各地方政府应严格执行隔离制度,切断传播途径,巩固成果,全力遏制2019-nCoV传播。
新型冠状病毒(2019-nCoV)是一种β属的冠状病毒,引起以肺部炎症病变为主的疾病,导致一系列与SARS相似的呼吸系统疾病症状,如呼吸困难、呼吸窘迫综合征或脓毒症休克,增加患者进入重症监护室(ICU)的概率和病死率[1]。自2019年12月在湖北省武汉市发现2019-nCoV以来,疫情快速蔓延。国家卫生健康委员会疫情通报显示,截至2020年2月8日,确诊人数33738人,密切接触者371905人[2]。为快速抑制病情蔓延,国家出台多种举措,如延长春节假期、企业延迟复工、学校延期开学、限制出行、居家隔离等。本文尝试从2019-nCoV的特点入手建立SEIR模型,分析目前防控手段的有效性。
1、对象与方法
1.1 数据来源
本文数据均取自国家卫生健康委员会官方网站2020年1月21日至2月8日公布的统计数据。随着2019-nCoV的快速传播,自1月25日开始,各地纷纷启动“重大突发公共卫生事件Ⅰ级响应”,自2月2日起更有多地限制社区居民外出,每家每天指派一人购买生活用品,因此我们认定从1月25日开始,病毒即进入有干预传播状态,自2月2日起,病毒进入强干预传播状态。
1.2 SEIR模型的建立
本文尝试使用SEIR模型揭示2019-nCov的传播过程,鉴于2020年1月26日国家卫生健康委员会宣布2019-nCoV肺炎潜伏期具有传染性,因此本文设计了包含潜伏期传染性的SEIR流行病模型。模型涉及4类人群:易感人群(S)、潜伏人群(E)、感染人群(I)、移除人群(R),假设康复概率是γ,易感者的感染概率是β,感染者接触的人数是r,潜伏者会传染易感者的概率为β2,潜伏者每天接触的人数是r1,另假设潜伏者的传染率与易感者的感染概率相同,则有如下微分方程组:
其中,i是2019-nCov传播的天数,N是总人数,N=Si+Ii+Ri+E,iSi,Ei,Ii,Ri分别表示第i天时易感人群人数、潜伏人群人数、感染人群人数、移除人群人数,α表示潜伏者转化为感染者的概率。
1.3 确定方程相关参数
总人口N:人口数据来源为中华人民共和国国家统计局公布的2018年年末武汉人口数据,共计883.73万人[5]。
感染概率β、β2以及潜伏者转化为感染者的概率α的计算:在实际生活中,因采取了干预措施,假设感染者平均每天接触人数为r=10,潜伏者接触人数r1=10,疫情期,易感人群数目等于总人数减去感染和潜伏人群,感染人群的变化率为dIi/di=βIiSi/N-γIi≈(β-γ)I,依据国家卫计委在疫情爆发期公布的数据进行拟合后估算β=β2≈0.045。依据α=Ii/Ei计算得到α≈0.10。基本再生数R0的估算:R0是指在发病初期,当所有人均为易感者时,一个病人在其平均患病期内所传染的人数[6]。基于SEIR模型,由公式R0=1+λTg+P(1-P)×(λTg)²计算R0值,其中λ是早期感染率,λ=lnY(i)/i,Y(i)是截至到i时刻的感染人数,Y(i)取值来自国家统计局实际公布人数,潜伏期用TL表示,感染期表示为T,i生成时间Tg的计算可近似于TL和Ti的序列间隔,计算得出Tg∈(7.0,8.4),P为潜伏期占生成时间的比值,即P=TL/Tg,由公式计算可得2019-nCoV的P值为0.70,将P值代入R0公式计算可得R0≈2.38-2.72[7-8]。本文预测的感染人数来自Python建立SEIR模型进行参数拟合后的预测结果。
2、结果
据方程确定相关参数后,拟合SEIR模型计算R0结果如表1所示,并将相关参数代入Python建立SEIR模型进行参数拟合后的预测结果如图1所示。
表1 基于SEIR模型的基本再生数及关键参数
2.1 停工停学但不限制出行措施实行后2019-nCoV的爆发趋势
在采取停工停学但不限制出行防控措施下2019-nCoV的爆发趋势如图1所示,疫情会在2月底附近迎来顶峰,并且感染人数将远超目前的情况。
图1 不限出行措施实行后2019-nCoV的爆发趋势
2.2 在限制出行措施下2019-nCoV的爆发趋势
图2显示了限制出行措施,即2月2日起限制社区居民外出,每家每天指派一人购买生活用品政策,假设感染者平均每天接触人数r=3,潜伏者接触人数r1=5,治愈率较之前有所提高,从国家卫生健康委员会疫情通报数据上显示,治愈率提高25%,即:γ从0.095提高了25个百分点至0.11875,还可以观察到,在限制社区出行后,感染曲线明显平滑减缓,潜伏和感染人群的峰值分别降低了45.71%和29.90%,我们还发现潜伏者的峰值和感染者的峰值相差一周左右,潜伏者的峰值更快到达,成为本次疫情拐点的标志。此外,模型的峰值出现在第107天,从12月12日武汉市卫健委进行的回顾式调查发现第一例患者发病,判断其感染病毒时间应大致在12月初甚至更早的时间,所以,推断出模型疫情拐点出现在3月初。
图2 限制出行措施下疫情变化趋势
图3显示了实际和预测确诊病例的变化趋势,实际确诊病例已经出现了下降,显示疫情已经得到明显的控制,真实疫情拐点或将提前。
图3 实际和预测确诊病例变化趋势
图4显示了严格限制出行措施下真实疫情的变化趋势,数据来源于国家卫生健康委员会官方网站[3],我们发现,疑似病例已出现平缓趋势,并且从2月4号以来不再增长并在2月8日出现了下降,表明疑似病例峰值已经出现,密切接触者仍是以较高比例增长,证明仍需执行强防控手段限制出行,静待密切接触者拐点来临,解除医学观察人群曲线平缓,证明医学观察方法日益谨慎,医学观察人员目前增长趋势减缓,证明疫情有逐渐被控制。
图4 严格限制出行措施下真实疫情变化趋势
图5显示了真实疫情中重症、死亡、出院人数比例的变化趋势,数据来源于国家卫生健康委员会官方网站[3],我们发现随着治疗经验和手段的不断丰富,全国医院紧急驰援武汉的情况下,治愈人数比例有了明显的上升,重症和死亡人数比例在不断下降,近期重症比例升高的原因是,感染人数不再上升,重症人数虽然减少,但是比例略有上升,综合表明病毒疫情正在被逐渐控制。
图5 重症、死亡、出院人数比例的变化趋势
3、讨论
本文应用Python建立潜伏期具有传染性的SEIR模型,根据国家卫生健康委员会疫情通报数据计算模型参数,得出本次疫情相关预测数据,并和实际数据进行了对比分析,得出以下结论:
模型显示减少潜伏和患病者接触的易感人群数量,更够减缓疫情发展速度,并使潜伏和感染人群的峰值将分别降低45.71%和29.90%。因此,疫情一旦爆发应及时采取相应等级的应急响应举措,以减缓疫情发展及减少疫情峰值发病人数[9]。相对SARS疫情从2002年12月底发现到2003年5月9日温家宝总理签署国务院第376号令,公布施行《突发公共卫生事件应急条例》,历时4个月[10],本次从发现2019-nCoV疫情到全国各地政府纷纷启动进入采取的措施仅历时一个月,阻止疫情蔓延非常及时高效。预测模型实际确诊病例人数的拐点出现在3月初,但国家及时出台强力的管控举措,及早限制出行、居家隔离、提示及强制出行戴口罩等防护措施切断病毒传染途径,大幅度减少潜伏和感染人群与易感人群的接触人数[11],使实际确诊病例比预测病例出现大幅度下降,预示疫情拐点或将提前。
模型显示,随着治愈率的上升,疫情也将尽早得到控制。从图5真实疫情中重症、死亡、出院人数比例图发现,在实际疫情中随着医疗经验和手段的增加,重症和死亡人数减少,出院人数增多,治愈率上升了25%,提示疫情可控,SARS病毒的毒株分离用了近5个月时间[12],而2019-nCoV病毒从2020年1月2日中国疾控中心病毒病所拿到武汉肺炎疫情第一份标本到1月24日,全球首发了第一株2019-nCov病毒毒株信息,毒株分离成功到向全世界公布仅用了一个月时间[13],这一重要信息公布使2019-nCov疫苗的研制成功更进了一步,进一步验证疫情可控。
应增加大众的医疗知识学习,提高群众传染病疫情防范意识[11],在疫情期间做好自我防护,主动拒绝食用野生动物,从源头切断病毒传播来源。2月7日上午华南农业大学针对新型冠状病毒肺炎疫情研究发现从穿山甲中分离出的毒株与新冠毒株相似度达99%[14],本次研究结果表明穿山甲是新型冠状病毒的潜在中间宿主。对本次疫情的源头防控具有重大意义,为相关政策从野生动物管控、重大疫情防控的角度出发,建立健全的监管机制,加强疫病防控知识教育和培训,有效防范人与野生动物交叉感染疫源、传播疫病等,提供了强有力的科学依据。
本文通过收集2019-nCoV在各级防控状态下的国家卫生健康委员会疫情通报数据代入SEIR模型得到2019-nCoV的基本再生数,应用Python建立潜伏期具有传染性的SEIR模型,得出本次疫情相关预测数据,发现预测结果的潜伏人群比感染人群更快达到峰值,这与实际的数据相符。本文建模过程中,因缺少二次感染人群数据,无法加入模型估算二次感染率,未来需更多微观感染数据对模型参数进行调整。
参考文献:
[1]中华人民共和国传染病防治法[J].中华老年口腔医学杂志,2005(01):56-58.
[3]中华人民共和国国家卫生健康委员会.截至2月8日24时新型冠状病毒肺炎疫情最新情况[EB/OL].2020-02-09]
[4]知乎.在家宅着也能抵抗肺炎!玩一玩SEIR传染病模型[EB/OL].
[5]中华人民共和国国家统计局.主要城市年度数据“武汉”2018[EB/OL].
[6]白振国.周期传染病模型的基本再生数[J].工程数学学报,2013,(2):175-183.
[7]周涛,刘权辉,杨紫陌,等.新型冠状病毒感染肺炎基本再生数的初步预测[J/OL].中国循证医学杂志:1-6.
[8]徐恭贤,冯恩民,王宗涛,等.SARS流行病的SEIR动力学模型及参数辨识[J].黑龙江大学学报(自然科学版),2005(04):459-462+467.
[9]王行兵,胡燕,吴满琳.隔离措施对北京SARS疫情控制影响的仿真分析[J].算机仿真,2005,22(11):299-302.
[10]刘长永.小汤山"非典"医院七昼夜建成[J].中国经济周刊,2009,(38):91.
[11]欧阳芬,吴荷玉,杨英,等.基于新型冠状病毒肺炎快速传播的应对措施[J/OL].全科护理:1-2.
[12]王立贵,张霞,郝荣章,等.传染病疫情现场预测预警系统的构建[J].医疗卫生装备,2012,33(7):10-11.
[13]吴少敏,钟哲.华农:从穿山甲中分离出的毒株与新冠毒株相似度达99%[EB/OL].
耿辉,徐安定,王晓艳,张勇,尹小妹,马茂,吕军. 基于SEIR模型分析相关干预措施在新型冠状病毒肺炎疫情中的作用[J]. 暨南大学学报(自然科学与医学版),:1-7.
[基金项目]国家社会科学基金一般项目(16BGL183);陕西省软科学研究计划一般项目(2017KRM211);西安交通大学第一附属医院院基金(2019RKX-05)
分享:
新型冠状病毒感染疫情是近年影响最大的一次传染病疫情,世界卫生组织将引起此次疫情的病毒命名为“新型冠状病毒”(SARS-CoV-2)。国家卫生健康委发布的各版新型冠状病毒感染(COVID-19)诊疗指南,均将病原学检测作为COVID-19的确诊指标,反转录聚合酶链反应(reverse transcription PCR,RT-PCR)作为检测人员日常工作中最广泛使用的检测方法。
2024-03-01药物性肝损伤(drug-induced liver injury,DILI)是大多药物的不良反应之一,是临床医师需要处理的最具挑战性的肝脏疾病之一,同时也是临床前研究及药物上市后导致其从市场撤出的最常见原因之一[10,11]。药物使用过程中,因药物本身和/或其代谢产物,或由于特殊体质对药物的超敏感性或耐受性降低导致DILI发生[10],临床上可表现为急性或慢性肝病。
2024-02-27新型冠状病毒感染(corona virus disease 2019,COVID-19)是病毒识别细胞受体血管紧张素转化酶2(angiotensin-converting enzyme 2,ACE2),ACE2表达在心脏、肾脏、脑、胃肠等器官组织[1],因此COVID-19侵犯呼吸系统外,还影响循环、神经、消化系统等各脏器。针灸是传统中医重要的特色技术,多版诊疗方案都有针灸干预指导意见。
2024-02-23目前,新型冠状病毒感染症(COVID-19)的诊断是通过RT-PCR进行的,但RT-PCR目前存在三个问题:(1)RT-PCR试剂包的供应跟不上需求;(2)城市以外的社区医院缺乏RT-PCR基础设施以适应高样本通量;(3)RT-PCR依赖于收集的样本中是否存在可检测到的SARS-CoV-2。
2024-02-22新型冠状病毒肺炎是由新型冠状病毒感染引发的呼吸道传染病,当前检测到的变异株有阿尔法(Alpha,α)、贝塔(Beta,β)、伽玛(Gamma,γ)、德尔塔(Delta,δ)和奥密克戎(Omicron)等[1]。2021年12月9日,我国首次检出奥密克戎变异株,其已成为主要流行株,目前已出现于多个省市;大多数感染奥密克戎者症状轻微,甚至是无症状感染,这使得该变异株隐匿性和感染性更强、传播速度更快。
2024-02-22新型冠状病毒感染发生以来,中国秉持生命优先的公共价值理念,构建以“人民至上、生命至上”为价值导向的抗疫卫生防控体系[1];形成群策群力、携手攻坚抗疫的良好局面。家庭是社会的微缩影,家庭环境卫生清洁消毒对打造社会环境清洁健康有重要的推动作用。因此,植根大健康理念,引导群众养成良好的卫生习惯[2],改善生存环境和生活质量,减少健康危害因素[3],有助于用文明素质与健康素养构筑起防控相关传染病的“铜墙铁壁”[4]。
2024-01-07新型冠状病毒感染(COVID-19)是由新型冠状病毒(severe acute respiratory syndrome Corona virus 2,SARS-CoV-2)感染引起的一种急性呼吸道传染病。自2019年新冠疫情席卷全球后,全世界范围内陆续发现了新型冠状病毒阿尔法(Alpha)、贝塔(Beta)、伽马(Gamma)、德尔塔(Delta)、奥密克戎(Omicron)等变异株[1]。2022年2月广州市报告中国首例奥密克戎亚变体BA.2毒株感染者[2],由于奥密克戎毒株潜伏期短、传染性强、传播速
2023-12-13新型冠状病毒肺炎(Corona Virus Disease 2019,COVID-19),简称新冠肺炎,是一种由新型冠状病毒(severe acute respiratory syndrome coronavirus 2,SARS-CoV-2)感染导致的急性呼吸系统传染性疾病,可表现为轻症至重症感染。自 2019 年末新冠肺炎疫情在全球肆意蔓延[1,2]。疫情期间肿瘤患者接受抗肿瘤治疗,存在感染新冠病毒的风险[3,4]。感染新冠病毒,影响抗肿瘤治疗时机和疗效,甚至发生严重病毒性肺炎引起呼吸衰竭。
2023-12-132020年中国认知障碍流行病学调查数据显示,我国60岁及以上老年人中痴呆症的患病率为6.04%,总人数超过1500万[1],位居世界首位。由于认知功能下降、高龄、多病共存、无法严格遵循疫情防控措施、脑组织中的血管紧张素转化酶2(ACE-2)表达增加等因素[2],痴呆病人发生新型冠状病毒感染(COVID-19)的风险增加,且病情进展快、死亡风险更高。
2023-12-12新型冠状病毒感染,又称新型冠状病毒肺炎(COVID-19),疫情爆发至今已有3年多时间,是全球仍然面临的一场异常艰巨的遭遇战。随着疫情的全球大流行,新型冠状病毒感染在我国境内传播与扩散的风险也将长时间持续存在。面对新型冠状病毒感染,中医药在疫情诊治中体现的独特优势受到业界的广泛关注。
2023-11-20人气:17986
人气:16674
人气:14049
人气:13856
人气:12909
我要评论
期刊名称:当代医学
期刊人气:8258
主管单位:中华人民共和国国家卫生和计划生育委员会
主办单位:中国医师协会
出版地方:北京
专业分类:医学
国际刊号:1009-4393
国内刊号:11-4449/R
邮发代号:82-829
创刊时间:1994年
发行周期:旬刊
期刊开本:大16开
见刊时间:7-9个月
影响因子:0.213
影响因子:0.000
影响因子:0.069
影响因子:0.345
影响因子:0.323
400-069-1609
您的论文已提交,我们会尽快联系您,请耐心等待!
你的密码已发送到您的邮箱,请查看!