Joint Inversion of Electromagnetic and Seismic Data Based on Structural Constraints Using Variational Born Iteration Method

An efficient 2-D joint full-waveform inversion method for electromagnetic and seismic data in a layered medium background is developed. The joint inversion method based on the integral equation (IE) method is first proposed in this paper. In forward computation, the IE method is employed, which usually has smaller discretized computation domain and less cumulative error compared with the finite-difference method. In addition, fast Fourier transform is used to accelerate the convolution between Green’s functions and induced sources due to the shift invariance property of the layered Green’s functions in the horizontal direction. In the inversion model, the cross-gradient function is incorporated into the cost function of the separate inversion to enforce the structure similarity between electric conductivity and seismic-wave velocity. We use the improved variational Born iteration method and two different iteration strategies to minimize the cost function and reconstruct the contrasts. Several typical models in geophysical applications are used to validate our joint inversion method, and the numerical simulation results show that joint inversion can improve the inversion results when compared with those from the separate inversion.

[1]  A. Roberts,et al.  A framework for 3-D joint inversion of MT, gravity and seismic refraction data , 2011 .

[2]  David L.B. Jupp,et al.  Joint Inversion of Geophysical Data , 2007 .

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

[4]  A. Abubakar,et al.  Joint inversion approaches for geophysical electromagnetic and elastic full-waveform data , 2012 .

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

[6]  Lin-Ping Song,et al.  A fast 2D volume integral‐equation solver for scattering from inhomogeneous objects in layered media , 2005 .

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

[8]  Alan W. Roberts,et al.  Verification of velocity‐resistivity relationships derived from structural joint inversion with borehole data , 2013 .

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

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

[11]  P. Dell’Aversana,et al.  Joint inversion of rock properties from sonic, resistivity and density well‐log measurements , 2011 .

[12]  Max A. Meju,et al.  Structure‐coupled multiphysics imaging in geophysical sciences , 2011 .

[13]  Zhang Yerong,et al.  Variational Born iteration method and its applications to hybrid inversion , 2000 .

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

[15]  A. Abubakar,et al.  Joint electromagnetic and seismic inversion using structural constraints , 2009 .

[16]  Jean Virieux,et al.  An overview of full-waveform inversion in exploration geophysics , 2009 .

[17]  Rongjiang Wang,et al.  Stability of rapid finite‐fault inversion for the 2014 Mw6.1 South Napa earthquake , 2015 .

[18]  G. E. Archie.  The electrical resistivity log as an aid in determining some reservoir characteristics , 1942 .

[19]  K. Michalski,et al.  Multilayered media Green's functions in integral equation formulations , 1997 .

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

[21]  Philip J. Morris,et al.  Formulas of Acoustics , 2004 .

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

[23]  J. Shewchuk.  An Introduction to the Conjugate Gradient Method Without the Agonizing Pain , 1994 .

[24]  L. MacGregor,et al.  Determination of reservoir properties from the integration of CSEM and seismic data , 2006 .

[25]  Aria Abubakar,et al.  Joint MT and CSEM data inversion using a multiplicative cost function approach , 2011 .

[26]  W. Chew.  Waves and Fields in Inhomogeneous Media , 1990 .

[27]  Charles L. Lawson,et al.  Solving least squares problems , 1976, Classics in applied mathematics.

[28]  S. Constable.  Ten years of marine CSEM for hydrocarbon exploration , 2010 .

[29]  Neil Genzlinger.  A. and Q , 2006 .

[30]  W. Marsden.  I and J , 2012 .

[31]  Weng Cho Chew,et al.  Integral Equation Methods for Electromagnetic and Elastic Waves , 2007, Synthesis Lectures on Computational Electromagnetics.

[32]  M. H. Waxman,et al.  Electrical Conductivities in Oil-Bearing Shaly Sands , 1968 .

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

[34]  Hicks,et al.  Gauss–Newton and full Newton methods in frequency–space seismic waveform inversion , 1998 .

[35]  Weng Cho Chew,et al.  An iterative solution of the two‐dimensional electromagnetic inverse scattering problem , 1989, Int. J. Imaging Syst. Technol..

[36]  Aria Abubakar,et al.  Joint petrophysical inversion of electromagnetic and full-waveform seismic data , 2012 .

[37]  Qing Huo Liu,et al.  Three-dimensional reconstruction of objects buried in layered media using Born and distorted Born iterative methods , 2004, IEEE Geosci. Remote. Sens. Lett..

[38]  W. Chew,et al.  Reconstruction of two-dimensional permittivity distribution using the distorted Born iterative method. , 1990, IEEE transactions on medical imaging.