2.5D controlled-source EM modeling with general 3D source geometries

Most 2.5D controlled-source electromagnetic (CSEM) modeling algorithms presented to date explicitly consider only sources that are point dipoles oriented parallel or perpendicular to the direction of constant conductivity. This makes simulations of complex source geometries expensive, requiring separate evaluations of many point dipole fields, and thus limits the practical applicability of such schemes for simulating and interpreting field data. We present a novel 2.5D CSEM modeling scheme that overcomes this limitation and permits efficient simulations of sources with general shape and orientation by evaluating fields for the entire source at once. We accommodate general sources by using a secondary field approach, in which primary fields are computed for the general source and a 1D background conductivity model. To carry out the required Fourier transforms between space and wavenumber domain using the same fast cosine and sine transform filters as in conventional algorithms, we split the primary and secondary fields into their symmetric and antisymmetric parts. For complex 3D source geometries, this approach is significantly more efficient than previous 2.5D algorithms. Our finite-difference algorithm also includes novel approaches for divergence correction at low frequencies and EM field interpolation across conductivity discontinuities. We describe the modeling scheme and demonstrate its accuracy and efficiency by comparisons of 2.5D-simulated data with 1D and 3D results.

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

[2]  Kerry Key,et al.  2D marine controlled-source electromagnetic modeling: Part 1 — An adaptive finite-element algorithm , 2007 .

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

[4]  W. A. Mulder,et al.  A multigrid solver for 3D electromagnetic diffusion , 2006 .

[5]  Y. Mitsuhata 2-D electromagnetic modeling by finite‐element method with a dipole source and topography , 2000 .

[6]  Christoph Schwarzbach,et al.  Stability of finite element solutions to Maxwell's equations in frequency domain , 2009 .

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

[8]  Michael Leppin Electromagnetic modeling of 3-D sources over 2-D inhomogeneities in the time domain , 1992 .

[9]  Stefan Lüth,et al.  Monitoring at the CO2 SINK site: A concept integrating geophysics, geochemistry and microbiology , 2009 .

[10]  Rita Streich,et al.  3D finite-difference frequency-domain modeling of controlled-source electromagnetic data: Direct solution and optimization for high accuracy , 2009 .

[11]  Rita Streich,et al.  Electromagnetic fields generated by finite‐length wire sources: comparison with point dipole solutions , 2011 .

[12]  R. Greenfield,et al.  Numerical solutions of the response of a two-dimensional earth to an oscillating magnetic dipole source , 1976 .

[13]  D. D. Snyder A method for modeling the resistivity and IP response of two-dimensional bodies , 1976 .

[14]  Patrick Amestoy,et al.  Hybrid scheduling for the parallel solution of linear systems , 2006, Parallel Comput..

[15]  Alan D. Chave,et al.  Controlled electromagnetic sources for measuring electrical conductivity beneath the oceans: 1. Forward problem and model study , 1982 .

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

[17]  David L. Alumbaugh,et al.  The practical application of 2D inversion to marine controlled-source electromagnetic data , 2010 .

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

[19]  R. N. Edwards,et al.  Transient marine electromagnetics: the 2.5-D forward problem , 1993 .

[20]  G. W. Hohmann,et al.  4. Electromagnetic Theory for Geophysical Applications , 1987 .

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

[22]  Alan D. Chave,et al.  Electromagnetic induction by a finite electric dipole source over a 2-D earth , 1993 .

[23]  Uri M. Ascher,et al.  Fast Finite Volume Simulation of 3D Electromagnetic Problems with Highly Discontinuous Coefficients , 2000, SIAM J. Sci. Comput..

[24]  Tage Røsten,et al.  A 2.5D finite-element-modeling difference method for marine CSEM modeling in stratified anisotropic media , 2008 .

[25]  Ki Ha Lee,et al.  A numerical solution for the electromagnetic scattering by a two-dimensional inhomogeneity , 1985 .

[26]  Walter L. Anderson Fourier cosine and sine transforms using lagged convolutions in double-precision (subprograms DLAGF0/DLAGF1) , 1983 .

[27]  Walter L. Anderson,et al.  Algorithm 588: Fast Hankel Transforms Using Related and Lagged Convolutions , 1982, TOMS.

[28]  Aria Abubakar,et al.  2.5D forward and inverse modeling for interpreting low-frequency electromagnetic measurements , 2008 .