关于研究纯净水安全监控问题
摘 要
水是人类生存的基本需求,生活饮用水的安全问题直接关系到广大人民群众的健康。为研究纯净水安全监控问题,本文建立4个数学模型对其进行实际分析。
模型I 对问题1,要求给出纯净水安全风险的评价方法,本文采用综合评价方法对批次进行评判排序,引进比较尺度进行归一化处理,对数据进行标准化处理,综合评价指标就是两者的乘积。综合评价指标临界值即为标准值为1的情况,其值为0.9375。两者进行比较,超出临界值则说明该批次的纯净水受到污染,反之未受到污染。通过表6.1-1可知有12个批次的纯净水受到了污染。
模型Ⅱ 对问题2,本文把霍顿水质指数法与层次分析法结合,通过建立层次分析模型,构造成对比较矩阵,求解出相应的权重向量,对假定的各个公司的水质风险度进行排序,对各个公司水质的主要危害因素进行分析。7月到9月各公司风险度排序为: > > > > > > > > 。12月到1月各公司的风险度排序为: > > > > > > > > 。7月到9月电导率主要危害者:I H A F;菌落总数主要危害者: B D H;大肠杆菌主要危害者:C F D;霉菌和酵母主要危害者:D。12月到1月A公司的危害因素为电导率,菌落总数;B公司的危害因素为菌落总数;D公司的危害因素为菌落总数;H公司的危害因素为菌落总数;I公司的危害因素为电导率,菌落总数。
模型Ⅲ 对问题3,通过对数据的分析得出各公司在纯净水生产流通环节(归为仓库和销售网点两类)的危害因素以及各个危害指标的分布规律,画出直方图。得到A公司的销售网点的危害因素为菌落总数,仓库的危害因素为电导率、菌落总数、霉菌和酵母;B公司的销售网点的危害因素为菌落总数,仓库无危害因素;D公司的销售网点的危害因素为菌落总数、霉菌和酵母,仓库的危害因素为电导率、菌落总数、霉菌和酵母;H公司的销售网点和仓库的危害因素均为电导率和菌落总数;I公司的销售网点无危害因素,仓库的危害因素为电导率和菌落总数。A、B、D三个公司中,A公司的管理状况最差,B公司的管理状况最好,D公司的管理状况处于两者中间。
模型Ⅳ 对问题4,在模型二的基础上对问题四进行求解,计算出各公司在2010的上半年和下半年各公司的Q值并分别对公司分类,依照比较尺度对四类公司进行较合理的加权并算出各公司所占比例,对各公司检验批次作出比较合理的安排。在2010年按上下半年分各公司检验的批次为:A公司上半年14次,下半年14次;B公司上半年10次,下半年6次;C公司上半年2次,下半年2次;D公司上半年6次,下半年10次;E公司上半年2次,下半年2次;F公司上半年24次,下半年2次;G公司上半年2次,下半年2次;H公司上半年6次,下半年10次;I公司上半年6次,下半年10次;
对问题5,本文通过对以上四个问题的探讨研究,及四个数学模型所得出的结果分析评价该城区的饮用水安全形势并给出相应的监控对策。
关键词:权重 归一化 风险度 层次分析 霍顿水质指数 综合指标临界值
一、问题重述
1.1问题背景
日趋加剧的水污染,已对人类的生存安全构成重大威胁,成为人类健康、经济和社会可持续发展的重大障碍。据世界权威机构调查,在发展中国家,各类疾病有8%是由于饮用了不卫生的水而传播的,每年因饮用不卫生水至少造成全球2000万人死亡,因此,水污染被称作"世界头号杀手"。
我国政府对纯净水安全问题十分重视,已将纯净水安全作为一项重要的公共管理目标,采取了一系列措施,强化纯净水安全的监管,并取得了初步成效。但纯净水安全问题的总体形势仍不容乐观,依然存在一系列隐忧,近年来食品安全方面的恶性、突发性事件屡屡发生。 2007年 07月 12日 ,南通一纯净水厂发生造假事件。2008年3月底,贵阳市发生数百人感染甲肝事件,经卫生部中国疾控中心专家组核查,确认“竹源牌”桶装水是造成疫情爆发的主因。 2009年 03月 25日 ,某大学B区学生饮用了“清清”牌桶装纯净水后,百余学生先后出现集体腹泻事件。 2009年 2月 26日 ,湖南师范某寝室在长沙爱高普纯净水有限公司订购的桶装纯净水中出现了黑色虫子事件。生物性和化学性污染对纯净水安全的影响愈来愈严重。
1.2相关情况
本问题主要考虑纯净水的以下危害因素: (按照危害的严重性依次给出)
“电导率”: 是纯净水的特征性指标,反映的是纯净水的纯净程度,以及生产工艺控制的好坏,“电导率”根本达不到国家卫生标准要求,与自来水无异,根本不能算做纯净水。
菌落总数: 是指纯净水检样经过处理,在一定条件下培养后所取1ml(g)检样中所含菌落的总数。它可以作为判定纯净水被污染程度的指标之一。
大肠菌群:反映纯净水加工过程中对大便污染程度的一个指标。数值越高证明污染越严重。
霉菌:食物霉变后产生,直接引起中毒,或产生致癌物质,毒害人体。
纯净水的安全危机的爆发,往往是日常的监控机制和管理长期存在漏洞的反映。完整、有效的纯净水安全风险分析监测预控,为政府及有关部门实施控制措施提供决策依据和技术支持,可以有效提高纯净水安全监管效率和管理水平,及时化解可能出现的安全危机。近年来,我国在从国家宏观层面探讨建立纯净水安全预警机制的研究方面,已取得了不少理论成果但由于我国地域辽阔,经济社会发展水平很不平衡,如何构建有效的预警机制并应用到饮用水安全监控过程还处于起步阶段。
1.3问题的提出
某城区共有九家生产并销售纯净水的公司,其中A公司和B公司规模较大,其余均为小公司。针对该城区提供的近年的关于各公司的纯净水检测报告(见附件),请你利用数学建模的方法回答以下问题:
1、结合本问题所给数据,给出纯净水安全风险分析的科学评价方法,确定评价的标准和评价的规则,对该城区所有批次的纯净水进行评判排序。
2、对该城区范围内的监控对象(各公司)按风险度进行排序评价,并对它们分类综合评价,指出各公司产品的主要可能的危害因素,并指出同类公司的实际特点。
3. 对检测出的不合格的样品成因分析:评价纯净水生产流通环节(归为仓库和销售网点两类)的危害因素以及各个危害指标的分布规律,并通过四类危害指标的分析,讨论A、B、D公司的管理状况。
4.国家相关部门每年要面对各种专项检验,对于纯净水专项检验的投入经费有限,已知该城区下一年度投在纯净水方面的检验总批次为100个批次,在现有历史数据的基础上,并考虑各公司的实际运行状况,如何设置各公司检验批次的分布,使得抽检方案的针对性最优(即检出的风险性为最大)。
5.结合你的工作,请你给该城区食品安全委员会写一篇短文,阐述你的观点,评价该城区的饮用水安全形势并给出监控对策。
二、问题分析
针对问题1,建立数学模型Ⅰ,首先进行模型的准备,沿用层次分析法中的比较尺度作为综合评价系数;其次,对数据按照统一尺度进行标准化处理;最后,得出综合评价指标为综合评价系数与数据标准化之后结果的乘积。统一标准化区间为[0,1],综合指标大于1说明该公司的纯净水超过安全监控指标,综合指标小于1则说明该公司的纯净水在安全监控指标之内。
针对问题2,本文建立了层次分析模型,应用霍顿水质指数方法对该城区范围内的监控对象(各公司)按风险度进行排序评价,并对它们分类综合评价,指出各公司产品的主要可能的危害因素。利用层次分析法求解各个水质危害指标的权重和各公司对该地区纯净水市场影响的权重,在应用霍顿水质指数表达式: ,得出7月到9月,12月到1月的各公司风险度排序。通过危害率表达式: ,得出7月到9月,12月到1月的各公司主要危害因素。
针对问题3,我们对检测出的不合格的样品进行成因分析,评价纯净水生产流通环节(归为仓库和销售网点两类)的危害因素以及各个危害指标的分布规律,并通过四类危害指标的分析,分析A、B、D公司的管理状况。通过对数据的分析处理,得出各公司销售网点和仓库的危害指标直方图分布。
针对问题4,此问在模型二的基础上对问题四进行求解,计算出相应的 ,并进行归一化处理后可得到各个公司纯净水水质对该城区的风险度,即各公司在2010的上半年和下半年各公司的Q值并分别对公司分类,其次依照比较尺度对四类公司进行较合理的加权并算出各公司所占比例,最后对各公司检验批次作出具体安排。
针对问题5,我们通过对以上四个问题的探讨研究,及四个数学模型所得出的结果分析评价该城区的饮用水安全形势并给出了相应的监控对策。
三、模型假设
1) 霉菌和酵母的标准值“不得检出”用0表示,结果“未检出”也用0表示
2) 不考虑四种因素以外的因素的影响
3) 受检方的综合评价指标与临界值相等说明纯净水未受到污染
4)所有的数据都客观准确
5)各批次有一个危害因素不合格则视为该批次不合格
6)采样地点里的不合格产品是随机分布的
7)在一定时期内,各公司的管理方式不变
8)模型IV中依据温度对四项指标的影响将春季和夏季归为一类考虑,秋季和冬季归为一类考虑
四、符号约定
1. :电导率的结果
2. :菌落的结果
3. :大肠杆菌的结果
4. :霉菌的结果
5. :比较尺度
6. :归一化结果
7. :综合评价指标
8. :综合评价指标临界值
9. :综合评价系数
10. :某项水质指标的实测质量浓度
11. :某项水质指标的标准测量浓度
12. :某项水质指标的权重值
13. :加权水质指数(也称霍顿指数)
14. : 危害率
15. :累计危害率
五、模型的建立
5.1模型Ⅰ
5.1.1 模型准备
当比较两个可能具有不同性质的因素 和 对于一个上层因素 的影响时,采用Saaty等人提出的1-9尺度,即 得取值范围是1,2,…,9及其互反数1,1/2,…,1/9。在进行定性的成对比较时,人们头脑中通常有5种明显的等级,但本文只需考虑4种,用1-9尺度可以方便的表示如下:
尺度 |
含义 |
1 |
与 的影响相同 |
3 |
比 的影响稍强 |
5 |
比 的影响强 |
7 |
比 的影响明显的强 |
2,4,6 |
与 的影响之比在上述两个相邻等级之间 |
1,1/2,…,1/7 |
与 的影响之比为上面 的互反数 |
表5.1-1 1-9尺度 的含义
5.1.2 数据处理
对电导率、菌落、大肠杆菌的结果进行标准化,所得数据如下:
受检方 |
电导率 |
菌落总数 |
大肠杆菌群 |
霉菌酵母 |
||||
标准值 |
结果 |
标准值 |
结果 |
标准值 |
结果 |
标准值 |
结果 |
|
C |
≤1 |
0.26 |
≤1 |
0.05 |
≤1 |
0 |
0 |
0 |
E |
≤1 |
0.26 |
≤1 |
0.05 |
≤1 |
0 |
0 |
0 |
F |
≤1 |
0.391 |
≤1 |
0 |
≤1 |
0.3333333 |
0 |
0 |
C |
≤1 |
0.172 |
≤1 |
0.25 |
≤1 |
0.6666667 |
0 |
0 |
G |
≤1 |
0.1 |
≤1 |
0.15 |
≤1 |
0 |
0 |
0 |
B |
≤1 |
0.314 |
≤1 |
40 |
≤1 |
0.3333333 |
0 |
0 |
A |
≤1 |
0.723 |
≤1 |
1.75 |
≤1 |
0.3333333 |
0 |
0 |
H |
≤1 |
0.802 |
≤1 |
7 |
≤1 |
0 |
0 |
0 |
A |
≤1 |
0.535 |
≤1 |
8 |
≤1 |
0 |
0 |
0 |
D |
≤1 |
0.521 |
≤1 |
9 |
≤1 |
0.3333333 |
0 |
3 |
B |
≤1 |
0.412 |
≤1 |
0 |
≤1 |
0.3333333 |
0 |
0 |
B |
≤1 |
0.389 |
≤1 |
0.05 |
≤1 |
0.3333333 |
0 |
0 |
A |
≤1 |
0.423 |
≤1 |
0.1 |
≤1 |
0 |
0 |
0 |
A |
≤1 |
0.523 |
≤1 |
0.1 |
≤1 |
0.3333333 |
0 |
0 |
A |
≤1 |
0.397 |
≤1 |
0 |
≤1 |
0 |
0 |
0 |
B |
≤1 |
0.553 |
≤1 |
0.05 |
≤1 |
0 |
0 |
0 |
B |
≤1 |
0.394 |
≤1 |
0.2 |
≤1 |
0 |
0 |
0 |
A |
≤1 |
0.337 |
≤1 |
0 |
≤1 |
0 |
0 |
0 |
H |
≤1 |
2.72 |
≤1 |
0.05 |
≤1 |
0 |
0 |
0 |
B |