水质评价问题(共24页)

1、精选优质文档-倾情为你奉上水质评价问题摘要本文主要讨论水质的综合评价问题，根据所给材料，结合所建立的数学模型及所给出的评价标准，综合比较了四种水源之间的差异，采用层次分析法对某村各相距500米以上的四口水井水质进行排序，并依据所得出的结论对四种井水进行水质等级的评价。对于问题一，采用层次分析法对四口井水质进行综合评价。首先建立递阶层次结构模型，构建出单层次中的所有矩阵和总层次矩阵，然后进行层次单排序及一致性检验和层次总排序及一致性检验，结果表明矩阵的排序具有较满意的一致性，可以接受该分析结果。四口水井水质由好到坏的排序为：东水井、南水井、北水井、西水井。对于问题二，运用多元分析法中的系统聚类分

2、析法，充分利用了所给的各项数据，依据四口井水的检测数据与各项指标标准化的差值做为矩阵的元素，形成距离矩阵，求出较为准确的结论，画出矩阵图，可得出结论：东井、南井、北井属于第类水质，西井属于第类水质。关键词：层次分析法、多元分析、系统聚类分析法。一、问题重述1. 问题背景水是地球上一切生物赖以生存也是人类生产生活不可缺少的最基本物质。不同用途的水质要求有不同的质量标准。我国规定，根据供水目的的不同，存在着、农用灌溉水水质标准等。各种工业生产对水质要求的标准也各不相同，依照地表水环境质量标准（GB3838-2002）中规定，地面水使用目的和保护目标，中国地面水分五大类：类水质，水质良好。地下水只需

3、消毒处理，地表水经简易净化处理（如过滤）、消毒后即可供生活饮用者。 类水质，水质受轻度污染。经常规净化处理（如絮凝、沉淀、过滤、消毒等），其水质即可供生活饮用者。 类水质，适用于集中式生活饮用水源地二级保护区、一般鱼类保护区及游泳区。 类水质，适用于一般工业保护区及人体非直接接触的娱乐用水区。 类水质，适用于农业用水区及一般景观要求水域。超过五类水质标准的水体基本上已无使用功能。2. 问题重述某村内有各相距500米以上的四口水井，分别位于村东、村西、村南和村北，由于农业和生活排放废物使地下浅表水遇到污染，水质监测资料如附录1所示. 请完成以下问题：（1）请用2种以上的数学方法对该村的四个井水的

4、水质进行排序，并比较是否由于方法的不同导致存在着异，以及差异产生的原因。（只用了一种方法做）（2）请对该村的四个井的地表水分别进行水质等级判断。（水质分级标准参考附录2，或自己查有关资料）二、问题分析该村四口水井是村民饮用水，根据以上提到的地表水环境质量标准（GB3838-2002）的标准，必须符合 、类水质。因此，只需考虑前三类水质。对于问题一，该问题要求我们对该村的四个井水水质进行排序，可根据层次分析法分析数据得到若干矩阵，进行整理综合，比较四口井水的水质，对四口井进行水质质量排序。对于问题二，综合利用所给出的四口井水的检测数据与附录2中所给的数据进行比较，多元分析建立数学模型。首先观察其

5、所给附录1中监测数据，与附录2水质等级分类中的数据相对照，附录1中PH、溶解氧、高锰酸盐指数、化学需氧量等数据均在附录2中有所对应的数据，可以直接对比，而SO42 、Cl- 、铁、锰、铜等数据却并未在附录2中有具体分类指标，于是我组通过网络浏览查阅了各种资料，根据地表水环境质量标准（GB 38382002）可知集中式生活饮用水地表水源地补充项目标准限值 单位：mg/L项 目标 准 极 限硫酸盐（以SO42-计） 250氧化物（以Cl-计） 250硝酸盐（以N计） 10铁 0.3锰 0.1利用以上所查标准对SO42 、Cl- 、铁、锰、铜等数据进行标准化处理，对四口井监测以上项目进行数据比较，作

6、出距离矩阵，对其进行求和，即可对四口井水水质等级进行等级分类。3、 模型假设1. 在所研究的时间内，四口井水的自然净化能力相似。2. 在所研究的时间内，不会出现较大的降水或干旱，四口井水的水质不会相互交换。3. 在所研究的时间内，四口井水被污染的程度是相似的，且周围的环境不会发生较大的变换。4、 符号说明O目标层：选择井水C准则层：所检测的各项指标P方案层：东井水、南井水、西井水、北井水 Ci准则层的i个指标（0i19）A=( Cij )19*19O-C之间的成对比较矩阵，即为判断矩阵P=(Pi)4*4P-C之间的矩阵max判断矩阵A的最大特征根W判断矩阵A的最大特征根max的特征向量CI 一

7、致性指标RI平均随机一致性指标CR 一致性比例样本集×-R+的一个距离函数5、 模型的建立与求解1. 问题一1.1数据的处理综合评价污染物对水质的影响数据标准化因为评价指标通常是一个介于01 之间的数，所以将四个项目这两年多的数据全部进行标准化处理。 （1）PH 依据附件 2 水质分级标准可知、级水的 PH 均为 6-9，对其进行标准化，取其平均值为 7.5 为其标准值。 （2）因类水质是比较好的饮用水，故其他数据均以类水作为标准对其进行标准化。 1.2建立结构模型O目标层井水水质排序C准则层所检测的各项指标P方案层东井水、南井水、西井水、北井水井水水质排序所检测的各项指标东井水北井

8、水西井水南井水建立井水优劣排序模型。现要比较19个准则层的因子Ci 对O的影响大小，故我们采用每次取两个因子Ci 和Cj ，以Aij表示Ci 和Cj对O的影响大小之比，全部比较结果以矩阵A=( Cij )19*19（见附录3）表示，称A为O-C之间的成对比较矩阵。不难看出，若Ci 和Cj对A的影响之比为Cij ，则Cj和Ci 对A的影响之比为。还要对4口井水的水质进行数据处理，及作出P-C之间的矩阵P=(Pi)4*4(）（见附录4）。1.3层次单排序及一致性检验判断矩阵A对应于最大特征值max的特征向量W，经归一化后即为同一层次相应因素对于上一层次某因素相对重要性的排序权值,这一过程即为层次单

9、排序。计算一致性指标。查找相应的平均随机一致性指标RI对n=1到9RI的值，如下表n123456789RI000.580.901.121.241.321.411.45用随机方法构造500个样本矩阵:随机地从1-9及其倒数中抽取数字构造正互反矩阵,求得最大特征根的平均值。计算平均随机一致性指标。则一致性比列。根据所给的数据可算出：成立，认为判断矩阵的一致性是可以接受的。1.4层次总排序及一致性检验我们最终要得到各元素,特别是最低层中各方案对于目标的排序权重,从而进行方案选择。总排序权重要自上而下地将单准则下的权重进行合成。根据所给的数据可算出：成立，认为层次总排序结果具有较满意的一致性并接受该分

10、析结果。1.5结算结果及结论（matlab（程序见附录5）东南西北总排序权重0.460.27 0.090.18故四口井水的水质排序是：东水井、南水井、北水井、西水井2. 问题二2.1数据处理（1）通过附件2可以发现，溶解氧越大，说明水质越好，所给四口井的溶解氧基本位于58(单位：)之间，并且水质的分类等级标准中只有溶解氧此项取“>”，其余各项均取“”，为了保持相同的比较规则，在标准化溶解氧时只需把溶解氧浓度取倒数即可。（2）SO42-、Cl-、铁、锰、硝酸盐（以N计）依据问题分析中表格对其进行标准化。（3）其他数据处理方式同问题一。2.2建立结构模型求解结果相似性度量是样本集，距离是&#

11、215;R+的一个函数，满足条件如下：（1） d(x,y)>=0， x,y（2） d(x,y)=0，当且仅当x=y（3） d(x,y)=d(x,y), x,y（4） d(x,y)<=d(x,z)+d(z,y), x,y,z 对于相似度的比较，我们选取“欧氏距离”作为标准，求值如下 ：得出距离矩阵如下表：表一 距离矩阵类水质标准类水质标准类水质标准东井水质西井水质南井水质北井水质类水质标准00.18060.98170.04060.19930.06830.1050类水质标准0.180600.80100.22070.05890.16880.1319类水质标准0.98170.801001.

12、02020.85870.96890.9320东井水质0.04060.22071.020200.19650.07960.1275西井水质0.19930.05890.85870.196500.13340.1080南井水质0.06830.16880.96890.07960.133400.0562北井水质0.10500.13190.93200.12750.10800.05620 同时用matlab（程序见附录6）得出聚类分析图，建立井水水质聚类分析模型，分类结果如下：图一 聚类方法图由图可知东井、南井、北井属于第类水质，西井属于第类水质。6、 模型评价1. 优点：在解题过程中，建立模型所应用层次分析法

13、和多元分析法均能对数据进行较为全面的处理，应用了题目中给予的所有数据，使结果较为准确，具有更高的可信度。问题一运用层次分析法，是系统性的分析方法，把研究对象作为一个系统， 按照分解、 比较判断、 综合的思维方式进行决策，使每个层次中的每个因素对结果的影响程度都是量化的，非常清晰、明确。而且层次分析法把定性方法与定量方法有机地结合起来，将复杂的系统分解，能将人们的思维过程数学化、系统化，便于人们接受， 将题中多目标、多准则又难以全部量化处理的决策比较问题化为多层次单目标问题，通过各个因子之间的两两比较确定同一层次元素相对上一层次元素的数量关系后，最后进行简单的数学运算。问题二运用多元分析法，它利

14、用降维技术用少数几个综合变量来代替原始多个变量，这些综合变量集中了原始变量的大部分信息。且它通过计算综合主成分函数得分，对问题进行了科学评价，得出较为准确的结论。并且聚类分析模型的优点就是直观，使结论形式简明。2. 缺点：对问题的处理过程中选择了较为简单的模型，因此对各方面因素的考虑可能不够充分。层次分析法是一种带有模拟人脑的决策方式的方法，因此必然带有较多的定性色彩，会因为标准的不同而导致结果不同，本题我们选用了类的标准计算，也可选用其他类别的标准来计算，这使问题一的结论带有过多的个人色彩。运用多元分析法建立的聚类分析模型并不是应用的题目所给的原始数据，而是对数据进行了整合，使一部分数据产生

15、丢失。3. 附录1某村井水水质监测数据报告编号：商水监/ SM089-2009 监测日期：2011.10.15 样品类别样品编号pH溶解氧总硬度SO42-Cl-铁锰铜锌挥发酚高锰酸盐指数化学需氧量mg/L东井DB8.10 5080.17<DL0.150.00623.851.4西井DB8.346.93271381600.110.19<DL0.150.00316.265.1 南井DB7.496.4524 232 162 0.070.20 <DL0.27<DL1.9 <10北井DB7.157.21020 209 312 1.080.65<

16、DL0.19<DL1.7 <10某村井水水质监测数据报告编号：商水监/ SM089-2009 监测日期：2011.10.15 样品类别样品编号总磷硝酸盐氮亚硝酸盐氮氨氮氟化物氰化物汞砷镉六价铬铅粪大肠菌群mg/L（个/L）东井DB0.781.350.02610.71.48<DL<DL0.0026<DL<DL<DL900西井DB0.982.630.01213.0 1.72<DL<DL0.0178<DL<DL<DL1805南井DB0.15 18.40.3210.050.53<DL<DL0.0025<DL<

17、;DL<DL600北井DB0.03 0.320.004<DL2.18<DL<DL0.0211<DL<DL<DL968附录2序号 分类标准值项目地表水环境质量标准GB3838-2002基本项目标准限值I类II类III类1pH值无量钢6-92溶解氧     （mg/L）饱和率90%（或7.5）653高锰酸盐指数  （mg/L）2464化学需氧量(COD) （mg/L）1515205五日生化需氧量(BOD5)（mg/L）3346氨氮（NH3-N）  （mg/L）0.150.51.07总磷（以P计）&

18、#160; （mg/L）0.02(湖、库0.01)0.1(湖、库0.025)0.2(湖、库0.05)8总氮（湖、库，以N计）（mg/L）0.20.51.09铜       （mg/L）0.011.01.010锌       （mg/L）0.051.01.011氟化物（以F计） （mg/L）1.01.01.012硒       （mg/L）0.010.010.0113砷    &

19、#160;  （mg/L）0.050.050.0514汞       （mg/L）0.000050.000050.000115镉       （mg/L）0.0010.0050.00516铬（六价）（mg/L）0.010.050.0517铅       （mg/L）0.010.010.0518氰化物     （mg/L）0.0050.050.219挥发酚 &

20、#160;   （mg/L）0.0020.0020.00520石油类     （mg/L）0.050.050.0521阴离子表面活性剂（mg/L）0.20.20.222硫化物0.050.10.223粪大肠菌群  （个/L）200200010000附录3pH值溶解氧 （mg/L）高锰酸盐指数 （mg/L）化学需氧量(COD) （mg/L）氨氮（NH3-N） （mg/L）总磷（以P计） （mg/L）总氮（湖、库，以N计）（mg/L）铜 （mg/L）锌 （mg/L）氟化物（以F计） （mg/L）砷 （mg/L）汞 （mg/L）镉

21、（mg/L）铬（六价）（mg/L）铅 （mg/L）氰化物 （mg/L）挥发酚 （mg/L）硫化物粪大肠菌群 （个/L）pH值113.750.55037537.5751507.51507500750750150037501500.0375溶解氧 （mg/L）113.750.55037537.5751507.51507500750750150037501500.0375高锰酸盐指数 （mg/L）0.0.10.13.10010200402404000020002002004001000400.01化学需氧量(COD) （mg/L）227.5110075075150030015300150001500

22、1500300075003000.075氨氮（NH3-N） （mg/L）0.020.020.0750.0117.50.751530.15330001501515307530.00075总磷（以P计） （mg/L）0.0.0.010.0.10.120.40.020.440020224100.40.0001总氮（湖、库，以N计）（mg/L）0.0.0.10.1.1012040.24400020020204010040.001铜 （mg/L）0.0.0.0050.0.0.50.0510.20.010.22001011250.20.00005锌 （mg/L）0.0.0.0250.0.2.50.2551

23、0.05110005055102510.00025氟化物（以F计） （mg/L）0.0.0.50.6.50510020120200001000100100200500200.005砷 （mg/L）0.0.0.0250.0.2.50.25510.05110005055102510.00025汞 （mg/L）6.66667E-066.66667E-060.3.33333E-060.0.00250.000250.0050.0010.000050.00110.050.0050.0050.010.0250.0010.镉 （mg/L）0.0.0.00056.66667E-050.0.050.0050.10

24、.020.0010.022010.10.10.20.50.020.铬（六价）（mg/L）0.0.0.0050.0.0.50.0510.20.010.22001011250.20.00005铅 （mg/L）0.0.0.0050.0.0.50.0510.20.010.22001011250.20.00005氰化物 （mg/L）0.0.0.00250.0.0.250.0250.50.10.0050.110050.50.512.50.10.挥发酚 （mg/L）0.0.0.0010.0.0.10.010.20.040.0020.044020.20.20.410.040.00001硫化物0.0.0.025

25、0.0.2.50.25510.05110005055102510.00025粪大肠菌群 （个/L）26.26.10013.1333.100001000200004000200400020000200004000040001附件4溶解氧 （mg/L）东西南北东1/5.111.1.1.西1/6.90.10.1.南1/6.40.1.11.125北1/7.20.0.0.1高锰酸盐指数 （mg/L）东西南北东23.811.12.14西16.20.18.9.南1.90.0.11.北1.70.0.0.1化学需氧量(COD) （mg/L）东西南北东51.410.3.3.西65.11.14.344.34南150

26、.0.11北150.0.11氨氮（NH3-N） （mg/L）东西南北东10.710.21471.西131.126086.南0.050.0.10.北0.150.0.31总磷（以P计） （mg/L）东西南北东0.7810.5.226西0.981.16.32.南0.150.0.15北0.030.0.0.21锌 （mg/L）东西南北东0.15110.0.西0.15110.0.南0.271.81.811.北0.191.1.0.1氟化物（以F计） （mg/L）东西南北东1.4810.2.0.西1.721.13.0.南0.530.0.10.北2.181.1.4.1砷 （mg/L）东西南北东0.002610.

27、1.040.西0.01786.17.120.南0.00250.0.10.北0.02118.1.8.441挥发酚 （mg/L）东西南北东0.0061233西0.0030.511.51.5南0.0020.0.11北0.0020.0.11粪大肠菌群 （个/L）东西南北东90010.1.50.西18051.13.1.南6000.0.10.北9681.0.1.1SO42-东西南北东18311.0.0.西1380.10.0.南2321.1.11.北2091.1.0.1附录5clc,clearfid=fopen('txt1.txt','r');n1=19;n2=4;a=;fo

28、r i=1:n1 tmp=str2num(fgetl(fid); a=a;tmp; % %读准则层判断矩阵endfor i=1:n1 str1=char('b',int2str(i),'=;'); str2=char('b',int2str(i),'=b',int2str(i),'tmp;'); eval(str1); for j=1:n2 tmp=str2num(fgetl(fid); eval(str2); %读方案层的判断矩阵 endendri=0,0,0.58,0.90,1.12,1.24,1.32,1.4

29、1,1.45; %一致性指标x,y=eig(a);lamda=max(diag(y);num=find(diag(y)=lamda);w0=x(:,num)/sum(x(:,num);cr0=(lamda-n1)/(n1-1)/ri(n1)for i=1:n1 x,y=eig(eval(char('b',int2str(i); lamda=max(diag(y); num=find(diag(y)=lamda); w1(:,i)=x(:,num)/sum(x(:,num); cr1(i)=(lamda-n2)/(n2-1)/ri(n2);endcr1, ts=w1*w0, cr

30、=cr1*w0txt1.txt格式如下：1 1 3.75 0.5 50 375 37.5 75 150 7.5 150 7500 750750 1500 3750 150 0.03751 1 3.75 0.5 50 375 37.5 75 150 7.5 150 7500 750750 1500 3750 150 0.03750. 0. 1 0. 13. 100 10 200 40 2 40 40000 2000 200 200 400 1000 40 0.012 2 7.5 1 100 750 75 1500 300 15 300 15000 1500 1500 3000 7500 300

31、0.0750.02 0.02 0.0750.01 1 7.5 0.75 15 3 0.15 3 3000 150 1515 30 75 3 0.000750. 0. 0.01 0. 0. 1 0.1 2 0.4 0.02 0.4 400 20 2 2 4 10 0.4 0.00010. 0.0.1 0.1. 10 1 20 4 0.2 4 4000 200 2020 40 100 4 0.0010. 0.0.005 0. 0. 0.5 0.05 1 0.2 0.01 0.2 200 10 1 1 2 5 0.2 0.000050. 0.0.025 0. 0. 2.5 0.25 5 1 0.05

32、1 1000 50 5 5 10 25 1 0.000250. 0.0.5 0.6. 505 100 20 1 20 20000 1000 100 100 200 500 20 0.0050. 0.0.025 0. 0. 2.5 0.25 5 1 0.051 1000 50 5 5 10 25 1 0.000250. 0. 0. 0.06 0. 0.00250.000250.005 0.001 0.00005 0.001 1 0.05 0.005 0.005 0.010.0250.0010.0. 0.0.0005 0.5 0. 0.05 0.005 0.1 0.02 0.001 0.02 20

33、 1 0.1 0.1 0.2 0.5 0.02 0.0. 0.0.005 0. 0. 0.5 0.05 1 0.2 0.01 0.2 200 10 1 1 2 5 0.2 0.000050. 0.0.005 0. 0. 0.5 0.05 1 0.2 0.01 0.2 200 10 1 1 2 5 0.2 0.000050. 0.0.0025 0. 0. 0.25 0.025 0.5 0.1 0.005 0.1 100 5 0.5 0.5 1 2.50.1 0.0. 0.0.001 0. 0. 0.1 0.01 0.2 0.040.002 0.04 40 20.2 0.20.4 1 0.04 0

34、.000010. 0.0.025 0. 0. 2.5 0.25 5 1 0.051 1000 50 5 5 10 25 1 0.0002526. 26.100 13. 1333. 10000 1000 20000 4000 200 4000 20000 20000 40000 400011 1. 1. 1.0. 1 0. 1.0. 1. 1 1.1250. 0. 0. 11 1. 12. 14 0. 1 8. 9.0. 0.1 1.0. 0.0. 11 0. 3. 3.1. 1 4.34 4.340. 0.1 10. 0.1 11 0. 214 71.1. 1 260 86.0. 0. 1 0.0. 0.3 11 0. 5.2 261. 1 6. 32.0