本文在复频域内,通过应用混合变量粘弹性波方程和线性常微分方程组的指数矩阵解法,给出了一种计算非均匀吸收介质中地震波传播的广义传播矩阵解法。该方法适用于各种粘弹性模型,可模拟任意震源及所产生的各种体波、面波,数值结果表明具有很高的计算精度。 In the complex frequency domain, a generalized propagation matrix method for seismic wave propagation in non-uniform absorption media is given by applying the mixed variable viscoelastic wave equation and the exponential matrix solution of linear ordinary differential equations. The method is applicable to all kinds of viscoelastic models and can simulate any seismic source and various kinds of body waves and surface waves. The numerical results show that the method has high computational accuracy.