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

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.
- Alford, R.M., Kelly, K.R. and Boore, D.M., 1974. Accuracy of finite-differencemodeling of the acoustic wave equation. Geophysics, 39: 834-842.
- Bérenger, J.P., 1994. A perfectly matched layer for the absorption of electromagneticwaves. J. Comput. Phys., 114: 185-200.
- 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.
- 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.
- 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.
- Chen, J.B., 2012. An average-derivative optimal scheme for frequency-domain scalarwave equation. Geophysics, 77(6): T201-T210.
- Chen, J.B., 2008. Variational integrators and the finite element method. Appl. Math.Comput., 196: 941-958.
- 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.
- Gauthier, O., Virieux, J. and Tarantola, A., 1986. Two-dimensional nonlinear inversionof seismic waveforms: Numerical results. Geophysics, 51: 1387-1403.
- Ha, W. and Shin, C., 2012. Laplace-domain full-waveform inversion of seismic datalacking low-frequency information. Geophysics, 77(5): R199-R206.
- Ha, W. and Shin, C., 2013. Why do Laplace-domain waveform inversions yield long-wavelength results? Geophysics, 78(4): R167-R173.
- 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.
- 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.
- 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.
- 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.
- 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.
- Lysmer, J. and Drake, L.A., 1972. A finite element method for seismology. MethodsComput. Phys. Adv. Res. Appl., 11: 181-216.
- Marfurt, K.J., 1984. Accuracy of finite-difference and finite-element modeling of thescalar and elastic wave equations. Geophysics, 49: 533-549.
- 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.
- 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.
- 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.
- Pratt, R.G., 1990. Frequency-domain elastic wave modeling by finite differences: A toolfor crosshole seismic imaging. Geophysics, 55: 626-632.
- Pratt, R.G., 1999. Seismic waveform inversion in the frequency domain; Part 1, Theoryand verification in a physical scale model. Geophysics, 64: 888-901.
- 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-
- 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.
- 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.
- Pyun, S., Shin, C., Lee, H. and Yang, D., 2008. 3D elastic full waveform inversion in the
- Laplace domain. Expand. Abstr., 78th Ann. Internat. SEG Mtg., Las Vegas:1976-1980.
- 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.
- Shin, C. and Cha, Y.H., 2009. Waveform inversion in the Laplace-Fourier domain.Geophys. J. Internat., 177: 1067-1079.
- Shin, C. and Cha, Y.H., 2008. Waveform inversion in the Laplace domain. Geophys. J.Internat., 173: 922-931.
- 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.
- Shin, C. and Sohn, H., 1998. A frequency-space 2-D scalar wave extrapolator usingextended 25-point finite-difference operator. Geophysics, 63: 289-296.
- 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.
- Stekl, I. and Pratt, R.G., 1998. Accurate viscoelastic modeling by frequency-domainfinite-differences using rotated operators. Geophysics, 63: 1779-1794.
- 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.
- Tarantola, A., 1984. Inversion of seismic reflection data in acoustic approximation.Geophysics, 49: 1259-1266.
- Virieux, J. and Operto, S., 2009. An overview of full-waveform inversion in explorationgeophysics. Geophysics, 74(6): WCC1-WCC26.
- Yang, Q., Hu, G. and Wang, L., 2014. Research status and development trend of fullwaveform inversion. Geophys. Prosp. Pet., 53: 77-83.
- 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.