Exploration Geophysics Exploration Geophysics Society
Journal of the Australian Society of Exploration Geophysicists
RESEARCH ARTICLE

An efficient waveform inversion using the common mid-point gather in the wavenumber-space-time domain

Yunhui Park 1 Sukjoon Pyun 1 2
+ Author Affiliations
- Author Affiliations

1 Department of Energy Resources Engineering, Inha University, 100 Inha-ro, Nam-gu, Incheon 22212, South Korea.

2 Corresponding author. Email: pyunsj@inha.ac.kr

Exploration Geophysics 48(3) 219-225 https://doi.org/10.1071/EG16019
Submitted: 18 February 2016  Accepted: 18 February 2016   Published: 16 March 2016
Originally submitted to KSEG 14 October 2015, accepted 15 February 2016  

Abstract

As full waveform inversion (FWI) requires large computation time, a variety of techniques have been suggested to reduce the computational burden. In this study, we use wavenumber-space-time domain modelling, which directly generates common mid-point (CMP) gathers, to implement the FWI algorithm. The modelling technique, which is suitable for laterally invariant velocity models, synthesises CMP gathers efficiently by using limited wavenumber components, and therefore allows reduced computation time for FWI. To consider the accuracy as well as the efficiency of FWI, the Gauss-Newton method using the approximate Hessian matrix is chosen in this study. Rather than using all of the wavenumber components, our FWI algorithm can be accelerated by using only a few components. The wavenumber components can be chosen through an analysis of the residual wavefields. To validate the usefulness of our method, we first use a one-dimensional (1D) velocity model. From the 1D model example, we note that our FWI algorithm can be successful if given a reliable initial velocity model and sufficient data with a long offset distance. Even though our algorithm is valid for only horizontally layered velocity models, we apply our algorithm to a two-dimensional (2D) velocity model with lateral velocity variations. Through the 2D velocity example, we confirm that our FWI can be used to estimate subsurface structures with dipping interfaces if the dips are moderate and the structures can thus be considered to be locally flat.

Key words: CMP gather, Fourier transform, full waveform inversion, Gauss-Newton method, wavenumber-space-time domain.


References

Ben-Hadj-Ali, H., Operto, S., and Virieux, J., 2008, 3D frequency-domain full-waveform tomography based on a domain decomposition forward problem: 78th Annual International Meeting, SEG, Expanded Abstracts, 1945–1949.

Ben-Hadj-Ali, H., Operto, S., and Virieux, J., 2011, An efficient frequency-domain full waveform inversion method using simultaneous encoded sources: Geophysics, 76, R109–R124
An efficient frequency-domain full waveform inversion method using simultaneous encoded sources:CrossRef |

Jeong, W., Pyun, S., Son, W., and Min, D. J., 2013, A numerical study of simultaneous-source full waveform inversion with l1-norm: Geophysical Journal International, 194, 1727–1737
A numerical study of simultaneous-source full waveform inversion with l1-norm:CrossRef |

Krebs, J. R., Anderson, J. E., Hinkley, D., Neelamani, R., Lee, S., Baumstein, A., and Lacasse, M. D., 2009, Fast full-wavefield seismic inversion using encoded sources: Geophysics, 74, WCC177–WCC188
Fast full-wavefield seismic inversion using encoded sources:CrossRef |

Liu, Z., and Zhang, J., 2013, Elastic full waveform inversion for near surface imaging in CMP domain: 83rd Annual International Meeting, SEG, Expanded Abstracts, 1904–1908.

Liu, Z., and Zhang, J., 2014, Pseudo 2D elastic full waveform inversion for near surface imaging: 84th Annual International Meeting, SEG, Expanded Abstracts, 1017–1021.

Liu, L., Ding, R., Liu, H., and Liu, H., 2015, 3D hybrid-domain full waveform inversion on GPU: Computers & Geosciences, 83, 27–36
3D hybrid-domain full waveform inversion on GPU:CrossRef |

Martin, G. S., Larsen, S., and Marfurt, K., 2002, Marmousi-2: an updated model for the investigation of AVO in structurally complex areas: 72nd Annual International Meeting, SEG, Expanded Abstracts, 1979–1982.

Moghaddam, P. P., Keers, H., Herrmann, F. J., and Mulder, W. A., 2013, A new optimization approach for source-encoding full-waveform inversion: Geophysics, 78, R125–R132
A new optimization approach for source-encoding full-waveform inversion:CrossRef |

Park, Y., and Pyun, S., 2015, Source wavelet estimation using common mid-point gathers: Geosystem Engineering, 18, 199–204
Source wavelet estimation using common mid-point gathers:CrossRef |

Park, Y., Pyun, S., and Ha, W., 2014, An efficient wavenumber-space-time domain finite-difference modeling of acoustic wave equation for synthesizing CMP gathers: Geosystem Engineering, 17, 287–293
An efficient wavenumber-space-time domain finite-difference modeling of acoustic wave equation for synthesizing CMP gathers:CrossRef |

Pratt, R. G., Shin, C., and Hick, G. J., 1998, Gauss–Newton and full Newton methods in frequency-space seismic waveform inversion: Geophysical Journal International, 133, 341–362
Gauss–Newton and full Newton methods in frequency-space seismic waveform inversion:CrossRef |

Pyun, S., Son, W., and Shin, C., 2011, 3D acoustic waveform inversion in the Laplace domain using an iterative solver: Geophysical Prospecting, 59, 386–399
3D acoustic waveform inversion in the Laplace domain using an iterative solver:CrossRef |

Reynolds, A. C., 1978, Boundary conditions for the numerical solutions of wave propagation problems: Geophysics, 43, 1099–1110
Boundary conditions for the numerical solutions of wave propagation problems:CrossRef |

Shin, C., Yoon, K., Marfurt, K. J., Park, K., Yang, D., Lim, H. Y., Chung, S., and Shin, S., 2001, Efficient calculation of a partial-derivative wavefield using reciprocity for seismic imaging and inversion: Geophysics, 66, 1856–1863
Efficient calculation of a partial-derivative wavefield using reciprocity for seismic imaging and inversion:CrossRef |

Shin, C., Pyun, S., and Bednar, J., 2007, Comparison of waveform inversion, part 1: conventional wavefield vs logarithmic wavefield: Geophysical Prospecting, 55, 449–464
Comparison of waveform inversion, part 1: conventional wavefield vs logarithmic wavefield:CrossRef |

Shin, J., Ha, W., Jun, H., Min, D. J., and Shin, C., 2014, 3D Laplace-domain full waveform inversion using a single GPU card: Computers & Geosciences, 67, 1–13
3D Laplace-domain full waveform inversion using a single GPU card:CrossRef |

Sirgue, L., and Pratt, R. G., 2004, Efficient waveform inversion and imaging: a strategy for selecting temporal frequencies: Geophysics, 69, 231–248
Efficient waveform inversion and imaging: a strategy for selecting temporal frequencies:CrossRef |

Son, W., Pyun, S., Shin, C., and Kim, H. J., 2014, An algorithm adapting encoded simultaneous-source full-waveform inversion to marine-streamer acquisition data: Geophysics, 79, R183–R193
An algorithm adapting encoded simultaneous-source full-waveform inversion to marine-streamer acquisition data:CrossRef |

Vigh, D., and Starr, E. W., 2008, 3D prestack plane-wave, full-waveform inversion: Geophysics, 73, VE135–VE144
3D prestack plane-wave, full-waveform inversion:CrossRef |

Yang, P., Gao, J., and Wang, B., 2015, A graphics processing unit implementation of time-domain full-waveform inversion: Geophysics, 80, F31–F39
A graphics processing unit implementation of time-domain full-waveform inversion:CrossRef |

Zhu, L., and Rivera, L. A., 2002, A note on the dynamic and static displacements from a point source in multilayered media: Geophysical Journal International, 148, 619–627
A note on the dynamic and static displacements from a point source in multilayered media:CrossRef |



Export Citation