1、研究背景:
　　布魯氏桿菌病(Brucellosis，以下簡稱布病)，亦被稱為“波狀熱”、“地中海熱”和“馬爾他熱”，是由布魯氏桿菌引起的高傳染性人畜共患病?；疾〉难颉⑴?、豬、犬等是人間布病的主要傳染源，通過直接接觸、消化道、呼吸道等途徑傳播給人類。該病是一種全身感染性疾病，可感染所有人，嚴(yán)重危害人類健康，并阻礙畜牧業(yè)與社會經(jīng)濟(jì)發(fā)展。
　　布病主要流行于亞洲，撒哈拉以南非洲地區(qū)，中東，地中海，拉丁美洲的一些國家和東南歐地區(qū)。

2、目前，布病仍然是一個嚴(yán)重的公共衛(wèi)生問題，并導(dǎo)致嚴(yán)重的經(jīng)濟(jì)損失和較高的全球健康負(fù)擔(dān)，尤其是在發(fā)展中國家。在我國，山東省報告病例數(shù)位于全國排名前10位。2005年以來，該病發(fā)病率快速增長且波及范圍不斷擴(kuò)大。迄今，山東所有17地（市）均有病例報道。
　　目前我國布病流行態(tài)勢復(fù)雜、疫區(qū)不斷蔓延，給人們的身體健康造成嚴(yán)重危害并導(dǎo)致嚴(yán)重的疾病負(fù)擔(dān)。山東作為我國北方地區(qū)典型疫區(qū)，其流行強(qiáng)度不斷增強(qiáng)，流行范圍不斷擴(kuò)大。本研究以山東省布病流行數(shù)據(jù)，

3、分析該病的流行特征，探討影響該病流行的因素，從而為該病的預(yù)防控制提供科學(xué)依據(jù)。
　　研究目的:
　　1.綜合應(yīng)用傳統(tǒng)流行病學(xué)、空間流行病學(xué)及時間序列分析方法描述該病的人間，時間及空間的流行分布;
　　2.通過1∶2匹配的病例對照研究探究布病的發(fā)病影響因素，為該病的預(yù)防控制提供科學(xué)依據(jù)。
　　資料來源與方法
　　1.資料來源
　　(1)描述性研究數(shù)據(jù)來源于山東省疾病報告信息系統(tǒng)(Shandong Dis

4、easesReporting Information System，SDRIS)2004至2014年布病監(jiān)測數(shù)據(jù)報告卡中監(jiān)測的病例資料。
　　(2)病例對照研究依據(jù)2013年山東省疾病報告信息系統(tǒng)（ShandongDiseases Reporting Information System，SDRIS）資料，選定2013年1月-12月間在7個研究地區(qū)利津、慶云、章丘、鄆城、無棣、青州、臨淄7個研究地區(qū)確診的布氏桿菌病患者，開展1∶2

5、匹配的病例對照研究。
　　2.研究方法:
　　(1)描述性研究采用SPSS16.0軟件，對布病監(jiān)測病例進(jìn)行包括年齡、性別、職業(yè)等人口特征描述，全省各縣（區(qū)）分布情況及時間季節(jié)分布特征。方法有x2檢驗、線性趨勢x2檢驗、Kruskal-WallisH秩和檢驗，以P＜0.05為差異具有統(tǒng)計學(xué)意義?；诮⒌牡乩硇畔⑾到y(tǒng)數(shù)據(jù)庫、地理信息系統(tǒng)軟件(ArcGIS9.3)及空間聚集分析軟件(SatScan7.0)，運用空間專題制圖、空間

6、聚集性分析及時間聚集性分析，探討其時空分布特征。基于SAS9.0進(jìn)行時間序列分析并建立ARIMA模型。
　　(2)病例對照研究EpiData3.0依據(jù)調(diào)查問卷建立數(shù)據(jù)庫。采用SPSS16.0進(jìn)行數(shù)據(jù)分析。病例組與對照組的可比性分析中，連續(xù)變量的比較使用t檢驗，方差不齊的使用t'檢驗;分類變量的比較使用卡方檢驗;有序變量的比較使用Wilcoxon秩和檢驗。以各暴露因素為自變量，以是否為病例為因變量，采用logistic回歸模型進(jìn)行單

7、因素分析，顯著性檢驗水平P＜0.10，結(jié)合布病的生物學(xué)意義，將符合的變量，利用逐步回歸法進(jìn)行多因素分析，計算各變量比值比(OR)及人群歸因危險度百分比（PAR％），顯著性檢驗水平P＜0.05。
　　研究結(jié)果:
　　1.布病流行特征
　　(1)布病“三間”分布:2005-2014年發(fā)病率呈逐年增長趨勢，年均增幅為35.56％。布病病例各月份均有發(fā)生，3-8月為發(fā)病高峰期，占67.41％(4160/6171)，發(fā)病峰值在5

8、月，占13.90％(858/6171)，春夏發(fā)病病例數(shù)所占比例逐年減少趨勢。男性病例占76.17％(4542/6171)，男女比為2.79∶1。發(fā)病年齡為最小為5個月，最大89歲，其中46-55歲最高占24.6％(1518/6171)。病例職業(yè)以農(nóng)民為主。全省17地市均有病例報道，地市前五位分別為菏澤、濰坊、濱州、濟(jì)寧、淄博，分別占本省病例的16.53％(1003/6067)、12.44％(755/6067)、11.44％(694/60

9、67)、10.25％(622/6067)、9.35％(567/6067)。2005-2014年病例波及縣（區(qū)）共計139個縣（區(qū)），分別為20、18、24、23、33、43、68、94、117、133個，呈逐年遞增趨勢（線性趨勢x2=421.57，P＜0.0001）。病例主要報自綜合性醫(yī)院，占41.94％(2588/6171)，病例主要經(jīng)臨床診斷10.71％(661/6171)及實驗室診斷88.64％(5470/6171)。
　　

10、(2)空間流行病學(xué)分析:2004-2012年山東省布病發(fā)病率由0.0382/10萬(35例)上升至0.6205/10萬（598例），年均發(fā)病率為0.2111/10萬。集中度分布M=0.3753，3-6月高發(fā)(56.27％，1041/1802)。布病發(fā)病率的全局指數(shù)值為0.198901（P=0.000120），表明布病發(fā)病率在全省范圍內(nèi)存在空間正相關(guān)性。2006、2007、2009-2012年山東省布病發(fā)病率為空間正相關(guān)。局部空間自相關(guān)分

11、析(LISA)得出8個縣（區(qū)）為“高高”(HH)聚集，經(jīng)LISA可視化顯示，山東省西南部和北部為高發(fā)區(qū)，與發(fā)病率較高地區(qū)空間分布一致。一類聚集區(qū)是以鄄城為中心點、輻射半徑為33.83km，RR=9.72;二類聚集區(qū)是以濱城為中心點、輻射半徑為62.78km，RR=4.99，8處HH區(qū)均在兩類聚集區(qū)域內(nèi)。
　　(3)時間序列分析:ARIMA(0，2，1)模型為最優(yōu)的模型(x2=5.58 P=0.35)。2014年的預(yù)測總發(fā)病率為1.

12、571/10萬人，而且每月發(fā)病率分別為0.101，0.118，0.143,0.166,0.160,0.172,0.169,0.133,0.122,0.105,0.103和0.079/10萬人，實際觀測值均在預(yù)測值的95％CI內(nèi)。MAPE值為58.79％。
　　2.布病流行因素
　　經(jīng)單因素及多因素分析與動物相關(guān)的職業(yè)(OR=1.846，95％CI=1.188-2.870)，家人患過布病(OR=2.572，95％CI=1.05

13、4-6.273)，與羊接觸(OR=4.871，95％CI=3.059-.7.756)，處理病畜難產(chǎn)、流產(chǎn)、正常產(chǎn)(OR=2.258，95％CI=1.257-4.054)，接觸牲畜時有外傷(OR=1.912，95％CI=1.043-3.504)，吃生拌肉或吃未熟的肉或者手不潔就拿食物吃(OR=6.402，95％CI=2.029-20.198)是布病流行的危險因素，與動物密切接觸時使用防護(hù)措施(OR=0.473，95％CI=0.293-0.

14、764)為其保護(hù)性因素。且與動物密切接觸時無防護(hù)措施與危險因素與動物相關(guān)的職業(yè)，家人患過布病，與羊接觸，處理病畜難產(chǎn)、流產(chǎn)、正常產(chǎn)，接觸牲畜時有外傷，吃生拌肉或吃未熟的肉或者手不潔就拿食物吃均有相乘交互作用。
　　結(jié)論:
　　(1)山東省布病發(fā)病率呈逐年增長趨勢疫情波及范圍逐年擴(kuò)大;布病發(fā)病主要集中在46-55歲男性農(nóng)民;布病全年均有病例在3-8月為流行高峰，即春夏季節(jié)。該病存在正空間相關(guān)性，包括8處‘H-H’空間聚集區(qū)，發(fā)

15、病主要位于山東西南部魯西平原和山東省北部的魯北平原，并存在存在兩類聚集區(qū)，發(fā)病趨勢由北部向中部、由西向東范圍逐漸擴(kuò)大，東部沿海威海煙臺出現(xiàn)病例，且范圍擴(kuò)大。布病可應(yīng)用ARIMA(0，2，1)較好的模擬并預(yù)測，并通過模型得出2014年的預(yù)測布病總發(fā)病率為1.571/10萬人，可見布病流行強(qiáng)度持續(xù)增強(qiáng)。
　　(2)與動物相關(guān)的職業(yè)，家人患過布病，與羊接觸，處理病畜難產(chǎn)、流產(chǎn)、正常產(chǎn)，接觸牲畜時有外傷，吃生拌肉或吃未熟的肉或者手不潔就拿