湖泊科学   2020, Vol. 32 Issue (5): 1539-1551.  DOI: 10.18307/2020.0525.
0

研究论文

引用本文 [复制中英文]

王昱, 卢世国, 冯起, 刘蔚, 刘娟娟, 刘开清, 左一锋, 梯级筑坝对黑河水质时空分布特征的影响. 湖泊科学, 2020, 32(5): 1539-1551. DOI: 10.18307/2020.0525.
[复制中文]
WANG Yu, LU Shiguo, FENG Qi, LIU Wei, LIU Juanjuan, LIU Kaiqing, ZUO Yifeng. Effect of cascade dam construction on the spatio-temporal variations of water quality in Heihe River. Journal of Lake Sciences, 2020, 32(5): 1539-1551. DOI: 10.18307/2020.0525.
[复制英文]

基金项目

国家自然科学基金项目(51669011)和甘肃省自然科学基金重大项目(18JR4RA002)联合资助

通信作者

王昱; E-mail:wangyu-mike@163.com

文章历史

2019-12-04 收稿
2020-02-29 收修改稿

码上扫一扫

梯级筑坝对黑河水质时空分布特征的影响
王昱1,2 , 卢世国1 , 冯起2 , 刘蔚2 , 刘娟娟1 , 刘开清3 , 左一锋1     
(1: 兰州理工大学能源与动力工程学院, 兰州 730050)
(2: 中国科学院西北生态环境资源研究院内陆河流域生态水文重点实验室, 兰州 730000)
(3: 甘肃省讨赖河流域水资源局, 酒泉 735000)
摘要:为探究梯级大坝建设对河流水质变化规律的影响,将黑河上中游划分为坝上河段、坝下河段及自然河段,于2017年12月-2018年8月选取了24个主要控制断面进行水质调查,并采用多元统计的方法对比分析了不同时空尺度上的水质分布特征.结果表明:黑河上中游水质时空变化的主要影响因子为水温(WT)、pH值、溶解氧(DO)、电导率(EC)、总氮(TN)、总磷(TP)和五日生化需氧量(BOD5).空间尺度上,WT、EC、BOD5、高锰酸盐指数(CODMn)、TN等指标具有显著性差异,其中坝上河段受BOD5、CODMn影响较大,自然河段WT、EC和TN为关键指标,而各个因子对坝下河段水质影响较小.时间尺度上,WT、EC、BOD5、氨氮与季节变化存在明显相关性,是不同河段水质时间变化的控制因子,且大多数水质因子在非汛期变化最明显.降水、温度、水文条件等季节性影响因素和梯级水库联合运用模式是该区域水质时间差异的主要原因;空间差异主要受祁连、张掖地区外源性污染物排放以及筑坝环境下水动力条件改变而产生的沉积滞留效应和沿程累积效应影响.研究表明,外源性污染源依然是导致水质变差的主要因素,梯级筑坝则是导致水质变差的间接因素.因此控制该区域人类活动所造成的外源性污染源,并针对不同种类污染物的季节变化特征实施合理的水库运行方式是改善水电梯级开发河段水质状况的关键.
关键词水质    多元统计    时空变化    梯级大坝    黑河    
Effect of cascade dam construction on the spatio-temporal variations of water quality in Heihe River
WANG Yu1,2 , LU Shiguo1 , FENG Qi2 , LIU Wei2 , LIU Juanjuan1 , LIU Kaiqing3 , ZUO Yifeng1     
(1: School of Energy and Power Engineering, Lanzhou University of Technology, Lanzhou 730050, P. R. China)
(2: Key Laboratory of Ecohydrology of Inland River Basin, Northwest Institute of Eco-environment and Resources, Chinese Academy of Sciences, Lanzhou 730000, P. R. China)
(3: Taolaihe Basin Water Resourse Authority of Gansu Province, Jiuquan 735000, P. R. China)
Abstract: Construction of casade dams in Heihe River had a certain influence on the river water quality. We divided the water into upstream and downstream reaches of dam and natural reaches. To explore the impact of casade dams on the water quality, multivariate statistical methods were used to analyze the spatio-temporal characteristics of water quality based on 24 different sampling sites from december 2017 to august 2018. The results showed that the main impact factors on water quality in Heihe River were water temperature (WT), pH, dissolved oxygen (DO), electrical conductivity (EC), total nitrogen (TN), total phosphorus (TP) and biochemical oxygen demand (BOD5). WT, EC, BOD5, chemical oxygen demand (CODMn) and TN varied signnificantly on the spatial scales. The water quality of the upstream reaches of dam were greatly affected by BOD5 and CODMn, while WT, EC and TN were the key drivers for the natural reaches and less environmental influence on the downstream water quality. In temporal terms, WT, EC, BOD5 and ammonia nitrogen excerted the greatest effects on the water quality, with most of the water quality factors changed obviously in non flood seasons. Factors such as temperature, rainfall, hydrological conditions and the joint operation mode of cascade reservoirs were the main factors affecting the seasonality of water quality in the region. The spatial differences were mainly influenced by the exogenous pollutants discharging in Qilian and Zhangye areas, and the influence of sediment retention and accumulation on the environment of dam construction due to the change of hydrodynamic conditions. The results show that exogenous pollution sources were still the main factors, followed by the cascade dam construction as the indirect factors, for the water quality deterioration. Therefore, controlling the exogenous pollution sources caused by human activities in the region, and implementing a reasonable reservoir operation mode according to the seasonal variation of pollutants were the key to improve the river water quality along with the cascade hydropower development.
Keywords: Water quality    multivariate analysis    spatio-temporal variation    cascade dams    Heihe River    

河流水质是流域的一个重要特征,通过研究水体中物理化学成分的变化,不仅可以获取地质、气象、水文等自然因素在河流中形态表征的信息,同时还能体现流域范围内社会经济、人类活动对河流水系作用的响应[1].近年来,随着流域梯级开发,原始河流的水力学连通性和纵向联通性发生改变,在很大程度上重新配置了上下游河道物质、能量和水流输送的时空分布特征,进而影响到河流的水文过程、水动力条件、介质运移机制等水文环境要素,这些将成为水体富营养效应的叠加因素[2-3];另外,河流水动力条件的变化势必会影响沉积物粒度的组成和分布[4],而河流沉积物作为水体污染物(如营养盐颗粒)的重要载体,其粒度特性的改变将会制约水体理化物质的迁移转化行为.因此,开展水电梯级开发情景下河流水质变化规律研究,对于流域生态修复和水库调度具有重要意义.目前,国内外学者针对水电开发对河流水质变化的影响进行了大量研究. Dynesius等[5]研究表明,水库建设将导致河流水动力条件减弱,阻碍了颗粒态营养要素的迁移方式,从而加重了库区水体富营养化趋势;Petts[6]认为水库水质发生变化的主要影响因素是水库运用方式、水位、分层、滞留时间、异重流等;刘丛强等[7]分析表明,太阳辐射和热量传输不平衡会导致季节性水体分层,主要体现在库区水温结构分层和下泄水温的延迟滞留效应上;此外,国内外不少学者[8-10]认为梯级水库运行对水质的影响必然会以某种形式叠加,产生累积效应.但上述研究大都集中在长江流域、黄河流域、澜沧江流域、Seine河流域等外流河及其支流上,对内陆河流域这样水质脆弱地区的案例性研究并不多见.

黑河作为我国第二大内陆河,是河西走廊绿洲赖以生存和社会经济可持续发展的重要水资源基地.近年来,黑河上游水电梯级开发极大地破坏了河流连续性,不仅使库区自身水体理化性质的时空分布发生变化,还对下游河道水环境系统的演替以及河流生态健康带来一系列影响[11-12].纵观以往研究,多是简单评价各区间河段的水体质量[13-14],而缺乏对筑坝河段水质大尺度、多断面的时空变化特征研究.为此,本文通过监测2017年12月—2018年8月汛期和非汛期的水质数据,运用主成分分析(PCA)和方差分析(ANOVA)的方法,识别水质因子的时空差异性,探讨导致该区域水质时空分异的原因以及水电梯级开发对河流水质的影响,以期为流域综合治理和水库调度提供科学依据.

1 材料与方法 1.1 研究区概况

黑河发源于祁连山北麓,流域范围介于37°45′~42°40′N,96°42′~102°04′E之间,流经青海、甘肃和内蒙古3个省区,全长821 km,流域总面积14.3×104 km2.黑河以莺落峡和正义峡为界,划分为上、中、下游,莺落峡以上为上游,主要由两大支流野牛沟河和八宝河及其汇合后形成的部分干流构成,其干流穿越祁连山谷后至莺落峡出山流入中游走廊平原,经正义峡流出,进入下游额济纳,最终汇入居延海.

本文选择研究黑河上中游地区,其中上游海拔在2600~4300 m之间,多年平均气温为-5~4℃,年降水量超过350 mm,部分区域可达600~700 mm,冰川融水量约为4.0×108 m3,气候阴湿寒冷,植被覆盖较好,是黑河流域的主要产流区.中游地区年降水量不足200 mm,蒸发量在1400 mm左右,区内光热资源充足,为黑河径流的主要利用区.黑河属于典型的内陆河,径流来源主要有降水形成的地表径流、冰川融冰、融雪和地下水等,其中降水补给约为95 %,冰川融水补给仅为4 %. 6—8月份是河源径流的丰水期,该时段径流量约占全年径流量的69.6 %,枯水期12—3月份来水量约占全年总量的15.9 %左右[15-16],研究区其他水文水力要素特征见表 1.

表 1 黑河上中游河道及水文要素特征* Tab. 1 Channel characteristics and hydrological features of the upper and middle Heihe River

黑河流域水电梯级开发主要分布在上游干流山区,其水能资源蕴藏量106×104 kW,可开发量为52.8×104 kW,年发电量为38.48×108 kW ·h.根据地形情况,干流从上到下共规划了8座梯级电站,依次为黄藏寺、宝瓶河、三道湾、二龙山(松木沟)、大孤山、小孤山(石羊岭)、龙首二级(西流水)、龙首一级以及龙渠电站;两支流末端分别建有地盘子和牛板筋电站(图 1).总体上该河段是以大孤山、宝瓶河两大电站为主,地盘子、牛板筋等小水电站为辅的开发格局.截止所有样品采集完成时,干流河段已建成7座梯级电站,支流河段建成2座(表 2).

图 1 研究区位置及电站分布 Fig.1 Study area and hydropower stations distribution
表 2 黑河上游水电开发情况 Tab. 2 Hydropower development on the mainstream of the upper Heihe River
1.2 样品采集与分析

水库蓄水运行后,一方面库区水位雍高,水深加大,流速减缓,加之受水库回水的顶托作用影响,阻碍了河流污染物向下流动,容易使坝上河段出现水体富营养化现象[17-18];另一方面,水库水位的调节以及来流水体理化性质的差异,造成水质因子处于频繁的交换状态,并且这种交换又会反过来进一步改变坝上和坝下河道水体的物理化学性质[19].另外,通常上游河道水力特性的改变会对下游水环境系统演变带来一系列影响[12].因此,为探明梯级水电开发对黑河水质变化的影响,本文在筑坝河流上布设坝上河段、坝下河段和下游自然河段3类采样区域,选择24个主要控制断面进行水样监测(图 1).自然河道上根据黑河干流水文站的分布情况、河流交汇位置及行政区的划分属性,设置6个采样断面(H19、H20、H21、H22、H23、H24);坝上河段依据梯级电站河流入库口位置及库容大小,设置9个采样断面(H1、H2、H5、H6、H9、H10、H12、H14、H16);坝下河段按照梯级水库泄流情况以及电站尾水出口位置,设置9个采样点(H3、H4、H7、H8、H11、H13、H15、H17、H18).分别于河流非汛期(2017年12月—翌年3月)、汛期(2018年6—8月)进行逐月采样,其中非汛期采样为每月下旬,汛期为每月上旬,取样过程按照从上到下的原则分3天进行,每个断面采样3次,分别采集1000 mL水样固定后低温保存运回实验室;另外,各采样时段并没有较大的降水.

监测指标包括水温(WT)、pH值、溶解氧(DO)、电导率(EC)、五日生化需氧量(BOD5)、高锰酸盐指数(CODMn)、氨氮(NH3-N)、总氮(TN)和总磷(TP),共计9项.其中WT、pH、DO、EC使用哈希便携式水质仪现场测定,剩余指标参照《水和废水监测分析方法》[20]中所示方法在实验室测定,并且每个样品皆平行测定3次,计算结果取其平均值.

1.3 数据分析方法

本文首先运用主成分分析法[21]提取了引起研究区域水质变化的主导因素并加以分析,再结合重复测量方差分析[21]判别黑河坝上河段、坝下河段与自然河段3种不同空间类别上的水质指标差异,最后利用单因素方差分析[21]识别不同空间类别上时间尺度的水质变化特征,从而在时空尺度上探讨梯级水库联合运行对河流水质的变化规律影响机制.主成分分析之前,先对原始数据进行KMO和Barlett球度检验,同时为了消除各指标间不同量纲的影响,将原数据进行标准化处理[22],即均值为0,方差为1.上述所有分析均在Microsoft Excel 2010、IBM SPSS Statistics 20.0和Origin 9.0软件中实现.

2 结果与分析 2.1 水质因子特征与评价

2017年12月—2018年8月研究区域水质因子的均值和实测范围见表 3.由表 3可知,监测时段内该区域水体pH值主要分布在8.24~8.99之间,属于弱碱性水;DO、BOD5、CODMn、NH3-N、TN和TP浓度均值分别为8.26、2.71、2.87、0.47、0.94和0.17 mg/L,根据国家《地表水环境质量标准》(GB 3838—2002)常规水质项目评价标准,TN和TP属于Ⅲ类水质标准,BOD5、CODMn和NH3-N属于Ⅱ类水质标准,DO属于Ⅰ类水质标准,基本能满足功能区目标限值.从实测范围来看,宝瓶至龙首河段的TP浓度、宝瓶水库和大孤山水库的BOD5浓度、平川和高台地区的TN浓度均超出了功能区Ⅲ类水标准限值,其中TN浓度最大值超出Ⅴ水标准的1.07倍,TP浓度最大值超出湖库Ⅴ类水标准的1.44倍,说明上述水域污染较重,甚至部分河段出现了劣Ⅴ类水现象.通过Pearson相关分析得出(表 4),EC与除pH、DO、TP外的其他因子呈极显著正相关,即EC随着NH3-N、TN、BOD5、CODMn浓度的增加而显著增加;BOD5、CODMn、NH3-N、TP之间有显著相关性,表明其来源具有一致性和相似性;DO与TP、NH3-N呈明显负相关,说明它们之间是相互影响的;WT与EC、TN、NH3-N、BOD5也表现出较强的相关性,体现了自然因素水温对河水的影响.

表 3 研究区水质因子特征 Tab. 3 Characteristics of water quality factors in the studied area
表 4 水质因子间的Pearson相关系数* Tab. 4 Pearson correlation coefficients of water quality indexes
2.2 水质主成分分析

经检验得出,本文用于分析的数据KMO值为0.67,Bartlett球度的显著性满足P < 0.001的置信度要求,表明原始数据之间有较强的相关关系,适合用于主成分分析.按照特征值λ>1的原则,提取了3个引起水质变化的主成分(表 5),累积方差贡献率为81.949 %,能够充分提取原始数据信息.

表 5 研究区监测指标主成分分析结果* Tab. 5 Principle component analysis of monitoring index in the studied area

表 5可知,第1主成分PC1的方差贡献率最大,为53.364 %,与之相关联的参数是WT、EC、BOD5、CODMn和TN,体现了气候条件(WT)变化对水质的影响,以及人类活动影响下有机污染物(BOD5、CODMn)和营养盐(TN)主导着该区域水质变化;且TN对PC1有最强的相关性,变化特征最明显,表明其受季节变换、人类活动等因素影响显著.第2主成分PC2的方差贡献率为16.671 %,与之相关联的参数是NH3-N、TP和DO,并且NH3-N与TP呈正相关,与DO呈负相关,同时Pearson相关分析也显示它们之间表现出显著负相关,说明随着营养盐(TP、NH3-N)的输入,DO浓度出现相反的变化.第3主成分PC3的方差贡献率最低,为11.931 %,与之关联的参数是pH,反映了水体的酸碱度,它对水体的氧化还原反应起着一定的控制作用,决定着河流水体系统的化学稳定性[23].

为了更清晰地反映研究区域水质状况,将各采样点得分与水质参数荷载矩阵进行PCA双序分析(图 2),可以直观反映3方面信息:①样点在排序空间上呈明显的季节分布,逆时针方向来看,汛期基本聚集在第2象限和第3象限,非汛期在第1象限和第4象限. ②水质因子的箭头越长表示因子荷载越大,对排序的贡献也越大,箭头间的夹角则表示水质因子间的相关性,夹角越小相关性越高;图中箭头较长的水质因子在排序空间分成WT和EC、TN、CODMn、NH3-H、BOD5、DO两簇排列,主导着该区域的水质变化,且每簇因子间夹角越小,相关性越高. pH、TP箭头短,对排序贡献小,样品间差异小. ③将采样点向水质因子箭头及其延长线作投影,投影位置在箭头指向方向越远数值越大,表明该水质因子对其正贡献越大,反之则越小.逆时针方向来看,1、3象限中,箭头指向了排序空间中EC、TN、CODMn、NH3-H、BOD5、DO变异最大的方向,反映了非汛期BOD5、TN、CODMn、NH3-H主导着该时段的水质变化,并引起了DO和EC的变化;2、4象限中WT和pH最高,说明WT和pH主导着汛期各采样点的水质变化.

图 2 研究区水质因子PCA分析双序图(X表示汛期样点,F表示非汛期样点) Fig.2 The PCA biplots of water quality factors in the research area
2.3 水质变化的时空分布特征 2.3.1 水质因子的空间差异性分析

本文将研究区域划分为坝上断面、坝下断面和自然断面3组,对2017年12月—2018年8月的水质数据进行重复测量方差分析,结果表明,WT、EC、BOD5、CODMn、TN等水质因子在不同组别上存在显著性差异(P < 0.05).从它们的空间变化规律可以看出(图 3),WT、EC、TN的最大平均值出现在自然断面,显著高于组间不存在明显差异的坝上断面和坝下断面. BOD5最大平均值出现在坝上断面,最小平均值出现在坝下断面,表现出组间显著性差异. CODMn最小平均值出现在坝下断面,显著低于不存在显著性差异的坝上断面和自然断面.其中,TN和EC浓度值在坝下断面较为均匀集中,而坝上断面的CODMn浓度值以及自然断面的WT、EC、BOD5浓度值则相对分散.

图 3 研究区水质因子的空间差异 Fig.3 Spatial variations of water quality factors in the research area
2.3.2 水质因子的时间差异性分析

由于黑河流域年降雨量的70 %都集中分布在6—8月[24].据此,本文将该区域的6—8月划分为汛期,其余月份(12月—翌年3月)归类为非汛期.水质参数的时间变化采用季节参数相关矩阵进行评价,通过单因素方差分析(One-way ANOVA)得出,除pH、DO、TP之外,WT、EC、BOD5、CODMn、TN和NH3-N与季节变化存在显著相关性(P < 0.05).这些显著性指标中,TN具有最大的相关系数(Spearman's R=0.701),其次是EC(R=0.700)、BOD5(R=0.692)、WT(R=0.656)和NH3-N(R=0.552).而与季节相关的参数可视为引起水质时间变化的主要指标.在不同空间河段上,随着季节的变化,EC、BOD5、CODMn、NH3-N、TN浓度整体呈上升趋势.其中EC在各河段变幅较小;BOD5和NH3-N浓度在各河段的变幅均较大,而坝上河段变幅最大;NH3-N和TN在自然河段变幅最大,坝上河段与坝下河段变幅较小;WT值无论是在坝上河段,或者是坝下河段,还是自然河段,汛期都要明显高于非汛期,体现了显著的气候区域差异.另外,本文研究发现DO浓度并没有呈现出显著的季节变化,且pH和TP浓度在汛期和非汛期差异也不明显.总体来看,研究区各水质因子在非汛期变化较汛期显著(图 4).

图 4 研究区水质因子的时间差异 Fig.4 Temporal variations of water quality factors in the research area
3 讨论 3.1 梯级大坝建设下库区河段水质变化特征分析

河流筑坝蓄水后,强水动力条件下的河流搬运作用将逐渐演变为弱水动力条件下的湖泊沉积作用,进而增强了水库对河流沉积物、营养盐颗粒等物质的沉积滞留效应,导致库区水体的物理化学特性发生改变[25-26].有研究表明[7],水库库容也可能是其影响因素之一.据此次调查结果显示,黑河库区河段上NH3-N、TP、EC、BOD5、CODMn浓度均随库容的增大而升高.主要原因是水库库容越大,与之相适应的调节性能越好,水力停留时间越长,因而营养盐、有机物等各种污染物的扩散和降解程度也就越缓慢.通常营养盐、有机物的滞留会消耗水中大量DO[27],本文中该河段水体DO浓度对BOD5、CODMn、TP浓度的变化也呈现出相似的响应规律(图 5),表明有机物和营养盐浓度是影响库区河段水体DO浓度变化的主要因素.有观点认为[28-29],梯级水库联合运行会将大坝阻隔对河流水环境系统的影响逐级“放大”,呈现出迥异于单一水库的变化特征,如水温的延迟滞留效应逐渐加大,各类离子、污染物浓度逐渐升高,水体富营养化趋势逐渐加重等.从本文研究结果来看,库区河段上NH3-N、TP、BOD5、EC都随距宝瓶坝前距离的减小而逐渐增大,且距坝前距离分别能解释其浓度变化的67.7 %、73.6 %、68.9 %和77.1 %,说明上游外源性污染源产生的EC、NH3-N、TP、BOD5浓度在下游筑坝环境下有较明显的沿程累积效应.但这些污染物浓度却与下游龙首一级水库坝前距离相关性很差(图 6),只能解释其浓度变化的6.9 %、9.0 %、8.4 %和28.0 %.表明河流污染物经梯级水库调节后(即沿程累积效应),各级河段水质污染程度均有所加重,且其污染程度类似.黑河梯级水库属于深水峡谷型水库,由于年际气温变化较小,水温分层变化不是特别明显,没有典型的温跃层,但正是这种不显著的温度分层结构,有效限制了水库上下水团的对流交换,阻断了下层水体中溶解氧的恢复补充途径[30],从而使水库下层形成了还原性环境,抑制了亚硝化作用和硝化作用,导致沉积物中氨化作用形成的NH3-N不能完全转化为硝态氮,进而造成下泄水体的NH3-N浓度升高,这是坝上河段和坝下河段NH3-N浓度没有差异的重要原因.筑坝蓄水对河流水温的另一重要影响则是滞温效应,具体表现为下泄水温出现平坦化、均一化特性[8].瞿群[31]、杨宏[32]分析发现黑河高坝水库建成后下泄水温常年维持在较稳定的低温状态,且秋、冬季有增温效应,春、夏季有降温效应,本文中相应河段水温的时空变化过程也有此类似现象.总体来看,黑河梯级水库建设对河流水质的影响具有三维结构的特征,即横向、纵向和垂向3个方向.横向上由于库区自身的沉积作用取代了河流搬运作用,营养盐、有机物等典型污染物易于发生沉积滞留效应,加大了库区水体富营养化的潜在风险,其程度跟库容大小有一定关系;纵向上受梯级水库联合运行影响,大坝阻隔将河流污染物的累积效应逐级放大,加重了各级库区的污染趋势;而垂向上水温分层改变了水团混合性质,干扰了水体中一些离子的迁移转换过程,并对下泄水体理化特性也产生了一定的影响.

图 5 水质因子与水库库容的关系 Fig.5 The relationship between water quality indexes and reservoir capacity
图 6 水质因子与坝前距离的关系 Fig.6 The relationship between water quality indexes and distance from dam
3.2 不同河段水质变化特征对梯级大坝建设的响应

从空间尺度上看,黑河梯级水库蓄水运行以后,坝上河段、坝下河段及下游自然河段水质的空间差异主要体现在WT、EC、BOD5、CODMn、TN等指标上.其中BOD5、CODMn浓度在坝上河段最高,坝下河段最低,而BOD5、CODMn是常被用来反映水体受工农业、生活等有机耗氧类物质污染程度的综合指标.考虑到上游祁连地区金属矿采选、冶炼等工矿业活动较强,因此上游来水中外源性污染源排放的有机物沿河累积于库区,并经水库调节后,上覆水体受污染较重,下层水体受影响较小,这虽然在一定程度上降低了下泄水体中污染物的本底值,但当遇到大规模降雨或洪水事件,储存于坝内的有机物被冲刷流出后,下游水体将受到严重污染. TN浓度在自然河段较高,主要原因是该河段途径张掖、临泽、高台等市县,人口稠密,各地区的工业废水以及村镇生活污水、农业径流等人类活动向河水排入了大量营养盐,导致TN质量浓度急剧升高,加之自然河段河床比降较小,水流缓慢,而且受上游筑坝影响,来流水动力条件减弱,故污染物质更容易聚集. EC跟TN情况类似,在自然河段最高. Pearson相关分析显示,随着水体盐类污染物的增加,EC也升高.可能是由于该地区造纸、农副产品加工企业产生的工业废水以及村镇生活污水等外源性点源污染向河水排放营养盐的同时,带入了大量溶解性离子,从而使EC升高[30],这也进一步反映了自然河段水质受人类活动(包括水土流失)影响较大.基于前文分析,各河段WT的变化主要与气候、海拔、人类活动等因素有关.值得一提的是,对于磷元素而言,水库的拦截效率较高,可高达80 %,是磷元素的一个有效“汇”[10, 33].有研究表明[34],猫跳河梯级水电开发对TP的滞留效果显著,约使磷酸盐浓度降低了90 %;Seine河上游水库滞留了60 %的磷元素,主要是沉降在沉积物中[31].本研究中水体TP浓度相对较高,但并没有明显的空间分布差异.一方面原因是水库大量拦沙拦磷破坏了水体纳污的缓冲机制,造成坝下河段的溶解性磷浓度增加[35];另一方面,上游支流区内强烈的放牧活动产生较多的磷元素,导致输入河流的TP浓度较高,而中游自然河段周围是农灌区,区内农业非点源污染、化肥和农药的使用以及雨水冲刷携带的颗粒磷使水体TP浓度增加.同时值得注意的是,宝瓶水库TP浓度远高于其他水库,主要原因是宝瓶电站位于黑河梯级水库的最上游,直接来源于上游支流区内的大量磷元素,被沉积滞留于库区底泥后,无论是在厌氧或好氧环境下,磷元素的释放都随温度的升高而增长[36].

时间尺度上,由于黑河梯级水库以季节性调度方式为主,故季节变化可能对水质变化影响更大.单因素方差分析表明,不同河段绝大多数水质因子随时间变化呈现出季节变化的特征,如非汛期EC、BOD5、CODMn、NH3-N、TN浓度显著高于汛期,这也与PCA双序分析结果基本一致.不同季节水质变化主要受降水、温度、水文条件、水库运用方式等因素影响[37].其中非汛期(12—3月)是黑河径流最小的季节,水库水量补给和换水频度减弱,水力停留时间延长,EC、BOD5、CODMn、NH3-N等典型污染物容易滞留;加之非汛期(冬季)气温较低,微生物活动减弱,水体自净能力低下,最终导致该时段河流污染负荷较高.而与此相反的是,汛期(6—8月)上游来水充足,水库弃水较多,污染物一起被排往下游的也多;另外,汛期降水的增加对污染物产生了较强的稀释净化作用,同时水温增加,流域及河道内微生物活动增强,能够迅速降低水体中污染物的浓度,因此汛期河流污染程度相对较轻.

纵向来看,杜有军等[38]的研究表明2001—2008年黑河自然河段水质有逐年加重的趋势,具体表现为CODMn、NH3-N浓度总体上升,其浓度年均值分别为5.04和0.80 mg/L;而2017—2018年监测结果为2.77和0.46 mg/L,呈现下降趋势,这从侧面反映出近年来流域范围内对企业废水、生活污水等点源污染的治理有一定效果,但本次监测数据也表明TN浓度有所上升,分析认为是自然河段水质仍受农业面源污染影响.杨永宇等[13]和李慧[39]研究指出2007—2011年黑河筑坝河段(札马什克至莺落峡河段)CODMn浓度变化不明显,TP、NH3-N浓度略有下降.但从2017—2018年监测结果来看,CODMn、TP、NH3-N浓度均有上升趋势,说明随着蓄水进程的推进,水质变化特征也发生了相应改变.

4 结论

1) 根据国家《地表水环境质量标准》(GB 3828—2002)常规水质项目评价标准,黑河上中游全年水质类别基本能满足功能区目标限值,但TP浓度及TN浓度在局部地区出现了超标现象.主成分分析表明该区域水质时空变化的主要影响因素为WT、pH、DO、EC、TN、TP和BOD5.

2) 水质因子在时空分布上有一定的差异性.空间尺度上,WT、EC、BOD5、CODMn、TN等指标具有显著性差异,其中坝上河段受BOD5、CODMn影响较大,自然河段WT、EC和TN为关键指标,而各因子对坝下河段水质影响较小.时间尺度上,WT、EC、BOD5、NH3-N与季节变化存在明显相关性,是不同河段水质时间变化的控制因子,且大多数水质因子在非汛期变化最明显.

3) 水质因子的空间差异主要受祁连、张掖地区外源性污染物排放以及筑坝环境下水动力条件改变而产生的沉积滞留效应和沿程累积效应的影响;时间差异主要受降水、温度、水文条件等季节性影响因素和梯级水库联合运用模式的影响.

4) 研究表明,外源性污染源依然是导致水质变差的主要因素,梯级筑坝是导致水质变差的间接因素.因此控制该区域人类活动所造成的外源性污染源,并针对不同种类污染物的季节变化特征实施合理的水库运行方式是改善水电梯级开发水域水质状况的关键.

5 参考文献

[1]
Li K, Li ZH, Chen HB et al. Spatial variations of the water quality in Jianhe River of Wudang Mountain in the upper reaches of Danjiangkou Reservoir. J Lake Sci, 2013, 25(5): 649-654. [李昆, 李兆华, 陈红兵等. 丹江口水库上游武当山剑河水质空间差异性分析. 湖泊科学, 2013, 25(5): 649-654. DOI:10.18307/2013.0505]
[2]
Mao ZP, Wang YC, Peng WQ et al. Advances in effects of dams on river ecosystem. Advances in Water Science, 2005, 16(1): 134-140. [毛战坡, 王雨春, 彭文启等. 筑坝对河流生态系统影响研究进展. 水科学进展, 2005, 16(1): 134-140.]
[3]
Yang K, Deng X, Li XL et al. Impacts of hydroelectric cascade exploitation on river ecosystem and landscape:a review. Chinese Journal of Applied Ecology, 2011, 22(5): 1359-1367. [杨昆, 邓熙, 李学灵等. 梯级开发对河流生态系统和景观影响研究进展. 应用生态学报, 2011, 22(5): 1359-1367.]
[4]
Wang Y, Lian YT, Feng Q et al. Effects of dam interception on the spatial distribution of sediment granularities in Heihe River. J Lake Sci, 2019, 31(5): 1459-1467. [王昱, 连运涛, 冯起等. 筑坝拦截对黑河河道沉积物粒度空间分布的影响. 湖泊科学, 2019, 31(5): 1459-1467. DOI:10.18307/2019.0505]
[5]
Dynesius M, Nilsson C. Fragmentation and flow regulation of river systems in the Northern third of the world. Science, 1994, 266(5186): 753-762. DOI:10.1126/science.266.5186.753
[6]
Petts G ed. Impounded rivers:perspectives for ecological management. New York: Wiley, Chichebster, 1984: 67-71.
[7]
Liu CQ, Wang FS, Wang CY et al. Responses of aquatic environment to river damming-from the geochem ical view. Resources and Environment in the Yangtze Basin, 2009, 18(4): 384-396. [刘丛强, 汪福顺, 王雨春等. 河流筑坝拦截的水环境响应——来自地球化学的视角. 长江流域资源与环境, 2009, 18(4): 384-396.]
[8]
Ji DB, Long LH, Xu H et al. Advances in study on cumulative effects of construction of cascaded reservoirs on water environment. Advances in Science and Technology of Water Resources, 2017, 37(3): 7-14. [纪道斌, 龙良红, 徐慧等. 梯级水库建设对水环境的累积影响研究进展. 水利水电科技进展, 2017, 37(3): 7-14.]
[9]
Ma ZL, Cheng GW. Progress in research on impacts of river hydro-power engineering on eco-environment. Advances in Water Science, 2006, 17(5): 748-753. [麻泽龙, 程根伟. 河流梯级开发对生态环境影响的研究进展. 水科学进展, 2006, 17(5): 748-753.]
[10]
Jossette G, Leporcq B, Sanchez N et al. Biogeochemical mass-balances (C, N, P, Si) in three large reservoirs of the Seine Basin(France). Biogeochemistry, 1999, 47(2): 119-146.
[11]
Wang Y, Lu SG, Feng Q et al. Spatio-temporal characteristics and source identification of water pollutants in the upper and middle reachers of Heihe River. China Environmental Science, 2019, 39(10): 4194-4204. [王昱, 卢世国, 冯起等. 黑河上中游水质时空分异特征及污染源解析. 中国环境科学, 2019, 39(10): 4194-4204.]
[12]
Li Z, Chen YB, Li C et al. Advances of eco-environmental effects and adaptive management in river cascading development. Advances in Earth Science, 2018, 33(7): 675-686. [李哲, 陈永柏, 李翀等. 河流梯级开发生态环境效应与适应性管理进展. 地球科学进展, 2018, 33(7): 675-686.]
[13]
Yang YY, Yin L, Liu C et al. Water quality assessment of Heihe River Basin based on grey correlation method and BP neural network. Yellow River, 2017, 39(6): 58-62. [杨永宇, 尹亮, 刘畅等. 基于灰关联和BP神经网络法评价黑河流域水质. 人民黄河, 2017, 39(6): 58-62.]
[14]
Li H, Zhou YC. Water quality evaluation of Heihe River Basin based on fuzzy matter-element model. Yellow River, 2015, 37(10): 78-80, 85. [李慧, 周轶成. 基于模糊物元模型的黑河流域水质评价. 人民黄河, 2015, 37(10): 78-80, 85.]
[15]
Guo QL, Yang YS, Chang XS et al. Annual variation of Heihe River runoff during 1957-2008. Progress in Geography, 2011, 30(5): 550-556. [郭巧玲, 杨云松, 畅祥生等. 1957-2008年黑河流域径流年内分配变化. 地理科学进展, 2011, 30(5): 550-556.]
[16]
Yu ZS, Bai F, Li WP. Characteristic and cause analysis of water resources evolution in the middle reaches of the Black River Basin. China Rural Water and Hydropower, 2007(6): 1-3, 7. [余志山, 白福, 李文鹏. 黑河流域中游地区水资源演变特征及成因分析. 中国农村水利水电, 2007(6): 1-3, 7.]
[17]
Cheng B, Wang X, Ma JC et al. Spatial Distribution characteristics of nutrients and chlorophyll-a in the Lancang River Basin under cascade Reservoirs. Environmental Science, 2019, 40(4): 1779-1787. [程豹, 望雪, 马金川等. 澜沧江流域梯级水库建设下水体营养盐和叶绿素a的空间分布特征. 环境科学, 2019, 40(4): 1779-1787.]
[18]
Ran GH, Ge JW, Miao WJ et al. Spatial-temporal differentiation of water quality in Gufu River of Three Gorges Reservoir. Acta Ecologica Science, 2013, 33(17): 5385-5396. [冉桂花, 葛继稳, 苗文杰等. 三峡库区古夫河水质时空分异特征. 生态学报, 2013, 33(17): 5385-5396.]
[19]
Ji DB, Liu DF, Yang ZJ et al. Hydrodynamic characteristics of Xiangxi Bay in Three Gorges Reservoir. Science China Press, 2010, 40(1): 101-112. [纪道斌, 刘德富, 杨正健等. 三峡水库香溪河库湾水动力特性分析. 中国科学:物理学力学天文学, 2010, 40(1): 101-112.]
[20]
Editorial Board of "Monitoring and analysis methods of water and wastewater", Ministry of Environmental Protection of the People's Republic of China ed. Monitoring and analysis methods of water and wastewater: fourth edition. Beijing: China Environmental Science Press, 2002. [国家环境保护总局《水和废水监测分析方法》编委会.水和废水监测分析方法: 第4版.北京: 中国环境科学出版社, 2002. ]
[21]
Li YG, Du HW, Huang YS et al eds. SPSS 19.0 introduction and improvement of statistical analysis. Beijng: Tsing hua University Press, 2014: 308. [李玉光, 杜宏巍, 黄永生等. SPSS19. 0统计分析入门与提高.北京: 清华大学出版社, 2014: 308. ]
[22]
Sun GY, Liu X, Li LH. Spatial distributions of water physicochemical properties in the Kaidu River. Arid Zone Research, 2017, 34(2): 259-265. [孙桂燕, 刘翔, 李兰海. 开都河水理化性质的空间分布特征. 干旱区研究, 2017, 34(2): 259-265.]
[23]
Ma XX, Wang LC, Liao LL. Spatio-temporal characteristics and source identification of water pollutants in Wenruitang River Watershed. Environmental Science, 2015, 36(1): 64-71. [马小雪, 王腊春, 廖玲玲. 温瑞塘河流域水体污染时空分异特征及污染源识别. 环境科学, 2015, 36(1): 64-71.]
[24]
Liu SY, Zhang GH, Zhang CY et al. Water resources system evolution and impact of human activities in Heihe River Basin. Journal of Jilin University:Earth Science Edition, 2008, 38(5): 806-812, 819. [刘少玉, 张光辉, 张翠云等. 黑河流域水资源系统演变和人类活动影响. 吉林大学学报:地球科学版, 2008, 38(5): 806-812, 819.]
[25]
Hart DD, Poff NL. A special section on dam removal and river restoration. BioScience, 2009, 52(8): 653-655.
[26]
Li GR, Liu CQ, Chen C et al. Dissolve inorganic carbon and its carbon isotope composition in cascade reservoir of the Maotiao River during summer and autumn. Environmental Science, 2009, 30(10): 2891-2897. [李干蓉, 刘丛强, 陈椽等. 猫跳河流域梯级水库夏-秋季节溶解无机碳(DIC)含量及其同位素组成的分布特征. 环境科学, 2009, 30(10): 2891-2897.]
[27]
Yang FX, Jian HM, Tian L et al. Distribution characteristics of COD and DO and its influencing factors in the daliaohe estuary. Environmental Science, 2014, 35(10): 3748-3754. [杨福霞, 简慧敏, 田琳等. 大辽河口COD与DO的分布特征及其影响因素. 环境科学, 2014, 35(10): 3748-3754.]
[28]
Stanford JA, Ward JV. Revisiting the serial discontinuity concept. River Research & Applications, 2001, 17(4): 303-310.
[29]
Zhang SJ, Yan JP, Li GQ. Cumulative effects of cascade development project adjustment on water temperature. Journal of Hydraulic Engineering, 2014, 45(11): 1336-1343. [张士杰, 闫俊平, 李国强. 流域梯级开发方案调整的水温累积影响研究. 水利学报, 2014, 45(11): 1336-1343.]
[30]
Wang CY, Zhu J, Ma M et al. Thermal stratification and paroxysmal deterioration of water quality in a canyon reservoir, Southwestern China. J Lake Sci, 2005, 17(1): 54-60. [王雨春, 朱俊, 马梅等. 西南峡谷型水库的季节性分层与水质的突发性恶化. 湖泊科学, 2005, 17(1): 54-60. DOI:10.18307/2005.0109]
[31]
Qu Q. Research on the river hydropower engineering PEIA-in the case of Heihe River (Huangzangsi reach to Yingluoxia) hydropower engineering PEIA[Dissertation]. Lanzhou: Lanzhou University, 2007. [瞿群.流域水电规划环境影响评价研究——以黑河干流(黄藏寺~莺落峡河段)水电梯级开发规划环境影响评价为例[学位论文].兰州: 兰州大学, 2007. ]
[32]
Yang H. Study on cumulative effects assessment of cascade hydropower development in drainage area[Dissertation]. Lanzhou: Lanzhou Uinverity, 2007. [杨宏.流域水电梯级开发累积环境影响评价研究[学位论文].兰州: 兰州大学, 2007. ]
[33]
Ran XB, Yu ZG, Yao QZ et al. Advances in nutrient retention of dams on river. J Lake Sci, 2009, 21(5): 614-622. [冉祥滨, 于志刚, 姚庆祯等. 水库对河流营养盐滞留效应研究进展. 湖泊科学, 2009, 21(5): 614-622. DOI:10.18307/2009.0502]
[34]
Li YQ, Lu CX, Deng O et al. Regional response of ecological effects of cascade reservoirs development within the Maotiaohe River Watershed. Resources Science, 2011, 33(8): 1454-1461. [李亦秋, 鲁春霞, 邓欧等. 猫跳河流域梯级开发的生态效应区域响应. 资源科学, 2011, 33(8): 1454-1461.]
[35]
Zhou JJ, Zhang M, Li Z. Dams altered Yangtze River phosphorus and restoration countermeasures. J Lake Sci, 2018, 30(4): 865-880. [周建军, 张曼, 李哲. 长江上游水库改变干流磷通量、效应与修复对策. 湖泊科学, 2018, 30(4): 865-880. DOI:10.18307/2018.0401]
[36]
Guo PC, Wang PF, Jia SB. Research progress of the environmental impact factors of internal phosphorus loading release in river system. Journal of Nanjing Forestry University: Natural Science Edition, 2008, 32(3): 117-121. [郭鹏程, 王沛芳, 贾锁宝. 河流内源磷释放环境影响因子研究进展. 南京林业大学学报:自然科学版, 2008, 32(3): 117-121.]
[37]
Zhang XY, Ma PM, Gao QH et al. Spatial-temporal variations of water quality in upstream and downstream of Three Gorges Dam. J Lake Sci, 2019, 31(3): 633-645. [张馨月, 马沛明, 高千红等. 三峡大坝上下游水质时空变化特征. 湖泊科学, 2019, 31(3): 633-645. DOI:10.18307/2019.0303]
[38]
Du YJ. Analysis of water pollution change and water quality evaluation in the middle reaches of Heihe River. Gansu Water Conservancy and Hydropower Technology, 2009, 45(1): 35-37. [杜有军. 黑河干流中游段水污染变化分析与水质评价. 甘肃水利水电技术, 2009, 45(1): 35-37.]
[39]
Li H. Water quality assessment and trend analysis in the Heihe River Basin[Dissertation]. Lanzhou: Lanzhou Uinverity of Technology, 2014. [李慧.黑河流域水质评价及变化趋势分析[学位论文].兰州: 兰州理工大学, 2014. ]