三维井间电磁场的正反演计算.pdf

第 45 卷 第 5 期 地 球 物 理 学 报 Vol . 45 , No . 5 　　　　　　 　　　　　　 　2002 年 9 月 CHINESE JOURNAL OF GEOPHYSICS 　 Sep . , 2 0 0 2 ( ) [文章编号] 　0001 - 5733 2002 05 - 0735 - 09 [ 中图分类号] 　P631 三维井间电磁场的正反演计算 魏宝君 　张庚骥 ( ( ) ) 石油大学 华东 应用物理系 , 山东东营 　257061 ( ) [摘 　要] 　提出了一种计算体积分方程的方法 ———改进型局域非线性迭代 MLNI ,并用其 对三维井间电磁场进行正反演计算. 该方法将井间大尺度散射体分为近场和远场区域两部 分 ,它们的位置和尺寸均随场点位置的变化而改变. 采用局域非线性近似计算近场区域的影 响 ,将远场区域的影响作为外部激励源 ,采用迭代方法计算. 该方法具有计算速度快 、所需内 存量少 、收敛性好的优点. 在反演中采用基于 MLNI 的非线性反演方法 ,利用这种方法可以将 成像区域集中于一定范围内而不是整个三维空间. 该反演方法由于考虑到了雅可比矩阵元 素的非线性项 ,因而与传统的Born 迭代反演方法相比更为精确和稳定. 考虑到信息量和计算 机内存的限制 ,第 1 次成像采用双重体元分割法进行较粗略的成像 ,然后缩小成像范围进行 第 2 次较精确的成像. 数值计算结果表明 ,MLNI 是一种有效的计算井间大尺度异常体散射 场的方法 ,将该方法用于反演过程能够得到较高分辨率的三维井间电导率图像. [ 关键词] 　改进型局域非线性迭代 ,三维 ,井间 , 电磁场. 1 　引　　言 利用体积分方程计算三维井间电磁场 , 由于数据点只分布在电导率异常区域而不是 整个三维空间 ,因而与有限元法和有限差分法相比数据量相对较小 ,是进行反演计算的有 效方法. 计算体积分方程的主要方法有 :逐次逼近解法 ,在异常体电导率与背景电导率对 比度较高时发散 ,适用范围很窄 ;体积分方程的直接解法[ 1] ,计算结果准确 ,但需进行直接 的大型复矩阵求逆运算 ,因而计算速度慢 ,所需内存量大 ,不适用于计算大尺度异常体的 散射场 ;系统迭代方法[2 ] ,较逐次逼近解法收敛性好 ,在井间有多个异常体时该方法较为 [3 ] 有效 ,但用于大尺度异常体在高电导率对比情况下也发散 ;局域非线性近似 和拟线性近 似[4 ] ,都是近似算法 ,无需迭代 ,具有所需内存量小 、速度快的优点 ,但在异常体尺寸较大 和高电导率对比情况下精度降低 ,在改变发射源的位置时 ,采用拟线性近似需重新计算张 [4 ] [5 ] λ 量因子 ^ ,降低了计算效率 ;改进型Born 序列 ,从理论上讲 ,无论电导率对比度有多高 均收敛 ,但由于实际计算中的误差 ,其收敛范围受到限制. 另外 ,该方法收敛速度很慢 ,在 相同模型中 ,其计算速度比