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

作者：马一行,胡俊华,刘祖鉴

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

发文机构：中国地质调查局发展研究中心 自然资源部矿产勘查技术指导中心 中国科学院地质与地球物理研究所 中国自然资源航空物探遥感中心

关键词：群表示C2v群体积分方程法电磁阻抗矩阵Group representationC2v groupVolume integral equation methodElectromagnetic impedance matrix

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