Electromagnetic induction in a fully 3‐D anisotropic earth

The bulk electrical anisotropy of sedimentary formations is a macroscopic phenomenon which can result from the presence of porosity variations, laminated shaly sands, and water saturation. Accounting for its effect on induction log responses is an ongoing research problem for the well-logging community since these types of sedimentary structures have long been correlated with productive hydrocarbon reservoirs such as the Jurassic Norphlet Sandstone and Permian Rotliegendes Sandstone. Presented here is a staggered-grid finite-difference method for simulating electromagnetic (EM) induction in a fully 3-D anisotropic medium. The electrical conductivity of the formation is represented as a full 3£ 3 tensor whose elements can vary arbitrarily with position throughout the formation. To demonstrate the validity of this approach, finite-difference results are compared against analytic and quasi-analytic solutions for tractable 1-D and 3-D model geometries. As a final example, we simulate 2C‐40 induction tool responses in a crossbedded aeolian sandstone to illustrate the magnitude of the challenge faced by interpreters when electrical anisotropy is neglected.

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

[2]  Roland W. Freund,et al.  Conjugate Gradient-Type Methods for Linear Systems with Complex Symmetric Coefficient Matrices , 1992, SIAM J. Sci. Comput..

[3]  J. Klein,et al.  The Petrophysics Of Electrically Anisotropic Reservoirs , 1995 .

[4]  Roland W. Freund,et al.  A Transpose-Free Quasi-Minimal Residual Algorithm for Non-Hermitian Linear Systems , 1993, SIAM J. Sci. Comput..

[5]  J. H. Moran,et al.  Basic theory of induction logging and application to study of two-coil sondes , 1962 .

[6]  D. L. Alumbaugh,et al.  Three-dimensional electromagnetic modeling and inversion on massively parallel computers , 1996 .

[7]  Mark E. Everett,et al.  Electric dipole fields over an anisotropic seafloor: theory and application to the structure of 40 ma pacific ocean lithosphere , 1999 .

[8]  Otto N. Fanini,et al.  Increased Oil-In-Place In Low-Resistivity Reservoirs From Multicomponent Induction Log Data , 2000 .

[9]  J. H. Moran,et al.  SOME EFFECTS OF FORMATION ANISOTROPY ON RESISTIVITY MEASUREMENTS IN BOREHOLES , 1958 .

[10]  G. W. Hohmann,et al.  A finite-difference, time-domain solution for three-dimensional electromagnetic modeling , 1993 .

[11]  Brian R. Spies,et al.  Three-Dimensional Electromagnetics , 1999 .

[12]  J. H. Moran,et al.  Effects of formation anisotropy on resistivity‐logging measurements , 1979 .

[13]  K. Yee.  Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media , 1966 .

[14]  Liang C. Shen,et al.  Simulation Of Induction And Mwd Resistivity Tools In Anisotropic Dipping Beds , 1997 .

[15]  Bülent Baygün,et al.  Oilfield Anisotropy : Its Origins and Electrical Characteristics , 1995 .

[16]  G. Newman,et al.  Frequency‐domain modelling of airborne electromagnetic responses using staggered finite differences , 1995 .

[17]  D. Georgi,et al.  Aspects of Multicomponent Resistivity Data and Macroscopic Resistivity Anisotropy , 2001 .

[18]  S. Constable,et al.  The electrical conductivity of an isotropic olivine mantle , 1992 .

[19]  J. Klein,et al.  Induction Log Anisotropy Corrections , 1993 .

[20]  T. D. Barber,et al.  The Effect Of Crossbedding Anisotropy On Induction Tool Response , 2001 .

[21]  L. Onsager.  Reciprocal Relations in Irreversible Processes. II. , 1931 .