基于频域有限差分的三维大地电磁的低感应数预条件求解方法研究.pdfVIP

学兔兔 石油天然气学报 (江汉石油学院学报)2014年9月第36卷第9期 Journal of Oil and Gas Technology(J．JPI)Sep．2014 Vo1．36 No．9 基于频域有限差分的三维大地电磁的 低感应数预条件求解方法研究 蒋涛( 薹 。 江大 ，湖北武汉 。刚。。) 胡文宝 (油气资源与勘探技术教育部重点实验室(长江大学)，湖北武汉4301 oo) [摘要]求解Maxwell微分方程经过离散后得到的大型、病态线性方程组的高机时制约了为电磁数据反演 而进行的三维正演算法的应用。一种新的基于电场的Helmholtz分解的低感应数预条件可以解决该问题。 通过对一个三维埋藏体进行有限差分正演，分析比较了低感应数预条件迭代求解过程中2种算法 (共轭 梯度法、最小残差法)的迭代次数与时间，结果表明，采用最小残差法构建低感应数预条件矩阵，线性 方程组的迭代求解次数和和求解时间都要比共轭梯度法少，具有一定的速度优势。 [关键词]频域有限差分；正演；预条件；Helmholtz分解；低感应数 [中图分类号]P631．325 [文献标志码]A [文章编号]1000—9752(2014l 09—0088—06 正演数值计算的算法研究是地球物理勘探的基础之一，只有充分了解和认识了已知场源所产生的地 球物理异常，才可能从获得的数据中反演引起异常的异常体电磁学性质。常用的数值计算方法主要有基 于积分方程的积分方程法、边界单元法和基于微分方程的有限差分法、有限元方法。随着计算机技术的 发展和新算法技术的引进，近十年来大地电磁三维正反演取得了丰硕的成果lI～】。有限差分法因其物理 意义明确、简单直观、计算量小、计算速度快被很多学者选用 _o。谭捍东等 ]采用双共轭梯度稳定 解法求解三维交错采样有限差分离散后得到的大型系数矩阵方程组，从迭代收敛稳定、计算精度高、速 度快等方面，通过2个理论模型检验了该算法的正确性和计算精度。肖骑彬等_7 从网格精度、形成系统 方程的方式、边界条件以及预条件线性算子等，对大地电磁有限差分数值解作了对比，结果表明，合适 的预条件再配合好的线性算子，不仅可以加速收敛，而且可以降低迭代次数。李焱等 为了更有效地提 高大地电磁三维正演的计算速度，根据各频率之间求取电磁场值的过程是相互独立的特点，引入并行处 理技术，对2个理论模型的大地电磁响应采用三维交错网格有限差分数值模拟，高效地实现了正演计 算。在上述研究中，如何快速、准确地求解经过有限网格离散后的大型稀疏线性方程组成为关键。 该次研究采用Yee[9 提出的交错网格有限差分法导出了频域 Maxwell方程的离散化关系式，将三维 电磁响应的模拟问题转化成求解大型线性代数方程组的问题。笔者将总场分离成一次场 (背景场)和二 次场 (散射场)，背景场用解析方法求取，二次场用有限差分实现。导出的线性代数方程组采用预条件 方法的Krylov子空间类迭代方法求解。 为了解决频域低于 1MHz时采用频域有限差分方法得到的线性方程组在迭代求解过程中收敛缓慢 的问题，Newman等[3“ 提出了一种能显著减少三维频域有限差分数值计算时间的方法，即基于电场的 Helmholtz分解定理的低感应数预条件方法。Maxwell微分方程经过频域有限差分网格离散后得到的线 [收稿日期]2O14 o1一o6 [基金项目]国家自然科学基金项目 (4127411 5；国家重点基础研究发展计划 (“973 )项目 (2013CB228605)；油气资 源与勘探技术教育部重点实验室 (长江大学)开放基金资助项目 (K2014-23)；长江大学工程技术学院科学研究发展基