基于吸收边界条件的瞬变电磁法三维矢量有限元快速正演  

Efficient 3D vector finite element modeling for TEM based on absorbing boundary condition

作　　者：张永超 王光杰[4,5] 李宏杰 廉玉广[1,2,3] 李文 邱浩[1,2,3] 牟义 ZHANG YongChao;WANG GuangJie;LI HongJie;LIAN YuGuang;LI Wen;QIU Hao;MOU Yi(Mine Safety Technology Branch of China Coal Research Institute,Beijing 100013,China;State Key Laboratory of Coal Mining and Clean Utilization(China Coal Research Institute),Beijing 100013,China;Beijing Coal Mine Safety Engineering Technology Research Center,Beijing 100013,China;University of Chinese Academy of Sciences,Beijing 100049,China;Key Laboratory of Shale Gas and Geoengineering,Institute of Geology and Geophysics,Chinese Academy of Sciences,Beijing 100029,China)

机构地区：[1]煤炭科学技术研究院有限公司安全分院,北京100013 [2]煤炭资源高效开采与洁净利用国家重点实验室(煤炭科学研究总院),北京100013 [3]北京市煤矿安全工程技术研究中心,北京100013 [4]中国科学院大学,北京100049 [5]中国科学院地质与地球物理研究所页岩气与地质工程重点实验室,北京100029

出　　处：《地球物理学报》2021年第3期1106-1118,共13页Chinese Journal of Geophysics

基　　金：煤炭科学技术研究院科技发展基金(2018CX06);国家自然科学基金青年科学基金(51704162);国家科技重大专项(2016ZX05045001-004)资助.

摘　　要：瞬变电磁法的三维有限元正演通常采用齐次边界条件,为满足该边界条件,需要构建较大尺寸的模型,这降低了正演问题的求解速度.针对该问题,本文采用吸收边界条件代替齐次边界条件,以缩小模型体积,加快正演速度:首先,从时间域麦克斯韦方程组出发,推导了基于库伦规范的矢量势的微分控制方程,结合一阶吸收边界条件推导了相应的的弱形式方程;在此基础上采用一阶四面体矢量单元进行单元分析、Newmark法进行时间离散,实现了瞬变电磁法的快速三维正演.通过均匀半空间模型的解析解,H型地电断面的CR1Dmod解和相应模型有限元解的对比,验证了本文算法的正确性.均匀半空间模型分别采用吸收边界条件和齐次边界条件的正演结果对比表明:吸收边界条件确实可以提高三维正演的精度或者缩小模型尺寸、加快计算速度.The homogeneous boundary condition is commonly used in the 3D FEM forward modeling of TEM.In order to satisfy the approximation of this boundary condition,it is necessary to construct a large-scale model,which reduces the solving speed of the forward problem.Aiming at the problem,in this paper,absorbing boundary condition was used to reduce the size of the model and accelerate the solution speed instead of homogeneous boundary condition.Firstly,the differential equation of magnetic vector potential based on Coulomb′s gauge was derived from Maxwell′s equations in time domain,combining the first-order absorbing boundary condition,the corresponding weak form equation was deduced and simplified.Then,by using the first-order tetrahedral vector element to carry out element analysis,Newmark method for time discretization,an efficient TEM 3D modeling method was realized.The proposed algorithm was validated by comparing the analytical solution of homogeneous half space,the CR1Dmod solution of H-type geoelectric section,with their corresponding finite element solution.The further comparison between the modeling results of absorbing boundary condition and homogeneous boundary condition showed that the absorbing boundary condition could improve the accuracy of 3D modeling or reduce the size of the model and accelerate the solution speed.

关 键 词：瞬变电磁 正演 矢量有限元 吸收边界条件 

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