地震地质 ›› 2008, Vol. 30 ›› Issue (2): 525-533.

• 研究论文 • 上一篇    下一篇

大地电磁测深二维正演中辅助场的新算法

马为, 陈小斌, 赵国泽   

  1. 中国地震局地质研究所, 北京, 100029
  • 收稿日期:2007-06-04 修回日期:2007-08-05 出版日期:2008-06-20 发布日期:2009-08-16
  • 作者简介:马为,男,1981年生,2004年毕业于吉林大学,2007年毕业于中国地震局地质研究所,获硕士学位,主要从事电磁测深方法研究与应用,电话:010-62009067,E-mail:xboy149@sohu.com.
  • 基金资助:
    国家自然科学基金(40674036);国家自然科学基金重点项目(40534023)共同资助

A NEW ALGORITHM FOR THE CALCULATION OF AUXILIARY FIELD IN MT 2D FORWARD MODELING

MA Wei, CHEN Xiao-bin, ZHAO Guo-ze   

  1. Institute of Geology, China Earthquake Administration, Beijing 100029, China
  • Received:2007-06-04 Revised:2007-08-05 Online:2008-06-20 Published:2009-08-16

摘要: 文中提出一种计算大地电磁二维正演辅助场的新方法。在主场满足线性插值基函数的前提下,通过构建二次插值基函数来提高辅助场计算中数值求导的精度。通过与线性插值函数法的计算结果进行对比,证明该算法在基本不增加计算量的同时可以大幅度提高大地电磁二维正演精度。文中还研究了网格对辅助场计算的影响,结果表明二次插值法较线性插值法对网格的要求更为宽松。此外还以均匀半空间模型为例讨论了地表网格纵向间距与模型之间的关系,这对正演计算中的网格设计有参考意义。

关键词: 大地电磁, 辅助场, 二维正演, 有限元直接迭代算法, 二次插值, 网格

Abstract: This paper focuses on the calculation of the auxiliary field-another important factor related the precision of MT 2D forward modeling.The purpose is to raise the precision of forward calculation and provide MT inversion with more sufficient supports.It suggests constructing a quadratic interpolation shape function in the surface grid for calculation of the auxiliary field,when employing the current available method of linear interpolation for calculating the principal field.In this work,the Direct Iteration Finite Element Algorithm(DIFE)is utilized to solve the calculations of interpolation nodes.This new algorithm is tested through model calculation.Comparing with the 2D forward modeling of the auxiliary field by the linear interpolation,the quadratic interpolation does not increase much amount of computation,but raises the precision of MD forward calculation.We also discussed the reason from theory,and on the basis of the result,the current design criterion for 2D forward grid has been improved.The element size can be enlarged to one skin effect depth.So the new algorithm has special meaning to the study of complex topography forward problem.We will further discuss this question and compare the results.Besides,the calculation of auxiliary field in 3D forward modeling is also very interesting.

Key words: magnetotelluric, auxiliary field, 2D forward modeling, DIFE, quadratic interpolation, grid

中图分类号: