4. 网格化方法选择与参数设置

构造等值线图的“底子”是网格。地震解释得到的层位数据是沿测线分布的离散点（解释层位、断层多边形边界、井分层校正点等），网格化就是把这些不规则散点插值成规则网格，再由网格追踪等值线。网格化方法选错，等值线会出现虚假闭合、锯齿、过冲或“牛眼”；参数设错，则会把真实的小幅度构造抹平，或把噪声放大成假构造。本节围绕反距离加权、克里金、最小曲率三类常用算法，说明其适用性与关键参数的设置原则。

三类算法的本质差异

反距离加权（IDW）以距离的幂次倒数作为权重，认为离待估点越近的已知点影响越大。它计算快、结果稳定、不会产生超出数据范围的极值，但本质是“加权平均”，天然具有平滑效应，对局部极值和陡坎不敏感，容易在数据稀疏区形成以已知点为中心的圆形“牛眼”假象。适合数据较密、构造平缓、只需快速成图的场景。

克里金基于变异函数（半变差函数）建模数据的空间相关性，分两步：先用实验变异函数拟合理论模型（球状、指数、高斯等），得到块金值、基台值、变程三个参数，再做最优无偏估计。它能给出估计方差，便于评估网格可靠性，对数据分布不均、存在各向异性的工区适应性强。代价是需要足够多的数据点才能稳定拟合变异函数，参数不当会出现“过度平滑”或“过拟合”。

最小曲率把网格面看作受数据点约束的薄板，求满足数据点且弯曲能量最小的曲面，等价于求解双调和方程。它生成连续光滑、无“牛眼”的曲面，特别适合区域构造趋势明显、需要平滑过渡的等值线，是构造图常用的“美观型”算法。但它对局部突变（小断层、小高点）响应较弱，且迭代求解对网格规模和松弛参数敏感。

网格间距：先定“分辨率预算”

网格间距决定等值线能表达的最小构造尺度。原则是：网格间距应小于目标构造最小宽度的一半，同时与地震测网密度匹配。若三维工区面元为 12.5 m×12.5 m、解释线距 25 m，网格取 10～12.5 m 通常足够；若二维测网线距 500 m，网格取 100～250 m 即可，再细只会放大插值噪声。判断标准很直接：把网格间距减半重做一次，若等值线形态与闭合圈数量基本不变，说明当前间距已足够；若出现大量新“小圈”，多为噪声被放大，应回退。

搜索半径：控制“借点”的范围

搜索半径决定每个网格节点从多远范围内取已知点参与插值。半径太小，节点周围无点，出现空洞或被迫外推；半径太大，会把断层另一盘的点“借”过来，造成等值线跨断层连续、构造被拉平。实用做法是统计工区有效解释点的平均点距 d，取 R≈2～3d，并保证圆内通常有 4～8 个点。对断层附近应单独处理：先建立断层多边形作为插值屏障，使搜索不跨越断层，再对断层两侧分别网格化，最后拼接。

平滑因子与迭代参数

最小曲率类算法通常有松弛因子（relaxation）与迭代次数。松弛因子过小收敛慢、易残留锯齿；过大则迭代震荡甚至发散。建议从 1.0 起，按 1.0→1.2→1.4 逐级试算，观察残差（网格值与原始点值之差）是否稳定下降并收敛。平滑因子用于压制短波长噪声，但会压低真实构造幅度，应遵循“先不平滑成图、确认构造真实存在、再加最小必要平滑”的顺序，避免用平滑掩盖解释质量问题。

方法选择与参数配置示例

下表为某二维工区（线距 500 m，有效解释点平均点距约 400 m，目标为识别幅度 20 ms 以上的低幅度圈闭）的合理假设配置，仅作教学示例，实际需按工区试验确定。

网格化完成后，应把网格回代到原始解释点，统计残差分布：残差应近似零均值、无系统性偏移；若断层附近残差偏大，说明搜索跨越了断层；若局部残差呈规律性正负交替，多为网格过粗或平滑过度。只有残差与等值线形态同时通过检查，网格化结果才可作为闭合检查的可靠输入。