1D+洪水淹没制图软件（洼地处理、流向计算、集水面积计算、绘制洪水淹没图）DEM GIS HecRAS

1. 软件原理简介

基于洪水频率计算得到的设计洪水成果，通过水位流量关系将常见年遇洪水的洪峰流量转化为水位，进一步结合洪泛平原DEM数据，采用GIS水文分析方法计算洪泛平原内任一栅格的水流流程，耦合一维水力学模拟计算，用于大范围洪水淹没制图。淹没制图算法近似于HecRAS等水力学模拟软件算法。

1.1 输入数据

1D+洪水淹没制图软件V1.0版的输入地形数据为数字高程模型DEM（Digital Elevation Model）。DEM是以数字形式表达地面高程空间分布的模型，包含了丰富的地表形态信息。DEM作为重要的基础地理信息数据，在水文、水力学研究中主要用于水文特征分析、水文过程建模和水力特征提取等。DEM规则格网模型将空间分割成规则的格网单元，以数字矩阵的形式表达高程属性。栅格型DEM简单直观、易于操作，在水文过程模拟中较为常用。1D+洪水淹没制图软件V1.0版采用栅格型DEM。1D+洪水淹没制图软件V1.0版输入的DEM数据格式为ESRI ASCII格式（图1），文件扩展名为“.asc”。所输入的样例研究区DEM样例数据展示如图2所示。

原文链接：1D+洪水淹没制图软件（洼地处理、流向计算、集水面积计算、绘制洪水淹没图）DEM GIS HecRAS

图1. 1D+洪水淹没制图软件V1.0版的输入DEM数据格式

图2. 1D+洪水淹没制图软件V1.0版的输入DEM样例数据展示

1.2 1D+洪水淹没制图计算流程

洼地处理

DEM中往往存在着一些凹陷区域，称之为洼地。栅格型DEM中存在洼地栅格，其相邻栅格的高度均高于中心栅格，形成“井”。洼地的存在会使得洼地栅格及其上游栅格的地表产流不能流出，造成提取的数字河网不连续。1D+洪水淹没制图软件V1.0版采用下挖的方法对洼地及平原进行处理，填洼结果使得平原、盆地等低洼地区的河网流域提取结果更加准确；

流向计算

1D+洪水淹没制图软件V1.0版采用经典且相对较为简单的D8流向算法。如图3所示，用0~7来表示中心栅格的8个方向，分别对应正南、东南、正东、东北、正北、西北、正西和西南，流域出口的流向记为9。

图3 数字表示流向

集水面积计算

栅格的集水面积是指水流入该栅格的上游面积。对于局部高点，除了自身产流外，无其他栅格水流入，那么局部高点的集水面积为1个栅格面积。从各局部高点出发，沿水流方向累计可得到汇流路径上各栅格的集水面积。样例研究区集水面积计算结果在ArcGIS的展示如图4所示。

图4集水面积计算结果展示

河网定义

产流在重力作用下从高处流向低处，经坡面汇流过程流入河道，再经河道汇流过程流出流域出口。显然流域内集水面积越大的地方越可能是河流，因此可通过将集水面积不低于某一阈值的栅格标记为河流来生成数字河网。1D+洪水淹没制图软件V1.0版可人为输入集水面积阈值，也可根据水面面积比例自动生成。样例研究区河网定义计算结果展示如图5所示。