目前针对水环境容量的总量控制大多根据90%保证率最枯月降雨条件下的设计流量进行计算^[1], 这样可以保证河流在枯水期水质达标, 也可避免水环境保护的“不足”。实际上河流水文条件具有动态变化特征, 当水量增大(如丰水期)时, 按最枯设计流量条件下的水环境容量进行污染物总量控制, 会造成水环境“过保护”, 浪费一部分水环境的纳污能力^[2]。流域降雨径流形成的非点源污染的季节性差别很大, 目前尚无公认合适的水文设计条件, 同时适用于丰水期和枯水期的水环境容量计算, 以便进行合理的容量总量控制。流域水环境容量在空间上也呈现一定差异^[3]。现有研究侧重于干流的环境容量分析, 由于水文资料的缺乏, 很多流域对支流研究少, 在支流较多的情况下, 容量分配易出现不合理的情况^[4]。

掌握水环境容量时空分布特征是合理分配流域水环境容量的基础。但是, 在水文资料缺乏的流域, 无法利用实测的径流资料开展水环境容量时空分布的分析。

本文以漠阳江流域为对象, 在现有水文资料的基础上, 利用HSPF水文模型模拟得到各子流域不同时间尺度的径流量, 然后利用一维稳态水质数学模型计算相应径流量下的漠阳江流域水环境容量, 进而开展水环境容量时空分布特征的分析。研究结果可为科学评价与合理利用漠阳江流域的水环境容量以及确定流域水环境承载力提供依据。

1 漠阳江流域概况

漠阳江(图 1)发源于广东省云浮市西南大云雾山南侧, 在阳江市的北津注入南海, 干流全长199 km, 流域总面积6091 km², 其中阳江市境内流域面积5604.4 km², 占总流域面积的93.1%, 境内流域面积占阳江市总面积的71.7%。集水面积超过100 km²的一级支流有11条, 其中潭水河、西山河和那龙河集水面积较大; 二级支流有6条。漠阳江流域雨量充沛, 年均降水量为2195 mm, 降水年内分配不均, 年际变化较大。年内分配具有干湿季分明、雨量集中于汛期的特点。

漠阳江是阳江市最主要的饮用水源, 也是流域内城乡工农业生产的最重要水源。漠阳江干流下游水质控制目标为Ⅲ类水, 上游与中游以及漠阳江各支流控制目标为Ⅱ类水。目前漠阳江全河段水质虽以Ⅱ类为主, 但水污染防治形势不容乐观。近年来, 社会经济的持续高速发展, 对漠阳江水质产生较大的影响, 水污染问题也逐渐显现, 水环境承载力的压力越来越大。如果没有科学系统的水环境综合管理方案, 流域水质保护的压力将越来越大, 也将对漠阳江供水水源水质安全构成严重威胁。以分区域水环境容量为基础, 单位面积的水环境容量可以表示水环境对区域社会经济发展的支持能力, 即水环境承载力的高低。

2 数据和方法

漠阳江地理信息数据库包括DEM (分辨率为30 m美国NASA的ASTER GDEM)、河网(广东省科学研究院提供, ArcGIS矢量数据)、土地利用矢量数据(广东省科学研究院提供, 包含水田、旱地等16类, 图 2)以及在此基础上划分的31个子流域(在ArcGIS中根据DEM图与河网图手动划分)。

气象数据来自国家气象信息中心, 包括1983-2012年阳春、阳江两个站(图 1)逐日(20-20时)降水量、日最高气温、日最低气温、露点温度、云量、风速、相对湿度等, 用于水文模型建立时气象数据的输入。

水质监测数据来自阳江市监测站, 包括2005-2011年7个监测断面(春湾、陂面、鱼皇石、中郎、河口、漠阳江水厂、江城)的数据(1-11月, 每月1至6个监测日), 指标包括COD (化学需氧量)、氨氮、重金属等24项, 用于水质模型的校准与验证。

2.1 水文模型

2.1.1 流域水文条件模拟

由于连续的实测水文资料缺乏, 考虑到各子流域空间异质性, 选用分布式水文模型HSPF模拟漠阳江各支流的径流量。在对HSPF模型完成验证的基础上, 将模拟得到的水文数据(各个子流域30年的逐月流量数据)用于分析研究漠阳江流域水环境容量时空分布。

HSPF模型中水文模拟以Stanford水文模型为基础, 根据土地利用特征, 分为透水地面和不透水地面, 针对不同下垫面的水文过程进行模拟。透水地面的模拟考虑降水、截留、地表填洼、蒸发、蒸腾、渗透、地表径流、壤中流和地下水流等水文过程; 不透水地面的模拟包括降水、截留、地表径流、蒸发。HSPF模型水文模块对于降雨径流复杂过程的描述, 能较好地再现流域实际水文过程^[5-6]。

国内外利用HSPF模型对降雨径流为主的水文过程、气候变化以及土地利用变化对水文过程的影响开展了很多研究, 表明HSPF模型在诸多地区水文过程的时空分布变化模拟中均有很好的表现, 石赟赟等^[7]、李燕等^[8]、Cryer等^[9]、Alarcon等^[10]和董延军等^[11]在不同时间与空间尺度上模拟所得结果的相对误差均小于15%, NASH系数都在0.8以上, 模拟效果良好。

2.1.2 漠阳江流域HSPF水文模型验证

漠阳江流域面积较大, 水文过程复杂, 水文敏感参数较多, 对河流流量影响较大, 水分平衡以及水文参数对整个流域非点源污染负荷模拟起着至关重要的作用。本文参考HSPF模型手册^[12]和HSPFParm^[13]以及相关研究成果设定参数初值, 然后对参数进行率定, 水文部分主要参数的最终取值见表 1。

利用2009年漠阳江流域水文观测资料对水文模块进行参数率定, 根据HSPF模型指定的参数范围, 反复调整参数, 使流域年出流总量的模拟值和实际出流总量趋于一致。以2010年实测资料验证模型, 水文模拟首先从年度水分平衡开始, 初步模拟年度水量低于测量值, 因此适当减小下层额定存储量(LZSN), 减少潜在蒸散, 增加产流, 从而增加降水引起的径流水量。年水量模拟结果中, 2009年模拟值为190.95 m³/s, 实测值为203.28 m³/s; 2010年模拟值为228.91 m³/s, 实测值为259.12 m³/s; 2009年相对误差为6.06%, 2010年的相对误差为11.7%, 相对误差均小于15%。

2009年和2010年的月水量模拟结果分析见图 3。除2009年5月模拟结果的相对误差为52.7%外(漠阳江地区有较多水库, 我们未得到相关数据或说明, 为了保持流域的完整性, 没有扣除水库所在子流域, 该月误差较大可能是由于水量激增引起的水库蓄水所致), 其他月份模拟值与观测值的相对误差均小于15%, Nash系数在0.9以上, 曲线走向一致, 较好地反映了漠阳江流域的水文过程。

2.2 水质模型与环境容量估算方法

2.2.1 水质模型的建立与验证

对于研究河流水污染迁移转化过程, 若主要关心污染物浓度的沿程变化, 对断面上的变化情况没有太多要求时, 可采用一维水质模型进行描述。水环境容量是在特定的水文设计条件下, 稳定的污染源和水文特征时, 河流能接受的最大污染物的负荷量。在河流流速、流量、水位都不随时间而变化, 污染源排放连续稳定情况下, 纵向分散作用远小于河流的移流和降解作用, 同时不同时刻的纵向分散作用在宏观上可以相互抵消。在这种特定情况下, 可以忽略纵向分散作用。考虑到国内, 尤其是漠阳江流域所在地区的水环境保护要求以及水环境容量计算的精度, 同时考虑到稳态模型发展的成熟性, 本文采用一维稳态模型计算水环境容量。

为保证水质模型能真实地反映污染物流入后的情形, 我们验证了所建立的水质模型。将199 km的漠阳江干流及主要支流划分为27个河段共81个计算单元(图 4), 每个计算单元的空间步长为1 km。以2010年漠阳江流域污染负荷和水文条件输入所建立的河流水质模型, 模拟出漠阳江干流沿程的COD和氨氮浓度值, 选取5个断面的实测值作为模型验证的对比资料, 选用相对误差(RE)和Nash-Sutcliffe标准来评估模型在校准和验证过程中的模拟效果。NASH公式计算如下:

$E = 1 - \sum\limits_{i = 1}^n {{{\left( {{Q_i} - {{Q'}_i}} \right)}^2}/} \sum\limits_{i = 1}^n {{{\left( {{Q_i} - \bar Q} \right)}^2}} $

(1)

Q_i为第i时刻观测流量(m³/s); Q'_i为第i时刻模拟流量(m³/s); n为时段总数; Q为平均观测流量。

NASH系数是描述计算值对目标值拟合精度的无量纲统计参数, 一般取值在-∞~1之间, 本研究中, NASH系数在0.8以上表明符合要求, 0.6~0.8为基本符合要求, 低于0.6则认为偏差过大, 结果不适宜。

参考中山大学环境科学研究所编制的《粤西主要河流水质保护规划修编(研究报告)》, 广东省粤西地区COD的衰减系数为0.08~0.1 d^-1, 氨氮的衰减系数为0.07 d^-1, 经过模型验证, COD的衰减系数K_COD取值为0.1 d^-1, 氨氮的衰减系数K_氨氮取值为0.07 d^-1。模拟结果见表 2。

将下游入海口作为距离基准, 为0 km, 上游发源地为199 km, 从春湾监测点附近开始模拟, 根据水质监测数据, 得到验证结果如图 5所示。

由图 5可以看出, 除中朗点氨氮实测值与模拟值误差较大(中朗位于阳春市中心的下游, 其实测数据具有很大的不确定性)外, 其余点的COD与氨氮的模拟值与实测值相对误差均在10%左右, Nash-Sutcliffe系数在0.8以上, 且实测值与模拟值的沿河变化趋势基本相同, 说明模拟结果符合漠阳江干流的水质变化特征。由此表明, 该一维水质模型可适用于漠阳江流域水环境容量的计算。

2.2.2 水环境容量估算方法

根据漠阳江流域的实际情况, 以COD和氨氮为研究对象。对于漠阳江干流的混合输移过程, 研究重点是污染物浓度的沿程变化, 因此采用一维水质模型模拟水质变化和计算水环境容量。

一维稳态水质数学模型假定各排污口连续、均匀地排污^[14]。根据旁侧入流和污染源位置, 将对象河段分成若干段, 旁侧入流和污染源位于段首, 每段保证流速、流量定常, 在同一段内, 采用《环境影响评价技术导则》(HJ/T 2.1-2.3-93)^[15]推荐的斯特里特-菲立浦(Streeter-Phelps)模式:

$C(x) = {C_0}\exp \left[ { - K\left( {\frac{{{x_{i - 1}}}}{{86400{u_{i - 1}}}}} \right)} \right]$

(2)

式中, C(x)和C₀分别为x=x和x=0处河水COD或氨氮浓度(mg/L), x为到排污口(x=0)的河水流动距离(m), u为河水平均流速(m/s), K为COD或氨氮的综合降解系数(d^-1)。

本研究采用段首控制方法, 即控制上游断面(即段首)的水质达到功能区段的要求, 这样, 由于有机物的降解, 则在该段内的下游水质处处达到或优于功能区段的控制指标。段首控制方法严格控制功能区段的水质不超标。

在功能区段的段首, 由于来水中污染物的浓度和功能区段水质要求的差别, 为来水提供稀释容量:

${E_0} = {Q_0}({C_{\rm{s}}} - {C_0})$

(3)

E₀为功能区段段首的稀释容量(t/d), C_s为功能区段水质标准(mg/L), Q₀为来水流量(m³/s), C₀为来水浓度(mg/L)。

第i个断面处的环境容量为

${E_i} = ({Q_i} + {q_i}){C_{\rm{s}}} - {Q_i}{C_{\rm{s}}}(x)$

(4)

E_i为第i个断面处的环境容量(t/d), q_i为第i个断面处的排污流量(m³/s), Q_i为混合后干流流量(m³/s), 其余符号意义同上。

功能区段内所具有的总环境容量为

$\begin{array}{l} E = {E_0} + \sum {{E_i} = {Q_0}} ({C_{\rm{s}}} - {C_0}) + \\ \;\;\;\;\;\;\sum\limits_1^n {\left[ {({Q_i} + {q_i}){C_{\rm{s}}} - {Q_i}{C_{\rm{s}}}(x)} \right]} \;。 \end{array}$

(5)

3 结果与讨论 3.1 径流变化特征

漠阳江流域径流量的年际变化主要受降水影响, 1983-2012年共30年的径流量结果统计得出:最大径流量为537.7 m³/s (2001年), 最小径流量为166.6 m³/s (1989年), 30年平均值为322.5 m³/s, 变差系数为0.29。

漠阳江流域年内丰水期与枯水期径流差异明显, 其中枯水期最大径流量为210.1 m³/s (1990年), 最小径流量为50.0 m³/s (2004年), 30年平均值为105.1 m³/s, 变差系数为0.4;丰水期最大径流量为990.2 m³/s (2001年), 最小径流量为235.2 m³/s (1999年), 30年平均值为540.0 m³/s, 变差系数为0.3。丰水期平均径流量是枯水期平均径流量的5倍, 且丰水期最大径流量与最小径流量与年际最大与最小径流量相同, 主要原因是丰水期降雨量大, 对漠阳江流域的径流量的产生量作用巨大。

漠阳江流域360个月中最小径流量为10.1 m³/s (2005年2月), 最大径流量为3042.3 m³/s (2001年6月), 月份差异巨大。

3.2 水环境容量的时间分布

3.2.1 年际变化

在完成HSPF模型验证之后, 模拟估算漠阳江流域各个子流域连续30年的径流量, 结合相关的河网、流速等物理数据, 利用验证后的水环境容量模型, 计算出漠阳江流域COD与氨氮的逐年环境容量变化曲线(图 6)。可以看出, COD与氨氮的容量走势与年际降雨量走势基本上一致, COD容量整体上是氨氮的25倍左右。1983-2012年共30年的COD容量最大值为27.42万t/a, 最小值为15.44万t/a, 平均容量为20.92万t/a, 变差系数为0.15;氨氮容量最大值为1.35万t/a, 最小值为0.45万t/a, 平均容量为0.83万t/a, 变差系数为0.27。

将30年COD与氨氮的容量数据对应的保证率进行统计, 结果如图 7所示。可以看出: 1)保证率随容量的增大而逐渐减小, 在3%左右的时候出现减缓拐点, 主要原因是高于25.5万t/a的容量(3%保证率)需要较为极端的气象情况(2001年降雨量较为极端)才会出现, 这个区间内容量的不确定性较大, 因此小于3%保证率的容量减缓趋势不明显; 2)在90% (传统的水文设计条件), 50%和10%保证率水文年下, 漠阳江流域COD的环境容量分别为16.45, 21.84和24.97万t/a, 氨氮的环境容量分别为0.51, 0.88和1.14万t/a; 3)利用该曲线可以在漠阳江流域科学地利用水环境容量, 以平水年2010年为例, 若全流域COD的排放量达到COD的容量(21.4万t/a), 则漠阳江流域超负荷排放的概率是50%。

3.2.2 季节变化

对漠阳江流域各个水文年枯水期(1-3月以及10-12月)与丰水期(4-9月)的环境容量进行统计, 如图 8所示。可以看出: 1)相同保证率的COD丰水期的容量都是枯水期容量的2倍左右, 氨氮丰水期的容量是枯水期容量的3倍左右, 季节性明显; 2)在90%(传统的水文设计条件), 50%和10%保证率水文年下, 漠阳江流域COD雨季的环境容量分别为9.3, 12.92和15万t/a, COD旱季的环境容量分别为4.69, 5.73和7.35万t/a, 氨氮雨季的环境容量分别为0.37, 0.69和0.9万t/a, 氨氮旱季的环境容量分别为0.08, 0.14和0.23万t/a; 3)保证率越小, COD与氨氮的容量在丰水期与枯水期的差异越大, 主要原因是漠阳江流域降丰水期节分明, 丰水期降雨量大, 导致丰水期比枯水期容量大, 而枯水期的降雨量的影响在于基流(基流在年际变化中保持相对稳定); 4) COD与氨氮在枯水期和丰水期的变化趋势基本上一致, 且容量与降雨量的关系并非简单的线性关系, 雨季的保证率所对应的容量变化幅度较大, 旱季较小。

3.2.3 逐月变化

对漠阳江全流域30年逐月的COD与氨氮环境容量进行统计分析, 结果见图 9。可以看出: 1)在逐月变化图中, 各月份平均数与中位数逐月变化走势一致, 以中位数为标准, COD与氨氮的容量在1月最小, 分别是0.7和101.84 t/a, 1至6月逐月升高, 并在6月达到最大容量, 分别是2.4和1385.27 t/a, 然后在6至12月逐月减少; 2)对于COD容量的12个箱体数据, 丰水期(4-9月)的容量较大, 枯水期(1-3月, 10-12月)的容量较小, 且丰水期容量数据之间的变幅大, 变差系数也大, 丰水期容量最大值为5.44万t/a (2001年6月), 最小值为0.65万t/a (1994年4月), 极差为4.79万t/a, 变差系数为0.44;枯水期容量最大值为1.97万t/a (1990年2月), 最小值为0.63万t/a (2005年2月), 极差为1.34万t/a, 变差系数为0.34。主要原因是丰水期雨量较大, 丰水期容量的影响因素更多, 更为复杂。

3.3 空间分布

选取平水年2010年为空间分布研究的基准年, 漠阳江各个子流域COD与氨氮单位面积的环境容量统计结果见图 10。可以看出以下几点。1)单位面积COD容量与单位面积氨氮容量的空间分布基本上相同。2)以COD为例, 漠阳江干流的沿河变化是由漠阳江上游的95.27 t/km²减小到漠阳江中游的51.62 t/km², 再增大到漠阳江下游的77.32 t/km²。上游到中游COD的减小主要是由于漠阳江上游水质较好, 而在上游及中游的各条支流汇入一定浓度的COD后, 中游能容纳的COD量变小, 而面积又增大, 因此单位面积容量减小。中游到下游COD的增大是因为下游的水质目标是三类水, 而上游与中游的水质目标是二类水, 所以下游能容纳更多的COD, 但因为其面积较大, 所以其单位面积的COD容量比中游大, 但没有上游大。3)除去水库以及漠阳江的干流流域, 漠阳江一级支流与二级支流单位面积的COD容量较小。一级支流中西山河下游流域容量(77.83 t/km²)最大, 主要原因是西山河下游流域有多条支流汇入, 流量大, 坡度大, 且流域面积小; 大八河流域容量(5.37 t/km²)最小, 主要原因是大八河流域流量小, 坡度小, 面积大, 且城镇及居民用地所形成的不透水地面面积大。对于漠阳江流域其他一级二级支流, 其承载力(以单位面积的环境容量表示)的大小由流域面积、流量、坡度以及土地利用等因素综合决定。

4 结论

根据漠阳江的实际情况, 在水文资料不足的情况下, 本文利用水文模型HSPF实现对各个子流域时间空间上水文数据计算, 并将水文模型与一维稳态水质数学模型相结合, 计算得到流域水环境容量总量的时空分布, 结论如下。

1) HSPF模型对于水文过程的模拟结果中, 2009年的相对误差为6.06%, 2010年的相对误差为11.7%, 月模拟值与观测值相对误差均小于15%, Nash系数在0.9以上, 曲线走向一致, 可较好地反映漠阳江流域的水文过程。一维稳态水质数学模型的对比资料显示模拟值与实测值相对误差均在10%左右, Nash-Sutcliffe系数在0.8以上, 且浓度沿河变化趋势基本上相同, 模拟结果符合漠阳江干流的水质变化情况。

2)漠阳江30年的容量统计分析显示:漠阳江流域雨量充沛, COD与氨氮环境容量很大, 30年容量变化曲线与降雨变化曲线走势基本上一致; 枯水期与丰水期季节差异明显, 丰水期容量远大于枯水期容量。对选取的5个代表性水文年的容量统计结果分析得出:漠阳江流域环境容量随保证率增大呈现增长减缓的趋势; 漠阳江流域丰水期容量变幅较大, 变差系数大; 漠阳江流域最枯月为1月, 最丰月为6月。

3)以单位面积的环境容量表示漠阳江流域各个子流域的承载力, 统计结果表明:漠阳江干流流域承载力较大, 漠阳江干流的一级与二级支流的承载力较小; 其水环境承载力由流域面积、流量、坡度以及土地利用等因素综合决定。