2020-2022年食用农产品安全监督抽检数据分析.pdf

40 食品安全导刊 2023年10月（下）食安管理20202022 年食用农产品安全监督抽检数据分析彭雯婧，罗 斌*，曹喜春，李越凡，熊 铮，卢素珍（江西省检验检测认证总院食品检验检测研究院，江西南昌 330001)摘 要：本研究采用统计学方法对食用农产品的抽检情况进行了多维度分析，并通过不合格率高位值分级、抽检情况二项分布检验、趋势分析等方法评估食品品种总体风险情况，利用卡方检验对抽检不合格的相关因素进行了验证，对食品品种进行风险画像，提高食品安全抽检计划的科学性和精准性。关键词：监督抽检；风险因素；不合格因素；食用农产品Analysis of Sampling Data of Edible Agricultural Products Safety Supervision in 20202022PENG Wenjing,LUO Bin*,CAO Xichun,LI Yuefan,XIONG Zheng,LU Suzhen(Jiangxi Provincial Institute of Food Inspection and Testing,Nanchang 330001,China)Abstract:This study used statistical methods to conduct multi-dimensional analysis of the sampling of edible agricultural products,and evaluated the overall risk of food varieties through methods such as high value grading of non conformities,binomial distribution test of sampling situation,and trend analysis.The chi-square test was used to verify the relevant factors of unqualified sampling inspection,and the risk portrait of food varieties was carried out to improve the scientificity and accuracy of food safety sampling inspection plan.Keywords:supervision sampling;risk factor;unsatisfactory rate;edible agricultural products食品安全与人们的身体健康息息相关，保障食品安全是社会稳定的重要举措。农产品作为人们日常生活中必不可少的食材，每年抽检量在40%以上，不合格率较高。本研究以 20202022 年食用农产品监督抽检数据为例，引入二项分布检验、卡方检验等多种统计学方法对食品农产品监督抽检数据进行分析，量化食品品种风险值，结合场所、购进/生产月份、不合格项目分析，为食品品种进行精准画像，以期为抽检计划执行提供科学支持。1 数据来源通过国家市场监督管理总局食品安全抽检公布结果查询系统收集 20202022 年食用农产品监督抽检数据（数据来源：https:/ 方法与步骤本研究用 Microsoft Excel、python、SPSS 等软件对下载的数据进行录入处理和统计分析。2.1 食品品种风险值的多维度分析与量化食品品种的不合格率被广泛用于食品安全风险评估，该指标是一个动态变量，在特定时间段内，不合格率存在较大的随机性，这给分析带来了不确定性。采用二项分布检验和基于食品品种不合格率的多维度分析方法，全面评估食品品种的安全风险。2.1.1 不合格率风险分级通常使用不合格率或不合格度等指标来表征食品安全风险。以全年总体不合格率作为衡量的标尺，基金项目：江西省市场监督管理局科技项目（GSJK202211）。作者简介：彭雯婧（1984），女，江西吉安人，本科，工程师。研究方向：食品监督抽检。通信作者：罗斌（1968），女，江西丰城人，本科，正高级工程师。研究方向：食品业务综合管理。E-mail:。Oct.2023 CHINA FOOD SAFETY 41食安管理将食品不合格率进行分级赋分，从而体现该食品品种在全年中不同的风险程度。当食品品种的不合格率所属级别越高，就可以认为该食品类别的风险较高，并给予相应的高分值。2.1.2 二项分布检验风险分级二项分布检验是一种适用于二分类问题的统计方法，用于检验二项分布中比例是否等于某个特定值。本文将合格和不合格作为二分类变量，然后使用二项分布检验来计算每个食品品种合格率的置信区间。通过分析置信区间的宽度和位置，可以将食品品种划分为低风险、中风险、高风险 3 个等级。如果置信区间较宽，表示合格率的不确定性较大，食品品种的风险也较高。2.1.3 不合格率的趋势分级将 20202022 年各类食品的不合格率通过最小二乘法进行线性拟合，找到一条与实际观测数据的残差平方和最小的直线，进而确定最佳拟合直线的参数，y kx b，y 是因变量，x 是自变量，k 是斜率，b 是截距。k 在一定程度上表达了该食品品种不合格率在近几年的整体发展趋势，一般来说 k 0 表示食品不合格率呈上升趋势，有风险增高的趋势，且k值越高，表示增长的趋势越高；k 0 表示该类食品不合格率整体呈下降趋势。2.1.4 食品品种风险量化将 3 个因素结果进行分级量化，综合评价各食品风险性。以比较范围数据的总体不合格率作为不合格率风险分级的标准，将品种不合格率与总体不合格率的比值以及不合格率 2%作为标准进行划分。二项分布检验可以将食品品种风险分为高、中、低风险，增长趋势按 k 值进行等级划分。等级赋分具体分值情况见表 1。综合上述风险因素，根据实际赋予不同权重对各品种风险值进行量化计算，为抽检计划的制定和执行提供指导。风险值计算公式为1niiiTv=（1）式中：T 为食品品种总风险值；vi为各项分值；i为各分项权重（根据工作重点确定，如对该品种分析的关注度等），i=1，2，。2.2 食品品种抽样月份、抽样区域、检验项目的风险点分析卡方检验是一种基于分布的假设检验方法，用于综合分析食品品种、生产月份、抽样区域的不合格率的显著差异性，并评估它们在食品品种中的不同风险表达方式，可结合不合格检验项目为食品品种进行预警。表 1 风险指标与分值对应表风险指标指标区间赋分值不合格率分级（a 为总体不合格率）x 2a5a x 2a42%x a30 x 2%2x 01二项分布检验高风险5中风险3低风险1不合格趋势分级k 0.0550.03 k 0.0540.01 k 0.0330 k 0.012k 013 结果与分析3.1 食品品种风险值根据表 1 对 20202022 年食用农产品的整体数据进行赋分，假定所有食品抽检的各分项权重 均为 1，计算得到各食品品种的风险值见表 2。表 2 中食荚豌豆和葱根据二项分布检验被归为低风险品种，不合格批次分别占总批次的 1/8、2/39，不合格率得分较高，可能是因为抽检量和不合格批次量较少导致二项分布检验与不合格率分级差别较大，后期监管中可以适当提高抽检批次进行风险验证。说明风险分析应避免单一指标的片面解读，需要考虑相关因素的综合影响。13 个高风险品种中，结合不合格率、斜率得分、二项分布检验分值分析，其他水产品（牛蛙）、海水蟹、豇豆和姜等几类产品在 3 个维度的结果都属于高风险分值品种，应重点加强监管；豆芽的不合格率得分和二项分布检验分值都属于高风险组别，斜率得分属于降低组别，可在重点监管范围内降低频次。3.2 食品品种时间、抽样场所、检验项目的风险点分析选取高风险品种进行风险点分析，包括豇豆、姜、海水蟹、其他水产品、香蕉、海水虾、韭菜、淡水鱼、辣椒、猕猴桃、生干籽类、食荚豌豆及豆芽。卡方42 食品安全导刊 2023年10月（下）食安管理检验结果及不合格项目统计见表 3。3.2.1 食品品种时间、抽样场所风险点分析食用农产品是非预包装食品，一般以购进/抽样日期统计。食用农产品抽样场所主要有餐饮（如餐馆和食堂等）和流通（如超市、菜市场等）场所。综合分析，韭菜、姜、豇豆、香蕉、其他水产品（牛蛙）和淡水鱼等食品品种绝大部分月份甚至全年的不合格率均具有显著差异性，应在全年加强抽检和监管频率，并在抽样场所对有差异性的品种提高抽检数量与监管力度；食荚豌豆为时间与场所没有显著性差异的高风险品种，可能是抽检量和不合格批次量偏少，应在全年对各场所适当提高抽检量与监管力度进行验证；对于豆芽、猕猴桃、海水虾、海水蟹和生干籽类等部分月份不合格率具有显著差异的食品品种，针对重点时间和场所适当提高抽检量与监管力度。3.2.2 食品品种检验项目的风险点分析农药残留不合格常检出于蔬菜水果中。由于农药以及作物特性不同，部分作物与农药间有紧密关联。4-氯苯氧乙酸钠、6-苄基腺嘌呤对种子发芽、细胞分裂有促进作用，且具有抑制胚根生长的生理作用，不合格样品只检出于豆芽1。氯吡脲对瓜果类植物有促进花芽分化、保花保果、促进果实膨大的作用，能改善猕猴桃的外观，并延长猕猴桃的货架期，只检出于猕猴桃2。腐霉利是一种微毒级的杀菌剂，表 2 各食品品种风险值食品细类不合格率得分二项分布检验斜率得分风险值食品细类不合格率得分二项分布检验斜率得分风险值豇豆55515贝类2125食荚豌豆51511淡水鱼55313豆芽45110淡水虾3115姜55515海水蟹55515芋3339海水虾55414甘薯3137海水鱼3328胡萝卜2125其他水产品55515马铃薯2125橙2327萝卜2114柑、橘3317黄瓜2125李子2125苦瓜2114桃2125辣椒45413枣2125甜椒3137猕猴桃45413茄子2327草莓3339番茄2125葡萄2125韭菜55313香蕉55515葱4138芒果3339大蒜3137火龙果2125茭白3137梨2125莲藕2125苹果2125芹菜3328猪肝2125菠菜3328其他畜副产品2125普通白菜3328羊肉2125油麦菜2125牛肉2114大白菜2114猪肉2114菜薹3137鸡肉2125结球甘蓝2125鸭肉2125花椰菜1124鸡蛋2114生干籽类45312其他禽蛋2114注：风险值总分在 10 分以上为高风险品种，6 10 分为较高风险品种，4 5 分为中等风险品种，3 分为低风险品种。其中风险值为 3 的数据表示不合格率为 0、二项分布检验结果为低风险、趋势持平或下降，因此未在表 2 中展示。Oct.2023 CHINA FOOD SAFETY 43食安管理被大量用于易被真菌和细菌感染的韭菜。豇豆与辣椒因成熟期较短且病虫害较多，不合格农药种类也较多，需要重点关注3。禁用农药甲胺磷和限用农药（毒死蜱、克百威、氧乐果等）也在辣椒和豇豆等产品中被检测出。以上情况主要是由农户在种植环节追求经济效益，过度依赖农药造成的。兽药残留主要检出于水产品中。恩诺沙星是一种用于治疗动物皮肤感染、呼吸道感染的抗生素，在其他水产品（牛蛙）、淡水鱼、海水虾中均有检出，可能是养殖场或者商户为了减少水产品在养殖过程或运输过程中的病死率，使用了含有大量恩诺沙星的饲料或暂养水，导致迁移性污染。此外，禁用的兽药硝基呋喃类药物和孔雀石绿也被检测出于其他水产品（牛蛙）、淡水鱼、海水虾。重金属铅、镉被检出于海水虾、海水蟹、辣椒、韭菜和姜。重金属污染主要是由于农产品在种植、养殖过程中受到环境污染造成的。如海水中的镉超标会造成甲壳类生物重金属的富集4，蔬菜受种植环境的污染造成重金属残留不合格。质量指标（酸价）和生物毒素（黄曲霉毒素 B1）不合格只检出于生干籽类（花生）。花生作为植物油脂的主要原料之一，在种植、储存、运输或加工环节的表 3 食品品种时间与抽样场所卡方结果以及不合格项目汇总表食品品种抽样/购进月份抽样场所不合格项目韭菜1（2=32.19）、2（2=68.07）、5（2=13.86）、6（2=17.43）、7（2=16.30）、8（2=5.61）、9（2=5.43）、10（2=18.76）、12（2=14.86）餐馆（2=4.70）腐霉利、镉、克百威、毒死蜱、铅、氯氟氰菊酯、氧乐果和氯氰菊酯姜1（2=198.72）、2（2=51.43）、3（2=200.55）、4（2=129.75）、5（2=187.92）、6（2=117.84）、7（2=133.66）、8（2=93.92）、9（2=98.12）、10（2=103.38）、11（2=8.13）、12（2=30.14）食堂（2=3.85）铅、噻虫胺、镉、水胺硫磷、吡虫啉、毒死蜱、噻虫嗪、氧乐果和甲拌磷豇豆1（2=27.01）、2（2=19.53）、3（2=45.68）、4（2=110.28）、5（2=383.51）、8（2=6.97）、9（2=65.79）、10（2=78.06）、11（2=41.33）、12（2=18.32）三唑磷、倍硫磷、甲基异柳磷、克百威、噻虫胺、甲氨基阿维菌素苯甲酸盐、灭蝇胺、水胺硫磷、噻虫嗪、氯唑磷、啶虫脒、氧乐果、灭多威、毒死蜱、甲胺磷和氟虫腈辣椒3（2=7.61）、5（2=13.11）、6（2=168.80）、7（2=26.60）、8（2=8.28）、10（2=18.76）、11（2=12.43）噻虫胺、毒死蜱、镉、啶虫脒、氧乐果、噻虫嗪、水胺硫磷、克百威、吡唑醚菌酯、铅、甲氨基阿维菌素苯甲酸盐、甲胺磷、吡虫啉和甲拌磷食荚豌豆多菌灵豆芽1（2=11.50）、5（2=33.26）、9（2=12.21）小食杂店（2=5.64）4-氯苯氧乙酸钠、6-苄基腺嘌呤猕猴桃9（2=27.35）、10（2=5.66）氯吡脲香蕉2（2=4.88）、3（2=10.42）、4（2=40.34）、5（2=27.85）、6（2=15.63）、7（2=11.64）、8（2=48.59）、9（2=106.13）、10（2=5.48）、12（2=11.29）吡虫啉、噻虫胺、噻虫嗪、腈苯唑、氟虫腈、苯醚甲环唑和吡唑醚菌酯其他 水产品2（2=8.72）、3（2=5.42）、4（2=127.18）、5（2=44.38）、6（2=59.93）、7（2=41.88）、8（2=176.60）、9（2=12.28）、10（2=141.94）、11（2=49.33）、12（2=5.62）恩诺沙星、西林代谢物、呋喃唑酮代谢物淡水鱼2（2=42.64）、4（2=12.69）、5（2=16.26）、6（2=20.86）、7（2=29.05）、8（2=16.88）、9（2=5.83）、10（2=60.55）、11（2=22.88）恩诺沙星、孔雀石绿、五氯酚酸钠、磺胺类（总量）、甲氧苄啶、呋喃唑酮代谢物、地西泮和氯霉素海水虾8（2=82.11）、9（22=13.95）、10（2=24.64）餐馆（2=55.49）恩诺沙星、镉、呋喃它酮代谢物、氯霉素和呋喃唑酮代谢物海水蟹5（2=14.40）、10（2=144.24）镉生干籽类4（2=10.85）、11（2=6.89）、12（2=12.14）批发市场（2=6.88）酸价、黄曲霉毒素 B注：2为卡方值，表内为 P 0.05 的月份与场所。44 食品安全导刊 2023年10月（下）食安管理不当操作都有可能使花生发生霉变或者腐败，造成花生内部油脂的酸败和黄曲霉毒素的产生与富集5。综上所述，13 个品种中农兽药残留和重金属污染问题较为普遍，应结合农兽药残留与品种的关联性以及农业农村部禁限用公告，对重点项目与品种进行重点监控。重金属污染的主要原因是环境污染和食品品种本身特性，应重点关注产生问题品种的种植、养殖环节，对整体过程加强监管。4 结语通过对 20202022 年食用农产品监督抽检的分析，将食用农产品分为 4 个风险等级，监管部门可以根据不同风险等级的具体情况分配任务。对于其他水产品（牛蛙）、豇豆等多维度分析都属于高风险的品种，应加大抽检与监管力度，如开展专项行动等。对于部分如葱、茭白等不同维度分析大部分属于中高风险情况但不同分析方法风险情况有差异，应适当提高抽检批次进行风险验证，并结合各品种在不同月份与环节场所不合格率显著性差异情况，合理执行抽检计划。根据不同的食品品种对农兽药残留、重金属污染等问题的不同敏感程度以及种养殖环境影响，应采取不同措施进行干预和解决，结合农兽药的禁限用要求，针对不合格率较高的品种、时间、区域分布，合理安排检查和抽检工作，避免潜在的食品安全风险，提高监管的效率和精确性。对重金属污染严重的种植、养殖地进行溯源分析。由于食用农产品的特殊性，需要各部门间加强合作与信息共享，农业、环境等部门进行联动监管，加强对农业生产者的管理和监督，确保农药合理使用和安全使用，包括正确选择以及使用农药、控制施药频次、遵守安全间隔期和对重金属污染源的分析治理等，同时加大对农产品的质量检测和监督力度，建立完善的农产品溯源制度，保障公众的食品安全。本文主要以 20202022 年食用农产品安全监督抽检数据为例进行了分析与预警，剖析品种、项目、时间等方面的风险因素，为高风险种类的风险画像，期望通过信息化手段，为后期的监管与抽检工作提供方向，为政府监管部门提供决策支持，有针对性地制定监管措施，最大限度地提高监管效率和精确性，确保食品安全。参考文献1 宋晓婉,宋旭凤,姚凯.豆芽中植物生长调节剂的残留及风险监测评估 J.现代农业科技,2020(6):229-230.2 董金磊.氯吡脲对猕猴桃理化参数的影响及猕猴桃膨大果识别方法研究 J.咸阳:西北农林科技大学,2016.3 王定军,胡容,曾强,等.余庆县辣椒主要病虫害发生种类及其防治技术 J.农技服务,2023,40(9):67-70.4 宋成成.水产品中重金属污染现状及检测技术 J.中国食品工业,2022(18):40-41.5 郑冬梅,姬婧,张文山,等.日照市花生黄曲霉毒素污染成因与防控技术J.现代农业科技,2022(10):181-183.效落实，全面保障食品安全。参考文献1 郝璟弢.食品检验检测过程中存在的问题及应对措施 J.现代食品,2023,29(6):193-195.2 韦双.浅议食品快检技术在食品安全监管中的应用 J.中国食品工业,2023(5):67-69.3 赵翠芳.完善我国食品检验检测体系的思考 J.中国食品工业,2022(21):62-65.4 郭国栋.流通环节食品检测问题及处理方法 J.中国食品,2023(12):73-75.5 韦双.浅议食品快检技术在食品安全监管中的应用 J.中国食品工业,2023(5):67-69.（上接第 39 页）