Full-Text Search:
Home|Journal Papers|About CNKI|User Service|FAQ|Contact Us|中文
《Chinese Journal of Geophysics》 2015-08
Add to Favorite Get Latest Update

3-D frequency-domain CSEM modeling using aparallel direct solver

HAN Bo;HU Xiang-Yun;HUANG Yi-Fan;PENG Rong-Hua;LI Jian-Hui;CAI Jian-Chao;Hubei Subsurface Multi-scale Imaging Key Laboratory,Institute of Geophysics and Geomatics,China University of Geosciences;  
Three-dimensional modeling of electromagnetic data is a computationally demanding problem.For frequently-used numerical techniques such as finite-element and finite-difference methods,solving the large linear systems arising from the discretization of Maxwell's equations is a key step which has a major impact on the applicability of the solution,and it has always been a research topic to solve the linear equations efficiently,robustly and accurately.A 3D modeling scheme based on direct solutions of the linear system is presented for frequency domain controlled-source electromagnetic(CSEM)surveys. The Helmholtz equation in terms of secondary electric fields is discretized using a finite-volume(FV)method over a staggered grid.Taking advantage of recent developments in numerical algorithms and the availability of computational resources,the resulting linear system of FV equations is solved directly using the massively parallel solver,namely MUMPS,instead of the most commonly used linear solvers,i.e.Krylov subspace iterative techniques.The direct solver carries out an LU(and possibly LDLT)decomposition of the system matrix and then computes solutions efficiently by applying forward and backward substitutions.To evaluate the computational performance of the direct solver,a series of numerical tests based on synthetic 1D models were conducted,and the results indicate that(1)Normalized residuals of solutions are almost independent of the conductivity value assigned to air layers but increase rapidly as the frequency value decreases.Nevertheless,the order of magnitude of the largest normalized residual is as small as 10-11.At the same time,although the matrix factorization time varies as either the air conductivity or the frequency changes,the variation is only a fraction of the total run time.(2)Both the execution time and required memory increase rapidly(more than linearly)with increasing grid sizes.(3)By executing MUMPS in parallel over multiple processors,not only the total run time but also the average memory used per processor can be reduced a lot.However,the total memory requirement increases with the number of processes.The scalability of MUMPS is limited.Additional numerical experiments considering specific survey settings were done to demonstrate the reliability and effectiveness of the code,and the results are as follows:(1)The FV numerical solutions show excellent agreement with semianalytic solutions for the 1D models.(2)The computation time of a multitransmitter problem is comparable to that of a single-transmitter problem.(3)Reasonable modeling results of 3D models can be obtained for both typical land and marine survey scenarios.In summary,compared with iterative linear solvers,the direct solvers generally benefit CSEM modeling in three aspects.The first is they often provide more accurate solutions.The second is that direct solvers are much more stable for ill-conditioned linear systems,which are almost inevitable because of large electrical conductivity contrasts and/or non-uniform grid.The last is in multitransmitter problems,only a single matrix factorization is necessary,and multiple solutions can be achieved very easily by reusing the factors.The presented 3D CSEM modeling scheme,which employs the MUMPS direct solver,possesses all these advantages.In addition,solving linear systems can be executed in parallel to speed up the computation and to reduce the average memory used per node,although the parallel scalability of MUMPS is limited.In spite of the fact that matrix factorizations for large models can entail tremendous computational cost,it can be anticipated that direct solvers will be used more and more widely as the development of both numerical algorithms and computers.
【Fund】: 国家自然科学基金项目(41274077;41474055);; 中国地质调查局项目(12120113101800)联合资助
【CateGory Index】: P631.325
Download(CAJ format) Download(PDF format)
CAJViewer7.0 supports all the CNKI file formats; AdobeReader only supports the PDF format.
【Citations】
Chinese Journal Full-text Database 9 Hits
1 WU Gui-ju,HU Xiang-yun,LIU Hui (Institute of Geophysics and Geometics,China University of Geosciences,Wuhan 430074,China);Progress in CSAMT three-dimensional forward numerical simulation[J];Progress in Geophysics;2010-05
2 Basin and Reservoir Research Center,Petroleum University, Beijing 102249,China;MODELING OF 3-D ELECTROMAGNETIC RESPONSES IN FREQUENCY DOMAIN BY USING STAGGERED-GRID FINITE DIFERENCE METHOD SHEN JINSONG[J];Chinese Journal of Geophysics;2003-02
3 WANG Ruo,DI Qing-Yun,WANG Miao-Yue,WANG Guang-Jie Institute of Geology and Geophysics,Chinese Academy of Sciences,Beijing 100029,China;Research on the effect of 3D body between transmitter and receivers on CSAMT response using Integral Equation method[J];Chinese Journal of Geophysics;2009-06
4 CHEN Gui-Bo~(1,2),WANG Hong-Nian~2,YAO Jing-Jin~3,YANG Shou-Wen~21 School of Physics,Changchun University of Science and Technology,Changchun 130022,China2 School of Physics,J ilin University,Changchun 130026,China3 Institute of Geophysical and Geochemical Exploration,Geological Academe of China,Hebei Lang fang 065000,China;Modeling of electromagnetic responses of 3-D electrical anomalous body in a layered anisotropic earth using integral equations[J];Chinese Journal of Geophysics;2009-08
5 XU Zhi-Feng,WU Xiao-Ping~* Mengcheng National Geophysical Observatory,School of Earth and Space Science,University of Science and Technology of China,Hefei 230026,China;Controlled source electromagnetic 3-D modeling in frequency domain by finite element method[J];Chinese Journal of Geophysics;2010-08
6 YANG Bo1,XU Yi-Xian1,2,HE Zhan-Xiang3,SUN Wei-Bin3 1 Institute of Geophysics and Geomatics,China University of Geosciences,Wuhan 430074,China 2 State Key Laboratory of Geological Processes and Mineral Resources,Wuhan 430074,China 3 BGP,CNPC,Hebei Zhuozhou 072751,China;3D frequency-domain modeling of marine controlled source electromagnetic responses with topography using finite volume method[J];Chinese Journal of Geophysics;2012-04
7 HAN Bo;HU Xiang-Yun;Adam SCHULTZ;ZUO Bo-Xin;LI Jian-Hui;CAI Jian-Chao;Institute of Geophysics and Geomatics,China University of Geosciences;College of Earth,Ocean,and Atmospheric Sciences,Oregon State University;;Three-dimensional forward modeling of the marine controlled-source electromagnetic field with complex source geometries[J];Chinese Journal of Geophysics;2015-03
8 Tan Han dong,Wei Wen bo,Deng Ming and Jin Sheng. Tan Han dong.Collage of Geophysics and Information Technique,China University of Geology,Beijing City,100083,China;General use formula in MT tensor impedance.[J];Oil Geophysical Prospecting;2004-01
9 Han Bo1,Hu Xiangyun1,He Zhan-xiang2 and Cai Jianchao1.1.Institute of Geophysics and Geomatics,China University of Geosciences(Wuhan),Wuhan,Hubei 430074,China 2.GME and Geochemical Survey Department,BGP Inc.,CNPC,Zhuozhou,Hebei 072751,China;Mathematical classification of magnetotelluric inversion methods[J];Oil Geophysical Prospecting;2012-01
【Co-citations】
Chinese Journal Full-text Database 10 Hits
1 Liu Changsheng~(12) Ren Zhengyong~(13) Tang Jingtian~1 Yan Yan~4 1.School of Info-physics and Geometrics Engineering,Central South University,Changsha 410083,China. 2.Changsha Aeronautical Vocational and Technical College,Changsha 410124,China 3.Institute of Geophysics,ETH-Swiss Federal Institute of Technology,Zurich 8092,Switzerland. 4.Institute of Higher Education,Central South University,Changsha 410083,China.;Three-dimensional magnetotellurics modeling using edgebased finite-element unstructured meshes[J];应用地球物理(英文版);2008-03
2 Lin Chang-Hong1,2,3, Tan Han-Dong1,2,3, and Tong Tuo1,2,31. State Key Laboratory of Geological Processes and Mineral Resources, China University of Geosciences, Beijing, 100083, China. 2. Key Laboratory of Geo-detection (China University of Geosciences), Ministry of Education, Beijing, 100083, China. 3. School of Geophysics and Information Technology, China University of Geosciences, Beijing, 100083, China.;Three-dimensional conjugate gradient inversion of magnetotelluric full information data[J];应用地球物理(英文版);2011-01
3 Wang Xindong Shen Jinsong.;Sensitivity Analysis and Imaging Technique of the 3D Crosswell Electromagnetic Measurements[J];Well Logging Technology;2003-05
4 SHEN Jin-song (Basin & Reservoir Research Centre, University of Petroluem, Beijing 102200, China);Modeling of the 3-D Electromagnetic Responses to the Anisotropic Medium by the Edge Finite Element Method[J];Well Logging Technology;2004-01
5 QIN Shi-yin 1,WANG Chang-xue 2,YANG Wei 2,TAO Guo 2,SHEN Jin-song 2 (1.China University of Geosciences, Beijing 100083,China; 2.School of Earth Resources & Information, University of Petroleum,Beijing 102249,China);The Computation of Electromagnetic Scattering Responses in Anisotropic Formation and Its Applications[J];Well Logging Technology;2004-03
6 Ling Mingyou(Linpan Oil Production Plant,Shengli Oilfield Ltd.Co.,SINOPEC,Shandong 251507,P.R.China);The Study of Remaining Oil by the Technic of Well-earth Electric Potential Image[J];Fault-block Oil & Gas Field;2003-04
7 ZHOU Lixin1,MI Suting2,HUANG Xiaobo1 (1.Air Force Engineering, Design & Research Institute Beijing 100077, China; 2.Center for analysis and prediction, CSB, Beijing 100036, China);AFE of sliding and bending deformation occurred in Daliushu Dam, Heishanxia Gorge of the Yellow River[J];Progress In Geophysics;2003-03
8 CHEN Shijun1,2, LIU Hong1, ZHOU Jianyu2,HE Xinghua3 (1.Institute of Geology and Geophysics ,Chinese Academy of Sciences, Beijing 100029, China; 2. Geophysical Research Institute of Shengli Oilfield Company Limited, SINOPEC, Dongying 257022 China;3. Shengli Oilfield Company Limited, SINOPEC, Dongying 257000, China);The current situation and prospect of crosswell seismic technique[J];Progress In Geophysics;2003-03
9 ZHOU Jianyu1,2, LI Jianping2,ZHANG Yonghui2 (1. Lanzhou Geology Institute, Chinese Academy of Sciences, Lanzhou 730000, China; 2. Shengli Oilfield, SINOPEC, Dongying 257000, China);The analysis and application of crosswell seismic survey[J];Progress In Geophysics;2003-03
10 SHEN Jin-song (Basin and Reservoir Research Center,Petroleum University, Beijing 102249,China);Modeling of the multi-component induction log in anisotropic medium by using finite difference method[J];Progress In Geophysics;2004-01
China Proceedings of conference Full-text Database 6 Hits
1 XUE Guo-Qiang~1,LI Xiu~2,DI Qing-Yun~1 (1.Institute of Geology and Geophysics,Chinese Academy of science,Beijing,100029,China; 2.School of Geology and Survey Engineering,Changan University,Xi'an 710054,China);Research progress in TEM forward modeling and inversion calculation[A];[C];2009
2 ZHAO Guo-ze, CHEN Xiao-bin, TNAG Ji (Institute of Geology,China Earthquake Administration, Beijing 100029,China);Advanced geo-electromagnetic methods in China[A];[C];2007
3 TANG Jing-tian, REN Zheng-yong, HUA Xi-rui(School of Info-physics and Geomatics Engineering, Changsha, 410083);The forward modeling and inversion in geophysical electromagnetic field[A];[C];2007
4 LI Xiu;XUE Guo-qiang;WANG Miao-yue;School of Geology and Survey Engineering,Chang,An University;Institute of Geology and Geophysics,Chinese Academy of Sciences;;Research progress of the transformation from TEM to pseudo-wave[A];[C];2011
5 WANG Ruo;WANG Miao-yue;DI Qing-yun;WANG Guang-jie;Institute of Geology and Geophysics,Chinese Academy of Sciences;;3D1C CSAMT modeling using finite element method[A];[C];2015
6 WANG Xian-xiang;DI Qing-yun;TANG Jing;GUO Hai-feng;Key Lab of Engineering Geomechanics,Institute of Geology and Geophysics,Chinese Academy of Sciences;Hebei Electric Power Design and Research Institute;;CSAMT resolution in 3D geoelectric structures[A];[C];2015
【Secondary Citations】
Chinese Journal Full-text Database 10 Hits
1 HUANG Lin-ping, DAI Shi-kun (Resources and Information Institute, University of Petroleum, Beijing102249, China);Finite Element Calculation Method of 3D Electromagnetic Field under Complex Condition[J];;2002-06
2 YANG Hui 1 DAI Shi kun 1 SONG Hai bin 2 Huang Lin ping 1 (1. Key Laboratory of Geophysics, University of Petroleum, Beijing 102200, China; 2. Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100101, China);Overview of Joint Inversion of Integrated Geophysics[J];Progress In Geophysics;2002-02
3 ZHOU Hou-fang, LIU Chuang,SHI Kun-fa(1. Institute of Geographical Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China; 2. The Graduate School of Chinese Academy of Sciences, Beijing 100039, China; 3. Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100101, China);A Review of study on geothermal resources exploration[J];Progress In Geophysics;2003-04
4 DI Qing-yun1,MARTYN Unsworth2,WANG Miao-yue1(1. Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100101,China;2. Department of Physics, University of Alberta, Edmonton Alberta Canaela T6G 2J1);2.5 D CSAMT modeling with finite element method[J];Progress In Geophysics;2004-02
5 SU Yong-jun,WANG Xu-ben,LUO Jian-qun(Chengdu University Of Technology,Chengdu 610059,China);The ahaeological application of high-density resistivity method to ditch exploration on Sanxingdui Site[J];Progress in Geophysics;2007-01
6 WANG Ruo1,WANG Miao-yue1,LU Yuan-lin2(1.Institute Geology of and Geophysics,Chinese Academy of Sciences,Beijing 100029,China2.China Geo-Equipment corporation,Beijing 100102,China);Preliminary study on 3D3C CSAMT method modeling using finite element method[J];Progress in Geophysics;2007-02
7 TANG Jing-tian,REN Zheng-yong,HUA Xi-rui(School of Info-physics and Geomatics Engineering,Changsha,410083);The forward modeling and inversion in geophysical electromagnetic field[J];Progress in Geophysics;2007-04
8 LIU Jian-Xin,TONG Xiao-Zhong,YANG Xiao-Hong,XIE Wei HU Hou-Ji(School of Info-physics and Geomatics Engineer,Central South University,Changsha 410083,China);Application of real coded genetic algorithm in two-dimensional magnetotelluric inversion[J];Progress in Geophysics;2008-06
9 YU LI-MING XU HONG-HAI (Beijing Graduate Department, Wuhan College of Geology);THE NUMERICAL SOLUTION FOR THE EM RESPONSE OF THREE-DIMENTIONAL CONDUCTING PERMEABLE BODIES[J];Chinese Journal of Geophysics;1986-02
10 CHEN JIU-PING CHEN LE-SHOU WANG GUANG-E(Department of Geophysics, China University of Geoscicrices. Beijing);MAGNETOTELLURIC MODELLING OF THREE-DIMENSIONAL BODIES IN LAYERED EARTHS[J];Chinese Journal of Geophysics;1990-04
©2006 Tsinghua Tongfang Knowledge Network Technology Co., Ltd.(Beijing)(TTKN) All rights reserved