摘要: 校園是一個特殊的人群密集場所,如果校園內(nèi)出現(xiàn)新冠病例往往有大面積傳染的風險,還可能成為家庭傳播的重要來源。為此提出了基于班級減少接觸策略,用于校園疫情防控,該策略可理解為依據(jù)班級結(jié)構對校園人群進行社團劃分后,減少社團間或社團內(nèi)的人際接觸。研究結(jié)果表明,在減少同等接觸人數(shù)或接觸時長的前提下,基于社團結(jié)構的有差別減少接觸策略可取得更好的疫情防控效果,使得校園人群的感染峰值和致病總?cè)藬?shù)進一步降低。
關鍵詞: COVID19;疾病傳播;社團結(jié)構;SEIR模型;校園疫情
中圖分類號: TP391;R183.3文獻標識碼: A
Suppression of COVID19 Campus Spreading Based on Community Structures
GAO Tian, XU Xiaoke
(College of Information and Communication Engineering, Dalian Minzu University, Dalian 116600, China)
Abstract:The campus is a special crowded place. If a patient with covid19 appears on campus, there is often the risk of large-scale infection, and it may also become an important source of family transmission. This paper proposes a class-based contact reduction strategy for campus epidemic prevention and control. This strategy can be understood as reducing human contact between and within communities after the campus population is divided into groups according to the class structure. The research results show that, on the premise of reducing the same number of contacts or contact time, differential contact reduction strategies based on community structure can achieve better epidemic prevention and control effects, making the infection peak and the total number of disease-causing population fell further on campus.
Keywords: COVID19; disease spread; community structure; SEIR model; campus outbreak
0 引言
自2019年底新冠疫情爆發(fā)以來,中國政府迅速行動,提出了“外防輸入,內(nèi)防反彈”總策略和“動態(tài)清零”總方針,取得了世界矚目的疫情防控成果。但隨著后期阿爾法、德爾塔、奧米克戎等傳染性更強變種的出現(xiàn)[1],使得國內(nèi)疫情態(tài)勢變得嚴峻起來,不但家庭聚集性疫情頻發(fā),而且校園聚集性疫情也多了起來[24]。在此形勢下,校園疫情防控越來越成為一個熱點話題,如何保證疫情期間中小學正常的教學活動成為亟待解決的社會問題。
當新冠肺炎病例在校園內(nèi)出現(xiàn)時,需要及時施加有效的防控措施。由于各種藥物性干預措施都無法完全阻止疫情傳播,因此各種非藥物性干預措施,如群體免疫[56]、數(shù)字跟蹤[79]、對所有人實行無差別的減少接觸[1013]等措施被提出。相比群體免疫策略和數(shù)字跟蹤策略,無差別減少接觸策略是目前非藥物干預措施中使用最為廣泛的防控策略,通過降低人群之間的接觸來延緩或者阻斷新冠傳播。如Laura等通過無差別減少25%、75%、95%的物理接觸措施,發(fā)現(xiàn)三種減少比例下總感染人數(shù)均會下降,且隨著減少比例的增大總感染人數(shù)也會快速下降[10]。Petrnio等提出了一種新的基于covid-SEIR的傳染病模型,減少30%~90%(以10%遞增)的接觸比例,結(jié)果表明減少接觸比例越大感染總?cè)藬?shù)也越少[11]。此外,Calistus等[12]和Saif等[13]的研究,通過社區(qū)接觸率參數(shù)基線值降低的方式來衡量無差別減少接觸措施的影響。雖然上述無差別減少接觸措施在防控疫情上有一定效果,但考慮到新冠病毒變體的高傳播性及在校園這種封閉環(huán)境下極快的傳播速度,很多時候?qū)W校只能選擇停課封校來緩解校園疫情,而停課隔離的大量學生返回家庭后會帶來極大的經(jīng)濟和社會成本,甚至不利于整個社會的穩(wěn)定[14]。針對此問題,需要制定更加科學有效的防控措施去應對新冠肺炎在校園的傳播。
在一些信息傳播研究中,研究人員會基于網(wǎng)絡的社團結(jié)構制定策略來控制社會輿論傳播[15],隱私信息傳播[16]等,受此啟發(fā),在校園新冠傳播研究中也可考慮基于社團結(jié)構來抑制傳播。由于校園一般存在明顯的班級結(jié)構,校園人群在同一班級內(nèi)的接觸較多,在不同班級間的接觸則較少,且校園的班級信息易于獲取,所以不同于信息傳播研究中使用社團檢測算法來檢測網(wǎng)絡的社團結(jié)構,在校園新冠傳播研究中可考慮依據(jù)班級結(jié)構來進行社團劃分?;诎嗉壗Y(jié)構進行社團劃分時,將同屬一個班級的人劃分到一個社團,然后減少社團間或社團內(nèi)的人際接觸來控制校園新冠傳播。
新冠通過真實的線下人際接觸進行傳播,而線下人際接觸網(wǎng)絡往往具有社團結(jié)構,基于此,本文首先通過帶有社區(qū)結(jié)構的真實線下人際接觸數(shù)據(jù)來構建真實校園新冠傳播環(huán)境即新冠傳播模型,然后將減少接觸分為減少接觸人數(shù)與減少接觸時長兩種手段在這個新冠傳播模型下仿真模擬進行研究。相對于無差別減少接觸策略,結(jié)果表明考慮減少接觸人數(shù)時,減少社團內(nèi)的接觸防控效果最好;而考慮減少接觸時長時,減少社團間的接觸防控效果最好。當社團間和社團內(nèi)減少接觸聯(lián)合使用時,針對減少接觸人數(shù)指標,先減少社團內(nèi)再減少社團間的接觸防控效果最好;而針對減少接觸時長指標,先減少社團間再減少社團內(nèi)的接觸防控效果最好。本研究考慮了班級具有操作性的社團結(jié)構在新冠防控中的作用,能夠為制定有效的校園疫情防控策略提供參考,更好的平衡校園疫情防控和學生正常學習秩序。
1 研究數(shù)據(jù)與方法
1.1 真實人際接觸數(shù)據(jù)和時效網(wǎng)絡的構建
由于技術和數(shù)據(jù)采集的限制,以前很多研究中使用網(wǎng)絡模型、線上社交網(wǎng)絡或者問卷調(diào)查數(shù)據(jù)構建了人工合成傳播環(huán)境進行新冠或其他疾病傳播研究[1719],并沒有使用人際之間的真實接觸數(shù)據(jù)。
近些年隨著信息科技和物聯(lián)網(wǎng)技術的快速發(fā)展,已經(jīng)能使用高科技手段采集適用于疾病傳播研究的真實線下接觸數(shù)據(jù),比如Wi-Fi熱點數(shù)據(jù)[20]、GPS數(shù)據(jù)[21]、定制射頻標簽(RFID)數(shù)據(jù)[2223]等。基于Wi-Fi熱點在采集時可能會因為系統(tǒng)設定或者電量問題不受控的關閉熱點從而導致線下接觸數(shù)據(jù)丟失,GPS設備則會面臨定位精度不準確而導致數(shù)據(jù)精確度較差。采集效果更好的是通過定制的射頻標簽(RFID),在采集人群接觸信息時,個體可通過佩戴定制的RFID傳感器以20s的時間分辨率檢測并記錄個體間的近距離接觸,最終收集的數(shù)據(jù)既包含接觸邊信息(可用于計算個體的接觸人數(shù)),也包含接觸次數(shù)信息(可用于計算個體間的接觸時長)。本文構建新冠傳播環(huán)境時使用了RFID設備收集的網(wǎng)絡規(guī)模較大的小學真實線下接觸數(shù)據(jù)集[22]和高中真實線下接觸數(shù)據(jù)集[23]。
真實小學數(shù)據(jù)集記錄了242名學生在校內(nèi)2天的線下接觸數(shù)據(jù),高中數(shù)據(jù)集記錄了329名學生在校內(nèi)5天的線下接觸信息。其中小學人群共有11個班級125 773條接觸數(shù)據(jù)(班級內(nèi)部107 599條,班級間18 174條),高中人群共有9個班級188 508條接觸數(shù)據(jù)(班級內(nèi)部173 558條,班級間14 950條),每條數(shù)據(jù)代表一次近距離接觸。這種近距離接觸會存在相同的人之間多次接觸的情況,通過去除這類重復的接觸數(shù)據(jù)獲得接觸邊總數(shù)的數(shù)據(jù),用于表示數(shù)據(jù)集網(wǎng)絡中人與人之間共有多少接觸邊,小學人群共有8 317條接觸邊(班級內(nèi)部3 769條,班級間4 548條),高中人群共有5 818條接觸邊(班級內(nèi)部3 578條,班級間2 240條)。
考慮到校園環(huán)境下,老師同學一般都周一到周五在學校工作學習,周末居家休息,因此本研究基于上述兩個數(shù)據(jù)集構建以周為周期的人際接觸時效網(wǎng)絡作為新冠傳播環(huán)境,而后在這個新冠傳播環(huán)境下進行仿真實驗。對于小學數(shù)據(jù)集來說,原始數(shù)據(jù)只有2天非雙休日的接觸數(shù)據(jù)(day1,day2),所以本文基于2天的原始數(shù)據(jù)將非雙休日接觸數(shù)據(jù)擴充為5天(周一至周五分別為day1,day2,day1,day2,day1),另外雙休日校園人群一般不會在校園內(nèi)聚集,故仿真實驗中不構建在雙休日的接觸網(wǎng)絡,認為個體在周末相互不接觸。高中數(shù)據(jù)集雖然有5天非工作日的接觸數(shù)據(jù),但第1天數(shù)據(jù)并不完整,只有5h的數(shù)據(jù),本文在構建網(wǎng)絡時使用第5天數(shù)據(jù)替代了第1天數(shù)據(jù)并在后續(xù)的仿真中驗證了替代數(shù)據(jù)并不會影響數(shù)據(jù)在時間上的關聯(lián),雙休日的處理同小學數(shù)據(jù)集一樣。
構建好的時效網(wǎng)絡中包含校園人群的接觸信息,這些接觸信息既可以算出每個人的接觸人數(shù),也可以根據(jù)記錄的接觸信息條數(shù)來計算校園人群人與人之間的接觸時長。由于原接觸數(shù)據(jù)集是每20s收集一次接觸信息的結(jié)果,其中會存在接觸不連續(xù)的情況,如果將40s或更大的時長定義為一次連續(xù)接觸的時長,那么將會出現(xiàn)沖突,因此本文將20s定義為一次接觸的接觸時長。以圖1為例,a與b存在兩次近距離接觸,與d存在一次近距離接觸,那么a的接觸人數(shù)為2,a與b的接觸時長為2*20s=40s,a與d的接觸時長為1*20s=20s。同理,也可以算出b、c、d的接觸人數(shù)以及與他人的接觸時長。
1.2 傳播模型
常用的傳染病模型有SI[24]、SIS[25]、SIR[26]、SEIR[27]等模型,每種模型的傳染病應用場景不同。SI模型適用于一旦患病就難以恢復的傳染病研究;SIS模型適用于患病后可治愈,但治愈后有再次感染疾病可能的傳染病研究;SIR模型適用于患病后可恢復為免疫態(tài),難以再次感染病毒的傳染病研究;SEIR模型可以看作增加了潛伏期(E)的SIR模型,潛伏期內(nèi)無任何感染癥狀,一旦過了潛伏期就會變?yōu)橛邪Y狀的感染態(tài),但最終可以恢復為免疫態(tài)。
新冠傳染屬于有潛伏期的傳染病,且治療好后在當前傳染期中再次感染的情況非常罕見[28],故可以視為在經(jīng)歷過感染態(tài)后,不會再次成為易感態(tài),只會成為免疫態(tài)的傳染病。因此針對新冠疫情的模擬仿真,使用SEIR模型最為合適。個體在SEIR模型中共有四種狀態(tài):易感態(tài)(S)、潛伏態(tài)或者暴露態(tài)(E)、感染態(tài)(I)、免疫態(tài)(R)。初始情況下,所有人都處于易感態(tài),然后隨機生成一個感染新冠的患者,感染態(tài)的患者通過人際接觸感染處于易感態(tài)的人;處于易感態(tài)的人一旦被感染,就會進入潛伏期,潛伏期內(nèi)沒有癥狀但是會有傳染性;在患者渡過潛伏期后,進入感染態(tài),感染態(tài)同樣有傳染性;最后在經(jīng)歷痊愈期后患者將會被治愈,成為免疫態(tài)。SEIR模型的感染率、潛伏期和痊愈期參數(shù)的具體設置如下:
1)感染率:基本繁殖數(shù)R0是傳染病學的一個重要指標,用于衡量一個群體中通過初始感染的個體直接感染的人數(shù)[29]。傳統(tǒng)的SEIR模型的感染率設置可依據(jù)傳染病的R0來設定平均感染率β-作為最終使用的仿真模擬參數(shù),初始時設置任意感染率,然后仿真模擬1 000次,通過得出的初始感染的個體直接感染的人數(shù)大小與文獻查到的基本繁殖數(shù)R0的比值乘上初始設置的任意感染率得出最終的平均感染率β-。雖然可以通過上述過程得出平均感染率β-,然而恒定的感染率參數(shù)并不適用于新冠病毒的傳播。根據(jù)Barry等的研究,動物在感染新冠后每日的病毒排放量En并不相等,而是隨著時間推移有先快速增加后緩慢下降的趨勢(day0-day14病毒排放量依次大致為0,10,158,501,1 995,1 500,1 000,658,316,160,3.98,3,2,1,0,單位:equivalent/ml),病毒排放量越大傳染性就越強[30]。本研究依據(jù)上述動物每天病毒排放量大小En來近似患者每一天的傳染率。首先計算感染期間病毒脫落量對數(shù)(以10為底)的平均值,計算公式:
E-=∑Nn=1lgEnN(1)
其中,N為被感染期間具有傳染能力的總天數(shù)(13天),En為患者在被感染后第n天時的病毒脫落量。通過E-可以計算每日傳染率,公式:
βn=β-·lgEnE-(2)
其中,βn為第n天的感染率,β-為被感染期間的平均感染率?,F(xiàn)有研究表明,新冠病毒傳播的R0值在2.24到5.71之間[3132],本文將小學和高中環(huán)境下的R0取值為4,之后在高中和小學環(huán)境下設置不同的感染率進行仿真模擬以計算R0值,最終仿真得出當平均感染率為2.8*10-4s-1時可滿足R0為4的條件,而后基于β-和病毒每日脫落量來計算每日的感染率。
2)潛伏期時長:據(jù)鐘南山團隊的研究,感染新冠后患者處于潛伏狀態(tài)天數(shù)的中位數(shù)為4天[33],故本文將仿真模型潛伏期時長參數(shù)設置為4天。
3)痊愈期時長:患者渡過潛伏期后進入感染態(tài),進入感染態(tài)的患者會被立即送入醫(yī)院隔離治療,不再接觸其他人員,在醫(yī)院治療的平均住院時間為10天[34]。
1.3 衡量指標
為了衡量不同減少接觸策略對于新冠疫情防控效果的影響,本文采用感染峰值和致病總?cè)藬?shù)這兩個指標,它們都是值越小說明對應的策略越具有優(yōu)勢。
1)感染峰值。感染峰值是指截至第i日處于潛伏或感染態(tài)的總?cè)藬?shù)Ni(i= …,n)中的最大值。使用感染峰值作為指標是衡量防控措施有效性的常見指標,感染峰值的降低會減輕醫(yī)療系統(tǒng)的瞬時壓力,可驗證防控策略的有效性。
2)致病總?cè)藬?shù)。致病總?cè)藬?shù)是指從疫情爆發(fā)到疫情結(jié)束總的感染人數(shù),致病總?cè)藬?shù)的降低可減少社會在治療患病人員的成本方面的投入。
2 基于社團結(jié)構減少接觸對校園疫情影響
2.1 減少接觸簡介
1)減少接觸人數(shù):此手段可理解為有些人能不見就不見,直接斷開與這些人的連接,比如減少接觸前a與b交互2次,減少接觸人數(shù)后兩者1次都不交互。對應到圖1的時效網(wǎng)絡上,減少a與b的接觸時,將會直接斷開a與b之間的接觸,相應地,a與b的接觸時長也降為0。
2)減少接觸時長:此手段可理解為減少個人之間接觸的總時長,比如減少接觸前a與b有2次交互,相應的接觸時長為2*20s=40s,現(xiàn)在減少為只有1次交互,接觸時長變?yōu)榱?*20s=20s。對應到圖1的時效網(wǎng)絡上,減少a與b的接觸時,將會使交互次數(shù)變?yōu)?次,但由于仍存在接觸,所以他們的接觸人數(shù)未發(fā)生變化。值得注意的是,當a與b的交互次數(shù)減為0時,接觸時長會變?yōu)?*20s=0s,此時相當于斷開了a與b之間的接觸,他們的接觸人數(shù)會各自減1。
2.2 基于社團結(jié)構減少接觸
新冠疫情在校園人群中通過人際接觸網(wǎng)絡進行傳播。一般情況下,校園環(huán)境下都存在明顯的班級結(jié)構,人們在班級內(nèi)的接觸人數(shù)較多、接觸時長較長,在班級間的接觸人數(shù)較少、接觸時長較短,因此可考慮校園人群的班級結(jié)構來改進無差別減少接觸策略。有可能減少班級間的接觸控制疫情傳播效果比較明顯,因為班級間的接觸斷開代表著大的學生群體間不再接觸。有可能減少班級內(nèi)的接觸控制疫情傳播效果更好,因為班級內(nèi)的學生接觸比較密切??紤]班級結(jié)構的減少接觸策略可看作首先依據(jù)班級結(jié)構對校園人群進行社團劃分,將同一班級的人看作一個社團,然后基于劃分出的社團結(jié)構來減少接觸。
基于社團結(jié)構的減少接觸策略可以分為減少社團內(nèi)接觸、減少社團間接觸兩種手段。實施策略首先要根據(jù)班級結(jié)構對校園人群進行社團劃分,而后將校園人群分為社團間人群和社團內(nèi)人群,分別用于減少社團內(nèi)接觸、減少社團間接觸的仿真實驗,同時在校園人群上進行無差別減少接觸的仿真實驗作為對比。考慮到社團間和社團內(nèi)人群的總接觸人數(shù)或總接觸時長不相同,為了對比三種減少接觸方式在同等減少接觸人數(shù)數(shù)量或接觸時長下的效果,所以還需設置一下減少接觸的上限。本節(jié)將減少上限設為社團內(nèi)人群和社團間人群的總接觸人數(shù)或總接觸時長的最小值,減少接觸人數(shù)上限用p表示,減少接觸時長上限用t表示。依照此規(guī)則,當減少接觸人數(shù)時,小學人群減少上限p設置為3 769人,高中人群減少上限p設置為2 240人,當減少接觸時長時,小學人群減少上限t設置為18 174*20s,高中人群減少上限t設置為14 950*20s。
通過上述介紹,首先確定了小學人群和高中人群減少接觸的上限值,并劃定了社團內(nèi)人群、社團間人群,之后將分別實施三種減少接觸方式的仿真實驗:1)社團內(nèi)減少:在社團內(nèi)人群中隨機選擇要減少接觸的人數(shù)或時長,進行基于社團結(jié)構的減少社團內(nèi)接觸的仿真實驗;2)社團間減少:在社團間人群中隨機選擇要減少接觸的人數(shù)或時長,進行基于社團結(jié)構的減少社團間接觸的仿真實驗;3)無差別減少:在校園人群中隨機選擇要減少接觸的人數(shù)或時長,進行無差別減少接觸的仿真實驗。
基于上述三種減少接觸方式,分別進行減少0.0p-1.0p(以0.1遞增)的接觸人數(shù)或減少0.0t-1.0t的(以0.1遞增)接觸時長的仿真實驗,減少0.0p或0.0t對應無減少接觸改動的情況。根據(jù)仿真實驗結(jié)果,本文比較兩種衡量指標的結(jié)果,一種表示當前感染峰值隨減少接觸人數(shù)或接觸時長增大而變化的情況,另一種表示每種減少不同接觸人數(shù)或接觸時長下致病總?cè)藬?shù)的變化情況。所有仿真實驗都在Windows10,python3.8的環(huán)境下進行了1 000次的重復模擬以排除偶然因素的影響。
2.3 減少接觸人數(shù)對校園疫情傳播影響
減少接觸是校園面臨疫情傳播威脅時可采取的防控措施,首先對減少校內(nèi)人群的接觸人數(shù)進行探究,設置小學人群減少上限p為3 769人,高中人群減少上限p為2 240人。圖2代表小學人群(左)和高中人群(右)在減少不同接觸人數(shù)下使用不同減少接觸方式的感染峰值變化結(jié)果。
由圖2可知,當減少接觸人數(shù)時,與無改動措施相比,雖然社團間、社團內(nèi)或者無差別減少接觸均會降低感染峰值,但減少社團內(nèi)接觸人數(shù)峰值降低幅度是最大的。分析產(chǎn)生這一結(jié)果的主要原因如下。當減少人群接觸人數(shù)的同時,也會相應減少人群的接觸時長,而人與人之間的接觸時長越短疫情越難以傳播。社團內(nèi)人群之間的平均接觸時長相比社團間人群或者整個校園人群的平均值都要長,所以當減少相同數(shù)量的接觸人數(shù)時,優(yōu)先社團內(nèi)減少更能抑制疫情的傳播,減低疫情傳播的峰值。
圖3代表小學人群和高中人群在減少不同接觸人數(shù)下使用不同減少接觸方式的致病總?cè)藬?shù)變化結(jié)果。從小學人群結(jié)果圖中可以看出,隨減少接觸人數(shù)的增大,減少社團內(nèi)接觸對降低致病總?cè)藬?shù)的效果會越來越顯著,另外兩種減少接觸方式未能有效降低致病總?cè)藬?shù)。高中人群的結(jié)果表明,3種減少接觸方式均會有效降低致病總?cè)藬?shù),但社團內(nèi)減少的效果最好。我們分析產(chǎn)生圖3這一結(jié)果的原因與產(chǎn)生圖2結(jié)果的原因一致,都是因為社團內(nèi)人群的平均接觸時長最大,所以當減少相同數(shù)量的接觸人數(shù)時,優(yōu)先社團內(nèi)減少的效果最優(yōu)。
2.4 減少接觸時長對校園疫情傳播影響
2.3節(jié)中分析了減少接觸人數(shù)對校園疫情的傳播影響,本小節(jié)將分析減少接觸時長對校園新冠疫情的傳播影響,在本小節(jié)中小學人群減少接觸時長上限t設置為18 174*20s,高中人群減少上限t設置為14 950*20s。小學人群和高中人群減少接觸時長的感染峰值變化結(jié)果如圖4所示。
從圖4中的曲線可以看出,無論是在小學人群還是高中人群,相比于其他兩種減少接觸方式,減少社團間接觸時長的仿真結(jié)果感染峰值最低,且隨著減少接觸時長的增大,這個效果越來越明顯。但同時也可以發(fā)現(xiàn)當使用減少接觸時長的手段時,社團內(nèi)減少或無差別減少的方式在感染峰值上變化不大,本文對這一現(xiàn)象分析的原因如下。社團內(nèi)人群或整個校園人群的平均接觸時長遠大于社團間人群的平均接觸時長,所以減少相同的接觸時長對社團間人群的影響比較大,對社團內(nèi)人群的影響比較小。由于社團間人群受影響較大,使得病毒難以從一個社團傳播到另一個社團,最終只能在一個或少數(shù)幾個受影響的社團內(nèi)進行傳播,從而感染峰值下降的效果明顯。與此同時,由于社團內(nèi)人群或整個校園人群的平均接觸時長較大,所以減少社團內(nèi)或無差別減少都未能對新冠在社團內(nèi)人群或整個校園人群的傳播產(chǎn)生顯著影響。
小學人群和高中人群減少不同接觸時長下致病總?cè)藬?shù)的變化如圖5所示。由5可知,對于社團內(nèi)減少或無差別減少方式來說,致病總?cè)藬?shù)隨減少接觸時長增大變化不大,這與圖4的結(jié)論相同。但不同的是,小學人群在減少社團間接觸方式的結(jié)果上存在一個0.8t的臨界點,減少接觸時長未超過0.8t時,致病總?cè)藬?shù)未有明顯變化,超過0.8t時,致病總?cè)藬?shù)開始快速下降。我們分析存在0.8t臨界點的原因如下。當小學人群減少接觸時長不超過0.8t時,疫情仍能輕易的從一個社團傳播到其他社團,使得減少社團間的效果并不明顯。當減少接觸時長超過0.8t時,社團間疫情傳播開始被有效抑制,從而使得致病總?cè)藬?shù)不斷下降。
3 基于社團結(jié)構減少接觸人數(shù)和接觸時長實驗擴展
上節(jié)在校園環(huán)境下的減少接觸仿真實驗可看作只減少社團間接觸或只減少社團內(nèi)接觸進行校園疫情防控進行研究,有時候并不一定能完全達到徹底消滅疫情傳播的效果。在疫情比較嚴重或者基本再生數(shù)比較大的情況下,需要將減少社團間接觸與減少社團內(nèi)接觸結(jié)合起來使用用來校園疫情防控,這種方式可以減少更多的接觸,同時達到更好的防止疾病傳播的效果。此時,我們探究是先減少社團間接觸再減少社團內(nèi)接觸效果好,還是先減少社團內(nèi)接觸再減少社團間接觸效果好。
擴展實驗將減少接觸上限提高到整個校園人群的總接觸人數(shù)(小學為8 317人,高中為5 818人)和總接觸時長(小學為125 773*20s,高中為188 508*20s)。仿真實驗的思路:1)先減少社團內(nèi)接觸,再減少社團間接觸:優(yōu)先減少社團內(nèi)人群的接觸人數(shù)或接觸時長,當社團內(nèi)的接觸人數(shù)或接觸p時長減為0時,再減少社團間的接觸人數(shù)或接觸時長。2)先減少社團間接觸,再減少社團內(nèi)接觸:優(yōu)先減少社團間人群的接觸人數(shù)或接觸時長,當社團間的接觸人數(shù)或接觸時長減為0時,再減少社團內(nèi)的接觸人數(shù)或接觸時長。3)整個校園人群無差別減少接觸:在整個校園人群中無差別減少接觸人數(shù)或接觸時長。
在本節(jié)的實驗中,使用減少接觸人數(shù)的手段時,設置小學人群減少人數(shù)上限p為8 317人,高中人群減少人數(shù)上限p為5 818人,致病總?cè)藬?shù)隨減少接觸人數(shù)增大的變化情況見圖6。使用減少接觸時長的手段時,設置小學人群減少時長上限t為125 773*20s,高中人群減少時長上限t為188 508*20s,致病總?cè)藬?shù)隨減少接觸時長增大的變化情況見圖7。兩圖中豎虛線對應X軸的位置為上節(jié)減少接觸人數(shù)或減少接觸時長對比實驗的上限值點,減少接觸人數(shù)或減少接觸時長大于這個上限值點時,可以觀察到先減少所有社團間的接觸再減少社團內(nèi)的接觸、先減少所有社團內(nèi)的接觸再減少社團間的接觸、無差別減少接觸三者的防控疫情傳播效果。另外因為減少1.0p或者1.0t的接觸時將會斷開網(wǎng)絡中所有的接觸,此時采取任何措施都不能使得疫情傳染開來,所以本文并未展示減少1.0p或者1.0t的結(jié)果。
由圖6可知,使用減少接觸人數(shù)的手段時,相比于另外兩種策略,先減少社團內(nèi)再減少社團間接觸效果最優(yōu)。對于小學人群,當減少接觸人數(shù)為0.5p時,先減少社團內(nèi)再減少社團間接觸人數(shù)的策略已經(jīng)可以將致病總?cè)藬?shù)控制在10人以下,有效控制了疫情在小學人群中的傳播。對于高中人群,當減少接觸人數(shù)為0.3p時,先減少社團內(nèi)再減少社團間接觸人數(shù)的策略同樣將致病總?cè)藬?shù)控制在10人以下,阻斷了疫情在高中人群中的傳播。另外從圖中的豎虛線可以看出,使用先減少社團內(nèi)再減少社團間的策略,當斷開所有社團內(nèi)的接觸時,已經(jīng)可以阻斷校園疫情的傳播,無需再進一步減少社團間的接觸。
由圖7可知,使用減少接觸時長的手段時,相比于另外兩種策略,先減少社團間再減少社團內(nèi)接觸效果最優(yōu)。對于小學人群和高中人群,當減少接觸時長分別為0.2t和0.1t時,先減少社團間再減少社團內(nèi)的策略已經(jīng)極大降低了致病總?cè)藬?shù),在此基礎上,如果繼續(xù)減小接觸,致病總?cè)藬?shù)會進一步下降。最后從圖中的豎虛線可以看出,當斷開社團間的所有接觸時,雖已經(jīng)極大降低了致病總?cè)藬?shù),但并沒有完全徹底消滅校園疫情,如果想達到更好的疫情防控效果可考慮進一步減少社團內(nèi)的接觸。
通過本節(jié)的仿真實驗可以得出結(jié)論,當同時減少社團內(nèi)和社團間的接觸時,單考慮減少接觸人數(shù),先減少社團內(nèi)的接觸再減少社團間的接觸控制校園疫情傳播效果最好;單考慮減少接觸時長,先減少社團間的接觸再減少社團內(nèi)的接觸控制校園疫情傳播效果最好。在實際應用中,當校園面臨疫情威脅,可以先采用單獨減少社團內(nèi)接觸人數(shù)或單獨減少社團間接觸時長的策略來預防校園疫情大規(guī)模爆發(fā)。當校園已經(jīng)發(fā)現(xiàn)疫情,應當立即采取更嚴格的防控措施,同時使用社團內(nèi)減少和社團間減少策略來控制新冠疫情傳播。
4 結(jié)語
本文基于真實校園中人際接觸數(shù)據(jù)構建了小學和高中新冠傳播環(huán)境,使用符合真實新冠傳染率規(guī)律的SEIR模型仿真模擬了無差別減少接觸策略與基于社團結(jié)構減少接觸策略在校園內(nèi)的新冠傳播情況,對各種減少接觸策略的仿真結(jié)果進行了對比分析。
結(jié)果表明,當使用減少校園人群接觸人數(shù)手段時,優(yōu)先減少社團內(nèi)的接觸在降低校園人群致病數(shù)和感染峰值上效果最好;當使用減少校園人群接觸時長手段時,優(yōu)先減少社團間的接觸效果最好。最后的擴展實驗表明,當社團間和社團內(nèi)減少接觸同時使用時,單考慮減少接觸人數(shù),先減少社團內(nèi)再減少社團間的接觸防控效果最好;單考慮減少接觸時長,先減少社團間再減少社團內(nèi)的接觸防控效果最好。
綜上所述,當校園人群面臨新冠傳播的威脅時,在防止校內(nèi)傳染源出現(xiàn)的同時,也應當在校內(nèi)實行合適的防控策略,較優(yōu)的策略是考慮校園人群的班級結(jié)構,先對校園人群進行社團劃分,然后根據(jù)實際情況優(yōu)先減少校園人群社團內(nèi)的接觸人數(shù)或社團間的接觸時長,實行基于社團結(jié)構的減少接觸策略。在此策略下,一旦病毒在校內(nèi)開始傳播,相比無差別減少接觸策略,基于社團結(jié)構的減少接觸策略可以取得更好的疫情防控效果。當然這些只是校園層面采取的防控措施,當面臨疫情威脅時,也要減少與社會人員的不必要接觸,做好校園人群的防護工作。
參考文獻:
[1]HEMMER C J, LBERMANN M, REISINGER E C. COVID19: epidemiologie und mutationen[J]. Der Radiologe, 202 61(10): 880887.
[2]溫競?cè)A, 宋晨. 當前疫情形勢如何?校園防疫如何做?[N]. 新華每日電訊, 20220507(4).
WEN J, SONG C. What is the current situation of the epidemic? How to do campus epidemic prevention?[N]. Xinhua Daily Telegraph, 20220507(4).
[3]陳晶. 落細防控舉措,抓好社區(qū)和校園疫情防控[N]. 人民政協(xié)報, 20220511(5).
CHEN J. Implement detailed prevention and control measures, and focus on community and campus epidemic prevention and control[N]. CPPCC DAILY, 20220511(5).
[4]郭麗菲. 太原學院疫情防控不松懈,筑牢校園“安全網(wǎng)”[N]. 山西青年報, 20220421(8).
GUO L. Taiyuan University remains vigilant in epidemic prevention and control, building a strong campus "safety net"[N]. Shanxi Youth Daily, 20220421(8).
[5]HABIB H. Has Sweden’s controversial covid19 strategy been successful?[J]. Bmj, 2020, 369: m2376.
[6]貝殼, 陸默. “群體免疫” 策略: 科學還是賭博[J]. 世界科學, 2020 (5): 1012.
BEI K, gdIuE1plgzkltuVEwuDg3w==LU M. "Group immunity" strategy: science or gambling[J]. World Science, 2020(5): 1012.
[7]FERRETTI L, WYMANT C, KENDALL M, et al. Quantifying SARS-CoV2 transmission suggests epidemic control with digital contact tracing[J]. Science, 2020, 368(6491): eabb6936.
[8]BAKER A, BIAZZO I, BRAUNSTEIN A, et al. Epidemic mitigation by statistical inference from contact tracing data[J]. Proceedings of the National Academy of Sciences, 202 118(32): e2106548118.
[9]SOLDANO G J, FRAIRE J A, FINOCHIETTO J M, et al. COVID19 mitigation by digital contact tracing and contact prevention (app-based social exposure warnings)[J]. Scientific Reports, 202 11(1): 18.
[10] MATRAJT L, LEUNG T. Evaluating the effectiveness of social distancing interventions to delay or flatten the epidemic curve of coronavirus disease[J]. Emerging Infectious Diseases, 2020, 26(8): 17401748.
[11] SILVA P C L, BATISTA P V C, LIMA H S, et al. COVID-ABS: an agent-based model of COVID19 epidemic to simulate health and economic effects of social distancing interventions[J]. Chaos, Solitons & Fractals, 2020, 139: 110088.
[12] NGONGHALA C N, IBOI E, EIKENBERRY S, et al. Mathematical assessment of the impact of non-pharmaceutical interventions on curtailing the 2019 novel Coronavirus[J]. Mathematical Biosciences, 2020, 325: 108364.
[13] ULLAH S, KHAN M A. Modeling the impact of non-pharmaceutical interventions on the dynamics of novel coronavirus with optimal control analysis with a case study[J]. Chaos, Solitons & Fractals, 2020, 139: 110075.
[14] EICHENBAUM M S, REBELO S, TRABANDT M. The macroeconomics of epidemics[J]. The Review of Financial Studies, 202 34(11): 51495187.
[15] SUN J, HAO Y. Hierarchicalsemantic community detection in information networks: a complete information graph approach[J]. Tehnicˇki Vjesnik, 2019, 26(6): 17001706.
[16] 張長倫, 寧楠, 金占勇. 基于雙重網(wǎng)絡社團結(jié)構的隱私信息傳播模型[J]. 信息技術, 2018(9): 3843.
ZHANG C, NING N, JIN Z. Privacy information spreading model based on community structure on dual netword[J]. Information Technology, 2018(9): 3843.
[17] MATRAJT L, LEUNG T. Evaluating the effectiveness of social distancing interventions to delay or flatten the epidemic curve of coronavirus disease[J]. Emerging Infectious Diseases, 2020, 26(8): 1740.
[18] SALATH M, JONES J H. Dynamics and control of diseases in networks with community structure[J]. PLoS computational biology, 2010, 6(4): e1000736
[19] ZHANG J, JIN Z. Epidemic spreading on complex networks with community structure[J]. Applied Mathematics and Computation, 2012, 219(6): 28292838.
[20] ZHOU M, MA M, ZHANG Y, et al. EDUM: classroom education measurements via large-scale WiFi networks[C]//Proceedings of the 2016 acm international joint conference on pervasive and ubiquitous computing. Heidelberg Germany, 2016: 316327.
[21] WANG R, CHEN F,CHEN Z, et al. StudentLife: assessing mental health, academic performance and behavioral trends of college students using smartphones[C]//Proceedings of the 2014 ACM International Joint Conference on Pervasive and Ubiquitous Computing. Seattle Washington, 2014: 314.
[22] STEHL J, VOIRIN N, BARRAT A, et al. High-resolution measurements of face-to-face contact patterns in a primary school[J].PloS One, 201 6(8): e23176.
[23] MASTRANDREA R, FOURNET J, BARRAT A. Contact patterns in a high school: a comparison between data collected using wearable sensors, contact diaries and friendship surveys[J].PloS One, 2015, 10(9): e0136497.
[24] BARTHLEMY M, BARRAT A, PASTOR-SATORRAS R, et al. Velocity and hierarchical spread of epidemic outbreaks in scale-free networks[J]. PhysicalReview Lletters, 2004, 92(17): 178701.
[25] QIN Y, ZHONG X, JIANG H, et al. An environment aware epidemic spreading model and immune strategy in complex networks[J]. Applied Mathematics and Computation, 2015, 261: 206215.
[26] SALATH M, JONES J H. Dynamics and control of diseases in networks with community structure[J].PLoS Computational Biology, 2010, 6(4): e1000736.
[27] 范如國, 王奕博, 羅明, 等. 基于SEIR的新冠肺炎傳播模型及拐點預測分析[J]. 電子科技大學學報, 2020, 49(3): 369374.
FAN R, WANG Y, LUO M, et al. SEIR-Based COVID19 Transmission Model and Inflection Point Prediction Analysis[J]. Journal of University of Electronic Science and Technology of China, 2020, 49(3): 369374.
[28] PEGHIN M, BOUZA E, FABRIS M, et al. Low risk of reinfections and relation with serological response after recovery from the first wave of COVID19[J]. European Journal of Clinical Microbiology & Infectious Diseases, 202 40(12): 25972604.
[29] WALLINGA J, TEUNIS P. Different epidemic curves for severe acute respiratory syndrome reveal similar impacts of control measures[J]. American Journal of Epidemiology,2004,160(6):50916.
[30] ROCKX B, KUIKEN T, HERFST S, et al. Comparative pathogenesis of COVID19, MERS, and SARS in a nonhuman primate model[J]. Science, 2020, 368(6494): 10121015.
[31] ZHAO S, LIN Q, RAN J, et al. Preliminary estimation of the basic reproduction number of novel coronavirus (2019-nCoV) in China, from 2019 to 2020: A data-driven analysis in the early phase of the outbreak[J]. InternationalJournal of Infectious Diseases, 2020, 92: 214217.
[32] SANCHE S, LIN Y T, XU C, et al. High contagiousness and rapid spread of severe acute respiratory syndrome coronavirus 2[J]. Emerging Infectious Diseases, 2020, 26(7): 14701477.
[33] GUAN W, NI Z, HU Y, et al. Clinical characteristics of coronavirus disease 2019 in China[J]. New EnglandJournal of Medicine, 2020, 382(18): 17081720.
[34] WANG D, HU B, HU C, et al. Clinical characteristics of 138 hospitalized patients with 2019 novel coronavirus-infected pneumonia in Wuhan, China[J]. Jama, 2020, 323(11): 10611069.
(責任編輯 李 進)