大地电磁偏移成像_宋维琪探究.pdfVIP

　第 27 卷　第 1 期 地 　震 　学 　报 Vol. 27 , No . 1 　 　2005 年 1 月　(102～108) ACTA SEISMOLO GICA SIN ICA Jan. , 2005 　 文章编号 : 025323782 (2005) 0120102207 大地电磁资料偏移成像 3 宋维琪1) 　邹文勇2) 1) 中国山东东营 257061 石油大学 (华东)地球资源与信息学院 2) 中国山东东营 255075 中石化胜利油田地质科学研究院　　 摘要 　在研究 Zhdanov 的二维电磁场偏移成像的基础上 , 针对它在某些方面的不足 , 研究并 提出了改进算法. 首先 , 用一维反演结果确定初始背景电阻率 , 然后在延拓计算过程中 , 利 用迭代算法进行逐层修正、延拓计算 , 从而获得较准确的电阻率 ; 第二 , 研究了有限差分方 程的改进算法 , 根据大地电磁场变化的特点 , 设计了纵向以等比级数变网格的迭代算法. 通 过以上两方面技术的改进 , 经理论模型和实际资料试算验证 , 取得了较好的结果. 关键词 　　大地电磁资料 　有限差分法 　数值计算 　偏移成像 中图分类号 : P318. 2 　　　文献标识码 : A 引言 电磁场偏移成像问题类似于地震勘探中的弹性波偏移成像问题. 但是二者又有本质的 不同 (Lee et al ,1989) , 弹性波场满足一般的波动方程 , 而电磁波限于我们讨论的方法满足 扩散方程. 不管是电磁波还是弹性波偏移成像的机理是类似的. 对于不同介质分界面 , 能 够体现此界面的像函数 map 为 (马在田 , 1989) map ( x , z) = u( x , z , tu) d ( x , z , td ) (1) 其中 , 函数 u和 d 分别为上行波函数和下行波函数 ; x , z 均为空间坐标 ; tu , td 分别为上下 行波的旅行时间 , 即入射波的初至和反射波刚刚开始传播的时间相等 , 体现在相位上即是 二者在界面上的相位相同. 我们在地面上观测的电磁场是向下传播的场和向上传播的场的 叠加的结果. 而向上传播的场是入射场经过地下各种媒质返回到地面上的场. 换句话说 , 这种场是地下各种介质源场的叠加场. 偏移成像就是把这种向上传播的场 , 以反时间方向 传播 , 依次归位到各个场源即不同介质的分界面. 在地面上观测的电磁场如何分离成向上传播的场和向下入射的场 , 这个问题在本文中 没有作过多的讨论 , 只是沿用了其他人的类似于地震勘探中的波场分离问题. 接下来的问 题就是用有限差分的方法重建这些场在地下的分布情况. 根据成像原理 , 对这些重建场进 行偏移成像 , 最后转换成电阻率成像问题. Zhdanov 等 (1996)经过多年的系统研究 , 借鉴了地震中的逆时偏移概念 , 将波场分析 的原理 , 即地震波场从地面向地下反向外推进行偏移成像的方法 , 用于对电磁场的解释 , 对大地电磁场进行逆时偏移成像 , 提出了偏移电磁场的概念 , 并在此基础上将偏移和反演3 2003202213 收到初稿 , 2004207208 收到最后修改稿 , 2004207226 决定采用. 问题结合起来. 但是其方法也存在不少问题 , 如对如何客观选取背景电导率以及电性复杂 界面的成像问题仍未能很好地解决 ; 另外 , 该方法直接用追赶法外推计算矩阵方程 , 其计 算精度不能满足高精度偏移成像要求. 根据这些方面的不足 , 在研究其方法的基础上 , 设 计了精度更高的有限差分算法. 通过迭代寻优方法选取电导率计算 , 并从理论和实际两方 面进行验证 , 取得了较原来算法更好的结果. 1 　二维电磁场有限差分偏移 对二维情况即设定 y 方向为地质体走向 (电性不变)方向 , 在介质的电阻率变化不是很 大的情况下 , 我们可以把电场 y 分量近似地表示为 ( Zhdanor et al , 1996) Ey ( x , z , ω) = QdE ( x , z , ω) exp (i kn z ) + QuE ( x , z , ω) exp ( - i kn z ) (2) 其中 , Qd , uE 是跟深度有关的待定系数 , kn ( x , z , ω) = iωσn ( x , z)μ0 , σn ( x , z) 是背景电 导率. 上述电场分量的上行、下行部分满足亥姆赫兹方程 , 即525 x2 + 525 z2 Ed , uy ( x , z , ω) + k2n ( x , z , ω) Ed , uy ( x , z , ω) = 0 (3) 将上述方程经过一系列推导得53 Qd , uE5 x2 5 z = o 2i kn 52 Qd , uE5