基于地理信息系统GIS的水流泥沙数学模型研讨.pdfVIP

十、河流模拟技术 ·1401· I 基于地理信息系统(GS)的水流泥沙数学模型研究 陈界仁1 杨 涛2 刘 云1 (1．河海大学水利水电工程学院南京 210098； 2．河海大学水资源环境学院 南京 210098) 摘要：地理信息系统(GIS)是一个具有创建、转化、编辑、显示、查询和管理数字化地图及其相连的属 性数据系统，因此GIS在水流泥沙数学模型研究应用与其集成将成为水流泥沙数学模型的方向之一。 本文针对河道形态特点，提出了河道数字高程模型(DEM)的生成方法，基于GIS的水流泥沙数学模型实 现步骤和方法，在赣江樟树至拖船埠河段初步实现了基于GIS的二维水流泥沙数学模型。 关键词：地理信息系统 数字高程模型 水流泥沙数学模型 集成 1前言 基于地理信息系统(GIS)的数字流域研究和运用已经十分广泛，但在数字河流研究和 应用方面还不很普遍。在数字流域中采用的大尺度单元，关注的是流域出口的总水沙过 程，对每个单元格的水沙过程一般不关注，而在研究河道水流泥沙运动及工程影响时， 所遇到的问题都是局部区域，其尺度较小，要求能较精细地反映各点的水流、泥沙运动。 天然河道长宽比悬殊，边界复杂，要求生成的单元尺度小，一定程度上限制了地理信息 系统在河道水流泥沙数学模型中的应用。 数学模型在计算中涉及到大量数据的输入、输出以及数据的显示、查询、分析，传 统的方法对数据要求的格式不同，没有较好的通用性，而且其效果不好。GIS是一个具 有创建、转化、编辑、显示、查询和管理数字化地图及其相连的属性数据系统，因此GIS 在水流泥沙数学模型研究应用与其集成将成为水流泥沙数学模型的方向之一。 2数字高程模型(DEM)构建 2．1 计算前处理中数字地面模型的应用 河道地形原始采样点可以是规则的，但多数情况下是不规则的。由于客观条件的限 制，地形测点的数量和密度都不能满足要求，一般通过内插可以对数字地面模型进行加 密，还可以对不同结果的数字地面模型进行变换，将离散型分布的数据点转化成规则网 格分布的数据点，即离散数据的网格化。 采用数字地面模型中的规则矩形网格。计算前处理需要将数字地图上的高程或和 水深数据离散到网格节点上，即前面所说的离散数据网格化。由于数字地图上高程和 作者简介：陈界仁，男，1964年出生，江苏靖江人，博士，教授。 ．1402· 第六届全国泥沙基本理论研究学术讨论会论文集 水深数据的分布不是规则的，应采用一定的算法从网格上提取高程和水深数据。目 前对于规则矩形网格中数据的提取主要有如下几种算法：①与距离成反比的加权法； (室)MQS(ModifiedQuadratic 的加权方法的基础上提出的；③径向基函数插值法：径向基函数中基函数是由单个变量 的函数构成，主要有Multiquadric方法和薄板样条法；④有限元法。 本文中采用与距离成反比的加权法，在网格中使用圈选，设置一个最小搜索半径、 一个最大搜索半径和增长步长，在网格中循环搜索，搜索步骤如下： (1)在最小半径圈内进行搜索，若圈选范围内有水深数据，则用与距离成反比的加权 平均法，将网格节点与各数据点之间的距离的倒数乘以水深数据进行相加，再除以网格 节点与各数据点之间的距离的倒数之和，即为节点处的水深值； (2)若在第一次搜索中，圈选范围内没有水深数据，则在最小半径上增加一个增长步 长，扩大范围，进行第(1)步的搜索工作。若在最小半径与最大半径范围内，都未搜索到 水深数据点，则系统进行报错处理； (3)循环搜索，直至所选区域内所有的网格节点都被赋以水深值。 当计算区域内所有网格节点都赋完高程和水深值后，统一提取所有节点上的高程和 水深值，按专业计算模型所需的格式，形成高程和水深数据参数文件，并进行保存。 2．2计算后处理中数字地面模型的生成 GRID建模方法有如下三种：从测点、从线数据集和从面数据集生成。在由点、线 和面数据集生成GRID的方法中，GRID网格角点的属性值由用户任意指定，这样用户 可根据自己的需要来构建GRID模型，如属性可以是土壤类型、密度值、高程、温度、 湿度等。属性甚至可以是对象的标识号、拓扑误差