可控源电磁测深数据场源效应及反演讨论.pdfVIP

物探化探计算技术 33卷 V盯()，，z)=矿()，，z)一IT。(y,Z)。将总场分离为主设为零，利用等式(13)求解初始场值E；③在等 场和二次场分量，这样解决了电磁场在电源和磁源 式(14)中带人求得的岔：计算B：；④最后使用碰来 附近的奇异性，也保证直接计算二次场的质量。主 估计近似场值罡。重复以上过程直到迭代收敛， 场值可以通过简单电导率结构计算得到。 收敛条件为地表最大和最小二次场值改变量小于 通过对x方向场分量的傅式变化可以将三维 给定的值。 计算减少为二维计算。傅里叶变换定义如下： 如果获得二次场E和或，那么就可以将二次 (5) 场和主场相加得到x方向的总场分量这样其它两 F(k，，，，，z)=f=。dxe：F(菇，Y，：) 个水平分量可以通过下式计算： F(茹，Y，z)=磊1『-。dk，e嘶P(后，，Y，：)(6) ～ 一托。J 其中 k。是沿着构造走向的波数；带上标的表示 B／(k,，Y，：)5而tWg』,O丽 复试变换域的值。在笛卡尔坐标系下方程(3)方 7 妒—■；■一¨七*—■广∥u∽掣“，掣)(15) 程(4)可以表示为：、 一ik，E—a应：一z彰=o (7) ’ L —B。J 岛(吃，y,z)2而tWl比士V西 (8) a,B：一a,B；一砌盯E=脚△盯醒 16， a。E+ik，￡：一iwB：=0 (9) (础。掣椭掣)c 最后，通过傅里叶逆变换(6)求得空间域(菇， 一以q—a，彰-muE：=PC△鸩(10) a,B：一a,B：一iwB：=0 (11)Y，：)的EM场值。 (12) a,B：+ik，台：一m盯台；=m△口彰 如果源在菇方向上不变(二维时k，=0)正如2模型正演计算结果 MT或者无限线源的情形。方程(7)到方程(12)可 以分解为两种模式。方程(7)到方程(9)定义为横 元正演算法。建立电性模型进行正演模拟。 电模式(TE)，只有在x方向的电场。方程(10)到 模型A：looQ·m均匀半空间下一个 方程(12)定义为横磁模式(TM)，只有x方向的磁 场分量。对于有限源，在傅里叶变换域中k，是非 10Q·m的低阻体模型，模型几何特征见图l频 零的。TE和TM模式不再是相互独立的模式，而点范围：2 上方均匀分布，测点间距100m。电偶源与着构造 是在x方向上耦合的。通过对(7)到方程(12)重 新整理后在傅式域沿着构造走向的二次场表示 走向一致。 为： V．严芒)-ttotrvE= 、 ’y 7 胁衍E：一溉(墨半)+ 、 ’y 7 ik,(咄：×V(专))％(13) 图1正演模型A