还剩2页未读,继续阅读
文本内容:
流域水文是一个复杂的过程,其动态变化不仅受气候条件、土壤因子等的影响,植被在其中也起着非常重要的作用,只有将三者有机的结合起来,才能客观的反映流域水文动态,才能有效地对森林的水文功能进行评估本研究拟采用的水文模型为TOPMODEL,气候模型为MT-CLIM,考虑到本研究所涉及的范围,植被模型只选取蒸散发和林冠截流两部分模型TOPMODELTOPMODEL由Beven和Kirkby在1979年提出,是一个基于地形的半分布式模型其理论基础为变动源面积Variable SourceAreas,根据这个概念坡面流只在整个流域的一小部分上发生,产生坡面流的陆地表面是那些在降雨事件中地下水位上升至地表的饱和区域O Loughlin1981年结合GIS和现代图形技术增强了其模拟和可视化能力TOPMODEL被发展主要是作为一种研究工具,并且已经广泛地运用于实际当中该模型的特点是结构简单,优选参数少,充分利用了容易获取的地形资料,而且与观测的物理水文过程密切联系的地形指数TOPMODEL地形指数的概念地形是影响流域对降雨响应的重要因子,它决定了重力作用下流域中水分运动的趋势,因此也影响着水文系统的许多其它方面TOPMODEL提出了地形指数Topographic Index的概念,并假定相同地形指数的水文单元上有着相似的水文学特性,而不考虑其所在的位置地形指数是水文模型的重要组成部分,它反映了土壤湿度、地表饱和度的空间分布,和径流生成的过程Zhang andMontgomery,1994,因此可以说地形指数是TOPMODEL的核心地形指数tan
5./T=lnA:特定集水区面积;P坡度角地形指数的计算地形指数通常用DEM计算,所选用的算法与所定义栅格的尺度或精度影响着计算结果计算地形指数时,首要的任务是确定特定集水区面积Specific catchmentarea,特定集水区面积由上坡面积Upslope area,a求得a被定义为一点point或等高线上总的集水区面积Moore等,1991年工A为每单位宽度L等高线的上坡面积A=a/L XA是一个有重要水文、地理意义的分布量Burges等,1994年,上坡面积和特定集水区面积的计算依靠于流向路径Flow direction的计算本文采用David G.Tarboton1997年提出的一种简明有效的单流向算法来计算特定集水区面积的基本方程TOPMODELTOPMODEL基本方程的数学出发点是连续性方程和Darcy定律,TOPMODEL应用中有以下几点假设1土壤饱和水力传导度G随着土壤厚度的增加服从指数递减;K=…式中KO:地表土壤的水力传导度;f:为参数,表示K随z的下降速度2饱和壤流受地形梯度驱动,被地下水位控制,状态稳定,并且地下水位线和饱和壤中流平行与地表;Qb=Q efl=阴$nMQ°=式中Qb饱和壤中流;Qo土壤饱和时的排水量;AB土壤可排水空隙率;S:流域平均;;土壤饱和缺水量To饱和导水率A:lna/tanp;3相同地形指数的单元,有着相似的水文学特性根据以上假设建立方程式中Sii处饱和缺水量;S流域平均饱和缺水;m水力传导度随深度递减的模型参数;一Si3-lnTja:排水面积;0:坡度;Ti:i处的饱和导水率;入流域Ina/tan0的平均值;6:流域In Ti的平均值通常Ti的空间分布是未知的,因此假定Ti在整个流域上为常数,方程变为,这个方程即为TOPMODEL的基本方程A-lnf—Si Itan
4.由于Si代表着局部饱和缺水量,当任一单元Si0时则表明土壤已经饱和,坡面流将发生从公式中还可看出地形指数的高值单元更容易达到饱和而产生坡面流的产流计算TOPMODELTOPMODEL产流计算包括壤中流Qb和饱和坡面流Qs两部分,其中壤中流计算如式,坡面流则需先计算流域内局部饱和缺水量的情况,模拟过程中,局部饱和缺水量是随时间变化的函数,可用下式表达式中St:t时刻土壤饱和缺水量;St;:t-l时刻土壤饱和缺水量;Qt-i t-l时刻径流量;r:净补给率;At:模拟的时间步长;坡面流的计算方程为式中At:整个流域的面积;As:流域中SiW0的面积。