关断电流校正的意义.PPT

* * ATTEM瞬变电磁系统 数据处理软件中关键技术的研究 报告人：嵇艳鞠 导师：林君教授 学术报告 一、时间域电磁场理论 二、瞬变电磁场关断电流的校正方法 三、感应电动势变磁场信号 四、一次场同二次场分离 在ATEM瞬变电磁系统数据软件研制的基础上，进一步开发适用于全程探测ATTEM系统的数据处理软件。为了解决探测浅部，提高浅层探测的分辨率。在发射机的电流关断不为零时，关断电流校正只是对电流延迟的影响的校正，关断电流的精确校正将提高早期阶段的精度和分辨率，提高ATTEM的浅层探测能力。而在电流衰减期间，记录的信号包含了一次场和二次场的叠加，如何消除一次场的影响变的至关重要。如果从记录的总场中消除一次场，在数据处理时，就可以使观测和解释扩展到全时域。 其中，E为电场强度（V/m）B为磁感应强度（Wb/m2），H为磁场强度（A/m）， D为电位移矢量（C/m2），J为电流密度（A/m2），J0为一次场电流密度（A/m2），ρ自由电荷密度（C/m2）ρ0w为一次电荷密度（C/m2）。如国考虑一次场源的作用，第一式和第四式右端分别应加上J0 和ρ0。 在各向同性均匀介质中麦克斯威尔方程组变为： 一、电磁场基本理论： 麦克斯威尔方程组描述了电磁场最基本的规律，在时间域中可表示为： 麦克斯威尔方程组的物理意义：建立了场强矢量、电流密度及电荷密度之间的关系。 物质条件：? 在导电介质中，物质条件可认为： 一、麦克斯威尔方程组 为了求解电场和磁场，引入矢量位函数，电场和磁场可由相应的矢量位函数的微分直接求出。矢量H是涡旋场，故可表示为矢量A的旋度，将下式代入麦克斯威尔方程中，利用矢量恒等式和罗伦兹条件化简。即得到电磁场的矢量位和标量位所满足的时间域波动方程。 在不同条件下求解波动方程，只需将相应的条件（有源还是无源，是电性源还是磁性源，位移电流和传导电流同时考虑还是忽略其中一项）代入方程并求解即可。 矢量位函数： 时间域波动方程： 罗伦兹条件： 麦克斯威尔方程组的边界条件： 磁感应强度B的法方向分量连续； 电位移矢量D法向分量的差等于电荷密度连续； 电场强度E的切方向分量连续； 磁感应强度H的切向分量关系； 电流密度的边界条件； 位函数的边界条件； 两媒质电导率有限时，传导电流和位移电流均为有限值。应用地球物理中的问题属于此中情况。 地球物理中的标量位描述的似稳场，界面时标量位连续，电流法向边界条件的右端为零。 为了求解波动方程，当在矩形回线（特别是正方形矩形回线）中心观测时，即发射、接收为中心回线排列方式下，矩形回线的响应可近似为与其面积相同的圆形回线的场，而圆形回线的场易于求解。 对于圆形回线的场的求解是将回线内一小块面积看成一垂直磁偶极子，以求得的垂直磁偶极子的场为基础。 即设 用R= |ρ－ρ’ |代替，将垂直磁偶极子计算结果代入上式，得回线元产生TE的波的矢量位： 在整个回线内求面积积分得到回线产生的总矢量电位： 这里有： 二、求解中心回线方式下的响应 对电测深方法来说，为了减少横向不均匀性的影响，我们常只测量发射回线中心处磁场的垂直分量，令圆形回线方式响应方程中ρ为零 ，可得到这种所谓“中心回线”排列的响应。 令ξ0==u0,kρ =λ； 得电磁场分量为： 由矢量电位经微分运算可得到电磁场： 当发射和接收装置均置于均匀大地表面，z=h=0; 园形回线方式响应： 在距离远小于自由空间波长的条件下，准静态近似都成立（k0 ≈0）。取u0=λ；ρ=0 中心回线排列的响应： 将u0=λ 代入，得到： 又因为一阶贝塞尔函数与零阶贝塞尔函数关系如下： 分子分母同乘以（ ）项 利用积李普希兹积分 ， 利用索墨菲尔德积分 对于发射和接收装置均置于地表时，可以将频率域转换时间域，做 s=iω代换，将方程除以s后，进行下列变换对方程进行反拉普拉斯变换，有： ?? ? 得到回线中心的磁场阶跃响应： ? ? 在晚延时段内， 很小时，有： 二、关断电流的校正方法 方法一、坐标移动法 斜阶跃波后沿的主要影响是响应值比垂直阶跃波的响应值低（如延时起点从后沿终点起算），且与延时 有关，故最简单的改正方法是坐标移动法。 M.Wasten和D.G.Price提出三种坐标移动法，如下： 关断电流校正的意义： 消除电流延迟的影响，得到精确的全程视电阻率 。 关断电流的精确校正将提高早期阶段的精度和分辨率，提高ATEM的浅层探测能力。 根据D.V.Fitterman 和 W.Anderson 从Duhamel积分出发对此效应的推导公式，对于均匀