GIS技术支持下的洪水模型建模_万洪涛资料.pdfVIP

　　文章编号:1000-0585(2001)04-0407-09 　　收稿日期:2001-03-19;修订日期:2001-07-10 　　基金项目:中国科学院知识创新工程领域前沿项目 (CXIOG -D00-05);STA fellowship(100025) f rom JST , Japan资助. 　　作者简介:万洪涛 (1969-), 男 , 江西吉安人 , 博士 , 中科院遥感所博士后 , 主要从事地理信息系统与水文模型 的集成 、地理信息系统及遥感在流域管理中的应用等方面的研究工作。 GIS 技术支持下的洪水模型建模 万洪涛1 , 2 , 周成虎2 , 万庆2 , 汪承义2 (1.中国科学院遥感应用研究所 , 北京　100101; 2.中国科学院北京地理科学与资源研究所资源与环境信息系统国家重点实验室 , 北京　100101) 摘要:在复杂区域建立洪水模型时 , 计算网格的手工生成方法容易出错甚至不可行 , 自动生 成算法则可大大节省计算网格生成的工作量。洪水模型中的计算网格与 GIS 栅格数据及不规 则三角网空间数据结构非常相似 , 因此 , GIS 中成熟的网格自动生成算法可用于生成洪水模 型计算网格。文章详细讨论了 GIS 支持下的洪水模型自动建立步骤 , 并以黄河下游花园口～ 夹河滩河段为例 , 利用地形图 、 土地利用图 、 水利工程设施分布 、 水文站点图等资料 , 通过 自动生成网格及其空间拓扑关系 , 建立了洪水过程数值模拟模型 , 并详细解释了计算网格数 据格式。 关　键　词:地理信息系统;洪水模型;计算网格;拓扑关系;模型建模 中图分类号:P208;TV122　　文献标识码:A 1　引　言 进行二维洪水模拟计算时 , 必须建立一个能够描述计算域的几何形状 、物理属性和边 界条件的计算模型。对于较简单的问题 , 输入的数据较少 , 所需要的数据手工处理即可。 但对于中到大尺度模拟计算问题 , 手工生成格网数据就变得容易出错甚至不可行。而且在 生成格网拓扑数据时 , 一旦无错误的模型生成 , 如果模型需要修改 , 如计算区域扩大或缩 小 , 模型则须从新生成 , 这必然造成大量的重复劳动[ 1] 。Haber估计[ 2] , 以前进行的水流 模拟 , 80%的费用花在采用传统手工方法获取数据。因此 , 自动生成计算格网以及为格网 编码或建立拓扑关系 , 有可能大大减少建立格网的工作量并减少模型计算量。格网自动生 成和编码算法各不相同 , 其取决于处理程序和分析程序的类型 。一般将计算域离散成规则 格网或不规则格网 , 这与 GIS栅格数据结构 (GRID)及不规则三角网 (T IN)管理空间 数据方式非常相似 , 因此 , GIS中有不少格网自动生成算法可用于生成水文模型中的计算 网格[ 3 , 4] , 能提高空间子单元的定义能力 , 可对计算域进行更好的离散 , 使计算域的数目 达到一定的要求;利用 GIS 工具还能够准确地内插得到子计算域的空间参数[ 5] 。建立洪 水模型时 , 需用到 GIS 数据库中的等高线 、 水系 、 水利工程数据 、土地利用图 、 居民地 分布图等空间分布数据;其次 , GIS 中的数据处理与空间分析功能模块 , 可用于洪水模型 第 20卷　第 4期 2001 年 9 月 地　　理　　研　　究 GEOGRAPHICAL　RESEARCH Vol.20 , No.4 Sept., 2001 自动建模 , 我们据此设计的 GIS 支持下的洪水模型建模步骤为:确定计算域范围 、 生成 计算网格 、设置边界与初始条件以及模型参数等[ 6] 。 2　二维洪水模型建模过程 2.1　计算域的生成 在洪水模型建立之前 , 应该确定模型计算域 。计算域既可用一个已存在的边界数据文 件确定 、 也可在背景图中勾绘得到 。即根据等高线图选择一个较历史最大淹没水深稍高的 等高线 , 根据该等高线包围的范围确定计算域的范围 。如果选取的等高线不连续 , 还应选 取其他的特征线 (如海岸线 、 湖岸线 、 河岸线 、 道路等), 使选取的线性要素组成一个封 闭的多边形 , 从而得到计算区域边界。 2.2　计算网格生成 2.2.1有结构网格和无结构网格 　　得到计算域边界之后 , 必须考虑如何对之进行离散 , 空间离散可采用有结构格网或无 结构格网 。有结构格网一般是一个逻辑长方形点阵 , 这些点阵为已定义好且与其相邻点具 有规则的关系。例如 , 一个典型的二维有结构格网 (图 1), 任意内部结点有四个相邻结 点 , 并可将函数 f(i , j)规则地表达在格网上 ,其右邻接点和上邻接点分别为 f(i +1 , j)和 f(i , j +1)。不规则结构网格可看作为规则正交网格 ,结点间的连线可