苯接触评估-BAYES统计学应用探究.docVIP

苯接触评估:BAYES统计学应用探究摘要： ［目的］ 用BAYES统计学方法估算车间空气中挥发性气体的区域浓度并检验其正确性。 ［方法］ 以某橡胶厂车间空气中苯浓度的历史监测资料为基础，结合专家主观概率及物理模型，进行BAYES统计推断，估算各岗位的区域浓度，并以现场测定数据对评估结果进行检验。 ［结果］ 估算的区域浓度最高的3个岗位为打浆、夹布机成型、搪浆；有7个岗位的区域浓度估算结果与8个岗位的现场采样结果无明显差别，而历史测量资料只有4个岗位的估算结果与现场采样结果无差别。 ［结论］ BAYES统计学方法能降低变异，更有效地评估车间空气中挥发性有机物的区域浓度水平。 关键词：应用；接触评估；BAYES统计学；苯 文章编号:1006-3617(2007)01-0016-05 中图分类号:R195.1 文献标识码:A 接触评估是职业流行病学的重要组成部分，准确的接触评估结果对得到合理解释的剂量-反应关系、采取有效预防措施、保护工人的安全和健康具有重要意义。 目前职业卫生的研究对接触水平大多是定性或半定量的评估，这样不能得到明确的剂量-反应关系。很多定量的接触评估研究也没有标准方法和一致的检验方法，其无偏性、有效性、一致性都有待提高［1］。 为此，本研究以苯为研究对象，用BAYES（非经典）统计学方法将物理模型和专家知识综合起来，估算工作场所苯的区域浓度［2］。 1 对象与方法 1.1 对象 以上海某大型国有橡胶制品企业为研究对象。该厂成立于1954年，现有炼胶、打浆、钢编等共11个生产车间。主要产品为各种橡胶管。新旧车间的自动化程度不同，基本工艺相差不大。生产过程中大量使用苯以增加橡胶的粘合性。车间的自然通风良好，各工房均有轴流风机、吊扇，无可运行的局部排、送风装置。 1.2 方法 1.2.1 现场劳动卫生基本情况调查 了解工厂历史、原料产品、工艺、职业卫生防护、主要工作岗位及其原料及用量、产品及产量、工作任务、操作方式、工作环境、工作频度，以及这些因素在什么时候发生过什么变化。 1.2.2 收集、整理浓度资料 到相关部门收集该厂历年监测资料，统一车间、工种、岗位及采样地点的命名并编码。 1.2.3 车间空气中苯浓度的现场测定 连续10 d对各操作岗位的苯浓度用碳管采样，每天1次。停产岗位则当天不能测定。 1.2.4 计算区域浓度 ①现场参数测定：无各车间的建筑设计、工业卫生设计资料，参数全部自行测定或推算。对打浆车间、夹布机车间的全部参数进行现场测定，包括车间大小、工作点与污染源的距离。测量污染物产生率时，先将车间密闭，车间内均匀布点，同时用12只采样泵在污染物产生前采样1次，在污染物开始产生后连续采样6次，每次10 min，持续2 d。室内污染物前后二次平均浓度之差乘以室内有效容积，除以采样时间，得到污染物产生率。用风速仪持续2 d分别于上午、中午、下午连续2 min测量近区界面、全部通风口风速，计算室内通风量、界面通风量。②收集外推信息：现场测量、询问污染物来源、挥发面积、通风面积、局部通风设施等，供专家推导主观概率。③选择专家：挑选3位从事多年职业卫生工作的专家作为主观概率的提供者。专家各自提出自己的主观概率。如不一致，则分别与作者协商。④编辑参数信息文本：主要内容为本研究的目的、物理模型参考文献、工艺及工作场所照片、模型需要的参数及相关信息、需要专家提供的结果。⑤专家主观概率（先验信息或先验输入）和物理模型：专家根据参数信息文本和专业知识、经验，根据可行性和有效性选择物理模型，给出各参数的数学分布形式、均值及该参数在某可信度下的变动范围。⑥编程计算：WinBugs编程物理模型，设置相应参数，计算区域浓度（后验或后验输出）。 1.3 质量控制 现场调查员均为从事多年职业卫生工作的本专业工作人员，经过多次培训和预调查；数据由双人录入、自动核对；实验室分析由专职人员完成。 1.4 统计分析 ACCESS建立数据库，SPSS11.0作正态性检验、χ2分析。 2 结果 2.1 各岗位历史区域浓度监测资料结果 1964～2003年，该厂有20年区域浓度资料，共有浓度数据515个。算术均数为127.4 mg/m3，几何均数为11.7 mg/m3。区域平均浓度最高的是夹布机成型岗位，浓度为216.85 mg/m3；最低的是东风吸引其他岗位，浓度为0.06 mg/m3。钢编成型岗位浓度最为分散，标准差为39.35 mg/m3，见表1。 2.2 现场各岗位区域采样结果 全厂浓度最高的岗位分别是打浆、夹布机成型和搪浆。除“钢编其他”岗位外，各岗位测量值均为对数正态分布［3］，但离散程度均较高。浓