ARTICLE

Time domain full waveform inversion of seismic data for the East Sea, Korea, using a pseudo-Hessian method with source estimation

JIHO HA1 SUNGRYUL SHIN2 CHANGSOO SHIN3 WOOKEEN CHUNG4
Show Less
4 Department of Energy and Resources Engineering, Korea Maritime and Ocean University, Busan, South Korea.,
3 Department of Energy Systems Engineering, Seoul National University, Seoul, South Korea.,
JSE 2015, 24(5), 419–437;
Submitted: 25 July 2014 | Accepted: 4 August 2015 | Published: 1 November 2015
© 2015 by the Authors. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution -Noncommercial 4.0 International License (CC-by the license) ( https://creativecommons.org/licenses/by-nc/4.0/ )
Abstract

Ha, J., Shin, S., Shin, C.S. and Chung, W., 2015. Time domain full waveform inversion of seismic data for the East Sea, Korea, using a pseudo-Hessian method with source estimation. Journal of Seismic Exploration, 24: 419-437. Waveform inversion is used to estimate subsurface velocity information using seismic datasets. To overcome the computational burden, the use of back-propagation algorithm and the pseudo-Hessian matrix are proposed. Many researchers using these algorithms have shown successful results with synthetic and field data tests. In particular, the computational efficiency of waveform inversion is improved by using a pseudo-Hessian with regularization using a virtual source vector. However, these theoretical concepts have been mainly applied to waveform inversion in the frequency or Laplace domains. We propose full waveform inversion using an estimated source wavelet with the pseudo-Hessian matrix and back-propagation in the time domain. We derive the virtual source vector for the first order hyperbolic equation based on 2D staggered-grid modeling. The updated gradient direction is obtained from both the virtual source and back-propagation wavefield vectors. To improve the availability for field data sets, we also perform the source estimation using deconvolution of the observed data based on the least-squares method. In a synthetic test with a modified Marmousi2 model, the inverted velocity model obtained by the proposed waveform inversion algorithm using estimated wavelets shows similarity to the true velocity. The estimated source wavelet shows good agreement with the true source wavelet. We also test the proposed waveform inversion with field data from the East Sea, Korea. The calculated traces with the estimated source wavelet and the inverted velocity model show direct and reflection events similar to those in the real seismic traces. This confirms that the proposed algorithm can be applied to field data.

Keywords
waveform inversion
pseudo-Hessian matrix
back-propagation algorithm
time domain
source estimation
seismic data
References
  1. Bae, H., Shin, C., Cha, Y., Choi, Y. and Min, D.-J., 2010. 2D acoustic-elastic coupled waveforminversion in the Laplace domain. Geophys. Prosp., 58: 997-1010.
  2. Bae, H.-S., Pyun, S., Shin, C., Marfurt, J.K. and Chung, W., 2012. Laplace-domain waveforminversion versus refraction-traveltime tomography. Geophys. J. Internat., 190: 595-606.
  3. Brossier, R., Operto, S. and Virieux, J., 2009. Seismic imaging of complex onshore structures by2D elastic frequency-domain full-waveform inversion. Geophysics, 74: WCC105-WCC118.
  4. Choi, Y. and Alkhalifah, T., 2011. Source-independent time-domain waveform inversion usingconvolved wavefield: Application to the encoded multisource waveform inversion.Geophysics, 76: R125- R1234.
  5. Claerbout, J.F., 1971. Toward a unified theory of reflector mapping. Geophysics, 36: 467-481.
  6. Fichtner, A., 2011. Full seismic waveform modeling and inversion. Springer Verlag, Heidelberg,Berlin.
  7. Gauthier, O., Virieux, J. and Tarantola, A., 1986. Two-dimensional nonlinear inversion of seismicwaveforms: numerical results. Geophysics, 51: 1387-1403.
  8. Geller, R.J. and Hara, T., 1993. Two efficient algorithms for iterative linearized inversion ofseismic waveform data. Geophys. J. Internat.,115: 699-710.
  9. Graves, R.W., 1996. Simulating seismic wave propagation in 3D elastic media, using staggered-gridfinite differences. Bull. Seismol. Soc. Am., 86: 1091-1106.
  10. Ha, T., Chung, W. and Shin, C., 2009. Waveform inversion using a back-propagation algorithmand a Huber function norm. Geophysics, 74: R15-R24.
  11. Ha, T., Pyun, S. and Shin, C., 2006. Efficient electric resistivity inversion using adjoint state ofmixed finite-element method for Poisson’s equation. J. Computat. Phys., 214: 171-186.
  12. Hampson, D.P., Russell, B.H. and Bankhead, B., 2005. Simultaneous inversion of pre-stack seismicdata. Expanded Abstr., 75th Ann. Internat. SEG Mtg., Houston: 1633-1637.
  13. Jin, S. and Madariaga, R., 1993. Background velocity inversion by a genetic algorithm. Geophys.Res. Lett., 20: 93-96.
  14. Jin, S. and Madariaga, R., 1994. Non-linear velocity inversion by a two-step Monte Carlo method.Geophysics, 59: 577-590.
  15. Kim, Y., Min, D.-J. and Shin, C., 2011. Frequency-domain reverse time migration with sourceestimation. Geophysics, 76: S41-S49.
  16. Koo, N.-H., Shin, C., Min, D.-J., Park, K.-P. and Lee H.-Y., 2011. Source estimation and directwave reconstruction in Laplace-domain waveform inversion for deep-sea seismic data.Geophys. J. Internat., 187: 861-870.
  17. Koren, Z., Mosegaard, K., Landa, E., Thore, P. and Tarantola, A., 1991. Monte Carlo estimationand resolution analysis of seismic background velocities. J. Geophys. Res. , 96: 20289-20299.
  18. Lailly, P., 1983. The seismic inverse problem as a sequence of before stack migrations: Conferenceon Inverse Scattering: Theory and Application. Expanded Abstr., Soc. Industr. Appl.Mathemat. Conf.: 206-220.
  19. Levander, A.R., 1988. Fourth-order finite-difference P-SV seismograms. Geophysics, 53:1425-1436.
  20. Métivier, L., Bretaudeau, F., Brossier, R., Operto, S., and Virieux, J., 2014. Full waveforminversion and the truncated newton-method: quantitative imaging of complex subsurfacestructures. Geophys. Prosp., 62: 1353-1375.TIME-DOMAIN FULL-WAVEFORM INVERSION 437
  21. Mosegaard, K. and Tarantola, A., 1995. Monte Carlo sampling of solutions to inverse problems.J. Geophys. Res., 100: 12431-12447.
  22. Pratt, P., 1999. Seismic waveform inversion in the frequency domain, part I: Theory andverification in a physical scale model. Geophysics, 64: 888-901.
  23. Pratt, R.G. and Shipp, R.M., 1999. Seismic waveform inversion in the frequency domain, Part 2:
  24. Fault delineation in sediments using cross hole data. Geophysics, 64: 902-914.
  25. Pratt, R.G., Shin, C. and Hicks, G.J., 1998. Gauss-Newton and full Newton methods infrequency-space seismic waveform inversion. Geophys. J. Internat., 133: 341-362.
  26. Sambridge, M. and Mosegaard, K., 2002. Monte Carlo methods in geophysical inverse problems.Rev. Geophys., 40: 1-29.
  27. Sheen, D.-H., Tuncay, K., Baag, C.-E. and Peter, J.0., 2006. Time domain Gauss-Newton seismicwaveform inversion in elastic media. Geophys. J. Internat., 167: 1373-1384.
  28. Shin, C. and Cha, Y., 2008. Waveform inversion in the Laplace domain. Geophys. J. Internat. , 173:922-931.
  29. Shin, C. and Min, D., 2006. Waveform inversion using a logarithm wavefield. Geophysics, 71:R31-R42.
  30. Shin, C., Jang, S. and Min, D.-J., 2001. Improved amplitude preservation for prestack depthmigration by inverse scattering theory. Geophys. Prosp., 49: 592-606.
  31. Shin, C., Pyun, S. and Bednar, J.B., 2007. Comparison of waveform inversion, part 1: conventionalwavefield vs. logarithmic wavefield. Geophys. Prosp., 55: 449-464.
  32. Shin, C., Yoon, K., Marfurt, K.J., Park, K., Yang, D., Lim, H.Y., Chung, S. and Shin, S., 2001.
  33. Efficient calculation of a partial-derivative wavefield using reciprocity for seismic imagingand inversion. Geophysics, 66: 1856-1863.
  34. Tarantola, A., 1984. Inversion of seismic reflection data in the acoustic approximation. Geophysics,49: 1259-1266.
  35. Toxopeus, G., Thorbecke, J., Wapenaar, C.P.A., Petersen, S., Slob, E. and Fokkema, J.T., 2008.
  36. Simulating migrated and inverted seismic data by filtering a geologic model. Geophysics,73: T1-T10.
  37. Virieux, J., 1986. P-SV wave propagation in heterogeneous media: velocity-stress finite-differencemethod. Geophysics, 51: 889-901.
  38. Virieux, J. and Operto, S., 2009. An overview of full-waveform inversion in exploration geophysics.Geophysics, 74: WCC1-WCC26.
  39. Zong, Z., Yin., X. and Wu, G., 2013. Multi-parameter nonlinear inversion with exact reflectioncoefficient equation. J. Appl. Geophys., 98: 21-32.
Share
Back to top
Journal of Seismic Exploration, Electronic ISSN: 0963-0651 Print ISSN: 0963-0651, Published by AccScience Publishing