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

2019/8/261食品安全风险评估统计模型及软件开发2019/8/262风险评估及统计模型CDEEM构建原理及数据库CDEEM软件及评估结果CDEEM验证与评价主要内容2019/8/263食品安全的2个策略零风险选择，即“黑与白”的方法“凡是对人和动物有致癌作用的化学物不得加入食品”（Delancy修正案）100%安全与零风险，人的生命与健康最重要、不惜代价风险分析“Everysubstanceisapoison–itisonlyamatterofdose”（Paracelsus）没有100%安全的食品，仅能将风险渐低到可接受水平2019/8/264风险评估（RiskAssessment）框架4.风险表征Riskcharacterisation1.危害鉴定Hazardidentification2.危害表征Hazardcharacterisation3.膳食暴露Dietaryexposure•国际标准的制定依赖暴露评估数据（发达国家主导）•不同国家食物消费量与食品污染水平不同•成员国有权制定与国际标准不一致的国家标准•国际标准不一定适合特定食物消费量与污染水平•标准的设立是保护大多数人，而不是平均人•暴露评估不是均数简单计算，需进行高暴露水平和特殊人群评估（大数据库）发展中国家必须开展并重视的研究2019/8/265现有技术标准现状膳食暴露评估方法对于获得的膳食污染物数据必须与食物消费量结合，暴露评估有三种方法：•以食品污染水平和食物消费量均为事前进行估计的简单点评估模型；•以污染物水平进行事前估计乘以食物消费量分布数据的分布点评估模型；•以污染物和消费量均为分布数据的随机概率评估模型。2019/8/266传统方法（非概率模型）点评估：固定食品消费量×固定污染水平=摄入量平均暴露量=食物消费均值×食物污染中位水平高端暴露量=高端食品消费量×高污染水平不能提供暴露量的可能范围，通常是保守的估计，作为筛选目的.分布模型：食品消费量分布数据X固定污染水平=摄入量范围（特定消费模式）比点估计更加精确，但因取决于假设仍属于保守。点评估和简单分布方法趋向用“最坏情况”假设，而不考虑化学物在食品中存在的概率，不同食品中化学物的污染水平不同，或者食物消费量不同。2019/8/267污染物数据库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.22.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暴露量=Σ消费量iX污染水平i点评估概率模型•简单到复杂•解决估计数据的不确定性问题•风险管理从意义不大到十分有意义过度•避免了点评估中的争论和对于实际摄入量过高估计概率模型分析•消费量数据库2019/8/268从点评估到概率分析2019/8/269国外研究进展欧盟：FP5，“随机化模拟人类对化学物和营养素暴露的开发、验证和应用”FP6，“发展全新整合的食品风险分析确保食品安全”•概率性暴露评估ProbabilisticExposureAssessment•急性暴露与慢性暴露风险AcuteandChronicrisks美国：一些主要概率评价模型ModelnameexpandednamedevelopermoduleSHEDSStochastichumanexposuredosesystemEPA/ORDdietary,nodietary,dermalandinhalationDEEM(Calender)DietaryexposureevaluationmodelExponentdietary,nodietary,dermalandinhalationCARESCumulativeandAggregateRiskEvaluationSystemDowagrosciencedietary,nodietary,dermalandinhalationLIFELINELifelinegroup,incdietary,nodietary,dermalandinhalation2019/8/2610结果0.010.1110100100010000020406080100PercentileDose•群体暴露量•群体剂量0.010.1110100020406080100PercentileExposure模型•计算个体暴露/剂量分布tot1TIMEEXP.orDOSEIngestiontot1TIMEEXP.orDOSEDermaltot1TIMEEXP.orDOSEInhalation数据库输入暴露分布•污染物数据库•食物消费量数据库•食品加工系数SHEDS模型StochasticHumanExposureandDoseSimulation(SHEDS)Model2019/8/2611中国开展的暴露评估•食物消费量数据：全国营养调查（1959,1982,1992,2002）•食品污染水平数据污染物监测网、总膳食研究（1990，1992，2000）、双份饭研究。人体组织/体液的直接监测，如母乳中有机氯化合物、二噁英的浓度。•2002年中国居民营养与健康状况调查分层多阶段等容整群随机抽样方法膳食调查：食品频数法、3天24小时回顾法、3天入户称重法。2019/8/2612国内研究现状•我国在卫生标准起草过程中，已经利用我国膳食营养调查提出的食物消费量模式和食品污染物现状调查，利用点评估技术进行暴露评估，提出我国的限量标准指标以及在采用国际标准时通过暴露评估评价限量标准指标的适用性，但这一过程属于筛选性质的多，进行定量暴露评估的较少，与国际水平存在一定差距。•1990年、1992年和2000年成功进行的3次“中国总膳食研究”，在研究内容的完整性和研究方法的科学性上已基本接近发达国家水平。但在将危险性评估理论和数理统计、计算机技术相结合，开展膳食暴露定量评估模型研究上却要落后和迟缓一些。2019/8/2613风险评估统计模型•点估计–IESTI•概率模型：–急性：–慢性：ipkijkijkijwcxy1ipkkjiijkijwcxy1**“理想”条件下的概率模型：实际应用的概率模型：ipkkijkijwcxy12019/8/2614IESTI(InternationalEstimateofShortTermIntake)•由世界粮农组织和世界卫生组织农药残留专家联席会议(JMPR)确定的欧盟及国际权威机构在制定农残最大限量标准（MRL）时计算农药急性暴露量所采用的方法。•它主要是针对一天内摄入可能引起急性反应的农药，以食品为对象，选取某食品消费人群的高端消费量和食品污染物监测样本中的高残留量计算24小时内膳食暴露量，与急性参考剂量比较进行评估。•为解决单位食品个体间的残留差异，引入了单位重量和变异因子两个概念。2019/8/2615•LP：高端消费量，摄食者消费量的P97.5，即能涵盖消费人群中97.5%的食用者每天的消费量；•HR：高残留量，食品污染物监测的各检测样品的最大残留浓度值；•U：食品单位重量，由污染物监测实验提供的单位食品重量的中位数；•V：变异因子，单位食品高端残留量（P97.5）与检测样品的平均残留量的比值；•P：加工因子，加工后食品中的污染物残留浓度与加工前原始农产品中污染物浓度的比值；•bw：消费人群的平均体重。参数解释2019/8/2616•情形1•食品单位重量小于25ｇ，如原始或经加工的农产品（谷物、小麦、油料种子及豆类等小粒农作物）；这种情况也适用于肉类食品、蛋类、肝脏、肾脏等可食动物内脏类食品。bwP)-HRorHR(*LPIESTI四种情形（情形1、2a、2b和3）2019/8/2617•食品单位重量大于25g，引入一个变异因子。–情形2a单位食品重量小于消费人群的每日高端消费量LP。如：桃、李等水果。假设个体一天内消费多于1个单位重量的某食品，且第一个单位重量的该食品残留水平为[HR×v]，其余为HR。–情形2b单位食品重量超过消费人群的每日高端消费量LP。如：大西瓜、大白菜等。假设个体一天内仅消费小于等于1个单位重量的某食品，且消费部分残留水平为[HR×v]。bwP)-HRorHP(*U)-LP(v*P)-HRor(HR*UIESTIbwv*P)-HRorHP(*LPIESTI情形22019/8/2618•对于散装或多种成分混合的加工食品，如果汁、牛奶等，以监测实验获得的各检测样品残留浓度的中位数STMR代表可能的最高残留浓度。bwP-STMR*LPIESTI情形32019/8/2619•WHO收集了一些国家部分水果、蔬菜等食品的单位重量信息（表1）。•食品的单位重量2019/8/2620•变异因子的引入是为了解决混合样品中食品个体间的残留差异。单个水果中的最高残留浓度可能要比混合样品高5-10倍。目前，关于变异因子的研究和讨论还在进行。变异因子2019/8/2621概率模型•概率法对所评价化学物在食品中存在概率与污染水平（浓度）及相关水平的消费量进行模拟。这种方法需要足够的食品中化学物浓度和食物消费量数据，评价才有意义。通常情况下，法定的市场监督检查并不能提供统计学上有意义的分布特征。如果所获资料质量有保证，概率方法可能是最合适的方法：通过科学的抽样，将食物中某化学物浓度与实际含有该物质的食品消费量结合起来，从而提供了一个真实的暴露评价基础，来估计某物质的暴露量是否超过预定的安全阈值。•是否对食物消费量数据和污染物监测数据进行分布拟合：–参数方法：用适当的理论分布对数据进行拟合–非参数方法：直接根据样本的经验分布（大样本）2019/8/2622参数法与非参数法的选择2019/8/2623以急性暴露评估为例建模•污染物浓度模型–非参：从现有数据随机抽样（经验分布）–参数：二项分布/对数正态分布•处理效应（加工因子）模型–非参：无（1）；固定（fk,nom，fk,upp）–参数：对数或logit转换后正态分布•变异性建模–非参：默认（按单位重量：１，３，５，７）–参数：贝塔模型、贝努力模型、对数正态模型2019/8/2624cdf1cdf2cdfmUncertaintyaboutexposureordoseforgivenpercentileUncertaintyaboutpercentileforgivenexposureordoseExposureorDosePercentileExposureorDosePercentile•变异Variability=在人群间真正存在的差异，是真实性的部分。•不确定度Uncertainty=缺乏知识或数据信息造成的测量值与真实值差别，原则上可能降低解决变异与不确定性难点：MonteCarlo方法量化变异性，Bootstrap方法量化不确定性2019/8/2625MonteCarlo方法•基本思想：以一个概率模型为基础，按照这个模型所描绘的过程，通过模拟实验结果，得出问题的近似解，即进行一种数字模拟实验。•三个主要步骤：构造或描述概率过程；实现从已知概率分布抽样；建立各种估计量。一般而言，对模拟次数M有一定要求，模拟次数较少，结果可能不稳定。2019/8/2626Bootstrap方法•以原始数据为基础进行有放回抽样的统计模拟方法,可用于研究一组数据某统计量的分布特征,适用于难以用常规方法解决的参数区间估计、假设检验等问题。•基本思想是:从样本含量为n的原始数据范围内有放回地随机抽取i(通常选取i≤n)个观察单位组成一个样本，每个观察单位每次被抽到的概率相等，所得样本称为Bootstrap样本。重复该过程若干次,记为B，所得B个Bootstrap样本均可看作是从相应经验分布中随机抽取而来。对于每个Bootstrap样本进行相关统计量（如p50、p95、p99）计算，就得到统计量的Bootstrap分布，它表示从总体中抽样获得的样本统计量的不确定性。2019