ARTICLE · 1245038
把标准动作装进插件——《地质灾害风险评估方法探索》·基础篇·第三部分·第2篇:实操① SEEP/W 渗流建模
上一篇我们拿 GeoStudio 给插件的边坡稳定性引擎做了一次全面对照——五处对账、五处修复，最终两边在同一种破坏模式上会师。尺子校准了，这一篇开始讲怎么用。先讲第一段：SEEP/W 渗流建模——把"降雨怎么进坡体"这件事，在插件界面上一次配齐。GeoStudio 里要十个步骤点出来的东西，插件收拢成一块参数面板、一个按钮。
一、一个页面，两段工作流
打开插件的"边坡稳定性"页，你会看到整块界面就是第 1 篇那张流程图的落地版，从上到下三大块：
GeoStudio 十步操作 ↔ 插件参数面板：左边一列是 GeoStudio 的动作，右边是插件里对应的输入位。
- ① SEEP/W 渗流模型建立
：剖面几何、岩土材料、降雨工况、初始条件——GeoStudio 的"新建工程→绘制几何→定义材料→建函数→赋值→设边界→建降雨函数→划分网格"八步，收拢成十几个输入框； - ② SLOPE/W 边坡稳定性分析
：计算方法下拉（七个条分法）、条分数、关联渗流结果——GeoStudio 的"关联→选方法→搜索→求解"四步； - 操作行
：三个按钮对应三种跑法——"① 渗流求解"、"② 稳定性求解"（须先完成①）、"一键全流程"。
两段式的好处和 GeoStudio 一样：渗流可以单独跑、单独看孔压场，满意了再做稳定性；换一种条分法重算②时，不用重跑几十分钟的渗流。
二、剖面几何：六个数字，一个都不能想当然
几何面板六个参数：坡面高差、坡顶宽、坡脚高、坡角、覆盖层厚、网格尺寸。默认值就是点位1：30 / 25 / 10 / 31° / 3 / 4 m。
第 1 篇对账时摔过的坑，这里展开讲清楚。头一个是"坡高"到底指什么：复现手册写"坡高 40 m"，我们第一版把它当成了坡面高差——结果界面输入框里填 40，生成的坡顶高程变成 50 m，整个坡高了三分之一。原工程文件里 40 其实是坡顶高程（坡脚平台 10 m，坡面高差只有 30 m）。所以插件把这一栏直接命名为"坡面高差"，不给歧义留门。
第二个是界面形态。下拉框里有两档：
- 等厚
：覆盖层界面 = 地表竖向下移"覆盖层厚"——一块厚度均匀的板； - 楔状（点位1原工程）
：界面倾角缓于坡面，覆盖层在坡肩处最厚（13 m）、往坡脚收窄到 3 m——原工程真实的几何。
别小看这个下拉框。第 1 篇说过：等厚 3 m 的薄板里滑面稍一切深就进了基岩，怎么搜 FoS 都是 2.3 往上；换成楔状，粘土层才有厚度容得下那条沿界面走的临界滑弧——同一个坡，两种几何，两个世界。
参数校验也埋在这一层：坡角限制在 5°~80°，覆盖层厚必须小于坡高+坡脚高，楔状模式要求坡脚高大于覆盖层厚——填错会直接拒绝并说明原因，而不是算出一个看似合法的错数。
三、材料面板：两行八字段
材料区两行——覆盖层·粘壤土、基岩·页岩，每行六个字段：
前四个是第 1 篇讲过的强度与渗透主角；后两个是非饱和渗流的"性格参数"。把插件的默认预设画成曲线：
插件默认预设下两种土的有效饱和度 Se 与相对渗透系数 kr 随吸力的变化。页岩的进气值（≈1/α）比粘壤土小一个量级，失水更快、更"娇气"。
读法三句话：α 是进气值的倒数——ψ 越过 1/α，土开始明显排水；n 控制曲线陡峭度——n 越大，土在进气附近"放手"放得越急；kr 决定非饱和区还能过多少水——粘壤土在吸力 1 m 时 kr 还有 0.2 量级，而页岩已经掉到 10⁻³ 以下，这正是"雨水在粘土层里走、页岩近似不透水"的曲线语言。插件默认取内置质地表的粘壤土/页岩预设——与 GeoStudio 原工程用数据点曲线拟合的方案是同一物理量的两种表达（第 1 篇对账时已确认这属于"可解释残差"）。
粘壤土与页岩的 k 差了 300 多倍，这个比数值得记住：它意味着界面天然是一道"毛细屏障"（第二部分讲过的老朋友）——雨水下渗到界面就走不动了，水在粘土层内积攒，孔压沿界面抬升。降雨型滑坡的很多戏，都发生在这一层薄界面上。
四、降雨工况、初始条件：三个框，一个单位陷阱
雨强列表是一个逗号分隔的文本框，默认 2, 20, 100（mm/d）——每个雨强一个独立工况，最多 5 个。时长默认 80 h（docx 口径）。初始水位默认 10 m——恰好与坡脚平台齐平，这意味着坡脚一带的覆盖层天然饱和。
单位陷阱值得单独一节。GeoStudio 的工程文件内部以英制存储：我们解包点位1的 .gsz 时，降雨边界写的是 7.59e-08——这是 ft/s，换算回来正是 2 mm/d；几何坐标 131.23 其实是 40 m。复现手册也专门警告过："原材料后部个别分析文字出现 mm/h，但工况表记录为 mm/d，必须结合工程单位核对通量换算。"插件内部统一 SI（m/h），界面上只认 mm/d——你在 GeoStudio 里抄参数时，抄数值之前先抄单位。
五、网格尺寸：精度与时间的天平
网格尺寸默认 4 m，范围 0.5~20 m。经验标尺（以点位1 80 h 单工况为准）：
- 4 m
：约 20~30 分钟，孔压场形态正确； - 2~3 m
：约 1 小时起； - 1 m
（GeoStudio 原工程口径）：小时级起步，低雨强可能数小时。
为什么低雨强反而慢？雨强越小，每一步水头变化越小，非线性迭代的步长压得越碎——2 mm/d 工况的耗时可以超过 100 mm/d。这是当前内置开源引擎与 GeoStudio 商用解算器差距最明显的地方（补充篇会放进"已知边界"清单）。
网格生成后有一道面积守恒校验：每种材料的三角网总面积必须对得上几何多边形，缺一块就报警拒算——这道闸门是对账时那个"覆盖层整除消失"bug 换来的，保证条分和渗流永远不会在"缺了一块的网格"上静默出错。
六、点下"① 渗流求解"之后
按钮按下，引擎在后台做五件事：生成结构化三角网（界面/地表逐列对齐，材料分区 1=覆盖层、2=基岩）→ 填充可复现的输入工作簿 model.xlsx → 组装 Richards 方程有限元 → 沿时间轴瞬态推进（每 1/10 时长存一帧孔压场）→ 导出末帧孔压与逐帧存档。
完成后看三样东西：
- 横幅
： 工况 3/3 收敛——不收敛的工况会在②里被跳过并注明； - 产物表
：每个工况目录下有四件套—— model.xlsx（完整输入工程）、model_mesh.json（网格）、model_seep.csv（末帧孔压）、seep_frames.npz（逐帧孔压场，②画 FoS-时间曲线要用）； - 输出目录
：插件数据目录下的 outputs/slope_stability/<运行号>/，一个运行一个文件夹，report.json 记录本次全部参数与结果。
到这里，①的活干完了：每个工况都攥着一沓"坡体内部孔隙水压力随时间演化"的快照。下一站，把这沓快照交给 SLOPE/W——一个一个安全系数是怎么被搜出来的。
（本篇是"基础篇·第三部分"第 2 篇；下一篇讲②稳定性求解与成果解读。）
