A parallel magnetotelluric (MT) 3D forward modeling algorithm was implemented on multi-core CPU, and was introduced into the 3D conjugate gradient (CG) inversion algorithm with MT data for parallel computation. The authors tested the parallel MT 3D forward modeling code with a theoretical model and the parallel 3D inversion code with one synthetic data. A comparison of the execution efficiency was made between serial computation and parallel computation. The results indicate that the parallel 3D forward modeling and inversion computation are correct, and the efficiency of computation is greatly improved.
[1] Newman G A,Alumbaugh D L.Three-Dimension Electromagnetic modeling and inversion on massive parallel computer[R].Sandia Report,1996:SAND96-0582,UC-400,403,405.[2] Newman G A,Alumbaugh D L.Three-Dimension massive parallel Electro- magnetic Inversion-I.Theory[J].Geophys.J.Int,1997:345-354.[3] Michael Commer,Newman G A.Large scale 3D EM inversion using optimized simulation grids nonconformal to the model space[R].SEG 2006,PAPR247.[4] Michael Commer,Newman G A.3D CSEM modeling and inversion for hydrocarbon reservoir mapping The bathymetry problem[R].SEG 2006:247.[5] Carazzone J J, Dickens T A, Green K E,et al. Inversion Study of a Large Marine CSEM Survey[R].SEG Las Vegas 2008 Annual Meeting,2008:644-647.[6] Yoshioka K, Zhdanov M S. Electromagnetic forward modeling based on the integral equation method using parallel computers[C]//75th Annual International Meeting, SEG, Expanded Abstracts, 2005:550-553.[7] Ken Yoshioka, Zhdanov M S.Modeling large-scale geoelectrical structures with inhomogeneous backgrounds using the integral equation method: applications to the bathymetry e.ect study in marine CSEM data[R].SEG 2006,PAPR211.[8] Endo M, Cuma M, Zhdanov M S.Multi-grid IE method for large-scale models with inhomogeneous background[R].SEG 2007, seg-annual,2007.[9] Zhdanov M S, Cuma M.Three-dimensional electromagnetic holographic imaging in offshore petroleum exploration[C]//SEG Las Vegas 2008 Annual Meeting,2008:629-633.[10] Tan H D,Tong T,Lin C H.The parallel 3D magnetotelluric forward modeling algorithm[J].Applied Geophysics,2006,3(4).[11] Lin C H,Tan H D,Tong T.Parallel rapid relaxation inversion of 3D magnetotelluric data[J].Applied Geophysics,2009,6(1).[12] Mitsuhata Yuji, Uchida Toshihiro.3D magnetotelluric modeling using the T-Ω finite-element method[J].Geophysics, 2004,69(1) :108-119.[13] Sadd Y. Iterative methods for sparse linear systems[R].Society of Inustrial and Applied Mathematics,1996.[14] Newman G A,Alumbaugh D L.Three-dimensional magnetotelluric inversion using non-linear conjugate gradients[J].Geophys J Int, 2000:140, 410-424.[15] Mackie R L, Madden T R.Three-dimensional magnetotelluric inversion using conjugate gradients[J].Geophysical J Int, 1993b:115, 215-229.[16] Zhdanov M S,Varentov I M,Weaver J T,et al.Methods for modeling electromagnetic fields: Results from COMMEMI-The international project on the comparison of modeling methods for electromagnetic induction[J].Journal of Applied Geophysics, 1997,37:133-271.