Register      Login
Exploration Geophysics Exploration Geophysics Society
Journal of the Australian Society of Exploration Geophysicists
RESEARCH ARTICLE

Least squares staggered-grid finite-difference for elastic wave modelling

Lei Yang 1 2 Hongyong Yan 1 3 Hong Liu 1
+ Author Affiliations
- Author Affiliations

1 Key Laboratory of Petroleum Resources Research, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing, 100029, China.

2 University of Chinese Academy of Sciences, Beijing, 100049, China.

3 Corresponding author. Email: yanhongyong@163.com

Exploration Geophysics 45(4) 255-260 https://doi.org/10.1071/EG13087
Submitted: 17 October 2013  Accepted: 14 January 2014   Published: 24 February 2014

Abstract

Staggered-grid finite-difference (SFD) methods have been used widely in seismic wave numerical modelling and migration. The conventional way to calculate the high-order SFD coefficients on spatial derivatives is the Taylor-series expansion method, which generally leads to great accuracy on just a small frequency zone. In this paper, we first derive the SFD coefficients of arbitrary even-order accuracy for the first-order spatial derivatives by the dispersion relation and the least-squares method, which can satisfy the specific numerical solution accuracy of the derivative on a wide frequency zone. Then we use the SFD coefficients based on the least-squares to solve the first-order spatial derivatives and analyse the accuracy of the numerical solution. Finally, we perform elastic wave numerical modelling with the least-squares staggered-grid finite-difference (LSSFD) method. Meanwhile, the numerical dispersion, the modelling accuracy and the computing costs of the new method are compared with that of the Taylor-series expansion staggered-grid finite-difference (TESFD) method. The numerical dispersion analysis and elastic wavefield modelling results demonstrate that the LSSFD method can efficiently suppress the numerical dispersion and has greater modelling accuracy than the conventional TESFD method under the same discretisations and without extra computing costs.

Key words: dispersion, elastic wave, least squares, numerical modelling, staggered-grid finite-difference (SFD).


References

Bérenger, J. P., 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 |

Chu, C., and Stoffa, P. L., 2012a, Implicit finite-difference simulations of seismic wave propagation: Geophysics, 77, T57–T67
Implicit finite-difference simulations of seismic wave propagation:Crossref | GoogleScholarGoogle Scholar |

Chu, C., and Stoffa, P. L., 2012b, Determination of finite-difference weights using scaled binomial windows: Geophysics, 77, W17–W26
Determination of finite-difference weights using scaled binomial windows:Crossref | GoogleScholarGoogle Scholar |

Dong, L. G., Ma, Z. T., Cao, J. Z., Wang, H. Z., Gao, J. H., Lei, B., and Xu, S. Y., 2000a, A staggered-grid high-order difference method of one-order elastic wave equation: Chinese Journal of Geophysics, 43, 411–419

Dong, L. G., Ma, Z. T., and Cao, J. Z., 2000b, A study on stability of the staggered-grid high-order difference method of first-order elastic wave equation: Chinese Journal of Geophysics, 43, 856–864

Etgen, J. T., 2007, A tutorial on optimizing time domain finite-difference schemes: ‘Beyond Holberg’: Stanford Exploration Project Report, 129, 33–43.

Graves, R. W., 1996, Simulating seismic wave propagation in 3D elastic media using staggered-grid finite difference: Bulletin of the Seismological Society of America, 86, 1091–1106

Holberg, O., 1987, Computational aspects of the choice of operator and sampling interval for numerical differentiation in large-scale simulation of wave phenomena: Geophysical Prospecting, 35, 629–655
Computational aspects of the choice of operator and sampling interval for numerical differentiation in large-scale simulation of wave phenomena:Crossref | GoogleScholarGoogle Scholar |

Igel, H., Riollet, B., and Mora, P., 1992, Accuracy of staggered 3-D finite-difference grids for anisotropic wave propagation: 62nd Annual International Meeting, SEG, Expanded Abstracts, 1244–1246.

Kim, J. W., and Lee, D. J., 1996, Optimized compact finite difference schemes with maximum resolution: AIAA Journal, 34, 887–893
Optimized compact finite difference schemes with maximum resolution:Crossref | GoogleScholarGoogle Scholar |

Kindelan, M., Kamel, A., and Sguazzero, P., 1990, On the construction and efficiency of staggered numerical differentiators for the wave equation: Geophysics, 55, 107–110
On the construction and efficiency of staggered numerical differentiators for the wave equation:Crossref | GoogleScholarGoogle Scholar |

Lee, C., and Seo, Y., 2002, A new compact spectral scheme for turbulence simulations: Journal of Computational Physics, 183, 438–469
A new compact spectral scheme for turbulence simulations:Crossref | GoogleScholarGoogle Scholar |

Liu, Y., 2013, Globally optimal finite-difference schemes based on least squares: Geophysics, 78, T113–T132
Globally optimal finite-difference schemes based on least squares:Crossref | GoogleScholarGoogle Scholar |

Liu, Y., and Sen, M. K., 2009a, A practical implicit finite-difference method: examples from seismic modeling: Journal of Geophysics and Engineering, 6, 231–249
A practical implicit finite-difference method: examples from seismic modeling:Crossref | GoogleScholarGoogle Scholar | 1:CAS:528:DC%2BD1MXivFeitbk%3D&md5=ba1dd9481927d9f8b7f4c928eea65ef5CAS |

Liu, Y., and Sen, M. K., 2009b, An implicit staggered-grid finite-difference method for seismic modeling: Geophysical Journal International, 179, 459–474
An implicit staggered-grid finite-difference method for seismic modeling:Crossref | GoogleScholarGoogle Scholar |

Liu, Y., and Sen, M. K., 2011, Scalar wave equation modeling with time-space domain dispersion-relation-based staggered-grid finite-difference schemes: Bulletin of the Seismological Society of America, 101, 141–159
Scalar wave equation modeling with time-space domain dispersion-relation-based staggered-grid finite-difference schemes:Crossref | GoogleScholarGoogle Scholar |

Madariaga, R., 1976, Dynamics of an expanding circular fault: Bulletin of the Seismological Society of America, 66, 639–666

Nakata, N., Tsuji, T., and Matsuoka, T., 2011, Acceleration of computation speed for elastic wave simulation using a Graphic Processing Unit: Exploration Geophysics, 42, 98–104
Acceleration of computation speed for elastic wave simulation using a Graphic Processing Unit:Crossref | GoogleScholarGoogle Scholar | 1:CAS:528:DC%2BC3MXisVKrtL4%3D&md5=5b847e5a7f11d813dcd9698790c77587CAS |

Pei, Z., 2004, Numerical modeling using staggered-grid high order finite difference of elastic wave equation on arbitrary relief surface: Oil Geophysical Prospecting, 39, 629–634

Tam, C. K. W., and Webb, J. C., 1993, Dispersion-relation-preserving finite difference schemes for computational acoustics: Journal of Computational Physics, 107, 262–281
Dispersion-relation-preserving finite difference schemes for computational acoustics:Crossref | GoogleScholarGoogle Scholar |

Tessmer, E., 2011, Using the rapid expansion method for accurate time-stepping in modeling and reverse-time migration: Geophysics, 76, S177–S185
Using the rapid expansion method for accurate time-stepping in modeling and reverse-time migration:Crossref | GoogleScholarGoogle Scholar |

Virieux, J., 1984, SH-wave propagation in heterogeneous media: velocity stress finite-difference method: Geophysics, 49, 1933–1942
SH-wave propagation in heterogeneous media: velocity stress finite-difference method:Crossref | GoogleScholarGoogle Scholar |

Virieux, J., 1986, P-SV wave propagation in heterogeneous media: velocity stress finite difference method: Geophysics, 51, 889–901
P-SV wave propagation in heterogeneous media: velocity stress finite difference method:Crossref | GoogleScholarGoogle Scholar |

Yan, H., and Liu, Y., 2012, High-order finite-difference numerical modeling of wave propagation in viscoelastic TTI media using rotated staggered grid: Chinese Journal of Geophysics, 55, 252–265
High-order finite-difference numerical modeling of wave propagation in viscoelastic TTI media using rotated staggered grid:Crossref | GoogleScholarGoogle Scholar |

Yan, H., Liu, Y., and Zhang, H., 2013, Prestack reverse-time migration with a time-space domain adaptive high-order staggered-grid finite-difference method: Exploration Geophysics, 44, 77–86
Prestack reverse-time migration with a time-space domain adaptive high-order staggered-grid finite-difference method:Crossref | GoogleScholarGoogle Scholar |

Zhang, J., and Yao, Z., 2013, Optimized finite-difference operator for broadband seismic wave modeling: Geophysics, 78, A13–A18
Optimized finite-difference operator for broadband seismic wave modeling:Crossref | GoogleScholarGoogle Scholar |

Zhou, H., and Zhang, G., 2011, Prefactored optimized compact finite-difference schemes for second spatial derivatives: Geophysics, 76, WB87–WB95
Prefactored optimized compact finite-difference schemes for second spatial derivatives:Crossref | GoogleScholarGoogle Scholar |