Numerical simulation of attenuated wavefields using a Padé approximant method

• Published in: Geophysical Journal International Vol. 78; no. 1; pp. 105 - 118
• Main Authors: Day, Steven M., Minster, J. Bernard
• Format: Journal Article
• Language: English
• Published: Oxford, UK Blackwell Publishing Ltd 01.07.1984
• Subjects: Approximants; Attenuation; Computer simulation; Exact solutions; Laws; Mathematical analysis; Mathematical models; Operators
• ISSN: 0956-540X; 0016-8009; 1365-246X
• DOI: 10.1111/j.1365-246X.1984.tb06474.x

Realistic anelastic attenuation laws are usually formulated as convolution operators, but this representation is intractable for time-domain synthetic seismogram methods such as the finite difference method. An approach based on Padé approximants provides a convenient, accurate reformulation of general anelastic laws in differential form. The resulting differential operators form a uniformly convergent sequence of increasing order in the time derivative, and all are shown to be causal, stable and dissipative. In the special case of frequency-independent Q, all required coefficients for the operators are obtained in closed form in terms of Legendre polynomials. Low-order approximants are surprisingly accurate. Finite-difference impulse responses for a plane wave in a constant-Q medium, calculated with the fifth-order convergent, are virtually indistinguishable from the exact solution. The formulation is easily generalized to non-scalar waves. Moreover, this method provides a framework for incorporating amplitude-dependent attenuation into numerical simulations.