4. 品位直方图与分布形态判读
品位直方图是把清洗后的单样品位按区间归并后统计频数得到的分布图,它回答的核心问题不是“平均品位是多少”,而是“品位在空间上如何铺开”。对铜矿而言,一个样品组合的品位分布往往不是理想的正态钟形:它可能右偏(低品位样品占多数、少数高品位样品拉长右尾)、可能尖峰厚尾,也可能出现两个甚至三个明显的峰。这些形态特征直接对应矿化期次、岩性组合与蚀变强度的差异,是分级方案能否贴合矿体实际的第一道判据。
一、从清洗数据到频数表:直方图的构造
绘制前先明确数据口径:使用清洗后的单样品位(已剔除超限、重复、空值与明显录入错误),单位统一为 %Cu,样品长度不等时宜按长度加权频数,避免长样品过度代表。区间划分建议采用等宽分箱,箱宽按数据跨度选取,常用 0.05%、0.1% 或 0.2% Cu;箱数控制在 12~25 个之间,过少掩盖多峰,过多则被随机噪声干扰。以某铜矿 1200 个样品、品位 0.05%~2.60% Cu 为例,取箱宽 0.1%,得到如下频数结构(示例数据,用于说明方法):
| 品位区间(%Cu) | 频数 | 累计频率 | 形态提示 |
|---|---|---|---|
| 0.0~0.3 | 96 | 8.0% | 低品位背景 |
| 0.3~0.6 | 288 | 32.0% | 第一峰主体 |
| 0.6~0.9 | 252 | 53.0% | 第一峰峰顶附近 |
| 0.9~1.2 | 168 | 67.0% | 谷区 |
| 1.2~1.5 | 180 | 82.0% | 第二峰上升 |
| 1.5~1.8 | 132 | 93.0% | 第二峰峰顶 |
| 1.8~2.1 | 60 | 98.0% | 右尾 |
| 2.1~2.6 | 24 | 100% | 高品位稀疏尾 |
把上表画成直方图,可以看到 0.6%~0.9% 与 1.5%~1.8% 两个频数高点,中间 0.9%~1.2% 出现明显低谷。这就是典型的双峰形态,而不是单峰右偏。
二、偏度与峰度的判读口径
偏度衡量分布的不对称方向:偏度为正表示右偏(长尾拖向高品位),为负表示左偏。峰度衡量尾部厚度与峰顶尖锐程度,常用超额峰度(正态为 0),大于 0 表示尖峰厚尾,小于 0 表示平顶。判读时不要只看一个数,应结合直方图形态与分位数一起看。上例数据计算得偏度约 +0.85、超额峰度约 +1.6,说明整体右偏且尖峰厚尾;但右偏主要来自 1.8% Cu 以上的高品位尾,而非主体分布本身。若误把整体右偏当作单一对数正态处理,就会把 1.2%~1.8% 这一段真实存在的第二峰“抹平”,导致高品位矿石被并入低级别。
三、多峰分布对应什么:矿化期次与岩性差异
多峰不是统计噪声,通常意味着样品总体由两个或多个地质子总体混合而成。常见解释有三类:其一,两期矿化叠加,早期浸染状黄铜矿化形成 0.5%~0.9% Cu 的低峰,晚期沿裂隙充填的脉状黄铜矿—斑铜矿形成 1.4%~1.8% Cu 的高峰;其二,岩性差异,细粒闪长岩中铜矿物分散、品位偏低,而接触带矽卡岩或角砾岩筒品位偏高;其三,蚀变分带,绢英岩化带与青磐岩化带品位系统不同。判别方法是把直方图按岩性、标高或矿化类型分组重绘:若分组后每组趋于单峰,则混合总体假设成立;若分组后仍为双峰,则需考虑采样尺度或分析误差。
四、据形态调整分级方案
分级边界应尽量落在谷区或分布的低频段,而不是切在峰顶。上例原方案以 1.0% Cu 为界划分低品位矿与高品位矿,恰好切在第二峰上升段,会把大量 1.0%~1.3% Cu 的脉状矿化样品划入低级别,造成高品位资源被低估。依据双峰形态,可把边界下移至谷区中部约 1.1%~1.2% Cu,或上移至第二峰起点约 1.3% Cu,具体取决于选矿回收试验对边界品位的经济论证。调整后建议形成三段式分级:边界品位 0.3% Cu、低品位矿 0.3%~1.3%、高品位矿 ≥1.3%,并对 ≥1.8% Cu 的右尾单独设高品位亚类,便于配矿。
调整后必须复核两件事:一是各分级吨位与平均品位是否随边界移动发生合理变化,可用累计频率曲线读取边界处的吨位占比;二是分组直方图是否支持“两个子总体”的判断,若按岩性分组后低品位组与高品位组各自单峰且峰位与总体双峰吻合,则分级方案有地质依据。若分组后仍重叠严重,应回到采样代表性与样品长度加权环节复查,而不是继续微调边界。
综上,直方图判读的落脚点始终是分级方案:峰的数量决定要不要分组分级,峰的位置决定边界往哪里放,谷区决定边界不能放在哪里。把这三条与地质分组、经济边界品位结合,品位分级才不是拍脑袋的数字游戏。