研究背景恙虫病(tsutsugamushi disease)又称丛林斑疹伤寒(scrub typhus),是由恙虫病东方体(Orientia tsutsugamushi)所致的急性自然疫源性传染病,临床上以特征性焦痂或溃疡、高热、淋巴结肿大、皮疹为特征,并可引起多种并发症甚至死亡恙虫病广泛流行于亚洲太平洋沿岸三角地区,北至日本北部,南至澳大利亚北部,西至巴基斯坦、阿富汗。我国属于重点疫区之一。1985年以前,恙虫病仅在北纬31。以南广泛流行,多于夏季高发,被称为夏季型恙虫病。1985年以后,我国长江以北该病新疫区不断被发现和证实,流行形式主要呈秋冬型。山东省为秋冬型恙虫病典型疫区,自1986年在蒙阴县发生该病流行后,其流行范围迅速扩大,至今,已波及全省约80%的广大地区。秋冬型恙虫病,作为北方的一种新发传染病,流行范围不断扩大,流行强度不断增强,给公众身体健康造成了严重危害。本研究以该病典型疫区——山东省为代表,分析了秋冬型恙虫病的流行特征,探讨了影响该病流行的因素,对该病防控对策与措施的制定具有重要意义。研究目的1.综合应用传统流行病学方法、空间流行病学方法及时间序列分析方法描述该病的流行特征；2.通过1：2匹配的病例对照研究分析探讨该病的影响因素,从而为该病防控对策与措施的制定提供科学依据。资料来源与方法1.资料来源(1)描述性研究研究数据来自山东省疾病报告信息系统(Shandong Diseases Reporting Information System, SDRIS)2006年至2011年秋冬型恙虫病监测资料,共计1746例。(2)病例对照研究基于山东省疾病报告信息系统(Shandong Diseases Reporting Information System, SDRIS)2010年秋冬型恙虫病监测资料,于2011年1-3月,在泰安新泰市、临沂沂南县等9县(市或区)对2010年恙虫病患者128例,及按1：2匹配的对照256例予以调查、收集资料。2.研究方法(1)描述性研究描述恙虫病患者的人口分布特征,包括性别、职业、年龄构成比,各地区分布以及时间分布特征,检验方法为x2检验、Kruskal-WallisH秩和检验,以P<0.05为差异有统计学意义。基于建立的地理信息系统数据库、地理信息系统软件(ArcGIS9.3)及空间聚集分析软件(SatScan7.0),运用空间专题制图、距离倒数插值、空间自相关分析、空间聚集分析、时间聚集分析及时间趋势下空间变异分析,探讨该病的时空分布特征。采用频谱分析、移动平均分析、季节效应分析以及建立ARIMA模型对该病进行时间序列分析。(2)病例对照研究采用问卷对研究对象进行调查；基于SPSS16.0软件,对资料进行整理、变量进行量化后,采用单因素、多因素条件Logistic回归进行分析,以OR及95%CI、人群归因危险度百分比(PAR%)表示疾病与暴露之间的关联强度。显著性检验水平α=0.05。研究结果1.恙虫病流行特征(1)恙虫病“三间”分布：1746例恙虫病患者,发病年龄分布为1-92岁,主要集中在55岁以上年龄组,占51.1%(892/1746)；男性占48.3%(843/1746),女性占51.7%(903/1746),男女比为1：1.07；农民、工人、学生及学龄前儿童均有发病,以农民为主,占84.6%(1477/1746)。该病主要流行于10-11月,以10月为高峰,占发病总数的71.0%(1240/1746)。全省17地市,以临沂、泰安、日照、莱芜、青岛发病数为多,分别占发病总数的22.9%(399/1746)、18.1%(316/1746)、14.7%(256／1746)、14.1%(246/1746)以及13.1%(229/1746)。(2)空间流行病学分析：全局Moran’sI指数值为0.324(P<0.01)。8处地点该病发病率的局部Moran’sI指数值有统计学意义(P<0.01),均为“高一高”(H-H)聚集,其中较为重要的地区为钢城、莱城、新泰、东港,发病数分别为38、154、160、105例,局部Moran’sI值分别为2.111、1.642、1.277、0.775。各年发病聚集时间分别为：2006年9月23日至11月20日(202例)、2007年10月2日至11月11日(197例)、2008年9月30日至11月7日(302例)、2009年9月25日至11月10日(204例)、2010年10月5日至11月13日(226例),RR值分别为45.55、34.60、50.64、53.09、79.84(P值均<0.01)。共发现以泰安新泰为中心的辐射半径为58.28km,及以日照东港为中心的辐射半径22.68km的两处恙虫病聚集区,病例数分别为542例、134例,RR分别为4.52、3.96(P值均<0.01),地貌分别为内陆低山丘陵和沿海丘陵；以临沂蒙阴县为中心的辐射半径为45.82km的区域,恙虫病年增长率最高,达45.04%,RR=3.68(P<0.01)。(3)时间序列分析：频谱分析结果显示秋冬型恙虫病流行具有3年的周期性；移动平均分析表明该病流行强度将持续增强,2012—2014年各年预测发病数分别为310、337、366例,年增长率约9%；乘法季节效应分析显示本病符合秋冬型流行特点,10、11月季节效应指数分别为8.454、2.230,其他月份该指数均小于1。构建本病发病ARIMA(0,1,1)(0,1,0)12模型为(1-B)(1-B12)X1=(1-0.811B)ut。模型残差序列为白噪声,Box-Ljung统计量为3.116,P=-0.999,模型拟合效果良好。应用该模型预测2011年11、12月该病发病数,预测值与实际发病相符,实际观测值均在预测值的95%CI内。2.恙虫病流行因素经采用单因素及多因素条件Logistic回归分析,结果显示住宅周围有杂草滋生或作物种植、居住村边、发病前感知被不明昆虫叮咬、院内堆放柴草或农作物秸秆、居住环境潮湿、个人卫生习惯欠佳、居住草房为该病的危险因素,OR及95%CI分别为2.815[1.493,5.307]、2.250[1.208,4.190]、4.500[2.206,9.180]、2.100[1.163,3.793]、5.788[1.819,18.418]、1.943[1.000,3.810]、4.049[1.062,15.434],PAR%分别为44.04%、41.45%、29.77%、27.44%、19.56%、13.95%、7.68%。鼠类传播疾病知识认知为该病的保护性因素,OR及95%CI为0.523[0.292,0.936],PART%为47.57%。结论1.秋冬型恙虫病以中老年农民为主,9月下旬到11月上旬为流行高峰；该病存在正空间相关性,均为“高-高”(H-H)聚集,其中较为重要的地区为钢城、莱城、新泰、东港；共存在内陆低山丘陵和沿海丘陵两个聚集区,其中,以临沂蒙阴为中心的区域该病年增长率最高。秋冬型恙虫病流行具有3年的周期性,总体流行强度将持续增强,2012—2014年该病预测发病数分别为310、337、366例,每年递增约9%。应用ARIMA(0,1,1)(0,1,0)12模型可较好地预测该病发病情况。2.秋冬型恙虫病的流行与居住环境状况、个人卫生状况及健康知识认知情况相关。其中,住宅周围有杂草滋生或作物种植、居住村边、发病前感知被不明昆虫叮咬、院内堆放柴草或农作物秸秆、居住环境潮湿、个人卫生习惯欠佳、居住草房为该病危险因素,鼠类传播疾病知识认知为该病的保护性因素。