基于非结构化网格的任意复杂2D RMT有限元模拟  

Two-dimensional complicated radio-magnetotelluric finite-element modeling using unstructured grids

作　　者：原源[1,2] 汤井田[1,2] 任政勇[1,2] 周聪[1,2] 肖晓[1,2] 张林成[1,2] 

机构地区：[1]中南大学有色金属成矿预测与地质环境监测教育部重点实验室,长沙410083 [2]中南大学地球科学与信息物理学院应用地球物理系,长沙410083

出　　处：《地球物理学报》2015年第12期4685-4695,共11页Chinese Journal of Geophysics

基　　金：国家深部探测专项第3项目(SinoProbe-03);"十二五"国家科技支撑计划课题(2011BAB04B01);国家自然科学基金(41574120);国家高技术研究发展计划(863计划)(2014AA06A602);国家重点基础研究发展计划(973计划)(2015CB060201);中南大学创新驱动计划(2015CX008)资助

摘　　要：射频大地电磁法(RMT)是以潜艇天线发射的射线源等作为场源的一种地球物理勘探方法,目前被广泛应用于近地表工程和环境地球物理勘探.RMT数据解释常采用基于静态假设的大地电磁法(MT)程序,往往会反演出不真实的浅层目标体.为解决这一长期困扰RMT资料解释的问题,本文实现了考虑位移电流的RMT有限元数值模拟.为了处理任意复杂模型,如起伏地形,非结构的三角形网格被用于离散RMT模型.首先通过算例对比,验证了程序的正确性和可靠性.通过Dike模型讨论了空气层厚度对RMT数值解的影响,结果表明当空气层厚度大于1/4波长即可满足精度要求.以山脊模型为例计算了位移电流对RMT响应的影响,表明位移电流的影响会随着山脊高程的增加而增大.最后通过舒家店实际模型进一步验证了位移电流在RMT中的重要性.As a newly developed geophysical exploration method,the radio-magnetotelluric（RMT） method is widely used in near-surface engineering and environment geophysical investigation.Since the interpretation of RMT data usually adopts the magnetotelluric（MT）forward modeling routine,in which displacement currents are generally neglected,the inverted conductivity model may not correctly reflect the true conductivity structure in the shallow subsurface.To solve this issue,we developed a finite-element forward modeling code for RMT data,in which displacement currents are considered.With this code,we analyze the dielectric effect of displacement currents on RMT responses over resistive subsurface models.First,we derived the boundary value problem（BVP）about the EM field components by theGalerkin method,in which the displacement currents were considered.Then we used the finite element method and PARDISO solver to calculate the electric and magnetic field components.To deal with complicated structure and surface topography,unstructured triangle grids were adopted for mesh generation.To improve the computation accuracy,the local refinement was used.At last,we developed a forward modeling code for RMT data with the consideration of displacement currents and analyzed the effect of displacement currents on 2DTM-mode,TE-mode data,which measured with the RMT method at frequencies between 10 and 300kHz.First,a synthetic model was used to verify the correction of our new code.The result shows that the response calculated by our code agrees well with other results.Utilizing the Dike model,the effect of the thickness of the air layer on accuracy of numerical solutions was investigated.The result shows that when the thickness of the air layer is greater than 1/4wavelength,highly accurate solutions can be guaranteed.Then the impact of displacement currents on RMT data with ridge terrain was studied on a hill model with complicated topography.From this model,the following results can be demonstrated：（1）The effect of displac

关 键 词：大地电磁 RMT 有限元 非结构网格 位移电流 

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