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

基于能量的激励波形电磁响应研究 瞬变的正演计算 空间电阻器，也称为半航空电阻器，通过地面传输和航空接收。因此，它具有航空和地面瞬变电阻器的优点。同地面瞬变电磁系统相比,它具有野外布线方便快捷、探测高效等优势;较航空瞬变电磁系统,具有信噪比更高、空间分辨率更好、勘探深度更深的优点,可快速实现较大面积内深部异常体的探测,较适用于我国地形复杂的山区资源探测。 最早的地-空电磁探测系统(TURAIR)出现在20世纪70年代的加拿大和俄罗斯,它是一种频率域的电磁探测系统。发射线圈直接由发电机供电,但是在文献中并没有明确说明激励场的波形。1991-1993年间,由澳大利亚研制的固定线圈地-空瞬变电磁系统(FLAIRTEM)和1997年加拿大人研制的地-空电磁探测系统(TerraAir)都是时间域电磁探测系统,激励波形均为方波。所不同的是TerraAir系统进行的是全波形的瞬变电磁响应计算而FLAIRTEM系统只进行了电流关断后的电磁响应计算。在1992年,日本基于地面长偏移距瞬变电磁法(LOTEM)提出了GREATEM系统,它采用接地长导线源发射,感应线圈在空中接收的工作方式,具有布线快速和高效的优点。Richard S. Smith在2001年对航空、半航空和地面电磁系统的数据进行了对比分析。 1998年,Guimin Liu研究了磁性源、自由空间下,不同电流发射波形对航空瞬变电磁响应的影响,但仅对电流完全关断后,即电流为零段(off-time段)的电磁响应进行了研究。2006年于生宝等分析了阶跃关断与线性关断早期瞬变电磁响应的区别。2007年嵇艳鞠等进行了浅层瞬变电磁法中全程瞬变场的畸变研究,激励信号是斜阶跃波。2010 年,许洋铖等研究了阶跃波关断的航空时域电磁法的初始场计算。2008年Yin, C.计算了磁性源,发射波形为半正弦波和梯形波的全波形航空感应电动势。 在前面已有研究的基础上,利用地-空系统阶跃波瞬变电磁响应与任意激励波形卷积的方法,研究了电性源,激励波形为方波、半正弦波、三角波和梯形波情况下电流不为零段(on-time段)和off-time段对应高阻异常和低阻异常的地-空瞬变电磁响应影响;并且提出了激励能量概念,通过激励能量大小分析了不同激励波形响应的大小。 1 计算参数说明 为了计算不同激励波形的地-空瞬变电磁响应,首先需要计算阶跃波关断后的地-空电磁响应。图1为时间域地-空电磁探测系统示意图。接收线圈距地面高30 m,接收线圈面积200 m2,导线长1 km,电流为200 A. 纳比吉安给出了有限长接地导线层状大地的频率域磁场垂直分量表达式 式中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 激励电流半周期 为了计算任意电流波形的瞬变电磁