1. 从“看电阻率”到“看时间”：瞬变电磁法的底层逻辑
做物探的人都知道，直流电法看的是稳定电流场，频率域电磁法看的是谐变场，而瞬变电磁测深（TEM）完全换了条思路——它不靠稳态，靠的是“断电后的那一瞬间”。这个区别，正是第三章整章的立论基础。
先说清楚它到底在干嘛。瞬变电磁测深的工作原理可以浓缩成一句话：在地面布置一个方形或圆形回线，通入恒定电流，让地下建立起一次磁场；然后突然切断电流，利用接收线圈记录地下感应涡流产生的二次场随时间衰减的响应。这个二次场的衰减曲线，携带了地下介质的电阻率分布信息。时间越靠后，对应探测深度越深，这就是“测深”二字的物理来源。
有个很重要也很容易混淆的概念：瞬变电磁法测的到底是什么？很多人以为它测的是电阻率本身，其实不对。仪器记录的是二次场电压随时间的变化，也就是衰减曲线，电阻率是通过反演计算反推出来的。早期早期仪器甚至直接用晚期段电压值做视电阻率估算，但实际工作中，衰减曲线的形态、早晚段的划分、关断时间的影响，任何一个环节处理不当，结果都可能偏离真实地电断面。
这章内容的定位，我理解下来是承上启下的。前面两章大概率讲的是电磁场基本理论、Maxwell方程组和扩散方程，从第三章开始才正式进入TEM的核心——瞬变过程怎么建立、怎么传播、怎么被观测，进而引出全期视电阻率、正演计算、反演约束这些实际处理问题。
在项目实践中，无论你用的是哪家国产或进口仪器，处理软件怎么换，底层物理模型都是这一章讲的东西。搞不清“关断时间”“一次场”“二次场”“扩散深度”这些基础概念，后面做一维反演和二维解释的时候会非常吃亏。
所以这篇总结，我打算把这章里几个最难啃、也最容易让新手翻车的知识点拆开揉碎讲清楚：瞬变响应是怎么产生的、扩散深度和探测深度的区别、全期视电阻率为什么比早期/晚期公式靠谱、以及正演计算中最容易被忽视的网格与时间道设置问题。
2. 从断电瞬间到涡流扩散：瞬变响应的物理过程拆解
2.1 一次场建立：发射回线通流后的稳态
在做瞬变测深之前，发射回线需要先通入恒定电流。这个过程通常持续几十毫秒到几百毫秒，视仪器型号和供电能力而定。电流稳定后，地下空间内会形成一个稳定的磁场分布，这就是一次场。
很多人会忽略一个细节：发射回线通流后，地下的磁场不是瞬间建立起来的。由于大地介质具有一定的导电性，磁场建立过程伴随着电磁能量的扩散，只是这个过程相对于观测窗口来说非常短暂，大多数情况下可以近似当作瞬态过程来忽略。但在某些高阻覆盖层发育的地区，一次场稳定时间可能会略微拉长，实际处理时如果发现早期道数据异常跳动，排除仪器因素后，可以考虑这个原因。
一次场的分布形态，在均匀半空间条件下可以用简单的偶极子场叠加来理解——回线每个边段都相当于一个磁偶极子，整个回线产生的总场是各段贡献的叠加。但在层状介质中，这个场会被界面反射和折射，变得复杂得多。这也就是为什么正演计算不能简单套均匀半空间公式的原因。
不过，第三章的重点并不是一次场怎么建立，而是断电后发生的事。因为TEM真正有用的信号，恰恰来自电流切断之后。
2.2 关断瞬间与感应电动势的产生
电源断开后，发射回线中的电流并非瞬间归零，而是经过一个短暂的关断过程。这个过程的时长叫“关断时间”（ramp time），在目前的仪器上通常为几微秒到几十微秒。关断时间越短，高频成分越丰富，对浅部分辨能力越强；关断时间越长，早期信号损失越严重。
也正是在这个关断过程中，地下介质中感生出涡旋电流。物理图像可以这样理解：一次场的变化，在地下导体中产生感应电动势，这个电动势驱动了二次涡流；而二次涡流又反过来产生二次磁场，被接收线圈记录下来。
要注意的是，感应电动势的大小与一次场的变化率成正比。所以在关断时间很短的情况下，感应电动势峰值非常大，如果接收线圈和发射回线之间的互感补偿做得不好，很容易造成接收放大电路饱和，早期道数据直接畸变。这种畸变不是简单滤波能消除的，需要在野外采集时通过“发射电流波形记录 + 同步对时”的方式处理，后期做关断时间校正才能修回来。
这里有一个非常关键的实践原则：野外必须记录发射电流的关断波形，不能只记录关断时间这一个参数。因为实际关断波形并非理想的线性斜坡，而是带有一定的弯曲，尤其在大电流、长导线条件下。数值模拟表明，用理想斜坡代替实际波形做校正，在早期道可引入百分之几到十几的误差。
2.3 地下涡流环的扩散：烟圈效应
断电以后，地下涡旋电流不是原地衰减的，而是会随着时间推移逐渐向深处、向外侧扩散。这个现象有个经典的类比——“烟圈效应”。你可以想象你往地下一口气吐出一个电流环，这个环会逐渐下沉、扩张、强度减弱，就像一枚烟圈在水中慢慢散开。随着扩散深度增加，所反映的介质体积越来越大，对应的电阻率信息也越来越趋向于深部平均效应。
从数学上说，涡流场的扩散遵循扩散方程，扩散深度与时间是平方根关系。均匀半空间条件下，扩散深度可写成海姆霍兹形式的近似表达式：
d ≈ √(2t/（μσ）)
其中t是延迟时间，μ是磁导率，σ是电导率。这个式子说明，同样一个时间道，导电性越好的介质中涡流扩散越慢，探测深度越小；而导电性差的介质中扩散快，同样时间能看更深。这一点对于解释TEM测深“在高阻区探测深度大、低阻区探测深度小”这一现象非常关键。
但也正是这种扩散特性，导致TEM的纵向分辨能力有限。涡流环的半径在随时间增长，携带的信息是一个越来越大的体积平均响应。很多新手拿到反演断面图的时候，总觉得电阻率分层平滑得像奶油蛋糕一样，这是正常的——TEM本身就是体积探测，不是点测。想要提高纵向分辨率，只能在观测系统中做文章，比如用更小的发射回线、更短的关断时间、更密的测点距，以及后续的中心回线/重叠回线装置选择。
2.4 接收信号：一次场与二次场的分离
接收线圈记录到的总磁场变化率，实际上既包含一次场的变化，也包含二次场的变化。在关断结束后的一次场残留期间，接收信号以一次场为主；随着一次场迅速衰减为零，后续记录段就是纯粹的二次场衰减。
问题在于，实际观测中很难做到将一次场完全分离干净。尤其在采用“中心回线装置”时，发射回线与接收线圈同心布置，关断瞬间一次场变化剧烈，接收线圈中感应出的电动势极大，如果接收系统动态范围不够，甚至会造成振荡，影响后续相当长一段时间的二次场记录。这就是“关断盲区”概念的实际来源——在关断后的一小段时间内，信号不可靠。
解决思路大致有三条：
- 硬件上尽量缩短关断时间，减小盲区宽度；
- 软件上做关断时间校正，从记录的电流波形反卷积还原纯二次场；
- 观测方式上采用“分离回线装置”，发射回线与接收线圈保持一定距离，从而降低一次场耦合。
三种方法各有利弊。硬件方案受限于发射功率和回路电感；软件校正对电流波形记录精度要求较高；分离回线装置则牺牲了浅部分辨率和信号强度，更适合探深型任务。实际项目里，我见得最多的还是中心回线+软件校正的组合，兼顾浅层分辨率和施工效率。
3. 时间域与频率域的对照：为什么瞬变测深更“省事”
要真正理解TEM，其实不能只盯着时间域，最好把它和频率域电磁法放在一起对照看。两者描述的是同一个物理过程，只是观察角度不同。
频率域电磁法（如可控源音频大地电磁法CSAMT、大地电磁法MT）用一系列不同频率的谐变场来分别探测不同深度——频率越高，集肤深度越小，看越浅；频率越低，看得越深。所以频率域仪器本质上是在扫频，每个频点对应一个深度范围。
瞬变电磁法则是一次发射包含所有频率成分的脉冲信号，然后在时间域里按不同延迟时间窗口分别采样。延迟时间越短，对应有效频率越高，反映浅部；延迟时间越长，有效频率越低，反映深部。也就是说，一次发射相当于同时完成多次不同频率的测深。
这样一来，TEM的最大优势就体现出来了：效率高。同样测一个剖面，频率域方法需要逐个频点观测，TEM只需要在每个测点发射几次、记录一条衰减曲线，后续通过时窗积分就能得到多深度信息。在野外施工效率和成本上，TEM的优势不是一星半点。
但也不能光说好处。频率域方法的优势在于它可以用相干积分、重复堆叠等手段非常有效地压制噪声，而TEM的时间域信号，特别是在晚期道，信号幅度已经非常小，往往被噪声淹没。这就是为什么TEM晚期数据质量普遍较差，也是为什么反演时通常要给晚期数据施加更大的权重衰减或直接截断。
另一个有趣的关联是换算关系。通过傅里叶变换，时间域响应和频率域响应可以相互转换。在实践中，某些正演算法是在频率域计算再变换回时间域，因为频率域解析表达式在层状介质中推导更成熟；也有一些算法直接在时间域做递推计算。两种路径得到的数值结果应该一致，这也是验证正演代码是否正确的一个有效手段。我早期做代码验证的时候，就喜欢用频率域结果做交叉校验，往往能发现一些不显眼的积分步长或边界条件错误。
4. 全期视电阻率：从“晚期近似”到“全程可用”
4.1 早期和晚期公式的适用条件
正如前面所说，仪器直接记录的是电压衰减曲线。要从这个曲线得到电阻率信息，最朴素的方法是套用均匀半空间的解析公式。
在时间域的早期段，电磁场还没有明显扩散到深层，可以假设半空间温度计，偏导简化后得到一个只与早期电压和时间相关的表达式，这是“早期视电阻率公式”。在晚期段，涡流扩散已经足够深，可以视为准静态极限，于是又得到“晚期视电阻率公式”。
问题在于，这两个公式的适用窗口在实际数据中常常衔接不上。中期的数据段，既不符合早期近似条件，也不符合晚期近似条件，用任何一个公式算出来的视电阻率都会失真。尤其在地下存在明显的低阻覆盖层时，早期段持续时间非常短，很快进入“过渡区”，这时无论套哪个公式都别扭。
我印象很深的一次经历，是在某矿区做数据整理时，用商业软件默认的晚期公式算全曲线视电阻率，浅层道算出来电阻率直接出现负值——这显然是公式失效的典型特征。当时我对参数做了半天检查，最后才意识到是公式选用的问题，换成全期视电阻率后，曲线才恢复正常形态。
4.2 全期视电阻率的求取思路
全期视电阻率的本质是：对于观测到的每一个时间道的电压值，反求一个均匀半空间模型的电阻率，使得该模型的理论响应恰好在相同时刻等于观测值。换句话说，不再区分早晚，全程都用一个“等效均匀半空间”来匹配观测响应。
数学上这是一个单变量求根问题。给定电阻率ρ，计算模型在该时刻的理论电压值V理论(ρ,t)，与实测V实测(t)比较，通过迭代调整ρ，使两者相等。常用的求解方法有二分法和牛顿迭代法。二分法实现简单、稳定性好，缺点是迭代次数略多；牛顿法收敛快，但要求初始值离真解不太远，否则容易震荡。
这里有个实际技巧：在用二分法求全期视电阻率时，初始搜索区间不能拍脑袋设成“0.001到100000Ω·m”这样的大区间。虽然二分法对区间没太多要求，但区间过大会导致前期无效迭代太多，浪费时间。一般可以先用晚期公式求一个粗估的电阻率值，然后把这个值上下各扩大两个数量级作为搜索区间，这样通常只需要十几二十次迭代就能收敛到位。
4.3 为什么全期公式在层状介质中仍有局限
全期视电阻率克服了早期/晚期近似局限于均匀半空间的缺陷，但它毕竟建立在一个“均匀半空间”模型假设上。当地下存在多层结构时，这个等效电阻率不是一个真实层的电阻率，而是一个复杂的加权平均值。
换句话说，全期视电阻率曲线仍然不是地下真实电阻率剖面，只是一个中间量。它的价值在于：作为反演的初始模型输入、用于快速判断地电断面的大致特征、用于质量控制（与理论曲线对比看拟合情况）。换句话说，它是“给人看的”，不是“给最终结果用的”。
这也是我反复跟同行强调的一个观点：视电阻率曲线再平滑、再漂亮，也不能直接拿来做地质解释。真正的分层和电阻率定量解释，必须走反演这一步。视电阻率的作用是指导反演初始模型和判断数据质量，不是替代反演。
5. 正演模拟的实操要点与常见踩坑
5.1 层状介质一维正演的算法选择
第三章必然会涉及一维正演计算。层状介质中，瞬变电磁响应的计算通常有两种路径：频率域汉克尔变换法和时间域直接递推法。前者通过将层状介质的反射系数表达式代入，计算频率域响应，再通过余弦/正弦变换转到时间域；后者则是直接在时间域内用传输矩阵表示各层之间的场耦合。
两种方法各有特点。频率域方法在计算速度和稳定性上总体更优，因为频率域反射系数的表达式是解析的，积分变换也有成熟的数值算法（如汉克尔变换的滤波器系数法）可用；时间域方法则更直观，容易做边界条件限定和复杂模型的扩展，但数值稳定性对时间步长和层厚度比较敏感。
在实操中，我一般优先建议频率域+数字滤波器法。原因很简单：实现路径成熟，参考资料多，且易于做精度对比验证。滤波器系数方面，业内用得比较多的是一组127点或201点的正弦/余弦变换滤波器系数，网上有公开数据集可下载。需要说明的是，滤波器系数的选择影响最终结果的精度，不同作者发布的系数之间存在细微差异，在深层晚期响应计算时，这种差异可能显现为几个百分点的偏差，必要时需要做两组系数的交叉验证。
5.2 发射回线参数的设置与离散化
正演计算中有个经常被忽略但影响很大的参数：发射回线的边长。很多初学者为了省事，直接把发射回线当作磁偶极子来近似处理，在小回线（边长小于20米）情况下勉强可用；但如果回线边长达到100米以上，这种近似就会带来明显偏差，尤其在浅层早期道响应中，误差可能超过10%。
正确处理方式是把回线离散成若干个小段，对每个小段的场贡献做数值积分。离散段数越多，精度越高，但计算量也越大。经验上，回线边长在100米以内时，每边离散成10段左右就能达到千分之几的精度；边长增大到500米以上，需要加密到每边20段以上。
另一个容易出错的点是发射电流波形的模拟。正演时如果简单假设阶跃关断（即电流瞬间变为零），会给早期道带来系统性偏差。实际观测中关断时间不为零，正演模型应该把梯形波或实测电流波形纳入计算，然后通过“同步观测”的方式与实测数据匹配。目前主流商业软件基本都支持自定义发射波形，使用方一定要确认自己导入的波形文件是否与实际采集参数一致，否则反演出的浅部电阻率会整体系统性偏离。
5.3 时间道的选择对计算精度的影响
时间道的对数等间隔分布是TEM正演和反演中最常见的设置。通常覆盖从几微秒到几十毫秒的范围，分成20到40个时间道。
有一个长期困惑新手的问题是：时间道设置得越密越好吗？并非如此。时间道过密会增加计算量，而又不增加独立信息量，因为相邻时间道的响应高度相关；时间道过疏则可能漏掉关键的衰减形态特征。实践中，每个数量级（以10为底）分布10个时间道左右通常就够了，一共设置30道左右已经是比较均衡的方案。
另外要注意时间道的覆盖范围应与仪器的实际观测窗口一致。正演算出的数据响应只在与观测窗口重合的区间内才有可比性，超出观测窗口的早期段或者晚期段，正演结果再精细也没有实际对应数据。做拟合时，必须只针对有效时间范围内的数据进行，不能拿全窗口正演曲线去跟截断后的实测数据对比。
5.4 正演结果常见异常及原因排查
我自己调试正演代码时积累了一份排错清单，这里分享出来，省得大家走弯路。
第一类异常：晚期响应出现震荡或负值。多半是数字滤波器系数选择不当，或者积分上限不够高、截断频率过低。解决办法是提高采样频率范围，或更换更长的滤波器系数。
第二类异常：早期响应出现跳变。这往往是时间步长太粗，或者发射回线离散段数不足。可以细划时间轴或增加回线分段数验证。有时也可能是因为发射波形未做同步延迟校正，导致脉冲位置偏移。
第三类异常：全期视电阻率曲线在浅层出现负值。这个问题在前面章节提过，多数是晚期公式误用导致的，但也可能是正演响应本身在早期段存在数值振荡。两者区分方法很简单：用均匀半空间模型做基准测试。均匀半空间下，全期视电阻率应该恒等于设定电阻率。如果这个最基本的标准都过不了，那代码里大概率有坐标符号或积分核函数写反的问题。
第四类异常：层状模型的响应与理论渐近线不一致。例如在极高阻基底上，晚期道响应会趋近于某个渐近值。此时可以解析估算晚期渐近值，与数值结果对比。如果偏差超过1%，说明深层响应的积分处理精度不足，需要对积分区间的末端做特殊处理（比如变步长加密）。
6. 野外数据质量控制的几个关键环节
6.1 发射回线与接收线圈的几何关系
野外采集时，我们最容易忽略的其实是几何位置误差。中心回线装置中，接收线圈偏离回线中心哪怕十几厘米，带来的浅层早期道偏差都不可小觑。尤其在复杂地形条件下，回线形状可能发生扭曲，不再是标准矩形，这时正演假设的几何形态就和实际情况产生了偏差。
解决办法：在布置回线的时候，尽量用GPS或全站仪把四个角点坐标精确测下来，回线边长误差控制在1%以内；使用通条固定回线，避免地形起伏造成回线悬空。如果地形条件确实不允许，至少要记录每个测点回线的实际形状，供后期正演校核。
6.2 关断时间记录与数据预处理
之前强调过的电流波形记录，在野外实施时有两种做法：一是仪器自动记录每次发射的电流波形回传；二是用采样率更高的示波器独立记录。前者操作简单，但采样率可能不足以及时捕捉关断沿的细节；后者精度高但对野外操作要求更高。
数据预处理阶段，最关键的步骤是同步校正。因为发射与接收存在一定的同步延迟，全部时间道都需要做整体平移。这个延迟如果校正不准，浅层早期道的电阻率会明显失真。常见做法是选择一个浅层已知地电模型（比如测区内的已知钻孔柱状图），用正演结果拟合实测早期道，反推时间零点位置。
6.3 噪声来源分析与叠加策略
TEM观测中的噪声来源主要有五种：一是仪器本身的电子噪声，二是接收线圈运动产生的摆动噪声，三是天电噪声（雷电等宽频干扰），四是人文电磁干扰（高压线、通信基站、工业游散电流），五是风致振动。
抑制手段对应如下：仪器噪声可以通过多次堆叠压制，堆叠次数通常取64到2048次；摆动噪声可以通过将接收线圈固定在三脚架或埋入浅坑减少；天电噪声往往集中在特定频段，可以用带通滤波滤除但要注意不伤及信号；人文干扰则主要通过选择干扰小的时段重复测量来规避；风致振动主要靠降低接收线圈重心和使用防震垫来缓解。
有一个容易被忽略的叠加技巧：叠加时不要简单求算术平均，建议采用中值叠加或剔除异常值后再叠加的策略。理由是偶然出现的尖峰脉冲干扰（比如附近有人使用对讲机）会在平均中保留残差，却能被中值过滤干净。实际使用中，20%的剔除率通常能带来明显的质量提升。
7. 从正演到反演：第三章内容的延伸思考
第三章的核心是理解瞬变响应本身，但它最终要服务于反演和解释。这里就涉及一个衔接问题：正演算出来的响应曲线，和实际观测数据拟合时，到底什么是“好的拟合”？
很多初学者拿着反演软件，习惯性地追求最小拟合差，把拟合差降低到零点几个百分点，然后认为反演结果一定正确。但TEM反演是个高度病态的问题，同样的观测数据可以由很多不同的地电断面拟合得到相似水平的响应。追求极低拟合差，往往逼迫反演算法生成锯齿状或极端值模型，反而失去了地质合理性。
更合理的思路是：在拟合差达到合理水平（通常3%到5%，视噪声水平而定）后，加入模型平滑约束或层厚度先验约束，使反演结果走向稳定且地质合理的解。反演的“多解性”不是靠提高拟合精度解决的，而是靠加入先验信息解决的。
第三章的正演计算能力，恰好为这种“约束反演”提供了基础——你能快速计算出某个修改后模型的理论响应，即可实时判断某个界面位置或电阻率调整方向是否合理。这也是我建议所有做TEM数据处理的人，哪怕不做反演开发，也要把正演计算吃透的原因。它就像是你的“理论标尺”，离开它，你根本无法辨别反演给出的模型是真实反映还是数值伪影。
在实际项目中，我通常的做法是：先用全期视电阻率曲线快速判断地电断面的大致层数和电阻率量级，然后据此设置反演的初始层数和电阻率范围，再让反演在带约束条件下迭代。这样既保证了收敛速度，又能避免反演掉进局部极小值的泥潭。
另外一个经验是，无论软件怎么自动反演，都要抽几条典型测点的数据，单独做手动正演拟合对比。这种“笨办法”虽然费时间，但能真正暴露反演结果中的不合理细节。比如某个看似完美的反演断面，在手动拟合时可能会发现某个测点的深层响应始终拟合不上，这就提示深层可能存在旁侧体或各向异性影响，而这是标准一维反演无法表达的信息。
8. 实操心得与进阶方向
写到这里，第三章的核心内容基本覆盖了。最后再分享几点我在长期项目实践中积累的操作体会，供同行参考。
第一，做TEM数据处理之前，务必先看一遍原始衰减曲线的形态。不要急着滤波、反演。如果曲线本身呈现明显的无规律跳变或台阶状畸变，那么任何后续处理都大概率是徒劳的。数据质量是解释可靠性的前提，这在TEM中比在直流电法中更为严苛。
第二，正演和反演计算的参数记录一定要完整。很多人处理数据时用了什么滤波器系数、什么时间道设置、什么网格剖分参数，做完就忘。等到数据到解释阶段需要回溯时，发现记录不全，只能重算——这浪费的时间足够让人崩溃。我的习惯是在每个测区的工作目录里建一个参数记录文本，每次处理前把软件参数导出存档，包括发射回线边长、电流大小、关断时间、时间道设置、滤波参数、反演约束条件等，全部写清楚。
第三，关于浅层和深层信息平衡的问题。TEM对浅层的分辨能力受关断时间和采样率限制，对深层则受信号幅度和噪声水平限制。如果想兼顾两者，在观测系统设计阶段就要做好权衡——大回线大电流有利于深层信号，但对浅层盲区加大；小回线小电流有利于浅层精细探测，但深部信号幅度不足。这是不可调和的矛盾，所以很多项目采用“大小回线组合”的方式，同一测区做两套观测，后期融合解释，效果比单套系统好得多。
第四，拓宽视野。第三章的内容虽然聚焦在层状介质一维响应上，但实际地质体没有严格的“层状”特征。当测区存在明显的构造破碎带、岩体接触带或地下采空区时，一维反演会出现明显偏差，这时需要考虑二维甚至三维正反演。建议在掌握了第三章的基础之后，尽早接触三维数值模拟软件，哪怕是简单模型的计算，也能帮你建立对多维效应的直观感知。
第五章讲到后期，我会继续围绕二维条件下的TEM响应特征展开。那部分的内容，跟实际工程和矿区应用的关联性更强，届时再详细展开聊。如果你正在做TEM相关的数据处理或项目设计，不妨先把第三章这些基础问题彻底弄通，后面再进阶时会顺畅得多。
