大地电磁三维正演聚集多重网格算法

doi:10.13278/j.cnki.jjuese.20160359

Abstract

Abstract: To speed up 3D magnetotelluric (MT) modelling, we introduce a novel algebraic multigrid-aggregation-based algebraic multigrid method (AGMG) into three dimensional forward modeling of magnetotelluric. We used the finite-volume algorithm based on the Yee's grids to discretize quasi-static Maxwell's equations with Dirichlet boundary conditions,and used the AGMG method to solve the final large sparse linear equation system in electric field. Through the coarsening of AGMG and the aggregating based on N-passes of a pairwise matching algorithm applied to the matrix graph,we proposed three different AGMG algorithms:1) the classic V-cycle AGMG algorithm;2) the AGMG preconditioned conjugate gradient algorithm (AGMG-CG); 3) the AGMG pretreated generalized conjugate residual method (AGMG-GCR). We performed 3D MT modelling for typical geo-electric models with different iteration,and analyzed the features of the AGMG techniques through comparing with ModEM algorithm. The results show that the AGMG methods are accurate and robust, the AGMG preconditioner improves the convergence of the classic V-cycle AGMG and Krylov subspace methods greatly. The AGMG-GCR method is the most effective one presented in this paper,which speeds up the modeling by ten times more than the ModEM codes for large-scale grids (144×152×104). The AGMG-GCR is especially suitable for large-scale 3D MT modeling because of its high precision, fast convergence, and robust iteration.

Key words: magnetotelluric (MT), three-domain modelling, multigrid method, finite volume method

CLC Number:

P631.3

Chen Hui, Yin Min, Yin Changchun, Deng Juzhi. Three-Dimensional Magnetotelluric Modelling Using Aggregation-Based Algebraic Multigrid Method[J].Journal of Jilin University(Earth Science Edition), 2018, 48(1): 261-270.

References

[1] Newman G A. A Review of High-Performance Com-putational Strategies for Modeling and Imaging of Electromagnetic Induction Data[J]. Surveys in Geophysics, 2014, 35(1): 85-100.
[2] Smith R. Electromagnetic Induction Methods in Mi-ning Geophysics from 2008 to 2012[J]. Surveys in Geophysics, 2014, 35(1): 123-156.
[3] Siripunvaraporn W. Three-Dimensional Magnetotellu-ric Inversion: An Introductory Guide for Developers and Users[J]. Surveys in Geophysics, 2012, 33(1): 5-27.
[4] 谭捍东,余钦范,Booker John,等. 大地电磁法三维交错采样有限差分数值模拟[J]. 地球物理学报,2003,46(5):706-711. Tan Handong, Yu Qinfan, Booker J, et al. Magnetotelluric Three-Dimensional Modelling Using the Staggered-Grid Fnite Difference Method[J]. Chinese Journal of Geophysics, 2003, 46(5): 706-711.
[5] 沈金松. 用交错网格有限差分法计算三维频率域电磁响应[J]. 地球物理学报,2003,46(2):281-289. Shen Jinsong. Modelling of 3-D Eectromagnetic Responses in Frequency Domain by Using Straggered-Grid Finite Difference Method[J]. Chinese Journal of Geophysics, 2003, 46(2): 281-289.
[6] Mackie R L, Madden T R, Wannamaker P E. Three-Dimensional Magnetotelluric Modeling Using Difference Equations: Theory and Comparisons to Integral Equation Solutions[J]. Geophysics, 1993, 58(2): 215-226.
[7] 李焱, 胡祥云, 杨文采, 等. 大地电磁三维交错网格有限差分数值模拟的并行计算研究[J]. 地球物理学报,2012,55(12):4036-4043. Li Yan, Hu Xiangyun, Yang Wencai, et al. A Study on Parallel Computation for 3D Magnetotelluric Modeling Using the Staggered-Grid Finite Difference Method[J]. Chinese Journal of Geophysics, 2012, 55(12): 4036-4043.
[8] Haber E, Ascher U M. Fast Finite Volume Simulation of 3D Electromagnetic Problems with Highly Discontinuous Coefficients[J]. SIAM Journal on Scientific Computing, 2000, 22(6): 1943-1961.
[9] Haber E, Ruthotto L. A Multiscale Finite Volume Method for Maxwell's Equations at Low Frequencies[J]. Geophysical Journal International, 2014, 199(2): 1268-1277.
[10] 陈辉,殷长春,邓居智. 基于Lorenz规范条件下磁矢势和标势耦合方程的频率域电磁法三维正演[J]. 地球物理学报,2016,59(8):3087-3097. Chen Hui, Yin Changchun, Deng Juzhi. A Finite-Volume Solution to 3D Frequency-Domain Electromagnetic Modelling Using Lorenz-Gauged Magnetic Vector and Scalar Potentials[J]. Chinese Journal of Geophysics, 2016, 59(8): 3087-3097.
[11] Ren Z, Kalscheuer T, Greenhalgh S, et al. A Finite-Element-Based Domain-Decomposition Approach for Plane Wave 3D Electromagnetic Modeling[J]. Geophysics, 2014, 79(6): E255-E268.
[12] 黄临平,戴世坤. 复杂条件下3D电磁场有限元计算方法[J]. 地球科学:中国地质大学学报,2002,27(6):775-779. Huang Linping, Dai Shikun. Finite Element Calculation Method of 3D Electromagnetic Field under Complex Condition[J]. Earth Science: Journal of China University of Geoscineces, 2002, 27(6): 775-779.
[13] Mitsuhata Y, Uchida T. 3D Magnetotelluric Mode-ling Using the T-Omega Finite-Element Method[J]. Geophysics, 2004, 69(1): 108-119.
[14] 李俊杰,严家斌,皇祥宇. 无单元Galerkin法大地电磁三维正演模拟[J]. 地质与勘探,2015,51(5):946-952. Li Junjie, Yan Jiabin, Huang Xiangyu. Three-Dimensional Forward Modeling of Magnetotellurics Using the Element-Free Galerkin Method[J]. Geology and Exploration, 2015, 51(5): 946-952.
[15] 严家斌,皇祥宇. 大地电磁三维矢量有限元正演[J]. 吉林大学学报(地球科学版),2016,46(5):1538-1549. Yan Jiabin, Huang Xiangyu. 3D Forward Modeling of Magnetotelluric Field by Vector Finite Element Method[J]. Journal of Jilin University (Earth Science Edition), 2016, 46(5): 1538-1549.
[16] Kruglyakov M, Geraskin A, Kuvshinov A. Novel Accurate and Scalable 3-D MT Forward Solver Based on a Contracting Integral Equation Method[J]. Computers & Geosciences, 2016, 96: 208-217.
[17] Wannamaker P E. Advances in Three-Dimensional Magnetotelluric Modeling Using Integral Equations[J]. Geophysics, 1991, 56(11): 1716-1728.
[18] 王书明,李德山,胡浩. 三维/三维构造下大地电磁相位张量数值模拟[J]. 地球物理学报,2013,56(5):1745-1752. Wang Shuming, Li Deshan, Hu Hao. Numerical Modeling of Magnetotelluric Phase Tensor in the Context of 3D/3D Formation[J]. Chinese Journal of Geophysics, 2013, 56(5): 1745-1752.
[19] de Groot-Hedlin C. Finite-Difference Modeling of Magnetotelluric Felds: Error Estimates for Uniform and Nonuniform Grids[J]. Geophysics, 2006, 71(3): G225-G233.
[20] Han N, Nam M J, Kim H J, et al. A Comparison of Accuracy and Computation Time of Three-Dimensional Magnetotelluric Modelling Algorithms[J]. Journal of Geophysics and Engineering, 2009, 6(2): 136.
[21] Smith J T. Conservative Mmodeling of 3-D Electro-magnetic Fields: Part Ⅱ: Biconjugate Gradient Solution and an Accelerator[J]. Geophysics, 1996, 61(5): 1319-1324.
[22] Siripunvaraporn W, Egbert G, Lenbury Y. Nume-rical Accuracy of Magnetotelluric Modelling: A Comparison of Finite Difference Approximations[J]. Earth Planets Space, 2002, 54: 721-725.
[23] Mackie R L, Madden T R. Conjugate Direction Relaxation Solutions for 3-D Magnetotelluric Modeling[J]. Geophysics, 1993, 58(7): 1052-1057.
[24] Weiss C J, Newman G A. Electromagnetic Induction in a Generalized 3D Anisotropic Earth: Part 2: The LIN Preconditioner[J]. Geophysics, 2003, 68(3): 922-930.
[25] 陈辉,邓居智,谭捍东,等. 大地电磁三维交错网格有限差分数值模拟中的散度校正方法研究[J]. 地球物理学报,2011,54(6):1649-1659. Chen Hui, Deng Juzhi, Tan Handong, et al. Study on Divergence Correction Method in Three-Dimensional Magnetotelluric Modeling with Staggered-Grid Finite Fifference Method[J]. Chinese Journal Geophysics, 2011, 54(6): 1649-1659.
[26] Streich R. 3D Finite-Difference Frequency-Domain Modeling of Controlled-Source Electromagnetic Data: Direct Solution and Optimization for High Accuracy[J]. Geophysics, 2009, 74(5): 95-105.
[27] Puzyrev V, Koric S, Wilkin S. Evaluation of Parallel Direct Sparse Linear Solvers in Electromagnetic Geophysical Problems[J]. Computers & Geosciences, 2016, 89: 79-87.
[28] Koldan J, Puzyrev V, de la Puente J, et al. Algebraic Multigrid Preconditioning Within Parallel Finite-Element Solvers for 3-D Electromagnetic Modelling Problems in Geophysics[J]. Geophysical Journal International, 2014, 197(3): 1442-1458.
[29] Mulder W A. Geophysical Modelling of 3D Electro-magnetic Diffusion with Multigrid[J]. Computing and Visualization in Science, 2008, 11(3): 129-138.
[30] Pan K, Tang J. 2.5-D and 3-D DC Resistivity Modelling Using an Extrapolation Cascadic Multigrid Method[J]. Geophysical Journal International, 2014, 197(3): 1459-1470.
[31] Trottenberg U, Clees T. Multigrid Software for Industrial Applications: From MG00 to SAMG[M]. Heidelberg: Springer, 2009: 423-436.
[32] Notay Y. An Aggregation-Based Algebraic Multigrid Method[J]. Electronic Transactions on Numerical Analysis, 2010, 37(6): 123-146.
[33] Henson V E, Yang U M. Boomer AMG: A Parallel Algebraic Multigrid Solver and Preconditioner[J]. Applied Numerical Mathematics, 2002, 41(1): 155-177.
[34] Saad Y. Iterative Methods for Sparse Linear Systems[M]. Philadelphia: SIAM, 2003.
[35] Pflaum C. A Multigrid Conjugate Gradient Method[J]. Applied Numerical Mathematics, 2008, 58: 1803-1817.
[36] MT 3D Inversion Workshop. Dublin Test Model 1(DTM1)[EB/OL].[2016-10-20] http://www.complete-mt-solutions.com/mtnet/workshops/mt3di-nv/2008_Dublin/Dublin/3dmodel.html.
[37] Kelbert A, Meqbel N, Egbert G D, et al. ModEM: A Modular System for Inversion of Electromagnetic Geophysical Data[J]. Computers & Geosciences, 2014, 66: 40-53.

Metrics

Viewed

Full text

Abstract

Cited

Shared

Discussed

Comments

Recommended 10

[1]	ZHANG Shi-guang, LIU Cheng-zhi, LEI Shuang-fang, ZHANG Yan, WU Gao-ping, LIU Qiu-hong. The Application of High-Resolution Sequence Stratigraphy in Multiplex Deposition System of the River, Lake and Delta--To Take the Development Block of Fuyu Oil Layer in Chaoyanggou Reservoir for Example[J]. J4, 2009, 39(3): 361 -368 .
[2]	Li Bi-le，Shen Xin，Chen Guang-jun，Yang Yan-qian，Li Yong-sheng. Geochemical Features of OreForming Fluids and Metallogenesis of Vein Ⅰ in Asiha Gold Ore Deposit，Eastern Kunlun，Qinghai Province[J]. Journal of Jilin University(Earth Science Edition), 2012, 42(6): 1676 -1687 .
[3]	YANG Xiao-ping, LI Yang-chun,LIU Zhen, WANG Yan, WANG Hong-jie. Classification of Tectonic Sequence and Dynamic Evolution of Jixi Basin, Eastern Heilongjiang Province[J]. J4, 2005, 35(05): 616 -621 .
[4]	LI Fa-wen, FENG Ping, ZHANG Chao. Vertically-Mixed Runoff Model and Application in Beisan River of Tianjin[J]. J4, 2011, 41(2): 459 -464 .
[5]	WANG Li-jun,LIU Guo-cai,HUANG Ji-guo,LI Yue-qian,CONG Ying,SHEN Zhao-li,ZHAO Xiao-bo. Examination and Study of Catalytic Oxidation in Sight Water Resumption[J]. J4, 2006, 36(03): 458 -461 .
[6]	Lei Ruxiong, Wu Changzhi, Qu Xun, Gu Lianxing, Chen Gang,Uerna Amy, Sun Hongtao, Liu Guoling. Geochronology, Geochemistry and Zircon Hf Isotope Compositions of the Ore-Bearing Gneiss Granite in the Tianhudong Iron-Molybdenum Ore Deposit in the Central Tianshan, West China: Implications for the Early Paleozoic Tectonic Evolution of Central Tianshan[J]. Journal of Jilin University(Earth Science Edition), 2014, 44(5): 1540 -1552 .
[7]	CHU Feng-you, HU Da-qian, YAO Jie. Mineral Composition and Element Geochemistry of Co-Rich Nodule from YJB Sea Mount in Central Pacific Ocean[J]. J4, 2007, 37(1): 8 -0014 .
[8]	ZHOU Yan，ZHENG Pei-xi，WANG Tie-fu，ZHANG Yan-jie. H-O Isotopes of Gold Deposits in the Hanging Wall of Zhao-Ping Fault Belt[J]. J4, 2007, 37(4): 668 -0671 .
[9]	TANG Hua-feng,WANG Pu-jun,JIANG Chuan-jin,LIU Jie,ZHANG Qing-chen,FENG You-liang. Seismic Characters of Volcanic Facies and Their Distribution Relation to Deep Faults in Songliao Basin[J]. J4, 2007, 37(1): 73 -0078 .
[10]	FAN Xiao-min, LI Zhou-bo. Analysis of Fluctuation and Amplitude Increase in Sonic Logs Caused by Fractures in Carbonate Rocks[J]. J4, 2007, 37(1): 168 -0173 .

Three-Dimensional Magnetotelluric Modelling Using Aggregation-Based Algebraic Multigrid Method

PDF (PC)

Abstract

Cite this article

share this article

References

Related Articles 4

Metrics

Comments

Recommended 10

[1]	Shao Yifei , Gong Qizhou, Xu Quan, Zhang Yue, Wang Yanzhang, Wang Shilong. 3D Forward Modeling of Time-Domain Airborne Electromagnetic Reponses in Deeply Incised Terrain Based on Octree Mesh [J]. Journal of Jilin University(Earth Science Edition), 2025, 55(6): 2120-2131.
[2]	Lu Yan, Liu Zongbin, Liao Xinwu, Li Chao, Wang Ya. Fluid Flow Characteristics in Low-Permeability Sandstone Reservoirs [J]. Journal of Jilin University(Earth Science Edition), 2025, 55(4): 1077-1090.
[3]	LIU Jian-xin, GUO Rong-wen, TONG Xiao-zhong, LIU Ying, LIU Peng-mao. General Fourier Analysis of Multi-Grid in Magnetotelluric Modeling [J]. J4, 2011, 41(5): 1587-1595.
[4]	MA Rui-jie, LI Xin. The Finite Element Volume Method and Application of Mathematic Model of Karst Groundwater Flow [J]. J4, 2005, 35(06): 762-0765.