观潮 · 模式预报与气象 AI 观察
河内内涝零资料建模：32.8 平方公里的 3Di 模型，8 个众包水深点里 6 个对到 2.2 厘米
> 城市水动力内涝模型通常要三样东西：实测地形、显式管网、雨量计。2025 年 10 月 7 日台风 Matmo 造成河内 Mỹ Đình 与 Tây Mỗ 内涝时，这三样一样都没有。作者的应对是：把 1 米 DEM 按河道湖泊地形、未测水体、下穿隧道和一处基坑系统修补，用按地表覆盖分配的入渗率顶替没有记录的地下管网，用 IMERG Late 半小时卫星降水做闭域二维强迫，再拿 8 个从新闻和社交媒体里挖出来的众包水深点做率定。结果是 6 个点对上了、误差在 2.2 厘米以内，另外 2 个点全程是干的。作者反复强调：8 个点全部用于调参、没有留出独立验证，所以这只是一个率定拟合，不是无约束的预测检验。
缺的是哪一块
河内位于红河三角洲，城区低地海拔只有 5 到 20 米，洪水是这个地理环境的固有特征。近二十年强降水事件频率显著上升，与东南亚城市在低地上的密集定居叠加。2024 年 9 月 7 日超强台风 Yagi 袭击越南北部，造成 21 人死亡，河内有 32 处水深超过 0.5 米，Chương Mỹ 县积水达 1.5 米并持续数日；2008 年 10 月 30 日至 11 月 1 日的洪水让主要街道积水达 1 米，交通瘫痪，越南北部与中国南方共 92 人死亡，经济损失超过 1.77 亿美元。Luo 等人用九个照片判读站点率定淹没模型，R² 只到 0.6——这说明用媒体影像做定性水深估计来约束河内洪水模拟有其先例，但也说明众包视觉证据能达到的拟合优度本来就有限。
2025 年 10 月这场是同类情形。台风 Matmo（越南编号 11 号）从中国广西内陆移入，10 月 7 日下午河内建设厅记录了 100 多个积水点，其中 23 个路口完全无法通行，Mỹ Đình 和 Tây Mỗ 是西部城区里受影响最重的两个。
水动力模型只能沿着它底层高程面里明确包含的拓扑路径走水，所以高程模型的选择和预处理给后续参数率定能恢复多少设定了上界。已有文献的做法基本都假定地表障碍是凸出地面、需要被剥离的东西（建筑、高架桥面、植被），却很少处理反向问题：未经处理的高程面把真实中过水的地下或下凹通道给堵死了。作者说，同时处理被遮挡的下穿隧道、人工路基、基坑和未测湖盆这四件事，此前没有研究在一个城市模拟域里一起做过。
三段数据缺口分别怎么补
缺口 · 替代做法
实测地形与水下地形 · 1 m 裸地 DEM 上修补河道与湖泊地形、未测水体、道路下穿、一处基坑
地下排水管网 · 按地表覆盖类型分配入渗率，作为未记录管网的集总替代
雨量计 · IMERG Late 半小时卫星降水，作为闭域二维强迫
模型用 3Di 亚网格求解器，在灵活四叉树网格上解二维浅水方程，并整合 1 米裸地 DEM 的亚网格地形。计算网格在 32.8 平方公里上共 20,838 个有效单元，单元边长 4 到 64 米，最细的 4 米和 8 米单元布置在下穿隧道、Nhuệ 河河道和地形陡变区。
计算域是水文闭合的，没有任何边界条件：既没有地下排水管渠、泵站，也没有外部边界水位和侧向入流。降水是唯一入口，地表入渗是唯一出口。这样能把严格落在流域边界内的降雨产流隔离出来，同时也意味着来自 Nhuệ 河上游的河道来水和市政管网的输移都无法复现。降水按二维面雨均匀分配到所有单元，床面摩擦用 Manning 公式，入渗用 3Di 简单入渗模块，两者都以空间栅格给出，模拟连续跑 65 小时。
水量平衡说明闭域把绝大部分降水留在了地表：总降水 7,989,142 立方米，参数化入渗带走 816,312 立方米（占降水的 10.2%），模拟结束时地表储存 7,935,230 立方米，恰好等于降水减入渗（7,172,830 立方米）加上前期储存 762,291 立方米，最大闭合误差只有 109 立方米。
结果数字
8 个率定位置里，6 个有水面，模拟水深与报告的中心值在每个湿点上都在 2.2 厘米以内，平均偏差 0.57 厘米，平均绝对误差 1.34 厘米，均方根误差 1.49 厘米。但 P01 和 P06 在全部 15 个输出时刻都是零水深，按 0 厘米计。
聚合误差指标主要由这两个干点决定，而不是那 6 个率定成功的点：
指标 · 8 个点 · 去掉 P01、P06
平均绝对误差 · 6.63 cm · 1.34 cm
均方根误差 · 11.39 cm · 1.49 cm
平均偏差 · 6.05 cm · 0.57 cm
均方根误差远大于平均绝对误差（11.39 对 6.63 厘米），说明残差集中在个别严重误差上，而不是均匀分布在模拟域内。两种口径下偏差都为负，即模型相对目击者的描述系统性偏干。时间敏感性也很现实：08:00 到 14:00 之间，模拟水深在 P03 从 25.5 退到 11.2 厘米、P04 从 40.2 退到 30.2 厘米、P05 从 39.9 退到 29.2 厘米，因此报告时刻差 30 分钟，模拟水深就差 2 到 4 厘米；而 P02、P07、P08 在整个 7 小时里波动不超过 4 厘米。
原文图
图：河内西部的研究区（红线）叠在卫星影像与坊（ward）界之上，范围约 32.8 平方公里，覆盖 Mỹ Đình、Tây Mỗ、Phương Canh、Nhân Chính 等城区，右侧标注比例尺。Nhuệ 河沿研究区东南边界穿过，Thang Long 大道的下穿隧道在极端暴雨时是主要的地形汇水点。
图：8 个众包水深观测点在研究区内的位置，观测日期为 2025 年 10 月 7 日。可见 P01、P06（Tây Mỗ Road 与 Road 70）位于研究区外围，其余各点分布在 Mỹ Đình 一带；两个位于高速公路走廊 Km4+944 与 Km8+200 的点意味着下穿隧道的水力是否算对，对局部率定和整体汇流都有影响。
图：(a) 各点观测水深范围与中心值（蓝）对照模拟值（橙），P01 与 P06 标为 dry——报告有 20 与 25 厘米积水、模拟全程为零；(b) 模拟对观测的散点，虚线为 1:1 线，P08 接近 1 米的深水点与 P07 都落在对角线附近，P03 则明显偏浅。
两个干点为什么比六个拟合点更有信息量
作者认为，P01（Tây Mỗ Road，20 厘米）和 P06（Road 70，25 厘米）是模型结构性缺陷唯一显形的地方，并列出三种都能解释这一失败、但在受控检验下会给出不同结果的物理假说：
- 如果 2017 年高程面里的农用池塘后来被回填、在模型里充当了人为蓄水点，那么换成更新的当代高程面，应该能在 P06 生成合理水深，而六个湿点不受影响。
- 如果 IMERG Late 的雨型系统性偏低（比 Mễ Trì 气象站低 47%），把降水整体放大会让两个干点变湿，但同时会抬高六个已率定点的水深，除非同时重新调整入渗能力。
- 如果存在未绘制的雨水管渠、涵洞或局地地表径流沿这两条道路走廊集中汇流，那么集总降水加均匀入渗根本分辨不了局地积水，调雨量和改地形都消除不了这个误差。
第 2 条和第 3 条之间还有一层等当量（equifinality）问题：建成区上率定出来的入渗率是 1 毫米/小时，远低于同类城市模拟文献里的 6 到 40 毫米/小时有效排水率。若卫星降水确实偏低 47%，这个极小的排水损失率正好可以作为一个等当量补偿——偏低的降水配极小的排水损失，与更高的雨量配真实的管网输移，会给出完全相同的地表积水量。六个湿点上的拟合与这种数学补偿相容，但也同样与「降水准确、真实基础蓄滞」相容，现有监测数据分不开这两种解释。
作者承认的限制
- 8 个点全部用于调参，没有留出独立验证点。 原文在摘要和结论里两次强调：这只是一个率定拟合，不能证明同一套参数化在其他未监测的城市集水区或下一场风暴中仍保持预测精度。
- 闭域边界排除了两类真实来源。 上游 Nhuệ 河的河道来水和市政管网的输移都无法复现，受这两个因素影响的积水深度模拟不出来。
- 率定后的损失率物理解释必须极其谨慎。 约 90% 的降水留在地表、入渗只占 10.2%，意味着模拟街道积水深度主要依赖总降水量和赋予不透水表面的最小损失能力，而降水与入渗之间的等当量交易使这些参数失去了唯一解释。
- 三项能分开假说的诊断实验都没做成。 用偏差订正后的雷达降水驱动、用未修补的原始地形做对比、以及对已回填农用池塘的当代地面测量，作者明说以现有野外观测都无法实现。
- 众包观测本身的不确定性不小。 原文在不同位置给出的报告不确定性区间分别为 10 到 25 厘米与 5 到 10 厘米，1.34 厘米的平均绝对误差落在观测噪声之内，无法与噪声统计区分。
- 相邻排水盆地是物理相连的。 运河、共用泵站和回水水位把河内各排水盆地连在一起，逐个盆地单独建模、逐时间步交换边界流量是作者提出的后续方向。
一句话判断
在资料稀缺的城市流域，用公开数据加众包水深点可以把二维水动力模型率定到厘米级，但那 2.2 厘米只能证明流程走得通；真正该看的是两个全程为干的点，它们把「卫星降水偏低、管网缺失、地形过期」三种可能一起摆在了桌面上。