1.5 资源量估算

1.5.1 地统计学基础

学习目标
  • 计算实验变异函数并拟合理论模型（球状、指数、高斯）
  • 解释变程（相关长度）、基台值、块金值——联系地质连续性与采样质量
  • 应用普通克里金（OK）估值，以克里金方差 σ²ₖ 作为局部不确定性指标
  • 诊断高块金效应源于采样误差还是真实短程变异
核心概念
术语 定义
变异函数 γ(h) 间距 h 样品对半方差均值，量化空间相关性随距离衰减
变程 a 达到基台值的距离，代表地质连续性尺度，与钻孔间距匹配检验
块金 C₀ 原点处截距，含微尺度变异与测量误差
基台 C₀+C 平台方差，近似样品总方差
普通克里金 OK 无偏线性估计，权重和为 1，最常用的品位估值法
克里金方差 局部估计不确定性，JORC 分类与风险量化输入
块金效应 短程变异远大于长程变异，常见于 Au、U 等矿种
指示克里金 对高于截止品位的概率估值，适合高次变异矿床
交叉验证 留一法或 K 折验证，评估估值无偏性与精度
参数与指标
参数 典型值 单位 说明
滞后距 钻孔间距 25–50% m 100 m 线距常用 25 m 滞后
滞后容差 滞后宽 50% m 方位容差 ±22.5°
每滞后最少对数 30–50 — 稳定实验变异函数
块金/基台比（中等） 0.25–0.50 — 斑岩 Cu 典型
块金/基台比（高） >0.75 — 检查采样与复合
搜索椭圆长轴 1.5–2.5 × 变程 m 各向异性沿走向
每块最多样品数 12–24 — 避免过度平滑
公式
  • 实验变异函数：γ∗(h)=(1/2N)×Σ[Z(xi)−Z(xi+h)]2\gamma *(h) = (1/2N) \times \Sigma [Z(x_{i}) - Z(x_{i} + h)]^{2}
  • 球状模型：γ(h)=C0+C×[1.5(h/a)−0.5(h/a)3]\gamma (h) = C_{0} + C \times [1.5(h/a) - 0.5(h/a)^{3}]（h≤a）；h>a 时为 C₀ + C
  • 指数模型：γ(h)=C0+C×[1−exp⁡(−3h/a)]\gamma (h) = C_{0} + C \times [1 - \exp(-3h/a)]
  • OK 估计：Z∗(x0)=Σλi×Z(xi)Z*(x_0) = \Sigma \lambda ᵢ \times Z(xᵢ)，约束 Σλi=1\Sigma \lambda ᵢ = 1
  • 克里金方差：σk2(x0)=C(0)−Σλi×C(x0,xi)−μ\sigma ^{2}ₖ(x_0) = C(0) - \Sigma \lambda ᵢ \times C(x_0, xᵢ) - \mu
工具与标准
  • 常用软件/仪器：Supervisor（Snowden）；Isatis.neo；Leapfrog Edge；Surpac
  • 相关标准：JORC Code 2012；GB/T 17766（固体矿产资源储量分类）；NI 43-101
操作步骤
  1. 数据检验：复合等长；去聚类；离群值审查（±3σ 或同行评议）；按岩性/蚀变分域
  2. 变异函数：全向与定向实验曲线 → 拟合块金、基台、变程、各向异性 → 变异函数云图验证
  3. 估值设置：块模型范围与 SMU → 搜索策略 → 选择 OK 或 SK（简单克里金）及漂移
  4. 验证：条带图、交叉验证（克里金效率）、剖面对比 → 系统偏差则迭代变异函数
知识延伸

地统计学不是黑箱插值，而是对「空间相关性假设」的显式建模。块金过高可能源于样品长度不匹配、化验 QA/QC 问题、或真实金块效应；需回到 1.4.3 追溯。变程小于钻孔间距时，外推块体不确定性极高，不应分类为 Indicated。

各向异性在斑岩矿床中常见：水平变程远大于垂向，反映层控富集或超基因带。使用全向变程会导致垂向过度连续、高估可采厚度。指示克里金、多指标克里金可在高次变异性矿床中改进分类，但 JORC 披露需说明方法选择与验证结果。

方法对比：ID² 反距离法简单但对数据聚类敏感；OK 在数据充足时更稳健且提供方差。行业惯例以 OK 为主估计，ID² 为基准对比；两者吨位差异 >5% 需调查域划分与离群值处理。验证环节应包含全局（吨位、平均品位）与局部（高品位壳、矿化边界）两个尺度。

变异函数模型选择
模型 形状特征 适用地质 注意
球状 近原点线性增至变程 均质浸染状 最常用，JORC 易接受
指数 近原点曲率大 渐变过渡带 变程定义为 95% 基台距离
高斯 近原点平缓 非常均质体 过平滑风险
空洞效应 原点凹陷 规则间距数据 提示数据聚类或复合问题

块金诊断流程：① 检查样品长度与地质界线对齐；② 审查 QA/QC（CRM/重复）；③ 分析离群值（直方图、箱线图）；④ 对比不同滞后距块金是否稳定；⑤ 若块金/基台 >0.75 且 QA/QC 合格，接受为高块金矿种特征。交叉验证指标：克里金效率 KE=1−σk2/σd2ataKE = 1 - \sigma ^{2}ₖ/\sigma ^{2}_data，KE >0.3 为可接受；平均误差应接近 0（无系统偏差）。

Top-cut（品位封顶）原则：对高于 P99 或地质师认定的异常值封顶，防止高品位孤岛扭曲变异函数与吨位。封顶值选择须文档化，敏感性分析展示封顶前后吨位变化。不封顶时 LG 优化与资源报告均可能显著高估。

常见误区
  • 未等长复合即计算变异函数
  • 强各向异性矿床使用全向变程
  • 忽视高块金仍分类为 Measured
关联章节
自测要点
  1. 块金/基台比偏高时，应依次排查哪些数据与地质因素？
  2. 垂向变程显著小于水平变程时，对块模型与露天台阶设计有何含义？

1.5.2 三维地质建模

学习目标
  • 在 Leapfrog Geo 建立隐式矿体模型，或在 Surpac/Datamine 手工线框固实
  • 应用域边界：截止品位、夹石、最小可采宽度、贫化壳
  • 块单元尺寸取钻孔间距 1/4–1/3；赋予密度与岩石类型码
  • 对比 ID² 与 OK 结果，JORC 报告披露差异
核心概念
术语 定义
隐式建模 Leapfrog 径向基函数面，由钻孔接触面自动生成，适合复杂接触
线框实体 Surpac 剖面手工圈矿传统方法，透明度高、劳动密集
估值域 地质与地统计均质的估算分区，跨域估值引入平滑偏差
截止品位 CoG 矿/岩分类最低品位，与 NSR 经济模型联动
SMU 最小可采单元，块模型父块尺寸参考
贫化壳 矿岩接触外扩薄壳，模拟采矿不可避免夹石
参数与指标
参数 典型值 单位 说明
块尺寸（水平） 25–50 m 对应 100 m 孔距
块尺寸（垂向） 5–15 m 匹配台阶高度
截止品位（斑岩 Cu） 0.20–0.30 %Cu 项目 NSR 特定
最小可采宽度（地下） 5–10 m 露天 15–30 m
内部夹石吞并阈值 3–5 m 薄夹石带并入矿体
矿石密度 2.7–3.2 t/m³ 实测为主
Leapfrog 面分辨率 1–5 m 接触带加密
公式
  • NSR 截止品位：CoG=(Cmining+Cproc+G&A)/(Price×Recovery×Payability)CoG = (C_{mining} + C_{proc} + G\&A) / (Price \times Recovery \times Payability)
  • 吨位：T=Σ(Vi×Di)T = \Sigma (V_{i} \times D_{i})，V_i 块体积，D_i 密度
  • 金属量：M=T×Gavg×RecoveryM = T \times G_{avg} \times Recovery
  • ID² 权重：λi=1/di2/Σ(1/dj2)\lambda ᵢ = 1/dᵢ^{2} / \Sigma (1/dⱼ^{2})
工具与标准
  • 常用软件/仪器：Leapfrog Geo / Edge；Surpac；Datamine Studio RM；Vulcan
  • 相关标准：JORC Code 2012 Table 1 估算方法；GB/T 17766
操作步骤
  1. 域建模：导入录井 → Leapfrog 岩性/蚀变面 → 品位壳或地质接触圈矿域
  2. 块模型：旋转网格至矿体主轴 → 设父块 → 接触带子块加密
  3. 估值：复合至块支持 → OK 主估、ID² 基准 → 应用 CoG 与最小宽度规则
  4. 验证：OK vs ID² 吨位/品位差异 <5%；剖面视觉检查；JORC Table 1 格式报告
知识延伸

地质模型与块模型应「同胚」：线框或隐式面定义域边界，块模型仅在域内估值。隐式建模对数据稀疏区外推敏感，需结合地质师解释限制曲面曲率。密度未实测而采用 2.7 t/m³ 默认值可导致 5–8% 吨位误差，对融资报告与剥采比优化均为重大风险。

截止品位不是静态常数：应随金属价格、汇率、选冶回收与成本变动做敏感性分析。内部夹石处理规则（薄夹石吞并厚度）显著影响可采吨位与平均品位，须在技术报告中明确并与采矿设计贫化假设一致。

与露天优化的接口：块模型导出至 Whittle 或 Surpac Pseudoflow 前，应完成品位封顶（top-cut）审查与密度赋值。未封顶的高品位异常会在 LG 算法中产生「高值孤岛」，人为抬高境界。矿体线框质量直接决定 2.1.1 境界圈定 可信度。

域划分决策依据
域类型 划分依据 估值策略 风险
岩性域 录井岩性码 分域独立变异函数 域界模糊时跨域泄漏
蚀变域 蚀变类型/强度 分域估值 蚀变渐变区难划分
矿化样式域 浸染/脉状/块状 不同变程 样式过渡带
氧化带域 氧化/过渡/硫化物 不同回收率假设 氧化带界线不确定

隐式建模（Leapfrog）优势：数据稀疏区可自动生成合理曲面，减少手工线框主观性。劣势：外推区曲面可能过度光滑，须地质师设置趋势约束与曲率限制。线框建模（Surpac）优势：地质师完全控制矿体边界，透明度高。劣势：劳动密集，复杂接触带耗时。行业趋势是 Leapfrog 建地质域 + Surpac/Edge 估值的混合工作流。

常见误区
  • 默认密度未实测 → 系统性吨位偏差
  • 块尺寸等于孔距 → 过度平滑
  • 未封顶异常值构建矿壳 → 高品位吨位膨胀
关联章节
自测要点
  1. 地质域划分错误如何导致 OK 与 ID² 系统性偏差？
  2. NSR 截止品位与地质截止品位在应用上有何区别？

1.5.3 资源分类

学习目标
  • 一致应用 JORC（2012）、NI 43-101、GB/T 17766 分类框架
  • 依据孔距、QA/QC、估计置信度将 Inferred 升级至 Indicated、Measured
  • 记录采矿损失（5–15%）与贫化（5–25%）及依据
  • 计算经济截止品位，整合金属价格、回收率与运营成本
核心概念
术语 定义
Inferred 低置信度，不确定性约 ±20–30%，不能单独支撑可行性研究
Indicated 中等置信度，约 ±10–15%，可支撑预可行性研究
Measured 高置信度，约 ±5–10%，可支撑可行性研究与储量转换
修正因素 采矿、选冶、经济、环境参数，Resource → Reserve 必备
矿石储量 Ore Reserve 经济可采的 Measured/Indicated 资源子集
经济截止品位 边际块 NPV=0 的品位，高于地质截止品位
参数与指标
参数 典型值 单位 说明
孔距 Inferred 100–200 m 斑岩 Cu
孔距 Indicated 50–100 m JORC 典型
孔距 Measured 25–50 m 生产钻探
采矿损失（露天） 5–10 % 设计未采到
贫化（露天） 5–15 % 接触夹石
贫化（地下崩落） 10–25 % 过冒与充填混入
价格假设（Cu） LME 三年均或共识 $/t JORC 要求合理依据
公式
  • 经济盈亏平衡品位：CoG=OPEX/(Price×Recovery×(1−Penalty))CoG = OPEX / (Price \times Recovery \times (1 - Penalty))
  • 储量吨位（简化）：TR=TResource×(1−Loss)×(1−Dilution)T_{\mathrm{R}} = T_{\mathrm{Resource}} \times (1 - Loss) \times (1 - Dilution)（贫化亦可增加废石吨）
  • 分类距离规则（典型）：25 m 内 ≥2 孔 → Measured；75 m 内 ≥2 孔 → Indicated；否则 Inferred
  • 价格披露：名义与实际价格须透明记录，JORC 禁止无依据现货价
工具与标准
  • 常用软件/仪器：Surpac 资源报告；Leapfrog Edge 分类壳；CoG 敏感性电子表
  • 相关标准：JORC Code 2012；GB/T 17766-2020；NI 43-101；SAMREC
操作步骤
  1. 分类编码：按孔距、σ²ₖ、地质连续性为块体赋 Measured/Indicated/Inferred（JORC 第 19 条）
  2. CoG 敏感性：CoG ±20% 跑境界或采场优化 → 记录吨位/品位敏感度
  3. 修正因素：损失、贫化、回收、计价条款 → 合格 Indicated/Measured 转 Probable/Proved Reserve
  4. 报告：按类别列表资源与储量；JORC Table 1 披露全部假设；CP 签字
知识延伸

资源分类是「地质知识+数据密度+估计不确定性」的综合判断，不能简化为孔距查表。JORC 强调地质连续性：Inferred 块体紧邻钻孔但若位于未验证的构造错断带，仍不应升级。克里金方差高区对应高估计风险，应与钻探工程优先级联动。

储量转换引入修正因素后，Probable/Proved 才是矿山融资与设备采购依据。贫化与损失假设须与 2.1 或 2.2 采矿方法匹配：崩落法贫化显著高于充填法。CoG 对运营支出敏感，在成本通胀周期应重新运行敏感性，避免经济矿体被错误剔除。

与项目经济的衔接：9.2.1 项目经济学 的 NPV 与 CoG 互为约束。分类升级不应受融资时间表压力而提前，CP 个人法律责任与项目声誉风险远高于短期资本成本。

JORC 分类决策要素
要素 Measured Indicated Inferred
孔距（斑岩 Cu） 25–50 m 50–100 m 100–200 m
地质连续性 已验证 合理推断 有限推断
估计方法 OK + 验证 OK/ID² + 验证 OK/ID²
σ²ₖ 相对水平 低 中 高
支持研究 可行性 预可行性 不可单独支撑 FS
钻孔数/块体 ≥3 孔穿过 ≥2 孔穿过 ≥1 孔

地质截止品位（CoG_geol）由矿体圈定需求决定，通常低于经济截止品位（CoG_econ）。JORC 资源量可用地质 CoG 圈定矿体壳，但储量转换必须用经济 CoG 并经过修正因素。CoG 敏感性分析应展示 ±20% 变动下吨位与金属量变化，是技术报告必备内容。

Resource → Reserve 转换须应用修正因素：采矿损失（露天 5–10%、地下 10–15%）、贫化（露天 5–15%、崩落 10–25%）、选冶回收率（须基于冶金试验）、计价条款（NSR 公式）。仅 Measured/Indicated 可转 Proved/Probable Reserve；Inferred 不得转换。

常见误区
  • 有邻近钻孔无地质连续性即标 Indicated
  • 用现货价定 CoG → JORC 不合规
  • 储量报表忽略贫化/损失 → 可采吨位高估
关联章节
自测要点
  1. JORC 中地质连续性与钻孔间距在分类决策中分别起什么作用？
  2. 资源量（Resource）与储量（Reserve）的本质区别是什么？