基于拟态有限体积法频率域可控源三维正演计算.pdf

第 59卷 第 1O期 地 球 物 理 学 报 Vo1．59。N0．10 2016年 1O月 CHINESE JOuRNAL OF GEOPHYSICS 0ct．，2016 彭荣华，胡祥云，韩波等．2016．基于拟态有限体积法的频率域可控源三维正演计算．地球物理学报，59(10)：3927—3939，doi： 10．6038／cj PengR H ，HuX Y ，HanB，eta1．2016．3D frequency-domainCSEM forwardmodelingbasedonthemimeticfinite—volume method．ChineseJ．Geophys．(inChinese)，59(10)：3927—3939，doi：10．6O38／cjg2O161O36． 基于拟态有限体积法的频率域可控源三维正演计算 彭荣华 ，胡祥云 ，韩波。，蔡建超 1中国地质大学(武汉)地球物理与空间信息学院，武汉 430074 2不列颠哥伦比亚大学地球 、海洋与大气科学学院，温哥华 V6T 1Z4，加拿大 3中国海洋大学海洋地球科学学 院，青岛 266100 摘要 大规模地球物理电磁数据的定量解释需要发展高效、稳定的三维正反演算法．本文通过求解离散化 的三维 电场矢量 Helmhohz方程 ，实现了基于有限体积法的频率域可控源 电磁 (CSEM)三维正演算法．为模拟具有强电性 差异的三维 电性介质 ，该算法采用拟态有限体积法 (MFV)对 Maxwell方程组进行离散化 ；另外，为获得稳定、高精 度 的正演数值结果 ，采用直接矩阵分解技术来求解离散所得到的大型稀疏线性方程组．对于具有多个发射源 的 CSEM测量来说，一次矩阵分解结果能够用于同频率下所有场源的正演计算．为降低场源奇异性及边界条件对数 值精度 的影响，采用虚拟场源校正技术 ，避免了散射场公式 中在构建场源项时所需的大量时间．对于具有多个频率 的CSEM 的模拟计算，采用分频并行策略来加快三维正演计算．最后，通过与一维层状模型及三维模型的数值结果 的对 比验证 了本文所开发的正演算法对频率域 CSEM 模拟计算的准确性及有效性，表 明该正演算法能够有效应用 于三维介质 的数值计算．另外，对于多频率 CSEM 的并行测试结果表 明基于分频并行策略的并行计算能够显著地 降低正演计算时间． 关键词 可控源 电磁法 ；有限体积法 ；虚拟场源校正技术 ；三维正演 ；直接分解法 doi：10．6O38／cjg2O161O36 中图分类号 P631 收稿 日期 2016—07—06，201609—12收修定稿 3D frequency—domainCSEM forwardmodelingbased on themimeticfinite-volumemethod PENG Rong—Hua ，HU Xiang—Yun卜，HAN Bo。，CAIJian—Chao 1InstituteofGeophysicsandGeomatics，ChinaUniversityofGeosciences，Wuhan430074，China 2DepartmentofEarth，OceanandAtmosphericSciences，UniversityofBritishColumbia，VancouverV6T1Z4，Canada 3CollegeofMarineGeosciences，OceanUniversityofChina，Qingdao266i00，China Abstract Quantitativeinterpretation of1arge—scalecontrolled—sourceelectromagnetic (CSEM ) data in the frequency domain requires effi