基于二次场算法的CSEM二维有限单元法正演

采用基于二次场的有限单元算法开展了二维可控场源电磁方法(CSEM)的数值模拟研究。研究工作从麦克斯韦方程组出发,导出了同时考虑电阻率与磁导率异常的可控源电磁法二次场方程;采用了基于二叉树结构的三角单元剖分,推导了双线性与双二次插值的单元钢度矩阵表达式;分别对电阻率异常及磁导率异常模型进行了试算。实践表明,二次场算法无需对场源区域剖分,可有效的减小计算区间,提高了计算效率;与总场算法相比二次场算法边界处理简单;二次场算法有利用提高计算精度,解决了低频计算不稳定问题;试算表明,较小的磁导率变化就可对可控源视电阻率产生明显影响,使电阻率变大。