进口鸡肉传带禽伤寒沙门氏菌的定量风险评估

李建军，高梦昭，邓 柯，李艺超，郭瀚民，方海萍，赵明刚

（1.海关总署标准法规研究中心，北京 100013；2.清华大学统计学研究中心，北京 100084；3.清华大学工业工程系，北京 100084）

禽伤寒沙门氏菌（Salmonella gallinarum）以禽类为专一性寄主，可引起禽败血症，造成重大经济损失，因此世界动物卫生组织（OIE）将其列为须通报动物疫病，我国在《一、二、三类动物疫病名录》及《中华人民共和国进境动物检疫疫病名录》中将其列为二类动物疫病。禽伤寒在亚洲、欧洲及南美洲时有发生，其中巴西、乌拉圭、阿根廷等南美国家多年来疫情持续存在，仅2017 年巴西就暴发禽伤寒疫情81 起，发病鸡75.2 万只[1]。鸡肉是我国重要的进口肉类产品之一，2018 年进口量达50.2 万t，2019 年达80 万t，其中南美是最主要的进口来源地[2-3]，因此进口鸡肉传带禽伤寒沙门氏菌的风险不容忽视。为此，本研究分析了进口鸡肉传带禽伤寒沙门氏菌的各类风险因子，以此建立数学模型，并基于2017 年巴西禽伤寒沙门氏菌流行数据，对进口鸡肉传带禽伤寒沙门氏菌的风险进行定量评估。

1 材料与方法

1.1 场景分析

鸡肉生产要经过鸡群选择、宰前检测、屠宰、宰后检测、分割、冷冻、包装、运输等过程。被屠宰的鸡带菌或鸡肉在生产过程中被病菌污染，可能导致进口鸡肉产品传带病菌。进口鸡肉传带禽伤寒沙门氏菌的风险场景见图1。

1.2 风险参数确定

以每批次进口鸡肉为评估单元。在特定批次鸡肉中，只要任何组分携带有禽伤寒沙门氏菌，则将该批次鸡肉传带疫病风险概率视为100%。各风险因子参数见表1。

1.3 计算模型与分析方法

（1）在本研究中，鸡群数量为n，鸡的平均质量为m，每批鸡肉质量为mt，则在每个鸡群选择屠宰鸡的数量为

式中k=1，2，…，n。

（2）鸡群疫病流行率为p1，从感染鸡群随机选择鸡制成的鸡肉产品传带禽伤寒沙门氏菌的概率为q0，则该批鸡肉携带疫病的概率为

（3）定义：

X1=1 表示从感染鸡群选择某只鸡携带病菌，0表示未携带病菌；

图1 进口鸡肉传带禽伤寒沙门氏菌的风险场景

表1 进口鸡肉传带禽伤寒沙门氏菌的风险参数

X2=1 表示该鸡在屠宰后携带病菌，0 表示未携带病菌；

X3=1 表示分割后的鸡肉携带病菌，0 表示未携带病菌；

X4=1 表示鸡肉在冷冻后携带病菌，0 表示未携带病菌；

Y1=1 表示该鸡在宰前检验中被认为是病鸡，0表示不被认为是病鸡；

Y1=1 表示该鸡在宰后检验中被认为是病鸡，0表示不被认为是病鸡。

根据全概率公式

将式（4）（5）（6）带入式（3），并化简，即得

利用上述数学模型及确定的风险参数值，采用固定参考值和蒙特卡洛方法，模拟分析进口鸡肉携带禽伤寒沙门氏菌的风险水平，以及不同风险因子对进口鸡肉携带禽伤寒沙门氏菌风险的影响。

2 结果

2.1 固定参考值分析

将表1 中的参考值带入式（2）和式（8），得q0=0.367 4，p=0.001 0，即在设定的参考值条件下，进口鸡肉携带禽伤寒沙门氏菌的概率为0.1%，也即每1 000 批次进口产品中携带该病菌的批次数期望值为1。

2.1.1 各风险因子对感染鸡群生产鸡肉带菌率的影响 通过固定其他参数值，分析表1 中相关参数对感染鸡群生产的鸡肉产品传带禽伤寒沙门氏菌概率[式（8）中q0]的影响。结果显示，在设定的风险参数取值范围内，q0值随病菌感染率和鸡肉带菌率的上升而线性上升（图2～3）。当鸡群感染率达100%时，产品携带病菌的概率为48.6%；当感染鸡鸡肉带菌率为100%时，产品携带病菌的概率达73.4%。其余风险参数对q0值影响很小。

图2 鸡群感染率对q0 值的影响

图3 感染鸡的鸡肉带菌率对q0 值的影响

2.1.2 各风险因子对随机挑选鸡群生产鸡肉产品带菌率的影响 利用式（2），进一步分析各参数对来自随机挑选鸡群每批次鸡肉产品携带禽伤寒沙门氏菌概率（p）的影响。结果显示：在本研究设定的条件下，仅来源鸡群数量、鸡群疫病流行率对p值有较大影响（图4～5），其余参数影响很小。当每批次鸡肉来源鸡群数量从1 增加至4 时，p值从0.000 5 增加至0.002 0。当鸡群疫病流行率从0增加至1 时，其p值从0 增加至1。

2.2 蒙特卡洛模拟

通过蒙特卡洛方法，从表1 各分布中抽样出参数组合，并利用全概率公式计算鸡肉携带禽伤寒沙门氏菌的风险概率。在进行1 000 次抽样后，得到的风险概率频数直方图如图6 所示。根据该模拟结果可以计算得到进口鸡肉携带禽伤寒沙门氏菌的风险概率均值为0.001 2，标准差为0.001 4。

图4 每批鸡肉来源鸡群数量对p 值的影响

图5 鸡群疫病流行率对p 值的影响趋势

图6 鸡肉携带禽伤寒沙门氏菌的风险概率直方图

从概率分布范围看，进口鸡肉携带禽伤寒沙门氏菌的风险概率25%在0.000 3 以下，50%在0.000 8 以下，95%在0.003 8 以下（表2）。

表2 进口鸡肉携带禽伤寒沙门氏菌的风险概率分位数

3 讨论

本研究建立了进口鸡肉传带禽伤寒沙门氏菌风险的定量评估模型，并以进口批次为分析单元，分析了各种风险因子对进口鸡肉传带禽伤寒沙门氏菌风险的影响，发现生产每批鸡肉所用的鸡只来源鸡群越多、鸡群疫病（禽伤寒）流行率越高，其传带禽伤寒沙门氏菌的概率就越高，特别是鸡群疫病流行率对传带概率的影响幅度非常大。其他风险因子如鸡群感染率、屠宰场感染率、冷冻灭活率、宰前抽检率、宰后抽检率等对进口鸡肉传带禽伤寒沙门氏菌的总体风险基本无影响。分析原因，主要在于生产每批次鸡肉所用的鸡只数量很大，即在风险概率计算模型p=1-{1-p1[1-（1-q0）nk]}n中，nk值很大。由于鸡群感染率、屠宰场感染率、冷冻灭活率、宰前抽检率、宰后抽检率等因素仅对q0值产生影响，因此不论这些参数取值如何，q0值都会在0～1之间，也即（1-q0）nk值均会非常接近于0。故进口鸡肉传带禽伤寒沙门氏菌的概率基本不受这些参数变化影响，而主要受鸡群疫病流行率（p1）和鸡群数量（n）影响。

当产品来源地鸡群禽伤寒沙门氏菌流行率保持在巴西2017 年的流行水平时，根据固定参考值法得到进口鸡肉传带禽伤寒沙门氏菌的风险概率值为0.1%，通过蒙特卡洛模拟得到的风险均值为0.12%，即总体风险在0.1%左右。我国进口鸡肉批次多、数量大，禽伤寒沙门氏菌输入风险不容忽视。根据本研究结果，鸡群一旦被禽伤寒沙门氏菌感染，通过降低鸡群感染率以及宰前检验、宰后检验、屠宰过程中交叉污染控制、储藏冷冻等后续措施降低鸡肉传带疫病风险的效果非常有限。因此，要降低进口鸡肉传带禽伤寒沙门氏菌风险，关键是控制输华鸡肉来源鸡群的疫病感染率。为此，建议采取以下两方面措施：一是尽量从禽伤寒非疫区进口鸡肉；二是如需从禽伤寒流行区进口鸡肉，应通过双边协议等形式要求输出国政府及企业加强疫病监测，保证输华鸡肉来源鸡群未受禽伤寒感染。

4 结论

禽伤寒沙门氏菌是较严重的禽类传染病病原，其随进口鸡肉的传入风险主要受来源鸡群数量及鸡群疫病流行率影响，尤其是鸡群疫病流行率。当进口鸡肉来源鸡群禽伤寒沙门氏菌流行率在0.05 的水平时，其传带病菌的风险概率在0.1%左右。鉴于该疫病在世界各地多发，禽伤寒沙门氏菌随进口鸡肉传入我国的风险不容忽视。为此，建议通过与输出方监管部门合作，加强对输华鸡肉来源养殖场的疫病防控与监测，避免输华鸡肉来源鸡群受到感染。