opmodel模型在枯水期日径流模拟中的应用.docxVIP

opmodel模型在枯水期日径流模拟中的应用 东江流域是珠江流域的主要支流之一，主要流经中国经济发达的亚热带。除了广东省惠州、深圳、东莞等市提供供水外，该流域还承担着向香港的供水任务。香港80%的供水来自东江流域。由于广东省主要受亚热带季风气候的影响，年流量在年内发生了很大变化。年，70-80%的年降水量集中在每年的5月至9月，旱季基本停止。近年来，中国南方干旱地区频繁发生旱旱，旱季遭受严重的供水严重的供水危机，以及沿海地区咸潮的支撑作用，严重影响了社会生活和经济活动。因此，预测江东流域水源中断流量的预测对于水资源管理在很大程度上是必要的。 对于水资源管理和生态环境保护而言,日径流过程的模拟是非常重要的,尤其是枯水期径流过程的模拟对水资源配置、生态环境需水量的研究更为必要.但一般情况下,仅利用降雨径流模型模拟出流域的日径流过程,对枯水期的径流过程很少给予关注,也很少对模拟结果的精度进行专门地评价.针对这一问题,Smakhtin等采用半分布式的降雨径流时变模型(VTI)对南非的五个研究流域进行模拟分析,分析结果认为该模型可以模拟出具有一定精度的枯水期径流过程,并建议采用其它模型对此问题进行分析.在此基础上,黄国如等以东江流域为例,较为深入地分析了新安江流域水文模型模拟流域枯水径流过程的可行性及其精度.上述研究均是将丰水期和枯水期作为一个整体进行考虑,先模拟出流域的整体日径流过程后,再对枯水径流进行分析.但由于枯水期和丰水期的模型参数并不完全相同,如果将两个阶段分开考虑,分别研究其径流模拟过程,应该会更为精确.TOPMODEL(a TOPography based hydrological MODEL)模型为基于地形特征的半分布式流域水文模型,该模型结构相对简单,且充分考虑了地形特征对流域径流模拟精度的影响.该模型的产流方式与蓄满产流类似,但与基于蓄满产流的新安江模型的基本原理和计算方法有着显著的差异.尽管TOPMODEL模型对日径流过程有着较高的模拟精度,但由于枯水径流有其固有的特征,利用该模型能否较好地模拟枯水径流过程值得深入探讨.因此,本文利用TOPMODEL模型对东江流域进行枯水季节的日径流过程模拟,在对径流模拟精度进行评估时,除了采用常规的评价指标外,还采用多种描述枯水径流的特征量来评估该模型的枯水模拟精度. 1 u3000地表径流及产汇流过程 1979年Beven和Kirkby提出了一个以地形为基础的半分布式流域水文模型TOPMODEL.该模型的主要特征是利用数字高程模型(DEM)来计算地形指数ln(a/tanβ),地形指数中的a为流经坡面任一点i处单位等高线长度的汇水面积,tanβ为该点处的坡度.用地形指数ln(a/tanβ)来反映流域下垫面的空间变化对产流区域形成和变化的影响. TOPMODEL模型主要借助地形指数ln(a/tanβ)来描述和解释径流趋势以及由于重力排水作用使得径流沿坡向的运动.该模型假设在流域任何一处的土壤上有三个不同的含水区:第一个是植被根系区,其次是土壤非饱和区,第三个是饱和地下水区(图1).水流运动规律假定如下:降水P首先下渗进入植被根系区,贮存在这里的水分部分蒸发,蒸发量为Ea,部分进入土壤非饱和区.在土壤非饱和区中的水分以一定速率Qv通过重力排水作用垂直进入饱和地下水区,然后通过侧向运动形成基流Qb.如果饱和地下水面不断升高,在流域某一山脚低洼汇合处(比如说下游河道)冒出,就会形成饱和坡面流Qs.因此,TOPMODEL模型中的流域总径流Q是基流Qb与饱和坡面流Qs之和.相应的产汇流过程见图1.TOPMODEL模型把全流域按DEM网格分块,每一个网格作为一个水文单元,对每一个水文单元流域进行产汇流计算,包括不饱和层水分运动、饱和层水分运动及饱和坡面流运动.饱和坡面流和壤中流均视为在空间上相等,可通过等流时线方法进行汇流演算,求出单元流域出口处的流量过程.采用近似运动波的常波速演算法对单元出口处的流量过程进行河道演算,得到流域出口断面处的流量过程. 在TOPMODEL模型中,除地形指数ln(a/tanβ)以外,所用的模型参数主要包括:m,为土壤下渗速率呈指数衰减的速率(m);T0,为土壤刚达到饱和时的有效下渗速率(m2/h);Td,为重力排水的时间滞时参数(d);SRmax,为植被根系区最大蓄水能力(m);SR0,为根带土壤的初始缺水量(m);RV,为地表坡面汇流的有效速度(m/h);CHV,为主河道汇流的有效速度(m/h). 尽管TOPMODEL模型的假定是对流域真实情况的近似,并非在任何地方都能适应,然而该模型结构相对简单,需要优选的“有效”参数较少,充分利用了容易获取的地形资料,另外,可以将TOPMODEL模型所得结果与野外实测资料进行比较,所以自从TOPMODEL模型被提出以来,被