江华 潘海霞 孙明伟 彭谨 杨浩 周舜泰 陈伟 周志远 刘义涛 何路 曾俊

DOI：10.3760/cma.j.issn.1671-0282.2014.09.006

基金项目：四川省科技厅科技支撑项目（SZ20120209、2012SZ0181）；成都市科技局科研基金资助（11PPYB099SF-289、12PPYB181SF-002）

作者单位：610101 成都，四川省医学科学院·四川省人民医院急诊医学与灾难医学研究所创伤代谢组多学科实验室 计算数学与生物统计学教研组（江华、潘海霞、孙明伟、彭谨、杨浩、周舜泰、陈伟、周志远、刘义涛、曾俊），急救中心急诊外科（江华、孙明伟、曾俊），肿瘤科（潘海霞）；北卡罗来纳大学教堂山分校医学院 Lineberger癌症研究中心（周舜泰）；中国医学科学院北京协和医院肠外与肠内营养科（陈伟）；泸州医学院生物化学教研室（周志远）；新西兰怀卡托大学管理学院（何路）

通信作者：曾俊，Email：cdjianghua@gmail.com

【摘要】目的 本研究旨在模拟外来输入性埃博拉出血热（Ebola virus disease， EVD）病例出现后，不同的防控策略对该病在我国传播及爆发流行模式的影响，并为相关决策部门提供一种循证决策依据。

方法 选用传染病SIR动力学模型建模。模拟病毒输入区域为我国南方某特大型城市群，该城市群人口5000万，所有人均为易感者。 出现输入性病例，模拟情况分别包含单个感染者入境I（0） = 1，某小型旅行者团体（全为感染者）入境I（0） = 5及I（0） = 10。I（0）我们称为“零号感染者”，其发病过程呈现典型的埃博拉出血热发病过程。其传播能力用β表示。对于输入性埃博拉感染者进入我国后，从其发病到得到控制（启动紧急预案）的时间t*设定了如下特征参数：t* = 6，12，18，24，30 h。传播模式设定如下：S0=5×108，I0[1， 20]；λ（0.1）；q（0.1]；t*{6，12，18，24}h。

结果 零号感染者输入后，随着应急预案启动的时刻t*变化，病毒暴发曲线为指数型，指数方程为y=1.053e1.801x ，模型擬合度R2=1，RSME=0.118 5。……