Arc Hydro模型提取流域水文信息及精度分析

Arc Hydro模型提取流域水文信息及精度分析―――以松花江流域为例
作者：李勇于宏兵艾丽娜张时佳耿丽娟
摘要: 流域水文模型已成为流域水文信息提取的重要手段。

运用Arc Hydro 模型以松花江流域为研究对象，以分辨率为90 m的SRTM DEM 为数据源，以松花江流域1∶ 25 万水系矢量数据为参照，分别在400，500和600 km2 三个不同集水面积阈值下进行河网提取。

运用图像比对软件Alike Duplicate Image Finder 对所提取河网进行精度分析。

结果显示: 集水面积阈值为500 km2 下提取的河网最为精确。

在500 km2 集水面积阈值基础上，运用ArcHydro 模型将松花江流域划分为642 个子流域。

基于Arc Hydro 模型提取的松花江流域河网与实际河网基本一致，说明利用该模型进行流域水文信息提取是合理可行的。

关键词: Arc Hydro 模型; 集水面积阈值; 水文信息; 子流域提取; 精度分析中图分类号: P334．92 文献标识码: A 文章编号: 1672-643X(2021) 06-0120-04 Application of Arc Hydro model to extract hydrological information and accuracy analysis: Songhua River basin as an example
LI Yong，YU Hongbing，AI Lina，ZHANG Shijia，GENG Lijuan
Abstract: Watershed hydrological model has become an important means of extraction hydrological information．The paper studied Songhua river basin based on SRTM DEM with 90 m resolution and 1: 250，000 river vector data．Arc Hydro model has been applied to extract river network on different catchment area threshold of 400，500 and 600 km2 ．The image alignment software such as Alike Duplicate Image Finder has been used to analyze extraction accuracy of river network．The results show that the river network is accurate when the catchment area threshold is 500 km2 ．Arc Hydro model is used to generate 642 sub-watersheds when the catchment area threshold is 500 km2 ．The river network of Songhua river basin extracted by Arc Hydro model is consistent with the actual one，which shows that using the model to carry out hydrological information extraction is reasonable．
Key words: Arc Hydro model; catchment area threshold; hydrological information; sub-watershed extraction;accuracy analysis
流域作为我国水资源管理的基本单元之一，其水文信息的获取十分重要。

传统的流域水文信息获取主要依靠人工手段，不仅消耗大量的人力物力，精度也难以保证。

随着地理信息系统以及水文模型的发展，依靠水文模型提取流域水文信息已成为流域水文数据获取的重要途径［1］。

数字高程模型(DEM) 是描述地面高程空间分布的有序数值阵列，是目
前GIS 中表现3D 地形的主要2．5 维手段［2］。

Arc Hydro 是ArcGIS 中用于水资源管
理的地理加时间数据模型，可以支持水文信息系统的搭建，把空间和时间的水资源数据结
合起来，辅助
完成水文分析和模拟［3］。

基于DEM 和Arc Hydro 模型提取流域水文信息，目前逐
渐成为流域水文分析的有效手段。

1 研究区域概况
松花江流域是我国七大流域之一，位于我国东北地区(东经119°52'～132°31'，北
纬41°42'～51°38')。

松花江全长2 328 km，流域面积约为56．12 万km2，径流总量759 亿m3。

松花江流域是我国的重工业基地之一，也是重要的商品粮基地和林业生产基地。

从某种意义上来说，松花江虽然是黑龙江的支流，但其对东北地区的工农业生产、内河航运、人民生活等方面的经济和社会意义都超过黑龙江和东北其他河流［4］。

另一方面，
松花江流域建立了大量的化工、冶金等企业，对水体造成严重污染，2021 年更是发生了
硝基苯污染重大事件［5］。

因此本文选取松花江流域作为研究区域，对Arc Hydro 模型
在不同集水面积阈值下提取的水文信息进行精度分析，确定最佳集水面积阈值，提取流域
河网及子流域，以期为松花江流域水体污染和生态保护等研究提供数据参考。

2 数据来源与Arc Hydro 模型 2．1 数据来源
DEM 数据源: 由美国太空总署(NASA) 和国防部国家测绘局(NIMA) 联合测量的SRTM DEM。

2000年2月11日，美国发射的“奋进”号航天飞机上搭载SRTM 系统采集数据，在
中国境内可用的数据为3″(分辨率约90 m) ，数据的大地平面基准为WGS84，高程基准为EGM96，垂直误差小于16 m(90% 置信度) ［6］。

本文所用的SRTM DEM 数据是从中科院
国际科学数据服务平台下载拼接后依据松花江流域边界裁剪获得(图1) 。

松花江流域水系图: 松辽委提供的1∶ 25 万松花江流域水系矢量数据(图2) 。

2．2 Arc Hydro 模型
Arc Hydro 模型是把GIS 和水文地理领域知识相结合的水文地理数据模型，由美国ESRI 联合德克萨斯大学的水资源研究中心(CRWR) 共同搭建。

Arc Hydro 模型在ArcGIS
里存储水文要素的空间、属性及时间数据，用来描述流域的地貌和地形特征，通过水的运
动路径来反映各要素类间的关系［3］。

数据模型包括5 个部分: Network(河流网络)、Drainage(集水区域)、Channel (河道描述)、Hydrography(水文地形) 和TimeSeries (时间序列) ［7-8］。

Arc Hydro 模型结构信息如表1。

3 不同集水面积阈值下水文信息提取
Arc Hydro 模型能够实现流域内河网和子流域的自动提取，其主要过程包括: 数据的
预处理、流向分析、汇流累积量计算、河网提取以及子流域划分。

本文在此基础上增加精
度分析环节，将确定最佳集水面积阈值问题转化为定量的数值比较。

3．1 DEM 数据预处理
DEM 是用一组有序数值阵列形式表示地面高程的一种实体地面模型，大多是由离散的
高程点或者等高线数据进行插值获得，由于数据采集错误或者真实存在的凹陷等，DEM 数
据可能存在一些洼地［9］。

这种栅格高程值的突然性变化，会导致水流流向的计算错误。

因此在提取水文信息之前，需要对DEM 数据进行填洼预处理，处理过程利用Arc Hydro
模型中的Fill Sinks 工具，结果如图3。

3．2 流向分析
利用无洼的DEM 进行流向分析是水文信息提取的关键，Arc Hydro 模型运用单流向
D8 算法进行流向分析。

D8 算法具体过程如下: 首先计算DEM数据中每个栅格与其相邻的
8 个栅格之间的高程差，然后按坡度最陡原则确定该栅格的水流方向，即高程差为正值，
则为流出，若为负值，则为流入［10］。

流向分析过程利用Arc Hydro 模型中的Flow Direction工具，结果如图4。

3．3 汇流累积量计算
在得到水流方向之后，可利用水流方向数据来计算汇流累积量。

汇流累积量表示区域
内经过每个栅格的流水累积量。

汇流累积量计算方法为: 假设数字高程模型每个栅格点处
有一个单位的水量，根据上一步流向分析结果计算流过每个栅格的水量大小，从而得到该
区域的汇流累积量。

栅格的汇流累积量越大，表示能够流入的栅格数目越多，其汇流能力
越强，该栅格所代表的地形就很可能是河谷; 反之，汇流累积量为零的栅格，则可能代表
流域的分水岭［11］。

汇流累积量的计算利用Arc Hydro模型中的Flow Accumulation 工具，结果如图5。

3．4 河网提取
当某个栅格点上的汇流累积量超过了设定的集水面积阈值，则可以认为该点属于某一
水系，将附近的各个相似栅格点连接，即可形成一条河道。

设定的集水面积阈值大小决定
了河网提取的精度和详细程度，换句话说，参照不同尺度的流域水系图，要达到相同的提
取精度，设定的集水面积阈值一般不同。

通过DEM 提取的河网能否与流域实际情况相吻合，取决于DEM 能否反映流域地形的
实际情况及所采用的集水面积阈值的大小［10］。

经过反复验证，选取400 ～600 km2 的集水面积阈值作为研究范围。

通过调整集水面积阈值，提取不同集水面积阈值下的河网。

河网提取过程利用Arc Hydro 模型中的stream definition和stream segmentation 工具，结果如图6。

3．5 精度分析
不同集水面积阈值下提取的河网精度一般不同，通常情况下集水面积阈值设定的越小，提取的河网分叉越多，河流长度越长。

Alike Duplicate ImageFinder 是一款图像相似度
比对软件，可以将不同图像之间的相似性以数值形式定量表现出来，相似度越接近1，则
两个图像越相似。

基于Alike DuplicateImage Finder 分别将400、500 和600 km2 的集水面积阈值下提取的河网与实际河网进行图像相似度比对，相似度分别为1．14、0．98
和0．82，显然集水面积阈值为500 km2 下提取的河网最为精确。

3．6 子流域划分
子流域与地理中实际子流域的概念基本一致，其划分的前提是提取的流域河网与实际
河网相符。

Arc Hydro模型根据流向分析结果确定子流域的拓扑关系，其划分边界通过DEM 的数字地形分析来确定，按照生成的起始河道和流域水系，从最低一级的河道交叉点
沿分水线勾绘出起始子流域，如此类推到高一级河道直至全流域［12］。

子流域划分利用Arc Hydro 模型中的catchment grid delineation 工具，结果如图7。

根据精度分析结果，在500 km2 集水面积阈值下对松花江流域进行子流域划分，共提取了642 个子流域。

4 结果与讨论
利用SRTM DEM 数据和Arc Hydro 模型，在400、500 和600 km2 集水面积阈值下进
行河网提取。

经过与实际河网相似度比对，相似度分别为1．14、0．98 和0．82，确定
最佳集水面积阈值为500 km2，在500 km2 集水面积阈值下进行松花江流域子流域划分，
共提取了642 个子流域。

基于Arc Hydro 模型提取的松花江流域河网与实际河网相似度
极高，说明利用该模型进行流域水文信息提取是合理可行的。

感谢您的阅读，祝您生活愉快。