复杂地形电导率线性变化二维大地电磁有限单元法正演模拟-地球探测与信息技术专业论文.docx

万方数据 万方数据 研究生学位论文独创性声明和版权使用授权书独 研究生学位论文独创性声明和版权使用授权书 独 创 性 声 明 本人声明：所呈交的论文是本人在导师指导下进行的研究工作及取得的研究成 果。据我所知，除了文中特别加以标注和致谢的地方外，论文中不包含他人已经发 表或撰写过的研究成果，也不包含为获得其它教育机构的学位或证书而使用过的材 料。对论文的完成提供过帮助的有关人员已在论文中作了明确的说明并表示谢意。 学位论文作者（签字）： 签字日期： 学位论文版权使用授权书 本学位论文作者完全了解(学校)有关保留、使用学位论文的规定，有权保留并 向国家有关部门或机构送交论文的印刷本和电子版本，允许论文被查阅和借阅。本 人授权(学校)可以将学位论文的全部或部分内容编入有关数据库进行检索，可以采 用影印、缩印或扫描等复制手段保存、汇编学位论文。同时授权中国科学技术信息 研究所将本学位论文收录到《中国学位论文全文数据库》，并通过网络向社会公众 提供信息服务。（保密的学位论文在解密后适用本授权书） 学位论文作者签名： 导师签字： 签字日期： 年 月 日 签字日期： 年 月 日 摘要 大地电磁测深法是基于岩石的电性差异为前提，利用天然交变电磁场作为场源 来研究地质体电性结构的勘探方法，它属于电磁法勘探中的频率域方法。目前，已 有大量的文献对二维大地电磁法正演方面进行了广泛的研究,比如用各种数值模拟 方法对二维大地电磁进行正演模拟。而有限差分法，有限单元法和边界单元法是比 较常用数值模拟方法。虽然这三种常规的方法在二维大地电磁正演研究方面上得到 较好的研究成果，但是存在一些不足之处，正如复杂地形模拟、剖分方式选取、电 性参数设定以及辅助场计算等。本文是在电导率线性变化情况下，推导出复杂地形 条件下二维大地电磁法有限单元法正演模拟，以提高正演模拟精度和效率为目的， 为今后反演的实现和资料解释提供更加充分的正演支持。 本文首先综合阐述了大地电磁法数值模拟方法的研究历史及其研究现状，然后 从麦克斯韦方程组出发，建立了大地电磁测深法边值问题和变分问题。利用三角单 元剖分模拟复杂地形及地质结构，而电导率在单元内线性变化，更加符合野外地质 条件，易于形成等值线图，同时给出了具有垂直断面地形的简单网格剖分方法。在 三角单元剖分电导率线性变化条件下，利用 Mathematica 计算出大地电磁 TE（横电） 模式有限单元法系数矩阵，同时，在 TM（横磁）模式下，利用高斯积分处理后计算 出有限单元法系数矩阵。在求取辅助场时，采用二次朗格朗日插值计算场值的偏导 数直接获得辅助场。最后利用 Fortran90 编写出复杂地形条件电导率线性变化的大 地电磁二维正程序，并详细介绍了该程序的实现过程。 通过一些模型来实现了 TE、TM 模式复杂地形电导率线性变化大地电磁二维有限 单元法正演模拟。首先设计相应的一维地电模型做正演模拟，将模拟结果与理论结 果进行比较，以其来检验正确性和精确性；然后用典型的二维地电模型做正演模拟， 将模拟结果与前人结果比较，以其来检验可靠性和准确性；最后设计复杂地电模型 （具有垂直间断面地形）做正演计算，对结果进行分析，并解释该地形地质条件下 大地电磁场的分布特征，同时检验本文方法简单适用性。 关键词：二维大地电磁，电导率线性变化，有限单元法，垂直间断面，正演模拟 I 万方数据 万方数据 Abstract Based on the difference between the rock for the foundation and precondition,the Magnetotelluric use the natural alternating electromagnetic field as a source to study geological exploration method of electrical structure. It is the electromagnetic exploration methods which belong to frequency domain. At present,Magnetotelluric method of two-dimensional forward modeling problem has been widely research in the number of literature.Such as all kinds of numerical simulation method is used to solve the problem of two-dimensional Magnetotelluric forward.And then,finite dif