二维有限元大地电磁正演模拟在Matlab上的实现.pdfVIP

地 震 地 磁 观 测 与 研 究 第 29卷 第5期 SEISM OLOGICALA D GEOMAGNETIC VoI．29 NO．5 2008年 10月 OBSERVATION AND RESEARCH Oct 2008 二维有限元大地 电磁正演模拟 在 Matlab上的实现 张 昆¨ 魏文博 ¨ 叶高峰 1 中国北京 100083中国地质大学地球物理与信息技术学院 2)中国北京 100083地下信息探测技术与仪器教育部重点实验室(中国地质大学 ．北京) 摘要 Matlab高级语言开发平 台对矩阵的计算简便快速，并且绘 图输出文件方便。以这些特点为 契机 ，在 Matlab平 台上 ，讨论利用有限单元法实现二维大地 电磁正演的问题 ，推导出相应 的公式 ， 计算地表或海底视电阻率及相位 ，并与解析解 以及主流方法结果对 比，证实 了本文方法的正确性 。 关键词 Matlab；二维正演 ；MT；有限元 中图分类号 ：P315．69 文献标识码 ：A 文章编号 ：1003—3246(2008)05—0083—06 引言 大地 电磁测深法是一种重要的地球物理勘探方法 ，具有勘探深度大 、不受高阻屏蔽 、对低 阻有较高分辨力等优点，因此广泛应用于地质构造研究以及资源勘查等方面。 当前 ，国内外很多学者对大地电磁测深法作了大量的研究工作，】971年 Coggon提出利用 有限元数值模拟方法计算大地电磁场。他详细阐述 了怎样用有限元法模拟大地电磁场 ，从能 量最小原理推导出亥姆霍兹方程的有限元格式。在 国内，很多学者曾在大地电磁的二维有限 元正演方面做过工作，如陈乐寿，周熙襄等。 目前 ，常规正演方法都是采用 Fortran或者 C语言为开发平台，用这样的开发平台可以实 现快速计算 ，在很短时间内得到正演模型。但是基于类似中级或低级语言开发程序代码 ，格式 和语法要求严格，代码 比较繁琐 ，例如对矩阵的计算就十分麻烦 ，需要编写很长的代码实现矩 阵计算 ，一段庞大的代码往往只能实现一项功能 ，不方便使用。而且常规方法的模型剖分数据 文件 比较复杂 ，不利于修改。 为克服以上不利因素，我们以现今发达的PC硬件及软件技术为基础，在 Matlab(矩阵工 作室)语言环境下实现了二维大地电磁正演。 1 大地电磁场 的二维问题 在进行海洋大地电磁场数值模拟研究时，我们首先从描述电磁场分布的麦克斯韦方程组 作者简介 ：张昆(1983 )，男，天津塘沽人 ，现于中国地质大 学 (北京)攻读地球物理学硕士学位 ，主要研 究方 向为地 下深部构造及大地 电磁正、反演方法研究 基金项 目：国家 自然科学基金重点项 目“华北地区现今岩石圈导 电性结构及其大陆动力学意义” 本文收到 日期 ：2008—04—23 地 震 地 磁 观 测 与 研 究 出发，推导出亥姆霍兹方程；然后利用二维有限元方法得到相应的 目标泛函，将解偏微分方程 问题转换为解泛函变分问题 ；再利用矩形单元剖分的有限元方法实现正演模拟。最后编写 出 正演代码并验证结果正确性。 1．1 方法原理 忽略位移电流 ，描述电磁场分布的麦克斯韦方程组是 ×E = “H X H = · H ： 0 (1) · D ： 0 为磁导率 ，为电导率。 为满足能量最小原理，利用不同模式的总能量表达式 ，推导出在二维边界给定条件下的泛 函表达式 ，将问题转化为解变分的问题 。 K — f (2) 建立数学模型，使用矩形单元剖