3-D direct current resistivity forward modeling by adaptive multigrid finite element method

Abstract

Based on the fact that 3-D model discretization by artificial could not always be successfully implemented especially for large-scaled problems when high accuracy and efficiency were required, a new adaptive multigrid finite element method was proposed. In this algorithm, a-posteriori error estimator was employed to generate adaptively refined mesh on a given initial mesh. On these iterative meshes, V-cycle based multigrid method was adopted to fast solve each linear equation with each initial iterative term interpolated from last mesh. With this error estimator, the unknowns were nearly optimally distributed on the final mesh which guaranteed the accuracy. The numerical results show that the multigrid solver is faster and more stable compared with ICCG solver. Meanwhile, the numerical results obtained from the final model discretization approximate the analytical solutions with maximal relative errors less than 1%, which remarkably validates this algorithm.

Keywords

adaptive multigrid / a-posteriori error estimator / unstructured mesh / V-cycle / finite element method

References

[1] LiY. G., SpitzerK.. Three-dimensional DC resistivity forward modelling using finite elements in comparison with finite-difference solutions [J]. Geophys J Int, 2002, 151: 924-934

[2] ZhouB., GreenhalghS. A.. Finite element three-dimensional direct current resistivity modelling: Accuracy and efficiency considerations [J]. Geophys J Int, 2001, 145: 679-688

[3] WuX.-ping.. A 3-D finite-element algorithm for DC resistivity modelling using the shifted incomplete Cholesky conjugate gradient method [J]. Geophys J Int, 2003, 154: 947-956

[4] ShenJ.-s., SunW.-bo.. 2.5-D modeling of cross-hole electromagnetic measurement by finite element method [J]. Pet Sci, 2008, 5: 126-134

[5] XiongB., LuoY.-zhong.. Finite element modeling of 2.5-D TEM with blocks homogeneous conductivity [J]. Chinese Journal Geophysics, 2006, 49(2): 590-597

[6] FengT., GullikssonM., LiuW. B.. Adaptive finite element methods for the identification of elastic constants [J]. Journal of Scientific Computing, 2006, 26(2): 217-235

[7] MollerM., KuzminD.. Adaptive mesh refinement for high-resolution finite element schemes [J]. Int J Numer Meth Fluids, 2006, 52: 545-569

[8] KerryK., ChesterW.. Adaptive finite-element modeling using unstructured grids: The 2D magnetotelluric example [J]. Geophysics, 2006, 71(6): G291-G299

[9] LiY. G., KerryK.. 2D marine controlled-source electromagnetic modeling: Part 1. An adaptive finite-element algorithm [J]. Geophysics, 2007, 72(2): WA51-WA62

[10] StephanP., MaischakM., LeydeckerF.. An hp-adaptive finite element/boundary element coupling method for electromagnetic problems [J]. Comput Mech, 2007, 39: 673-680

[11] TangJ.-t., RenZ.-y., HuaX.-rui.. Theoretical analysis of geo-electromagnetic modeling on Coulomb gauged potentials by adaptive finite element method [J]. Chinese Journal Geophysics, 2007, 50(5): 1584-1597

[12] WilliamL. B.A multigrid tutorial [M], 20002nd edPhiladelphia, Society for Industrial and Applied Mathematics

[13] ChenZ. M., WangL., ZhengW. Y.. An adaptive multigrid method for time-harmonic Maxwell equations with singularities [J]. SIAM Journal Science Computation, 2007, 20(1): 118-139

[14] MulderW. A.. A multigrid solver for 3D electromagnetic diffusion [J]. Geophysical Prospecting, 2006, 54: 633-649

[15] ZienkiewiczO. C., ZhuJ. Z.. Super-convergent patch recovery and a posteriori error estimates: Part 1. The recovery technique [J]. International Journal in Numerical Method Engineering, 1992, 33(7): 1331-1364

[16] ZienkiewiczO. C., ZhuJ. Z.. Super-convergent patch recovery and a posteriori error estimates: Part 2. Error estimates and adaptivity [J]. International Journal in Numerical Method Engineering, 1992, 33(7): 1365-1382

[17] ZienkiewiczO. C., TaylorR. L.The finite-element method [M], 2000, Butterworth-Heinemann, Elsevier

[18] JinJ.The finite element method in electromagnetics [M], 2001, New York, John Wiley & Sons

[19] RenZ.-y., TangJ.-tian.. Object-oriented philosophy in designing adaptive finite-element package for 3D elliptic deferential equations [J]. EOS on Transactions, 2007, 88(52): NS11A-0164

[20] [20] HYPRETEAM. Hypre: A high performance preconditioner [EB/OL]. [2007-01-04]. http://www.llnl.gov/CASC/hypre.

[21] [21] PARAVIEWTEAM. Paraview: An open-source multi-platform parallel visualization application [EB/OL]. [2008-06-09]. http://www.paraview.org/New/index.html.

[22] RichardB., MichaelB., TonyF. C., DemmelJ., DonatoJ. M., DongarraJ., EijkhoutV., PozoR., RomineC., van der VorstH.Templates for the solution of linear systems, building blocks for iterative methods [M], 2006, New York, Society for Industrial and Applied Mathematics