S U M M A R YThis paper presents a generalized wave equation which unifies viscoelastic and pure elastic cases into a single wave equation. In the generalized wave equation, the degree of viscoelasticity varies between zero and unity, and is defined by a controlling parameter. When this viscoelastic controlling parameter equals to 0, the viscous property vanishes and the generalized wave equation becomes a pure elastic wave equation. When this viscoelastic controlling parameter equals to 1, it is the Stokes equation made up of a stack of pure elastic and Newtonian viscous models. Given this generalized wave equation, an analytical solution is derived explicitly in terms of the attenuation and the velocity dispersion. It is proved that, for any given value of the viscoelastic controlling parameter, the attenuation component of this generalized wave equation perfectly satisfies the power laws of frequency. Since the power laws are the fundamental characteristics in physical observations, this generalized wave equation can well represent seismic wave propagation through subsurface media.Key words: Non-linear differential equation; Seismic attenuation; Wave propagation.
I N T RO D U C T I O NFor seismic wave simulation, there are two basic types of wave equations: pure elastic wave equation and viscoelastic wave equation. This paper proposes to unify these two types of wave equations into a single wave equation.The seismic property of subsurface media should generally be considered to be a superposition of Hook's pure elastic model and a viscoelastic model. Customarily, the Kelvin-Voigt model assumes the seismic property to be a stack of the elastic model and a Newtonian viscous model. The latter is a linear viscosity in which the stress-strain relationship is presented by a first-order temporal derivative. Thus, the Kelvin-Voigt model is given as the following:where σ (t) is the time-dependent tensile stress, ε(t) is the corresponding strain and E 0 and E 1 are two parameters of the model. In eq.(1), E 0 ε is the elastic term and E 1 dε/dt is the linear viscous term. This leads to a well-known wave equation called the Stokes equation. However, subsurface media behave in a much generalized viscoelastic manner. To reflect this characteristic, a stress-strain relationship can be described by a fractional order, instead of the first order, temporal derivative. Then, the superposition model can be expressed aswhere β is the fractional order of the time differential. The last term E 1 d β ε/dt β uses a compact form of fractional derivative to describe the frequency and time dependency of a viscoelastic system (Scott-Blair 1947).The mathematical-physics history of this empirical method, using a fractional derivative to represent the viscoelastic property, begins with Nutting (1921) who reported the observation that the stress relaxation could be modelled by fractional powers of time. Such a model is in sharp contrast to the traditional view that stress relaxation is best described by decaying exponentials. Gemant (1936) f...