1. 引言
水质评价是水环境质量综合评价的简称,是一种在环境质量标准下对水质样本进行分类识别的问题 [1] 。随着水质评价研究的深入,水质评价工作在积极吸纳各相关学科和领域的研究成果的同时,自身也在不断发展和成熟。
将秦皇岛市桃林口水库作为研究区域,以2008~2017年十年水源站和出库站两个监测站监测数据的PH、总磷、氨氮、高锰酸盐指数、溶解氧和五日生化需氧量实测值数据作为水质评价的水质指标对桃林口水库水质进行评价,为秦皇岛市桃林口水库水质预测预警研究及工作提供研究基础与依据。
桃林口水库位于秦皇岛市西北部,滦河支流青龙河上,于1998年建成,控制流域面积5060 km2,总库容8.59亿m3,每年可为秦皇岛市提供1.82亿m3,为唐山、秦皇岛地区补充农业水源5.
2 亿立方米
。桃林口水库是一座具有城市供水、农业灌溉、防洪以及发电等综合效益的大型水利枢纽工程 [2] ,其水质对青龙县乃至秦皇岛市的社会稳定和经济发展起着举足轻重的作用。选取桃林口水库水源站和出库站两个监测断面(水源站和出库站)的水质指标数据对桃林口水库水质进行评价。
2. 研究方法
2.1. 单因子评价法
单因子评价是指将监测的实测值与各自对应的评价标准进行比较,其中最不好的一个被用作最后的水质类别 [3] 。单因子评价法是我国水质评价中应用最为广泛的水质评价方法 [4] ,其方法简单容易操作,可以直观的对水质监测指标予以描述,但该法是对于最坏情形的考量,对水质环境不能较为客观的反馈,导致最终对评价水体不能较为适度地保护。郭蓉 [5] 提到单因子评价法和综合评价法中的指数评价法是使用较多的水质评价方法,具有简单明了、直接反映水质状况与评价标准间的关系的特点。郑琨 [6] 等通过单因子指数法和综合污染指数法,对同水质类别进行洁净程度判别。刘荣荣等 [7] 采用单因子标准指数法对其进行水环境质量评价,研究结果表明浏河水质未达到地表水环境质量的要求,已受到一定程度的污染。
2.2. 综合污染指数法
综合污染指数法判断水质量情况的出发点是通过根据水质组分浓度相较于其环境质量标准的大小对比进而做出判断的 [8] 。计算公式为:
(1)
式中,Pj为j监测站的综合污染指数;n为参加评价的水质指标;Pi为i项水质指标的污染分指数;Ci为第i项水质指标的浓度值,Si为第i项水质指标的水质标准值,参照《地表水环境质量标准》(GB3838-2002)表1中的II类水质标准 [8] 。
郭晶 [8] 等采用综合污染指数法对洞庭湖水质进行综合评价,研究结果表示综合污染指数法虽不能对各断面的水质等级进行判定,但能对水质进行定量分析,可以清楚地反映洞庭湖的污染程度。
2.3. 神经网络评价法
神经网络法是20世纪80年代迅速崛起的新兴的非线性科学,是至今理论发展得较为完善的智能学习算法,因神经网络在水质模型应用中具有高速计算能、记忆能力、容错能力及学习能力等优点,使得其在水质评价中被广泛应用 [9] [10] ,其中较为常用的方法为反相传播(BP)网络模型和Hopfield网络模型等 [11] [12] [13] 。
曹阳阳 [14] 等应用RBF神经网络构建水质评价模型,对燕山南麓水库库区水质予以量化评价。结果表明库区水质属于II、III类水,存在一定程度的化学污染。何钰 [15] 等采用BP神经网络模型和综合指数评价法对嘉陵江南充段干流的水质进行评价,研究结果表示:与综合指数评价法相比,BP神经网络具有评价客观和应用方便等优点。
3. 水质指标特征
图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类水质标准限值),五日生化需氧量的浓度值总体表现为丰水期 > 平水期 > 枯水期,且五日生化需氧量的浓度值呈现出上升的趋势。
通过桃林口水库监测断面水源站和出库站十年监测数据结果可以看出,在一般情况下,在水库丰水期,总磷、氨氮、高锰酸盐指数和五日生化需氧量的浓度值较高,而在平水期和枯水期的浓度值较低。溶解氧浓度值随丰水期、平水期和枯水期的变化情况与总磷、氨氮、高锰酸盐指数和五日生化需氧量浓度值等的变化情况相反,溶解氧浓度值在枯水期较高,平水期和丰水期较低。可能是因为桃林口水库丰水期的降雨等将水库上游及水库周边富含氮、磷的土壤汇入水库,桃林口水库中有机物的含量得以升高,为微生物繁殖、生长提供创造条
![](//html.hanspub.org/file/12-2410773x11_hanspub.png)
Figure 1. Monitoring values of water quality indicators at water source and outbound stations
图1. 水源站和出库站各水质指标监测值
件。还有可能在丰水期暴雨形成径流,携带大量泥沙入库,造成淤积,影响水库库区的水质。并且桃林口水库丰水期集中在5~9月,这段时期水库的水温较于其他时间较高,有利于微生物的繁殖生长,微生物的繁殖生长过程会消耗水中的溶解氧,致使溶解氧含量在丰水期较低,总磷、氨氮、高锰酸盐指数和五日生化需氧量等浓度较高,因此桃林口水库丰水期的水质较于其他时期较差。
氨氮、高锰酸盐指数、五日生化需氧量等浓度值呈现一定的上升趋势,造成这个状况的原因可能是:青龙县的经济、社会发展造成污染源的增多、水底淤积的累积、丰水期降雨量及温度升高为微生物的繁殖增长提供机会等。
4. 评价结果分析
4.1. 单因子评价法结果
由表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类水质。可见,平水期和枯水期的水质相对较好,丰水期的水质相对较差。
4.2. 综合污染指数法结果
采用综合污染指数法对桃林口水库水源站和监测站两个监测断面2008~2017年十年全年、丰水期、平水期和枯水期等时期的综合污染指数分别进行计算,计算结果如表3所示,水源站全年的综合污染指数在0.437~0.578,出库站全年的综合污染指数在0.439~0.569,水质级别由表4可知,均为轻污染,其中2012年水源站综合污染指数为0.578,出库站综合污染指数为0.563,污染最为严重,2008年、2009年水源站和出库站综合污染指数均小于0.5,污染最轻。
![](Images/Table_Tmp.jpg)
Table 1. Water quality category of single source evaluation method for water source and outbound stations
表1. 单因子评价法水源站和水库站水质类别
![](Images/Table_Tmp.jpg)
Table 2. Single factor method do not count total phosphorus in the water source and outbound stations
表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,均为轻污染。一般情况下,水源站和出库站均是丰水期综合污染指数 > 平水期综合污染指数 > 枯水期综合污染指数,可见枯水期水质整体上优于丰水期和平水期,丰水期水质较差。
![](Images/Table_Tmp.jpg)
Table 3. Comprehensive pollution index of water source and outbound stations of Taolinkou Reservoir
表3. 桃林口水库水源站和出库站综合污染指数
![](Images/Table_Tmp.jpg)
Table 4. Classification criteria for water environmental quality pollution index
表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%,表明这三种水质指标对桃林口水库水质的影响较大,应加以关注。
![](//html.hanspub.org/file/12-2410773x12_hanspub.png)
Figure 2. Pollution share rate of water source (left) and outbound (right) stations in 2008~2017
图2. 水源站(左)和出库站(右)2008~2017年污染分担率
4.3. BP神经网络模型评价结果
4.3.1
. 训练样本的选择
本研究训练样本数据集选取《地表水环境质量标准》(GB3838-2002)中的总磷、氨氮、高锰酸盐指数、溶解氧和五日生化需氧量等的水质指标浓度限值作为训练集。根据水质标准将水质分为5类,I、II、III、IV和V类,形成五组样本区间,本研究运用rand(·)函数对五组样本范围随机插值,每类水质生成50个样本,共生成250个样本,样本生成规则如表5所示。
![](Images/Table_Tmp.jpg)
Table 5. Sample generation rule table
表5. 样本生成规则表
将总磷、氨氮、高锰酸盐指数、溶解氧和五日生化需氧量的监测数据作为输入变量,水质类别作为输出变量,对水质进行评价。取1、2、3、4和5分别作为I、II、III、IV和V类等5类水质类别的期望输出。
4.3.2
. 网络层数和隐含层的确定
根据Kolmogorov定理“对于具有一个隐层的3层BP神经网络,有足够多的隐节点数,就可以实现对任意非线性函数的逼近” [16] ,并且三层BP神经网络结构可以提高网络学习速度,而过于复杂的网络结构意味着过多的参数和模型精度的降低,所以本研究选择隐含层数为1的3层BP神经网络。
本研究中输入层为五项指标值(总磷、氨氮、高锰酸盐指数、溶解氧和五日生化需氧量值),输入层神经元数为5,输出层为水质类别,输出层神经元数为1。根据Komogorov定理和Hecht-Nielsen理论确定隐含层节点数范围为3~12,通过试错法,选择训练次数最少并且误差最小对应的神经元数,最终确定的隐含层节点数为7。故本研究的BP神经网络的拓扑结构最终确定为5:7:1,结构如图3所示。
![](//html.hanspub.org/file/12-2410773x13_hanspub.png)
Figure 3. BP neural network water quality prediction model structure diagram
图3. BP神经网络水质预测模型结构图
4.3.3
. 选择算法
程序的实现基于MATLAB环境,对生成的样本进行学习训练。建立的BP神经网络模型的参数设定如下,训练函数采用的是自适应学习率动量因子梯度训练法,输入层与隐含层之间的传递函数为对数S型函数logsig,隐含层与输出层之间的传递函数是正切S型函数tansig。期待误差为0.0001,学习速率设定0.0001,学习目标为0.00004,最大迭代次数为5000。
![](Images/Table_Tmp.jpg)
Table 6. Water quality evaluation results of water source and outbound stations based on BP neural network model
表6. 基于BP神经网络模型的水源站和出库站水质评价结果
BP神经网络模型对桃林口水库水源站和出库站水质评价结果水质类别基本一致。如表6所示,枯水期水质优于平水期和丰水期,水源站十年中有五年的枯水期水质为I类水质,出库站十年中有七年的枯水期水质为I类水质;水源站十年中有两年的枯水期水质为I类水质,出库站十年中有三年的枯水期水质为I类水质;2012年水源站和出库站的丰水期水质类别都为III类水质。2008~2017年十年全年平均水质都为II类水质,丰水期除2012年外,都为II类水质,满足集中式饮用水水源地水质指标要求。BP神经网络模型对桃林口水库水源站和出库站水质评价结果表明,平水期和枯水期水期时,出库站的水质优于水源站的水质。
4.4. 三种方法的比较
单因子评价法作为我国最常用的水质评价方法,简单易行,在桃林口水库这种污染浓度不大,污染状况简单的环境下,单因子评价法可以对桃林口水库水质做出定性的评价,判定出水源站和出库站2008~2017年十年的丰、平、枯水期的水质等级,结果表明,桃林口水库水质除总磷外,2008~2017年各项指标全年平均、丰、平、枯水期均达到II类水质目标要求,全因评价达到III类水质目标要求,满足集中式饮用水水源地水质指标要求。环境管理部门是可以将单因子评价方法运用于环境影响评价过程的,单因子评价法有利于水环境的保护,但是会在一定程度上降低水域的功能,尤其对于污染情况复杂的水域,单因子评价法将无法发挥出水环境的综合作用。
综合污染指数法计算方法相对来说较为复杂,可以对水质进行定量的分析,对比出水源站和出库站两个监测站点的水质状况,并能分析出污染分担率,但无法判定水库水质的等级。通过综合污染指数法可以看到,总磷、高锰酸盐有指数、五日生化需氧量这三项水质指标对桃林口水库的水质影响较大。
BP神经网络模型只通过水质标准进行学习,摒弃了人为干预,获得较为复杂的非线性处理的能力,BP神经网络模型是基于多个水质指标的综合评价,在一定程度上对单因子评价法的不足进行弥补,更具有客观性。丰水期时水源站和出库站在两种方法评价下的评价结果等级差别不大,只有2015、2016年时单因子评价法评价结果为III类水质,BP神经网络模型评价结果为Ⅱ类水质;平水期时水源站和出库站在两种方法评价下的评价结果等级差别也不大;但枯水期时水源站和出库站在两种方法评价下的评价结果等级差别较大,水源站评价结果完全一致的有5次,出库站评价结果完全一致的只有3次,BP神经网络模型评价下枯水期水库水质多为Ⅰ类水质;水源站和出库站的全年平均在两种方法评价下的评价结果等级几乎一致。可以看出BP人工神经网络模型评价得出水质类别结果是明显优于单因子评价法得出的水质类别结果,这是由于在BP人工神经网络模型评价水质类别的过程中,采用的是一种更为客观的方式,基于多个水质指标的综合评价,因此得到了一个与单因子评价法相较下较为客观的评价结果。
5. 结论
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)。
参考文献