基于非均匀网格有限差分法的大地电磁静位移模拟
Modeling of Magnetotelluric Static Shift Based on Non-Uniform Meshes Finite Difference Method
DOI: 10.12677/AG.2018.88148, PDF,    科研立项经费支持
作者: 邓小康:中南勘测设计研究院水能资源利用关键技术湖南省重点实验室,湖南 长沙;湖南工业大学土木工程学院,湖南 株洲;童孝忠*:中南大学地球科学与信息物理学院,湖南 长沙
关键词: 大地电磁二维模型静位移有限差分法非均匀网格Magnetotelluric Two-Dimensional Model Static Shift Finite Difference Method Non-Uniform Meshes
摘要: 在大地电磁勘探中,视电阻率曲线总会存在一些静态现象,而静态问题又一直是电磁测深领域一个棘手的问题,它直接影响到电磁测量结果的解释,能否正确识别静态是实际勘探工作能否顺利进行的一个关键,同时也关系到测深工作的成败。为了计算二维静位移地电模型的大地电磁响应,本文采用非均匀网格有限差分法进行了数值模拟。从电磁场满足的微分方程出发,利用非均匀网格有限差分法导出了二维大地电磁正演计算的线性方程组。通过对浅部不均匀体模型的模拟计算,得到了大地电磁的静位移响应规律,同时也验证了非均匀网格有限差分正演算法的准确性。
Abstract: Because of the existence of shallow inhomogeneous bodies, static shift occurs in magnetotelluric sounding curves, which makes the interpretation of magnetotelluric data very difficult, or even impossible. It is very important to identify static shift in the interpretation of magnetotelluric data. In this paper, in order to compute the two-dimensional magnetotelluric responses of static shift model, non-uniform grids finite difference method was adopted for numerical results. From the electric field and magnetic field based on the differential equation, the linear equation for magnetotelluric numerical modeling was derived by non-uniform meshes finite difference method. By the numerical simulation for the geo-electric model with shallow inhomogeneous body, the character of static shift can be obtained and the magnetotelluric responses computed by non-uniform girds finite difference method are effective.
文章引用:邓小康, 童孝忠. 基于非均匀网格有限差分法的大地电磁静位移模拟[J]. 地球科学前沿, 2018, 8(8): 1353-1361. https://doi.org/10.12677/AG.2018.88148

参考文献

[1] 周邵民, 柳建新, 孙欢乐. 基于非负最小二乘法的一维MT正则化反演研究[J]. 工程地球物理学报, 2017, 14(3): 253-261.
[2] 伍亮, 李桐林, 朱成, 等. 大地电磁测深法中静态效应及其反演[J]. 地球物理进展, 2015, 30(2): 840-846.
[3] 李国瑞, 席振铢, 龙霞. 基于实测积累电荷静电场消除静态效应的方法[J]. 物探与化探, 2017, 41(4): 730-735.
[4] 陈辉, 王春庆, 孟小红, 等. 均匀大地中低阻静态模型及其静校正效果[J]. 工程地球物理学报, 2007, 4(4): 282-286.
[5] 冯德山,王珣. 大地电磁双二次插值FEM正演及最小二乘正则化联合反演[J]. 中国有色金属学报,2013,23(9):2524-2531.
[6] 田红军, 佟铁钢. 带地形的二维各向异性大地电磁测深数值模拟[J]. 工程地球物理学报, 2015, 12(2): 139-144.
[7] 叶益信, 胡祥云, 何梅兴, 等. 电导率分块均匀大地电磁场二维有限元数值模拟[J]. 工程地球物理学报, 2009, 6(1): 9-14.
[8] de Groot-Hedlin, C. (2006) Finite-Difference Modeling of Magnetotelluric Fields: Error Estimates for Uniform and Nonuniform Grids. Geophysics, 71, 97-106. [Google Scholar] [CrossRef
[9] Xiao, Q.B. and Zhao, G.Z. (2010) Comparison of Finite Difference Numerical Solutions in Magnetotelluric Modeling. Chinese Journal of Geophysics, 53, 622-630.
[10] Rao, K.P. and Babu, G.A. (2006) EMOD2D—A Program in C++ for Finite Difference Modeling of Magnetotelluric TM Mode Responses over 2D Earth. Computer & Geosciences, 32, 1499-1511. [Google Scholar] [CrossRef
[11] 王涛, 柳建新, 童孝忠, 等. 有限单元法与有限差分法在MT一维正演模拟中的对比分析[J]. 物探化探计算技术, 2013, 35(5): 538-543.
[12] Zhdanov, M. and Tolstaya, E. (2004) Minimum Support Nonlinear Parametrization in the Solution of a 3D Magnetotelluric Inverse Problem. Inverse Problem, 20, 937-952.
[13] Wannamaker, P.E. (1991) Advances in Three-Dimensional Magnetotelluric Modeling Using Integral Equations. Geophysics, 56, 1716-1728. [Google Scholar] [CrossRef
[14] 张辉, 唐新功. 一种新型有限差分网格剖分方法在大地电磁一维正演中的应用[J]. 物探与化探, 2015, 39(3): 562-566.
[15] 柳建新, 童孝忠, 郭荣文, 等. 大地电磁测深法勘探[M]. 北京: 科学出版社, 2012.