基于srmdem的大凌河流域地表水文模拟模型构建.docxVIP

基于srmdem的大凌河流域地表水文模拟模型构建 dem（数字高度模型）用于表示地形的三维向量有限序列（xi，yi，zi）（i.1，2，…，n），表示区域d中的实际地形起伏，其中（xi，yi）d是平面坐标，zi是（xi，yi）的相应高度。随着遥感技术的发展,高精度DEM的获取十分便捷。同时,地理信息系统以及计算机技术的发展也使得DEM应用范围越来越广泛,如基于DEM的分布式流域水文模型。用DEM模拟地表水文的方法有多种,ArcGIS中的水文分析模块采用比较传统的D8算法。笔者利用大凌河流域SRTMDEM数据,运用ArcGIS ModelBuilder构建大凌河流域地表水文模拟模型,实现DEM填洼、水流方向的计算、流域的提取、汇流累积量的计算和河网的提取等步骤的自动执行,从而实现大凌河地表水文的自动模拟。 1 数据和方法 1.1 资料来源 SRTM(Shuttle Radar Topography Mission)是美国太空总署(NASA)和国防部国家测绘局(NIMA)联合测量的,SRTM DEM对外公开,获取方便,目前广泛应用到各种地表水文模拟中。采用4幅分辨率为90 m的STRM DEM,其DEM的行号和列号分别为(60,5)、(60,4)、(61,5)和(61,4),拼接后的4幅DEM(113.7°～125.3° E,34.6°～45.3° N)覆盖了整个大凌河流域。此外,所用资料还包括大凌河流域部分区域2000年10月TM影像以及网上部分水文数据。 1.2 地表反应和真实地表的计算 DEM中的洼地区域是指水流方向不合理的地方,在进行水文分析前需要进行洼地区域的判断。洼地有两种类型:一是由数据采集或数据内插造成的错误地表描述;二是类似湖泊、陷穴和喀什特地貌等造成的真实地表反应。对于第一种洼地需要在进行地表水文模拟前进行填平,对于第二种洼地要给予保留。这样才能将DEM 中的所有洼地填平,使之成为一个具有“水文学意义”的DEM,从而保证从DEM 中提取的河网是连续的。洼地深度是最低高程值与最低出水口高程之差。 1.3 中心栅格的水流方向 水流方向是指水流离开每一个栅格的方向,其计算是流域范围界定和水系提取的前提。D8算法的思想是:在3×3的窗口中,对中心栅格与周围栅格高程差值进行对比,利用最大坡降的方法来判断中心栅格的水流方向。水流方向由最大坡降ΔH决定,其计算公式为: ΔH=ΔZ/d*100 (1) 式中,ΔZ是与周围单元格的高程差;d是与周围各单元格的距离,它与栅格的方向有关。处于1、4、16、64这4个方向时(图1),d=1;处于2、8、32、128这4个方向时,d=2√d=2。 比较中心栅格与周围相邻的8栅格的坡降ΔH的大小,最大坡降ΔH的方向就是该中心栅格的水流流向。ArcGIS中的水流方向分析工具Flow Direction的水流方向编码规则(图1),1、2、4、8、16、32、64、128分别代表东、东南、南、西南、西、西北、北、东北方向。图2中,左边是DEM,右边是用Flow Direction进行水流方向分析后水流方向的编码。如果中心栅格到周围栅格的坡降ΔH相当,将扩大扫描窗口,直到确定方向为止;如果中心单元格比周围的8个栅格的高程值都要小,水流无法流出,则认为该中心单元格为“噪声”,将水流方向赋值为周围最低的高程值。第一行的第二列高程为27的栅格,周围栅格高程值分别为92、33、99、40、36,在右边的水流方向编码为33(图2)。 1.4 dem累积量的计算 汇流累积量是每点的汇流累积总量,它是河网提取的基础,其基本思想是:以规则格网表示的DEM每点处有一个单位的水量,按照自然水流由高到低的自然规律,按照水流方向的数字矩阵计算经过每个栅格的累积量。一个栅格的汇流累积量是其上游全部经过该栅格汇流累积量总和。 1.5 计算流域边界 由于拼接后的4副DEM有相当大的一部分在大凌河流域之外,为方便流域的提取,需要对其进行裁剪处理。先将拼接后的DEM投影到大地坐标系,再根据已知的大凌河信息,对DEM进行大概的裁剪,得到比较小的又覆盖大凌河流域的DEM,从而缩小运算区域,大大减少模型的运行时间。基于裁剪后的DEM,先计算洼地深度,判断DEM的填洼阈值。然后按照选定的阈值进行DEM填洼,得到无洼DEM。基于无洼DEM,得到水流方向,然后进行盆地流域计算,可以得到流域边界。再根据水流方向和流域边界,可以提取大凌河流域内的水系(图3)。 2 结果与分析 2.1 形成地表模型 采用ArcGIS的水文分析工具集下的Fill工具进行填洼时,如果不对填洼阈值项Z limit进行选择,默认为将所有的洼地填充。为了避免将真实地形中的洼地填充,需要计算洼地深度分布图,然后结合实际地形,可以得到一个合理的填洼阈值。而计算洼地深