3. IDW关键参数体系

反距离加权（IDW）的估值结果并不是“算法自动算出来的客观真值”，而是由一组参数共同塑造的模型输出。同一套样品数据，换一组参数，矿体边界、品位高低和储量数字都可能明显变化。因此，理解IDW的参数体系，本质上是在理解“你向算法表达了怎样的地质假设”。本节围绕五个关键参数展开：幂指数 p、搜索椭球体（各向异性与主轴方向）、最少/最多样品数、距离计算方式，以及如何在上机前对参数敏感性做出预判。

幂指数 p：控制权重的集中程度

IDW 的估值公式为：待估点品位等于各已知样品品位按其距离的 −p 次方加权平均。p 越大，近处样品权重越占绝对优势，估值曲面越“贴合”样品点，但样品之间的区域容易出现“牛眼”状尖峰和台阶；p 越小，权重趋于平均，估值更平滑，但会过度抹平高品位矿囊和低品位夹层。工程上常用 p=2（距离平方反比），这也是多数储量估算软件的默认值。

选择 p 时应结合样品密度与矿化变化性：样品稀疏、矿化连续性好时，p 可取 1.5～2；样品密集、矿化局部突变明显时，可提高到 2.5～3，但必须检查是否出现不合理的孤立高值圈闭。一个实用的检验方法是沿一条穿过若干见矿孔的勘探线做剖面估值，观察估值曲线在样品点附近是否出现明显“尖刺”，以及两孔之间是否被压得过于平坦。

搜索椭球体：各向异性与主轴方向

搜索椭球体决定“用哪些样品参与估值”。它由三个半轴长度（主轴、次轴、短轴）和三个方向角（方位角、倾角、倾伏角）定义。各向异性比（主轴/次轴、主轴/短轴）应反映矿体的真实几何形态：层状、似层状矿体沿走向和倾向延伸远、厚度方向短，各向异性比常取 3～8；等轴状矿体则接近 1。若把各向异性比设错，算法会把沿走向远处的样品与垂直方向近处的样品等同对待，导致矿体边界被“拉平”或“拖尾”。

主轴方向必须与矿体最大连续方向一致，通常由地质剖面、矿体三维模型或变差函数分析确定。方向角定义需在全文和软件中保持一致：本文约定方位角为从正北顺时针旋转，倾角为主轴向下为正。方向设反 90°，等价于把矿体走向与厚度方向互换，估值结果会出现系统性偏差，而这种错误在平面图上往往不易察觉，必须用垂直矿体走向的剖面来检查。

最少/最多样品数与距离计算方式

最少样品数保证估值有足够约束，一般设为 3～5；若某待估块体在椭球体内样品不足，应将其标记为“未估值”而不是强行外推。最多样品数用于控制计算量并避免远距离样品稀释局部信息，常取 8～16。两者配合“每象限/每八分圆最多样品数”使用，可防止样品在某一方向过度聚集而其他方向缺失。

距离计算方式直接影响权重：三维直线距离最简单，但对倾斜矿体，沿层距离更符合地质连续性；各向异性距离则先对坐标做方向缩放再求距离，使椭球搜索与距离权重保持一致。若搜索用椭球而距离用直线，会出现“椭球内但权重不合理”的矛盾，建议二者统一采用各向异性距离。

参数敏感性预判：上机前先想清楚

参数敏感性预判不是算完再解释，而是在建模前就判断“哪些参数会主导结果”。可按以下顺序做单因素扰动：先固定椭球与样品数，只变 p；再固定 p，只变各向异性比与主轴方位角；最后变最少/最多样品数。每次记录矿体平均品位、高品位单元体积占比、未估值块体比例三个指标。若 p 从 2 变到 2.5 使平均品位变化超过约 5%，或主轴方位角偏转 15° 就使矿体边界明显移动，这些参数就是本项目的敏感因子，必须用地质剖面、变差函数或交叉验证加以约束。

需要强调的是，上述数值均为教学示例中的合理假设，实际项目应结合样品数量、矿化类型和勘探网度确定。参数没有脱离地质背景的“最优值”，只有与地质认识一致、且经得起剖面检查和交叉验证的“可接受值”。把参数体系当作一套需要解释和验证的地质假设来管理，IDW 的结果才具备用于储量估算的可信度。