王晶囡 逯蘭芬 高旭 許云中 丁悅航

摘 要：為了研究疫苗免疫和群外個體遷入對SIR型傳染病傳播行為的影響，在復(fù)雜網(wǎng)絡(luò)上建立了一類含有免疫作用的SIR傳染病模型，討論了具有齊次性質(zhì)的SIR復(fù)雜網(wǎng)絡(luò)平衡點穩(wěn)定性和具有非齊次性質(zhì)的無標度網(wǎng)絡(luò)（度分布關(guān)聯(lián)和度分布不關(guān)聯(lián)）染病平衡點的存在性，從中得到了3種情況的染病臨界值，根據(jù)臨界值討論了比例免疫和目標免疫策略優(yōu)缺點，為控制病毒傳染提供了的理論參考。

關(guān)鍵詞：無標度網(wǎng)絡(luò)；穩(wěn)定性；傳染病模型；臨界值；免疫

DOI：10.15938/j.jhust.2018.02.026

中圖分類號： O175；O193

文獻標志碼： A

文章編號： 1007-2683（2018）02-0144-05

Abstract：In order to study the effect of the vaccine immunization and the individual immigration on the spread behaviors of the SIR epidemic， the SIR epidemic model for complex networks is constructed. The equilibrium stability of complex networks with homogeneous characteristics and the existence of epidemic equilibriums of scale-free network with non-homogeneous （when the node degrees distribution are related and when the node degrees distribution）are discussed. The epidemic critical threshold values under three different conditions are obtained. Based on the epidemic critical thresholds values， the advantages and disadvantages of uniform ratio immunization and targeted ratio immunization are compared，which presents the theoretical reference for controlling viral infection.

Keywords：scale-free networks； stability； epidemic model； threshold value； immunity

0 引 言

有些病毒性傳染病的傳播范圍廣，速度快又無法被有效治療，對人類生命造成嚴重威脅，有效控制傳染病的傳播現(xiàn)成為急需解決的問題之一。據(jù)世界衛(wèi)生組織（WHO）在2015年1月27日發(fā)布的最新報告顯示，在全球范圍內(nèi)，感染埃博拉病毒病例為22057例，死亡人數(shù)上升至8795人；感染沙特冠狀病毒導(dǎo)致死亡人數(shù)也已達到400多人。因此，對這些病毒性傳染病的發(fā)病機理，傳染規(guī)律及防治策略的研究就顯得尤為重要。

病毒的傳播可以看作是服從某種規(guī)律的網(wǎng)絡(luò)傳播行為。21世紀初，人們利用節(jié)點和連接節(jié)點的連線，建立了可以描述傳染性疾病傳播的復(fù)雜網(wǎng)絡(luò)。最常見的均勻網(wǎng)絡(luò)是小世界網(wǎng)絡(luò)[1]，該網(wǎng)絡(luò)節(jié)點數(shù)不變，可以增加或減少少量的隨機捷徑，即每個節(jié)點的度都圍繞平均值上下波動；而無標度網(wǎng)絡(luò)[2]模型揭示了原有網(wǎng)絡(luò)的連線不變，增加新節(jié)點數(shù)目以及新節(jié)點與原有網(wǎng)絡(luò)結(jié)點的連線來改變原有網(wǎng)絡(luò)的拓撲結(jié)構(gòu)?？茖W(xué)家們發(fā)現(xiàn)有些病毒性傳染病的傳播是由于其周圍環(huán)境的不同（即能接觸到傳染病人的數(shù)量可能會相差很大）會導(dǎo)致被感染的概率不同，這種傳播特點符合無標度網(wǎng)絡(luò)的范疇。

近些年，將復(fù)雜網(wǎng)絡(luò)應(yīng)用到傳染病的研究已受到了人們的極大關(guān)注，例如文[3]在研究具有齊次性質(zhì)的susceptible-infected- removed （SIR）傳染病網(wǎng)絡(luò)系統(tǒng)時，發(fā)現(xiàn)當(dāng)傳染率比傳染臨界值大時，染病者會在人群中占有一定比例，當(dāng)傳染率小于傳染臨界值時，疾病最終會消失；文[4]在文[3]的基礎(chǔ)上考慮了通過接種疫苗來降低被感染的幾率，建立了具有免疫作用齊次性質(zhì)SIR模型的復(fù)雜網(wǎng)絡(luò)（所有結(jié)點的度都服從一個常數(shù)k的指數(shù)分布），該模型如下：

此外，文[4]還利用非齊次性質(zhì)（度分布不服從指數(shù)分布， 如冪分布）的復(fù)雜網(wǎng)絡(luò)（模型如（2）所示），討論了免疫作用以及非齊次性質(zhì)對疾病傳播的影響，說明了通過接種疫苗確實可以起到預(yù)防和控制疾病傳播的作用。

文[5]在無標度網(wǎng)絡(luò)上建立了一類SEI1I2R的傳染病模型，研究了無限維無標度網(wǎng)絡(luò)中的傳播閾值，證明了無病平衡點的穩(wěn)定性與疾病的一致持續(xù)生存性，得到了目標免疫好于比例免疫的結(jié)論。

文[6]研究了一類現(xiàn)實復(fù)雜情形下SIRS型傳染病模型的無病平衡點與地方病平衡點的唯一存在性，漸近穩(wěn)定性及傳染病綜合控制策略，本文在文[4-6]的基礎(chǔ)上，將種群分為易感個體S，染病個體I，康復(fù)個體R，考慮到非線性傳染率、群外個體遷入、生育與死亡以及疾病可水平和垂直傳播等因素，來研究病毒性傳染病傳播的規(guī)律。

1 具有免疫功能齊次性質(zhì)的無標度網(wǎng)絡(luò)

本文考慮會有新的易感個體輸入到此種群中，引入易感個體的輸入率As，種群內(nèi)會有出生率b∈[0，1]和死亡率d∈[0，1]，（d>b），慮考遺傳因素，假設(shè)康復(fù)個體R僅生育康復(fù)個體R，易感個體S僅生育易感個體S，感染病毒的人群在患病期間痛苦不堪，不會考慮生育后代，即假設(shè)染病者不會生育孩子（如麻風(fēng)病）。用λ表示傳染病通過有效接觸進行傳播的傳播系數(shù)，k為單位時間內(nèi)接觸的人數(shù)。γ為治愈率，γ∈[0，1]?？紤]到染病者與易感者、康復(fù)者的死亡率可能不同，因此用d-∈[0，1]表示染病個體的死亡率。a為免疫接種的比例，則可以得到如下的動力學(xué)模型：

2 具有免疫功能非齊次性質(zhì)的無標度網(wǎng)絡(luò)

在復(fù)雜網(wǎng)絡(luò)中，無標度網(wǎng)絡(luò)可以極好地演繹病毒性疾病的傳播情況，本文根據(jù)文[4]和文[5]的研究方法，建立一個SIR免疫無標度網(wǎng)絡(luò)模型來研究疫苗免疫對病毒性疾病在非齊次性質(zhì)的無標度網(wǎng)絡(luò)中傳播的影響。

現(xiàn)令Sk（t），Ik（t），Rk（t）分別為在t時刻易感染者、感染者、移出者在度為k的個體中占有的比例，Θ（t）表示任選一條連線指向感染者的概率。病毒傳播流程圖如圖1所示。

根據(jù)流程圖1，可得到如下所示的無標度網(wǎng)絡(luò)SIR免疫模型：

2.1 度分布不關(guān)聯(lián)

度分布不關(guān)聯(lián)時，任選網(wǎng)絡(luò)中一條連線指向感染者的概率如下：

3 傳染病控制方法與免疫策略

病毒性傳染病爆發(fā)的特點是來勢兇、危害大、傳播廣、社會影響大。人類一直與各種傳染性疾病斗爭。隨著交通運輸水平的大幅提升，染病者可以在短時間內(nèi)到達許多地方，若沒有做好預(yù)防措施，危害不可小視。因此，傳染病的預(yù)防與控制在現(xiàn)代生活中顯得尤為重要。

通過免疫方法來抑制病毒傳播是人們所關(guān)心的話題，采用的免疫策略不同，得到的免疫效果也不同。下面根據(jù)前面所得到的理論分析，來比較一下均勻比例免疫策略和目標比例免疫策略的效果。

3.1 免疫策略的比較分析

本文第1部分所研究的齊次性質(zhì)復(fù)雜網(wǎng)絡(luò)模型（3）所描述的是在均勻網(wǎng)絡(luò)上采用均勻比例免疫的數(shù)學(xué)模型，即指從網(wǎng)絡(luò)中隨機的選取一定比例的節(jié)點進行免疫，對網(wǎng)絡(luò)中度大的節(jié)點和度小的節(jié)點平等對待，沒有考慮網(wǎng)絡(luò)中節(jié)點間的差異。

本文第2部分研究的非齊次性質(zhì)的無標度網(wǎng)絡(luò)模型（7）所描述的是目標比例免疫，即網(wǎng)絡(luò)中各節(jié)點的度不均勻，可以對這些節(jié)點采取有目標的免疫，選取少量度較大的節(jié)點進行免疫，而在復(fù)雜網(wǎng)絡(luò)中，一旦該節(jié)點被免疫，則可以把與這些節(jié)點連接的邊刪除，其余節(jié)點不會傳染免疫節(jié)點，被免疫節(jié)點更不會傳染其他節(jié)點。被刪去的節(jié)點都是網(wǎng)絡(luò)中度較大的，刪去的邊也多，故大大減少了病毒的傳播途徑。在無標度網(wǎng)絡(luò)中，目標比例免疫的免疫臨界值與無標度網(wǎng)絡(luò)中新節(jié)點加入新連接的邊數(shù)相關(guān)。傳播率在很大范圍內(nèi)取不同值時，與均勻比例免疫相比，可以得到較小的免疫臨界值。相對于均勻比例免疫，目標比例免疫在現(xiàn)實實行起來也更容易[8]。

為更直觀描述這兩種免疫策略的免疫效果，根據(jù)本文對傳染臨界值λc，λ-c，λ～c算出的結(jié)果及在文[9-10]的基礎(chǔ)上，得到的結(jié)果如表1所示。

3.2 傳染病控制其他方法

考慮到現(xiàn)實世界中的傳染病傳播[11-16]和控制過程[17-20]中存在下列情形：①疾病可通過水平和垂直傳播；②群內(nèi)個體會生育；③群內(nèi)個體有死亡（自然和非自然）；④對染病個體采取隔離措施但不完全到位；⑤群外易感個體遷入。根據(jù)本文模型中的主要參數(shù)，針對流行于人類的傳染病的特點，經(jīng)過分析可以從如下幾個方面來控制病毒性傳染病的傳播：

1）接種疫苗：從本文計算出的傳染臨界值λ-c的表達式可以看出，若接種比例a變大，傳染臨界值變大。此時疾病需要跨越更大的傳染臨界值才能流行，這說明通過接種疫苗可以預(yù)防和控制疾病的傳播。

2）提高治愈率：從本文計算出的傳染臨界值λ-c的表達式還可以看出，治愈率γ變大，傳染臨界值也可以變大，所以可以通過提高醫(yī)學(xué)水平可以提高治愈率，從而控制疾病的傳播。

3）增大隔離率：同樣從本文計算出的傳染臨界值λ-c的表達式可以看出，引入易感個體的輸入率As越小，傳染臨值就越大，也就是說通過控制人口流動，且對染病者進行有效隔離，可以減少易感者與感染者的接觸率，減少易感者感染病毒的幾率。

另外，傳染病的流行需要合適的媒介，對生態(tài)環(huán)境進行改善，消除細菌賴以生存的因素，這些都對控制與預(yù)防流行性傳染性疾病具有重要意義。

4 結(jié) 論

本文討論了無標度網(wǎng)絡(luò)上SIR傳染病模型平衡點的存在性及其穩(wěn)定性的，得到了一些在無標度網(wǎng)絡(luò)上含有免疫SIR模型的動力學(xué)行為。在研究無標度網(wǎng)絡(luò)具有非齊次性質(zhì)時，分別討論度分布關(guān)聯(lián)和度分布不關(guān)聯(lián)時的傳染臨界值。并根據(jù)傳染臨界值提出了控制傳染病的方法，比如注射疫苗，提高治愈率，控制流動人口數(shù)量等。為有效抑制病毒傳染提供了一定的理論依據(jù)。從均勻比例免疫和目標比例免疫的差異比較中，發(fā)現(xiàn)目標比例免疫策略優(yōu)于均勻比例免疫。并根據(jù)得到的染病臨界值提出了控制病毒傳染綜合策略。

參 考 文 獻：

[1]史定華.無標度網(wǎng)絡(luò)基礎(chǔ)理論和應(yīng)用研究[J].電子科技大學(xué)學(xué)報，2010， 39（5）：644-650.

[2]黃萍，張許杰，劉剛.小世界網(wǎng)絡(luò)的研究現(xiàn)狀與展望[J].情報雜志， 2007（4）：66-68.

[3]MORENO Y， PASTOR-SATORRAS R， VESPIGANI A. Epidemic out Breaks in Complex Heterogeneous Networks [J]. European Phys J B， 2002（26）： 521-529.

[4]張海峰，傅新楚.含有免疫作用的SIR傳染病模型在復(fù)雜網(wǎng)絡(luò)上的動力學(xué)行為[J].上海大學(xué)學(xué)報（自然科學(xué)版）， 2007，13（2）：189-192.

[5]杜雯，黃堯，婁潔.生殖器疣在無標度網(wǎng)絡(luò)上的傳播與免疫策略研究[J].生物數(shù)學(xué)學(xué)報，2014，29（2）：361-368.

[6]孫有發(fā)，郭旭沖，梁肖肖等.現(xiàn)實復(fù)雜情形下的SIRS型傳染病模型及其控制策略[J].系統(tǒng)仿真學(xué)報，2010，22（1）：195-200.

[7]王文娟.無疾病潛伏期的傳染病動力學(xué)模型的閾值和再生數(shù)分析[J].運城學(xué)院學(xué)報，2005，23（5）：12-14.