留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

冠状病毒SARS-CoV-2、SARS-CoV和MERS-CoV的传染动力学分析

梁凯豪 张文峰 张小花 吴卓葵 刘芹 张超龙 李梓龙

梁凯豪, 张文峰, 张小花, 吴卓葵, 刘芹, 张超龙, 李梓龙. 冠状病毒SARS-CoV-2、SARS-CoV和MERS-CoV的传染动力学分析[J]. 电子科技大学学报, 2020, 49(3): 349-356. doi: 10.12178/1001-0548.2020067
引用本文: 梁凯豪, 张文峰, 张小花, 吴卓葵, 刘芹, 张超龙, 李梓龙. 冠状病毒SARS-CoV-2、SARS-CoV和MERS-CoV的传染动力学分析[J]. 电子科技大学学报, 2020, 49(3): 349-356. doi: 10.12178/1001-0548.2020067
LIANG Kai-hao, ZHANG Wen-feng, ZHANG Xiao-hua, WU Zhuo-kui, LIU Qin, ZHANG Chao-long, LI Zi-long. Analysis of Spread Dynamics of Coronavirus SARS-CoV-2, SARS-CoV and MERS-CoV[J]. Journal of University of Electronic Science and Technology of China, 2020, 49(3): 349-356. doi: 10.12178/1001-0548.2020067
Citation: LIANG Kai-hao, ZHANG Wen-feng, ZHANG Xiao-hua, WU Zhuo-kui, LIU Qin, ZHANG Chao-long, LI Zi-long. Analysis of Spread Dynamics of Coronavirus SARS-CoV-2, SARS-CoV and MERS-CoV[J]. Journal of University of Electronic Science and Technology of China, 2020, 49(3): 349-356. doi: 10.12178/1001-0548.2020067

冠状病毒SARS-CoV-2、SARS-CoV和MERS-CoV的传染动力学分析

doi: 10.12178/1001-0548.2020067
基金项目: 广东省自然科学基金(2018A0303130136);广东省科技计划(2015A030402008, 2015A070704059);广东省教育厅育苗工程(2012WYM_0081);广东省乡村振兴战略专项省级资金(粤财农[2018]125号)
详细信息
    作者简介:

    梁凯豪(1983-),男,博士,主要从事压缩感知、数学建模、农业金融等方面的研究

    通讯作者: 张文峰, E-mail:zhwf999@163.com
  • 中图分类号: R181.2

Analysis of Spread Dynamics of Coronavirus SARS-CoV-2, SARS-CoV and MERS-CoV

  • 摘要: 为了分析比较COVID-19、SARS和MERS这3种传染性疾病的传染动力学,该文通过传染病增长率和传染抑制常数,建立了传染病的传播增长模型,然后通过非线性拟合得到SARS-CoV-2、SARS-CoV和MERS-CoV这3种冠状病毒传播增长模型的参数。分析表明,SARS-CoV-2的增长率约为SARS-CoV和MERS-CoV病毒的两倍,SARS-CoV-2的倍增周期为2~3天。湖北的传染抑制常数比其他地区低两个数量级,与湖北地区的现状吻合。
  • 图  1  沙特阿拉伯的MERS病例

    图  2  2020年COVID-19的真实病例数和拟合病例数

    图  3  2003年SARS的真实病例数和拟合病例数

    图  4  MERS四个爆发周期的真实病例数和拟合病例数

    表  1  拟合结果

    肺炎地区相关系数(R)决定系数(DC)增长率r0传染抑制常数$s$$( \times {10^{ - 4}})$${N_{\max }}$拟合值/例真实病例数/例平均相对误差/%
    COVID-19HB0.994 80.989 60.2110.02778 920.01059 98931.56
    GD0.999 40.998 70.2982.2701 314.4161 328−1.02
    ZJ0.998 40.996 40.3432.9571 159.9261 172−1.03
    HN0.999 30.998 40.2932.3441 249.6591 257−0.58
    SARSNat0.998 40.996 70.1480.2765 355.9915 3270.54
    GD0.982 20.964 70.1070.7031 516.1511 5120.27
    BJ0.999 70.999 40.1760.6972 520.5222 521−0.02
    HK0.998 60.997 20.0850.4821 755.7881 7550.04
    MERSC#10.999 60.999 20.0911.434632.7626290.60
    C#20.999 70.999 40.0794.899161.8231573.07
    C#30.998 10.995 70.1485.482269.4752670.93
    C#40.999 70.999 40.1084.736228.4992270.66
    下载: 导出CSV
  • [1] DONNELLY C A, GHANI A, LEUNG G, et al. Epidemiological determinants of spread of causal agent of severe acute respiratory syndrome in Hong Kong[J]. Lancet, 2003, 361(9371): 1761-1766. doi:  10.1016/S0140-6736(03)13410-1
    [2] LIPSITCH M, TED C, BEN C, et al. Transmission dynamics and control of severe acute respiratory syndrome[J]. Science, 2003, 300: 1966-1970. doi:  10.1126/science.1086616
    [3] RILEY S, CHRISTOPHE F, DONNELLY C A, et al. Transmission dynamics of the etiological agent of SARS in Hong Kong: Impact of public health interventions[J]. Science, 2003, 300: 1961-1966. doi:  10.1126/science.1086478
    [4] 姚玉华, 孙丽华. SARS流行病传染动力学模型[J]. 数学的实践与认识, 2004, 34(3): 1-5. doi:  10.3969/j.issn.1000-0984.2004.03.001

    YAO Yu-hua, SUN Li. Dynamics model of infection of SARS epidemic[J]. Journal of Mathematics in Practice and Theory, 2004, 34(3): 1-5. doi:  10.3969/j.issn.1000-0984.2004.03.001
    [5] 李铮, 陈曦, 滕虎, 等. SARS流行病传染动力学研究[J]. 生物化学与生物物理进展, 2004, 31(2): 167-171. doi:  10.3321/j.issn:1000-3282.2004.02.013

    LI Zheng, CHEN Xi, TENG Hu, et al. Infectious kinetics of SARS epidemic[J]. Progress in Biochemistry and Biophysics, 2004, 31(2): 167-171. doi:  10.3321/j.issn:1000-3282.2004.02.013
    [6] 徐恭贤, 冯恩民, 王宗涛, 等. SARS流行病的SEIR动力学模型及其参数辨识[J]. 黑龙江大学自然科学学报, 2005, 22(4): 459-467. doi:  10.3969/j.issn.1001-7011.2005.04.009

    XU Gong-xian, FENG En-ming, WANG Zong-tao, et al. SEIR dynamic model of SARS epidemic and its parameter identification[J]. Journal of Natrual Science of Heilongjiang University, 2005, 22(4): 459-467. doi:  10.3969/j.issn.1001-7011.2005.04.009
    [7] 谭欣欣, 冯恩民, 徐恭贤, 等. SARS流行病动力学建模及其参数控制系统的研究[J]. 工程数学学报, 2003, 20(8): 39-44.

    TAN Xin-xin, FENG En-ming, XU Gong-xian, et al. SARS Epidemic modeling and the study on its parameter control system[J]. Chinese Journal of Engineering Mathmatics, 2003, 20(8): 39-44.
    [8] 曹亚丹, 魏顺行, 徐恭贤. 埃博拉流行病的离散logistic动力学模型及其参数辨识[J]. 渤海大学学报(自然科学版), 2019, 40(3): 257-262.

    CAO Ya-zhou, WEI Shun-xing, XU Gong-xian. Discrete logistic dynamic model and its parameter identification for the Ebola epidemic[J]. Journal of Bohai University (Natural Science Edition), 2019, 40(3): 257-262.
    [9] 李睿琪, 唐明, 许伯铭. 多关系网络上的流行病传播动力学研究[J]. 物理学报, 2013, 62(16): 504-510.

    LI Rui-qi, TANG Ming, XU Bo-min. Epidemic spreading on multi-relational networks[J]. Acta Physica Sinica, 2013, 62(16): 504-510.
    [10] 李睿琪, 王伟, 舒盼盼, 等. 复杂网络上流行病传播动力学的爆发阈值解析综述[J]. 复杂系统与复杂性科学, 2016, 13(1): 1-39.

    LI Rui-qi, WANG Wei, SHU Pan-pan, et al. Review of threshold theoretical analysis about epidemic spreading dynamics on complex networks[J]. Complex Systems and Complexity Science, 2016, 13(1): 1-39.
    [11] KIM K, JUNG K. Dynamics of inter organizational Public health emergency management networks: Following the 2015 MERS response in South Korea[J]. Asia-Pacific Journal of Public Health, 2018, 30(3): 893-931.
    [12] NISHIURA H, ENDO A, SAITOH M, et al. Identifying determinants of heterogeneous transmission dynamics of the Middle East respiratory syndrome (MERS) outbreak in the Republic of Korea, 2015: A retrospective epidemiological analysis[J]. BMJ Open, 2016, 6(2): 1012-1033.
    [13] CHEN X, TAN M, ZHAO J, et al. Identifying influential nodes in complex networks based on a spreading influence related centrality[J]. Physica A: Statistical Mechanics and its Applications, 2019, 22(5): 536-557.
    [14] LU Y, LIU J. The impact of information dissemination strategies to epidemic spreading on complex networks[J]. Physica A: Statistical Mechanics and its Applications, 2019, 21(3): 379-396.
    [15] NING Y, LIU X, CHENG H, et al. Effects of social network structures and behavioral responses on the spread of infectious diseases[J]. Physica A: Statistical Mechanics and its Applications, 2020, 1(2): 38-57.
    [16] DING Q, LI W, HU X, et al. The SIS diffusion process in complex networks with independent spreaders[J]. Physica A: Statistical Mechanics and its Applications, 2019, 16(5): 1102-1136.
    [17] HUANG S, CHEN F, ZHANG Y. Global analysis of epidemic spreading with a general feedback mechanism on complex networks[J]. Advances in Difference Equations, 2019(1): 1-20. doi:  10.1186/s13662-018-1939-6
    [18] ZHANG X, WU J, ZHAO P, et al. Epidemic spreading on a complex network with partial immunization[J]. Soft Computing, 2018, 22(14): 4525-4533. doi:  10.1007/s00500-017-2903-1
    [19] LIU Q, AJELLI M, ALETA A, et al. Measurability of the epidemic reproduction number in data-driven contact networks[J]. PNAS, 2018, 115(50): 12680-12685. doi:  10.1073/pnas.1811115115
    [20] LU L, CHEN D, REN X, et al. Vital nodes identification in complex networks[J]. Physics Reports, 2016, 650(1): 1-63.
    [21] WANG W, LIU Q, LIANG J, et al. Coevolution spreading in complex networks[J]. Physics Reports, 2019, 820(1): 1-51.
    [22] 范如国, 王奕博, 罗明, 等. 基于SEIR的新型肺炎传播模型及拐点预测分析[J]. 电子科技大学学报, 2020, DOI: 10.12178/1001-0548.2020029.

    FAN Ru-guo, WANG Yi-bo, LUO Ming, et al. SEIR-based novel pneumonia transmission model and inflection point prediction analysis[J]. Journal of University of Electronic Science and Technology of China, 2020, DOI: 10.12178/1001-0548.2020029.
    [23] 杨政, 原子霞, 贾祖瑶. 基于迁徙数据估计武汉感染新型冠状病毒的人员数量[J]. 电子科技大学学报, 2020, DOI: 10.12178/1001-0548.2020030.

    YANG Zheng, YUAN Zi-xia, JIA Zu-yao. Estimating the number of people infected with COVID-19 in Wuhan based on migration data[J]. Journal of University of Electronic Science and Technology of China, 2020, DOI: 10.12178/1001-0548.2020030.
  • [1] 刘益安, 马瑞辰, 李国, 于奇, 刘洋, 胡绍刚.  负阻态忆阻Hopfield神经网络动力学 . 电子科技大学学报, 2023, 52(1): 38-43. doi: 10.12178/1001-0548.2022294
    [2] 吴宗柠, 狄增如, 樊瑛.  多层网络的结构与功能研究进展 . 电子科技大学学报, 2021, 50(1): 106-120. doi: 10.12178/1001-0548.2020068
    [3] 杨晨, 文成, 许小可.  新冠肺炎疫情对民族地区人口流动的影响 . 电子科技大学学报, 2021, 50(6): 938-946. doi: 10.12178/1001-0548.2021147
    [4] 林自展, 肖井华, 周金连, 吴晔.  基于观点动力学的在线点评研究 . 电子科技大学学报, 2020, 49(1): 155-160. doi: 10.12178/1001-0548.2018320
    [5] 杨政, 原子霞, 贾祖瑶.  基于迁徙数据估计武汉感染新型冠状病毒的人员数量 . 电子科技大学学报, 2020, 49(3): 330-338. doi: 10.12178/1001-0548.2020030
    [6] 陈端兵, 白薇, 王岩, 王敏, 俞伍平, 周涛.  新型冠状病毒肺炎防控效果的定量评估 . 电子科技大学学报, 2020, 49(3): 339-344. doi: 10.12178/1001-0548.2020144
    [7] 郑文, 赵偲, 李泽堃, 武啸泽, 胡涛.  基于Web数据挖掘的COVID-19流行病学特征分析 . 电子科技大学学报, 2020, 49(3): 408-414. doi: 10.12178/1001-0548.2020079
    [8] 谭索怡, 曹自强, 秦烁, 陈洒然, 赛斌, 郭淑慧, 刘楚楚, 蔡梦思, 周涛, 张伟, 吕欣.  基于密切接触者人数推断新冠肺炎疫情发展趋势 . 电子科技大学学报, 2020, 49(5): 788-794. doi: 10.12178/1001-0548.2020263
    [9] 许小可, 文成, 张光耀, 孙皓宸, 刘波, 王贤文.  新冠肺炎爆发前期武汉外流人口的地理去向分布及影响 . 电子科技大学学报, 2020, 49(3): 324-329. doi: 10.12178/1001-0548.2020033
    [10] 林俊锋.  基于引入隐形传播者的SEIR模型的COVID-19疫情分析和预测 . 电子科技大学学报, 2020, 49(3): 375-382. doi: 10.12178/1001-0548.2020083
    [11] 傅家旗, 刘敏, 邓春燕, 黄娟, 江明珠, 郭强, 刘建国.  复杂人流网络下的COVID-19传播模型 . 电子科技大学学报, 2020, 49(3): 383-391. doi: 10.12178/1001-0548.2020109
    [12] 喻孜, 张贵清, 刘庆珍, 吕忠全.  基于时变参数-SIR模型的COVID-19疫情评估和预测 . 电子科技大学学报, 2020, 49(3): 357-361. doi: 10.12178/1001-0548.2020027
    [13] 乔晓华, 徐毅, 孙玉霞, 武花干.  忆阻超混沌Lü系统的隐藏动力学特性研究 . 电子科技大学学报, 2018, 47(3): 402-409. doi: 10.3969/j.issn.1001-0548.2018.03.013
    [14] 楼凤丹, 周银座, 庄晓丹, 张新荣.  时效网络结构及动力学研究进展综述 . 电子科技大学学报, 2017, 46(1): 109-125. doi: 10.3969/j.issn.1001-0548.2017.01.017
    [15] 王伟, 舒盼盼, 唐明, 高辉.  网络传播动力学模拟方法评述 . 电子科技大学学报, 2016, 45(2): 288-294.
    [16] 史定华.  网络科学中的互联网加——理念、评述和展望 . 电子科技大学学报, 2016, 45(4): 616-624. doi: 10.3969/j.issn.1001-0548.2016.04.014
    [17] 符丁, 李明江, 黎路.  基于价值驱动的人类行为动力学实证研究和建模 . 电子科技大学学报, 2015, 44(5): 652-656. doi: 10.3969/j.issn.1001-0548.2015.05.002
    [18] 荣智海, 吴枝喜, 王文旭.  共演博弈下网络合作动力学研究进展 . 电子科技大学学报, 2013, 42(1): 10-22. doi: 10.3969/j.issn.1001-0548.2013.01.005
    [19] 王昱青, 陈华富, 尧德中.  一种扩展的功能磁共振BOLD动力学模型研究 . 电子科技大学学报, 2007, 36(2): 291-293.
    [20] 徐红兵, 吕炳朝, 陈光.  一类非线性动力学系统的变结构混沌控制 . 电子科技大学学报, 1999, 28(3): 283-285.
  • 加载中
图(4) / 表(1)
计量
  • 文章访问数:  15284
  • HTML全文浏览量:  4801
  • PDF下载量:  278
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-02-20
  • 修回日期:  2020-03-24
  • 网络出版日期:  2020-05-14
  • 刊出日期:  2020-05-01

冠状病毒SARS-CoV-2、SARS-CoV和MERS-CoV的传染动力学分析

doi: 10.12178/1001-0548.2020067
    基金项目:  广东省自然科学基金(2018A0303130136);广东省科技计划(2015A030402008, 2015A070704059);广东省教育厅育苗工程(2012WYM_0081);广东省乡村振兴战略专项省级资金(粤财农[2018]125号)
    作者简介:

    梁凯豪(1983-),男,博士,主要从事压缩感知、数学建模、农业金融等方面的研究

    通讯作者: 张文峰, E-mail:zhwf999@163.com
  • 中图分类号: R181.2

摘要: 为了分析比较COVID-19、SARS和MERS这3种传染性疾病的传染动力学,该文通过传染病增长率和传染抑制常数,建立了传染病的传播增长模型,然后通过非线性拟合得到SARS-CoV-2、SARS-CoV和MERS-CoV这3种冠状病毒传播增长模型的参数。分析表明,SARS-CoV-2的增长率约为SARS-CoV和MERS-CoV病毒的两倍,SARS-CoV-2的倍增周期为2~3天。湖北的传染抑制常数比其他地区低两个数量级,与湖北地区的现状吻合。

English Abstract

梁凯豪, 张文峰, 张小花, 吴卓葵, 刘芹, 张超龙, 李梓龙. 冠状病毒SARS-CoV-2、SARS-CoV和MERS-CoV的传染动力学分析[J]. 电子科技大学学报, 2020, 49(3): 349-356. doi: 10.12178/1001-0548.2020067
引用本文: 梁凯豪, 张文峰, 张小花, 吴卓葵, 刘芹, 张超龙, 李梓龙. 冠状病毒SARS-CoV-2、SARS-CoV和MERS-CoV的传染动力学分析[J]. 电子科技大学学报, 2020, 49(3): 349-356. doi: 10.12178/1001-0548.2020067
LIANG Kai-hao, ZHANG Wen-feng, ZHANG Xiao-hua, WU Zhuo-kui, LIU Qin, ZHANG Chao-long, LI Zi-long. Analysis of Spread Dynamics of Coronavirus SARS-CoV-2, SARS-CoV and MERS-CoV[J]. Journal of University of Electronic Science and Technology of China, 2020, 49(3): 349-356. doi: 10.12178/1001-0548.2020067
Citation: LIANG Kai-hao, ZHANG Wen-feng, ZHANG Xiao-hua, WU Zhuo-kui, LIU Qin, ZHANG Chao-long, LI Zi-long. Analysis of Spread Dynamics of Coronavirus SARS-CoV-2, SARS-CoV and MERS-CoV[J]. Journal of University of Electronic Science and Technology of China, 2020, 49(3): 349-356. doi: 10.12178/1001-0548.2020067
  • 新型冠状病毒SARS-CoV-2是一种与SARS-CoV和MERS-CoV病毒类似的冠状病毒。SARS疫情在2002~2003年对世界多个国家地区造成严重影响。MERS-CoV病毒于2012年在沙特出现,疫情主要发生在中东国家和韩国,其他国家地区偶有发生。

    SARS和MERS的动力学研究已有一些进展。目前对流行病传染动力学研究的数学模型主要有确定性模型和随机模型。文献[1]得到SARS的平均潜伏期约为6.4天,从出现临床症状到入院的平均时间在3~5天。在病死率方面,60岁以下的患者的病死率为13.2%,而60岁以上患者的病死率为43.3%。潜伏期和病死率会影响SARS的病例变化。文献[2]研究得出在未采取控制措施的情况下,一个传染性非典型肺炎病例会感染约2~3个继发病例,采取有力的公共卫生措施能降低疫情规模。文献[3]对香港SARS疫情进行分析,发现SARS疫情的特征包括单独的超级传播事件(SSEs)和持续的社区传播。在疫情初期,平均每个病例产生2.7次继发性感染,其中很大一部分来自医院传播。传播率在流行期间下降,主要是由于人口接触率下降和医院感染控制得到改善。文献[4]建立了SARS流行病的数学模型,计算出模型参数,分析了SARS流行特征,并对疫情发展进行了预测。文献[5]利用部分国家、地区及中国内地部分城市的数据,通过Logistic确定型增长模型进行拟合,揭示了各个地区SARS传染力不均匀的现象以及防控措施上的差异。文献[6]建立了带有潜伏期及终身免疫的SARS流行病的SEIR动力学模型,对模型中的参数进行了参数辨识,得到了一个SARS流行病的控制区域。文献[7]建立了带有潜伏期及终身免疫的SARS流行病SEIR模型及参数辩识系统,论证了该类控制模型的主要数学性质以及系统的流不变性和弱不变性。文献[8]研究了埃博拉流行病的传染动力学模型及其参数辨识问题,构建了埃博拉流行病的离散logistic动力学模型,给出了参数辨识优化问题。文献[9]研究了多关系网络上的流行病传播动力学问题,提出一种双关系网络模型(工作−朋友关系网),研究了多关系对流行病传播动力学行为的影响和复杂网络上流行病传播动力学的爆发阈值[10],梳理总结了SIS和SIR模型爆发阈值的异同。文献[11]研究了MERS在组织间公共卫生应急管理网络的动态发展,绘制了MERS爆发期间的沟通和响应网络模式。文献[12]研究了2015年韩国中东呼吸综合征(MERS)疫情异质性传播动态的决定因素,利用传播树识别超级传播者,估计不同类型宿主的再生产数,即单个主病例产生的平均继发病例数。文献[13]为了识别复杂网络中的超级传播者,提出了一种融合指数−扩展影响相关中心度,通过提取和综合传统中心度指标和扩展影响拓扑特征信息来识别节点的影响。经典的SIR(易感染−恢复)模型描述了感染的传播过程。人们对流行病的反应影响传播,因此许多研究也考虑了信息传播。文献[14]分析了信息传播的操作如何影响受感染的个体以及流行病的传播条件,并基于社区规模和个人意识可能对感染率有影响的假设,提出了一个SIR-A(易感−恢复−活动)模型,将感染和信息传播映射到双层网络。文献[15]使用复杂网络的演化博弈模型研究了流行病传播过程,通过假设个人在疾病传播过程中将自己的收益与相邻个体的收益进行比较来选择策略,使用4种不同的仿真原则来更新策略,分析了违反直觉的现象,发现不同的仿真原则会影响违反直觉的现象的范围和严重性。文献[16]通过独立传播者参与的流行病SIS(易感−感染−易感)调查,解释了在发生扩散动力学的复杂系统中,通常具有许多潜在关系可以促进传染病的传播过程。文献[17]研究人类行为和接触异质性对传染病传播的影响,在考虑与他人潜在接触次数不同的个体的不同恐惧程度基础上,提出了具有一般反馈机制的基于网络的SIRS流行病模型。文献[18]提出了一种基于复杂网络中的部分免疫和免疫无效性的新型病毒传播模型,即易感性−感染性−易感性−恢复性−易感性,并基于均值场理论研究了该模型在统一网络和无标度网络中的流行动力学行为。对于复杂网络中的传染动力过程,可以从社会人口数据中模拟感染传播过程[19],此外,关键节点识别和共同进化传播机制在流行病传播的控制中至关重要[20-21]。文献[22]基于复杂网络理论建立了带有潜伏期的2019-nCoV流行病SEIR动力学模型,对疫情拐点进行了预测。以上工作从确定性模型、随机性模型、复杂网络传播动力学等方面研究了SARS-CoV和MERS-CoV的传播规律,但是对于SARS-CoV-2的传播研究仍然不够。

    目前,COVID-19作为一种新型传染病,其传播机制尚不明确。虽然,SARS-CoV-2与SARS-CoV和MERS-CoV同为冠状病毒,但是其传染特性与后两者有何区别,仍然有很多问题需要研究。

    • 对传染病感染人数进行研究分析,假设t时刻的感染人数为N(t)。在新型冠状病毒发生初期,人们对病原体、传播途径等方面缺乏充分认识,导致对传染病的防控不够重视,在新型冠状病毒发生初期感染人数的增长率可以看做常数${r_0}$。但当感染病例达到一定数量时,感染人数的增长率随着感染人数N的增加而减少,增长率函数$r(N)$$N$的单调递减函数。简单起见,不妨假设增长率为感染人数的线性函数,$r(N) = {r_0} - sN$表示人们没有采取任何措施并且任由病毒随意传播的增长率,$s$为传染抑制常数,反映人们所采取的防控措施对传染病的抑制效果。$s$越大,表示人们所采取的防控措施效果越明显。假设在病毒传染力和防控措施的共同作用下,感染的最大人数为${N_{\max }}$,即可以得到当$N = {N_{\max }}$时,增长率$r({N_{\max }}) = 0$。此时,可以得到传染抑制常数$s = {r_0}/{N_{\max }}$

    • 根据模型假设,在连续模型下,单位时间$\Delta t$内增加的感染者人数$\Delta N$可以表示为增长率$r(N)$$N$的乘积,并且${t_0}$时刻的感染人数为${N_0}$,因此可以得到感染人数微分方程的初值问题:

      $$\left\{ \begin{aligned} & \frac{{{\rm{d}}N}}{{{\rm{d}}t}} = {r_0}\left( {1 - \frac{N}{{{N_{\max }}}}} \right)N \\ & N({t_0}) = {N_0} \\ \end{aligned} \right.$$ (1)

      式中,等式第一项${r_0}N$表示在无防控措施情况下传染病的自然流行趋势;第二项$ - {r_0}{N^2}/{N_{\max }}$表示防控措施对传染病流行产生的效果。在传染病流行初期,由于缺乏有力的传染病防控措施,疫情的发展规律主要受病毒本身的传染特性影响,第一项起到主导作用,这时可以认为$ - {r_0}{N^2}/{N_{\max }} = 0$,此时微分方程的解为:

      $$N = {N_0}{\operatorname{e} ^{{r_0}(t - {t_0})}}$$ (2)

      这说明在无人为干预的情况下任病毒随意传播,感染者人数呈多代指数级增长。在传染病爆发期,第二项的作用开始逐渐占据主导地位,新增发病人数逐渐减少并趋向于零,传染病得到控制。

      式(1)是一个可分离变量的微分方程初值问题,可以解得:

      $$N(t) = \frac{{{N_{\max }}}}{{1 + \left( {\dfrac{{{N_{\max }}}}{{{N_0}}} - 1} \right){{\rm{e}}^{ - {r_0}(t - {t_0})}}}}$$ (3)

      可得$N(t)$的二阶导数为:

      $$\frac{{{{\rm{d}}^2}N}}{{{\rm{d}}{t^2}}} = r_0^2\left( {1 - \frac{N}{{{N_{\max }}}}} \right)\left( {1 - \frac{{2N}}{{{N_{\max }}}}} \right)N$$ (4)

      根据式(3)和式(4)可以得到感染人数的两个规律:

      1) $\mathop {\lim }\limits_{t \to + \infty } N(t) = {N_{\max }}$,这说明无论在${t_0}$时刻感染人数的初始值${N_0}$如何,感染人数的极限总是${N_{\max }}$

      2)由$N(t)$的二阶导数可以知道,感染者变化率${\rm{d}}N/{\rm{d}}t$$N = {N_{\max }}/2$时取得最大值,即感染者人数达到极限值一半时是加速增长时期。这一点之后,新增感染人数会逐渐变小,并最终变为零。因此,在$N = {N_{\max }}/2$所对应的时刻可以称之为新增感染人数的拐点。

    • 传染病所感染的人数从$N$变化到$2N$的时间跨度称之为倍增周期。如果在无人为干预的情况下,根据式(2),可以得到倍增周期${T_0} = \ln 2/{r_0}$。倍增周期${T_0}$仅依赖于无人为干预的感染人数增长率${r_0}$

      在采取防护措施等人为因素干预下,倍增周期将不仅仅依赖于增长率${r_0}$。由式(3)可以得到$N(t)$的反函数:

      $$t = \frac{1}{{ - {r_0}}}\ln \left( {{{\frac{{{N_{\max }} - N}}{N}} \bigg/ {\frac{{{N_{\max }} - {N_0}}}{{{N_0}}}}}} \right) + {t_0}$$ (5)

      假设${t_1}$时刻的感染人数为$N$${t_2}$时刻的感染人数为$2N$,则可以得到倍增周期:

      $$T = {t_2} - {t_1} = \frac{1}{{ - {r_0}}}\ln \left[ {\frac{{{N_{\max }} - 2N}}{{2({N_{\max }} - N)}}} \right]$$ (6)

      式中,函数$T$是关于感染人数$N$的函数,且要求${N_{\max }} - 2N > 0$,其导数为:

      $$\frac{{{\rm{d}}T}}{{{\rm{d}}N}} = \frac{{{N_{\max }}}}{{{r_0}({N_{\max }} - N)({N_{\max }} - 2N)}}$$

      可知$T$为单调递增函数,当$N \to 0$时,$T \to {T_0}$,而当$N \to {N_{\max }}/2$时,$T \to \infty $。这说明,当出现拐点时,感染人数的倍增现象则不再出现。

    • 本文采用的数据均来自世界卫生组织(WHO)(http://www.who.int)和国家卫生健康委员会官方网站(http://www.nhc.gov.cn/)。国家卫生健康委员会自2020年1月20日起公布新型冠状病毒SARS-CoV-2的确诊病例和新增确诊病例数据,本文选择SARS-CoV-2的数据起止日期为2020年1月20日−2020年2月16日,选取的地区分别是湖北(HB)、广东(GD)、浙江(ZJ)、河南(HN)。SARS-CoV病毒确诊病例数据选择的起止日期为2003年4月21日−2003年6月30日,选取的国家地区分别是全国(Nat)(不包括香港、澳门和台湾地区)、广东(GD)、北京(BJ)、香港(HK)。图1给出了沙特阿拉伯从2012年第12周~2019年第24周的每周MERS-CoV感染人数的柱状图。可见,沙特阿拉伯经历了4次明显的MERS疫情爆发,分别是2014年第7周~第23周、2015年第1周~第12周、2015年第15周~第26周、2015年第28周~第39周,这4个爆发周期分别记为C#1、C#2、C#3、C#4。选取这4个周期的MERS发生病例进行分析。本文应用1stOpt软件进行方程拟合,拟合过程采用Levenberg-Marquardt优化算法。

      图  1  沙特阿拉伯的MERS病例

    • 利用式(3)的模型对COVID-19、SARS和MERS在各个地区发生的病例数进行拟合,拟合效果显著。表1给出了拟合的结果,其中相关系数是度量所研究变量之间线性相关的程度,其取值范围是0~1,越接近1表示变量之间的相关程度越高;决定系数定义为回归方程中因变量的已解释离差与总变差的比值,反映回归方程的因变量变化可靠程度,其取值范围也是0~1,越接近1表示方程拟合效果越好。从相关系数和决定系数看,他们的最小值分别是0.982 2和0.964 7。除了湖北的COVID-19病例数外,${N_{\max }}$与真实病例数都相当接近,平均相对误差的绝对值在0.03%以下。湖北的COVID-19的真实病例数59 989例是截至2020年2月18日的确诊累计病例,目前在湖北的确诊新增病例仍然在每日1 000~2 000例之间,预测疫情结束之后,湖北的确诊累计病例将达7万多例。

      表 1  拟合结果

      肺炎地区相关系数(R)决定系数(DC)增长率r0传染抑制常数$s$$( \times {10^{ - 4}})$${N_{\max }}$拟合值/例真实病例数/例平均相对误差/%
      COVID-19HB0.994 80.989 60.2110.02778 920.01059 98931.56
      GD0.999 40.998 70.2982.2701 314.4161 328−1.02
      ZJ0.998 40.996 40.3432.9571 159.9261 172−1.03
      HN0.999 30.998 40.2932.3441 249.6591 257−0.58
      SARSNat0.998 40.996 70.1480.2765 355.9915 3270.54
      GD0.982 20.964 70.1070.7031 516.1511 5120.27
      BJ0.999 70.999 40.1760.6972 520.5222 521−0.02
      HK0.998 60.997 20.0850.4821 755.7881 7550.04
      MERSC#10.999 60.999 20.0911.434632.7626290.60
      C#20.999 70.999 40.0794.899161.8231573.07
      C#30.998 10.995 70.1485.482269.4752670.93
      C#40.999 70.999 40.1084.736228.4992270.66

      对于COVID-19,选取了疫情较为严重的4个省份进行模型拟合,分别是湖北、广东、浙江、河南。湖北在2020年2月13日统计确诊病例时,统计方式发生改变,将SARS-CoV-2核酸检测确诊改为临床确诊,导致确诊病例出现跳跃性递增。2月13日当天湖北新增确诊病例14 840。因此在图2中,累计病例在2020年2月13日有一个跳跃过程。总体来说,4个省份的模型拟合效果都较好。

      SARS病例在中国各省份发生的情况与COVID-19不同。截至2020年2月18日,中国内地(不包括香港、澳门和台湾地区)累计发生COVID-19共72 436例,其中湖北59 989例,湖北的累计病例占据了中国内地病例的绝大部分。2003年,SARS在中国多个主要省份均有发生,疫情较为严重的几个主要省份所发生的病例数量之间并没有明显的差异。因此,在地区病例数上,选取了中国内地(不包括香港、澳门和台湾地区)的病例数以及广东、北京、香港这几个地区的病例数。拟合结果如图3所示,拟合效果较好。

      图  2  2020年COVID-19的真实病例数和拟合病例数

      图  3  2003年SARS的真实病例数和拟合病例数

      对于MERS,沙特阿拉伯经历了4次明显的MERS疫情爆发周期。对这4个周期的MERS病例的模型拟合如图4所示。这4个周期的平均相对误差最大的是第2个周期,为3.07%,总体拟合效果较好。

    • 增长率${r_0}$主要反映传染病自然传染的能力,它受多个因素的影响,包括病毒本身的传染性强弱、人口流动情况和民众的卫生素质。人口密度虽然也会影响传染病的传播,但对传播的影响并不十分明显。对于SARS-CoV-2的增长率${r_0}$,由表1可知,湖北、广东、浙江、河南4个地区的增长率介于0.2~0.35之间,平均值为0.286,相互之间的差异并不十分明显。就SARS-CoV来说,北京的增长率约为香港的两倍。这可能与北京的人口基数较大、人口流动较快有关。MERS-CoV病毒在沙特发生的4个周期的增长率相差不明显,平均值为0.106。将3种冠状病毒的传播能力进行比较发现,SARS-CoV-2的增长率明显比SARS-CoV和MERS-CoV的增长率大得多。SARS-CoV-2的增长率约为SARS-CoV增长率的两倍,说明SARS-CoV-2的传染力比SARS-CoV要大的多。这就解释了在中国COVID-19疫情比2003年的SARS疫情病例增长快得多的原因。

      图  4  MERS四个爆发周期的真实病例数和拟合病例数

      倍增周期与增长率的倒数有关,SARS-CoV和MERS-CoV的倍增周期相当,都是5~10天,而SARS-CoV-2的倍增周期只有2~3天,与文献[23]的分析吻合,在指数增长作用之下,SARS-CoV-2的患病人数会迅速增加。

    • 传染抑制常数$s$与一些人为干预因素有关,主要包括将染病者与正常人隔离,将确诊病例与疑似病例隔离,民众对传染病的重视程度,对疾病流行区的清洁消毒。这些人为干预措施与疫情发生后当地政府、医疗机构所采取的措施有密切关系,反映政府管理部门和医疗卫生部门对疫情的应急管理能力。原则上来说,疫苗也属于人为干预因素之一,但是这3种冠状病毒目前都没有疫苗。传染抑制常数$s$越大,则反映所采取的人为干预措施越及时有效,疫情越容易在短时间内被控制。反之,则表示疫情越容易失控。

      3种冠状病毒SARS-CoV-2、SARS-CoV和MERS-CoV的传染抑制常数表现出显著差异,见表1。对于SARS-CoV,全国内地、广东、北京、香港的差异不大,都在一个数量级${10^{ - 5}}$,广东和北京相当,约为全国总体情况的两倍。MERS-CoV的传染抑制常数比SARS-CoV高一个数量级,为${10^{ - 4}}$,第一至第四周期几乎是逐渐递增,这说明沙特在MERS疫情后续爆发周期采取了更强的应急干预措施,阻断MERS-CoV病毒传染。在广东、浙江、河南,SARS-CoV-2的传染抑制常数与MERS-CoV同一个数量级。但是,在湖北,SARS-CoV-2的传染抑制常数为$2.7 \times {10^{ - 6}}$,比广东、浙江和河南低两个数量级。这就解释了COVID-19在湖北失控的原因。当然,广东、浙江和河南的传染抑制常数比湖北高的多,也得益于湖北所采取的强有力的“封城”措施。2020年初,广东、浙江和河南也较早启动突发公共卫生事件一级响应,使得政府管理部门、卫生管理部门和民众都采取了较为积极的应对措施。中国绝大部分省市最后都启动突发公共卫生事件一级响应,要求民众在家隔离、外出戴口罩,提高个人卫生要求。而针对湖北省内来说,由于染病者众多,早期而言,湖北的医疗机构应付不暇,医院的病床、设备、药品、口罩、防护服等物资严重不足,导致染病者不能有效隔离,因此导致湖北的传染抑制常数显著偏高。

    • 本文建立了传染病的动力学数学模型,分别利用该模型分析了SARS-CoV-2、SARS-CoV和MERS-CoV的流行特点。传染病增长率决定了传染病初期的流行程度,而传染抑制常数取决于各地采取的防控措施。3种冠状病毒动力学模型的参数分析合理解析了3种冠状病毒传播的特点以及疫情爆发过程中各地所采取的措施。SARS-CoV-2的增长率比SARS-CoV和MERS-CoV的增长率大得多,倍增周期约为SARS-CoV和MERS-CoV的1/2。在湖北,SARS-CoV-2的传染抑制常数比其他地区低两个数量级,使得湖北地区的疫情比其他地区严重得多。

参考文献 (23)

目录

    /

    返回文章
    返回