Abstract
Using staggered-grid finite difference method to solve seismic wave equation, large spatial grid and high dominant frequency of source cause numerical dispersion, staggered-grid finite difference method, which can reduce the step spatial size and increase the order of difference, will multiply the calculation amount and reduce the efficiency of solving wave equationThe optimal nearly analytic discrete (ONAD) method can accurately solve the wave equation by using the combination of displacement and gradient of spatial nodes to approach the spatial partial derivative under rough grid and high-frequency condition. In this study, the ONAD method is introduced into the field of reverse-time migration (RTM) for performing forward- and reverse-time extrapolation of a two-dimensional acoustic equation, and the RTM based on ONAD method is realized via normalized cross-correlation imaging condition, effectively suppressed the numerical dispersion and improved the imaging accuracy. Using ONAD method to image the groove model and SEG/EAGE salt dome model by RTM, and comparing with the migration sections obtained by staggered-grid finite difference method with the same time order 2nd and space order 4th, results show that the RTM based on ONAD method can effectively suppress numerical dispersion caused by the high frequency components in source and shot records, and archive accurate imaging of complex geological structures especially the fine structure, and the migration sections of the measured data show that ONAD method has practical application value.
Similar content being viewed by others
References
Alterman, Z., and Karal, F. C., 1968, Propagation of elastic wave in layered media by finite difference methods: Bulletin of the Seismological Society of America, 58(1), 67–398.
Chen, K. Y., 2009, Wave equation pre-stack reversetime migration method based on high-order finite-difference approach: Geophysical Prospecting for Petroleum, 73(5), 475–478.
Chen, K. Y., 2010, Study on absorbing boundary conditions of perfectly matched layer: Geophysical Prospecting for Petroleum, 49(5), 472–477.
Chen, K. Y., 2015, Numerical dispersion mechanism analysis and its wave field characteristic of seismic wave numerical propagation: Petroleum Geophysics, 13(3), 7–10.
Dong, L. G., Ma, Z. T., Cao, J. Z., et al., 2000, A Staggered-grid high-order difference method of one-order elastic wave equation: Chinese Journal of Geophysics, 43(3), 411–419.
Dong, L. G., and Li, P. M., 2004, Dispersive Problem in seismic wave propagation numerical modeling: Natural Gas Industry, 24(6), 53–56.
D’Hooge, J., and Mcaleavey, S. A., 2011, Medical Imaging 2010: Ultrasonic Imaging, Tomography, and Therapy: Proceedings of SPIE - The International Society for Optical Engineering, 18(2), S180.
Fei, T., and Larner, K., 1995, Elimination of numerical dispersion in finite-difference modeling and migration by flux-corrected transport: Geophysics, 60(6), 1830–1842.
Kaelin, A., and Biondi, B., 2006, Least-square attenuation of reverse time migration artifacts: 76th Annual International Meeting, SEG, Expanded Abstracts, 2348–2352.
Kondoh, Y., Hosaka, Y., and Ishii, K., 2010, Kernel optimum nearly-analytical discretization (KOND) algorithm applied to parabolic and hyperbolic equations: Computers & Mathematics with Applications, 27(27), 59–90.
Liu, Y., 2013, Noise reduction by vector median filtering: Geophysics, 78(3), V79–V86.
Li, J. S., Yang, D. H., and Liu, F. Q., 2013, An efficient reverse-time migration method using local nearly analytic discrete operator: Geophysics, 78(1): s15–s23.
Li, J. S., Michael, F., Yang, D. H., et al., 2015, 3D weak-dispersion reverse time migration using a stereo-modeling operator: Geophysics, 80(1), s19–s30.
Li, J. S., Yang, D. H., Michael, F., et al., 2018, Using a stereo-modeling operator to reduce aliasing during migration, Geophysics, 83(5), 1–41.
Mou, Y. G., and Pei, Z. L., 2005, Seismic numerical simulation of 3D complex media: Petroleum Industry Press, China, 14–15.
Tang, H. G., and He, B. S., 2017, Pseudo spectrum method of first-order acoustic wave equation finite-difference schemes with fourth-order time difference accuracy: Oil Geophysical Prospecting, 52(1), 1–80.
Virieux, J., 1986, P-SV wave propagation in heterogeneous media: Velocity-stress finite-difference method: Geophysics, 51(4), 889–901.
Wang, Y., Duan, Y. W., Wang, T., et al., 2017, Simulation of 2D elastic wave field separation by optimizing approximate analytical discretization method: Oil Geophysical Prospecting, 52(3), 458–467.
Xie, W., Yang, D. H., Liu, F. Q., et al., 2014, Reverse-time migration in acoustic VTI media using a high-order stereo operator: Geophysics, 79(3), WA3–WA11.
Yang, D. H., Lu, M., Wu, R. S., et al., 2004. An optimal nearly analytic discrete method for 2D acoustic and elastic wave equations: Bulletin of the Seismological Society of America, 94(5), 1982–1991.
Yang, D. H., Peng, J. M., Lu, M., et al., 2006, Optimal Nearly Analytic Discrete Approximation to the Scalar Wave Equation: Bulletin of the Seismological Society of America, 96(3), 1114–1130.
Yang, D. H., Song, G., and Lu, M., 2007, Optimally Accurate Nearly Analytic Discrete Scheme for Wave-Field Simulation in 3D Anisotropic Media: Bulletin of the Seismological Society of America, 97(5), 1557–1569.
Yang, D., Liu, E., Song, G., and Wang, N., 2009, Elastic wave modelling method based on the displacement–velocity fields: an improving nearly analytic discrete approximation: Journal of Seismology, 13(2), 209–217.
Yang, D., Tong, P., and Deng, X., 2012, A central difference method with low numerical dispersion for solving the scalar wave equation: Geophysical Prospecting, 60(5), 885–905.
Yao, G., Wu, D., and Debens, H. A., 2016, Adaptive finite difference for seismic wavefield modelling in acoustic media: Scientific Report, 6, 30302.
Zhang, H. X., and Ning, S. N., 2002, Pre-stack reverse time migration of elastic wave equation: Journal of China University of Mining & Technology, 31(5), 371–375.
Zhang, J. H., and Yao, Z. X., 2013a, Optimized explicit finite-difference schemes for spatial derivatives using maximum norm: Journal of Computational Physics, 250(10):511–526.
Zhang, J. H., and Yao, Z. X., 2013b, Optimized finite-difference operator for broadband seismic wave modeling: Geophysics, 78(1), A13–A18.
Acknowledgement
Thanks to Li Li-Ping, master of Ocean University of China, for his help in this research.
Author information
Authors and Affiliations
Corresponding author
Additional information
This work was financially supported by the National Key R&D Program of China (No. 2018YFC1405900), the National Natural Science Foundation of China (No. 41674118), the Fundamental Research Funds for the Central Universities (No. 201822011) and the National Science and Technology Major Project (No. 2016ZX05027-002).
Liu Ming-Zhu, graduated from Anhui University of Science and Technology in 2017 with a bachelor’s degree in exploration technology and engineering. She is currently a master’s degree student at Ocean University of China. Her main research interest is the study of seismic wave RTM.
Rights and permissions
About this article
Cite this article
Liu, MZ., He, BS. Suppress numerical dispersion in reverse-time migration of acoustic wave equation using optimal nearly analytic discrete method. Appl. Geophys. 17, 133–142 (2020). https://doi.org/10.1007/s11770-019-0790-1
Received:
Revised:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11770-019-0790-1