A graphics processing unit implementation of time-domain full-waveform inversion

ABSTRACTThe graphics processing unit (GPU) has become a popular device for seismic imaging and inversion due to its superior speed-up performance. We implemented GPU-based full-waveform inversion using the wavefield reconstruction strategy. Because computation on the GPU was much faster than CPU-GPU data communication, in our implementation, the boundaries of the forward modeling were saved on the device to avert the issue of data transfer between the host and device. We adopted the Clayton-Enquist absorbing boundary to maintain the efficiency of the GPU computation. A hybrid nonlinear conjugate gradient algorithm combined with the parallel reduction scheme was used to do computation in GPU blocks. The numerical results confirmed the validity of our implementation.

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

[2]  Markus Hadwiger,et al.  Multisource Reverse-time Migration and Full-waveform Inversion on a GPGPU , 2010 .

[3]  A. Tarantola,et al.  Nonlinear Inversion of Seismic Reflection Data , 1984 .

[4]  Yaxun Tang,et al.  Joint Preconditioned Least-squares Inversion of Simultaneous Source Time-lapse Seismic Data Sets , 2009 .

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

[6]  Heiner Igel,et al.  Accelerated 3-D full-waveform inversion using simultaneously encoded sources in the time domain: application to Valhall ocean-bottom cable data , 2013 .

[7]  Jungkyun Shin,et al.  3D Laplace-domain full waveform inversion using a single GPU card , 2014, Comput. Geosci..

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

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

[10]  W. Hager,et al.  A SURVEY OF NONLINEAR CONJUGATE GRADIENT METHODS , 2005 .

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

[12]  William W. Symes,et al.  Migration velocity analysis and waveform inversion , 2008 .

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

[14]  Paulius Micikevicius,et al.  3D finite difference computation on GPUs using CUDA , 2009, GPGPU-2.

[15]  Paul Sava,et al.  Madagascar: open-source software project for multidimensional data analysis and reproducible computational experiments , 2013 .

[16]  David Yingst,et al.  Viscoacoustic waveform inversion of velocity structures in the time domain , 2014 .

[17]  F. Herrmann,et al.  A new optimization approach for source-encoding full-waveform inversion , 2013 .

[18]  Hongwei Liu,et al.  Wavefield reconstruction methods for reverse time migration , 2013 .

[19]  C. Shin,et al.  Waveform inversion in the Laplace domain , 2008 .

[20]  Jinghuai Gao,et al.  CUDA-based Acceleration of full waveform inversion on GPU , 2011 .

[21]  William W. Symes,et al.  Computational Strategies For Reverse-time Migration , 2008 .

[22]  Antoine Guitton,et al.  Constrained full-waveform inversion by model reparameterization1 , 2012 .

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

[24]  A. Pica,et al.  Nonliner inversion of seismic reflection data in a laterally invariant medium , 1990 .

[25]  신창수.  Apparatus, methods and recording medium for imaging a subsurface using the waveform inversion in the laplace-fourier domain , 2008 .

[26]  Baoli Wang,et al.  RTM using effective boundary saving: A staggered grid GPU implementation , 2014, Comput. Geosci..

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

[28]  A. Majda,et al.  Absorbing boundary conditions for the numerical simulation of waves , 1977 .