4. 地质数据采集与整理实战
剖面法与断面法储量计算的精度,首先取决于进入计算流程的原始地质数据是否规范、可比、可追溯。本节以一条假想的铁铜多金属矿区勘探线整理任务为例,把勘探线布设、样品组合、品位统计和特高品位处理串成一条完整的数据整理链,并给出可直接套用的校验表。示例中的孔号、坐标、品位均为教学假设值,不代表任何真实矿区实测结果。
一、勘探线布设:先定参照系,再定线距
勘探线布设的核心是建立统一的坐标与方位参照。工程上通常以矿区主矿体走向为基准,令勘探线垂直矿体走向布置,线方位角记为 α(从正北顺时针量取)。若矿体走向为 NE35°,则勘探线方位角约为 125°(35°+90°)。全文及图示均采用“正北为 0°、顺时针为正”的约定,后续方位角、倾向描述不得混用其他参照系。
线距的确定要兼顾矿体规模与储量类别要求。以本例为例,矿体走向长约 1200 m,倾向延深约 400 m,采用 100 m 线距布设 13 条勘探线(编号 0~12 线),其中 0、4、8、12 线为首批施工的主干线,用于控制矿体总体形态;其余为加密线。每条勘探线上布置 2~4 个钻孔,孔位偏离勘探线的距离一般不应大于线距的 1/10,本例控制在不大于 10 m,超出时须在孔斜测量中记录偏斜方位与偏距,供剖面投影时归位。
布设完成后应立即形成“勘探线—钻孔—坐标”对照关系,这是后续所有数据挂接的主键。建议字段包括:线号、孔号、孔口坐标 X/Y/Z、终孔深度、测斜成果、开孔方位与倾角。缺少测斜数据的钻孔,在剖面投影时只能按直孔处理,会带来系统性偏差,因此测斜必须作为布设阶段的强制项。
二、样品组合:把不等长样品变成等长样段
岩芯取样常按岩性、矿化程度分段,样长从 0.5 m 到 2.0 m 不等。剖面法计算需要沿孔深方向按固定样长(本例取 2.0 m)组合,才能与剖面面积、块段体积匹配。组合的基本原则是:
- 组合样长与工业指标、矿体圈定精度相匹配,本例取 2.0 m;
- 组合只能在矿体内或同一工业矿层内进行,不得跨越矿体边界、夹石剔除边界;
- 组合品位按样长加权,公式为 C = Σ(Ci·Li) / ΣLi;
- 组合后样段必须保留原始样号清单,保证可追溯。
以 4 线 ZK401 孔为例,某矿段原始样品为:0.8 m@0.42%、1.2 m@0.55%、2.0 m@0.61%、1.5 m@0.48%、0.5 m@0.30%,合计 6.0 m,恰好组合为 3 个 2.0 m 样段。第一段由 0.8 m 与 1.2 m 两个样组成,加权品位为 (0.42×0.8+0.55×1.2)/2.0 = 0.498%;第二段为单一 2.0 m 样,品位 0.61%;第三段由 1.5 m 与 0.5 m 组成,加权品位为 (0.48×1.5+0.30×0.5)/2.0 = 0.435%。若某矿段总长不能被 2.0 m 整除,剩余尾段可单独成样并标注“非标准样长”,在储量计算中按实际长度参与,不得强行并入相邻样段。
三、品位统计:先看分布形态,再选统计量
组合完成后,对全矿区组合样做品位统计。本例共获得 186 个组合样,铜品位算术平均 0.52%,标准差 0.31%,变异系数 0.60,属中等变异。统计时至少应输出:样品数、最小值、最大值、算术平均、标准差、变异系数、分位数(P10/P50/P90)。变异系数大于 1.0 时,说明品位分布明显偏态,算术平均易受高值拉动,应同时给出中位数与截尾均值作为对照。
统计还必须分单元进行。本例按 0~4 线、4~8 线、8~12 线三个块段分别统计,发现 8~12 线平均品位 0.71%,显著高于其他块段,进一步核查发现该块段存在 3 个异常高值样,进入特高品位处理环节。若不分单元统计,这一局部富集会被全区平均值掩盖。
四、特高品位处理:识别、判定、替代、复核
特高品位指明显偏离总体分布、对平均值产生过度拉动的个别品位值。常用判定方法有两种:一是经验倍数法,以平均品位的 6~8 倍作为上限,本例取 6 倍,即 0.52%×6 ≈ 3.1%;二是分位数法,以 P95 或 P97.5 作为阈值。本例 3 个异常样品位分别为 5.8%、4.9%、3.6%,均超过 3.1% 阈值,判定为特高品位。
处理方法常用“替代法”:用特高品位所在单元(或全矿区)的平均品位替代原值,或按“特高品位+相邻若干样”的平均值替代。本例对 5.8% 样采用所在块段平均品位 0.71% 替代,对 4.9% 与 3.6% 样采用相邻 3 样平均品位替代。处理后重新统计,8~12 线块段平均品位由 0.71% 降至 0.58%,全区平均由 0.52% 降至 0.49%。处理前后必须保留原始值、替代值、替代依据三列,任何一次替代都可回溯。
五、整理与校验:一张可填写的校验表
数据整理的最后一步是校验。建议按下表逐项核对,任一项不通过则退回上一步修正,不得带病进入剖面投影与面积计算。
| 校验项 | 校验内容 | 本例结果 | 判定 |
|---|---|---|---|
| 坐标与方位 | 孔口坐标完整,测斜齐全,偏距≤10 m | 13 线 41 孔,测斜齐全,最大偏距 8.6 m | 通过 |
| 样长组合 | 组合样长 2.0 m,未跨矿体边界 | 186 个组合样,尾段 4 个已标注 | 通过 |
| 品位统计 | 分块段统计,变异系数已计算 | 三块段 CV 分别为 0.54、0.57、0.78 | 通过 |
| 特高品位 | 阈值明确,替代有依据,原始值保留 | 3 个样已替代并记录 | 通过 |
| 数据挂接 | 线号—孔号—样号三级可追溯 | 全部挂接,无孤立样 | 通过 |
完成校验后,数据即可进入剖面投影:按勘探线方位将钻孔投影到剖面,组合样按孔深落位,为后续圈定矿体边界、量算剖面面积做好准备。整个整理过程的关键不是某一个公式,而是“参照系统一、组合有边界、统计分单元、替代有记录、校验可回溯”这五条纪律。