摘要：介绍了模糊层次法评价水环境质量的基本原理和方法步骤，建立了水环境质量综合评价模型。利用该评价方法就无锡市主要河流上控制断面的水质监测结果进行了具体分析计算，确定水质优劣等级，为后续污染优先防控治理工作奠定基础。
关键词 ：水环境质量 综合评价 模糊 AHP MATLAB
1建立水质综合评价的指标体系
1.1水质综合评价的基本思路是
用层次分析法确定各指标的权重，在此基础上用模糊综合评测理论进行单指标和综合评价。
1.2水质综合评价指标体系的建立
进行水质综合评价首先要设立相应的综合评价指标体系。一般应遵循完全性、简单性、客观性等几项原则。据此建立能全面反映水质状况的综合评价指标体系，归纳如下
表1 水质综合评价指标体系
| 项目 | 内容 |
| PH值（无量纲） | 6～9 |
| 感官性因素 | 色、味、嗅、透明度、浑浊度、悬浮物、总固体 |
| 氧平衡因子类 | 溶解氧、化学耗氧量、生化需氧量 |
| 营养盐因子类 | 盐、氨盐、硫酸盐 |
| 毒物因子类 | 挥发酸、氰化物、汞、铬、砷、镉、铅、有机氯 |
| 微生物因子类 | 大肠杆菌 |
模糊评价即在评价过程中引入模糊性概念，运用模糊数学来处理水质评价的一些问题，以反映水质质量状况或水体污染程度的不确定性；层次分析法是将要评价系统的各种要素进行两两比较，并计算出各要素的权重，根据综合权重按最大权重原则确定最优方案。将这两种方法相结合对水质进行综合评价的方法称为模糊层次综合评价。
2.1建立评价因素集
建立影响评价对象的n个因素组成的集合，称因素集
2.2建立决策集
建立有m个评语结果组成的评价集，亦称决策集
2.3建立隶属函数确定单因素评价矩阵R
应用模糊数学的基本概念，评价因素集中每一个指标隶属于评语集中不同评语分级的程度，称为隶属度。隶属函数一般以uV(x)表示，其含义是：因素集U中的某一指标隶属于模糊集合的程度，其取值在[0 ，1]区间，并以rji表示。此值即为uj的隶属度。评价因素集U中全部指标隶属度的合成，即为单因素评价矩阵。
（1）根据各指标的特征，拟定各隶属函数为线性函数，其数学表达式为：
（1）
将实测的各指标值带入隶属函数，即可求得第j个因素的单因素模糊评价为V上的模糊子集
（2）
于是单因素评价矩阵R为
（3）
2.4用层次分析法建立权重集
2.4.1两两比较构造判断矩阵
根据上面建立的水质综合评价指标体系，通过两两比较，构造判断矩阵
（4）
表示相对于水质状况而言，指标与的比较，采用1～9的标度得到的比重要的程度。
2.4.2计算重要性排序
借助MATLAB软件计算判断矩阵的最大特征值及与之相对应特征向量。
2.4.3一致性检验
判断矩阵的一致性检验指标为
（5）
其中；
——判断矩阵的随机一致性指标，如下
表2 随机一致性指标
| n | 1 | 2 | 3 | 4 | 5 | 6 | 7 |
| RI | 0 | 0 | 0.58 | 0.9 | 1.12 | 1.24 | 1.32 |
| n | 8 | 9 | 10 | 11 | 12 | 13 | 14 |
| RI | 1.41 | 1.45 | 1.49 | 1.51 | 1.54 | 1.56 | 1.58 |
2.4.4模糊层次综合评价
由模糊评判综合原理得：
（6）
对的所得结果进行归一化，并取，即得综合评判结果。
3水环境模糊层次综合评价实例
将上述模型应用于对无锡市周边数条河流的水环境评价。
3.1评价指标及基本数据
对无锡市周边数条河流进行监测，监测内容和结果（数据来源于《太湖水资源水环境研究》第133页）如下
表3 无锡市主要河流水质监测结果 （单位：mg/L）
| 河流 | 断面 | 挥发酚 | COD | TP | 石油类 | NH3-N | Cr6+ |
| 转水河 | W1 | 0.007 | 31.6 | 0.62 | 0.21 | 11.4 | 0.002 |
| W2 | 0.21 | 42.0 | 0.28 | 0.49 | 4.41 | 0.002 | |
| 新兴塘 | W3 | 0.784 | 38.3 | 0.27 | 0.57 | 5.60 | 0 |
| 九里河 | W4 | 0.432 | 27.4 | 0.21 | 0.41 | 4.50 | 0 |
| W5 | 0.004 | 12.5 | 5.51 | 0.09 | 3.83 | 0 | |
| 北兴塘 | W6 | 0.014 | 34.6 | 0.30 | 0.16 | 8.23 | 0.002 |
| W7 | 0.004 | 30.2 | 0.58 | 0.39 | 3.98 | 0.002 | |
| 东亭港 | W8 | 0.025 | 40.2 | 0.62 | 0. | 10.1 | 0.002 |
| 吼山港 | W9 | 0.004 | 26.1 | 0.24 | 0.15 | 6.10 | 0.025 |
| 三坝港 | W10 | 0.004 | 23.9 | 0.22 | 0.15 | 5.45 | 0.002 |
表4 据《地表水环境质量标准》（GB3838-2002） （单位：mg/L）
| 指标 | 一类 | 二类 | 三类 | 四类 | 五类 |
| 挥发酚 | 0.002 | 0.002 | 0.005 | 0.01 | 0.1 |
| COD | 15 | 15 | 20 | 30 | 40 |
| TP | 0.02 | 0.1 | 0.2 | 0.3 | 0.4 |
| 石油类 | 0.05 | 0.05 | 0.05 | 0.5 | 1.0 |
| NH3-N | 0.15 | 0.5 | 1.0 | 1.5 | 2.0 |
| Cr6+ | 0.01 | 0.05 | 0.05 | 0.05 | 0.1 |
依《地表水环境质量标准》建立决策集
3.4建立单因素评价矩阵
将断面W1的各实测值代入公式（1）得该断面的单因素评价矩阵R1
同理，可得其他各断面的单因素评价矩阵。
3.5确定各因素的权重
3.5.1通过对该体系中的指标两两比较得判断矩阵为：
3.5.2借助MATLAB计算判断矩阵的最大特征值及相应特征向量。
表5 计算所得各特征值如下表
| 6.1849 | -0.0365+1.0309i | -0.0365-1.0309i | -0.0366+0.2588i | -0.0366-0.2588i | -0.0386 |
对上面的判断矩阵归一化，可得
即为各指标的权重。
3.5.3进行一致性检验
由（4）式知
所以判断矩阵具有满意一致性，可以用做综合评价模型各指标的权重。
3.5.4监测断面的综合评价结果为：
据此判断断面W1的水质为五类，即重污染。同理可得各条河流各个断面的水质评价结果如下
对以上计算结果进行统计可得下表
表6 各条河流各个断面的水质评价结果
| 河流 | 转水河 | 新兴塘 | 九里河 | ||
| 断面 | W1 | W2 | W3 | W4 | W5 |
| b* | b5 | b5 | b5 | b5 | b5 |
| 检测结果 | 重污染 | 重污染 | 重污染 | 重污染 | 重污染 |
| 河流 | 北兴塘 | 东亭港 | 吼山港 | 三坝港 | |
| 断面 | W6 | W7 | W8 | W9 | W10 |
| b* | b4 | b5 | b5 | b3 | b3 |
| 检测结果 | 中污染 | 重污染 | 重污染 | 轻污染 | 轻污染 |
（1）无锡市周围数条河流的水质状况不容乐观：转水河、新兴塘、九里河、东亭港都已遭受了严重的污染，其他几条河流如北兴塘、吼山港、三坝港也在不同程度上存在着污染问题。
（2）水质评价牵涉因素众多，问题复杂，本文将模糊数学和层次分析法相结合建立数学模型，解决了由于不同污染因子对水质状况的影响程度不同而造成的评价难题。
5参考文献
黄贤金，王腊春.太湖水资源水环境研究[M].科学出版社，2008.
仇雁翎，陈玲，赵建夫.饮用水水质监测与分析[M].化学工业出版社，2006.
薛定宇，陈阳泉.高等应用数学问题的MATLAB求解[M].清华大学出版社，2004.
韩中庚.数学建模方法及其应用[M].高等教育出版社，2005.
焦瑞峰，吴泽宁，刘进国.水环境质量模糊层次综合评价模型及应用[J].黄河水资源保护科技成果与论文选编[M].黄河水利出版社，2005.
李士勇.工程模糊数学及其应用[M].哈尔滨工业大学出版社，2004.下载本文
