食品安全评估

1食品安全问题的定量评估摘要本题是一个食品安全定量评估的问题。我们小组根据网上公布的数据，先选出几个比较突出的化学污染物（如镉、汞、砷、有机氯、农药DDT），找出中国人均食物日摄入量及种类，之后选取深圳的网上公布的数据分析研究深圳的食品安全问题。我们小组提出两个问题：一是如何确定污染物各自之间的比重及人均日摄入总量，二是化学污染物的摄入量对人体的危害程度的评价。针对问题一建立模型一：运用层次分析，结合CAC的食物分类标准将化学污染物分成3类，运用熵权法求权重求出各个污染物各自之间的比重及在日食入量中占的比重；模型二是运用概率评估模型进行分析的，我们知道对总人群来说，其每日摄入食物量和自身体重都是服从正态分布的。我们用MATLAB软件中normrnd函数分别模拟出一组服从正态分布的日摄入和体重的随机数，由某种化学污染物摄入量和安全摄入量得到某种化学污染物对人体安全危害指数；然后调用深圳市2001年全市食品安全检查数据，并对其进行检验，得出当地食品危害指数表如下图，并对其做出合理的分析。化学污染物名称浓度危害指数中毒程度镉0.010241.2安全无毒汞0.008128.9正常砷0.03556.2安全无毒有机氯0.0196.5安全无毒农药DDT0.0211.4正常最后我们对模型优缺点进行分析，并对模型中存在的问题进一步改进。关键词：CAC层次分析，熵权法，概率评估，危害指数2一、问题重述俗话说的好“安居乐业”，其实社会也是如此，如果必需品都没有安全保证，那社会又怎么会和谐！近年来，地沟油、硫磺姜、回炉面包、三鹿牛奶等食品安全事件频发，食品安全问题依然是当下中国最热门的一幕丑剧。当瘦肉精事件尚未平息，上海染色馒头又像一重磅炸弹引爆全国；随后接棒的则有温州毒馒头、广东墨汁粉条、辽宁毒豆芽等。常年风起云涌的食品安全事件，像一双双粗重的钢琴手，一次次肆虐地撩拨人们敏感的神经。为此，人大常委会建议：把食品安全作为国家安全组成部分，其重要性不亚于金融安全、粮食安全、能源安全、生态安全；加大对违反食品安全行为的惩处力度，以法治的方式维护公民的“食品安全权”；推进食品，尤其是农产品的规模化、产业化生产经验，从根本上保证食品质量；参照科学依据，结合中国国情，制定出全国统一的食品安全检验和认定标准；充分发挥社会组织作用，强化食品生产行业协会的作用，加强行业自律；加强宣传力度，将涉及到食品安全方面的法律法规汇编成册，针对从业者流动性强等特点，加大宣传力度并从严考核；充分发挥社会监督的作用，鼓励全社会共同参与食品安全监督等。食品安全问题内容广泛、复杂，为了突出主要削弱次要问题，我们小组选择食品中污染物的分布作为研究食品安全的一个突破口，建立数学模型，利用互联网数据，定量评估食品安全问题。二、问题分析建立数学模型即是抽象实际问题，由于关于食品安全问题的方面又太过繁杂，我们组仅选择其中一个方面即食品中污染物的分布进行研究，为了定量评估题目中涉及到的食品安全问题。1）问题的背景“民以食为天，食以安为先”，食品卫生安全关系国计民生，关系经济繁荣和社会稳定。近年来，随着人民生活水平的不断提高，食品安全问题更加受到人民的关注。我国是人口大国，食品消费量大，而粮食生产加工分散，因而食品卫生安全监控难度较大。学习借鉴欧美国家食品卫生安全监控的做法，建立化学污染在食品中的安全模型，用以处理调查或检测数据，对公共食品卫生安全做出评估应该是一条可行的途径。欧美已经建立了有关模型并且开始使用，用时间证明了其可行性。但是模型的建立和使用还存在着很多困难和问题。在建立食品安全模型时需要设计合理的周阳调查方案或者数据的搜集，并由抽样数据建立比较准确的人群食品摄入量模3型。化学污染物分布模型的建立面临的问题包括：怎样利用抽样率很低的随机抽样数据建立模型；怎么由随机变量取值大于某一数值的部分样本数据再加上其他可以利用的信息估计出这个随机变量的整体分布。危害评估模型就是利用模型一的结果对全国、某地区、某类食品的安全状况作出评价。2）问题分析通过查阅相关资料，我们了解到食品中对人体健康造成威胁的因素大致有三类：物理因素，化学因素，生物因素。其中，物理因素可以通过良好的生产操作加以避免，因此基本不做讨论；生物因素涉及到外部条件，以及该生物危害人体的繁殖和变异过程，具有较高的复杂性，也不易做研究；所以，我们主要分析化学因素。化学因素主要包括重金属、有机农药、化肥、食品添加剂等，他们会通过各种复杂的途径先污染食物，进而危害人体健康，根据所污染所产生的危害的严重性程度划分等级，确定食品中化学污染物的危害指数，进而完成对食品安全的定量评估。为此，我们需要解决一下两个方面的问题：问题一：查阅相关资料，根据所得数据求出各种食品中不同污染物的污染水平，及人均日摄入量中各化学污染物的含量。问题二：计算出各类化学元素的人均摄入量,结合问题一得到的分类后各类食物中所含的化学污染物的含量，怎样计算出各类化学污染物的人均摄入量和人体每日安全摄入量，如何根据这些数据确定出该化学污染物人均摄入量与危害指数的关系；并最终确定某种化学污染物实际测定值与食品安全的关系。三、模型假设3.1假设各化学污染物之间不会相互作用；3.2假设被评估人群在调查期内的体重是一个定值；3.3假设被评估人群的体重服从正态分布；3.4假设被评估人群每天摄入量是服从正态分布的；3.5假设化学残留物在食物中的各项指标不会随时间改变；3.6被评估人群中不同个体对化学污染物的耐受剂量与体重成正比。四、符号说明CAC……国际食品法典委员会食品与饲料分类标准ijX……指标组成的矩阵ijK……矩阵列向量归一化后的矩阵g……人均每日消耗食物量jw……第j个指标的熵定义4je……第j个指标的熵权值iA……第i种食物的摄入量ia……第i种食物中镉的浓度ib……第i种中食物汞的浓度ic……第i种食物中砷的浓度id……第i种食物中有机氯的浓度ie……第i种食物中农药DDT的浓度W……个体的体重mEXP……食品中某种化学污染物的人群摄入量(/)gkgmFSRI……食品安全危害指数FQPA……安全系数五、模型的建立与求解5.1模型的建立与求解首先，建立对食品中某种化学污染物摄入量的评估模型，得到食品中该化学污染物的摄入量评估值；其次，建立食品安全危害指数模型，将实际食品中化学污染物浓度的测定量与通过模型得到的危害指数等级表中的各级浓度范围进行比较，最终可定量评估出食品的安全程度。5.2问题一模型的建立我们参考国际食品法典委员会食品与饲料分类标准对食物进行分类。食品法典分类的标准是按4个层次来进行食品和动物饲料分类的。首先，按照原料来源及是否经过加工分为A、B、C、D、E5个等级。这5个等级分别是A级植物来源的初级食品、B级动物来源的初级食品、C级初级饲料、D级植物来源的加工食品、E级动物来源的加工食品。各大类型中根据污染物的相似性再分为不同的组，分类标准的每个组中都包含各种推荐指标，如潜在污染物残留量、作物和动物在正常情况下的消费量以及应用农药最大残留限量的商品比例等。所以我们结合CAC分类标准将部分食物进行分类，并分别研究各类食品中的不同化学污染物。由于该污染物的种类繁多，所以我们挑选对人体危害性较大的砷、汞、镉，有机氯和农药DDT进行研究。通过运用层次分析法将各种食品分成三大类，即植物性食物、动物性食物、其他类食物，得到如下层次结构图：5依据层次分析法确定指标分层体系，然后采用熵权法计算食品中各化学污染物所占的比重。熵权法的基本原理及模型：熵是系统无序程度的度量，它表示从一组不确定事物中提取信息量的多少，可用于度量已知数据所包含的有效信息量和确定权重，在食品安全评价中可以采用。通过熵的计算确定权重，即根据各项监测指标值的差异程度u，确定各指标的权重。1）数据矩阵归一化。设有m个评价指标，n个评价对象，其原始数据矩阵为：111212122212.....................nnijmmmnxxxxxxXxxx式中ijx（i=1，2,…，m；j=1,2，…，n）为第i个指标的第j个采样值。矩阵归一化：()ijijmnKk2）在m个评价指标n个评价对象的评价估计问题中，第j个指标的熵的定义为：1njijijiwkKLNK()kLNn（i=1，2，…，m）当ijK=0时，令0ijijkLNk定义第j个指标的熵后，则其熵权为：食物中化学污染物的比重植物性动物性其他粮食蔬菜水果肉类鱼类蛋奶类奶粉零食饮料、茶叶奶类果蔬罐头水果鱼肉罐头鱼类皮蛋蛋类61(1)/()njjjiewmw其中101;1mjjiww5.12问题一模型的求解：层次分析法中对于第四层的分类做一些说明，由于对于不同的化学污染物，有的调查数据不全包含所有指标，因此在计算比重时采用化学污染物含量与分类总量的比值作为指标的权重(平均值与样本数的乘积再求和与总样本数的比)。上述层次结构图中，我们以汞、镉、砷、有机氯和DDT共有四层，其中第四层对第三层的权重为其加权平均值，根据附表中的一、二、三、四、五表中数据，整理化简，然后运用熵权法，求出各层次每个成员的熵权值。Xij是各种化学污染物在三类食物中所占的比重组成的矩阵如下：785414.0709724.0635796.0042291.0306973.0187482.0125367.0199966.0345861.0225455.0027104.016491.0164239.0611848.0467571.0ijX第j个指标的熵值矩阵：547385.0729012.0825124.0729615.0959236.0jw五种化学污染物的权重矩阵：374177.0224026.014457.0223527.0033699.0je用熵权法求得的镉、汞、砷、有机氯、农药在植物性食物、动物性食物、其他食物中的比重如下表：化学污染物的比重表（1）化学污染物名称植物性食物动物性食物其他食物镉的含量比0.4675710.2254550.306973汞的含量比0.6118480.3458610.042291砷的含量比0.1642390.1999660.635796有机氯的含量比0、164910.1253670.709724农药DDT的含量比0、0271040.1874820.785414再次运用熵权法求镉、汞、砷、有机氯、农药DDT的权重比，结果如下：食物中各污染物的比重表（2）7名称镉汞砷有机氯农药DDT食物中比值0.0336990.2235270.144570.2240260.374177从表（1）中我们可以看到化学污染物镉在各类食物中比较均衡；汞在植食性食物中含量比较大，而在其他食物中含量比很小；砷在植食性和动食性含量比相当，都不大，而在其他食物中含量比很大；有机氯和农药DDT在其他食物中含量比远高于植食性和动食性食物。从表（2）中我们可以看出，就整体食用食物来说，有机氯和农药DDT占比重较大，汞和砷所占比重一般，而镉占最小。根据我们分析的结果，希望有关部门加强对这方面的注意力度5.2问题二模型的建立与求解5.21问题二模型的建立：由于模型一仅仅是对食品中化学污染物分布情况作分析，未能得出该类化学污染物对人体的危害程度，因此我们用概率评估模型对此做一分析。由于食品中化学污染物的危害作用与其进入人体的绝对量有关，因此评价食品安全以人体对化学污染物的实际摄入量与其安全摄入量比较更为科学合理，在这样一种理论背景下，导出可以用来评价食品中某种化学污染物浓度对人体影响的食品安全危害指数：mmmEXPFSRIRSI其中，mSI是某种化学污染物的安全摄入量，我们经查有关资料得出。由于以上模型主要针对于成人人群，考虑到不同人群的生理特性存在差异，引入一个安全系数FQPA，且FQPA可以由美国食品质量安全法确定，进而得到危害指数的最终模型，如下：mmmEXPFSRIRFQPASI其中我们经查资料得出安全系数FQPA=100，校正系数R=1。参照美国环保署(EPA)模式，即摄