二次插值的时间域激电2.5维有限元数值模拟

摘要

通过研究基于三角单元二次插值的有限元2.5维时间域激电正演,引入Cole-Cole模型推导了2.5维复电阻率满足的边值问题,采用Guptasarma提出的滤波算法,实现了频率域激电到时域的转换;从刚度矩阵的特点出发、结合基于符号分析的线性方程组直接解法,设计实现了适合时间域2.5维激电正演的快速算法.该算法通过对计算区边界的近似处理,使处理后的刚度矩阵只与供电频率及波数有关,而与供电点位置无关,从而有效减少线性方程组直接解法中矩阵分解的工作量;设计了基于图论理论的矩阵重排与填入元分析算法,实现了高效的矩阵LDLT分解;优化了计算流程,在正演计算过程中对于同一剖分结构只需进行一次符号分析,相同频率和波数条件下所有供电点只需一次LDLT分解.最后,利用快速算法,研究了异常体的视电阻率、视极化率、视频散率等参数的异常特征.