4. 阻抗张量估算与测深曲线生成

大地电磁测深的数据处理链条中，从时间序列到最终可解释的视电阻率—相位曲线，中间最关键的一步是阻抗张量的估算。这一步决定了后续反演所依赖的“观测数据”是否可靠。本节围绕谱估计、Robust 阻抗估算、张量旋转与坐标系统、以及视电阻率和相位曲线的计算流程展开，并给出可操作的实现思路。

一、从时间序列到互功率谱

MT 的天然场源信号微弱，且常被工频、人文噪声、仪器漂移污染。阻抗估算的第一步是把两个水平电场分量 (E_x, E_y) 与两个水平磁场分量 (H_x, H_y) 的时间序列转换到频率域，并计算互功率谱。常用做法是分段（如每段 4096 或 8192 个采样点）、加窗（Hanning 或 Prolate 窗）、做 FFT，再对多段做平均得到谱矩阵。

对每个频点，谱矩阵包含 E 的自功率、H 的自功率以及 E–H 互功率。设 S_ij 表示分量 i 与 j 的互功率谱，则阻抗满足线性关系：

E_x = Z_xx·H_x + Z_xy·H_y
E_y = Z_yx·H_x + Z_yy·H_y

用最小二乘可解出四个阻抗元素。但普通最小二乘假设噪声只存在于 E 道，而实际中 H 道同样受噪声影响，且噪声常不服从高斯分布。这正是需要 Robust 估算的原因。

二、Robust 阻抗估算

Robust 方法的核心思想是：不追求让所有数据点都参与拟合，而是通过加权迭代降低离群点（强噪声段）的影响。工程上常用的是基于 Huber 权或 Tukey 双权（biweight）的 M 估计，配合“远参考”技术效果更佳。

一个可操作的迭代流程如下（示例参数，实际需按数据调整）：

1. 用普通最小二乘得到初始阻抗 Z⁽⁰⁾；
2. 计算每个频点、每段的预测残差 r = E − Z·H；
3. 按残差的稳健尺度（如 MAD，中位数绝对偏差）归一化，得到标准化残差 u；
4. 用 Tukey 双权 w(u) = (1 − u²)²（|u|<1 时）计算权重，|u|≥1 时权重置 0；
5. 以权重重新加权求解 Z^(k+1)，重复 2–4 步直至阻抗变化小于阈值（如 1%）。

需要注意：Robust 只能压制“少数”离群段。若某频段整体被强噪声占据（如 50 Hz 工频及其谐波、近场人文干扰），权重会大量置零导致该频点自由度不足，此时应结合远参考道或直接剔除坏频点，而不是盲目相信 Robust 结果。

三、张量旋转与坐标系统

阻抗张量是二阶张量，其元素依赖于测量坐标系的取向。设测量坐标系相对地质走向旋转了角度 θ，则旋转后的阻抗为：

Z'(θ) = R(θ) · Z · Rᵀ(θ)
R(θ) = [cosθ  sinθ; −sinθ  cosθ]

这里约定 θ 为从测量 x 轴逆时针旋转到新 x 轴的角度，全文与图示保持一致。旋转的目的一是寻找使张量非对角元素最大、对角元素最小的“主轴”方向（二维近似下的构造走向），二是把不同测点统一到同一构造坐标系下便于对比。

常用判据是阻抗张量旋转不变量：二维偏离度（skew）和椭率。当 skew 接近 0 时，数据近似一维/二维，旋转到主轴后 Z_xx、Z_yy 应接近零。实际中可扫描 0°–180°，绘制各元素随 θ 的变化曲线，取 |Z_xy−Z_yx| 最大或对角元素最小的角度作为走向估计。

四、视电阻率与相位曲线的计算

得到主轴方向阻抗后，对 TE、TM 两种极化模式分别计算视电阻率和相位。以非对角元素为例：

ρ_xy = (1 / (ω·μ₀)) · |Z_xy|²
φ_xy = atan2( Im(Z_xy), Re(Z_xy) )   （单位：度）

其中 ω = 2πf，μ₀ = 4π×10⁻⁷ H/m。若用 Z_yx 则得到 ρ_yx、φ_yx。相位定义需与阻抗实部、虚部符号约定一致，通常取 0°–90° 为正常感应范围，出现负相位或超过 90° 往往提示噪声、静态位移或三维效应。

实操中还应同步输出误差棒：由互功率谱的相干性和多段自由度估计阻抗方差，再传播到 ρ 和 φ。相干性低的频点（如 <0.5）其曲线可信度差，应在成图时标注或降权。

五、实操要点与结果判读

把上述流程落到一次实际处理中，建议按以下顺序检查结果：

• 相干性：逐频点查看 E–H 相干，低于 0.5 的频点应标记，必要时用远参考改善；
• 相位范围：正常感应相位应在 0°–90°；出现负相位或 >90° 需排查噪声、静态位移或三维效应；
• 二维偏离度：skew 较大（如 >0.3）说明三维性明显，此时旋转到“主轴”的意义有限，应谨慎解释；
• 曲线连续性：ρ 曲线出现孤立跳点多为坏频点，应结合误差棒判断是否剔除；
• TE/TM 分离：两条曲线在中高频分离常反映浅部二维/三维不均匀性，低频趋于重合则指示深部一维背景。

需要强调：视电阻率和相位曲线是“带误差的观测估计”，不是真值。成图时应同时给出误差棒和相干性信息，反演时也应把误差作为权重输入，而不是把曲线当作无误差的确定数据。只有把谱估计、Robust 加权、坐标旋转和误差传播作为一个整体来理解，测深曲线才具备可靠的地质解释基础。