Abstract：Based on Biot theory, the wave propagation in saturated porous media is coupling issue. Because of the complexity of the model, it's difficult to give high-efficient numeric solution. For that reason, the paper gave a practical solution that the seismic wave propagates in viscoelastic saturated layered and porous media, which is to use Laplace transform in time domain and Fourier transform in space domain for studied 3-D layered media and to transform the original issue into ordinary differential equation set with six independent variables, and each variable is function only having depth coordinate and one horizontal coordinate, and the transmission matrix of wave propagation issue in each layer can be further got, using numeric approach to compute the eigenvalues and the eigenvectors of transmission matrix. The kernel of the algorithm is first to get the wavefield in post-transform space, and then to use FFT to carry out inverse Fourier and Laplace transforms and to get wavefield in time domain. The paper gave the detailed deduce of algorithm and numeric cases.
