亓庆新，席振铢，徐昱，木仁

(1.湖南五维地质科技有限公司，湖南 长沙 410205; 2.中南大学 地球科学与信息物理学院，湖南 长沙 410083；3.湖南省计量检测研究院，湖南 长沙 410014)

0 引言

瞬变电磁法[1-3](TEM)，也称时间域电磁法，在地球物理勘探领域中是一种非常重要的人工源电磁法[4-5]。TEM受地形影响小、施工效率高，广泛应用于地热、煤炭、水资源、矿产勘查、工程勘查、地质灾害调查、地质环境调查、考古探测、军事探测、城市管线检测等诸多方面[6-8]，近几年已逐渐向大深度矿产[9]和油气资源勘探[10-11]方向发展。此方法以电磁感应原理为基础，通过不接地回线向地下发射一次脉冲磁场，在一次脉冲磁场的激励下，地下地质体将产生感应涡流，感应涡流产生的二次磁场不会随一次磁场消失而立即消失，在一次脉冲磁场间歇期间利用天线或接地电极观测二次磁场,通过研究感应二次场响应幅值和衰减速度变化规律，从而达到探测目标体空间展布的目的[12-13]。

在实际应用中，发射线圈与接收线圈作为相邻线圈，两者之间会发生互感现象，在关断一次场后，接收线圈受一次场变化影响会产生感应电动势，而且此一次场电动势通常比二次场在接收线圈中产生的电动势大得多，并难以区分、无法剔除，导致瞬变电磁法在浅层探测中出现盲区[14-16]，严重干扰早期接收信号精度。Smith 等[17]、嵇艳鞠等[18]、Walker等[19]通过对瞬变响应和发射电流波形的实时监测，再根据数值计算方法把接收信号中受一次场影响产生的感应信号剔除，获得纯二次场响应，但是实际勘探中受地质条件等原因影响，线圈大小和位置与理论值存在偏差，导致计算值与实际感应值不符；范涛等[20]研究了多匝小线圈的自感和互感引起的接收信号失真问题，采用曲线偏移的方法修正信号；杨海燕等[21]根据场源提出了一种新型圆锥型瞬变电磁场源装置，可减弱线圈间的互感耦合，提高浅部探测能力。总之，在浅层瞬变电磁法中，收发线圈的互感尚没有较理想的解决方案，因此导致的早期信号失真而产生浅层探测盲区的问题依然存在。

为了消除一次场对接收线圈的影响，实现瞬变电磁早期纯二次场响应观测，提出了等值反磁通瞬变电磁法(opposing coils TEM，简称OCTEM)。这是一种新型的瞬变电磁法，其原理与传统的瞬变电磁法一致，不同的是：OCTEM以等值反磁通的规律作理论依据，采用通以大小相同、方向相反电流的大小相同、平行共轴的天线作为发射源发射一次场，在反向对偶磁源的正中间平面接收二次场信号，测量地下纯二次场的响应[22]。根据矢量叠加原理，正负反向磁源产生的一次场始终存在至少一个磁场的垂直分量为零的平面，该位置称为等值反磁通的零磁通面。由于需要将接收天线设置在所述的等值反磁通零磁通面，因此对于等值反磁通瞬变电磁法的零磁通面的位置确认就显得尤为重要。目前，等值反磁通瞬变电磁法的零磁通面的位置确认一般采用机械调节接收天线的位置，通过目测二次场衰减曲线早期没有比较明显的过渡过程(明显畸变)，来定性判断零磁通面的位置。但是，这种方法主观性较大，无法科学、准确地确认零磁通面的位置。为此，本文提供了一种能够快速、准确对零磁通面位置进行确定的方法。

1 零磁通调零方法原理

根据上述频率的取值规则，本次实验以20 kHz频率为例，在等值反磁通调零模式发送20 kHz双极性方波，使用式(1)对采集的数据进行离散傅里叶变换，计算出0～40 kHz的频谱数据:

。

(1)

式中：xn为采集的离散时间序列;N为变换区间长度;k为谐波序号,k=0,1，…，N-1。

信号的平均功率Pavg为

(2)

根据帕塞瓦尔定理(Parseval's theorem)有：

(3)

可以得出：

。

(4)

考虑到功率谱密度(PSD)，平均功率有：

，

(5)

式中：Δf=fs/N，因此可得：

。

(6)

式中，PSD的单位为V2/Hz。

对PSD取对数，使振幅较低的成分相对高振幅成分得以拉高，以便观察掩盖在低幅噪声中的周期信号：

。

(7)

2 零磁通调零方法应用

本次测试使用自主研发的等值反磁通瞬变电磁仪。此电磁仪分为主机和天线两部分，天线设计采用微线圈发射(线圈直径不超过2 m)、中心回线接收的一体机装置，即将2个多匝小线圈分为2段串联安装在接收骨架上作为发射天线，正中间平面放置接收天线，上下线圈导线材料、绕制方法、匝数、大小都相同。采用微线圈发射可以在浅层得到高聚集度能量的一次场，能保证浅层探测的激励场，而收发一体机装置能保障每个测点发射一次场以及发射、接收相对位置的一致性，不仅大大减小了野外施工布线误差，还降低了由记录点位置不同导致的二次场误差。除此之外，线圈一体机装置的设计使在实际勘探工作中无需野外布线，既节省时间、降低成本，又能在高密度扫面测量工作中实现快速移动，提高了效率，有利于实现高精度的浅层勘探。

本次实验选取编号为1、2的2台等值反磁通瞬变电磁仪进行实验，发射50%占空比双极性方波，实际发射波形见图1，近似线性关断。

图1 发射电流波形

通过天线上的调零旋钮实时调整接收天线的位置，直至实时接收数据的功率谱密度A位于功率谱密度曲线的第N个极小值时(N为自然数)，此时接收天线的位置即为等值反磁通瞬变电磁法中零磁通面的位置。首先大幅调节接收天线的上下位置，同时观测采集数据的功率谱密度曲线，如果功率谱值逐渐增大，则向相反方向调节；当功率谱密度明显下降至一个极小值A1后上升时，则继续向同一方向调节天线位置，开始细调；当功率谱密度再次出现一个极小值A2时，此时改变调节方向，进入微调，直至对应的极小值为细调调零系数A3，应A3

对2台仪器进行上述的调零操作，实时调整天线上的调零旋钮，得到如图2所示的仪器随时间变化的功率谱曲线。在调零前后分别进行一次数据采集，并绘制成如图3所示的2台仪器衰减曲线对比。计算出2台仪器每个取样道调整前与调整后的数值差异的百分比，绘制随时间变化的百分比图(图4)。

图2 功率谱曲线

根据图3和图4可以看出，调零前的数据曲线在早期过渡现象明显，经过幅度起伏较大的过渡过程后数据曲线才变得平滑。经过调零后，早期的过渡过程得到了抑制，采集信号比较平滑，振荡减弱，对于中、晚期则无明显影响，能够符合野外工作条件，可以用于实际勘探任务。

(a)1号机 (b)2号机

(a)1号机 (b)2号机

3 结论

通过实时采集接收信号，并进行离散傅里叶变换得到频谱数据，获得功率谱密度随时间变化的曲线，并根据其变化确定零磁通面的位置。通过对实际研制的等值反磁通仪器进行测试，得出如下结论：

1)使用此方法确定零磁通面位置，精度高，测试方法简便可行。

2)使用此方法可准确找到零磁通面位置，抑制早期的过渡过程，减弱早期信号的振荡幅度。

3)此零磁通面位置确定方法可为高性能等值反磁通瞬变电磁仪器的设计和制造提供参考。