2. 工频干扰对电磁物探数据质量的影响评估
工频干扰之所以必须严肃对待，不是因为它“看起来脏”，而是因为它会沿着一条清晰的链条，把电压或磁场测量误差逐级放大为信噪比下降、视电阻率畸变，最终污染反演模型。评估的核心任务，就是把这条链条上的每一环量化出来，并据此判断数据处于哪一级干扰状态、还能不能用于后续处理与解释。
一、从信噪比开始量化危害
电磁物探的有效信号通常是宽频带弱信号，天然电场（如AMT、CSAMT）在几百赫兹以下的有效信号幅度常在微伏到毫伏量级，而50 Hz工频及其谐波通过大地、电缆、电极和仪器接地回路耦合进来，幅度可能达到有效信号的数倍。信噪比的量化最直接的方式是频域功率比：对一段时窗数据做FFT，分别统计50 Hz整数倍频点附近的窄带功率与相邻背景频带的功率，得到工频频点的信噪比恶化量。
设某测点原始数据在1～1000 Hz频段内背景噪声功率谱密度约为 P_bg，50 Hz处谱峰功率为 P_50，则工频干扰对总噪声功率的贡献可近似为 ΔP = P_50 − P_bg。若有效信号在该频点的功率为 P_s，则干扰前后的信噪比变化为：
SNR_before = 10·log10(P_s / P_bg) （dB）
SNR_after = 10·log10(P_s / (P_bg + ΔP)) （dB）
SNR_loss = SNR_before − SNR_after （dB）举例说明（示例参数，非实测）：某AMT测点有效信号在50 Hz附近功率为1×10⁻¹² V²，背景噪声为2×10⁻¹³ V²，则SNR_before≈7.0 dB；若50 Hz谱峰功率达4×10⁻¹² V²，ΔP=3.8×10⁻¹² V²，SNR_after≈−5.8 dB，损失约12.8 dB。这意味着原本勉强可用的频点直接跌入不可用区间。工程上常把单频点SNR损失超过6 dB视为“显著污染”，超过12 dB视为“严重污染”。
二、视电阻率计算中的畸变机制
视电阻率由阻抗张量或电场磁场比值计算，公式为 ρ_a = (1/ωμ)|E/H|²（或对应标量形式）。工频干扰对E和H通道的污染并不对称：电场通道通过电极接地电阻和电缆天线效应耦合较强，磁场通道相对较弱，因此E/H比值被系统性抬高或压低，导致视电阻率曲线在50 Hz及其谐波处出现尖峰或凹陷。更麻烦的是，若干扰幅度随时间漂移（如负荷变化引起的工频幅值波动），相邻频点的相位也会被拉扯，使相位曲线在干扰频段附近出现非物理的跳变。
量化这一畸变，可用“相对视电阻率偏差”指标：以干扰频段两侧的干净频点做平滑插值，得到该频点的期望视电阻率 ρ_exp，实测值为 ρ_obs，则偏差率 δ = |ρ_obs − ρ_exp| / ρ_exp × 100%。δ在5%以内通常可接受，5%～20%需标记，超过20%则该频点应剔除或降权。实际处理中，还要检查同一测点相邻频点的连续性——单点尖峰往往是工频干扰，而整段抬升则可能是静态位移或系统性问题，二者不能混为一谈。
三、反演结果层面的危害
反演对噪声的响应不是线性的。少量强干扰频点若未剔除，会在反演中形成“虚假高导层”或“虚假高阻体”：因为反演算法倾向于用模型结构去拟合这些异常频点，把干扰误认为地下真实电性变化。更隐蔽的是，若干扰频点被错误地当作有效数据参与加权，反演的目标函数会被拉偏，导致整体模型的分辨率下降、深部界面模糊。评估反演层面的危害，可对比“含干扰数据反演模型”与“剔除干扰后反演模型”的差异：统计两者在目标深度上的电阻率相对偏差，以及数据拟合残差中50 Hz频点的贡献占比。若残差中工频频点贡献超过总残差的30%，说明干扰已实质性主导了反演过程。
四、建立干扰严重程度分级意识
把上述指标整合起来，可以形成一套可操作的分级判断。下图给出从原始数据到分级结论的评估链路，以及各级对应的处理策略。
分级的意义在于把“感觉干扰大”变成“可比较、可决策”的判断。Ⅰ级数据经常规陷波与平滑即可使用；Ⅱ级必须做频点剔除或降权，并在反演报告中记录处理方式；Ⅲ级数据原则上不应进入反演，应优先考虑复测、改善接地与布线，或改用抗干扰能力更强的采集参数。需要强调的是，三项指标应取最严重者定级，因为反演对局部强干扰的敏感度远高于对整体噪声水平的敏感度。
五、评估时的实操要点
第一，评估窗口要覆盖完整观测时段，不能只看一段“看起来干净”的数据，工频干扰常随用电负荷呈日变化。第二，E、H通道要分别评估，不能只算合成阻抗，否则会掩盖单通道的严重污染。第三，谐波不可忽略，150 Hz、250 Hz等奇次谐波在AMT高频段常与有效信号重叠，评估时应把50 Hz整数倍频点一并纳入统计。第四，保留评估记录：测点号、时段、三项指标值、定级结论、处理动作，形成可追溯的质量档案。只有把量化指标、分级结论和处理策略绑定在一起，干扰评估才真正服务于数据质量控制和反演可靠性。