复杂地形条件下可控源音频大地电磁法2D正演数值模拟研究-地球探测与信息技术专业论文.docx

万方数据 万方数据 研究生学位论文独创性声明和版权使用授权书 独 创 性 声 明 学位论文作者（签字）：签字日期：本人声明：所呈交的论文是本人在导师指导下进行的研究工作及取得的研究 成果。据我所知，除了文中特别加以标注和致谢的地方外，论文中不包含他人已 经发表或撰写过的研究成果，也不包含为获得其它教育机构的学位或证书而使用 过的材料。对论文的完成提供过帮助的有关人员已在论文中作了明确的说明并表 学位论文作者（签字）： 签字日期： 学位论文版权使用授权书 本学位论文作者完全了解(学校)有关保留、使用学位论文的规定，有权保留 并向国家有关部门或机构送交论文的印刷本和电子版本，允许论文被查阅和借 阅。本人授权(学校)可以将学位论文的全部或部分内容编入有关数据库进行检 索，可以采用影印、缩印或扫描等复制手段保存、汇编学位论文。同时授权中国 科学技术信息研究所将本学位论文收录到《中国学位论文全文数据库》，并通过 网络向社会公众提供信息服务。（保密的学位论文在解密后适用本授权书） 学位论文作者签名：签字日期： 年 月 日导师签字：签字日期： 学位论文作者签名： 签字日期： 年 月 日 导师签字： 签字日期： 年 月 日 否 摘 要 目前，可控源音频大地电磁法的二维数值模拟正演大多数都是以矩形单元或 者三角形单元为网格剖分，并且在模拟信号场源时大都是利用一个节点上的电流 函数来表示。但是，矩形单元网格剖分在模拟复杂地形时候必然会造成精确度不 够，同时也失去一定的灵活度，三角形网格剖分单元虽然能够解决矩形这一问题， 但是三角形剖分可能会增加节点的数量，从而增加了计算时间；场源函数分布的 节点只有一个就会造成解方程组的不稳定，从而影响数值模拟结果（即电阻率值 和视相位值）的稳定性和可靠性。 为了解决以上问题，本文在推导出正演的边值问题和变分问题所满足的微分 方程基础上，选择 FORTRAN 90 语言对方程进行有限单元法编译，在有限单元法 中利用矩形中的三角形网格剖分离散求解区域和双线性插值进行模拟计算，这样 不增加节点节省了计算时间又能灵活而精确地模拟变化复杂的地形。同时，根据 模型的变化而选择不均匀网格进行剖分模拟，在异常体和场源部分采用密网格， 其他采用稀疏网格。在解决解方程组的不稳定性的问题上，则引入视狄拉克函数 对场源进行模拟，使场源被分布在多个结点上，增加了方程右端项的非零值的个 数，从而降低这种不稳定性。在验证程序可行性和精确性时，本文通过模拟计算 均匀大地模型和“k”型模型到得了相应的视电阻率曲线图和相位图曲线图，再 对比两组模型的结果图的解析解和数值解进行了验证。 最后针对具有实际生产意义的几组二维模型(纯地形山脊、纯地形山谷和纯 地形斜坡、背斜、向斜、以及含异常体的水平和含异常体起伏地形)进行正演数 值模拟分析。结果表明纯地形的山谷会在底部出现低阻高相位异常而地表空气引 起的高阻区域向深部拉伸；不含异常体的山脊地形引起底部高阻低相位异常；斜 坡模型会引起地表高阻向底部延伸；向斜模型的异常与山谷相似，并在凹陷下面 出现一处低阻低相位的异常；背斜与山脊的异常相似，在凸出部分底下出现一处 低阻高相位异常；含异常体水平地形模型的异常区域在纵向上有拉伸现象，并对 异常体的响应十分明显；山脊会凸出低阻体的异常而掩盖高阻体异常，山谷则相 反。本文通过对所设置的各个模拟结果进行了总结，这给后期的反演解译提供了 一定的指导和借鉴。 关键字：可控源音频大地电磁法，二维数值模拟，有限单元法，起伏地形 I Abstract Today, HYPERLINK /w/rectangular/ rectangular or triangle grid subdivision is widely used in Two dimensional numerical simulation of forward modeling of The Controlled Sourceacoustic Magnetotelluric (CSAMT), so is the current function HYPERLINK /w/express/ expressed the Controlled HYPERLINK /w/signal/ signal HYPERLINK /w/source/ source HYPERLINK /w/distribute/ distributed on one node . But in the matter of fact, it will cause less HYPERLINK /w/accuracy/ accuracy and loss HYPERLINK /w/flexibility/ flexibility