Boundary element solutions for broad-band 3-D geo-electromagnetic problems accelerated by an adaptive multilevel fast multipole method

Geophysical Journal International

Volumn 192, Issue 2, 2013, Pages 473-499

  Ren, Zhengyong ^a   Kalscheuer, Thomas ^a   Greenhalgh, Stewart ^a   Maurer, Hansruedi ^a  

^a ETH ZURICH   (Switzerland)

Author keywords

Electromagnetic theory; Magnetic and electrical properties.; Magnetotelluric; Numerical solutions

Indexed keywords

BOUNDARY ELEMENT METHOD; ELECTRIC FIELDS; ELECTROMAGNETIC FIELDS; ELECTROMAGNETIC INDUCTION; ELECTROMAGNETS; GEOMETRY; MAGNETOTELLURICS; NUMERICAL METHODS; NUMERICAL MODELS; SAILING VESSELS; TOPOGRAPHY;

ELECTROMAGNETIC THEORIES; HIGH FREQUENCY SIMULATION; MAGNETIC AND ELECTRICAL PROPERTIES; MAGNETOTELLURIC; MULTILEVEL FAST MULTIPOLE METHOD; NUMERICAL SOLUTION; PERMITTIVITY DISTRIBUTIONS; SYSTEM OF LINEAR EQUATIONS;

FINITE ELEMENT METHOD;

BOUNDARY ELEMENT METHOD; COMPLEXITY; ELECTROMAGNETIC FIELD; FINITE ELEMENT METHOD; MAGNETOTELLURIC METHOD; NUMERICAL MODEL; THREE-DIMENSIONAL MODELING; TOPOGRAPHY;

NORWAY;

EID: 84875772060     PISSN: 0956540X     EISSN: 1365246X     Source Type: Journal    
DOI: 10.1093/gji/ggs043     Document Type: Article

References (78)

1. Abramowitz, M. & Stegun, I., 1964. Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, Dover publications, New York.
2. Ahmed, M., Lavers, J. & Burke, P., 1988. A comparision of boundary integral formulations for 3D axisymmetric eddy-current problems, IEEE Trans. Magn., 24, 138-141.
3. Amestoy, P., Duff, I., L'Excellent, J.&Koster, J., 2002.Afully asynchronous multifrontal solver using distributed dynamic scheduling, SIAM J. Matrix Anal. Appl., 23, 15-41.
4. Amos, D., 1986. Algorithm 644: a portable package for Bessel functions of a complex argument and nonnegative order, ACM Trans. Math. Softw., 12, 265-273.
5. Aprea, C., Booker, J. & Smith, J., 1997. The forward problem of electromagnetic induction: accurate finite-difference approximations for twodimensional discrete boundaries with arbitrary geometry, Geophys. J. Int., 129, 29-40.
6. Avdeev, D., Kuvshinov, A., Pankratov, O. & Newman, G., 2002. Threedimensional induction logging problems, Part I: an integral equation solution and model comparisons, Geophysics, 67, 413-426.
7. Badea, E., Everett, M., Newman, G. & Biro, O., 2001. Finite-element analysis of controlled-source electromagnetic induction using Coulombgauged potentials, Geophysics, 66, 786-799.
8. Balay, S. et al., 2010. PETSc Users Manual, Tech. rep., Argonne National Laboratory.
9. Bapat, M. & Liu, Y., 2010. A new adaptive algorithm for the fast multipole boundary element method, Comput. Model. Eng. Sci., 58, 161-184.
10. Berdichevsky, M.N. & Dmitriev, V.I., 2008. Models and Methods of Magnetotellurics, Springer, Berlin and Heidelberg.
11. Bleszynski, E., Bleszynski, M. & Jaroszewicz, T., 1996. Aim: adaptive integral method for solving large-scale electromagnetic scattering and radiation problems, Radio Sci., 31, 1225-1252.
12. Blome, M., Maurer, H.R. & Schmidt, K., 2009. Advances in threedimensional geoelectric forward solver techniques, Geophys. J. Int., 176, 740-752.
13. Borup, D. & Gandhi, O., 1984. Fast-Fourier-transform method for calculation of SAR distributions in finely discretized inhomogeneous models of biological bodies, IEEE Trans. Microw. Theory Tech., 32, 355-360.
14. Buchau, A., Huber, C., Rieger, W. & Rucker, W., 2000. Fast BEM computations with the adaptive multilevel fast multipole method, IEEE Trans. Magn., 36, 680-684.
15. Çakir, Ö ., 2006. Themultilevel fastmultipolemethod for forwardmodelling the multiply scattered seismic surface waves, Geophys. J. Int., 167(2), 663-678.
16. Çakir, Ö ., 2009. Forward modelling the multiply scattered 2.5-D teleseismic P waves accelerated by the multilevel fast multipole method, Geophys. J. Int., 176(2), 505-517.
17. Candansayar, M.E. & Tezkan, B., 2008. Two-dimensional joint inversion of radiomagnetotelluric and direct current resistivity data, Geophys. Prospect., 56, 737-749.
18. Colton, D. & Kress, R., 1983. Integral Equation Methods in Scattering Theory, John Wiley & Sons, New York.
19. Doherty, J., 1988. EM modelling using surface integral equations, Geophys. Prospect., 36, 644-668.
20. Epton, M. & Dembart, B., 1995. Multipole translation theory for the 3- dimensional Laplace and Helmholtz equations, SIAM J. Sci. Comput., 16, 865-897.
21. Franke, A., B̈orner, R.-U. & Spitzer, K., 2007. Adaptive unstructured grid finite element simulation of two-dimensional magnetotelluric fields for arbitrary surface and seafloor topography, Geophys. J. Int., 171, 71-86.
22. Fujiwara, H., 2000. The fast multipole method for solving integral equations of three-dimensional topography and basin problems, Geophys. J. Int., 140(1), 198-210.
23. Galassi, M., Gough, B., Jungman, G., Theiler, J., Davies, J., Booth, M. & Rossi, F., 2009. GNU Scientific Library Reference Manual, Network Theory Ltd., United Kingdom.
24. Gumerov, N. & Duraiswami, R., 2004a. Fast Multipole Methods for the Helmholtz Equation in Three Dimensions, Elsevier Science, Amsterdam.
25. Gumerov, N. & Duraiswami, R., 2004b. Recursions for the computation of multipole translation and rotation coefficients for the 3-D Helmholtz equation, SIAM J. Sci. Comput., 25, 1344-1381.
26. Haber, E., Ascher, U., Aruliah, D. & Oldenburg, D., 2000. Fast simulation of 3D electromagnetic problems using potentials, J. Comput. Phys., 163, 150-171.
27. Hackbusch, W. & Nowak, Z., 1989. On the fast matrix multiplication in the boundary element method by panel clustering, Numer. Math., 54, 463-491.
28. Harrington, R., 1993. Field Computation by Moment Methods, Wiley-IEEE Press, New York.
29. Hohmann, G., 1975. Three-dimensional induced polarization and electromagnetic modeling, Geophysics, 40, 309-324.
30. Homentcovschi, D., 2002. Minimum-order regular boundary integral equations for three-dimensional eddy-current problem, IEEE Trans. Magn., 38, 3433-3438.
31. Hou, J., Mallan, R.K. & Torres-Verdin, C., 2006. Finite-difference simulation of borehole EMmeasurements in 3D anisotropicmedia using coupled scalar-vector potentials, Geophysics, 71, G225-G233.
32. Ismail, N. & Pedersen, L., 2011. The electrical conductivity distribution of the Hallandsas Horst, Sweden: a controlled source radiomagnetotelluric study, Near Surf. Geophys., 9, 45-54.
33. Jin, J., 2002. The Finite Element Method in Electromagnetics, Wiley-IEEE Press, New York.
34. Jiracek, G., 1990. Near-surface and topographic distortions in electromagnetic induction, Surv. Geophys., 11(2), 163-203.
35. Kalscheuer, T., Bastani, M., Donohue, S., Persson, L., Pfaffhuber, A.A., Reiser, F. & Ren, Z., 2012. Delineation of a quick clay structure at Smørgrav, Norway, with electromagnetic methods, J. appl. Geophys., in revision.
36. Kalscheuer, T., Pedersen, L.B. & Siripunvaraporn, W., 2008. Radiomagnetotelluric two-dimensional forward and inverse modelling accounting for displacement currents, Geophys. J. Int., 175, 486-514.
37. Kalscheuer, T., Juanatey, G., Meqbel, N. &Pedersen, L.B., 2010. Non-linear model error and resolution properties from two-dimensional single and joint inversions of direct current resistivity and radiomagnetotelluric data, Geophys. J. Int., 182, 1174-1188.
38. Key, K. & Weiss, C., 2006. Adaptive finite-element modeling using unstructured grids: the 2D magnetotelluric example, Geophysics, 71, G291- G299.
39. Lee, J.,Zhang, J.&Lu, C., 2004. Sparse inverse preconditioning of multilevel fast multipole algorithm for hybrid integral equations in electromagnetics, IEEE Trans. Antennas Propag., 52, 2277-2287.
40. Li, Y.&Key, K., 2007. 2D marine controlled-source electromagnetic modeling, Part 1: an adaptive finite-element algorithm, Geophysics, 72,WA51- WA62.
41. Liu, E. & Lamontagne, Y., 1998. Geophysical application of a new surface integral equation method for EM modeling, Geophysics, 63, 411-423.
42. Liu, G. & Becker, A., 1992. Evaluation of terrain effects in AEM surveys using the boundary element method, Geophysics, 57, 272-278.
43. Liu, Y., 2009. Fast Multipole Boundary Element Method-Theory and Applications in Engineering, Cambridge University Press, Cambridge.
44. Mackie, R., Smith, J. & Madden, T., 1994. 3D electromagnetic modeling using finite-difference equations- the magnetotelluric example, Radio Sci., 29, 923-935.
45. May, D. & Knepley, M., 2011. Optimal, scalable forward models for computing gravity anomalies, Geophys. J. Int., 187, 161-177.
46. Mayergoyz, I., 1982. Boundary integral equation of minimum order for the calculation of three-dimensional eddy current problems, IEEE Trans. Magn., 18, 536-539.
47. Mitsuhata, Y. & Uchida, T., 2004. 3D magnetotelluric modeling using the T-Omega finite-element method, Geophysics, 69, 108-119.
48. Mukherjee, S. & Everett, M.E., 2011. 3D controlled-source electromagnetic edge-based finite element modeling of conductive and permeable heterogeneities, Geophysics, 76, F215-F226.
49. Nam, M., Kim, H., Song, Y., Lee, T., Son, J. & Suh, J., 2007. 3D magnetotelluric modelling including surface topography, Geophys. Prospect., 55, 277-287.
50. Newman, G. & Alumbaugh, D., 2002. Three-dimensional induction logging problems, Part 2: a finite-difference solution, Geophysics, 67, 484-491.
51. Operto, S., Virieux, J., Amestoy, P., L'Excellent, J.-Y., Giraud, L. & Ali, H. B.H., 2007. 3d finite-difference frequency-domain modeling of viscoacoustic wave propagation using a massively parallel direct solver: a feasibility study, Geophysics, 72(5), SM195-SM211.
52. Parry, J. &Ward, S., 1971. Electromagnetic scattering from cylinders of arbitrary cross section in a conductive half space, Geophysics, 36, 67-100.
53. Pedersen, L. & Engels, M., 2005. Routine 2D inversion of magnetotelluric data using the determinant of the impedance tensor, Geophysics, 70, G33-G41.
54. Pedersen, L., Bastani, M. & Dynesius, L., 2005. Groundwater exploration using combined controlled-source and radiomagnetotelluric techniques, Geophysics, 70, G8-G15.
55. Pedersen, L.B., Bastani, M. & Dynesius, L., 2006. Some characteristics of the electromagnetic field from radio transmitters in Europe, Geophysics, 71, G279-G284.
56. Pérez-Gavilán, J.&Aliabadi, M., 2001. Symmetric Galerkin BEMformulticonnected bodies, Commun. Numer. Methods Eng., 17, 761-770.
57. Pusch, D. & Ostrowski, J., 2010. Robust FEM/BEM coupling for magnetostatics on multiconnected domains, IEEE Trans. Magn., 46, 3177-3180.
58. Ren, Z., Kalscheuer, T. & Maurer, H., 2010. Accurate and robust evaluation of singular boundary element integrals for computing 3D radiomagnetotelluric fields, in Proceedings of the Expanded Abstracts EAGE Near Surface Meeting Zurich, EAGE, Switzerland.
59. Rokhlin, V., 1992. Diagonal forms of translation operators for Helmholtz equation in three dimensions, Tech. rep., DTIC Document.
60. Rucker, W., Hoschek, R.&Richter, K., 1995.VariousBEMformulations for calculating eddy currents in terms of field variables, IEEE Trans. Magn., 31, 1336-1341.
61. Seher, T. & Tezkan, B., 2007. Radiomagnetotelluric and Direct Current Resistivity measurements for the characterization of conducting soils, J. appl. Geophys., 63, 35-45.
62. Shen, J., 1995. Computational Electromagnetics Using Boundary Elements: Advances in Modelling Eddy Currents, Computational Mechanics Publications, Southampton.
63. Shewchuk, J., 1996. Triangle: engineering a 2D quality mesh generator and Delaunay triangulator, Lecture Notes Comput. Sci., 1148, 203-222.
64. Streich, R., 2009. 3D finite-difference frequency-domain modeling of controlled-source electromagnetic data: direct solution and optimization for high accuracy, Geophysics, 74(5), F95-F105.
65. Taskinen, M. & Yl̈a-Oijala, P., 2006. Current and charge integral equation formulation, IEEE Trans. Antennas Propag., 54, 58-67.
66. Tezkan, B. & Saraev, A., 2008. A new broadband radiomagnetotelluric instrument: applications to near surface investigations, Near Surf. Geophys., 6, 245-252.
67. Tezkan, B., Goldman, M., Greinwald, S., Ḧordt, A., M̈uller, I., Neubauer, F. &Zacher, G., 1996. A joint application of radiomagnetotellurics and transient electromagnetics to the investigation of a waste deposit in Cologne (Germany), J. appl. Geophys., 34, 199-212.
68. Tezkan, B., Georgescu, P. & Fauzi, U., 2005. A radiomagnetotelluric survey on an oil-contaminated area near the Brazi Refinery, Romania, Geophys. Prospect., 53, 311-323.
69. Tsuboi, H. & Misaki, T., 1987. Three-dimensional analysis of eddy-current distributions by the boundary element method using vector variables, IEEE Trans. Magn., 23, 3044-3046.
70. Ward, S. & Hohmann, G., 1987. Electromagnetic theory for geophysical applications: electromagnetic methods in applied geophysics, Tulsa: Soc. Explor. Geophys., Nabighian MN (Ed.), 1, 131-311.
71. Weidelt, P., 1975. Electromagnetic induction in three-dimensional structures, J. Geophys., 41, 85-109.
72. Xu, S., Ruan, B., Zhou, H., Chen, L. & Xu, S., 1997. Numerical modeling of 3-D terrain effect on MT field, Sci. China Ser. D-Earth Sci., 40, 269-275.
73. Yl̈a-Oijala, P., 2008. Numerical analysis of combined field integral equation formulations for electromagnetic scattering by dielectric and composite objects, Prog. Electromag. Res., 3, 19-43.
74. Yl̈a-Oijala, P. & Taskinen, M., 2003. Calculation of CFIE impedance matrix elements with RWG and n x RWG functions, IEEE Trans. Antennas Propag., 51, 1837-1846.
75. Yl̈a-Oijala, P. & Taskinen, M., 2007. Improving conditioning of electromagnetic surface integral equations using normalized field quantities, IEEE Trans. Antennas Propag., 55, 178-185.
76. Zhao, J. & Chew, W., 2000. Integral equation solution of Maxwell's equations from zero frequency to microwave frequencies, IEEE Trans. Antennas Propag., 48(10), 1635-1645.
77. Zhdanov, M.S., Lee, S.K. & Yoshioka, K., 2006. Integral equation method for 3D modeling of electromagnetic fields in complex structures with inhomogeneous background conductivity, Geophysics, 71, G333-G345.
78. Zheng, D., 1997. Three-dimensional eddy current analysis by the boundary element method, IEEE Trans. Magn., 33, 1354-1357.