2005年12月重庆师范大学学报(自然科学版)Dec.2005星丝查星兰塑!型翌型堕曼皇!!盟!坚盟!竺生旦!堕坐f翌翌堡!竺竺兰!!生竺2∑!!:丝翌!:兰用复杂网络理论研究疾病的传播+孙胜秋(北京交通大学理学院,北京100044)摘要:概述了用复杂网络的理论来研究两个传染病模型,即SIR模型和SIS模型,吡及在这两个模型中如何计算疾病的传播闯值。由复杂网络的局部渐近性质——在图中节点的总个数趋于∞时,图具有局部树形结构,可知SIR模型可以用分支过程理论来研究,从而根据母函数的方法可得到一些很好的性质=而SIS模型可以通过平均场的理论得毗理解,且知不论在度相关Scale—Free网络中还是度不相关Sesle-Free网络中都不存在非零传播阈值。关键词:复杂网络;度分布;渗流;传播阉值中图分类号:m1112;R181.2+5文献标识码:A文章编号:1672—6693(2005)04—0001—05StudyofDiseaseSpreading-UsingtheTheoryofComplexNetworksSUNSheng—qiu(SchoolofScience,BeijingJiaoTongUniversitytBeijing100044,China)ILece/tttoAbstract:InthispaperwereviewbrieflythemodelandSISmodelandtheproblemofhowthelocalasymptoticpropertiesofdiscussioninliteratureaboutthetwoinfectiousdiseasesmodels:SIRtwocalculatethetransmissionthresheldinvaluleofthemodelsAccordingtacomplexnetworks,whenthepitchpointnutitberinthegraphgoesSIRmodeltoirffinlty。thegraphislo-tocallyshapedintotree-likestuetttrebestudiedusingthetheoryofbranchingprocess.Itleadsgetsomegoodpropertiesbyusingthemethodofgeneratingfunctions.SISmodelwhetherindegreecorrelatedScale—FreenetworkstransmissionthresbohtorbeundemmodbythemesJflfieldtheory.andnotindegreeuncon'elatedScale—Freenetworks。theredoesexistKeywords:complexnetworks;degreedistribution;percolation;transmissionthreshold疾病的传播可看作是一个随机复杂网络,在这里用点代表可能通过传染或感染影响疾病传播行为的个体,如果两个个体之间通过某种方式可能直接发生传染与被传染的关系,就在他们之间连一条边。l模型把人划分为3类:易感类(susceptible,此类人并未染病但会被传染类个体传染),传染类(infective,此类人已患有疾病并会传染给其他人),康复类(re-covered,此类人已从疾病中康复过来,并具有免疫力)。SIR模型用于研究这3类人的数量随着时间的变化规律。用S(t),,(t),R(t)表示这3类人在某时刻t的数量,s(£)=S(t)/N,i(£)=,(t)/Ⅳ,r(t)=R(f)/N,用口表示单位时间内任意易感类个体被一传染类个体传染的一致概率,y表示单位时间内传染类个体康复并具有免疫力的比率,则此模型可用下列微分方程来描述:两个传统的传染病模型先介绍两个传统的传染病模型,此模型被称为充分混合模型。在此模型中,假定在网络中任意两个节点之间都可能存在传播疾病的边,相当于在完全图上讨论问题”“o。1.1SIR模型19世纪20年代LowellReed和WadeHamptonFrost…最初构建了SIR模型(充分混合模型)。该收稿日期:2005·09—05资助项目:国家自然科学基金重点项目(No10531070)作者简介.剥、胜秋(1982一),女,山东临沂』、,硕{:研究生,研究方向为复杂网络。万 方数据2重庆师范大学学报(自然科学版)第22卷一dt2一胁,面2肛一Tt,面2T。1.2SIS模型并不是任何疾病其康复者都会对这种病产生免疫力,一些不能自我控制但能通过药物治疗的疾病,康复的人通常会被病人再次传染上。结核病和淋病是两个大量研究的例子。计算机病毒也属于这种类型(计算机可以通过杀毒软件清除其中的计算机病毒,却抵御不了病毒的再次袭击)。所谓的SIS模型就是携带这种病的人治愈后不是从传染类转到康复类,而是转回到易感类。此模型可用微分方程描述:老=巾+'yi,面di=膨一Ti在上述2个模型中,虽然给出了疾病传播的基本的动力系统,但在实际生活中是不现实的,因为疾病只能在有实际身体接触的人之间传播,而在上述模型中任意人都有可能被所有病人传染,亦即疾病传播网络中所有结点具有相同的度(结点的度是指它连接边的个数)。在现实世界的疾病传播中,有必要考虑度分布不均匀的情况”1(度分布是指具有某个度的结点出现的概率分布)。近年来文献中有不少关于疾病在任意度分布网络中传播的研究‘8“3。笔者已经严格证明了在度序列给定的随机图中,有一个渐近局部性质:在图中节点的总个数趋于*时,随机选择一个点y,距离它小于等于m的所有点形成的导出子图L(y)无圈。在巨大分支还没出现的情况下,分支的大小都是有限的,结果表明有限分支包含圈的概率是趋于0的,是一个树形结构,所以可用G—W分支过程的理论来讨论疾病在任意度分布网络中传播的性质。本文简要介绍如何把传统的疾病传播模型推广为任意度分布网络上的传染过程。把传统的SIR模型推广到任意度分布网络之上的传染过程,在数学上可以把它看作在同一个网络上的边渗流过程”o。在疾病传播网络中,边以一定概率丁传播疾病相当于边渗流理论中边以一定概率r被占据,从而可以根据渗流模型得到有关传染病的预测情况。比如疾病爆发规模的分布,疾病的传播闽值,以及相变发生后传染病规模的大小。万 方数据下面所介绍的这个SIR模型,是讨论随机出现一个病人,在给定度分布的网络上传播疾病的过程,而源于很多实际情况。它描述某疾病有一个传染源,再在人群中传播开来的过程,例如SARS疾病的传播。作为初步近似,假定个体对之间的传播概率。与传染类个体的患病周期n均是同分布的。2.1.1度不相关网络在本小节简要介绍度不相关网络中几个重要的量。以P(k’lk)记以从网络中k度点出发到达t’度点的概率。所谓度不相关网络,是指P(k7Ik)与k无关。假设网络中具有某个度k的节点出现的概率(度分布)是P。。对x∈[0,1],记G。(m)=三p。x‘,称Go(x)为{PI}的母函数(母函数是研究非负整数值随机变量的有力数学工具)。在度不相关网络中,从任意点出发,到达一个度为t的邻居点的概率为黾。记吼一t=妻ZJ≥Pi或吼=…乙jp。坐j:;:;睑,叽表示以任意点出发到达一个邻居点,厶强j而这一邻居点另外还有k个邻居(称k为出度)的概率(在这里邻居是指有边相连的结点)。记{吼}的母函数为G。(z)=∑qkz‘。以H,(z)记在网络中随机选择一条边,该边所到达的有限分支的节点个数分布所对应的母函数。以凰(x)表示在网络中随机选择一点,该点所在分支的节点个数所对应的母函数。可以得到日.(z)=xqo+xql日1(z)+xq2[H1(Ⅲ)]2+…=xGl(Hl(z))140(z)=xpo+xplⅣl(z)+xp2[Ⅳl(x)]2+…=xGo(Hl(¥))笔者已证明网络中渐近地无圈,因此上述公式渐近地成立。由上述二公式可得随机选择点所在有限分支的平均大小(s)=1t'o(1)=1+#蒜-一L,1LJ当G’.(1)=1时,(s)=*,相变发生,出现巨分支。巨分支所占比例大小S=1一有限分支所占比例大小,因此皿(1)=G。(哦(1)),S=1一峨(1)=1一Go(H.(1))设u是G。(“)=u的最小非负实数解,由上式知日,(1)=“,故S=1一c0(u)。而u=tt,(1)是随机边到达有限分支的概率。2.1.2任意度分布网络中疾病的传播具体描述这个传播过程:随机出现的一个病人,通过与其邻居2现实世界中的传染病模型2.1任意度分布的SIR模型第4期孙胜秋:用复杂网络理论研究疾病的传播3接触,把病传给这些邻居,然后每个邻居再各自进行接触传播,如此进行下去,就可知道在整个网络中疾病的传播情况。假设i是传染类个体』是易感类个体,i把病传播给,的概率是o,个体i的患病周期为f.,个体康复后不会再传染他人,记t,为i传播疾病给J的概率。则i与J之间不传播疾病的概率为当T>T时,巨分支产生,在这里称之为有瘟疫发生,即疾病爆发。瘟疫影响人所占比例大小S(r)=1一‰(1,,)=1一c0(“,r)这里“;H.(1,T)是u;G。(u,T)的最小非负实数解。u表示随机边到达一有限分支的概率。一个一度点不受瘟疫影响的概率”=1一T+儿,这是因为与这点相连的边没有被占据的概率为1一T,与这点相连的边被占据但又指向一有限分支的概率为m,故一度为k的点不受瘟疫影响的概率为。‘,所以随机选择一度为k的点且没受瘟疫影响l一巧2舰(1一r。釉叫“2e”“,因此t.=1一e”n。设o,t均是同分布的随机变量,设其分布分别为P(r),P(r),则n也是同分布的随机变量,由此可得任意两个人之间传播疾病的概率为的概率为茄,其母函数为哿’故不受瘟疫影响的点(巨分支以外的点)的平均度为口G’。(")屯mT=<%)=1一【P(r)P(7)e~drd,",(0≤r≤1)。从而在此网络中,每条边传播疾病的概率为r,对应到渗流理论中,每条边被占据的概率为r。从而若随机选择一个点的度为k,与它相连的m条边被占据的概率为c?r(1一T)“…,故随机出现的一个病人,传染疾病给其邻居的个数分布的母函数为2i五FS2梨z:皿二!型瓣(卜s≤H)1一………G(")。G0(*,T)=Z∑仇c?r(一)…x。=,却^(一+o*r)‘=(1+(T1TGo1*一1)T)因为在任意度分布的网络中,随机选择点的邻点的出度分布的母函数为G。(x)=∑吼*‘,所以随机出现的一个病人,其每个邻居传播疾病人数分布的母函数为皇丢苦三击,母函数为鱼导掣,所以受瘟疫g—in:糍:生坠S旦墼巩影响的点(巨分支以内的点)的平均度为一1一繇(")4一随机选择一度为k的点且受瘟疫影响的概率为…。上式表明受瘟疫影响的点的平均度比不受瘟疫影响的点的平均度要大或至少相等。这跟实际情况吻合:在瘟疫发生时受瘟疫影响的人平均接触的人数要比不受瘟疫影响的人平均接触的人数要多。2.1.3度相关网络中的讨论在度相关网络中两邻点之间有某种相关性,即P(k’lk)与k和^’都有关。以联合概率分布ew=p(k’Ik)q;表示随机选择的边的两端顶点的出度分别为k和k7的概率∑p(k’lk)=1,e“,=e“。G。(z,T)=z∑吼c?r(一)…x‘=∑^一同理日,(x,T)=xG。(H.(x,T),T),风(x,T)=xG。(H,(z,T),T)。Ho(z,T)表示由随机出现的一个病人所在有限分支中被传染上疾病人数的分布的母函数,所以疾病爆发规模的平均大小为1TxTG1—1)T)从k度点出发到达有限分支大小的母函数为吧:=赤=黜=志时发蜘用(s)表示随机变量s的均值,当1一TO;(1)=O时,变。即该疾病在此网络的传播阈值为∑和。1。一∑☆(%一1)P^(s):=日’。c·,r,=-+器=,+i;;;;;;;i万方数据 风(*)=*旦—弋‘一,随机点所在整个分支大寺钆。∑e。【风,(x)]。小的母函数为日(x)=p。x+*王Pd砜一,(x)]‘。当有巨分支出现时,巨分支所占比例大小S=1一P。一G。(‰,),‰=巩(1)=上i_。这里‰表示从k寺%7度点出发的随机边到达有限分支的概率。类似1.2节其它讨论,在此不再详细表述。2.2现实世界的SIS模型∑e从,u∥”4重庆师范大学学报(自然科学版)第22卷因为康复的病人不具有免疫力,有可能再次被传染上,疾病的传播过程对应到网络中就相当于存在圈,因此不能用分支过程的思想来刻画疾病的传播过程。文献中用另外一种统计物理的平均场方法来研究在任意度分布网络中SIS模型疾病的传播。Pastor.Satorras和Vespi目lanl_6’针对SIS传染病给出了一个详细的平均场解法。用i,5。分别表示度为k的点中传染类个体和易感类个体各占比例。对于一个k度的点,用kAo(A)来代替传统模型中的口.(见第二节)。它表示的是单位时间内间易感类个体被传染所占的比例,这里A表示经接触一个传染类个体而被传染上疾病的概率。p(A)表示随机边到达的点属于传染类个体的概率。一般而言,日(A)与顶点度数有关,但在后面的讨论p(A)中对所有结点都相同(这也是统计物理中平均场近似的思想)。作为初步近似,在下面的讨论中还假设传染类个体的康复概率,=1。2.2.1度不相关情形在度不相关网络中,因为随在scale—free网络中不存在非零传播闽值。因为在sF网络中(k2)一*,所以A。=0,表明无论传播概率有多小,只要大于0,总会用瘟疫发生,Inter—net网上疾病的传播就是这种情况。2.2.2度相关情形本小节在度相关的网络中运用度的平衡性来讨论是否存在传播阚值。所谓网络度的平衡性是指,从k度点出发到达&7度点的总边数等于从A7度点出发到达k度点的总边数,用公式表述为如俨(k’lk)=&_^,P(klk’)。在度相关网络中,从一k度点出发的边到达一%7度点的概率为P(k’Ik),而该&7度点为传染类个体的概率为i。故从k度点出发的边到达度为k’的传染类个体的概率为P(k’lk)i。所以从k度点出发的边到达的点为传染类个体的平均概率为口(A)=∑P(k’1k)i∽从而在康复概率T=1时,有微分方程机边到达一k度点的概率为磐,而此点为传染类Lyypj些dt=ka0Ut)(1一i。)一i。即个体的概率为i。,从而随机边到达的点是传染类个∑印。i。掣叫m)+Ak[1屯㈤】pⅢtk)i“f)成立。若i。(t)的值很小,则(2)式的一阶近似为(2)体的平均概率为口(^)2专忑■。从而在721时有以下微分方程成立。掣=从日(A)(1一“(t))一i(f)考虑系统在一段时间以后的稳定状态,此时掣,乳““叭伊一“。…7’山其中£Ⅲ=一6“.+AkP(k7lk)为Jacobian矩阵(6n.是Kroneeker业dt-0'扣拦1A‰kOM吣)=豢得~+‘一+(A)’““”“’一∑如k”(1)delta函数)。当Jacobian矩阵£“至少存在一个非负特征值时,“=0不是上述方程的稳定解。而如果L。的特征值都严格小于0,则易知i。=O是上述方程的稳定解。设以。是矩阵c。.=细(%’1%)的最大特征值,则£w的特征值为一1+AA。,所以当一1+AA。<0时,方程的解i;=0是稳定解。记A。=A。~,由上述讨论知A,是疾病的传播阈值。根据Frobenius定理,对任意的非负不可约矩阵(A。),令/L。为其最大特征值,有口(^)2南莩概——口(^)=击莩概—Ak—01+AkO(&)表示随机变量的均值。有下面求i。的非零稳定解。对(1)式两端求导,未(南军",盅)。≥t即学:静巩趴:%,抓删,i。(t)在有限时间内减小到0,表明系统处于稳定状态时,疾病不会传播开来;而当A≥^。时,“(t)有非零解,表明系统处于稳定状态时,网络中疾病以大于0的概率持续传播,因此有疾病爆发。称^,为传播阈值。Am≥叩高善Aa,妒(6妒(k)=k,A=c2得7),(3)其中妒为任意的m维向量。在上述讨论中特别令A:≥n÷n善军^’扣(1I^)P(k’11)根据网络度的平衡性,有万方数据 第4期孙胜秋:用复杂网络理论研究疾病的传播5如伊(k’lk)=k'pIrP(klk’),(4)对(4)式两边同乘以^,再对k,∥求和,得到(k2)=善☆_F;如(kIk7)令;。(k’)=∑如(kl☆)表示&’度点的最近邻的平均度,故(^2)=∑&’‘.k。(k’)。在度相关scale-free网络中,可以证明如果幂率指数y<3时,(k2>随着节点增加而趋于*,此时k。(k’)一*,故有A:≥nin王如(1Ik)k。。(2)—+*,因此由上述的讨论知^,=0,即此时疾病传播的阈值为0。3结束语通过复杂网络的理论来研究疾病的传播,能深入理解到网络的拓扑结构对疾病传播有重要影响,从而找到控制疾病传播的更有效的方法,这对于社会上的很多流行病提供了一些理论研究的基础。本文介绍的两个模型只是真实情况的简单抽象,能够说明一些问题,当然也存在一些与实际不符的地方,比如缺少与实际传染病学和计算机病毒学所对应的参数等。因此建立更加真实的模型,并在此基础上研究疾病的传播是目前一个很重要的课题。参考文献:[1]ANDERSONRM,MAYRM.InfectiousDiseasesofHu—l'lftan¥[M]Oxford:OxfordUniversityPress,1991.[2]BARABA’SALBERTR,JEONGH.MemttleldTheoryforScale-freeRandomNetworks[J]PhysicalA,1999,272.173.187万 方数据[3]CALLWAYDS,NEWMANMEJ,STROGATZSH,eta1.NetworkRobustnessandFragility:PercolationRan-domGraphs[J].PhysRevlett,2000,85:5468-5471.[4]COHENR,EREZK,BEN·AVRAHAMD,etalResili—eneeoftheIntemettoRandomBreakDowns[J].PhysRev1ett,2000.86:3682-3685.[5]NEWMANMEJ,STROGATZSH,WATTSDJ.RandomGraphsWithArbitaryDegreeDistributionsand11lehAppli·cations[j].PhysRevE,2001,64:026118.[6]PASTOR.SATORRASR,VESPIGNANIA.Epide叫cSpreadinginScale-freeNetworks[J],PhysRevlett,2001,86:3200-3213.[7]BOGUNA’M,PASTOR—SATORRASR,VESPIGNAN[A.AbsenceofEpidemicThresholdinScale—freeNetworksWithConnectivityCorrelations.[EB/OL].http://www.arxiv.org/abs/eond-maI/0208163,2002[8]PASTOR.SATORRASR,VESPIGNANIAImmunizationofComplexNetworks[J].PhysRevE,2002,65:036104.[9]COHENR,BEN—AVRAHAMD,HAVLINS.Efficientlm—munizationofThePopulationandComputers[EB/OL].ht-tp://www.arxiv.org/abs/eond—mat/0207387,2002.r10]NEWMANMEJSpreadofEpidemicDiseaseonNet.works[J]PhysRevlett,2002,89:208701[11]VAZQUEZA,MORENOYResiliencetoDamdaeofGraphsWithDegreeGa'relariom[J].PhysRevE,2003,67:015101『12]NEWMANMEJ.TheStraetureandFunctionofComplexNetworks[EB/OL].http://wwwarxlv.org/absdcond—mai/0303516.2003.(责任编辑游中胜)用复杂网络理论研究疾病的传播
作者:作者单位:刊名:英文刊名:年,卷(期):被引用次数:
孙胜秋, SUN Sheng-qiu
北京交通大学,理学院,北京,100044
重庆师范大学学报(自然科学版)
JOURNAL OF CHONGQING NORMAL UNIVERSITY(NATURAL SCIENCE EDITION)2005,22(4)6次
1.ANDERSON R M;MAY R M Infectious Diseases of Humans 1991
2.BARAB A'S;ALBERT R;JEONG H Meanfield Theory for Scale-free Random Networks 1999
3.CALLWAY D S;NEWMAN M E J;STROGATZ S H Network Robustness and Fragility:Percolation on RandomGraphs[外文期刊] 2000(25)
4.COHEN R;EREZ K;BEN-AVRAHAM D Resilience of the Internet to Random Break Downs 2000
5.NEWMAN M E J;STROGATZ S H;WATTS D J Random Graphs With Arbitary Degree Distributions and TheirApplications[外文期刊] 2001
6.PASTOR-SATORRAS R;VESPIGNANI A Epidemic Spreading in Scale-free Networks[外文期刊] 2001(14)7.BOGU NA' M;PASTOR-SATORRAS R;VESPIGNANI A Absence of Epidemic Threshold in Scale-free NetworksWith Connectivity Correlations 2002
8.PASTOR-SATORRAS R;VESPIGNANI A Immunization of Complex Networks 2002
9.COHEN R;BEN-AVRAHAM D;HAVLIN S Efficient Immunization of The Population and Computers 200210.NEWMAN M E J Spread of Epidemic Disease on Networks[外文期刊] 2002
11.VAZQUEZ A;MORENO Y Resilience to Damdae of Graphs With Degree Correlations[外文期刊] 200312.NEWMAN M E J The Structure and Function of Complex Networks 2003
1. 余雷.薛惠锋.YU Lei.XUE Hui-feng 基于复杂网络的传染病传播模型研究[期刊论文]-陕西科技大学学报(自然科学版)2007,25(3)
2. 杨钢 传染病研究的几种方法[期刊论文]-焦作大学学报2007,21(3)
3. 张海峰.傅新楚.ZHANG Hai-feng.FU Xin-chu 含有免疫作用的SIR传染病模型在复杂网络上的动力学行为[期刊论文]-上海大学学报(自然科学版)2007,13(2)
4. 尹礼寿.闫喜红.YIN Li-shou.YAN Xi-hong 复杂网络中具有媒介传播SIS模型的稳定性分析[期刊论文]-长春大学学报(自然科学版)2010,20(2)
5. 许鹏远 复杂网络上的传染病模型研究[学位论文]2007
6. 陈进良.CHEN Jin-liang 基于小世界网络的甲型流感传播的研究[期刊论文]-计算机时代2010(5)7. 林文祖.陶军 复杂网络上的知识传播模型的研究[期刊论文]-中国科技信息2008(24)8. 王剑凌 复杂网络上传染病防控的研究[期刊论文]-福建教育学院学报2011(2)
9. 邵英英.刘孟.傅新楚.Shao Yingying.Liu Meng.Fu Xinchu 无标度网络上分片线性传染力与免疫作用下的流行病动力学[期刊论文]-动力学与控制学报2009,7(3)
1.杨钢 传染病研究的几种方法[期刊论文]-焦作大学学报 2007(3)
2.苟清明 一个具有阶段结构的SIS流行病模型[期刊论文]-四川师范大学学报(自然科学版) 2007(5)3.余雷.薛惠锋 基于复杂网络的传染病传播模型研究[期刊论文]-陕西科技大学学报(自然科学版) 2007(3)4.范旭.马军海.修妍 复杂供应链网络中的不确定性分析[期刊论文]-复杂系统与复杂性科学 2006(3)
5.黄淑琼.魏晟.岳丽.向浩.刘丽.刘莉.甘秀敏.许川.聂绍发 复杂网络模型在传染病预防控制中的应用[期刊论文]-中华流行病学杂志 2009(4)
6.贾秀丽.蔡绍洪.张芙蓉 一种动态的无标度网络模型[期刊论文]-四川师范大学学报(自然科学版) 2009(6)
本文链接:http://d.g.wanfangdata.com.cn/Periodical_cqsfxyxb200504001.aspx