1. 引言
流域污染负荷包括点源和非点源两部分。目前,我国几乎没有系统的长系列入河点源和非点源污染监测资料。河流水质污染监测则主要在可控制的河道断面上,监测数据中的污染物既有来自点源污染物, 也有非点源污染物 [1] 。如何在现有监测资料的基础上确定点源和非点源污染负荷各自比重,明确河流主要污染负荷的贡献类型,是一项具有重要现实意义的工作,将对污染源监控和管理有重要作用,并直接决定污染治理投资的方向及强度。
河流污染物通量可为水污染控制提供重要信息和决策依据。目前,关于河流污染物通量的估算方法大致有四种类型,分时段通量和 [2] [3] [4] [5] 、时段平均浓度与时段水量之积 [5] [6] 、通量频率分布之和 [7] 、对流–扩散模式 [4] [8] 等。前2种类型较为常用,后者比前者粗略。Webb等 [9] (1997)采用实测断面瞬时浓度、时段平均浓度、断面瞬时流量、采样期间平均流量、时段平均流量等构造了5种时段通量的计算公式,并进行了评估。结果显示各种方法估算的结果相差较大。富国 [10] (2002)通过对5种估算方法的误差比较分析,指出时均离散通量对估算误差具有一定的贡献,在选择时段通量估算方法时应考虑污染源类别、污染物特性等因素影响。表1列出了5种时段通量估算方法、特点及应用取向分析。
根据河流断面污染物通量可以确定实际点源和非点源对通量的贡献,对此国内外学者已开展一些研究 [11] [12] [13] [14] 。常规点源和非点源分割方法是枯水期平均通量法:采用枯水期河道平均通量作为点源贡献进而推广到全年。该方法基本原理是认为非点源污染主要是由汛期地表径流引起的,枯水期河水水质主要反映点源的污染情况,可以忽略非点源贡献。其不足之处在于,由于枯水期点源的流达率最低,因此将枯水期平均通量作为点源平均通量估算到全年,则可能导致点源年通量的低估。袁宇等 [14] (2008)从实用的角度提出了以月径流量与月通量相关系数确定丰水期通量增量的非点源比例系数,建立了河流年通量的点源和非点源分割公式,但该方法在将枯水期平均月通量作为点源通量时未考虑背景负荷(可视为非点源部分)贡献,可能一定程度上导致点源年通量的高估。
针对上诉情况本研究总结前人工作的基础上,建立了一个简单的适合点源和非点源多种情况的时段通量估算方法,改进了传统的点源和非点源通量分割方法,提出了水污染负荷贡献类型的判别标准,并选取岷江流域主要控制断面进行了实例分析。
2. 研究方法
2.1. 河流断面污染物时段通量估算方法
本研究建立了一个简单的适合点源和非点源多种情况的时段通量估算方法,见式1。
(1)
式中,
为瞬时浓度;
为瞬时流量;
为时段平均流量;
为估算时间段内的样品数量;
为时段长度;
为权重系,根据
和
的相关系数
确定,若
和
相关系数
,则
;若
和
相关系数
,则
;若
和
相关系数
,
。
2.2. 河流断面通量点源和非点源的分割方法
本研究采用的点源和非点源通量分割方法如下:
枯水期月平均通量
:
(2)
点源年通量
:
(3)
非点源年通量
:
(4)
Table 1. The estimation methods, features and applied range analysis of the river period-fluxes
表1. 河流时段通量的估算方法、特点及应用取向分析 [10]
备注:
代表估算时间段内的样品数量;
为估算时间段转换系数(取时段长度)。
以上各式中,
为枯水期第
月的月通量;
为丰水期或平水期第
月的月通量;
为枯水期月数;
为丰、平水期月数;
为扣除背景负荷的调整系数,根据枯水期平均浓度
和背景浓度
确定,即
;
为丰水期点源通量增量调整系数,根据丰水期通量
和径流量
的相关系数
确定,即
。
2.3. 河流断面水污染负荷类型的判别标准
在上述断面通量估算及点源和非点源分割方法的基础上,本文将河流污染负荷贡献类型划分成5种类型,分别为点源为主型、点源占优型、混合型、非点源占优型、非点源为主型。表2给出了以点源贡献比例为判据的水污染负荷贡献类型的判别标准。
3. 应用实例
3.1. 岷江主要断面污染物月通量
本研究方法应用于长江主要一级支流岷江,表3为2013年岷江主要控制断面高锰酸盐指数和氨氮
Table 2. Description of water pollution load types
表2. 水污染负荷类型分类
Table 3. Estimation of fluxes per month of Min River in 2007
表3. 2013年岷江主要控制断面污染物月通量估值
污染物月通量估值。可以看出,在2013年6~9月的丰水期内,各断面高锰酸盐指数的月通量较高,与径流量有较高的相关性。表4为各断面高锰酸盐指数和氨氮指标月均浓度和月通量分别与径流量的相关系数,由表4可知,各断面高锰酸盐指数的月通量与径流量的相关系数分别为0.94,0.94,0.89和0.92,说明岷江高锰酸盐指数通量与径流量关系密切。高场和大渡河断面氨氮丰水期(6~9月)月通量较高,与径流量的关相系数为0.90和0.76。而彭山和五通桥断面氨氮月通量和径流量的相关系数较低,为0.62和0.65,这一差别应该是这两个断面氨氮指标的点源贡献率较高所致。这一结论可由表4中显示的彭山和五通桥断面氨氮月均浓度和径流量呈现较好的负相关性得到进一步证实。
3.2. 污染负荷贡献类型分析
根据前文确定的河流断面水污染负荷类型的判别标准,我们对四个断面高锰酸盐指数和氨氮的污染负荷贡献类型进行分析,结果表明,就高锰酸盐指数而言,彭山断面属于混合型,为点源和非点源共同作用。五通桥、高场及大渡河断面都显示以非点源占优。氨氮在彭山断面显示出明显的点源贡献为主,五通桥断面也显示为点源占优,高场断面属于混合型,大渡河断面表现出非点源贡献较高。总的来说,流域高锰酸盐以非点源贡献为主,氨氮以点源贡献为主(表5)。
彭山断面上游的成都市是岷江流域污染负荷最大的输入区,是流域高锰酸盐指数和氨氮点源污染负荷的主要输入区,致使流域成都以下江段水质受到污染。五通桥和高场断面在大渡河支流汇入岷江干流以下,大渡河水质以非点源来源为主,其水量在一定程度上对岷江干流彭山断面下来的物污染物进行了稀释。
4. 结论和建议
1) 本文在前人工作的基础上建立了一个简单的通量计算方法,该方法适合点源和非点源多种情况的时段。该方法的下一步工作是利用自动监测站数据,验证其准确性,并与其他方法进行对比分析。
2) 本文在概念上明确点源和非点源分割改进的方法,从实用角度提出以月径流量与月通量相关系数确定丰水期点源通量增量调整系数。该方法的主要改进方向在于丰水期点源通量增量的调整系数更精确的估算。
Table 4. Linear correlations of runoff volume per month and concentration and flux per month of Min River
表4. 岷江2013年月径流量与各污染物月均浓度和月通量的线性相关系数
Table 5. Distribution of point souse and non-point source
表5. 点源与非点源对通量的贡献率及类型分析
3) 岷江流域主要断面实例计算分析表明本文提出的负荷估算及贡献类型识别方法的实际效果。流域高锰酸盐污染主要以非点源贡献为主,氨氮以点源贡献为主,结果可为流域水污染控制规划提供科学依据。
项目基金
国家水体污染控制与治理科技重大专项2017ZX07301。
NOTES
*通讯作者。