3. 地质数据库构建与质量评价

资源量估算的可靠性,首先取决于进入估算流程的那张“数据底表”是否干净、一致、可追溯。地质数据库并不是把钻孔编录表简单堆进一个文件,而是把孔口信息、测斜、岩性、蚀变、矿化、采样、化验、密度、工程坐标等多源数据,按统一的空间参照和编码体系组织起来,并保留每一次修改的痕迹。本节围绕钻孔数据管理、样品组合、特高品位处理、完整性检查与误差识别五个环节展开,重点说明“实际怎么做”以及“做到什么程度算合格”。

一、钻孔数据管理:以孔为骨架的多表关联

一个规范的地质数据库通常以“孔号”为主键,拆分为孔口表(collar)、测斜表(survey)、岩性/蚀变表(geology)、采样表(sample)、化验表(assay)、密度表(density)等。孔口表记录开孔坐标(X、Y、Z)、终孔深度、孔向、施工日期;测斜表记录不同深度的倾角与方位角;采样表记录样号、起止深度、样长、采样方法与岩性归属;化验表记录样号对应的元素品位。多表之间通过孔号与深度区间建立关联,任何一条化验记录都必须能回溯到具体的采样段和空间位置。

实际管理中最容易出问题的是坐标与深度参照。必须事先明确:X、Y采用何种坐标系与投影带,Z是相对标高还是绝对标高,深度是沿孔轴测量还是垂直深度。以某斑岩铜矿示例(以下数据均为教学假设)为例,孔口坐标为X=500 000 m、Y=3 200 000 m、Z=1 250 m,孔深600 m,开孔方位角90°、倾角-60°。若测斜显示孔深300 m处倾角变为-55°、方位角92°,则后续样品的空间位置必须按测斜分段计算,而不能用开孔方向一插到底。数据库应保留原始测斜记录,并由软件按最小曲率法或线性插值生成三维孔迹,供后续块体赋值使用。

二、样品组合:把不等长样变成等长样

野外采样受岩性边界、矿化边界和采样规范限制,样长往往不等,常见0.5~2.0 m。而品位估值通常要求等长组合样,以消除长样对品位的过度加权。组合的核心是“长度加权平均”:在给定组合长度(如2.0 m)下,沿孔轴把连续样品切分为等长区间,区间内品位按样长加权。

示例:某孔0~6 m内样品为0~1.0 m品位0.42%、1.0~2.5 m品位0.55%、2.5~4.0 m品位0.38%、4.0~6.0 m品位0.61%(均为假设)。若组合长度为2.0 m,则0~2.0 m组合品位为(1.0×0.42+1.0×0.55)/2.0=0.485%;2.0~4.0 m为(0.5×0.55+1.5×0.38)/2.0=0.4225%;4.0~6.0 m为0.61%。组合时必须处理“跨岩性/跨矿化边界”的问题:通常要求在矿化域或岩性域内部组合,不让组合样跨越域边界,否则会把废石品位混入矿石。组合后还应保留组合样长度、原始样号列表和加权权重,便于审计。

三、特高品位处理:识别、截切与敏感性

特高品位(极端值)会显著拉高局部估值,尤其在样品数量少、矿化不均匀的矿床中。处理不是简单“删掉”,而是先识别再决策。常用识别方法包括:按矿化域统计品位分布,以均值加若干倍标准差、或上四分位数加1.5倍四分位距作为阈值;也可用累积频率曲线取99%分位。示例中某金矿矿化域样品均值1.2 g/t、标准差2.8 g/t,取均值加3倍标准差约9.6 g/t作为初选阈值,超过者列为待处理样品。

处理方式主要有三种:截切(capping,把超过阈值的品位替换为阈值)、限制影响范围(限制高品位样品的搜索半径或赋值权重)、以及保留原值但单独做敏感性分析。截切阈值不能拍脑袋,应通过对比不同阈值(如8、10、12 g/t)下资源量变化来判断:若阈值从12降到8导致金属量下降超过约10%,说明估值对极端值敏感,需要在报告中披露并谨慎取值。所有被截切的样品应记录原值、截切值、阈值依据和操作人。

四、数据完整性检查与误差来源识别

完整性检查应形成固定清单并逐项留痕:孔口坐标是否落在矿区范围内、孔深是否与测斜和采样最大深度一致、样品区间是否重叠或存在未采样空段、化验样号是否与采样表一一对应、品位是否出现负值或超出仪器检出限的异常、密度数据是否覆盖主要岩性。下表给出一份可直接套用的检查示例(数值为教学假设)。

检查项检查规则示例发现处理建议
样品区间重叠同孔相邻样品起止深度不得重叠ZK1203在45.0~46.0 m与45.5~47.0 m重叠核对原始采样记录,修正边界
未采样空段矿化域内不应存在大于组合长度的空段ZK1207在88~92 m无样品补采或标注为废石并说明
样号缺失化验表样号必须存在于采样表3个化验样号无对应采样记录追溯送样单,补录或剔除
品位异常品位为负或超过合理上限1个样品Cu品位为-0.02%判定录入错误,回查原始化验单
测斜缺失孔深大于100 m应有不少于3条测斜ZK1211仅开孔1条测斜补测或按直孔假设并标注风险

误差来源可分为四类:一是测量误差,如坐标、孔深、测斜的系统偏差;二是采样与化验误差,包括样长不均、采样污染、化验精度与检出限;三是录入与关联误差,如样号错配、深度单位混用(米与英尺)、区间重叠;四是地质解释误差,如岩性边界与矿化域圈定不一致。识别误差的关键是“交叉验证”:用测斜计算的孔底坐标与相邻孔见矿标高对比,用组合前后品位分布对比,用不同阈值下的资源量对比。只有把误差来源写清楚,后续的估算方案对比和不确定性评价才有依据。

钻孔数据库:从原始表到可估算数据底表 孔口表 collar 测斜表 survey 岩性/蚀变表 采样表 sample 化验表 assay 密度表 density 关联主键 孔号 + 深度区间 统一坐标系 统一深度参照 样号唯一 ① 完整性检查 重叠/空段/样号/负值 ② 样品组合 2.0 m 长度加权 ③ 特高品位处理 识别→截切→敏感性 ④ 误差来源识别 测量/采样/录入/解释 可估算数据底表 等长组合样 + 截切记录 三维孔迹 + 域编码 检查日志 + 修改留痕 示例(教学假设):ZK1203 组合与截切 0–1.0 m 0.42% | 1.0–2.5 m 0.55% | 2.5–4.0 m 0.38% | 4.0–6.0 m 0.61% 2.0 m 组合:0–2.0 m = 0.485% | 2.0–4.0 m = 0.4225% | 4.0–6.0 m = 0.61% 特高品位阈值示例:均值1.2 + 3×标准差2.8 ≈ 9.6 g/t,超阈值样品记录原值与截切值
建议把上述检查表固化为数据库的“入库前校验脚本”:每次导入新钻孔数据自动运行,输出异常清单并阻断不合格记录进入估值流程。检查日志与截切记录应随资源量报告一并归档,作为不确定性评价的证据链。

总体而言,地质数据库构建与质量评价的目标不是“数据看起来齐全”,而是每一条进入估算的样品都能说明来源、经过何种处理、处理依据是什么。钻孔数据管理保证空间与关联正确,样品组合保证权重合理,特高品位处理控制极端值影响,完整性检查与误差识别则把潜在偏差显性化。这四步做扎实,后续的估算方案对比才有共同的可信基础。