ArcGIS水文分析+HEC-RAS水动力模拟：从DEM到洪水风险图完整链路

搞水文风险分析这么多年，我最大的感触是：单独用GIS做水文分析，容易停留在“纸上谈兵”的程度；单独用HEC-RAS做水力学模拟，又容易卡在“没有靠谱边界条件”的尴尬上。把ArcGIS的水文分析模块和HEC-RAS的水动力学模拟串成一条完整链路，才真正把洪水危险性评估从“画几根等值线”推进到“能落地的风险图”。

这套组合拳我前后跑了不下几十个项目，从山区小流域到城市内涝片区都试过。今天就把从DEM处理、河网提取、子流域划分，到HEC-RAS建模、参数率定、危险性分级，再到风险评估制图的完整流程，按我的实操经验捋一遍。这篇文章适合正在做洪水风险图、山洪灾害调查评价、城市内涝风险评估项目的同行，以及想从“只会点ArcGIS工具栏”进阶到“能独立跑通一个完整水文水动力项目”的新手。

先说清楚一个底层认知：ArcGIS水文分析解决的是“水从哪里来、往哪里去”的问题，HEC-RAS解决的是“水来了以后淹多深、流多快”的问题。前者输出的是汇流路径、集水面积、河道分级这些“骨架”，后者输出的是水面线、淹没范围、流速分布这些“血肉”。两者结合，才能回答洪水风险评价里最核心的三个问题：淹哪儿、淹多深、多危险。

1. 项目整体思路与方案选型

1.1 为什么把GIS水文分析和HEC-RAS搭在一起用

很多人问我，为什么不用HEC-HMS直接做水文模型，再用HEC-RAS做水力学？或者干脆全用HEC-RAS自带的地形处理功能？我做下来觉得，ArcGIS和HEC-RAS的搭配，在中小流域尺度上性价比最高。

ArcGIS的水文分析工具集（Hydrology）基于地表径流模型，原理很简单：水流永远沿着最陡坡向流动。基于这个朴素的假设，我们可以从DEM里提取出完整的河网、流域边界、汇流累积量，这些是后续一切分析的底子。这套工具跑得快、参数少、结果直观，特别适合做初筛和前期方案设计。

HEC-RAS则是真正的杀手锏。它是美国陆军工程兵团开发的开源水力学模型，行业认可度极高，支持一维/二维水动力模拟。一维模型把河道简化成断面序列，用能量方程和动量方程求解水面线；二维模型则用全二维浅水方程模拟漫滩水流。ArcGIS负责把地形“翻译”成模型能吃的几何数据，HEC-RAS负责把水文成果“翻译”成洪水风险指标。

我用一个生活化的类比：ArcGIS水文分析就像给整个流域做了一次“地形体检”，摸清了血管（河网）和器官（子流域）的分布；HEC-RAS则是在这些血管上做“血流动力学检测”，看血液（洪水）流过时哪里会溢出、溢出多深。

1.2 这套方案的适用边界和预期产出

这套组合不是万能的，我踩过坑之后明确它的适用条件：

• 流域面积：最适合中小流域（几十到几千平方公里）。太大流域需要分布式水文模型耦合，太小流域直接用地貌学方法更简单。
• 数据条件：必须有质量过得去的DEM（分辨率至少30米，最好10米或5米），否则提取的河网就是“鬼画符”。
• 洪水类型：适合暴雨洪水、溃坝洪水、河道漫溢型洪水。不适合风暴潮、内涝（除非做精细化二维耦合）。

一套标准流程跑完，你能得到这些产出：

我做完这套流程最深的体会：ArcGIS的“栅格思维”和HEC-RAS的“断面思维”是完全不同的两种语言，转录这一步是项目成败的关键。后面会详细说怎么处理这个转录。

2. 数据准备与前期处理

2.1 基础数据清单与质量要求

别一上来就开跑模型，数据这块偷懒的话坑你坑得死死的。列个清单给你：

• 数字高程模型（DEM）：核心数据，决定一切分析的下限。我一般要求至少30米分辨率，优选12.5米或5米（比如某公开的全球高精度地形数据）。山区用30米会漏掉很多细小沟道，平缓地区30米甚至会提取出乱七八糟的伪河网。
• 土地利用/覆被数据（LU/LC）：用于确定曼宁糙率系数（n值）。不同地表覆被糙率差别巨大：河道0.025-0.05，农田0.03-0.05，城镇不透水面0.013-0.02，林地0.05-0.15。
• 土壤类型数据：估算下渗率，用于产流计算。至少要有土壤质地分类（黏土、壤土、砂土等）。
• 雨量/水文数据：设计暴雨、实测流量过程线。拿不到实测流量就得用推理公式或区域水文手册推算。
• 河道断面数据：如果你有实测断面（测量成果或水深数据），精度直接上一个档次；没有的话就用DEM提取的断面，下游需要标注，精度会打折。
• 桥梁、涵洞、堰坝等涉水建筑物资料：这些对HEC-RAS计算有决定行影响，尤其是山区桥梁壅水，能抬高上游水位几十厘米。

用星号标注一下优先级：DEM > 降雨/水文数据 > 实测断面 > 土地利用 > 土壤数据 > 建筑物资料。项目时间紧的话，按这个顺序取舍。

2.2 DEM预处理里的隐藏细节

DEM拿到手不要直接扔进填洼（Fill）工具。我一开始也犯过这毛病，后来发现一堆伪沟道就这么被“填”出来的。

第一步：投影转换。先确认DEM的高程单位是米，坐标系统是投影坐标系（最好是UTM或者高斯-克吕格）。如果你的数据是经纬度（度分秒），算坡向、流向时会把距离单位搞混，流速失真。我习惯用投影坐标系，因为后面HEC-RAS的曼宁公式对单位极其敏感。

第二步：异常值处理。DEM里经常有负值、空心值（NODATA）、或者错误的“尖刺”高程点。用【栅格计算器】条件判断把这些异常值替换成邻域均值，或者在ArcGIS Pro里用【焦点统计】做一次平滑。这一步不做，流向栅格会出现一堆莫名其妙的局部洼点。

第三步：填洼（Fill）。填洼是所有水文分析的“起手式”。原理是不断地把局部低洼点抬高到和周围最低的溢出点一样高，这样水流才能连续地从上游流向下游。要注意填洼深度阈值：默认填到和周边最小高程齐平，但如果你的DEM有真实可信的洼地（比如喀斯特地貌、人工水库），这个步骤会把它们抹掉，结果就是河网偏离实际。我的做法是先做一次【条件函数】把已知水域掩膜掉再填洼，保护真实地貌。

第四步（进阶）：人工修正地形。如果你在平原地区，路堤、桥墩、涵洞这些人工地物在DEM里往往体现不出来。有经验的做法是用【欧几里得距离】和【栅格计算器】把道路中心线两侧的高程抬升，模拟路堤的挡水效应。这个操作对后续淹没范围影响极大，我在城市内涝项目中试过，做与不做路堤修正，淹没范围能差出20%。

3. ArcGIS水文分析实操：从DEM到流域参数

3.1 流向计算与河网提取的参数选择

填完洼，进入水文分析的核心主线。流程是：计算流向 → 计算汇流累积量 → 提取河网 → 河网分级 → 矢量化和栅格化。

流向计算用【流向（Flow Direction）】工具。它会给每个栅格单元分配一个方向编码（1、2、4、8、16、32、64、128），指向最陡下坡方向。ArcGIS默认用的是D8单流向算法，把水流全部导向一个邻域单元。算法简单、稳定，但对平坦地区和发散地形（如山前冲积扇）表现不佳。如果你的研究区很平缓，考虑用D-infinity（多流向）算法，不过HEC-RAS建模时河网还是要回到单流向提取的主河道上。我大部分项目用D8就足够了，关键是不要对流向结果反复“滤波”或平滑，让数据保持“原生性格”。

汇流累积量用【汇（Fill Accumulation）】工具，原理是把流向当作“水管网络”，统计每个单元上游有多少个单元的水流汇入。这个值就是“上游贡献面积”（单位：像元数）。有个经验：流量累积量和实际流量是近似正比的，所以它可以直接用来划定河道起点。

河网提取最简单的是【栅格计算器】条件判断：

con("flowacc" > 1000, 1)

阈值1000意味着“上游至少1000个栅格才算是河道”。这个阈值很关键，我整理过经验值：

确定阈值时，我通常会做一次“目视收敛”：把不同阈值提取的河网和实际水系（比如卫星影像上的河道）叠加对比，调到两者重合度最高的那个值。别指望一次成功，多试几个值对比，这个操作叫“阈值敏感性分析”，论文里能写，项目里更是保命。

3.2 子流域划分与特征参数计算

提取河网之后，下一个步骤是划分集水区。ArcGIS里可以用【分水岭（Watershed）】工具，把每个河网交汇点以上的区域都划成一个子流域。实际操作是先用【河流链接（Stream Link）】给每条河段分配唯一编号，然后再用【分水岭】基于流向栅格和链接栅格划分。

子流域数量是个平衡问题：划太细，后面HEC-RAS每个子流域都要输入洪峰流量，工作量爆炸；划太粗，落不了地图，风险评价的空间精度不够。我的经验是：根据项目面积和防洪需求，控制在10-30个子流域之间。比如300平方公里的流域，20个子流域比较合适，每个子流域面积10-20平方公里，既有统计意义又不至于面目模糊。

每个子流域都要计算这些特征参数：

• 面积（km²）——用栅格面积和投影面积换算
• 主沟长度（km）——河网中最长路径长度，影响汇流时间
• 流域平均坡度——计算【坡度（Slope）】后做【分区统计】
• 流域平均高程——同上，【分区统计】得到
• 形状系数——周长和面积的关系

这些参数汇总后，下一步要用它们推求设计洪水。

3.3 从ArcGIS导出HEC-RAS需要的几何文件

这一步是整个流程里最容易被低估、却最决定模型结果的环节。我管它叫“语言的翻译”，因为ArcGIS用的是矢量面（河岸、断面线），HEC-RAS用的是中心线+断面链（剖面）的语言。

具体操作流程：

1. 清理河网：把提取的河网矢量转成单线（【Feature to Line】），删除短小的支流毛刺，保留主河道。河道中心线必须连续、无分叉，因为HEC-RAS是沿水流方向逐断面计算的。
2. 生成断面线：在【HEC-GeoRAS】工具栏中用【断面线】工具，在中心线上每隔一段距离画横断面线。断面间距我一般采用河道宽度的1-2倍，或者100-500米。急弯、桥梁、汇流口附近加密（50-100米），平直河段可以放宽。每条断面线要完全覆盖可能淹没的范围，一般从一侧河漫滩到另一侧河漫滩，宁长勿短。
3. 赋予断面属性：用【HEC-GeoRAS】的【属性】工具给每条断面线赋予：「断面编号」（从上游往下游依次递增）、「河岸位置」（左岸/右岸）、「曼宁糙率」（按岸别和植被分开赋值）、「桥梁/涵洞标记」（有涉水建筑物时）。
4. 导出HEC-RAS几何文件：用【HEC-GeoRAS】菜单“Export to HEC-RAS”，生成一个几何数据文件（.g01格式），里面包含了河流中心线坐标、断面剖面的高程点序列、糙率分段、建筑物信息。
5. 在HEC-RAS中导入：打开HEC-RAS，新建项目，选择几何数据编辑器，从文件导入.g01格式的几何文件。

这步最容易出的三个问题：断面线方向不一致（必须都是左岸指向右岸，否则HEC-RAS报错）、断面与中心线交点不唯一（断面线和中心线“打结”了）、高程点间距太密（剖面上几百个点，HEC-RAS计算效率大崩）。前两个要在ArcGIS里修正，第三个可以用HEC-RAS的“断面点过滤”功能适当抽稀。

4. HEC-RAS建模与模拟实操

4.1 构建几何数据：断面、糙率、桥梁的处理

HEC-RAS的一维水动力模型基于三个方程：连续方程、动量方程、能量方程。人话翻译就是：水流量守恒、动量守恒、能量沿程损失。它沿河道按断面逐段求解，给定上游入流流量和下游边界条件，计算出每个断面的水位和流速。

在HEC-RAS几何数据编辑器里，你要做几件事：

断面数据检查。导入的断面如果发现“交叉断面”报错，多半是断面线方向或者轮廓点顺序有问题。HEC-RAS对断面的核心要求是：断面必须从左岸到右岸连续，且从上游看必须是逆时针方向（默认约定）。我一般会在导入后逐个断面快速浏览一遍，重点看有没有明显跳变的高程点。

糙率设置。曼宁糙率系数n值的选择，直接影响水位计算结果。用ArcGIS的水文分析提取了土地利用图之后，把每个断面的糙率按河槽和河漫滩分开赋值：河槽用主河道土质或植被对应的n值（常见0.03-0.045），河漫滩用与土地利用类型对应的n值（农田0.035、草地0.04、灌木0.05、城市不透水面0.015）。在HEC-RAS里打开“断面/河漫滩”窗口，分别输入左右河漫滩和主河槽的n值。糙率这参数敏感性非常高，推荐做一次“糙率敏感性分析”，看水位对n值的响应范围，报告里能说明稳健性。

桥梁/涵洞处理。在ArcGIS中标注了桥梁位置后，HEC-RAS几何编辑器里通过“桥梁/涵洞”工具添加。桥梁的几何参数包括：桥面最低点高程、桥孔宽度、桥墩形状和尺寸。别小看桥墩的形状系数，它直接进入局部水头损失计算，处理不好会让桥梁上下游水位差离谱。常见问题内存：桥孔净宽算错或者把桥墩忘掉，计算水面线会高出实际好几十厘米。

4.2 恒定流与非恒定流模拟的参数设置

HEC-RAS有两种模拟模式：恒定流（Steady Flow）和非恒定流（Unsteady Flow）。

恒定流是入门首选，也适合设计洪水评估。它假设流量不随时间变化，输入每个断面的设计洪峰流量，算出对应的沿程水面线。操作流程：

1. 编辑“稳态流数据”，输入各子流域入口的洪峰流量（就是你用推理公式或水文比拟法推求的设计洪水）。
2. 下游边界条件设置：常见选项有“正常水深”（需要输入河道比降）或“已知水位”。我用正常水深比较多，但河道比降不能直接取整体坡度，要取下游拟合断面的局部比降。
3. 计算模式选“水面线”，逐个断面运行。

非恒定流则更接近真实洪水过程。你需要给上游边界输入流量过程线（时间-流量数据），而不是一个固定流量。它能模拟洪水的涨落过程、洪峰传播时间、溃坝波的演进。设置要点：

• 计算时间步长一般取秒级（比如30-60秒），太小计算太慢，太大结果发散。
• 初始条件可以是“恒定流热身”或“给定初始水深”。
• 下游边界建议用“正常水深”或“水位流量关系曲线”。

我在做溃坝洪水分析时几乎只用非恒定流，因为它能还原溃口逐渐扩大、洪水波向下游推进的完整物理过程。但要注意：非恒定流计算参数多，调试成本高，新手建议先跑通恒定流，确认几何数据无误后再切换非恒定流。

4.3 结果导出与水面线提取

模拟完成后，HEC-RAS可以输出每个断面的水位、流速、过水面积、弗劳德数等成果。要导出到ArcGIS做后续风险分析：

1. 在HEC-RAS中打开“结果”菜单，选择“水位/流速”等参数，右键“导出为GIS格式”。
2. 用HEC-GeoRAS导入ArcGIS，把每个断面的计算结果（水位、流速）关联到断面线图层。
3. 用【栅格插值】工具（如TIN插值、克里金插值）把断面上的水位和流速插值成连续栅格。
4. 用【栅格计算器】计算淹没深度 = 水面高程栅格（插值结果）− 地形高程栅格（DEM）。注意：水面高程是通过断面水位插值出来的，如果断面间距太大或太少，插值结果和实际水面形态会有偏差，建议加密断面。
5. 用【栅格计算器】裁剪淹没深度>0.05米或>0.3米的区域为淹没范围（阈值根据项目精度要求定）。

插值那一步我要多叮嘱一句：不要直接用“反距离权重”就把水位插值完事，要先做一次合理性检查。因为在山区，河道水位差很大，断面间的插值必须考虑沿河方向的变化趋势，最好用“样条函数”或“带障碍的插值”。我见过一个项目拿IDW插值，结果在山坡上插出一片“空中的水”，就是因为没设置断线障碍。

5. 洪水危险性评价与风险综合评估

5.1 危险性等级划分：水深、流速、历时

有了淹没水深和流速栅格，就可以做危险性分级了。危险性（Hazard）指数常用水深×流速的组合来定义，行业里最常见的分级标准是：

具体分级阈值要结合当地的人口、建筑抗淹能力和行业标准调整，但整体逻辑不变：水深和流速的耦合决定危险性等级。

在ArcGIS里实现：用【栅格计算器】把水深和流速栅格重分类，再通过多条件判断合成一个危险性等级栅格。

con(("depth" < 0.5) & ("velocity" < 1), 1, con(("depth" < 1) & ("velocity" < 2), 2, con(("depth" < 2) & ("velocity" < 3), 3, 4)))

危险历时（淹没时间）一般只能靠非恒定流模拟的时间序列结果统计，做起来工作量大。如果项目不做，至少要把水深和流速两个指标做扎实。

5.2 易损性与暴露性分析

危险性只是“洪水本身有多猛”，风险评估还要看“淹到的东西值不值钱、多不多人”。这就要做易损性（Vulnerability）和暴露性（Exposure）分析。

暴露性回答“风险区内有什么”：把土地利用图、人口分布图、道路、学校、医院、重要基础设施等叠加到淹没范围上，统计各类承灾体的面积、数量、人口。在ArcGIS里用【叠加分析】（Intersect/Clip）和【分区统计】（Zonal Statistics）来实现，简单粗暴但极其有效。

易损性回答“承灾体有多脆弱”：不同用地类型对洪水的敏感性不同。农田淹几天可以恢复，但医院、变电站一淹就损失惨重。常见易损性系数表（0-1之间）：

易损性系数不是死的，你可以根据当地建筑结构、防汛能力做修正。用ArcGIS把土地利用栅格【重分类】成易损性系数栅格，得到易损性栅格。

5.3 风险矩阵与结果制图

风险（Risk）= 危险性（Hazard）× 易损性（Vulnerability）× 暴露量（Exposure）的乘积关系。最常用的做法是风险矩阵：

在ArcGIS里，把危险性栅格和易损性栅格用【栅格计算器】相乘，得到一个连续风险值栅格，再按自然间断点法（Jenks）或分位数法分成四级：低风险、中风险、高风险、极高风险。注意，相乘之前两个栅格的坐标系、像元大小、范围要完全一致，否则结果会出现错位。

最后制图时，我习惯的输出是：

• 底图：影像或地形图
• 主图层：危险性等级（红-橙-黄-绿渐变）
• 附加图层：风险等级（分层设色，叠加透明度）
• 必要元素：指北针、比例尺、图例、子流域编号、网格线

风险图的注记别偷懒，要写清数据来源、模拟工况（比如“100年一遇设计洪水”）、模型版本和参数。工程验收时这一页给评审人员看，权威性直接拉满。

6. 常见问题与排查技巧实录

跑这套流程这些年，我收集了一堆“翻车现场”的素材，整理成速查表：

有几个“独家避坑技巧”我单独拎出来说：

第一，DEM坐标系单位坑。印象最深的一次，我师弟用了经纬度坐标的DEM，流向计算时所有水平距离都按度算，结果汇流累积量在平原区爆炸式增长，提取的河网成了一张蜘蛛网。排查半天，发现就是投影没做。所以每次拿到数据第一步永远都是看坐标系。

第二，ArcGIS和HEC-RAS之间的“编号错位”。在ArcGIS里给断面编号如果是从1到100，但HEC-RAS导入后识别成从100逆序到1，整个河道的上下游方向就反了。我的习惯是在导出GeoRAS前用【排序工具】把断面按从上游到下游重新排序，并在HEC-RAS的几何编辑器里用“交叉断面”显示确认一下河流流向箭头方向。

第三，非恒定流模拟的“热启动”。如果从零水深开始跑非恒定流，前几个时间步往往会出现水位剧烈震荡，看上去就像“洪水倒流”。解决办法是先用恒定流计算出一个初始水面线，然后在非恒定流初始条件里选择“使用恒定流结果作为初始条件”。这个操作我每次必做，能省下大量调试时间。

第四，碎屑和局部洼地导致的“孤立淹没区”。做完淹没水深栅格后，你会发现有些淹没像元孤零零地出现在远离河道的地方，这是DEM噪声或插值误差造成的。不要急着下结论，先看看这些“斑块”是否和真实地形（比如局部低洼地、死水塘）对应。如果是伪斑块，用【焦点统计】的多数滤波抹掉；如果是真实地形，保留并单独分析。这条对我做城市内涝特别重要，有些“孤立淹没区”其实是下穿隧道或地下空间入口，反而需要重点警示。

第五，数据版本管理。跑项目的时候模型参数、DEM版本、河网阈值经常要反复试，如果没有版本管理的意识，后期想复现某个结果会非常痛苦。我现在每个方案的中间结果都按“日期+工况”命名存档，比如“20250101_100y_q100_g01”，模型参数单独记一个Excel表。这条看着朴素，几次项目验收前回看数据时救过我的命。

结尾

最后分享一个我反复和团队强调的工作习惯：一直做下去，直到你手里同时握着ArcGIS的水文分析图、HEC-RAS的水面线成果和最终风险等级图，再回头检查一遍“数据流链路”是否闭环。我遇到太多人做到中间某一步被工具卡住，结果东改西改，最后的成果图和最初DEM的坐标系都不一样。

实操中还有个小技巧值得试：把你的流域出口设一个“虚拟监测断面”，把HEC-RAS模拟出的水位过程线和实测（或水文手册发布的）洪水水位做对比。如果偏差在合理范围内，说明整条链路是通的，后续的风险评价结果才敢拿出去用。

这套ArcGIS+HEC-RAS的组合只覆盖了洪水风险分析的“主干道”，再往下走还能接HEC-HMS做全分布式水文模拟，或用HEC-RAS 2D做更精细的漫滩演进，甚至叠加FEMA的生命安全模型做人员疏散分析。先把这条主干道跑顺，后面拓展起来才有底气。遇到什么问题，欢迎同行有机会一起交流跑模型踩过的坑。