3. 实验变异函数计算原理与参数选择

实验变异函数是地质统计学建模的起点，它把离散采样点之间的空间相似性转化为一条可量化的曲线。计算过程本身并不复杂：对每一对样本点，求其距离 h 与半方差 γ(h)=½[Z(x)−Z(x+h)]²，再按距离分组求平均。真正决定这条曲线“长什么样”的，是四个参数——滞后距组距、滞后容差、角度容差和最大距离。它们不是软件里的默认值，而是需要根据采样密度、矿化尺度和各向异性方向主动设定的分析决策。

四个参数的含义与作用

滞后距组距（lag spacing）是相邻两个距离分组中心之间的距离，决定了曲线在距离轴上的采样密度。组距过小，每个分组内的点对数量不足，半方差估计噪声大，曲线呈锯齿状；组距过大，曲线被过度平滑，变程和块金效应被掩盖。经验上，组距应使每个分组至少包含 30 对以上的样本点，同时保证在变程范围内有 8～12 个有效分组。

滞后容差（lag tolerance）是每个分组允许的距离波动范围，通常取组距的一半。若组距为 20 m、容差为 10 m，则中心 100 m 的分组实际接收 90～110 m 的点对。容差过小会漏掉大量点对，容差过大则相邻分组重叠，造成曲线人为平滑。当采样在距离上分布不均时，可适当放大容差以“接住”稀疏距离段的数据。

角度容差（angular tolerance）用于方向性变异函数计算。以某一方向为中心轴，只有连线方向落在中心轴±容差范围内的点对才被计入。容差通常取 22.5°（即 45° 扇形）或 30°。容差越小，方向选择性越强，但可用点对越少；容差过大，方向性被稀释，各向异性结构难以显现。判断各向异性时，一般先以 45° 扇形在 0°、45°、90°、135° 四个方向计算，再对变异最大的方向细化。

最大距离（max distance）是参与计算的最远点对距离，一般取研究区最大连续尺度的 1/2～2/3，或取变程的 1.5～2 倍。超过该距离的点对空间相关性已很弱，纳入计算只会增加噪声并拉平曲线尾部。最大距离过小则看不到基台，过大则尾部点对稀少、方差虚高。

参数敏感性：从曲线形态反推参数是否合理

参数选择是否得当，最直接的检验方式是观察实验变异函数曲线的形态。以下判断依据来自教学示例的合理假设，实际项目中应以真实点对统计为准。

当组距过小时，曲线在近距离段剧烈跳动，相邻分组半方差忽高忽低，甚至出现负值或非单调上升。这说明每个分组内点对太少，平均化不足。处理办法是增大组距，或放宽滞后容差把邻近分组的点对并入。

当组距过大时，曲线过于光滑，块金效应被“吞掉”，变程被低估。若理论模型拟合时块金值接近零而实际采样存在明显短程变异，应回头检查组距是否超过了最小采样间距的 2～3 倍。

当角度容差过小时，方向性变异函数在部分距离分组上点对数为零，曲线断裂或缺失。此时应放宽容差至 30°～45°，或合并相邻方向。反之，若四个方向的曲线几乎重合，而地质上已知存在条带状或层状结构，则可能是容差过大稀释了方向差异，应缩小容差重新计算。

当最大距离过大时，曲线尾部出现上翘或剧烈波动，这是因为远距离点对数量骤减、方差估计不稳定。此时应截断最大距离，使尾部每个分组仍保有足够点对。

一组可复现的参数设置示例

假设某沉积型矿床采样间距约 25 m，研究区东西向长约 1200 m，地质上推测矿化沿 90° 方向连续性更强。一组合理的初始参数如下（示例假设，非实测）：

lag_spacing      = 25 m      # 约等于最小采样间距
lag_tolerance    = 12.5 m    # 组距的一半
max_distance     = 500 m     # 约研究区长度的 40%
angular_tolerance= 22.5°     # 45°扇形
directions       = 0°, 45°, 90°, 135°
min_pairs_per_lag= 30        # 低于此值的分组不参与拟合

按此设置先计算全向（omnidirectional）实验变异函数，确认曲线在 500 m 内是否出现基台。若 90° 方向曲线明显高于 0° 方向，则各向异性成立，再以 90° 为主轴、22.5° 容差细化，并逐步调整组距与最大距离，观察变程估计是否稳定。参数敏感性分析的核心不是找一组“唯一正确”的值，而是确认在合理扰动范围内，变程、基台和块金比的估计不发生数量级变化——曲线稳健，后续拟合才有意义。