基于二次场的二维大地电磁有限元法数值模拟.pdfVIP

( ) 第 35 卷第 8 期 同 济 大 学 学 报 自 然 科 学 版 Vol. 35 No. 8 2007 年 8 月 ( ) 　Aug. 2007 JOURNAL OF TON GJ I UN IV ERSITY NA TURAL SCIENCE 基于二次场的二维大地电磁有限元法数值模拟 刘小军 ,王家林 ,于　鹏 ( 同济大学 海洋地质国家重点实验室,上海　200092) 摘要 : 引入基于计算二次场响应的方法 ,将其用于二维大地电磁数值模拟. 导出了二维大地电磁二次场的表达式 , 在使用有限元法解偏微分方程时 ,利用二次场偏微分方程与总场偏微分方程的相似性 ,直接采用总场法中的系数 矩阵作为二次场的合成矩阵. 将根据此原理所编制的二次场有限元法模拟软件 ,对几个典型模型进行了试算 ,并与 传统总场模拟法作了比较分析. 结果表明 ,直接计算二次场的结果精度更高 ,更接近真实解. 关键词 : 大地电磁 ; 有限元法 ; 总场 ; 二次场 ; 数值模拟 ( ) 中图分类号: P 631 　　　　　　　文献标识码 : A 　　　　　　　　文章编号 : 0253 - 374X 2007 08 - 1113 - 05 Secondary FieldBased TwoDimensional Ma gnetotelluric Numerical Modeling by Finite Element Method L IU Xiaojun , W A N G Jialin , Y U Pen g ( State Key Laboratory of Marine Geology , Tongji Universit y , Shanghai 200092 , China) Abstract : This paper presents an efficient method named secondary finite element method ( FEM) for twodimensional magnetotelluric numerical modeling , which provides solution for secondary variations in the field. To deduce the rationale of the new procedure ,the total field is divided into primary and secondary field. Only the secondary field equations which need to be solved for the primary field can be obtained easily. The equations of the secondary field can be solved effortlessly by finite element m