瞬变电磁法浅层探测技术
第21卷 第2期
2006年4月
CHINESE JOURNAL O F RADIO SCIENC E
电 波 科 学 学 报
ol. 21, No. 2 V Ap ril, 2006
文章编号 1005-0388(2006) 02-0284-04
瞬变电磁法浅层探测技术
于生宝 王 忠 嵇艳鞠 林 君
①
(吉林大学智能仪器与测控技术研究所, yusheng bao 2000@163. co m, 吉林长春130026)
摘 要 从理论上推导发射电流关断时间对瞬变电磁法勘探结果的影响, 给出均匀半空间大地模型重叠回线互阻抗计算方法, 分析阶跃关断与线性关断早期瞬变响应的区别, 讨论了线性非零关断瞬变电磁法发射信号对浅层探测结果的影响。利用
IGBT(绝缘栅双极型晶体管) 实现瞬变电磁法发射电路, 减小发射机关断时间, 使之接近理想的阶跃函数, 消除非零关断的影响, 解决了加大发射电流与快关断的矛盾。实验结果表明, 减小发射电流关断时间, 可以提高瞬变电磁浅层勘探结果分辨率, 减小瞬变电磁法勘探盲区。
关键词 瞬变电磁, 关断时间, 探测盲区中图分类号 TH 762 文献标识码 A
Time domain transient electromagnetic shallow prospecting
YU Sheng -bao W ANG Zhong JI Yan -ju LIN Jun
(I nstitute of I ntelligent I nstr ument and M easur ement T echnology , J ilin Univers ity , y us heng bao 2000@163. com, Changchun J il in 130026, China)
Abstract The effect of turn -o ff tim e of tr ansm itting cur rent on transient electro -magnetic ex ploratio n is derived theo retically , the mutua-l im pedance of coincident loo p in uniform half -space is evaluated, the difference effect of step turn -off and lin -ear turn -off on early tr ansient signal is analyzed, and the effect of linear nonzero turn -off of transmitting current on shallow subsur face sounding is discussed. A in -sulated -g ate bipo lar transistor (IGBT ) is used in tr ansient electromagnetic trans -mitter sy stem to sho rten tur n -off time. The turn -o ff w avefor m is desig ned to ap -pr oach ideal step function to elim inate the effect of nonzero turn -off and solve the conflict betw een streng thening tr ansm itter and rapid sw itching. T he ex perimental results indicate that sho rtening turn -off time can enhance reso lution o f transient e -lectrom ag netic for shallow subsurface sounding and reduce the blind zone o f tr ans-i ent electrom agnetic so unding.
Key words transient electro magnetic, turn -off tim e, blind zo ne
1 引 言
电磁方法是浅部环境与工程探测的重要手段, 就目前来说, 所有浅部探测方法均是由相应的深部
探测方法发展而来的[1]。电磁方法分两大类:频率域电磁方法(如MT ) 和时间域电磁方法(如TEM ) 。近地表目标探测对于电磁方法来说是一个难题, 而这一范围却与人类活动极为密切, 地下2~100m 左
①收稿日期:2004-10-09
第2期 于生宝等:瞬变电磁法浅层探测技术右的深度范围是频率域方法的一个弱视区(或半盲区) 。对于时间域来说, 探测深度是由观测时间的早晚决定的, 目前时间域瞬变电磁法(TEM ) 仪器记录的几乎全部是晚期信号, 或者在数据处理时仅采用了晚期信号, 这种情况将产生两种后果:第一, 它损失了T EM 方法探测浅部结构的能力, 因为浅部结构的信息主要由早期信号携带; 第二, 它降低了TEM 方法的分辨能力, 因为关断电流的影响将使瞬变响应发生畸变[2], 因此, 从地表到地下20m 左右也是瞬变电磁法探测的半盲区。瞬变电磁法发射机输出信号的关断沿近似为抛物线或指数函数, 用实际的关断沿进行解释非常困难, 实际中都是用理想的阶跃函数近似解释, 但结果会与实际情况有很大差异。要探测浅层目标体, 有效的办法就是减小发射电流的关断时间。
[2]
-1
285
普拉斯变换, L p 为关于p 的拉普拉斯逆变换, J 为贝赛尔函数[3]。
当发射线圈与接收线圈用相同的回线时(半径为a ) , 即采用重叠回线, 引入几何因子:
Q( , K ) =J 1( ) =J 1( a) J 2( b ) K = /a, = a 代入方程(1) 并简化为:Z(t) =- a
2
1
L -p [I (p ) p A 0]Q( , K ) d (2)
线性关断电流如图1(b) 所示, 其表达式为:1t
I (t) =
2t 2-t 1
t 1t 2
2 非零关断的影响
瞬变电磁法中一般假定电流关断是一个理想的阶跃函数, 如图1(a) 所示, 但实际上不可实现, 由于发射线圈的非零闭环响应和发射机特性决定, 电流的下降需要一定时间, 下降方式可能是线性, 也可能是按抛物线或指数规律, 这里只就线性关断进行讨论, 如图1(b) 所示。通常瞬变电磁法的算法都是与发射电流和关断时间有关, 且假定电流关断如图1(a)
所示。
I (p ) =L p [I (t) ]
-pt 1-pt 2[1-+=
p (t 2-t 1) p (t 2-t 1) p
代入方程(2) 中并设:
1
G 0(t) =L -p (A 0)
G 1(t) =L p [p A 0e 1](t 2-t 1)
1-1
G 2(t) =L -p [p A 0e -pt 2](t 2-t 1) -1方程(2) 可以简化为:Z(t) =- a
-1-1-pt -1
[G (t) -G (t) +G (t) ]
1
2
Q( , K ) d (3)
在实际地球物理勘探中, 相对于t 0来说t 1是足
够大的, 初始时间一般选择t 1=0, 其中G 0(t) 是稳态电流I 0响应, 与关断时间无关, 可以忽略。平移电流波形坐标原点至(t 1, I 0) 点, 方程(3) 可以化为:
- p
) -1A 0(1-e Z(t) = a 0L p []Q( , K ) d
(4)
其中: =t 2-t 1
1
将F(t, ) =-L -p [
]代入方程(4) 可以简化为:p
图1 阶跃关断与线性关断波形
Z(t) =[F(t - , ) -0
F(t, ) ]Q( , K ) d
(5)
在层状大地介质中, 半径为a 和b 同心圆形线圈瞬变电磁互阻抗Z (t) 计算公式为:
Z(t) =- ab
[3]
1
L -P [I (p ) p A 0(P, p , ) ]
J 1( a) J 2( b) d (1)
(据A P Raiche. 1984. ) 。式中A 0(P , p , ) 为层状结构大地阻抗函数, P 为厚度与其电阻率的乘积, p 为与-j 对应的拉普拉斯算子, 为角频率, 为汉克尔逆变换的积变量, I (p ) 为归一化的电流波形拉
若q = a p / =t/ a
2
-12
由L p [f (p ) ]=) 知:2F( a
21010F( ) =-L -q []=-L -p []
q p
[3]
根据A P Raiche 的计算结果, 使用同步线圈, 并进行拉普拉斯变换、展开贝赛尔函数, 方程(5) 最终可以简化成如下形式:
Z(t) =t
222
286
电 波 科 学 学 报 第21卷
k=0
k k
T z ( , k)
k !(k +1) !(k +2) !(2k +5) !
(6)
3 测试与应用
在地矿部 九 五 重点项目《同步测量技术及其在电法中的应用》中采用IGBT 及双桥式逆变电路实现了瞬变电磁发射机, 关断时间达1. 2 s(纯电阻负载) 。关断时间测试结果如表1(测试仪器:TEK -TRONIX 100M 示波器; 输出电压20V; 负载4 100m 线圈; 输出频率12. 5Hz; 测试次数20 6次) 。
表1 关断时间测试
负载电阻( ) 输出电流(A) 关断时间( s)
7336
5448
4564
2. 8780
1. 513160
120400
其中:
T z ( , k) =
-1k+
(1- ) 2
(k +)
2
2
= /t X = a /4t 当 →0时, T z ( , k) →1, Z (t) 即为阶跃函数关断的响应。
从方程(6) 中可以看出关断时间t 0=t 2-t 1对早期互阻抗有较大影响, 如图2所示, 视电阻率为1000 m 的均匀半空间模型在不同关断时间的响应, t 0=0的曲线反映了视电阻率早期和中期内的下降特性, 最后接近均匀半空间的电阻率, 随着关断时间t 0的增加, 早期视电阻率误差也逐渐增大, 然后随时间增加逐渐接近均匀半空间电阻率, 图3给出同一均匀半空间模型线性关断与阶跃关断视电阻率比值(K t ) 随时间变化曲线, 更直观地给出线性关断电流的关断时间对早期视电阻率的影响
[4]
为了验证系统在浅层工程探测中的有效性, 在长春市郊某地进行了城市活断层探测。活断层是导致多种地质灾害, 特别是地震、崩塌、滑坡、地表错断的重要根源, 对地面设施、建筑破坏力极大, 因此, 城市活断层探测对于工程安全性评价、减灾防灾具有重要意义。为验证瞬变电磁法发射电流关断时间对浅层勘探结果的影响, 采用快关断浅层瞬变电磁探测系统对该活断层进行了探测, 由于测区宽度有限, 只能铺设10m 的线框, 采用了重叠回线装置, 发射回线为边长10m 的正方形线圈, 接收回线的面积100m 2, 点距为10m, 发射电流为6. 2A, 发射频率为25H z, 关断时间t 0=18 s 。并与常规瞬变电磁探测结果进行了对比。图4给出常规瞬变电磁法解释结果与浅层瞬变电磁法解释结果对比图。从视电阻率的断面图中可以看出浅层瞬变电磁探测方法对近地表的电性变化信息有了一定反映, 从测线680m 处高程为-10m 开始至测线640m 处高程为-50m
之
。
图2
均匀半空间模型关断时间对视电阻率的影响
图4 蔡家村城市活断层常规瞬变电磁结果(上, 据曾昭发2003) 与浅层瞬变电磁法
解释结果剖面(下) 对比
图3 均匀半空间模型线性关断与阶跃关断视电阻率比
第2期 于生宝等:瞬变电磁法浅层探测技术间, 断层反映清楚。勘探结果在距地面10m 以下准确反映了活断层的信息, 比常规瞬变电磁法提高了浅层分辨能力。
287
[8] 嵆艳鞠. 浅层高分辨率全程瞬变电磁系统中全程二次
场提取技术研究[D].长春:吉林大学博士论文, 2004. 40~55.
[9] Hesham M. E-l K aliouby and Essam A Eldiw any z.
T ransient electro magnetic responses of 3D po lar izable body[J]. G eophy sics, 2004, 69(2) :426~430. [10] 白登海等. 时间域瞬变电磁法中心方式全程视电阻率
的数值计算[J]. 地球物理学报, 2003, 46(5) :697~704.
[11] 唐新功, 胡文宝. 层状地层中三维薄板的瞬变电磁响
应[J]. 石油地球物理勘探, 2000, 35(5) :628~633.
[12] 李 貅. 瞬变电磁测深的理论与应用[M ]. 西安:陕西
科学技术出版社, 2002, 75~90.
[13] 闫述等. 瞬变电磁场的直接时域数值分析[J].地球
物理学报, 2002, 45(2) :276~
282.
4 结论
传统电磁法一般假定关断电流波形是理想的阶跃函数, 电磁法发射机电流关断时间大约为50 s~300 s, 如果忽略了关断时间, 解释结果中早期视电阻率就会产生很大误差, 如图3推导结果及图4实际应用所证明, 关断时间越长, 影响越大, 被测对象电阻率越大, 影响也越大, 关断后, 影响随时间增加而减小[4]。为了使解释结果准确可靠, 一种方法是在处理瞬变电磁数据时进行校正, 从数据中消除关断时间的影响, 通过曲线匹配进行数据解释。另一种方法尽量减小发射机关断时间, 使之接近理想的阶跃函数, 消除非零关断的影响。参考文献
[1] M N Nabig hian. Elect romag net ic M ethods in A pplied
Geophysics[M ].T ulsa:So ciety o f Explorat ion Geo -physicists, 1987, 402~431.
[2] 白登海. 瞬变电磁法中两种关断电流对响应函数的影
响及其应对策略[J]. 地震地质, 2001, 23(2) :245~251.
[3] A P Raiche. T he effect o f ramp function turn -off on
the T EM response of lay ered earth [J].Ex plor ation Geophysics, 1984, 15(1) :37~41.
[4] D V F itter man, W L Anderson. Effect o f tr ansmitter
tur n -off t ime o n tr ansient soundings [J].G eoex plor a -t ion, 1987, 24(2) :131~146.
[5] 朱凯光等. 基于M A T L A B 的高频谐变电磁场的数值
计算[J]. 电波科学学报, 2002, 17(3) :224~228. K G Zhu, et al . . N umer ical co mputation o f hig h fr e -quency elect romag net ic fields under M A T L AB [J ]. Chinese Journal o f Radio Science, 2002, 17(3) :224~228.
[6] 于生宝等. 瞬变电磁法中发射机关断时间的影响研究
[J].石油仪器, 1999, 12(6) :15~17.
[7] 程德福. 近区磁源瞬变电磁法信号检测技术研究[D].
长春:吉林大学博士论文, 2002, 17~35.
于生宝 (1963-) , 男, 吉林人, 吉林大学智能仪器与测控技术研究所, 教授, 博士毕业于俄罗斯远东国立技术大学, 主要参与完成
863计划、国家十五科技攻关、自然科学基金项目5项, 主持部级项目
1项, 学校项目1项; 获国家发明专利2项, 实用新型专利1项; 获教育部科技进步2等奖1项、吉林省
科技进步1等奖1项。主要研究方向:地球物理仪器。
王 忠 (1979-) , 男, 山东人, 博士研究生, 吉林大学智能仪器与测控技术研究所, 助教。
嵇艳鞠 (1972-) , 女, 黑龙江人, 博士, 吉林大学智能仪器与测控技术研究所, 副教授。