流域单次降雨侵蚀过程的动态模拟计算.docxVIP

流域单次降雨侵蚀过程的动态模拟计算 一、 与gis的结合 土壤侵蚀模型是定量评估土壤侵蚀影响、土地利用和土壤规划的科学工具。近年来, 土壤侵蚀模型由经验模型向基于侵蚀机理的过程模型发展, 和由单坡面模型向流域模型发展。GIS以其对空间数据的输入、存储、查询分析、管理和显示输出能力, 日益成为土壤侵蚀模型应用不可或缺的工具, 并出现了大量的侵蚀模型与GIS结合的实例。将流域侵蚀过程模型与GIS结合, 研究流域径流、侵蚀、产沙的时空变化规律, 是土壤侵蚀模型发展的一个趋势。 GIS与侵蚀模型的结合主要有三种方式:松散结合、交互界面结合和完全结合。松散结合主要指模型与GIS是各自独立开发的情况。首先用GIS对空间数据进行查询和预处理, 然后按模型所需的格式将数据输入模型, 最后将模型计算结果再转为GIS文件格式, 进行显示和分析。松散结合的优点是无须改变模型的代码, 而不足之处是有大量的数据管理和转换工作。交互方式是开发一个交互界面, 为模型提供输入数据, 以及对模型计算结果进行处理和演示, 所有数据转换是通过交互界面自动进行。此类结合实例最多, 以美国普度大学集成GRASS与ANSWERS和AGNPS模型的工作最为典型 PCRaster 二、 pcrast空间动态模型语言 1. 个图层与点数据文件 PCRaster系统采用四种文件来表达空间动态数据。栅格图层文件是PCRaster最重要的文件, 每个图层文件存储了地表景观的一个属性特征, 各属性特征均与某种数据类型相联系 (表1列出了PCRaster的数据类型) ;点数据文件包含点数据的地理坐标和一个或多个属性数据, 用于表达点与图层文件的关系;交叉表文件用于表示各图层之间的关系;时间序列文件则与图层文件结合, 用于表达空间实体的时-空变化特征。 2. gis空间分析函数 PCRaster提供了120多个空间和时间命令函数。常用的命令包括单点运算 (代数、逻辑、比较、关系和函数运算) 、邻域运算 (高通滤波, 移动平均等) 和面域运算 (类别统计、耗费距离等) 。PCRaster最初的发展是为了解决地貌和水文问题, 所以提供了一系列专门用于解决地貌和水文问题的GIS空间分析函数, 包括坡面和集水区分析、物质输移 (如径流-泥沙输移) 分析等。 PCRaster系统的动态模型语言在地图代数和制图模型语言的基础上, 扩展了动态迭代模型所需要的动态输入输出命令函数, 它大大简化了动态模型与GIS数据库之间的数据交换工作, 在形式上也类似于静态的制图模型命令。PCRaster空间动态模型语言与制图模型语言在概念上有着本质区别。制图模型语言是一种自然语言, 而PCRaster的动态模型语言命令则类似于数学方程, 每个命令函数计算一个表达式, 得到一个结果, 这对熟悉数学函数的研究人员建立动态模型是很方便的。如泥沙输移能力方程: G=C 对于栅格图层CoverFactor (C) , FlowVolume (Q) , Slope (S) 和指数1.7, 写成PCRaster的空间动态模型语言的命令函数: TransportCapacity=CoverFactor 3. 动态模拟的程序 用PCRaster空间动态语言编写的动态模型按逻辑关系可划分为5个模块 (如图1所示) :数据联结 (binding) 、空间位置属性定义 (areamap) 、时间控制 (timer) 、初始化 (initial) 和动态模拟 (dynamic) 。 数据联结模块是模型与GIS数据库输入/输出的接口。由于模型使用的变量与GIS数据库文件相对应, 数据联结模块可以为这些文件定义变量别名, 使得数据库的文件与模型的变量联结在一起;空间位置属性定义模块 (areamap) 决定了模型模拟区域的面积和空间分辨率。它用一个全局图层定义模型模拟区域的空间位置属性, 所有的图层文件均有与全局图层文件一致的地理坐标、范围和网格大小;时间控制模块 (timer) 通过设置模型的开始时间、结束时间和时间步长来决定模型的历时和模拟次数;初始化模块 (initial) 用于设置模型的初始条件, 空间属性和非空间属性数据均可在此赋初值。它由一系列pcrcalc操作命令组成, 可以看作是静态制图模型语言为动态模型的第一步循环赋初值;动态模拟模块 (dynamic) 定义了每一时间步长内的所执行的一系列命令操作, 它们由一系列pcrcalc操作命令组成, 在每一时间步长内从顶向下顺序执行, 计算的结果作为下一时间步长的初始值。第一步计算使用初始化模块中定义初始值, 第二步用第一步的结果, ..., 循环迭代直到结束时间为止。在每一时间步长内既可以从GIS数据库读取数据, 也可将计算结果存储到数据库中。 动态输入输出语句 (timeinpu