恙虫病是由恙虫病东方体感染引起的一种自然疫源性疾病，鼠类为主要储存宿主，经恙螨幼虫叮咬而传播。恙虫病广泛流行于亚太地区，全球每年估计发生100万病例^[1]。我国恙虫病原仅在长江以南地区流行，为夏季型。自1986年山东省蒙阴县发现恙虫病暴发以来，该病在北纬31°以北地区迅速蔓延。迄今，除山东省外，天津、北京、山西、河南、江西、江苏等地均已有该病流行的报告，流行特征为秋冬型。近年来，我国北方地区秋冬型恙虫病流行区域不断扩大，流行强度不断增加，但目前与该病流行相关的许多问题尚不清楚。为此本文应用制图、空间自相关分析、空间及时空扫描分析等技术探讨山东省秋冬型恙虫病的时空分布特征和流行趋势，为其预防控制提供参考。

1. 资料来源：恙虫病监测数据源自2006－2014年山东省疾病报告信息系统。人口学资料来自2006－2014年山东省各市统计年鉴。病例诊断标准按中国CDC恙虫病监测标准^[2]。

2. 分析方法：

（1）描述性分析：根据监测数据的发病率资料，运用地理信息系统（GIS）平台绘制发病率地区分布图，对恙虫病的空间分布进行可视化描述。发病率地区分布图用ArcGIS 9.3软件（ESRI Inc.，Redlands，CA，USA）制作。

（2）空间自相关分析：以全局空间自相关分析Moran's Ⅰ指标探测整个研究区域内的空间聚集模式。数据采用GeoDa 0.9.5-i软件处理，显著性检验定为α＝0.05。Moran's Ⅰ的显著性由Monte Carlo随机化检验评估，检验次数设置为999次。Moran's Ⅰ的范围介于-1～1之间。Moran's Ⅰ＞0时，表示疾病存在正向空间自相关；Moran's Ⅰ＜0时，表示疾病存在负向空间自相关。Moran's Ⅰ的值越接近1或-1，表明空间自相关性越强^[3]。

局域空间自相关分析通过局域Moran's Ⅰ与局域型空间自相关（local indicators of spatial autocorrelation，LISA）聚集图实现。LISA聚集图上共有4种局部空间聚集模式：“高-高”（H-H）聚集，“低-低”（L-L）聚集，“低-高”（L-H）聚集，“高-低”（H-L）聚集。其中，“高-高”聚集，“低-低”聚集对应于全局自相关指标中的正向空间自相关，“低-高”聚集，“高-低”聚集对应于全局自相关指标中的负向空间自相关^[4]。利用Z检验对LISA统计量进行假设检验，若P＜0.05，即推断恙虫病病例的空间分布具有自相关性。

（3）空间扫描及时空扫描分析：根据病例发生地的经纬度编码，建立GIS数据库。空间扫描及时空扫描分析以空间动态窗口扫描统计为基础，基于Poisson模型，应用SatScan 9.1.1软件，采用1∶100 000的多边形地图在县级水平上分析恙虫病的时空聚集性^[5]。以圆形或椭圆形窗口表示研究区域的地理范围，以高度表示时间。零假设为在没有协变量存在的情况下，期望病例数与人口数成正比，即扫描窗口内外RR值大小相同。随着扫描窗口的变化，采用对数似然比（LLR）检验计算所有窗口内外区域之间的发病率差异，寻找其中LLR最大的窗口为最有可能存在聚集性的区域。计算该区域的RR值，并检验差异有无统计学意义，推断确定聚集区域^{[6, 7]}。本研究最大空间集群设置为总人口的50%，Monte Carlo随机化检验评估检验次数设置为999次，统计学检验水平设为0.05。

1. 描述性分析：2006－2014年山东省共报告恙虫病4 453例，发病率由2006年的0.23/10万上升至2014年的1.53/10万，发病率呈逐年上升趋势，发病地区逐年增多，高发地区有向东移动并逐渐扩大的趋势（图 1）。全省17个地（市）中有13个发现恙虫病病例，其中枣庄市自2013年开始出现报告病例，并由2013年的10例迅速增加为2014年的124例。

图 1 2006－2014年山东省恙虫病发病率地区分布

图选项

2. 空间自相关分析：2006－2014年山东省恙虫病平均发病率的全局Moran's Ⅰ＝0.501 5，Z＝8.776 4（P＜0.01），表明其在全省范围内存在正向空间自相关性。LISA聚集图显示H-H聚集区主要地区为日照市五莲县、岚山区和莒县，泰安新泰市，莱芜市莱城区和钢城区，淄博市沂源县与临沂市蒙阴县（图 2）。2006－2014年各年发病率的自相关分析结果见表 1 ，全局Moran's Ⅰ在0.05的检验水准下均有统计学意义。

图 2 2006－2014年山东省恙虫病年均发病率LISA聚集图

图选项

表 1 2006－2014年山东省恙虫病年均发病率自相关分析

表选项

3. 空间扫描分析及时空扫描分析：2006－2014年山东省恙虫病病例呈非随机分布，各年份有统计学意义的高危聚集区列于表 2。2006－2014年山东省恙虫病空间扫描聚集显示（图 3），一级聚集区在2006－2010年及2012年主要集中在中南部山区，2011、2013及2014年主要集中在东部丘陵地区，聚集区的范围在2008、2011、2013和2014年呈不断扩大趋势，并于2014年首次在枣庄市发现空间聚集区台儿庄区与市中区。2006年及2013年发现二级聚集区。

表 2 2006－2014年山东省恙虫病空间扫描分析

表选项

图 3 2006－2014年山东省恙虫病空间扫描聚集图

图选项

时空扫描结果发现2014年10月1日至2014年11月30日以日照东港区为中心、辐射半径为222.34 km的聚集区，此聚集区包含64个区（县），RR＝38.72（P＜0.01）。

本文采用空间流行病学分析方法探讨山东省2006－2014年秋冬型恙虫病的时空分布特征，结果发现恙虫病流行地区呈逐年增多。山东省自1986年首次在蒙阴县发现该病以来，流行地区蔓延迅速，17个地（市）中有13个已发现该病流行，且流行强度不断增加，是秋冬型恙虫病的典型流行区域^[8]。近年研究发现，在山东省检测出恙虫病东方体（Ot）的基因型呈多样性特点^[9]，秋冬型恙虫病病原体基因变异活跃是否与该病流行特征有关还有待深入探讨。

本文结果显示，山东省恙虫病的分布呈非随机分布，Moran's Ⅰ值均为正值，说明2006－2014年该病存在正向空间自相关性，即某一空间区域的发病率高（或低），其相邻区域相应的发病率也高（或低），LISA聚集图显示大多数发病聚集区位于山区和丘陵，提示该病在山东省具有空间聚集性。山区和丘陵地区植被等自然环境，有利于鼠类和恙螨的孳生繁殖^{[10, 11]}。从空间扫描聚集图可以看出空间聚集区从中南部山区向东部丘陵地区移动，聚集区面积大约2～3年扩大一次。这种现象可能与降水量等气象因素有关；此外，农业、林业基本建设与发展也可以对鼠类和恙螨的栖息环境产生影响^{[12, 13]}。

空间扫描结果发现枣庄市为秋冬型恙虫病新流行地区，2013年之前枣庄地区没有报告病例，最近两年开始出现并迅速增加，空间扫描聚集图发现有2个空间聚集区在该地区形成。枣庄地区是山东省降水量最丰富的地区之一，日照时间充足，这些条件可为虫媒传染病的发生和传播提供有利条件，从而促使恙虫病新疫源地的形成和扩散。

利益冲突 无