口岸媒介生物性传染病传入风险分析系统的建立与评价

徐 军，肖云霞，2，王安东，2，张江国，莫善明，徐新龙，戴 莉

（1. 阿拉山口出入境检验检疫局综合技术服务中心，新疆阿拉山口 833418；2.石河子大学动物科技学院，新疆石河子 832000）

口岸媒介生物性传染病传入风险分析系统的建立与评价

徐 军1，肖云霞1，2，王安东1，2，张江国1，莫善明1，徐新龙1，戴 莉1

（1. 阿拉山口出入境检验检疫局综合技术服务中心，新疆阿拉山口 833418；2.石河子大学动物科技学院，新疆石河子 832000）

本研究参照澳大利亚和新西兰联合开发的风险管理标准AS/NZS4360，对比国内多个风险分析系统，结合调研、统计分析等方法，建立了口岸地区媒介生物性传染病风险分析系统，并应用本系统对阿拉山口口岸地区鼠类携带棘球蚴病的传入风险进行了分析。与其他多个系统对比评价发现，本分析系统结果较为可信，且分析过程较其他相对比系统更严谨科学。但实际情况复杂多变，风险分析系统需要实际检验过程中不断优化改善。

风险分析；口岸；媒介生物性传染病；鼠类携带棘球蚴；阿拉山口

随着国际贸易的不断增加，出入境人员、交通工具、集装箱及货物等在国际的流动也日益频繁。国家质检总局的统计，2014年上半年的出入境货物共计487.34万批，货值6 108.8亿美元的产品物流通过我国口岸，每年的出入境人数也超过2亿人次[1]。国际贸易因其打破了疾病暴发流行的生态地域限制，为媒介生物及媒介生物性传染病的传播，提供了人为的便利[2，3]，并已经成为传播传染病及其媒介生物的最主要途径[4，5]。因此，媒介生物性传染病通过口岸贸易往来以及媒介生物的分布流行，使其传入我国的风险也进一步加大。

目前，国内在实际检验检疫工作中，对于媒介生物性传染病仍没有建立完备的风险评估体系，对风险的评估多依靠定性主观判断等方法，存在较大的误差[6，7]；另外，虽也建立了一些评估系统，但因存在系统性缺失、客观性不足等缺点，而造成缺

乏全面性、操作烦琐、普及率低等问题[8]。这一方面可能导致低风险传染病被高估，而造成人力、物力等资源的浪费；另一方面也可能导致因高风险传染病被低估，而使得因防控力度不够使传染病传入国内，造成重大经济损失甚至严重危害人畜健康。因此建立国境口岸媒介生物性传染病传入风险分析系统显得尤为必要[6，9]。

本研究在参考澳大利亚和新西兰联合开发的风险管理标准AS/NZS4360以及澳大利亚通用监控系统的基础上，借鉴相关风险分析参考文献[3-5，10-12]，并结合我国口岸地区目前所面临的输入性媒介生物性传染病的问题，建立了一种半定量多指标综合评估的方法：通过专家咨询法（Delphi）确定对媒介生物性传染病传入可能性与危害程度两方面的影响因素后，运用调查问卷结合AHP层次分析法确定各影响因素权重系数，最终通过风险概率的计算结果得到风险评估等级。运用此方法对阿拉山口口岸地区鼠类携带棘球蚴病的传入风险进行分析评估，对建立的方法进行实际操作展示，以指导实际应用；同时与国内多部分析系统进行对比评价，对本研究的优势进行分析。

1 风险分析系统的建立

1.1 建立影响因素综合评估体系

媒介生物性传染病作为传染病的一部分同样被传染源、传播途径和易感动物3个流行条件制约，同时也受自然因素以及社会因素的影响[6，12]。本研究从媒介生物性传染病传入的可能性及其危害程度两方面分别进行影响因素评估：媒介生物性传染病传入可能性影响因素即对病原体通过媒介生物传入国内的可能性造成影响的因素，一般由传染源的状况、传播途径和防控措施及国家信息三个方面确定；危害程度影响因素即对病原体传入国内后造成的危害具有影响作用的因素，一般由易感动物、国际影响和传染源的危害三个方面确定。而后通过专家咨询法（Delphi）确定其对媒介生物性传染病的影响因素[10-11，13-15]。其遵循的原则如下：

1.1.1 重点选择。影响虫媒性传染病传播流行的因素很多，有针对性地选择主要的影响因素进行分析显得尤为重要。若列出所有因素，对于非主要及某些不确定性的因素过多将会冲淡主要影响因素对风险的评判，从而影响风险评估的准确性。因此，应重点选择起决定作用的因素[10，12]。

1.1.2 相对稳定。某些影响因素，如所有生物种类、自然空气因素等，存在大量的不确定因素，无法统计搜集。则应选取相对稳定的因素进行评价[15]。

1.1.3 易于评价。对于不易进行分类评估、统计分析的影响因素，可不必细分类别[11-12]。

1.1.4 相对独立。确定影响因素时应尽量保证其相对独立性。若选择的因素在内涵上有交叉，会加重或削弱该因素的权重，从而影响评估结果的准确性。因此在确定影响因素时应尽量避免交叉，或加以特别规定进行处理，保证其相对的独立性[10，16]。

1.1.5 合理协调。虫媒传染病作为一类传染病有其特性，但由于背景的不尽相同，可能造成其复杂多样性，本指标体系的建立适用于宏观的评价。对于具体病种、地区可对其作必要调整，设置合理缺项[11]。

1.2 建立权重系数综合评估体系

由于媒介生物性传染病风险评估过程中存在着不确定因素，因此需要取得相关数据资料，以避免个人主观因素的影响。本研究通过专家咨询法（Delphi）与层次分析法（AHP）进行权重系数[17]的确定以增加此风险分析系统的科学、合理性[13]。步骤如下：

1.2.1 确定各影响因素比重。将上述通过咨询法确立的影响因素分不同层面列出，并对口岸媒介生物监测实验室专家和一线医学媒介生物监控工作人员分别填写调查问卷，利用Delphi法定量统计分析后以确定系统中可能性和危害严重程度的各影响因素比重。

1.2.2 确定各影响因素权重。本研究运用由北京欣晟允软件技术有限公司开发的层次分析法辅助软件yaahp开展分析。其分析结果可准确算得各影响因素的权重系数。

1.3 建立风险水平综合评估体系

本研究通过对多标准及评估模型进行选择对比，将可能性和危害严重程度两方面影响因素定性分析进行风险概率评分，结合权重系数计算，确立媒介生物性传染病评估模型风险水平[18]，最终利用矩阵分析指数表确立风险水平综合评估值。

1.3.1 确定可能场景。首先确立第二准则层每个影响因素的不同可能场景，进行定性描述后对每个可能场景分配定性估计值如“低”、“中等”、“极高”等（表1）[18]。

表1 风险评分等级和概率范围

1.3.2 概率计算。对第二准则层每个影响因素设立多个可能影响途径，分配定性后的估计值在相应表2中可确立其概率范围，将概率范围中值相乘后为该影响因素可能概率PEij。以此类推可得第二准则层所有影响因素可能概率，将所得该因素的概率与其相应的权重系数Wij相乘（式（1））可得该因素概率值PPEij。

即：PPEij=PEij×Wij（1）

第一准则层影响因素可能概率PEi为与其相应的第二准则层因素概率值PPEij之和（式（2））。分别求得第一准则层影响因素可能概率PEi后并将其与相应权重系数Wi结合相乘可得第一准则层因素概率值PPEi（式（3））。将第一准则层因素概率值PPEi累加最终得出该风险概率值（式（4）），根据概率分级范围表（表2、3）可确定可能性和危害程度两风险概率。

即：PEi=∑PPEij（2）

PPEi=PEi×Wi（3）

P=∑PPEi（4）

1.3.3 风险评估确定。通过以上方法求得媒介生物性传染病传入可能性及其危害程度两方面的风险概率并取得相应风险概率等级（表2、3）。参照AS/NZS4360:2004矩阵评估指数表[19]（表4），找到相对应风险概率最终确定风险水平。

表2 媒介生物性传染病传入可能性概率分级

表3 媒介生物性传染病传入危害程度的界定

表4 风险水平矩阵评估指数表

2 风险分析系统的应用与评价

2.1 风险分析系统的应用

运用前文建立的风险分析系统对阿拉山口口岸地区鼠类携带棘球蚴病的传播风险进行分析评估。经调查，该口岸鼠类媒介生物本底情况及驻地状况基本清晰，一是与其同类自然疫源地之间没有自然的屏障；二是随着经济贸易的不断发展、人员的频繁流动、现代化的交通工具等因素使口岸地区的病媒生物种类及密度趋于一致。通过专家咨询法确立了该病传入可能性及其危害程度两方面影响因素，并通过调查问卷及yaahp软件计算得各因素权重系数（表5）。

对第二准则层每个影响因素进行定性描述，同时对其描述的不同概率途径进行概率分配。例：病原体阳性情况（图1）。

由图1可知该影响因素病原体阳性情况可能由3个可能路径（L1～L3）的概率影响。我们通过描述分配定性，分配后得估计值：L1为“中等”；L2为“低”；L3为“高”。根据概率范围表（表2）可取其中值，即：L1=0.50、L2=0.155、L3=0.80；PE11=6.2×10-2。以此类推算得第二准则层各影响因素的可能概率PEij，表6示范了这一评估步骤：

表5 媒介生物性传染病发生风险评估影响因素体系

图1 病原体阳性情况概率范围计算示例

表6 计算第一准则层可能性影响因素传染源的可能概率

表7 媒介生物性传染病传入可能性风险概率

算得第一准则层的可能概率PEi后，同理根据式（3）和式（4）可得最终风险概率P。如表7：

根据相同步骤可得鼠类媒介生物性传染病传入危害程度概率P=7.9×10-2=“较小”。将两结果根据矩阵指数表可能性“罕见”与危害程度“较小”相对应风险概率最终确定风险水平为L：低危险度风险。

2.2 风险分析系统的评价

本文对应用案例经国内多个不同模板[11，12，19]进行风险分析评价，并分别对分析过程以及分析结果进行对比，以通过分析过程及结果对比来评价本研究所建立的系统。对该示例进行分析，最终评估分析结果相似，证明了其本分析系统结果的合理性。另外，本系统分析过程同王育兵[12]等建立的模型分析过程比较可发现，王育兵等建立的模型分析在最终分析计算过程中概率运算部分，其评估模型最终计算表达式为：y = ΠXiWi。由于其评估指标通过连乘最终计算得评估指标。忽略了可能由于某一低指标因素判断为0，从而造成整个系统得到结果为0，最后对最终风险评判结果造成很大误差。本系统的考虑各因素的影响，参照Nigel Perkins与Mark Stevenson[18]提出的半定量评分概率计算方法，避免了连乘所造成的因某一因素偏差而对结果影响过大的缺陷。

同时与曹鑫等[11]和曾晓芃等[19]所述方法相比较，其方法在风险界定过程当中多运用定性分析，虽简单易于评价，但对可变性和不确定性的无法全面表述。本模型对各因素概率进行概率分配，增加定量分析过程扩大了信息量，对变量的阐述更加全面。在权重系数的确立及概率运算部分，本研究系统运用了AHP层次分析法软件，更加准确的得到每个影响因素的权重系数。减少了主观因素的影响，分析结果更具有科学、客观性。因此本文中建立的方法，科学有效，具有推广价值。

3 结论

媒介生物性传染病风险评估多见于口岸等媒介生物性传染病输入风险的评估，由于相关国内资料有限，本研究系统参考AS/NZS4360等多个国际标准，对比结合多部国内风险分析模型。采取Delphi、AHP等方法，并借助层次分析法辅助软件yaahp计算得到各影响因素权重系数。同时借鉴新西兰麦思大学Niegel Perkins和Mark Stevenson博士著风险分析培训手册[18]中半定量评分方法进行概率范围的判定。最终建立了比较完善的风险评估分析模型。各口岸可根据自身情况，对于媒介生物性传染病进行风险评估，其结果对预警、控制有一定的理论意义和现实意义，对制定病媒生物防治对策和卫生资源的分配提供理论依据。

同时借助口岸媒介生物性传染病传入风险开展分析对本系统进行实践和评估，经过对阿拉山口口岸地区鼠类携带棘球蚴病的风险分析，验证了本分析系统可行性。通过与其他风险分析方法的对比，风险分析结果具有一致性，即本系统分析结果可信，同时对比分析过程，也凸显了本系统的客观、严谨性。实际运用过程中因情况复杂多变，可能会存在某些不足，需在应用中根据反馈不断优化。

[1]国家质检总局.全国货物检验检疫概况[Z].国家质检总局信息公开[2014].北京：国务院，2014.

[2]Lines J. Chikungunya in Italy：globalisation is to blame，not climate change [J]. BMJ：British Medical Journal，2007，335（7620）：576.

[3]Tatem AJ，Rogers DJ，Hay SI. Estimating the Malaria Risk of African Mosquito Movement by Air Travel[J]. Malaria journal， 2006，5（1）：57.

[4]Tatem AJ，Hay SI，Rogers DJ. Global Traffic and Disease Vector Dispersal[J]. Proceedings of the National Academy of Sciences of the United States of America，2006，103（16）：6242-6247.

[5]Tatem AJ，Rogers DJ，Hay SI. Global Transport Networks and Infectious Disease Spread[J]. Advances in parasitology，2006，62：293-343.

[6]马红梅，陈海婴，熊志伟，等. 媒介生物及媒介生物性疾病风险评估技术研究进展[J]. 中国媒介生物学及控制杂志，2012，23（6）：592-595.

[7]李翔，戚立松，戴明，等. 船舶外来医学媒介生物传入风险分析初探[J]. 中国国境卫生检疫杂志，2012，（4）：261-265.

[8]柳小青，陈海婴，马红梅，等. 大型活动重点行业蜚蠊侵害与风险评估研究及应用[J]. 中国媒介生物学及控制杂志，2014，25（1）：51-54.

[9]郭齐，孙向东，刘拥军，等.国内动物卫生风险分析研究评述[J].中国动物检疫，2011，28（9）：42-44.

[10]East IJ，Wicks RM，Martin PAJ，et al. Use of a Multi-Criteria Analysis Framework to Inform the Design of Risk Based General Surveillance Systems for Animal Disease in Australia[J]. Preventive veterinary medicine，2013，112（3）：230-247.

[11]曹鑫. 国境口岸输入性虫媒传染病及媒介生物风险分析和风险管理系统研究[D].北京：中国人民解放军军事医学科学院，2010.

[12]王育兵，张丽芬，窦跃龙，等. 某部媒介生物性传染病风险评估体系的建立与应用[J]. 中华卫生杀虫药械，2014，（2）：119-121，125.

[13]徐蓓蓓，李辉，鹿文军，等. 进出口商品风险评价中风险因素权重的确定[J]. 检验检疫科学，2008，18 （3）：27-30.

[14]Standrads Australia，Standrads New Zealand. AS/NZS 4360.Australian/New Zealand Standard Risk Management[S]. Sydney：Standrads Australia International Ltd，2006-06-14.

[15]闰俊平，魏伟，陈溥言，等.动物卫生风险评估方法研究[J].上海畜牧兽医通讯，2009，12（2）：57-59.

[16]张凡建. 动物和动物产品进口风险分析及其要素研究[D]. 北京：中国农业大学动物医学院，2005.

[17]范春彦，韩晓明，汤伟华.AHP中专家判断信息的提取及指标权重的综合确定法[J].空军工程大学学报（自然科学版），2003，4（1）：65-67.

[18]Nigel Perkins，Mark Stevenson. 动物及动物产品风险分析培训手册[M]. 王承芳译. 北京：中国农业出版社，2004.

[19]曾晓芃，付学锋，张勇，等. 2008年北京奥运会病媒生物风险识别与评估方法研究[J]. 中国媒介生物学及控制杂志，2009，20（2）：98-101.

Establishment of Risk Analysis System for Introduction of Vector-borne Infectious Diseases via Border Ports and its Application and Evaluation

Xu Jun1，Xiao Yunxia1，2，Wang Andong1，2，Zhang Jiangguo1，Mo Shanming1，Xu Xinlong1，Dai Li1
（1. Ala-shankou entry-exit inspection and quarantine bureau of the P.R. China，Alashankou 833418; 2. College of Animal Science and Technology，Shihezi University，Shihezi 832000）

A risk analysis system on introduction of vector-borne infectious diseases via border ports was established in the study on the basis of Australian/New Zealand Standard Risk Management（AS/NZS4360） and several domestic risk analysis systems in combination with investigation and statistical analysis. The risk analysis system was used to assess the risk of introduction of echinococcosis in rodents via Ala-shankou port area. Compared with other risk analysis systems，the result of the system described in this paper was genuine and believable，and the parsing process was more rigorous and reasonable. However，the risk analysis system that established in the study needs further optimization to face the complicated practical situation.

risk analysis；border port；vector-borne infectious disease；echinococcosis in rodents；Ala-shankou port

S852.65

C

1005-944X（2015）05-0013-05

新疆维吾尔自治区自然科学基金（201442137-1）

徐军