膳食暴露评估数学模型及模型改进若干方法哈尔滨工业大学祝园园桑俊俊李国华.docVIP

全国第四届研究生数学建模竞赛 题 目 中国膳食暴露评估数学模型及模型改进的若干方法 摘 要： 本文建立由膳食模型、污染物分布模型、风险评估模型组成的中国膳食暴露评估数学模型。首先，参考全国总膳食研究设计合理抽样方案，采用三日称重法加询问法加记录法的调查方法，并根据污染物分布模型对食物进行合理分类，设计调查表，建立基于多元回归的全国膳食模型。其次，综合考虑污染物排放数据、日常食物和进出口监测数据及各自权重，建立基于对数正态分布的污染物分布模型来计算各类食物中各污染物的含量，并用已有的调查数据对模型参数进行估计和验证。然后，利用前两个模型的输出结果，建立基于支持向量机方法的风险评估模型，对全国、某个地区、某类食品的安全状况做出评价并给出预警，而且采用Monte Carlo方法，给出居民某项污染物摄入量的99.999%的右分位点。最后，针对符合性检验只能获得部分样本数据，前两个模型的数据可能不匹配，不同统计分类标准的转化以及利用省市监测数据估计全国情况等问题，提出一些建议和解决方法。 参赛队号 1021313 参赛学校 哈尔滨工业大学 参赛队员姓名 祝园园 桑俊俊 李国华 目录 1 模型假设和符号说明 2 2 人群食物摄入量模型 2 2.1 抽样方法的设计 3 2.2 称重+记录+回顾的膳食调查方法 4 2.3 抽样调查中食物分类办法 6 2.4 人群食品摄入量模型的建立 7 3 污染物分布模型 8 3.1 问题分析 8 3.2 模型假设 10 3.3 模型建立 10 3.4 参数的估计 12 3.5 对部分数据的建模结果 13 4 风险评估模型的建立 15 4.1 风险评估概论 15 4.2 模型的建立和求解 16 4.3 模型分析 20 4.4 模型改进 20 4.5 对右分位点的估计 21 5 改进模型的若干理论研究 21 5.1 污染物分布模型的修正 21 5.2 不配套的抽样调查数据的衔接方法 22 5.3 调查数据中不同统计分类标准不同的转化 22 5.4 部分省市与全国的概率密度函数关系 23 6 结论及后续工作 23 模型假设和符号说明 假设全国总膳食调查数据大部分真实可信，失真数据比重极少，可以忽略。 假设全国膳食调查数据能够真实反映全国平均膳食分布情况。 假设全国膳食调查数据异常数据比例很少，可以剔除。例如某些具有特殊体质的人具有与平常人差异很大的饮食数据。 假设食品卫生检测部门日常对市场上食物的监测数据是真实可信的。 假设食品的生产、销售、流通过程对食品中污染物的分布没有影响。 假设在人群食品摄入模型中的调查对象和污染物分布模型中被调查食品的消费者是在同一个城市，且他们所从事的职业的性质比较接近。这样由于两个模型中的数据不匹配造成的误差较小。 假设在人群食品摄入模型和污染物分布模型中被调查对象的男女比例与该地区或该省的男女比例非常接近。 ：满足一定概率分布的随机变量 ：对数正态分布的对数方差 ：对数正态分布的对数平均值 ：对数正态分布的对数方差的估计 ：对数正态分布的对数平均值估计 ：不同来源食品污染物含量分布的加权值 人群食物摄入量模型 人群食物摄入量模型的主要目的是用于估计不同地区、不同性别、不同年龄、不同季节、不同劳动强度、不同经济收入的人群各类食品的一天摄入量。 要建立人群食物摄入量模型，需获得我国总膳食数据。我国总膳食数据应该由调查人员入户调查获得。但是，由于这项调查工作量太大，如果实行普查，其工作量甚至超过全国人口普查，故而只可能在全国几亿户家庭中随机抽取几千户，至多几万户进行一次性调查。另外，中国居民消费的食品种类比其他国家居民消费的食品种类复杂得多，包括：主食、肉类、蔬菜、水果、水、饮料、各种调味剂和各种经过加工的食品，而主食又将包括大米、粳米、糯米、早釉、晚釉、米粉、面粉、馒头、花卷、挂面、面条、方便面、油条、烧饼、大饼、面包、桃酥、蛋糕、小米、玉米面、玉米棒等各种食品。不难看出，如果以此细分，食物将达数千种以上，在实际调查过程中进行如此详细地分类，其调查工作量太大，而如果随意粗糙进行分类，则将影响调查的精度。 因此，如何根据我们国家的实际情况和污染物分布模型数据，设计合理的抽样调查方案和食物分类办法，使调查结果能尽量反映全国的实际情况，调查结果的数据使用起来效果比较理想，同时使调查的全部工作量在可以承受的范围内，是我们首要解决的问题。 抽样方法的设计 我国是一个拥有13亿人口，3.7