分布式TOPKAPI模型的数值解法及其应用.pdfVIP

第 37卷第2期 河 海 大 学 学 报 (自然 科 学 版 ) Vo1．37No．2 2009年 3月 JournalofHohaiUniversity(Natural~iences) Mar．2009 DOI：10．3876／j．issn．1000—1980．2009．02．007 分布式 TOPKAPI模型的数值解法及其应用 邓 鹏 ，李致家 ，刘志雨2 (1．河海大学水文水资源学院，江苏 南京 210098；2．水利部水文局，北京 100053) 摘要：为了提高分布式TOPKAPI模型的计算精度，用数值分析的方法对模型非线性水库方程进行 求解．通过分析得出非线性水库方程的一般格式，运用四阶龙格库塔算法建立求解数值解的递推关 系式，并且通过变步长链提高算法的效率，保证递推关系式的收敛和稳定．利用变步长的四阶龙格 库塔算法在 1km网格精度下对布柳河流域进行洪水模拟，模型率定期和检验期的确定性系数分别 达到了0．908和0．912． 关键词：TOPKAPI模型；龙格库塔 ；布柳河；洪水模拟 中图分类号：P338．2 文献标识码 ：A 文章编号：1000-1980(2009)02—0158—04 TOPKAPI模型是一个以物理概念为基础的分布式流域水文模型．模型 自1995年提出至今经过了一系列 的发展．最新的TOPKAPI模型对壤中流计算进行了调整，增加了土壤水深层渗漏以及其他模块[1]．TOPKAPI 模型可以利用数字高程图(DEM)、土壤分布图、土地利用图等信息识别模型结构和参数．模型假定土壤及地 表网格内侧向水流运动可用运动波模型来模拟，将建立在空间点上的假设在一定的空间尺度上进行积分，从 而把初始的微分方程转变为非线性水库方程 ． 最新的TOPKAPI模型中非线性水库方程采用基于适当近似的解析法求解 ．该方法对指数项用二次方程 近似，二次方程的系数通过最小二乘法得到，从而将复杂的指数方程转化为低阶方程并得到解析解．但该方 法在对指数项的近似中有明显的舍入误差，并且二次方程系数也难 以确定．为了提高模型的计算精度，本文 采用变步长的四阶龙格库塔算法求解TOPKAPI模型方程，并将该算法应用于布柳河流域． 1 模型基本原理 TOPKAPI模型的流域特性参数、降水和水文响应的空间分布在水平方 向用正交网格系统 (DEM的方网 格)，在垂直方向用各网格所对应的水平土柱进行模拟．每个网格均是体现模型物理特性(质量平衡和动量平 衡)的1个计算节点．根据相邻网格的最小能量原理，将活动网格与周边4个相邻网格联系起来，并假定活动 网格只能与 1个下游网格相联系，但可以从其他 3个网格接受上游来水．在上层非饱和区有： oJv 1 = [(fax+q。+qus)一fbXl一 s (1) 式中：Vl——网格单元内的上层蓄水量；f——时间；-厂a——降水人渗；——网格大小；qno——地面径流输 入；qus——壤中侧向人流；^——深层渗漏；C——局地传导系数；a。——反映上层土壤透水特性的指数．与 壤中径流的模拟相似，应用曼宁阻力定律，近似用运动波方程来描述水在地面上和河道中任一处水平方向上 的运动，将偏微分运动波方程积分到网格尺度上，就能分别得到结构相似的非线性水库方程： OVo =r0 一 (2) u dvc ： XW+Quc)一 5 (3) 收稿 日期：2008—04—23 基金项目：国家 自然科学基金 (5407~17) 作者简介：邓鹏(1983一)，男，江苏兴化人，博士研究生，主要从事水