[发明专利]花生黄曲霉毒素B1污染对人体健康风险的定量评估方法

摘要

本发明涉及花生黄曲霉毒素B1污染对人体健康风险的定量评估方法。它包括S1：整理得到符合样本数据分布特征的花生黄曲霉毒素B1污染数据分布集；S2：根据评估需要，依据相应的人群特征对待评估的区域人群进行划分，然后对划分而得的人群的相应数据整合备用，用于后续评估使用；S3：构建人群花生黄曲霉毒素B1暴露量的评估模型和摄入花生黄曲霉毒素B1致原发性肝癌超额病例数的评估模型：S4：整合得到花生黄曲霉毒素B1污染对人体健康风险的评估结果，包括变异性评估结果和不确定性评估结果。该评估方法用于花生黄曲霉毒素B1污染对人体健康风险定量评估，评估结果科学、准确，可给政府、企业等实施风险管理提供关键技术与数据支撑。

主权项

花生黄曲霉毒素B1污染对人体健康风险的定量评估方法，其特征在于：它包括以下步骤：S1：根据待评估的区域，选取与待评估区域相应的花生黄曲霉毒素B1污染原始数据，决定未检出值至检出限值或定量限值的数据取值，整理得到符合样本数据分布特征的花生黄曲霉毒素B1污染数据分布集；S2：根据评估需要，依据相应的人群特征对待评估区域人群进行划分，然后对划分而得的人群的相应数据包括该人群的花生消费量数据、体重数据进行整合、备用，用于后续评估使用；S3：构建评估模型，所述评估模型有两个，分别用于人群花生黄曲霉毒素B1暴露量的评估，和在前述评估结果基础上的花生黄曲霉毒素B1污染对人体健康风险定量评估即膳食摄入花生黄曲霉毒素致原发性肝癌超额风险即原发性肝癌超额病例数的评估：评估模型1为： EXP = C × X × F BW 其中：EXP‑人群花生黄曲霉毒素B1暴露量，单位为纳克每千克体重每天(ng/kgbw/day)；C‑花生中黄曲霉毒素B1污染浓度，单位为微克每千克(μg/kg)；X‑人体每日花生消费量，单位为克每天(g/day)；F‑加工因子，0＜F≤1；BW‑人体体重，单位为千克(kg)；评估模型2为：R＝EXP×[M×P+N×(1‑P)]其中：R‑原发性肝癌年超额病例数，单位为例每十万人每年(例/十万人/年)；EXP‑人群花生黄曲霉毒素B1暴露量，单位为纳克每千克体重每天(ng/kgbw/day)；P‑人群中乙型肝炎病毒携带率，单位％；M‑花生黄曲霉毒素B1对乙型肝炎病毒携带者的平均致癌强度，单位为例每十万人每年(例/十万人/年)；N‑花生黄曲霉毒素B1对非乙型肝炎病毒携带者的平均致癌强度，单位为例每十万人每年(例/十万人/年)；S4：分析选用上述评估模型并结合适合的取样方法进行评估，整合得到花生黄曲霉毒素B1污染对人体健康风险的定量评估结果，包括变异性评估结果和不确定性评估结果：变异性评估结果：将S2中经划分而得的待评估人群的花生消费量X、体重数据BW、加工因子F相应作为固定参数值，将花生黄曲霉毒素B1污染值C作为评估模型1的输入量，选择Monte Carlo抽样方法从S1中的花生黄曲霉毒素B1污染数据分布集中重复随机抽样，将由此随机抽取的黄曲霉毒素B1浓度对应输入评估模型1，计算获得人群花生黄曲霉毒素B1暴露量，然后通过评估模型2，将人群中乙型肝炎病毒携带率P、花生黄曲霉毒素B1对乙型肝炎病毒携带者的平均致癌强度M和花生黄曲霉毒素B1对非乙型肝炎病毒携带者的平均致癌强度N固定参数以及由评估模型1计算得到的人群花生黄曲霉毒素B1暴露量输入评估模型2，进行多次重复抽样模拟运算，分别得到待评估的特定区域人群人体中花生黄曲霉毒素B1暴露量概率分布和花生黄曲霉毒素B1污染对特定区人群人体健康风险即原发性肝癌年超额病例数的变异性评估结果；不确定性评估结果：基于S1的花生黄曲霉毒素B1污染数据集，样本容量为m，采用Bootstrap抽样方法从中有放回地随机抽取多个Bootstrap小样本，每个小样本的样本容量为i，i≤m，然后对每个小样本通过Monte Carlo模拟抽样，重复上述变异性评价过程即选择Monte Carlo抽样方法从各Bootstrap小样本中重复随机抽样，将由此随机抽取的黄曲霉毒素B1浓度对应输入评估模型1，计算人群花生黄曲霉毒素B1暴露量，然后通过评估模型2，将人群中乙型肝炎病毒携带率P、花生黄曲霉毒素B1对乙型肝炎病毒携带者的平均致癌强度M和花生黄曲霉毒素B1对非乙型肝炎病毒携带者的平均致癌强度N固定参数以及由评估模型1计算获得的人群花生黄曲霉毒素B1暴露量输入评估模型2，进行多次重复抽样模拟运算，得到各Bootstrap小样本的花生黄曲霉毒素B1暴露量概率分布，并得到各Bootstrap小样本的原发性肝癌年发生率即原发性肝癌年超额病例数结果，然后进行统计分析，分别得到待评估的区域人群人体花生黄曲霉毒素B1暴露分布中相应统计量的概率分布和花生黄曲霉毒素B1污染对特定区域人群人体健康风险即原发性肝癌年超额病例数的不确定性评估结果。