基于SWAT模型的流域径流模拟.docx

??

?

??

基于SWAT模型的流域径流模拟

?

??

?

?

?

?

?

?

?

???

?

?

?

?

?

摘要：研究基于GIS平台，建立了冷水河流域地理信息空间数据库及SWAT模型基本属性数据库，将流域划分为31个子流域，223个水文响应单元（HRU）。利用SWAT-CUP模型的SUFI-2算法，进行了模型的校准与验证，表明模型适用于冷水河流域。模型设定了不同土地利用状况研究土地利用方式的转变对流域内径流的影响，模拟了不同气候条件下区域径流的变化，为区域水资源开发利用、防洪减灾等提供了科学支撑。

关键词：冷水河流域；SWAT模型；SUFI-2算法；防洪减灾

1.研究区基本情况

南郑县冷水河流域地处米仓山主脊以北，全长53.5km，平均比降29.9‰，流域面积641.41km2。研究区年平均气温14.9℃，年平均降水970.4mm，降水主要集中在7～9月份。近年来由于汉中地区暴雨中心南移，研究区南部山区碑坝、牟家坝等地，降雨强度大，致灾频率较高。

2.模型建立

冷水河流域模型子流域个数为31，对土地利用类型、土壤类型、坡度类型进行重分类，共划分223个水文响应单元[1]。

本次径流的模拟主要采用日降水数据，通过SCS径流曲线法进行模拟，径流演算主要以日为时间步长；模型主要采用偏正态分布模拟产生日降水量；基于前期气象资料，本文选择Penman—monteith法来模拟潜在蒸发量；由于本次河道较顺直且宽浅，河道演算方法选取马斯京根法[2]。

3.模型校准与验证

3.1模型敏感性分析

选取冷水河流域2000～2014年月径流数据对模型进行参数敏感性分析，对相关参数进行不确定性分析，通过SWAT-CUP模型SUFI-2算法调参计算，设定调参调用SWAT模型运行2000次，得到模型建立的7个主要影响参数[3]，CN2、GW_DELAY、ESCO、CH_N2、ALPHA_BF、SURLAG、SFTMP。

3.2模型校准及验证

选取冷水河流域第10个子流域的逐月流量来进行模型的校准及验证。校准及验证资料为三花石水文站2008～2014年实测逐月流量资料。

根据现有资料，模型校准选取三花石水文站2008～2011年4年逐月流量资料，进行模型校准。校准期月径流的效率系数Ens为0.85，模型相关系数R2为0.97，P因子和R因子分别为0.49和0.29；研究采用2012～2014年3年的逐月数据进行径流验证，验证期模型效率系数Ens为0.89，模型相关系数R2为0.82，P因子和R因子分别为0.45和0.51。通过模拟图形及验证系数表明SWAT模型对冷水河流域比较适合，适宜在此参数基础上建立模型。

4.环境变化条件下的径流模拟

4.1土地利用状况变化

土地利用状况变化条件模拟主要设定3种情况。以现有的土地利用状况为基础，模拟土地利用方式恶化及水保措施等相对较好的情况，研究土地利用变化对流域水循环的影响[4]。

从冷水河流域多年平均流量的模拟结果来看，土地利用方式的转变对流域内径流的年际径流变化影响较小。当水保措施较好的情况下，流域内大面积林地及草地对水源的涵养作用明显，遏制了水土流失，有效储存了水量，年径流量相对较小为10.95m3/s；而当林地大面积消亡，植被覆盖率降低，植被储存水量不断减少，水土流失严重，且居民建筑用地及耕地不断增加的情况下，导致部分裸地显现，加剧了水土流失，径流增大为11.53m3/s。

从冷水河流域多年月平均流量的模拟状况来看，土地利用方式的转变对流域径流月变化影响较为明显。雨季（4月～10月）较枯季（11月～次年3月）径流的变化显著。在枯水期土地利用方式的转变对径流的影响程度相对较小，相反，在雨季时土地利用方式的转变对径流影响程度较大。雨季的初期土地利用方式的转变相对雨季后期对径流的影响程度大，径流的增幅更为显著。

由土地利用方式转变情景模拟结果的年际及年内的分布来看，植被覆盖率，对流域径流的影响明显。林地的涵养水源作用较为明显，森林植被的大面积分布能够有效的储存水量，减缓坡面径流流速，通过增加土壤的下渗量，减少水土流失。林地的增加也可消减汛期的洪峰，减少流域洪水风险，同时汛期的截存水量在枯水季节可补给径流，解决枯水季节的干旱灾害。林地的减少会加剧径流的年内分配不均，更易造成区域旱涝灾害。因此对于流域内植被的保护及土地资源的合理利用，不但关系着区域水资源含量的多少，同时也是减少区内洪水及干旱灾害的一项重要的措施。

4.2气候变化

结合南郑县近50年气候变化规律，按照气温在现状情况下分别+2℃、+1℃、不变、-1℃、-2℃，降水在现状情况下分别+20