1. 引言
寄主型传染病受到多种异质性因素的影响,如遗传、年龄、空间等等 [1] [2] 。对于像性传播疾病这样的直接性的传染病来讲,寄主群体的性活跃度(性伴侣的个数)是影响其传播的主要的异质因素 [1] [3] 。
从Lajmanovich和Yorke的一篇会议论文开始 [4] ,学者们建立了许多数学模型来研究性活跃度的异质性对于传染病传播的影响 [5] 。近来,用网络的拓扑结构描述传染病传播的数学方法很是流行,具体来说就是用网络中的节点表示个体,用边来表示个体之间的接触。有趣的是,许多刻画性传播疾病的社交网络都是服从幂律分布: [6] 。按照网络表示论,不同个体性活跃度的异质性可以用寄主群体的度分布来刻画 [7] 。
众所周知,包括HIV在内的一些性传播疾病属于易感–感染–移除(SIR)型传染病,以往大家关注的焦点是该型传染病的阈值现象 [8] 。但逐渐地有人认识到SIR型传染病可以映射到无穷维的键逾渗问题上来研究 [9] [10] 。相应的,传染病模型中的一些量可以用逾渗理论中的量来代替,如传播强度可以用键占据率来替换,用逾渗阈值来代替爆发阈值。
另外,Newman等人利用概率母函数的重整化方法 [11] ,得到了含有传播率的自组织方程。这一方法是求解Golton-Waterson分支树阈值问题的著名方法。正如Durrett [12] 所指出的,该方法用在随机图中并不是严格的,其中的问题是当传染病速率在阈值之上时,传染病的传播就会产生回路,即:被传染的个体会再次被传染。此时再用分支树模型会有误差,但他指出这种误差是很小的,所以经过一个小的变动,该自组织方程仍然适用 [12] 。
基于上述结果,本文研究当时,即:时的标度律问题。并发现在不同的底层结构下,如同质结构、泊松分布、几何分布等,SIR型传染病爆发流行的标度律系数都是1,而对于无标度网络,它的标度律系数是依赖于参数的。
本文的框架如下:第二节我们介绍SIR型传染病的键逾渗理论;第三部分给出不同情形下的爆发流行的标度律。
2. 具有任意度分布的随机图的键逾渗理论
在文献 [11] 中,Newman把概率母函数与键逾渗理论相结合,来研究具有任意度分布的随机图上的SIR型传染病模型。首先,从节点出发,节点的度分布(任选一个节点度为的概率)的概率母函数定义如下:
(2.1)
如果从边出发,即:从网络中任选一条边并沿着该边到达的节点,此时该节点除了原始入边之外,剩余边的总数叫做该节点的超出度,其相应度分布的母函数为:
(2.2)
当传染病在网络上传播时,它并不是通过所有的边进行传播而是以平均传播速率在某些边上进行传播,显然。在流行病中边被传染了又叫做被占领,正如Newman所说,网络中一旦有某些边被占领,那么原来网络的结构就会发生变化。此时要想得到传染病的爆发流行(流行病中感染人数的比例),就必须考虑节点的被占领的边的分布,并且他用表示该分布的母函数,另外还给出了与之间的关系:
(2.3)
类似的,定义为节点的超出度中被占领边的分布(任选一条边所到达的节点的剩余边中被占领的边数的分布)的母函数,它与之间的关系如下:
(2.4)
设表示的是从节点出发(随机选择一个顶点)得到的小规模爆发的传染病中个体的数目,并令为爆发规模分布的母函数,它可以表示为:
(2.5)
其中表示的是由单个感染个体以平均传播速率引发的爆发规模为的传染病的概率。令为任选一条边,并沿着它到达的节点与规模为的簇连接的概率的母函数,或者按照分支过程的观点,它表示的是第0代的一个感染个体,经过几轮感染之后所有的感染个体的总数目的母函数。关于的Dyson方程的自组织形式如下 [11] :
(2.6)
该方程可以简单粗略的解释为:随机选择一条边所得到的大集团的规模等于,沿着边到达节点所引发的通过被占领的边连接的所有簇(小集团)规模之后再加上节点本身 [10] 。类似的,随机选择一个节点得到的大集团规模母函数的重整化形式为:
(2.7)
在文献 [11] 中,Newman等人得到了传播速率的临界值,并给出它的表达式:
(2.8)
称为流行病的阈值 [10] 。
当时,一场流行病或者用逾渗理论里的“大集团”就会出现,此时用来表示大集团里感染的人数占总人数的比例,并且它等于:
(2.9)
此时的意义发生变化,它不再是所有的簇的规模分布的母函数,而是除了大集团以外所有小规模的簇的母函数。利用方程(2.7)我们又可以得到流行病规模的另一个表达式:
(2.10)
其中是自组织方程(2.6)的解,即:
(2.11)
上面所提到的表示的是随机选择一条边,并沿着它到达的某个节点仍然是健康的易感节点的概率,且它的范围是在0到1之间。
基于上述结论,我们可以研究关于的爆发流行的标度律,即:
(2.12)
其中表示不同结构下的传染病爆发流行的标度律系数。
3. 标度律
把式(2.4)和式(2.11)结合起来得到:
(3.1)
这时令,而,;它的导数、都是恒大于0的,也就是说,函数在区间上是单调递增且是凸的。并且
(3.2)
而式(2.11)可以变形为。当()时,该方程只有唯一的正解;而当 ()时,方程除了1以外还有一个之间的解。从图1也可以直观地看出函数的性质以及解的存在。
既然方程的存在性已经得到了证明了,那么很自然地下一步就想解出该方程并得到解的具体形式。接下来将考虑在阈值附近的解的具体形式。
度分布的异质性将会导致不同的疾病传播模式。为了研究传染病的爆发流行对于传染速率的敏感度,我们必须得到自组织方程(2.11)的解。事实上,当(在的附近)时,利用泰勒展开式,可以得到几种具体度分布下的方程的解的具体形式。下面以泊松分布为例给出该情形下的具体推导过程。
当节点的度分布为泊松分布时,,相应的母函数以及传播阈值。此时方程(2.11)就变形为:
(3.3)
显然是方程(3.3)的一个根,另外利用泰勒展式得到单位根以外的解:
(3.4)
Figure 1. Two curves of. When, there exists unique intersection in
图1. 两种情形下的非平凡解的存在唯一性
其中,并且当时,是一个无穷小量。此时方程(2.12)的左边
(3.5)
容易得到泊松分布的标度律系数近似为1。同质情形、两种节点和几何分布 [13] 标度律系数如表1所示()。
我们得到以上几种情形下的标度律系数都是1。
这里再考虑无标度网络的情形,即网络的底层结构服从幂律分布: (是归一化常数),文献 [14] 中对此情形作了分析计算。此时:
(3.6)
令和,再利用近似:
(3.7)
可以得到:
, (3.8)
当时,,并得到标度律系数
Table 1. The values of the scaling law coefficient corresponding to different degree distributions
表1. 不同度分布下的标度律值
, (3.9)
但当时,,得到当时标度律系数为。
总之,当网络的度分布服从幂律分布时,此时SIR型传染病爆发流行的标度律系数是依赖于参数的取值的,当时,标度律系数和上述几种结构情形一样是为1的,此时也符合无穷维键逾渗的标度律情形。但当和时,标度律系数是大于1的且等于。
4. 结论
数学模型有助于我们更好地理解底层接触模式是如何影响传染病传播的,除此之外,理论分析的结果也有利于我们设计统计方案来收集和分析数据资料。与其它研究SIR传染病模型相比,本文的方法有下面两个优点:一是把随机图理论和无穷维的逾渗理论结合起来,这就使得对SIR型传染病模型的研究变得更简单;另一方面,用标度律来刻画对于传播速率爆发流行的影响,这就从宏观角度更好的说明了R. May所提到的现象 [1] 。
不同于Newman,我们给出了当时,自组织方程(2.11)的解析解。利用那个不等于1的,我们得到了包括同质、泊松分布和几何分布情形下SIR型传染病爆发流行的标度律系数,它们的值都是1。这个结果证实了逾渗理论中的结论,即:无穷维的逾渗模型的标度律都是一样的。另外我们发现当底层的网络结构是幂律分布时,此时爆发流行的标度律系数是与参数相关的。
基金项目
国家自然科学基金(11001004)。
*通讯作者。
参考文献