观潮 · 模式预报与气象 AI 观察
门头沟23·7泥石流：出沟口5.8米每秒，进村只剩1到3米，最大堆积4.06米
> 2023 年 7 月 29 日至 8 月 2 日北京遭遇历史罕见特大暴雨，门头沟区斋堂镇白虎头村发生滑坡型泥石流。作者用无人机正射影像建出 0.266 米格网的 DEM、室内土工实验定出容重与强度参数，再用 DAN3D 的 Voellmy 模型反演整场运动，把反演出的流速与流深直接当作荷载送进结构计算，逐构件核算距沟口约 60 米的一栋砖混—框架混合结构民居。结果给出了一条完整的链条：泥石流出沟口时流速 5.8 米每秒，冲进居民点时降到 1—3 米每秒，最大堆积厚度 4.06 米；砖砌墙体因抗剪承载力不足而倒塌，钢筋混凝土梁柱因抗弯抗剪储备足够而保留。作者同时说明了三处没做的事：梯田拦挡造成的 2 号沟堆积只作现象解释、未做定量推断；因缺乏高精度灾前影像，侵蚀率只能参照前人方法反算；二层泥石流等效荷载按一层的 30% 简化取值。
缺的是哪一块
按动力学条件，泥石流分土力类和水力类两种。土力类以高位斜坡土石体失稳后的滑动、错落、崩塌、坠落为动力，运动中因震动和液化导致土石体力学状态突变、结构丧失、流态化加剧；滑坡型泥石流是其中典型，隐蔽性强、预测难度大、运动速度快，是我国半干旱半湿润山区防灾的难点。
滑坡型泥石流的破坏对象是山区的房子，但「成灾动力学」和「建筑物结构损毁」这两段长期是分开做的：基于光滑粒子流体动力学（SPH）的 DAN3D 之类工具在复杂地形、大变形和侵蚀挟带上有优势，已被广泛用于泥石流的动力过程反演；而北方土石山区已有的风险研究大多只做单一泥石流动力过程，或者只做风险区划。动力模拟给出的流速、流深，和结构工程师习惯用的荷载与构件承载力之间，缺一座桥。
白虎头村这栋房子的灾后形态本身就是线索：墙体倒了，梁柱没倒。这种破坏模式恰好可以用来验证「把模拟流速流深当荷载输入」这条路走不走得通。
方法怎么搭的
整条链路分四段，每段的输入都是上一段的输出。
环节 · 做法
地形 · 野外调查加无人机航摄，用 Agisoft PhotoScan 合成像元 0.266 米 × 0.266 米的 DEM
物源参数 · 现场采样做室内土工实验：粒径小于 0.005 毫米的占 3.0%，最大内摩擦角 24.8°，黏聚力 18.2 千帕
动力反演 · DAN3D，选用 Voellmy 模型，DEM 划分成 1 米 × 1 米的均匀网格，设定两条沟域同步启动
结构计算 · 从模拟结果里取建筑所在位置的流速、流深，按 GB 50003—2011 和 GB/T 50010—2010 逐构件验算
参数取值上作者做了几处保守处理。最大侵蚀深度按 5 米输入（实测沟道最深多在 3—3.5 米），理由是「更好地覆盖不利地形」；室内试验给的容重 18.03 千牛每立方米，摩擦系数取 μ = 0.10、湍流系数 ξ = 500 米每二次方秒，取值参照实测、前人研究和敏感性分析。
物源量用沟道断面反算：先按断面间距估算体积，再反推平均侵蚀率。初步计算的物源体积约 920 立方米，最终体积约 19171 立方米——也就是说运动中挟带进来的物源比滑源本身大一个量级，这决定了后面为什么必须逐段看侵蚀深度。
结果数字
雨量与地形条件
项 · 数值
斋堂站（距研究区 3 公里）激发泥石流总雨量 · 460.4 毫米
最大 1 小时降雨量 · 45.1 毫米
最大 6 小时降雨量 · 100.4 毫米
最大 24 小时降雨量 · 327.4 毫米
灾害发生前平均雨量 / 发生时刻雨量 · 8.9 毫米每小时 / 45.1 毫米每小时
流动距离 / 高差 · 约 1000 米 / 约 383 米
沟道断面面积 · 12—40 平方米
泥石流容重 · 1.84 吨每立方米（黏性泥石流）
附近三村发生时刻 · 2023 年 8 月 1 日 1 时
运动过程分五个阶段（模拟给出的时间分段）
时段 · 过程
0—30 秒 · 两条沟域同时启动，斜坡中上部浅表层物源活化并沿主沟道运移
30—70 秒 · 两股泥石流分别沿各自支沟运移
70—90 秒 · 双沟汇流，形成复合型泥石流
90—120 秒 · 突破沟口进入堆积区
120—300 秒 · 完成扩散堆积，形成最大厚度 4.06 米、展宽约 130 米的堆积扇
速度
位置 · 数值
1 号沟 50 秒后稳定速率 · 4—5 米每秒
2 号沟峰值（10—30 秒，受地形落差驱动） · 14.6 米每秒
2 号沟 30—70 秒衰减后 · 6—7 米每秒
汇流后（70—80 秒） · 4—5 米每秒
距沟口 100 米束窄段（80—90 秒，束窄压缩效应） · 反弹至 5—10 米每秒
冲出沟口时残余速度 · 5.8 米每秒
沟口外居民点一带 · 1—3 米每秒
沟口外流速虽然只有 1—3 米每秒，但高密度流体携带的巨量物源仍对居民区建筑造成显著冲击，作者把这种「低速高冲量」特征单独点了出来。
堆积厚度（7 个监测点，7 号位于受损建筑附近）
监测点 · 位置 · 最大堆积厚度
2 · 2 号沟落差显著段 · 1.25 米
4 · 2 号沟前缘、梯田一带 · 2.5 米
3 · 1 号沟汇流前主沟道 · 0.5 米
5 · 汇合后距沟口 100 米束窄段 · 1.25 米
扇核心区 · 沟口外 · 4.06 米
侵蚀深度
位置 · 侵蚀深度
1 号沟上游 · 0.1—0.5 米
1 号沟汇聚点前 300 米内 · 0.5—1 米
2 号沟上游 · 0.5—1 米
2 号沟前缘 300 米内 · 0.5—2 米
两股汇流后主沟道 · 1—2 米
流经狭窄沟段 · 2—3 米
冲出沟口后 · 锐减至 0.1—1 米，转为堆积主导
验证：现场最大堆积厚度约 4 米，模拟 4.06 米，误差小于 2.0%；模拟最大侵蚀深度 3 米。
结构验算：建筑位于距沟口 60 米处，承重体系是 MU10 烧结砖砌体墙配 C30 钢筋混凝土构造柱与圈梁。荷载输入取建筑处流速上限 3 米每秒、流深 1.5 米、入射角 45°、爬高约 8 米（房屋总高约 8.4 米）。来流弗劳德数算出来是 0.782，小于 3，所以泥石流以反射波形式冲击结构。
项 · 数值
3 米跨墙体抗剪承载力 · 39.17 千牛
9 米跨墙体抗剪承载力 · 155.93 千牛
6 米跨墙体抗剪承载力 · 100.43 千牛
泥石流对垂直墙面的单位宽度冲击力 · 92.625 千牛每米
六面墙的抗剪承载力都小于所受泥石流冲击力，因此倒塌；梁柱因为有纵向钢筋与箍筋，抗弯、抗剪储备足够而保留——这与灾后现场「墙倒梁柱完整」的形态一致。
原文图
图1 研究区位置：紫线为白虎头村边界，红星是泥石流研究点，底图为高程着色。村落沿大独山西麓分布，背靠西山隆起带，周边发育密集沟谷，形成 25°—45° 陡坡与阶地交错的地貌，山地占比 98.5%，最高海拔 906.6 米。
图2 泥石流运动速率的时空演化，六个时刻依次为 t = 10、50、90、120、150、300 秒，色标 0—15 米每秒，黄色数字 1—7 为监测点。2 号沟在 10—50 秒率先出现高值区，90 秒后两沟汇流形成连续高速带，150 秒后堆积扇上速度整体回落。
图3 堆积厚度时间曲线，左（a）为汇聚前各沟监测点，右（b）为汇聚后监测点。2 号沟监测点 2 在约 50 秒冲到 3.9 米后快速回落，说明该处在汇流前就形成了临时堆积；汇聚后的监测点 6 在 200 秒后稳定在 3.5 米左右，对应扇面核心区持续堆积。
图4 结构计算的对象与荷载：（a）房屋示意，三跨分别为 3、9、6 米，总高 8.4 米；（b）泥石流冲击示意，冲击前流深 h₀ = 1.5 米、来流速度 U₀ = 3 米每秒、爬高 h₁ = 8 米。这两个数字就是后面墙体抗剪验算的输入。
作者承认的限制
限制 · 原文口径
梯田拦挡只作现象解释 · 2 号泥石流在梯田一带形成堆积「可能由多因素共同造成」，包括体积相对较小、局地流速变化导致能量快速衰减、梯田边界的拦挡与粗糙效应，作者明确写「本文在此仅作现象解释，未作定量推断」
侵蚀率靠反算 · 因缺乏高精度灾前影像，侵蚀区面积参照 Ellett 等、Jaboyedoff 等的计算方法，结合 0.266 米分辨率 DEM 的断面估算；侵蚀深度按 5 米保守输入，而实测沟道最深多在 3—3.5 米
二层荷载简化 · 爬高约 8 米、房屋总高约 8.4 米，二层梁可能承受局部横向泥石流作用；为与一层口径一致，二层等效荷载按一层的 30% 取值做简化复核
物源不全部施加于单体建筑 · 建筑处的流速 3 米每秒、流深 1.5 米是局地模拟结果，两条沟汇合后的物源仅用于过程重现，并非全部施加于该单体建筑
沟道展宽未出现的归因 · 2 号沟保持高坡度、势能较高，致使「未能出现沟道展宽现象」——这是作者的机理解释，未给出定量判别
单沟试算不足以重现 · 单沟情形的试算难以再现束窄段堆积与扇端轮廓，因此设定两条沟域同步启动
一句话判断
这篇文章的价值在于把「模拟出的流速流深」真的变成了结构验算的输入，并用「墙倒梁柱不倒」这个灾后事实校验了整条链路；但它同时也说明，梯田这类小型拦挡工程在极端降雨面前防护能力有限——真正管用的还是选址避让，或者在迎流面加一道钢筋混凝土挡墙。