High-order finite-difference simulations of marine CSEM surveys using a correspondence principle for wave and diffusion fields

The computer time required to solve a typical 3D marine controlled-source electromagnetic surveying (CSEM) simulation can be reduced by more than one order of magnitude by transforming low-frequency Maxwell equations in the quasi-static or diffusive limit to a hyperbolic set of partial differential equations that give a representation of electromagnetic fields in a fictitious wave domain. The dispersion and stability analysis can be made equivalent to that of other types of wave simulation problems such as seismic acoustic and elastic modeling. Second-order to eighth-order spatial derivative operators are implemented for flexibility. Fourth-order and sixth-order methods are the most numerically efficient implementations for this particular scheme. An implementation with high-order operators requires that both electric and magnetic fields are extrapolated simultaneously into the air layer. The stability condition given for high-order staggered-derivative operators here should be equally valid for seismic-wave simulation. The bandwidth of recovered fields in the diffusive domain is independent of the bandwidth of the fields in the fictitious wave domain. The fields in the fictitious wave domain do not represent observable fields. Propagation paths and interaction/reflection amplitudes are not altered by the transform from the fictitious wave domain to the diffusive frequency domain; however, the transform contains an exponential decay factor that damps down late arrivals in the fictitious wave domain. The propagation paths that contribute most to the diffusive domain fields are airwave (shallow water) plus typically postcritical events such as refracted and guided waves. The transform from the diffusive frequency domain to the fictitious wave domain is an ill-posed problem. The transform is nonunique. This gives a large degree of freedom in postulating temporal waveforms for boundary conditions in the fictitious wave domain that reproduce correct diffusive frequency-domain fields.

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

[2]  M. A. Dablain,et al.  The application of high-order differencing to the scalar wave equation , 1986 .

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

[4]  Olav Holberg,et al.  COMPUTATIONAL ASPECTS OF THE CHOICE OF OPERATOR AND SAMPLING INTERVAL FOR NUMERICAL DIFFERENTIATION IN LARGE-SCALE SIMULATION OF WAVE PHENOMENA* , 1987 .

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

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

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

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

[9]  R. Kosloff,et al.  Absorbing boundaries for wave propagation problems , 1986 .

[10]  M. Gershenson.  Simple interpretation of time-domain electromagnetic sounding using similarities between wave and diffusion propagation , 1997 .

[11]  A. Tarantola.  Inversion of seismic reflection data in the acoustic approximation , 1984 .

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

[13]  Moshe Reshef,et al.  A nonreflecting boundary condition for discrete acoustic and elastic wave equations , 1985 .

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

[15]  B. Ursin,et al.  Electromagnetic fields in planarly layered anisotropic media , 2007 .

[16]  I. Brevik,et al.  Investigating the exploration potential for 3D CSEM using a calibration survey over the Troll Field , 2009 .

[17]  Jean-Pierre Berenger,et al.  A perfectly matched layer for the absorption of electromagnetic waves , 1994 .

[18]  Michael Commer,et al.  An accelerated time domain finite difference simulation scheme for three‐dimensional transient electromagnetic modeling using geometric multigrid concepts , 2006 .

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

[20]  Lucy MacGregor,et al.  Sea Bed Logging (SBL), a new method for remote and direct identification of hydrocarbon filled layers in deepwater areas , 2002 .

[21]  Jenö Gazdag,et al.  Modeling of the acoustic wave equation with transform methods , 1981 .

[22]  Svein Ellingsrud,et al.  The Meter Reader—Remote sensing of hydrocarbon layers by seabed logging (SBL): Results from a cruise offshore Angola , 2002 .