Joint electromagnetic and seismic inversion using structural constraints

We have developed a frequency-domain joint electromagnetic (EM) and seismic inversion algorithm for reservoir evaluation and exploration applications. EM and seismic data are jointly inverted using a cross-gradient constraint that enforces structural similarity between the conductivity image and the compressional wave (P-wave) velocity image. The inversion algorithm is based on a Gauss-Newton optimization approach. Because of the ill-posed nature of the inverse problem, regularization is used to constrain the solution. The multiplicative regularization technique selects the regularization parameters automatically, improving the robustness of the algorithm. A multifrequency data-weighting scheme prevents the high-frequency data from dominating the inversion process. When the joint-inversion algorithm is applied in integrating marine controlled-source electromagnetic data with surface seismic data for subsea reservoir exploration applications and in integrating crosswell EM and sonic data for reservoir monitoring and evaluation applications, results improve significantly over those obtained from separate EM or seismic inversions.

[1]  Philip E. Gill,et al.  Practical optimization , 1981 .

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

[3]  A. Tarantola A strategy for nonlinear elastic inversion of seismic reflection data , 1986 .

[4]  R. Pratt,et al.  INVERSE THEORY APPLIED TO MULTI‐SOURCE CROSS‐HOLE TOMOGRAPHY.: PART 1: ACOUSTIC WAVE‐EQUATION METHOD1 , 1990 .

[5]  R. N. Edwards,et al.  12. Electrical Exploration Methods for the Seafloor , 1991 .

[6]  Douglas W. Oldenburg,et al.  Inversion of time-domain electromagnetic data for a horizontally layered Earth , 1993 .

[7]  P. M. van den Berg,et al.  Seismic applications of acoustic reciprocity , 1993 .

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

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

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

[11]  Z. M. Song,et al.  Frequency-domain acoustic-wave modeling and inversion of crosshole data; Part II, Inversion method, synthetic experiments and real-data results , 1995 .

[12]  C. S. Cox,et al.  Marine controlled-source electromagnetic sounding. 2. The PEGASUS experiment , 1996 .

[13]  Gene H. Golub,et al.  Matrix computations (3rd ed.) , 1996 .

[14]  Fadil Santosa,et al.  Recovery of Blocky Images from Noisy and Blurred Data , 1996, SIAM J. Appl. Math..

[15]  Curtis R. Vogel,et al.  Iterative Methods for Total Variation Denoising , 1996, SIAM J. Sci. Comput..

[16]  Steven Constable,et al.  Marine controlled‐source electromagnetic sounding: 1. Modeling and experimental design , 1996 .

[17]  D. Oldenburg,et al.  NON-LINEAR INVERSION USING GENERAL MEASURES OF DATA MISFIT AND MODEL STRUCTURE , 1998 .

[18]  P. M. Berg,et al.  Extended contrast source inversion , 1999 .

[19]  Vladimir Druskin,et al.  Optimal finite difference grids and rational approximations of the square root I. Elliptic problems , 2000 .

[20]  Oszkar Biro,et al.  Finite-element analysis of controlled-source electromagnetic induction using Coulomb-gauged potentials , 2001 .

[21]  Steven Constable,et al.  Electrical resistivity structure of the Valu Fa Ridge, Lau Basin, from marine controlled‐source electromagnetic sounding , 2001 .

[22]  P. M. Berg,et al.  Contrast Source Inversion Method: State of Art , 2001 .

[23]  A. Abubakar,et al.  Multiplicative regularization for contrast profile inversion , 2003 .

[24]  M. Meju,et al.  Characterization of heterogeneous near‐surface materials by joint 2D inversion of dc resistivity and seismic data , 2003 .

[25]  Luis A. Gallardo,et al.  Evidence for correlation of electrical resistivity and seismic velocity in heterogeneous near‐surface materials , 2003 .

[26]  M. Meju,et al.  Joint two-dimensional DC resistivity and seismic travel time inversion with cross-gradients constraints , 2004 .

[27]  Aria Abubakar,et al.  A multiplicative regularization approach for deblurring problems , 2004, IEEE Transactions on Image Processing.

[28]  A. Abubakar,et al.  A General Framework for Constraint Minimization for the Inversion of Electromagnetic Measurements , 2004 .

[29]  R. Gerhard Pratt,et al.  Efficient waveform inversion and imaging: A strategy for selecting temporal frequencies , 2004 .

[30]  Laura M. Wallace,et al.  Subduction zone coupling and tectonic block rotations in the North Island, New Zealand , 2004 .

[31]  Y. Rubin,et al.  Direct reservoir parameter estimation using joint inversion of marine seismic AVA and CSEM data , 2005 .

[32]  N. Linde,et al.  Local earthquake (LE) tomography with joint inversion for P‐ and S‐wave velocities using structural constraints , 2006 .

[33]  Changsoo Shin,et al.  Waveform inversion using a logarithmic wavefield , 2006 .

[34]  A. Binley,et al.  Improved hydrogeophysical characterization using joint inversion of cross‐hole electrical resistance and ground‐penetrating radar traveltime data , 2006 .

[35]  David R. Davies,et al.  Petroleum industry and renewables , 2007 .

[36]  S. Hubbard,et al.  Joint inversion of crosshole radar and seismic traveltimes , 2008 .

[37]  L. Gallardo Multiple cross‐gradient joint inversion for geospectral imaging , 2007 .

[38]  Daniele Colombo,et al.  Geophysical modeling via simultaneous joint inversion of seismic, gravity, and electromagnetic data: Application to prestack depth imaging , 2007 .

[39]  Max A. Meju,et al.  Joint two‐dimensional cross‐gradient imaging of magnetotelluric and seismic traveltime data for structural and lithological classification , 2007 .

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

[41]  A. Abubakar,et al.  Simultaneous multifrequency inversion of full-waveform seismic data , 2009 .