做洪水风险图这些年我最大的感触是HEC-RAS 和 ArcGIS 这对组合单拿出来任何一个都不稀奇但真正把它们耦合在一起跑通全流程才是让风险图从“画得出来”变成“算得出来、拿得出手”的关键。这篇就把我从数字高程模型处理、产汇流分析、水动力模拟到风险图出图的完整技术路线整理出来给正准备从头做洪水风险图的同行一个可以直接照抄的作业。先说清楚这个东西是干什么的。洪水风险图不是简单把水位线画在地图上它要回答三个问题哪些地方会被淹、淹到多深、淹了以后后果多严重。前两个靠水动力模型算第三个靠 GIS 做空间叠加分析两件事合在一起就绕不开 HEC-RAS 和 ArcGIS 的耦合。适合谁来参考刚接手洪水风险图项目的工程师、做流域防洪规划的技术人员以及想把手头零散数据整理成规范成果的研究生都能从这条流程里找到自己的切入点。1. 整体思路与耦合方案选型1.1 为什么必须让 HEC-RAS 和 ArcGIS 联动很多人一开始会问HEC-RAS 自己不是已经能算水位、出淹没范围了吗ArcGIS 在这里到底扮演什么角色这个问题得从洪水风险图的数据流说起。HEC-RAS 的核心是水动力计算它需要河道断面、糙率、流量、边界水位这些输入算完以后输出沿程水位、流速、水深。但它有个致命短板——几何数据的来源。河道的中心线怎么提、横断面在哪里切、断面地形从哪来HEC-RAS 自己基本干不了或者干得很吃力。这些恰恰是 ArcGIS 最擅长的事情。反过来ArcGIS 擅长画图和分析但它没有水动力计算能力如果只靠等高线或者经验法画一个淹没范围既没有物理依据也经不起审查。只有把两者耦合起来让 ArcGIS 准备模型、HEC-RAS 做计算、再回到 ArcGIS 出图分析整条链路才闭环。我见过不少项目在这上面吃过亏。有人只用了 HEC-RAS 自带的 RAS Mapper地形、断面全部在里面处理结果遇到复杂河网或者需要精细控制断面位置的时候完全使不上劲也有人只靠 ArcGIS 手工勾淹没范围评审的时候拿不出计算依据直接被专家问住。这两种偏科做法都不可取正确思路是“用 ArcGIS 管好模型输入用 HEC-RAS 保证水力计算精度再用 ArcGIS 完成成果表达”。1.2 两种耦合路线怎么选耦合路线目前主流有两种我分别踩过把各自的适用场景说清楚。第一种是经典路线ArcMap HEC-GeoRAS 扩展。HEC-GeoRAS 是 USACE 官方提供的工具条安装在 ArcMap 里专门用来做流域提取、河道中心线生成、横断面切割、地形数据提取最后直接生成 HEC-RAS 能读的几何文件。这条路线成熟稳定网上教程多国内大多数老工程师都是这么做的。第二种是 HEC-RAS 6.x 自带的 RAS Mapper 路线。RAS Mapper 内置了地形处理、河网生成、断面切割等功能理论上不需要 ArcGIS 也能完成建模。但实际用下来它在精细控制断面位置、复杂河网编辑、大批量地块处理上跟 ArcGIS 的编辑能力还是有差距尤其是做成果图的时候RAS Mapper 的出图效果远不如 ArcGIS 灵活。我做项目时用的是一条混合路线也是我想重点分享的在 ArcGIS 中完成 DEM 预处理、河网提取、流域边界划定在 ArcGIS 中手工编辑河道中心线、横断面切割线导出几何数据在 HEC-RAS 中完成水动力计算通过 RAS Mapper 导出淹没水深栅格回到 ArcGIS做淹没范围矢量化、风险等级划分、专题图出图。这套路线的优势在于把“模型前处理”和“成果后处理”两个最需要精细操作的环节放在 ArcGIS 里把“计算”这个黑盒子交给 HEC-RAS两边各干各擅长的事效率和可靠性都最高。1.3 版本搭配与安装建议版本搭配是很多新手一上来就卡住的地方。HEC-GeoRAS 这个工具条对版本非常敏感ArcMap 10.2 配 HEC-GeoRAS 10.2ArcMap 10.8 配 HEC-GeoRAS 10.8混搭基本装不上或者功能异常。我自己现在的环境是ArcMap 10.8 HEC-GeoRAS 10.8 HEC-RAS 6.3这套组合稳得很。如果你用的是 ArcGIS Pro那需要注意HEC-GeoRAS 的传统工具条在 Pro 里已经不支持了Pro 用户更推荐直接用 RAS Mapper或者通过 ArcGIS Pro 的 Interoperability 工具把 HEC-RAS 结果转进来。安装顺序也有讲究先装 ArcGIS再装 HEC-GeoRAS最后装 HEC-RAS。装完以后第一次打开 ArcMap要在自定义菜单里勾选 HEC-GeoRAS 工具条如果没有出现多半是扩展模块没启用到“自定义 扩展模块”里把 HEC-GeoRAS 勾上。还有个小坑安装路径里尽量不要带中文和空格HEC-GeoRAS 在中文路径下容易读不出工程文件这个我踩过不止一次。2. 前期数据准备ArcGIS 端的水文分析2.1 DEM 选型、裁剪与填洼处理整个流程从 DEM 开始。DEM 的精度直接决定后面所有成果的上限用什么数据要看项目范围和资金。大范围流域尺度的分析用公开的 30 米 SRTM 或者 ALOS 12.5 米数据就够了可以在地理空间数据云或者 USGS EarthExplorer 下载。重点防洪城区或者河谷地带最好用 5 米或者 1 米的激光雷达数据这种数据一般需要找测绘部门购买或者申请。拿到 DEM 以后第一件事不是急着做水文分析而是先裁剪出研究区范围。在 ArcGIS 里用“按掩膜提取”工具把 DEM 裁剪到你关心的流域边界或者行政区范围。这里有个细节裁剪的时候环境设置里的“捕捉栅格”一定要选上否则输出栅格的对齐方式会和原始 DEM 不一致后面叠加分析会出现半个像元的偏移。裁剪完以后做填洼。这个步骤很多人不理解其实是因为真实地形里有不少洼地雨水流到洼地就停了没法继续往下游流D8 算法就断了。填洼工具会把这些洼地填平保证水流能从头流到尾。填洼参数一般用默认就行但如果你的区域有真实的坑塘、采石场这类需要保留的地形填洼前要先把它们手工处理掉否则水流会失真。填洼完成以后我建议顺便做一次投影检查。如果原始 DEM 是经纬度坐标GCS_WGS_1984 或者 CGCS2000 经纬度一定要用“投影栅格”工具转到高斯投影或者 UTM 投影。经纬度坐标下做水文分析单位是度不是米流量累积和面积计算的结果全部是错的这个坑能让人排查到怀疑人生。2.2 河网提取与流域边界划定填洼结束以后水文分析的完整链路是流向 → 流量累积 → 河网提取 → 流域边界。流向用 D8 算法每个像元的水只流向周围 8 个像元里坡度最陡的那个ArcGIS 里对应“流向”工具。做完流向以后做“流量累积”这个工具统计的是每个像元汇集了多少个上游像元的水量值越大说明越可能是河道。接下来是河网提取。流量累积栅格需要设定一个阈值超过阈值的像元才被认为是河网。阈值怎么定小了河网密得跟毛细血管一样大了只剩干流。我的经验是先用一个比较大的阈值比如总像元数的 1% 左右看干流位置是否正确再逐步调小观察支流是否合理找到适合你研究区的那一档。ArcGIS 里把河网变成矢量线有三条路栅格计算器里写“Con(FlowAcc 阈值, 1)”再用“栅格河网矢量化”用“河网链接”工具先生成河段编码再矢量化直接用“Map Algebra”表达式配合“Stream to Feature”。我一般用第三种原因是河网矢量化工具直接输出的线往往有很多细碎分叉需要在编辑状态下手工整理。整理河网线有一个很重要的操作——加密折点。默认提取出来的河网线折点很稀疏直接拿去生成断面切割线切割线方向会乱。我会用“编辑工具 构造折线 平滑”或者“增密”工具把中心线变成平滑连续的折线。流域边界用“分水岭”工具输入流向栅格和出水口点就能把整个汇水区圈出来。注意出水口点要准确落在河道上否则分出来的流域会偏。2.3 河道中心线与横断面切割线的制作这是整个流程里最费功夫、也最考验耐心的一步。先说河道中心线。虽然可以用河网提取的结果直接转但实际项目中我强烈建议人工检查一遍。把提取的河网线和遥感影像叠加一段一段确认河道位置尤其是弯曲河段、分汊河段、人工裁弯取直的地方电脑提取经常不靠谱。ArcGIS 里画中心线的操作很简单开始编辑 → 选择中心线图层 → 构造工具选“折线” → 从上游往下游用鼠标点击放置折点 → 双击结束。但有几个经验要记住线段尽量平滑不要有尖锐的折角否则后续提取断面会出问题中心线方向必须统一从上游指向下游这个直接影响 HEC-RAS 里的桩号方向遇到江心洲、分汊河段要画两条中心线分别代表左右汊中间用“连接线”在分汊点衔接避免线段自相交和互相压盖压盖会导致断面切割时重复取点。断面切割线的生成是另一个技术点。在 HEC-GeoRAS 里有一个工具专门生成断面线它会在中心线上按设定的间距垂直生成切割线。但实际项目里这个自动生成的结果往往不满足要求因为河道的走向是弯曲的一个垂直方向切过去在高滩地或者堤防位置经常切得乱七八糟。我的做法是先自动生成一批再手工微调。切割线间距怎么设定平原城市河段我习惯取 200 到 300 米一条山区陡峭河段加密到 100 米桥梁、闸坝等水工建筑物上下游 50 米范围内还要单独加密。切割线的长度要覆盖整个可能淹没的范围一般取河道两侧各到堤防外 200 米以上宁可多切不能少切。还有一个关键操作——切割线必须和中心线垂直。这个听起来简单但在弯道处特别容易出问题。HEC-GeoRAS 生成断面线的时候是以中心线上的点为基准的如果中心线折点太稀疏生成的断面线角度就会偏。解决方法是在生成前用“增密”工具把中心线折点加密到 10 米一个这样每个位置的切线方向都准确断面线自然就垂直了。2.4 断面地形信息导出与几何文件生成断面切割线画好以后就要从 DEM 上提取地形信息了。这一步的核心逻辑是沿着每条切割线把 DEM 上的高程点按顺序取出来形成断面的地形剖面同时提取中心线与切割线交点的河床高程。HEC-GeoRAS 里的工具链是Stream Centerline Attributes生成河流代码和流向→ Cross-Section Cut Lines Attributes生成断面代码和上下游关系→ 然后用“Extract GIS Data”工具自动把断面地形数据提取到一个 GIS 图层里。这一步做完以后用 HEC-GeoRAS 的“Export GIS Data”功能可以直接导出 HEC-RAS 能识别的几何文件.g01 格式。导出的时候有两个选项一定要选对Horizontal Units 选 meters 或者你自己用的投影单位Vertical Units 选 meters千万别搞混了。导出成功以后在 HEC-RAS 里新建工程用“Import Geometry Data”把几何文件导入就能看到完整的河道断面了。我第一次做的时候兴奋了半天因为这一步意味着你从 GIS 到水力模型的桥终于通了。这里插一句糙率的事。糙率曼宁系数虽然是在 HEC-RAS 里设置的但它的取值依据跟土地利用类型、河道形态直接相关。在做断面切割的时候可以顺手把河道的糙率分区一并考虑进去把主槽、滩地、居民区分别标记出来这样后面在 HEC-RAS 里赋值就有了依据。3. HEC-RAS 水动力模拟核心环节3.1 几何文件导入与模型检查几何文件导入 HEC-RAS 以后第一件事不是急着跑计算而是做模型检查。打开几何数据编辑器逐条断面检查以下几个方面断面桩号顺序是否正确从上游到下游递增左右岸标记是否与实际河道一致左岸指面向下游时的左手边断面地形是否合理有没有异常跳点比如突然冒出一个 500 米高程的孤立点那多半是 DEM 没填洼干净或者提取的时候取到了噪声点断面之间的中心线连接是否连续有没有交叉、遗漏。HEC-RAS 里有个很实用的功能叫“River Schematic”就是那张河网示意图。几何文件导入后河网示意图上会显示河道和水流方向如果箭头指向不对说明中心线方向有问题赶紧回到 ArcGIS 里改。还有一个我踩过坑的地方——断面范围。HEC-RAS 计算的时候断面两端分别是左岸和右岸的终点如果断面范围没有覆盖到可能淹没的范围水位一高水就会从断面外侧漏出去算出来的淹没范围偏小。这种问题在界面上还不容易发现要到“Plot Cross Section”里看水位线是否越过了断面端点。3.2 设计洪水过程与边界条件构建HEC-RAS 支持稳态Steady Flow和非稳态Unsteady Flow两种计算模式做洪水风险图两种都会用到。稳态计算相对简单适合推求某一频率洪水比如 20 年一遇、50 年一遇、100 年一遇的最高水位和最大淹没范围。只需要给定上游边界流量和下游边界水位模型就会沿河算出每个断面的水面线。非稳态计算要复杂一些需要给定流量过程线时间序列适合做洪水演进、淹没过程分析成果包括不同时刻的淹没范围、最大淹没水深、淹没历时。做风险图一般建议至少用非稳态算一遍因为很多地区的淹没风险不是由洪峰流量单独决定的洪水的涨落过程、下游顶托都对淹没范围有影响。上游流量怎么定这个要靠水文计算不是 GIS 或水力模型能解决的。常用的方法是查本地区暴雨洪水手册用推理公式或者瞬时单位线法推求设计洪峰流量和洪水过程线。如果项目区附近有水文站直接用实测洪水资料进行频率分析用 P-III 型频率曲线求各种重现期的设计洪水这是最可靠的办法。下游边界条件的选择也有讲究。感潮河段或者靠近大江大河干流的支流河口一般用给定水位或水位流量关系如果下游没有控制断面用“正常水深”加上一个能坡一般取河床坡降 0.001 到 0.005也能凑合。但能坡取多少一定要看实际的河道纵坡取大了水位偏低取小了水位偏高两个方向的偏差在评审会上都容易被质疑。3.3 桥梁涵洞在水动力模型里的处理河道上的桥梁、涵洞、堰坝对洪水水位的影响非常明显忽略它们算出来的水位可能低一到两米。我在项目里遇到过最夸张的情况一座老桥的桥孔过流能力不足壅水达到 1.5 米直接让上游一片原本不淹的农田变成了淹没区。HEC-RAS 里处理桥梁有两种常用方式第一种是直接用 Bridge 工具建模。需要输入桥面的顶部高程、桥孔宽度、桥墩尺寸和数量、桥面底部高程等信息。HEC-RAS 会根据水位高低自动判断是桥孔过流还是桥面溢流计算相对精细缺点是桥梁构造数据不好拿老桥根本没有设计图纸。第二种是用 Inline Structure线状建筑物模拟。适用于堰坝、水闸等规则建筑物设定堰顶高程、宽度、流量系数就行操作比 Bridge 简单很多。实际项目里我的建议是重点桥梁必须用 Bridge 工具建普通小桥如果资料不全可以简化成 Inline Structure甚至在有足够数据支撑的情况下忽略。但要注意凡是做了简化的桥梁都要在技术报告里写清楚评审专家问起来要能说出依据。涵洞的处理思路类似HEC-RAS 里有 Culvert 工具需要输入涵洞类型圆管涵、箱涵、进口形式、尺寸、进出口高程、糙率等参数。涵洞计算里最容易出错的是淹没出流和自由出流的状态判断软件会自动处理但你要确认输入的进出口高程和实际地形一致。3.4 模拟计算、结果判定与导出几何、边界、建筑物都设置好以后就可以开始计算了。稳态计算直接运行 Steady Flow Analysis一般几分钟就能算完。算完以后一定要看“Errors, Warnings, Notes”窗口很多新手忽略这个结果模型算了也出图了但水位是错的还不知道错在哪。非稳态计算要重点设置计算时间步长和输出间隔。时间步长一般取 1 分钟到 5 分钟地形变化剧烈的山区河道取小值平原河道取大值。输出间隔影响结果文件的大小做风险图需要最大淹没水深可以设成 1 小时输出一次计算完用“Max Depth”来出图。计算完以后第一个要检查的是水位合理性。把算出来的水面线和实际调查的洪水痕迹点对比如果出入大于 0.5 米就要回头检查糙率、断面、边界条件。第二个要检查的是纵剖面线是否光滑如果出现锯齿状跳动说明有断面局部壅水或者糙率突变需要调整。结果确认没问题以后就用 RAS Mapper 导出成果。RAS Mapper 里可以加载地形、显示水位、水深、流速分布操作很直观。导出的时候选择“Depth”栅格范围选整个模型范围分辨率尽量和原始 DEM 一致这样回到 ArcGIS 里叠加分析最省心。这一节最后提一个非常关键的细节导出栅格前要在 RAS Mapper 里把坐标系设置正确。HEC-RAS 6.x 默认以地形 DEM 的坐标系为基准但有时候你导入的地形没有正确带投影信息导出的栅格就是个无坐标系的裸数据回到 ArcGIS 里还得手动定义投影一不小心定义错了就全部偏移。4. ArcGIS 端洪水风险图制作与出图4.1 模拟结果栅格导入与坐标统一从 RAS Mapper 导出的水深栅格回到 ArcGIS 里第一步是检查坐标系。我一般会先在图层属性里看栅格的坐标系如果是 Unknown直接右键“定义投影”选成和 DEM 一致的坐标系。这里必须说清楚一个概念定义投影Define Projection和投影转换Project Raster是完全不同的操作——定义投影只是告诉软件数据是哪个坐标系不会改变坐标值投影转换才是真正的坐标变换。很多人搞混导致数据跟底图对不上。如果栅格坐标系正确但和底图对不上那就需要做投影转换了。例如你导出的栅格是 CGCS2000_3_Degree_GK_CM_114E而底图是 WGS_1984_UTM_Zone_49N就要用“投影栅格”工具转换过去。栅格导入以后第一眼看到的往往是黑乎乎一片因为水深栅格的范围里没有淹没的区域是 NoData显示为黑色。这个不是数据坏了只是显示设置问题。在符号系统里把 NoData 显示为白色或者透明就能看到清晰的淹没范围了。4.2 淹没范围矢量化与水深分级水深栅格看清楚以后做两件事。第一把淹没范围提取出来。用“条件函数”工具Con把大于 0 的像元赋值 1小于等于 0 的赋值 NoData这样就能得到一个二分栅格再用“栅格转面”工具转成矢量面。转出来的面可能有大量细碎小图斑用“消除”工具按面积阈值合并掉一般小于 1000 平方米的碎图斑直接消除地图整洁度会大幅提升。第二做水深分级。洪水风险图不是只看一个“淹了还是没淹”关键要看淹多深。水深分级国内各行业规范略有差异但大致可以参考水深区间风险等级适用说明0 0.5 米低风险一般影响步行和车辆通行0.5 1.5 米中风险影响一层建筑人员需转移1.5 3 米高风险低层建筑进水救援困难大于 3 米极高风险严重威胁生命财产安全操作上用“重分类”工具把水深栅格按上述区间分成 1、2、3、4 级然后“栅格转面”得到分级面。或者用“栅格计算器”不断 Con 嵌套也行但层级多了容易出错我推荐先重分类再转换逻辑更清晰。转成矢量面以后每个面要素携带一个水深等级属性用不同颜色渲染区分这就是洪水风险图的核心图层了。4.3 风险等级划分与空间叠加分析单纯的淹没水深还不等于风险真正完整的洪水风险图还要考虑承灾体的脆弱性。ArcGIS 里做这一步用的是空间叠加分析。你需要准备土地利用现状、人口分布、道路网络、居民点等数据然后做以下几种典型分析淹没范围内人口统计用“面转栅格”把人口分布数据转成栅格再用“分区统计”Zonal Statistics按水深等级分区统计人口数量淹没范围内土地利用损失把淹没范围面和土地利用面做“相交”Intersect按地类统计淹没面积转移道路可达性把淹没水深大于 0.5 米的区域叠加到路网上标注出不可通行路段为防汛转移预案提供依据。这些分析的产出比单纯一张淹没水深图要实用得多。我做过的几个项目最后评审专家问得最多的就是“你这个范围里大概影响多少人”“有没有妨碍转移的路段”如果没有提前做叠加分析现场根本答不上来。风险等级图可以在水深分级的基础上再叠加人口密度和土地利用的权重得到一个综合风险等级——高风险区、中风险区、低风险区。具体权重可以根据当地实际情况调整每个区域标准不完全一样做之前先和项目负责人确认清楚别自己想当然。4.4 专题图整饰与统计表输出风险图最终要变成一张规范的成果图ArcGIS 的制图功能这时候派上用场。进入布局视图设置纸张大小一般用 A3 或者 A0 出图插入以下要素指北针、比例尺、图例、标题栏、坐标网。字体建议用宋体或者黑体地图要素的图层顺序按“影像底图 → 面状风险图层 → 线状道路水系 → 点状居民点”排列。制图里面有几个容易忽略的细节图例里的水深分级色带要跟正文描述一致别换个配色就完全对不上了比例尺设置要跟当前地图范围的投影单位一致用“动态比例尺”模式缩放的时候比例尺自动更新标题栏里要写清楚洪水频率如“100 年一遇洪水风险图”和坐标系信息风险图有多个洪水频率的版本出图前一定要在文件名和标题里写清楚不然后期整理报告的时候分不清哪个是哪个。统计表输出可以用“属性表导出”功能把淹没面积、影响人口、影响道路长度等指标整理成 Excel 表格作为风险报告的附表。导出的时候字段编码要是中文乱码用“字段计算器”把字段名先改成英文导出完再映射成中文可以绕开乱码问题。做完这些一套标准的洪水风险图成果就算齐了风险等级图、水深分布图、淹没范围图、统计报表一个都不少。5. 常见问题与排查技巧实录5.1 河网提取断裂或流向混乱这是 ArcGIS 水文分析里最高频的问题之一症状是提取的河网断断续续支流接不上干流或者流向箭头乱指。排查思路按顺序来检查 DEM 是否填洼了。没填洼直接做流向是最常见的断裂原因检查流向栅格的赋值范围D8 流向值只有 1、2、4、8、16、32、64、128 这 8 个值如果出现 0 或者负数说明 DEM 有 NoData 区域先做“填洼”或“条件函数”补值检查河网提取阈值阈值太大河网会断在细支流上调小阈值试试检查 DEM 投影单位经纬度坐标做流量累积会出现大面积异常值必须改成投影坐标。5.2 断面切割线方向不统一断面线在 HEC-GeoRAS 里每个断面都有左右岸属性方向不一致意味着左右岸属性错乱导入 HEC-RAS 以后模型会报错。这个问题的根源几乎都在中心线。中心线折点方向不统一——上半段上游到下游下半段反向生成断面的时候方向就全乱了。解决方法是画中心线之前用“编辑工具”里的“反向”功能统一线的方向确认所有中心线都是从上游指到下游。还有一个检查技巧在 ArcGIS 里给断面线加上箭头符号运行“符号系统 → 高级 → 箭头”看所有箭头是否指向河流一侧。如果有几根箭头指向反了单独选中“反向”那根线就可以了。5.3 HEC-RAS 迭代不收敛怎么调非稳态计算最常见的问题是迭代不收敛报错信息五花八门但本质都是水位或流量在相邻时间步之间振荡。我的排查顺序是缩短计算步长从 5 分钟改成 1 分钟检查断面间距如果相邻断面距离超过 500 米在 HEC-RAS 中插入中间断面检查下游边界下游水位给得太低或者能坡取得太陡容易导致计算发散把计算选项里的θ权重从默认值调到 1全隐式稳定性显著提升如果还不行退回检查几何文件——很可能是断面范围不足水位顶到断面端点导致水量泄漏。5.4 栅格导入 ArcGIS 背景全黑或位置偏移水深栅格导入以后全黑先别慌。右键图层 → 符号系统 → 拉伸换一个色带NoData 设成透明绝大多数情况就好了。如果还是黑的用“属性”查看栅格最小值和最大值有时候栅格是浮点型范围显示为 0 到 2.147e9其实是某个异常的像元值在“符号系统”里把拉伸范围手动改小就行。位置偏移的问题90% 是坐标系不统一。检查栅格和底图的坐标系不一致就投影转换千万别手动平移硬凑。5.5 中文路径、坐标系混乱引发的麻烦事这是我最想强调的最后一件事。ArcGIS 和 HEC-RAS 的中文路径兼容性都一般尤其是 HEC-GeoRAS 导出几何文件的环节路径里一有中文就容易报“无法读取数据”或者生成空文件。我的习惯是给每个项目单独建一个纯英文路径的工程文件夹比如 D:\FloodRisk\Project_Zhanghe所有中间数据、模型文件都放里面从源头避免这种低级问题。坐标系混乱我前面也提过几次。定义投影和投影转换是两个工具很多人随手用“定义投影”强行把矢量数据定义到另一个坐标系结果位置不变但数值全错。记住一句话不知道数据原本用什么坐标系先打开属性看能查到最好查不到问提供数据的人千万别瞎猜。最后说几句体会做洪水风险图这几年我最深的体会是模型本身不难学难的是把数据管清楚、把每一步的原理想明白。曾经为了一个断面的高程对不上反复折腾了两天才发现是中心线折点方向导致切割线取错位置这类问题工具书里不会写全靠经验积累。如果你准备上手这类项目我的建议是不要一上来就追求非稳态、高精度先从一条小流域的稳态计算跑通全流程把 ArcGIS 和 HEC-RAS 的数据流打通再逐步加桥梁、加非稳态、加风险分析。工具永远在更新但“GIS 准备数据 — 水动力计算 — GIS 表达成果”这条主线不会变。最后分享一个小技巧每次做完整套流程以后把中间的所有关键参数DEM 来源、填洼阈值、河网阈值、糙率取值、边界条件整理成一张配置表存到项目文件夹里。这不是给甲方看的是给三个月后的自己看的——到那时候你绝对会感谢当初记了这张表。 HEC-RAS与ArcGIS耦合：从DEM到洪水风险图的完整技术路线
做洪水风险图这些年我最大的感触是HEC-RAS 和 ArcGIS 这对组合单拿出来任何一个都不稀奇但真正把它们耦合在一起跑通全流程才是让风险图从“画得出来”变成“算得出来、拿得出手”的关键。这篇就把我从数字高程模型处理、产汇流分析、水动力模拟到风险图出图的完整技术路线整理出来给正准备从头做洪水风险图的同行一个可以直接照抄的作业。先说清楚这个东西是干什么的。洪水风险图不是简单把水位线画在地图上它要回答三个问题哪些地方会被淹、淹到多深、淹了以后后果多严重。前两个靠水动力模型算第三个靠 GIS 做空间叠加分析两件事合在一起就绕不开 HEC-RAS 和 ArcGIS 的耦合。适合谁来参考刚接手洪水风险图项目的工程师、做流域防洪规划的技术人员以及想把手头零散数据整理成规范成果的研究生都能从这条流程里找到自己的切入点。1. 整体思路与耦合方案选型1.1 为什么必须让 HEC-RAS 和 ArcGIS 联动很多人一开始会问HEC-RAS 自己不是已经能算水位、出淹没范围了吗ArcGIS 在这里到底扮演什么角色这个问题得从洪水风险图的数据流说起。HEC-RAS 的核心是水动力计算它需要河道断面、糙率、流量、边界水位这些输入算完以后输出沿程水位、流速、水深。但它有个致命短板——几何数据的来源。河道的中心线怎么提、横断面在哪里切、断面地形从哪来HEC-RAS 自己基本干不了或者干得很吃力。这些恰恰是 ArcGIS 最擅长的事情。反过来ArcGIS 擅长画图和分析但它没有水动力计算能力如果只靠等高线或者经验法画一个淹没范围既没有物理依据也经不起审查。只有把两者耦合起来让 ArcGIS 准备模型、HEC-RAS 做计算、再回到 ArcGIS 出图分析整条链路才闭环。我见过不少项目在这上面吃过亏。有人只用了 HEC-RAS 自带的 RAS Mapper地形、断面全部在里面处理结果遇到复杂河网或者需要精细控制断面位置的时候完全使不上劲也有人只靠 ArcGIS 手工勾淹没范围评审的时候拿不出计算依据直接被专家问住。这两种偏科做法都不可取正确思路是“用 ArcGIS 管好模型输入用 HEC-RAS 保证水力计算精度再用 ArcGIS 完成成果表达”。1.2 两种耦合路线怎么选耦合路线目前主流有两种我分别踩过把各自的适用场景说清楚。第一种是经典路线ArcMap HEC-GeoRAS 扩展。HEC-GeoRAS 是 USACE 官方提供的工具条安装在 ArcMap 里专门用来做流域提取、河道中心线生成、横断面切割、地形数据提取最后直接生成 HEC-RAS 能读的几何文件。这条路线成熟稳定网上教程多国内大多数老工程师都是这么做的。第二种是 HEC-RAS 6.x 自带的 RAS Mapper 路线。RAS Mapper 内置了地形处理、河网生成、断面切割等功能理论上不需要 ArcGIS 也能完成建模。但实际用下来它在精细控制断面位置、复杂河网编辑、大批量地块处理上跟 ArcGIS 的编辑能力还是有差距尤其是做成果图的时候RAS Mapper 的出图效果远不如 ArcGIS 灵活。我做项目时用的是一条混合路线也是我想重点分享的在 ArcGIS 中完成 DEM 预处理、河网提取、流域边界划定在 ArcGIS 中手工编辑河道中心线、横断面切割线导出几何数据在 HEC-RAS 中完成水动力计算通过 RAS Mapper 导出淹没水深栅格回到 ArcGIS做淹没范围矢量化、风险等级划分、专题图出图。这套路线的优势在于把“模型前处理”和“成果后处理”两个最需要精细操作的环节放在 ArcGIS 里把“计算”这个黑盒子交给 HEC-RAS两边各干各擅长的事效率和可靠性都最高。1.3 版本搭配与安装建议版本搭配是很多新手一上来就卡住的地方。HEC-GeoRAS 这个工具条对版本非常敏感ArcMap 10.2 配 HEC-GeoRAS 10.2ArcMap 10.8 配 HEC-GeoRAS 10.8混搭基本装不上或者功能异常。我自己现在的环境是ArcMap 10.8 HEC-GeoRAS 10.8 HEC-RAS 6.3这套组合稳得很。如果你用的是 ArcGIS Pro那需要注意HEC-GeoRAS 的传统工具条在 Pro 里已经不支持了Pro 用户更推荐直接用 RAS Mapper或者通过 ArcGIS Pro 的 Interoperability 工具把 HEC-RAS 结果转进来。安装顺序也有讲究先装 ArcGIS再装 HEC-GeoRAS最后装 HEC-RAS。装完以后第一次打开 ArcMap要在自定义菜单里勾选 HEC-GeoRAS 工具条如果没有出现多半是扩展模块没启用到“自定义 扩展模块”里把 HEC-GeoRAS 勾上。还有个小坑安装路径里尽量不要带中文和空格HEC-GeoRAS 在中文路径下容易读不出工程文件这个我踩过不止一次。2. 前期数据准备ArcGIS 端的水文分析2.1 DEM 选型、裁剪与填洼处理整个流程从 DEM 开始。DEM 的精度直接决定后面所有成果的上限用什么数据要看项目范围和资金。大范围流域尺度的分析用公开的 30 米 SRTM 或者 ALOS 12.5 米数据就够了可以在地理空间数据云或者 USGS EarthExplorer 下载。重点防洪城区或者河谷地带最好用 5 米或者 1 米的激光雷达数据这种数据一般需要找测绘部门购买或者申请。拿到 DEM 以后第一件事不是急着做水文分析而是先裁剪出研究区范围。在 ArcGIS 里用“按掩膜提取”工具把 DEM 裁剪到你关心的流域边界或者行政区范围。这里有个细节裁剪的时候环境设置里的“捕捉栅格”一定要选上否则输出栅格的对齐方式会和原始 DEM 不一致后面叠加分析会出现半个像元的偏移。裁剪完以后做填洼。这个步骤很多人不理解其实是因为真实地形里有不少洼地雨水流到洼地就停了没法继续往下游流D8 算法就断了。填洼工具会把这些洼地填平保证水流能从头流到尾。填洼参数一般用默认就行但如果你的区域有真实的坑塘、采石场这类需要保留的地形填洼前要先把它们手工处理掉否则水流会失真。填洼完成以后我建议顺便做一次投影检查。如果原始 DEM 是经纬度坐标GCS_WGS_1984 或者 CGCS2000 经纬度一定要用“投影栅格”工具转到高斯投影或者 UTM 投影。经纬度坐标下做水文分析单位是度不是米流量累积和面积计算的结果全部是错的这个坑能让人排查到怀疑人生。2.2 河网提取与流域边界划定填洼结束以后水文分析的完整链路是流向 → 流量累积 → 河网提取 → 流域边界。流向用 D8 算法每个像元的水只流向周围 8 个像元里坡度最陡的那个ArcGIS 里对应“流向”工具。做完流向以后做“流量累积”这个工具统计的是每个像元汇集了多少个上游像元的水量值越大说明越可能是河道。接下来是河网提取。流量累积栅格需要设定一个阈值超过阈值的像元才被认为是河网。阈值怎么定小了河网密得跟毛细血管一样大了只剩干流。我的经验是先用一个比较大的阈值比如总像元数的 1% 左右看干流位置是否正确再逐步调小观察支流是否合理找到适合你研究区的那一档。ArcGIS 里把河网变成矢量线有三条路栅格计算器里写“Con(FlowAcc 阈值, 1)”再用“栅格河网矢量化”用“河网链接”工具先生成河段编码再矢量化直接用“Map Algebra”表达式配合“Stream to Feature”。我一般用第三种原因是河网矢量化工具直接输出的线往往有很多细碎分叉需要在编辑状态下手工整理。整理河网线有一个很重要的操作——加密折点。默认提取出来的河网线折点很稀疏直接拿去生成断面切割线切割线方向会乱。我会用“编辑工具 构造折线 平滑”或者“增密”工具把中心线变成平滑连续的折线。流域边界用“分水岭”工具输入流向栅格和出水口点就能把整个汇水区圈出来。注意出水口点要准确落在河道上否则分出来的流域会偏。2.3 河道中心线与横断面切割线的制作这是整个流程里最费功夫、也最考验耐心的一步。先说河道中心线。虽然可以用河网提取的结果直接转但实际项目中我强烈建议人工检查一遍。把提取的河网线和遥感影像叠加一段一段确认河道位置尤其是弯曲河段、分汊河段、人工裁弯取直的地方电脑提取经常不靠谱。ArcGIS 里画中心线的操作很简单开始编辑 → 选择中心线图层 → 构造工具选“折线” → 从上游往下游用鼠标点击放置折点 → 双击结束。但有几个经验要记住线段尽量平滑不要有尖锐的折角否则后续提取断面会出问题中心线方向必须统一从上游指向下游这个直接影响 HEC-RAS 里的桩号方向遇到江心洲、分汊河段要画两条中心线分别代表左右汊中间用“连接线”在分汊点衔接避免线段自相交和互相压盖压盖会导致断面切割时重复取点。断面切割线的生成是另一个技术点。在 HEC-GeoRAS 里有一个工具专门生成断面线它会在中心线上按设定的间距垂直生成切割线。但实际项目里这个自动生成的结果往往不满足要求因为河道的走向是弯曲的一个垂直方向切过去在高滩地或者堤防位置经常切得乱七八糟。我的做法是先自动生成一批再手工微调。切割线间距怎么设定平原城市河段我习惯取 200 到 300 米一条山区陡峭河段加密到 100 米桥梁、闸坝等水工建筑物上下游 50 米范围内还要单独加密。切割线的长度要覆盖整个可能淹没的范围一般取河道两侧各到堤防外 200 米以上宁可多切不能少切。还有一个关键操作——切割线必须和中心线垂直。这个听起来简单但在弯道处特别容易出问题。HEC-GeoRAS 生成断面线的时候是以中心线上的点为基准的如果中心线折点太稀疏生成的断面线角度就会偏。解决方法是在生成前用“增密”工具把中心线折点加密到 10 米一个这样每个位置的切线方向都准确断面线自然就垂直了。2.4 断面地形信息导出与几何文件生成断面切割线画好以后就要从 DEM 上提取地形信息了。这一步的核心逻辑是沿着每条切割线把 DEM 上的高程点按顺序取出来形成断面的地形剖面同时提取中心线与切割线交点的河床高程。HEC-GeoRAS 里的工具链是Stream Centerline Attributes生成河流代码和流向→ Cross-Section Cut Lines Attributes生成断面代码和上下游关系→ 然后用“Extract GIS Data”工具自动把断面地形数据提取到一个 GIS 图层里。这一步做完以后用 HEC-GeoRAS 的“Export GIS Data”功能可以直接导出 HEC-RAS 能识别的几何文件.g01 格式。导出的时候有两个选项一定要选对Horizontal Units 选 meters 或者你自己用的投影单位Vertical Units 选 meters千万别搞混了。导出成功以后在 HEC-RAS 里新建工程用“Import Geometry Data”把几何文件导入就能看到完整的河道断面了。我第一次做的时候兴奋了半天因为这一步意味着你从 GIS 到水力模型的桥终于通了。这里插一句糙率的事。糙率曼宁系数虽然是在 HEC-RAS 里设置的但它的取值依据跟土地利用类型、河道形态直接相关。在做断面切割的时候可以顺手把河道的糙率分区一并考虑进去把主槽、滩地、居民区分别标记出来这样后面在 HEC-RAS 里赋值就有了依据。3. HEC-RAS 水动力模拟核心环节3.1 几何文件导入与模型检查几何文件导入 HEC-RAS 以后第一件事不是急着跑计算而是做模型检查。打开几何数据编辑器逐条断面检查以下几个方面断面桩号顺序是否正确从上游到下游递增左右岸标记是否与实际河道一致左岸指面向下游时的左手边断面地形是否合理有没有异常跳点比如突然冒出一个 500 米高程的孤立点那多半是 DEM 没填洼干净或者提取的时候取到了噪声点断面之间的中心线连接是否连续有没有交叉、遗漏。HEC-RAS 里有个很实用的功能叫“River Schematic”就是那张河网示意图。几何文件导入后河网示意图上会显示河道和水流方向如果箭头指向不对说明中心线方向有问题赶紧回到 ArcGIS 里改。还有一个我踩过坑的地方——断面范围。HEC-RAS 计算的时候断面两端分别是左岸和右岸的终点如果断面范围没有覆盖到可能淹没的范围水位一高水就会从断面外侧漏出去算出来的淹没范围偏小。这种问题在界面上还不容易发现要到“Plot Cross Section”里看水位线是否越过了断面端点。3.2 设计洪水过程与边界条件构建HEC-RAS 支持稳态Steady Flow和非稳态Unsteady Flow两种计算模式做洪水风险图两种都会用到。稳态计算相对简单适合推求某一频率洪水比如 20 年一遇、50 年一遇、100 年一遇的最高水位和最大淹没范围。只需要给定上游边界流量和下游边界水位模型就会沿河算出每个断面的水面线。非稳态计算要复杂一些需要给定流量过程线时间序列适合做洪水演进、淹没过程分析成果包括不同时刻的淹没范围、最大淹没水深、淹没历时。做风险图一般建议至少用非稳态算一遍因为很多地区的淹没风险不是由洪峰流量单独决定的洪水的涨落过程、下游顶托都对淹没范围有影响。上游流量怎么定这个要靠水文计算不是 GIS 或水力模型能解决的。常用的方法是查本地区暴雨洪水手册用推理公式或者瞬时单位线法推求设计洪峰流量和洪水过程线。如果项目区附近有水文站直接用实测洪水资料进行频率分析用 P-III 型频率曲线求各种重现期的设计洪水这是最可靠的办法。下游边界条件的选择也有讲究。感潮河段或者靠近大江大河干流的支流河口一般用给定水位或水位流量关系如果下游没有控制断面用“正常水深”加上一个能坡一般取河床坡降 0.001 到 0.005也能凑合。但能坡取多少一定要看实际的河道纵坡取大了水位偏低取小了水位偏高两个方向的偏差在评审会上都容易被质疑。3.3 桥梁涵洞在水动力模型里的处理河道上的桥梁、涵洞、堰坝对洪水水位的影响非常明显忽略它们算出来的水位可能低一到两米。我在项目里遇到过最夸张的情况一座老桥的桥孔过流能力不足壅水达到 1.5 米直接让上游一片原本不淹的农田变成了淹没区。HEC-RAS 里处理桥梁有两种常用方式第一种是直接用 Bridge 工具建模。需要输入桥面的顶部高程、桥孔宽度、桥墩尺寸和数量、桥面底部高程等信息。HEC-RAS 会根据水位高低自动判断是桥孔过流还是桥面溢流计算相对精细缺点是桥梁构造数据不好拿老桥根本没有设计图纸。第二种是用 Inline Structure线状建筑物模拟。适用于堰坝、水闸等规则建筑物设定堰顶高程、宽度、流量系数就行操作比 Bridge 简单很多。实际项目里我的建议是重点桥梁必须用 Bridge 工具建普通小桥如果资料不全可以简化成 Inline Structure甚至在有足够数据支撑的情况下忽略。但要注意凡是做了简化的桥梁都要在技术报告里写清楚评审专家问起来要能说出依据。涵洞的处理思路类似HEC-RAS 里有 Culvert 工具需要输入涵洞类型圆管涵、箱涵、进口形式、尺寸、进出口高程、糙率等参数。涵洞计算里最容易出错的是淹没出流和自由出流的状态判断软件会自动处理但你要确认输入的进出口高程和实际地形一致。3.4 模拟计算、结果判定与导出几何、边界、建筑物都设置好以后就可以开始计算了。稳态计算直接运行 Steady Flow Analysis一般几分钟就能算完。算完以后一定要看“Errors, Warnings, Notes”窗口很多新手忽略这个结果模型算了也出图了但水位是错的还不知道错在哪。非稳态计算要重点设置计算时间步长和输出间隔。时间步长一般取 1 分钟到 5 分钟地形变化剧烈的山区河道取小值平原河道取大值。输出间隔影响结果文件的大小做风险图需要最大淹没水深可以设成 1 小时输出一次计算完用“Max Depth”来出图。计算完以后第一个要检查的是水位合理性。把算出来的水面线和实际调查的洪水痕迹点对比如果出入大于 0.5 米就要回头检查糙率、断面、边界条件。第二个要检查的是纵剖面线是否光滑如果出现锯齿状跳动说明有断面局部壅水或者糙率突变需要调整。结果确认没问题以后就用 RAS Mapper 导出成果。RAS Mapper 里可以加载地形、显示水位、水深、流速分布操作很直观。导出的时候选择“Depth”栅格范围选整个模型范围分辨率尽量和原始 DEM 一致这样回到 ArcGIS 里叠加分析最省心。这一节最后提一个非常关键的细节导出栅格前要在 RAS Mapper 里把坐标系设置正确。HEC-RAS 6.x 默认以地形 DEM 的坐标系为基准但有时候你导入的地形没有正确带投影信息导出的栅格就是个无坐标系的裸数据回到 ArcGIS 里还得手动定义投影一不小心定义错了就全部偏移。4. ArcGIS 端洪水风险图制作与出图4.1 模拟结果栅格导入与坐标统一从 RAS Mapper 导出的水深栅格回到 ArcGIS 里第一步是检查坐标系。我一般会先在图层属性里看栅格的坐标系如果是 Unknown直接右键“定义投影”选成和 DEM 一致的坐标系。这里必须说清楚一个概念定义投影Define Projection和投影转换Project Raster是完全不同的操作——定义投影只是告诉软件数据是哪个坐标系不会改变坐标值投影转换才是真正的坐标变换。很多人搞混导致数据跟底图对不上。如果栅格坐标系正确但和底图对不上那就需要做投影转换了。例如你导出的栅格是 CGCS2000_3_Degree_GK_CM_114E而底图是 WGS_1984_UTM_Zone_49N就要用“投影栅格”工具转换过去。栅格导入以后第一眼看到的往往是黑乎乎一片因为水深栅格的范围里没有淹没的区域是 NoData显示为黑色。这个不是数据坏了只是显示设置问题。在符号系统里把 NoData 显示为白色或者透明就能看到清晰的淹没范围了。4.2 淹没范围矢量化与水深分级水深栅格看清楚以后做两件事。第一把淹没范围提取出来。用“条件函数”工具Con把大于 0 的像元赋值 1小于等于 0 的赋值 NoData这样就能得到一个二分栅格再用“栅格转面”工具转成矢量面。转出来的面可能有大量细碎小图斑用“消除”工具按面积阈值合并掉一般小于 1000 平方米的碎图斑直接消除地图整洁度会大幅提升。第二做水深分级。洪水风险图不是只看一个“淹了还是没淹”关键要看淹多深。水深分级国内各行业规范略有差异但大致可以参考水深区间风险等级适用说明0 0.5 米低风险一般影响步行和车辆通行0.5 1.5 米中风险影响一层建筑人员需转移1.5 3 米高风险低层建筑进水救援困难大于 3 米极高风险严重威胁生命财产安全操作上用“重分类”工具把水深栅格按上述区间分成 1、2、3、4 级然后“栅格转面”得到分级面。或者用“栅格计算器”不断 Con 嵌套也行但层级多了容易出错我推荐先重分类再转换逻辑更清晰。转成矢量面以后每个面要素携带一个水深等级属性用不同颜色渲染区分这就是洪水风险图的核心图层了。4.3 风险等级划分与空间叠加分析单纯的淹没水深还不等于风险真正完整的洪水风险图还要考虑承灾体的脆弱性。ArcGIS 里做这一步用的是空间叠加分析。你需要准备土地利用现状、人口分布、道路网络、居民点等数据然后做以下几种典型分析淹没范围内人口统计用“面转栅格”把人口分布数据转成栅格再用“分区统计”Zonal Statistics按水深等级分区统计人口数量淹没范围内土地利用损失把淹没范围面和土地利用面做“相交”Intersect按地类统计淹没面积转移道路可达性把淹没水深大于 0.5 米的区域叠加到路网上标注出不可通行路段为防汛转移预案提供依据。这些分析的产出比单纯一张淹没水深图要实用得多。我做过的几个项目最后评审专家问得最多的就是“你这个范围里大概影响多少人”“有没有妨碍转移的路段”如果没有提前做叠加分析现场根本答不上来。风险等级图可以在水深分级的基础上再叠加人口密度和土地利用的权重得到一个综合风险等级——高风险区、中风险区、低风险区。具体权重可以根据当地实际情况调整每个区域标准不完全一样做之前先和项目负责人确认清楚别自己想当然。4.4 专题图整饰与统计表输出风险图最终要变成一张规范的成果图ArcGIS 的制图功能这时候派上用场。进入布局视图设置纸张大小一般用 A3 或者 A0 出图插入以下要素指北针、比例尺、图例、标题栏、坐标网。字体建议用宋体或者黑体地图要素的图层顺序按“影像底图 → 面状风险图层 → 线状道路水系 → 点状居民点”排列。制图里面有几个容易忽略的细节图例里的水深分级色带要跟正文描述一致别换个配色就完全对不上了比例尺设置要跟当前地图范围的投影单位一致用“动态比例尺”模式缩放的时候比例尺自动更新标题栏里要写清楚洪水频率如“100 年一遇洪水风险图”和坐标系信息风险图有多个洪水频率的版本出图前一定要在文件名和标题里写清楚不然后期整理报告的时候分不清哪个是哪个。统计表输出可以用“属性表导出”功能把淹没面积、影响人口、影响道路长度等指标整理成 Excel 表格作为风险报告的附表。导出的时候字段编码要是中文乱码用“字段计算器”把字段名先改成英文导出完再映射成中文可以绕开乱码问题。做完这些一套标准的洪水风险图成果就算齐了风险等级图、水深分布图、淹没范围图、统计报表一个都不少。5. 常见问题与排查技巧实录5.1 河网提取断裂或流向混乱这是 ArcGIS 水文分析里最高频的问题之一症状是提取的河网断断续续支流接不上干流或者流向箭头乱指。排查思路按顺序来检查 DEM 是否填洼了。没填洼直接做流向是最常见的断裂原因检查流向栅格的赋值范围D8 流向值只有 1、2、4、8、16、32、64、128 这 8 个值如果出现 0 或者负数说明 DEM 有 NoData 区域先做“填洼”或“条件函数”补值检查河网提取阈值阈值太大河网会断在细支流上调小阈值试试检查 DEM 投影单位经纬度坐标做流量累积会出现大面积异常值必须改成投影坐标。5.2 断面切割线方向不统一断面线在 HEC-GeoRAS 里每个断面都有左右岸属性方向不一致意味着左右岸属性错乱导入 HEC-RAS 以后模型会报错。这个问题的根源几乎都在中心线。中心线折点方向不统一——上半段上游到下游下半段反向生成断面的时候方向就全乱了。解决方法是画中心线之前用“编辑工具”里的“反向”功能统一线的方向确认所有中心线都是从上游指到下游。还有一个检查技巧在 ArcGIS 里给断面线加上箭头符号运行“符号系统 → 高级 → 箭头”看所有箭头是否指向河流一侧。如果有几根箭头指向反了单独选中“反向”那根线就可以了。5.3 HEC-RAS 迭代不收敛怎么调非稳态计算最常见的问题是迭代不收敛报错信息五花八门但本质都是水位或流量在相邻时间步之间振荡。我的排查顺序是缩短计算步长从 5 分钟改成 1 分钟检查断面间距如果相邻断面距离超过 500 米在 HEC-RAS 中插入中间断面检查下游边界下游水位给得太低或者能坡取得太陡容易导致计算发散把计算选项里的θ权重从默认值调到 1全隐式稳定性显著提升如果还不行退回检查几何文件——很可能是断面范围不足水位顶到断面端点导致水量泄漏。5.4 栅格导入 ArcGIS 背景全黑或位置偏移水深栅格导入以后全黑先别慌。右键图层 → 符号系统 → 拉伸换一个色带NoData 设成透明绝大多数情况就好了。如果还是黑的用“属性”查看栅格最小值和最大值有时候栅格是浮点型范围显示为 0 到 2.147e9其实是某个异常的像元值在“符号系统”里把拉伸范围手动改小就行。位置偏移的问题90% 是坐标系不统一。检查栅格和底图的坐标系不一致就投影转换千万别手动平移硬凑。5.5 中文路径、坐标系混乱引发的麻烦事这是我最想强调的最后一件事。ArcGIS 和 HEC-RAS 的中文路径兼容性都一般尤其是 HEC-GeoRAS 导出几何文件的环节路径里一有中文就容易报“无法读取数据”或者生成空文件。我的习惯是给每个项目单独建一个纯英文路径的工程文件夹比如 D:\FloodRisk\Project_Zhanghe所有中间数据、模型文件都放里面从源头避免这种低级问题。坐标系混乱我前面也提过几次。定义投影和投影转换是两个工具很多人随手用“定义投影”强行把矢量数据定义到另一个坐标系结果位置不变但数值全错。记住一句话不知道数据原本用什么坐标系先打开属性看能查到最好查不到问提供数据的人千万别瞎猜。最后说几句体会做洪水风险图这几年我最深的体会是模型本身不难学难的是把数据管清楚、把每一步的原理想明白。曾经为了一个断面的高程对不上反复折腾了两天才发现是中心线折点方向导致切割线取错位置这类问题工具书里不会写全靠经验积累。如果你准备上手这类项目我的建议是不要一上来就追求非稳态、高精度先从一条小流域的稳态计算跑通全流程把 ArcGIS 和 HEC-RAS 的数据流打通再逐步加桥梁、加非稳态、加风险分析。工具永远在更新但“GIS 准备数据 — 水动力计算 — GIS 表达成果”这条主线不会变。最后分享一个小技巧每次做完整套流程以后把中间的所有关键参数DEM 来源、填洼阈值、河网阈值、糙率取值、边界条件整理成一张配置表存到项目文件夹里。这不是给甲方看的是给三个月后的自己看的——到那时候你绝对会感谢当初记了这张表。