基于能量的激励波形电磁响应研究

基于能量的激励波形电磁响应研究

瞬变的正演计算空间电阻器，也称为半航空电阻器，通过地面传输和航空接收。因此，它具有航空和地面瞬变电阻器的优点。同地面瞬变电磁系统相比,它具有野外布线方便快捷、探测高效等优势;较航空瞬变电磁系统,具有信噪比更高、空间分辨率更好、勘探深度更深的优点,可快速实现较大面积内深部异常体的探测,较适用于我国地形复杂的山区资源探测。最早的地-空电磁探测系统(TURAIR)出现在20世纪70年代的加拿大和俄罗斯,它是一种频率域的电磁探测系统。发射线圈直接由发电机供电,但是在文献中并没有明确说明激励场的波形。1991-1993年间,由澳大利亚研制的固定线圈地-空瞬变电磁系统(FLAIRTEM)和1997年加拿大人研制的地-空电磁探测系统(TerraAir)都是时间域电磁探测系统,激励波形均为方波。所不同的是TerraAir系统进行的是全波形的瞬变电磁响应计算而FLAIRTEM系统只进行了电流关断后的电磁响应计算。在1992年,日本基于地面长偏移距瞬变电磁法(LOTEM)提出了GREATEM系统,它采用接地长导线源发射,感应线圈在空中接收的工作方式,具有布线快速和高效的优点。RichardS.Smith在2001年对航空、半航空和地面电磁系统的数据进行了对比分析。1998年,GuiminLiu研究了磁性源、自由空间下,不同电流发射波形对航空瞬变电磁响应的影响,但仅对电流完全关断后,即电流为零段(off-time段)的电磁响应进行了研究。2006年于生宝等分析了阶跃关断与线性关断早期瞬变电磁响应的区别。2007年嵇艳鞠等进行了浅层瞬变电磁法中全程瞬变场的畸变研究,激励信号是斜阶跃波。2010年,许洋铖等研究了阶跃波关断的航空时域电磁法的初始场计算。2008年Yin,C.计算了磁性源,发射波形为半正弦波和梯形波的全波形航空感应电动势。在前面已有研究的基础上,利用地-空系统阶跃波瞬变电磁响应与任意激励波形卷积的方法,研究了电性源,激励波形为方波、半正弦波、三角波和梯形波情况下电流不为零段(on-time段)和off-time段对应高阻异常和低阻异常的地-空瞬变电磁响应影响;并且提出了激励能量概念,通过激励能量大小分析了不同激励波形响应的大小。1计算参数说明为了计算不同激励波形的地-空瞬变电磁响应,首先需要计算阶跃波关断后的地-空电磁响应。图1为时间域地-空电磁探测系统示意图。接收线圈距地面高30m,接收线圈面积200m2,导线长1km,电流为200A.纳比吉安给出了有限长接地导线层状大地的频率域磁场垂直分量表达式式中u0=(λ2-k2002)1/2,在准静态条件下,u0=λ,因此式(1)可写为式中:R1=[(x-x′)+y2]1/2均匀半空间反射系数rTE=λ−λ2+iωμ0σ√λ+λ2+iωμ0σ√rΤE=λ-λ2+iωμ0σλ+λ2+iωμ0σ三层大地反射系数计算参数说明见表1.式(2)中有两重积分,内层积分存在振荡的Bessel函数,一般只能求出均匀半空间的解析式,如果求解两层及以上层状大地的电磁响应需采用数值积分方法。使用D.Guptasarma和B.Singh提出的数字滤波法解决内层积分的hankel变换问题,而外层积分利用辛普森积分法就可以实现。频率域到时间域的转换是利用D.Guptasarma提出的线性滤波方法,将-iωμ0sHz变化到时间域就得到了感应电动势−μ0s∂hz∂t-μ0s∂hz∂t.Nabighian,M.N.和Oristaglio,M.L.给出了均匀半空间接地导线赤道轴线上感应电动势的解析表达式,将计算结果与式(3)进行对比,以验证算法的正确性。Vz=−μ0s∂hz∂t=2sIπσy3{(1+θ2y2)e−θ2y2erf(θL)−LR2(1+y22R22)erf(θR)+θLy2π√R22e−θ2R22}(3)Vz=-μ0s∂hz∂t=2sΙπσy3{(1+θ2y2)e-θ2y2erf(θL)-LR2(1+y22R22)erf(θR)+θLy2πR22e-θ2R22}(3)式中:θ=(μ0σ4t)1/2θ=(μ0σ4t)1/2;R2=(y2+L2)1/2,计算参数如表2所示。从图2可以看出,应用解析式法和数值积分法计算得出的感应电动势衰减曲线重合的很好。2瞬变电磁响应任意激励波形的瞬变电磁响应为阶跃波瞬变电磁响应与电流一阶导数的卷积,如公式(4)为式中:V任意(t)为任意波形的瞬变电磁响应;I任意(t)为任意瞬时电流;V阶跃(t)为阶跃波的瞬变电磁响应。结合公式(4)及阶跃波地-空瞬变电磁响应即可求得层状大地任意激励波形的地-空瞬变电磁响应,并通过激励能量对瞬变响应进行分析。2.1激励电流半周期为了计算任意电流波形的瞬变电磁响应,图3给出所需要的四种发射电流波形,其中Ti为激励电流半周期,i=1,2,3,4分别对应于半正弦波、方波、三角波和梯形波。t1,t2,t4,t7为各激励电流的脉冲持续时间,u(t)为阶跃函数,δ(t)为冲激函数。对应的波形解析式和一阶导数表达式如下:1)半正弦波2)方波3)三角波4)梯形波2.2瞬变电磁响应计算不论发射机的有限长导线中通入的是哪一种电流,本质上都是一种能量,激励能量的通用表达式为E=∫t000t0I2(t)dt(13)式中:t0为脉冲持续时间;I(t)为发射电流。半正弦波能量方波能量三角波能量梯形波能量其中脉冲宽度取1/4倍半周期长,脉冲频率f=25Hz,因此,四种波形的半周期长均为0.02s.I1、I2、I3、I4分别是正弦波、方波、三角波和梯形波的电流最大值,均取200A.ω为正弦波的角频率。经计算可得各激励能量间关系为E方波=1.5E梯形波=2E正弦波=3E三角波(18)将等式(6)、(8)、(10)、(12)分别与阶跃波的瞬变电磁响应进行卷积,就可得出这四种波形的层状大地瞬变电磁响应,如图4所示。接收位置坐标为(0,100,30),接收线圈匝数为1,第一层电导率σ1=1/100S/m,第二层电导率σ2=1/5S/m,第三层电导率σ3=1/100S/m;第一层大地厚度h1=100m,第二层厚度h2=50m.图4给出了off-time段的瞬变电磁响应。在整个的off-time段,方波的响应是最大的,其次是梯形波和半正弦波,三角波的瞬变电磁响应最小,这一大小关系与不同激励能量间的关系是完全吻合的,见式(18),也就是说,激励能量大对应的响应就大。因此如果只想了解哪个激励波形所产生的响应大,就不需要依次对响应大小进行计算,而只需计算它们的激励能量大小。而在on-time段,当t=0.01μs时,方波的响应最大,为7.096×107nV,在0.06ms至1.24ms梯形波的响应最大,为7.471×105nV,在1.24ms到3.75ms三角波响应最大,为3.767×105nV,而从3.75ms至完全关断时刻t=5ms之间,响应的最大值变为梯形波最大,为7.46×105nV,因此激励能量与瞬变电磁响应间的大小关系只适用于off-time段。3对各波的响应以三层大地为模型,将地-空瞬变电磁响应分为on-time段和off-time段进行分析。根据电导率值的不同分为高阻模型、低阻模型和均匀半空间模型进行计算,除电导率外其余参数均与图4所取参数相同,电导率参数见表3.图5～图8分别为半正弦波、方波、三角波和梯形波的on-time与off-time段电磁响应。分析并得出如下结论:1)在on-time段,无论是低阻模型还是高阻模型,半正弦波与三角波的响应无明显变化。而对于方波和梯形波,低阻模型响应变化明显。在关断早期,也就是(0.005s,0.0053s)之间,高阻模型与低阻模型的响应混叠在一起,不利于进一步的反演。2)分析图6和8,地-空电磁探测系统的响应在off-time段,高阻模型比低阻模型衰减的快,如果只记录off-time段响应,当中间层电导率逐渐减小时,对于实测数据而言,高阻异常很快进入噪声区,因此无法得到足够多的高阻段有用信息。解决的办法是发射波形采用梯形波或方波,接收机同时记录off-time和on-time段的响应。3)半正弦波和三角波在on-time段的响应随电导率的变化不明显。在t=3.3ms时刻,半正弦波以均匀半空间模型计算的响应结果为比较基准,高阻模型和低阻模型的平均响应幅值差比仅为0.823%,三角波为0.839%,而方波和梯形波分别为139.41%和135.17%.因此,采用方波和梯形波激励时,记录on-time段和off-time段数据进行异常分析,可使反演分辨率更高。4瞬变电磁响应计算以阶跃波层状大地地-空瞬变电磁响应为基础,利用卷积方法计算了方波、梯形波、半正弦波和三角波的on-time和o