Modeling and Visualizing Seismic Wave Propagation in Elastic Medium Using Multi-Dimension Wave Digital Filtering Approach

Authors: Jason Chien-Hsun Tseng, Nguyen Dong-Thai Dao, Chong-Ching Chang

Abstract:

A novel PDE solver using the multidimensional wave digital filtering (MDWDF) technique to achieve the solution of a 2D seismic wave system is presented. In essence, the continuous physical system served by a linear Kirchhoff circuit is transformed to an equivalent discrete dynamic system implemented by a MD wave digital filtering (MDWDF) circuit. This amounts to numerically approximating the differential equations used to describe elements of a MD passive electronic circuit by a grid-based difference equations implemented by the so-called state quantities within the passive MDWDF circuit. So the digital model can track the wave field on a dense 3D grid of points. Details about how to transform the continuous system into a desired discrete passive system are addressed. In addition, initial and boundary conditions are properly embedded into the MDWDF circuit in terms of state quantities. Graphic results have clearly demonstrated some physical effects of seismic wave (P-wave and S–wave) propagation including radiation, reflection, and refraction from and across the hard boundaries. Comparison between the MDWDF technique and the finite difference time domain (FDTD) approach is also made in terms of the computational efficiency.

Keywords: Seismic Wave Propagation, Multi-dimension WaveDigital Filters, Partial Differential Equations.

Digital Object Identifier (DOI): doi.org/10.5281/zenodo.1333190

References:

[1] O. Julius Smith III's, A Basic Introduction to Digital Waveguide Synthesis, CCRMA, Stanford University, 2006.
[2] A. Fettweis and G. Nitsche, "Numerical integration of partial differential equations using principles of multidimensional wave digital filters", J. of VLSI Sig. Proc., vol. 3, pp. 7-24, 1991.
[3] A. Fettweis, "Wave digital filter: theory and practice", Proc. IEEE, vol. 74, no. 2, pp. 270-327, 1986.
[4] A. Fettweis, ``Robust numerical integration using wave-digital concepts," Multidim. Syst. Sig. Process, Vol. 17, pp. 7-25, 2006.
[5] A. Udias, Principles of seismology, Cambridge University Press, 1999.
[6] L. Dresen "The influence of an asymmetry in the sequence rock-coalrock on the propagation of Rayleigh seam waves", Geophysical Prospecting, 1985.
[7] P. K. Kundu, Fluid Mechanics, Academic Press Limited, London, 1990.
[8] J. P. Berenger, "A perfectly matched layer for the absorption of electromagnetic waves-. J. Computational Physics, vol 114, 185-200, 1994.