A 3D full waveform inversion method based on reconstructed source wavelet with wave equation

LIANG Zhanyuan1 WU Guochen1 WANG Yumei2

(1.School of Geosciences, China University of Petroleum (East China), Qingdao, Shandong, China 266580)
(2.Research Institute of Geophysics, Shengli Oilfield Branch Co., SINOPEC, Dongying, Shandong, China 257022)

【Abstract】We propose in this paper a 3D full waveform inversion method based on reconstructed source wavelet with wave equation. We resolve the wave equation with the direct waves at near offsets of prestack shot gathers as boundary conditions and take the time series at the source point as source wavelet. This method uses the principle of source wavefield reconstruction to simulate the inverse process of the wave propagation along the surface, and extracts direct waves at near offsets, which avoids the interference from the waves such as refraction wave, reflection wave and diving wave. So the method has less dependence on the initial model. For tests, we select the 3D SEG/EAGE overthrust model and field data to analyze the phase difference of 3D seismic wavelet based on the source wavefield reconstruction theory of wave equation. The test results show that there is only phase difference between the reconstructed wavelet and the input wavelet. After the phase adjustment, the reconstructed wavelet is used for the 3D full waveform inversion; the inversion velocity is close to the real one. The results are more reliable at channels. The waveform and phase between the observed records and simulated records have also good correspondence, which validates the feasibility and applicability of the proposed method.

【Keywords】 seismic wavelet; full waveform inversion; source wavefield reconstruction; wave equation; SEG/EAGE overthrust model;


【Funds】 National Key Basic Research Program of China (973 Program) (2013CB228604) Joint Funds of the National Natural Science Foundation of China (U1562215)

Download this article


    [1] Lailly P. The seismic inverse problem as a sequence of before stack migrations. Conference on Inverse Scattering: Theory and Application. Society for Industrial and Applied Mathematics, Philadelphia, PA, 1983, 206–220.

    [2] Tarantola A. Inversion of seismic reflection data in the acoustic approximation. Geophysics, 1984, 49 (8): 1259–1266.

    [3] Wang Y W, Dong L G, Huang C et al. A multi-step strategy for mitigating severe nonlinearity in elastic full-waveform inversion. OGP, 2016, 51 (2): 288–294 (in Chinese).

    [4] Li Z Y, Li Z C, Liu Y Y et al. Tomographic waveform inversion with wave-field decomposition. OGP, 2016, 51 (2): 295–300 (in Chinese).

    [5] Virieux J, Operto S. An overview of full-waveform inversion in exploration geophysics. Geophysics, 2009, 74 (6): WCC1–WCC26.

    [6] Fu J Y, Li Z C, Yang G Q et al. Wave equation reflection travel-time inversion in acoustic media. OGP, 2015, 50 (6): 1134–1140 (in Chinese).

    [7] Maurer H, Greenhalgh S, Latzel S. Frequency and spatial sampling strategies for crosshole seismic waveform spectral inversion experiments. Geophysics, 2009, 74 (6): C79–C89.

    [8] Krebs J R, Anderson J E, Hinkley D et al. Fast fullwavefield seismic inversion using encoded sources. Geophysics, 2009, 74 (6): C177–C188.

    [9] Habashy T, Abubakar A, Pan G et al. Source-receiver compression scheme for full-waveform seismic inversion. Geophysics, 2011, 76 (4): R95–R108.

    [10] Ben-Hadj-Ali H, Operto S, Virieux J. An efficient frequency-domain full waveform inversion method using simultaneous encoded sources. Geophysics, 2011, 76 (4): R109–R124.

    [11] Li X, Aravkin A, van Leeuwen T et al. Fast randomized full-waveform inversion with compressive sensing. Geophysics, 2012, 77 (3): A13–A17.

    [12] Gauthier O, Virieux J, Tarantola A. Two-dimensional nonlinear inversion of seismic waveforms: Numerical results. Geophysics, 1986, 51 (7): 1387–1403.

    [13] Feng o, Wang H Z, Tian L X et al. A strategy for source wavefield reconstruction in reverse time migration. SEG Technical Program Expanded Abstracts, 2011, 20: 3164–3168.

    [14] Chen L S, Wang G E. 大地电磁测深法. Beijing: Geological Publishing House, 1984 (in Chinese).

    [15] Yang P J, Yin X Y. Summary of seismic wavelet pick-up. OGP, 2008, 43 (1): 123–128 (in Chinese).

    [16] Lee K H, Kim H J. Source-independent full-waveform inversion of seismic data. Geophysics, 2003, 68 (6): 2010–2015.

    [17] Xu K, Greenhalgh S A, Wang M Y. Comparison of source-independent methods of elastic waveform inversion. Geophysics, 2006, 71 (6): R91–R100.

    [18] Choi Y, Min D J. Source-independent elastic waveform inversion using a logarithmic wavefield. Journal of Applied Geophysics, 2012, 76: 13–22.

    [19] Zhou B, Greenhalgh S A. Crosshole seismic inversion with normalized full-waveform amplitude data. Geophysics, 2003, 68 (4): 1320–1330.

    [20] Choi Y, Alkhalifah T. Source-independent time-domain waveform inversion using convolved wavefields: Application to the encoded multisource waveform inversion. Geophysics, 2011, 76 (5): R125–R134.

    [21] Song Z M, Williamson P R, Pratt R G. Frequency-domain acoustic-wave modeling and inversion of crosshole data: Part Ⅱ—Inversion method, synthetic experiments and real-data results. Geophysics, 1995, 60 (3): 796–809.

    [22] Ao R D, Dong L G, Chi B X. Source-independent envelope-based FWI to build an initial model. Chinese Journal of Geophysics, 2015, 58 (6): 1998–2010 (in Chinese).

    [23] Zhou H, Amundsen L, Zhang G. Fundamental issues in full waveform inversion. SEG Technical Program Expanded Abstracts, 2012, 31.

    [24] Warner M, Ratcliffe A, Nangoo T et al. Anisotropic 3D full-waveform inversion. Geophysics, 2013, 78 (2): R59–R80.

    [25] Ma Y, Hale D. Quasi-Newton full-waveform inversion with a projected Hessian matrix. Geophysics, 2012, 77 (5): R207–R216.

This Article


CN: 13-1095/TE

Vol 52, No. 06, Pages 1200-1207+1119-1120

December 2017


Article Outline


  • 1 Introduction
  • 2 Methods and principles
  • 3 Numerical experiment
  • 4 Field application
  • 5 Conclusions
  • References