食品安全风险评估统计模型及软件开发课件.ppt

1、食品安全风险评估统计模型及软件开发演示文稿风险评估及统计模型风险评估及统计模型 CDEEM CDEEM 构建原理及数据库构建原理及数据库 CDEEM CDEEM 软件及评估结果软件及评估结果CDEEM CDEEM 验证与评价验证与评价食品安全的食品安全的2 2个策略个策略零风险选择，即零风险选择，即“黑与白黑与白”的方法的方法“凡是对人和动物有致癌作用的化学物不得加入食品凡是对人和动物有致癌作用的化学物不得加入食品”（DelancyDelancy修正案）修正案）100%100%安全与零风险，人的生命与健康最重要、不惜代价安全与零风险，人的生命与健康最重要、不惜代价风险分析风险分析“Every

2、substance is a poison“Every substance is a poison it is only a matter of dose”it is only a matter of dose”（ParacelsusParacelsus ）没有没有100%100%安全的食品，仅能将风险渐低到安全的食品，仅能将风险渐低到大多数人支持大多数人支持难以理解！难以理解！风险评估风险评估（Risk Assessment）框架框架4.风险表征风险表征 Risk Risk characterisationcharacterisation危害鉴定危害鉴定 Hazard identificat

3、ionHazard identification2.危害表征危害表征Hazard characterisationHazard characterisation3.3.膳食暴露膳食暴露Dietary exposureDietary exposure国际标准的制定依赖暴露评估数据国际标准的制定依赖暴露评估数据 （发达国家主导）（发达国家主导）不同国家食物消费量与食品污染水平不同不同国家食物消费量与食品污染水平不同成员国有权制定与国际标准不一致的国家标准成员国有权制定与国际标准不一致的国家标准国际标准不一定适合特定食物消费量与污染水平国际标准不一定适合特定食物消费量与污染水平标准的设立是保护大多数

4、人，而不是平均人标准的设立是保护大多数人，而不是平均人暴露评估不是均数简单计算，需进行高暴露暴露评估不是均数简单计算，需进行高暴露 水平和特殊人群评估（大数据库）水平和特殊人群评估（大数据库）现有技术标准现状现有技术标准现状膳食暴露评估方法对于获得的膳食污染物数据必须膳食暴露评估方法对于获得的膳食污染物数据必须与食物消费量结合，暴露评估有三种方法：与食物消费量结合，暴露评估有三种方法：以食品污染水平和食物消费量均为事前进行估计的简单点以食品污染水平和食物消费量均为事前进行估计的简单点评估模型；评估模型；以污染物水平进行事前估计乘以食物消费量分布数据的分以污染物水平进行事前估计乘以食物消费量分布

5、数据的分布点评估模型；布点评估模型；以污染物和消费量均为分布数据的随机概率评估模型。以污染物和消费量均为分布数据的随机概率评估模型。传统方法（非概率模型）传统方法（非概率模型）点评估：点评估：固定食品消费量固定食品消费量 固定污染水平固定污染水平 =摄入量摄入量 平均暴露量平均暴露量=食物消费均值食物消费均值 食物污染中位水平食物污染中位水平 高端暴露量高端暴露量=高端食品消费量高端食品消费量 高污染水平高污染水平不能提供暴露量的可能范围，通常是保守的估计，作不能提供暴露量的可能范围，通常是保守的估计，作为筛选目的为筛选目的.分布模型：分布模型：食品消费量分布数据食品消费量分布数据X X固定污

6、染水平固定污染水平=摄入量摄入量范围（特定消费模式）范围（特定消费模式）比点估计更加精确，但因取决于假设仍属于保守。比点估计更加精确，但因取决于假设仍属于保守。点评估和简单分布方法趋向用点评估和简单分布方法趋向用“最坏情况最坏情况”假设，假设，而不考虑化学物在食品中存在的概率，不同食品而不考虑化学物在食品中存在的概率，不同食品中化学物的污染水平不同，或者食物消费量不同。中化学物的污染水平不同，或者食物消费量不同。污染物污染物数据库数据库0.00.20.40.00.71.42.22.93.6-0.7-1.4-2.2-2.9-3.6消费量和残留数据随机采样0.00.20.40.00.71.42.2

7、2.93.6-0.7-1.4-2.2-2.9-3.60.00.20.40.00.71.42.22.93.6-0.7-1.4-2.2-2.9-3.6暴露量=消费量i X 污染水平i点评估点评估 概率模型概率模型 简单到复杂简单到复杂 解决估计数据的不确定性问题解决估计数据的不确定性问题 风险管理从意义不大到十分有风险管理从意义不大到十分有意义过度意义过度 避免了点评估中的争论和对于避免了点评估中的争论和对于实际摄入量过高估计实际摄入量过高估计概率模型分析概率模型分析 消费量消费量数据库数据库“理想”条件下的概率模型：食品组别 groupdietary,nodietary,dermal and（发

8、达国家主导）食品消费量数据描述性分析结果腌制食品：蜜制和盐制，烟熏制品数据随机匹配，计算污染物暴露量，要求每种农产国际标准不一定适合特定食物消费量与污染水平Monte Carlo方法量化变异性，膳食调查数据中所调查的食物均采用膳食五位码。HR：高残留量，食品污染物监测的各检测样品的最大残留浓度值；Dietary exposure evaluation modelHR：高残留量，食品污染物监测的各检测样品的最大残留浓度值；概率性暴露评估Probabilistic Exposure Assessment以污染物和消费量均为分布数据的随机概率评估模型。实现从已知概率分布抽样；WHO收集了一些国家部分

9、水果、蔬菜等食品的单位重量信息（表1）。Monte Carlo方法量化变异性，膳食调查数据的描述性分析从点评估到概率分析从点评估到概率分析国外研究进展 欧盟：欧盟：FP5FP5，“随机化模拟人类对化学物和营养素暴露的开随机化模拟人类对化学物和营养素暴露的开 发、验证和应用发、验证和应用”FP6FP6，“发展全新整合的食品风险分析确保食品安全发展全新整合的食品风险分析确保食品安全”概率性暴露评估概率性暴露评估Probabilistic Exposure Assessment急性暴露与慢性暴露风险急性暴露与慢性暴露风险Acute and Chronic risks美国：一些主要概率评价模型美国：一

10、些主要概率评价模型Model Model namenameexpanded nameexpanded namedeveloperdevelopermodulemoduleSHEDSSHEDSStochastic human exposure dose Stochastic human exposure dose systemsystemEPA/ORDEPA/ORDdietary,nodietary,dermal and dietary,nodietary,dermal and inhalationinhalationDEEMDEEM(Calender)(Calender)Dietary exp

11、osure evaluation modelDietary exposure evaluation modelExponentExponentdietary,nodietary,dermal anddietary,nodietary,dermal and inhalation inhalationCARESCARESCumulative and Aggregate Risk Cumulative and Aggregate Risk Evaluation SystemEvaluation SystemDow Dow agroscienceagrosciencedietary,nodietary

12、,dermal and dietary,nodietary,dermal and inhalationinhalationLIFELINELIFELINELifeline group,Lifeline group,incincdietary,nodietary,dermal and dietary,nodietary,dermal and inhalationinhalation结果结果0.010.1110100100010000020406080100PercentileDose 群体暴露量群体暴露量 群体剂量群体剂量0.010.1110100020406080100PercentileEx

13、posure模型模型 计算个体暴露计算个体暴露/剂量分剂量分布布 tot1TIMEEXP.orDOSEIngestiontot1TIMEEXP.orDOSEDermaltot1TIMEEXP.orDOSEInhalation数据库输入数据库输入 暴露分布暴露分布 污染物数据库污染物数据库 食物消费量数据库食物消费量数据库 食品加工系数食品加工系数 SHEDS SHEDS 模型模型 S Stochastic tochastic H Human uman E Exposure and xposure and D Dose ose S Simulation(SHEDS)Modelimulation(

14、SHEDS)Model Bootstrap方法量化不确定性（01、02、03、04）WHO收集了一些国家部分水果、蔬菜等食品的单位重量信息（表1）。人口学数据库：2002年中国居民营养与健康状况调查家庭成员基本情况登记表收集的数据，主要提供年龄、性别、体重等信息。发、验证和应用”Bootstrap方法解决变异与不确定性难点：污染物总数据库合计监测了510种食物，135种污染物,487875条记录。假设个体一天内消费多于1个单位重量的某食品，且第一个单位重量的该食品残留水平为 HR v，其余为HR。概率模型与点估计模型估计结果比较（小麦粉，ug/kg/day）对于每个Bootstrap样本进行相

15、关统计量（如p50、p95、p99）计算，就得到统计量的Bootstrap分布，它表示从总体中抽样获得的样本统计量的不确定性。HR：高残留量，食品污染物监测的各检测样品的最大残留浓度值；（3）当中国的食物是具有中国特色的而CAC编码中没有的食物以同类食物两位字母代码加6000开始进行独立编码；perofARFD:急性参考剂量的百分数（%）它主要是针对一天内摄入可能引起急性反应的农药，以食品为对象，选取某食品消费人群的高端消费量和食品污染物监测样本中的高残留量计算24小时内膳食暴露量，与急性参考剂量比较进行评估。Evaluation System2005-2006年海关农产品进出口监测数据，有效

16、数据包括监测了44种食物，104种污染物，88843条记录。成员国有权制定与国际标准不一致的国家标准参数：二项分布/对数正态分布加工产品,熟食,加工方法：烘焙烤加工产品,熟食,加工方法：蒸dietary,nodietary,dermal and中国开展的暴露评估中国开展的暴露评估 食物消费量数据：食物消费量数据：全国营养调查全国营养调查 （1959,1982,1992,20021959,1982,1992,2002）食品污染水平数据食品污染水平数据污染物监测网、总膳食研究（污染物监测网、总膳食研究（19901990，19921992，20002000）、双）、双份饭研究。份饭研究。人体组织人体

17、组织/体液的直接监测，如母乳中有机氯化合物、体液的直接监测，如母乳中有机氯化合物、二噁英的浓度。二噁英的浓度。20022002年中国居民营养与健康状况调查年中国居民营养与健康状况调查分层多阶段等容整群随机抽样方法分层多阶段等容整群随机抽样方法膳食调查：膳食调查：食品频数法、食品频数法、3 3天天2424小时回顾法、小时回顾法、3 3天入户称重法。天入户称重法。国内研究现状 我国在卫生标准起草过程中，已经利用我国膳食我国在卫生标准起草过程中，已经利用我国膳食营养调查提出的食物消费量模式和食品污染物现营养调查提出的食物消费量模式和食品污染物现状调查，利用点评估技术进行暴露评估，提出我状调查，利用点

18、评估技术进行暴露评估，提出我国的限量标准指标以及在采用国际标准时通过暴国的限量标准指标以及在采用国际标准时通过暴露评估评价限量标准指标的适用性，但这一过程露评估评价限量标准指标的适用性，但这一过程属于筛选性质的多，进行定量暴露评估的较少，属于筛选性质的多，进行定量暴露评估的较少，与国际水平存在一定差距。与国际水平存在一定差距。19901990年、年、19921992年和年和20002000年成功进行的年成功进行的3 3次次“中国中国总膳食研究总膳食研究”，在研究内容的完整性和研究方法，在研究内容的完整性和研究方法的科学性上已基本接近发达国家水平。但在将危的科学性上已基本接近发达国家水平。但在将

19、危险性评估理论和数理统计、计算机技术相结合，险性评估理论和数理统计、计算机技术相结合，开展膳食暴露定量评估模型研究上却要落后和迟开展膳食暴露定量评估模型研究上却要落后和迟缓一些。缓一些。风险评估统计模型风险评估统计模型 点估计 IESTI 概率模型：急性：急性：慢性：慢性：ipkijkijkijwcxy 1ipkkjiijkijwcxy 1*“理想理想”条件下的概率模型：条件下的概率模型：实际应用的概率模型：实际应用的概率模型：ipkkijkijwcxy1IESTI(International Estimate of Short Term Intake)由世界粮农组织和世界卫生组织农药残留专家

20、联席会议(JMPR)确定的欧盟及国际权威机构在制定农残最大限量标准（MRL）时计算农药急性暴露量所采用的方法。它主要是针对一天内摄入可能引起急性反应的农药，以食品为对象，选取某食品消费人群的高端消费量和食品污染物监测样本中的高残留量计算24小时内膳食暴露量，与急性参考剂量比较进行评估。为解决单位食品个体间的残留差异，引入了单位重量和变异因子两个概念。LPLP：高端消费量，摄食者消费量的：高端消费量，摄食者消费量的P97.5P97.5，即能涵盖消费人，即能涵盖消费人群中群中97.5%97.5%的食用者每天的消费量；的食用者每天的消费量；HRHR：高残留量，食品污染物监测的各检测样品的最大残留：高

21、残留量，食品污染物监测的各检测样品的最大残留浓度值；浓度值；U U：食品单位重量，由污染物监测实验提供的单位食品重量：食品单位重量，由污染物监测实验提供的单位食品重量的中位数；的中位数；V V：变异因子，单位食品高端残留量（：变异因子，单位食品高端残留量（P97.5P97.5）与检测样品）与检测样品的平均残留量的比值；的平均残留量的比值；P P：加工因子，加工后食品中的污染物残留浓度与加工前原：加工因子，加工后食品中的污染物残留浓度与加工前原始农产品中污染物浓度的比值；始农产品中污染物浓度的比值；bwbw：消费人群的平均体重。：消费人群的平均体重。参数解释参数解释 情形情形1 1 食品单位重量

22、小于25，如原始或经加工的农产品（谷物、小麦、油料种子及豆类等小粒农作物）；这种情况也适用于肉类食品、蛋类、肝脏、肾脏等可食动物内脏类食品。bwP)-HRor HR(*LPIESTI 四种情形（情形四种情形（情形1 1、2a2a、2b 2b 和和3 3）食品单位重量大于食品单位重量大于25g25g，引入一个变异因子。，引入一个变异因子。情形情形2a2a 单位食品重量小于消费人群的每日高端消费量单位食品重量小于消费人群的每日高端消费量LPLP。如：。如：桃、李等水果。假设个体一天内消费多于桃、李等水果。假设个体一天内消费多于1 1个单位重量的某食品，个单位重量的某食品，且第一个单位重量的该食品残

23、留水平为且第一个单位重量的该食品残留水平为 HR HR v v，其余为，其余为HRHR。情形情形2b2b 单位食品重量超过消费人群的每日高端消费量单位食品重量超过消费人群的每日高端消费量LPLP。如：。如：大西瓜、大白菜等。假设个体一天内仅消费小于等于大西瓜、大白菜等。假设个体一天内仅消费小于等于1 1个单位重量个单位重量的某食品，且消费部分残留水平为的某食品，且消费部分残留水平为 HR HR v v。bwP)-HRor HP(*U)-LP(v*P)-HRor (HR*UIESTIbwv*P)-HRor HP(*LPIESTI 情形情形2 2 对于散装或多种成分混合的加工食品，如果汁、牛奶等，

24、以监测实验获得的各检测样品残留浓度的中位数STMR代表可能的最高残留浓度。bwP-STMR*LPIESTI 情形情形3 3 WHOWHO收集了一些国家部分水果、蔬菜等食品的单位重量收集了一些国家部分水果、蔬菜等食品的单位重量信息（表信息（表1 1）。）。食品的单位重量食品的单位重量 变异因子的引入是为了解决混合样品中食品个体间的残留变异因子的引入是为了解决混合样品中食品个体间的残留差异。单个水果中的最高残留浓度可能要比混合样品高差异。单个水果中的最高残留浓度可能要比混合样品高5-5-1010倍。目前，关于变异因子的研究和讨论还在进行。倍。目前，关于变异因子的研究和讨论还在进行。变异因子变异因子

25、概率模型概率模型 概率法对所评价化学物在食品中存在概率与污染水平（浓概率法对所评价化学物在食品中存在概率与污染水平（浓度）及相关水平的消费量进行模拟。度）及相关水平的消费量进行模拟。这种方法需要足够这种方法需要足够的食品中化学物浓度和食物消费量数据，评价才有意义。的食品中化学物浓度和食物消费量数据，评价才有意义。通常情况下，法定的市场监督检查并不能提供统计学上有通常情况下，法定的市场监督检查并不能提供统计学上有意义的分布特征。如果所获资料质量有保证，概率方法可意义的分布特征。如果所获资料质量有保证，概率方法可能是最合适的方法：通过科学的抽样，将食物中某化学物能是最合适的方法：通过科学的抽样，将

26、食物中某化学物浓度与实际含有该物质的食品消费量结合起来，从而提供浓度与实际含有该物质的食品消费量结合起来，从而提供了一个真实的暴露评价基础，来估计某物质的暴露量是否了一个真实的暴露评价基础，来估计某物质的暴露量是否超过预定的安全阈值。超过预定的安全阈值。是否对食物消费量数据和污染物监测数据进行分布拟合：是否对食物消费量数据和污染物监测数据进行分布拟合：参数方法：用适当的理论分布对数据进行拟合参数方法：用适当的理论分布对数据进行拟合 非参数方法：直接根据样本的经验分布（大样本）非参数方法：直接根据样本的经验分布（大样本）中国食品污染物暴露评估与国际接轨的基础为解决单位食品个体间的残留差异，引入了

27、单位重量和变异因子两个概念。与食物消费量结合，暴露评估有三种方法：急性概率评估模型分析结果Exposure or Dose桥梁数据库包含了连接膳食调查数据和污染物点评估和简单分布方法趋向用“最坏情况”假设，而不考虑化学物在食品中存在的概率，不同食品中化学物的污染水平不同，或者食物消费量不同。不能提供暴露量的可能范围，通常是保守的估计，作为筛选目的.污染物数据监测数据中农产品具有CAC编码（1）当中国的食物在CAC的食物编码中能找到相应的编码时，全部采用CAC的编码，编码类型为两位字母代码加四位数字；假设个体一天内消费多于1个单位重量的某食品，且第一个单位重量的该食品残留水平为 HR v，其余为

28、HR。中国食物成份表2004 北京大学医学院出版社膳食暴露评估方法对于获得的膳食污染物数据必须Hazard characterisation2005-2006年海关农产品进出口监测数据，有效数据包括监测了44种食物，104种污染物，88843条记录。罐装食品：各类罐装食品加工产品,熟食,加工方法：蒸平均暴露量=食物消费均值 食物污染中位水平参数方法：用适当的理论分布对数据进行拟合加工产品,熟食,加工方法：烘焙烤perofARFD:急性参考剂量的百分数（%）膳食暴露评估方法对于获得的膳食污染物数据必须参数法与非参数法的选择参数法与非参数法的选择以急性暴露评估为例建模以急性暴露评估为例建模 污染物

29、浓度模型污染物浓度模型 非参：非参：从现有数据随机抽样（经验分布）参数：参数：二项分布/对数正态分布 处理效应（加工因子）模型处理效应（加工因子）模型 非参：非参：无（1）；固定（fk,nom，fk,upp）参数：参数：对数或logit转换后正态分布 变异性建模变异性建模 非参：非参：默认（按单位重量：，）参数：参数：贝塔模型、贝努力模型、对数正态模型 cdf1cdf2cdfmUncertainty about exposure or dose for given percentile Uncertainty about percentilefor given exposure or dose

30、Exposure or DosePercentileExposure or DosePercentile变异Variability=在人群间真正存在的差异，是真实性的部分。不确定度Uncertainty=缺乏知识或数据信息造成的测量值与真实值差别，原则上可能降低 解决变异与不确定性难点：解决变异与不确定性难点：Monte CarloMonte Carlo方法量化变异性，方法量化变异性，BootstrapBootstrap方法量化不确定性方法量化不确定性 Monte CarloMonte Carlo方法方法 基本思想：以一个概率模型为基础，按照这个模型所描绘的过程，通过模拟实验结果，得出问题的近

31、似解，即进行一种数字模拟实验。三个主要步骤：构造或描述概率过程；实现从已知概率分布抽样；建立各种估计量。一般而言，对模拟次数一般而言，对模拟次数M M有一定要求，模拟次数较少，有一定要求，模拟次数较少，结果可能不稳定。结果可能不稳定。BootstrapBootstrap方法方法 以原始数据为基础进行有放回抽样的统计模拟方法,可用于研究一组数据某统计量的分布特征,适用于难以用常规方法解决的参数区间估计、假设检验等问题。基本思想是:从样本含量为n的原始数据范围内有放回地随机抽取i(通常选取i n)个观察单位组成一个样本，每个观察单位每次被抽到的概率相等，所得样本称为Bootstrap样本。重复该过

32、程若干次,记为B，所得B个Bootstrap样本均可看作是从相应经验分布中随机抽取而来。对于每个Bootstrap样本进行相关统计量（如p50、p95、p99）计算，就得到统计量的Bootstrap分布，它表示从总体中抽样获得的样本统计量的不确定性。模型原理模型原理定量模型的变异性定量模型的变异性（图（图1 1）估计模型的不确定估计模型的不确定度（图度（图2 2）1234膳食调查数据库及结构膳食调查数据库及结构 污染物监测数据库及结构污染物监测数据库及结构 农产品的编码及分类农产品的编码及分类 膳食数据与污染物数据整合膳食数据与污染物数据整合 模型数据基础模型数据基础 膳食量数据库：2002年

33、中国居民营养与健康状况调查中通过24小时膳食回顾法收集的有效数据，包括22567个家庭66172人连续3天共计193814人日、1810703条膳食量记录；人口学数据库：2002年中国居民营养与健康状况调查家庭成员基本情况登记表收集的数据，主要提供年龄、性别、体重等信息。膳食调查数据库膳食调查数据库标识变量人天数记录条数0（正常值）19381418092381（奇异值）119914062（删除值）5659flag flag 说明说明:2000-2006年全国污染物监测数据中，有效数据包括监测了497种食物，75种污染物，399032条记录。2005-2006年海关农产品进出口监测数据，有效数据

34、包括监测了44种食物，104种污染物，88843条记录。污染物总数据库合计监测了510种食物，135种污染物,487875条记录。污染物监测数据库污染物监测数据库一般而言，对模拟次数M有一定要求，模拟次数较少，结果可能不稳定。单个水果中的最高残留浓度可能要比混合样品高5-10倍。核心模块：bootstrap&Monte Carlo、暴露分布、模型验证；水平和特殊人群评估（大数据库）以污染物和消费量均为分布数据的随机概率评估模型。情形2b 单位食品重量超过消费人群的每日高端消费量LP。Bootstrap方法量化不确定性高暴露人群（P95）摄入各食物贡献度（ug/kg体重/天）Uncertaint

35、y about exposure单个水果中的最高残留浓度可能要比混合样品高5-10倍。数据随机匹配，计算污染物暴露量，要求每种农产假设个体一天内消费多于1个单位重量的某食品，且第一个单位重量的该食品残留水平为 HR v，其余为HR。合理的污染物监测农产品合并归组是桥梁数据库的建立的关键；假设个体一天内消费多于1个单位重量的某食品，且第一个单位重量的该食品残留水平为 HR v，其余为HR。Uncertainty about percentile2005-2006年海关农产品进出口监测数据，有效数据包括监测了44种食物，104种污染物，88843条记录。（2）当中国的食物是具体加工的食物时，在其原

36、料食物编码后面加两位数字01、02等；中国食品污染物暴露评估与国际接轨的基础点评估：固定食品消费量 固定污染水平=摄入量（01、02、03、04）食品频数法、3天24小时回顾法、3天入户称重法。中国短期摄入量估计值-NESTI（ug/kg体重/天）化学物名称：毒死蜱 中国食品污染物暴露评估与国际接轨中国食品污染物暴露评估与国际接轨的基础的基础 实现污染物数据库与膳食调查数据库实现污染物数据库与膳食调查数据库关联的重要途径关联的重要途径 识别同一食品不同表达的唯一符号，识别同一食品不同表达的唯一符号，便于后续计算机处理。便于后续计算机处理。农产品编码及分类农产品编码及分类食物编码主要参考四个国内

37、外文件：食物编码主要参考四个国内外文件：Pesticides Residues in Food,Second Edition,Section 2,Codex Alimentarius Volume 2,1993.GEMS/FOOD Consumption Cluster Diets.杨月欣杨月欣 王光亚王光亚 潘兴昌潘兴昌.中国食物成份表中国食物成份表2002 2002 北京大学北京大学医学院出版社医学院出版社 杨月欣杨月欣 王光亚王光亚 潘兴昌潘兴昌.中国食物成份表中国食物成份表2004 2004 北京大学北京大学医学院出版社医学院出版社编码总则编码总则 删除具有明显的地区性和国家性的食删除

38、具有明显的地区性和国家性的食物编码，保留普遍通用的食物编码；物编码，保留普遍通用的食物编码；删除删除CACCAC食物编码中关于饲料的食物食物编码中关于饲料的食物分类编码；分类编码；删除删除CACCAC食物编码中食物编码中4000-50004000-5000同物同物异名的重复编码；异名的重复编码；编码筛选原则编码筛选原则（1 1）当中国的食物在）当中国的食物在CACCAC的食物编码中能找到相的食物编码中能找到相应的编码时，全部采用应的编码时，全部采用CACCAC的编码，编码类型的编码，编码类型为两位字母代码加四位数字；为两位字母代码加四位数字；FB0269FB0269葡萄葡萄FB0275FB02

39、75草莓草莓FC0001FC0001柑桔柑桔FC0004FC0004橙子橙子FC0005FC0005柚子柚子FI0327FI0327香蕉香蕉FI0341FI0341猕猴桃猕猴桃FI0342FI0342桂圆桂圆FI0353FI0353菠萝菠萝FP0226FP0226苹果苹果FP0230FP0230梨梨编码方法编码方法（2 2）当中国的食物是具体加工的食物时，在其原）当中国的食物是具体加工的食物时，在其原料食物编码后面加两位数字料食物编码后面加两位数字0101、0202等；等；CF1211CF1211面粉面粉CF121101CF121101面条面条 *CF121102CF121102方便面方便面

40、*CF121103CF121103三明治三明治 *CF121104CF121104馒头馒头 *CF121105CF121105月饼月饼 *CF121106CF121106饼干饼干 *CF121107CF121107酥酥 *CF121108CF121108蛋糕蛋糕 *CF121109CF121109其他面点其他面点 *编码方法编码方法（3 3）当中国的食物是具有中国特色的而当中国的食物是具有中国特色的而CACCAC编编码中没有的食物以同类食物两位字母代码加码中没有的食物以同类食物两位字母代码加60006000开始进行独立编码；开始进行独立编码；IM1000IM1000蛤蛤IM1002IM1002

41、乌贼乌贼IM1003IM1003贻贝类贻贝类IM1005IM1005扇贝类扇贝类IM6000IM6000鲍鱼鲍鱼 *IM6001IM6001蚬蚬 *IM6002IM6002蛏子蛏子 *IM6003IM6003蚶蚶 *IM6004IM6004章鱼章鱼 *IM6005IM6005桃花蛸桃花蛸 *IM6006IM6006海参海参 *IM6007IM6007海蜇海蜇 *编码方法编码方法（4 4）对于一些食用频率较少又不常见的食物）对于一些食用频率较少又不常见的食物以同类食物代码加以同类食物代码加99999999表示。表示。VR0075VR0075根茎蔬菜根茎蔬菜FB0269FB0269葡萄葡萄VR05

42、05VR0505芋头芋头FB0275FB0275草莓草莓VR0506VR0506萝卜萝卜FB9999FB9999其他浆果其他浆果 *VR0508VR0508红薯红薯VC0428VC0428丝瓜丝瓜VR0577VR0577胡萝卜胡萝卜VC0429VC0429南瓜南瓜VR0589VR0589马铃薯马铃薯VC0432VC0432西瓜西瓜VR0600VR0600山药山药VC6000VC6000哈密瓜哈密瓜 *VR6000VR6000水生蔬菜水生蔬菜 *VC6001VC6001甜瓜甜瓜 *VR9999VR9999其他根茎蔬菜其他根茎蔬菜 *VC9999VC9999其他瓜其他瓜 *编码方法编码方法比点估计

43、更加精确，但因取决于假设仍属于保守。点评估 概率模型我国在卫生标准起草过程中，已经利用我国膳食营养调查提出的食物消费量模式和食品污染物现状调查，利用点评估技术进行暴露评估，提出我国的限量标准指标以及在采用国际标准时通过暴露评估评价限量标准指标的适用性，但这一过程属于筛选性质的多，进行定量暴露评估的较少，与国际水平存在一定差距。实际应用的概率模型：食品来源 category5，即能涵盖消费人群中97.Bootstrap方法兴奋中枢神经系统的食品污染物监测农产品的编码是桥梁数据库构建的基础；发展中国家必须开展并重视的研究以污染物水平进行事前估计乘以食物消费量分布数据的分布点评估模型；假设个体一天内

44、消费多于1个单位重量的某食品，且第一个单位重量的该食品残留水平为 HR v，其余为HR。急性概率评估模型分析结果污染物名称：铅Hazard characterisation合理的污染物监测农产品合并归组是桥梁数据库的建立的关键；加工产品,熟食,加工方法：煎炸假设个体一天内消费多于1个单位重量的某食品，且第一个单位重量的该食品残留水平为 HR v，其余为HR。Bootstrap方法急性概率评估模型分析结果污染物名称：铅概率模型与点估计模型估计结果比较（小麦粉，ug/kg/day）点评估和简单分布方法趋向用“最坏情况”假设，而不考虑化学物在食品中存在的概率，不同食品中化学物的污染水平不同，或者食物

45、消费量不同。平均暴露量=食物消费均值 食物污染中位水平食品来源食品来源 categorycategory食品类别食品类别 typetype食品组别食品组别 groupgroup食品序列号食品序列号 numbernumber食品加工方法食品加工方法 methodsmethods*每个食品赋予二个编码：国内六位编码和CAC食品编码（01、02、03、04）（A、B、C、D）（001、002、003）（FC0001、FC0004）（00、01、02）五级分类编码（参照五级分类编码（参照CACCAC）CDEEM 食物分类食物分类A级级 初级初级植物产品植物产品B级级 初级初级动物产品动物产品G级级 其

46、他加其他加工食品工食品*D级级 植物来植物来源加工食品源加工食品E级级 动物来动物来源加工食品源加工食品F级级 多组分多组分加工食品加工食品*01水果水果02蔬菜蔬菜03禾本禾本04坚果种子坚果种子05香料香料001柑橘类柑橘类002仁果类仁果类003核果类核果类004浆果类浆果类005热带水果热带水果FC0001柑桔柑桔FC0002柠檬柠檬FC0003蜜橘蜜橘FC0004橙子橙子FC0005柚子柚子CategoryTypeGroupCAC code食物谱系图食物谱系图 CAC编码编码植物类加工食品植物类加工食品植物来源单植物来源单组分产品组分产品 谷物碾磨部分谷物碾磨部分杨月欣的六位码杨月欣

47、的六位码加工方法加工方法加工方法编码（加工方法编码（2121种）种）加工代码加工代码加工方法加工方法0000原始状态原始状态0101碾压碾压(粉粉,汁汁)0202刨光刨光,削皮削皮0303加工产品加工产品,需再加工才可食用：面需再加工才可食用：面条条,粉丝粉丝,饺子等饺子等0404加工产品加工产品,熟食熟食,加工方法：蒸加工方法：蒸0505加工产品加工产品,熟食熟食,加工方法：烘加工方法：烘 焙焙 烤烤0606加工产品加工产品,熟食熟食,加工方法：煎加工方法：煎 炸炸0707加工产品加工产品,熟食熟食,加工方法：炒加工方法：炒0808加工产品加工产品,熟食熟食,加工方法：煮加工方法：煮 卤卤0

48、909加工产品加工产品,熟食熟食,加工方法：粥加工方法：粥 汤汤 羹羹加工代码加工代码加工方法加工方法1010加工产品加工产品,熟食熟食,加工方法：加工方法：冷冻冷冻1111罐装食品：各类罐装食品罐装食品：各类罐装食品1212腌制食品：蜜制和盐制，烟熏腌制食品：蜜制和盐制，烟熏制品制品1313发酵食品发酵食品1414脱水植物产品脱水植物产品1515脱水乳制品：奶粉脱水乳制品：奶粉1616脱水肉制品：肉松，牛肉干脱水肉制品：肉松，牛肉干1717海鲜干燥制品：虾皮，鱼干海鲜干燥制品：虾皮，鱼干1818婴儿配方食品婴儿配方食品1919冷冻食品类：蔬菜等冷冻食品类：蔬菜等9999加工方法不明加工方法不

49、明type_idtype_idtypetypetype_idtype_idtypetype1 1谷物谷物9 9蔬菜类蔬菜类2 2根茎类根茎类1010鱼类及海产食品鱼类及海产食品3 3豆类豆类1111蛋类蛋类4 4蔗糖和蜂蜜蔗糖和蜂蜜1212水果类水果类5 5坚果与油料种子坚果与油料种子1313奶类及其制品奶类及其制品6 6植物油和油脂植物油和油脂1414肉类及其脏器类肉类及其脏器类7 7兴奋中枢神经系统的食品兴奋中枢神经系统的食品1515动物油脂类动物油脂类8 8香料香料1616饮料饮料食品分类食品分类(GEMS/FOOD consumption cluster diets)(GEMS/FOO

50、D consumption cluster diets)食物分类食物分类参照参照20022002年食物成分表分类标准、国际食品法典委员会（年食物成分表分类标准、国际食品法典委员会（CACCAC）通用分类标准和通用分类标准和GEMS/FOOD Consumption Cluster Diets要求进要求进行分类：行分类：type_idtype_idtypetypetype_idtype_idtypetype1 1谷类及制品谷类及制品1212鱼虾蟹贝类鱼虾蟹贝类2 2薯类、淀粉及制品薯类、淀粉及制品1313婴幼儿食品婴幼儿食品3 3干豆类及制品干豆类及制品1414小吃、甜饼小吃、甜饼4 4蔬菜类及