1、层次分析法在水环境质量评价的方法X学院毕业论文中文题目:层次分析法在水环境质量评价的应用学生姓名丄 系 别赵 专业班级丄 指导教师丄 成绩评定 2011 年6月1引言 12水环境评价国内外研究现状及发展趋势 12.1水环境质量评价国内外研究现状 12.2水环境质量评价发展趋势 43水环境质量评价应用 53.1玄武湖水质概况 53.2层次分析法简介及基本步骤 53.3玄武湖质量评价计算 53.3.1西北湖质量评价计算 63.3.2东南湖质量评价计算 93.3.3东北湖质量评价计算 113.3.4西南湖质量评价计算 133.4计算结果 154 结论 15参考文献 17致谢 18附表 19摘要玄武湖
2、是南京城区内最大的湖泊，分为西北，东南，西南和东北四个湖区。 随着人口的不断增长， 大量的生活污水被排入湖中， 使得湖水逐年水质污染严重， 近年来大有加重趋势。水质将变浑发臭，严重影响了城市环境和市民身体健康。 因此，对玄武湖整治是势在必行的。本次研究是根据玄武湖实测水质指标资料， 应用层次分析法， 通过建立层次分析模型体系， 经计算最终得出各指标权重， 从 而对玄武湖四个湖区水资源质量状况进行相关分析与评价， 并提出保护水资源措 施，为更好的开发和利用水资源提供参考。此次研究结果表明，湖水中 BOD超标严重，应采取针对性治理。层次分析法在玄武湖水环境评价的应用X1引言玄武湖位于南京城东北，为
3、全国五大城市湖泊之一。湖泊总面积 4.4 km 2 ，其中湖水面积 3.91 km ，湖中陆地面积 0.49 km 。从 80 年代开始， 玄武湖流域 人口迅速增加， 大量的生活污水排入湖内， 使湖水中的有机物， 营养盐的负荷不 断增加。 2003 年的水质监测报告表明：玄武湖水质的总氮超标率为 8.3%。另一 项针对玄武湖底层泥巴的调查实验表明， 表面上一汪清水的玄武湖， 底层的那些 淤泥重金属含量，已经超过了南京城市的土壤的重金属含量 1 。本研究的目的是 运用层次分析法来研究玄武湖水环境问题， 通过分析计算确定其污染程度及主要 污染物，并对玄武湖水质进行相关评价， 以期为提高和改善玄武湖
4、的水环境提供 有效途径和理论依据。 水环境评价能为决策者提供有效的辅助决策信息， 对于水 环境保护和实现可持续发展具有重要意义。 其研究成果将大大提高玄武湖水环境 管理与决策的水平， 对于更好地利用和保护水资源、 控制水污染， 都具有重要的 现实意义和深远的历史意义。2国内外研究现状及发展趋势2.1国内外水环境评价研究现状国内外环境质量评价方法多种多样， 但目前国内还没有制定出统一的评价方 法标准供环保工作者使用 2 。其中水环境质量评价方法较多如：布朗水质指数、 普拉特水质指数、 罗斯水质指数、 内梅罗水质指数、 综合污染指数、 模糊数学法 和地图叠加法等， 最后一种方法是国内目前普遍采用的
5、方法， 简单且实用。 以上 种种评价方法都要首先确定断面单项指标代表值， 大多用平均值作为代表值， 而 内梅罗水质指数法则既考虑到平均值， 同时考虑端值对评价结果的影响， 但评价 工作中有很多具体问题不易解决， 很少采用。 其他方法还牵涉到标准问题或评价 指标的权重问题，也很少采用 3 。水质评价方法有两大类， 一类是以水质的物理化学参数的实测值为依据的评 价方法；另一类是以水生物种群与水质的关系为依据的生物学评价方法。 较多采 用的是物理化学参数评价方法， 其中又分： 单项参数评价法即用某一参数的实 测浓度代表值与水质标准对比， 判断水质的优劣或适用程度。 多项参数综合评 价法即把选用的若干
6、参数综合成一个概括的指数来评价水质，又称指数评价法。指数评价法用两种指数即参数权重评分叠加型指数和参数相对质量叠加型指数 两种。 参数权重评分叠加型指数的计算方法是， 选定若干评价参数， 按各项参数 对水质影响的程度定出权系数，然后将各参数分成若干等级，按质量优劣评分， 最后将各参数的评分相加， 求出综合水质指数。 数值大表示水质好， 数值小表示 水质差。用这种指数表示水质， 方法简明，计算方便。 参数相对质量叠加型指数 的计算方法是， 选定若干评价参数， 把各参数的实际浓度与其相应的评价标准浓 度相比， 求出各参数的相对质量指数， 然后求总和值。 根据生物与环境条件相适 应的原理建立起来的生
7、物学评价方法，通过观测水生物的受害症状或种群组成， 可以反映出水环境质量的综合状况， 因而既可对水环境质量作回顾评价， 又可对 拟建工程的生态效应作影响评价， 是物理化学参数评价方法的补充。 缺点是难确 定水污染物的性质和含量。在选择城市水质评价指标时，一般可按下列几个原则进行选择：（1）据水环境监测和污染源监测所掌握的水中主要污染物种类、数量以及 迁移、转化规律及危害大小进行选择。（2）选择指标应有相应的水质标准及其它明确的规定。（3）所选择指标标应有明确的物理意义，能够表征所受污染的类型。（4）应注意到各指标之间的可比性和指标的概括面及代表性，应尽可能选 择较少的指标用以反映水体污染的特征
8、。以下列举几种国内外较为常用的水质分析方法：水质指数法是用一个综合指数表示水质的物理和化学参数， 特点是便于比较。 分类评价 法是根据水质分类标准， 先确定单项因子所属的等级。 然后进行叠加并综合评价， 特点是参数选择范围广。 数理统计法是把每一次监测数据看成一个随机变量， 利 用数理统计的方法推求出各种概率下的水质状况， 特点是能较完整地提供监测水 体质量的时空变化。 模糊数学法是根据水体质量好与差的界面不明确， 具有模糊 性特点， 运用模糊数学中的隶属函数剖析水质分类界限， 综合评判水体质量。 灰 色系统法是利用水质信息系统中具有信息不完备的特点， 把灰色系统理论中的方 法应用于水质评价。
9、模糊数学法水体污染及其轻重的程度属于模糊概念， 而根据水质标准进行水质评价是一 个典型的模糊模式识别问题， 所以采用模糊集理论方法， 能使水质评价的理论与 方法建立在较为严谨的数学模型上， 能较为准确的反映出水质的状况。 采用模糊 向量加权计分法， 是用隶属度函数来描述水质分级的界限， 计算权重， 并进行模 糊矩阵的复合运算，取隶属度最大者所对应的水质级别作为评价水体的评价级 别，最后以模糊向量为权， 乘以相应级别的分值， 计算出水质的分值。 这种方法 能够较为直观准确的表达出目前水质情况， 因此采用模糊向量加权计分法做为评 价的方法。层次分析法层次分析法（AHp是将决策总是有关的兀素分解成目
10、标、准则、方案等层 次，在此基础之上进行定性和定量分析的决策方法。 该方法是美国运筹学家匹茨 堡大学教授萨蒂于本世纪 70 年代初，在为美国国防部研究 根据各个工业部门对 国家福利的贡献大小而进行电力分配 课题时，应用网络系统理论和多目标综合 评价方法，提出的一种层次权重决策分析方法。 这种方法的特点是在对复杂的决 策问题的本质、 影响因素及其内在关系等进行深入分析的基础上， 利用较少的定 量信息使决策的思维过程数学化， 从而为多目标、 多准则或无结构特性的复杂决 策问题提供简便的决策方法。尤其适合于对决策结果难于直接准确计量的场合 4 。层次分析法的步骤如下：（1） 通过对系统深刻认识，确定
11、该系统的总目标，弄清规划决策所涉及的 范围、所要采取的措施方案和政策、实现目标的准则、策略和各种约束条件等， 广泛收集信息。（2） 建立一个多层次的递阶结构，按目标的不同、实现功能的差异，将系 统分为几个等级层次。（3） 确定以上递阶结构中相邻层次兀素间相关程度。通过构造两比较判断 矩阵及矩阵运算的数学方法， 确定对于上一层次的某个兀素而言， 本层次中与其 相关兀素的重要性排序 - 相对权值。（4） 计算各层兀素对系统目标的合成权重，进行总排序，以确定递阶结构 图中最底层各个兀素的总目标中的重要程度。（5） 根据分析计算结果，考虑相应的决策 5。层次分析法整个过程体现了人的决策思维的基本特征，
12、 即分解、判断与综合， 易学易用，而且定性与定量相结合，便于决策者之间彼此沟通，是一种有你 效 的系统分析方法， 广泛地应用在经济管理规划、 能源开发利用与资源分析、 城市 产业规划、人才预测、交通运输、水资源分析利用等方面 6 。层次分析法也是本文所采用的分析方法。2.2国内外水环境评价发展趋势水环境质量评价是认识和研究水环境的一项重要内容， 其目的是准确反映环 境的质量和污染状况， 预测未来的发展趋势， 是水环境管理保护和治理的一项重 要基础性工作。目前在进行水环境质量评价时，主要存在以下问题：(1) 迄今没有一个被大家公认通用的具有可比性的水环境质量评价数学模 型，各部门进行评价时， 选
13、用数学模型的任意性很大， 常常不能反映本地区污染(2) 多因子综合评价中确定环境因子权重存在任意性， 缺乏比较客观可靠的 确定环境因子权重的方法 7因此，未来的水质评价模型需要解决好上述两项问题。 目前科研人员研究出了一种新型的水质评价模型人工神经网络( ANN)，该方法是复杂非线性科学和人工智能科学的前沿， 其在水质评价的应用研究在国 内外尚处于初创阶段。在详细分析LM算法的基础上，提出了基于LM算法的水质综 合评价BP莫型，并将该模型应用于实例，进行效果检验。结果表明： LM-BP莫型用于水质综合评价是可行的，与其它评价方法相比，评价结果更加客观、合理，此外，随着科学技术的飞速发展， 各种
14、高科技手段也将应用到水质评价中来。 例如目前计算机已经被广泛应用到水环境评价中， 在全国水质数据集中管理实现 后，可以开发基于GIS和WE的水质信息管理系统，能够对海量数据进行直观、高 效的分析与访问，为水环境管理提供辅助工具， 为水资源评价提供可靠数据支持。 另外随着自动化技术的应用，大大提高了工作效率并降低了各种人为误差的可 能，也为水资源评价的准确性提供了支持。 目前还可以利用数学莫型莫拟氮磷营 养物的流失过程， 是研究营养物污染来源和扩散途径的有效手段， 从而能够有针 对性地选择流域中的部分污染区进行控制。 莫拟氮磷流失的数学莫型有负荷莫型 和迁移转化模型两类。负荷模型中，CREAM用
15、于研究土地管理对营养物和杀虫剂 的影响；GLEAM用于模拟地下水中的营养物负荷；EPIC用来计算侵蚀对农作物产 量的影响；AGNP用于计算氮磷的流失量；SWRF可用于决定农作物轮作、种植和 收获日期、化肥的施用日期和用量等管理因子的影响； Ketelsen 的水文地质模型， 可用于解释岛上含水层氮的空间分布 8； Qian 的组合模型可预测营养物的截获 量9 ；土壤测试值法是利用磷的吸附与解析动力学模拟计算磷的土壤测试值，可 用于指示污染风险。输移模型中，NTT-Watershed模型用于模拟氮硝态氮、铵态氮、有机氮的输移和最后的存在形式， 可供集水区规划人员开发最佳的管理措施 方案，以便最大
16、程度地减少氮的输出量，从而减少氮的污染； LPM-N于评价从农业用地渗出的氮，适合于管理肥料的应用和估计进入地下水的氮的渗透浓度 10 ； 土壤水评价工具SWA能帮助水资源管理者评价营养物等非点源污染和管理措 施，进行规划和决策。因此未来的水资源评价与高科技的结合是必然趋势。3水环境质量评价应用3.1玄武湖水质概况本文所要实例分析的玄武湖是城市浅水型湖泊， 随着南京城工业的发展和人 口的增长，这种类型的湖泊十分容易受到大量污染物的污染， 而水体遭受污染后， 一些有毒的重金属最容易沉淀在底部的淤泥里， 释放有害健康的物质。 近年来的 玄武湖富营养化严重， 蓝藻等各类水藻不断生长， 研究表明， 大
17、多数蓝藻会使水 质变腥变臭， 许多种的藻类含有毒素， 可富集在水生生物体内， 通过食物链影响 人类健康。长期下去会使湖泊变浅，甚至成为沼泽。此次通过层次分析法的应用， 通过质量评价计算分析出玄武湖的水质状况， 并 通过权重直观分析出主要污染物，从而找到合理的方法对玄武湖水进行综合整 治，改善湖水质量。3.2层次分析法简介及基本步骤层次分析法简介：层次分析法本质上是一种决策思维方式， 它把复杂的问题分解为各个组成因 素，将这些因素按支配关系分组形成有序的递阶层次结构， 通过两两比较的方式 确定层次中诸因素的相对重要性， 然后综合人的判断以决定决策诸因素相对重要 性总的顺序。APH现了人们决策思维
18、的这些基本特征，即分解、判断和综合。 层次分析法基本步骤：层次分析法经常可以作为一种确定指标权重的方法加以应用。在项目评估 中，运用层次分析法确定指标权重一般分为五个步骤： （1）建立问题的递阶层次 结构； （2）构造两两比较判断矩阵； （3）进行层次因素（指标）单排序数值的计 算；（4）进行APH判断矩阵一致性检验；（5）计算各层指标的组合权重。本次计算将通过上述几个基本步骤对玄武湖水质指标进行相关计算并分析。3.3玄武湖质量评价计算玄武湖分为四个湖区： 西北湖， 西南湖， 东南湖和东北湖， 本次研究将通过 实际的水质评价计算，对玄武湖水质质量进行相关分析。1建立问题的递阶层次结构，其结构如
19、图建立问题的递阶层次结构，其结构如图1:水资源评价指标体系 A2构造判断矩阵根据有关部门提供的玄武湖的四个湖区水质评价监测的主要6种因素 2001-2002年的平均值，将其列入下列矩阵中，即得实测样本矩阵，其单位为mg/L。 另外，为了研究玄武湖整体水质情况，把求得的各测点监测因素的平均值也列入 样本矩阵中，以反映所有测点总体平均状况。玄武湖水体水质监测数据见表 1.表1玄武湖水体水质监测数据表 单位：mg/LNH3- NDOBOD5Cr6+挥发酚石油类西北湖0.1210.05.80.0020.0010.18东南湖0.289.48.90.0020.0020.18东北湖0.1510.35.20.
20、0020.0010.17西南湖0.219.47.00.0020.0010.18全湖0.209.67.40.0020.0010.18通过查阅表1结合地表水环境质量标准基本项目标准限值 (附表1)，计算B-C层六项水质指标的各自权重值。3.3.1西北湖质量评价计算水质质量评价计算(DO,BOD5氨氮，Cr6+，挥发酚，石油类)，构造D(权重 如表3所示：表2 构造DO权重判断矩阵C1C2C3C4C5C111.622.832C20.62511.251.752C30.50.811.41.6C40.3750.5710.71411.143C50.3130.50.6280.8751根据表 2 计算,得到：
21、W(G:Bi)=0.358,W(C2：Bi)=0.224,W(C3：Bi)=0.179,W(C:Bi)=0.128,W(C5：Bi)=0.112。.判断矩阵一致性检验，并计算矩阵最大特征值。由公式（1）得，max =6 ,由公式（2）得，ci =0。所以CR =CI R =0 :：: 0.1。满足一致性，矩阵满足要求。表3平均随即致性指标（RI)n123 456789RI000.52 0.891.121.261.361.411.46表4构造BOD权重判断矩阵C1C2C3C4C5C1110.6430.0711.5C2110.6430.0711.5C31.5561.55610.1112.333C4
22、14149121C50.6670.6670.4290.0481跟据表 4计算，得到： W（C:B2）=0.055,W（C 2：B2）=0.055,W（C 3：B2）=0.085, W（C:B2）=0.768,W（C 5：B2）=0.037。进行判断矩阵一致性检验，CR = CI . RI = 0 : 0.1 ,满足一致性，矩阵满足要 求。表5构造氨氮权重判断矩阵cGGC4GG112.66729.3334662.667Cb0.07912.3163.6324.947G0.0340.43211.5682.316G0.0220.2750.63811.362C50.0160.2020.4680.7341
23、根据表 5计算，得到： W（CB3）=0.869，W（C:B3）=0.069，W（B3）=0.030，W(C:B3)=0.019，W(C:B3)=0.014。进行判断矩阵一致性检验，CR = CI . RI = 0 - 0.1 ,满足一致性，矩阵满足要求。表6构造Cr6+权重判断矩阵CiC2C3C4C5Ci166612.25C20.1671112.042C30.1671112.042C40.1671 COD112.042C50.0820.4900.4900.4901根据表6计算，得到：W(CB4)=0.633，W(G:B4)=0.105，W(GB4)=0.105，W(C:B4)=0.105，W
24、(B4)=0.052。进行判断矩阵致性检验，CR = CI RI =:0 ：: 0.1 ,满足致性,矩阵满足要求。表7构造挥发酚权重判断矩阵C1C2C3C4C5C1114949C2114949C30.250.2512.2524.75C40.1110.1110.444111C50.0100.0100.0400.0901根据表7计算，得到：W(CB5)=0.422，W(c:B5)=0.422，W(B5)=0.105，W(c:B5)=0.047，W(B5)=0.004。进行判断矩阵致性检验，CR = Cl RI =0 ：: 0.1 ,满足致性,矩阵满足要求。表8构造石油类权重判断矩阵C1C2C3C4
25、C5C11112.4626.308C21112.4626.308C31112.4626.308C40.4060.4060.40612.563C50.1590.1590.1590.3901根据表8计算，得到：W(CB6)=0.281，W(G:B6)=0.281 ,W(GB6)=0.281，W(C:B6)=0.114 , W(B6)=0.045。进行判断矩阵一致性检验，CR =CI RI =0 ：0.1,满足一致性，矩阵满足要求。计算 B-C层组合权重，W=2.415，W=3.677, W=1.242，W=1.838，W=1.789，W=2.367。表9 建立B层元素之间相对重要度比较表示两个因素
26、相比，具有相同重要性3表示两个因素相比，前者比后者稍重要5表示两个因素相比，前者比后者明显重要7表示两个因素相比，前者比后者强烈重要9表示两个因素相比，前者比后者极端重要2，4，6, 8表示上述相邻判断的中间值根据表10构造A-B层6项水质指标重要度判断矩阵表10构造A-B层重要度判断矩阵B1B2B3B4B5B6B110.1430.250.3330.3330.5B2711.752.3332.3333.5B340.57111.3331.3332B430.4290.75111.5B530.4290.75111.5B650.7140.50.6670.6671根据表 10计算，得到： W(B:A)=0
27、.050 , W(B:A)=0.350 , W(B:A)=0.200 ,W(B:A)=0.150，W(B:A)=0.150，W(B:A)=0.100。进行判断矩阵一致性检验，CR二ClRI = 0 ： 0.1 ,满足一致性，矩阵满足要 求。计算西北湖水质体系水质指标综合权重W=2.43 &3.3.2东南湖质量评价计算水质质量评价计算(DO,BOD氨氮，C6+,挥发酚，石油类)表11构造DC权重判断矩阵C1 C2C3C4C5C11 1.7892.3163.3683.895C20.559 11.2941.8822.176C30.432 0.72311.4551.682C40.297 0.5310.
28、68811.156C50.257 0.4590.5950.8651根据表11计算，得到：W(CB1)=0.394，W(G:B1)=0.220，W(G:B1)=0.168，W(CB)=0.117，W(C:B1)=0.101。进行判断矩阵一致性检验，CR =CI . RI=0 ：: 0.1 ,满足一致性，矩阵满足要求。表12构造BOD权重判断矩阵C1C2C3C4C5C1110.8310.4920.186C2110.8310.4920.186C31.2041.20410.5920.224C42.0342.0341.69010.397C55.3645.3644.4552.6361根据表 12计算，得到： W(CB2)=0.094 , W(G:B2)=0.094 , W(G:B2)=0.113 ,W(C:B2)=0.193 , W(B 2)=0.505。进行判断矩阵一致性检验，CR.CI .RI 0.1 ,满足一致性，矩阵满足要求。表13构造氨氮权重判断矩阵C1C2C3C4C5C111.6925.5389.38513.231C20
