一座 6.2 km² 的铜矿区，68 个钻孔、19 条剖面，模型交付那天在屏幕上转得很顺，评审说漂亮。第二年坑体往下采了 40 m，安全科问东帮最终边坡在雨季会不会失稳，工程部问涌水量会不会翻倍，打开同一套模型，没人答得上来。三维地质模型数值模拟要解决的就是这一步：把能看的模型，变成能算的模型。
一、模型的三种用法：记录、判断、推演
同一套块体有三种用法，要的数据不是一套。当记录看，只要形态对：剖面切得上、块体能闭合、图能出。当判断用，要的是估值：块段尺寸、品位插值的方式与分类边界，回答多少吨、什么级别。当推演跑，多出一个时间维度：应力、位移、水位随时间怎么变，边界条件怎么给，岩体用什么本构。前两层做不好，推演一定做不了；前两层做得好，推演仍然可能一步都跑不动，缺的是参数与边界。
多数项目停在第一层，不是技术不够，是需求没写清。判断线上的第一句话该是工程问题，不是软件名。写成“东帮最终边坡在雨季高水位下安全系数会不会低于 1.25”，网格尺寸、参数表、边界条件就都有了着落；写成“做个模拟看看”，做出来的多半是一段能转的动画。问题句里有会不会、多久、多少，才需要推演；只问长什么样、多少吨，静态模型足够。
需求方说不出一个能用数字验收的工程问题，先停下来，把那句话写出来再开工。这句话写不出来，后面每一步都会变成试错，网格划三遍也定不下来。
二、三类推演场景，各自要什么
三类场景对应三套方程，混着做是第一处浪费。采动应力与地表沉陷走连续介质力学，把岩体当等效连续体，给开挖步序，算应力重分布与沉陷盆，常用 FLAC3D 这类显式差分程序。边坡与排土场走极限平衡或强度折减，算安全系数与滑面位置，MIDAS GTS NX 这类有限元程序能直接做折减。含水层与矿井涌水走渗流，给渗透系数与补给边界，算降深与涌水量，MODFLOW 这类有限差分程序是主力。
表1 三类推演场景的输入、输出与常见误用
| 场景 | 求解类型 | 进场前必须齐全的输入 | 交付的输出 | 常见误用 |
|---|---|---|---|---|
| 采动应力与地表沉陷 | 连续介质力学（FLAC3D 类） | 岩体力学参数、开挖步序、原岩应力 | 应力场、位移场、沉陷等值线 | 把三年开采压成一步算 |
| 边坡与排土场稳定 | 极限平衡或强度折减（MIDAS GTS NX 类） | 结构面产状、地下水位、堆载 | 安全系数、潜在滑面 | 只交一个安全系数 |
| 含水层与矿井涌水 | 渗流（MODFLOW 类） | 渗透系数、给水度、补给与排泄边界 | 降深曲线、涌水量 | 四周边界设成不透水 |
表里最后一列是三处方向性错误。采动沉陷最容易漏步序，三年开采压成一步，位移一次到位，应力路径全丢，沉陷盆比实测深三成。边坡最容易只交一个数，评审要的是滑面位置与治理建议，一个安全系数撑不住答辩，评审潜规则那篇写过评审到底看什么。涌水最容易把边界设死，四周全不透水，水只能从坑里出，算出来偏小一半，实测一对比就露。
一个模型只回答一个问题。两个问题要两块网格，或者分步耦合；硬塞进一个模型，第一次改参数就得全部重划。场景没定就动手划网格，返工不在网格上，在参数表与边界条件上。
三、从面到体：三条接口路径
建模软件里做出来的是面与体，求解器要的是单元。中间这段路只有三条走法，选错的代价通常第三周才显出来。
直接导出原生网格，适合地层不超过五层、单元规整、只做一次的项目。控制权在导出按钮上：块段 45 m，导出的单元就是 45 m，坑体附近想加密得回建模软件改块体模型，改一次重导一次。中间格式 STEP 或 DXF 几何干净，多数前处理软件能读，属性会丢，层号、岩性、参数全变成无名曲面，得在求解器一侧重新赋；缝合容差要试，给小了缝补不上，给大了薄层被吃掉。脚本重建读库、按层号映射本构与参数、生成网格、写命令流，前期投入大，之后改模型只改参数表，适合每季度都要更新的项目。
三条路径都要交同一张卡片：模型原点、单位、正方向、坐标带号。这张卡片缺一项，接手的人要拿两天去猜。接口跑通的标准不是能导入，是导进去能算出一个已知答案：拿一块规则六面体试算，位移与解析解差超过 5%，回去查单位与原点，别在网格上找原因。
四、网格划分的四个坑
三维模型网格划分是模型与求解器之间唯一一次不可逆的转换，坑都在几何上，看不出来，算出来才知道。非流形边与自相交排第一：2.1 万个面片里 37 条边被三个面共用，体积计算出错，求解器直接报错。缝合容差先小后大，从 0.05 m 起步，超过 0.1 m 就会吃掉薄层；缝合后重算体积，与块体模型差超过 3%，说明还有洞没补上。
单位与坐标、地层归并这两处坑，动的都是参数。毫米与米混用，一份 FLAC3D 建模文件能从 180 m 高变成 180 000 m；原点离大地坐标 400 多万米，浮点精度先丢在后两位，位移能比监测值大一个量级，做法是局部坐标建模，平移量单独记在卡片上，回大地坐标再补。11 层并成 5 层看着清爽，1.2 m 的软弱夹层被并进 8 m 的厚层，加权平均后黏聚力偏高近三成，安全系数跟着偏高；单层厚度不足网格尺寸三分之一的层，要么保留并局部加密，要么把归并规则写进报告，参数按厚度加权。
网格尺寸与算量失衡排第四。全局 15 m 是 214 万单元，单轮 9 小时；全局 45 m 是 76 万单元，单轮 3.2 小时。坑体与边坡附近 15 m、远场 45 m，才算得动；长宽比超过 8 的单元先处理，否则收敛慢，结果还在这些单元上失真。单元数先定上限，按机器 8 小时能跑完一轮为准，超了就在远场放粗；网格划不动多半是数据没过关，建模前必做的 7 类数据体检里列的孔深闭合与剖面偏差，最后都会变成几何缺口。
五、参数反演：用监测数据校准
手册给的力学参数区间很宽，黏聚力 20–40 kPa 都算合理，取中值算出来的位移可能差一个量级。反演做的是把区间收窄：拿现场监测到的位移、水位、应力，反过来找那组能重现实测的参数。
先挑敏感参数，一般三到五个。逐个做敏感性试算，参数上下动 30%，看输出动多少；动不动的参数不参与反演，动得多的才值得调。挑完再谈算法，试算与最小二乘都能用，关键是留验证点。
位移反演黏聚力与内摩擦角。边坡布 12 个监测点，9 个参与反演，3 个留作验证；月位移 3–8 mm 的区段，反演后黏聚力从 26 kPa 落到 21 kPa，内摩擦角从 34° 到 31°，验证点残差 2 mm。水位反演渗透系数，3 个观测孔埋深 8–14 m，渗透系数从 0.9 m/d 调到 1.5 m/d，涌水量预测从 1 800 m³/d 修到 2 300 m³/d，与实测对上。
验证点残差超过监测精度两倍，别再调参数，回去查本构与网格。把参数调到能对上全部监测点，通常是把模型调成了曲线拟合，换一个工况就崩。反演结果写成区间，不写单值，位移给范围，安全系数给最不利组合。
六、动态更新：模型不是一次性交付物
一次性交付的模型，第二年就只剩展示价值。坑体推进、补勘加孔与水位变化，任何一项都能让上一版的边界条件失效。可行的做法是把更新做成流程，而不是出事再返工。
三个触发条件写进合同：新增钻孔超过原有孔数 10%，坑体推进超过一个台阶 12 m，观测孔水位变幅超过 2 m。任一项触发就排一次重构，不等年度评审。
流程五步：库更新、块体重算，网格重划之后复用参数，再与上一版出差异图。参数复用是省时间的关键，反演过的参数存成参数表，重构时只改几何不改参数，除非监测显示参数也变了。版本要留基线，旧版不覆盖，位移累积量不重置，两版之间才能算出增量。
一次季度重构 3–5 天，含出图；从零重建要两周以上。连续两个季度不更新，模型在评审里只能当示意图用，别指望它支撑治理方案。
七、精度与算量的取舍
三维数值模拟的算量由三个圈决定：影响圈按坑深的 3 倍取，边界放到坑深的 5 倍，网格尺寸取台阶高的三分之一。边界离坑底不足两倍坑深，结果会被边界效应污染，位移偏大还是偏小，取决于边界约束怎么给。
精度不是越细越好。网格加密一倍，单元数涨四倍，误差可能只降几个百分点；把力气花在参数与边界上，回报比加密网格高。算量与费用的关系，三维地质建模到底多少钱那篇拆过四因子，网格规模是其中之一。
一个可用的判据：同一工况用粗细差一倍的两套网格各算一次，关键量差超过 10%，说明还没收敛，继续加密或改局部尺寸；差在 10% 以内，这套网格可以定稿。
八、边界与风险
连续介质假设有适用边界。节理密集、结构面控制的岩体，等效参数能把平均行为算对，把破坏路径算错；这类边坡要做离散元或块体模型，不能拿连续体结果直接定治理方案。
参数不确定性消不掉。反演能把区间收窄，收不没；同一组监测数据换一种本构，参数能差两成。推演不能替代勘察，不能替代评审，也不承诺位移绝对值，它能做的是把方案之间的差别算出来，让治理顺序有依据。
多场耦合只在有明确水力梯度时值得做。渗流与应力耦合的计算量是单向的两三倍，没有水位实测支撑，耦合出来的孔压场只是假设。结论用于设计之前，至少两个独立监测源互相印证；只有一个源，把结论写成趋势判断，不写绝对值。
九、结语
模型的精度不由网格决定，由要回答的问题决定。先写出那句能用数字验收的工程问题，再决定模型做到什么份上；顺序反过来，做出来的模型一定很漂亮，也一定没人敢用。三维地质模型数值模拟的门槛不在求解器，在那句问题、那张参数表、和每年肯不肯重算一次。
