目录一 问题提出...2 二 问题分析...3 三 基本假设...4 四 符号说明...4 五 数据收集分析 蔬菜数据收集及预处理 蔬菜各季节消费种类及比例分析 四季蔬菜市场分析 深圳市蔬菜输入渠道比例 安全

Transcription

1 深圳蔬菜安全风险评估与抽样检验方案研究 摘要 : 民以食为天, 蔬菜食品安全检测关系到每一个人的饮食健康, 是民生大事 本文根据题目所给的信息, 收集了相关数据, 并做了描述性统计分析, 然后建立的蔬菜风险评估模型, 并设计了合理的蔬菜抽检方案 首先, 根据近两年抽检蔬菜的种类和数量分析 得到深圳市各个季节的全部消费蔬菜均有 种 统计分析显示芸薹类 叶菜类 瓜类蔬菜合占了市场蔬菜总量的 65-80%, 菜心 白菜 上海青 番茄均为各季节的主流消费蔬菜 深圳市蔬菜输入渠道中农贸市场 批发市场 超市分别为 35% 32% 26%, 占了总体蔬菜市场的 93%, 本文还对主要蔬菜类的输入渠道进行统计分析, 得到各类蔬菜的输入渠道分布比例 然后, 分析抽检忽略的蔬菜输入渠道, 包括网购蔬菜 流动摊贩 城中村蔬菜消费分析等 得到通过网购渠道购买蔬菜对深圳市蔬菜安全风险影响为 1%, 但生鲜电商平台方兴未艾, 需要严格对网购的蔬菜进行安全检查 流动摊贩的蔬菜多为本地散户种植的蔬菜, 抽检难度大, 占深圳市蔬菜比例约为 2%, 给深圳市蔬菜安全带来较大的风险影响 其次, 针对深圳市蔬菜安全风险评估问题, 本文建立 IFS 蔬菜质量安全指数模型 模糊风险评估理论模型 农药危险系数模型 基于蔬菜抽检结果,IFS 模型求解得到深圳市平均安全指数 IFS = , 模糊理论模型求解得到蔬菜质量安全风险级别为中度风险, 其安全风险来源于氟虫腈 毒死蜱 敌敌畏, 氧乐果等 13 种超标农药, 其中氟虫腈的 IFS 值达 求解农药风险系数模型得到 13 种超标农药风险级别, 敌敌畏等 5 种农药属于低度风险, 氧乐果等 7 种农药属于中度风险, 毒死蜱为高度风险 最后, 根据深圳市历史抽检数据, 运用随机抽样估量模型估计抽样检验方案的样本数为 批次, 结合深圳市各季度蔬菜合格率对各季度进行不等概率抽检, 再根据深圳市消费蔬菜数量分布 输入渠道等指标, 设计详细的蔬菜抽检方案, 方案中将农药残留分为重点农药检测项目和常规农药检测项目, 还增加了重金属污染 食品添加剂项目的监测 基于对蔬菜抽检数据的多维度统计分析及模型求解的结果, 本文从各辖区蔬菜抽检的比例 各类蔬菜抽检权重 蔬菜输入渠道 农药知识普及的方面提出蔬菜质量安全监测的建议 关键词 : 蔬菜安全检测 IFS 模型 模糊数学 抽样检验 ~ 1 ~

2 目录一 问题提出...2 二 问题分析...3 三 基本假设...4 四 符号说明...4 五 数据收集分析 蔬菜数据收集及预处理 蔬菜各季节消费种类及比例分析 四季蔬菜市场分析 深圳市蔬菜输入渠道比例 安全监测忽略的蔬菜输入渠道 深圳市蔬菜监测结果分析 年全国抽检蔬菜质量检测结果分析 六 模型建立与求解 卡方检验 基于 IFS 的深圳市蔬菜质量安全风险评估模型 基于模糊理论的风险评价模型 农药风险系数模型 蔬菜抽样检验方案设计 深圳市新鲜蔬菜抽检具体方案 多目标优化模型 七 模型评价与推广 模型的优点 模型的缺点 模型的推广方向 八 推动深圳市蔬菜食品安全质量提高的建议信 九 参考文献 十 附录 ~ 2 ~

3 一 问题提出 民以食为天 食品安全问题是国家和民族安全问题 深圳市目前新鲜蔬菜消费属性是输入型的 除了批发市场 农贸市场 超市 配送中心等常规渠道外, 还有网络购物 流动摊贩 大单位 ( 企业饭店 ) 的点对点直接配送等输入方式 蔬菜品种繁多, 即使同一品名, 也存在产地的差别 ; 即使同一地域来源的同一品名蔬菜, 不同生产管理模式也会存在较大的质量差异 ; 此外还有反季节生产方式等 ; 这诸多因素都对传统的蔬菜安全监管工作方法构成新的挑战 根据附件中蔬菜抽检的流程及检测的数据, 完成以下任务 (1) 收集数据, 调查, 分析深圳市各季节新鲜蔬菜消费种类及不同输入渠道的比例 ; (2) 评估检验中忽略某些输入渠道的抽样导致的食品安全风险 ; (3) 基于附件评估深圳新鲜蔬菜食品安全风险 ; (4) 根据科学的统计推断原理, 设计一套深圳市新鲜蔬菜抽样检验方案, 以满足不同层面监测工作的需要, 提高食品安全防控效率与质量 ; (5) 基于研究结果给深圳市政府写一封建议信 ; 二 问题分析 在蔬菜种植过程中, 需要考虑农产品的收益, 因此会使用一些杀虫剂来防止害虫灾害, 在销售过程中会对一些蔬菜使用防腐添加剂以延长蔬菜保质期, 或者蔬菜受到周围空气的农药污染 而这些措施都会使蔬菜含一些有对人体有危害作用的化学成分 初步问题是在统计分析方面, 基于题目给出的附件数据, 本文结合深圳市市场和质量监督管理委员会官网 ( 上的数据, 收集了 2014 年 年各期的农产品质量安全报告 ( 详细数据请看附件 1) 年全国生鲜蔬菜检测结果 根据深圳市食用农产品质量安全检测抽检方案可知 : 每月的蔬菜质量安全抽检均选择当月的主销产品和时令产品进行安全检测 本文将季节按月份可划分为 : 春季 :3 月 5 月, 夏季 :6 月 -8 月, 秋季 :9 月 -11 月, 冬季 :12 月 -2 月, 并以此分析各季节新鲜蔬菜的消费种类 各蔬菜所占市场比例及输入渠道 接下来的重要问题是寻找蔬菜安全风险指标, 建立蔬菜食品安全风险评估体系, 本文认为蔬菜的安全风险指标包括蔬菜抽检的合格率 不合格食品有害物质的超标率 蔬菜产地 季节性等 在建立食品安全风险评估体系方面, 相关领域学者提出了基于 Delphi Method 和 AHP 的区域蔬菜质量风险评价方法 基于 IFS 法的蔬菜质量风险评价方法等 本文认为需要综合考虑实际数据特征及农药风险系数才能建立较好的安全风险评估模型 最后的问题是需要设计深圳市新鲜蔬菜检验方案, 通过抽验样本精确地推断总体样本情况的关键是运用科学合理的抽验方法 本文认为难点在于需要找出蔬菜食品安全监测的全部核心指标并协调分配到各层面监测工作中, 既要做到全面切实可行的安全检测也需避免重复检测 并根据在研究的结果, 从蔬菜食品安全检 普及农药残留基本知识 生产业农药施放管理制度方面予以参考性建议 ~ 3 ~

4 三 基本假设 (1) 抽检数据真实可靠 (2) 深圳市蔬菜市场供求关系处于平衡状态 (3) 蔬菜中农药残留量服从正态分布 (4) 蔬菜中农药残留的法定限制范围不发生更改 (5) 蔬菜检验过程不会出现技术性失误 (6) 短期内蔬菜产地 生产模式等不发生变化 四 符号说明 符号说明符号说明 P c 农药 c 的超标率 IFS c 各种农药人体健康的综合影响 F c 农药 c 的施检频率 M c 农药 c 的检出值 S c 农药 c 的敏感因子 N 蔬菜抽检年样本总量 q n 种不同蔬菜的购买指数向量 T 深圳市蔬菜年输入总量 IFS c 农药 c 对人体健康的影响 uα 2 标准正态分布面积之和等于 a 时的标准值 EDI c 农药 c 的实际摄入量估算值 p 以往蔬菜抽查的不合格率 f 安全摄入量矫正因子 δ 总体标准差 SIc 安全摄入量 B i 第 i 种蔬菜的合格率 m b 人体平均质量 pe k R c 蔬菜中农药 c 的最大检出值 nx ik 深圳市第 i 个辖区的人口比例深圳 k 月份第 i 各辖区应抽检的蔬菜数量 ~ 4 ~

5 五 数据收集分析 5.1 蔬菜数据收集及预处理 (1) 深圳市蔬菜抽检数据的采集本文结合深圳市市场和质量监督管理委员会官网 ( 上的数据, 收集了题目附件数据 深圳市 2014 年 年各个季度的农产品质量安全报告 深圳市 2016 年 年每月的农产品质量安全报告以及各期的蔬菜质量安全监测结果, 数据类型如表 1 表 2 所示 ( 详细数据请看附件 ) 表 年 1 月份全市食用农产品监测总体情况 区域类型抽样数量 不合格样品数量 合格率 (%) 总合格率 (%) 种植业产品 盐田畜禽产品 区水产品 种植业产品 大鹏畜禽产品 新区水产品 排名 区域 表 年 1 月份各区种植业产品监测情况汇总表 抽样数量 ( 份 ) 不合格样品数量合格率 (%) 排名 光明新区 盐田区 福田区 坪山区 罗湖区 龙华区 大鹏新区 龙岗区 南山区 宝安区 合计 (2) 全国蔬菜抽检数据的采集本文结合附件 3 提供的数据, 并在食品资讯中心 ( 收集了 2015 年 年全国抽检蔬菜制品检验结果, 如表 3 所示 ( 详细数据见附件 ) ~ 5 ~

6 表 3 蔬菜及监督抽检产品合格信息 序号 标识生产企业名称上海乐购物流有限公司闵行分公司 上海乐购物流有限公司闵行分公司 上海乐购物流有限公司闵行分公司 标识生产企业地址 上海市闵行区都会路 2385 号 上海市闵行区都会路 2385 号 上海市闵行区都会路 2385 号 被抽样单位名特易购商业 ( 上海 ) 有特易购商业 ( 上海 ) 有特易购商业 ( 上海 ) 有 被抽样单位所 上海 上海 上海 食品名称 小白菜 卷心菜 黄瓜 规格型号 500g/ 袋 500g/ 袋 500g/ 袋 生产日期 / 批 ( 进货 ) ( 进货 ) ( 进货 ) 分类 食用农产品 食用农产品 食用农产品 公告号 2016 第 23 期 2016 第 23 期 2016 第 23 期 公告日期 任务来源 / 项 上海 / 省抽 上海 / 省抽 上海 / 省抽 (3) 蔬菜名称规范化根据收集的数据, 发现在蔬菜抽检记录数据中有些是同一种蔬菜, 但记录的别名不同, 例如 : 结球甘蓝 (Brassica oleracea L.var.capitata L.) 是十字花科植物, 别名有 : 卷心菜 洋白菜 疙瘩白 包菜等 因此, 本文根据 中华人民共和国国家标准 : 蔬菜名称 1(GB ) 及 2017 年深圳市食用农产品质量安全监测抽检方案 对抽检蔬菜进行分类, 制定了蔬菜分类标准化表, 分类标准如表 4 所示 ( 详细表格请看附录 4), 并以此对抽检记录的蔬菜名称进行规范整理 叶菜类 菠菜莴苣芹菜蕹菜苋菜叶菾菜菊苣冬寒菜落葵 表 4 叶菜类包含的蔬菜品种 菠稜菜 赤根菜 角菜 波斯草千金菜 莴笋 生菜 青笋 莴苣笋 莴菜 油麦菜 莜麦菜芹 药芹 苦堇 堇葵 堇菜 旱芹空心菜 竹叶菜 通菜 藤菜 蓊菜米苋 赤苋 刺苋 青香苋 苋莙荙菜 牛皮菜 厚皮菜 光菜 叶甜菜欧洲菊苣 苞菜 吉康菜 法国莒荬菜冬苋菜 冬葵 葵菜 滑肠菜木耳菜 软浆叶 软姜子 染浆叶 胭脂豆 豆腐菜 藤菜 紫果菜 (4) 各月份适合种植的蔬菜类型蔬菜按气候类型可分为喜热型 喜寒型 耐寒型, 按成长类型又可分为长得快型 长得慢型等, 因此不同蔬菜的种植环境有差异 本文收集了蔬菜种植时间表, 记录了每个月适宜种植的蔬菜种类, 如表 5 为 1 月份适合种植的蔬菜 ( 所有月份数据请见附录 5) ~ 6 ~

7 表 5 1 月份适合种植的蔬菜类型 月份适合播种的蔬菜 1 月油菜 四月曼 菠菜 芥蓝 生菜 马铃薯 葱 茄子 番茄 辣椒 芋头 芫荽 5.2 蔬菜各季节消费种类及比例分析根据 深圳市食用农产品质量安全检测抽检方案 可知 : 每月的蔬菜质量安全抽检均选择当月的主销产品和时令产品进行安全检测 因此, 本文将季节按月份可划分为 : 春季 :3 月 5 月, 夏季 :6 月 -8 月, 秋季 :9 月 -11 月, 冬季 : 12 月 -2 月, 并根据 2016 年 2017 年期间的各期蔬菜质量安全检测的数据分析深圳市各季节新鲜蔬菜的消费种类 比例及输入渠道等信息, 所得结果如下 : (1) 春季 (3 月 -5 月 ) 表 6 深圳市春季所有销售的新鲜蔬菜及其比例 蔬菜名称所占比例 蔬菜名称 所占比例 蔬菜名称 所占比例 菜心 土豆 包芥菜 生菜 萝卜 香麦菜 青瓜 番茄 茼蒿 西红柿 胡萝卜 平包菜 上海青 白菜心 甜豆 结球甘蓝 红萝卜 蚕菜 黄瓜 荬菜 毛瓜 苦瓜 西葫芦 白萝卜 小白菜 洋葱 奶白菜 四季豆 菠菜 韭黄 油麦菜 麦菜 甘蓝 芥蓝 莴笋 南瓜 大白菜 花椰菜 西兰花 春菜 圆椒 大葱 豇豆 皇帝菜 黄皮椒 茄子 尖椒 红菜心 莲藕 本地白 苋菜 西洋菜 姜 苦麦菜 荷兰豆 红菜苔 芦笋 娃娃菜 通菜 蒜苗 青椒 番薯叶 芥菜 芹菜 白苋菜 丝瓜 ~ 7 ~

8 水生蔬菜类 4% 茄果类 13% 豆类 7% 瓜类 14% 其它蔬菜 3% 叶菜类 20% 芸薹类 32% 根菜类 4% 薯芋类 3% 图表标题 根菜类薯芋类芸薹类叶菜类瓜类茄果类豆类水生蔬菜类其它蔬菜 图 7 深圳市春季新鲜蔬菜种类及其比例 数据分析 -- 春季蔬菜 : 由表 6 和图 7 可知, 春季深圳市约有 66 种可消费的新鲜蔬菜, 其中菜心 青瓜 生菜等蔬菜销售比例较大, 分别占蔬菜市场的 7.0% 4.3% 4.8% 的比例 而在各大类蔬菜中, 芸薹类 叶菜类和瓜类的蔬菜消费的比例最大, 分别为 32% 20% 14%, 占总体蔬菜市场 66% 的比例, 其原因在于芸薹类 叶菜类和瓜类蔬菜品种比较丰富 有营养价值且易于种植, 是主流的消费蔬菜食品, 因此供应也比较多 (2) 夏季 (6 月 -8 月 ) 表 8 深圳市夏季所有销售的新鲜蔬菜及其比例 蔬菜名称所占比例蔬菜名称所占比例蔬菜名称所占比例 菜心 蕹菜 红薯叶 青瓜 西葫芦 麦菜 西红柿 结球甘蓝 红苋菜 生菜 通菜 秋葵 豇豆 白萝卜 花椰菜 苦瓜 菠菜 蒜苔 大白菜 芹菜 洋葱 芥蓝 娃娃菜 黄叶白 上海青 韭菜 皇帝菜 春菜 萝卜 潺菜 小白菜 荷兰豆 五彩椒 番茄 红萝卜 奶白菜心 奶白菜 尖椒 包生菜 芥菜 本地白菜 豇豆 ( 红 ) 荬菜 西兰花 莲藕 茄子 番薯叶 蒲瓜 油麦菜 苋菜 南瓜 黄瓜 白苋菜 丝瓜 圆椒 银丝王 冬瓜 四季豆 苦麦菜 土豆 ~ 8 ~

9 瓜类 20% 茄果类 13% 根菜类 4% 薯芋类 1% 葱蒜类 2% 芸薹类 27% 豆类 10% 叶菜类 23% 水生蔬菜类 0% 根菜类薯芋类葱蒜类芸薹类豆类水生蔬菜类叶菜类瓜类茄果类 图 9 深圳市夏季新鲜蔬菜种类及其比例 数据分析 -- 夏季蔬菜 : 由表 8 和图 9 可知, 夏季深圳市约有 60 种可消费的新鲜蔬菜, 其中菜心 青瓜 西红柿等蔬菜销售比例较大, 分别占蔬菜市场的 9.4% % 的比例 而在各大类蔬菜中, 芸薹类 叶菜类和瓜类的蔬菜消费的比例最大, 分别为 27% 23% 20%, 占总体蔬菜市场 70% 的比例 (3) 秋季 (9 月 -11 月 ) 表 10 深圳市秋季所有销售的新鲜蔬菜及其比例 蔬菜名称所占比例 蔬菜名称 所占比例 蔬菜名称 所占比例 菜心 结球甘蓝 苋菜 青瓜 番薯叶 小瓜 生菜 荷兰豆 银丝王 上海青 萝卜 大葱 西红柿 芹菜 红菜苔 大白菜 娃娃菜 红菜心 苦瓜 西芹 红豇豆 小白菜 西洋菜 红辣椒 春菜 包心芥菜 淮山 豇豆 白菜 黄皮椒 芥蓝 空心菜 尖椒 油麦菜 包生菜 韭菜 奶白菜 红苋菜 苦麦菜 荬菜 黄白菜 莲藕 茄子 土豆 毛瓜 麦菜 莴笋 南瓜 菠菜 茼蒿 生姜 黄瓜 水白菜 水瓜 芥菜 蕹菜 笋 白萝卜 长白菜 西兰花 ~ 9 ~

10 蒲瓜 白菜苔 香麦菜 圆椒 白菜心 小南瓜 四季豆 胡萝卜 雪里红 西葫芦 花椰菜 洋葱 红萝卜 青椒 长豆角 茄果类 10% 豆类 5% 水生蔬菜类 1% 根菜类 4% 薯芋类 2% 瓜类 17% 芸薹类 41% 叶菜类 20% 图表标题 根菜类薯芋类葱蒜类芸薹类叶菜类瓜类茄果类豆类水生蔬菜类 图 11 深圳市秋季新鲜蔬菜种类及其比例 数据分析 --- 秋季蔬菜 : 由表 10 和图 11 可知, 秋季深圳市约有 75 种可消费的新鲜蔬菜, 其中菜心 青瓜 生菜等蔬菜销售比例较大, 分别占蔬菜市场的 8.5% 8.3% 6.7% 的比例 而在各大类蔬菜中, 芸薹类 叶菜类和瓜类的蔬菜消费的比例最大, 分别为 41% 20% 17%, 占总体蔬菜市场 78% 的比例 (4) 冬季 (12 月 -2 月 ) 表 12 深圳市冬季所有销售的新鲜蔬菜及其比例 样品名称所占比例样品名称所占比例样品名称所占比例 菜心 土豆 葱 青瓜 荬菜 番薯叶 生菜 胡萝卜 佛手瓜 上海青 莲藕 皇帝菜 小白菜 西葫芦 青椒 油麦菜 白菜 蒜苔 西红柿 红萝卜 豌豆 大白菜 茼蒿 小南瓜 豇豆 圆椒 长白菜 苦瓜 奶白菜 包生菜 菠菜 芹菜 春菜心 春菜 红薯叶 大葱 番茄 韭菜 红苋菜 白萝卜 莴笋 黄皮椒 结球甘蓝 尖椒 黄芽白菜 麦菜 姜 空心菜 ~ 10 ~

11 黄瓜 萝卜 毛瓜 四季豆 南瓜 泡椒 芥蓝 西兰花 沙葛 娃娃菜 包心芥菜 山药 西洋菜 红菜苔 生姜 茄子 蒲瓜 水白菜 荷兰豆 通心菜 水瓜 红菜心 白菜苔 蕹菜 芥菜 白菜心 西生菜 茄果类 9% 豆类 8% 图表标题 根菜类 5% 瓜类 16% 芸薹类 35% 水生蔬菜 1% 叶菜类 23% 葱蒜类 1% 薯芋类 2% 根菜类芸薹类薯芋类葱蒜类叶菜类水生蔬菜瓜类茄果类豆类 图 13 深圳市冬季新鲜蔬菜种类及其比例 数据分析 - 冬季蔬菜 : 由表 12 和图 13 可知, 冬季深圳市约有 75 种可消费的新鲜蔬菜, 其中菜心 青瓜 生菜等蔬菜销售比例较大, 分别占蔬菜市场的 7.5% 6.8% 6.3% 的比例 而在各大类蔬菜中, 芸薹类类 叶菜类和瓜果类的蔬菜消费的比例最大, 分别为 35% 23% 16%, 占总体蔬菜市场的 71%, 这也反映出市民购买蔬菜更偏向于白菜类 叶菜类和瓜类的蔬菜 5.3 四季蔬菜市场分析本文统计发现, 一些主流蔬菜例如青瓜 菜心 菠菜 西红柿等, 且这些蔬菜所占的市场比例最大, 是每个季节主销食品, 其原因在于这些类型蔬菜产量大 能够四季种植 蔬菜口感较好 5.4 深圳市蔬菜输入渠道比例 (1) 深圳市蔬菜输入渠道分析 : 根据 2017 深圳市食用农产品质量安全检测抽检方案 可知, 目前常用的蔬菜渠道包括 : 批发市场 农贸市场 超市 蔬菜生产基地 本文根据 2014 年 2017 年各期的 食用农产品质量安全报告 的数据统计, 得到深圳市蔬菜输入渠道比例 ~ 11 ~

12 批发市场 32% 生产基地 3% 商场超市 30% 农贸市场 35% 商场超市农贸市场批发市场生产基地图 14 深圳市蔬菜输入渠道比例 蔬菜输入渠道结果分析 : 根据图 14 可知, 现阶段深圳市蔬菜输入渠道中, 农贸市场 批发市场 商场超市 生产基地占得比例分别为 35% 32% 30%,3% (2) 各大类蔬菜的输入渠道分析根据蔬菜消费种类的分析发现芸薹类 叶菜类 瓜类在四季当中所占的比重都很大, 因此本文根据数据分别分析这三类蔬菜的输入渠道比例情况 芸薹类输入渠道比例 : 如图 15 所示为各渠道输入叶菜类蔬菜的比例 深圳市芸薹菜类的蔬菜输入渠道主要包括农贸市场 批发市场 超市 生产基地和商场超市, 其中农贸市场占了总体蔬菜市场的 42%, 如图 15 所示为各渠道输入芸薹类蔬菜的比例 商场超市批发市场 8% 22% 生产基地 3% 超市 25% 农贸市场 42% 超市农贸市场批发市场商场超市生产基地 图 15 深圳市芸薹菜类输入渠道比例 叶菜类输入渠道比例 : 如图 16 所示为各渠道输入叶菜类蔬菜的比例 深圳市叶菜类的蔬菜输入渠道主要包括农贸市场 批发市场 超市 生产基地和商场超市, 其中农贸市场占了总体蔬菜的 40%, 原因在于叶菜类蔬菜存放周期较短, 且为普通市民主要消费的蔬菜 商场超市 8% 生产基地 8% 超市 21% 批发市场 23% 农贸市场 40% 超市农贸市场批发市场商场超市生产基地 图 16 深圳市叶菜类输入渠道比例 ~ 12 ~

13 瓜类输入渠道比例 : 深圳市瓜类的蔬菜输入渠道主要包括农贸市场 批发市场 超市 生产基地和商场超市, 其中商场市场输入渠道的占了总体蔬菜 39%, 如图 17 所示为各渠道输入瓜类蔬菜的比例 商场超市 39% 生产基地 1% 批发市场 19% 超市 16% 农贸市场 25% 超市农贸市场批发市场商场超市生产基地 图 17 深圳市瓜类输入渠道比例 5.5 安全监测忽略的蔬菜输入渠道目前深圳市蔬菜质量安全检查场所包括各区所有的批发市场 农贸市场 超市 蔬菜生产基地, 覆盖了深圳市大多数蔬菜的输入渠道, 但缺少了对网购蔬菜 流动摊贩 大单位点对点运输蔬菜的监测 (1) 网购蔬菜分析随着市民消费支出不断提高, 对蔬菜食品品质 安全 便捷性方面提出了更高的需求, 因此有的消费者会直接在互联网上购买蔬菜 本文根据 2016 年上半年深圳居民购物行为调查报告 的数据统计 2, 得到深圳市居民日常购买生鲜类产品的主要渠道比例分布, 如图 18 所示 根据 2016 年中国生鲜电商市场现状分析及发展趋势预测 的数据统计 3, 得到用户网购生鲜类产品的比例分布 中国生鲜电商市场发展趋势, 如表 19 图 20 所示 网购, 4% 线下, 96% 网购线下图 18 深圳市居民日常购买生鲜产品主要渠道来源 : 2016 年上半年深圳居民购物行为调查报告 表 年网购用户网购的生鲜类产品分布 排名 产品 排名 网购品类分布 1 蔬菜类 1 水果类 72.0% 2 水果类 2 乳品 49.0% 3 肉类 3 肉类 46.7% ~ 13 ~

14 4 水产类 4 水产 4.63% 5 乳品 5 蛋类 24.9% 6 蛋类 6 蔬菜类 23.7% 令 Q sz 为深圳市民网购生鲜产品的比例,U sz 为 2015 年网购生鲜类产品的用户中网蔬菜类的比例, P sz 为深圳市新鲜蔬菜中网购输入渠道占的比例, 则可得到 : P sz = Q sz U sz Q sz = 4%,U sz = 23.7% 解得 :P sz = 0.948% 图 20 中国生鲜电商市场交易规模及增长率来源 : 2016 年中国生鲜电商市 场现状分析及发展趋势预测 结果分析 : 在深圳市总体蔬菜市场中, 通过网络渠道购买蔬菜的比例约为 0.95% 即忽略网购渠道的蔬菜抽样导致的食品安全风险占深圳市食品安全风险的比例约为 1%, 但未来生鲜电商交易预计会以 53%~80% 的年增长率迅速升高, 而网购的蔬菜来源于全国各地的蔬菜电商, 食品监督局需要与生鲜电商平台进行协作, 严格对网购蔬菜的质量安全进行监测 (2) 流动摊贩售卖蔬菜分析流动商贩的蔬菜大多来源于个体户种植, 属于深圳市本地种植蔬菜, 还有部分来源于批发市场 根据深圳市统计局发布的数据显示, 深圳现有农业用地约 10 万亩, 其中包括基本农田 3.23 万亩, 蔬菜年产量 S=10.18 万吨 设深圳市种植蔬菜占深圳市年消费蔬菜的比例为 P l, 深圳市年消费蔬菜总量为 W, 深圳市拥有的人口为 L 据 2016 年全国平均身高体重, 中国营养学会 2007 年中国居民平衡膳食宝塔建议, 综合计算得到 : 深圳市人均摄入的蔬菜量 F = 275g/( 人 天 ) 深圳市蔬菜年消费总量 W = F L 365 ~ 14 ~

15 深圳市本地种植蔬菜比例 P l = S/W 求解可得 P l = , 因此深圳市本地种植的蔬菜占总蔬菜市场的比例为 4.02%, 深圳蔬菜生产基地的蔬菜占的市场总量 P s = 3%, 因此计算深圳市流动摊贩销售的蔬菜占比 P t: P t = P l P s, 解得 :P t = 1.02% 流动摊贩 = 本地种植 + 常规渠道 结果分析 : 深圳市通过流动摊贩渠道消费的蔬菜仍存在安全风险问题, 而抽检流动摊贩的蔬菜难度大且样本代表性不强, 应加强对深圳市本地种植蔬菜的统一化管理, 从源头把控本地种植蔬菜的安全质量 (3) 大单位点对点蔬菜分分析当下抽检部门着重对蔬菜的重点输入渠道进行抽检, 而缺少对大型连锁餐饮单位点对点运输蔬菜的质量抽检, 可能导致蔬菜质量在餐饮这一环节出现安全问题 根据 深圳市各区餐饮服务食品安全量化 A 级单位名单通报, 就 2010 年 年每年深圳市各区被评为食品安全量 A 级的连锁餐饮单位数量情况, 从侧面反映大型连锁餐厅蔬菜质量安全情况 结果以图 21 展示如下 : 福田罗湖南山宝安龙岗 盐田光明坪山龙华大鹏 图 年 年各区 A 级连锁餐饮单位数量情况 结果分析 : 深圳市餐饮行业 A 级单位数量呈现在逐年上升趋势, 进而说明了蔬菜安全质量监管呈现良好态势 而盐田区 大鹏新区 坪山新区和光明新区被评为安全量 A 级的餐饮单位较少, 从侧面反映出此四个辖区餐饮食品质量有待提高, 蔬菜质量同样需要提高 根据 2016 年深圳市食品及食用农产品抽检情况报告 数据, 对餐饮行业蔬菜质量安全进行评估 ~ 15 ~

16 表 21-2 餐饮环节蔬菜制品合格率 食品类别食品品种批次 不合格批 次 不合格项目 ( 项目数 ) 合格率 蔬菜制品 酱腌菜 二氧化硫残留量 (10) 其他蔬菜制品 蔬菜干制品 二氧化硫残留量 (11) 食用菌制品 % 结果分析 : 餐饮环节蔬菜质量安全抽检中, 蔬菜合格率较高, 达 96.72% 而在不合格批次中, 酱腌菜 蔬菜干制品分别约占 50%, 且不合格原因均为二氧化硫残留超标, 警示餐饮单位需要严格把控添加剂使用量 另外,2016 年深圳市市内确诊了 11 起餐饮环节食物中毒事件, 其中因致病微生物造成事故占比为 63.64%( 图 21-3): 图 起食物中毒事件中各原因 结果分析 : 深圳市餐饮环节在蔬菜抽检合格率高的情况下, 致病微生物 添加剂等抽检项目占不合格 有问题蔬菜较大比例, 示意着餐饮单位仍存在蔬菜二次污染问题, 餐饮单位仍需致力于改善此种状况 (4) 城中村小型市场抽检渠道忽略分析城中村是深圳市居民主要生活区, 也是流动摊 - 贩高度集中的地方, 因此市场食品安全影响范围较广, 对当中食品的抽检便显得尤为必要 深圳市外来人口聚居地城中村扮演着重要角色, 本文分析得到深圳市外来人口居住在城中村房屋所占比例为 58.7%, 城中村是外来人口生活地集中地 ~ 16 ~

17 图 21-4 深圳市 2010 年外来人口空间聚集度分布图 外来人口分布 : 多数集中在罗湖区一带, 而罗湖区为深圳市旧城区, 当中城中村分布较多, 此结果可以进一步阐述上述观点 : 城中村中人口数量所占比多, 外来人口聚集现象突出 本文从两个角度分析城中村小型市场蔬菜安全风险 : (1) 流动摊贩平均通勤时间考虑到广州市与深圳市同为一线发达城市, 城中村现象也相似, 所以广州市市场的一些相关情况可以酌情在深圳市予以考虑 城中村小型市场的聚集通常都是流动摊贩的聚集形成, 而通过分析流动摊贩的平均通勤时间, 便可大致推断城 [4] 中村小型市场流动摊贩集中情况 在黄耿志等人基于广州市 200 个摊贩访谈的实证研究中发现, 得出摊贩从居住地到摆摊地的交通方式和通勤时间如下表所示 : 表 21-5 摊贩从居住地到摆摊地的交通方式和通勤时间表 结果分析 : 由以上分析可知平均通勤时间为 23.6 分钟, 可见, 距离因素致使流动摊贩更趋于集中在人口密度高 租金成本低的地区, 而城中村便符合这样的要求 事实证明, 城中村小型市场便是流动摊贩的密集集中场所 流动摊贩蔬菜抽检难度大, 质量安全难以得到保证, 而且其影响范围和人口较多, 故而风险不可忽视 (2) 居民经济收入水平外来人口收入不稳定, 故而居住情况受影响因素较多, 深圳市 2010 年外来 L 人口达到 798 万, 记为 wl P, 而租金选择因素所占比记为 x, 因经济收入较低人 ~ 17 ~

18 L 群数量为 d, 可得 : 解得 : L 798, P 25% L wl d L wl x * P 低收入人群较为关注日常开支, 而城中村里低收入人群趋向选择小型市场解决日常蔬菜需求, 而若忽略小型市场蔬菜质量抽检,25% 的外来人口可能受到影响, 而这个数字较大, 不可忽视 随着深圳市经济发展, 外来人口数量仍在增加, 目前深圳市外来人口已远超以上数字, 所以小型市场蔬菜质量抽检任务便显得尤为必要 (3) 居民蔬菜消费流通假设市内蔬菜从农田到消费者之间构成一个网络, 蔬菜最终全部进入消费者终端, 各单位蔬菜量处于动态平衡状态, 蔬菜供求关系平衡 ( 图 21-6) x 外源地输入 本地种植 农贸市场 批发市场生产基地散户种植 a b c 超市 d 城中村小型市场 消费者 表 21-6 深圳市蔬菜流通渠道简图 按照蔬菜从 农田到消费者 这条输送链分析可知, 可将此输送阶段大致分成三部分 : 第一部分是外源地输入和本地种植, 第二部分是农贸市场 批发市场 生产基地和散户种植, 第三部分是超市 城中村小型市场和消费者 上文分析已得输入渠道之比, 即为蔬菜流通量之比, 比例如下表 : 表 21-7 深圳市蔬菜输入渠道比例 流通节点生产基地超市批发市场农贸市场 流通比例 3% 30% 32% 35% 从输送链第一部分分析, 外源地输入占比 95%, 所以批发市场占比 =[32/ (32+35)]*95%=45%, 农贸市场占比 =95%-45%=50%; 设生产基地占全市蔬菜流通 量为 X s, 则 ~ 18 ~

19 35% 50% 3% X s 解得 : X 4.29% s 第二部分分析, 因为上表流通量之比是根据全市流通量占比得到的, 故求超市流通量在第二部分流通量占比时, 应予以还原 设蔬菜流经第二部分时超市蔬 菜量占第二部分蔬菜总量为 X c, 所以有 X c 1 X c 30% 求得 X 43% c 第三部分, 因为蔬菜流通量最后流通到消费者处, 故有以下线性方程组, 设 当蔬菜在第二部分流通时城中村蔬菜流通量占第二部分蔬菜总量为 X cz, 有 a b 1 c d 1 45%* a 43% X 1 cz 50%* c 4.29% 43% 0 a 1 0 b 1 0 c 1 0 d 1 解得 : X a c 77.42% d 22.58% cz 第三部分中, 消费者得到蔬菜方式大致可分为农贸市场 超市购买和城中村中小型市场购买 因农贸市场 大型超市多位于城中村外, 故分析深圳市城中村内外常住人口比大致可求出 a,b 两值 深圳市城中村居住约 700 万人口, 而 2016 年统计局数据显示, 深圳市常住人口达 1191 万, 则城中村外人口约 491 万, 以此 a b 解得 : a 41% b 59% 结果分析 : 城中村蔬菜流通量占第三部分蔬菜总量约 39%, 折算成全市占比约 24%, 若不经过抽检, 则风险较大 城中村小型市场是流动摊贩高密度集中地, 当中卫生情况较差, 蔬菜二次污染概率较大, 不合格率可能较蔬菜流通链第一部分低, 即风险性可能更高 所以需加强对城中村小型市场的摊位登记和抽检, 以降低风险 ~ 19 ~

20 5.6 深圳市蔬菜检测结果分析 (1) 深圳市蔬菜抽检总体合格率对深圳市 2017 年 1 月 -3 月食用农产品监测抽检结果进行分析, 如图 22 所示, 第一期蔬菜检测合格率为最低值 88.80%, 第十期的合格率达到最高 98.2%, 蔬菜质量安全监测合格率呈波动上涨的态势, 合格率的总体水平在 94% 以上 % 98.00% 96.00% 94.00% 92.00% 90.00% 88.00% 86.00% 84.00% 94.40% 94.00% 98.00% 97.70% 98.20% 96.40% 94.90% 94.30% 92.50% 88.80% 第一期 第二期 第三期 第四期 第五期 第六期 第七期 第八期 第九期 第十期 图 年 1 月 -3 月深圳各期蔬菜检测项目合格率 (2) 各抽检场所蔬菜的合格率本文选取了深圳市 2016 年 9 月 年 3 月的各期蔬菜食品抽检结果, 统计分析各大蔬菜场所的抽检合格情况, 如图 23 所示 : % 98.00% 96.00% 94.00% 92.00% 90.00% 88.00% 86.00% 84.00% 98.90% 94.90% 96.80% 95.70% 90% 生产基地 农贸市场 批发市场 商场超市 小超市 图 23 深圳市 2016 年 9 月 年 3 月各场所蔬菜抽检合格率 结果分析 : 根据图 23 可知, 生产基地蔬菜抽检的合格率最高, 达到 98.9%, 农贸市场 批发市场 商场超市的合格率也处于较高值, 达到 94.5% 以上 但小超市场所的合格率为 90%, 相对其他场所较低 本文分析认为, 生产基地的蔬菜受到运送过程中的二次污染较少, 因此合格率较高, 农贸市场 批发市场 商场超市有较为严格的市场监管制度, 而小超市的蔬菜因管理较为疏松且容易受周围环境污染或储存不规范导致合格率较低 ~ 20 ~

21 (3) 各期抽检不合格的蔬菜数量深圳市 1-10 期的蔬菜中, 检出农药超标的蔬菜数量为 43 批, 各期不合格蔬菜数量如图 24 所示 第一期第二期第三期第四期第五期第六期第七期第八期第九期第十期 图 年 1 月 -3 月深圳各期蔬菜检测不合格数量 (4) 超标农药比例分布深圳市 1-10 期抽检蔬菜需要检测的农药 58 种, 其中有检出 13 种不同农药超标, 主要包括敌敌畏 毒死蜱 氟虫腈等, 超标农药的检出量在 (mg/kg) 之间, 其中的高毒农药包括 : 氧化乐果 克百威 甲拌磷 甲胺磷 甲基异柳磷 5 种, 同时因毒死蜱超标的蔬菜占不合格蔬菜的 33%, 如图 25 所示 2% 2% 4% 8% 2%2% 6% 2% 4% 12% 33% 19% 4% 敌敌畏 啶虫脒 毒死蜱 多菌灵 氟虫腈 甲胺磷 甲拌磷 甲基异柳磷 克百威 氯氰菊酯 图 年 1 月 -3 月深圳市不合作蔬菜中的农药分布 (5) 农药的超标倍数统计各种农药的超标倍数在 之间, 如表 26 所示, 其中甲基异柳磷的含量是合格指标的 53 倍, 不合格蔬菜中农药的平均超标倍数达 8.51, 因此不合格蔬菜食品极易对身体健康造成损害, 需制订有关农药施放限制政策严格控制在农产品中使用农药的剂量 ~ 21 ~

22 超标倍数 检测项目 表 26 深圳市 1-10 期蔬菜检测中不合格蔬菜超标情况 超标检出平均值 (mg/kg) 法定标准 (mg/kg) 超标倍数 敌敌畏 啶虫脒 毒死蜱 多菌灵 氟虫腈 甲胺磷 甲拌磷 甲基异柳磷 百克威 氯氰菊酯 辛硫磷 溴氰菊酯 氧化乐果 图 年 1 月 -3 月深圳市不合格蔬菜中农药超标倍数 (6) 深圳市蔬菜合格率按季度统计根据深圳市食品药品监督局公布的 年的蔬菜抽检结果, 统计各个季度蔬菜检测的合格率, 如表 28 所示 表 28 深圳市 年各季度蔬菜抽检合格率 季度 年份 抽检总数量 4*360 4*360 12* *3 第一季度 (1-3) 91.7% 89.2% 91.7% 95.6% 95.4% 95.0% ~ 22 ~

23 94.0% 95.2% 95.6% 92.8% 98.1% -- 第二季度 (4-6) 97.1% 97.9% 93.9% 93.3% 98.8% -- 第三季度 (7-9) 97.7% 98.3% 88.6% 92.2% 96.3% -- 第四季度 (10-12) 96.9% 97.5% 一年抽检合格率 92.45% 91.9% 96.64% % 98.00% 96.00% 94.00% 92.00% 90.00% 88.00% 86.00% 84.00% 82.00% 图 深圳不同季度蔬菜中检测项目合格率 表 年深圳市各季度合格率平均值 季度 第一季度 第二季度 第三季度 第四季度 合格率 92.46% 95.35% 95.15% 92.57% 结果分析 : 由图表 (28-30) 可知, 深圳市 年各季度呈现 双峰梯形型 分布, 出现 2 个高峰,2014 年第二季度,2016 年第二季度 第四季度的蔬菜合格率有明显的升高 第一季度 第四季度蔬菜抽检合格率低于第二季度 第三季度的蔬菜抽检合格率, 本文分析认为第一季度南方地区多为潮湿雨水季节害虫增多, 第四季度蔬菜种植种类较少, 因此反季节蔬菜较多, 需要施放各种植物生长剂及灭虫剂, 导致蔬菜农药残留量增加 (7)2016 年深圳市各类蔬菜抽检合格率根据 2016 年深圳市蔬菜抽检结果, 统计分析各类蔬菜的不合格率, 如表 31 所示 表 年深圳市八大类蔬菜抽检合格率 月份 不同蔬菜类型的合格率 (%) 茄果类瓜类根茎类水生类 叶菜类 鳞茎类芸薹属类豆类 ~ 23 ~

24 平均值 结果分析 : 由表 31 可知,2016 年深圳市各类型蔬菜合格率都达到 91% 以上, 其中, 主流消费蔬菜叶菜类 瓜类 茄果类的全年平均合格率都达到了 95% 以上, 水生类蔬菜合格率达 100%, 原因在于水生蔬菜的生长环境水量多, 能够起去除蔬菜中多数农药的作用 年全国抽检蔬菜质量检测结果分析对题目附件 3 的数据进行分析, 提取出鲜菜类不合格信息情况表, 并根据链接, 收集了 年全国抽检蔬菜的合格信息 (1) 全国生鲜蔬菜抽检合格率根据附件 3 中的 310 批不合格生鲜蔬菜的抽检来源分析, 得到抽检合格的生鲜蔬菜数量为 7159, 近两年全国蔬菜合格率为 95.85%, 如图 32 所示 深圳市蔬菜属于输入型, 因此全国生鲜蔬菜抽检合格率, 可作为深圳市蔬菜抽检误差检验的参考标准 不合格蔬菜, 310, 4.15% 合格蔬菜, 7159, 95.85% 不合格蔬菜 合格蔬菜 图 32 全国蔬菜抽检合格情况 (2) 不合格生鲜蔬菜分析统计不合格新鲜蔬菜的类别, 如表 33-1 所示, 导致蔬菜食品不合格的主要包括添加剂 农药残留 重金属污染等 8 种类型, 其中因农药超标残留导致蔬菜不合格的比例占了 60%, 其次是食物添加剂占了 22%, 如图 33-2 所示, 出 ~ 24 ~

25 农药残留外的因素导致蔬菜不合格的比例占了 40%, 因此在蔬菜检测中, 也需要对重金属 添加剂进行的严格的监控才能保证蔬菜食品的安全性 表 33-1 全国蔬菜抽检不合格类别分布 不合格项目分类 数量 不合格项目分类 数量 标签 1 重金属污染 34 非法添加剂 14 质量指标 1 农药残留 186 微生物 4 生物毒素 1 食物添加剂 69 食物添加剂 22% 微生物 1% 标签 1% 非法添加剂 5% 质量指标 0% 重金属污染 11% 生物毒素 0% 农药残留 60% 标签非法添加剂农药残留生物毒素 图 33-2 蔬菜不合格类别的比例 (3) 不合格蔬菜类型统计根据近两年全国抽检不合格生鲜蔬菜的结果, 统计分析不合格蔬菜的类型分布, 如图 34 所示 水生蔬菜 1% 豆类 11% 长年生长的其它蔬菜 18% 根菜类 3% 薯芋类 8% 葱蒜类 10% 芸薹类 6% 茄果类 8% 瓜类 3% 叶菜类 32% 图 34 全国抽检不合格蔬菜类型比例 结果分析 : 根据图 34 可知, 在不合格蔬菜中, 市场主销蔬菜叶菜类的比例达到 32%, 而同为主流蔬菜瓜类占的比例为 3%, 原因是叶菜类的蔬菜品种最多, 消费需求高, 而且叶菜类蔬菜容易受到虫害, 使得农药施放比例较高 因此需要增加叶菜类蔬菜抽检的批次, 提高抽检要求 (4) 农药超标类型统计 ~ 25 ~

26 苯醚甲环唑 代森锰锌 滴滴涕 敌敌畏 啶虫脒 毒死蜱 多菌灵 氟虫腈 腐霉利 高效氯氟氰菊酯 甲氨基阿维菌素苯甲酸盐 甲拌磷 甲基异柳磷 克百威 乐果 硫丹 氯氟氰菊酯 氯氰菊酯 灭蝇胺 水胺硫磷 涕灭威 辛硫磷 氧乐果 乙酰甲胺磷 甲胺磷 根据近两年全国抽检不合格生鲜蔬菜的结果, 统计分析全国抽检不合格蔬菜农药残留类型及其数量, 如图 35 所示 图 35 全国抽检不合格蔬菜农药残留类型 结果分析 : 根据图 35 可知, 在不合格蔬菜中包括了 25 种不同类型的农药, 包括氧化乐果 克百威 甲拌磷 甲胺磷 甲基异柳磷等高危或不得检出的农药 其中毒死蜱 克百威 氧乐果检出数量分别为 , 明显高于其它农药检出数量 (5) 不合格的蔬菜中农药的超标倍数农药残留量是衡量蔬菜安全风险最重要的指标之一, 统计分析近两年全国抽检不合格蔬菜农药的超标倍数, 如表 36 图 37 所示 根据统计结果可知, 全国蔬菜抽检农药超标平均倍数为 14.37, 其中甲拌磷超标倍数达 139, 为所有农药残留超标倍数的最高值, 不可检出的超标农药中硫丹 水胺硫磷含量分别为 0.122mg/kg 0.535mg/kg, 因此不合格蔬菜食品仍会引发较大的安全问题, 严格限制甲拌磷 氧乐果 水胺硫磷等超标严重的农药的施放剂量 农药类型 平均值 (mg/kg) 图 36 全国抽检不合格蔬菜农药超标倍数统计表 法定限制 超标倍数 农药类型 抽检平均值 (mg/kg) 法定限制 ~ 26 ~ 超标倍数 苯醚甲环唑 甲基异柳磷 代森锰锌 克百威 滴滴涕 乐果 敌敌畏 硫丹 啶虫脒 氯氟氰菊酯 毒死蜱 氯氰菊酯

27 安徽 北京 福建 甘肃 广西 广东 贵州 海南 河北 河南 吉林 江苏 江西 内蒙古 青海 山东 上海 四川 天津 新疆 云南 重庆 多菌灵 灭蝇胺 氟虫腈 水胺硫磷 腐霉利 涕灭威 氯氟氰菊酯 辛硫磷 甲氨基 阿维菌素苯甲酸盐 氧乐果 甲胺磷 乙酰甲胺磷 甲拌磷 平均超 标倍数 图 37 全国抽检不合格蔬菜农药超标倍数 (6) 不合格蔬菜的地区分布蔬菜种植会受到当地环境 气候的影响, 即同一种蔬菜统也会因不同地区种植产生质量安全的差异, 统计全国抽检不合格生鲜蔬菜的生产地区, 如表 38 所示 表 38 不合格生鲜蔬菜地区分布 ~ 27 ~

28 结果分析 : 根据表 38 可知, 全国蔬菜抽检中, 山东 上海 贵州地区生产的蔬菜不合格数量较高分别为 , 深圳市在选择蔬菜来源时, 应对同一蔬菜种类的不同地区加以分析, 选择合格率较高的生产地区的蔬菜 六 模型建立与求解 6.1 相关性分析 (1) 蔬菜抽检渠道的卡方检验表 40-1 深圳市蔬菜抽检统计表 批发市场农贸市场总计 合格 不合格 总计 根据表 40-1 的蔬菜统计数据, 本文计算得出卡方值, 取自由度 df = 1, 查找卡方检验临近值表, 解得 : K 2 = < K 2 < 因此, 有 95% 的把握认为在不同抽检场所与蔬菜合格情况有关, 本文分析认为批发市场与农贸市场的蔬菜管理模式不同, 蔬菜种类数量有差异, 从而对蔬菜抽检指标有影响 (2) 蔬菜抽检季度的卡方检验表 年深圳市各季度蔬菜抽检结果 第一季度 第二季度 总计 合格 不合格 总计 根据表 40-2 的蔬菜统计数据, 本文计算得出卡方值, 取自由度 df = 1, 查找卡方检验临近值表, 解得 : K 2 = K 2 > 因此, 有 99.9% 的把握认为蔬菜合格率与不同季度抽检有关, 本文分析认为第一季度 (1-3 月 ) 为冬季过渡到春季, 适宜自然种植蔬菜较少, 反季节蔬菜供应比例较多 而反季节蔬菜需要施放多种植物激素促进生长且大棚蔬菜温度 湿度较高, 不利于农药的分解, 导致农药残留升高 6.2 基于 IFS 的深圳市蔬菜质量安全风险评估模型基于蔬菜抽检不合格率的评价模型具有可靠性 直观性的特点, 但也存在许多不足之处, 例如没有结合蔬菜具体的食用部分 人体摄入量等指标, 使得评估风险的能力较低 为进一步评估蔬菜安全风险, 本文基于食品安全指数法 (IFS) 来评估深圳市新鲜蔬菜食品的安全风险 (1) 模型建立 ~ 28 ~

29 IFS c = EDI c f SI c mb IFS = n i=1 IFS c n EDI c = R F E P, SI c = ADI c 考虑忽略渠道的蔬菜产生风险, 对深圳市蔬菜安全风险进行如下计算 : R=IFS *(1+2%) (2)IFS 模型求解根据 2016 年全国平均身高体重 的报告, 全国 18 岁以上居民平均体重约为 61kg, 据中国营养学会 2007 年中国居民平衡膳食宝塔建议 显示, 体重为 60kg 的成人每天摄入的蔬菜量为 275g, 而芸薹类 叶菜类 瓜类蔬菜可食用部分可认为是 100%,R 取不合格蔬菜检测中农药最大检出值, 因此得到以下指标 : mb = 61kg, F = 278g/( 人 天 ) E = 1, f = 1, P = 1 R c = Max(M c ) 本文根据 GB 查询到蔬菜抽检项目中 58 种农药的 ADI 值, 如表 41 所示为农药相对应的 ADI 值 (58 种农药的 ADI 值请看附录 6) 表 41 检测项目中的 14 种农药的 ADI 值 农药名称 ADI/(mg kg 1 ) 农药名称 ADI/(mg kg 1 ) 敌敌畏 甲胺磷 啶虫脒 0.07 甲拌磷 毒死蜱 0.01 甲基异柳磷 多菌灵 0.03 百克威 氟虫腈 氯氰菊酯 0.02 氧化乐果 辛硫磷 溴氰菊酯 0.01 六六六 根据深圳市 1-10 期蔬菜抽检数据, 在 58 中农药检测中, 超标的有 13 种, 因此本文分析 13 种超标农药的安全风险, 运用 Matlab 设计基于 IFS 的蔬菜质量安全风险评估算法 ( 详细代码见附录 2), 得到深圳市 2017 年 1-3 月份蔬菜食品安全指数, 如表 42 所示 表 年 1 月 -3 月深圳市蔬菜检测 IFS 农药名称 IFS ADI(mg/kg) 最大检出值 R(mg/kg) 氟虫腈 毒死蜱 敌敌畏 啶虫脒 多菌灵 甲胺磷 甲拌磷 ~ 29 ~

30 甲基异柳磷 克百威 氯氰菊酯 辛硫磷 溴氰菊酯 氧乐果 IFS 结果分析 : 风险系数大于 1 时, 表示该类蔬菜含有的农药有损害人体健康的风险 由表 28 可知, 风险系数大于 1 的农药有氟虫腈 毒死蜱 敌敌畏, 氧乐果, 其中氟虫腈的 IFS 值达 12.99, 存在的安全隐患比较大 蔬菜检测不合格的 13 种农药给深圳市带来的安全风险为 : IFS = 蔬菜安全风险系数大于 1, 说明深圳市 2017 年 1 月 -3 月的蔬菜食品的农药残留高于蔬菜食品安全标准, 考虑忽略渠道的蔬菜比例, 则风险评估值为 R=1.61, 不合格蔬菜食品引起的安全风险较大, 存在有较大的蔬菜食品安全风险, 深圳市仍需采取相关措施加强蔬菜质量安全检测 农药名称 龙岗区 宝安区 表 43 深圳市各辖区蔬菜检测 IFS 福田区 深圳市各辖区 IFS 南山区 罗湖区 龙华新区 光明新区 坪山新区 大鹏新区 氟虫腈 毒死蜱 敌敌畏 啶虫脒 多菌灵 甲胺磷 甲拌磷 甲基异柳磷 克百威 氯氰菊酯 辛硫磷 溴氰菊酯 氧乐果 IFS 盐田区 结果分析 : 由表 43 可知, 深圳市其中的 9 个辖区的蔬菜食品安全风险小, 处于安全状态, 而光明新区的 IFS =1.153, 大于 1, 存在蔬菜食品安全风险 大鹏新区和盐田区均未检出蔬菜农药残留超标, 因此食品安全级别较高 6.3 基于模糊理论的蔬菜安全风险评估 ~ 30 ~

31 本文发现, 基于 IFS 模型评价蔬菜的安全系数, 只考虑了超标农药最大剂量, 若农药检测出现异常值, 对模型影响较大 因此, 为更好衡量蔬菜安全风险, 结合各种蔬菜输入渠道比例, 蔬菜种类及合格率, 建立了基于模糊理论的蔬菜评估模型 模型建立令 U 表示影响因素集, 将蔬菜输入渠道作为一级指标, 各渠道的输入类型作为二级指标, 令 U i 表示第 i 种输入渠道蔬菜合格率,u ij 表示第 i 种渠道的第 j 种蔬菜合格率,ω i 表示第 i 种渠道蔬菜输入比例,α ij 表示第 i 种渠道中第 j 种 蔬菜输入比例, 则可表示为 U 0 = {U 1, U 2, U 3, U 4 } U 1 = {u 1,1, u 1,2,, u 1,8 } U 2 = {u 2,1, u 2,2,, u 2,8 } U 3 = {u 3,1, u 3,2,, u 3,8 } { U 4 = {u 4,1, u 4,2,, u 4,8 } 一级指标权重 : W = [ω 1, ω 2, ω 3, ω 4 ] 二级指标权重 : α 11 α 1j T = ( ) α i1 α ij (1) 指标集合和评价集合 U = { 鳞茎类, 叶菜类, 芸薹属类, 瓜类, 茄果类, 豆类, 根茎类, 水生类 } V = { 非常低, 低风险, 中等风险, 高风险, 非常高 } (2) 设置隶属函数采用梯形隶属函数, 记 M = (a, b, c, d), a b c d, 则 0, x a x a, a x b b a μ(x) = 1, b x c d x, c x d d c { 0, d x (3) 建立 U,V 的模糊评价矩阵 R r 11 r 1m R = ( ) r n1 r nm (4) 设置指标的权数分配 A (5) 矩阵合成 A = (a 1, a 2,, a n ), a i = 1 n i=1 ~ 31 ~

32 (6) 对 B 进行归一化处理 Y = A R b 2 C = ( b 1,,, b m ) b i b i b i 模型求解 (1) 模糊风险级别评价 : 根据各类蔬菜抽检的合格率, 将蔬菜风险设为 5 个等级, 如下表所示表 44 蔬菜风险等级表 风险等级 取值区间 非常高 (0.07,0.08,0.09,0.10) 高 (0.05,0.06,0.07,0.08) 中等 (0.03,0.04,0.05,0.06) 低 (0.01,0.02,0.03,0.04) 非常低 (0,0,0.01,0.02) (2) 构建模糊矩阵 R 农贸 = ( ) R 商场超市 = ( 图 45 风险值隶属度函数 ) R 批发市场 = ( ) R 生产基地 = ( ) ~ 32 ~

33 (3) 设置指标的权数分配 A A 农贸市场 = (0.009,0.3163,0.253,0.1446,0.1054,0.0904,0.0422,0.0392) A 批发市场 = (0.0314,0.2199,0.2618,0.1728,0.1623,0.0681,0.0681,0.0157) A 商品超市 = (0.0267,0.3316,0.2406,0.123,0.1016,0.0856,0.0588,0.0321) A 生产基地 = (0,0.6,0.4,0,0,0,0,0) (4) 蔬菜风险评价矩阵合成 Y = A R 解得 : Y 农贸市场 = (0.1928,0.1655,0.2984,0,0.3434) Y 批发市场 = (0.3435,0.1266,0.2556,0,0.2618) Y 商品超市 = (0.2406,0,0.0752,0.4436,0.2400) Y 生产基地 = (0.6,0,0,0,0.4) (5) 风险综合评价指标的合成 C = [0,1,2,3,4], B = Y C 解得 :B 农贸市场 = , B 批发市场 = B 商品超市 = , B 生产基地 = 1.6 (6) 深圳市蔬菜质量安全风险评估一级指标 : U = {B 农贸市场, B 批发市场, B 商场超市, B 生产基地 } 一级指标权重 : 根据深圳市 2017 年 1-3 月抽检蔬菜渠道的比例分配权重 : W = [0.4442,0.2602,0.2562,0.0394] 风险值计算 :F = U W T 解得 :F = 结果分析根据深圳市 1-10 期蔬菜抽检数据, 结合模型分析可得, 深圳市蔬菜安全风险值为 :F =2.0757, 属于中度风险, 四大蔬菜输入渠道的安全风险分别为 : B 农贸 = , B 批发 = , B 超市 = , B 生产基地 = 1.6 模型求解考虑忽略渠道的蔬菜, 即 : 蔬菜安全风险 = 常规渠道安全风险评估 + 忽略渠道导致的安全风险, 我们进行如下计算 R(F)=F*(1+2%) 解得 :R(F)=2.117 ~ 33 ~

34 R(B 农贸 ) = 2.179, R(B 批发 ) = 1.719, R(B 超市 ) = 2.490, R(B 基地 ) = 可见, 深圳市输入渠道安全风险在低风险 ~ 中等风险之间, 仍需在蔬菜抽检 风险监管等各层面提高蔬菜质量安全管理 6.4 农药风险系数模型蔬菜中农药的残留会对人体带来不同程度的健康危害 根据深圳市蔬菜抽检结果的数据指标, 为直观反映农药残留产生的安全风险, 考虑农药的超标率 农药毒性等因素, 定义农药的风险系数 R c (1) 模型建立 R c = ap c + b F c + S c P c = N c /N 0 a,b 为权重因子, 风险系数 R 与农药超标率 P 危害程度 S 成正比, 与农药的施检频率 F 成反比,P S F 均为一段时间的计算值, 根据风险系数随时间的动态变化 将风险 R 划分为 3 个级别 : 表 46 不同类型农药敏感因子 风险级别 R 低度风险 R [0,1.5) 中度风险 R (1.5,2.5] 高度风险 R (2.5, + ] (2) 模型求解建立中期风险系数, 时间段为 3 个月, 根据深圳市 2017 年 1 月 -3 月蔬菜抽检结果统计, 合格率为 95% 以上, 综合风险系数的历史运用, 选取权重因子 a=100,b=0.1 进行模型求解计算 对 58 种农药进行敏感因子 S 的估算, 将近期蔬菜检出的高危农药的敏感因子 S 设为 2, 近期检出的常规农药的 S 值设为 1, 在蔬菜种植中比较少施放的农药的敏感因子 S 设为 0.5, 如表 47 所示 表 47 不同类型农药敏感因子 敏感因子 S 农药类型 2 甲拌磷 甲基异柳磷 克百威 氧乐果 硫丹 涕灭威 水胺硫磷 百草枯 乙酰甲胺磷 乐果 甲胺磷 1 苯醚甲环唑 敌敌畏 啶虫脒 毒死蜱 多菌灵 氟虫腈 腐霉利 氯氟氰菊酯 灭蝇胺 甲氨基阿维菌素苯甲酸盐 氯氟氰菊酯 0.5 丙溴磷 阿维菌素 吡虫啉 咪鲜胺 二嗪磷等 (3) 结果分析根据深圳市 1-10 期蔬菜抽检数据, 得到蔬菜各种农药残留中期危险系数, 如表 48 所示,13 种超标农药中, 毒死蜱风险系数为 3.18 属于高风险 敌敌畏 多菌灵 氯氰菊酯 辛硫磷 溴氰菊酯危险系数小于 1.5 属于低度风险 其余 7 种农药危险系数属于中度风险 整体而言, 深圳市蔬菜农药残留危险系数为中度风险, 原因是在蔬菜中毒死蜱超标率较高且还含有其他高危农药 监管部门应加 ~ 34 ~

35 强蔬菜种植中毒死蜱 甲胺磷 氧乐果 克百威 甲拌磷 甲基异柳磷农药的使用管理规范, 严格限制使用剂量 表 48 深圳市 2017 年 1 月 -3 月蔬菜农药残留中期危险系数 农药类型 N 0 N 1 F P S R 风险级别 氟虫腈 中度风险 毒死蜱 高度风险 敌敌畏 低度风险 啶虫脒 中度风险 多菌灵 低度风险 甲胺磷 中度风险 甲拌磷 中度风险 甲基异柳磷 中度风险 克百威 中度风险 氯氰菊酯 低度风险 辛硫磷 低度风险 溴氰菊酯 低度风险 氧乐果 中度风险 6.5 蔬菜抽样检验设计 (1) 蔬菜抽检方案基本要求评估蔬菜质量安全风险的前提是抽检蔬菜的样本要有代表性, 因此, 选择样本遵循以下原则 : 1. 覆盖深圳市各个辖区, 统计分析各辖区消费蔬菜的输入渠道 数量, 是抽检蔬菜样本覆盖批发市场 农贸市场 大小型超市等 2. 抽样检验方案需切实可行, 需要考虑抽检所需的成本 时间及抽检数据的及时处理 (2) 蔬菜抽样总数量计算本文根据简单随机抽样的原则对蔬菜抽检总量进行估算, 令 N 表示样本总量,T 表示为 2017 年深圳市蔬菜流通总量 : T uα 2 p (1 p ) 2 N = T δ 2 + uα 2 p (1 p ) 2 深圳市拥有人口 2000 万, 一年蔬菜消费量极大, 可认为 : T δ 2 uα 2 p (1 p ) 2 解得 :N = uα 2 p (1 p ) 2 δ 2 深圳市蔬菜来源于全国各地, 分析深圳市蔬菜抽检结果误差可根据全国蔬菜抽检合格率作为标准, 令 K i 表示第 i 种蔬菜的全国抽检合格率, B i 表示深圳市第 i 种蔬菜抽检的合格率,N 0 表示深圳 1-10 期抽检蔬菜数量,N 1 表示不合格蔬菜 ~ 35 ~

36 数量, 根据前面对深圳市蔬菜抽检和全国蔬菜抽检结果的统计分析, 求解得到以下指标 : p = N 1 N 0,N 1 = 42, N 0 = 761 解得 :p = 5.52% n δ 2 = N i (B N i K i ) 2 0 i=1 δ = 0.232% 取样本置信概率为 95%, 结合标准正态分布函数, 可得 : α = 0.05 uα = 解得 : N = (3) 样本数量按季节分配将一年分为四个季度, 按季度分配抽样数量, 根据前面卡方检验结果蔬菜合格率与季度有关, 因此, 为了使抽检资源得到合理分配, 以深圳市各季度蔬 菜抽检平均不合格率的比例分配样本数量 令 u m 表示各季度不合格率,P j 表示 季度抽样比例,N j 表示每个季度抽样数量, 得到一年四个季度的抽检样本数量, 如表 49 所示 P j = u m /(u 1 + u 2 + u 3 + u 4 ) N j = N P j m = 1,2,3,4 表 49 深圳市各季度抽检蔬菜分配数量 季度 第一季度 第二季度 第三季度第四季度 历史不合格率 u m 7.54% 4.65% 4.95% 7.43% 样本分配比例 P j 30.69% 18.93% 20.15% 30.23% 样本分配数量 N*P j (4) 样本数量按月分配将蔬菜抽检分为 12 个月, 根据各个季度的样本数量, 等比例分配到季度包含的月份中得到个月抽检数量 N k N k = N j 3, k = 1,2,3 12, j = 1,2,3,4 (5) 按分配蔬菜种类比例分配样本数量根据深圳市抽检蔬菜类型比例, 推断深圳市市民消费蔬菜的类型及其数量比例, 如表 50 所示, 将蔬菜抽检的样品量按比例分配 表 50 深圳市市民消费蔬菜类型比例 蔬菜类型 比例 (%) 月抽检数量 (k=1,2,3 12) 叶菜类 N k 芸薹类 N k ~ 36 ~

37 瓜类 N k 茄果类 N k 豆类 N k 根菜类 N k 葱蒜类 N k 其他蔬菜 N k (6) 深圳市辖区抽样数量分配根据深圳市各辖区的居民人口分布, 推断蔬菜数量的地区分布, 令 pe i 表示深圳市第 i 个辖区的人口比例, nx ik 表示 k 月份第 i 各辖区应抽检的蔬菜数量, 可得 nx ik = pe i N k 用矩阵表示为 NX = PE N PE = (pe 1 pe 2 pe 9 pe 10) nx 1,1 nx 1,12 解得 :NX = ( ) nx 10,1 nx 10,12 (7) 深圳市蔬菜输入渠道抽检分配令 ns ij 表示第 i 个辖区, 第 j 类渠道的蔬菜抽检样品数量,bc i 表示第 i 类渠道蔬菜消费比例 ns ij = bc i nx jk 矩阵表示形式为 NS = BC NX 6.6 深圳市新鲜蔬菜具体抽样检验方案 (1) 抽检数量及时间安排抽检总量为全年 批次, 分 12 个月完成, 根据季度不同, 每月抽检数量不一, 如表 51 所示 表 51 深圳市各月份抽检蔬菜分配数量 第一季度第二季度第三季度第四季度 月份 分配蔬菜 样本数量 每月抽检批次由各辖区共同完成, 每月分配具体的抽检任务数量, 各辖区抽检批次如表 52 所示 ~ 37 ~

38 表 52 深圳市各月份各辖区抽检蔬菜分配数量 月份 各辖区抽检批次 南山区福田区宝安区罗湖区龙华新区光明新区盐田区坪山新区大鹏新区龙岗新区 (2) 辖区蔬菜抽检点分布为保证抽检地区分布合理, 根据深圳市地图, 各辖区蔬菜抽检先采用 5 点取样法, 建立各辖区的主要蔬菜检测点, 再从这 5 个主监测点向周围蔬菜场所地区发散, 直到覆盖辖区内所有蔬菜场所 图 53 五点抽样设计图 (2) 蔬菜抽检场所各辖区抽检蔬菜需覆盖该地区内的所有经营场所, 包括农贸市场 批发市场 商场超市以及小型的零售超市,, 每个辖区按比例分配具体的抽检样本数量到各蔬菜场所, 分配比例如表 54 所示 表 54 深圳市各辖区抽检蔬菜场所分配数量 抽检场所农贸市场批发市场商场超市生产基地零售超市 抽检样本比例 40% 23% 21% 8% 8% (3) 抽检蔬菜抽检品种 ~ 38 ~

39 每个抽检场所的蔬菜抽检品种均以叶菜类 芸薹类 瓜类 茄果类 豆类 根菜类 葱蒜类蔬菜为主, 其他类蔬菜为辅, 样本中各类型蔬菜所占比例如表 55 所示 表 55 深圳市各蔬菜场所抽检蔬菜类型 蔬菜类型叶菜类芸薹类瓜类茄果类豆类根菜类葱蒜类其他类抽检比例 23% 27% 16% 12% 10% 4% 2% 6% (4) 其他蔬菜输入渠道抽检对策低收入人群较为关注日常开支, 而城中村里低收入人群趋向选择小型市场解决日常蔬菜需求 深圳市以行政村为单位的城中村有 241 个, 各辖区根据每月新鲜蔬菜抽检的批次, 先在辖区内的各城中村分配 1 个批次, 再均衡分配到辖区内各蔬菜渠道 (5) 监测项目 检测依据和判断依据及原则将检测的 58 种农药项目划分为重点检测项目和常规检测项目, 其中重点检测项目是近两年全国蔬菜抽检 深圳市蔬菜抽检中出现的超标农药项目, 其余为常规检测项目, 如表 56 所示 抽检蔬菜可能会存在污泥 杂草的其它影响抽检的物质, 对于蔬菜样品的处理方法参照 农药残留分析样本的采样方法 (NY/T ) 表 56 农药检测项目和检测依据 重点检测项目 常规检测项目 农药项目氟虫腈 毒死蜱 敌敌畏 啶虫脒 多菌灵 甲胺磷 甲拌磷 甲基异柳磷 克百威 氯氰菊酯 辛硫磷 溴氰菊酯 氧乐果 苯醚甲环唑 代森锰锌 滴滴涕 甲氨基阿维菌素苯甲酸盐 腐霉利 水胺硫磷 乐果 21 种农药 对硫磷 甲基对硫磷 乙酰甲胺磷 三唑磷 丙溴磷 杀螟硫磷 二嗪磷 马拉硫磷 亚胺硫磷 伏杀硫磷 六六六 氰戊菊酯 甲氰菊酯 氯氟氰菊酯 氟氯氰菊酯 联苯菊酯 氟胺氰菊酯 氟氰戊菊酯 三唑酮 百菌清 异菌脲 涕灭威 ( 包括涕灭威砜 涕灭威亚砜 ) 灭多威 甲奈威 三氯杀螨醇 腐霉利 五氯硝基苯 乙烯菌核利 哒螨灵 苯醚甲环唑 嘧霉胺 阿维菌素 除虫脲 灭幼脲 吡虫啉 烯酰吗啉 虫螨腈 咪鲜胺 嘧菌酯 二甲戊乐灵 噻虫嗪 氟啶脲等农药残留 检测依据 NY/T GB/T GB/T SN/T GB/T GB/T NY/T NY/T GB/T NY/T 等 6.7 多目标优化模型在蔬菜质量安全监测中, 需要分析蔬菜从种植到进入消费市场中的一系列过程中农药的效用函数, 抽检方案的设计不仅需要考虑蔬菜对消费者健康影响, 还需要考虑蔬菜生产提供方的经济效益 即在达到蔬菜农药使用安全标准下, 需要 ~ 39 ~

40 考虑如何提高蔬菜生产量及销售量, 提高蔬菜的经济利益 (1) 农药来源途径分析在蔬菜种植 加工 运输 储存 销售等过程中, 都会存在加入农药或其他添加剂的可能, 以达到提高整体收益的目的 图 59 蔬菜污染来源及降解方向 图 60 农药对蔬菜的影响逻辑图 (2) 农药对蔬菜种植 销售的效用分析令 xi 表示蔬菜种植第 i 种农药或激素的施放剂量 在施放农药或植物激素后 : Ns 表示种植蔬菜 a 的损失量,Ts 表示种植蔬菜 a 节省成本,Zs 表示种植蔬菜 a 节省成本, 则 : 蔬菜 a 损失函数 : 种植蔬菜 a 节省时间函数 k N s = p i (x i ) i=1 ~ 40 ~

41 蔬菜 a 节省成本函数 蔬菜 a 保鲜时间函数 (3) 抽检结果的误差值检验 k T s = q i (x i ) i=1 m a Z a = z i (x i ) i=1 m a B a = z 0 + z i i=1 设 x ji 表示深圳市第 j 个辖区的第 i 类蔬菜抽检的合格率, 则深圳市各辖区的 各类蔬菜抽检合格率 X 可表示为 : x 11 x 1 i X = ( ) x j1 xj i 抽检结果误差值计算 S: n S = (B i K i ) 2 i=1 B i = x ji (4) 蔬菜抽检成本令 F 表示开展蔬菜抽检活动所需的成本,G 表示检测一批蔬菜所需的成本,M 表示蔬菜抽检的批次数量, Z 表示抽查需要的资金成本,t 表示抽检一批蔬菜所需要的时间,Te 表示蔬菜所有批次数量抽检完成需要的时间,ZC 表示抽检所需的总成本 : Z = F + G M Te = M t 则抽检过程所需的总成本 10 j=1 得到 :ZC = f(z, Te) 对于 f 的表示, 根据抽检难度的不同, 需要实际中衡量时间成本和资金花费成本的占比, 本文假设其为线性关系, 则抽检成本函数 f(z, Te) = a Z + b Te + C 0 (5) 多目标优化抽检方案设计需要考虑花费成本 抽检结果的误差值 蔬菜经济效益等因素, 令 QZ 表示抽检方案的可执行度, 则 QZ 为多目标影响的指标, 令 S 0 表示抽检允许的误差,ZC 0 表示抽检的最高成本限制, 则 QZ 约束条件为 ~ 41 ~

42 { S S 0 ZC ZC 0 ZC = a Z + b Te + C 0 n s = (B i K i ) 2 i=1 max (N s ) max (T s ) max (Z a ) max (B a ) 基于以上抽样检验模型, 根据深圳市历史蔬菜抽检数据, 求解可得到抽检流程的样本数量及其具体分布, 并以此为基础设计深圳市新鲜蔬菜抽检方案, 能够在达到抽检标准条件下, 更好的平衡蔬菜经济利益, 抽检成本, 减少对蔬菜市场带来的影响 七 模型评价与推广 7.1 模型的优点 (1) 模型 1 建立合格率 超标倍数等指标对蔬菜抽检数据进行了详细的描述性统计分析, 精确地展示了各个数据的统计结果 (2) 模型 2 基于统计分析的数据, 建立 IFS 模型, 考虑多种影响因素, 将蔬菜安全风险进行量化, 使深圳市蔬菜质量安全状况得到清晰的展示 (3) 模型 3 建立农药风险系数, 根据蔬菜中农药的超标率及是否为高危农药指标, 量化农药风险, 直观展示深圳市蔬菜农药残留的风险 (5) 基于随机抽样及不等概率抽样方法设计深圳市蔬菜抽检方案, 对抽检资源配置进行了优化 (4) 基于模糊理论蔬菜风险评价模型, 综合了蔬菜合格率 各类蔬菜比例 输入渠道比例, 较好分析了深圳市蔬菜安全风险状况 7.2 模型的缺点 (1) 模型未考虑蔬菜农药超标原因是种植污染, 还是运输管理过程中的二次污染, 对风险来源分析不够深入 (2) 基于 IFS 模型的安全系数只考虑了超标农药最大剂量, 若农药检测出现异常值, 对模型影响较大 (3) 抽样检测方案的设计中没有考虑各季度的反季节蔬菜 蔬菜来源地区参数指标, 仍存在优化抽检方案的空间 7.3 模型的推广方向模型可应用于多因素导致的安全风险评估, 如应用于医药检验安全性 区域空气环境评估 水源品质评估等 八 推动深圳市蔬菜食品安全质量提高的建议信 根据本文研究结果, 从蔬菜食品安全检查 蔬菜抽检制度 蔬菜种植管理方 ~ 42 ~

43 面予以参考性建议建议 (1) 增加蔬菜安全风险高的辖区抽检比例 : 根据本文建立的模型求解, 发现光明新区 龙岗区 宝安区的蔬菜安全风险系数 IFS 分别为 , 风险系数较大, 盐田区 大鹏新区在 1-10 期的蔬菜抽检合格率为 100%, 蔬菜安全风险低, 因此在蔬菜抽检工作安排上, 可适当减少盐田区 大鹏新区抽检工作, 并加大在光明新区 龙岗区 宝安区的抽检比例 (2) 多方协作跨区域蔬菜安全监测 : 网购蔬菜的安全性必须从源头得到保证, 根据本文对生鲜农产品电商平台的发展分析, 目前网购蔬菜的比例约为 0.94%, 但生鲜蔬菜电商平台规模以 65%-80% 的年增长速率扩大, 随人们对蔬菜的种类质量要求的不断提高, 未来网购生鲜蔬菜会成为深圳市不可或缺的蔬菜输入渠道, 需要食品监督机构监控各大电商平台蔬菜安全检测情况, 及时了解电商平台蔬菜来源, 并与来源区域的蔬菜检测机构协作, 共同保证蔬菜食品的安全 (3) 把控蔬菜来源地区 : 深圳市蔬菜属于输入型, 来源于全国各个地区, 根据本文对近两年全国蔬菜抽检结果的统计分析发现, 不同地区的蔬菜生产不合格情况有差异, 在选择蔬菜输入时, 应综合考虑蔬菜种植地区, 从源头上减少不合格蔬菜的出现 (4) 提高叶菜类抽查比重 : 根据全国蔬菜抽检结果分析, 不合格蔬菜中 32% 属于叶菜类蔬菜, 而深圳市四季消费的叶菜类蔬菜均达到 20%, 叶菜类蔬菜消费需求高且不合格数量较多, 应在蔬菜抽检中提高叶菜类蔬菜的抽检比重, 能够提高蔬菜的安全防控能力 (5) 实行散户种植蔬菜合并制度 : 根据本文的模型分析, 深圳市流动摊贩的蔬菜 78% 属于散户种植的蔬菜, 这给蔬菜抽检带来很大难度, 本文建议对散户密集种植蔬菜的区域进行合并起来, 建立批量生产的蔬菜种植基地, 既方便了散户种植和销售蔬菜, 也能更好地对蔬菜进行质量安全抽检 (6) 加强对易超标农药残留的检测 : 根据对深圳市蔬菜抽检结果及全国蔬菜农残抽检的统计分析, 蔬菜农药残留中毒死蜱 百克威 氧乐果的超标率最高, 且超标倍数较大, 易引起食品安全问题 需要在蔬菜安全抽检过程中着重检测这些农药在蔬菜中的残留量, 同时也应告诫蔬菜种植企业减少这些农药的使用量 九 参考文献 [1] 董秀金. 区域农业环境与蔬菜质量安全风险评价 - 以浙江省为例 [D]. 浙江大学, [2] 深圳万人市场调查 年上半年深圳居民购物行为调查报告 [J], 连锁业共享论坛,2016. [3] 中国产业调研网 年中国生鲜电商市场现状分析及发展趋势预测, [4] 黄耿志, 薛德升等. 城市流动的微区位选择机制 [J]. 人文地理.2016,1(147): [5] 吴建国. 数学建模案例精篇 [M]. 中国水利水电出版社 [6] 贾敬敬, 陈春明. 中国食品安全态势分析 [M]. 北京 : 中国科学技术出版社, 2003:16. [7] 何雨帆, 左方华等. 广西蔬菜农药残留动态及质量安全风险评估 [J]. 安徽农业科学.2011,39(1): ~ 43 ~

44 十 附录 附录 1: 统计各季节新鲜蔬菜种类及比例的 Matlab 代码 [data1,textdata1,aum1]=xlsread('d\shenzhen \veg_winter'); [data2,textdata2,aum2] = xlsread('d\ ShenZhen \veg_fall'); [data3,textdata3,aum3] = xlsread('d \ ShenZhen \veg_summer'); [data4,textdata4,aum4] = xlsread('d\ ShenZhen \veg_spring'); veg_winter=cell(1000,2); veg_fall=cell(1000,2); veg_spring=cell(1000,2); veg_summer=cell(1000,2); for t=1:1:1000 veg_winter{t,2}=1; veg_fall{t,2}=1; veg_spring{t,2}=1; veg_summer{t,2}=1; end Nj=1; C=1; k=1; veg_winter{k,1}=aum1{2,3}; for i=3:1:length(aum1) for j=1:1:k if(strcmp(aum1{i,3},veg_winter{j,1})) veg_winter{j,2}=veg_winter{j,2}+1; C=0; end end if C==0 C=1; else veg_winter{k+1,1}= aum1{i,3}; k=k+1; end end C=1; k=1; veg_fall{k,1}=aum2{2,3}; for i=3:1:length(aum2) ~ 44 ~

45 for j=1:1:k if(strcmp(aum2{i,3},veg_fall{j,1})) veg_fall{j,2}=veg_fall{j,2}+1; C=0; end end if C==0 C=1; else veg_fall{k+1,1}= aum2{i,3}; k=k+1; end end C=1; k=1; veg_summer{k,1}=aum3{2,3}; for i=3:1:length(aum3) for j=1:1:k if(strcmp(aum3{i,3},veg_summer{j,1})) veg_summer{j,2}=veg_summer{j,2}+1; C=0; end end if C==0 C=1; else veg_summer{k+1,1}= aum3{i,3}; k=k+1; end end C=1; k=1; veg_spring{k,1}=aum4{2,3}; for i=3:1:length(aum4) for j=1:1:k if(strcmp(aum4{i,3},veg_spring{j,1})) veg_spring{j,2}=veg_spring{j,2}+1; C=0; end ~ 45 ~

46 end if C==0 C=1; else veg_spring{k+1,1}= aum4{i,3}; k=k+1; end end xlswrite('c:\users\administrator\desktop\veg_winter.xls', veg_winter); xlswrite('c:\users\administrator\desktop\ veg_fall.xls',veg_fall); xlswrite('c:\users\administrator\desktop\veg_summer.xls', veg_summer); xlswrite('c:\users\administrator\desktop\veg_spring.xls', veg_spring); clc; clear; [numdata1,textdata1,aum1]=xlsread('c:\users\administrator \Desktop\ veg_spring.xls '); [a,b]=size(aum1); k=0; for i=1:1:length(aum1) k=k+aum1{i,2}; end for j=1:1:length(aum1) aum1{j,2}=aum1{j,2}/k; end xlswrite('c:\users\administrator\desktop\ com_spring_rite.xls',aum1); 附录 2: 蔬菜超标率统计的 Matlab 代码 [numdata1,textdata1,aum1]=xlsread('c:\users\administrator \Desktop\rg\ thy'); [a,b]=size(aum1); cj=cell(length(aum1),4); for i=1:1:length(aum1) cj{i,1}=aum1{i,1}; cj{i,2}=aum1{i,2}; cj{i,3}=aum1{i,3}; ~ 46 ~

47 cj{i,4}=(aum1{i,2}-aum1{i,3})/aum1{i,3}; end xlswrite('c:\users\administrator\desktop\rg\ cbl.xls',cj); 附录 3: 基于 IFS 指标评估蔬菜质量安全风险的 Matlab 代码 [numdata1,textdata1,aum1]=xlsread('c:\users\administrator \Desktop\sz\yu'); E=1; f=1; P=1; F=278; mb=61; ADI=cell(13,1); ifsc=cell(14,1); R=cell(13,2); for i=1:1:13 R{i,1}=aum1{i+1,1}; R{i,2}= aum1{i+1,2}; ADI{i,1}=aum1{i+1,3}; end for j=1:1:13 EDI=(R{j,2}*F*E*P); af=edi/1000; SI=ADI{j}; ifsc{j,1}=(af*f)/(si*mb); end xlswrite('c:\users\administrator\desktop\sz\ert.xls',ifsc ); 附录 4: 蔬菜名称规范化表 类别 主要品种 别名 萝卜 葖 芦菔 莱菔 胡萝卜 红萝卜 黄萝卜 丁香萝卜 药性萝卜 番萝卜 根菜类芜菁 蔓菁 圆根 盘菜 芜菁甘蓝 洋蔓菁 洋疙瘩 洋大头菜 根菾菜 红菜头 紫菜头 美洲防风 欧防风芹菜萝卜 蒲芹萝卜 牛蒡 东洋萝卜 蝙蝠刺 大力子 ~ 47 ~

48 ~ 48 ~ 根芹菜根洋芹菜 球根塘蒿婆罗门参西洋牛蒡 蒜叶婆罗门参山葵黑婆罗门参菊牛蒡 鸦葱薯芋类马铃薯土豆 山药蛋 洋芋 地蛋 荷兰薯 爪哇薯姜生姜 黄姜芋芋头 芋艿 毛芋魔芋磨芋 蒟蒻山药薯蓣 白苕 脚板苕 山薯 大薯 佛掌薯甘薯红薯 白薯 地瓜 番薯 红薯豆薯地瓜 凉薯 沙葛 新罗葛葛菊芋洋姜 鬼子姜菜用土圞儿香芋 美洲土圞儿蕉芋蕉藕 姜芋 食用美人蕉 食用莲蕉草石蚕螺丝菜 宝塔菜 甘露儿 地蚕葱蒜类韭草钟乳 起阳草 懒人菜 韭黄大葱木葱 汉葱洋葱葱头 圆葱大蒜胡蒜 蒜蒜薹蒜苔蒜苗蒜黄 青蒜分葱四季葱 菜葱 冬葱胡葱蒜头葱 瓣子葱 火葱 肉葱细香葱四季葱 香葱 蝦夷葱韭葱扁葱 扁叶葱 洋蒜苗楼葱龙爪葱 龙角葱蕹藠头 藠子 莱芝白菜类大白菜结球白菜 黄芽菜 包心白菜普通白菜白菜 小白菜 青菜 油菜乌塌菜榻菜 塌棵菜 榻地菘 黑菜菜薹菜心 绿菜薹 菜尖薹菜紫菜薹红菜薹芥菜类茎芥茎瘤芥 青菜头 菜头 包包菜 羊角菜 菱角菜抱子芥 儿菜 娃娃菜笋子芥 棒菜叶芥青菜 苦菜 春菜 辣菜 雪里蕻根芥辣疙瘩 冲菜 芥头 大头菜 疙瘩菜薹芥

49 ~ 49 ~ 甘蓝类结球甘蓝洋白菜 卷心菜 包心菜 椰菜 莲花白 包包白 圆白菜 茴子白球茎甘蓝苤蓝 菘 玉蔓菁 芥蓝头 擘蓝花椰菜花菜 菜花青花菜木立花椰菜 意大利花椰菜 嫩茎花椰菜 绿菜花 西兰花芥蓝白花芥蓝抱子甘蓝芽甘蓝 子持甘蓝羽衣甘蓝绿叶甘蓝 菜用羽衣甘蓝 叶牡丹 花包菜叶菜类菠菜菠稜菜 赤根菜 角菜 波斯草莴苣千金菜 莴笋 生菜 青笋 莴苣笋 莴菜 油麦菜 莜麦菜芹菜芹 药芹 苦堇 堇葵 堇菜 旱芹蕹菜空心菜 竹叶菜 通菜 藤菜 蓊菜苋菜米苋 赤苋 刺苋 青香苋 苋叶菾菜莙荙菜 牛皮菜 厚皮菜 光菜 叶甜菜菊苣欧洲菊苣 苞菜 吉康菜 法国莒荬菜冬寒菜冬苋菜 冬葵 葵菜 滑肠菜落葵木耳菜 软浆叶 软姜子 染浆叶 胭脂豆 豆腐菜 藤菜 紫果菜茼蒿蒿子秆 大叶茼蒿 蓬蒿 春菊芫荽香菜 胡荽 香荽茴香小茴香 香丝菜 结球茴香 鲜茎茴香 甜茴香菊花脑菊花叶 黄菊籽 路边黄 黄菊仔荠菜护生草 菱角菜 地米草 扇子草菜苜蓿草头 金花菜 黄花苜蓿 刺苜蓿 南苜蓿 黄花草子番杏新西兰菠菜 洋菠菜 夏菠菜 白番苋 海滨莴苣 宾菜 蔓菜苦苣花叶生菜 花苣紫背天葵血皮菜 紫背菜 红凤菜 观音苋 双色三七草罗勒毛罗勒 九层塔 零陵香 兰香草 光明子 省头草马齿苋马齿菜 长命菜 五行草 瓜子菜 马蛇子菜紫苏荏 赤苏 白苏 香苏 苏叶 桂荏 回回苏榆钱菠菜食用滨藜 洋菠菜 山菠菜 山菠薐草薄荷山野薄荷 蕃荷菜莳萝土茴香 草茴香 小茴香鸭儿芹鸭脚板 三叶芹 山芹菜 野蜀葵 三蜀葵 水芹菜蕺菜鱼腥草 蕺儿根 侧耳根 狗帖耳 鱼鳞草 菹菜蒲公英黄花苗 黄花地丁 婆婆丁 蒲公草马兰马兰头 红梗菜 紫菊 田边菊 马兰菊 鸡儿肠 竹节草

50 ~ 50 ~ 香芹菜荷兰芹 洋芫荽 欧芹 法国香菜 旱芹菜珍珠菜角菜 白苞蒿 山芹菜 珍珠花菜 甜菜子 鸭脚艾 乳白艾瓜类黄瓜王瓜 胡瓜冬瓜东瓜节瓜毛瓜南瓜中国南瓜 倭瓜 番瓜 饭瓜笋瓜印度南瓜 玉瓜 北瓜西葫芦美洲南瓜 蔓瓜 白瓜 香瓜越瓜白瓜 脆瓜 酥瓜 梢瓜菜瓜蛇甜瓜 老羊瓜 酱瓜丝瓜圆筒丝瓜 蛮瓜 水瓜 棱角丝瓜 胜瓜苦瓜凉瓜瓠瓜扁蒲 葫芦 蒲瓜 夜开花 瓠子佛手瓜洋丝瓜 合掌瓜 菜肴梨 瓦瓜 万年瓜 拳头瓜蛇瓜蛇豆 蛇丝瓜 长豆角茄果类番茄西红柿 番柿 柿子 洋柿子茄子落苏辣椒海椒 辣子 辣角 番椒青椒大椒 灯笼椒 柿子椒 彩色甜椒 甜椒酸浆红姑娘 洋姑娘 灯笼草 洛神珠豆类菜豆四季豆 芸豆 芸扁豆 豆角 刀豆 敏豆 玉豆 油豆长豇豆长豆角 豆角 带豆 裙带豆菜用大豆毛豆 枝豆豌豆青斑豆 麻豆 青小豆 荷兰豆 淮豆 留豆 金豆 麦豆 回回豆 甜豌豆蚕豆胡豆 佛豆 寒豆 罗汉豆扁豆峨眉豆 沿篱豆 眉豆 肉豆 龙爪豆莱豆金甲豆 科马豆 荷豆 玉豆 雪豆 洋扁豆 白豆 状元豆 棉豆 荷包豆刀豆大刀豆 关刀豆 洋刀豆多花菜豆红花菜豆 大白芸豆 大花芸豆 看花豆四棱豆翼豆 翅豆 四角豆 杨桃豆 热带大豆 四稔豆藜豆黧豆 黎豆 猫猫豆 毛毛豆 毛胡豆 毛狗豆 小狗豆 狸豆 八升豆 狗爪豆水生蔬菜莲藕莲 藕 荷茭白茭瓜 茭笋 菰首慈姑剪刀草 燕尾草水芹刀芹 楚葵 蜀芹 紫堇 蕲荸荠马蹄 地栗 乌芋 凫茈

51 ~ 51 ~ 菱菱角 龙角 水栗豆瓣菜西洋菜 水田芥 水蔊菜芡实鸡头米 鸡头 水底黄蜂莼菜马蹄草 水葵 水荷叶 湖菜 露葵蒲菜香蒲 甘蒲海带江白菜 昆布紫菜多年生及杂类蔬菜笋用竹竹笋芦笋石刁柏 龙须菜黄花菜萱草 金针菜百合夜合 中篷花香椿香椿树 红椿 椿花 椿甜树枸杞枸杞菜 枸杞头 枸杞芽蘘荷阳霍 野姜 蘘草 茗荷 苴蓴菜蓟朝鲜蓟 洋蓟 荷兰百合 法国百合辣根西洋山萮菜 山葵萝卜食用大黄原叶大黄 圆叶大黄黄秋葵秋葵 羊角豆桔梗地参 四叶菜 绿花根 铃铛花 沙油菜 梗草 道拉基蕨蕨菜 蕨苔 龙头菜 蕨儿菜鹿蕨菜乾苔发菜 头发菜 石发蒌蒿芦蒿 水蒿 香艾蒿 小艾 水艾薇菜野豌豆, 大巢菜, 斑矛架, 野苕子车前草车轮菜 牛舌菜 蛤蟆衣食用菊甘菊 臭菊玉米笋玉笋 多穗玉米 珍珠笋 番麦笋嫩玉米菜玉米 菜苞谷 青玉谷 御麦糯玉米中国玉米 糯苞谷甜玉米甜苞谷 甜玉蜀黍 菜玉米食用菌香菇香蕈 冬菇 花菇双孢蘑菇蘑菇 白蘑菇 双孢菇 洋菇 褐蘑菇 棕色蘑菇糙皮侧耳平菇 北风菌 青蘑 桐子菌草菇兰花菇 美味包脚菇 秆菇 麻菇 中国蘑菇金针菇毛柄金钱菇 冬菇 朴菇 朴菰黑木耳木耳 光木耳 云耳银耳白木耳 雪耳猴头菇猴头蘑 刺猬菌毛头鬼伞鸡腿蘑 鸡腿菇姬松茸巴西蘑菇 巴氏蘑菇茶薪菇杨树菇 柱状田头菇 柳环菌 茶树菇

52 ~ 52 ~ 真姬菇玉蕈 斑玉蕈 蟹味菇 胶玉蘑 鸿喜菇 海鲜菇灰树花贝叶多孔菌 云蕈 栗蘑 舞茸 莲花菌 千佛菌滑菇珍珠菇 光帽鳞伞 滑子蘑刺芹侧耳雪茸 干贝菇 杏鲍菇白灵侧耳白灵菇阿魏侧耳阿魏菇盖襄侧耳台湾平菇 鲍鱼菇毛木耳竹荪长裙竹荪 短裙竹荪 棘托竹荪肺形侧耳姬菇 秀珍菇 小平菇金顶侧耳榆黄蘑 玉皇菇大球盖菇长根菇大杯蕈猪肚菇洛巴伊口蘑金福菇北冬虫夏草蛹虫草牛肝菌松茸鸡枞羊肚菌榛蘑蜜环菌鸡油菌红菇口蘑蒙古口蘑青冈菌离褶伞一窝鸡芽苗菜绿豆芽绿豆芽幼芽黄豆芽黄豆幼芽黑豆芽黑豆幼芽青豆芽青豆幼芽红豆芽红豆幼芽蚕豆芽蚕豆幼芽红小豆苗红小豆幼芽豌豆苗豌豆幼芽花生芽花生幼芽苜蓿芽苜蓿幼芽或幼苗小扁豆芽小扁豆幼芽或幼苗萝卜芽萝卜芽幼苗菘蓝芽菘蓝幼芽或幼苗沙芥芽沙芥幼芽或幼苗芥菜芽芥菜幼芽或幼苗

53 芥蓝芽白菜芽独行菜芽香椿苗向日葵芽荞麦芽胡椒芽紫苏芽水芹芽小麦苗胡麻芽蕹菜芽芝麻芽黄秋葵芽花椒脑芽球菊苣苦苣芽佛手瓜稍辣椒尖豌豆尖草芽碧玉笋 芥蓝幼芽或幼苗白菜幼芽或幼苗独行菜幼苗香椿幼苗向日葵幼芽荞麦幼苗胡椒幼芽或幼苗紫苏幼芽或幼苗水芹幼苗小麦幼苗胡麻幼芽或幼苗蕹菜幼苗芝麻幼芽或幼苗黄秋葵幼苗花椒嫩芽菊苣芽球苦苣幼芽或幼苗佛手瓜幼稍辣椒幼稍豌豆幼稍草芽幼嫩假茎黄花菜幼嫩假茎 附录 5: 各月份适宜种植的蔬菜类型 月份适合播种的蔬菜 1 月油菜 四月曼 菠菜 芥蓝 生菜 马铃薯 葱 茄子 番茄 辣椒 芋头 芫荽 2 四月曼 菠菜 芥蓝 生菜 马铃薯 葱 芫荽 芋瓠 黄瓜 四季豆 茄子 番茄 青花菜 辣椒 芋头 白菜 萝卜 甘蓝 3 四月曼 菠菜 芥蓝 白菜 萝卜 芋瓠 黄瓜 四季豆 茄子 番茄 丝瓜 冬瓜 南瓜 苦瓜 辣椒 芋头 葱 芫荽 豇豆 毛豆 空心菜 苋菜 甘蓝 4 白菜 萝卜 芋瓠 黄瓜 四季豆 茄子 番茄 丝瓜 冬瓜 南瓜 苦瓜 辣椒 芋头 葱 芫荽 豇豆 毛豆 空心菜 甘蓝 油菜 苋菜 韭菜 芹菜 5 白菜 萝卜 芋瓠 黄瓜 四季豆 茄子 番茄 丝瓜 冬瓜 南瓜 苦瓜 葱 芫荽 豇豆 毛豆 空心菜 油菜 苋菜 甘蓝 芹菜 韭菜 6 空心菜 甘蓝 油菜 苋菜 韭菜 白菜 黄瓜 豇豆 毛豆 苦瓜 葱 萝卜 芹菜 花椰菜 7 空心菜 甘蓝 油菜 苋菜 韭菜 白菜 黄瓜 豇豆 毛豆 苦 ~ 53 ~

54 瓜 葱 萝卜 花椰菜 四季豆 茄子 番茄 芹菜 8 空心菜 甘蓝 油菜 苋菜 韭菜 白菜 黄瓜 豇豆 毛豆 苦瓜 葱 萝卜 花椰菜 四季豆 芹菜 青花菜 豌豆 胡萝卜 大蒜 9 空心菜 甘蓝 油菜 苋菜 韭菜 白菜 菠菜 生菜 葱 芫荽 青菜 青花菜 豌豆 胡萝卜 大蒜 萝卜 花椰菜 茼蒿 10 甘蓝 油菜 韭菜 白菜 菠菜 芥蓝 生菜 葱 芫荽 芹菜 青花菜 豌豆 胡萝卜 大蒜 萝卜 花椰菜 茼蒿 11 甘蓝 油菜 韭菜 白菜 菠菜 芥蓝 生菜 葱 芫荽 芹菜 青花菜 豌豆 胡萝卜 大蒜 萝卜 花椰菜 茼蒿 12 甘蓝 油菜 韭菜 白菜 菠菜 芥蓝 生菜 葱 芫荽 芹菜 青花菜 豌豆 胡萝卜 大蒜 萝卜 花椰菜 茼蒿 茄子 番茄 辣椒 附录 6: 深圳市抽检蔬菜指标 58 种农药的 ADI 值 ( 根据 GB 查询 ) 农药名称 ADI/(mg kg 1 ) 农药名称 ADI/(mg kg 1 ) 敌敌畏 氯氟氰菊酯 0.02 啶虫脒 0.07 氟氯氰菊酯 0.04 毒死蜱 0.01 溴氰菊酯 0.04 多菌灵 0.03 联苯菊酯 0.01 氟虫腈 氟胺氰菊酯 甲胺磷 氟氰戊菊酯 0.02 甲拌磷 三唑酮 0.03 甲基异柳磷 百菌清 0.02 百克威 异菌脲 0.06 氯氰菊酯 0.02 涕灭威 辛硫磷 灭多威 0.02 溴氰菊酯 0.01 甲奈威 氧化乐果 三氯杀螨醇 对硫磷 腐霉利 0.1 甲基对硫磷 五氯硝基苯 0.01 水胺硫磷 乙烯菌核利 0.01 乐果 氟虫腈 乙酰甲胺磷 0.03 哒螨灵 三唑磷 0.03 苯醚甲环唑 0.01 丙溴磷 0.03 嘧霉胺 0.2 杀螟硫磷 阿维菌素 二嗪磷 除虫脲 0.02 氟啶脲 灭幼脲 1.25 亚胺硫磷 0.01 吡虫啉 0.06 伏杀硫磷 0.02 甲氨基阿维菌素 - - 苯甲酸盐 ~ 54 ~

55 六六六 烯酰吗啉 0.2 氯氰菊酯 0.02 虫螨腈 0.03 氰戊菊酯 0.02 咪鲜胺 0.01 甲氰菊酯 0.03 嘧菌酯 0.2 二甲戊乐灵 0.03 噻虫嗪 ~ 55 ~