多重网格方法大地电磁正演模拟的分析.pdf

摘 要 计算大型线性方程组是大地电磁正演问题的重要部分。求解线性方程 组的一般迭代算法都可以有效光滑高频误差分量，但对低频误差分量的光 滑效率比较差，严重制约了正演模拟的效率。为了有效提高二维大地电磁 正演模拟的效率，本文从求解大型线性方程组入手，初步引入了多重网格 方法。通过分析多重网格方法在一维问题上的应用，介绍了多重网格方法 的基本原理，并证明了该方法可以高效地光滑高频误差分量和低频误差分 且 亘。 本文采用有限单元法进行二维正演模拟，求得每层网格上的刚度矩阵 和方程右端项，并根据模拟特点设计了多重网格方法的循环策略。为了证 明多重网格方法在大地电磁二维正演模拟中优越性，本文选用了三种高效 的迭代方法，即广义最小残余法、QR分解的最小二乘法和双共轭梯度稳 定法，通过比较这些算法在多重网格方法处理前后的计算结果，得到多重 网格方法在大地电磁二维正演模拟中的特点。 首先对水平层状模型进行了模拟，定量分析了多重网格方法在计算精 度上的优越性；其次分别对地下不均匀体模型和起伏地形模型进行了模 拟，定性分析了多重网格方法在迭代次数、计算速度和精度上的优越性。 关键词大地电磁，多重网格方法，正演模拟，有限单元法，迭代 ABSTRACT linear isall ofthe importantpart Solvinglarge equations forward work．For iterative modeling general algorithms，the error canbesmoothed thesmooth of components effectively,butefficiency the error are seriously low-frequencycomponentsrelativelypoor，SO，this theeffectsofforwardsimulation．Inorderto the hampered improve two-dimensionalforward efficiency,themulti-grid magnetotelluricmodeling tOthetwo-dimensionalforward methodshasbeenintroduced magnetotelluric the of methodin modeling．Byanalyzingapplicationmulti—grid one-dimensionalbasic of methodhasbeen principlesmulti-鲥d problem，the thismethodcall smoothboth itis that introduced，andproved efficiently error andlow error componentfequencycomponent． hi曲-frequency isusedint