人间布鲁氏菌病（布病）是世界上最常见的人畜共患病之一^[1-3]。我国监测数据显示，布病的发病率急剧上升，在法定传染病排名中从2000年的第16位上升至2014年的第6位。宁夏回族自治区（宁夏）在20世纪80年代到90年代初为布病的低发地区，年发病率＜0.100/10万^[4]，从1995年开始发病率逐渐上升，在2008年后，布病的发病人数呈百位数增长，发病率位居全国第3位，发病增速位居全国首位。既往研究显示，布病的再次“复燃”在流行病学特征方面出现了新的发展趋势^[5]，布病的疫区、受侵人群和流行形式均发生了变化，而这种变化存在明显的时空共变现象^[6]。为探讨宁夏地区布病疫情的时空分布特征及其与主要家畜存栏数的相关性，本研究对2012－2018年宁夏各县（区）布病发病率进行时空分析，从而为采取相应的防制措施提供参考依据。

资料与方法

1.资料来源：2012－2018年宁夏各县（区）布病发病数资料来自宁夏回族自治区传染病报告信息管理系统（数据实时上传中国疾病预防控制中心法定传染性疾病网络直报系统），各县（区）人口数和牛、羊、猪存栏数均来自《宁夏统计年鉴》。

2.统计学分析：用动态数列描述2012－2018年宁夏布病发病率的发展动态，用ArcGIS 10.6软件^[7]分析布病空间自相关性，应用SPSS 23.0软件进行Spearman相关分析，分析2015年宁夏22个县（区）布病发病率与牛、羊、猪存栏数的相关性。

（1）动态数列是一系列按时间顺序排列起来的指标，用以观察和比较事物在时间上的变化和发展趋势。绝对增长量是说明事物在一定时期增长的绝对值，可分为累积增长量：即报告期指标与基期指标之差；逐年增长量：即报告期指标与前一期指标之差。发展速度和增长速度是说明事物在一定时期的速度变化。定基比发展速度：即报告期指标与基期指标之比；环比发展速度：报告期指标与前一期指标之比；定基比增长速度：定基比发展速度－1；环比增长速度：环比发展速度－1。

（2）空间自相关：空间自相关是指疾病在空间上的分布存在一定的相互关系，相邻地区的变量值会相互影响。可分为全局空间自相关（global spatial autocorrelation）和局部空间自相关（local spatial autocorrelation）^[8]。使用GIS软件生成布病发病率可视化地图，并对发病率进行全局自相关和局部自相关聚集性分析，获得的地图既可以对比同年各县（区）的发病情况，也可以纵向比较不同年份各县（区）的发病情况。

① 全局空间自相关分析：采用全局Moran’s I值来度量，反应空间相邻区域属性值的相似度，其取值范围为[-1, 1]。I＞0表示空间正相关；I＜0表示空间负相关；I值接近0，则不存在空间相关关系。I的绝对值越大，空间相关性越强。当|Z|＞1.96，P＜0.05为差异有统计学意义，存在空间自相关^[9-10]。其公式如下：

式中，I为Moran’s I值，n为空间要素的个数，χ_i与χ_j是在单元格i与j处的观察值，χ为所有空间单元属性值的均值，ω_ij指要素i与j之间的空间权重。

② 局部空间自相关分析：采用局部Moran’s I值来度量，用于探测局部区域是否属于高值或者低值聚集，解释为中心（高/低）发病率地区被周围（高/低）发病率地区所包围^[11]。局部空间自相关分析聚集图共有4种类型：高-高聚集（H-H）、低-低聚集（L-L）、低-高聚集（L-H）、高-低聚集（H-L）。其公式如下：

式中字母表示与全局Moran’s I系数相同。

结果

1.时间分布：计算各年发病率，经趋势性检验后，未发现差异有统计学意义（χ²＝0.781，P＝0.377），但根据累积增长量却发现从2012年以后，各年发病率较2012年都有所上升。2012－2015年，宁夏地区布病的发病率逐年增高，2015年达到最高，发病率为43.656/10万。自2015年以后，发病率显示出较小幅度的下降。见表 1。

2.全局自相关分析：以县级单位做全局自相关分析（表 2），2012、2013、2016年Moran’s I值均为正值，说明呈空间正相关分布，存在聚集性。而2014、2015、2017和2018年不存在空间自相关，差异均无统计学意义（均Z＜1.96，P＞0.05）。

3.局部自相关分析：全局空间自相关只能说明研究对象存在自相关，但却无法说明具体的聚集区，只有通过局部自相关来得知某一区域与其临近区域的关系，因此对全局Moran’s I值具有统计学意义的年份进行局部空间自相关分析并生成局部空间自相关聚集分布图（图 1）。从图中可以看出，2012年宁夏“H-H”聚集区包括红寺堡区和同心县，其平均发病率为41.964/10万；2013年宁夏“H-H”聚集区在同心县，其发病率为51.867/10万；2016年不存在“H-H”聚集区，这可能与2016年同心县和红寺堡区发病率降低有关。总体来看，宁夏地区2012、2013以及2016年的布病聚集区相对稳定。

4.布病与牛、羊、猪存栏数的相关关系：因2015年宁夏地区布病的发病率为历年最高，故进一步分析发病率和牛、羊、猪存栏数的相互关系。宁夏各县（区）布病发病率与牛、羊、猪存栏数的分布关系见图 2。Spearman相关性分析结果显示，布病的分布与羊存栏数存在正相关关系（r＝0.692，P＝0.000），说明随着羊存栏数的增加，布病发病率也相应增加；而牛和猪存栏数与布病发病率之间尚未发现相关关系（r＝0.047、0.223，均P＜0.05）。图 2的颜色越深表明布病发病率越高，圆点越大表明牛、羊、猪的存栏数越多。在图 2B中，颜色较深的发病率高的盐池县，其羊存栏数的圆点也同样较大，而颜色较浅的发病率低的惠农区、大武口区以及兴庆区等，其羊存栏数的圆点也相应较小；而在图 2A中，尚未看到颜色越深的布病发病率高的地区，其牛存栏数的圆点越大，在图 2C中，也未看到颜色越深的布病发病率高的地区，其猪存栏数的圆点越大。

讨论

本研究基于2012－2018年布病的监测数据，研究了宁夏布病发病率的时空分布情况，发现宁夏布病具有独特的时空聚类特征。对2012－2018年宁夏布病发病率的发展动态分析发现，2012－2015年宁夏布病呈逐年上升，这与国内其他地区基本一致^[12-14]，在2015年发病率处于历年最高，随后从2015－2018年宁夏布病的发病率呈现一个较缓的下降态势。全局空间自相关研究结果表明宁夏布病在2012、2013和2016年具有空间正相关性，呈现聚集性，而2014、2015、2017和2018年呈随机分布。进一步局部自相关分析表明，宁夏在2012、2013和2016年布病的“H-H”“L-L”聚集区基本稳定。

了解宁夏布病的空间分布对于控制布病发病，减少布病带来的疾病负担具有重要意义，同时也有利于资源的最优化配置^[15]。有助于识别布病发病高风险区域，从而对“H-H”聚集区域采取重点监测的预防方针，而对“L-L”聚集区域及时发现控制疾病发病的有利因素，为布病监测提供指导性作用，并且在一定程度上为相关部门的决策制定提供依据^[16]。本研究表明，2012、2013和2016年宁夏布病“L-L”发病均在相对固定的地方聚集，因此这种聚集状况的有利影响因素有待于进一步研究。

在宁夏布病与牛、羊、猪存栏数关系的分析中，发现羊存栏数与布病发病率存在正相关关系，羊存栏数的增加将引起布病发病的增加。畜间布病的存在是引起布病的主要危险因素，原因是病畜会排出大量病原菌，如果处理不当，会严重污染周围环境，从而给其他动物、接触人员的感染创造条件。但畜间布病却是布病感染的必要不充分条件，因为布病是由多种原因引起的，即使畜间存在布病，若处理得当也可避免引起布病，另外，人们的卫生防护意识、接触人员的身体素质和机体免疫力也是引起布病的相关因素。

基于引起布病的危险因素，应当加强控制传染源、切断传播途径和保护易感人群的防控措施。对于病畜，要及时进行扑杀，且要将其运至远离居民区的地方进行填埋或者焚烧，并要及时处理病畜的排泄物及分泌物，对污染的圈舍也要进行严格的消毒。对于布病的职业相关人员应加强自身防护工作，在接触牲畜时需佩戴眼罩、口罩、长筒橡胶手套以及专门的防护服，且在接触牲畜后及时洗手和消毒。

综上所述，随着宁夏地区畜牧业的发展，加强宁夏布病的预防检测刻不容缓。针对本次研究的布病疫情分布特征，建议应进一步加强宁夏中部地区疫情的检测，对于高发病地区，适当增加防控资源，以遏制布病疫情的蔓延，而对于低发病地区，要及时寻找有利影响因素来减少布病的发生。

利益冲突  所有作者均声明不存在利益冲突所有作者均声明不存在利益冲突