Efficient pre-conditioned iterative solution strategies for the electromagnetic diffusion in the Earth: finite-element frequency-domain approach
Open Access
- 12 March 2013
- journal article
- Published by Oxford University Press (OUP) in Geophysical Journal International
- Vol. 193 (3), 1460-1473
- https://doi.org/10.1093/gji/ggt071
Abstract
We formulate a 3-D finite-element frequency-domain (FEFD) solution for electromagnetic (EM) diffusion and present efficient solution strategies. A system of FEFD equations is pre-conditioned by incomplete LU (ILU) and subsequently solved by the quasi-minimal residual (QMR) method. A rule of thumb for choosing an effective drop tolerance of ILU is proposed. When multiple sources are simulated in a given model, ILU is computed only once and is reused as a pre-conditioner for multiple QMR computations with different source vectors. Resulting solution vectors are also bootstrapped to reduce the number of QMR iterations required for the convergence. We demonstrate that when conductivity structures of an earth model and source frequencies are updated/perturbed, ILU that is computed from the previous model is still an effective and useful pre-conditioner for new forward modelling problems. Using the reusability of ILU, we also propose a new efficient way to overcome the slow convergence rate of the iterative FEFD solution in the static limit. We show that the reuse of ILU and solution bootstrapping serve as effective strategies for improving the computational efficiency of the iterative FEFD solution. Finally, we apply the proposed efficient solution strategies to marine EM survey scenarios in complex offshore models and further demonstrate their effectiveness.Keywords
This publication has 35 references indexed in Scilit:
- Three-dimensional finite-element modelling of magnetotelluric data with a divergence correctionJournal of Applied Geophysics, 2011
- Three‐dimensional induction logging problems, Part 2: A finite‐difference solutionGeophysics, 2002
- Interpretation of long-offset transient electromagnetic data from the Odenwald area, Germany, using two-dimensional modellingGeophysical Journal International, 2000
- Developments and trends in the parallel solution of linear systemsParallel Computing, 1999
- A fast vector-potential method using tangentially continuous vector finite elementsIEEE Transactions on Microwave Theory and Techniques, 1998
- A Fast and High Quality Multilevel Scheme for Partitioning Irregular GraphsSIAM Journal on Scientific Computing, 1998
- Conservative modeling of 3-D electromagnetic fields, Part II: Biconjugate gradient solution and an acceleratorGeophysics, 1996
- Frequency‐domain modelling of airborne electromagnetic responses using staggered finite differences1Geophysical Prospecting, 1995
- Conjugate Gradient-Type Methods for Linear Systems with Complex Symmetric Coefficient MatricesSIAM Journal on Scientific and Statistical Computing, 1992
- Mixed finite elements in ?3Numerische Mathematik, 1980