3D marine magnetotelluric modeling and inversion with the finite-difference time-domain method

ABSTRACTFrequency-domain methods, which are typically applied to 3D magnetotelluric (MT) modeling, require solving a system of linear equations for every frequency of interest. This is memory and computationally intensive. We developed a finite-difference time-domain algorithm to perform 3D MT modeling in a marine environment in which Maxwell’s equations are solved in a so-called fictitious-wave domain. Boundary conditions are efficiently treated via convolutional perfectly matched layers, for which we evaluated optimized parameter values obtained by testing over a large number of models. In comparison to the typically applied frequency-domain methods, two advantages of the finite-difference time-domain method are (1) that it is an explicit, low-memory method that entirely avoids the solution of systems of linear equations and (2) that it allows the computation of the electromagnetic field unknowns at all frequencies of interest in a single simulation. We derive a design criterion for vertical node spacin...

[1]  E. C. Titchmarsh.  Introduction to the Theory of Fourier Integrals , 1938 .

[2]  Yongwimon Lenbury,et al.  Numerical accuracy of magnetotelluric modeling: A comparison of finite difference approximations , 2002 .

[3]  J. T. Smith.  Conservative modeling of 3-D electromagnetic fields, Part I: Properties and error analysis , 1996 .

[4]  G. Newman,et al.  Three-dimensional magnetotelluric inversion using non-linear conjugate gradients , 2000 .

[5]  D. Avdeev,et al.  Detecting a salt dome overhang with magnetotellurics: 3D inversion methodology and synthetic model studies , 2012 .

[6]  S. Constable,et al.  Marine magnetotellurics for petroleum exploration Part I: A sea-floor equipment system , 1998 .

[7]  Y. Sasaki.  Full 3-D inversion of electromagnetic data on PC , 2001 .

[8]  Colin Farquharson,et al.  Three-dimensional finite-element modelling of magnetotelluric data with a divergence correction , 2011 .

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

[10]  Kerry Key,et al.  Marine Electromagnetic Studies of Seafloor Resources and Tectonics , 2011, Surveys in Geophysics.

[11]  N. Georgieva,et al.  Problem-independent enhancement of PML ABC for the FDTD method , 2003 .

[12]  Michael S. Zhdanov,et al.  Contraction integral equation method in three‐dimensional electromagnetic modeling , 2002 .

[13]  Catherine de Groot-Hedlin,et al.  Finite-difference modeling of magnetotelluric fields: Error estimates for uniform and nonuniform grids , 2006 .

[14]  L. Knizhnerman,et al.  Spectral approach to solving three-dimensional Maxwell's diffusion equations in the time and frequency domains , 1994 .

[15]  Douglas W. Oldenburg,et al.  An Algorithm For the Three-dimensional Inversion of Magnetotelluric Data , 2002 .

[16]  D. R. Fokkema,et al.  BICGSTAB( L ) FOR LINEAR EQUATIONS INVOLVING UNSYMMETRIC MATRICES WITH COMPLEX , 1993 .

[17]  A general correspondence principle for time-domain electromagnetic wave and diffusion fields , 1996 .

[18]  Weng Cho Chew,et al.  A 3D perfectly matched medium from modified maxwell's equations with stretched coordinates , 1994 .

[19]  Bruce Archambeault,et al.  The Finite-Difference Time-Domain Method , 1998 .

[20]  A. Orange,et al.  Magnetotelluric exploration for hydrocarbons , 1989, Proc. IEEE.

[21]  K. Vozoff,et al.  8. The Magnetotelluric Method , 1991 .

[22]  J. Carcione.  A spectral numerical method for electromagnetic diffusion , 2006 .

[23]  Rolf Schuhmann,et al.  A general stability analysis of FIT/FDTD applied to lossy dielectrics and lumped elements , 2004 .

[24]  Irene A. Stegun,et al.  Handbook of Mathematical Functions. , 1966 .

[25]  Alan G. Jones,et al.  Magnetotelluric 3-D inversion—a review of two successful workshops on forward and inversion code testing and comparison , 2013 .

[26]  Toshihiro Uchida,et al.  3D magnetotelluric modeling using the T‐Ω finite‐element method , 2004 .

[27]  G. Michael Hoversten,et al.  Marine magnetotellurics for petroleum exploration, Part II: Numerical analysis of subsalt resolution , 1998 .

[28]  Chester J. Weiss,et al.  Mapping 3D salt using the 2D marine magnetotelluric method: Case study from Gemini Prospect, Gulf of Mexico , 2006 .

[29]  Olivier Cessenat,et al.  Sophie, an FDTD code on the way to multicore, getting rid of the memory bandwidth bottleneck better using cache , 2013, ArXiv.

[30]  Ralph-Uwe Börner,et al.  Numerical Modelling in Geo-Electromagnetics: Advances and Challenges , 2010 .

[31]  R. L. Mackie,et al.  Three-dimensional magnetotelluric modelling and inversion , 1989, Proc. IEEE.

[32]  Jorge Nocedal,et al.  A Limited Memory Algorithm for Bound Constrained Optimization , 1995, SIAM J. Sci. Comput..

[33]  J. J. Zach,et al.  Gradient Calculations for 3D Inversion of CSEM Data Using a Fast Finite-difference Time-domain Modelling Code , 2008 .

[34]  Chester J. Weiss,et al.  Mapping thin resistors and hydrocarbons with marine EM methods, Part II -Modeling and analysis in 3D , 2006 .

[35]  Bengt Fornberg,et al.  The pseudospectral method; accurate representation of interfaces in elastic wave calculations , 1988 .

[36]  P. Wannamaker,et al.  Three-dimensional magnetotelluric modeling using difference equations­ Theory and comparisons to integral equation solutions , 1993 .

[37]  Frank A. Maaø,et al.  Fast finite-difference time-domain modeling for marine-subsurface electromagnetic problems , 2007 .

[38]  J. Wait.  On the relation between telluric currents and the Earth's magnetic field , 1954 .

[39]  Dmitry B. Avdeev,et al.  Three-Dimensional Electromagnetic Modelling and Inversion from Theory to Application , 2005 .

[40]  Ki Ha Lee,et al.  A new approach to modeling the electromagnetic response of conductive media , 1989 .

[41]  Stephen D. Gedney,et al.  Convolution PML (CPML): An efficient FDTD implementation of the CFS–PML for arbitrary media , 2000 .

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

[44]  Philip E. Wannamaker,et al.  Advances in three-dimensional magnetotelluric modeling using integral equations , 1991 .

[45]  Michael S. Zhdanov,et al.  Fast numerical modeling of multitransmitter electromagnetic data using multigrid quasi-linear approximation , 2006, IEEE Transactions on Geoscience and Remote Sensing.

[46]  E. Süli,et al.  A convergence analysis of Yee's scheme on nonuniform grids , 1994 .

[47]  E. C. Titchmarsh,et al.  The theory of functions , 1933 .

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

[49]  Rune Mittet,et al.  High-order finite-difference simulations of marine CSEM surveys using a correspondence principle for wave and diffusion fields , 2010 .

[50]  I. S. Gradshteyn,et al.  Table of Integrals, Series, and Products , 1976 .

[51]  Gerald W. Hohmann,et al.  Diffusion of electromagnetic fields into a two-dimensional earth; a finite-difference approach , 1984 .