張夏男，張洪偉，紀 穎，王秋芳，謝立璟，葉 萌

布魯氏菌病(簡稱布病)是由布魯氏菌引起的人獸共患的一種慢性傳染病[1]，是我國的乙類傳染病。布病的主要臨床表現(xiàn)是發(fā)熱、多汗、乏力，肌肉和關(guān)節(jié)疼痛等[2]，布病在急性期如治療不徹底，極易轉(zhuǎn)為慢性，病程長達數(shù)年甚至幾十年，甚至喪失勞動力[3]。我國大部分省市報告人間布魯氏菌病的病例數(shù)有升高趨勢[4-6]，本研究擬采用GeoDa-1.14.0.0軟件、集中度法及圓形分布法、構(gòu)成比對全國31個省市布病流行特征的變化趨勢進行分析。

1 資料與方法

1.1資料來源 2006-2017年布病疫情資料來自公共衛(wèi)生科學數(shù)據(jù)共享中心網(wǎng)站(http://www.phsciencedata.cn)，該網(wǎng)站提供的數(shù)據(jù)庫涵蓋了39種傳染病、慢性非傳染性疾病、健康危險因素、生命登記等內(nèi)容，其中傳染病收集的主要是自2004年傳染病網(wǎng)絡(luò)直報以來報告的全部傳染病數(shù)據(jù)，主要包括分地區(qū)、年齡組的發(fā)病人數(shù)、發(fā)病率、死亡人數(shù)和死亡率等。為了分析2006-2017年布病病例的空間和時間屬性，將每3年的新發(fā)病例數(shù)進行合并形成4個時間段的數(shù)據(jù)，第1時間段為2006-2008年，第2時間段為2009-2011年，第3時間段為2012-2014年，第4時間段為2015-2017年。每個時間段各省市的累積發(fā)病率以3年新發(fā)病例數(shù)的總和作為分子，各省3年平均人口數(shù)為分母計算。

1.2 分析方法

1.2.1地區(qū)分布及時空分布 將4個時間段的累積發(fā)病率(1/10萬)數(shù)據(jù)導(dǎo)入GeoDa-1.14.0.0，以全國基礎(chǔ)地理數(shù)據(jù)為地圖，建立Rook權(quán)重矩陣文件，呈現(xiàn)全國31個省市累積發(fā)病率的箱地圖，同時進行全局空間自相關(guān)分析和局部空間自相關(guān)分析。全局空間自相關(guān)分析使用全局Moran’sⅠ指數(shù)反映空間要素屬性值的聚集程度；局部空間自相關(guān)分析反映空間單元與相鄰單元疫情的相關(guān)程度，可用LISA集聚圖識別四類聚集區(qū)域：高-高聚集、高-低聚集、低-高聚集和高-低聚集。將局部Moran’sⅠ指數(shù)通過10%顯著性水平檢驗的區(qū)域分類顯示于LISA集聚圖中。

1.2.2季節(jié)性分析 采用集中度法的M值和圓形分布法的R值對4個時間段布病的季節(jié)性特點進行分析。M值和R值的計算公式及評定標準參見文獻[7-8]。評定標準如下：越接近1季節(jié)性越強，越接近0季節(jié)性越弱，0.9以上代表嚴格季節(jié)性，0.7～0.9代表很強的季節(jié)性，0.5～0.7代表明顯的季節(jié)性，0.3～0.5代表一定的季節(jié)性，0.3以下代表季節(jié)性較弱。

1.2.3年齡分布 利用構(gòu)成比分析各年齡段分布的變化情況，并用趨勢性χ2檢驗分析老少組(14歲及以下和65歲及以上)和非老少組(15～64歲組)在4個時間段的構(gòu)成比有無統(tǒng)計學差異，檢驗水準α=0.05。

2 結(jié) 果

2.1全國布病疫情現(xiàn)狀 2006-2017年全國布病疫情總體呈上升趨勢，年平均發(fā)病率為2.836/10萬，發(fā)病率最低為2006年(1.454/10萬)，最高為2014年(4.223/10萬)，4個時間段的年平均發(fā)病率分別是1.686/10萬、2.691/10萬、3.457/10萬和3.469/10萬。如圖1所示，布病病例數(shù)與發(fā)病率曲線完全吻合。

圖1 2006-2017年全國布病發(fā)病曲線Fig.1 Brucellosis incidence curve in China from 2006 to 2017

2.2各省市布病累積發(fā)病率 北方省份和南方省份的界定參照文獻[9]。16個北方省份和15個南方省份中，第1時間段江西省、海南省、貴州省均沒有報告病例，之后的3個時間段31省市均有報告病例。高發(fā)省份(內(nèi)蒙古自治區(qū)、山西省、黑龍江省、吉林省、河北省)布病的累積發(fā)病率在第4個時間段有所下降，但寧夏回族自治區(qū)和新疆維吾爾自治區(qū)的累積發(fā)病率出現(xiàn)明顯提高，如圖2所示。

圖2 各省市布病在4個時間段的累積發(fā)病率Fig.2 Cumulative incidence rate of brucellosis in four time periods in various provinces

2.3布病累積發(fā)病率的空間分布 4個時間段的箱式地圖顯示全國布病累積發(fā)病率的空間分布。第1時間段的累積發(fā)病率在0至113.41/10萬之間，離群省為內(nèi)蒙古、黑龍江、吉林、河北、山西和陜西。第2時間段的累積發(fā)病率在0.002 3/10萬至205.69/10萬之間，離群省為內(nèi)蒙古、黑龍江、吉林、山西、河北和寧夏。第3時間段的累積發(fā)病率在0.029 4/10萬至124.76/10萬之間，離群省為內(nèi)蒙古、黑龍江、山西、寧夏和新疆。第4時間段的累積發(fā)病率在0.08/10萬至100.94/10萬之間，離群省為內(nèi)蒙古、黑龍江、山西、寧夏和新疆。

2.4全局空間自相關(guān)分析 以4個時間段31個省市布病累積發(fā)病率為單位進行全局Moran’s Ⅰ空間自相關(guān)分析，Moran’s Ⅰ指數(shù)分別為0.212、0.190、0.377和0.264(Z值分別為3.800、4.368、4.108和2.861，均有P<0.05)，表明4個時間段31個省市布病累積發(fā)病率存在空間正相關(guān)，即相鄰省市之間有相似的布病發(fā)病風險，隨著時間的推移，發(fā)病風險有上升趨勢。

2.5局部空間自相關(guān)分析 4個時間段的LISA聚集圖主要表現(xiàn)為“高高”“低低”聚集，且只有“高高”“低低”“低高”3種局域相關(guān)類型，無“高低”聚集區(qū)。前3個時間段“高高”聚集發(fā)生在中國的北方區(qū)域主要包括內(nèi)蒙古、黑龍江、吉林、河北、山西等5省市；第4個時間段“高高”聚集發(fā)生在中國的北方區(qū)域內(nèi)蒙古、寧夏和甘肅，并且隨著時間的推移寧夏從“低高”聚集轉(zhuǎn)變?yōu)椤案吒摺本奂?，提示寧夏布病發(fā)病率有顯著的升高。甘肅先從“非顯著性區(qū)域”轉(zhuǎn)變?yōu)椤暗透摺本奂?，最后轉(zhuǎn)變?yōu)椤案吒摺本奂８拭C聚集類型的轉(zhuǎn)變除了與甘肅本身布病發(fā)病率提高有關(guān)外，還與寧夏和新疆發(fā)病率的提高有關(guān)?！暗透摺本奂霈F(xiàn)在“高高”聚集周圍，陜西省在4個時間段均為“低高”聚集?！暗偷汀本奂l(fā)生在中國的南方區(qū)域主要包括云南、貴州、廣西、廣東、湖南、江西和浙江。這種空間聚集現(xiàn)象顯示南北方布病發(fā)病率的差異依然存在。

2.6時間分布 隨著時間的推移，M值和R值均有減小的趨勢，逐漸低于0.3，見表1，提示布病從具有一定季節(jié)性逐漸轉(zhuǎn)變?yōu)榧竟?jié)性不明顯。

表1 全國布病在4個時間段的集中度和圓形分布法分析Tab.1 Analysis of concentration and circular distribution of brucellosis in four time periods in China

2.7年齡分布 因未獲得2017年分年齡組的病例數(shù)，所以第4個時間段只統(tǒng)計2015和2016年的分年齡組的病例數(shù)。4個時間段內(nèi)，均以35～54歲組的發(fā)病人數(shù)居多，老少組的構(gòu)成比分別是7.3%、7.4%、9.9%和12.7%。趨勢性χ2檢驗的結(jié)果表明，隨時間變化，老少組和非老少組構(gòu)成比的趨勢變化差異有統(tǒng)計學意義(χ2=1 946.95，P<0.001)。

3 討 論

本研究以每3年作為一個時間段研究全國布病在2006-2017年流行特征的變化情況，既反映了布病的長期變化趨勢，也避免了布病發(fā)病的逐年波動性，所以具有較強的實際意義。結(jié)果表明2006-2017年全國布病疫情總體呈上升趨勢，但自2015年來布病的年發(fā)病率出現(xiàn)下降趨勢，尤其是布病的高發(fā)省份，如內(nèi)蒙古、黑龍江、山西?！秶抑虚L期動物疫病防治規(guī)劃(2012-2020年)》[10]和《國家布魯氏菌病防治計劃(2016-2020年)》[11]的頒布和實施對遏制布病疫情上升態(tài)勢起到了積極的作用。而南方省份布病累積發(fā)病率明顯低于北方省份，但處于持續(xù)的上升態(tài)勢。這可能與南方省份作為新發(fā)疫區(qū)，缺乏經(jīng)驗有一定的關(guān)系。

空間分析方法可用來分析肺結(jié)核、病毒性肝炎、手足口病等多種傳染病的流行病學特征[12-13]。不同研究者在不同省份根據(jù)該方法均發(fā)現(xiàn)布病分布存在空間自相關(guān)[5-6, 14]，本研究顯示，2006-2017年全國布病在4個時間段的Moran’s Ⅰ指數(shù)為0.190～0.377，說明存在明顯聚集性，主要集中在內(nèi)蒙古附近的北方省份，與2006-2012年的研究結(jié)果一致[15]。LISA集聚圖顯示除新疆維吾爾自治區(qū)外，高高聚集區(qū)與高發(fā)病率區(qū)基本一致；另外寧夏回族自治區(qū)和甘肅省隨著時間的變化逐漸演變?yōu)楦吒呔奂瘏^(qū)，陜西省在4個時間段均為低-高聚集區(qū)。低-高聚集區(qū)要采取積極的措施控制牲畜及其產(chǎn)品從高發(fā)區(qū)域流入，從而預(yù)防疫情蔓延。我國南方省份雖屬于低-低聚集區(qū)，但病例數(shù)卻逐年遞增，建議這些區(qū)域也不能放松警惕，加強與北方省份的交流，積極探索布病的防控措施。

集中度法的M值和圓形分布的R值在4個時間段呈現(xiàn)遞減趨勢，說明布病的季節(jié)性越來越弱。本研究中，2012-2014年的R值為0.273 5，與2012-2016年的R值接近[16]。2004-2010年濟南市的M值和R值為0.44左右[7]，2004-2012年山東省的M值為0.375 3[17]，2004-2015年北京市M值和R值為0.3左右[18]，提示不同省份和不同時間段的M值和R值存在差異。

4個時間段內(nèi)，均以35～54歲組的發(fā)病人數(shù)居多，但隨著時間的推移，老少組的構(gòu)成比升高。2011-2016年福建省的病例中有11.91%是散居兒童、學生、教師、餐飲食品業(yè)等，14起聚集性疫情[6]。2013-2018年甘肅省學生占比為2.42%[14]。這些證據(jù)提示布病從職業(yè)人群向一般人群蔓延，可能與感染方式的多樣化有關(guān)，所以在布病防控過程中既要關(guān)注職業(yè)人群也要關(guān)注一般人群。

綜上所述，全國布病的南北差異依然存在，季節(jié)性逐漸趨于不明顯，老少人數(shù)有增多趨勢。布病防控形勢依然比較嚴峻，既要關(guān)注職業(yè)人群也要關(guān)注一般人群，高發(fā)省份依然要認真落實各項防控措施，低發(fā)省份加強監(jiān)測預(yù)防。

利益沖突：無