大地电磁测深二维正演中辅助场的新算法-地震地质.pdf

第 ３０卷　第２期 地　震　地　质 Ｖｏｌ．３０，Ｎｏ．２ ２００８年６月 ＳＥＩＳＭＯＬＯＧＹＡＮＤＧＥＯＬＯＧＹ Ｊｕｎｅ，２００８ 大地电磁测深二维正演中辅助场的新算法 马　为　陈小斌　赵国泽 （中国地震局地质研究所，北京　１０００２９） 摘　要　　文中提出一种计算大地电磁二维正演辅助场的新方法。在主场满足线性插值基函数的前 提下，通过构建二次插值基函数来提高辅助场计算中数值求导的精度。通过与线性插值函数法的计 算结果进行对比，证明该算法在基本不增加计算量的同时可以大幅度提高大地电磁二维正演精度。 文中还研究了网格对辅助场计算的影响，结果表明二次插值法较线性插值法对网格的要求更为宽松。 此外还以均匀半空间模型为例讨论了地表网格纵向间距与模型之间的关系，这对正演计算中的网格 设计有参考意义。 关键词　　大地电磁　辅助场　二维正演　有限元直接迭代算法　二次插值　网格 ＋ 中图分类号：Ｐ