基于区域径流的swa分布式水文模型.docxVIP

基于区域径流的swa分布式水文模型 摘要：通过构建适用于闽江流域的SWAT分布式水文模型,结合情景设置法,分别模拟研究区内不同土地利用情景下的径流过程,以定量分析土地利用变化对流域径流影响。结果表明:(1)经过参数率定的SWAT模型,率定期和验证期的评价指标 土地利用作为人类对自然利用的一种途径闽江流域作为福建省重要的生态单元,流域内生态保护区和水源保护区众多,水资源变化直接关系到全省1/3居民的生产生活1 研究区概况闽江流域位于东经116°23′—119°35′,北纬25°23′—28°16′。闽江作为福建省最大的河流,发源于武夷山脉,流经省内36个县(市),流域面积59 922 km2 研究方法2.1 基础数据本研究所采用的数据有地形数据、土地利用类型数据、土壤类型数据和气象数据等。数字高程模型(DEM)数据是由地理空间数据云平台提供下载的30 m分辨率ASTER GDEMV2数据,后期通过ENVI图像处理软件校正拼接而成(图1)。土地利用数据是清华地学系Gong等2.2 模型构建首先基于闽江流域DEM数据进行河网提取,考虑到模型运行效率,确定最小集水面积阈值为300 km2.3 模型参数率定与验证SWAT模型进行径流模拟时以大量相关参数为基础,为提高模型的适用性,本文借助SWAT-CUP软件在参数敏感性分析的基础上采用SUFI-2算法 式中:2.4 土地利用情景设置为定量研究闽江流域不同土地利用类型的径流响应情况,在模型模拟精度达标的基础上,分别设置5种土地利用情景,各情景土地利用面积变化统计见表1。(1)以闽江流域土地利用类型现状作为基准情景S0,方便与其他情景进行对比分析。(2)根据研究区丘陵山地面积广大,耕地后备资源相对匮乏的情况,为保障粮食安全,设置耕地储备情景S1:将基准情景中坡度15°的裸地、草地、林地类型转化为耕地类型。(3)根据福建省海峡西岸经济区城镇化、工业化迅速发展导致人地矛盾突出等情况,设置建设开发情景S2:将基准情景中坡度25°的林地、裸地、草地类型转化为建设用地类型。(4)根据研究区内生态保护区和水源保护区众多的特点,为保障当地生态安全,设置植被恢复情景S3:将基准情景中坡度≥25°的裸地、耕地、建设用地类型转化为林地类型,坡度25°的裸地类型和坡度介于15°～25°的耕地、建设用地类型转化为草地类型。(5)根据《福建省土地利用总体规划(2006—2020年)》设计分区调控情景S4:闽东北区鼓励推广坡度工业,将坡度25°的裸地类型转化为建设用地类型;闽西北区要加大森林资源保护,增加耕地后备资源,将坡度25°的裸地、草地类型转化为林地类型,将坡度15°的裸地类型转化为耕地类型。将以上各个情景分别代入校准后的SWAT模型,模拟闽江流域2008—2018年间不同土地利用情景下的水文径流过程,得到流域径流量,并结合年径流变差系数 式中: 式中:Cv为年径流变差系数;3 结果与分析3.1 模型参数率定与验证SWAT模型参数众多,需要在率定与验证过程中对其进行筛选。根据模型运算原理及各参数意义由表2可知,主要影响闽江流域径流模拟的参数有CN2、CANMX、ESCO、CH_K2、ALPHA_BF等。其中第1敏感参数CN2综合反映了研究区土壤渗透性、前期土壤水分条件以及土地利用类型等特征,是影响流域径流模拟的主要参数。降水过程中,植被冠层可减小雨滴的侵蚀能力并截留部分降雨,CANMX是植被冠层完全成长时的最大截留量,对下渗、地表径流及蒸散发都具有显著影响。ESCO主要作用于土壤蒸发过程,其值越小,模型可以从下层获得更多的蒸发需水量。CH_K2代表主河道冲积物的有效渗透系数,对于持续得到地下水补给的常年性河流,其有效渗透系数为0。基流图2为闽江流域各水文站点逐月径流模拟值与实测值的比较结果,可见模型模拟的逐月径流趋势与实测径流波动趋势大体一致,峰值拟合表现良好。目前在SWAT模型模拟精度的评价标准方面,有研究3.2 情景模拟结果分析根据本文设置的5种土地利用情景模式,通过率定好的SWAT模型,模拟闽江流域不同土地利用情景下的2008-2018年平均径流量,并与基准情景模拟结果进行对比分析(表3)。由表3可知,闽江流域不同土地利用情景下,多年平均径流量与年径流变差系数Cv值变化趋势一致,由大到小依次为:S2(建设开发情景)、S1(耕地储备情景)、S0(基准情景)、S4(分区调控情景)、S3(植被恢复情景)。植被恢复情景下,流域年均径