科研之友
微信
新浪微博
Facebook
分享链接
摘要
快速的电磁法正演模拟是反演的核心要素之一.本文针对陆地可控源电磁(CSEM)问题，实现了三维多频的快速正演.首先基于电场的微分控制方程，利用非结构化网格和伽辽金方法建立关于电场矢量的有限元线性方程组.为快速求解大规模复系数有限元线性方程组，本文引入了基于有理Krylov子空间的单极点模型降阶算法.该方法只需求解实系数矩阵方程，然后使用Lanczos算法快速构建有理Krylov正交基子空间，达到系数矩阵维数降阶的目的，从而加快正演求解速度.层状模型和低阻块体模型的计算模拟表明：在精度满足条件的情况下，同常规有限元方法相比，加速比达到10～20多倍.对于百万未知数级、几十至上百个频率的三维可控源电磁正演问题，在个人PC机上只需几分钟便可快速求解.
- 单位
