基于并行化直接解法的频率域可控源电磁三维正演  

3-D frequency-domain CSEM modeling using aparallel direct solver

作　　者：韩波[1] 胡祥云[1] 黄一凡[1] 彭荣华[1] 李建慧[1] 蔡建超[1] 

机构地区：[1]中国地质大学(武汉)地球物理与空间信息学院,地球内部多尺度成像湖北省重点实验室,武汉430074

出　　处：《地球物理学报》2015年第8期2812-2826,共15页Chinese Journal of Geophysics

基　　金：国家自然科学基金项目(41274077;41474055);中国地质调查局项目(12120113101800)联合资助

摘　　要：电磁法的三维数值模拟是一个对数值算法和计算机硬件要求都非常高的问题.对常用的微分类方法如有限单元法和有限差分法而言,求解最后所得的大型线性方程组是至关重要的一步,直接影响到正演算法的实用性.如何高效、稳定且准确地解线性方程长期以来一直是被探讨的问题.本文实现了基于线性系统直接求解技术的频率域可控源电磁(CSEM)三维正演.使用交错网格有限体积法(FV)来离散化关于二次电场的Helmholtz方程;使用直接解法取代传统的迭代解法来求解离散线性系统,即对系统矩阵进行完全LU分解,具体通过调用大规模并行矩阵直接求解器(MUMPS)来实现.基于理论模型做了一系列数值实验,首先证明了直接解法的高精度和稳定性,并考察了其内存需求、计算时间和并行可伸缩性等主要计算性能,最后检验了所开发的算法快速模拟多场源CSEM问题的能力以及对常规海洋和陆地CSEM模拟的有效性.Three-dimensional modeling of electromagnetic data is a computationally demanding problem.For frequently-used numerical techniques such as finite-element and finite-difference methods,solving the large linear systems arising from the discretization of Maxwell＇s equations is a key step which has a major impact on the applicability of the solution,and it has always been a research topic to solve the linear equations efficiently,robustly and accurately.A 3D modeling scheme based on direct solutions of the linear system is presented for frequency domain controlled-source electromagnetic（CSEM）surveys. The Helmholtz equation in terms of secondary electric fields is discretized using a finite-volume（FV）method over a staggered grid.Taking advantage of recent developments in numerical algorithms and the availability of computational resources,the resulting linear system of FV equations is solved directly using the massively parallel solver,namely MUMPS,instead of the most commonly used linear solvers,i.e.Krylov subspace iterative techniques.The direct solver carries out an LU（and possibly LDLT）decomposition of the system matrix and then computes solutions efficiently by applying forward and backward substitutions.To evaluate the computational performance of the direct solver,a series of numerical tests based on synthetic 1D models were conducted,and the results indicate that（1）Normalized residuals of solutions are almost independent of the conductivity value assigned to air layers but increase rapidly as the frequency value decreases.Nevertheless,the order of magnitude of the largest normalized residual is as small as 10-11.At the same time,although the matrix factorization time varies as either the air conductivity or the frequency changes,the variation is only a fraction of the total run time.（2）Both the execution time and required memory increase rapidly（more than linearly）with increasing grid sizes.（3）By executing MUMPS in parallel over multiple processors,not only the total run time but

关 键 词：可控源电磁法 三维模拟 LU分解 直接解法 

分 类 号：P631[天文地球—地质矿产勘探]