实验16 地表水文分析.docxVIP

实验16 地表水文分析 实验目的 了解单流向水文模型中的基本概念和建立地表径流模型的基本方法； 熟悉水文学上出水口，河网，集水流域，径流强度指数，地形湿度指数等基本概念； 掌握利用ArcGIS中的水文分析工具来提取流域和地表径流特征；注意区别栅格和矢量两种特征表达的方法；径流强度指数和地形湿度指数的计算方法。 实验背景 地表水文分析应用于与地表水流有关的诸多领域，对地球表面形状的理解也具有十分重要的意义。利用ArcGIS的水文分析工具，建立地表水的运动模型，辅助分析地表水流从哪里产生以及要流向何处，再现水流的流动过程，有助于了解排水系统和地表水流过程的一些基本的概念和关键的过程。 无洼地DEM 由于数据噪音、内插方法的影响，DEM数据中常常包含一些“洼地”（被较高高程区域围绕的低地）。有些洼地是在DEM生成过程中带来的数据错误，有些则表示真实的地形如采石场或岩洞等。“洼地”将导致流域水流不畅，不能形成完整的流域网络，因此在利用模拟法进行流域地形分析时，要首先对DEM数据中的洼地进行处理。 去除洼地最常用的方法是对DEM做平滑处理，然后再填充。其方法是利用一个小窗口在整个DEM上进行遍历。首先判断窗口中心点是否比四周点的高程都低，如果低，即认为这个点是洼地点。若是洼地点，泽将洼地点的高程值改为它4个方向中高程最小的值。则将洼地点的高程值改为它4个方向中高程最小的值。 单流向模型 水流方向的计算是水文分析的关键步骤，用它来进一步确定汇水面积，水流长度，汇水流域等重要的水文参数。水流方向编码如图16-1所示。 水流的流向是通过计算中心栅格与邻域栅格的最大距离权落差来确定的。距离权落差是指中心栅格与邻域栅格的高程差除以两栅格间的距离，栅格间的距离与方向有关，如果邻域栅格对中心栅格的方向值为2、8、32、128，则栅格间的距离约为1.414，否则距离为1。如果高程差为正值，则为流出。为负值则为流入。 汇流积累量 在地表径流模拟过程中，汇流累积量是基于水流方向数据计算而来的。对每一个栅格来说，其汇流累积量的大小代表着其上游有多少个栅格的水流方向最终汇流经过该栅格，汇流累积的数值越大，该区域越易形成地表径流。由水流方向数据到汇流累积量计算的过程如图16-1所示。汇水面积为汇流积累量乘以每个像元的面积。 图16-1 流水累积量的计算 河网 基于地表径流漫流模型，当汇流累积量达到一定值的时候，就会产生地表水流，那么所有那些汇流量大于那个临界数值的栅格就是潜在的水流路径，由这些水流路径构成的网络，就是河网。 出水口 确定小级别的流域的出水口的位置（汇流点）是进行流域分割的前提。结合数字高程模型和高分辨率的航空像片来目视确定小级别的流域出水口的位置；接下来，利用ArcToolbox中spatial analysis tools\hydrology\snap pour point捕捉节点工具使出水口的位置符合DEM模型。snap pour point的思想是利用一个记录着point点栅格数据，在这个数据层中，那些属性值存在的点作为潜在的出水点，在该点位置上在指定距离内在汇流累积量的数据层上搜索那些具有较高汇流累积量栅格点的位置，这些搜索到的栅格点就是小级别的流域的出水点。 集水流域 集水流域（集水区域；流域盆地）是地表水与地下水分水线所包围的集水区或汇水区，习惯上以地表水的分水线来确定，可以根据地形图勾勒出流域的分水线。 径流强度指数 径流强度指数（Stream Power Index，SPI）用来测量水流侵蚀能力，可以被定义： 其中：As是汇水面积；β是坡度。 小提示： 函数tan要求坡度的单位是弧度，所以需要将提供的坡度单位（度数，degree）转为弧度（radians），具体方法是: radians = degrees *π / 180 地形湿度指数 地形湿度指数(topographic wetness index)可定量模拟流域内土壤水分的干湿状况,在流域的土壤及分布式水文模型等研究中具有重要的意义。地形湿度指数可以被定义为： 其中：As是汇水面积；β是坡度。 实验内容 基于单水流模型，利用DEM数据生成流域的基本结构：河网，出水口，集水流域，集水盆地。 利用DEM数据计算径流强度指数和地形湿度指数的方法。 实验数据 数据路径：E16 路径 文件名称 格式 说明 \ \mn_ws_clip grid 30米分辨率的数字高程模型 \ \orthoimage image 研究区1米 实验步骤 空间分析环境的设置 分别在ArcMap和ArcCatalog的主菜单上点击Tools\Extensions,打开Extensions对话框中，点击选中Spatial Analysis，3D Analyst和Network Analysis前的复选框，添加拓