复杂条件下3d 电磁场有限元计算方法 - 科学网.pdfVIP

第 27 卷 第 6 期 地球科学 ———中国地质大学学报 Vol . 27 　No . 6 2 002 年 11 月 Earth Science —Journal of China University of Geosciences Nov . 　2 0 0 2 复杂条件下 3D 电磁场有限元计算方法 黄临平 ,戴世坤 (石油大学资源与信息学院 ,北京 102249) 摘要 : 从电磁场的 Maxwell 方程出发 ,根据电磁场的边值问题及变分公式建立了有限元方程 组. 采用可以模拟较为复杂的空间地质构造和地形起伏的四面体单元离散计算区域 ;单元中 的插值函数选择了精度较高的十点双二次多项式;并采用连续的双二次多项式插值函数来模 拟计算区域内单元中电导率 σ的空间变化. 推导出了地下变电导率 σ条件下计算三维电磁 场的有限元单元方程的解析表达式 ;采用伽辽金方法推导出了散度效正有限元方程组. 根据 所推导的公式 ,编制了三维有限元的计算程序. 数值计算结果表明 ,上述公式推导正确 ,为三 维电磁场的数值计算提供了一条有效的新途径. 关键词 : 三维电磁场 ;有限元 ; 电导率. 中图分类号: O441 　　　　文献标识码 : A ( ) 文章编号 : 1000 - 2383 2002 06 - 0775 - 05 ( ) 作者简介 : 黄临平 1964 - ,男 ,副教授 ,现在石油大学做博士后研究 ,主要研究方向为地球 物理正、反演方法理论. 　　3D 电磁法的数值模拟计算方法,通过国内外不 确 ,为三维电磁场的数值计算提供了一条有效的新 少学者多年的研究[1～4 ] ,取得了一些有意义的成果. 途径. 但由于其数值计算过于繁复并受计算机内存、计算 时间及计算量的限制 , 除了电导率分界面和某些简 1 　三维介质中的电磁场与变分原理 ( ) 单几何形体 如球体或柱体 外 ,三维模型响应的计 算较为困难. 所以 3D 电磁场的正、反演计算在实际 1. 1 　三维介质中电磁场的边值问题 资料处理中的应用仍然不多见. 本文从电磁场的 由电磁场的Maxwell 方程 ,角频率为 ω( 时间因 ω Maxwell 方程出发 ,根据电磁场的边值问题及变分公 子为 e - i t ) 的电磁场的基本方程为: 式建立有限元方程组 ,首次推导出了地下变电导率 ωμ ( ) ×E = - i H , 1 σ条件下采用四面体单元计算三维电磁场的有限元 ( σ ωε) ( ) ×H = - i E = 0 . 2 单元方程的解析表达式 ; 并利用伽辽金方法推导出 ( ) ( ) 将式 2 代入式 1 导出导电媒质中的 3D 正演 了散度效正有限元方程组. 数值计算时采用可以模 计算的 E 的微分方程为: 拟较为复杂的空间地质构造和地形起伏的四面体单 1