具有脉冲接种SIQR模型的危险源扩散动力学分析.docx

-1-

具有脉冲接种SIQR模型的危险源扩散动力学分析1

戴东甫赵林度

东南大学系统工程研究所，南京（210096）

摘要：在应对生物恐怖袭击事件中，如果能够准确地了解危险源扩散的规律，将有助于更好地规避或降低危险源扩散造成的风险。本文采用传染病动力学的方法，借鉴带有脉冲接种的SIQR模型来分析危险源扩散机制，得到了控制危险源扩散的脉冲接种周期和接种率的关系，并提出了一种在一定脉冲周期下，实现消除危险源最小接种率的具体算法。通过实际数据分析得到，可以在较小的脉冲接种率的条件下，控制危险源的扩散，比应用一次性接种的方法控制危险源扩散更有效。

关键词：SIQR模型；生物危险源扩散；脉冲接种；扩散动力学；生物恐怖

1、引言

生物恐怖威胁是现实存在的，它具有容易实施、施放突然、心理威胁大、后果严重等特点，以及生物恐怖威胁的跨国性、不确定性和高难度防扩散性。生物反恐的关键挑战来自不知何时、何地和何种方式遭受生物恐怖袭击[1]，只能通过疫苗、抗生素、药物和某些手段控制和预防。然而，如果能了解危险源扩散的规律，将有助于更好地规避或降低危险源扩散造成的风险。为此，国内外学者开展了大量的研究工作，主要研究了具有时滞的传染病动力学模型、具有年龄结构的传染病模型、在多种群中传播的传染病模型、非自治的传染病动力学模型、具有迁移的传染病模型和具有脉冲效应的传染病模型[2]。

相对于传统的致使疾病消除的方法，脉冲疫苗接种方法能够在很低的疫苗价值下消除疾病[3]。目前的研究还是起步阶段，国内外研究不是很多，主要有柳合龙，徐后宝等建立了带有脉冲免疫和传染年龄的SIV传染病模型，证明了此类解的存在唯一性[4]；勒祯，马知恩考虑了连续预防接种且传染率是标准的SIRS传染病模型，在连续预防接种和脉冲接种下，分别给出模型的基本再生数[5]；勒祯在脉冲预防接种下，给出了SIRS传染病模型基本再生数，证明无病周期解的存在性和全局渐进稳定性[6]。MGRoberts和RRKao（1998）研究了具有脉冲出生的SI模型，给出了无病平衡点、地方病平衡点、平凡平衡点和种群灭绝平衡点，并讨论了这些点的稳定性[7]；Borisshulgin，Lewistone和Zviaagur（1998）通过构建和分析具有脉冲接种的SIR模型，考虑了非线性和混沌情况下的SIR模型，给出了无病平衡点，并且讨论了其稳定性，给出消除疾病的条件[8]；他们后来（2000）进一步考虑了在此模型下消除疾病的最大周期的情况，给出了最大周期的确定方法[9]；LuZ，ChiX和ChenL（2002）在此基础上进一步讨论了具有脉冲接种的水平和垂直传染的SIR模型[10]；AlbertoD’Onofrio

1国家自然科学基金资助项目“生物反恐体系中应急物流网络优化与仿真研究”

-2-

（2002）通过分析带有脉冲接种的SEIR模型，给出了在几种不同情况下的无病周期解，并且对其稳定性进行了分析[11]，随后又有学者于2006年和2007年在此基础上进一步研究了具有脉冲接种的SEIR模型，并且给出疾病是否传播的阀值[12,13]。我国应付传染病扩散，常使用的措施之一就是对已感染人群进行隔离，但是国内外对具有脉冲接种隔离的传染病模型研究还很少。

本文假设如果一个国家或地区遭到生物恐怖袭击，区域内的每一个人都有可能成为感染者，每个已感染者都可以看成一个危险源，这样可以应用传染病动力学的方法研究危险源的扩散机制。

2、模型建立

2.1传统的模型

传统的SIQR模型：

要用一次性疫苗接种消除传染病，则要0，即：

βIS?(γ+δ+m+α1)I0（2）

从而，得到S，假设m，得到：

（3）

则最小接种率：

in=1?（4）

2.2具有脉冲接种的SIQR模型

2.2.1模型假设

（1）出生率等于自然死亡率；

（2）一旦一个人遭到危险源的传染，这个人将会体现出来，即危险源没有潜伏期；

（3）对危险源进行隔离处理时，危险源就不再扩散；

-3-

（4）一旦从危险源中恢复，就不会再次被危险源传播；

（5）除了进行特殊处理，否则人人都有可能遭到危险源的袭击；

（6）按一定周期对未遭到危险源传播的人群进行处理，使得不再遭到危险源的威胁，接种