3. 样品组合与品位代表长度
在进入变异函数拟合与克里金权重求解之前，必须先回答一个基础问题：参与计算的每一个“样品点”，其品位值究竟代表多长一段岩芯或多少个刻槽米？原始采样往往长度不一——有的样长0.8 m，有的1.5 m，有的因矿化界线卡样而只有0.4 m。若直接把不同长度的样品当作等权点投入变异函数计算，短样与长样被赋予相同权重，会系统性放大短样的局部波动，使实验变异函数在近距离被虚假抬高。样品组合（compositing）正是为解决这一问题而设的预处理环节：把不等长原始样品按规则合并成等长的组合样品，使每个组合品位代表一个固定的“代表长度”，从而让后续统计与插值在同一尺度上进行。
组合的原理：长度加权而非简单平均
组合的核心是长度加权平均。设某钻孔在区间内有一组原始样品，样长为 li、品位为 gi，则组合品位为：
G = Σ(l_i × g_i) / Σ(l_i)注意分母是参与组合的总长度，而不是样品个数。举例说明：某孔在 0～3.0 m 内有三段样品，分别为 0.5 m@1.2 g/t、1.0 m@2.0 g/t、1.5 m@1.6 g/t。若简单算术平均得 (1.2+2.0+1.6)/3 = 1.60 g/t；而长度加权为 (0.5×1.2 + 1.0×2.0 + 1.5×1.6)/3.0 = (0.6+2.0+2.4)/3.0 = 1.667 g/t。两者相差约 4%，在品位分布偏斜时差异会更明显。这就是必须坚持长度加权的原因。
组合长度的确定
组合长度不是随意选取的，它需要与三个因素匹配：一是采样长度，组合长度通常取原始样长的整数倍或接近其统计中位数，避免频繁切割；二是矿体厚度与开采单元，组合长度应小于最小可采厚度，且与选冶、采矿的台阶高度或回采进尺相协调；三是变异函数的结构尺度，组合长度应明显小于变异函数的变程，否则会把短尺度结构“抹平”。工程实践中，若原始样长中位数为 1.0 m，常见组合长度取 1.0 m 或 2.0 m；对于厚大矿体、样长 2 m 的钻孔，可取 2 m 或 4 m。本例后续统一采用 2.0 m 作为组合长度进行演示。
残留样品的处理规则
当孔内某一段的总长度不是组合长度的整数倍时，末端会剩下一个不足 2.0 m 的“残留样”。处理方式有三种，需在项目中明确并全程一致：
- 并入前一组：把残留段与相邻的最后一个完整组合段合并，形成一个超长组合样。适用于残留段较长（如≥1.0 m）且矿化连续的情况。
- 独立保留：残留段单独作为一个组合样，长度记为其实际长度，并在数据库中标注“非标准长度”。适用于残留段接近组合长度、或位于矿体边界需保留的情况。
- 剔除：残留段过短（如<0.3 m）且位于矿体外围时舍弃。必须记录剔除长度与位置，防止在品位统计中造成偏差。
组合的操作流程与结果校验
实际操作通常在地质统计软件或数据库中完成，关键步骤是：按孔号与深度排序原始样品→按组合长度从孔口向下累加→判断每个组合区间覆盖了哪些原始样→按长度加权求品位→处理末端残留→输出组合品位表。组合完成后必须做三项校验：组合样总长度与原始总长度是否闭合；组合品位是否落在参与组合的原始品位最小最大值之间（加权平均的必然性质，越界说明计算有误）；组合前后品位均值与方差的变化是否在合理范围——组合会降低方差，但均值应基本不变。
上图以一段 8.4 m 的钻孔为例演示了组合过程：原始五段样品长度分别为 0.5、1.0、1.5、0.8、1.2 m，按 2.0 m 组合后得到四个标准组合样和一个 0.2 m 残留样。可以看到，组合品位 1.73、1.70、1.56、1.56 均落在其参与组合的原始品位极值之间，符合加权平均的性质；组合后总长 8.2 m 加残留 0.2 m 等于原始 8.4 m，长度闭合。需要说明，上述数值为教学示例，实际项目应以真实样品数据计算。
组合对后续变异函数分析的影响
组合后的数据才适合进入变异函数计算。原因在于：实验变异函数依赖“距离—品位差平方”的配对统计，若样品长度不一，则同一距离上的品位差既包含真实的空间变化，也混入了样长差异带来的噪声。等长组合样消除了这一混杂，使变异函数的块金效应、基台与变程更能反映真实的地质结构。同时，组合长度也决定了变异函数能分辨的最小尺度——组合长度为 2.0 m 时，小于 2.0 m 的结构细节已被平滑，滞后距（lag）的起点通常取组合长度或其整数倍，避免在无法分辨的尺度上强行拟合。
完成组合并校验通过后，得到的等长组合品位表就是后续变异函数分析、实验变异函数计算与理论模型拟合的标准输入。此时每个样品点都代表相同的 2.0 m 长度，样品点坐标通常取组合区间的中点深度，这一坐标约定也需在后续所有计算中保持一致。