Journal of Jilin University(Earth Science Edition) ›› 2017, Vol. 47 ›› Issue (1): 215-223.doi: 10.13278/j.cnki.jjuese.201701301

Previous Articles     Next Articles

Magnetotelluric Responses of Two-Dimensional Complex Conductivity Structures

Luo Tianya1, Xiong Bin1, Cai Hongzhu2, Chen Xin1, Liu Yunlong1, Lan Huaikang1, Li Zuqiang1, Liang Zhuo1   

  1. 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: 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

CLC Number: 

  • P631.3
[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.
[1] Zhang Bo, Cao Hongkai, Sun Jianmeng, Zhang Pengyun, Yan Weichao. Numerical Simulation of Response Characteristics of Array Induction Logging in Heavy Oil Thermal Recovery Formation [J]. Journal of Jilin University(Earth Science Edition), 2018, 48(4): 1277-1286.
[2] Yin Changchun, Yang Zhilong, Liu Yunhe, Zhang Bo, Qi Yanfu, Cao Xiaoyue, Qiu Changkai, Cai Jing. Characteristics of 3D DC Resistivity Response for Arbitrary Anisotropic Models Using Circular Scanning Measurement [J]. Journal of Jilin University(Earth Science Edition), 2018, 48(3): 872-880.
[3] Li Jianping, Weng Aihua, Li Shiwen, Li Dajun, Li Sirui, Yang Yue, Tang Yu, Zhang Yanhui. 3-D Forward Method for Geomagnetic Depth Sounding Based on Finite Difference Method in Spherical Coordinate [J]. Journal of Jilin University(Earth Science Edition), 2018, 48(2): 411-419.
[4] Zeng Zhaofa, Huo Zhijun, Li Wenben, Li Jing, Zhao Xueyu, He Rongqin. 3D Vector Finite-Element Airborne Electromagnetic Modelling in an Arbitrary Anisotropic Medium [J]. Journal of Jilin University(Earth Science Edition), 2018, 48(2): 433-444.
[5] An Zhenfang, Zhang Jin, Zhang Jianzhong. Geometry Optimization Design of Three Dimensional Marine Vertical Cable [J]. Journal of Jilin University(Earth Science Edition), 2018, 48(1): 271-284.
[6] Weng Aihua, Li Sirui, Yang Yue, Li Dajun, Li Jianping, Li Shiwen. Basic Principle, Current Status and Prospect of Magnetometric Resistivity [J]. Journal of Jilin University(Earth Science Edition), 2017, 47(6): 1838-1854.
[7] Li Dajun, Weng Aihua, Yang Yue, Li Sirui, Li Jianping, Li Shiwen. Three-Dimension Forward Modeling and Characteristics for Surface-Borehole Transient Electromagnetic by Using Staggered-Grid Finite Difference Method [J]. Journal of Jilin University(Earth Science Edition), 2017, 47(5): 1552-1561.
[8] Gao Jianshen, Sun Jianmeng, Jiang Yanjiao, Cui Likai. Effect of Electrode Array Structures in Laterolog and a New Array Measurement Method [J]. Journal of Jilin University(Earth Science Edition), 2016, 46(6): 1874-1883.
[9] Yan Jiabin, Huang Xiangyu. Vector Finite Element Method [J]. Journal of Jilin University(Earth Science Edition), 2016, 46(5): 1538-1549.
[10] Liu Haifei, Liu Jie, Gao Han, Guo Rongwen, Tong Xiaozhong, Ma Changying. 3D Finite Element Forward Modeling of Five-Pole Longitude Induced Polarization Sounding [J]. Journal of Jilin University(Earth Science Edition), 2016, 46(3): 884-892.
[11] Li Qingyang, Li Zhenchun, Huang Jianping, Li Na, Su Zairong. Factor Analysis of Seismic Modeling with Topography Based on a Fully Staggered Body-Fitted Grids [J]. Journal of Jilin University(Earth Science Edition), 2016, 46(3): 920-929.
[12] Yin Changchun, Qiu Changkai, Liu Yunhe, Cai Jing. Weighted Laterally-Constrained Inversion of Time-Domain Airborne Electromagnetic Data [J]. Journal of Jilin University(Earth Science Edition), 2016, 46(1): 254-261.
[13] Jiang Fuyu, Xie Leilei, Chang Wenkai, Huang Yan, Zhang Zuohong. Forward Calculation of Three Dimensional Gravity Vector Using Finite Element Method [J]. Journal of Jilin University(Earth Science Edition), 2015, 45(4): 1217-1226.
[14] Zhu Chao, Xia Zhiyuan, Wang Chuanwu, Song Guangyong, Wei Xuebin, Wang Peng, Wang Haifeng, Wang Bo. Seismic Prediction for Sweet Spot Reservoir of Tight Oil [J]. Journal of Jilin University(Earth Science Edition), 2015, 45(2): 602-610.
[15] Tang Wenwu, Liu Jianxin,Tong Xiaozhong. Finite Element Forward Modeling on Line Source of FCSEM for Continuous Conductivity [J]. Journal of Jilin University(Earth Science Edition), 2013, 43(5): 1646-1654.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
No Suggested Reading articles found!