2005—2018年玉林市流行性腮腺炎流行特征和時(shí)空分析

(廣西玉林市疾病預(yù)防控制中心，廣西 玉林 537000)

流行性腮腺炎簡稱流腮，是由腮腺炎病毒感染引起的急性呼吸道傳染病，以腮腺非化膿性腫脹、疼痛伴發(fā)熱為主要癥狀，可導(dǎo)致多種腺體組織、神經(jīng)系統(tǒng)和肝臟等器官損傷，可引起腦膜炎、睪丸炎、卵巢炎等并發(fā)癥，對兒童危險(xiǎn)極大[1]。因此，各地及時(shí)分析性腮腺炎流行特征，掌握疫情動態(tài)，對開展防病宣教，減少發(fā)病非常重要。當(dāng)前對流腮的分析主要集中在三間分布，而三間分布分析不僅容易忽視不同時(shí)期和地區(qū)傳染病的高發(fā)具有偶然性，而且還存在人為選擇偏倚[2]。因此，有必要在行三間分布分析的同時(shí)，結(jié)合地理信息系統(tǒng)(geographic information system，GIS)分析傳染病的流行特征，以解決三間分布分析的不足。目前GIS在公共衛(wèi)生領(lǐng)域的應(yīng)用十分廣泛，通過GIS可分析疾病流行的空間特征和分布模式，探索可能的影響因素及病因，為疾病流行的監(jiān)控、預(yù)警、防制效果評價(jià)等提供參考依據(jù)[3]。因此，本研究在對玉林市流腮流行特征行三間分布分析的同時(shí)，利用GIS探索其時(shí)空聚集性，客觀定量地對疫情進(jìn)行分析，為玉林市各縣(市、區(qū))科學(xué)地制定流腮的防控策略提供參考依據(jù)。

1 資料與方法

1.1 資料來源 2005—2018年玉林市流腮發(fā)病數(shù)據(jù)來自中國疾病預(yù)防控制信息系統(tǒng)中的傳染病報(bào)告信息管理系統(tǒng)，人口資料來源于基本信息系統(tǒng)；地圖來源于2015年1∶100萬全國基礎(chǔ)地理信息數(shù)據(jù)。

1.2 分析方法

1.2.1 流腮流行特征分析 采用描述流行病學(xué)方法對玉林市2005—2018年流腮病例資料進(jìn)行三間分布分析，其中年齡分組采用傳染病報(bào)告信息管理系統(tǒng)的年齡分組形式。利用Microsoft Excel 2010對病例數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析。

1.2.2 空間自相關(guān)分析 采用ArcGIS 10.2軟件對中國矢量地圖進(jìn)行裁剪得到玉林市以縣(市、區(qū))為界的矢量地圖，通過ArcMap將發(fā)病數(shù)據(jù)進(jìn)行空間可視化展示。采用GeoDa 1.14軟件對2005—2018年玉林市流腮病例資料進(jìn)行空間自相關(guān)分析，分別以全局Moran’s I與局部空間關(guān)聯(lián)指數(shù)(local indicators of spatial association，LISA)為指標(biāo)進(jìn)行全局空間自相關(guān)分析和局部空間自相關(guān)分析[4]。全局空間自相關(guān)分析從全市角度分析是否存在流腮空間聚集，局部空間自相關(guān)分析則進(jìn)一步分析具體聚集縣(市、區(qū))及聚集模式。Moran’s I的取值范圍為 [-1，1]，若P<0.05且Moran’s I值>0，表示研究變量存在空間正相關(guān)；若P<0.05且Moran’s I值<0，表示研究變量存在空間負(fù)相關(guān)；Moran’s I值接近于0表示空間隨機(jī)分布，不存在相關(guān)性。采用蒙特卡羅模擬方法對 Moran’s I值進(jìn)行檢驗(yàn)，計(jì)算得到Z值和P值，以P<0.05為差異有統(tǒng)計(jì)學(xué)意義。通過LISA分布圖可獲得5種局部空間自相關(guān)分析模式：不顯著、高-高聚集模式、低-低聚集模式、低-高聚集模式和高-低聚集模式。

1.2.3 時(shí)空掃描分析 采用SaTScan 9.5軟件進(jìn)行時(shí)空掃描統(tǒng)計(jì)分析。選用離散變量Poisson分布模型，用圓形移動窗口法創(chuàng)建掃描窗口，掃描時(shí)間范圍設(shè)置為“2005—2018 年”，時(shí)間間隔“1年”，地區(qū)重疊設(shè)置為“無”，最大聚集范圍窗口為總?cè)丝诘?0%[5]。選擇蒙特卡羅模擬次數(shù)“999”，計(jì)算在不同圓心和不同半徑下的動態(tài)窗口區(qū)域內(nèi)外空間單元屬性的對數(shù)似然比 (log likelihood ratio，LLR)進(jìn)行統(tǒng)計(jì)推斷，以探測聚集中心和聚集半徑。根據(jù)LLR的大小確定聚集區(qū)類型，若P<0.05，則認(rèn)為該區(qū)域存在聚集性[2，6]。

2 結(jié)果

2.1 2005—2018年玉林市流腮發(fā)病的時(shí)間分布特征 2005—2018年玉林市共報(bào)告流腮10 798例，年均報(bào)告發(fā)病率為13.80/10萬，發(fā)病率最高的年份是2012年(23.12/10萬，1278例)，最低的是2017年(5.47/10萬，315例)；2005—2009年發(fā)病率呈下降趨勢，2010—2012年發(fā)病率回升并達(dá)到最高峰，2013—2015年發(fā)病率又大幅下降，2016—2018年發(fā)病率波動較大，但總體是上升趨勢。因此，2005—2018年玉林市流腮發(fā)病呈現(xiàn)“降-升-降-升”模式(見表1、圖1)。按季節(jié)和月份分析，發(fā)現(xiàn)每年全年各月均有病例報(bào)告，但存在明顯的季節(jié)性，4～7月為發(fā)病主高峰，11月至次年1月為發(fā)病次高峰，呈雙峰型分布，其中以2010—2013年期間發(fā)病高峰最密集。見圖2。

表1 2005—2018年玉林市流腮發(fā)病情況

圖1 2005—2018年玉林市流行性腮腺炎報(bào)告病例數(shù)和年均發(fā)病率

圖2 2005—2018年玉林市流行性腮腺炎月發(fā)病曲線

2.2 2005—2018年玉林市流腮發(fā)病的地區(qū)分布特征 玉林市各縣(市、區(qū))均有病例報(bào)告，玉州區(qū)年均發(fā)病率最高且報(bào)告病例數(shù)最多(28.77/10萬，2596例)，福綿區(qū)報(bào)告病例數(shù)最少(641例)，但博白縣年均發(fā)病率最低(5.08/10萬)，見表2。

表2 2005—2018年玉林市各縣(市、區(qū))流腮發(fā)病順位統(tǒng)計(jì)表

2.3 2005—2018年玉林市流腮發(fā)病的人群分布特征 按年齡分析，病例主要集中在3～14歲(占81.77%)，其中5～7歲年齡組發(fā)病率最高，分別為87.37/10萬、96.71/10萬、93.34/10萬，見表3。按性別分析，發(fā)現(xiàn)共報(bào)告男性病例7181例(占66.50%)，女性病例3617例(占33.50%)，男女年均發(fā)病率分別為17.51/10萬、9.72/10萬，男性發(fā)病率高于女性，各年男女發(fā)病率比在(1.43～2.14)∶1之間，見表1。按職業(yè)分析，發(fā)病以中小學(xué)生為主(6088例，占56.38%)，其次是幼托兒童(2151例，占19.92%)，散居兒童(1612例，占14.93%)，農(nóng)民(591例，占5.47%)，其他職業(yè)者(356例，占3.30%)。

表3 2005—2018年玉林市不同年齡組流腮發(fā)病情況

2.4 2005—2018年玉林市流腮疫情的GIS分析

2.4.1 空間自相關(guān)分析 采用GeoDa軟件，生成Rook鄰接的空間矩陣，分析2005—2018年玉林市流腮的Global Moran’s I與Local Moran’s I指數(shù)，可見除2017年(Moran’s I=0.336，Z=2.255，P=0.020)存在空間正相關(guān)外，其余年份均不存在空間自相關(guān)性，即無空間聚集性；LISA聚集分布圖顯示2017年流腮發(fā)病率呈高-高聚集模式，聚集區(qū)為北流市，見表4。

表4 2005—2018年玉林市流腮發(fā)病率空間自相關(guān)分析結(jié)果

2.4.2 時(shí)空掃描分析 時(shí)空掃描分析顯示，2005—2018年玉林市流腮的空間分布是隨機(jī)的，存在1個(gè)主要聚集區(qū)和4個(gè)次要聚集區(qū)，2016年北流市存在一個(gè)主要時(shí)空聚集區(qū)(LLR=331.22，RR=3.55，P<0.001)；次要聚集區(qū)分別為2018年博白縣、2010—2013年興業(yè)縣和福綿區(qū)、2008—2009年容縣、2005—2007年玉州區(qū)，結(jié)果見表5、圖3。

表5 2005—2018年玉林市流腮發(fā)病率時(shí)空掃描統(tǒng)計(jì)結(jié)果

圖3 2005—2018年玉林市流腮時(shí)空掃描聚集性可視化

3 討論

2005—2018年玉林市流腮年均報(bào)告發(fā)病率(13.80/10萬)低于同期全國(21.44/10 萬)[7]和廣西(37.92/10萬)[8]年均發(fā)病率水平，其中2012年為玉林市流腮發(fā)病高峰年(同期，全國亦出現(xiàn)流腮高發(fā)[7，9])。玉林市流腮發(fā)病曲線與謝群等[10]報(bào)道的流腮每隔2～5年會出現(xiàn)一次較大規(guī)模的流行的結(jié)論相一致，當(dāng)易感者積累到一定數(shù)量時(shí)，新的發(fā)病高峰便會再次出現(xiàn)。玉林市流腮發(fā)病呈現(xiàn)明顯的季節(jié)性，與相關(guān)報(bào)道一致[8，10]。流行曲線與學(xué)校開學(xué)和放假時(shí)間安排相吻合，每個(gè)學(xué)期開學(xué)后就會逐漸出現(xiàn)一次發(fā)病高峰。提示在學(xué)校和托幼機(jī)構(gòu)出現(xiàn)病例后若未能及時(shí)采取有效的隔離和治療措施，則很容易發(fā)生交叉感染而導(dǎo)致暴發(fā)。

從地區(qū)分布來看，玉林市玉州區(qū)年均發(fā)病率最高，可能原因是玉州區(qū)作為玉林市政治、經(jīng)濟(jì)、文化的中心，流動人口多且流動頻繁，進(jìn)城務(wù)工人員的子女免疫接種不規(guī)范，易感人群聚集，感染和傳播傳染病風(fēng)險(xiǎn)高，因此更易使城區(qū)流腮高發(fā)。分析發(fā)現(xiàn)，玉林市流腮發(fā)病年齡主要集中在3～14歲，尤以5～7歲兒童發(fā)病率最高，與全國流腮發(fā)病情況相似[7]，男性發(fā)病率高于女性，與既往文獻(xiàn)結(jié)果一致[2，8，10]。也說明流腮的易感人群為幼兒和中小學(xué)生[11]，其中又以男生為甚。這可能與男生衛(wèi)生意識差，天性好動，感染機(jī)會大等有關(guān)。

通過采用GIS進(jìn)行空間自相關(guān)分析和時(shí)空掃描分析，可從時(shí)間和空間兩個(gè)方面同時(shí)探討疾病在地理空間上聚集區(qū)域隨時(shí)間變化的規(guī)律性，并獲取聚集區(qū)域的相對危險(xiǎn)度，更準(zhǔn)確地定位聚集區(qū)域的位置[6]?？臻g自相關(guān)分析結(jié)果顯示，2005—2016年和2018年玉林市流腮發(fā)病不存在空間聚集性，可能是流腮在玉林市是隨機(jī)發(fā)病的，不存在空間自相關(guān)；在實(shí)際研究中，空間自相關(guān)是不均勻的，而Moran’s I指數(shù)僅反映空間整體自相關(guān)，且個(gè)別局部對象的屬性取值對全局分析影響顯著[12]。時(shí)空掃描分析結(jié)果顯示存在5個(gè)聚集區(qū)，聚集時(shí)間與流腮月發(fā)病曲線高峰相同，結(jié)合地區(qū)分布可進(jìn)一步認(rèn)為玉林市流腮的發(fā)生是隨機(jī)分布的，考慮可能的原因是流腮自然流行周期引起隨機(jī)散發(fā)、人口流動密集使易感者增多引起爆發(fā)、疫情監(jiān)測和疫苗接種工作的不斷完善使各地發(fā)病率得到有效控制，具體原因有待進(jìn)一步探索。

根據(jù)玉林市流腮流行特征分析和時(shí)空分析結(jié)果。建議：一是對兒童要做好兩劑次含腮腺炎成分疫苗(MuCV)的接種工作，兩劑接種時(shí)間間隔13～36個(gè)月，可最大限度減低轄區(qū)流腮發(fā)病率[13]。二是要做好學(xué)校和托幼機(jī)構(gòu)入托入學(xué)兒童腮腺炎疫苗接種證查驗(yàn)，及時(shí)發(fā)現(xiàn)漏種兒童，做好補(bǔ)種。三是要學(xué)校和托幼機(jī)構(gòu)做好晨午檢及因病缺勤登記，及時(shí)發(fā)現(xiàn)和處置流腮疫情，防止擴(kuò)散蔓延。四是要加大流腮防控知識宣教力度，特別是在流腮漏報(bào)率高的地區(qū)，更應(yīng)提高人群預(yù)防意識[14]；同時(shí)應(yīng)加強(qiáng)醫(yī)療衛(wèi)生機(jī)構(gòu)包括CDC的專業(yè)技術(shù)人員的繼續(xù)醫(yī)學(xué)教育，以提高他們的傳染病防控知識，因?yàn)閭魅静》揽刂R欠缺，對傳染病認(rèn)識不足，缺乏應(yīng)有的警惕性等都與醫(yī)學(xué)教育密切相關(guān)[15]。