利用皇岛桃林口水库2008~2017年十年的实测数据,选取PH、总磷(TP)、氨氮(NH3-N)、高锰酸盐指数(CODMn)、溶解氧(DO)和五日生化需氧量(BOD5) 6项指标,采用单因子评价法、综合污染指数法和BP神经网络模型等3种方法进行水质评价。研究结果表明:1) 单因子水质评价结果表明,桃林口水库水质全因评价达到III类水质目标要求。2) 综合污染指数法表明,水源站全年的综合污染指数在0.437~0.578,出库站全年的综合污染指数在0.439~0.569,均为轻污染。3) BP神经网络模型表明,桃林口水库十年全年平均水质都为II类水质。该水质评价研究经验为桃林口水库的管理提供较为客观的决策依据,也可为其他地区的水质评价工作提供参考借鉴。 Based on measured data during 2008~2017 at Taolinkou Reservoir in Qinhuangdao, six water quality indicators (PH, total phosphorus (TP), ammonia nitrogen (NH3-N), permanganate index (CODMn), dis-solved oxygen (DO) and five-day BOD5) were evaluated by three different evaluation methods (single factor evaluation method, comprehensive pollution index method and BP neural network model). The results showed that: 1) Single-factor water quality evaluation indicated that the water quality all-caliber evaluation met the requirements of Class III water quality. 2) The comprehensive pollution index method showed that the comprehensive pollution index of the water source station was 0.437~0.578, and the comprehensive pollution index of the outbound station was 0.439~0.569, which was light pollution. 3) The BP neural network model showed that the average annual water quality was Grade II water quality. It provided a more objective decision-making basis for the management of Taolinkou Reservoir, and could also provide reference for water quality assessment in other areas.
戴天骄1,赵林1,陈亮1,刘琦1,冯琛雅1,赵明杰2
1天津大学环境科学与工程学院,天津
2秦皇岛天大环保研究院有限公司,河北 秦皇岛
收稿日期:2019年6月27日;录用日期:2019年7月17日;发布日期:2019年7月29日
利用皇岛桃林口水库2008~2017年十年的实测数据,选取PH、总磷(TP)、氨氮(NH3-N)、高锰酸盐指数(CODMn)、溶解氧(DO)和五日生化需氧量(BOD5) 6项指标,采用单因子评价法、综合污染指数法和BP神经网络模型等3种方法进行水质评价。研究结果表明:1) 单因子水质评价结果表明,桃林口水库水质全因评价达到III类水质目标要求。2) 综合污染指数法表明,水源站全年的综合污染指数在0.437~0.578,出库站全年的综合污染指数在0.439~0.569,均为轻污染。3) BP神经网络模型表明,桃林口水库十年全年平均水质都为II类水质。该水质评价研究经验为桃林口水库的管理提供较为客观的决策依据,也可为其他地区的水质评价工作提供参考借鉴。
关键词 :水质,评价,BP神经网络,桃林口水库
Copyright © 2019 by author(s) and Wuhan University.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
水质评价是水环境质量综合评价的简称,是一种在环境质量标准下对水质样本进行分类识别的问题 [
将秦皇岛市桃林口水库作为研究区域,以2008~2017年十年水源站和出库站两个监测站监测数据的PH、总磷、氨氮、高锰酸盐指数、溶解氧和五日生化需氧量实测值数据作为水质评价的水质指标对桃林口水库水质进行评价,为秦皇岛市桃林口水库水质预测预警研究及工作提供研究基础与依据。
桃林口水库位于秦皇岛市西北部,滦河支流青龙河上,于1998年建成,控制流域面积5060 km2,总库容8.59亿m3,每年可为秦皇岛市提供1.82亿m3,为唐山、秦皇岛地区补充农业水源5. 2 亿立方米 。桃林口水库是一座具有城市供水、农业灌溉、防洪以及发电等综合效益的大型水利枢纽工程 [
单因子评价是指将监测的实测值与各自对应的评价标准进行比较,其中最不好的一个被用作最后的水质类别 [
综合污染指数法判断水质量情况的出发点是通过根据水质组分浓度相较于其环境质量标准的大小对比进而做出判断的 [
P j = 1 / n ∑ i = 1 n P i , P i = C i / S i (1)
式中,Pj为j监测站的综合污染指数;n为参加评价的水质指标;Pi为i项水质指标的污染分指数;Ci为第i项水质指标的浓度值,Si为第i项水质指标的水质标准值,参照《地表水环境质量标准》(GB3838-2002)表1中的II类水质标准 [
郭晶 [
神经网络法是20世纪80年代迅速崛起的新兴的非线性科学,是至今理论发展得较为完善的智能学习算法,因神经网络在水质模型应用中具有高速计算能、记忆能力、容错能力及学习能力等优点,使得其在水质评价中被广泛应用 [
曹阳阳 [
图1显示桃林口水库监测断面水源站和出库站十年监测数据结果,出库站PH、总磷、溶解氧和五日生化需氧量浓度均高于水源站各项指标浓度,水源站氨氮、高锰酸盐指数浓度高于出库站氨氮、高锰酸盐指数浓度。PH值范围为7.02~8.41,呈弱碱性。水源站和出库站监测断面总磷浓度十年全年平均值范围分别在0.0103~0.0276 mg·L−1和0.013~0.027 mg·L−1之间,大部分低于0.025 mg·L−1 (II类水质标准限值),均低于0.05 mg·L−1 (III类水质标准限值),总磷浓度值总体表现为丰水期 > 平水期 > 枯水期。水源站和出库站监测断面氨氮浓度十年全年平均值范围分别在0.061~0.4 mg·L−1和0.054~0.285 mg·L−1之间,均低于0.5 mg·L−1 (II类水质标准限值),氨氮浓度总体表现为丰水期 > 平水期 > 枯水期,且随着时间的发展,氨氮的浓度值呈现出上升的趋势。水源站和出库站监测断面高锰酸盐指数浓度十年全年平均值范围分别在1.897~2.805 mg·L−1和1.95~2.54 mg·L−1之间,均低于4 mg·L−1 (II类水质标准限值),高锰酸盐指数浓度总体表现为丰水期 > 平水期 > 枯水期,高锰酸盐指数浓度略微呈现上升的趋势。水源站和出库站监测断面溶解氧浓度十年全年平均值范围分别在8.821~11.08 mg·L−1和9.93~11.468 mg·L−1之间,均高于7.5 mg·L−1 (I类水质标准限值),溶解氧浓度总体表现为枯水期 > 平水期 > 丰水期。水源站和出库站监测断面五日生化需氧量浓度十年全年平均值范围分别在2.107~2.656 mg·L−1和2.057~2.704 mg·L−1之间,均低于3 mg·L−1 (I类水质标准限值),五日生化需氧量的浓度值总体表现为丰水期 > 平水期 > 枯水期,且五日生化需氧量的浓度值呈现出上升的趋势。
通过桃林口水库监测断面水源站和出库站十年监测数据结果可以看出,在一般情况下,在水库丰水期,总磷、氨氮、高锰酸盐指数和五日生化需氧量的浓度值较高,而在平水期和枯水期的浓度值较低。溶解氧浓度值随丰水期、平水期和枯水期的变化情况与总磷、氨氮、高锰酸盐指数和五日生化需氧量浓度值等的变化情况相反,溶解氧浓度值在枯水期较高,平水期和丰水期较低。可能是因为桃林口水库丰水期的降雨等将水库上游及水库周边富含氮、磷的土壤汇入水库,桃林口水库中有机物的含量得以升高,为微生物繁殖、生长提供创造条
图1. 水源站和出库站各水质指标监测值
件。还有可能在丰水期暴雨形成径流,携带大量泥沙入库,造成淤积,影响水库库区的水质。并且桃林口水库丰水期集中在5~9月,这段时期水库的水温较于其他时间较高,有利于微生物的繁殖生长,微生物的繁殖生长过程会消耗水中的溶解氧,致使溶解氧含量在丰水期较低,总磷、氨氮、高锰酸盐指数和五日生化需氧量等浓度较高,因此桃林口水库丰水期的水质较于其他时期较差。
氨氮、高锰酸盐指数、五日生化需氧量等浓度值呈现一定的上升趋势,造成这个状况的原因可能是:青龙县的经济、社会发展造成污染源的增多、水底淤积的累积、丰水期降雨量及温度升高为微生物的繁殖增长提供机会等。
由表1和表2可知,在2008~2017年十年的水质监测时间段内,桃林口水库水源站和出库站在单因子评价法下水质类别评价结果基本一致,除2016年水源站平水期的评价结果是I类水,出库站平水期评价结果是II类水外,其余水质类别评价结果完全一致。2012年全年平均和丰水期的水质评价结果是III类水,2015、2016年丰水期水质评价结果是III类水。桃林口水库水质除总磷外,2008~2017年各项指标全年平均、丰、平、枯水期均达到II类水质目标要求,全因评价达到III类水质目标要求,满足集中式饮用水水源地水质指标要求。从丰水期、平水期和枯水期的水质评价就结果来看,总磷不参与水质评价时,水源站和出库站2008~2017年十年的水质100%都达到II类水质标准。总磷不参与评价时,水源站和出库站均是10%丰水期、20%平水期和20%枯水期为I类水质。全因评价时,水源站丰水期30%为III类水质,10%平水期和10%枯水期为I类水质;出库站丰水期30%为III类水质,100%平水期为II类水质,10%枯水期为I类水质。可见,平水期和枯水期的水质相对较好,丰水期的水质相对较差。
采用综合污染指数法对桃林口水库水源站和监测站两个监测断面2008~2017年十年全年、丰水期、平水期和枯水期等时期的综合污染指数分别进行计算,计算结果如表3所示,水源站全年的综合污染指数在0.437~0.578,出库站全年的综合污染指数在0.439~0.569,水质级别由表4可知,均为轻污染,其中2012年水源站综合污染指数为0.578,出库站综合污染指数为0.563,污染最为严重,2008年、2009年水源站和出库站综合污染指数均小于0.5,污染最轻。
时期 | 水源站水质类别 | |||||||||
---|---|---|---|---|---|---|---|---|---|---|
2008年 | 2009年 | 2010年 | 2011年 | 2012年 | 2013年 | 2014年 | 2015年 | 2016年 | 2017年 | |
丰水期 | II类 | II类 | II类 | II类 | III类 | II类 | II类 | III类 | II类 | II类 |
平水期 | II类 | II类 | II类 | II类 | II类 | II类 | II类 | II类 | I类 | II类 |
枯水期 | II类 | II类 | II类 | II类 | II类 | II类 | II类 | II类 | II类 | I类 |
全年平均 | II类 | II类 | II类 | II类 | III类 | II类 | II类 | II类 | II类 | II类 |
时期 | 出库站水质类别 | |||||||||
丰水期 | II类 | II类 | II类 | II类 | III类 | II类 | II类 | III类 | III类 | II类 |
平水期 | II类 | II类 | II类 | II类 | II类 | II类 | II类 | II类 | II类 | II类 |
枯水期 | II类 | II类 | II类 | II类 | II类 | II类 | II类 | II类 | II类 | I类 |
全年平均 | II类 | II类 | II类 | II类 | III类 | II类 | II类 | II类 | II类 | II类 |
表1. 单因子评价法水源站和水库站水质类别
时期 | 水源站水质类别 | |||||||||
---|---|---|---|---|---|---|---|---|---|---|
2008年 | 2009年 | 2010年 | 2011年 | 2012年 | 2013年 | 2014年 | 2015年 | 2016年 | 2017年 | |
丰水期 | I类 | II类 | II类 | II类 | II类 | II类 | II类 | II类 | II类 | II类 |
平水期 | I类 | II类 | II类 | II类 | II类 | II类 | II类 | II类 | I类 | II类 |
枯水期 | II类 | II类 | II类 | II类 | II类 | II类 | II类 | II类 | I类 | I类 |
全年平均 | I类 | II类 | II类 | II类 | II类 | II类 | II类 | II类 | II类 | II类 |
时期 | 出库站水质类别 | |||||||||
丰水期 | II类 | I类 | II类 | II类 | II类 | II类 | II类 | II类 | II类 | II类 |
平水期 | I类 | I类 | II类 | II类 | II类 | II类 | II类 | II类 | II类 | II类 |
枯水期 | II类 | II类 | II类 | II类 | I类 | II类 | II类 | II类 | II类 | I类 |
全年平均 | II类 | II类 | II类 | II类 | II类 | II类 | II类 | II类 | II类 | II类 |
表2. 单因子总磷不计评价水源和出库站水质类别
水源站丰水期综合污染指数在0.435~0.688,均为轻污染;出库站丰水期综合污染指数在0.506~0.918,其中2012、2015、2016年综合污染指数大于0.7,为中污染。水源站平水期综合污染指数在0.457~0.564,出库站平水期综合污染指数在0.407~0.507,均为轻污染。水源站枯水期综合污染指数在0.391~0.555,出库站枯水期综合污染指数在0.395~0.579,均为轻污染。一般情况下,水源站和出库站均是丰水期综合污染指数 > 平水期综合污染指数 > 枯水期综合污染指数,可见枯水期水质整体上优于丰水期和平水期,丰水期水质较差。
水源站 | 出库站 | |||||||
---|---|---|---|---|---|---|---|---|
全年 | 丰水期 | 平水期 | 枯水期 | 全年 | 丰水期 | 平水期 | 枯水期 | |
2008年 | 0.437 | 0.435 | 0.462 | 0.429 | 0.461 | 0.506 | 0.484 | 0.452 |
2009年 | 0.459 | 0.495 | 0.487 | 0.483 | 0.439 | 0.575 | 0.441 | 0.396 |
2010年 | 0.533 | 0.486 | 0.529 | 0.478 | 0.529 | 0.535 | 0.407 | 0.453 |
2011年 | 0.566 | 0.607 | 0.492 | 0.490 | 0.484 | 0.594 | 0.465 | 0.395 |
2012年 | 0.578 | 0.634 | 0.564 | 0.466 | 0.563 | 0.784 | 0.486 | 0.398 |
2013年 | 0.527 | 0.688 | 0.493 | 0.465 | 0.514 | 0.530 | 0.507 | 0.579 |
2014年 | 0.551 | 0.546 | 0.551 | 0.555 | 0.519 | 0.560 | 0.489 | 0.449 |
2015年 | 0.536 | 0.589 | 0.457 | 0.538 | 0.569 | 0.875 | 0.438 | 0.479 |
2016年 | 0.553 | 0.682 | 0.488 | 0.450 | 0.547 | 0.918 | 0.496 | 0.401 |
2017年 | 0.517 | 0.549 | 0.428 | 0.391 | 0.500 | 0.682 | 0.498 | 0.439 |
表3. 桃林口水库水源站和出库站综合污染指数
污染指数 | ≤0.20 | 0.21~0.40 | 0.41~0.70 | 0.71~1.00 | 1.01~1.99 | ≥2.00 |
---|---|---|---|---|---|---|
污染级别 | 清洁 | 尚清洁 | 轻污染 | 中污染 | 重污染 | 严重污染 |
表4. 水环境质量污染指数分级标准
水质指标的污染分担率如图2所示,总磷、高锰酸盐指数和五日生化需氧量污染分担率高于PH、氨氮、溶解氧。水源站总磷十年单项污染指数范围为0.455~1.105,污染分担率在14%~32%之间,2008年~2014年均大于25%,2015年后总磷的污染分担率有所下降;出库站总磷十年单项污染指数范围为0.535~1.083,污染分担率在20%左右浮动。水源站和出库站高锰酸盐指数十年单项污染指数范围分别为0.474~0.633和0.488~0.663,污染分担率均在20%左右。水源站和出库站五日生化需氧量十年单项污染指数范围分别为0.694~0.885和0.686~0.901,污染分担率均在25%左右,近两年五日生化需氧量污染分担率有升高的趋势。总磷、高锰酸盐有指数、五日生化需氧量的污染分担率之和基本达到70%~80%,表明这三种水质指标对桃林口水库水质的影响较大,应加以关注。
图2. 水源站(左)和出库站(右)2008~2017年污染分担率
本研究训练样本数据集选取《地表水环境质量标准》(GB3838-2002)中的总磷、氨氮、高锰酸盐指数、溶解氧和五日生化需氧量等的水质指标浓度限值作为训练集。根据水质标准将水质分为5类,I、II、III、IV和V类,形成五组样本区间,本研究运用rand(·)函数对五组样本范围随机插值,每类水质生成50个样本,共生成250个样本,样本生成规则如表5所示。
样本编号 | 样本区间范围 | 期望输出 | 识别结果 |
---|---|---|---|
1 | (0,0,0,18,0)~(0.01,0.15,2,7.5,3) | 1 | I |
2 | (0.01,0.15,2,7.5,3)~(0.025,0.5,4,6,3) | 2 | II |
3 | (0.025,0.5,4,6,3)~(0.05,1,6,5,4) | 3 | III |
4 | (0.05,1,6,5,4)~(0.1,1.5,10,3,6) | 4 | IV |
5 | (0.1,1.5,10,3,6)~(0.2,2,15,2,10) | 5 | V |
表5. 样本生成规则表
将总磷、氨氮、高锰酸盐指数、溶解氧和五日生化需氧量的监测数据作为输入变量,水质类别作为输出变量,对水质进行评价。取1、2、3、4和5分别作为I、II、III、IV和V类等5类水质类别的期望输出。
根据Kolmogorov定理“对于具有一个隐层的3层BP神经网络,有足够多的隐节点数,就可以实现对任意非线性函数的逼近” [
本研究中输入层为五项指标值(总磷、氨氮、高锰酸盐指数、溶解氧和五日生化需氧量值),输入层神经元数为5,输出层为水质类别,输出层神经元数为1。根据Komogorov定理和Hecht-Nielsen理论确定隐含层节点数范围为3~12,通过试错法,选择训练次数最少并且误差最小对应的神经元数,最终确定的隐含层节点数为7。故本研究的BP神经网络的拓扑结构最终确定为5:7:1,结构如图3所示。
图3. BP神经网络水质预测模型结构图
程序的实现基于MATLAB环境,对生成的样本进行学习训练。建立的BP神经网络模型的参数设定如下,训练函数采用的是自适应学习率动量因子梯度训练法,输入层与隐含层之间的传递函数为对数S型函数logsig,隐含层与输出层之间的传递函数是正切S型函数tansig。期待误差为0.0001,学习速率设定0.0001,学习目标为0.00004,最大迭代次数为5000。
时期 | 水源站水质类别 | |||||||||
---|---|---|---|---|---|---|---|---|---|---|
2008年 | 2009年 | 2010年 | 2011年 | 2012年 | 2013年 | 2014年 | 2015年 | 2016年 | 2017年 | |
丰水期 | II类 | II类 | II类 | II类 | III类 | II类 | II类 | II类 | II类 | II类 |
平水期 | I类 | II类 | II类 | II类 | II类 | II类 | II类 | I类 | I类 | II类 |
枯水期 | I类 | II类 | II类 | II类 | I类 | I类 | II类 | II类 | I类 | I类 |
全年平均 | II类 | II类 | II类 | II类 | II类 | II类 | II类 | II类 | II类 | II类 |
时期 | 出库站水质类别 | |||||||||
丰水期 | II类 | II类 | II类 | II类 | III类 | II类 | II类 | II类 | II类 | II类 |
平水期 | I类 | II类 | II类 | II类 | II类 | II类 | II类 | I类 | I类 | II类 |
枯水期 | I类 | I类 | II类 | I类 | I类 | II类 | I类 | II类 | I类 | I类 |
全年平均 | II类 | II类 | II类 | II类 | II类 | II类 | II类 | II类 | II类 | II类 |
表6. 基于BP神经网络模型的水源站和出库站水质评价结果
BP神经网络模型对桃林口水库水源站和出库站水质评价结果水质类别基本一致。如表6所示,枯水期水质优于平水期和丰水期,水源站十年中有五年的枯水期水质为I类水质,出库站十年中有七年的枯水期水质为I类水质;水源站十年中有两年的枯水期水质为I类水质,出库站十年中有三年的枯水期水质为I类水质;2012年水源站和出库站的丰水期水质类别都为III类水质。2008~2017年十年全年平均水质都为II类水质,丰水期除2012年外,都为II类水质,满足集中式饮用水水源地水质指标要求。BP神经网络模型对桃林口水库水源站和出库站水质评价结果表明,平水期和枯水期水期时,出库站的水质优于水源站的水质。
单因子评价法作为我国最常用的水质评价方法,简单易行,在桃林口水库这种污染浓度不大,污染状况简单的环境下,单因子评价法可以对桃林口水库水质做出定性的评价,判定出水源站和出库站2008~2017年十年的丰、平、枯水期的水质等级,结果表明,桃林口水库水质除总磷外,2008~2017年各项指标全年平均、丰、平、枯水期均达到II类水质目标要求,全因评价达到III类水质目标要求,满足集中式饮用水水源地水质指标要求。环境管理部门是可以将单因子评价方法运用于环境影响评价过程的,单因子评价法有利于水环境的保护,但是会在一定程度上降低水域的功能,尤其对于污染情况复杂的水域,单因子评价法将无法发挥出水环境的综合作用。
综合污染指数法计算方法相对来说较为复杂,可以对水质进行定量的分析,对比出水源站和出库站两个监测站点的水质状况,并能分析出污染分担率,但无法判定水库水质的等级。通过综合污染指数法可以看到,总磷、高锰酸盐有指数、五日生化需氧量这三项水质指标对桃林口水库的水质影响较大。
BP神经网络模型只通过水质标准进行学习,摒弃了人为干预,获得较为复杂的非线性处理的能力,BP神经网络模型是基于多个水质指标的综合评价,在一定程度上对单因子评价法的不足进行弥补,更具有客观性。丰水期时水源站和出库站在两种方法评价下的评价结果等级差别不大,只有2015、2016年时单因子评价法评价结果为III类水质,BP神经网络模型评价结果为Ⅱ类水质;平水期时水源站和出库站在两种方法评价下的评价结果等级差别也不大;但枯水期时水源站和出库站在两种方法评价下的评价结果等级差别较大,水源站评价结果完全一致的有5次,出库站评价结果完全一致的只有3次,BP神经网络模型评价下枯水期水库水质多为Ⅰ类水质;水源站和出库站的全年平均在两种方法评价下的评价结果等级几乎一致。可以看出BP人工神经网络模型评价得出水质类别结果是明显优于单因子评价法得出的水质类别结果,这是由于在BP人工神经网络模型评价水质类别的过程中,采用的是一种更为客观的方式,基于多个水质指标的综合评价,因此得到了一个与单因子评价法相较下较为客观的评价结果。
1) 桃林口水库水质呈弱碱性,水库丰水期,总磷、氨氮、高锰酸盐指数和五日生化需氧量的浓度值较高,而在平水期和枯水期的浓度值较低;溶解氧浓度值在枯水期较高,平水期和丰水期较低。
2) 单因子水质评价结果显示,桃林口水库水质除总磷外,2008~2017年各项指标全年平均、丰、平、枯水期均达到II类水质目标要求,全因评价达到III类水质目标要求,满足集中式饮用水水源地水质指标要求。
3) 水源站全年的综合污染指数在0.437~0.578,出库站全年的综合污染指数在0.439~0.569,均为轻污染。枯水期水质整体上优于丰水期和平水期,丰水期水质较差。且总磷、高锰酸盐有指数、五日生化需氧量这三项水质指标的污染分担率之和基本达到70%~80%,对桃林口水库的水质影响较大。
4) BP神经网络模型对桃林口水库水质评价结果显示,2008~2017年十年全年平均水质都为II类水质,丰水期除2012年外,都为II类水质,满足集中式饮用水水源地水质指标要求。
感谢我的导师赵林老师和陈亮老师还有我的同学们刘琦和冯琛雅的无私帮助,无论是生活中还是科研过程中,因为他们的陪伴让我的研究生生涯很是圆满且快乐,望大家以后的日子里仍能保持今日的温暖和善良。
水质自动监测系统在地表水水库生态安全预警中的应用研究(项目编码:17273905D)。
戴天骄,赵 林,陈 亮,刘 琦,冯琛雅,赵明杰. 桃林口水库水污染特征及水质评价Pollution Characterization and Water Quality Assessment of Taolinkou Reservoir[J]. 水资源研究, 2019, 08(04): 412-422. https://doi.org/10.12677/JWRR.2019.84048