2. 搜索半径的基本概念与地质统计基础

在矿产资源量估算中，“搜索半径”并不是一个随手设定的数字，而是把地质统计模型与工程可采性连接起来的关键参数。它决定了估值时一个待估块段会“看”到哪些样品、以多大权重看待这些样品，以及最终估值结果在空间上的平滑程度。理解搜索半径，必须同时理解三个层次：它在几何上如何定义、在三维空间中如何被各向异性椭球体塑形、以及它与变异函数变程之间的定量关系。

一、搜索半径的定义与几何含义

搜索半径（search radius）指以待估块段中心（或块段质心）为原点，在给定方向上能够纳入参与估值的样品所允许的最大距离。它通常以“主方向半径 + 次方向半径 + 垂向半径”三个分量给出，而不是单一标量。原因在于矿体在空间上极少各向同性：沿走向可能连续数百米，沿倾向数十米，垂直方向上仅数米至十余米。若用单一球体半径，要么在走向上纳入过多样品造成过度平滑，要么在垂向上纳入不到足够样品导致估值失败。

从几何上看，搜索半径定义的是一个“搜索邻域”（search neighborhood）。邻域内样品数量、样品与块段的距离、样品之间的相对位置共同决定克里金权重。搜索半径过小，邻域内样品不足，估值方差增大、块段出现“空洞”；搜索半径过大，邻域跨越不同矿化域，权重被远处不相关样品稀释，估值趋于区域平均值，局部变化被抹平。因此搜索半径本质上是“偏差—方差”权衡的工程化表达。

二、各向异性椭球体：把半径放进真实地质空间

实际工作中，搜索邻域通常实现为一个各向异性椭球体。椭球体由三个半轴长度（a、b、c）和一组旋转角（绕 Z 轴方位角、绕 X 轴倾角、绕 Y 轴侧伏角）确定。三个半轴分别对应矿化的最大连续方向、中等连续方向和最短连续方向。建模时通常先由变异函数在各方向的变程比确定半轴比例，再结合勘探网度与矿体规模确定绝对尺度。

椭球体的旋转必须与地质坐标系一致。例如某斑岩铜矿体走向北东 40°、倾向南东、倾角约 55°，则椭球体长轴应旋转到走向方向，中轴沿倾向，短轴近似垂直矿体。若旋转角设置错误，即使半轴长度合理，邻域也会“斜切”矿体，把不同矿化带的样品混入同一估值块段。这是搜索半径实践中最常见、也最隐蔽的错误之一。

三、搜索邻域：不只是“半径”，还有数量与分布约束

工程实现中，搜索邻域通常由“椭球体 + 样品数量上下限 + 每象限/每扇区最大样品数”共同约束。仅给半径而不限制数量，会在样品密集区纳入过多近邻、在稀疏区纳入过少远邻。常见做法是设定最少样品数（如 4～8 个）和最多样品数（如 12～24 个），并把椭球体划分为若干扇区（如 8 个水平扇区 × 3 个垂向分层），限制每个扇区最多取若干样品，以保证权重在空间上分布均衡，避免“一侧样品扎堆、另一侧无样品”的偏倚。

当第一次搜索未达到最少样品数时，系统通常按比例放大椭球体（如 1.5 倍、2 倍）进行第二轮、第三轮搜索，直至满足条件或达到最大搜索距离。这一“多轮搜索”机制必须在报告中说明，因为它直接影响估值块段的样品来源与可靠性等级。

四、与变异函数变程的定量关系

变异函数变程 a 表示空间相关性消失的距离：滞后距小于 a 时，样品间存在可建模的相关性；大于 a 时，协方差趋近于零。搜索半径若显著大于变程，纳入的远处样品权重接近零，对估值贡献很小，却会增加计算量并可能引入不同矿化域的样品；若显著小于变程，则丢弃了本可利用的相关信息，估值方差偏大。

实践中常用的经验区间是：搜索半径取变程的 0.6～1.0 倍，主方向可接近 1.0a，次方向约 0.7a～0.8a，垂向约 0.5a～0.7a，具体比例由各方向变程比确定。例如某金矿主方向变程 120 m、倾向变程 80 m、垂向变程 30 m，则可设椭球半轴 a=110 m、b=75 m、c=28 m，比例约 1 : 0.68 : 0.25，与变程比 1 : 0.67 : 0.25 基本一致。这一设置既尊重了地质统计模型，又保留了工程上可接受的局部分辨能力。

五、参数设置示例与判断依据

下表给出一个教学示例（参数为合理假设，非实测结果），用于说明搜索半径与变程、样品数量约束之间的配套关系：

设置完成后，应通过交叉验证观察：估值与实测值的平均误差是否接近零、误差方差是否与克里金方差匹配、是否存在系统性高估或低估。若误差随空间位置呈明显趋势，往往提示椭球体旋转角或半轴比例与矿化方向不一致，需要回到变异函数模型重新核对，而不是简单调大或调小搜索半径。

综上，搜索半径是变异函数模型在估值空间中的“落地形态”：变程给出相关性的尺度，各向异性椭球体把这一尺度旋转、拉伸到真实地质方向，样品数量与扇区约束则保证邻域内的信息分布合理。三者配套，才能使搜索半径既反映地质统计规律，又满足资源量估算的工程可靠性要求。