考慮個體警覺行為的雙層網(wǎng)絡(luò)流行病傳播模型

王 帥，宋玉蓉，2，宋 波

（1.南京郵電大學(xué)自動化學(xué)院，南京 210003；2.江蘇省物聯(lián)網(wǎng)智能機器人工程實驗室，南京 210003；3.南京郵電大學(xué)現(xiàn)代郵政學(xué)院，南京 210003）

0 概述

近年來，流行病傳播在復(fù)雜網(wǎng)絡(luò)領(lǐng)域受到廣泛關(guān)注［1］。人類社會中流行病傳播過程常伴隨信息傳播等其他傳播過程，流行病信息通過媒體在人群中迅速擴散，進而改變大眾對流行病的認識。人們在了解到流行病相關(guān)信息后，為避免被感染其行為也會發(fā)生相應(yīng)改變，如戴口罩、避免接觸以及服用藥物等［2］，這些行為在一定程度上可有效抑制流行病傳播［3-5］。因此，對流行病傳播與信息擴散關(guān)系的研究具有重要意義。

目前，流行病傳播與信息擴散的交互模型為雙層網(wǎng)絡(luò)［6-7］，一層為模擬流行病在物理接觸網(wǎng)絡(luò)中傳播過程的物理接觸網(wǎng)絡(luò)，另一層為描述流行病傳播所引發(fā)信息擴散的信息網(wǎng)絡(luò)。在物理接觸網(wǎng)絡(luò)中，通常采用易感-感染-易感（Susceptible-Infected-Susceptible，SIS）流行病傳播模型或易感-感染-恢復(fù)（Susceptible-Infected-Recovered，SIR）流行病傳播模型［8-10］模擬流行病的傳播。在SIS 模型中，個體分為易感（S）個體和感染（I）個體，易感個體在和感染個體接觸后以概率β變?yōu)楦腥緜€體，而感染個體則以概率μ自我恢復(fù)，重新成為易感個體。在SIR 模型中，個體除了SIS 模型中的S 狀態(tài)和I狀態(tài)外，還存在恢復(fù)（R）狀態(tài)。I狀態(tài)的個體不再重新變?yōu)镾 狀態(tài)，而是以概率γ恢復(fù)或死亡成為R 狀態(tài)。在信息網(wǎng)絡(luò)中，通常采用未知-已知-未知（Unaware-Aware-Unaware，UAU）信息傳播模型描述流行病傳播引起的信息擴散。其中，個體根據(jù)對信息的了解程度分為不了解信息的無意識（U）狀態(tài)和已了解信息的有意識（A）狀態(tài)，其傳播動力學(xué)與SIS 模型相同。

現(xiàn)有研究通常將流行病傳播與信息擴散之間的相互作用建模為多層網(wǎng)絡(luò)中兩個相互競爭的傳播過程。與物理接觸網(wǎng)絡(luò)中流行病傳播過程相同，意識信息從有意識個體傳給其在信息網(wǎng)絡(luò)上的無意識鄰居，通過促使更多個體在信息網(wǎng)絡(luò)中產(chǎn)生意識，從而使其能夠采取預(yù)防行為以抑制流行病在物理接觸網(wǎng)絡(luò)中傳播。文獻［6］提出UAU-SIS 網(wǎng)絡(luò)傳播模型以描述網(wǎng)絡(luò)中意識和流行病之間的動態(tài)相互作用，指出存在意識傳播速率亞臨界點，當(dāng)意識傳播速率小于該臨界點時，流行病閾值不受意識傳播影響。文獻［11］發(fā)現(xiàn)易感個體在了解信息后采取的自我意識行為對流行病傳播的抑制作用明顯優(yōu)于感染個體的自我意識行為，其原因在于易感個體的自我意識行為能降低其被感染的可能性。文獻［12］研究多層網(wǎng)絡(luò)UAU-SIR 模型后發(fā)現(xiàn)，隨著自我感知率v的增加，個體自我感知會極大降低流行病感染規(guī)模。文獻［13］在UAU-SIS 傳播模型的基礎(chǔ)上調(diào)整UAU 信息傳播過程，使個體在其具有意識的鄰居比例超過給定值時會產(chǎn)生自我意識，并指出流行病爆發(fā)與意識擴散及接觸網(wǎng)絡(luò)的拓撲結(jié)構(gòu)有關(guān)。文獻［14］在接觸網(wǎng)絡(luò)的SIR 傳播模型基礎(chǔ)上引入免疫狀態(tài)，個體可通過接種疫苗對流行病免疫，且信息網(wǎng)絡(luò)的同質(zhì)性能促進疫苗接種規(guī)模擴大，從而有效抑制流行病傳播。

以上研究均假設(shè)網(wǎng)絡(luò)中個體具有相同特性，而現(xiàn)實中個體行為通常受鄰居數(shù)量等異質(zhì)性影響。考慮到不同個體在網(wǎng)絡(luò)中鄰居關(guān)系的異質(zhì)性，文獻［15］提出一種非線性耦合信息—流行病模型，指出流行病閾值由耦合網(wǎng)絡(luò)的拓撲結(jié)構(gòu)決定，其中節(jié)點的度分布的不均勻性會降低流行病閾值，控制流行病傳播需增強對已了解該流行病信息個體的保護。文獻［16］通過研究具有異質(zhì)感染率的多重網(wǎng)絡(luò)流行病傳播過程，構(gòu)建可調(diào)冪指數(shù)與節(jié)點的度相關(guān)的異質(zhì)感染率函數(shù)，證明負冪指數(shù)會使流行病感染率降低，即異質(zhì)感染率與節(jié)點的度負相關(guān)，免疫樞紐節(jié)點有利于抑制流行病傳播。雖然上述研究考慮到個體在了解流行病傳播產(chǎn)生的信息后都會產(chǎn)生警覺意識，卻忽視了不同個體在接觸行為選擇上的異質(zhì)性，即不同個體在了解到流行病相關(guān)信息后會采取不同行為以避免感染流行病或者維持網(wǎng)絡(luò)功能正常運行。例如，在流行病傳播期間，一部分具有流行病信息意識的個體可能會提高自身警覺性，注意基本的健康防護或采取其他措施來保證自身健康。此外，其還可能對其他鄰居產(chǎn)生警覺性，通過主動切斷與部分鄰居的接觸來避免自身被感染。而另一部分個體雖然具有流行病信息的意識，但是由于其存在日常工作與社交需求，不會對周圍個體產(chǎn)生警覺性，無法切斷與部分鄰居的接觸。

考慮不同個體在接觸行為上的異質(zhì)性，文獻［17］在雙層網(wǎng)絡(luò)UAU-SIS 傳播模型中引入個體的積極狀態(tài)與非積極狀態(tài)，在流行病傳播過程中，積極狀態(tài)的個體在信息網(wǎng)絡(luò)和接觸網(wǎng)絡(luò)中會斷開與同為積極狀態(tài)個體的連接，并證明這種基于個體活躍性狀態(tài)的引入可有效降低流行病傳播規(guī)模。此外，積極狀態(tài)與非積極狀態(tài)的切換是隨機發(fā)生的，其與信息和流行病的傳播相互獨立，由于現(xiàn)實中意識信息的傳播會對個體行為產(chǎn)生影響，因此由意識信息傳播引發(fā)個體行為的變化更符合現(xiàn)實情況。本文受文獻［17］啟發(fā)，考慮到流行病傳播時個體接觸行為存在異質(zhì)性，即個體了解流行病信息后會以一定概率成為警覺個體，進而改變個體接觸行為以規(guī)避被流行病所感染，提出一種基于個體警覺狀態(tài)的UAU-SIS 雙層網(wǎng)絡(luò)傳播模型。建立物理接觸網(wǎng)絡(luò)和信息擴散網(wǎng)絡(luò)，設(shè)計兩種個體警覺性接觸行為策略，并研究其在3 種雙層網(wǎng)絡(luò)中對流行病傳播的影響。

1 流行病傳播與信息擴散交互模型

1.1 傳播模型

由于社交網(wǎng)絡(luò)中信息傳播和個人行為的改變是一種復(fù)雜的社會心理過程，因此本文使用UAU-SIS雙層網(wǎng)絡(luò)建立模型并簡化傳播機制，如圖1 所示。定義下層SIS 網(wǎng)絡(luò)為物理接觸網(wǎng)絡(luò)，用以描述流行病在人群中的傳播；節(jié)點表示個體，連邊表示人與人之間的物理接觸。定義上層UAU 網(wǎng)絡(luò)為信息擴散網(wǎng)絡(luò)，用以描述隨流行病傳播發(fā)生的信息擴散；節(jié)點表示個體，連邊表示個體之間的信息交流。兩層網(wǎng)絡(luò)的對應(yīng)節(jié)點為同一個體，由于相同個體在進行物理接觸和信息交流時其鄰居集合通常不相同，因此上下兩層網(wǎng)絡(luò)結(jié)構(gòu)不同。個體在信息網(wǎng)絡(luò)中成為已了解信息的有意識個體后會以概率a變成警覺性個體（警覺節(jié)點），否則為非警覺性個體（非警覺節(jié)點）。

圖1 UAU-SIS 雙層網(wǎng)絡(luò)傳播模型Fig.1 UAU-SIS bilayer network spreading model

在物理接觸網(wǎng)絡(luò)中用SIS 模型模擬流行病在網(wǎng)絡(luò)中的傳播，個體在每個時刻為易感狀態(tài)或感染狀態(tài)。在流行病傳播過程中，當(dāng)感染個體與易感個體接觸時，易感個體會以概率β被感染，同時感染個體也會以概率μ恢復(fù)為易感個體。在信息網(wǎng)絡(luò)中用UAU 模型模擬流行病傳播時的信息擴散，與SIS模型類似，個體在每個時刻均為不了解信息的無意識狀態(tài)或了解信息的有意識狀態(tài)。假設(shè)信息層有意識個體來源于兩個方面：1）接觸網(wǎng)絡(luò)中無意識的感染個體在信息網(wǎng)絡(luò)中以概率κ轉(zhuǎn)變?yōu)橛幸庾R個體，并自發(fā)地向鄰居個體傳播流行病相關(guān)信息；2）信息網(wǎng)絡(luò)中了解信息的有意識個體在每個時刻會向信息層中的鄰居傳播流行病相關(guān)信息，在與無意識鄰居接觸后以概率λ使其變?yōu)橛幸庾R個體。同時，因為大部分流行病的傳播都具有一定季節(jié)性或周期性，使得散播信息的有意識個體遺忘或者不再傳播流行病相關(guān)信息，所以有意識個體會以概率δ重新變?yōu)闊o意識個體。

在流行病相關(guān)信息的傳播過程中，由于不了解流行病相關(guān)信息的無意識個體不會考慮如何預(yù)防被流行病感染的問題，而了解流行病相關(guān)信息的有意識個體會采取適當(dāng)?shù)姆雷o措施以減少被感染的風(fēng)險，因此有意識的易感個體被流行病感染的概率小于無意識個體。在本文模型中，假設(shè)有意識的易感個體被流行病感染的概率為無意識易感個體的σA∈[0，1]倍，若無意識的易感個體被流行病感染的概率為β，則有意識的易感個體被流行病感染的概率為σAβ?？紤]個體接觸行為異質(zhì)性對流行病傳播的影響，本文在模型中引入基于個體警覺性的接觸行為，并提出兩種不同的個體接觸行為策略。

1.2 個體警覺行為

本文引入個體警覺狀態(tài)表示不同個體在接觸行為選擇上的異質(zhì)性，信息層中了解流行病相關(guān)信息的有意識個體會以概率a變?yōu)榫X狀態(tài)（V）。由于警覺性個體充分了解流行病的嚴重性，因此其愿意改變自身社交行為以避免感染。然而隨著對流行病相關(guān)信息的遺忘，警覺性個體會停止散播信息，并以概率b重新恢復(fù)為非警覺狀態(tài)（N）。非警覺性個體因為未認識到流行病的嚴重性或存在工作與社交需求，所以會與全部鄰居接觸。在流行病傳播期間，若個體更偏向于進入警覺狀態(tài)，則其通常會希望保持警覺以最大程度的避免感染風(fēng)險，因此假設(shè)a+b=1?；谏鲜瞿Ｐ?，本文根據(jù)個體產(chǎn)生警覺性時接觸行為的不同偏好，提出個體警覺行為1 和個體警覺行為2 兩種相反的接觸行為策略。

1.2.1 個體警覺行為1

具有警覺性的個體會和同樣具有警覺性的鄰居個體保持物理接觸，同時斷開與不具有警覺性鄰居個體的物理接觸，而不具有警覺性的個體則會和所有鄰居保持接觸。例如，在流行病傳播期間，了解相關(guān)信息的個體比平時更加謹慎并采取帶口罩、勤洗手等防護措施，少數(shù)充分了解流行病嚴重性的警覺性個體為最大程度地避免被流行病感染，還會切斷與不了解流行病嚴重性的非警覺性鄰居個體之間的物理接觸，僅與同樣具有警覺性的鄰居接觸，由于這些鄰居已了解流行病相關(guān)信息并采取防護措施，因此與其接觸帶來的感染風(fēng)險更小。

1.2.2 個體警覺行為2

具有警覺性的個體會和不具有警覺性的鄰居個體保持物理接觸，同時斷開與警覺性鄰居個體的物理接觸，而不具有警覺性的個體則會和所有鄰居保持接觸。例如，在流行病傳播期間有意識的個體來源于兩個方面：1）由已感染流行病的個體自發(fā)轉(zhuǎn)化而來；2）在信息傳播網(wǎng)絡(luò)中接收到鄰居個體散播的消息從而獲得意識。因此，當(dāng)一部分有意識的個體產(chǎn)生警覺性后，警覺性個體中會有一定比例的感染態(tài)個體，如果具有警覺性的易感個體和不具備警覺性的鄰居個體進行物理接觸，則被流行病感染的風(fēng)險會降低。

考慮到個體警覺行為1 和個體警覺行為2 兩種行為策略都有其合理性，本文采用這兩種策略建立模型，通過參數(shù)控制不同個體選取相應(yīng)警覺行為策略并分析其對流行病傳播的抑制效果。

1.3 數(shù)學(xué)描述

在傳播過程中，每個時刻都被分為意識傳播（UAU 過程）、警覺狀態(tài)的改變、流行病傳播（SIS 過程）以及自我意識的產(chǎn)生（信息上傳過程）4 個階段。結(jié)合式（1）～式（5），通過微觀馬爾科夫分析［17］建立每個節(jié)點i分別從8 種狀態(tài)演化的方程為：

2 仿真與結(jié)果分析

本文提出個體在產(chǎn)生警覺性后的兩種個體接觸行為策略：處于警覺狀態(tài)的個體傾向于與同樣處于警覺狀態(tài)的鄰居個體保持物理接觸并和處于非警覺狀態(tài)的鄰居個體斷開物理接觸；處于警覺狀態(tài)的個體傾向于與非警覺狀態(tài)的鄰居個體保持物理接觸并和處于警覺狀態(tài)的鄰居個體斷開物理接觸。由于處于非警覺狀態(tài)的個體不具有警覺性，因此其會和所有鄰居個體保持接觸。以下針對不同的多層網(wǎng)絡(luò)類型進行仿真實驗，對比分析兩種個體警覺行為對流行病傳播的影響。其中，流行病傳播以及信息擴散參數(shù)設(shè)置為β=0.2、μ=0.4、λ=0.3、δ=0.6、α=0.4 以及σS=0.5。每個仿真實驗均進行200 次蒙特卡洛循環(huán)［7-9］。

由于現(xiàn)實網(wǎng)絡(luò)結(jié)構(gòu)具有小世界特性［20］和無標(biāo)度特性［21］，因此本文實驗選取具有這兩種特性的小世界網(wǎng)絡(luò)（WS）和無標(biāo)度網(wǎng)絡(luò)（BA）分別表示物理接觸層和信息傳播層。小世界網(wǎng)絡(luò)又稱均勻網(wǎng)絡(luò)或指數(shù)網(wǎng)絡(luò)，其網(wǎng)絡(luò)連接度分布可近似用Poisson 分布表示，該網(wǎng)絡(luò)具有明顯的聚類特性和較短的平均路徑長度；無標(biāo)度網(wǎng)絡(luò)的連接度分布函數(shù)具有冪律形式，其無明顯特征長度且具有增長特性和優(yōu)先連接特性2 個重要特性。本文從物理接觸層為無標(biāo)度網(wǎng)絡(luò)且信息傳播層為無標(biāo)度網(wǎng)絡(luò)的BA-BA 雙層網(wǎng)絡(luò)、物理接觸層為無標(biāo)度網(wǎng)絡(luò)且信息傳播層為小世界網(wǎng)絡(luò)的BA-WS 雙層網(wǎng)絡(luò)、物理接觸層為小世界網(wǎng)絡(luò)且信息傳播層為小世界網(wǎng)絡(luò)的WS-WS 雙層網(wǎng)絡(luò)，分別研究考慮個體警覺行為的雙層網(wǎng)絡(luò)流行病傳播與信息擴散交互模型。由于物理接觸網(wǎng)絡(luò)結(jié)構(gòu)較符合BA 網(wǎng)絡(luò)模型的增長特性與優(yōu)先連接特性，而信息傳播網(wǎng)絡(luò)結(jié)構(gòu)可能具有無標(biāo)度特性，也可能隨著時間的演化變得相對均勻，即具有小世界特性，因此在上下兩層網(wǎng)絡(luò)不同的情況下，本文暫未考慮物理接觸層為小世界網(wǎng)絡(luò)且信息傳播層為無標(biāo)度網(wǎng)絡(luò)的WS-BA 雙層網(wǎng)絡(luò)。

本文實驗中BA 網(wǎng)絡(luò)的節(jié)點數(shù)為2 000，每個節(jié)點在加入網(wǎng)絡(luò)時與3 個節(jié)點建立連接，WS 網(wǎng)絡(luò)的節(jié)點數(shù)為2 000，平均度＜k＞=4，網(wǎng)絡(luò)在生成過程中邊的重連概率p=0.3。在信息傳播層，現(xiàn)實社會中相同個體在接觸關(guān)系網(wǎng)絡(luò)和在線社交關(guān)系網(wǎng)絡(luò)中雖然無相同結(jié)構(gòu)，但也不會產(chǎn)生差異化較大的人際關(guān)系。例如，在接觸關(guān)系網(wǎng)絡(luò)中擁有較多好友的個體，其在信息傳播層中也同樣可能擁有較多好友，即使在這兩層網(wǎng)絡(luò)中擁有的好友并不完全相同。因此，為保證這種相關(guān)性，在BA-BA 雙層網(wǎng)絡(luò)和WS-WS 雙層網(wǎng)絡(luò)中，上層信息傳播網(wǎng)絡(luò)以下層接觸網(wǎng)絡(luò)為基礎(chǔ)隨機增加400 條邊［7-9］，在BA-WS 雙層網(wǎng)絡(luò)中，對上層信息傳播WS 網(wǎng)絡(luò)也進行隨機加邊處理。

當(dāng)采取個體警覺行為1 時，BA-BA、BA-WS 和WS-WS 3 種雙層網(wǎng)絡(luò)中流行病傳播規(guī)模i（t）在不同個體警覺性概率a下隨模擬的步驟序數(shù)t的變化如圖2 所示。可以看出，隨著個體警覺性概率a的增大，網(wǎng)絡(luò)中流行病傳播規(guī)模也增大，即個體警覺行為1對于流行病抑制作用減弱。從易感個體角度考慮，具有警覺性的易感個體需受到保護，當(dāng)其產(chǎn)生警覺性后可通過調(diào)整接觸行為策略減小被感染的概率。例如，通過切斷與部分非警覺性感染鄰居個體的接觸，在一定程度上可降低易感個體被感染的風(fēng)險。

圖2 3 種雙層網(wǎng)絡(luò)采取個體警覺行為1 時流行病傳播規(guī)模隨步驟序數(shù)的變化Fig.2 The change of epidemic spreading scales with ordinal number of steps when three bilayer networks adopt individual alert behavior 1

當(dāng)a較小時，網(wǎng)絡(luò)中只有少數(shù)有意識個體處于警覺狀態(tài)，其中具有警覺性的易感個體會切斷與其他非警覺性感染態(tài)鄰居的接觸以降低自身感染概率。由于一部分個體意識來自信息擴散，由與有意識的鄰居通信所獲得，另一部分個體意識由個體自身感染流行病而自發(fā)產(chǎn)生，因此有意識的個體中患病個體的比例會比無意識個體中患病個體的比例更高。當(dāng)a偏小時，由于不具備警覺性的易感個體比例高，采取個體警覺行為1 可在一定程度上實現(xiàn)對具有警覺性個體的隔離，因此可保護不具備警覺性的易感個體。隨著a的增大，非警覺性概率b減小，即網(wǎng)絡(luò)中個體由警覺性恢復(fù)為非警覺性的概率降低，最終網(wǎng)絡(luò)中警覺性個體比例增大，警覺個體之間形成的子網(wǎng)絡(luò)擴大，在其子網(wǎng)內(nèi)部各個體都會與其全部鄰居進行物理接觸，造成個體警覺行為1 的效果減弱。當(dāng)a=1 時，網(wǎng)絡(luò)中所有具有意識的個體永久轉(zhuǎn)變?yōu)榫X性個體，此時警覺性個體形成的子網(wǎng)絡(luò)擴展到幾乎整個網(wǎng)絡(luò)，因此個體警覺行為1 基本不起作用。特別地，當(dāng)a較小時，在流行病傳播初期會出現(xiàn)一個較高峰值（見圖2中a=0.1 時的曲線峰值），這是由于網(wǎng)絡(luò)中最初警覺性個體比例很小，隨著流行病在網(wǎng)絡(luò)中不斷傳播，警覺性個體比例隨有意識個體數(shù)量的增加而提高，此時個體警覺行為1 逐漸產(chǎn)生作用，使流行病感染規(guī)模逐漸下降，且隨著時間延長傳播規(guī)模逐漸穩(wěn)定。

當(dāng)采取個體警覺行為2 時，BA-BA、BA-WS 和WS-WS 3 種雙層網(wǎng)絡(luò)中流行病傳播規(guī)模在不同個體警覺性概率a下隨模擬的步驟序數(shù)的變化如圖3 所示?？梢钥闯觯S著個體警覺性概率a增大，網(wǎng)絡(luò)中流行病傳播規(guī)模減小，即個體警覺行為2 對流行病抑制作用增強。從易感個體角度考慮，警覺性易感個體會選擇與非警覺性鄰居個體接觸，因為有意識個體中的部分個體產(chǎn)生意識是因為其在物理接觸網(wǎng)絡(luò)中感染流行病，所以與無意識個體相比，有意識個體感染比例略高。警覺性個體由有意識個體根據(jù)概率轉(zhuǎn)化而來，通過避免與同樣具有警覺性的個體物理接觸，在一定程度上可降低個體被感染的風(fēng)險。當(dāng)a較小時，網(wǎng)絡(luò)中只有少數(shù)有意識個體處于警覺狀態(tài)，其中警覺性易感個體會保持與非警覺性鄰居個體（包括大部分感染個體）的物理接觸，此時個體警覺行為2 對于流行病傳播的影響不顯著。隨著a的增大，有意識個體產(chǎn)生警覺性的概率增大，由于有意識個體中感染個體的比例較高，因此警覺性易感個體通過改變接觸行為，切斷與其他警覺性感染鄰居個體的物理接觸可有效降低感染風(fēng)險。當(dāng)a=1時，網(wǎng)絡(luò)中很多個體隨著信息擴散永久成為警覺性個體，而警覺性個體不會主動與警覺性鄰居個體接觸，因此，當(dāng)網(wǎng)絡(luò)中警覺性個體比例很高時，網(wǎng)絡(luò)中大部分連邊均被切斷，警覺個體之間相互孤立，此時網(wǎng)絡(luò)基本不連通，流行病無法有效傳播，隨著已有感染個體逐漸康復(fù)，流行病最終在網(wǎng)絡(luò)中消失。

圖3 3 種雙層網(wǎng)絡(luò)采取個體警覺行為2 時流行病傳播規(guī)模隨步驟序數(shù)的變化Fig.3 The change of epidemic spreading scales with ordinal number of steps when three bilayer networks adopt individual alert behavior 2

特別地，與個體警覺行為1 時情況類似，當(dāng)采取個體警覺行為2 時，若a較大則在流行病傳播初期會出現(xiàn)一個較高峰值（見圖3中a=1.0 時的曲線峰值），此時網(wǎng)絡(luò)中處于感染態(tài)的有意識個體比例偏高，在流行病傳播初期這些個體在產(chǎn)生警覺性后通過與大量非警覺性個體接觸傳播流行病，而這些非警覺性個體并無流行病相關(guān)意識且防護能力較弱，隨著流行病在網(wǎng)絡(luò)中不斷傳播，警覺性個體比例隨著有意識個體數(shù)量的增加而提高，流行病的傳播規(guī)模逐漸穩(wěn)定。

因為BA 網(wǎng)絡(luò)具有無標(biāo)度特性，所以在網(wǎng)絡(luò)中存在少數(shù)度較大的個體，當(dāng)其被流行病感染時會通過與大量鄰居個體接觸而加快流行病在網(wǎng)絡(luò)中的傳播，體現(xiàn)出BA 網(wǎng)絡(luò)的脆弱性。WS 網(wǎng)絡(luò)由于節(jié)點的度分布較均勻，其流行病傳播速度比同等規(guī)模的BA網(wǎng)絡(luò)更慢。因此，當(dāng)采用WS-WS 網(wǎng)絡(luò)時，個體警覺行為對流行病的抑制作用比采用BA-BA 網(wǎng)絡(luò)更明顯。當(dāng)采用BA-WS 網(wǎng)絡(luò)時，由于僅信息傳播層為WS 網(wǎng)絡(luò)，信息意識傳播速度較同等規(guī)模的BA 網(wǎng)絡(luò)更慢，使得上層信息傳播網(wǎng)絡(luò)中具有信息意識的個體比例較少，造成具有警覺性意識的個體比例也較少，因此，當(dāng)采取個體警覺行為1 時，若a較小，則BA-WS 網(wǎng)絡(luò)在傳播初期的峰值較BA-BA 雙層網(wǎng)絡(luò)更高（見圖2（a）和圖2（b））；當(dāng)采取個體警覺行為2時，若a較大，則BA-WS 網(wǎng)絡(luò)在傳播初期的峰值較BA-BA 雙層網(wǎng)絡(luò)更高（見圖3（a）和圖3（b））。

圖4～圖6 分別為BA-BA、BA-WS 和WS-WS 3 種雙層網(wǎng)絡(luò)中不同個體警覺行為對傳播閾值的影響。

圖4 BA-BA 網(wǎng)絡(luò)采用不同個體警覺行為對流行病閾值的影響Fig.4 The influences of BA-BA network with different individual alert behaviors on epidemic threshold

圖5 BA-WS網(wǎng)絡(luò)采用不同個體警覺行為對流行病閾值的影響Fig.5 The influences of BA-WS network with different individual alert behaviors on epidemic threshold

圖6 WS-WS網(wǎng)絡(luò)采用不同個體警覺行為對流行病閾值的影響Fig.6 The influences of WS-WS network with different individual alert behaviors on epidemic threshold

可以看出：在3 種雙層網(wǎng)絡(luò)中，當(dāng)采取個體警覺行為1 時，隨著有意識個體警覺性概率a的增大，流行病爆發(fā)閾值會降低；當(dāng)采取個體警覺行為2 時，隨著有意識個體警覺性概率a的減小，流行病爆發(fā)閾值會降低。本文提出的兩種基于個體警覺行為的接觸策略均可通過提高流行病爆發(fā)閾值有效抑制流行病爆發(fā)。由于BA 為非均勻網(wǎng)絡(luò)，其傳播閾值遠小于WS 網(wǎng)絡(luò)［22］，因此當(dāng)物理接觸層為BA 網(wǎng)絡(luò)時，其流行病傳播閾值較物理接觸層為WS 網(wǎng)絡(luò)更小。由圖4～圖6 還可以看出：當(dāng)采取個體警覺行為1 時，若a減小則流行病傳播規(guī)模減小，流行病傳播閾值升高，而隨著a增大，個體警覺行為1 對流行病傳播規(guī)模的影響減小，流行病閾值降低；當(dāng)采取個體行為2時，若a增大則流行病傳播規(guī)模減小，流行病傳播閾值升高，而隨著a減小，個體警覺行為2對流行病傳播規(guī)模的影響減小，流行病閾值降低。由上述分析可知，基于個體警覺行為的接觸策略能有效抑制流行病傳播。

3 結(jié)束語

考慮到流行病傳播時個體警覺性對其接觸行為的影響，本文提出一種基于個體警覺行為的雙層網(wǎng)絡(luò)傳播模型。分別建立下層物理接觸網(wǎng)絡(luò)和上層信息擴散網(wǎng)絡(luò)，設(shè)計警覺個體避免與非警覺個體接觸以及警覺個體避免與警覺個體接觸兩種接觸行為策略，并在3 種雙層網(wǎng)絡(luò)模型中研究上述行為策略對流行病傳播的影響。仿真結(jié)果表明，該模型中兩種基于個體警覺狀態(tài)的接觸行為策略能有效降低流行病傳播規(guī)模，提高流行病爆發(fā)閾值。后續(xù)將研究自適應(yīng)雙層網(wǎng)絡(luò)模型，通過引入節(jié)點自適應(yīng)斷邊重連策略，進一步探索信息擴散對流行病傳播的抑制作用。