基于预处理共轭梯度的大地电磁快速正演 - 中南大学学报.pdf

第 41 卷第 5 期 中南大学学报( 自然科学版) Vol.41 No.5 2010 年 10 月 Journal of Central South University (Science and Technology) Oct. 2010 基于预处理共轭梯度的大地电磁快速正演 1, 2 2 2 2 张继锋 ，汤井田 ，王烨 ，肖晓 (1. 长安大学 地质工程与测绘学院，陕西 西安，710054； 2. 中南大学 信息物理工程学院，湖南 长沙，410083) 摘 要：针对大地电磁法有限元模拟中计算量大的特点，采用按行压缩存储方式的不完全 LDLT 预处理共轭梯度 T T −1 法快速求解大型复系数方程组。引入不完全 LDL 预处理，提出快速求解(LDL ) r 的方法，以加快预处理共轭梯 度法的收敛速度。研究结果表明：当网格节点自由度超过 1 万时，压缩率达到 99.9%，求解方程组时间在 1 s 以 内，为进一步快速反演奠定了基础。 关键词：按行压缩存储；预处理；共轭梯度；大地电磁 中图分类号：P631.3 文献标志码：A 文章编号：1672−7207(2010)05−1877−06 Magnetotellurics fast forward based on preconditioning conjugate gradient 1, 2 2 2 2 ZHANG Ji-feng , TANG Jing-tian , WANG Ye , XIAO Xiao (1. College of Geology Engineering and Geomatics, Chang’an University, Xi’an 710054, China; 2. School of Info-physics and Geomatics Engineering, Central South University, Changsha 410083, China) Abstract: Based on the characteristics of large computing cost in finite element modeling about magnetotelluric method, preconditioning conjugate gradient method based on compressed storage was proposed to solve large and complex T T −1 systems. Incomplete preconditioner LDL was introduced and fast solving (LDL