ARTICLE

Frequency-domain weighted-averaging finite-difference numerical simulation of qP wave propagation in TTI media

GUOCHEN WU KAI LIANG XINGYAO YIN
Show Less
College of Geo-Resources and Information, China University of Petroleum (East China), Beier Road #271, Dongying, Shan Dong 257061, P.R. China.,
JSE 2010, 19(3), 207–229;
Submitted: 6 August 2009 | Accepted: 10 March 2010 | Published: 1 July 2010
© 2010 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

Wu, G., Liang, K. and Yin, X., 2010. Frequency-domain weighted-averaging finite-difference numerical simulation of qP wave propagation in TTI media. Journal of Seismic Exploration, 19: 207- 229. The finite-difference method is widely used in numerical simulation of the propagation of seismic waves, but has the limitation that numerical dispersion reduces the accuracy and resolution of seismic wavefield simulation. In order to decrease the numerical dispersion of conventional finite-difference operators, this paper presents a frequency-domain weighted-averaging finite-difference operators defined on a 25-point stencil for numerical simulation of qP waves propagating in transversely isotropic media with a tilted symmetry axis (TTI media). We first approximate the differential operators using finite-difference analogues defined on 25-point stencils and then calculate the weighted average of the difference operators with weighting coefficients. The weighting coefficients are determined by the Gauss-Newton method of optimization theory. Using the weighted-averaging finite-difference analogues and combined boundary conditions, we successfully simulate qP wave propagation in homogeneous TTI media, layered TTI media and VTI Salt model. The seismic wavefields in the time and frequency domains are obtained and used to generate single shot records. The result of numerical simulation indicates that 25-point weighted-averaging finite-difference analogues can improve the accuracy of the numerical simulation of wavefields and efficiently suppress the numerical dispersion of conventional difference operators. This method may be employed in the foundation of qP migration and inversion in TTI media.

Keywords
TTI media
weighted-averaging
finite-difference
weighting coefficients
numerical dispersion.
References
  1. Alkhalifah, T., 1998. Acoustic approximation for processing in transversely isotropic media.Geophysics, 63: 623-631.
  2. Alkhalifah, T., 2000. An acoustic wave equation for anisotropic media. Geophysics, 65: 1239-1250.
  3. Alkhalifah, T., 2003. An acoustic wave equation for orthorhombic media. Geophysics, 68:1169-1172.
  4. Han, Q.Y. and Wu, R.S., 2004. Depth imaging in strongly heterogeneous VTI media. Expanded
  5. Abstr., 74th Ann. Internat. SEG Mtg., Denver: P1061-1064.
  6. Jo, C.H., Shin, C. and Suh, J.H., 1996. An optimal 9-point, finite-difference, frequency-space, 2-Dscalar wave extrapolator. Geophysics, 61: 529-537.
  7. Komatitsch, D. and Tromp, J., 2003. A perfectly matched layer absorbing boundary condition forthe second-order seismic wave equation. Geophys. J. Internat., 154: 146-153.
  8. Liang, K., Wu, G.C. and Yin, X.Y., 2008. Born approximation migration method of qP wave inVTI media. J. Seismic Explor., 17: 169-180.
  9. Lines, L.R. and Treitel, S., 1984. A review of least-squares inversion and its application togeophysical problems. Geophys. Prosp., 32: 159-186.
  10. Min, D.J., Shin, C., Kwon, B.D. and Chung, S., 2000. Improved frequency-domain elastic wavemodeling using weighted- averaging difference operators. Geophysics, 65: 884-895.
  11. Min, D.J., Yoo, H.S., Shin, C., Hyun, H.J. and Suh, J.H., 2002. Weighted-averagingfinite-element method for scalar wave equation in the frequency domain. J. Seismic Explor.,11: 197-222.
  12. Min, D.J., Shin, C., Pratt, R.G. and Yoo, H.S., 2003. Weighted-averaging finite-element methodfor 2D elastic wave equations in the frequency domain. Bull. Seismol. Soc. Am., 93:904-921.
  13. Pratt, R.G., 1990. Frequency-domain elastic wave modeling by finite differences: A tool forcrosshole seismic imaging. Geophysics, 55: 626-632.
  14. Shin, C. and Sohn, H., 1998. A frequency-space 2-D scalar wave extrapolator using extend 25-pointfinite-difference operator. Geophysics, 63: 289-296.
  15. Stekl, I. and Pratt, R.G., 1998. Accurate viscoelastic modeling by frequency-domain finitedifferences using rotated operators. Geophysics, 63: 1779-1794.226 WU, LIANG & YIN
  16. Thomsen, L,. 1986. Weak elastic anisotropy. Geophysics, 51: 1954-1966.
  17. Wu, G.C. and Liang, K., 2005. Quasi P-wave forward modeling in frequency-space domain in VTImedia. Oil Geophys. Prosp., 40: 535-545.
  18. Wu, G.C. and Liang, K., 2005. Combined boundary conditions of quasi-P wave withinfrequency-space domain in VTI media. Geophys. Prosp. for Petrol., 44: 301-307.
  19. Zeng, Y.Q., He, J.Q. and Liu, Q.H., 2001. The application of the perfectly matched layer innumerical modeling of wave propagation in poroelastic media. Geophysics, 66: 1258-1266.
  20. Zhang, L.B., Rector III, J.W. and Hoversten, G.M., 2003. An acoustic wave equation for modelingin tilted TI media. Expanded Abstr., 73rd Ann. Internat. SEG Mtg., Dallas: P153-156.
Share
Back to top
Journal of Seismic Exploration, Electronic ISSN: 0963-0651 Print ISSN: 0963-0651, Published by AccScience Publishing