1、纯净水安全监控问题何挺 严茂花 叶宁内容摘要水是人类赖以生存的资源，日趋加剧的水污染，已对人类的生存安全构成重大威胁，成为人类健康、经济和社会可持续发展的重大障碍，对水污染的研究，在现实社会中有着重大的的意义。本文是针对九家纯净水公司的具体检查情况建立模型对九家纯净水公司进行综合评价。模型 I:针对问题 1,本文建立了综合评价模型。首先，对题中所给的数据进行无量纲化处理, 建立纯净水风险评价指标体系，采用专家给出的各个项目（危害因素）的权值 分别为 0.17,0.25,0.28,0.30,最后给出纯净水1234,w安全风险分析的科学评价方法，即建立综合评价模型 12340.7.50.8.yxx
2、并对城区所有批次的纯净水进行了评价排序。 针对问题 2，本文在模型 I 的基础上进行求解，由九家公司每个项目（危害因素）的均值得到各公司纯净水中四个项目的得分权重，即各公司的风险度。风险度排序评价为： DHBIACFEG公 司 公 司 公 司 公 司 公 司 公 司 公 司 公 司 公 司然后对他们分类综合评价，指出了各公司产品的主要可能的危害因素，并给出了同类公司的实际特点。模型 II：针对问题 3，本文建立了熵权模型，通过熵值法求出四个危害因素在不同环节的权重，据此分析得出纯净水在生产流通环节的危害因素以及各个危害指标的分布规律，并在此基础上，讨论得出 公司的管理状况。ABD、 、模型 I
3、II:针对问题 4,我们在现有的历史数据下,求出各个危害因素所占的百分比,据此得到每家公司在下一年度 100 个批次的检验总数中的检验分配数,在此基础上,我们建立以检测出不合格品的风险度最高为目标函数的整数线性规划模型,利用 求解,可得到各公司检验批次的分布。LINGO针对问题 5，结合上述结果，本文评价了该城区的饮用水安全形势并给出了监控对策（见正文） 。关键词：综合评价 危害因素 熵权模型 熵值法 线形规划 软件LINGO一 问题重述1.1 问题背景水是人类赖以生存的资源，日趋加剧的水污染，已对人类的生存安全构成重大威胁，成为人类健康、经济和社会可持续发展的重大障碍。据世界权威机构调查，在
4、发展中国家，各类疾病有 8%是由于饮用了不卫生的水而传播的，每年因饮用不卫生水至少造成全球 2000 万人死亡，因此，水污染被称作“世界头号杀手“。纯净水的安全危机的爆发，往往是日常的监控机制和管理长期存在漏洞的反映。完整、有效的纯净水安全风险分析监测预控，为政府及有关部门实施控制措施提供决策依据和技术支持，可以有效提高纯净水安全监管效率和管理水平，及时化解可能出现的安全危机。近年来，我国在从国家宏观层面探讨建立纯净水安全预警机制的研究方面，已取得了不少理论成果但由于我国地域辽阔，经济社会发展水平很不平衡，如何构建有效的预警机制并应用到饮用水安全监控过程还处于起步阶段。本问题主要考虑纯净水的以
5、下危害因素: （按照危害的严重性依次给出）“电导率”: 是纯净水的特征性指标，反映的是纯净水的纯净程度，以及生产工艺控制的好坏， “电导率”根本达不到国家卫生标准要求，与自来水无异，根本不能算做纯净水。菌落总数: 是指纯净水检样经过处理，在一定条件下培养后所取 1ml（g）检样中所含菌落的总数。它可以作为判定纯净水被污染程度的指标之一。大肠菌群：反映纯净水加工过程中对大便污染程度的一个指标。数值越高证明污染越严重。霉菌：食物霉变后产生，直接引起中毒，或产生致癌物质，毒害人体。 1.2 问题的提出某城区共有九家生产并销售纯净水的公司，其中 A 公司和 B 公司规模较大，其余均为小公司。针对该城区
6、提供的近年的关于各公司的纯净水检测报告（见附件） ，请你利用数学建模的方法回答以下问题：1、结合本问题所给数据，给出纯净水安全风险分析的科学评价方法，确定评价的标准和评价的规则，对该城区所有批次的纯净水进行评判排序。2、对该城区范围内的监控对象（各公司）按风险度进行排序评价，并对它们分类综合评价，指出各公司产品的主要可能的危害因素，并指出同类公司的实际特点。3. 对检测出的不合格的样品成因分析：评价纯净水生产流通环节（归为仓库和销售网点两类）的危害因素以及各个危害指标的分布规律，并通过四类危害指标的分析，讨论 A、B、D 公司的管理状况。 4.国家相关部门每年要面对各种专项检验，对于纯净水专项
7、检验的投入经费有限，已知该城区下一年度投在纯净水方面的检验总批次为 100 个批次，在现有历史数据的基础上，并考虑各公司的实际运行状况，如何设置各公司检验批次的分布，使得抽检方案的针对性最优（即检出的风险性为最大） 。5.结合你的工作，请你给该城区食品安全委员会写一篇短文，阐述你的观点，评价该城区的饮用水安全形势并给出监控对策。二 问题分析2.1 问题 1 的分析首先，对题中所给的数据进行分析，同各个项目的标准值进行比较，得到每家公司各个批次的纯净水中不符合要求的纯净水的批次。其次，建立纯净水风险评价指标体系，由专家打分给出的四个项目的权重。最后，根据四个项目的权重建立综合评价模型，并利用此模
8、型对该城区所有批次的纯净水进行评判排序。2.2 问题 2 的分析问题要求对城区范围内的监控公司按风险度进行排序评价，我们在找出九家公司每个指标的平均值的基础上，利用问题 1 中所建立的综合评价模型（模型 I） ，求出九家公司的纯净水中四个指标加起来的风险度，即为各公司的风险度；接着分别按四个危害因素将九家公司进行排序评价，由此评价，可以看出各公司的主要危害因素。根据主要危害因素，我们指出了同类公司的实际特点。2.3 问题 3 的分析根据题中所给数据找出的不合格的样品所对应的环节，找出四个危害因素在各个环节的权值，因为抽取的样品中，存在一定的不确定度，所以我们建立了熵权模型，利用熵值法解决这一问
9、题，并且给出各个危害指标的分布规律。在此基础上，我们通过分析，讨论得出 公司的管理状况。ABD、 、2.4 问题 4 的分析本问要求我们设置各公司检验批次的分布，使得抽样方案的针对性最优，也即根据历史数据，使得在下一年度抽样检验中各公司生产的纯净水被检验到不合格的风险度最大，由于批次的分配数目不仅仅要落实到各家公司，还要具体到各家公司的两个环节（仓库和销售网点） ，所以我们先根据现有的历史数据，在算出各个危害因素所占的百分比的前提下，确定分配到各家公司的检验批次，再以检测出不合格品的风险度最高为目标函数，建立线形整数规划，在软件中求解，得到每家公司在各个环节的检验次数，也即各公司检验批LING
10、O次的分布。三 模型假设1.假设项目 4（霉菌和酵菌）中，标准值不得检出记为 0 值；2.各批次有一个危害因素不合格则视为该批次不合格；3.假设同一批次的纯净水各项检测指标都一样； 4.在一年内没有被检查，则视为该批次在这一年为合格；5.所有的数据都客观准确；6.不考虑四种因素以外的因素的影响；7.在一定时期内，各公司的管理方式不变。四 符号约定ixi i： 表 示 第 个 项 目 （ 危 害 因 素 ） （ =1.4） ；, i： 标 准 化 后 的 数 据 ；:iiMxm表 示 的 最 大 值 ；表 示 的 最 小 值 ；iw各 个 项 目 （ 危 害 因 素 ） 的 权 重 (i=1,2
11、34);:i， 第 个 危 害 因 素 在 各 个 环 节 的 权 重 ；9ipi第 公 司 的 不 合 格 数 （ ， ） ;: ;ijyj第 个 危 害 因 素 下 第 个 时 间 检 测 样 品 数 量ijz第 个 危 害 因 素 下 第 个 时 间 检 测 的 值 ；:Q抽 样 方 案 风 险 度 ；ijqj： 第 个 公 司 在 第 个 环 节 的 风 险 度 ；:ijR第 个 公 司 在 第 个 环 节 抽 样 分 配 数 ；iN： 第 个 公 司 的 抽 样 分 配 个 数 ；jS： 第 个 环 节 抽 样 分 配 数 。五 模型的建立与求解5.1 模型 I 的建立与求解（针对问
12、题 1）5.1.1 模型 I 的建针对问题 1，本文建立了综合评价模型。假设这九个公司为评价对象，共有四个评价指标ABCDEFGHI、 、 、 、 、 、 、 、（项目）电导率、菌落总数、大肠菌群、霉菌和酵母，分别记为 和 。123x、 、 4x首先：对四个指标进行无量纲化处理，也即标准化处理：因为这四个指标都是极小型指标，所以对数据直接进行极差变换将数据标准化，即 ，各个指标对应的数据变为如下表：iiixmM，表 1：标准化后的数据批号 项目 1(电导率0);data:N=22 24 4 32 1 2 1 9 5;x= 0.1505073 0.0597484 0.0026983 0 0.0022905 00.0496635 0.1525125 0.0306815 0.0043798 0 00 0 0 0 0.00423814 00.154619 0.1772801 0 0 0 00 0 0 0 0 0.00193110 0 0 0 0.0218883 0.00117150.0460465 0.0284019 0 0 0 00 0 0.0186512 0 0 00.0550745 0 0 0 0 0;q=45 42 5 1 6 1;enddataend
