Finite‐difference solution of the three‐dimensional electromagnetic problem using divergence‐free preconditioners

Author(s):  
Mike Zaslavsky ◽  
Sofia Davydycheva ◽  
Vladimir Druskin ◽  
Aria Abubakar ◽  
Tarek Habashy ◽  
...  
Geophysics ◽  
2002 ◽  
Vol 67 (2) ◽  
pp. 484-491 ◽  
Author(s):  
Gregory A. Newman ◽  
David L. Alumbaugh

A 3‐D finite‐difference solution is implemented for simulating induction log responses in the quasi‐static limit that include the wellbore and bedding that exhibits transverse anisotropy. The finite‐difference code uses a staggered grid to approximate a vector equation for the electric field. The resulting linear system of equations is solved to a predetermined error level using iterative Krylov subspace methods. To accelerate the solution at low induction numbers (LINs), a new preconditioner is developed. This new preconditioner splits the electric field into curl‐free and divergence‐free projections, which allows for the construction of an approximate inverse operator. Test examples show up to an order of magnitude increase in speed compared to a simple Jacobi preconditioner. Comparisons with analytical and mode matching solutions demonstrate the accuracy of the algorithm.


Sign in / Sign up

Export Citation Format

Share Document