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

• 作者: 周峰 张志勇 李泽林
• 来源期刊: 石油地球物理勘探
• 关键词: 激发极化法 稀疏矩阵 有限元 二次插值
• 年，卷（期）: 2015,（5）
• 所属期刊栏目: 非地震
• 页码范围: 999-1006
• 分类号: P631
• 语种: 中文
• DOI: 10.13810/j.cnki.issn.1000-7210.2015.05.025

摘要

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

作者信息

序号	姓名	单位	发文数	被引次数	H指数	G指数
1	张志勇	中南大学地球科学与信息物理学院	13	79	5.0	8.0
5	周峰	东华理工大学核工程与地球物理学院	12	16	2.0	3.0