1887

Abstract

Summary

Real-time 3D imaging of the induction log is essential for improving the decision-making process in geosteering. To fulfill this need, we investigated various strategies for reducing the computational cost for 3D modelling of induction logs using the integral equation (IE) method, including the use of iterative Krylov solver, convolution with FFT algorithm, contraction IE formulation, computation acceleration with GPUs, and domain decomposition. We present two cases example to demonstrate the implementation of IE with these strategies. In the first case, we show that the application of domain decomposition allows one to only discretize the inhomogeneous domain and save the computation cost in the case of isolated domains. We present a logging while drilling scenario on a complex model for the second case. Our implementation of the efficient IE on GPUs enables significant acceleration and allows the computation of 3D forward modelling within less than two minutes for each local 3D simulation domain with approximately two million grid cells on a laptop. The implementation of domain decomposition formulation shows a different arrangement of solving IE by decomposing the domain.

Loading

Article metrics loading...

/content/papers/10.3997/2214-4609.202310505
2023-06-05
2026-01-14
Loading full text...

Full text loading...

References

  1. Dyatlov, G., Onegova, E., and Dashevsky, Y. [2015]. Efficient 2.5D electromagnetic modeling using boundary integral equations. GEOPHYSICS, 80(3), E163–E173.
    [Google Scholar]
  2. Endo, M., Čuma, M., and Zhdanov, M. S. [2009]. Large-scale Electromagnetic Modeling for Multiple Inhomogeneous Domains. Communications in Computational Physics, 6(2), 269–289.
    [Google Scholar]
  3. Fang, S., Gao, G., and Verdin, C. T. [2006]. Efficient 3D Electromagnetic Modelling in the Presence of Anisotropic Conductive Media, Using Integral Equations. Exploration Geophysics, 37(3), 239–244.
    [Google Scholar]
  4. Hursán, G., and Zhdanov, M. S. [2002]. Contraction integral equation method in three-dimensional electromagnetic modeling. Radio Science, 37(6), 1-1–1-13.
    [Google Scholar]
  5. Jakobsen, M., and Tveit, S. [2018]. Distorted Born iterative T-matrix method for inversion of CSEM data in anisotropic media. Geophysical Journal International, 214(3), 1524–1537.
    [Google Scholar]
  6. Puzyrev, V., Pardo, D., Calo, V., and Torres-Verdín, C. [2019]. Recent advances on the inversion of deep directional borehole resistivity measurements. ASEG Extended Abstracts, 2019(1), 1–3.
    [Google Scholar]
  7. SaputeraD.H., JakobsenM., JahaniN., AlyaevS., EikremK.S., van=DongenK.W.A. [2022]. GPU-accelerated integral equation method for 3D modelling of induction logs. [Poster]. Geosteering and Formation Evaluation Workshop by NORCE, Stavanger, Norway, 1–2 November 2022. https://nfes.org/2022/10/14/gpu-accelerated-workshop-abstract-23.html
    [Google Scholar]
  8. van der Vorst, H. A. [1992]. Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems. SIAM Journal on Scientific and Statistical Computing, 13(2), 631–644.
    [Google Scholar]
  9. Wilson, G., Marchant, D., Haber, E., Clegg, N., Zurcher, D., Rawsthorne, L., and Kunnas, J. [2019]. Real-Time 3D Inversion of Ultra-Deep Resistivity Logging-While-Drilling Data. SPE Annual Technical Conference and Exhibition. OnePetro.
    [Google Scholar]
  10. Zhdanov, M. S. [2009]. Geophysical electromagnetic theory and methods. Elsevier.
    [Google Scholar]
/content/papers/10.3997/2214-4609.202310505
Loading
/content/papers/10.3997/2214-4609.202310505
Loading

Data & Media loading...

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error