8.1 矿山软件体系
现代矿山规划依赖地质建模、资源估算、矿山设计、生产调度与岩土分析等集成软件。本节梳理软件格局、数据流架构及各专业域工作流。
软件对比矩阵
| 软件 | 领域 | 优势 | 适用规模 |
|---|---|---|---|
| Leapfrog Geo | 地质建模 | 隐式 RBF 曲面,模型更新快 | 中大型 |
| Surpac | 地质 + 采矿 | 块体模型与矿山设计;国内应用广 | 中大型 |
| Datamine Studio | 地统计 | 变异函数、Kriging、资源报告 | 大型 |
| Whittle | 露天优化 | 嵌套矿坑壳;战略生产规划 | 可研–运营 |
| Vulcan | 地质 + 采矿 | 露天井下一体化工作流 | 大型 |
| Deswik | 调度 | 甘特图调度、设备仿真 | 生产计划 |
| Slide / RS2 | 岩土 | 边坡 LEM(Slide)、2D FEM(RS2) | 设计阶段 |
| FLAC3D | 岩土 | 3D 显式 FEM;大变形、地压 | 科研与设计 |
数据流架构
flowchart LR
drill[钻孔数据库] --> geo[地质模型]
geo --> block[块体模型]
block --> whittle[Whittle优化]
block --> fms[FMS调度]
fms --> erp[ERP/MES]
geo --> grade[品位控制]
grade --> fms
8.1.1 地质建模软件
8.1.1 地质建模软件
学习目标
- 按矿床复杂度与团队技能选择 Leapfrog、Surpac、Datamine
- 执行标准流程:钻孔导入 → 域建模 → 变异函数 → 块体模型 → 资源报告
- 比较隐式(RBF)与显式(剖面线框)建模
- 导出兼容 Whittle、Vulcan、FMS 品位控制的块体模型
核心概念
| 术语 | 定义 |
|---|---|
| 钻孔数据库 | 孔口、测斜、化验、岩性、密度表 |
| 域(domain) | 品位分布独立的地质单元(矿体、上盘等) |
| 隐式建模 | RBF 插值由数据点生成曲面 |
| 显式建模 | 剖面数字化线框三角化实体 |
| 变异函数 | 块金值、基台值、变程等空间连续性模型 |
| 块体模型 | 三维网格块体估算品位、密度、岩性 |
| Kriging | 地统计估算,遵循变异函数 |
| 资源分类 | JORC/NI 43-101 测定/控制/推断 |
参数与指标
| 参数 | 典型值 | 单位 | 说明 |
|---|---|---|---|
| 块体尺寸(露天) | 5×5×5 至 15×15×15 | m | 匹配 SMU |
| 块体尺寸(井下) | 2×2×2 至 5×5×5 | m | 窄脉需更小 |
| 组合长度 | 1–2 | m | 孔内化验组合 |
| 变异函数变程(Cu) | 50–200 | m | 矿床特定 |
| 搜索椭圆长轴 | 1.5–2.0 × 变程 | m | Kriging 邻域 |
| 每块最小样品数 | 8–16 | — | 控制资源 |
| 每块最大样品数 | 24–32 | — | 限制平滑 |
| 模型更新周期 | 每季度 | — | 每轮钻探后 |
公式
- 块体体积:
(m³) - 吨位:
- 金属量:
(% 品位金属) - 块体估算方差:
— Kriging 系统输出,驱动分类
工具与标准
- 常用软件:Leapfrog Geo、Surpac(GeoVIA)、Datamine Studio RM、Micromine、QGIS(QC)
- 相关标准:JORC 2012、NI 43-101、GB/T 17766(矿产资源估算)、ISO/IEC 17025(化验 QA)
操作步骤
- 导入钻孔:验证孔口、测斜闭合、化验 QA/QC(标准样、空白、重复样)。
- 定义域:岩性 + 矿化包络;复杂褶皱用隐式建模器。
- 化验组合:域内 1–2 m;95 百分位封顶异常值。
- 变异函数:每域每元素实验变异函数;拟合块金值、基台值、变程。
- 块体设置:父块体;脉状子分块(若支持)。
- 品位估算:普通 Kriging(或概算 IDW);条带图与交叉验证。
- 资源分类:至样品距离、Kriging 方差、地质连续性 → 测定/控制/推断。
- 报告:吨位、品位、金属量;分域表;导出 Whittle(.blk)或 CSV。
- 版本控制:每轮钻探后标注日期更新。
知识延伸
地质建模误差向下游放大:未封顶的高品位化验值可将品位晕染至数百块体;块体尺寸远大于 SMU 会使 Whittle 高估可选性与 NPV 5–15%。域间渗漏(搜索跨地质边界)是贫化的模型根源——搜索椭圆须按连续性设硬边界或软限制。
隐式与显式:褶皱 BIF 或复杂断层网络适合 Leapfrog 快速迭代;简单脉状矿体用 Surpac 剖面仍高效。密度默认 2.7 t/m³ 对所有岩性可导致 5–10% 吨位误差,须实验室测定分岩性填入。
与 1.5 资源估算 分类规则一致是合规基础;块体导出至 Whittle 前须 QC 字段名与单位。
常见误区
- 未封顶高品位化验 → 品位晕染
- 域界搜索渗漏 → 矿体贫化
- 块体过大 → Whittle NPV 不现实
- 忽视密度 → 吨位误差
- 无条带图验证 → 变异函数错误未识别
关联章节
自测要点
- Cu 矿变异函数变程 120 m,指定 Kriging 搜索椭圆与最小/最大样品数。
- 褶皱 BIF 矿床选隐式还是显式工作流的决策依据。
- 块体 45 Mt @ 0.65% Cu 与组合校验 38 Mt @ 0.72% Cu,差异诊断路径。
8.1.2 矿山设计与调度软件
8.1.2 矿山设计与调度软件
学习目标
- 以符合实际的成本、回收率、边坡参数运行 Whittle 矿坑优化
- 在 Vulcan/Surpac 设计露天台阶、井下采场与开拓
- 在 Deswik 构建衔接采矿速率与选厂产能的生产计划
- 接口块体模型与 FMS 现场品位控制
核心概念
| 术语 | 定义 |
|---|---|
| Whittle 矿坑壳 | Lerchs-Grossmann 优化矿坑轮廓 |
| 嵌套壳 | 收益因子系列矿坑壳 |
| 收益因子(RF) | 金属价缩放;RF 1.0 = 基准情形 |
| 边界品位 | 入选与排废最低品位 |
| 采矿计划 | 分时段采剥计划:吨数、品位、剥采比 |
| SMU | 最小采矿可选性单元 |
| 品位控制模型 | 爆孔化验更新的短程模型 |
| ERP/MES 集成 | 企业计划与矿山执行衔接 |
参数与指标
| 参数 | 典型值 | 单位 | 说明 |
|---|---|---|---|
| Whittle 边坡角 | 38–55 | ° | 分 sector 岩性 |
| 采矿成本(露天) | 2–8 | USD/t moved | 铲运全含 |
| 选矿成本 | 8–25 | USD/t ore | 选矿 + 管理费用 |
| 回收率 | 85–95 | % | 冶金试验 |
| 折现率 | 8–12 | % | Whittle NPV |
| 计划步长 | 月至季 | — | Deswik 甘特图 |
| SMU(露天) | 5³ 至 10³ | m | 匹配铲斗、爆破孔网 |
| 品位控制更新 | 每周 | — | 爆孔 → FMS |
公式
- 块体价值(Whittle):
- 剥采比:
- NPV(矿坑壳):
- 边际边界品位:
— % Cu 转 USD/lb
工具与标准
- 常用软件:Whittle、Vulcan(Maptek)、Surpac、Deswik Suite、Minemax Scheduler
- 相关标准:JORC 修正因素、公司技术报告规范
操作步骤
- 导入块体模型:校验块体尺寸、密度、品位字段。
- Whittle 参数:成本、回收率、价格、分区边坡角、采矿/选矿能力限制。
- 嵌套壳:RF 0.3–1.3;选 RF 0.85–1.0 实用矿坑(通道 + 岩土)。
- 矿坑分期:推帮 ≤ 30 m 台阶;运输道路 25 m 宽、10% 坡度。
- 井下:采场、开拓、通风节点(Vulcan/Surpac UG)。
- Deswik 甘特图:工序先后;选厂产能与车队规模约束。
- 导出运营数据:矿坑 DXF → FMS;计划吨数/品位 → ERP。
- 对账:月度模型与实际吨数、品位、剥采比。
知识延伸
Whittle 优化矿坑必须与 Slide/FLAC 岩土安全系数交叉校验——过陡边坡的 NPV 壳不可直接采用。运输深度限制与卡车车队匹配是计划可行性的常见约束;忽视检修停机可高估年产量 10–15%。
单一价格情景不是战略——嵌套壳 + 敏感性分析是投资委员会最低要求。品位控制模型与地质块体模型的漂移在 6 个月内可使入选品位超容差,须每周爆孔更新短程模型。
常见误区
- Whittle 边坡过陡 → 岩土不满足
- 忽视运输能力 → 计划不可行
- 块体未更新 → 入选品位偏离
- 单一价格情景
- 计划零停机时间
关联章节
自测要点
- 0.55% Cu、8000 USD/t、采矿 3 USD/t、选矿 12 USD/t ore、90% 回收率,估算边际边界品位。
- 计划 20 Mt/a 矿石 vs 选厂 18 Mt/a,Deswik 约束设置。
- 块体模型导出 → FMS 品位控制导入的数据交接清单。
8.1.3 岩土分析软件
8.1.3 岩土分析软件
学习目标
- Slide 2D LEM(Bishop、Spencer)边坡稳定分析
- RS2(2D FEM)、FLAC3D(3D FEM)模拟开挖应力/变形
- 解读安全系数与设计验收准则
- 将边坡参数反馈 Whittle 与露天设计
核心概念
| 术语 | 定义 |
|---|---|
| 安全系数 FoS | 抗滑力/滑动力;静力目标 > 1.3 |
| LEM | 条分法稳定分析;Bishop、Janbu、Spencer |
| FEM | 连续介质应力应变 + 本构模型 |
| Mohr-Coulomb | c + σ tan φ 线性弹塑性 |
| Hoek-Brown | 脆性岩体非线性强度 |
| 孔隙水压力 | 地下水降低有效应力、安全系数 |
| Unwedge | 井下顶板楔体解析稳定 |
| 地层反应曲线 | 隧道收敛与支护压力关系(FEM) |
参数与指标
| 参数 | 典型值 | 单位 | 说明 |
|---|---|---|---|
| FoS(静力,运营) | ≥ 1.3 | — | 验收标准 |
| FoS(静力,设计) | ≥ 1.5 | — | 长期 |
| FoS(拟静力) | ≥ 1.1 | — | 地震 |
| Bishop FoS(典型) | 1.2–1.8 | — | 最终矿坑边坡 |
| 岩体粘聚力 | 0–500 | kPa | 取决于 GSI |
| 内摩擦角 | 25–45 | ° | 岩体 |
| K₀ | 0.5–1.5 | — | σ_h/σ_v |
| FLAC3D 网格(近开挖面) | 2–5 | m | 边界更粗 |
公式
- Hoek-Brown(岩体):
- Bishop FoS:圆弧/非圆弧滑面力矩平衡迭代解
- 孔隙压力(ru):
- 支护压力:地层反应曲线与支护刚度交点
工具与标准
- 常用软件:Slide 6、RS2、FLAC3D、Unwedge、DIPS、Slope Radar(SSR)
- 相关标准:ANCOLD、McMahon 1985 边坡验收、GB 50021(岩土勘察)
操作步骤
- 数据:构造编录、实验室剪切、地下水、边坡 DXF。
- 材料:各岩性 Mohr-Coulomb 或 Hoek-Brown;孔隙水压力模型。
- Slide LEM:自动搜索临界滑面;Bishop + Spencer;报告最小 FoS。
- 敏感性:φ ± 5°、c ± 50%、水位 + 5 m。
- RS2/FLAC3D:分步开挖;塑性区范围。
- Unwedge:隧道交叉口锚杆长度/间距。
- 建议:可接受边坡角或支护 → Whittle/矿坑设计。
- 监测:FoS < 1.5 边坡布设 SSR/棱镜;位移速率触发。
知识延伸
LEM 使用**完整岩石 φ = 55°**是最危险错误——岩体 GSI 折减后 FoS 可从 3.0 降至 1.1。未模拟地下水的干坡 FoS = 1.5 在饱和条件下可至 1.1;测压管监测与设计须联动。
二维剖面不具代表性时须多剖面分析——最弱剖面控制。复杂 FEM 模型无现场监测验证则仅为学术练习。地震、暴雨渐进破坏须在风险台账中明确分析。
常见误区
- LEM 用完整岩石参数
- 忽视地下水
- 单一二维剖面代表全场
- FEM 无现场验证
- 仅静力 FoS
关联章节
自测要点
- 50° 边坡风化花岗岩(c=25 kPa,φ=32°)FoS=1.22,设计变更建议。
- LEM 与 FEM 对 30 m 高帮、10 m 台阶 + 5 m 平台比较框架。
- 三组节理交汇,Unwedge 锚杆布置要点。