基于SWAT模型的那棱格勒河流域径流模拟及水量平衡分析.pdf
《基于SWAT模型的那棱格勒河流域径流模拟及水量平衡分析.pdf》由会员分享,可在线阅读,更多相关《基于SWAT模型的那棱格勒河流域径流模拟及水量平衡分析.pdf(9页珍藏版)》请在咨信网上搜索。
1、刘兆峰,王建萍,韩进军,等.基于SWAT模型的那棱格勒河流域径流模拟及水量平衡分析 J.盐湖研究,2024,32(1):48-56.Liu Z F,Wang J P,Han J J,et al.Runoff simulation and water balance analysis in the Nalenggele River Basin using the SWAT Model J.Journal ofSalt Lake Research,2024,32(1):48-56.DOI:10.12119/j.yhyj.202401006基于SWAT模型的那棱格勒河流域径流模拟及水量平衡分析刘兆峰
2、1,2,3,王建萍1,2*,韩进军1,2,陈元军4(1.中国科学院青海盐湖研究所,中国科学院盐湖资源综合高效利用重点实验室,青海 西宁810008;2.青海省盐湖地质与环境重点实验室,青海 西宁810008;3.中国科学院大学,北京100049;4.青海中信国安锂业发展有限公司,青海 格尔木816000)摘要:那棱格勒河是柴达木盆地内流域面积和流量最大的内陆河流。本文在那棱格勒河上游构建了SWAT分布式水文模型,采用SUFI-2算法进行参数敏感度分析、率定和验证,对19792018年的月值径流量进行模拟分析。结果表明,敏感度最强的前3个参数依次为融雪基温、降雪气温及气温直减率,在率定期和验证期
3、决定系数R和Nash-Sutcliffe效率系数分别高于0.81和0.65。水量平衡分析结果显示19792018年流域的蒸发量、地表径流量、地下水补给量、地下侧流量分别占降水的92.30%、17.76%、16.18%和5.84%。表明流域内径流量主要来自冰雪融水,各水量平衡分量受温度影响较大。降水量变化趋势与地表径流量变化趋势较一致,说明降水转化为径流的时间较短。本研究对那棱格勒河流域水文研究及水资源调控具有参考意义。关键词:径流模拟;SWAT模型;径流成分;那棱格勒河流域中图分类号:P333.1文献标识码:A文章编号:1008-858X(2024)01-0048-09深入研究径流量模拟有利于
4、了解流域内水资源分布状况,规划和管理流域水资源,为水资源时空合理调配提供参考依据1。研究径流的发生发展过程对于水质和水量的定量预报至关重要,并为流域内防洪防旱做出预警。水文模拟是水文科学研究重要的方法之一,水文模型作为研究径流过程的有效工具,在水文领域的作用日益凸显。目前全世界已开发许多基于集总式、半分布式、分布式水文模型,其在不同的流域取得了较好的应用效果。如邵成国等2在乌鲁木齐河流域运用有融雪结构的新安江模型进行了径流模拟;高红凯等3对长江源区有冰川覆盖的冬克玛底河流域进行日径流模拟试验,得到流域冰川物质平衡变化;吕栋等4在黄泥庄流域应用了 VIC 模型,探讨了该模型在中小尺度流域上的适用
5、性并对径流进行模拟;李抗彬等5通过对 TOP-MODEL 模型的蒸发产流模块以及汇流模块进行改进,使其适用于半湿润地区降雨径流模拟;Gomez-landesa 等6在西班牙的庇里牛斯山地区应用 SRM 模型模拟了流域雪水当量的变化,预测了多个盆地的每日流量。与这些水文模型相比,SWAT模型能够与 GIS 软件集成,在国内外应用广泛,能够利用遥感影像提供的空间信息模拟多种不同的水文物理化学过程7。SWAT 模型可以很好的反应不同水文要素和下垫面分布的不均匀性8,在流域水文平衡、长期的地表径流模拟等方面得到了广泛的应用9。国内外对该模型展开了广泛的研究,孟现勇等研究了 CMADS 数据集在东亚地区
6、的适用性10,并在黑河流域所属高寒山区进行应用验证11;李帅等12在宁夏清水河流域应用 SWAT 模型发现气候变化主要通过降雨量和温度变化调控径流量;张菡等13在海拔高差悬殊的岷江上游地区应用 SWAT 模型模拟了日径流变化并基于模型计算气象干旱指收稿日期:2023-08-22;修回日期:2023-09-18基金项目:中国科学院重点部署项目(ZDRW-ZS-2020-3);第二次青藏高原综合科学考察研究(2019QZKK0805-02);青海省自然科学基金团队项目(2022-ZJ-903)作者简介:刘兆峰(1996-),男,硕士研究生,主要研究方向为水文地质方向。Email:。通讯作者:王建萍
7、(1972-),女,博士,研究员,主要从事自然资源管理与评价。Email:。第 32 卷第 1 期2 0 2 4年 2 月JOURNAL OF SALT LAKE RESEARCH盐湖研究Vol.32 No.1Feb.2024刘兆峰,等:基于SWAT模型的那棱格勒河流域径流模拟及水量平衡分析第1期数;Kassem 等在伊拉克库尔德斯坦地区使用 SWAT和人工神经网络 ANN 耦合模型进行日径流量的预测,结果显示其精度均优于单一模型14;Li 等使用了卫星遥感降雨产品作为降水量输入 SWAT 模型进行径流量模拟,并对比了气象站实测数据的模拟结果,为无资料地区水文模拟提供了借鉴15。北方高纬度寒区
8、由于特殊的气候和下垫面条件,径流过程往往存在较大的非线性16。受到地理条件的限制,很多高山地区都缺少观测资料,目前现有的无资料地区的水文模拟方法主要是利用水文系统响应、数理统计以及水文模型等研究方法。水文系统响应和数理统计的方法缺乏水文物理基础,很难获取水文时空变异特征以及反映不同尺度的气候、下垫面因素对水文过程控制作用的机理17。为了能更好的在严重稀缺资料地区开展水文模拟,选取 CFSR 数据集(Climate Forecast System Reanalysic)作为模型输入18,订正后的降雨数据能较好地实现模拟。那棱格勒河上游流域属于典型的冰雪融水补给山区,且资料较为匮乏,水文方面的研究
9、缺少,王涛在那棱格勒河流域建立了 TOPKAPI 模型进行该流域 1993 年月流量水文模拟19,由于气象水文等资料匮乏,不能对模型进行有效调试。本文在结合多种研究方法的前提下,建立适合该流域产汇流的SWAT 模型,为验证模型适用性,优化模拟效果,尤其在数据资料缺失或者精度不够时,采用顺序不确定拟 合 算 法(Sequential Uncertainty Fitting Version2,SUFI-2)进行参数率定20,使用拉丁超立方体随机采样法(Latin hypercube sampling),通过较少迭代次数抽样来准确重建输入分布,在多次模拟后得到一组最优模拟数值21,模拟结果是由一组模
10、型参数决定而非单一参数,能通过不断模拟减少参数不确定性的范围,使得模拟结果不断接近实测值。将 SWAT模型应用于那棱格勒河流域,在分析流域水文要素的基础上,对该流域的径流进行模拟,研究那棱格勒河流域山区水量要素的变化规律,为今后水资源的开发利用及调配提供参考依据。1研究区域概况那陵格勒河(9240E9315E,3640N3715N)位于青藏高原柴达木盆地西南部,流域发源于昆仑山脉的布喀达坂峰,该河是山前冲洪积平原上的一条常年性河流,全长 435 km,是柴达木盆地流域面积和水量最大的内陆河22(图 1)。那棱格勒河源区位于东昆仑阿尔金断裂碰撞交汇地带,补给水源主要是高山冰雪融水。高海拔的山间盆
11、地发育有很多盐湖,如勒斜武旦湖、西金乌兰湖、太阳湖、可可西里湖、库塞湖等,这些盐湖的矿产资源是支撑当地经济发展的命脉23。该流域气候具有多大风、蒸发强烈、气温低、降水极少、冬长夏短、昼夜温差大等特点,属于典型高原内陆高寒干旱气候。上游地段河谷狭窄,水流清澈,接受地下水、降水和夏季冰川融水的补给,全年有水;中游为季节性河流,流出山口后迅速渗入地下,地表水补给地下水;下游又重新溢出地表,形成集泉河,其中最大的有台吉乃尔河与乌图美仁河,最终注入台吉乃尔湖24。河流径流量呈现明显的季节性变化规律。河流径流量最枯时期一般为每年 10 月至翌年 3 月,河流仅接受山区内地下水补给,径流量小且较为稳定;4
12、月份后,随着气温升高,冰雪融化,河流接受冰雪融水及山区地下水补给,径流量逐渐增大;59 月份为河流径流量最丰时期,也是降水最集中的雨季,降水直接补给河水,河流径流量显著增大,7 月中旬至 8 月末达到最大;后续气温降低,降水减少,河流径流量也不断减少,10月份以后进入河流枯水期25。2数据来源及SWAT模型构建2.1数据来源及处理在那棱格勒河流域建立的分布式水文模型SWAT 模型需要输入 DEM、土壤数据、土地利用数据以及气象数据,模型所需数据和用途详见表 1。通过GIS 和数理统计等手段将数据处理并添加到模型所需的属性及空间数据库中以便使用。910 0 东920 0 东930 0 东370
13、0 北360 0 北910 0 东920 0 东930 0 东370 0 北360 0 北040 km高程59805322496045894228385734953134图例水文站河流边界N图1研究区地理位置图Fig.1Location of the study area49盐湖研究第 32卷研究区 DEM 数据使用地理空间数据云分辨率为 30 m30 m 的 GDEM 格式的数字高程数据,反映地形、河网和子流域等水文相关信息,基于 GIS 平台,对该数据进行拼接、剪裁、投影、提取等操作,制作模型所需的数据资料。土地利用数据来源为清华大学土地利用监测数据,分辨率为 10m,按照 SWAT模型里
14、面的土地利用分类标准并将其重新进行分类,划分为草地、水域、裸地 3 类(见图 2),其中草地占 43.1%,水域占 3.2%,裸地占 53.7%。土壤数据采用全球土壤数据库(HWSD)中的中国土壤数据集,空间分布率为 1 km1 km。在SWAT 模型中构建土壤数据库需要的土壤参数如土壤湿密度、有效持水量、饱和导水率等由 SPAW 软件计算得到26。研究区经过重新分类后的土壤类型有15 种(见图 2),得到该区域的土壤类型主要为薄层土、雏形土、砂性土、黑土和潜育土。2.2子流域划分及水文单元确定SWAT 模型可以在水文响应单元的尺度上进行地表径流、入渗、地下水、融雪径流、蒸散发等多种过程的模拟
15、27,在输入研究区域水文气象资料、土地利用及土壤数据、数字高程数据之后,利用 SWAT 模型流域划分模块将那棱格勒河流域上游流域划分为 45个子流域(见图 3)。水文响应单元(HRU)是模型汇流演算及计算水量平衡的最小单元,是指下垫面单一均匀的区域,可以用来反映不同土地利用类型、土壤类型和坡度组合的空间异质性。本次研究在划分水文单元时,通过设置阈值的方式,将小于阈值的部分按照比重分配到大于该阈值的土地利用状况、土壤特性和地形坡度数据当中,其中土地利用阈值设为 10%,土壤类型阈值设为 10%,坡度阈值设为10%28,将研究流域划分为 613 个水文响应单元,计算各个子流域面积、坡长、坡度变化等
16、参数。地表径流的计算在 SWAT 模型中采用美国农业部水土保持局(Soil Conservation Service)研制的SCS 径流曲线方法计算29,是用于估算无资料地区径流量的经验模型。SCS 法计算过程简单,需要参910 0 东370 0 北920 0 东930 0 东360 0 北370 0 北360 0 北910 0 东920 0 东930 0 东910 0 东920 0 东930 0 东370 0 北360 0 北370 0 北360 0 北N03060 km03060 kmN图例草地水体雪地裸地图例简育砂性土过渡性红砂土黑色石灰薄层士石灰性疏松岩性士石灰性黑土潜育黑土松软潜育土
17、永冻薄层土永冻雏形土永冻薄层土永冻薄层土永冻层土黑色石灰薄层土水体图2土地利用数据分布图,土壤属性数据分布图Fig.2Land use of the study area,soil type of the study area表1模型数据类型及来源Table 1Model data type and source数据类型数字高程DEM土壤数据土地利用数据气象数据水文数据分辨率30 m1 000 m10 m日值月值数据来源地理空间数据云HWSD全球土壤数据库清华大学土地监测数据CFSR数据集H19测站数据用途生成河网、划分流域构建土壤数据库制作HRU构建属性数据库径流模拟子流域出口水文站子流域流
18、域边界Watershed1235467910131617191815353637301211143132292026272822242921233443404244413345图3那棱格勒河流域子流域划分Fig.3Nalenggele river basin sub-watershed division50刘兆峰,等:基于SWAT模型的那棱格勒河流域径流模拟及水量平衡分析第1期数较少,资料易获取,对观测数据的要求不严格,并有效地考虑了不同流域下垫面的特点。CN 值(Curve number)是 SCS 法的主要参数,可将前期土壤湿润程度、坡度、土壤类型和土地利用现状等综合在一起后用量的指标来反
19、映下垫面对产汇流过程的一个影响30。计算原理为S=(25 400/CN)-254In=0.2SQsurf=()Rday-In2/()Rday-In+S,Rday 0.2S;Qsurf=0,Rday0.7,NSE0.5,就认为模拟效果较好36,从模型的适应性评价指标可知,P 因子和 R 因子在率定期和验证期也均符合要求,证明 SWAT 模型能较好地模拟那棱格勒河流域的水文径流过程。3.2径流恢复重建及水量平衡分析流域的水循环过程对区域水量水质有着重要的影响,对于区域生态环境的可持续发展起到重要作用37,流域的水量平衡中的各个分量对于资料缺少地区来说观测非常困难,通过借助水文模型的手段,可以探求各
20、分量的分配。流域的水量平衡是指在任意时间段内流域获得水量与支出水量之间的差额为该时间段内蓄水的变化量,通过 SWAT-CUP 对那棱表3月径流模拟结果精度评价Table 3Evaluation of monthly runoff simulation results时间段率定期验证期评价因子P-因子0.570.47R-因子0.500.42R0.830.71NSE0.810.652011/012011/032011/052011/072011/092011/112012/012012/112012/032012/052012/072012/09时间模拟径流/(m3/s)2019/012019/0
21、32019/052019/072019/092019/112020/012020/112020/032020/052020/072020/09时间实测径流模拟径流实测径流模拟径流水量250200150100500水量25020015010050018016014012010080604020200175150125100755025y=0.87x+5.95,R2=0.83模拟径流/(m3/s)实测径流/(m3/s)y=0.79x+14.53,R2=0.710255075100125150175实测径流/(m3/s)255075100125150 175200a、c分别为率定期、验证月径流量实测
22、值和模拟值;b、d分别为率定期、验证期实测径流拟合效果图图4那棱格勒河流域SWAT模型径流模拟情况Fig.4Rate regular runoff simulation of SWAT model in Nalenggele River52刘兆峰,等:基于SWAT模型的那棱格勒河流域径流模拟及水量平衡分析第1期格勒河上游流域模型参数校准和率定后,将验证识别后的参数带入到 SWAT 模型种再次模拟,借此分析那棱格勒河上游流域的水量平衡38,其水量平衡关系式如下:S=PREC+SNOMELT-SURQ-LATQ-PERC-ET(3)式中 S 为土壤储存水变化量(mm),PREC 为流域降水量(mm
23、),SNOMELT 为融雪量,SURQ 为地表径流量(mm),PERC 为土壤对地下水补给量(mm),LATQ为地下侧流量(mm),ET 为实际蒸散发量(mm)。通过模拟得到那棱格勒河流域 19792018 年年平均水量平衡分析结果,结合表 4 得知:各分量所占比例大小依次为蒸散发量、融雪量、地表径流量、地下水补给量、地下侧流量,而那棱格勒河流域上游山区可以忽略人为影响,可认为土壤平均含水量保持不变,土壤水变化量为 0;那棱格勒河流域的上游区域中,实际的蒸散发是流域水量中主要的输出项,占年均降水量的 92.30%;融雪量占降水量的 78.96%,可知流域的径流组成中冰雪融水所占比例较大,同湿润
24、区有很明显的差异39,地表径流略大于地下径流。表4那棱格勒河流域年平均水量平衡表Table 4 Average annual water balance table for the Nalenggele river basin197920181979198019811982198319841985198619871988198919901991199219931994199519961997199819992000200120022003200420052006200720082009201020112012201320142015201620172018平均值比例/%降水量409.88405
25、.83420.18367.65395.08414.00367.29376.71361.82368.06383.71358.65411.65387.49412.77436.86423.56389.18438.51407.59438.50453.19413.17443.82438.81395.78421.29466.69396.06429.77467.94341.32410.34381.51377.65426.12361.03360.85364.19382.30402.67-融雪量317.84330.56337.45275.36306.75323.15278.13285.75283.58277.9
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 SWAT 模型 棱格勒 河流 径流 模拟 水量 平衡 分析
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。