刘杰，段炜，王俊，刘成，戴国强

(1.中国电建集团 昆明勘测设计研究院有限公司，云南 昆明650000；2.云南省地震局，云南 昆明650000)

0 引言

瞬变电磁法是先利用接地导线或不接地线圈，向地下发射一次磁场，然后突然关断一次场以激励地下导电地质体产生涡流，并在一次场关断间隙观测该涡流产生的二次磁场随时间变化的响应，通过分析响应曲线提取地下地质体的电阻率、埋深等信息，从而解决相关地质问题[1-3]。传统大定源装置在城市、公路、铁路工程勘察中难以快速布线，快速勘探，应用局限性大，近年来国内很多专家将小线圈中心回线、小线圈重叠回线应用于工程勘探、隧道超前预报等[4-12]。然而非固定小回线形态变化对一次磁场影响很大[13]，此外，由于接收天线与发射天线之间的强耦合关系，造成接收二次场畸变严重，反演电阻率失真，反演深度难以把控[14-17]。席振铢等[18-19]提出的等值反磁通瞬变电磁法(OCTEM)采用反向对偶磁源中心回线装置，消除了收发线圈互感影响，提高了浅层探测分辨率，因而在岩溶探测、铁路选线、公路病害探测等工程勘察中应用日益广泛。

本文主要研究OCTEM在某公路隧道塌陷区的抢险探测应用。针对待探测隧道的探测深度浅、围岩背景电阻率相对低、地形地貌复杂等特点以及快速勘察的需求，选择采用OCTEM进行探测。

1 工程概况及物探方法选择

1.1 工程概况

在建公路隧道位于云南山区，隧道进口端覆盖层以第四系残积的硬塑状粉质黏土为主，电阻率约100 Ω·m，平均厚约10 m；基岩主要为砂岩夹黏土岩，局部夹有褐煤层，属半成岩，岩质极软，围岩电阻率200～300 Ω·m，相对低。隧道长约300 m，隧道轴线距离地面高程20～70 m，埋深较浅。在隧道进口左线开挖到里程21 554 m附近时地面发生塌陷，隧道进口左线内出现涌水、涌泥现象。为避免再次出现塌陷和隧道涌水、涌泥，需要对隧道轴线上可能的塌陷区进行探测，为后期隧道的掘进施工提供依据。

1.2 物探方法选择

针对探测深度在几十米至百米深的工程勘探，常用的物探方法有地震波法、高密度电阻率法、瞬变电磁法等。由于该隧道地面地形差异大，而地震波法在复杂地形测区施工相对困难，数据解释比较复杂，因而未采用。高密度电阻率法在接地条件良好的地区适应性强，适用于本次勘察，然其不足在于探测深度与横向测点密度(横向分辨率)有时难以兼顾；为了提高高密度电阻率法的探测效果，考虑采用瞬变电磁法作为补充勘察。传统瞬变电磁法在山区布线困难，且由于收发电磁互感影响严重，早期信号受一次场影响严重失真，导致浅层探测盲区较大，而等值反磁通瞬变电磁法(OCTEM)克服了收发线圈互感影响，可大大减小浅层探测盲区，采用收发一体固定小回线装置不仅便携方便施工，且保证了各测点激发场源的一致性以及收发相对位置的一致性，可避免地形差异以及回线形态变化对一次场的影响，进一步减小观测噪声，提高浅层分辨率，因此OCTEM适用于本次勘察。

2 OCTEM方法原理

OCTEM方法原理与传统回线源瞬变电磁法相同，即先由人工源向发射回线供入电流产生一次场，然后瞬间关断电流，激励大地感应出二次涡流场，再通过感应线圈接收随时间衰减变化的二次场，根据其幅值和变化规律，获得地下电阻率分布信息。

OCTEM与传统瞬变电磁法不同之处在于其采用反向对偶磁源小型中心回线装置(图1)，上下两个发射线圈相同且平行共轴，二者电流同步且等值反向，接收线圈在距上下线圈相等的中间平面，在该平面上下两磁性源产生的垂向磁场大小相等方向相反，矢量叠加后垂向磁场恒为0，磁力线呈水平状，在该平面一次场总磁通始终为0，根据法拉第电磁感应定律，在发射电流关断前后接收线圈中将不产生一次场感应电动势，等效于消除了一次场的影响。在接收面以外的大地区域，关断电流前后一次场磁通发生变化激发感应涡流，被感应线圈接收。可见，OCTEM在早期不受一次场影响，将接收到大地纯二次场响应，可减小浅层探测盲区。

图1 OCTEM线圈及其径切面一次场磁力线Fig.1 OCTEM coils and its primary magnetic field line on the diametric plane

OCTEM测深是基于瞬态电磁场扩散原理。在非磁性(磁导率μ0)导电均匀大地中，当地表上方发射线圈中的电流突然关断时，首先在距离发射线圈最近的地表附近产生感应涡流以维持在断开电流以前存在的磁场。接着由于大地电阻引起欧姆损耗，感应涡流将随时间衰减变化，随着延时t增加，感应涡流极大值向下向外形如“烟圈”扩散。扩散速度v和扩散深度δ都与电阻率的平方根成正比：

而在地面接收到的二次磁场时间导数响应曲线幅值与大地电阻率平方成反比[20]，当大地电阻率减小至原来的1/10，其响应增加至原来的100倍，对低阻体反映较灵敏。综上，结合扩散深度和接收到的二次场衰减曲线幅值及其衰减规律可以分析得到测点下方电阻率的变化规律，达到测深目的。

3 理论计算

先通过正演计算[21]，为方法可行性提供理论依据。根据隧道实际情况建立简化的三维地电模型(图2)：围岩背景设为两层，第一层电阻率为100 Ω·m，厚10 m，代表浅表覆盖层；第二层电阻率为200 Ω·m，代表基岩层；塌陷地质体电阻率为50 Ω·m，设为规则长方体，其走向长度1 500 m，远大于其宽度W、厚度H及顶部埋深D，塌陷体中心在地表投影点对应X0位置，垂直塌陷体走向(沿着X方向)布置测线，点距25 m，剖面范围-200～200 m;发射电流1 A,发射有效面积1 m2，接收有效面积1 m2，激发场源为单位阶跃脉冲。

图2 塌陷体地电模型Fig.2 Geoelectric model of underground collapse

本文计算了大小不同、埋深不同的4组低阻塌陷体模型的OCTEM响应，塌陷体模型的宽度×厚度(即边长)分别为：25 m×25 m和50 m×50 m，顶部埋深分别为50 m和25 m，计算结果见图3，上图为磁场时间导数响应剖面曲线，下图为相对异常剖面曲线。

图3 不同厚度、不同埋深塌陷体OCTEM响应Fig.3 OCTEM responses of different size target collapse under different depths

从响应剖面曲线来看，除了边长25 m、顶深50 m塌陷体模型响应曲线(图3a)看不出明显隆起异常，其他3个模型响应曲线隆起异常明显,塌陷体越大，埋深越浅，异常越明显,若取相对异常值0.2为临界可分辨相对异常，则理论上4个塌陷体异常都可以被分辨。边长25 m的塌陷体，顶深25 m时的相对异常极大值为87.5%(图3a)，顶深50 m时相对异常极大值30.05%(图3b)；边长50 m的塌陷体，顶深25 m时的相对异常极大值为341%(图3c)，顶深50 m时相对异常极大值116%(图3d)；粗略测算，塌陷体规模增加至4倍，相对异常极大值也增加至约4倍;塌陷体埋深增至2倍，相对异常极大值减至约1/3。据此判断，在本文所述地电条件下，若以20%相对异常为可分辨阈值，则理论上某塌陷体的最大可探测深度大概为其横向尺寸的2倍，当围岩电阻率增加或塌陷体电阻率减小时，相对异常幅度也将增加，探测能力将提升。

4 应用分析

4.1 仪器设备

本次探测采用湖南五维地质科技有限公司研制的HPTEM18系统。设备由收发一体天线、主机、天线与主机连接电缆及移动PC组成，主机与移动PC通过WIFI连接。该系统运用等值反磁通原理消除收发线圈之间的感应耦合；利用对偶中心耦合原理提高横向分辨率；采用统一标准的微线圈对偶磁源、高灵敏度磁感应传感器、高速24位数据采集卡以及高密度测量技术实现浅层高精度瞬变电磁勘探。

4.2 数据采集与处理

本次数据采集采用定点测量模式，测点距2 m，发射频率6.25 Hz，发射电流10 A，关断时间50 μs，叠加周期400次。

采用HPTEM系统自带处理软件进行数据处理，主要包含信号校正、平滑滤波和电阻率成像三个步骤：①信号校正：收发一体小回线瞬变电磁系统具有施工灵活方便的特征，但是其振荡耦合噪声增强，包含了信号线缆耦合、接收与发射之间的残余感应耦合等噪声，需要对实测信号进行校正。②平滑滤波：瞬变电磁法在实际观测数据时，由于某一周期性的干扰或气候等其他因素的影响，无法通过多次叠加统计的方法将其剔除，导致衰减曲线在晚期出现跳变，出现不正常衰减，不利于数据的进一步处理和解释，因此需要对数据进行平滑滤波。③电阻率成像：类似于传统瞬变电磁法，HPTEM系统针对等值反磁通瞬变电磁法响应曲线特点提供了电阻率快速成像方法。该方法基于瞬变电磁“烟圈”扩散理论，综合考虑了实测电压幅值及其随时间衰减速率特征，能很好地反映局部电阻率变化特征，适合本次勘察目的。

4.3 成果图分析

图4为在该公路隧道轴线上里程21 502～21 810 m段的OCTEM探测结果，剖面长308 m。实测感应电压剖面(图4a)呈现多处大大小小的局部隆起异常，对应在电阻率断面(图4b)中出现电阻率等值线下凹或局部低阻异常区，推测为岩体塌陷或节理裂隙集中发育区。将OCTEM探测结果与高密度电法(电极距5 m，单个排列长295 m)探测结果(图5)进行对比分析：OCTEM异常S-1与高密度电法异常G-1对应，反映了已知塌陷，证明了两者方法的有效性；另外OCTEM异常S-2与高密度电法异常G-2接近；此外高密度电法异常G-3在OCTEM电阻率断面上呈现独立的3个异常S-4、S-5、S-6，OCTEM异常S-3在高密度电法图上没有明显对应，原因一方面是OCTEM反映电阻率变化灵敏度更高，另一方面因为其点距更小，横向分辨率更高。所以OCTEM能够弥补高密度电法分辨率与探测深度难以兼顾的问题。

图4 OCTEM感应电压剖面(a)及反演电阻率断面解释成果(b)Fig.4 Induced voltage profile (a) and resistivity section (b) of OCTEM

图5 高密度电阻率法反演电阻率断面及解释Fig.5 Resistivity section of high density electrical resistivity

5 结论

OCTEM接收发射固定一体，设备轻便，野外施工便捷，作业效率高，其采用反向对偶磁源中心耦合装置从物理上直接消除了一次场影响，大大减小了浅层探测盲区。另外正演计算表明：均匀导电大地中OCTEM对地电阻率变化反应灵敏，针对该隧道建立的良导围岩下的局部低阻塌陷体模型，当其埋深不大于其横向尺寸2倍时OCTEM相对异常大于20%。在本次公路隧道塌陷区的探测应用中，OCTEM对已知塌陷反应明显，与高密度电法探测结果异常吻合，证明了方法的有效性。

致谢：向中南大学瞬变电磁法研究团队成员提供的帮助表示感谢。