基于Delphi和AHP的羊布鲁菌病预警指标体系研究

于丽敏 ,兰邹然 ，黄 芳 ，王 翠 ，杨 晶 ，刁有江

（1.山东农业工程学院，山东 济南 250100；2.山东省动物疫病预防与控制中心，山东 济南 250100）

布鲁菌病（Brucellosis）是由布鲁菌感染引起一种人兽共患慢性传染病，广泛分布于世界各地，常发生于羊、牛、猪等多种家畜及野生动物，也传染给人，可以引起哺乳动物发热、关节疼痛、流产、睾丸炎和疲劳等症状[1]。布鲁菌中有强致病力的主要有6种，以羊布鲁菌对人和动物的致病力最强，猪、牛布鲁杆菌次之[2]。布鲁菌可以经消化道、呼吸道、结膜、皮肤损伤、胎盘等多种途径穿透动物机体并快速进入血液，通过血液循环侵入到肝脏、脾脏、淋巴结等器官的细胞内并定居，由于其寄生特性有助于逃避宿主免疫反应（先天性或适应性）和抗生素的作用，持续对动物机体造成损害[3,4]。

布鲁菌感染作为一种人兽共患疫病，全球每年约有50万人被确诊患该病，但实际比报告病例要高出一个量级以上，严重危害人类健康和畜牧业经济的发展[5-7]。在上世纪八十年代，我国布鲁菌病疫情曾得到很好的控制，但进入本世纪以来，我国多个省份对布鲁菌病进行了详细报道，其中动物源性布鲁菌病的发生主要集中在全国各农牧区域[4]，且发病率呈逐年上升的趋势，若不加以防范极有可能存在潜在流行的风险[5-7]。

近年来我国养羊业发展迅速，截至2021年底我国肉羊存栏量已远超3亿只，占世界总存栏量的16.67%，羊肉产量占世界的33.33%[8]。养羊业的快速发展给羊场疫病防控带来了很大的压力。加强羊布鲁菌病预警体系建设是防控羊布鲁菌病蔓延的有效措施，但动物布鲁菌病的预警系统在我国尚处于初始的探索阶段，布鲁菌病的防控依旧面临着严峻的挑战[9]。构建完善的羊布鲁菌病预警系统，对人民生命健康和养羊业的可持续发展具有重要意义。本文利用德尔菲法和层次分析法，从羊群、饲养人员、环境因素和既往发病等方面，对羊布鲁菌病预警指标体系进行初步探讨，以期对我国布鲁菌病的防控提供参考依据。

1 方法

1.1 影响因素分析

通过文献调研和现场调查，分析影响羊布鲁菌病预警指标的因素。首先以 “羊”“布鲁菌病”(brucellosis)为关键词在“中国知网”“维普”“万方数据医药信息镜像系统”和“PubMed生物医学数据库”检索有关羊布鲁菌病研究现状的文献，阅读筛选出有关布鲁菌病流行病学、预警指标体系的文献，对获取的到文献整理，筛选出与羊布鲁菌病密切相关的文献[1,2,9]。同时，在山东峰祥畜牧种业科技有限公司、郓城县大鹏农牧科技有限公司、嘉祥县种羊场开展现场调研，通过文献调研和现场座谈考察，确定饲养动物羊[10-12]、饲养人群[13,14]、环境因素[10-12]、既往发病水平[10-12]等 4 个方面为影响羊布鲁菌病预警的主要因素。

1.2 德尔菲法

德尔菲法（Delphi）为设计问卷发送给专家，从而征询专家意见，其特点是匿名反馈、反复进行[15]。其实施步骤为：首先，根据既有研究内容，确定专家小组，专家人数一般为8～20人，10人以上的专家评审可以取得较为理想的精准度；其次，根据影响因素分析结果，设计调查问卷并向各位专家发送问卷，设计问卷时要说明研究背景材料及填写问卷注意事项；再次，按时回收专家问卷，初步整理所有专家意见，开展第二轮问卷发放[16]；最后，回收二轮专家问卷，整理分析专家意见[16,17]。

本研究专家选取原则为熟悉羊布鲁菌病知识，在该领域工作10年以上，有一定的研究基础。根据选取原则，共邀请13名专家，分别来自畜牧养殖、兽医、动物防疫等专业领域。根据已有研究和影响羊布鲁菌病的因素分析结果，设计专家问卷并发放问卷，确定羊布鲁菌病暴发早期预警指标的专家依据。

1.3 层次分析法

为了获取羊布鲁菌病各指标的权重，采用层次分析法(AHP)，各指标的权重分析基于AHP的1-9标度相对重要性打分，将相对重要性评估分为1分至9分的9个等级，将已构建的指标发给专家打分[15]。过程如下：

基于判断矩阵构造的AHP(如公式(1))，对指标两两逐对比较，比较结果写成矩阵形式。

ui：评价指标，ui∈U(i=1,2,…n)。 uij∶ui对 uj的相对重要性数值（j=1,2,…m），其取值按1-9标度法表示，1到9分表示重要性上升，1表示相等，分数越高，重要性越高。定义uji=1/uij。

根据判断矩阵，计算因子i(Wi1)的权重，如公式(2)-(6)。m是矩阵第i行或第j列的因子个数。λmax是矩阵的最大特征值。Wi是综合评价指标水平中因素i的权重。Wi1为因子i在整体水平中的最终权重，计算为各因子的权重与相应元素的权重的乘积。

基于层次分析法的主观性，需要测试一致性的判断矩阵(如公式(7)所示)。随机一致性比(CR)低于0.1表示一致性满意，否则必须调整矩阵直到符合。CI是判断矩阵一致性的指标 (如公式(8))。当CI=0时，它是完全一致的，当CI！=0时，它是不一致的。

2 指标体系构建结果

2.1 布鲁菌病预警指标体系

通过分析影响羊布鲁菌病预警指标的因素，专家咨询意见结果的整理汇总，获得羊布鲁菌病预警指标体系的4个二级指标、13个三级指标，如表1所示。

表1 羊布鲁菌病预警指标体系

2.2 布鲁菌病预警指标权重结果

专家打分表发放13份，回收13份，回收率100%。对发给专家的指标打分结果通过层次分析法运算后指标权重结果如图1所示，总目标为羊布鲁菌病预警指标体系结果，总权重为1.0000。羊、饲养员、环境因素、既往发病水平4项二级指标的权重分别为：0.610，0.146，0.102，0.141。

图1 羊布鲁菌病预警指标体系指标权重

3 讨论与结论

3.1 布鲁菌病预警指标体系权重分析

研究结果表明，13个三级指标中指标 “未免疫羊只的血清学阳性率大于10%”，指标“羊出现疑似布鲁菌病的症状和体征（长期发热、多汗、口渴、关节炎和跛行、流产或睾丸炎等）”2个指标在指标体系中的分值相对较高，均超过0.100，这2个指标均为我国目前羊布鲁菌病防控中应用较多、价值较高的指标，说明本研究构建的羊布鲁菌病指标体系及权重分布符合我国羊布鲁菌病的实际情况。而有些指标如“布鲁菌病疫区人员进入羊场工作”“从布鲁菌病疫区调入牧草以及其他饲料原料”等指标对羊布鲁菌病的预警价值在指标权重中分值较低，等于或低于0.03，表明这2个指标对羊布鲁菌病预警的价值较低。

从4个二级指标来看，饲养动物指标“羊”占的权重最大，为0.610，本研究的对象羊是羊布鲁菌病的最主要传染源，是羊布鲁菌病的自然宿主，因此在羊上体现的6个方面的指标决定了羊布鲁菌病预警的主要因素；其他3个二级指标“饲养员、环境因素、既往发病水平”权重分布在0.102-0.146之间，说明该指标均为一般指标，提示这三个指标在布鲁菌病的预警工作中起辅助作用，可以作为补充因素考虑。

3.2 基于Delphi和AHP构建指标体系的优势

德尔菲法是一种体现专家经验知识的主观定性方法，需要参与的专家具备代表性、专业性和权威性，这是保证测评结果有效的基础，否则会带来一定偏倚[18,19]。本研究遴选的13位专家分别来自畜牧养殖、兽医、动物防疫等专业领域，熟悉羊布鲁菌病的流行病学和防控知识，具有丰富的现场经验和独到的理论见解。在专家经验的基础上，本研究采用层次分析法定量计算方法为指标赋权重，获取到羊、饲养员、环境因素、既往发病水平4项二级指标的权重分别为：0.610，0.146，0.102，0.141， 实现对专家知识经验的定量化输出和定量平衡。

3.3 结论

本研究将定性和定量方法相结合，初步建立了适用于我国羊场羊布鲁菌病防控的预警指标体系，体现了专家知识经验性和定量客观性，对羊布鲁菌病预警方案的制定，具有一定的参考意义。在实际监测过程中，各地应对羊布鲁菌病预警指标体系进行适当的调整，制定适合当地实际情况的预警方案。