论文部分内容阅读
随着点源污染得到有效治理,非点源污染已成为许多地区水质超标的主要原因。与点源污染不同,非点源污染受诸多环境因素的影响,发生和发展过程具有变异性、随机性、滞后性等特点,对非点源污染的定量研究和有效治理仍面临较大挑战。定量识别流域非点源污染氮磷输出量和入河量是有效控制非点源污染的关键基础。本研究以浙江省绍兴市长乐江流域为研究对象,基于2003~2016年间的水文水质监测数据,建立了长乐江氮磷污染的LOADEST模型,模拟了河流氮磷污染的连续变化趋势;采用小波分析方法,分析了河流氮磷负荷量的周期性变化规律;定量分析了驱动流域非点源污染发生的降雨因子和影响流域非点源污染发展过程的地形因子、地表径流因子、地下蓄渗因子和植物截留因子;最后,采用基于河流一维水质方程-输出系数耦合方程和影响污染物入河过程的环境影响因子,估算了长乐江流域氮磷污染物的入河量和输出量。主要研究结果如下:(1)在2003~2016年间,长乐江平均流量为16.39m3·s-1,河流水质受到氮、磷等营养物质的明显污染。河流出口监测断面总氮(TN)多年监测浓度变化范围为0.20~13.29mg·L-1,多年平均浓度为3.75±0.91mg·L-1,远超该地区参照基准值(1.70 mg·L-1);总磷(TP)多年监测浓度变化范围为0.01~0.69 mg·L-1,多年平均浓度为0.14±0.09mg·L-1,处于地表水III类水水平;LOADEST模型估算结果表明长乐江TN的年均负荷量为2208±862t·yr-1,TP年均负荷为71±32 t·yr-1。河流的TN和TP负荷量都与径流量高度相关,但与TN和TP浓度的相关性较弱。(2)2003~2016年期间,长乐江TN和TP负荷量存在明显的周期性变化,TN和TP负荷量都以18个月的时间尺度为主发生着多时频的周期性变化,但TP负荷量的变化强度要比TN负荷量的变化强度弱,且两者的小波系数频谱等值线图与河流流量的相似性较大,都存在相同的时间尺度和震荡发生日期;进一步对TN和TP负荷量的小波方差展开分析,揭示了TN和TP负荷量周期性的变化是由河流流量的周期变化引起的,而与TN和TP的浓度变化相关性较小。(3)驱动流域氮磷非点源污染发生的降雨因子,以及影响流域非点源污染发展过程的地形因子、地表径流因子、地下蓄渗因子和植物截留因子均表现出一定的空间差异性:TN的降雨驱动因子变化范围为1.15~1.24,均值为1.20,TP的降雨驱动因子变化范围为1.22~1.33,均值为1.28;地形因子变化范围为0~2.22,均值为0.77;SCS-CN模型方法计算得到的TN地表径流值在1207~1597mm之间,均值为1463mm,以土壤流失通用方程计算得到的TP年侵蚀量在0~1433t·hm-2·a-1,均值为140t·hm-2·a-1;地下渗蓄因子的变化范围为1.00~1.04,均值为1.01;TN、TP植物截留因子均值分别为3.66和3.37;尽管不同环境影响因子之间存在强烈的共线性,但是考虑到这些环境影响因子与污染物运移关系的直接性,认为地表径流和植物截留是长乐江流域非点源污染发展过程的主要环境影响因子。(4)2003~2016年期间,长乐江流域水田、旱地、人居地、林地和荒地的TN输出系数分别为26.34±0.41kg·ha-1·yr-1、64.79±0.40 kg·ha-1·yr-1、42.97±0.38 kg·ha-1·yr-1、9.00±0.78 kg·ha-1·yr-1和6.38±3.47kg·ha-1·yr-1,TP输出系数分别为2.29±0.44kg·ha-1·yr-1、2.28±0.44kg·ha-1·yr-1、4.57±0.29 kg·ha-1·yr-1、1.38±0.88 kg·ha-1·yr-1和0.57±0.89kg·ha-1·yr-1;流域TN、TP入河量分别为1745t·yr-1和128t·yr-1。对不同土地利用类型的TN、TP入河量进行来源解析,其中对TN入河量贡献率从大到小排序为:旱地(51%)、水田(20%)、林地(16%)、人居地(11%)和荒地(1%),对TP入河量贡献率从大到小排序为:林地(34%)、旱地(25%)、水田(24%)、人居地(16%)和荒地(1%)。