3. 组合样长度不统一的成因与影响

组合样(Composite Sample)是把相邻若干单个样品的品位按长度加权合并后得到的代表值,其基本前提是“长度可加、品位按长度加权”。一旦参与组合的样品长度参差不齐,这个前提就被破坏,误差会从单样品层面一路传导到变异函数和资源量估算。理解成因与影响链条,是后续“长度统一”和“异常值处理”两步操作的必要铺垫。

一、长度不统一的四类来源

实际矿山数据库中,组合样长度不一致很少是单一原因造成的,常见来源可归为四类:

  • 采样方法差异。刻槽样常按固定长度(如1 m或2 m)连续采取;岩芯样受回次长度控制,回次长度由钻探工艺和岩层破碎程度决定,常见0.3~3 m不等;拣块样、全巷样则根本没有“长度”概念,只能折算为等效长度。
  • 矿体边界截断。当样品跨越矿体顶底板时,只能保留矿体内的那一段,形成“残样”,例如设计2 m的样在边界处只剩0.6 m。
  • 岩芯采取率与破碎带。破碎带、断层泥、溶蚀带岩芯采取率低,回次被拆分,长度被切碎;同一钻孔内软硬岩层交替,回次长度天然不均。
  • 数据录入与历史数据拼接。不同勘查阶段、不同施工单位、不同规范版本的数据合并入库,长度单位(m/ft)、起止深度记录方式不一致,也会造成名义长度不统一。

这四类来源中,前两类属于“客观必然”,后两类属于“可管理误差”。区分来源很重要:前者只能通过组合规则消化,后者应优先在数据库层面清洗。

二、对品位加权平均的偏差

设某组合区间内各样品的长度为 Li、品位为 gi,正确的组合品位应为:

G = Σ(L_i · g_i) / Σ L_i

若误用等权算术平均 G' = Σgi / n,则偏差为:

G' - G = Σ g_i/n - Σ(L_i·g_i)/ΣL_i

偏差方向取决于“长样与短样的品位谁高”。举例(示例数据,非实测):某组合段含三个样品,长度分别为1.0 m、0.5 m、2.5 m,品位分别为1.2 g/t、3.0 g/t、0.8 g/t。长度加权品位为 (1.0×1.2+0.5×3.0+2.5×0.8)/4.0 = (1.2+1.5+2.0)/4.0 = 1.175 g/t;等权平均为 (1.2+3.0+0.8)/3 = 1.667 g/t,偏高约42%。原因很清楚:高品位的短样被赋予了与长样相同的权重。

更隐蔽的问题是“组合样长度本身不一致”。若组合样A长2 m、组合样B长5 m,在后续统计中若对A、B再作等权平均,等于默认2 m与5 m代表相同体积,同样引入偏差。因此长度统一不仅是单样品层面的事,也是组合样层面的事。

注意:长度加权只对“沿矿体厚度方向连续采样”成立。对于沿走向或倾向的采样,长度应替换为对应的线密度权重或面积/体积权重,不能直接套用长度公式。

三、对变异函数与资源量估算的连锁影响

长度不统一的影响不是孤立的,它会沿“样品→组合样→变异函数→估值→资源量”逐级放大。下图展示了这条传导链及关键控制点。

组合样长度不统一的误差传导链 采样方法差异 刻槽/岩芯/拣块 矿体边界截断 残样、短样 采取率与破碎带 回次拆分 数据拼接 单位/深度记录 组合样长度不一致 等权平均 → 品位系统性偏高/偏低 变异函数失真 块金效应增大、变程缩短、各向异性扭曲 资源量估算偏差 克里格权重失衡、局部高估/低估 误差反馈放大 控制点① 统一组合长度 控制点② 异常值识别与处理

具体来看,变异函数失真有三个典型表现。第一,块金效应被抬高:短样与长样代表不同体积,其品位波动中混入了“体积差异”这一非地质因素,实验变异函数在原点附近迅速上升,块金值C₀被人为放大。第二,变程被压缩:长度不一的样品在空间上分布不均,短样密集区与长样稀疏区的信息密度不同,拟合出的变程往往偏短,导致估值过于“局部化”。第三,各向异性被扭曲:若沿矿体厚度方向长度差异大、沿走向差异小,各向异性椭球的长短轴比例会被错误估计。

这些失真最终反映到资源量上。以普通克里格为例,权重由变异函数决定;块金增大使克里格权重趋于平均化,估值被“拉平”,高品位段被低估、低品位段被高估,总体资源量可能变化不大但品位-吨位曲线形态改变;若变程被压缩,则估值过度依赖近邻样品,局部出现“牛眼”效应,高品位样品的影響被放大,资源量可能系统性偏高。对于品位分布偏态的矿床(如金、铜),这种偏差尤为敏感。

四、实践中的判断与处理思路

面对长度不统一,第一步不是急着组合,而是先做长度分布诊断:统计单样品长度的频数分布、均值、变异系数,并分钻孔、分岩性、分勘查阶段对比。若长度变异系数小于约0.3,等权与长度加权差异通常有限;若大于0.5,则必须严格按长度加权,并考虑统一组合长度。

第二步是确定组合长度。常用做法是取样品长度的中位数或众数作为组合长度(如2 m),对超长样进行劈分、对残样按实际长度保留并单独标记。劈分时按长度比例分配品位,残样不参与等权统计,避免短样被过度加权。

第三步是在变异函数拟合前复核:检查实验变异函数在短距离处的形态是否异常陡升,对比长度加权前后块金值的变化幅度。若块金值变化超过30%,应回到组合环节排查,而不是直接调整拟合参数掩盖问题。

判断口诀:先看长度分布,再定组合长度;先加权再拟合,先复核再估值。长度不统一是“输入问题”,不应靠“输出端调参”来补偿。