分布式水文模型中栅格汇流演算顺序的确定.pdfVIP

三、水力学、防洪与泥沙 ·305· 分布式水文模型中栅格汇流演算顺序的确定* 熊立华 (武汉大学水资源与水电工程科学国家重点实验室 武汉430072) 【摘要】在基于DEM的分布式流域水文模型中，确定流域DEM中每个栅格的汇流演算顺序是准确描述水 流时空分布的关键之一。本文比较研究了两种确定栅格汇流演算顺序的方法：随机抽取法和分级确定法。 作为例子，本文采用分级确定法来确定一个小型流域的栅格汇流演算顺序，并给出了结果。对于栅格数目很 大的情况，建议采用分级确定法来确定流域DEM中栅格汇流演算顺序。 【关键词】数字高程模型分布式水文模型汇流演算 在基于DEM的分布式流域水文模型中，流量演算包括两个步骤：先进行栅格汇流演 算；再进行河网汇流演算。同河网汇流演算一样，确定流域DEM中每个栅格的汇流演算顺 序是准确描述水流时空分布的关键之一。河网汇流演算主要基于河道节点之间的上下游关 系；而栅格汇流演算主要基于每个栅格和周围相邻栅格之间的水流联系。由于一个流域上 的栅格数目远远大于流域水系中河道节点的个数，因此确定流域DEM中每个栅格的汇流 演算顺序比确定流域水系中河网汇流演算顺序要复杂一些。本文主要详细研究两种确定栅 格汇流演算顺序的方法：随机抽取法和分级确定法。在讨论这两种方法之前，假定一个流域 的DEM上的每个栅格的流向已被唯一确定(通常采用最陡坡度法)。 1随机抽取法 通过随机抽取来确定栅格汇流演算顺序的步骤如下。 格的汇流演算顺序，咒为整数。 (2)赋初值Seq(1：N)=一1，以及7z=0。 (3)从演算顺序未定的所有栅格中随机抽取一个编号为．j的栅格。 (4)对编号为J的栅格，找到其所有的相邻上游栅格，然后考虑两种情况。如果其所有 的相邻上游栅格的演算顺序已定，那么，z的值加1，即矩一雄-I-1；并且确定编号为歹的栅格 的演算顺序为Seq(j)=咒。如果在其相邻上游栅格中，至少有一个的演算顺序还未定，表明 此次随机搜索无效。 (5)如果咒N，返回步骤(3)和(4)；否则，终止程序。 随机抽取法确定栅格汇流演算顺序的思路清楚简单，缺点是工作量比较大，效率不高。 师基金(101077)资助。 中国水利学会第三届青年科技论坛论文集 2分级确定法 分级确定法的思路是先定出各个栅格的汇流级别，然后再确定栅格的演算顺序。。 流级别，用数组Seq(1：N)来寄存每个栅格的汇流演算顺序。k、，z为整数。 (2)赋初值Order(1：N)=一1，Seq(1：N)=一1，以及k=0，7z=0。 (3)栅格分级的方法如下。①首先找到流域出口栅格，将k的值加1，即志一是+1，然后 流域出口栅格的级别确定为k。②对级别为k的一个或者多个栅格，找到它们所有的相邻 上游栅格。令忌一惫+1，然后将所有这些相邻上游栅格的级别确定为k。③重复步骤②，直 到数组Order(1：N)中每个元素的值都大于0。所得到的最大的栅格级数记为K。记同一 级别的栅格个数为t(k)，满足条件∑”-￡(s)=N。栅格分级的示意图参见文献[2]。 (4)对k=1而言，这种级别的栅格只有一个，是流域出口，故t，=1。假设其对应的栅 格编号为z，则S叼(z)=N。 ． (5)对忌=2，…，K中的每一个取值，重复以下步骤。找出级别为k的所有栅格，总个数 为t(忌)，所有对应栅格的编号用Z[1：t(”]来表示。该级别的所有栅格的演算顺序由下式 确定：Seq[1(m)]=(N+1)一m一∑：：：一￡(s)，m=1，2，…，￡(惫)。‘ 该方法虽然步骤较随机抽取法要多，但是效率较高。 3计算实例 ． 作为例子且便于图形显示，所选择流域的DEM包含的栅格数目不是很大。该DEM矩 此流域的面积大约为440km2。该流域的数字高程见图1。图2用颜色变化来更清楚地反 映流域高程的变化。在确定栅格流向时，一般采用不同的整数值来表示不同的栅格流向。 图3表示的是比较常见的栅格流向标注方法，比如正左流向用16来表示，右上角流向用 128来表示。采用图3的流向标注法后，采用最陡坡度法确定出流域DEM栅格的流向分布 图(见图4)。有了流域DEM栅格的流向分布图，