大地电磁数据处理中稳健回归估计量的性能与评估

《Minerals》：Robust Regression Estimators in Magnetotelluric Data Processing: Performance and Evaluation

　　

大地电磁（MT）传递函数通常采用统计回归方法进行估计。可靠的MT派生地电结构对于研究深部地球动力学过程与成矿系统至关重要；然而，噪声污染会扭曲阻抗估计并导致误导性的地质解释。在本研究中，研究人员将MM估计量引入MT数据处理，并利用合成算例与中国东北一条MT测线的数据，将其与M估计量和S估计量的性能进行了评估。结果表明，三种估计方法均可有效抑制复杂背景噪声，提高MT阻抗估计的稳定性和可靠性。各方法给出的视电阻率与相位响应总体一致，其中MM估计量得到相对最低的均方根（RMS）值。但在实际应用中应谨慎评估偏差—不确定性权衡，尤其在使用MM估计量处理长周期MT数据时。最后，研究人员建议将三种稳健估计方法之一与远程参考（remote reference）技术结合，以同时抑制预测变量（磁场数据）与因变量（电场数据）中的噪声。

研究背景与动机

大地电磁（MT）法作为常规电磁地球物理方法，广泛用于地壳及岩石圈结构、深部地球动力学、资源勘探与造山变形系统研究。地下导电异常常与岩浆储库、流体迁移通道、部分熔融及控矿断裂弱带相关，因此MT派生电阻率是约束成矿系统的重要依据。然而矿区人文电磁噪声复杂，易扭曲MT传递函数、引入阻抗估计不确定性并降低后续反演质量。传统普通最小二乘（OLS）在残差非高斯时对离群值敏感；尽管M估计量、S估计量等稳健方法已被采用，且在预测变量（磁场）受污染时仍存在杠杆点（leverage point）问题，需结合远程参考技术，但不同稳健估计量在典型统计场景下的表现差异仍需系统评估。为此，研究人员将MM估计量引入MT处理，并在合成与东北野外数据中对比M、S、MM三类估计量。

论文发表于《Minerals》，研究意义在于为MT阻抗估计提供兼具高崩溃点（breakdown point）与高统计效率的实用方案，并明确长周期数据处理中的偏差—不确定性权衡。

主要关键技术方法

研究人员采用0.5 Hz采样约10天的LEMI-423实测磁场驱动均匀半空间（电阻率100 Ω·m）合成电场，叠加高斯、方波、尖峰、锯齿混合噪声（SNR=30 dB与5 dB，污染比例30%/50%/70%）构建合成数据集；野外数据取自2016年4月内蒙古通辽至吉林白山NW—SE向宽频MT剖面中站点520、700及远程参考站500、690（Metronix ADU-07e，电偶极100 m，MFS06磁探头）。处理流程统一为：去均值与预白化、汉明窗分段（512样点，25%重叠）、七级降采样（每级/8）、FFT求谱密度矩阵、相干预选，随后分别用M估计量（Tukey双平方权、MAD尺度、IRLS）、S估计量（多起点最小化残差尺度）、MM估计量（S初估尺度+M细化，并均加入基于帽子矩阵对角元的杠杆权重）估计阻抗张量Z_ij，并以远程参考+M估计量作对照。

研究结果

3.1 合成实验

3.1.1 合成数据：以均匀半空间理论阻抗由实测磁场生成无噪电场，再逆推磁场与垂直场，构造可控干净基数据。

3.1.2 稳健估计结果：SNR=30 dB时，三种方法在周期<80 s短频段均偏离理论值，中长周期视电阻率与相位随污染比30%→70%仍稳定；S估计量高频视电阻率略偏，MM估计量RMS最低但长周期不确定性增大。SNR=5 dB时各方法波动增强、短周期受瞬态噪声影响更显著，但中长周期结构仍被保留，RMS在三者间接近。权重极坐标图显示高权值点聚集于稳定阻抗幅相区，说明三类估计均能识别异常谱窗。偏差—不确定性合成总误差（Total Error=偏差²+不确定性）表明MM在最长周期总误差偏高，印证偏差—不确定性权衡。

3.2 东北野外应用：站点520中长周期M估计量局部起伏大，S较平滑，MM长周期误差棒增大；站点700四法整体吻合。远程参考+M估计量在多场景中仍具优势。

讨论与结论翻译

讨论指出：自然电磁信号非平稳且偏离正态，M最小化残差目标、S最小化尺度目标致不同信噪比下表现分化；仅靠因变量残差权无法完全消除磁场（预测变量）杠杆污染，故需远程参考或有界影响估计量辅助。Tukey双平方权中调节常数依赖经验，需自适应。合成噪声未完全覆盖人文噪声持续性，结果偏乐观。

结论部分翻译：研究人员提出一种基于迭代MM估计量的MT频域线性回归算法，并与M、S估计量在合成与野外数据中系统评估。结果表明中等噪声下三者均可给出可靠阻抗张量，但不同SNR下稳健性与崩溃点之权衡需深究；单电场或单磁场污染实验显示杠杆权重不能完全抑制预测变量离群，故建议三种稳健法任选其一与远程参考技术联用，以同时压制本地电场与磁场噪声。