二维大地电磁正演中无网格算法研究.pdfVIP

石油天然气学报 2012年5月 第34卷 第5期 JournalofOilandGasTechnology May2012 Vo1．34 No．5 二维大地 电磁正演 中无网格算法研究 苏 l，胡文宝，朱 毅 (油气资源与勘探技术教育部重点实验室 (长江大学)，湖北荆州434023) [摘要]无网格法是近几年来发展的一种新的基于变分原理的数值计算方法，由于在计算形函数中不需要 划分网格 ，在力学、电磁学等领域得到 了广泛的研究。基于无网格法在大地 电磁勘探正演 中的应用进行 了研 究，首先对无网格法 的基本原理进行 了阐述 ，并利用广义变分原理推导 出了相应 的离散方程 ，编制 了相应 的程序 。最后通过理论模型的计算结果检验 了该算法的正确性。 [关键词]无网格法 ；大地 电磁法 ；移动最小二乘拟合 (MLS)；有限元法 [中图分类号]P631．325 [文献标识码]A [文章编号]1000—9752(2012)05—0087一O4 大地电磁测深法 (MT)是一种以岩石 的电性差异为基础和前提 ，利用天然交变电磁场研究地球电 性结构的勘探方法_1]。在解决 MT正演问题时 ，有限差分法 和有限元法 捌是主要的数值计算方法； 尽管这些方法取得了巨大成功 ，但是这些方法都是基于网格的数值计算方法 ，也有其 自身 的局限性_7]， 如 网格部分计算成本高、计算精度依赖于单元剖分的形状和大小等。产生上述 问题 的一个主要原 因是这 些方法都是基于网格剖分来建立离散系统方程的，为克服这些困难 ，无网格方法 (Meshless)作为有限 元法 (FEM)的一种重要补充具有重要的研究意义。 无网格法 (Meshless)是近 1O年来兴起的一种数值计算方法口]，它是在有限元法的基础上发展起 来的，但是不同于有限元 ；无网格法 的近似函数是建立在一系 B 列离散节点上的，然后在这些节点上利用形 函数求 出其近似 空气 值 ，代人相应的变分问题得到系统方程 。 地下／／／1J 无 网格法 已在力学 ]、油藏渗流 问题lg]、电磁学[1。。等领 y l 域得到了广泛的研究 ，然而在求解地球物理正演问题 中的理论 z 和应用研究鲜少 出现 。为此 ，笔者 利用无 网格法对二 维 MT 正演做 了相关研究 。首先从麦克斯韦方程组出发 ，推导出二维 介质中边值所对应 的广义变分 问题 ；然后用节点离散求解 区 Q 域 ，得到建立在离散节点上的无网格形 函数 ，求解变分问题得 C D 到边值问题对应的离散方程 ，解方程后代入视 电阻率计算公式 图 1 理论模型 计算地面上待求点的视电阻率 。 1 边值问题及其变分问题 大地 电磁场可以看作是从高空垂直入射到地面的平面电磁波Ⅲ ，基于这种假设可以推导二维情况下 大地 电磁测深法所对应的边值 问题[2]，如图 1所示 ： {f ·(r “)+Au=0E Ul一“ r 一r。 ∈r (1) l“一1EAB 二==0EAc，BD +k一0ECD ． tJ H o 7t 式中， 是哈密顿算子 ；为图1所示 的求解区域ABCD；为法线方向；r为内边界；对于电场极化 (TE极 [收稿 日期]2012—02—10 [基金项目]国家 “973 规划项目 (20o7cB2o9607)；国家自然科学基金项 目 。 [作者简介]苏洲 (