Abstract: Starting from the governing equations for axisymmetric consolidation of poroelastic materials with anisotropic permeability and compressible fluid, and using the technique of Laplace-Hankel transform, the generalized expressions for the displacements and stresses of materials in transformed domain are derived. Based on the theory of the analytical layer-element method, the global stiffness is assembled by imposing the boundary conditions and continuity between layers on the problem. Applying the inversion of the Laplace-Hankel transform, the solution to the consolidation problem of multi-layered foundation due to well point sink is obtained. In order to verify the proposed theory, numerical results are carried out by means of the compiled procedure and compared with the existing references. The numerical results demonstrate that the bigger the velocity of well point sink, the more quickly the surface displacement, and the feature of compressible fluid can decease the development of surface displacement under the same time factor.