1. 引言
在目前所应用的分布式水文模型中,美国农业部(USDA)开发研制的SWAT模型比较成熟 [1] [2] [3] 。SWAT模型是基于连续时段和流域尺度的综合水文模型。在国外,Srinivasan R.和Arnold [4] 在德克萨斯州的Seco Creek盆地上游应用SWAT模型进行了生态水文模拟;Andrea [5] 和Tripathi [6] 分别应用SWAT模型在德国中部Aar流域和印度进行了径流过程以及非点源污染负荷模拟,均取得了满意的结果。近些年,随着地理信息系统(GIS)的发展SWAT模型在国内一些流域也得到了较好的应用,比如洱海流域 [7] 、长江流域 [8] [9] 、淮河流域 [10] 等都有SWAT应用的研究案例。
洞庭湖区垸内资源性缺水、水质性缺水日益严重,开展垸内水体生态需水的研究迫在眉睫。本研究以大通湖垸为研究对象,利用大通湖区的气象数据、水文数据、土地利用类型图以及土壤类型图等数据资料,建立SWAT模型,通过模拟计算得出不同水期的生态补水量。
2. 数据与方法
2.1. 研究区概况
大通湖垸 [11] 为湖南省洞庭湖区11个重点垸之一,位于洞庭湖腹地,行政区域涉及大通湖区、沅江市以及南县。大通湖垸人口较为稠密且分布不均,湖区周边人口密集,低山丘陵地区人口密度相对较小。地形为典型的洞庭湖冲积平原,区内的湖泊、水面、沟渠纵横交错。土壤成土母质系河湖沉积物,地势低洼平坦,地面海拨高程24~32 m之间。区境属中亚热带向北亚热带过渡的大陆性季风湿润气候区,年平均气温16.5℃,无霜期264天;年均降水1240.8 mm。
2.2. 基础数据
SWAT水文模型的构建需要输入空间数据和属性数据,其中空间数据主要包DEM、土壤类型和土地利用数据等;属性数据主要包括水文数据和气象数据(见表1) [12] 。
2.3. SWAT分布式水文模型
利用大通湖垸所处区域2003年~2012年径流以及水质数据对模型进行率定,采用垸内实测的水文、水质数
Table 1. Basic data of ecological runoff model in Datong Lake
表1. 大通湖垸生态径流模拟基础数据
据进行验证。模型将大通湖流域划分为多个子流域,每个子流域由不同的水文响应单元(HRU)组成,利用空间数据和属性数据信息模拟多种不同的水文物理过程,包括径流、泥沙、化学物质的转移和转化过程。
模型的率定与验证采用专门针对SWAT模型开发的提供敏感性分析、率验证和不确定性分析的SWAT-CUP计算程序,评价指标包括相对误差(Re)、决定系数(R2)和Nash-Suttcliffe系数(Ens) [13] 。其中,Re的绝对值越小,模拟效果越好,反之,模拟效果则较差;R2值越接近于1,说明模拟值与实测值越接近,即模拟值对实测值的解释度越高;Ens的取值范围为0~1,值越大,则模拟值越接近观测值,模拟效果越好,反之,则越差。
2.4. 生态补水量计算
大通湖垸地区水环境污染物包括总磷、总氮、氨氮、氯化物以及化学需氧量,污染负荷比在10%以上,其中大部分来源于氨氮和总磷 [14] ,在此次研究中污染物总量控制主要考虑氨氮和总磷2种污染源。
生态补水量除了受垸内氨氮和总磷污染物量的影响,还与垸内径流和补给水源中污染物浓度有关,由于不同时段内垸内径流量和补给水源中污染物的浓度差异较大,因此,本研究我们分别计算了枯水期(9~12月份)、平水期(1~4月份)和丰水期(5~8月份)的生态补水量。
3. 结果与讨论
3.1. SWAT模型模拟
利用SWAT水文模型自有的敏感性数据分析模块,对大通湖垸内径流模拟的参数进行敏感性分析,遵循先调整地表径流,再调整土壤水、蒸发、地下径流以及降雪融雪参数的原则发现:基流α系数、径流曲线系数、土壤可利用水量、土壤蒸发补偿系数和浅层地下水径流系数的变化对研究区径流模拟结果的影响较大,因此在率定过程中要对这些参数进行调参。敏感性排序前5位的参数如表2所示。
SWAT分布式水文模型和氨氮、总磷负荷模型的塞阳河典型断面的相关系数(R2)和纳什效率系数(Ens)计算结果见表3和图1。
采用纳什效率系数Ens和相关系数R2来分别反映模型观测值和模拟值的拟合度和相关程度。径流、氨氮和总磷模拟验证期间的相关系数R2都接近了0.6,说明模拟值与实测值总体趋势一致,模拟效果较好。总磷负荷模拟的相关系数则相对较低,可能是由于大通湖垸区域面积小,垸内缺乏气象站和水文、水质站,率定模型相关监测数据来源于接近的周边地区,验证数据根据典型气象条件在垸内测量得到。尽管如此,总体而言模拟结果基本都在许可的范围之内,这表明SWAT模型对于平原地区小流域的径流过程和非点源氮磷负荷模拟也有一定的适用性。
3.2. 生态补水量
为探讨大通湖垸内维持生态环境降低水体中氨氮和总磷浓度的生态需水量,此研究在SWAT模型计算的基
Table 2. Definition and Range of parameters
表2. 参数率定及取值变化范围
Table 3. Results of calibration and verification for regional flow and non-point source nitrogen and phosphorus model
表3. 研究区域流量与非点源氮磷模拟校准验证结果表
Figure 1. Measured and predicted results in the cross section at Saiyang River
图1. 塞阳河断面实测与计算结果
础上,计算了大通湖垸内维持水体正常功能所需的最小生态需水量,污染物控制需水量可直接从南洞庭湖引水。根据物质守恒原理,要保持输出的污染物总量与输入的污染物总量相同,则有:
(1)
式中:W为生态补水量,W0为大通湖地区的径流量,P为达到地表水环境质量标准中三类标准的污染物浓度,M为大通湖地区目标污染物的总量(包括点源和非点源),Pe为从南洞庭湖补进来的水中污染物的背景浓度。
各不同水文时期的径流量W0和氨氮、总磷的总量M通过SWAT模型计算获得。通过式(1)可以计算大通湖垸氨氮和总磷达到III类水标准的生态补水量:枯水期生态补水量为0.68亿m3,平水期生态补水量为1.02亿m3,丰水期生态补水量为1.71亿m3。
从结果中可以看出,尽管丰水期大通湖垸内径流较大,但生态补水量比枯水期和平水期大,这主要是因为丰水期(5~8月份)汇水区的氨氮和总磷产生量较大,这和当地的农业种植结构有关,化肥的大量使用导致水体中氨氮和总磷的量增加。枯水期虽然氨氮和总磷的总量较小,但是由于9~12月份垸内的径流量较小,导致水体污染物浓度达不到标准,同样需要外源补水。通过计算枯水期、平水期和丰水期不同时段内大通湖垸的生态需水量,可以为补水及垸内水环境治理提供参考。
4. 结论
1) SWAT模型在大通湖垸的径流过程以及非点源氮磷污染负荷模拟相关系数R2分别接近0.80和0.70,这表明SWAT模型也可以应用于平原地区小流域的水文水质模拟。
2) 计算得出枯水期生态补水量为0.68亿m3,平水期生态补水量为1.02亿m3,丰水期生态补水量为1.71亿m3。研究可以为垸内水环境改善提供支持。
基金项目
湖南省水利科技项目,国家自然科学基金(51579094, 51679082, 51521006, 51479072)。