山谷线容易发生洪水?汇水范围怎么划?
“山谷线容易发生洪水?汇水范围怎么划?”这个问题,本质上是在问:地形中的水会往哪里流、哪些区域会把雨水汇集到同一条沟谷或河道里,以及这些汇水区域在 GIS 中如何可靠地划出来。对 GIS 学生、规划分析人员和水文制图人员来说,理解山谷线、汇流路径和汇水范围,比单纯点一个工具更重要。
本文以 DEM 数字高程模型为核心,说明为什么山谷线更容易形成洪水通道,并给出在 QGIS、ArcGIS Pro 或类似 GIS 软件中划分汇水范围的通用流程。你可以把它用于山洪风险初判、小流域边界提取、排水沟选线、道路涉水点分析等场景。

引言:山谷线为什么要和汇水范围一起看
山谷线通常是地表径流的汇集路径。降雨发生后,坡面水会沿地势较低的方向流动,并逐渐向沟谷、河槽、低洼地集中。因此,山谷线本身不等于洪水一定发生的位置,但它往往是洪水、泥石流、冲沟发育和积水风险分析中最需要重点检查的线状地形。
判断山谷线是否容易发生洪水,不能只看一条线的位置,还要看它上游有多大的汇水范围。汇水范围越大,表示有更多地表区域的雨水可能汇入同一出口;在相同降雨条件下,汇流量通常也越大。GIS 中常用 DEM 水文分析来回答这个问题。
背景:山谷线容易发生洪水的几个地形原因
在自然地形中,水往低处流是最基本的水文逻辑。山谷线之所以更容易成为洪水通道,主要有以下原因:
- 地形汇聚:山谷两侧坡面向谷底倾斜,降雨径流会向谷底集中。
- 上游汇水面积大:一条山谷线可能承接多个支沟、坡面和小流域的来水。
- 沟道输水能力有限:当降雨强度超过沟道排水能力时,容易形成漫溢或冲刷。
- 下垫面影响明显:裸地、硬化地面、陡坡区域会增加地表径流,减少下渗。
- 局部障碍改变水流:道路、涵洞、挡墙、堆土等可能阻断原有汇流路径,形成新的积水点。
因此,洪水风险分析不能只依赖等高线上的“沟谷形态”。更稳妥的做法是用 DEM 计算流向、汇流累积,再结合出水口划分汇水范围。
原理:汇水范围怎么划的核心逻辑
汇水范围,也常叫集水区或流域,是指所有能够通过地表径流汇入某个指定出口点的区域。在 GIS 中,汇水范围划分通常依赖 DEM 栅格数据。
常见计算逻辑如下:
- 准备 DEM:DEM 记录每个栅格单元的地面高程。
- 填洼:修正 DEM 中不合理的洼地,避免水流被错误困在单个像元或小坑中。
- 计算流向:判断每个像元的水应该流向相邻哪个像元。
- 计算汇流累积:统计上游有多少像元会流到当前像元。
- 提取山谷线:汇流累积值较高的像元通常对应沟谷、河道或潜在径流通道。
- 设置出水口:在山谷线、河道断面、涵洞口或分析断面上放置控制点。
- 划分汇水范围:反向追踪所有流向该出水口的上游像元,得到闭合区域。
这里最关键的是出水口位置。出水口放在不同位置,得到的汇水范围会完全不同。比如在沟谷上游放一点,划出的只是小汇水区;在下游村庄入口处放一点,可能包含整个上游山谷系统。
步骤:用 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 制图提供可靠基础。