4. 网格化实战

最小曲率法（Minimum Curvature）是 Oasis montaj 中处理磁测数据网格化时最常用的插值算法之一。它的核心思想是：在满足所有测点数值约束的前提下，寻找一个“弯曲能量”最小的光滑曲面。直观地说，它既不像最近邻法那样产生阶梯状块斑，也不像反距离加权那样在测点稀疏处出现“牛眼”状同心圆，而是在数据点之间生成连续、平滑且尽量贴合原始观测的曲面。对于磁异常这种既有区域背景又有局部尖峰（如矿体引起的正负伴生磁异常）的数据，最小曲率法能在保真与平滑之间取得较好的平衡。

在 Oasis montaj 中启动网格化，通常通过 Grid and Image → Gridding → Minimum Curvature 进入。整个实战过程可以拆成四个相互牵制的环节：选择输入通道与坐标、设定网格间距、配置张量（Tension）与迭代参数、确定边界外推策略。下面逐一展开，并给出一个可复现的教学示例。

一、输入准备与网格间距的确定

网格化的输入是数据库中的 X、Y 坐标通道和磁异常通道（例如经日变改正、IGRF 去除后的 TMI 或化极后的 RTP）。在 Oasis montaj 中，坐标通常以地图单位（米）存储，需确认 X 为东向、Y 为北向，且投影一致。网格间距（Cell Size）不是随意填写的数字，它直接决定网格文件的大小与信息保真度。经验规则是：网格间距取测点平均点距的 1/3～1/2。若测网为 100 m × 50 m（线距 × 点距），则点距 50 m，网格间距可取 20～25 m。间距过大，局部磁异常被抹平；间距过小，则产生大量插值“伪细节”，文件体积膨胀却无真实信息增益。

教学示例：假设某测区范围 X 从 500000 到 505000 m，Y 从 3000000 到 3004000 m，共 5000 m × 4000 m，点距约 50 m。取网格间距 20 m，则网格规模为 250 × 200 个节点，约 5 万个单元，属于轻量级、便于反复试算的规模。

二、张量与迭代参数的含义

最小曲率法在 Oasis montaj 中通过两个关键参数控制曲面行为：Tension（张量/张力） 和 Iterations（迭代次数）。Tension 取值范围通常为 0～1，它控制曲面在数据点之间的“绷紧”程度。Tension 接近 0 时，曲面趋向于最小弯曲的薄板样条，非常平滑，但可能在陡变处过冲；Tension 接近 1 时，曲面被拉向数据点之间的线性连接，更贴合局部变化，但可能出现轻微棱角。对磁异常数据，常用 0.25～0.5 作为起点：区域场平滑、局部异常又不过度振荡。

Iterations 控制松弛迭代的收敛次数。次数太少，曲面未充分收敛，残差偏大；次数过多，收益递减且耗时增加。一般 100～300 次即可满足教学与常规生产需求。实际判断收敛与否，应观察软件输出的最大残差（Max Residual）是否已降到可接受量级（例如小于数据动态范围的 1%），而不是盲目加大迭代次数。

三、边界外推策略

测区边缘往往没有数据约束，最小曲率法必须对边界外节点做外推，否则网格边缘会出现剧烈下坠或翘起。Oasis montaj 提供若干边界处理选项，常见的有：None（不外推）、Constant（常数延拓）、Linear（线性外推） 以及基于边缘趋势的延拓。对磁异常网格，推荐使用线性外推或基于边缘梯度的延拓，使边界外的场值沿边缘趋势平缓延伸，避免成图时出现“边框效应”。若后续要做延拓、导数等频率域处理，边界外推的质量尤其重要，因为不合理的边界会在 FFT 中引入虚假高频。

四、一次完整的试算与结果判读

按上述示例参数执行：输入 TMI 通道，网格间距 20 m，Tension = 0.3，Iterations = 150，边界采用 Linear 外推。运行后软件会输出网格文件（如 TMI_MC.grd）及残差统计。判读时重点看三处：其一，最大残差是否在合理范围，若残差过大说明迭代不足或 Tension 过高导致曲面无法贴合数据；其二，在 Map 中叠加测点位置，检查局部磁异常峰谷是否被“削平”或出现虚假旁瓣；其三，观察测区四边，边界外推段是否与内部场值平滑衔接，有无明显折角或跳变。

若发现局部异常被过度平滑，可将 Tension 从 0.3 调到 0.45 再试；若边缘出现翘曲，则检查外推方式是否误设为 None。每次只改一个参数并保留网格文件命名（如加后缀 _t30、_t45），便于对比。最终交付的网格应满足：异常形态与原始剖面一致、边界平缓、残差可接受，三者兼顾后再进入成图与后续处理环节。