2005—2018年玉林市流行性腮腺炎流行特征和时空分析

(广西玉林市疾病预防控制中心，广西 玉林 537000)

流行性腮腺炎简称流腮，是由腮腺炎病毒感染引起的急性呼吸道传染病，以腮腺非化脓性肿胀、疼痛伴发热为主要症状，可导致多种腺体组织、神经系统和肝脏等器官损伤，可引起脑膜炎、睾丸炎、卵巢炎等并发症，对儿童危险极大[1]。因此，各地及时分析性腮腺炎流行特征，掌握疫情动态，对开展防病宣教，减少发病非常重要。当前对流腮的分析主要集中在三间分布，而三间分布分析不仅容易忽视不同时期和地区传染病的高发具有偶然性，而且还存在人为选择偏倚[2]。因此，有必要在行三间分布分析的同时，结合地理信息系统(geographic information system，GIS)分析传染病的流行特征，以解决三间分布分析的不足。目前GIS在公共卫生领域的应用十分广泛，通过GIS可分析疾病流行的空间特征和分布模式，探索可能的影响因素及病因，为疾病流行的监控、预警、防制效果评价等提供参考依据[3]。因此，本研究在对玉林市流腮流行特征行三间分布分析的同时，利用GIS探索其时空聚集性，客观定量地对疫情进行分析，为玉林市各县(市、区)科学地制定流腮的防控策略提供参考依据。

1 资料与方法

1.1 资料来源 2005—2018年玉林市流腮发病数据来自中国疾病预防控制信息系统中的传染病报告信息管理系统，人口资料来源于基本信息系统；地图来源于2015年1∶100万全国基础地理信息数据。

1.2 分析方法

1.2.1 流腮流行特征分析 采用描述流行病学方法对玉林市2005—2018年流腮病例资料进行三间分布分析，其中年龄分组采用传染病报告信息管理系统的年龄分组形式。利用Microsoft Excel 2010对病例数据进行统计分析。

1.2.2 空间自相关分析 采用ArcGIS 10.2软件对中国矢量地图进行裁剪得到玉林市以县(市、区)为界的矢量地图，通过ArcMap将发病数据进行空间可视化展示。采用GeoDa 1.14软件对2005—2018年玉林市流腮病例资料进行空间自相关分析，分别以全局Moran’s I与局部空间关联指数(local indicators of spatial association，LISA)为指标进行全局空间自相关分析和局部空间自相关分析[4]。全局空间自相关分析从全市角度分析是否存在流腮空间聚集，局部空间自相关分析则进一步分析具体聚集县(市、区)及聚集模式。Moran’s I的取值范围为 [-1，1]，若P<0.05且Moran’s I值>0，表示研究变量存在空间正相关；若P<0.05且Moran’s I值<0，表示研究变量存在空间负相关；Moran’s I值接近于0表示空间随机分布，不存在相关性。采用蒙特卡罗模拟方法对 Moran’s I值进行检验，计算得到Z值和P值，以P<0.05为差异有统计学意义。通过LISA分布图可获得5种局部空间自相关分析模式：不显著、高-高聚集模式、低-低聚集模式、低-高聚集模式和高-低聚集模式。

1.2.3 时空扫描分析 采用SaTScan 9.5软件进行时空扫描统计分析。选用离散变量Poisson分布模型，用圆形移动窗口法创建扫描窗口，扫描时间范围设置为“2005—2018 年”，时间间隔“1年”，地区重叠设置为“无”，最大聚集范围窗口为总人口的30%[5]。选择蒙特卡罗模拟次数“999”，计算在不同圆心和不同半径下的动态窗口区域内外空间单元属性的对数似然比 (log likelihood ratio，LLR)进行统计推断，以探测聚集中心和聚集半径。根据LLR的大小确定聚集区类型，若P<0.05，则认为该区域存在聚集性[2，6]。

2 结果

2.1 2005—2018年玉林市流腮发病的时间分布特征 2005—2018年玉林市共报告流腮10 798例，年均报告发病率为13.80/10万，发病率最高的年份是2012年(23.12/10万，1278例)，最低的是2017年(5.47/10万，315例)；2005—2009年发病率呈下降趋势，2010—2012年发病率回升并达到最高峰，2013—2015年发病率又大幅下降，2016—2018年发病率波动较大，但总体是上升趋势。因此，2005—2018年玉林市流腮发病呈现“降-升-降-升”模式(见表1、图1)。按季节和月份分析，发现每年全年各月均有病例报告，但存在明显的季节性，4～7月为发病主高峰，11月至次年1月为发病次高峰，呈双峰型分布，其中以2010—2013年期间发病高峰最密集。见图2。

表1 2005—2018年玉林市流腮发病情况

图1 2005—2018年玉林市流行性腮腺炎报告病例数和年均发病率

图2 2005—2018年玉林市流行性腮腺炎月发病曲线

2.2 2005—2018年玉林市流腮发病的地区分布特征 玉林市各县(市、区)均有病例报告，玉州区年均发病率最高且报告病例数最多(28.77/10万，2596例)，福绵区报告病例数最少(641例)，但博白县年均发病率最低(5.08/10万)，见表2。

表2 2005—2018年玉林市各县(市、区)流腮发病顺位统计表

2.3 2005—2018年玉林市流腮发病的人群分布特征 按年龄分析，病例主要集中在3～14岁(占81.77%)，其中5～7岁年龄组发病率最高，分别为87.37/10万、96.71/10万、93.34/10万，见表3。按性别分析，发现共报告男性病例7181例(占66.50%)，女性病例3617例(占33.50%)，男女年均发病率分别为17.51/10万、9.72/10万，男性发病率高于女性，各年男女发病率比在(1.43～2.14)∶1之间，见表1。按职业分析，发病以中小学生为主(6088例，占56.38%)，其次是幼托儿童(2151例，占19.92%)，散居儿童(1612例，占14.93%)，农民(591例，占5.47%)，其他职业者(356例，占3.30%)。

表3 2005—2018年玉林市不同年龄组流腮发病情况

2.4 2005—2018年玉林市流腮疫情的GIS分析

2.4.1 空间自相关分析 采用GeoDa软件，生成Rook邻接的空间矩阵，分析2005—2018年玉林市流腮的Global Moran’s I与Local Moran’s I指数，可见除2017年(Moran’s I=0.336，Z=2.255，P=0.020)存在空间正相关外，其余年份均不存在空间自相关性，即无空间聚集性；LISA聚集分布图显示2017年流腮发病率呈高-高聚集模式，聚集区为北流市，见表4。

表4 2005—2018年玉林市流腮发病率空间自相关分析结果

2.4.2 时空扫描分析 时空扫描分析显示，2005—2018年玉林市流腮的空间分布是随机的，存在1个主要聚集区和4个次要聚集区，2016年北流市存在一个主要时空聚集区(LLR=331.22，RR=3.55，P<0.001)；次要聚集区分别为2018年博白县、2010—2013年兴业县和福绵区、2008—2009年容县、2005—2007年玉州区，结果见表5、图3。

表5 2005—2018年玉林市流腮发病率时空扫描统计结果

图3 2005—2018年玉林市流腮时空扫描聚集性可视化

3 讨论

2005—2018年玉林市流腮年均报告发病率(13.80/10万)低于同期全国(21.44/10 万)[7]和广西(37.92/10万)[8]年均发病率水平，其中2012年为玉林市流腮发病高峰年(同期，全国亦出现流腮高发[7，9])。玉林市流腮发病曲线与谢群等[10]报道的流腮每隔2～5年会出现一次较大规模的流行的结论相一致，当易感者积累到一定数量时，新的发病高峰便会再次出现。玉林市流腮发病呈现明显的季节性，与相关报道一致[8，10]。流行曲线与学校开学和放假时间安排相吻合，每个学期开学后就会逐渐出现一次发病高峰。提示在学校和托幼机构出现病例后若未能及时采取有效的隔离和治疗措施，则很容易发生交叉感染而导致暴发。

从地区分布来看，玉林市玉州区年均发病率最高，可能原因是玉州区作为玉林市政治、经济、文化的中心，流动人口多且流动频繁，进城务工人员的子女免疫接种不规范，易感人群聚集，感染和传播传染病风险高，因此更易使城区流腮高发。分析发现，玉林市流腮发病年龄主要集中在3～14岁，尤以5～7岁儿童发病率最高，与全国流腮发病情况相似[7]，男性发病率高于女性，与既往文献结果一致[2，8，10]。也说明流腮的易感人群为幼儿和中小学生[11]，其中又以男生为甚。这可能与男生卫生意识差，天性好动，感染机会大等有关。

通过采用GIS进行空间自相关分析和时空扫描分析，可从时间和空间两个方面同时探讨疾病在地理空间上聚集区域随时间变化的规律性，并获取聚集区域的相对危险度，更准确地定位聚集区域的位置[6]。空间自相关分析结果显示，2005—2016年和2018年玉林市流腮发病不存在空间聚集性，可能是流腮在玉林市是随机发病的，不存在空间自相关；在实际研究中，空间自相关是不均匀的，而Moran’s I指数仅反映空间整体自相关，且个别局部对象的属性取值对全局分析影响显著[12]。时空扫描分析结果显示存在5个聚集区，聚集时间与流腮月发病曲线高峰相同，结合地区分布可进一步认为玉林市流腮的发生是随机分布的，考虑可能的原因是流腮自然流行周期引起随机散发、人口流动密集使易感者增多引起爆发、疫情监测和疫苗接种工作的不断完善使各地发病率得到有效控制，具体原因有待进一步探索。

根据玉林市流腮流行特征分析和时空分析结果。建议：一是对儿童要做好两剂次含腮腺炎成分疫苗(MuCV)的接种工作，两剂接种时间间隔13～36个月，可最大限度减低辖区流腮发病率[13]。二是要做好学校和托幼机构入托入学儿童腮腺炎疫苗接种证查验，及时发现漏种儿童，做好补种。三是要学校和托幼机构做好晨午检及因病缺勤登记，及时发现和处置流腮疫情，防止扩散蔓延。四是要加大流腮防控知识宣教力度，特别是在流腮漏报率高的地区，更应提高人群预防意识[14]；同时应加强医疗卫生机构包括CDC的专业技术人员的继续医学教育，以提高他们的传染病防控知识，因为传染病防控知识欠缺，对传染病认识不足，缺乏应有的警惕性等都与医学教育密切相关[15]。