瞬变电磁在实际中的应用设计.docx

《瞬变电磁在实际中的应用设计.docx》由会员分享，可在线阅读，更多相关《瞬变电磁在实际中的应用设计.docx（21页珍藏版）》请在第一文库网上搜索。

1、引 言瞬变电磁法由于其具有抗干扰能力强、信息量丰富、勘探深度大，能克服低阻覆盖影响等特点，而受到国内外地球物理学家的重视。尤其近几年来，无论在仪器的改进与完善,理论研究及应用效果等方面都取得了很大进展。目前其应用范围也更加广泛。除了用于金属矿勘探外，还被应用于地质构造填图，石油、地热、煤田及地下水的探测。但瞬变电磁法的数据量大，测量结果受工作方式、回线大小、发射波形等因素的影响，解释工作比较困难。因此实际工作需要一种比较有力的解释方法，使人们能够从测得的瞬变特性曲线了解地电结构情况。传统的方法都是根据均匀半空间晚延时的瞬变特性,利用简单的公式把测得的电位(03/次)一一时间曲线转换成视电阻率一

2、一时间曲线。以此分析，推断地电断面的结构分布情况。这种方法从某种意义上说是对测量的电位用时间归一，它虽然有助于人们对地下地质情况的了解，但由于结果没深度概念，所以要据此来具体解释地下地电结构还是很困难的，最后得出的结论是很粗糙的，只是定性的分析,而没有定量的结果。尤其在实际工作中，由于受人文噪声、仪器功率等方面因素的影响,有时采样区间不够宽，那么计算的视电阻率曲线就不能反映地下电阻率的变化特征。因此,随着瞬变电磁法在各个领域的普及，工作量的大量增加，近年来人们都在寻求一种既快速又比较准确地解释方法，这种方法，这种方法能使人们很快地测量瞬变特性曲线了解到地下地电断面的变化情况，这就是把实测的电位

3、一一时间曲线转换成视电阻率一一深度曲线。当然，这里的电阻率值及深度大小不一定非常确切，但由于有了量的概念，对解释工作来说方便多了，它能使解释水平有一个很大的提高。最早从事这项工作研究的是Raicho,他在1985年导出了均匀半空间中感应电流线垂直散射速度及穿透深度，然后换算成视勘探深度，但结果不是非常理想。James和Macnae等人在1987年提出了用镜像技术处理TEM数据，把电位一一时间曲线转换成了视电阻率视深度曲线。同年A . G . Nekut基于同样基本思想，结合大地电磁法中的Bostic反演理论提出了时域电磁数据的直接反演方法。紧接着P . A . Eaton和Hohmann根据地

4、下感应电流线的散射速度及电阻率变化比较敏感的特点提出了瞬变电磁测深的近似反演方法，三者提出的方法的基本原理都是基于Nabighian (1979)的“烟圈”理论，只不过在具体的处理方法上有所不同，所以他们给出的结果基本相同的。本文就“烟圈”理论进行介绍，进而在介绍在实际应用中讨论瞬变电磁信号在地下与地面的响应特征。摘要：“烟圈”理论最早由Nabiphian提出（1979）。在地面观测的瞬变电磁场是大地感应涡流产生的，Nabiphian, Hoversten和Morrison, Reid和Macnae的研究表明，这个涡流可以近似地用圆形电流环等效。本文将较为详细介绍其理论原理，进而模拟计算不同时

5、刻瞬态电场在地下的分布形态基地上面感生电动势相应的变化，揭示了低阻异常体对感应涡流的聚集作用，地租覆盖层对瞬变场扩散减速作用，及瞬变场的延时效应。因此，瞬变电磁法对低阻体式敏感的，有上覆低阻层是探测同样的深度需较长时间，而延时效应瞬变长的晚期时段可反映埋藏较浅的异常体。关键词：“烟圈”概念 瞬变电磁场 散射速度 穿透深度第一章：电磁场基本理论一、麦克斯威尔方程组麦克斯威尔方程组的物理意义：建立了场强矢量、电流密度及电荷密度之间的关系。麦克斯威尔方程组描述了电磁场最基本的规律，在时间域中可表示为:V x HQJ） =J） + 00了）（安培环路定律）dtV X Eg）= -2（法拉第感应定律）d

6、tV Bg） = 0（ 5是涡旋场高斯磁场定律）V 。（乙，）=小心，）（库伦定律即高斯电烧律）其中，E为电场强度（V/m） B为磁感应强度（Wb/m2）, H为磁场强度（A/m）, D为电位移矢量（C/m2）, J为电流密度（A/m2）, J0为一次场电流密度（A/m2）, P自由电荷密度（C/m2） POw为一次电荷密度（C/m2）。如果考虑一次场源的作用，第一式和第四式右端分别应加上J0和P 0o物质条件：j = bE,B = JUH,D = E在导电介质中，物质条件可认为：y = o,VD = O在各向同性均匀介质中麦克斯威尔方程组变为:乂“(几/)二龙+迎8dtdBg)V x =dt

7、VeB(r,0 = 0V D（r,r） = 0为了求解电场和磁场，引入矢量位函数，电场和磁场可由相应的矢量位函数的微分直接求出。矢量H是涡旋场，故可表示为矢量A的旋度，将下式代入麦克斯威尔方程中，利用矢量恒等式和罗伦兹条件化简。即得到电磁场的矢量位和标量位所满足的时间域波动方程。矢量位函数：n VxAv2a =耳id2A dA-+ U(Tdt2 dtv2o =耳i+ U(J dt时间域波动方程:罗伦兹条件：+ 日电dt在不同条件下求解波动方程，只需将相应的条件（有源还是无源，是电性源还是磁性源,位移电流和传导电流同时考虑还是忽略其中一项）代入方程并求解即可。麦克斯威尔方程组的边界条件：1 .磁

8、感应强度B的法方向分量连续；2 .电位移矢量D法向分量的差等于电荷密度连续； （/1一，2）二P5n(E1 -E2) = 0x(H”2)= 0n _ 12)= 一冬dt3 .电场强度E的切方向分量连续;4磁感应强度H的切向分量关系;5 .电流密度的边界条件；6 .位函数的边界条件；地球物理中的标量位描述的似稳场，界面时标量位连续，电流法向边界条件的右端为零。二、中心回线方式下的响应为了求解波动方程，当在矩形同线（特别是正方形矩形同线）中心观测时，即发射、接收为中心回线排列方式下，矩形回线的响应可近似为与其面积相同的圆形回线的场，而圆形回线的场易于求解。对于圆形回线的场的求解是将回线内一小块面积

9、看成一垂直磁偶极子，以求得的垂直磁偶极子的场为基础。即设：dm = Ipdpdcp9用R=|P 一P”代替，将垂直磁偶极子计算结果代入上式，得回线元产生TE的波的矢量位：R PXdp ,P八adFlp, z)=侬吸迎卜 f *R, z) J(kpR)dk诙里有户应,z） =3efi+/s）这里有：5在整个回线内求面积积分得到回线产生的总矢量电位： j ooa 2 兀方（p,Z）=詈户（勺,z）J J JkpR）pdpd（pdkp4 oo o由矢量电位经微分运算可得到电磁场：纥=-Hm = （cr + icos）F + -VVFicOjLl令 o=uo,k p =入；得电磁场分量为:纭=竺出生j

10、 /“山+）一X. 4 （而）J（九夕）32o一o（z+）-%/(z _)卬 I (即)T 隼22H_ =_6。+)+ %。一力一 J(Aa)J0(Ap)dA2 ioEp=Ez=E=2对电测深方法来说，为了减少横向不均匀性的影响，我们常只测量发射回线中心处磁场的垂直分量，令圆形回线方式响应方程中P为零，可得到这种所谓“中心回线”排列的响应。园形同线方式响应：Hz 4（%）/。（2夕）“2 io在距离远小于自由空间波长的条件下，准静态近似都成立（kO 0）o取u0=入；P=o中心回线排列的响应：一“。q+力）+ 方 e（zi） J故 %/。(O) = 1当发射和接收装置均置于均匀大地表面，z=h

11、=O;Jx_/若口I下=K J 1 + 方2 ou0将uo=入Hz =代入，得到：:22la fJ, (Aa)dA+ w 1又因为一阶贝塞尔函数与零阶贝塞尔函数关系如下:AZ（助）= g（J0 （曲）daHz - la- J0（xzz）aX da A + u丁一/ 二左2（左2 _ _jo3R 0）分子分母同乘以（A-/）项Hz=Ia茄,丁丁 （）Hz=Ia22Hz =la d00CO 1- f , J()初u* 2 - w0077 J 尤，。(2a)d2 - jk2 da利用积李普希兹积分利用索墨菲尔德积分八一ikcieOO QJe-uz JA()u升=a/ 夕？ + z?la d d2 1

12、52 eika/瓦左点（一T）屋。3-(3 + 3ika-k2a2 )e-ikak q对于发射和接收装置均置于地表时，可以将频率域转换时间域，做S=i3代换，将方程除以S后，进行下列变换对方程进行反拉普拉斯变换，有：a = (/cr)1ik3 a3 = a3s32lea1 = -a2s ika = asX2乩=广2=/-高 P-(3+心+)厂1sk a s s sr 1_as2， CC 、L -e 、 = erfcQ-ya/2=1,2/4,S，?J (Ry厂(p_a、2 y =aa/.2万1/2产20 0cl erfcOa = 1 一 erf (6心=1 - J eu du得到回线中心的磁场阶

13、跃响应：/2-H.= 3erf(6a) - -%(3 + 3% 202a2 )夕”2。a 2万 cV = .icoRo0 0Hz(a, z)ad(pda12V = -iCDI f )+ 方，。5) 土 词0“000 12V = 一,斗0m2/ fj (曲)” + u 1在晚延时段内，e =4。很小时，有:vt = 0aJo/a1 3/2(vC)r= u(3 + 3u 22)e 11J不ottt 3/2 3/2 2T 7dH -I(y o a s/2y=_w_.w(_r -)j 3/25/2 2V Jb4a /5/2V J1 / C i20-2第二章烟圈的概念及其基本参数概念的引出敷设在均匀大地

14、上的发送回线中的电流产生磁力线，当该回线中电流突然切断的瞬间,在地下导电半空间中产生感应涡流场，开始瞬间的电流集中于地表，随后向下及向外扩散,等效电流:扩散半径:所在深度:v/d _2向下扩散速度：v 0 =式中：t为时间，s；。为电导率，S/ni；为空气的磁导率，1X4 71 X10 7H/m； c2= 8 / ( ji 2 )由以上各式看出，随着时间的增加，电流线的深度随之增大，同时其半径也相应增大。换句话说，导电半空间上的瞬变响应可近似地用圆电流线系统表示，这些电流线就像由发射回线“吹出”的“烟圈”。这些“烟圈”虽时间的增加向下向外运动(图1)。这就提示我们Z图1_ 4 t当计算均匀半空间的瞬变电磁响应时，可以用一个位于z =