An approximate 3D computational method for real-time computation of induction logging responses

Over many years, induction logging systems have been used to create well formation logs. The major drawback for the utilization of these tools is the long simulation time for a single forward computation. We proposed an efficient computational method based on a contrast-type of integral-equation formulation, in which we applied an approximation for the 3D electromagnetic field. We assumed that the dominant contribution in the integral equation is obtained by the contribution around the singularity of Green’s kernel. It is expected that the approximation yields reliable results when the (homogeneous) background conductivity around the logging tool is close to the actual conductivity at the location of the tool. We have developed a data-driven method to determine this background conductivity from the dominant part of the measured coaxial magnetic fields, which are mainly influenced by the conductivity at the tool sensors. For a synthetic model, the results were compared to the ones of a rigorous solution of the integral equation and show a good simulation response to small-scale variations in the medium. Further, the method was used to simulate the response of a realistic reservoir model. Such a model is created by a geological modeling program. We concluded that our approximate method was able to improve the approximation results in highly heterogeneous structures compared to the Born approximation and provide an effective medium-gradient around the tool. Our method, based on the wavefield approximation, also estimates the error, and hence yields a warning when the method becomes unreliable.

[1]  Sofia Davydycheva,et al.  Triaxial induction tool with electrode sleeve : FD modeling in 3D geometries , 2009 .

[2]  Gregory A. Newman,et al.  Electromagnetic induction in a fully 3‐D anisotropic earth , 2000 .

[3]  Ki Ha Lee,et al.  Electromagnetic fields in a non-uniform steel-cased borehole , 2006 .

[4]  S. A. Petersen,et al.  Optimization Strategy for Shared Earth Modeling , 2004 .

[5]  A. A. Kaufman,et al.  Principles of Induction Logging , 2013 .

[6]  Gregory A. Newman,et al.  Three-dimensional induction logging problems, Part 2: A finite-difference solution , 2002 .

[7]  G. Deschamps,et al.  Singularity in Green's function and its numerical evaluation , 1980 .

[8]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[9]  Ki Ha Lee,et al.  3D interpretation of electromagnetic data using a modified extended Born approximation , 2003 .

[10]  J. Bladel,et al.  Electromagnetic Fields , 1985 .

[11]  Jean G. Van Bladel,et al.  Electromagnetic Fields (IEEE Press Series on Electromagnetic Wave Theory) , 2007 .

[12]  Misac N. Nabighian,et al.  Electromagnetic Methods in Applied Geophysics , 1988 .

[13]  Qing Huo Liu,et al.  Applications of the BCGS-FFT method to 3-D induction well logging problems , 2003, IEEE Trans. Geosci. Remote. Sens..

[14]  Sergei Knizhnik,et al.  3D integral equation modeling with a linear dependence on dimensions , 2009 .

[15]  V. Dmitriev,et al.  Quasi-analytical approximations and series in electromagnetic modeling , 2000 .

[16]  Sofia Davydycheva,et al.  An efficient finite‐difference scheme for electromagnetic logging in 3D anisotropic inhomogeneous media , 2003 .

[17]  Gregory A. Newman,et al.  Three‐dimensional induction logging problems, Part I: An integral equation solution and model comparisons , 2002 .

[18]  David Pardo,et al.  Simulation of borehole-eccentered triaxial induction measurements using a Fourier hp finite-element method , 2013 .

[19]  A. Kaufman,et al.  Chapter 12 - Mathematical Modeling of the Response of Induction Logging Tools in 3D Geometries , 2003 .

[20]  B. I. Anderson,et al.  Modeling and Inversion Methods for the Interpretation of Resistivity Logging Tool Response , 2001 .

[21]  Weng Cho Chew,et al.  Effect of tool eccentricity on some electrical well-logging tools , 1990 .

[22]  A. Abubakar,et al.  A Green function formulation of the Extended Born approximation for three-dimensional electromagnetic modelling , 2005 .

[23]  Richard Liu,et al.  A fast integral equation solver for 3D induction well logging in formations with large conductivity contrasts , 2013 .