f Implementation of anelasticity to an integral equation method for seismic modelling

The introduction of attenuation to direct methods for seismic modelling such as finite differences or the Fourier method is still quite difficult. Recently an integral equation method has been presented which allows an easy incorporation of attenuation mechanisms (Freter and Kummer, 1989, Freter, 1990). This method is based on the calculation of the Green's function associated to a heterogenous medium by solving a Fredholm integral equation of the second kind. It is sufficiently general as it allows the implementation of various models of viscoacoustic media with an arbitrary dependence of the quality factor on both space and frequency (Freter, Schrader and Behle, 1991).