1. 课程导入与CSAMT方法概述
可控源音频大地电磁法(Controlled-Source Audio-frequency Magnetotellurics,简称CSAMT)是在天然场大地电磁法(MT/AMT)基础上发展起来的一种人工源频率域电磁测深方法。它通过接地电偶极子或水平磁偶极子向地下发送可控的音频频段电磁信号,在收发距适中的远区观测相互正交的电场与磁场分量,进而计算卡尼亚视电阻率和阻抗相位,反演地下介质的电阻率分布。与依赖天然场的方法相比,CSAMT具有信号强度可控、信噪比高、频带与测点布置灵活、生产效率较高等优点,因而在矿产资源勘查、地热资源调查、地下水与工程环境探测、深部构造研究等领域得到广泛应用。
一、CSAMT的基本原理
CSAMT的物理基础是电磁感应与趋肤效应。发射源供以交变电流后,在周围空间建立交变电磁场;电磁场向地下传播时,高频成分衰减快、穿透浅,低频成分衰减慢、穿透深。因此在同一测点上,通过改变发射频率并观测场值,就可以获得随频率变化的响应,从而推断不同深度的电性结构。频率越低,勘探深度越大,这就是频率测深的基本思想。
在测量方式上,通常沿测线布设接收点,测量两个水平正交的电场分量 Ex、Ey 和两个水平磁场分量 Hx、Hy(实际工作中常测一组正交的 E 与 H)。在源为电偶极子、测点位于远区(收发距 r 远大于趋肤深度)的条件下,场近似满足平面波假设,此时可由正交的电场与磁场分量计算波阻抗:
Zxy = Ex / Hy (卡尼亚阻抗,单位 Ω)
ρa = |Z|² / (ωμ₀) = |Ex/Hy|² / (2πfμ₀)
φ = arctan( Im(Z) / Re(Z) ) (阻抗相位,单位 mrad 或度)
其中 ω = 2πf 为角频率,μ₀ = 4π×10⁻⁷ H/m 为真空磁导率。由上式可见,视电阻率 ρa 是频率 f 与波阻抗模值的函数,它并不是某一深度的真实电阻率,而是该频率所对应“有效勘探体积”内介质电阻率的综合等效反映。这正是“视电阻率”一词的含义,也是后续曲线处理与解释工作的出发点。
二、视电阻率的定义与曲线特征
需要特别强调的是,CSAMT中计算得到的视电阻率受多种因素影响:一是地下介质的真实电阻率分布;二是收发距与场区条件,若测点不在远区,近场效应和过渡场效应会使视电阻率曲线发生畸变;三是场源方向与测线方向之间的几何关系,会产生所谓的场源效应与阴影效应;四是地形、静态位移、电磁干扰等非地质因素。因此,原始视电阻率曲线往往并不是一条光滑、单调的曲线,而可能出现飞点、跳变、高频段发散、低频段下掉等异常形态。
从曲线形态看,常见的有:均匀半空间对应的平直曲线;两层断面(如低阻覆盖层下伏高阻基岩)对应的“G型”或“D型”曲线;三层断面可能出现的“H型”“K型”“A型”“Q型”曲线。识别这些基本形态,是判断曲线是否合理、是否需要处理的第一道关口。下图给出了CSAMT从发射、传播、观测到视电阻率计算与异常解释的整体信号链与流程位置关系。
三、曲线处理与异常解释在勘探流程中的位置
从上面的信号链可以看出,CSAMT勘探并不是“测完即得结论”的线性过程,而是一个“采集—处理—反演—解释—验证”的闭环。曲线处理处在承上启下的关键位置:向上,它要消化野外采集质量、场源几何、干扰水平等信息;向下,它直接决定反演输入数据的可靠性。如果原始视电阻率曲线中的飞点、近场畸变、静态位移没有被识别和压制,反演就会把这些假异常“固化”成地下结构,导致解释结论偏离真实地质情况。
具体而言,曲线处理通常包括几个层次:第一是数据质量检查,逐点逐频查看视电阻率与相位曲线,识别明显偏离总体趋势的飞点,并结合原始时间序列、相干度、信噪比判断其成因;第二是畸变校正,包括近场与过渡场校正(当收发距不足时,远区公式失效,需用全场公式或校正曲线修正)、静态位移校正(由浅部小尺度电性不均匀引起的曲线整体平移,常用空间滤波或参考模型法处理);第三是曲线圆滑与频点取舍,在保留真实异常特征的前提下压制随机噪声。异常解释则是在反演断面基础上,结合地质、钻探、物探多源信息,对低阻/高阻异常的地质属性作出判断,并给出验证建议。
四、课程目标与学习路径
本课程围绕“视电阻率曲线处理”与“异常解释”两条主线展开,目标是使学习者能够独立完成从原始CSAMT数据到地质解释结论的主要环节。具体包括:理解CSAMT的场源、远区条件与视电阻率定义,能判断数据是否满足方法前提;掌握曲线质量检查与常见畸变识别方法,能针对飞点、近场效应、静态位移等采取合理处理措施;理解反演的基本流程与参数选择,能对反演结果进行合理性评价;掌握异常解释的基本逻辑与多源信息互证方法,能编写结构清晰的解释成果。