AccScience Publishing / JSE / Article
Submit to JSE
Cite this article
75
Views
Journal Browser
Volume | Year
Issue
Search
Forthcoming Issue
Current Issue
View All
News and Announcements
View All
ARTICLE

A generalized 17-point scheme based on the directional derivative method for highly accurate finite-difference simulation

WEI LIU1,2 ,  YANMIN HE1,2 ,  SHU LI1,3 ,  HAO WU1,2 ,  LIFENG YANG1,2 ,  ZHENMING PENG1,2
Show Less
1 School of Information and Communication Engineering, University of Electronic Science and Technology of China, Chengdu 610054, P.R. China.,
2 Information Geoscience Research Center, University of Electronic Science and Technology of China, Chengdu 610054, P.R. China.,
3 School of Information Science and Engineering, Jishou University, Jishou 416000, P.R. China.,
JSE 2019, 28(1), 4–34;
Submitted: 30 March 2018 | Accepted: 8 November 2018 | Published: 1 February 2019
© 2019 by the Authors. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution 4.0 International License ( https://creativecommons.org/licenses/by/4.0/ )
Download PDF
Cite
XML
Abstract

Liu, W., He, Y.M., Li, S, Wu, H., Yang, L.F. and Peng, Z.M., 2019. A generalized 17- point scheme based on the directional derivative method for highly accurate finite- difference simulation of the frequency-domain 2D scalar wave equation. Journal of Seismic Exploration, 28: 41-71. Forward modeling of the frequency-domain wave equation represents an essential foundation for full waveform inversion in the frequency domain, the accuracy and efficiency of which rely heavily on the forward modeling method employed. To reduce the numerical dispersion, anisotropy, and number of grids per the shortest wavelength in forward modeling methods, rotating coordinate systems have been successfully applied to establish finite-difference (FD) schemes for the forward modeling of the frequency- domain wave equation. However, rotated optimal FD schemes are incapable of handling rectangular sampling grids, which are ubiquitous in practice. Fortunately, optimal FD schemes based on the average-derivative method (ADM) overcome this restriction on different directional sampling intervals. However, the ADM itself is merely an algebraic approach and therefore does not inherit the geometrical properties of the rotating coordinate system. Based on the principle of a rotating coordinate system, a novel optimal directional derivative method (DDM)-based 4th-order, 17-point FD scheme is developed in this paper for the forward modeling of the frequency-domain, two- dimension scalar wave equation to approximate the spatial derivatives. The conventional 0963-065 1/19/$5.00 © 2019 Geophysical Press Ltd. 42 4th-order, 9-point scheme and rotated optimal 17-point FD scheme can be derived as special cases of the proposed scheme. Compared with the rotated optimal 17-point FD scheme, the proposed scheme is capable of addressing arbitrary rectangular sampling grids, including equal and unequal directional sampling intervals; moreover, the optimized weighted coefficients can reduce the number of grids per the shortest wavelength from 2.56 to less than 2.4 with maximum phase velocity errors of 1%. Furthermore, the proposed scheme is superior to the ADM-based optimal 17-point FD scheme in suppressing numerical dispersion due to the inherited geometrical properties of the rotating coordinate system. A perfectly matched layer boundary condition is applied to the final FD equation to attenuate boundary reflections. Numerical examples demonstrate the validity and adaptability of our 17-point FD scheme.

Keywords
seismic forward modeling
acoustic wave equation
frequency domain
finite difference
numerical dispersion analysis
directional derivative
References
  1. Alford, R.M., Kelly, K.R. and Boore, D.M., 1974. Accuracy of finite-differencemodeling of the acoustic wave equation. Geophysics, 39: 834-842.
  2. Bérenger, J.P., 1994. A perfectly matched layer for the absorption of electromagneticwaves. J. Comput. Phys., 114: 185-200.
  3. Brossier, R., Operto, S. and Virieux, J., 2009. Seismic imaging of complex onshorestructures by 2D elastic frequency-domain full-waveform inversion. Geophysics,74(6): WCC105-118.
  4. Cao, S.H. and Chen, J.B., 2012. A 17-point scheme and its numerical implementationfor high-accuracy modeling of frequency-domain acoustic equation. Chin. J.Geophys. (in Chinese), 55: 3440-3449.
  5. Chen, J.B., 2013. A generalized optimal 9-point scheme for frequency-domain scalarwave equation. J. Appl. Geophys., 92: 1-7. doi: 10.1016/j.jappgeo.2013.02.008.
  6. Chen, J.B., 2012. An average-derivative optimal scheme for frequency-domain scalarwave equation. Geophysics, 77(6): T201-T210.
  7. Chen, J.B., 2008. Variational integrators and the finite element method. Appl. Math.Comput., 196: 941-958.
  8. Fan, N., Zhao, L.-F. and Xie, X.-B., 2017. A general optimal method for a 2Dfrequency-domain finite-difference solution of scalar wave equation. Geophysics,82(3): T121-T132.
  9. Gauthier, O., Virieux, J. and Tarantola, A., 1986. Two-dimensional nonlinear inversionof seismic waveforms: Numerical results. Geophysics, 51: 1387-1403.
  10. Ha, W. and Shin, C., 2012. Laplace-domain full-waveform inversion of seismic datalacking low-frequency information. Geophysics, 77(5): R199-R206.
  11. Ha, W. and Shin, C., 2013. Why do Laplace-domain waveform inversions yield long-wavelength results? Geophysics, 78(4): R167-R173.
  12. Hustedt, B., Operto, S. and Virieux, J., 2004. Mixed-grid and staggered-grid finite-difference methods for frequency-domain acoustic wave modelling. Geophys. J.Internat., 157: 1269-1296.
  13. Jo, C.H., Suh, J.-H. and Shin, C.S., 1996. An optimal 9-point, finite-difference,frequency-space, 2-D scalar wave extrapolator. Geophysics, 61: 529-537.
  14. Kamei, R., Pratt, R.G. and Tsuji, T., 2015. Misfit functionals in Laplace-Fourier domainwaveform inversion, with application to wide angle ocean bottom seismographdata. Geophys. Prosp., 62: 1054-1074.
  15. Lee, D., Cha, Y.H. and Shin, C., 2008. The direct-removal method of waveforminversion in the Laplace inversion for deep-sea environments. Expanded Abstr.,78th Ann. Internat. SEG Mtg., Las Vegas: 1981-1985.
  16. Liu, L., Liu, H. and Liu, H., 2013. Optimal 15-point finite difference forward modelingin frequency-space domain. Chin. J. Geophys. (in Chinese), 56: 644-652.
  17. Lysmer, J. and Drake, L.A., 1972. A finite element method for seismology. MethodsComput. Phys. Adv. Res. Appl., 11: 181-216.
  18. Marfurt, K.J., 1984. Accuracy of finite-difference and finite-element modeling of thescalar and elastic wave equations. Geophysics, 49: 533-549.
  19. Min, D.-J., Shin, C., Pratt, R.G. and Yoo, H.S., 2003. Weighted-averaging finite-elementmethod for 2D elastic wave equations in the frequency domain. Bull. Seismol.Soc. Am., 93: 904-921.
  20. Operto, S., Virieux, J., Amestoy, P.R., Excellent, J.L., Giraud, L. and Ali, H.B.H., 2007.3D finite-difference frequency-domain modeling of visco-acoustic wavepropagation using a massively parallel direct solver: A feasibility study.Geophysics, 72(5), SM195-SM211.
  21. Operto, S., Virieux, J., Ribodetti, A. and Anderson, J.E., 2009. Finite-differencefrequency-domain modeling of viscoacoustic wave propagation in 2D tiltedtransversely isotropic (TTI) media. Geophysics, 74(5), T75—T95.
  22. Pratt, R.G., 1990. Frequency-domain elastic wave modeling by finite differences: A toolfor crosshole seismic imaging. Geophysics, 55: 626-632.
  23. Pratt, R.G., 1999. Seismic waveform inversion in the frequency domain; Part 1, Theoryand verification in a physical scale model. Geophysics, 64: 888-901.
  24. Pratt, R.G., Shin, C. and Hick, G.J., 1998. Gauss-Newton and full Newton methods infrequency-space seismic waveform inversion. Geophys. J. Internat., 133: 341-
  25. Pratt, R.G. and Worthington, M.H., 1990a. Inverse theory applied to multi-source cross-hole tomography. Part 1: acoustic wave-equation method. Geophys. Prosp., 38:287-310.
  26. Pratt, R.G. and Worthington, M.H., 1990b. Inverse theory applied to multi-source cross-hole tomography. Part 2: elastic wave-equation method. Geophys. Prosp., 38:311-329.
  27. Pyun, S., Shin, C., Lee, H. and Yang, D., 2008. 3D elastic full waveform inversion in the
  28. Laplace domain. Expand. Abstr., 78th Ann. Internat. SEG Mtg., Las Vegas:1976-1980.
  29. Saenger, E.H., Gold, N. and Shapiro, S.A., 2000. Modeling the propagation of elasticwaves using a modified finite-difference grid. Wave Motion, 31: 77-92.
  30. Shin, C. and Cha, Y.H., 2009. Waveform inversion in the Laplace-Fourier domain.Geophys. J. Internat., 177: 1067-1079.
  31. Shin, C. and Cha, Y.H., 2008. Waveform inversion in the Laplace domain. Geophys. J.Internat., 173: 922-931.
  32. Shin, C., Koo, N.H., Cha, H.Y. and Park, K.P., 2010. Sequentially ordered single-frequency 2-D acoustic waveform inversion in the Laplace-Fourier domain.Geophys. J. Internat., 181: 935-950.
  33. Shin, C. and Sohn, H., 1998. A frequency-space 2-D scalar wave extrapolator usingextended 25-point finite-difference operator. Geophysics, 63: 289-296.
  34. Shin, J. Kim, Y., Shin, C. and Calandra, H., 2013. Laplace-domain full waveforminversion using irregular finite elements for complex foothill environments. J._ Appl. Geophys., 96(9): 67 一 76.
  35. Stekl, I. and Pratt, R.G., 1998. Accurate viscoelastic modeling by frequency-domainfinite-differences using rotated operators. Geophysics, 63: 1779-1794.
  36. Tang, X. De, Liu, H., Zhang, H., Liu, L. and Wang, Z.Y., 2015. An adaptable 17-pointscheme for high-accuracy frequency-domain acoustic wave modeling in 2Dconstant density media. Geophysics, 80(6), T211-T221.
  37. Tarantola, A., 1984. Inversion of seismic reflection data in acoustic approximation.Geophysics, 49: 1259-1266.
  38. Virieux, J. and Operto, S., 2009. An overview of full-waveform inversion in explorationgeophysics. Geophysics, 74(6): WCC1-WCC26.
  39. Yang, Q., Hu, G. and Wang, L., 2014. Research status and development trend of fullwaveform inversion. Geophys. Prosp. Pet., 53: 77-83.
  40. Zhang, H., Liu, H., Liu, L., Jin, W.J. and Shi, X.D., 2014. Frequency domain acousticequation high-order modeling based on an average-derivative method. Chin. J.Geophys. (in Chinese), 57: 1599-1612.
Share
Back to top
Journal of Seismic Exploration, Print ISSN: 0963-0651, Published by AccScience Publishing
We use cookies on our website to ensure you get the best experience.
Read more about our cookies here
Accept