Finite Element Resistivity Forward Modeling Using Algebraic Multigrid Preconditioned Conjugate Gradient Method

Abstract:

An algebraic multigrid by smoothed aggregation preconditioned conjugate gradient method is developed to solve the liner system arising from 3-D direct current finite element resistivity forward modeling. The algorithm combines the efficiency of algebraic multigrid method and the stability of conjugate gradient method. Algebraic multigrid by smoothed aggregation keep in high-efficiency while simulation using local quasi-uniform mesh and its convergence effect will reduce while numerical modeling using anisotropic stretched grids. However tensor product non-equidistant mesh, a kind of anisotropic stretched grids, is often used in 3-D direct current resistivity forward modeling. In order to improve this situation, a factor is added to guide correct aggregation. Consequently, a typical example is used to prove that the improvement is the right. Finally, it is natural to conclude that the algorithm suggested in this paper is efficient and robust whether simulation using local quasi-uniform mesh or tensor product non-equidistant mesh

• Periodical: Applied Mechanics and Materials (Volumes 249-250)
• Pages: 792-797
• DOI: https://doi.org/10.4028/www.scientific.net/AMM.249-250.792
• Online since: December 2012
• Authors: Gui Hong Zou, Hua Qing Liang
• Keywords: 3D Resistivity Forward, Algebraic Multigrid Method, Finite Element Method (FEM), Preconditioned Conjugate Gradient Method, Smooothed Aggregation

References

[1] X. Wu,G. Xu, S. Li, The Caculation of Three-dimensional Geoeletric Field of Point Source by incomplete Cholesky conjugate gradient, Chinese Journal of Geophysics, 41(1998) 848-854.

[2] X. Wu, G. Xu, S. Li, A 3-D Finite Element Resistivity Forward Modeling Using Conjugate Gradient Algorithm, Chinese Journal of Geophysics, 46(2003) 428-432.

DOI: 10.1002/cjg2.3379

[3] J. LU, X. Wu, K. Spitzer, Algebraic Multigrid Method for 3D dc Resistivity modeling, Chinese Journal of Geophysics, 53(2010) 700-707.

[4] J. Tang, F. Wang, Y. Ren, 3-D direct current resistivity forward modeling by adaptive multigrid finite element method, Journal of Central South University of Technology, 17(2010) 587-592.

DOI: 10.1007/s11771-010-0527-z

[5] W.L. Briggs, S.F. McCormick, A multigrid tutorial, second ed., Tsinghua Press, Beijing, (2000).

[6] A. Dey, H.F. Morrison, Resistivity modeling for arbitrarily shaped three-dimensional structures, Geophysics, 44(1979a) 753-780.

DOI: 10.1190/1.1440975

[7] S. Xu, The finite element method in geophysics, People's Education Press, Beijing, (1995).

[8] S.K. Zhao, M.J. Yedlin, Some refinements on the finite-difference method for 3-D dc resistivity modeling, Geophysics, 61(1996) 1301-1307.

DOI: 10.1190/1.1444053

[9] T. Mifune, T. Iwashita, and M. Shimasaki, New Algebraic Multigrid Preconditioning for Iterative Solvers in Electromagnetic Finite Edge-Element Analyses, IEEE TRANSACTIONS ON MAGNETICS, 39(2003)1677-1680.

DOI: 10.1109/tmag.2003.810350

[10] P. Vanek, J. Mandel, M. Brezina, Algebraic multigrid based on smoothed aggregation for second and fourth order problems. Computing 56(1996) 179–196.

DOI: 10.1007/bf02238511

[11] M. W. Gee, J.J. Hu, R.S. Tuminaro, A new smoothed aggregation multigrid method for anisotropic problems, Numerical Linear Algebra with Applications, 16(2009) 19–37.

DOI: 10.1002/nla.593