未来30年金华江流域的径流与枯水指数计算,水文学论文

未来30年金华江流域的径流与枯水指数计算,水文学论文气候变化日渐成为世界关注的问题。气候变化背景下，各个水文要素，如降水、蒸发、径流乃至整个水文循环都发生了不同趋势的变化。很多专家就气候变化对水文循环的影响做出了预测和评估，这些评估牵涉到降水、温度、径流和水资源等各个方面[1,2].枯水水文是水文学的重要组成部分，与人类社会生活、经济活动以及自然生态密切相关[3].气候变化的加剧和工农业用水的急剧增加，使得枯水期水资源的问题影响到供水、农业和生态安全。但与传统水文相比拟，枯水水文研究还处在相对被忽视的阶段。我们国家的枯水水文从50年代开场发展，到80年代才得到大面积开展。近年来，一些学者对枯水开展了一系列研究，并获得了一些成果。然而，在已经知道的研究中，甚少有人结合统计降尺度模型和水文模型来分析将来枯水径流的变化趋势。本文使用统计降尺度中的SDSM模型来模拟钱塘江支流-金华江流域五个主要站点将来2020s〔2018-2040年〕的日降水和金华站的日最高和最低温度，并结合萨克拉门托水文模型，计算了将来30年〔2018-2040年〕金华江流域的径流。在这里基础上，分析其枯水指数7Q10和30Q10的变化趋势。1研究方式方法1.1数据来源金华江是钱塘江的支流之一，主源是东阳江北江，起源于磐安县龙乌尖，汇于衢江，全长195km,比降3.1。流域面积5990平方公里。本文研究对象为金华江金华站点以上的流域范围。研究所用的原始数据为1961-1990年的逐日实测气象和水文资料，包括降水、温度和径流等。各气象站点信息如表1所示。【1】本文采用的GCM是英国哈德利中心的HadCM3,采用的温室气体排放情景是IPCCAR4中的SRES情景中的A2和B2情景。在情景A2中，人口持续增长，人均经济增长和技术变化有明显地方性，全球化不明显；B2情景是人口增长低于A2,经济中等发展，技术更多样化，侧重于从局地解决经济、社会和环境可持续性问题的情景[4].1.2SDSM模型本文中，预测期2020s的降水和最高最低温度使用SDSM模型[5]进行计算。华而不实最高和最低温度用于计算流域潜蒸发。SDSM模型是统计降尺度的一种。统计降尺度的统计关系是建立于大尺度气候状况和下垫面各种特征之间的，依靠多年观测资料与地域特征建立。统计关系需要有独立的观测资料检验，然后把这种关系应用于GCMs输出的大尺度气候信息，来预估区域将来的气候变化情景[6].统计降尺度的目的是建立大尺度气候预报因子和预报量的函数关系，进而在三个假设条件下，即大尺度气候场和区域气候要素之间有显着的统计关系，大尺度气候场能很好的被GCMs模拟，建立的统计关系不由于气候情景变化而失效[7],实现大尺度气候信息的降尺度。SDSM方式方法是统计降尺度的一种详细方式方法，已经在我们国家多个地区应用，它结合了天气发生器和多元回归来进行降尺度计算，其核心部分是预报因子的挑选，具有易于操作，计算快速，结果可靠等优点。1.3水文模型萨克拉门托模型是一个概念性水文模型，已应用于多个流域，在中国南方地区也有应用。在模型中，流域面积被分为三个部分，分别是透水、不透水和可变不透水三个部分。模型要求的输入数据为实测的日降水、潜蒸发和径流以及将来的降水和潜蒸发。输出的结果是模拟的径流。模拟后的径流分四个部分，分别是来自透水面积和可变透水面积的直接径流，上层土壤水容量知足后过剩降水构成地表径流，来自上层自由水的壤中流，以及来自下层自由水的地下径流。本文的输入数据为金华江流域基准期的降水，潜蒸发和径流的日序列，而将来的降水和温度则由SDSM模拟计算而得。水文模型率定和验证的优化方式方法采用遗传算法，并以Nash效率系数系数作为第一目的函数，以Grayson〔1996〕定义的基流函数[8]作为辅助目的函数，以提高对枯水的模拟能力。潜蒸发量通过HARGREAVES方式方法计算，公式如下[9]:【2】2计算结果2.1SDSM模拟检验为了验证SDSM模型的可靠性，将基准期的实测数据与模拟数据相比拟。美国国家环境预报中心〔NCEP〕的全球再分析资料通常在气候模拟研究中作为驱动来检验模拟的结果。因而，本文选用NCEP再分析数据作为SDSM的驱动，来生成模拟数据，检验其可靠性。由于降雨站点较多，这里故选取八达站为例。温度的比拟选取金华站的日最高温度为例。比拟的目的分别为日降雨量最大值，日降雨量均值，日最高温度的最大值以及日最高温度的最小值。比拟结果如此图1所示。分析图表可知，日降雨量的平均值、日最高温度的最大值和最小值实测值与模拟值拟合程度较高，而日降雨量的最大值在夏季的拟合程度较差。这表示清楚，SDSM能够较好模拟温度和降雨量的时间序列，但在夏季，由于有多个影响因素的存在，日降雨量的最大值模拟的准确性会遭到十扰。2.2将来降雨量预测以八达站为例，选取p5vCSOOhPa高度径向风速〕，p8_zC800hPa高度涡度〕，r500CSOOhPa高度相对湿度〕，r850C800hPa高度相对湿度〕，rhum〔表而相对湿度〕作为预报因了。预报因了与刘吕明[7]关于SDSM模型预报因了的总结相吻合。实测数据和两种不同情景下的模拟数据在不同重新期下的值如此图2所示。结果表示清楚，除了夏季的五月、六月和九月，两种情景下的将来平均降雨量都有所增加，华而不实七月份的增LLJ隔度最大，而情景A2和情景B2的区别较小，目不具备明显的规律。降雨极值方而，则出现了较大变化。除了六月、八月和九月外，模拟值比预测值都有较大的增长。这讲明在气候变化影响下，金华江流域在将来30年在降水方面更容易出现极端天气。2.3潜蒸发计算潜蒸发的计算由日最高温度和日最低温度计算而得，华而不实将来30年的日最高温度和日最低温度由SDSM模型模拟而得。以日最高温度为例，选取malpC海平面气压〕，p500CSOiOhpa的位势高度〕，shum〔地表比湿〕，temp〔地表平均温度〕作为预报因子。模拟温度与实测温度以及潜蒸发比拟如此图3所示。分析图表能够发现，预测期与基准期相比拟，最大位有明显增大的趋势，情景A2比情景B2增加的幅度更大；均值增长幅度较小，两个情景没有明显区别；极小值方面，除了七月、八月和十一月，都呈现出减小的趋势；从一月到九月的潜蒸发均值也相应呈现出增长的趋势。这表示清楚在气候变化影响下，将来二十年的金华江流域将出现更多的极端天气，高温将更高层次，低温则更低。这和诸多学者关于气候变化的预测相符。2.4径流模拟结果预测期的径流用萨克拉门托模型模拟而得。模型使用遗传算法按日率定，率定结果表示清楚，枯水的模拟值与实际值比拟略微偏小，但模型的Nash效率系数到达了0.84,即总体能较好的模拟将来的径流。模型率定效果见图4.2.5枯水指数计算==为了分析预测期2020s径流的枯水变化趋势，本文采用了两个重要的枯水指标7Q10和30Q1007Qlp是指10年一遇的连续七天径流均值的最小值，30Q10则是10年一遇的连续30天径流均值的最小值。在枯水水文中，7Q10是最重要的枯水指标之一。使用Weibull0]分布分析实测径流序列和模拟径流序列，能够得出重要的枯水指标7Q10III30Q10.在实际应用中，7Q10是设计枯水的重要指标之一，在保卫水质、水资源调度、水生生物和湿地生态的保卫、十旱预测以及水利设施建设参考等各方面都发挥了重要作用}a}.而30Q10相比拟7Q10,应用较少，但在水质十分是饮用水，生态等方面也发挥了重要的作用。基准期的实测枯水指标和预测期2020s的模拟枯水指标见表2和表3.比拟不同重现期的连续平均流量能够发现，将来30年和基准期相比拟，都有不同幅度的增长。考虑到水文模型对枯水的模拟值略微偏小，将来30年实际增大的幅度可能会更大。7天连续平均值中，情景A2增长更为明显。而30天连续平均值中，则B2增长更为明显。能够推断，情景B2更容易出现长时间的极端气候现象。【3】3结论通过SDSM方式方法模拟得出的降雨量和温度表示清楚，该方式方法能有效地模拟将来的降水和温度。总体而言，温度的模拟效果优于降水，平均值的模拟优于极端值的模拟。某些月份的模拟值出现反常现象，造成反常现象的因素还需要进一步分析。在枯水模拟的经过中，出现了模拟值略微偏小的情况。引起这一现象的不确定因素，尚需要进一步研究。枯水指数的分析结果表示清楚，在两种情景下，枯水指数均有所增长。对7d连续指数系列而言，A2情景增长幅度较大，而对30d连续系列指数而言，B2情景增长幅度较大。这讲明在将来的30年中，干旱的极端气候可能会有所缓和，但A2情景比B2情景更容易出现长期的干旱情况。以下为参考文献：[1]夏军，左其亭。国际水文科学研究的新进展[J].地球科学进展，2006,21〔3〕：256-261.XIAJun,ZUOQiting.Advancesininternationalhydrologicalscienceresearch[J].AdvancesinEarthScience,2006,21〔3〕：256-261.〔inChinese〕[2]张徐杰，许月萍，高希超，等。CCSM3形式下汉江流域设计暴雨计算[J].水力发电学报，2020,31〔4〕：49-53.ZHANGXujie,XUYueping,GAOXichao,etal.EstimationifdesignstorminHanriverbasinwithCCSM3model[J].JournalofHydroelectricEngineering,2020,31〔4〕：49-53.〔inChinese〕[3]SmakhtinV.U.Lowflowhydrology:areview[J].JournalofHydrology.240〔2001〕：147186.[4]XuYP,ZhangXujie,TianYe.2020.Impactofclimatechangeon24-hdesignrainfalldepthestimationinQiantangRiverBasin,EastChina.HydrologicalProcesses,26〔26〕：4067-4077.[5]WilbyRL.TheStatisticalDownscalingModel:insightsfromonedecadeofapplication[J].InternationalJournalofClimatology.inPress.[6]范丽军，符淙斌，陈德亮。统计