山谷线容易发生洪水?汇水范围怎么划?

GIS基础理论
Dr.GIS
wowwwai GIS研习社 · 工具流程与项目排障

“山谷线容易发生洪水?汇水范围怎么划?”这个问题,本质上是在问:地形中的水会往哪里流、哪些区域会把雨水汇集到同一条沟谷或河道里,以及这些汇水区域在 GIS 中如何可靠地划出来。对 GIS 学生、规划分析人员和水文制图人员来说,理解山谷线、汇流路径和汇水范围,比单纯点一个工具更重要。

本文以 DEM 数字高程模型为核心,说明为什么山谷线更容易形成洪水通道,并给出在 QGIS、ArcGIS Pro 或类似 GIS 软件中划分汇水范围的通用流程。你可以把它用于山洪风险初判、小流域边界提取、排水沟选线、道路涉水点分析等场景。

山谷线容易发生洪水 汇水范围划分 DEM水文分析流程
基于 DEM 提取山谷线、确定出水口并划分汇水范围的典型 GIS 工作流。

引言:山谷线为什么要和汇水范围一起看

山谷线通常是地表径流的汇集路径。降雨发生后,坡面水会沿地势较低的方向流动,并逐渐向沟谷、河槽、低洼地集中。因此,山谷线本身不等于洪水一定发生的位置,但它往往是洪水、泥石流、冲沟发育和积水风险分析中最需要重点检查的线状地形。

判断山谷线是否容易发生洪水,不能只看一条线的位置,还要看它上游有多大的汇水范围。汇水范围越大,表示有更多地表区域的雨水可能汇入同一出口;在相同降雨条件下,汇流量通常也越大。GIS 中常用 DEM 水文分析来回答这个问题。

背景:山谷线容易发生洪水的几个地形原因

在自然地形中,水往低处流是最基本的水文逻辑。山谷线之所以更容易成为洪水通道,主要有以下原因:

  • 地形汇聚:山谷两侧坡面向谷底倾斜,降雨径流会向谷底集中。
  • 上游汇水面积大:一条山谷线可能承接多个支沟、坡面和小流域的来水。
  • 沟道输水能力有限:当降雨强度超过沟道排水能力时,容易形成漫溢或冲刷。
  • 下垫面影响明显:裸地、硬化地面、陡坡区域会增加地表径流,减少下渗。
  • 局部障碍改变水流:道路、涵洞、挡墙、堆土等可能阻断原有汇流路径,形成新的积水点。

因此,洪水风险分析不能只依赖等高线上的“沟谷形态”。更稳妥的做法是用 DEM 计算流向、汇流累积,再结合出水口划分汇水范围。

原理:汇水范围怎么划的核心逻辑

汇水范围,也常叫集水区或流域,是指所有能够通过地表径流汇入某个指定出口点的区域。在 GIS 中,汇水范围划分通常依赖 DEM 栅格数据。

常见计算逻辑如下:

  1. 准备 DEM:DEM 记录每个栅格单元的地面高程。
  2. 填洼:修正 DEM 中不合理的洼地,避免水流被错误困在单个像元或小坑中。
  3. 计算流向:判断每个像元的水应该流向相邻哪个像元。
  4. 计算汇流累积:统计上游有多少像元会流到当前像元。
  5. 提取山谷线:汇流累积值较高的像元通常对应沟谷、河道或潜在径流通道。
  6. 设置出水口:在山谷线、河道断面、涵洞口或分析断面上放置控制点。
  7. 划分汇水范围:反向追踪所有流向该出水口的上游像元,得到闭合区域。

这里最关键的是出水口位置。出水口放在不同位置,得到的汇水范围会完全不同。比如在沟谷上游放一点,划出的只是小汇水区;在下游村庄入口处放一点,可能包含整个上游山谷系统。

步骤:用 DEM 划分汇水范围的实操流程

步骤一:准备 DEM 和研究范围

首先准备覆盖研究区及其上游区域的 DEM。不要只裁剪到行政边界或项目红线,因为真实水流不会按行政界线停止。如果研究区边缘刚好在山谷上游,裁剪过小会导致汇水范围被截断。

建议检查以下内容:

  • DEM 是否完整覆盖上游地形。
  • 坐标系是否为投影坐标系,单位最好是米。
  • DEM 分辨率是否适合分析尺度,例如小流域或山洪沟分析不宜使用过粗数据。
  • 高程是否存在明显空洞、异常尖峰或条带噪声。

步骤二:统一坐标系和裁剪范围

如果 DEM 是经纬度坐标系,建议先投影到适合本区域的投影坐标系,例如 CGCS2000 高斯克吕格分带、UTM 分带或当地常用工程坐标系。这样后续面积、长度、坡度计算更可靠。

裁剪 DEM 时,应在研究区外保留足够缓冲区,尤其是上游方向。对于山谷线洪水分析,宁可多保留一些上游地形,也不要只留下沟谷下游的一小段。

步骤三:填洼处理

DEM 中经常存在局部洼地。洼地可能是真实湖泊、采坑,也可能是数据误差。如果不处理,流向计算时水会停在洼地中,导致汇流路径中断。

在 QGIS 中可以使用处理工具箱中的水文分析工具,例如 GRASS、SAGA 或 WhiteboxTools 提供的 Fill sinks、Fill depressions 类工具。在 ArcGIS Pro 中可以使用 Spatial Analyst 的 Fill 工具。

填洼并不是越强越好。过度填洼可能抹平真实地形,尤其是在城市、矿区、梯田和湖泊区域。对于真实水体或人工排水系统复杂区域,需要结合影像、河网和现场资料判断。

步骤四:计算流向

流向计算用于确定每个 DEM 像元的水流去向。常见算法是 D8,即每个像元只流向周围 8 个方向中坡降最大的一个邻近像元。它简单、稳定,适合很多基础流域划分任务。

在 ArcGIS Pro 中,一般使用 Flow Direction 工具。在 QGIS 中,可通过 GRASS r.watershed、SAGA Flow Accumulation 或 WhiteboxTools D8 Pointer 等工具实现。

如果研究区地形平缓,或水流可能出现分散流,D8 结果可能过于“单线化”。这时可以考虑多流向算法,但对于入门级汇水范围划分,D8 仍然是最常用的起点。

步骤五:计算汇流累积并提取山谷线

汇流累积表示每个像元上游有多少像元汇入它。数值越大,说明上游汇水面积越大。将汇流累积结果按阈值筛选,就可以提取潜在山谷线或河网。

例如,如果 DEM 分辨率为 10 米,一个像元面积为 100 平方米。若某像元汇流累积值为 10000,则其上游贡献面积约为 10000 × 100 = 1000000 平方米,即约 1 平方公里。

提取山谷线时,阈值没有固定答案。你需要根据地形尺度和分析目的调整:

  • 阈值太小:会提取出大量细碎沟线,图面杂乱。
  • 阈值太大:只保留主沟,可能漏掉小支沟和山洪沟。
  • 适合做法:先用多个阈值试算,再叠加影像、等高线、已有河网进行判断。

步骤六:确定出水口位置

出水口是划分汇水范围的控制点。常见出水口包括:

  • 沟谷下游与道路交汇处。
  • 涵洞、桥梁、排水口位置。
  • 村庄、厂区或建设用地上游边界处。
  • 山洪沟进入主河道的位置。
  • 需要计算来水面积的任意断面。

实际操作中,出水口点必须落在正确的水流路径上。如果点偏离山谷线几个像元,软件可能划出很小或错误的汇水范围。很多工具提供 Snap Pour Point 或类似功能,可以将出水口吸附到附近汇流累积值最大的像元上。

步骤七:划分汇水范围

完成流向栅格和出水口点后,就可以运行 Watershed、Basin、Catchment Area 等工具划分汇水范围。输出结果通常是栅格,再转换为矢量面用于制图、面积统计和叠加分析。

划分后建议做三项检查:

  • 汇水范围边界是否大致沿山脊线闭合。
  • 出水口是否位于汇水范围边界的最低出口处。
  • 汇水范围内的山谷线是否能够连续汇入出水口。

步骤八:计算汇水面积并评估洪水关注程度

将汇水范围转为面要素后,可以计算面积。面积单位建议使用平方米、公顷或平方公里,前提是数据已经处于合适的投影坐标系。

需要注意,汇水面积大并不等于一定发生洪水。洪水风险还与降雨强度、坡度、土壤入渗、植被覆盖、土地利用、沟道断面、排水设施等因素有关。但在缺少详细水文模型时,汇水范围和山谷线是非常重要的初判依据。

常见坑:汇水范围划错通常不是软件问题

坑一:DEM 裁剪范围太小

如果 DEM 只裁剪到项目区,软件只能在项目区内部寻找上游来水,结果必然偏小。山谷线分析尤其要保留上游地形,否则汇水范围会被人为截断。

坑二:出水口没有吸附到河道或山谷线

出水口偏离汇流路径是最常见错误之一。表现为输出汇水范围很小、形状奇怪,或者完全不包含真实上游沟谷。解决办法是叠加汇流累积图,把出水口吸附到高汇流累积像元上。

坑三:坐标系单位不对导致面积错误

如果数据仍是经纬度坐标系,直接计算面积会产生明显误差。汇水范围面积统计前,应使用投影坐标系,并确认单位是米。

坑四:汇流累积阈值随意设置

阈值决定山谷线提取的密度。不同 DEM 分辨率下,同一个阈值代表的实际汇水面积不同。设置阈值时,不要只看数值,要换算成实际面积,并结合影像和地形判断。

坑五:把山谷线等同于正式河道

DEM 提取的山谷线表示潜在汇流路径,不一定等于有常年水流的河道。它可能是季节性冲沟、干沟,也可能受道路和人工排水影响发生偏移。用于工程判断时,应结合遥感影像、实测河网和现场踏勘。

方法比较:QGIS、ArcGIS Pro 和命令行工具怎么选

方法 适合场景 优点 注意事项
QGIS 处理工具箱 学习、教学、小流域分析、开源工作流 免费,集成 GRASS、SAGA、WhiteboxTools 等水文工具 不同算法参数名称不同,需要理解流向和汇流累积原理
ArcGIS Pro Spatial Analyst 工程项目、规范化制图、已有 Esri 工作环境 工具链成熟,Fill、Flow Direction、Flow Accumulation、Watershed 流程清晰 需要相应扩展许可,数据管理要注意投影和像元大小
GRASS GIS r.watershed 较大区域、自动化批处理、水文地形分析 功能强,适合复杂水文处理 参数较多,初学者需要先做小范围测试
WhiteboxTools 自动化脚本、Python 工作流、批量处理 DEM 工具丰富,适合与 Python 结合 要注意工具版本、输入输出格式和无数据值设置

如果你刚开始学习,建议先用 QGIS 或 ArcGIS Pro 完整跑通一次流程。等理解了填洼、流向、汇流累积、出水口和汇水范围之间的关系,再考虑用脚本批处理。

检查清单:划分汇水范围前后要核对什么

在提交图件或分析结果前,可以按下面清单逐项检查。

  • DEM 是否覆盖出水口上游完整地形。
  • DEM 是否已经处理无数据值、异常高程和明显空洞。
  • 数据是否使用合适的投影坐标系。
  • 填洼参数是否合理,是否误改真实水体或大型洼地。
  • 流向栅格是否连续,没有异常断裂。
  • 汇流累积高值区是否与影像中的沟谷、河道大体一致。
  • 山谷线提取阈值是否换算过实际汇水面积。
  • 出水口是否吸附到正确的汇流路径上。
  • 汇水范围边界是否沿山脊线闭合。
  • 汇水面积计算是否使用米制投影坐标系。
  • 结果是否与已有河网、等高线、遥感影像或现场资料交叉验证。

FAQ:山谷线洪水和汇水范围划分常见问题

山谷线一定容易发生洪水吗?

不一定。山谷线是水流汇集的地形通道,所以更值得关注,但是否发生洪水还取决于降雨、汇水面积、坡度、土地利用、沟道断面和排水设施。GIS 分析可以帮助识别潜在风险区,不能替代完整的水文水力计算。

汇水范围怎么划才算准确?

准确的汇水范围需要满足几个条件:DEM 覆盖完整上游区域,流向计算合理,出水口位置正确,并且结果边界与山脊线、沟谷线、影像和已有水系基本一致。只要出水口放错,汇水范围就会明显偏差。

汇流累积阈值应该设多少?

没有统一数值。阈值应结合 DEM 分辨率和目标沟谷尺度确定。建议把阈值换算成实际上游面积,再通过多次试算选择最符合影像、等高线和已有河网的结果。

为什么我划出来的汇水范围特别小?

常见原因是出水口没有落在流向路径上,或者 DEM 裁剪范围太小。先检查出水口是否在汇流累积高值像元上,再检查 DEM 是否包含完整上游区域。必要时使用出水口吸附工具。

城市区域也能用 DEM 划汇水范围吗?

可以做初步地形分析,但城市区域的水流受雨水管网、道路边沟、下穿通道、泵站和建筑物影响很大。单纯 DEM 划分的汇水范围可能与真实排水分区不一致,需要结合管网资料和精细地形数据。

山谷线和分水岭有什么区别?

山谷线是水流汇集的位置,通常位于地形低处;分水岭是相邻汇水范围之间的边界,通常位于山脊或高地。划分汇水范围时,边界一般沿分水岭闭合,内部水流沿山谷线汇向出水口。

结论:先确定出水口,再谈山谷线洪水风险

判断山谷线容易不容易发生洪水,不能只看沟谷线本身,而要把它放在完整的汇水范围中理解。DEM 水文分析的核心流程是:准备 DEM、填洼、计算流向、计算汇流累积、识别山谷线、确定出水口、划分汇水范围、计算面积并核查结果。

实际工作中,最容易出错的是 DEM 范围太小、出水口位置不准、坐标系单位错误和汇流累积阈值随意设置。只要把这些环节控制好,汇水范围划分就能为山洪风险初判、工程选址、排水分析和 GIS 制图提供可靠基础。