Abstract: This paper presents the solutions for three dimensional Biot's consolidation problem of multilayered soils, and the compressibility of pore fluid is considered. Based on the equations of three dimensional Biot's consolidation with compressible constituents, by introducing displacement functions, intermediate variables and the Laplace-Fourier transform, two sets of transfer matrix of 6 order and 2 order are obtained. The solutions for three dimensional Biot's consolidation problem of multilayered soils with compressible constituents in the transform domain are obtained by using the transfer matrix method and combined with the continuity conditions and the boundary conditions. The actual solutions can be acquired by inverting the Laplace-Fourier transform. Corresponding program is carried out to analyze the influence of the compressibility of pore fluid on the consolidation behavior of the soils.