2. 重磁异常数据特征与预处理实战

向上延拓和边界增强算法对输入数据的“干净程度”极其敏感：一条未压制的测线噪声，经延拓后会被放大成虚假的深部异常；一个错误的网格间距，会让导数类边界增强结果整体偏移。因此，在进入位场变换之前，必须完成一条完整的数据预处理链：坐标与投影统一 → 网格化 → 去噪 → 区域场与局部场分离 → 质量检查 → 异常幅值标定。本节按这条链的实际操作顺序展开，并给出可复现的示例参数。

一、重磁异常的数据特征：先看清数据再动手

重磁观测数据有三个区别于一般地球物理数据的特征，直接决定了预处理策略。第一是测点分布不规则：重力测点常沿道路、村镇布设，磁测可能为航空测网，点距从几十米到几公里不等，直接做傅里叶变换会因采样不均产生假频。第二是量级跨度大：布格重力异常可达数百 mGal，而局部构造引起的异常可能只有 0.5～2 mGal；磁异常从几 nT 到上万 nT 都有，弱异常极易被区域背景淹没。第三是噪声来源多样：重力有仪器漂移、固体潮、地形改正残差；磁测有日变、飞行高度起伏、人文电磁干扰。这些噪声在空间域表现为高频抖动或单点跳变，在波数域表现为高频能量抬升。

因此，预处理的目标不是“把数据变好看”，而是在保留有效异常频谱的前提下压制干扰，并为后续延拓提供等间距、无空值、量纲统一的规则网格。

二、网格化：从离散测点到规则网格

网格化是把不规则测点插值到规则网格的过程，核心参数是网格间距和插值方法。网格间距的经验取法是取平均点距的 1/2～1/4；间距过大丢失局部异常，过小则引入大量插值假象。例如某测区平均点距 200 m，可选 50 m 网格。插值方法上，位场数据推荐最小曲率法或克里金法，前者对平滑场保形好，后者可给出估计方差用于质量评估；反距离加权在点稀疏区易产生“牛眼”假异常，应慎用。

网格化后必须检查两项：一是空值区，测区边缘外推区域应做掩膜标记，避免延拓时边界效应向内部扩散；二是网格方向与投影，所有数据统一到同一投影坐标系（如 CGCS2000 高斯投影），并记录网格原点、间距、行列数，供后续波数域计算使用。

三、去噪：空间域与波数域配合

去噪分两步走。第一步在空间域做单点跳变剔除：对每条测线计算相邻点差值，超过 3 倍中误差的点判为可疑，结合日变改正记录和飞行高度数据人工复核后剔除或替换。第二步在网格域做低通或正则化滤波：对网格数据做滑动平均（窗口 3×3 或 5×5）压制高频噪声，或用维纳滤波在波数域按信噪比加权。需要注意，低通滤波的截止波数不能低于后续边界增强所需的最短波长，否则会把要提取的局部构造信息一并抹掉。

四、区域场与局部场分离

区域场是深部或大范围地质体产生的宽缓背景，局部场是浅部目标引起的相对短波长异常。分离方法常用三种：多项式拟合（1～3 阶趋势面，适合区域场形态简单的情形）、滑动平均法（窗口大于目标异常波长 2～3 倍）、匹配滤波或小波多尺度分解（按波数或尺度分离，适合复杂叠加场）。以多项式拟合为例，对网格节点做 2 阶趋势面最小二乘拟合，残差即为局部场；阶数过高会把局部异常“吃掉”，阶数过低则残留区域背景。

分离结果必须做闭合检查：区域场 + 局部场应能逐点还原原始网格，残差仅为数值误差量级。这一步是后续向上延拓和边界增强可信度的前提。

五、数据质量检查：用可量化的指标说话

质量检查不能停留在“看起来平滑”。建议至少做四项量化检查：①空值率，网格内空值节点占比应低于 5%，否则延拓边界效应显著；②功率谱形态，去噪后高频段能量应明显下降而中低频保持连续，若出现异常尖峰说明存在周期性干扰；③闭合差，区域场与局部场之和与原始网格的均方根误差应小于原始异常幅值的 0.1%；④剖面抽查，沿主构造走向和垂直走向各取一条剖面，人工核对异常形态是否与已知地质、钻孔或已知矿点吻合。四项中任一项不通过，都应回到对应环节调整参数，而不是直接进入延拓。

六、异常幅值标定：让不同批次数据可比

幅值标定解决的是“同一测区不同时间、不同仪器、不同高度观测的数据能否拼接使用”的问题。重力数据以国家重力基准点为基准做绝对点值校正和仪器格值校正，并统一到同一高程基准；磁力数据以日变站记录做日变改正，以国际地磁参考场（IGRF）做正常场改正，航空磁测还需做飞行高度归算。标定完成后，所有网格必须使用同一单位（重力 mGal、磁力 nT）、同一基准面、同一投影，并在元数据中记录标定日期、基准点编号和改正模型版本。

完成上述六步后，得到的规则网格才具备进入向上延拓与边界增强计算的条件。判断预处理是否合格的最终标准不是某一步的参数多漂亮，而是延拓结果是否稳定、边界增强结果是否与已知地质边界可对比——预处理的质量，最终要在位场变换结果中接受检验。