A 21-point finite difference scheme for 2D frequency-domain elastic wave modelling
Bingluo Gu 1 2 Guanghe Liang 1 3 Zhiyuan Li 1 2
+ Author Affiliations
- Author Affiliations
1 Key Laboratory of Mineral Resources, Institute of Geology and Geophysics, The Chinese Academy of Sciences, Beijing 100029, China.
2 College of Earth Sciences, University of The Chinese Academy of Sciences, Beijing 100049, China.
3 Corresponding author. Email: lgh@mail.iggcas.ac.cn
Exploration Geophysics 44(3) 156-166 https://doi.org/10.1071/EG12064
Submitted: 22 October 2012 Accepted: 7 April 2013 Published: 8 May 2013
Abstract
The 21-point finite difference scheme for the frequency-space domain elastic wave forward modelling is designed through optimising the impedance matrix, especially calculating the spatial derivative terms and the mass acceleration terms of the elastic wave displacement equation as accurately as possible. Comparative tests show that the 21-point finite difference scheme is much better in grid dispersion, memory requirement, and computation time than the 9-point scheme and slightly better than the 25-point scheme. The 21-point finite difference scheme is ~15% lower in memory consumption and computing time than the 25-point scheme. The numerical examples show that the 21-point finite difference scheme is valid in the sense of the numerical simulation of ideal models.
Key words: 21-point finite difference scheme, computing time, dispersion, frequency domain, memory consumption.
References
Berenger, J., 1994, A perfectly matched layer for the absorption of electromagnetic waves: Journal of Computational Physics, 114, 185–200| A perfectly matched layer for the absorption of electromagnetic waves:Crossref | GoogleScholarGoogle Scholar |
Hustedt, B., Operto, S., and Virieus, J., 2004, Mixed-grid and staggered-grid finite-difference methods for frequency-domain acoustic wave modelling: Geophysical Journal International, 157, 1269–1296
| Mixed-grid and staggered-grid finite-difference methods for frequency-domain acoustic wave modelling:Crossref | GoogleScholarGoogle Scholar |
Jo, C., Shin, C., and Suh, J., 1996, An optimal 9-point finite-difference frequency-space 2-D scalar wave extrapolator: Geophysics, 61, 529–537
| An optimal 9-point finite-difference frequency-space 2-D scalar wave extrapolator:Crossref | GoogleScholarGoogle Scholar |
Li, G., Feng, J., and Zhu, G., 2011, Quasi-P wave forward modelling in viscoelastic VTI media in frequency-space domain: Chinese Journal of Geophysics, 54, 200–207
Liao, J. P., Wang, H. Z., and Ma, Z. T., 2009a, Compression storage for frequency-space domain 2-D elastic wave forward modelling: CPS/SEG Beijing 2009 International Geophysical Conference & Exposition.
Liao, J. P., Wang, H. Z., and Ma, Z. T., 2009b, 2-D elastic wave modelling with frequency-space 25-point finite-difference operators: Applied Geophysics, 6, 259–266
| 2-D elastic wave modelling with frequency-space 25-point finite-difference operators:Crossref | GoogleScholarGoogle Scholar |
Liao, J. P., Liu, H. X., Wang, H. Z., Yang, T. C., Wang, Q. R., Liu, X. H., and Ma, Z. T., 2011, Study on rapid highly accurate acoustic wave numerical simulation in frequency space domain: Progress in Geophysics, 26, 1359–1363
Luo, Y., and Schuster, G., 1990, Parsimonious staggered grid finite differencing of the wave equation: Geophysical Research Letters, 17, 155–158
| Parsimonious staggered grid finite differencing of the wave equation:Crossref | GoogleScholarGoogle Scholar |
Lysmer, J., and Drake, L. A., 1972, A finite-element method for seismology, in B. A. Bolt, ed., Methods in computational physics: Academic Press Inc. 11, 181–216.
Marfurt, K., 1984, Accuracy of finite-difference and finite-element modelling of the scalar and elastic wave equations: Geophysics, 49, 533–549
| Accuracy of finite-difference and finite-element modelling of the scalar and elastic wave equations:Crossref | GoogleScholarGoogle Scholar |
Marfurt, K., and Shin, C., 1989, The future of iterative modelling in geophysical exploration, in E. Eisner, ed., Handbook of geophysical exploration: I - Seismic exploration, 21 - Supercomputers in seismic exploration: Pergamon Press, 203–228.
Min, D., Shin, C., Kwon, B., and Chung, S., 2000, Improved frequency-domain elastic wave modelling using weighted-averaging difference operators: Geophysics, 65, 884–895
| Improved frequency-domain elastic wave modelling using weighted-averaging difference operators:Crossref | GoogleScholarGoogle Scholar |
Min, D., Shin, C., Pratt, R., and Yoo, H., 2003, Weighted-averaging finite-element method for 2D elastic wave equations in the frequency domain: Bulletin of the Seismological Society of America, 93, 904–921
| Weighted-averaging finite-element method for 2D elastic wave equations in the frequency domain:Crossref | GoogleScholarGoogle Scholar |
Min, D., Shin, C., and Yoo, H., 2004, Free-surface boundary condition in finite-difference elastic wave modelling: Bulletin of the Seismological Society of America, 94, 237–250
| Free-surface boundary condition in finite-difference elastic wave modelling:Crossref | GoogleScholarGoogle Scholar |
Operto, S., Virieux, J., Ribodetti, A., and Anderson, J., 2009, Finite-difference frequency-domain modelling of viscoacoustic wave propagation in 2D tilted transversely isotropic (TTI) media: Geophysics, 74, T75–T95
| Finite-difference frequency-domain modelling of viscoacoustic wave propagation in 2D tilted transversely isotropic (TTI) media:Crossref | GoogleScholarGoogle Scholar |
Paffenholz, J., Mclain, B., Zaske, J., and Keliher, P., 2002, Subsalt multiple attenuation and imaging: observations from the Sigsbee2B synthetic data set: 72nd Annual International Meeting, SEG, Expanded Abstracts, 2122–2125.
Pratt, R., 1990, Frequency-domain elastic wave modelling by finite differences: a tool for crosshole seismic imaging: Geophysics, 55, 626–632
| Frequency-domain elastic wave modelling by finite differences: a tool for crosshole seismic imaging:Crossref | GoogleScholarGoogle Scholar |
Pratt, R., and Worthington, M., 1990, Inverse theory applied to multi-source cross-hole tomography, Part 1: acoustic wave-equation method: Geophysical Prospecting, 38, 287–310
| Inverse theory applied to multi-source cross-hole tomography, Part 1: acoustic wave-equation method:Crossref | GoogleScholarGoogle Scholar |
Qian, J., Wu, S. G., and Cui, R. F., 2013, Extension of split perfectly matched absorbing layer for 2D wave propagation in porous transversely isotropic media: Exploration Geophysics, 44, 25–30
| Extension of split perfectly matched absorbing layer for 2D wave propagation in porous transversely isotropic media:Crossref | GoogleScholarGoogle Scholar |
Ren, H., Wang, H., and Gong, T., 2009, Seismic modelling of scalar seismic wave propagation with finite-difference scheme in frequency-space domain: Geophysical Prospecting for Petroleum, 48, 20–26
Saenger, E., Gold, N., and Shapiro, S., 2000, Modelling the propagation of elastic waves using a modified finite-difference grid: Wave Motion, 31, 77–92
| Modelling the propagation of elastic waves using a modified finite-difference grid:Crossref | GoogleScholarGoogle Scholar |
Shin, C., and Sohn, H., 1998, A frequency-space 2-D scalar wave extrapolator using extended 25-point finite-difference operators: Geophysics, 63, 289–296
| A frequency-space 2-D scalar wave extrapolator using extended 25-point finite-difference operators:Crossref | GoogleScholarGoogle Scholar |
Štekl, I., and Pratt, R., 1998, Accurate viscoelastic modelling by frequency-domain finite differences using rotated operators: Geophysics, 63, 1779–1794
| Accurate viscoelastic modelling by frequency-domain finite differences using rotated operators:Crossref | GoogleScholarGoogle Scholar |
Yin, W., Yin, X. Y., and Wu, G. C., 2006, The method of finite difference of high precision elastic wave equations in the frequency domain and wave-field simulation: Chinese Journal of Geophysics, 49, 561–568