Direct Method Solution of 3-D Magnetotelluric Modeling Using Vector Finite Element Method
Keywords:direct solver, electromagnetic, magnetotelluric modeling, rectilinear mesh, vector finite element
As exploration is forced into more difficult areas with complex three-dimensional (3-D) structural environments, the importance of 3-D magnetotelluric (MT) modeling is essential for the correct interpretation of MT data. To reduce the complexity of the calculations introduced by 3-D models, iterative forward calculation of MT response functions is used as basis for inversion of 3-D MT data. This paper proposes an alternative procedure for making reliable 3-D MT modeling codes for forward calculation that is not only effective but also accurate. This is accomplished by using a direct method to solve the linear systems arising from the discretization process in the vector finite element approach. The vector finite element method is known for its capability of overcoming difficulties in modeling caused by possible jumps of normal components across discontinuities of material properties. Meanwhile, by using a direct method rather than an iterative method, the process of solving the linear equations does not suffer from slow convergence. Here, we present a comparison between our modeling codes and codes based on a different approach. In the resulting 3-D MT responses it was found that the proposed method has high accuracy.
Tietze, K. & Ritter. O., 3D Magnetotelluric Inversion in Practice - The Electrical Conductivity Structure of the San Andreas Fault in Central California, Geophysical Journal International, 195(1), pp. 130-147, 2013.
Wei. W., Pape, F.L., Jones, A.G., Vozar, J., Dong, H., Unsworth, M.J, Jin, S., Ye, G., Jing, J., Zhang, L. & Xie, C., Northward Channel Flow in Northern Tibet Revealed From 3D Magnetotelluric Modeling, Physics of the Earth and Planetary Interiors, 235, pp. 13-24, 2014.
Ting, S.C. & Hohmann, G.W., Integral Equation Modeling of Three-Dimensional Magnetotelluric Response, Geophysics, 46(2), pp. 182-197, 1981.
Wannamaker, P.E., Advances in Three-Dimensional Magnetotelluric Modeling using Integral Equations, Geophysics, 56(11), pp. 1716-1728, 1991.
Mackie, R.L., Madden, T.R. & Wannamaker, P.E., Three-Dimensional Magnetotelluric Modeling using Difference Equations - Theory and Comparisons to Integral Equations Solutions, Geophysics, 58(2), pp. 215-226, 1993.
Smith, J.T., Conservative Modeling of 3-D Electromagnetic Fields Part I: Properties and Error Analysis, Geophysics, 61(5), pp. 1308-1318, 1996.
Smith, J.T., Conservative Modeling Of 3-D Electromagnetic Fields Part II: Biconjugate Gradient Solution and an Accelerator, Geophysics, 61(5), pp. 1319-1324, 1996.
Fomenko, Y.E. & Mogi, T., A New Computation Method for A Staggered Grid of 3D EM Field Conservative Modeling, Earth Planets Space, 54(5), pp. 499-509, 2002.
Livelybrooks, D., Program 3Dfeem: A Multidimensional Electromagnetic Finite Element Model, Geophysical Journal International, 144, pp. 443-458, 1993.
Mogi, T., Three-Dimensional Modeling of Magnetotelluric Data using Finite Element Method, Journal of Applied Geophysics, 35(2), pp. 185-189, 1996.
Coggon, J.H., Electromagnetic and Electrical Modeling by the Finite Element Method, Geophysics, 36(1), pp. 132-155, 1971.
Nam, M.J., Kim, H.J., Song, Y., Lee, T.J., Son, J.S. & Suh, J.H., 3D Magnetotelluric Modeling Including Surface Topography, Geophysical Prospecting, 55(2), pp. 277-287, 2007.
Ren, Z., Kalscheuer, T., Greenhalgh, S., & Maurer, H., A Goal-Oriented Adaptive Finite-Element Approach for Plane Wave 3-D Electromagnetic Modeling, Geophysical Journal International, 194(2), pp. 700-718, 2013.
Brner, R.U., Numerical Modeling in Geo-Electromagnetics: Advances and Challenges, Surveys in Geophysics, 31(2), pp. 225-245, 2010.
Schwarzbach, C., Brner, R.U. & Spitzer, K., Three-Dimensional Adaptive Higher Order Finite Element Simulation for Geo-Electromagnetics - A Marine CSEM Example, Geophysical Journal International, 187(1), pp. 63-74, 2011.
Ansari, S. & Farquharson, C.G., 3D Finite-Element forward Modeling of Electromagnetic Data using Vector and Scalar Potentials and Unstructured Grids, Geophysics, 79(4), pp. E149-E165, 2014.
Mitsuhata, Y. & Uchida, T., 3D Magnetotelluric Modeling using the T- Finite-Element Method, Geophysics, 69(1), pp. 108-119, 2004.
Farquharson, C.G. & Miensopust, M.P., Three-Dimensional Finite-Element Modeling of Magnetotelluric Data with A Divergence Correction, Journal of Applied Geophysics, 75(4), pp. 699-710, 2011.
Siripunvaraporn, W., Egbert, G. & Lenbury, Y., Numerical Accuracy of Magnetotelluric Modeling: A Comparison of Finite Difference Approximations, Earth Planets Space, 54(6), pp. 721-725, 2002.
Jin, J.M., The Finite Element Method in Electromagnetics, 3rd edition, Wiley-IEEE Press, 2014.
Aruliah, D.A., Ascher, U.M., Haber, E. & Oldenburg, D., A Method for the Forward Modeling of 3D Electromagnetic Quasi-Static Problems, Mathematical Models and Methods in Applied Sciences, 11(1), pp. 1-21, 2001.
Gould, N.I.M., Scott, J.A. & Hu, Y., A Numerical Evaluation of Sparse Direct Solvers for the Solution of Large Sparse Symmetric Linear Systems Of Equations, ACM Transactions on Mathematical Software, 33(2), pp. 1-32, 2007.
Schenk, O. & Gartner, K., Solving Unsymmetric Sparse Systems of Linear Equations with PARDISO, Future Generation Computer Systems, 20(3), pp. 475-487, 2004.
Miensopust, M.P, Queralt, P., Jones, A.G. & the 3D MT modellers, Magnetotelluric 3-D Inversion - A Review of Two Successful Workshops on Forward and Inversion Code Testing and Comparison, Geophysical Journal International, 193(1), pp. 1216-1238, 2013.