为掌握秦皇岛桃林口水库未来水质变化的状况,本文选取总磷(TP)、硝酸盐氮(NO 3 -N)、氨氮(NH 3 -N)、高锰酸盐指数(COD Mn )和溶解氧(DO)等5项水质指标,基于2014~2015年水源站和出库站两个监测站的5项水质指标的监测数据,确定WASP模型的参数;再根据水源站2016、2017年5项水质指标来预测出库站的水质情况。研究结果表明:该WASP模型对5项水质指标的模拟值的判定系数R2均大于0.7,相关系数r均大于0.83;各项水质指标的实测值与模拟值的相对误差一般低于10%。WASP模拟效果良好,可以运用到桃林口水库的预测预警系统中和水库水质预测与管理的工作中,同时为秦皇岛市青龙县水环境保护与管理工作提供决策性参考与依据。 In order to grasp the future water quality changes of Taolinkou Reservoir in Qinhuangdao, the total phosphorus (TP), nitrate nitrogen (NO 3 -N), ammonia nitrogen (NH 3 -N), permanganate index (COD Mn ) and dissolved oxygen (DO) were selected. Based on these five water quality monitoring data in 2014-2015 at the water source and the outbound stations, the parameters of WASP model were determined and then used to predict water quality of Taolinkou Reservoir in 2016 and 2017. The results show that the R2 of the simulated values of the five water quality indicators of the WASP model is greater than 0.7, and the correlation coefficient is greater than 0.83; the relative error between the measured values and the simulated values of the water quality indicators is generally less than 10%. It is shown that the established WASP can simulate water quality well and apply to the prediction and early warning system of Taolinkou Reservoir. At the same time, it provides decision-making reference and basis for Qinhuangdao water environmental protection and management work.
戴天骄1,赵林1,陈亮1,冯琛雅1,刘琦1,赵明杰2
1天津大学环境科学与工程学院,天津
2秦皇岛天大环保研究院有限公司,河北 秦皇岛
收稿日期:2019年8月7日;录用日期:2019年8月29日;发布日期:2019年9月17日
为掌握秦皇岛桃林口水库未来水质变化的状况,本文选取总磷(TP)、硝酸盐氮(NO3-N)、氨氮(NH3-N)、高锰酸盐指数(CODMn)和溶解氧(DO)等5项水质指标,基于2014~2015年水源站和出库站两个监测站的5项水质指标的监测数据,确定WASP模型的参数;再根据水源站2016、2017年5项水质指标来预测出库站的水质情况。研究结果表明:该WASP模型对5项水质指标的模拟值的判定系数R2均大于0.7,相关系数r均大于0.83;各项水质指标的实测值与模拟值的相对误差一般低于10%。WASP模拟效果良好,可以运用到桃林口水库的预测预警系统中和水库水质预测与管理的工作中,同时为秦皇岛市青龙县水环境保护与管理工作提供决策性参考与依据。
关键词 :水质,模拟,WASP模型,桃林口水库
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/
水质预测是对研究的水体的质量做出未来的预测,是水污染综合防治的基础工作。自上世纪中期开始,随着经济的快速发展、城市化进程的加快,水体污染问题受到广泛的关注,对水环境质量变化趋势进行的预测分析愈加受到重视。常用的水质预测方法主要有水质模型、灰色系统理论、神经网络等 [
水质模型是较早引入水质预测的方法之一,是可以随时间和空间的变化,定量地描述环境污染物在水体中迁移转化规律以及其他因素之间相互影响的数学表达方式 [
水质模型现今处于广泛应用及深入研究的时期,在于完善QUAL、WASP等多介质的水质模型,在我国WASP水质模型针对不同水体都取得了进一步的研究与发展。张佩芳 [
本文主要针对秦皇岛市桃林口水库水源地保护的相关问题,通过对水库水质模型的构建与模拟,来对桃林口水库的水质模拟和保护产生实际意义,同时,对于我国北方的大型水库水质模型的研究和水源地保护的研究具有借鉴作用。
桃林口水库位于秦皇岛市西北部,滦河支流青龙河上,于1998年建成,控制流域面积5060平方公里,总库容8.59亿立方米,每年可为秦皇岛市提供1.82亿立方米,为唐山、秦皇岛地区补充农业水源5. 2 亿立方米 。桃林口水库是一座具有城市供水、农业灌溉、防洪以及发电等综合效益的大型水利枢纽工程 [
本研究选取桃林口水库水源站和出库站两个监测断面水质指标数据建立WASP模型,利用水源站的水质指标监测值模拟出库站水质指标值。水源站和出库站的地理位置图如图1所示。
图1. 桃林口水库各监测点地理位置图
样本集数据来源于2014~2015年桃林口水库水源站和出库站两个监测站的水质监测的监测资料,利用2016~2017年水源站的监测数据对出库站的总磷(TP)、硝酸盐氮(NO3-N)、氨氮(NH3-N)、高锰酸盐指数(CODMn)和溶解氧(DO)等5项水质指标进行模拟研究。
WASP模型本质上是利用各个物质的质量平衡方程,通过平移扩散的公式,来描述水体中各种物质的迁移和转化。WASP模型的质量平衡方程 [
∂ C ∂ t = − ∂ ( U x C ) ∂ t − ∂ ( U y C ) ∂ t − ∂ ( U z C ) ∂ t + ∂ ( E x ∂ C ∂ x ) ∂ t + ∂ ( E z ∂ C ∂ z ) ∂ t + S L + S B + S K (1)
式中:C为水质组分的浓度,mg/L;T为时间步长,d;Ux,Uy,Uz为纵向、横向和垂向速度,m/d;Ex,Ey,Ez为纵向、横向和垂向扩散系数,m2/d;SL项为点源和面源负荷,g/m3;SB项为边界负荷,g/m3;SK项为水质组分的总转化率,g/m3·d。
假设水体的垂向和横向是均匀的,将可以得到以下的一维水质组分的运动方程:
∂ ( A C ) ∂ t = ∂ ( − U x A C + E x A ∂ C ∂ x ) ∂ x + A S L + A S B + A S K (2)
式中:A为水体的截面面积,m2;其余方程符号意义同上。
WASP模型的建立需要输入水源站和出库站两个监测断面的基本信息,包括一维河流的物理模型,长度、宽度、水深、流量等,还有气象条件如水温、光照等,边界参数、水化学参数及相应的动力学参数。本研究基于2014~2015年水源站和出库站两个监测站的监测数据来建立WASP模型,通过试算法确定参数的调节和选定,其中的动力学参数的设定参考了用户手册和通过相关资料的参数取值确定,如表1所示,而其他相关数据从秦皇岛市青龙县环境监测站相关部门获得。
参数名称 | 取值 | 单位 | 范围 |
---|---|---|---|
20 ℃ 时COD衰减速率 | 0.06 | 0~5.6 | |
COD温度系数 | 1.05 | d−1 | 0~1.07 |
20 ℃ 时氨氮硝化速率 | 0.02 | 0~10 | |
氨氮温度系数 | 1.043 | 0~1.07 | |
20 ℃ 时溶解性有机磷的矿化率 | 0.1 | d−1 | 0~0.22 |
溶解性有机磷的温度系数 | 1.06 | 0~1.08 | |
20 ℃ 时的复氧系数 | 0.88 | d−1 | 0~10 |
氧碳比 | 1.22 | 0~2.67 | |
20 ℃ 反硝化速率常数 | 0.02 | d−1 | 0~0.09 |
反硝化的温度系数 | 1.04 | 0~1.04 |
表1. WASP模型主要动力学参数
基于2014~2015年水源站和出库站两个监测站的监测数据建立的WASP模型,根据2016~2017年水源站5项指标模拟出的出库站各项指标模拟值和实测值对比如图2所示,相对误差计算结果见表2、表3。由图2可知WASP模型对于各项水质指标的模拟值和实测值的重合度是较高的,趋势一致,偶有偏差,但是从整体上的趋势可以看到重合度很高,可见WASP模型对各项水质指标的拟合结果较为理想。由图2仍可见,桃林口水库6~9月的总磷、氨氮和高锰酸盐指数含量与其他月份相比相对较高,但6~9月时溶解氧含量与其他月份相比相对较低。可能是因为桃林口水库丰水期在6~9月,这段时期水库的水温较于其他时间较高,有利于微生物的繁殖生长,微生物的繁殖生长过程会消耗水中的溶解氧,致使溶解氧含量在丰水期较低,总磷、氨氮和高锰酸盐指数等浓度较高,因此桃林口水库6~9月的水质较于其他时期较差,这与王佳 [
由表2可知,总磷实测值与模拟值的平均相对误差为9.4%,其中相对误差在10%以下的有8组(共13组数据),占61.5%;硝酸盐氮实测值与模拟值的平均相对误差为3.84%,其中相对误差在5%以下的有8组(共13组数据),占61.5%。由表3可知,氨氮实测值与模拟值的平均相对误差为29.57%,其中相对误差在30%以下的有8组(共13组数据),占61.5%;高锰酸盐指数实测值与模拟值的平均相对误差为7.31%,其中相对误差在10%以下的有11组(共13组数据),占84.62%;溶解氧实测值与模拟值的平均相对误差为5.35%,其中相对误差在10%以下的有10组(共13组数据),占76.92%。由于在对桃林口水库水质进行监测的过程中可能会有人为误差的存在,所以从整体上来看WASP模型的模拟结果良好,该模型可用于桃林口水库水质预测与管理工作。且氨氮2016~2017年实测值的平均值为0.22 mg·L−1,其中小于0.2 mg·L−1的监测值有10组(共13组数据),占76.92%,可见氨氮浓度较低,对于浓度较低的指标监测时,监测仪误差较大 [
图2. 各水质指标实测值和模拟值对比图
时间 | 总磷 | 硝酸盐氮 | 氨氮 | 高锰酸盐指数 | 溶解氧 | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
年.月 | 实测值 | 模拟值 | 相对误差 | 实测值 | 模拟值 | 相对误差 | 实测值 | 模拟值 | 相对误差 | 实测值 | 模拟值 | 相对误差 | 实测值 | 模拟值 | 相对误差 |
2016.3 | 0.01 | 0.009 | −10.00% | 1.93 | 1.95 | 1.04% | 0.056 | 0.049 | −12.50% | 2.16 | 2.07 | −4.17% | 13.8 | 14.23 | 3.12% |
2016.5 | 0.01 | 0.01 | 0.00% | 1.71 | 1.75 | 2.34% | 0.17 | 0.125 | −26.47% | 2.38 | 2.321 | −2.48% | 12.8 | 11.08 | −13.44% |
2016.7 | 0.02 | 0.02 | 0.00% | 1.82 | 1.98 | 8.79% | 0.92 | 0.785 | −14.67% | 3.17 | 2.92 | −7.89% | 8.8 | 9.25 | 5.11% |
2016.9 | 0.028 | 0.03 | 7.14% | 3.85 | 3.65 | −5.19% | 0.117 | 0.101 | −13.68% | 2.69 | 2.751 | 2.27% | 6.45 | 6.9 | 6.98% |
2016.11 | 0.011 | 0.01 | −9.09% | 3.86 | 3.95 | 2.33% | 0.16 | 0.146 | −8.75% | 2.3 | 2.096 | −8.87% | 12 | 12.24 | 2.00% |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
2017.1 | 0.009 | 0.01 | 11.11% | 3.81 | 4.05 | 6.30% | 0.11 | 0.15 | 36.36% | 1.93 | 2.435 | 26.17% | 10.98 | 10.21 | −7.01% |
2017.3 | 0.011 | 0.01 | −9.09% | 4.08 | 4.2 | 2.94% | 0.17 | 0.1399 | −17.71% | 2.46 | 2.443 | −0.69% | 14.72 | 13.551 | −7.94% |
2017.5 | 0.009 | 0.01 | 11.11% | 4.03 | 3.82 | −5.21% | 0.24 | 0.2938 | 22.42% | 2.3 | 2.762 | 20.09% | 10.7 | 10.4758 | −2.10% |
2017.7 | 0.009 | 0.01 | 11.11% | 3.97 | 4.12 | 3.78% | 0.13 | 0.223 | 71.54% | 2.26 | 2.416 | 6.90% | 10.28 | 9.85 | −4.18% |
2017.8 | 0.012 | 0.01 | −16.67% | 3.97 | 3.82 | −3.78% | 0.43 | 0.2087 | −51.47% | 2.83 | 2.999 | 5.97% | 8.05 | 8.04 | −0.12% |
2017.9 | 0.012 | 0.01 | −16.67% | 3.77 | 3.84 | 1.86% | 0.19 | 0.1283 | −32.47% | 3.15 | 3.151 | 0.03% | 12.85 | 12.21 | −4.98% |
2017.1 | 0.011 | 0.01 | −9.09% | 3.49 | 3.67 | 5.16% | 0.16 | 0.0742 | −53.63% | 2.61 | 2.723 | 4.33% | 6.92 | 6.8432 | −1.11% |
2017.11 | 0.011 | 0.01 | −9.09% | 3.49 | 3.53 | 1.15% | 0.05 | 0.0614 | 22.80% | 1.9 | 1.999 | 5.21% | 8.3 | 9.25 | 11.45% |
表2. 总磷、硝酸盐、氨氮、高锰酸盐指数和溶解氧模拟结果误差分析(单位:mg∙L−1)
机制十分复杂,相较于模型公式而言实际水体更为复杂,使得实测值与模型公式产生差别所致的误差 [
r X Y = r X Y = ∑ i = 1 N ( X i − X ¯ ) ( Y i − Y ¯ ) ∑ i = 1 N ( X i − X ¯ ) 2 ∑ i = 1 N ( Y i − Y ¯ ) 2 (3)
式中:rXY为某水质指标实测值和模拟值的相关系数;Xi为该水质指标的实测值,mg·L−1; X ¯ 为该水质指标实测值的平均值,mg·L−1;Yi为该水质指标的模拟值,mg·L−1; Y ¯ 为该水质指标模拟值的平均值,mg·L−1。
由表3、图3可知,各项水质指标的判定系数R2均大于0.7,相关系数r均大于0.839;除高锰酸盐指数外,总磷、硝酸盐氮、氨氮和溶解氧的判定系数R2均大于0.883,相关系数r均大于0.9,其中总磷、硝酸盐氮和溶解氧的相关系数r均大于0.95,可见所建立的WASP模型模拟出的模拟值与实测值吻合度很高,具有较好的相关性。
总磷 | 硝酸盐 | 氨氮 | 高锰酸盐指数 | 溶解氧 | |
---|---|---|---|---|---|
判定系数R2 | 0.9632 | 0.9298 | 0.883 | 0.7042 | 0.93 |
相关系数r | 0.98142753 | 0.964261 | 0.939681 | 0.839166 | 0.964365 |
表3. 2016~2017年模拟结果分析
图3. 各水质指标的线性回归方程
运用WASP模型对桃林口水库的总磷、硝酸盐氮、氨氮、高锰酸盐指数和溶解氧等水质指标进行模拟研究。研究结果表明,该五项水质指标的模拟值与实测值吻合度很高,具有较好的相关性,判定系数R2均大于0.7,相关系数r均大于0.83,证明建立的WASP模型适用于桃林口水库,可以运用于桃林口水库水质预测与管理的工作中,同时为秦皇岛市青龙县水环境保护与管理工作提供决策性参考与依据。
感谢我的导师赵林老师、陈亮老师和我的朋友们冯琛雅和刘琦的无私帮助,希望我们能依旧开心乐观地度过每一天、越过每一个困难,做最善良舒服的自己。
项目编码:17273905D,项目名称:水质自动监测系统在地表水水库生态安全预警中的应用研究。
戴天骄,赵 林,陈 亮,冯琛雅,刘 琦,赵明杰. 基于WASP模型的桃林口水库水质模拟研究 Water Quality Simulation of Taolinkou Reservoir Based on WASP Model[J]. 水资源研究, 2019, 08(05): 499-507. https://doi.org/10.12677/JWRR.2019.85057