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 限制平滑
模型更新周期 每季度 每轮钻探后
公式
  • 块体体积V=X×Y×ZV = X \times Y \times Z(m³)
  • 吨位T=V×density(t/m3)T = V \times density (t/m^{3})
  • 金属量M=T×grade/100M = T \times grade / 100(% 品位金属)
  • 块体估算方差σk2riging\sigma ^{2}_kriging — Kriging 系统输出,驱动分类
工具与标准
  • 常用软件:Leapfrog Geo、Surpac(GeoVIA)、Datamine Studio RM、Micromine、QGIS(QC)
  • 相关标准:JORC 2012、NI 43-101、GB/T 17766(矿产资源估算)、ISO/IEC 17025(化验 QA)
操作步骤
  1. 导入钻孔:验证孔口、测斜闭合、化验 QA/QC(标准样、空白、重复样)。
  2. 定义域:岩性 + 矿化包络;复杂褶皱用隐式建模器。
  3. 化验组合:域内 1–2 m;95 百分位封顶异常值。
  4. 变异函数:每域每元素实验变异函数;拟合块金值、基台值、变程。
  5. 块体设置:父块体;脉状子分块(若支持)。
  6. 品位估算:普通 Kriging(或概算 IDW);条带图与交叉验证。
  7. 资源分类:至样品距离、Kriging 方差、地质连续性 → 测定/控制/推断。
  8. 报告:吨位、品位、金属量;分域表;导出 Whittle(.blk)或 CSV。
  9. 版本控制:每轮钻探后标注日期更新。
知识延伸

地质建模误差向下游放大:未封顶的高品位化验值可将品位晕染至数百块体;块体尺寸远大于 SMU 会使 Whittle 高估可选性与 NPV 5–15%。域间渗漏(搜索跨地质边界)是贫化的模型根源——搜索椭圆须按连续性设硬边界或软限制。

隐式与显式:褶皱 BIF 或复杂断层网络适合 Leapfrog 快速迭代;简单脉状矿体用 Surpac 剖面仍高效。密度默认 2.7 t/m³ 对所有岩性可导致 5–10% 吨位误差,须实验室测定分岩性填入。

1.5 资源估算 分类规则一致是合规基础;块体导出至 Whittle 前须 QC 字段名与单位。

常见误区
  • 未封顶高品位化验 → 品位晕染
  • 域界搜索渗漏 → 矿体贫化
  • 块体过大 → Whittle NPV 不现实
  • 忽视密度 → 吨位误差
  • 无条带图验证 → 变异函数错误未识别
关联章节
自测要点
  1. Cu 矿变异函数变程 120 m,指定 Kriging 搜索椭圆与最小/最大样品数。
  2. 褶皱 BIF 矿床选隐式还是显式工作流的决策依据。
  3. 块体 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)BV=(grade×recovery×priceprocessingcost)×tonnesminingcost×tonnesBV = (grade \times recovery \times price - processing_{cost}) \times tonnes - mining_{cost} \times tonnes
  • 剥采比SR=Wastetonnes/OretonnesSR = Waste_{tonnes} / Ore_{tonnes}
  • NPV(矿坑壳)NPV=Σ(Cashflowt/(1+r)t)NPV = \Sigma (Cashflow_{t} / (1 + r)^t)
  • 边际边界品位CoG=Processingcost/(Recovery×Price×22.046)CoG = Processing_{cost} / (Recovery \times Price \times 22.046) — % Cu 转 USD/lb
工具与标准
  • 常用软件:Whittle、Vulcan(Maptek)、Surpac、Deswik Suite、Minemax Scheduler
  • 相关标准:JORC 修正因素、公司技术报告规范
操作步骤
  1. 导入块体模型:校验块体尺寸、密度、品位字段。
  2. Whittle 参数:成本、回收率、价格、分区边坡角、采矿/选矿能力限制。
  3. 嵌套壳:RF 0.3–1.3;选 RF 0.85–1.0 实用矿坑(通道 + 岩土)。
  4. 矿坑分期:推帮 ≤ 30 m 台阶;运输道路 25 m 宽、10% 坡度。
  5. 井下:采场、开拓、通风节点(Vulcan/Surpac UG)。
  6. Deswik 甘特图:工序先后;选厂产能与车队规模约束。
  7. 导出运营数据:矿坑 DXF → FMS;计划吨数/品位 → ERP。
  8. 对账:月度模型与实际吨数、品位、剥采比。
知识延伸

Whittle 优化矿坑必须与 Slide/FLAC 岩土安全系数交叉校验——过陡边坡的 NPV 壳不可直接采用。运输深度限制与卡车车队匹配是计划可行性的常见约束;忽视检修停机可高估年产量 10–15%。

单一价格情景不是战略——嵌套壳 + 敏感性分析是投资委员会最低要求。品位控制模型与地质块体模型的漂移在 6 个月内可使入选品位超容差,须每周爆孔更新短程模型。

常见误区
  • Whittle 边坡过陡 → 岩土不满足
  • 忽视运输能力 → 计划不可行
  • 块体未更新 → 入选品位偏离
  • 单一价格情景
  • 计划零停机时间
关联章节
自测要点
  1. 0.55% Cu、8000 USD/t、采矿 3 USD/t、选矿 12 USD/t ore、90% 回收率,估算边际边界品位。
  2. 计划 20 Mt/a 矿石 vs 选厂 18 Mt/a,Deswik 约束设置。
  3. 块体模型导出 → 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(岩体)σ1=σ3+σci(mbσ3/σci+s)a\sigma_{\mathrm{1}} = \sigma_{\mathrm{3}} + \sigma_{\mathrm{ci}} (m_{b} \sigma_{\mathrm{3}}/\sigma_{\mathrm{ci}} + s)^a
  • Bishop FoS:圆弧/非圆弧滑面力矩平衡迭代解
  • 孔隙压力(ru)ru=u/(γ×h)ru = u / (\gamma \times h)
  • 支护压力:地层反应曲线与支护刚度交点
工具与标准
  • 常用软件:Slide 6、RS2、FLAC3D、Unwedge、DIPS、Slope Radar(SSR)
  • 相关标准:ANCOLD、McMahon 1985 边坡验收、GB 50021(岩土勘察)
操作步骤
  1. 数据:构造编录、实验室剪切、地下水、边坡 DXF。
  2. 材料:各岩性 Mohr-Coulomb 或 Hoek-Brown;孔隙水压力模型。
  3. Slide LEM:自动搜索临界滑面;Bishop + Spencer;报告最小 FoS。
  4. 敏感性:φ ± 5°、c ± 50%、水位 + 5 m。
  5. RS2/FLAC3D:分步开挖;塑性区范围。
  6. Unwedge:隧道交叉口锚杆长度/间距。
  7. 建议:可接受边坡角或支护 → Whittle/矿坑设计。
  8. 监测:FoS < 1.5 边坡布设 SSR/棱镜;位移速率触发。
知识延伸

LEM 使用**完整岩石 φ = 55°**是最危险错误——岩体 GSI 折减后 FoS 可从 3.0 降至 1.1。未模拟地下水的干坡 FoS = 1.5 在饱和条件下可至 1.1;测压管监测与设计须联动。

二维剖面不具代表性时须多剖面分析——最弱剖面控制。复杂 FEM 模型无现场监测验证则仅为学术练习。地震、暴雨渐进破坏须在风险台账中明确分析。

常见误区
  • LEM 用完整岩石参数
  • 忽视地下水
  • 单一二维剖面代表全场
  • FEM 无现场验证
  • 仅静力 FoS
关联章节
自测要点
  1. 50° 边坡风化花岗岩(c=25 kPa,φ=32°)FoS=1.22,设计变更建议。
  2. LEM 与 FEM 对 30 m 高帮、10 m 台阶 + 5 m 平台比较框架。
  3. 三组节理交汇,Unwedge 锚杆布置要点。