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

来源期刊:中南大学学报(自然科学版)2009年第2期

论文作者:柳建新 蒋鹏飞 童孝忠 徐凌华 谢维 王浩

文章页码:484 - 484

关键词:大地电磁;有限单元法;二维正演;不完全LU分解;BICGSTAB算法

Key words:magnetotelluric; finite element method; two-dimensional forward modeling; incomplete LU decomposition; BICGSTAB algorithm

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

Abstract: Based on the finite element method of quadratic interpolation in a rectangular element to solve the two-dimensional magnetotelluric forward problem, no-uniform mesh based magnetotelluric responses calculation formula was derived. Based on the fact that the finite element method was used to form a linear equation, and the coefficient matrix was a large sparse, banded, symmetric, conditioned and complex matrix, its condition number is far larger than 1, and it was an ill-conditioned matrix, solving large scale ill-conditioned linear equation was very difficult, and so the BICGSTAB (Bi-conjugate gradient stabilized method) algorithm with incomplete LU decomposition for preconditioning was used to solve this system linear equation. This approach was verified through the calculations of layered earth model and two-dimensional earth model. The results show that the BICGSTAB algorithm has high speed, high precision and stability. The apparent resistivity curves and the impedance phase curves are very accurate to synthetic magnetotelluric responses.

基金信息:国家自然科学基金资助项目
教育部博士点基金资助项目

相关论文

  • 暂无!

相关知识点

  • 暂无!

有色金属在线官网  |   会议  |   在线投稿  |   购买纸书  |   科技图书馆

中南大学出版社 技术支持 版权声明   电话:0731-88830515 88830516   传真:0731-88710482   Email:administrator@cnnmol.com

互联网出版许可证:(署)网出证(京)字第342号   京ICP备17050991号-6      京公网安备11010802042557号