群表示理论在三维电磁体积分方程高效求解中的应用

• 作者: 马一行 胡俊华 刘祖鉴
• 来源期刊: 地球物理学进展
• 关键词: 群表示 Cv2群 体积分方程法 电磁阻抗矩阵
• 年，卷（期）: 2020,（5）
• 所属期刊栏目: 应用地球物理学Ⅰ(油气及金属矿产地球物理勘探)
• 页码范围: 1837-1844
• 页数: 8页
• 分类号: P631
• 语种: 中文
• DOI: 10.6038/pg2020DD0203

摘要

电磁法是矿产资源勘探和地壳结构研究的重要手段.传统电磁法理论局限于一维情况,无法处理复杂地质条件.积分方程法只需对异常体进行剖分,极大减小了内存需求,在三维数值模拟中存在较大优势.但对于大型异常体,积分方程离散涉及大型稠密电磁阻抗矩阵求逆和异常体散射电流线性方程组求解.阻抗矩阵需要计算大量的张量格林函数体积分,显著制约了积分方程法的模拟效率.群表示是对称性的一种有效数学表达.本文的目的是利用异常体对称性对散射电流的影响,减少阻抗矩阵元素的计算量,从而提高三维数值模拟效率.本文首先简要介绍了体积分方程和群表示的基本原理.然后,以Cv2群为例给出了基函数和检验函数的对称群表示.在此基础上,推导了电磁阻抗矩阵引入Cv2群的简化方程.最后,我们采用改进的积分方程法模拟了均匀半空间中具有两个垂直对称面的三维异常体的电磁响应.结果表明,Cv2的群表示可将复阻抗矩阵转化为块对角阵,使矩阵计算时间减少至1/4,总存储量和内存量分别减少至1/4和1/16,线性方程组求解时间缩短至1/12.改进算法适用于非均匀地单元离散和电导率分布.此外,本方法与观测系统的几何参数及矩量法的基函数和权函数集无关,极具灵活性,方便程序的编制和植入.总之,群表示对提高三维体积分方程法的模拟效率具有巨大的应用前景.