直流电阻率法2.5维正演的外推瀑布式多重网格法  

Extrapolation cascadic multigrid method for 2.5D direct current resistivity modeling

作　　者：潘克家[1,2,4] 汤井田[1] 胡宏伶[3,4] 陈传淼[3,4] 

机构地区：[1]中南大学有色金属成矿预测教育部重点实验室,地球科学与信息物理学院,长沙410083 [2]中南大学数学与统计学院,长沙410075 [3]湖南师范大学数学与计算机科学学院,长沙410081 [4]高性能计算与随机信息处理教育部重点实验室,长沙410081

出　　处：《地球物理学报》2012年第8期2769-2778,共10页Chinese Journal of Geophysics

基　　金：国家科技专项"深部探测技术与实验研究专项(SinoProbe-03-05)";国家自然科学基金(41204082,41174105);中国博士后科学基金(2011M501295);中南大学博士后科学基金资助

摘　　要：引入外推瀑布式多重网格法(EXCMG)求解2.5维直流电阻率有限元计算形成的大型稀疏线性方程组,结合基于地址矩阵的压缩存贮方式以及最优化离散波数,使得2.5维电阻率正演程序的计算速度大大提高而内存需求大大减小.研究结果表明:EXCMG法的收敛速度与网格尺寸无关,计算速度明显优于不完全Cholesky共轭梯度(ICCG)法.并且,随着问题规模的增大,EXCMG法的效率优势更加明显.对1600×1600网格的2.5维电阻率法模拟问题,正演程序仅耗时28s,视电阻率平均相对误差控制在0.22%以内,为进一步研究快速反演奠定了基础.Combining the compression storage mode based on address matrix with the optimized discrete wave-number,an extrapolation cascadic multigrid（EXCMG） method was introduced to solve the large sparse systems of linear equations resulting from 2.5D finite-element direct current（DC） resistivity modeling,which greatly increased the computing speed and significantly reduced the memory requirements of 2.5D forward modeling program.The research results showed that EXCMG method,which converged at a rate independent of the mesh size,was much more efficient than the incomplete Cholesky conjugate gradient（ICCG） method.Moreover,as the size of the problem increases,the efficiency advantage of EXCMG over other methods becomes more obvious.The EXCMG program for 2.5D DC resistivity modeling with a grid of 1600×1600 only took 28 seconds,and the average relative error of apparent resistivity was less than 0.22%,which laid the foundation for the latter studies of rapid inversion.

关 键 词：直流电阻率 外推瀑布式多网格法 FOURIER变换 共轭梯度法 

分 类 号：P631[天文地球—地质矿产勘探]