金沙江中上游流域未来径流模拟研究_杜智毅.pdf

第 卷 第 期 年 月水 资 源 与 水 工 程 学 报 ，收稿日期：；修回日期：基金项目：四川省教育厅农村水安全工程研究中心项目（）作者简介：杜智毅（），男，四川眉山人，硕士研究生，主要从事农业水土环境研究。通讯作者：倪福全（），男，四川成都人，博士，教授，硕士生导师，主要从事水安全和农业水土环境研究。：金沙江中上游流域未来径流模拟研究杜智毅，倪福全，邓 玉，周家豪（四川农业大学 水利水电学院，四川 雅安）摘 要：为了探究金沙江中上游流域未来径流变化趋势，为流域防洪规划提供依据，基于 水文模型，选用 数据集建立未来时段的全球气候模式，从时间和空间尺度解析研究区 年径流变化趋势。结果表明：流域 年降水量和平均气温均高于基准期，并且呈现上升趋势，其中流域南部降水量增幅较大，流域北部气温增幅较大。在 、种气候情景下，年年径流量均呈现增大趋势，变化率分别为 、。相较于基准期，未来春季和秋季径流量呈现减少趋势，夏季和冬季径流量呈现增加趋势，冬季径流量增幅达到了。流域产流量呈现从西北到东南依次增加的特点，相较于基准期，流域南部产流量均呈现增加趋势。未来径流量呈现增加趋势，冬季径流量增幅较大，可能会发生冬汛等极端水文事件，流域南部受洪水威胁的可能性进一步增大。关键词：未来径流模拟；产流量空间分布；模型；金沙江中上游流域中图分类号：；文献标识码：文章编号：（），（，）：，；，；，：；研究背景全球气候变化将改变流域的水文过程，直接或间接地影响流域的产流量和时空分布规律。随着全球变暖趋势的加剧，世纪 年代以来，极端水文事件的强度和频率发生了明显变化。洪水、干旱等极端水文事件往往会造成巨大损失，因此未来径流模拟研究引起了中外学者的高度关注，成为水文与水资源领域的研究焦点。（）模型由于其效率高、精度高、模拟周期较长和易于操作，已被广泛应用于不同流域的水资源评价中，在研究变化环境下的流域径流过程及响应规律、变化趋势以及归因分析等方面均有很好的应用。等使用 模型对洛博流域进行了水量分析，发现 模型能够很好地适用于径流模拟；贾何佳等利用 模型对黄河源区未来径流进行了模拟；马月运用 模型研究了金沙江流域气候变化对径流的影响；陈华等对考虑在水库调蓄下的金沙江径流进行了模拟研究。近年来，模型在金沙江流域得到了一定应用，但仍缺乏对金沙江流域的未来径流变化和未来产流量空间分布的研究。金沙江地处长江上游，流域水能资源丰富，年径流变化幅度较大，极端水文事件频发。该流域上游通天河段是南水北调西线工程的引水水源地，中游多发局域性暴雨，洪水灾害等极端事件时有发生。曾小凡等研究发现，从 世纪 年代以来，金沙江中上游流域平均气温上升显著。第六次联合国政府间气候变化专门委员会（，）报告指出，水安全风险随着全球升温水平的提高而增加，干旱和洪水等极端水文事件发生的概率也会增大。金沙江中上游流域的气温持续升高可能会导致洪水等极端水文事件的发生频率增加。因此，本文建立金沙江中上游流域 模型，结合未来气候模式对金沙江中上游流域未来径流进行模拟研究，评估年尺度和季节尺度下的径流量变化趋势，探讨研究区未来产流量的时空变化特征，以期为金沙江中上游流域防洪规划提供科学依据。数据来源与研究方法 研究区概况金沙江中上游流域（，）跨青海省、四川省、云南省和西藏自治区共 个省份（自治区），总面积 。金沙江全长约 ，自西北向东南穿过青藏高原、横断山脉、云贵高原、四川盆地西南部，地势上西北高、东南低，地形复杂多变，地貌差异大。其中金沙江在云南省丽江市石鼓镇以上为上游、石鼓镇至四川省攀枝花市雅砻江口为中游。金沙江中上游流域概况见图。图 金沙江中上游流域概况金沙江全流域多年平均气温为 ；降水量时空分布不均，上游和中游的多年平均降水量分别为 和 ，降水量年内分布不均匀，月份的降水量占全年的。金沙江中游水能资源丰富，是整个长江上游水电开发体系中不可或缺的一环。此外，金沙江中上游流域是重要的旅游区，旅游资源十分丰富，在当地社会经济的发展中占有重要地位。数据来源及用途 模型数据包括数字高程数据（，）、土地利用数据、土壤数据、气象数据和水文数据等。气象数据为流域 年的日降水量、最高气温、最低气温、湿度、风速和太阳辐射数据；土地利用分为耕地、林地、草地、水域、城镇用地和未利用土地 类；土壤数据来源于全球土 壤 数 据 库（，），可直接用于 模型土壤数据构建；水文数据为 年攀枝花水文站月径流观测数据，构建金沙江中上游流域 模型所需数据见表。水 资 源 与 水 工 程 学 报 年表 模型数据一览表数据名称数据来源分辨率及处理方式数据用途数字高程数据地理空间数据云，裁剪、投影变换生成流域河网、划分子流域土地利用数据地理空间数据云，遥感解译、重分类、投影变换 划分土壤类型数据 土壤数据库，裁剪、投影变换 划分，构建土壤数据库气象数据研究区 个气象站日值，建立天气发生器、按一定格式应用 软件处理构建气象数据库径流数据攀枝花水文站月值率定及验证 未来气候模式未来气候模式已经成为研究未来径流变化的重要工具之一。世界气候研究计划（，）组织的（）为国际耦合模式的评估和后续发展提供了重要的平台。目前有 个大气环流模式，被广泛应用于气候变化研究领域。本文选择在中国表现较好的大气环流模式 下的 、共 种不同排放情景进行径流模拟研究，其中 为高排放情景，为中等稳定化情景，为低排放情景，即分别代表在 世纪末辐射强迫达到、和 。水文模型 模型是美国农业部（，）开发的流域尺度的分布式水文物理模型，径流模拟是 模型最基本、最重要的功能，运用 模型进行径流模拟研究是 研究的焦点。模型中有众多参数，本文结合金沙江中上游流域的区域特点选取了 个参数对构建好的 模型模拟结果进行了校准与验证，并选取决定系数、纳什系数、误差百分数、均方根误差与标准误差的比率 用于评价模型的模拟精度。决定系数 表示的是模型模拟值与径流实测值之间的贴合程度，越接近于，说明模拟值与实测值越贴合，模型模拟效果越好；是指模型模拟值与径流实测值之间偏离程度，越接近于，说明两者偏差越小；是指模型模拟值与径流实测值之间的变化趋势，说明模型模拟值偏大，反之，模型模拟值偏小；是指均方根误差与实测径流值标准差的比值，绝对值越小，则模型模拟效果越好。一般情况下，当 ，时，可认为模拟结果与实测数据拟合情况较好，说明模型可以基本准确地反映该流域的水文循环过程。参数敏感性分析参数与模型运行的不确定性与模型的运行效率密切相关，但是 模型的参数众多，因此需要进行敏感性分析找出敏感性较强的参数，本研究使用 （）算法对参数进行敏感性分析。算法结合了 抽样法和 敏感性分析法，该算法以 抽样思想将整个参数空间进行分层，在各分层中依次进行随机抽样。现将参数空间划分为层，每层包含 个参数，即生成了 个 抽样参数组。在此基础上根据 敏感性方法思路，将 个抽样参数组中的 个参数进行微小调整改变，即改变参数，计算目标函数随每次微小改变的情况。本文以 和 检验值参数作为敏感程度的评判标准，值越接近于 且绝对值越大，则说明参数越敏感。结果与分析 模型的适用性分析采用攀枝花水文站月径流观测值对构建的 模型进行校准验证，选取 年为模型模拟的率定期，年为模型模拟的验证期，率定期和验证期流量模拟值与观测值对比见图，模型率定及验证结果各评价指标值见表。由图 和表 可见，模型参数经校准和验证后，流域出口攀枝花水文观测站月径流模拟效果较好，、和 指标均符合模型评价标准，表明 模型在金沙江中上游流域径流模拟中具有较好的适用性。表 模型率定及验证结果各评价指标值水文站阶段年份 攀枝花率定期 验证期 第 期 杜智毅，等：金沙江中上游流域未来径流模拟研究 参数敏感性分析结果本文结合金沙江中上游流域的特点，选取 个参数对模型进行校准验证，采用 算法对所选取的 个模型参数进行敏感性分析，个参数的名称、定义和敏感性排名以及敏感性分析结果见表 及图。图 模型率定期和验证期流量模拟值与观测值对比表 模型参数敏感性排名排名参数名称参数定义排名参数名称参数定义降雪基准温度植物吸收补偿因子降雪最高融化速率湿润土壤反照率最大冠层截留量基流储存因子饱和水力传导系数深蓄水层蒸发系数平均坡长地下水滞后系数土层表层至底层深度年中降雪最低速率土壤蒸发系数浅层地下水蒸发系数融雪基准温度平均坡度渠道曼宁系数 径流曲线数土壤湿容重温度递减率图 模型参数敏感性分析结果 水 资 源 与 水 工 程 学 报 年 由表 和图 可见，对金沙江中上游流域径流模拟结果影响最大的参数是决定降雪基准温度的，越大，代表流域降雪的基准温度越高；其次是代表年中降雪最高融化速率的 和代表流域植被最大冠层截留量的。另外，饱和水力传导系数、平均坡长 和土层表层至底层深度 也均对金沙江中上游流域径流模拟结果较为敏感。未来气候变化分析 未 来 降 水 情 况 分 析 利 用 、数据，将 年作为基准期，分析 年金沙江中上游流域的年降水量变化趋势，结果见图。由图 可知，相对于基准期，种气候情景下流域年降水量均有增大的趋势，由多年平均降水量对比来看，、和 未来气候情景下降水量较基准期分别增加了 、和。在 种不同情景下，年年降水量随时间整体均呈现增大的趋势，年变化率分别为、和 。图 基准期及不同气候变化情景下 年流域年降水量变化 为了研究空间尺度下的降水变化规律，将金沙江中上游流域划分为 个子流域，计算各个子流域的多年平均降水量。图 为金沙江中上游流域各个子流域在不同情景下较基准期的降水量变化情况。由图 可以看出，在、和 未来气候情景下，除 号子流域外，其他各个子流域的降水量均呈现增加趋势；不同气候情景下的降水量变化量空间分布无显著差异，研究区南部降水量增幅最大。图 不同气候变化情景下各个子流域较基准期的降水量变化量第 期 杜智毅，等：金沙江中上游流域未来径流模拟研究 未来气温情况分析利用未来气候情景数据（、），将 年作为基准期，分析金沙江中上游流域 年的最高、最低气温变化趋势以及平均气温变化量的空间分布，如图 所示。由图、可知，相对于基准期，种气候情景下流域 年多年平均最高气温和平均最低气温均有升高的趋势。、和 情景下的多年平均最高气温较基准期分别升高、和 。、和 情景下的多年平均最低气温较基准期分别升高、和 。年 种气候情景下流域各年最高气温和最低气温均呈波动式上升趋势，最高气温变化率分别为 、；最低气 温 变 化 率 分 别 为 、。其中 情景下的年最高气温和年最低气温上升率最大，情景下的年最高气温和年最低气温上升率最小。由图 可见，在 种情景下，各子流域多年平均气温几乎均呈现上升趋势，其中研究区北部升温幅度最大。图 基准期及不同气候情景下流域 年最高气温变化图 基准期及不同气候情景下流域 年最低气温变化 水 资 源 与 水 工 程 学 报 年图 不同气候变化情景下各个子流域较基准期的多年平均气温变化量 结合以上分析，与基准期相比，年金沙江中上游流域降水量和气温均呈现增加和上升趋势。张英娟等研究发现中国西部地区气温、降水量及湿度均呈现显著的增加趋势，与本文的结果一致。气温的显著上升可能会导致金沙江上游冰川、积雪的加速消融，从而影响金沙江中上游流域的产流机制。未来径流量模拟利用构建好的 模型及全球气候模式数据，对金沙江中上游流域 年年径流量变化趋势进行模拟预测，基准期及不同气候变化情景下流域年径流量变化见图。图 基准期及不同气候变化情景下 年流域年径流量变化 由图 可知，种情景下，未来多年平均径流量分别为 、和 ，分别较基准期增加了 、和减少了 。不同气候情景下未来年径流量均呈现增大趋势，情景下年径流量变化最明显，年径流量增长率为 ；情景下年径流量增长趋势最小，增长率为 。金沙江中游洪水遭遇明显的年份为 和 年，年径流量均高达 以上（异常值），引发了严重的洪涝灾害。情景下，年年径流量超过 ，、年的年径流量接近 ，因此、及 年很可能再次出现大的洪水，引发洪涝灾害；情景下，年年径流量与 年相近，因此在 情景下 年流域可能发生类似于 年的洪涝灾害；在 情景下，只有、年的年径流量比较接近于，其他年份均未达到 年的径流量值。