文献类型：期刊文章

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

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

出　　处：《地球物理学报》

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

年　　份：2015

卷　　号：58

期　　号：8

起止页码：2812-2826

语　　种：中文

收录情况：AJ、BDHX、BDHX2014、CAS、CSA、CSA-PROQEUST、CSCD、CSCD2015_2016、EI、GEOBASE、GEOREFPREVIEWDATABASE、IC、INSPEC、JST、PA、RCCSE、SCIE、SCOPUS、WOS、ZGKJHX、核心刊

摘　　要：电磁法的三维数值模拟是一个对数值算法和计算机硬件要求都非常高的问题.对常用的微分类方法如有限单元法和有限差分法而言,求解最后所得的大型线性方程组是至关重要的一步,直接影响到正演算法的实用性.如何高效、稳定且准确地解线性方程长期以来一直是被探讨的问题.本文实现了基于线性系统直接求解技术的频率域可控源电磁(CSEM)三维正演.使用交错网格有限体积法(FV)来离散化关于二次电场的Helmholtz方程;使用直接解法取代传统的迭代解法来求解离散线性系统,即对系统矩阵进行完全LU分解,具体通过调用大规模并行矩阵直接求解器(MUMPS)来实现.基于理论模型做了一系列数值实验,首先证明了直接解法的高精度和稳定性,并考察了其内存需求、计算时间和并行可伸缩性等主要计算性能,最后检验了所开发的算法快速模拟多场源CSEM问题的能力以及对常规海洋和陆地CSEM模拟的有效性.

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

分 类 号：P631]