张永超^1,2,3,,
王光杰^4,5,,,
李宏杰^1,2,3,
廉玉广^1,2,3,
李文^1,2,3,
邱浩^1,2,3,
牟义^1,2,3
1. 煤炭科学技术研究院有限公司安全分院, 北京 100013
2. 煤炭资源高效开采与洁净利用国家重点实验室(煤炭科学研究总院), 北京 100013
3. 北京市煤矿安全工程技术研究中心, 北京 100013
4. 中国科学院大学, 北京 100049
5. 中国科学院地质与地球物理研究所 页岩气与地质工程重点实验室, 北京 100029

基金项目: 煤炭科学技术研究院科技发展基金（2018CX06），国家自然科学基金青年科学基金（51704162），国家科技重大专项（2016ZX05045001-004）资助


详细信息

作者简介: 张永超, 男, 1983年生, 助理研究员, 主要从事瞬变电磁法正反演研究及应用. E-mail: 120460817@qq.com

通讯作者: 王光杰, 男, 1966年生, 中国科学院地质与地球物理研究所副研究员, 主要从事电磁法勘探研究和应用工作. E-mail: gjwang@mail.iggcas.ac.cn

中图分类号: P631

收稿日期:2019-04-04
修回日期:2021-01-27
上线日期:2021-03-10

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

ZHANG YongChao^1,2,3,,
WANG GuangJie^4,5,,,
LI HongJie^1,2,3,
LIAN YuGuang^1,2,3,
LI Wen^1,2,3,
QIU Hao^1,2,3,
MOU Yi^1,2,3
1. Mine Safety Technology Branch of China Coal Research Institute, Beijing 100013, China
2. State Key Laboratory of Coal Mining and Clean Utilization(China Coal Research Institute), Beijing 100013, China
3. Beijing Coal Mine Safety Engineering Technology Research Center, Beijing 100013, China
4. University of Chinese Academy of Sciences, Beijing 100049, China
5. Key Laboratory of Shale Gas and Geoengineering, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China


More Information

Corresponding author: WANG GuangJie,E-mail:gjwang@mail.iggcas.ac.cn

MSC: P631

--> Received Date: 04 April 2019
Revised Date: 27 January 2021
Available Online: 10 March 2021


摘要
摘要:瞬变电磁法的三维有限元正演通常采用齐次边界条件，为满足该边界条件，需要构建较大尺寸的模型，这降低了正演问题的求解速度.针对该问题，本文采用吸收边界条件代替齐次边界条件，以缩小模型体积，加快正演速度：首先，从时间域麦克斯韦方程组出发，推导了基于库伦规范的矢量势的微分控制方程，结合一阶吸收边界条件推导了相应的的弱形式方程；在此基础上采用一阶四面体矢量单元进行单元分析、Newmark法进行时间离散，实现了瞬变电磁法的快速三维正演.通过均匀半空间模型的解析解，H型地电断面的CR1Dmod解和相应模型有限元解的对比，验证了本文算法的正确性.均匀半空间模型分别采用吸收边界条件和齐次边界条件的正演结果对比表明：吸收边界条件确实可以提高三维正演的精度或者缩小模型尺寸、加快计算速度.
关键词: 瞬变电磁/
正演/
矢量有限元/
吸收边界条件

Abstract: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.
Key words:TEM/
Forward modeling/
Vector finite element/
Absorbing boundary condition



PDF全文下载地址:

http://www.geophy.cn/data/article/export-pdf?id=608009f7c23e67607ae103d6