An efficient waveform inversion using the common mid-point gather in the wavenumber-space-time domain
Yunhui Park 1 Sukjoon Pyun 1 21 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
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 | GoogleScholarGoogle Scholar |
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 | GoogleScholarGoogle Scholar |
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 | GoogleScholarGoogle Scholar |
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 | GoogleScholarGoogle Scholar |
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 | GoogleScholarGoogle Scholar |
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 | GoogleScholarGoogle Scholar |
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 | GoogleScholarGoogle Scholar |
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 | GoogleScholarGoogle Scholar |
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 | GoogleScholarGoogle Scholar |
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 | GoogleScholarGoogle Scholar |
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 | GoogleScholarGoogle Scholar |
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 | GoogleScholarGoogle Scholar |
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 | GoogleScholarGoogle Scholar |
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 | GoogleScholarGoogle Scholar |
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 | GoogleScholarGoogle Scholar |
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 | GoogleScholarGoogle Scholar |
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 | GoogleScholarGoogle Scholar |
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 | GoogleScholarGoogle Scholar |