摘要
瞬变电磁法的三维有限元正演通常采用齐次边界条件,为满足该边界条件,需要构建较大尺寸的模型,这降低了正演问题的求解速度.针对该问题,本文采用吸收边界条件代替齐次边界条件,以缩小模型体积,加快正演速度:首先,从时间域麦克斯韦方程组出发,推导了基于库伦规范的矢量势的微分控制方程,结合一阶吸收边界条件推导了相应的的弱形式方程;在此基础上采用一阶四面体矢量单元进行单元分析、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.
作者
张永超
王光杰
李宏杰
廉玉广
李文
邱浩
牟义
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)
出处
《地球物理学报》
SCIE
EI
CAS
CSCD
北大核心
2021年第3期1106-1118,共13页
Chinese Journal of Geophysics
基金
煤炭科学技术研究院科技发展基金(2018CX06)
国家自然科学基金青年科学基金(51704162)
国家科技重大专项(2016ZX05045001-004)资助.
关键词
瞬变电磁
正演
矢量有限元
吸收边界条件
TEM
Forward modeling
Vector finite element
Absorbing boundary condition