2016-2020年陜西省榆林市人間布魯氏菌病流行特征與時空聚集性

羅波艷，聶守民，范鎖平，郭少華，安翠紅，王文靜，孫養信

榆林市位于陜西省最北部，東臨黃河與山西相望，西連寧夏、甘肅，北鄰內蒙古，轄1市2區9縣，20世紀90年代以后羊產業成為了榆林市農民經濟收入的重要來源，據統計2019年榆林市羊只存欄量為5 491 233只，占全省羊只總量的60%左右[1]。布魯氏菌病(簡稱布病)是由布魯氏菌引起的一種人獸共患病，易感動物以羊、牛、豬為主。人患病后臨床癥狀不典型，容易轉為慢性，造成不同程度的勞動能力的喪失[2]。從1951年以來，榆林市一直是陜西省布病流行的重點區域，布病發病率水平也一直位于全省地市首位，且遠高于全省布病平均發病率水平[3]，病例感染來源以病羊為主。近年來，在陜西省政府加快肉牛肉羊產業發展意見和榆林市百億元級羊產業政策的推動下，羊產業發展加速，羊只存欄量迅速增加，人間布病發病水平也一直居高不下，嚴重影響了榆林市畜牧業發展和人民群眾生命健康，但關于榆林市布病防控相關研究甚少。本研究通過對2016-2020年榆林市布病發病情況進行時空聚集性分析研究，探討布病在榆林市的時空分布特征，確定相關的高風險地區，為制訂有效的防控策略提供科學依據。

1 材料與方法

1.1 資料來源 2016-2020年陜西省榆林市布病發病數據和人口數據分別來源于《中國疾病預防控制信息系統》和《中國疾病預防控制信息系統》的《基本信息系統》，榆林市縣區電子地圖矢量數據來源于全國地理信息資源目錄服務系統(https://www.webmap.cn)。布病診斷依據《布魯氏菌病診斷標準》(WS269-2019)。

1.2 方 法

1.2.1 描述性分析 用Microsoft Excel 2010建立2016-2020年榆林市布病發病信息數據庫，對數據進行整理分析，運用描述流行病學方法對榆林市布病流行特征進行分析，年均發病率(/10萬)=發病總例數/總人口數×100 000。

1.2.2 三維趨勢分析 將2016-2020年榆林市各縣(區、直轄市)布病發病率視為三維空間中的一個點(X，Y，Z)，其中X和Y表示各縣(區、直轄市)的經緯度坐標，Z表示年發病率水平，運用ArcGIS 10.8軟件對布病發病水平的三維趨勢進行分析。

1.2.3 空間自相關分析 空間自相關分析可以反映疾病在空間上的分布信息及疾病發生與周圍疾病發生的關聯強度，參照寧文艷[4]的方法，采用GeoDa 1.6軟件對2016-2020年榆林市布病進行空間自相關分析。1)全局空間自相關分析，計算全局Moran'sI值，Moran'sI值取值在(-1～1)，計算值越接近1，說明空間正相關性越強；計算值越接近-1，說明空間負相關性越強。當Z值的絕對值大于1.96，P<0.05，認為差異有統計學意義，即疾病存在空間自相關性。2)局部自相關分析：用于進一步探索聚集類型和聚集區域位置，局部空間自相關分析聚集圖共有“高-高(H-H)”聚集、“低-低(L-L)”聚集、“低-高(L-H)”聚集和“高-低(H-L)”聚集。

1.2.4 時空掃描分析 運用SaTScan 9.4.1軟件對2016-2020年榆林市布病發病數據進行分析。建立Possion模型，以縣(區)為空間單位，年為時間單位進行掃描，聚類范圍的最大空間掃描區域設定為30%的風險人口[5]，最大掃描時間設定為30%的研究期限，用對數似然比(Log likehood ratio，LLR，為掃描計算窗口內與窗口外區域之間的死亡率差異)來評估掃描窗口內布病發病的異常程度，概率P值由蒙特卡羅隨機化模擬得到，當P<0.05時，差異有統計學意義，認為該區域存在聚集性發病，且相對危險度(relativerisk，RR)值越大，說明聚集程度越高[6-7]。

2 結 果

2.1 流行基本情況 2016-2020年陜西省榆林市共報告布病病例1 773例，年均發病率為10.40/10萬，無死亡病例，其中2016年報告發病率最高(12.11/10萬)，2018年報告發病率最低(8.47/10萬)(見表1)，且各年報告發病率水平均高于陜西省全省報告發病率水平。

2.1.1 時間分布 2016-2020年榆林市1-12月均有布病病例報告，7月份報告病例數最多(243例，占13.71%)，總體以5-8月為主，占總報告病例數的49.41%(876/1 773)，詳見表2。

2.1.2 地區分布 2016-2020年榆林市10個縣區和省直轄神木市和府谷縣均有布病病例報告，累計報告病例數排在前3位的縣區有綏德縣(303例)、定邊縣(232例)和靖邊縣(206例)；年均發病率水平≥10.00/10萬的縣區有8個，占榆林市縣(區、市)總數的75%(8/12),平均發病率排在前3位的縣區為吳堡縣(29.17/10萬)、子洲縣(21.05/10萬)和清澗縣(21.03/10萬)，詳見表1。

表1 2016-2020年陜西省榆林市布病發病數(例)和發病率(/10萬)情況

2.1.3 人群分布 在2016-2019年報告的布病病例中，男女性別比為3.53∶1(1 382/391)，平均發病率分別為15.34/10萬、4.87/10萬；0～14歲、15～44歲、45～64歲、≥65歲的布病病例分別占0.90%(16/1 773)、27.75%(492/1 773)、60.46%(1 072/1 773)、10.89%(193/1 773)。見圖1。職業以農民為主，占89.57%(1 588/1 773)。

圖1 2016-2020年陜西省榆林市布病病例的年齡構成情況

2.2 時空分布特征

2.2.1 三維趨勢分析 圖2中，藍色曲線代表南北方向，綠色曲線代表東西方向。通過三維趨勢分析，我們發現2016-2020年榆林市布病發病率呈現由北向南、由西向東逐漸升高。具體來看，除了2016-2020年榆林市布病發病率在東西方向上表現為東部高于西部、中部最低，其中2018年發病率呈U型的發病特點，其余年份呈S型的發病特點；在南北方向上，榆林市布病發病率呈現出南部高于北部，其中2016-2018年、2019年、2020年榆林市布病發病率分別呈U形、倒U形、S形的發病特點。

圖2 2016-2020年陜西省榆林市布病發病率三維趨勢分析

2.2.2 全局自相關分析 2016-2020年榆林市12個縣(區、市)人間布病報告發病率的全局Moran's I統計量結果顯示Moran'sI系數均為正值，但除2020年外，其余年份Moran'sI系數均有統計學意義(P<0.05)，表明2016-2019年榆林市各縣(區、市)人間布病報告發病率呈現明顯的空間正相關性(詳見表2)。

表2 2016-2020年陜西省榆林市布病病例月報告情況

2.2.3 局部自相關分析 2016-2020年對榆林市人間布病報告發病率的局部Moran'sI系數進行假設檢驗有統計學意義上的區域有：1)高—高聚集區共有2個，2016年為清澗縣，2017、2018、2020年均為綏德縣，2019年為清澗縣和綏德縣；2)低-低聚集區共有3個，2016年為榆陽區和神木市，2017、2019年均為神木市，2018年為神木市、榆陽區和橫山區，見圖3。

A: 局部Moran's I顯著水平圖；B: 局部Moran's I聚類圖。

表3 2016-2020年陜西省榆林市人間布病報告發病率全局自相關分析結果

2.2.4 時空掃描分析 2016-2020年榆林市布病發病共掃描到1個一類時空聚集區和1個二類時空聚集區。一類時空聚集區以米脂縣為中心，半徑為51.62 km的圓形區域，覆蓋于佳縣、綏德縣、吳堡縣、子洲縣，聚集時間為2019年，該地區在高危時期共報告了198例布病病例。二類時空聚集區為靖邊縣為中心，因半徑為0 km，所以未顯示，見表4和圖4。

表4 2016-2020年陜西省榆林市布病發病時空掃描結果

圖4 2016-2020年陜西省榆林市布病時空聚集區

3 討 論

榆林市是陜西省第一個報告布病病例的地市，1951-1981年榆林市布病發病率達到了100/10萬以上，在人、畜間堅持開展了菌苗大面積連續性預防接種工作之后，發病率有所降低，但1996年榆林市綏德縣再次發生布病暴發流行，造成了布病在榆林地區的再度活躍[8]。經過人、畜間共同努力，近年來榆林市布病發病率水平總體維持在10/10萬以上，在全省范圍內仍屬于布病高流行區域。

本研究通過對2016-2020年榆林市布病病例月報告情況分析，發現每年基本上從3月病例數開始上升，7月達到最高，9月開始下降，但2019年從5月開始上升，7月達到最高，之后高峰一直維持到12月份，基本上呈全年流行狀態。這可能與近年來榆林市不斷加強封山禁牧監管力度有關，封山禁牧導致牲畜飼養模式由放養轉為圈養，加大了高危人群與牲畜接觸的頻率，也增加了經呼吸道感染布病的可能，導致月報告病例數居高不下。因此，榆林市應該在布病流行期來臨之前，開展高危人群摸底，加強布病防控知識宣傳，評估布病發病風險，根據風險評估結果采取有效的防控措施，切實達到防病關口前移的目的。

趨勢分析可以排除不穩定因素的干擾，以三維圖的形式直觀展示出發病率在大尺度上的變化趨勢[9]。本研究三維趨勢分析結果顯示，2016-2020年榆林市布病發病率呈現北低南高、東部略高于西部的特點，這可能與榆林市北部的神木市和府谷縣大力發展能源產業，羊只養殖業發展滯后，羊只存欄量少有關?？臻g自相關分析結果顯示榆林市布病存在明顯的聚集性，且呈現出 “高-高”和“低-低”聚集，聚集類型與其發病率高低分布基本一致。有研究指出，布病是一種自然疫源性疾病，受相鄰地區的自然和社會環境影響較大，導致發病地區多表現為“高-高”或“低-低”聚集[10]，所以說榆林市東、西部發病率高可能是因為其東部與山西省相鄰，西部與寧夏、甘肅省相鄰，這3個省份均為我國布病高發區域[11-12]，此外，榆林市榆陽區、神木市和橫山區雖屬于“低-低”聚集區，但從2017年開始，報告布病病例數卻呈逐年遞增的趨勢，建議榆林市重視這些區域的布病防控工作，加大布病防控知識宣傳力度。

空間自相關分析對疾病的分析局限在空間維度上，而時空掃描分析從時間和空間2個維度開展疾病分析，可以進一步提高分析結果的準確性和實用性[13]。本研究結果顯示，2016-2020年榆林市發病聚集時間為2019年，發現的聚集區與空間自相關分析掃描結果基本一致，聚集區域還包括了發病率高的子洲縣和吳堡縣。據人間摸底調查，靖邊縣羊只存欄量已達到130萬只左右，交通便利，羊只引進頻繁，畜間檢疫難度大；綏德縣、佳縣、子洲縣羊只存欄量均在20萬左右，養殖量在全市屬于中等水平，但其形成的規?；B殖集散市場導致布病傳播風險升高[14]；吳堡縣、米脂縣羊只存欄量不多，以散養戶為主，但山大溝深、交通不便，養殖方式落后，牲畜檢疫難度大，這些都給當地布病防控帶來很大難度。

近年來，隨著榆林市湖羊養殖項目的大力推進，養殖場規模和數量擴大，羊只存欄量迅速擴增，散養戶數量也只增不減，同時《榆林市封山禁牧管理辦法》規定的出臺使得家庭散養戶對羊只的飼養模式由原來的放養轉為圈養，很大程度上增加了人畜接觸機會和經呼吸道感染的可能性，這些共同導致榆林市布病防控形勢異常嚴峻。因此，建議榆林市應結合布病發病高峰時間段、發病人群、時空聚集區、羊只存欄量等情況，開展轄區布病風險研判，采取分區域、有重點、人畜同步聯防的布病防控策略，堅持“消滅傳染源、切斷傳播途徑、保護易感染人群”的傳染病防控方針，加強人畜間布病聯防聯控工作，同時加強養殖戶布病防控知識和個人防護行為重要性的宣傳，遏制布病發病水平。

局限性：本研究未得到2016-2020年榆林市各縣(區、市)的羊只存欄數據和陽性病畜數據，這造成無法對各縣(區、市)人間布病發病情況進行進一步探討分析。

利益沖突：無

引用本文格式：羅波艷，聶守民，范鎖平，等. 2016-2020年陜西省榆林市人間布魯氏菌病流行特征與時空聚集性[J]. 中國人獸共患病學報，2022，38(2)：122-127. DOI:10.3969/j.issn.1002-2694.2022.00.025