Multiscale Full Waveform Inversion Using GPU

We develop an efficient multi-scale time domain full wavefrom inversion (FWI) method using Graphics Processing Unit (GPU). Benefited from the acceleration of GPU, the intensive computational cost of time domain FWI becomes affordable on our desktop. As the geological targets in subsurface area are usually multi-scale, we propose a multi-scale time domain FWI method using GPU to invert the model more efficiently and effectively. The strategy for multi-scale FWI is very similar to frequency domain FWI. Therefore we make a comparison between the results of our multi-scale method and the frequency domain FWI results. We use 2D Overthrust model and 2D Marmousi model to demonstrate the validity of our method.

[1]  Romain Brossier,et al.  A multi-scale strategy for handling broadband seismic data , 2010 .

[2]  Gerard T. Schuster,et al.  Applications of multiscale waveform inversion to marine data using a flooding technique and dynamic early-arrival windows , 2010 .

[3]  G. A Theory for Multiresolution Signal Decomposition : The Wavelet Representation , 2004 .

[4]  A. Abubakar,et al.  Preconditioned non‐linear conjugate gradient method for frequency domain full‐waveform seismic inversion , 2011 .

[5]  R. Clapp Reverse time migration with random boundaries , 2009 .

[6]  O. Savasta,et al.  Wave-equation Traveltime Inversion Using 2D Cross-correlation Function , 1997 .

[7]  Stéphane Mallat,et al.  A Theory for Multiresolution Signal Decomposition: The Wavelet Representation , 1989, IEEE Trans. Pattern Anal. Mach. Intell..

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

[9]  Saulo P. Oliveira,et al.  Linearized seismic waveform inversion using a multiple re‐weighted least‐squares method with QR preconditioning , 2007 .

[10]  Jinghuai Gao,et al.  Fast Full Waveform Inversion of Multi-shot Seismic Data , 2010 .

[11]  Bjørn Ursin,et al.  Preconditioning of full‐waveform inversion in viscoacoustic media , 1999 .

[12]  Changsoo Shin,et al.  2D acoustic‐elastic coupled waveform inversion in the Laplace domain , 2010 .

[13]  Gerard T. Schuster,et al.  Wavepath eikonal traveltime inversion: Theory , 1993 .

[14]  Tariq Alkhalifah,et al.  Source-independent time-domain waveform inversion using convolved wavefields: Application to the encoded multisource waveform inversion , 2011 .

[15]  Xiang Li,et al.  Efficient least‐squares imaging with sparsity promotion and compressive sensing , 2012 .

[16]  R. Pratt Seismic waveform inversion in the frequency domain; Part 1, Theory and verification in a physical scale model , 1999 .

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

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

[19]  P. Mora Nonlinear two-dimensional elastic inversion of multioffset seismic data , 1987 .

[20]  A. Guitton,et al.  Attenuating crosstalk noise with simultaneous source full waveform inversion ★ , 2012 .

[21]  Changsoo Shin,et al.  A comparison between the behavior of objective functions for waveform inversion in the frequency and Laplace domains , 2008 .

[22]  Denes Vigh,et al.  3D prestack plane-wave, full-waveform inversion , 2008 .

[23]  C. Bunks,et al.  Multiscale seismic waveform inversion , 1995 .

[24]  A. Tarantola,et al.  Two‐dimensional nonlinear inversion of seismic waveforms: Numerical results , 1986 .

[25]  B. Engquist,et al.  Absorbing boundary conditions for acoustic and elastic wave equations , 1977, Bulletin of the Seismological Society of America.

[26]  S. Operto,et al.  Which data residual norm for robust elastic frequency-domain full waveform inversion? , 2010 .

[27]  S. Operto,et al.  3D finite-difference frequency-domain modeling of visco-acoustic wave propagation using a massively parallel direct solver: A feasibility study , 2007 .

[28]  Jean Virieux,et al.  Velocity model building by 3D frequency-domain, full-waveform inversion of wide-aperture seismic data , 2008 .

[29]  J. Sochacki Absorbing boundary conditions for the elastic wave equations , 1988 .

[30]  Jean Virieux,et al.  Quantitative imaging of complex structures from dense wide‐aperture seismic data by multiscale traveltime and waveform inversions: a case study , 2004 .

[31]  D. Komatitsch,et al.  An unsplit convolutional perfectly matched layer improved at grazing incidence for the seismic wave equation , 2007 .

[32]  D. Vigh,et al.  3D Pre-Stack Plane Wave Full Waveform Inversion , 2008 .

[33]  Gerard T. Schuster,et al.  An efficient multiscale method for time-domain waveform tomography , 2009 .

[34]  Jean Virieux,et al.  An efficient frequency-domain full waveform inversion method using simultaneous encoded sources , 2011 .

[35]  Yaxun Tang Target-oriented wave-equation least-squares migration/inversion with phase-encoded Hessian , 2009 .

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

[37]  William W. Symes,et al.  Reverse time migration with optimal checkpointing , 2007 .

[38]  R. Shipp,et al.  Seismic waveform inversion in the frequency domain, Part 2: Fault delineation in sediments using crosshole data , 1999 .