湖泊科学   2019, Vol. 31 Issue (3): 633-645.  DOI: 10.18307/2019.0303.
0

研究论文

引用本文 [复制中英文]

张馨月, 马沛明, 高千红, 严海涛, 钱宝, 三峡大坝上下游水质时空变化特征. 湖泊科学, 2019, 31(3): 633-645. DOI: 10.18307/2019.0303.
[复制中文]
ZHANG Xinyue, MA Peiming, GAO Qianhong, YAN Haitao, QIAN Bao. Spatial-temporal variations of water quality in upstream and downstream of Three Gorges Dam. Journal of Lake Sciences, 2019, 31(3): 633-645. DOI: 10.18307/2019.0303.
[复制英文]

基金项目

国家重点研发计划项目(2016YFA0600901)和湖北省自然科学基金项目(2017CFB312)联合资助

作者简介

张馨月(1986~), 女, 工程师; E-mail:157130765@qq.com

通信作者

钱宝, E-mail:jacber@163.com

文章历史

2018-06-07 收稿
2018-09-29 收修改稿

码上扫一扫

三峡大坝上下游水质时空变化特征
张馨月1 , 马沛明2 , 高千红1 , 严海涛1 , 钱宝3     
(1: 长江水利委员会长江三峡水文水资源勘测局, 宜昌 443000)
(2: 水利部中国科学院水工程生态研究所水利部水工程生态效应与生态修复重点实验室, 武汉 430079)
(3: 长江水利委员会水文局, 武汉 430010)
摘要:为探索三峡大坝上下游(坝上99.9 km、坝下63.0 km、全长162.9 km)水质时空变化特征,运用主成分分析和方差分析对2016年近坝段水质时空变化特征进行了分析.主成分分析表明,水文因子流量(Q)、气温(T)、水位(Z)和水质因子(水温(WT)、pH、电导率(EC)、溶解氧(DO)、悬浮物(SS)、高锰酸盐指数(CODMn)、硫酸盐(SO42-)、氟化物(F-)、总硬度(T-Hard)、硝态氮(NO3--N)、总氮(TN)和硒(Se))的变化主导着研究区域水质变化;各采样点主成分得分和双因素方差分析结果显示研究区域水质因子时间变化主要呈现出季节和不同水库运行时期的差异.消落期(2-5月),T-Hard、F-、SO42-和EC是影响河流水质变化的主导因子;汛期(7-8月),Q、SS、CODMn、NO3--N、TN和Se是影响河流水质变化的主导因子;T和WT主导着汛末(9月)河流水质变化,并引起了DO等理化特性的变化;高水位运行期(12月),Cl-是影响河流水质变化的主导因子.现阶段,DO、有机污染物(CODMn)、无机盐(SO42-和F-)、营养盐类(NO3--N和TN)、类金属元素(Se)和水体的矿化程度(T-Hard)的变化主导着区域水质的变化,是三峡大坝近坝段水域水质的控制因子.方差分析表明,河流的理化特性(DO、pH和SS)、营养盐组分构成(NH3-N和NO3--N)、无机盐类(EC和Cl-)、石油类有机污染物及粪大肠菌群(FC)等指标在坝上与坝下断面存在显著性差异.气温、水温、降雨、含沙量的季节性影响因素和水库调度运行模式是影响近坝段水质时间差异的主要因子;空间差异主要受城区污染排放和三峡水库调度引起的坝上和坝下水文和水动力学条件差异影响.因此控制研究区域因人类活动等造成的外源性污染,并针对不同类污染物质的季节变化特征实施合理的水库运行方式是近坝段水质提升的关键.
关键词水质    多元分析    时空变化特征    三峡大坝    上下游    
Spatial-temporal variations of water quality in upstream and downstream of Three Gorges Dam
ZHANG Xinyue1 , MA Peiming2 , GAO Qianhong1 , YAN Haitao1 , QIAN Bao3     
(1: Three Gorges Hydrology and Water Resources Survey Bureau, Bureau of Hydrology, Changjiang Water Resources Commission, Yichang 443000, P.R.China)
(2: Key Laboratory of Ecological Impacts of Hydraulic-Projects and Restoration of Aquatic Ecosystem of Ministry of Water Resources, Institute of Hydroecology, Ministry of Water Resources and Chinese Academy of Sciences, Wuhan 430079, P.R.China)
(3: Bureau of Hydrology, Changjiang Water Resources Commission, Wuhan 430010, P.R.China)
Abstract: In order to explore the spatio-temporal variations of water quality in the upstream and downstream (99.9 km upstream and 63.0 km downstream of the dam, total length 162.9 km) of the Three Gorges Dam, principal component analysis and variance analysis were used to analyze the vicinity dam water quality in 2016. Principal component analysis showed that the main factors dominating the water quality changes in the study area were hydrological factors, such as discharge (Q), temperature (T), water level (Z) and water quality factors, such as water temperature (WT), dissolved oxygen (DO), suspended substance (SS), permanganate index (CODMn), sulfate (SO42-), fluoride (F-), total hardness (T-Hard), nitrate nitrogen (NO3--N), total nitrogen (TN) and selenium (Se). The principal component scores of each sampling point and two-way ANOVA showed that the time variation of water quality factors in the study area were mainly different in season and operation period of different reservoirs, and the water quality changed in reservoir drawdown period (February to May) was mainly dominated by T-Hard, F-, SO42- and electronic conductivity (EC), while the main leading factors were Q, SS, CODMn, NO3--N, TN and Se in flood season (July to August), and the water quality in the end of flood season (September) was dominated by T and WT, and caused changes in the physical and chemical properties such as DO. In addition, chloride (Cl-) was the leading factor in the period of high water level (December). At the present stage, changes in DO, organic pollutants (CODMn), inorganic salts (SO42-, F-), nutrient salts (NO3--N and TN), Se and mineralization degree of water (T-hard) dominate the changes in regional water quality and make the greatest effect to the water quality among all parameters. The analysis of variance showed that the physical and chemical characteristics (DO, pH and SS), nutrient composition (NH3-N and NO3--N), inorganic salts (EC and Cl-), petroleum organic pollutants and fecal coliform (FC) of rivers were varied significant in the upstream and downstream of the dam. The seasonal influence factors such as temperature, water temperature, rainfall, sediment concentration and reservoir operation mode were the main factors affecting the time difference of water quality near the dam. The spatial difference were mainly affected by the pollution of urban and the change of hydrological and hydrodynamical conditions by Three Gorges Reservoir operation. Therefore, controlling the exogenous pollution caused by human activities in the study area and implementing a reasonable reservoir operation mode for the seasonal variation characteristics of types of pollutants are the key to water quality improvement in the near dam section.
Keywords: Water quality    multivariate analysis    spatio-temporal distribution characteristics    Three Gorges Dam    upstream and downstream    

河流的水质主要受到地质、水文、气象等自然因素周期性波动以及受污染排放、兴水利设施等人类活动而导致的随机波动两方面的影响[1].随着社会经济的发展,人类活动加剧,长江水质状况面临巨大压力.习近平2018年4月在深入推动长江经济带发展座谈会上再一次强调坚持共抓长江大保护、不搞大开发, 湖北省也提出了九大行动助推长江大保护, 长江水质状况日趋受到国家和地方政府的重视.三峡大坝蓄水以来,形成了三峡库区河道型水库与坝下天然河道两种形态,其水质因子也呈现出不同的变化特征, 三峡工程为防洪、发电、航运等带来经济效益的同时,带来的水环境问题也引起了社会各界学者的关注[2].近年来关于三峡水质变化的研究大多集中在库区支流[3-4],对三峡大坝近坝段长江干流的水质缺乏大尺度和多断面时空变化特征研究.

水质监测数据具有长期连续、多站点、多指标的特点,指标之间的信息通常存在相关性和重叠现象[5],如何浓缩这些信息,关系到准确和高效地对水质进行分析与评价.多元统计方法中主成分分析(PCA)和方差分析(ANOVA)被广泛应用于水质的时空特性分析,如运用PCA对湖泊[6]、湿地[7]、水库[8]、海洋[9]、河流[10-11]水质主导因子进行提取,定性识别污染源,运用ANOVA分析水质因子的时空差异性[12-13].这些研究中大多利用主成分提取水质变化的主导因子, 进而通过主成分得分评价水质好坏, 或单一地对因子载荷进行分析评判, 而将环境因子载荷系数与样点得分结合起来进行时空特性分析的并不多.本研究对2016年三峡大坝上下游162.9 km(培石-虎牙滩断面)的水环境因子特征进行评价和相关性分析,而后进一步结合PCA和ANOVA,识别主导研究河段水质变化的主要因子,判别水质时空差异性,并对形成时空差异的原因进行探讨,旨在为长江干流水质监测数据提供合理的分析评价,为三峡水库调度、长江持续健康发展提供决策参考.

1 研究区域与研究方法 1.1 研究区域概况

研究区域河段位于湖北省西南部,西起长江上游培石(重庆市与湖北省交界的鳊鱼溪口下300 m),东至长江中游虎牙滩(宜昌市猇亭区虎牙滩水文断面),坝上99.9 km,坝下63.0 km.该区域处于中纬度亚热带地区,属亚热带季风气候,雨热同季,四季分明.

1.2 研究方法 1.2.1 样品采集与分析

本研究共布置了8个采样断面(图 1),三峡库区坝上断面分别为培石、官渡口、巴东、庙河,三峡大坝下游断面分别为黄陵庙、南津关、宜昌、虎牙滩.每个断面分别在左岸、中泓、右岸布设3条采样垂线,2016年每月上旬采样1次,采样方法参考《水环境监测规范》(SL219-2013).样点名称用断面字母缩写加月份表示,培石断面1-12月即为P1~P12,以此类推,官渡口为G1~G12, 巴东为B1~B12, 庙河为M1~M12, 黄陵庙为HLM1~HLM12, 南津关为N1~N12, 宜昌为Y1~Y12, 虎牙滩为HYT1~HYT12.

图 1 三峡大坝上下游水质采样断面 Fig.1 Distribution of water sampling sections in upstream and downstream of Three Gorges Dam

监测指标包括气温(T)、水位(Z)、流量(Q)、水温(WT)、pH值、电导率(EC)、溶解氧(DO)、悬浮物(SS)、高锰酸盐指数(CODMn)、五日生化需氧量(BOD5)、氯化物(Cl-)、硫酸盐(SO42-)、氟化物(F-)、总硬度(T-Hard)、氨氮(NH3-N)、硝态氮(NO3--N)、总磷(TP)、总氮(TN)、硒(Se)、砷(As)、石油类有机污染物(POP)、粪大肠菌群(fecal coliforms,FC),共计22项.其中气温、水温、pH值、电导率、溶解氧采用便携式水质分析仪HQ40d(美国)现场测定,悬浮物采用重量法(GB/T 11901-1989)测定,水位、流量摘自水文三峡局水情报汛数据,其余指标测定参照国家《地表水环境质量标准》(GB 3838-2002).

1.2.2 基于多元统计方法的水质特征数据分析

本文首先运用主成分分析方法[14]提取引起该流域水质变化的主导因素并加以分析;再利用重复测量方差分析[14]判别三峡大坝上游与下游不同空间类别上有显著性差异的水质指标,再在不同空间类别上采用双因素方差分析量化季节变化与水库调度对水质的影响机制,识别季节变化和水库不同运行时期影响下水质状况的时间差异;并对主导研究区域水质状况的水质因子时空差异的潜在因素进行探索分析.对于低于检出限的参数,在统计分析时按其方法检出限的半值处理.主成分分析前,首先对原始数据进行KMO和Bartlett球度检验,一般认为KMO>0.7且P < 0.05时适合做主成分分析,为了消除各指标间不同量纲的影响,需将原始数据进行Z-Score标准化处理[15](标准化数据=(原数据-均值)/标准差;溶解氧数据进行线性变换[16], a=12,b=-1,然后再标准化).对空间差异显著的水质指标,分为坝上和坝下两组分别求年均值,并对平均值做箱式图进行空间差异直观分析比较.双因素方差分析季节划分为春(3-5月)、夏(6-8月)、秋(9-11月)、冬(12-2月),分别用A、B、C、D表示;三峡水库于每年9月上旬开始蓄水,11月初蓄至175 m, 11月至次年4月,根据生态、航运等需要逐步消落至155 m,5月份加大出力运行,6月上旬消落至汛限水位145 m,根据水库周期蓄水方式和2016年水库调度水位运行情况(图 2),三峡水库运行可以明显划分为消落期(1-4月)、集中消落期(5-6月)、汛期(7-9月)、蓄水期(10月)和蓄水后高水位运行期(11-12月)5个时期,分别用T1、T2、T3、T4和T5表示,由于2016年9月10日开始蓄水,而水质取样在开始蓄水之前,故把9月归为汛期.上述所有分析均在IBM SPSS Statistics 19.0和Origin 8.0软件下进行.

图 2 2016年三峡水库坝前水位变化及调度节点 Fig.2 Water-level fluctuation and scheduling node behind the dam of the Three Gorges Reservoir in 2016
2 结果与分析 2.1 水环境因子特征与评价

2016年研究区域水质因子均值和实测范围见表 1,研究区域水体呈弱碱性,根据国家《地表水环境质量标准》(GB 3838-2002)基本项目(FC除外)评价标准,研究区域水体全年均达到Ⅲ类水标准,其中汛期(5-10月)为Ⅲ类水,非汛期(1-4月,11-12月)为Ⅱ类水,整体水质状况较好.从实测范围来看,巴东、庙河、黄陵庙、虎牙滩断面POP出现超Ⅲ类水标准限值的情况,研究区域全年FC均超Ⅲ类水标准限值(图 3).水环境因子Pearson相关系数矩阵见附表 1,其中,Q与除F-、T-Hard外的其他因子呈显著相关,与CODMn的正相关性最强,与DO、SO42-、EC呈负相关;DO、SO42-与WT表现出极强的显著负相关;EC与SO42-、T-Hard呈显著正相关;CODMn与SS呈极强显著正相关;F-与T-Hard呈显著正相关;NO3--N与TN、Se呈极显著正相关,TN与Se呈极显著正相关;Se与As呈极显著正相关.

图 3 石油类有机污染物(POP)和粪大肠菌群(FC)的实测范围 Fig.3 Actual ranges of petroleum organic pollutants (POP) and fecal coliforms (FC)
表 1 研究区域水质因子特征 Tab. 1 Characteristics of water quality factors in the studied area
2.2 水质主成分分析

经检验得出, KMO值为0.795且Bartlett检验的P值小于0.0001,说明原始数据之间有较强的相关关系,可以进行主成分分析.本文按照特征值λ>1的原则,提取了引起水质变化的5个主成分(表 2),其累积方差贡献率为81.511 %,能够充分反映原始数据信息.

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

表 2可知,第1主成分PC1的方差贡献率最大,为33.182 %,与之关联的主要参数为Q、T、WT、DO、SS、CODMn和SO42-,反映了水文水动力(Q)和气候条件(T、WT)的变化对水质的影响,并引起了河流的理化特性(DO、SS)变化;人类活动影响下有机污染物(CODMn)和无机盐(SO42)变化主导着研究区域水质的变化. CODMn对PC1有最强的正相关性,变化特征明显.第2主成分PC2的方差贡献率为26.859 %,与之相关联的主要参数为F-、T-Hard、NO3--N、TN和Se,反映了水体营养盐(NO3--N、TN)、无机污染物(F-)和类金属元素(Se)的变化主导着区域水质变化以及水体矿化程度的变化(T-Hard);其中Se对PC2有最强的正相关性,表明其随季节、水文等因素变化明显.第3主成分贡献率为9.725 %,与之关联的参数为Z,反映了水位涨落对水质的影响.第4、5主成分没有与之相关联大于0.7的载荷因子.

为了更清晰地反映研究区域的水质状况,将各采样点得分和环境因子的载荷系数代替原始数据对水质进行分析(图 4),反映了3个方面的信息:(1)样点在排序空间呈明显的季节分布,逆时针方向来看,其中1-5月基本聚集在第3象限,6-8月在第4象限,9和10月在第1象限,11和12月在第2象限;(2)环境因子的箭头越长表示其因子载荷越大,对排序的贡献也越大,箭头间的夹角表示环境因子间的相关性,夹角越小表示相关性越高;图中箭头较长的环境因子在排序空间中分成WT、DO、T,SS、CODMn、Q,NO3--N、TN、Se、As、TP和T-Hard、F-、SO42-、EC四簇排列,主导着研究区域水体水质的变化,且每簇环境因子间夹角小,相关性高. Z、POP、pH、和BOD5的箭头短,其对排序的贡献小,样品间的差异小;(3)将样点向环境因子箭头及其延长线做投影,投影位置在箭头指向方向越远则数值越大,表明该环境因子对其正贡献越大,反之则越小.逆时针方向看,第3象限中,箭头指向了排序空间中T-Hard、F-、SO42-、EC变异最大的方向,可看到3-5月T-Hard、F-、SO42-、EC数值高,主导着水质变化,7-10月相对较低;7月Q、SS、CODMn、NO3--N、TN、Se达到最高值,主导着7-8月各样点的水质变化;9月WT最高,2、3月份较低,而DO浓度相反(因DO浓度越低反映水质越差,分析前对DO浓度做了线性转换,故在图中显示其与水温呈正相关),WT、DO主导着9月各样点水质变化;12月Cl-浓度全年最高,主导各样点水质变化.

图 4 研究区域水文水环境因子PCA分析双序(载荷系数和得分值)图(圆圈表示样点,箭头表示环境因子) Fig.4 The PCA biplots (scores & loadings) of hydrology and water quality factors in the research area
2.3 水质因子时空变化特征 2.3.1 水质因子空间差异特征分析

本文将坝上和坝下断面分成2组对其1-12月数据进行重复测量方差分析,结果表明,NH3-N、NO3--N、EC、SS、DO、Cl-、FC、POP、pH值等水质因子在坝上与坝下断面存在显著性差异(P < 0.05).

图 5可以看出,坝下断面NH3-N、NO3--N、DO、FC、POP、pH值的平均水平要高于坝上断面,坝上断面EC、SS、Cl-的平均水平显著高于坝下断面.其中,NH3-N、DO数值在坝上断面较坝下断面均匀集中,而EC、Cl-等离子水平的数值在坝上断面比坝下断面要分散得多.

图 5 水质因子的空间差异 Fig.5 Spatial variations of water quality factors
2.3.2 水质因子时间差异特征分析

以季节和三峡水库不同运行时期为因子,水质因子为因变量,对坝上和坝下断面分别进行双因素方差分析.结果表明(表 3),研究区域除POP以及坝下断面BOD5以外,其他水质因子均随季节或水库运行时期差异显著(P < 0.01),且季节与水库运行交互作用对水质因子的影响在坝上和坝下断面基本一致(EC除外).

表 3 双因素方差分析显著性检验结果 Tab. 3 Two-way ANOVA significance test results

根据主成分分析的结果,在时间差异显著的指标中挑选DO、T-Hard、CODMn、SO42-、Se、TN和NO3--N 7个代表性指标进一步做季节、水库运行时期及季节和水库运行时期组合平均值间的多重比较,除NO3--N外,其他指标坝上、坝下断面的多重比较结果基本一致,故仅选取坝上断面的多重比较结果为代表进行分析, 结果见表 4.

表 4 季节、水库运行时期及各水平组合的平均值多重比较* Tab. 4 Multiple comparisons of seasons, reservoir operation periods and average combinations of different levels

表 4可知,春、冬季(A、D)和消落期(T1)对应的DO浓度最高,交互作用使得DO浓度最高值出现在消落期和春季的组合(A:T1),最低值出现在夏季和集中消落期组合(B:T2);春季(A)和整个消落期(T1、T2)T-Hard显著高于其他时期,春季和消落期组合(A:T1)T-Hard最高,与春季和集中消落期组合(A:T2)差异不显著,但显著高于其他水平组合,最低值出现在冬季高水位运行期组合(D:T5);CODMn主要受季节和水库不同运行时期的主效应影响,两者的交互作用影响并不显著,最高值出现在夏季和汛期(B、T3),最低值出现在春、冬季节;春季和整个消落期SO42-浓度显著高于其他时期,最高值出现在春季和消落期组合(A:T1),最低值出现在夏季汛期组合(B:T3);Se浓度受季节和水库运行交互作用影响显著,最高值出现在夏季汛期组合(B:T3), 最低值出现在冬季高水位运行期(D:T5)和秋季的蓄水期(C:T4);夏季TN浓度显著高于其他季节,汛期(T3)与整个消落期(T2、T3)差异不显著,但显著高于蓄水期和高水位运行期,交互作用使得夏季和汛期组合(B:T3)TN浓度最高,最低值出现在秋季汛期组合,即9月份;坝上、坝下断面NO3--N浓度最高值均出现在夏季,坝下断面集中消落期(T2)显著高于其他时期,受季节与水库运行的交互影响,坝上断面NO3--N最高值出现在夏季汛期组合(B:T3),但差异与春、夏季的消落期并不显著,而坝下断面最大值出现在夏季与集中消落期组合(B:T2),与夏季汛期(B:T3)无显著差异.与TN浓度一样,坝上和坝下断面NO3--N浓度最低值均出现在9月份.

3 讨论 3.1 大坝上下游水质空间变化特征

三峡大坝蓄水以来近坝段上下游水质空间差异主要表现为河流的理化特性(DO、pH和SS)、营养盐组分构成(NH3-N和NO3--N)、无机盐类(EC和Cl-)、石油类有机污染物及微生物(FC)等指标的差异.大坝建设后改变了河流的水文、水动力学条件,对水质也产生了一定影响,主要表现为水库富营养化、藻类繁殖加剧、水库盐度增高、水库水温分层的次生影响[16].据报道,韩国锦江建坝后加速了水体的富营养化[17],南非的Bospoot大坝,水库坝前水域出现了明显的水温分层、盐度和电导率随时间升高等问题[18].本研究表明,坝上断面的无机盐(EC、Cl-)浓度显著高于坝下断面,这说明坝上断面的盐度较坝下游天然河道出现升高的现象.研究数据显示[19],三峡水库蓄水后,入库沙量远大于出库沙量,排沙比不到50 %,而泥沙对Cl-等离子具有较强的吸附作用,从而导致了坝上和坝下断面SS、EC、Cl-浓度的差异.其中坝上断面EC、Cl-浓度较坝下断面波动大,考虑由于坝上断面受上游来水、三峡水库调度等水文和水动力学条件变化的直接影响,且区间较长,泥沙存在区间变化,导致其EC、Cl-浓度变化要大于坝下断面,经三峡水库调节后,这些因素对坝下断面的影响要小得多.双因素方差分析结果也进一步说明坝上断面EC受季节和水库调度交互影响显著而坝下断面并未受到交互影响.坝下断面NH3-N、NO3--N、POP、FC显著高于坝上断面.有研究认为[20],长江干流主要污染物来源为城镇生活污水排放和农业污染排放.坝下断面主要位于城市生活主城区,且在宜昌和虎牙滩断面之间有一大型以NH3-N为主要污染的混合废污水入河排污口,考虑受城区生活等混合废污水影响,坝下游NH3-N、NO3--N等营养盐浓度高于坝上断面.值得一提的是,虽然水库中底泥可与磷结合,对河流TP具有较高的拦截率[21],但在此次调查结果坝上断面与坝下断面TP、TN浓度并未出现显著差异,推测与坝上SS浓度高,颗粒态N、P浓度高有关,而坝下SS浓度相对较低,水中颗粒态N、P浓度低,但坝下断面受城市污染较坝上断面高,溶解态N、P浓度高,故TP、TN未出现显著差异,而表现为营养盐组份构成的差异.周建军等[22]最新研究指出水库大量拦沙拦磷,水体纳污的缓冲机制被破坏,使得坝下游溶解磷浓度增加.另外,双因素方差分析表明,坝上断面NH3-N和NO3--N浓度仅受季节影响极显著,而坝下断面受季节和水库不同运行时期双因素影响显著(P < 0.01),因此坝上、坝下营养盐空间差异还受水库调度的影响;坝下断面FC丰度显著高于坝上断面,这与江磊等[23]研究的粪大肠菌群丰度受人类干扰影响显著的结论相符.三峡水库调度,不同水位运行下对水质产生了一定的影响[24],但近坝段水质状况的空间差异并不如时间差异显著.

3.2 大坝上下游水质时间变化特征

主成分与双因素方差分析均表明研究区域大多数水质因子时间变化呈现出季节和不同水库运行时期的变化特征,但由于三峡水库以季节周期性的调度方式为主,故其季节变化对水质因子的影响更大.其中理化指标(DO、T-Hard)、无机盐(Cl-)、水体营养盐(NO3--N、TN)、类金属(Se、As)等受季节与水库运行交互作用的影响极显著(P < 0.01),表明水温、降雨、含沙量等季节因素与水库调度引起的水文水动力因素之间的交互作用产生的新效应对这几类水质指标的影响是极其显著的.主成分分析表明消落期(2-5月)研究区域水体无机盐类浓度(F-和SO42-)和T-Hard较高,是影响水质变化的主导因子;流量、泥沙、有机污染物(Q、SS和CODMn)、水体营养盐(NO3--N和TN)和类金属(Se)主导汛期(7-8月)水质变化;气温、水温等气象因子是汛末(9月)水质变化的主导因子,并引起了DO等理化特性的变化;高水位运行期(12月)Cl-浓度较高,是影响水质变化的主导因子.坝上、坝下断面POP以及坝下断面的BOD5季节和水库不同运行时期均差异不显著,BOD5反映水体有机耗氧类污染,其浓度受水体滞留时间[25]和水体浑浊度[26]等多方面影响,坝上断面BOD5浓度受水库调度与季节变化直接影响,而坝下断面为天然河道的水流形态,单一的季节或水库运行时期难以识别其变化规律;POP未表现出明显的时间差异,可能是因为研究区域位于长江上游与中游的航运黄金水道,来往船只较多,石油类污染存在较大的随机性.

多重比较分析结果显示,春季消落期理化指标DO、T-Hard均显著高于其他时期,冬季高水位运行期T-Hard显著低于其他时期,与主成分分析结果一致. T-Hard与SS有极强的负相关性(P < 0.01),冬季高水位运行期,大部分粗颗粒泥沙沉入河底,可能由于蓄水刚完成,水中具有较强吸附能力的极细颗粒泥沙还未完全沉入水底,钙镁离子被大量吸附,从而使得水体硬度降低,而经过消落期长江水体的进一步澄清,长江自然背景的硬度达到最大;有研究发现,SO42-浓度受降雨和蓄、弃水过程等因素影响[27],本次监测发现汛期SO42-浓度最低,并随着水库蓄水逐渐升高,在消落期的春季达到最高,随后急剧降低. SO42-浓度受季节影响极为显著,春季雨水较少,水体离子浓度较高,同时在高水位运行期之后水中细颗粒泥沙逐渐减少,水中被泥沙吸附的SO42-减少,使得峰值出现在春季,随着夏季降雨量增加的稀释作用和泥沙的吸附作用,SO42-浓度迅速降低,这也从侧面反映了研究区域SO42-并非主要来自陆源污染,另外由于夏季泥沙含量大,采用离子色谱仪分析过滤后的样品,也可能是导致SO42-夏季浓度低的原因之一;CODMn常被作为反映受工农业、生活有机污染物和还原性无机物质污染程度的综合指标[28],主成分分析双序图与Pearson相关性分析结果一致,反映CODMn与SS、Q之间具有极强的相关性(P < 0.01),表明水体SS浓度对CODMn有直接的影响,而流量等水力学条件变化对SS和CODMn又有着直接或间接的影响.有研究表明[29],CODMn与泥沙有着显著的相关性,在水体未受污染的情况下,泥沙及其中有机质浓度决定了CODMn浓度大小,因此夏季长江水体较大的含沙量是导致汛期CODMn浓度高于其他时期的关键因素;TN、NO3--N营养盐受季节和水库调度交互作用的影响极显著,其最大值和最小值分别出现在夏季汛期和汛期初秋(9月),由于汛期来水来沙量大,水库调度汛期阶段一方面使汛期TN浓度升高,均衡了与消落期、集中消落期的差异性,另一方面由于9月蓄水前期含沙量迅速降低,TN浓度随之迅速下降,随着蓄水进程的推进,库区高水位运行,库区营养盐物质随泥沙沉降,清水下泄,水库运行使秋、冬季TN浓度降低.而坝下断面集中消落期的NO3--N浓度显著高于其他时期,这也从侧面反映了水库的集中消落,下泄流量增加导致坝下NO3--N增加.另外由宜昌水文站2016年降雨量观测数据发现,出现TN、NO3--N浓度低值的9月份降雨量全年最低,而最高值出现的集中消落期和汛期其降雨量也是全年最高的,因此除了水库调度的影响外,因降雨径流的淋溶、冲刷作用带入水体的陆源性营养盐污染也不容忽视. Se等类金属在夏季汛期显著高于其他时期,而最低值出现在蓄水期和高水位运行期,可能是由于夏季降雨较多,使得陆地上Se、As等类金属污染随降雨进入河流,再加之夏季流量大,水体紊动强烈,河流底质扰动,底质中的Se、As重新进入水中.

纵向来看,陈静生等[30]的研究表明1960s-1980s长江中上游水质呈现酸化趋势,水质变化趋势具体表现为碱度降低,SO42-、Ca2+、Cl-和氮化合物浓度升高,1962-1966年宜昌断面Cl-和SO42-浓度年均值分别为6.86和16.03 mg/L[31], 本次2016年监测数据分别为13.23和40.36 mg/L,呈现明显的上升趋势,长江中上游TN浓度较1960s-1980s也表现了较明显的上升趋势,总碱度表现为降低的趋势.长江干流1960s与2000s(2000-2007年)[32]离子总量全年均值均在180~230 mg/L之间,这反映了随着时间的推移,长江流域的地质条件并未发生大的变化,长江水体未受到严重污染.戴润泉等[33]研究指出三峡水库蓄水前(1996-2001年)库区江段主要污染因子为CODMn、总铅、NH3-N,本次主成分分析显示CODMn仍是主导近坝段水质的变化的主要因子,而NH3-N不再是主要污染因子,这也从侧面反映出由于国家保护长江、污染治理力度加大,工业废污水的入江排放得到了有效控制.三峡水库蓄水初期(2003-2004年),近坝水域DO、CODMn、BOD5等指标在坝后略有改善,Cu、As等重金属和TP、TN等营养盐浓度在坝上、坝下断面均降低[34], 但TP、TN浓度偏高.随着蓄水进程的推进,最高蓄水位从135 m升高到175 m,水质状况变化趋势也发生了改变,如坝上断面(培石)DO和NH3-N浓度呈下降趋势[35],本次监测数据较蓄水后至2008年稳中有升,TN浓度在库首、库中呈现明显的上升趋势[36].此次分析表明现阶段(2016年),DO、有机污染物(CODMn)、无机盐(SO42-、F-)、营养盐类(NO3--N、TN)、类金属元素(Se)和水体的矿化程度(T-hard)主导着研究区域的水质变化,是三峡大坝近坝段水域水质的控制因子.

4 结论

1) 研究区域水体全年均达到《地表水水环境质量标准》(GB 3838-2002)Ⅲ类水标准,其中汛期(5-10月)为Ⅲ类水,非汛期(1-4月,11-12月)为Ⅱ类水,整体水质状况较好.主成分分析表明研究区域水质主导因素为水文因子Q、T、Z和水质因子WT、DO、SS、CODMn、SO42-、F-、T-Hard、NO3--N、TN、Se.

2) 水质因子空间差异主要表现为河流的理化特性(DO、pH和SS)、营养盐组分构成(NH3-N和NO3--N)、无机盐类(EC和Cl-)、石油类有机污染物及微生物(FC)等指标的差异.

3) 水质因子时间差异主要呈现出明显的季节和不同水库运行时期的差异.现阶段(2016年)DO、有机污染物(CODMn)、无机盐(SO42-和F-)、营养盐类(NO3--N和TN)、类金属元素(Se)和水体的矿化程度(T-Hard)的变化主导着区域水质的变化,是三峡大坝近坝段水域水质的控制因子,且绝大多数控制因子(有机污染物、无机盐除外)受季节与水库运行调度交互影响极显著,可见,水库调度与季节的交互效应对近坝段水质有着不可忽视的作用.

4) 水质因子空间差异主要受城区污染排放和三峡水库调度引起的坝上和坝下断面之间的水文和水动力学条件差异的影响;时间差异主要受气温、水温、降雨、含沙量的季节性影响因素和水库调度运行模式影响.因此控制研究区域因人类活动等造成的外源性污染,并针对不同类污染物质的季节变化特征实施合理的水库运行方式是近坝段水质提升的关键.

5 附录

附录Ⅰ见电子版(DOI:10.18307/2019.0303).

附录Ⅰ 水环境因子Pearson相关系数矩阵(n=96) Appendix Ⅰ Pearson correlation coefficients matrix of water environmental factors

致谢: 感谢长江三峡水文水资源勘测局教授级高级工程师李云中对写作的帮助!

6 参考文献

[1]
Hong YP. Impact of the construction of the Yangtze Three Gorges Water Control Project on water quality. Water Resources Protection, 1986(1): 29-34. [洪一平. 长江三峡水利枢纽兴建对水质的影响. 水资源保护, 1986(1): 29-34.]
[2]
Wu JG, Huang JH, Han XG et al. Three-Gorges Dam-experiment in habitat fragmentation. Science, 2003, 300(5623): 1239-1240. DOI:10.1126/science.1083312
[3]
Ye L, Li DF, Tang T et al. Spatial distribution of water quality in Xiangxi River, China. Chinese Journal of Applied Ecology, 2003, 14(11): 1959-1962. [叶麟, 黎道丰, 唐涛等. 香溪河水质空间分布特性研究. 应用生态学报, 2003, 14(11): 1959-1962. DOI:10.3321/j.issn:1001-9332.2003.11.034]
[4]
Ran GH, Ge JW, Miao WJ et al. Spatial-temporal differentiation of water quality in Gufu River of Three Gorges Reservoir. Acta Ecologica Sinica, 2013, 33(17): 5385-5396. [冉桂花, 葛继稳, 苗文杰等. 三峡库区古夫河水质时空分异特征. 生态学报, 2013, 33(17): 5385-5396.]
[5]
Wang LJ, Wang X, Liu LS et al. Multivariate analysis of water factors in Dongting Lake. Research of Environmental Sciences, 2013, 26(1): 1-7. [王丽婧, 汪星, 刘录三等. 洞庭湖水质因子的多元分析. 环境科学研究, 2013, 26(1): 1-7. DOI:10.13198/j.res.2013.01.5.wanglj.005]
[6]
Yu Y, Zhang M, Qian SQ et al. Current status and development of water quality of lakes in Yunnan-Guizhou plateau. J Lake Sci, 2010, 22(6): 820-828. [于洋, 张民, 钱善勤等. 云贵高原湖泊水质现状及演变. 湖泊科学, 2010, 22(6): 820-828. DOI:10.18307/2010.0603]
[7]
Li J, Qi PS, Ma Y et al. Phytoplankton community characteristics of the icebound season and its relationship to the environmental variables in the Zhalong Wetland, China. African Journal of Biotechnology, 2011, 10(38): 7416-7426.
[8]
Shi LD, Zhu WJ, Hu JX et al. Characteristics and influencing factors analysis of spatio-temporal variations of water quality in Tanpu Reservoir and its tributaries. Journal of Hydroeclology, 2013, 34(5): 9-15. [施练东, 竺维佳, 胡菊香等. 汤浦水库及入库支流水质时空变化特征与影响因素分析. 生态学杂志, 2013, 34(5): 9-15. DOI:10.3969/j.issn.1674-3075.2013.05.002]
[9]
Wu ML, Wang YS, Lin L et al. Study on changes characterization of water quality by principal component analysis in Daya Bay. Marine Environmental Science, 2009, 28(3): 279-282. [吴梅林, 王友紹, 林立等. 基于主成分分析研究大亚湾水质时空变化特征. 海洋环境科学, 2009, 28(3): 279-282. DOI:10.3969/j.issn.1007-6336.2009.03.012]
[10]
Boyacioglu H. Surface water quality assessment using factor analysis. Water SA, 2006, 32(3): 389-393.
[11]
Liu DL, Liu XZ. Application of principal component analysis to the comprehensive evaluation of water quality in river. Reserch of Soil and Water Conservaion, 2006, 13(3): 124-128. [刘德林, 刘贤赵. 主成分分析在河流水质综合评价中的应用. 水土保持研究, 2006, 13(3): 124-128. DOI:10.3969/j.issn.1005-3409.2006.03.040]
[12]
Xu HS, Zhao L, Sun HS et al. Water quality analysis of Beijing segment of south-to-north water diversion middle route project. Environmental Science, 2017, 38(4): 1357-1364. [徐华山, 赵磊, 孙昊苏等. 南水北调中线北京段水质状况分析. 环境科学, 2017, 38(4): 1357-1364.]
[13]
Chen HY, Liu MS, Yuan F et al. Analysis of relationship between microhabitat characteristics and major water quality indicators in littoral zone of Gonghu Bay, Tai Lake. Wetland Science, 2016, 19(4): 580-585. [陈怀艳, 刘茂松, 袁芳等. 太湖贡湖湾滨岸带微生境特征与主要水质指标的关系分析. 湿地科学, 2016, 19(4): 580-585.]
[14]
Li YG, Du HW, Huang YS et al. SPSS19.0 introduction and improvement of statistical analysis. Beijing: Tsinghua University Press, 2014: 308. [李玉光, 杜宏巍, 黄永生等. SPSS19.0统计分析入门与提高. 北京: 清华大学出版社, 2014: 308.]
[15]
Qi JY, Ruan XH. Dam construction-induced environmental impact on riverine ecosystem. Journal of Hohai University: Natural Sciences, 2005, 33(1): 37-40. [祁继英, 阮晓红. 大坝对河流生态系统的环境影响分析. 河海大学学报(自然科学版), 2005, 33(1): 37-40. DOI:10.3321/j.issn:1000-1980.2005.01.009]
[16]
Yao HM, Huang RT, Gan FX et al. Principal component analysis of the water quality evaluation in East Lake. Geomatics and Information Science of Wuhan University, 2005, 30(82): 732-735. [姚焕玫, 黄仁涛, 甘复兴等. 用改进的主成分分析法对东湖的水质污染进行评价. 武汉大学学报, 2005, 30(82): 732-735.]
[17]
Jeong YH, Yang JS, Park K. Changes in water quality after the construction of an estuary dam in the Geum River estuary dam system, Korea. J Coastal Res, 2014, 298: 1278-1286. DOI:10.2112/JCOASTRES-D-13-00081.1
[18]
Mogakabe DE, Ginkel CV ed. The water quality of Bospoort Dam. Department of water affairs and forestry, South Africa, 2008.
[19]
Lu JY, Huang Y, Wang J. The analysis on reservoir sediment deposition and downstream river channel scouring after impoundment and operation of TGP. Engineering Science, 2011, 13(7): 129-136. [卢金友, 黄悦, 王军. 三峡工程蓄水运用后水库泥沙淤积及坝下游河道冲刷分析. 中国工程科学, 2011, 13(7): 129-136. DOI:10.3969/j.issn.1009-1742.2011.07.021]
[20]
Ding XY, Zhou HD, Lei XH et al. Hydrological and associated pollution load simulation and estimation for the Three Gorges Reservoir of China. Stoch Environment Resource & Risk Assessment, 2013, 27(3): 617-628.
[21]
Maavara T, Parsons CT, Ridenour C et al. Global phosphorus retention by river damming. Pro Natl Acad Sci USA, 2015, 112: 15603-15608.
[22]
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]
[23]
Jiang L, Zhu DJ, Chen YC et al. Analysis on fecal coliform pollution in surface water of China. Advances in Science and Technology of Water Resources, 2015, 35(3): 11-18. [江磊, 朱德军, 陈永灿等. 我国地表水体粪大肠菌群污染现状分析. 水利水电科技进展, 2015, 35(3): 11-18.]
[24]
Huang QC, Liu GL, Wang CY et al. Water quality variation characteristics of the main stream of the Three Gorges Reservoir under different water levels operation. Yangtze River, 2015, 46(s1): 132-136. [黄庆超, 刘广龙, 王雨春等. 不同水位运行下三峡库区干流水质变化特征. 人民长江, 2015, 46(s1): 132-136.]
[25]
Zhang S, Li CM, Lv PY et al. Spatial variation of BOD5 and CODMn in the Three Gorges Reservoir. J Lake Sci, 2007, 19(1): 70-76. [张晟, 李崇明, 吕平毓等. 三峡水库成库后水体中CODMn、BOD5空间变化. 湖泊科学, 2007, 19(1): 70-76. DOI:10.18307/2007.0111三峡水库成库后水体中CO]
[26]
Zhang J, Yang ZY, Yang XX. A study on the relation between the turbidity degree of sewage and the concentration of CODCr, BOD5 and coliform bacteria. Modern Preventive Medicine, 2004, 31(4): 505-512. [张静, 杨正炎, 杨湘霞. 生活污水浑浊度与CODCr、BOD5及粪大肠菌群浓度的关系. 现代预防医学, 2004, 31(4): 505-512. DOI:10.3969/j.issn.1003-8507.2004.04.014]
[27]
Xing CX. Analysis of sulfate change for Huangbizhuang Reservoir and Gangnan Reservoir. South-to-North Water Transfers and Water Science & Technology, 2013, 11(1): 66-67. [邢晨曦. 黄壁庄和岗南水库硫酸盐变化规律分析. 南水北调与水利科技, 2013, 11(1): 66-67.]
[28]
Editorial board of "water and wastewater monitoring and analysis method", Ministry of Environmental Protection of the People's Republic of China. Monitoring and analysis methods of water and wastewater: fourth edition. Beijing: China Environmental Science Press, 2002. [国家环境保护总局《水和废水监测分析方法》编委会. 水和废水监测分析方法:第4版. 北京: 中国环境科学出版社, 2002.]
[29]
Zhang JS, Zhang Y, Yu T et al. Influences of the suspended matter on the water quality parameters including COD, Potassium Permanganate Index and BOD5 in the Yellow River, China. Acta Scientiae Circumstantiae, 2004, 24(3): 70-76. [张静生, 张宇, 于涛等. 泥沙对黄河水质参数COD、高锰酸钾指数和BOD5的影响——三论黄河的COD与高锰酸盐指数不能真实反映其污染状况. 环境科学学报, 2004, 24(3): 70-76.]
[30]
Chen JS, Xia XH, Zhang LT et al. Relationship between water quality changes in the Yangtze, Yellow and Songhua Rivers and the economic development in the river basins. Acta Scientiae Circumstantiae, 1999, 19(5): 500-505. [陈静生, 夏星辉, 张利田等. 长江、黄河、松花江60-80年代水质变化趋势与社会经济发展的关系. 环境科学学报, 1999, 19(5): 500-505. DOI:10.3321/j.issn:0253-2468.1999.05.007]
[31]
Chen JS, Guan WR, Xia XH et al. A probe into several problems of water-quality trends in the mainstream of Yangtze river from 1960's to1980s. Environmental Chemistry, 1998, 17(1): 8-13. [陈静生, 关文荣, 夏星辉等. 长江干流近三十年来水质变化探析. 环境化学, 1998, 17(1): 8-13.]
[32]
Wang JS, Bian JJ, Chen XG. Analysis on changing tendency of water quality in trunk stream of Yangtze River since 2000. Hubei Water Power, 2009(2): 1-3. [汪金成, 卞俊杰, 陈新国. 2000年以来长江干流水质变化趋势分析. 湖北水利发电, 2009(2): 1-3.]
[33]
Dai RQ, Zang XP, Qiu GS. Water quality in the Three Gorges Reservoir before impoundment. Resources and Environment in the Yangtze Basin, 2004, 13(2): 124-127. [戴润泉, 臧小平, 邱光胜. 三峡水库蓄水前库区水质状况研究. 长江流域资源与环境, 2004, 13(2): 124-127. DOI:10.3969/j.issn.1004-8227.2004.02.005]
[34]
Chen YC, Fu J, Liu ZW et al. Evaluation of water quality near the dam before and after the storage of Three Gorges Reservoir. Journal of Hydroelectric Engineering, 2007, 26(4): 51-55. [陈永灿, 付健, 刘昭伟. 三峡水库蓄水前后近坝水域的水质评价与分析. 水利发电学报, 2007, 26(4): 51-55.]
[35]
Liu XB, Liu XX, Fu DL. Analysis on water quality variation of the mainstream and tributaries in Wushan section, Three Gorges Reservoir region. Environmental Science and Management, 2010, 35(5): 124-125. [刘学斌, 刘晓霞, 傅道林. 三峡工程库区巫山段干支流水质变化分析研究. 水利发电学报, 2010, 35(5): 124-125.]
[36]
Peng FL, He LH, Yu Y et al. Studies on the total nitrogen, total phosphorus and chlorophyll a variations in the mainstream and main tributaries of the Yangtze River before and after the impoundment in the Three Gorges Project area. Scientia Sinica Technologica, 2017, 47(8): 845-853. [彭福利, 何立环, 于洋等. 三峡库区长江干流及主要支流氮磷叶绿素变化趋势研究. 中国科学:科技科学, 2017, 47(8): 845-853.]