不完全LU分解预处理的BICGSTAB算法在大地电磁二维正演模拟中的应用

• 作者: 徐凌华 柳建新 王浩 童孝忠 蒋鹏飞 谢维
• 来源期刊: 中南大学学报（自然科学版）
• 学科: 地球科学
• 关键词: 大地电磁 有限单元法 二维正演 不完全LU分解 BICGSTAB算法
• 年，卷（期）: 2009,（2）
• 所属期刊栏目: 地质工程·土本工程
• 页码范围: 484-491
• 页数: 8页
• 分类号: P631.3
• 字数: 3629字
• 语种: 中文

摘要

基于双二次插值的有限单元法求解大地电磁二维正演问题,以不均匀网格剖分为基础,推导出大地电磁响应的计算公式.针对有限单元法最后形成一个线性方程组,系数矩阵是大型稀疏的带状对称正定复系数矩阵,并且其条件数远大于1,为严重病态矩阵,求解其对应方程组会遇到很多困难等问题,采用不完全LU分解(即上三角与下三角分解)处理的稳定双共轭梯度算法(BICGSTAB算法)求解该线性方程组,通过对层状介质和二维模型电磁响应进行计算,获得二维大地电磁的视电阻率曲线和阻抗相位曲线.研究结果表明,BICGSTAB算法具有速度快、精度高和稳定性好等优点.

作者信息

序号	姓名	单位	发文数	被引次数	H指数	G指数
1	柳建新	中南大学信息物理工程学院	157	1584	20.0	32.0
2	蒋鹏飞	中南大学信息物理工程学院	3	58	1.0	3.0
3	童孝忠	中南大学信息物理工程学院	48	391	11.0	17.0
4	徐凌华	中南大学信息物理工程学院	2	73	2.0	2.0
5	谢维	中南大学信息物理工程学院	25	285	8.0	16.0
6	王浩	中南大学信息物理工程学院	66	553	13.0	22.0