1887
3D Electromagnetics
  • ISSN: 2202-0586
  • E-ISSN:

Abstract

We present a novel technique to simulate numerically the measurements performed by a borehole induction-logging tool in 3D anisotropic rock formations. The simulations are based on an integral equation formulation. Previously, such a formulation was considered impractical for solving large-scale problems due to the resulting large full matrix. To avoid this difficulty, we assume a uniform background model and make use of a uniform grid whereupon there is no need to construct explicitly all of the entries of the full Green’s function matrix.

Using a uniform background model, the entries of the corresponding electric and magnetic Green’s tensors are relatively easy to calculate. In the presence of a uniform grid (not necessarily cubic), it is only necessary to calculate the first row of the resulting electric Green’s function matrix. Further, since the matrix is block Toeplitz, it can be rewritten into a block circulant form, and therefore matrix-vector multiplication can be efficiently performed with two FFTs and one inverse FFT. This strategy reduces the computation cost from 0(N*N) to 0(N*logN). In addition to the substantial computer savings, the FFT technique also substantially reduces memory storage requirements because only the first row and the first column in the block Toeplitz matrix are needed to perform the computations of the remaining entries of the matrix.

Numerical simulations of the measurements performed with an induction tool in dipping and anisotropic rock formations are benchmarked against accurate 3D finite-difference and ID codes. These benchmark exercises show that the newly developed integral-equation algorithm produces accurate and efficient simulations for a variety of borehole and formation conditions.

Loading

Article metrics loading...

/content/journals/10.1071/ASEG2003_3DEMab006
2003-04-01
2026-01-21
Loading full text...

Full text loading...

References

  1. Avdeev, D. B., Kuvshinov, A. V., Pankratov, O. V., and Newman, G. A., 2002, Three-dimensional induction logging problems, Part 1: An integral equation solution and model comparisons, Geophysics, Vol. 67, No. 2, pp.413-426.
  2. Born, M., 1933, Optics: Springer-Verlag, New York.
  3. Catedra, M.F., Torres, R.P., Basterrechea, J., and Gago, E., 1995, The CG-FFT method: Application of signal processing techniques to electromagnetics. Boston: Artech House.
  4. Fang S., and Wang, T., 2000, Accurate Born simulation of induction response using an optimal background, SEG Expanded Abstract, Calgary, pp. 1806-1809.
  5. Gan, H. and Chew, W. C, 1994, A discrete BCG-FFT algorithm for solving 3D inhomogeneous scatterer problems, Electromagnetic Laboratory, Department of Eletrical and Computer Engineering, University of Illinois, Urbana, IL, 61801.
  6. Gao, G-Z., Fang, S., and Torres-Verdin, C, 2003, A new approximation for 3D electromagnetic scattering in the presence of anisotropic conductive media, 3DEM-III symposium, Feb. 20-21, Adelaide, Australia.
  7. Gianzero, S., 1999, The paradox of anisotropy revisited, The log analyst, Vol. 40, No. 6, pp.485-491
  8. Golub, G. H. and Van Loan, C. F., 1996, Matrix Computations, Third edition, The Johns Hopkins University Press.
  9. Habashy, T.M., Groom, R.W. and Spies, B., 1993, Beyond the Born and Rytov approximations: A nonlinear approach to electromagnetic scattering, J. Geophys. Res., 98, pp. 1759-1775.
  10. Hohmann, G.W., 1983, Three-dimensional EM modelling, Geophys. Surv., 6, pp.27-53.
  11. Hursan, G., and M. S. Zhdanov, 2002, Contraction integral equation method in three-dimensional electromagnetic modeling: Radio Science, Vol. 37, No 6. pp.1-13.
  12. Klein, J.D., Martin, P.R., and Allen, D.F., 1997, the petrophysics of electrically anisotropic reservoirs: The Log Analyst, May-June, Vol. 38, No. 3, pp.25-36.
  13. Kriegshauser B., Fanini, O., Yu, L., 2000, Advanced inversion techniques for multicomponent induction data: Presented at the 70th Ann. Mtg., Soc. Expl. Geophys.
  14. Liu, Q.H., Zhang, Z.Q., Xu, X.M., 2001, The hybrid extended Born approximation and CG-FFT method for electromagnetic induction problems, IEEE Transactions on Geoscience and Remote sensing, Vol. 39, No. 2 .
  15. Moran, J.H. and Gianzero, S., 1979, Effects of formation anisotropy on resistivity-logging measurements: Geophysics, Vol. 44, No. 7, pp.1266-1286.
  16. Newman, G. A. and Alumbaugh D. L., 2002, Three-dimensional induction logging problems, Part 2: A finite-difference solution for 3-D induction logging problems, Geophysics, Vol. 67, No. 2, pp.484-491.
  17. Singer, B. Sh., 1995, Method for solution of Maxwell’s equations in non-uniform media, Geophys. J. Int., Vol 120, pp.590-598.
  18. Singer, B. Sh., Mezzatesta, A., and Wang, T., 2002, 3D modelling of electromagnetic field in subsurface and borehole environments, 72nd SEG Annual Meeting, Expanded Abstracts, Soc. Of Explor. Geophys., Salt Lake City, Utah, pp.312-315.
  19. Sleijpen, Gerard L.G. and Fokkema, Diederik R., 1993, BiCGSTAB(l) for Linear Equations Involving Unsymmetric Matrices with Complex Spectrum, Electronic Transactions on Numerical Analysis, Vol. 1, pp.11-32.
  20. Wang T. and Fang S., 2001, 3-D electromagnetic anisotropy modeling using finite differences, Geophysics Vol. 66, No.5 pp.1386-1398.
  21. Wannamaker, P.E., Hohmann, G.W., SanFilipo, W.A., 1983, Electromagnetic modelling of three-dimensional bodies in layered earths using integral equations, Geophysics, Vol. 49, No. 1, pp. 60-74.
  22. Xiong, Z., 1992, Electromagnetic modeling of 3-Dstructures by the method of system iteration using integral equations: Geophysics,Vol. 57, No. 12, pp.1556-1561.
  23. Xiong, Z., and Tripp, A. C, 1995, Electromagnetic scattering of large structures in a layered earth using integral equations: Radio Sci., 30, pp.921-929.
  24. Zhdanov M. S. and Fang S., 1996, Quasi-linear approximation in 3-D electromagnetic modeling: Geophysics, Vol. 61, No. 3, pp.646-665.
  25. Zhdanov M. S. and Fang S., 1997, Quasi-linear series in three-dimensional electromagnetic modeling: Radio Science, Vol. 32, No. 6, pp.2167-2188.
/content/journals/10.1071/ASEG2003_3DEMab006
Loading
  • Article Type: Research Article
Keyword(s): anisotropy; electromagnetic; induction tool; integral equation; numerical simulation
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