黄 明
(安徽建筑工业学院 环境与能源工程学院,安徽 合肥230022)
巢湖属浅水、人工控制型湖泊,在正常蓄水位8.0m时,其水域面积约755.0km2,入湖河流共33条,分别属杭埠-丰乐河、派河、南淝河-店埠河、柘皋河、白石山河、裕溪河7条水系[1-3]。近年来,巢湖水质富营养化问题仍然很突出,因此,对巢湖水质富营养化进行综合评价是一项重要的工作。本文依据2010年巢湖水质监测数据建立评价模型,并利用Matlab软件尝试对巢湖水质富营养化进行综合评价。
巢湖12个水质监测点位置见图1。选取南淝河入湖区、西半湖湖心、忠庙、东半湖湖心和巢湖船厂5个监测点数据对巢湖水质富营养化进行综合评价,见表1。
图1 巢湖水质监测点位置示意图
巢湖水质富营养化评价标准见表2。
(1)建立无量纲化实测数据矩阵和评价标准矩阵[4-5]
根据表1和表2,可得到实测数据矩阵X=(χij)5×4和评价标准矩阵Y=(ykt)4×5,
表1 2010年巢湖水质实测数据
表2 巢湖水质富营养化评价标准[1]
再建立无量纲化实测数据矩阵:A=(aij)(i=1,2,3,4,5;j=1,2,3,4)和无量纲化评价标准矩阵B=(bkt)(k=1,2,3,4;t=1,2,3,4,5),其中
利用Matlab可以得到
(2)计算各评价指标的权重
首先计算矩阵B的各向量的均值与标准差:
再计算变异系数:
得到各指标权重向量为
根据权重大小,总磷所起的作用最大,耗氧量次之,透明度与总氮权重大小接近,作用相对较小。
(3)建立水质综合评价模型
利用向量之间的距离来衡量2个向量之间的接近程度,可以调用Matlab软件中的欧氏距离函数dist(w,p)和绝对值距离函数 mandist(w,p)来计算向量间的距离[4-5]。
计算A中各行向量到B中各列向量的欧氏距离:
计算A中各行向量到B中各列向量的绝对值距离:
计算结果如表3和表4。
从表3和和表4可以看出,尽管欧氏距离与绝对值距离意义不同,但对各监测点的富营养化的评价等级是一样的。
表3 欧氏距离判别表
表4 绝对值距离判别表
对湖泊富营养化程度评价的方法较多,有的是将各指标分别与营养状态对照进行判定,但各指标间往往缺乏有机联系,尤其是当各指标反映的营养状态不一致时,综合评价就存在不少困难。
通过运用Matlab软件中的欧氏距离函数dist(w,p)和绝对值距离函数 mandist(w,p)来计算向量间的距离,并建立水质综合评价模型对各监测点实测数据进行评价,其评价结果一致,表明该评价方法具有稳定性。
[1]屠清瑛,顾丁锡,徐卓然,等.全国湖泊富营养化研究之二——巢湖富营养化研究[M].合肥:中国科技大学出版社,1991.
[2]安徽水利厅,安徽省环境保护局.安徽省水功能区划[M].北京:中国水利水电出版社,2004.
[3]张之源,王培华,张崇岱.巢湖营养化状况评价及水质恢复探讨[J].环境科学研究,1999,12(5):45-48.
[4]胡守信,李柏年.基于 MATLAB的数学实验[M].北京:科学出版社,2004.
[5]苏金明,阮沈勇.MATLAB6.1实用指南(上册)[M].北京:电子工业出版社,2002.