Frequency domain elastic waveform inversion using the Gauss-Newton method

Chung, W., Shin, J., Bae, H.S., Yang, D. and Shin, C.S., 2012. Frequency domain elastic waveform inversion using the Gauss-Newton method. Journal of Seismic Exploration, 21: 29-48. In spite of several advantages of using the Gauss-Newton method for waveform inversion, calculating a Hessian-matrix is problematic in the current computing environment. Consequently, many researchers have suggested a Conjugate Gradient Least Square (CGLS) algorithm to circumvent the Hessian matrix problem in the Gauss-Newton method. There has been no attempt, however, at using the CGLS method to perform elastic media waveform inversion. It is necessary to consider the characteristics of elastic media because the acoustic wave equation does not accurately simulate waveforms propagated through elastic media. Our objective is to develop a frequency domain waveform inversion algorithm using the CGLS method for elastic media. Numerical experiments with the Marmousi-2 model and the SEG/EAGE salt model verify enhancement of the inversion results. We confirm that a greatly improved waveform inversion can be carried out over a complicated layer structure with the Marmousi-2 model. We also verify that structures containing a high velocity salt can be reproduced correctly through an inversion experiment using the SEG/EAGE salt model. However, the improved inversion result of the SEG/EAGE salt model is limited to use with the long wavelength initial velocity, necessitating a method for obtaining the long wavelength velocity model.
- Aminzadeh, F., Brac, J. and Kunz, T., 1997. 3-D Modeling Series No.1: 3-D Salt and Overthrust
- Models. SEG, Tulsa, OK and EAGE, Houten.
- Chen, P., Jordan, T.H. and Zhao, L., 2007. Full three-dimensional tomography: a comparison
- between the scattering-integral and adjoint-wavefield methods. Geophys. J. Internat., 170:
- 175-181.
- Epanomeritakis, I., Akcelik, V., Ghattas, O. and Bielak, J., 2008. A Newton-CG method for
- large-scale three dimensional elastic full-waveform seismic inversion. Inverse Probl., 24:
- Golub, G.H. and Van Loan, C.F., 1996. Matrix Computation, 3rd Ed. The Johns Hopkins
- University Press, Baltimore.
- Hu, W., Abubakar, A. and Habashy, T.M., 2009. Simultaneous multifrequency inversion of
- full-waveform seismic data. Geophysics, 74: R1-R14.
- Kolb, P., Collino, F, and Lailly, P., 1986. Pre-stack inversion of a 1-D medium. Proc. IEEE, 74:
- 498-508.
- Lailly, P., 1983. The seismic inverse problem as a sequence of before stack migrations. In: Bednar,
- J.B. (Ed.), Conference on Inverse Scattering: Theory and Application. University of Tulsa,
- OK: 206-220.
- Marfurt, K.J., 1984. Accuracy of finite-difference and finite-element modeling of the scalar and
- elastic wave equations. Geophysics, 49: 533-549.
- Martin, G.S., Wiley, R. and Marfurt, K.J., 2006. Marmousi-2: An elastic upgrade for Marmousi.
- The Leading Edge, 25: 156-166.
- Mora. P., 1987. Nonlinear two-dimensional elastic inversion of multioffset seismic data. Geophysics,
- 52: 1211-1228.
- Pratt, R.G., Shin, C.S. and Hicks, G.J., 1998. Gauss-Newton and full Newton methods in
- frequency-space seismic waveform inversion. Geophys. J. Internat., 133: 341-362.
- Shin, C.S., Pyun, S. and Bednar, J.B., 2007. Comparison of waveform inversion, Part 1:
- Conventional wavefield vs. logarithmic wavefield. Geophys. Prosp., 55: 449-464.
- Shin, C.S. and Min, D.J., 2006. Waveform inversion using a logarithmic wavefield. Geophysics,
- 71: R31-R42.
- Tarantola, A., 1984. Inversion of seismic reflection data in the acoustic approximation. Geophysics,
- 49: 1259-1266.
- Versteeg, R., 1994. The Marmousi experience: Velocity model determination on a synthetic complex
- data set. The Leading Edge, 13: 927-936.