基于SQMR方法的三维CSAMT有限差分法数值模拟

秦策 王绪本 赵宁 徐玉聪

在三维可控源音频大地电磁(CSAMT)正演中,为消除源点处场值的奇异性,将电场分为一次场和二次场分别进行计算。一次场通过解析方法进行求解,二次场的二阶偏微分方程通过Yee氏交错网格进行有限差分离散。为节省内存空间,最终形成的稀疏对称复系数线性方程组采用行压缩格式(CSR)进行存储,并使用对称拟最小残差法(SQMR)求解。最后通过几个算例验证了程序在求解方程时,计算速度快、收敛稳定的特点。