地球电磁三维数值模拟多重网格方法及其应用的论文的研究

《地球电磁三维数值模拟多重网格方法及其应用的论文的研究》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、f，·’～；一}乙。’。j，．：二

2、UniversityofScienceandTechnologyofChinaAdissertationfordoctor’SdegreeStudiesofMultigridAlgorithmfor3DGeo-·ElectromagneticModelingandItsApplicationsAuthor’SName：JingjinLuSpeciality：GeophysicsSupervisor：Prof．XiaopingWuProf．KlausSpitzerFinishedtime：Novem

3、ber14m，2010中国科学技术大学学位论文原创性声明，本人声明所呈交的学位论文，是本人在导师指导下进行研究工作所取得的成果。除己特别加以标注和致谢的地方外，论文中不包含任何他人己经发表或撰写过的研究成果。与我一同工作的同志对本研究所做的贡献均已在论文中作了明确的说明。作者签名：签字日期：卫么：么2：丛中国科学技术大学学位论文授权使用声明作为申请学位的条件之一，学位论文著作权拥有者授权中国科学技术大学拥有学位论文的部分使用权，即：学校有权按有关规定向国家有关部门或机构送交论文的复印件和电子版，允许论文被查阅和借阅，可以将学位论文编

4、入《中国学位论文全文数据库》等有关数据库进行检索，可以采用影印、缩印或扫描等复制手段保存、汇编学位论文。本人提交的电子文档的内容和纸质论文的内容相一致。保密的学位论文在解密后也遵守此规定。作者签名：签字日期：摘要地球电磁勘探方法在能源及矿产资源勘探开发、水文工程地质勘察以及地壳上地幔深部研究中均发挥着重要作用。随着地质勘探目标日益复杂，三维精细结构探测成为迫切需求，导致复杂地下结构的地球物理三维数值模拟在计算效率上面临着越来越大的挑战。就目前地球电磁三维数值模拟的有限差分、有限元法而言，大规模非均匀网格、大的电性差异以及复杂的电性结

5、构等因素，使得最后形成的大型线性方程组的系数矩阵条件数很差，严重影响着现有迭代解法的收敛速度(如不完全Cholesky共轭梯度法，即ICCG法)。多重网格(MG)方法为近年来计算数学领域发展的新的高效算法，在数学上被证明是求解线性椭圆型偏微分方程的最优解法，其收敛速度与网格节点数无关。但由于诸多困难，在地球电磁三维数值模拟中的应用尚极少见。论文首先针对地球物理中常见的Poisson方程及Helmholtz方程进行了几何多重网格法数值模拟，并与目前应用于电磁三维数值模拟中较好的ICCG迭代算法进行比较。结果表明，多重网格算法收敛速度均

6、与网格节点数或网格尺度无关；而ICCG法的收敛速度则随网格节点数增加迅速减慢。几何多重网格法的计算速度比ICCG法要快数倍，证明了多重网格法三维数值模拟的高效性。然而以上几何多重网格(GMG)法应用于直流电阻率三维数值模拟遇到很大困难。点源直流电阻率法的数值模拟是一个开域的问题，三维正演中的网格高度非均匀，同时大的电性差异(不连续电性界面)、点电源的奇异性等问题都使得几何多重网格法无法获得其固有的高效性。Zhebel(2006)在其数学系博士论文中改进了传统的几何多重网格法，应用基于矩阵的多重网格法求解带电性差异的直流电阻率三维正演

7、问题，但收效甚微。代数多重网格(AMG)法近10年发展非常快，在计算流体力学领域的应用表明该方法在非均匀网格和物性系数不连续的情况下依旧能保持其固有的收敛性。论文对直流电阻率三维数值模拟的有限差分格式进行求解，针对大规模非均匀网格、大的电性差异等主要问题进行了研究，实现了电阻率三维数值模拟的代数多重网格算法，并与ICCG法进行系统的比较。大量的模型计算结果表明，AMG算法收敛速度与网格节点数、计算区域、网格精细度、电性差异及模型尺寸等因素基本无关，迭代收敛次数均在7至8次左右，其计算时间随网格节点数增加基本呈较低的线性增长：ICCG

8、法收敛曲线则呈振荡型，受上述各种因素的影响，其迭代收敛次数从数百次变化至数干次，计算时间随网格节点数增加呈非线性快速增长。一般大规模网格情况下，单个点源的三维数值摘要模拟，AMG法的计算速度约为ICCG法的3至4倍；对多个点源进行求解时，如21个点电源，AMG方法的计算速度约为ICCG法的9倍，优势将更加明显。在电阻率三维正演的代数多重网格快速算法基础上，将其应用于巷道超前探测的三维数值模拟研究中。近年来，国内煤矿屡次发生透水矿难，造成巨大的人员伤亡和财产损失，直流电阻率巷道超前探测面临着迫切的生产需求，然而由于可利用的观测空间有限

9、、巷道掘进面前方异常信号弱以及巷道空腔和旁侧异常体的影响等诸多问题，目前仍是国际上的研究难点。本文在大量地下巷道模型的电阻率三维数值模拟基础上，总结分析了巷道掘进面前方异常体真实位置与观测得到的视电阻率曲线极值点位置之间的关系，首次给