基于第三类边界条件的大地电磁三维矢量有限元任意各向异性正演  

3D vector finite element modeling of magnetotelluric in general anisotropic media with the third type boundary condition

作　　者：余益敏 黄松[1,2] 肖调杰 徐亚 YU YiMin;HUANG Song;XIAO TiaoJie;XU Ya(Key Laboratory of Petroleum Resource Research,Institute of Geology and Geophysics,Chinese Academy of Sciences,Beijing 100029,China;Innovation Academy for Earth Science,Chinese Academy of Sciences,Beijing 100029,China;University of Chinese Academy of Sciences,Beijing 100049,China;School of Computer,National University of Defense Technology,Changsha 410073,China)

机构地区：[1]中国科学院地质与地球物理研究所,中国科学院油气资源研究重点实验室,北京100029 [2]中国科学院地球科学研究院,北京100029 [3]中国科学院大学,北京100049 [4]国防科技大学计算机学院,长沙410073

出　　处：《地球物理学报》2023年第5期2181-2200,共20页Chinese Journal of Geophysics

基　　金：国家自然科学基金项目(91858214和42104078)资助。

摘　　要：地球内部介质的各向异性对地球物理场解译有很大影响,研究各向异性介质中大地电磁响应具有重要的意义.边界条件是影响电磁场正演精度的一个关键因素,其中第一类边界条件需要将底面边界设置在离异常体足够远的地方,面临着计算规模大、求解速度慢的问题.相比第一类边界条件,第三类边界条件具有计算规模更小、计算精度更高的优点,在三维各向同性正演中被广泛使用.然而,各向异性使得第三类边界条件理论变得更为复杂,目前尚未见到采用第三类边界条件的大地电磁三维各向异性正演.本文推导出各向异性介质中三维矢量有限元底界面的第三类边界条件,通过模型计算对算法的可靠性和精度进行了验证.结果表明,当精度相同时,第三类边界条件的计算规模可以更小;而当底界面离异常体较近时,第三类边界条件有着更高的精度.The anisotropy of the earth's internal medium has a great influence on the geophysical exploration.Therefore,it is important to study the responses of magnetotelluric(MT)method in anisotropic media.The boundary conditions play an important role in the accuracy of electromagnetic modeling.The first type boundary condition requires the bottom boundary to be set far enough away from the anomalies,which leads to a larger computational domain and longer computational time.The third type boundary condition has smaller computational scale and higher computational accuracy than the first type boundary condition,which is widely used in the modeling of 3D isotropic media.However,the theory of the third type boundary condition is more complex in anisotropic media.At present,there are no published papers that achieved 3D MT modeling with the third type boundary condition in general anisotropic media.In this paper,for 3D vector finite element the third type boundary condition of the bottom interface in anisotropic media is derived.The reliability and accuracy of the algorithm are verified by several models.The results show that the computational domain with the third type boundary condition can be smaller when the accuracy is the same;and the third type boundary condition has a higher accuracy when the bottom boundary is closer to the anomaly.

关 键 词：大地电磁 矢量有限元法 各向异性 第三类边界条件 

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