Abstract: Based on the Biot's theory of porous media,the calculation formula of reverberation ray matrix method is established for one-dimensional transient response of fluid-saturated soil,where the non-homogeneous,inertial,viscous and the compressible of solid particles and fluid are taken into account.The present methodology is validated by comparing solutions with some known result.As numerical examples,assuming that the material properties of the saturated soil have an exponential law distribution along the thickness-coordinate,the transient response,in terms of displacement,stress and pore press,are examined using numerical inverse Laplace transform.The effect of non-homogeneous on transient responses of gradient saturated soil is demonstrated and discussed.