复杂电性结构大地电磁二维响应特征

doi:10.13278/j.cnki.jjuese.201701301

吉林大学学报(地球科学版) ›› 2017, Vol. 47 ›› Issue (1): 215-223.doi: 10.13278/j.cnki.jjuese.201701301

复杂电性结构大地电磁二维响应特征

罗天涯¹, 熊彬¹, 蔡红柱², 陈欣¹, 刘云龙¹, 兰怀慷¹, 李祖强¹, 梁卓¹

1. 桂林理工大学地球科学学院, 广西桂林 541006;
2. 犹他大学矿产和地球科学学院, 美国犹他州盐湖城 84112

收稿日期:2016-05-08 出版日期:2017-01-26 发布日期:2017-01-26
通讯作者: 熊彬(1974-),男,教授,博士,主要从事电磁场理论及反演成像方面的教学与科研工作,E-mail:xiongbin@msn.com E-mail:xiongbin@msn.com
作者简介:罗天涯(1989-),男,博士研究生,主要从事电磁场数值模拟研究,E-mail:tianyaluo@aliyun.com
基金资助:
国家自然科学基金项目（40974077,41164004,41674075）；广西自然科学基金项目（2016GXNSFA380004,2013GXNSFAA019277）；桂林市“漓江学者”专项经费资助项目（2013005）；桂林理工大学研究生创新项目（BS201601）

Magnetotelluric Responses of Two-Dimensional Complex Conductivity Structures

Luo Tianya¹, Xiong Bin¹, Cai Hongzhu², Chen Xin¹, Liu Yunlong¹, Lan Huaikang¹, Li Zuqiang¹, Liang Zhuo¹

1. College of Earth Sciences, Guilin University of Technology, Guilin 541006, Guangxi, China;
2. College of Mines & Earth Sciences, University of Utah, Salt Lake City 84112, UT, USA

Received:2016-05-08 Online:2017-01-26 Published:2017-01-26
Supported by:
Supported byNational Natural Science Foundation of China (40974077,41164004, 41674075), Natural Science Foundation of Guangxi Province (2016GXNSFA380004, 2013GXNSFAA019277), Lijiang Scholars Special Expenditure Grant (2013005)and Innovation Project of GLUT Graduate Education (BS201601)

摘要/Abstract

摘要： 为了改进计算区域离散化问题,本文利用自适应非结构化网格有限单元法求解二维地电结构下大地电磁场满足的加权余量表达式。在有限元求解电磁场的过程中,网格剖分越精细、计算精度越高,计算量也会越大。此外,结构化网格难以适应任意地形以及复杂地质构造。而自适应非结构化网格在电性变化剧烈的区域会自动加密,在电性缓变的区域则生成粗疏的网格,从而优化网格质量与数量。因此,文中引入COMSOL Multiphysics软件,以实现若干地电模型的构建及非结构化自由四边形单元网格化。将网格数据信息导入本文算法,计算大地电磁场响应,并与解析解及数值解对比。结果表明,基于非结构化网格的正演模拟精度高、适应性强,为计算区域网格化提供了新的方法。

关键词: 非结构化网格, 大地电磁, 有限单元法, 正演

Abstract: To improve the discretization technique, adaptive unstructured quadrilateral grid finite element (FE) procedure is presented for solving weighted residual formulation of magnetotelluric (MT) fields in 2-D conductivity structures. The finer the meshes are, the more accurate the numerical solutions are. However, this is at the cost of time-consuming and large memory requirement. Adaptive mesh refinement generally provides an optimization of both run time and accuracy since in the method meshes are only refined where the electrical parameters vary intensely, while coarsen in the slow electrical resistivity change region. Due to the fact that regular grids have no potential to incorporate arbitrary topography and complex geologic structures, we introduced COMSOL Multiphysics software package to build arbitrary geometries and subdivide them with free quadrilaterals, then apply meshes information to our algorithm. Comparison with analytical solutions and previous numerical solutions demonstrate a high-accuracy and flexibility of numerical simulations based on unstructured grid, a new technology to discretize the computating domain.

Key words: unstructured grid, magnetotelluric(MT), finite element method, forward

中图分类号:

P631.3

罗天涯, 熊彬, 蔡红柱, 陈欣, 刘云龙, 兰怀慷, 李祖强, 梁卓. 复杂电性结构大地电磁二维响应特征[J]. 吉林大学学报(地球科学版), 2017, 47(1): 215-223.

Luo Tianya, Xiong Bin, Cai Hongzhu, Chen Xin, Liu Yunlong, Lan Huaikang, Li Zuqiang, Liang Zhuo. Magnetotelluric Responses of Two-Dimensional Complex Conductivity Structures[J]. Journal of Jilin University(Earth Science Edition), 2017, 47(1): 215-223.

参考文献

[1] 赵广茂,李桐林,王大勇,等. 基于二次场二维起伏地形MT有限元数值模拟[J]. 吉林大学学报(地球科学版),2008,38(6):1055-1059. Zhao Guangmao, Li Tonglin, Wang Dayong, et al. Secondary Field-Based Two-Dimensional Topographic Numerical Simulation in Magnetotellurics by Finite Element Method[J]. Journal of Jilin University (Earth Science Edition), 2008, 38(6):1055-1059.
[2] 顾观文,吴文鹂,李桐林. 大地电磁场三维地形影响的矢量有限元数值模拟[J]. 吉林大学学报(地球科学版),2014,44(5):1678-1686. Gu Guanwen, Wu Wenli, Li Tonglin. Modeling for the Effect of Magnetotelluric 3D Topography Based on the Vector Finite-Element Method[J]. Journal of Jilin University (Earth Science Edition), 2014, 44(5):1678-1686.
[3] 李桐林,张镕哲,朴英哲.大地电磁测深与地震初至波走时交叉梯度反演[J].吉林大学学报(地球科学版),2015,45(3):952-961. Li Tonglin, Zhang Rongzhe, Pak Yongchol. Joint Inversion of Magnetotelluric and First-Arrival Wave Seismic Traveltime with Cross-Gradient Constraints[J]. Journal of Jilin University (Earth Science Edition), 2015, 45(3):952-961.
[4] Zhdanov M S, Lee S K, Yoshioka K. Integral Equa-tion Method for 3D Modeling of Electromagnetic Fields in Complex Structures with Inhomogeneous Background Conductivity[J]. Geophysics, 2006, 71(6):333-345.
[5] Newman G A, Alumbaugh D L. Frequency-Domain Modelling of Airborne Electromagnetic Responses Using Staggered Finite Differences[J]. Geophysical Prospecting, 1995, 43(8):1021-1042.
[6] Aprea C, Booker J R, Smith J T. The Forward Pro-blem of Electromagnetic Induction:Accurate Finite-Difference Approximations for Two Dimensional Discrete Boundaries with Arbitrary Geometry[J]. Geophys J Int, 1997, 129(1):29-40.
[7] Cai H Z, Xiong B, Han M, et al. 3D Controlled-Source Electromagnetic Modeling in Anisotropic Medium Using Edge-Based Finite Element Method[J]. Computers & Geosciences, 2014, 73:164-176.
[8] Key K, Weiss C. Adaptive Finite-Element Modeling Using Unstructured Grids:The 2D Magnetotelluric Example[J]. Geophysics, 2006, 71(6):G291-G299.
[9] Fox R C, Hohmann G W, Killpack T J, et al.Topo-graphic Effects in Resistivity and Induced-Polarization Surveys[J]. Geophysics, 1980, 45(1):75-93.
[10] Wannamaker P E, Stodt J A, Rijo L. A Stable Finite Element Solution for Two-Dimensional Magneto-telluric Modelling[J]. Geophys J Int, 1987, 88(1):277-296.
[11] Reddy I K, Rankin D. Magnetotelluric Response of Laterally Inhomogeneous and Anisotropic Media[J]. Geophysics, 1975, 40(6):1035-1045.
[12] Franke A, B rner R U, Spitzer K. Adaptive Unstru-ctured Grid Finite Element Simulation of Two-Dimensional Magnetotelluric Fields for Arbitrary Surface and Seafloor Topography[J]. Geophys J Int, 2007, 171(1):71-86.
[13] McFee S, Giannacopoulos D. Optimal Discretizations in Adaptive Finite Element Electromagnetics[J]. Int J Numer Methods Eng, 2001, 52(9):939-978.
[14] Badea E A, Everett M E, Newman G A, et al.Finite-Element Analysis of Controlled-Source Electro-magnetic Induction Using Coulomb-Gauged Potentials[J]. Geophysics, 2001, 66(3):786-799.
[15] Li Y, Pek J.Adaptive Finite Element Modelling of Two-Dimensional Magnetotelluric Fields in General Anisotropic Media[J]. Geophys J Int, 2008, 175(3):942-954.
[16] Puzyrev V, Koldan J, Puente J, et al.A Parallel Finite-Element Method for Three-Dimensional Controlled-Source Electromagnetic Forward Modelling[J]. Geophys J Int, 2013, 193(2):678-693.
[17] Lee S K, Kim H J, Song Y, et al. MT2D Inv Matlab:A Program in MATLAB and FORTRAN for Two-Dimensional Magnetotelluric Inversion[J]. Computers & Geosciences, 2009, 35(8):1722-1734.
[18] Rodi W L. A Technique for Improving the Accuracy of Finite Element Solutions for Magnetotelluric Data[J]. Geophys J Int, 1976, 44(2):483-506.
[19] 陈乐寿. 有限元法在大地电磁场正演计算中的应用及改进[J]. 石油勘探,1981(3):84-103. Chen Leshou. Application and Improvement of Finite Element Method in Forward Calculation of Geo-Electromagnetic Field[J]. Geophysical Prospecting for Petroleum, 1981(3):84-103.
[20] 邓建辉,熊文林,葛修润. 复杂区域非结构化四边形网格全自动生成方法[J]. 计算结构力学及应用,1995,12(2):196-204. Deng Jianhui, Xiong Wenlin, Ge Xiurun. Automatic Generation of Unstructured Quadrilateral Mesh for Complex Domain[J]. Computational Structural Mechanics and Application, 1995, 12(2):196-204.
[21] 陈建军,郑耀,陈立岗,等. 非结构化四边形网格生成新算法[J]. 中国图像图形学报,2008,13(9):1796-1803. Chen Jianjun, Zheng Yao,Chen Ligang, et al. A New Unstructured Quadrilateral Mesh Generation Algorithm[J]. Journal of Image and Graphics, 2008, 13(9):1796-1803.
[22] Jin J M.The Finite Element Method in Electroma-gnetics[M]. New Jersey:Wiley-IEEE Press, 2002:1-780.
[23] Wannamaker P E, Stodt J A, Rijo L. Two-Dimen-sional Topographic Responses in Magnetotellurics Modeled Using Finite Elements[J]. Geophysics, 1986, 51(11):2131-2144.
[24] 李磊. 湘南骑田岭锡铅锌多金属矿区岩矿石电性研究[J]. 物探与化探,2007,31(增刊1):77-80. Li Lei. Researches on Rock Electrical Properties in the Qitianling Tin,Lead and Zinc Polymetallic Ore Deposit,Southern Hunan[J]. Geophysical and Geochemical Exploration,2007, 31(Sup. 1):77-80.
[25] Chouteau M, Bouchard K. Two-Dimensional Terrain Correction in Magnetotelluric Surveys[J]. Geophysics, 1988, 53(6):854-862.
[26] Jiracek G R. Near-Surface and Topographic Distor-tions in Electromagnetic Induction[J]. Surveys in Geophysics, 1990, 11:163-203.

Metrics

Viewed

Full text

Abstract

Cited

Shared

Discussed

复杂电性结构大地电磁二维响应特征

Magnetotelluric Responses of Two-Dimensional Complex Conductivity Structures

PDF (PC)

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

Metrics

本文评价

推荐阅读 10

[1]	张阳阳, 杜威. 基于超收敛单元片梯度恢复的自适应有限元法音频大地电磁正演[J]. 吉林大学学报(地球科学版), 2026, 56(2): 673-683.
[2]	黄兴国, 翁央央, 韩丽. 基于广义递归卷积的孔隙黏弹地震波正演模拟[J]. 吉林大学学报(地球科学版), 2026, 56(1): 377-385.
[3]	邵逸飞, 龚启舟, 许权张越, 王言章, 王世隆. 基于八叉树网格的时间域航空电磁深切割地形三维正演模拟[J]. 吉林大学学报(地球科学版), 2025, 55(6): 2120-2131.
[4]	龙刚, 沈金松, 苏朝阳, 冉尚. 电各向异性海底多金属硫化矿地层时域电磁响应模拟与特征分析[J]. 吉林大学学报(地球科学版), 2025, 55(1): 274-288.
[5]	雷松达, 王显祥, 刘遂明. 复杂地质条件下浅海区电性源瞬变电磁法三维响应特征[J]. 吉林大学学报(地球科学版), 2024, 54(3): 1016-1030.
[6]	余年, 张学健, 孔文新, 韦昌. 地磁感应电流多学科交叉研究进展[J]. 吉林大学学报(地球科学版), 2024, 54(1): 54-67.
[7]	张鑫崇, 殷长春. 基于Debye模型的时间域航空电磁激电效应正演模拟 [J]. 吉林大学学报(地球科学版), 2023, 53(5): 1573-1581.
[8]	胡向阳, 秦敏, 张恒荣, 杨冬, 袁伟. 基于SVR的定向井随钻电磁波测井层厚影响快速计算[J]. 吉林大学学报(地球科学版), 2023, 53(2): 589-.
[9]	于明浩, 巩向博, 万晓杰. 复杂地表相控震源正演模拟与照明分析[J]. 吉林大学学报(地球科学版), 2022, 52(6): 2034-2041.
[10]	骆祖江, 成磊, 张兴旺, 于丹丹, 徐成华. 悬挂式止水帷幕深基坑降水方案模拟优化[J]. 吉林大学学报(地球科学版), 2022, 52(6): 1946-1956.
[11]	罗明, 何叶, 周江江, 贾培蒙, 吕华星. 砾砂岩不均匀展布引起的构造畸变[J]. 吉林大学学报(地球科学版), 2022, 52(4): 1360-.
[12]	高科宁, 邓居智, 刘文玉, 余辉, 陈辉, 周聪, 王彦国, 刘宝山. 黑龙江省三合屯金矿床深部电性结构与成矿模式 [J]. 吉林大学学报(地球科学版), 2022, 52(3): 684-700.
[13]	王新宇, 严良俊, 毛玉蓉, 黄鑫, 谢兴兵, 周磊 . 起伏地形条件下长偏移距瞬变电磁三维正演 [J]. 吉林大学学报(地球科学版), 2022, 52(3): 754-765.
[14]	赵东东, 王旭龙, 周印明, 戴世坤. 基于泊松方程的空间波数混合域二度体磁异常数值模拟[J]. 吉林大学学报(地球科学版), 2022, 52(2): 592-601.
[15]	张志立, 韩复兴, 孙文艳, 王怡, 杨安琪, 焦艳艳, 薛诗桂. 利用正演模拟实现面波衰减[J]. 吉林大学学报(地球科学版), 2021, 51(6): 1890-1896.