有限单元法与有限差分法在MT一维正演模拟中的对比分析  

The comparative analysis of the finite element method and finite difference method in the MT one-dimensional forward modeling

作　　者：王涛[1,2] 柳建新[2] 童孝忠[2] 曹创华[2] 谭神湘[2] 

机构地区：[1]中国海洋大学海洋地球科学学院,青岛266100 [2]中南大学地球科学与信息物理学院,湖南长沙410083

出　　处：《物探化探计算技术》2013年第5期538-543,共6页Computing Techniques For Geophysical and Geochemical Exploration

基　　金：国家自然科学基金(41174103);教育部博士点基金(20110162130008);国家科技支撑计划项目(2011BAB04B08);中国地质调查局科研项目(资[2011]03-01-64);有色资源与地质灾害探查湖南省重点实验室项目(2010TP4012-6)

摘　　要：首先从电磁场所满足的麦克斯韦方程组出发,介绍了大地电磁测深正演的基本理论,并针对一维大地电磁模型加以讨论。运用有限单元法及有限差分法分别推导了大地电磁测深一维正演算法,并运用Matlab 7.0软件编写了相应的程序。为了检验这两种一维正演算法的准确性,设计了均匀半空间模型和层状介质模型,并给出了由本文程序、解析解得到的相应结果和图件,从而对正演结果进行对比分析。结果表明:两种方法的正演结果均真实地反映了模型的地电参数。In this paper, firstly, based on satisfaction of the Maxwell＇s equations for the electromagnetic field,the basic theory about the forward problem of magnetotelluric sounding is introduced and one-dimensional magnetotellurie modeling is discussed. The finite element method and finite difference method one-dimensional magnetotelluric forward algorithms are, then, derived. And the programs are written by using Matlab 7.0 software. In order to prove correctness of these magneto- tellurie forward algorithms, we design the homogeneous half -space model and the horizontal layered medium model and make a comparison of results which are from the forward modeling and analytic results in one-dimensional magnetotelluric. The re sults showed that： two ways of the forward results are a true refleetion of the ground electrical parameters of the model.

关 键 词：大地电磁测深 正演 数值模拟 有限单元法 有限差分法 

分 类 号：P631.325[天文地球—地质矿产勘探]