[发明专利]一种水量动态分配的有源积水扩散方法

说明书

本发明公开了一种水量动态分配的有源积水扩散算法，在计算出每个扩散源的积水量后，以扩散栅格为中心，根据水量平衡原理，在中心栅格及其周围八个栅格中扩散，每扩散完一圈进行扩散结果的合理性判断，判断已形成的积涝水位是否大于其外圈某个栅格的高程值，如果存在，则将扩散范围增加一圈再次进行扩散，直到积涝水位小于其外圈的任何一个栅格的高程值。在所有节点扩散完成之后，判断是否存在栅格参与了多次节点扩散过程，若有，则对这些节点进行合并扩散处理。在扩散过程中，从地域连通性及地形阻碍作用出发，不断判断积水路径的合理性，动态地进行水量分配，避免了人为选择多流向影响因子的主观性，提高了暴雨积水模拟的准确性。

技术领域

本发明涉及一种水量动态分配的有源积水扩散算法，属于城市内涝监测预警技术领域。

背景技术

城市暴雨积涝模拟过程主要分为降雨过程模拟、产汇流计算、暴雨淹没分析。城市暴雨积水即指降雨下落城市区域形成径流后由于地势低洼、排水不及时等原因造成的一定径流深度的地面积水。城市暴雨淹没分析是指在求得城市暴雨积水总量后，根据城市地形将积水合理分布到区域空间中，得到城市可能的积水深度以及积水的淹没情况。对于这一问题的研究基本从两个角度进行：(1)给定水位下的积水淹没分析，这些算法都是基于递归算法改进的角度进行的，需要给定积涝水位，但是在实际降雨中积涝水位获取存在一定难度；(2)已知洪水体积的积涝扩散：有些算法未突破汇水区边界限制，实际积水扩散时并不存在任何边界，会导致汇水区边界处的扩散结果不合理。有些算法虽然不存在边界限制问题，但其扩散过程可能会使得一个节点形成的扩散区域的积涝水位不一致，这与实际水塘形成的情景不一致。

数字高程模型(Digital Elevation Model，简称DEM)，是以数字的形式按一定结构组织在一起，表示实际地形特征空间分布的模型，也是地形形状大小起伏的数字描述。最基本的DEM由一系列地面点(x,y)(x为广义的该点的经度坐标值，y为广义该点的纬度坐标值)位置及其相联系的高程Z所组成，用数学函数式表达是：Z＝(x,y)。DEM 通常用地表规则网格单元构成的高程矩阵表示，广义的DEM还包括等高线、三角网等所有表达地面高程的数字表示。在本发明中，DEM是指覆盖研究区域的每个规则网格的高程值的集合。

发明内容

本发明的目的在于克服现有技术存在的不足，提供一种计算速度快准确度高的水量动态分配的有源积水扩散算法。

本发明采用的技术方案为：一种水量动态分配的有源积水扩散算法，包括以下步骤：

(1)读取DEM数据，获得区域每个栅格的高程值；

(2)生成与DEM栅格一一对应的网格编号数据，该数据唯一地确定了每个栅格的编号，读取编号数据，并根据经纬度坐标与DEM数据构成一一匹配的关系，并将其作为DEM 栅格的一个属性；

(3)根据地理位置，确定扩散节点(即扩散源)的栅格编号及其积涝水量，同样将节点的编号与积水量作为DEM栅格属性。

(4)从第一个节点开始，依次进行积水扩散，当前节点的溢流量为V_i，假设积涝高程为H，每个DEM栅格的面积为area，节点所在栅格及其他周围的栅格与积涝高程的高程依次为X₁,X₂…X₉，则根据水量平衡原理有((X-X₁)+(X-X₂)+…+(X-X₉))*area＝V_i，得到X。

(5)依次判断X与相邻8个栅格相邻的16个栅格的高程值差，当所有的差值都小于设定的阈值s时则扩散结束。

(6)当出现一个差值大于阈值s时，将扩散范围扩大到为以节点为中心的5*5栅格，然后重复步骤(4)得到所有节点扩散完成。(当一个栅格参与了两个及以上的节点扩散时，其每次的高程以原始高程参与计算。)