3-D induction log modelling with integral equation method and domain decomposition pre-conditioning

3-D induction log modelling with integral equation method and domain decomposition pre-conditioning

Saputera, D. H. (University of Bergen)
Jakobsen, M. (University of Bergen)
van Dongen, K.W.A. (TU Delft ImPhys/Medical Imaging; TU Delft ImPhys/Van Dongen goup)
Jahani, N. (NORCE Norwegian Research Centre AS)
Eikrem, K. S. (NORCE Norwegian Research Centre AS)
Alyaev, S. (NORCE Norwegian Research Centre AS)

2024

The deployment of electromagnetic (EM) induction tools while drilling is one of the standard routines for assisting the geosteering decision-making process. The conductivity distribution obtained through the inversion of the EM induction log can provide important information about the geological structure around the borehole. To image the 3-D geological structure in the subsurface, 3-D inversion of the EM induction log is required. Because the inversion process is mainly dependent on forward modelling, the use of a fast and accurate forward modelling tool is essential. In this paper, we present an improved version of the integral equation (IE) based modelling technique for general anisotropic media with domain decomposition pre- conditioning. The discretized IE after domain decomposition equals a fixed-point equation that is solv ed iterativ ely with either the block Gauss-Seidel or Jacobi pre-conditioning. Within each iteration, the inverse of the block matrix is computed using a Krylov subspace method instead of a direct solver. An additional reduction in computational time is obtained by using an adaptive relative residual stopping criterion in the iterativ e solv er. Using this domain decomposition scheme, numerical experiments show computation time reductions by factors of 1.97-2.84 compared to solving the full-domain IE with a GMRES solver and a contraction IE pre- conditioner. Additionally, the reduction of memory requirement for covering a large area of the induction tool sensitivity enables acceleration with limited GPU memory. Hence, we conclude that the domain decomposition method is improving the efficiency of the IE method by reducing the computation time and memory requirement.

Electromagnetic theory
Numerical modelling
Numerical solutions

http://resolver.tudelft.nl/uuid:8455eb64-fe10-4196-8c60-14f16c292381

https://doi.org/10.1093/gji/ggad454

0956-540X

Geophysical Journal International, 236 (2), 834-848

Institutional Repository

journal article

© 2024 D. H. Saputera, M. Jakobsen, K.W.A. van Dongen, N. Jahani, K. S. Eikrem, S. Alyaev