The Biot linearized quasi-static theory of fluid-infiltrated porous materials is used to formulate the problem of the two-dimensional plane strain deformation of a multi-layered poroelastic half-space by surface loads. The Fourier-Laplace transforms of the stresses, displacements, pore pressure and fluid flux in each homogeneous layer of the multi-layered half-space are expressed in terms of six arbitrary constants. Generalized Thomson-Haskell matrix method is used to obtain the deformation field. Simplified explicit expressions for the elements of the 6 × 6 propagator matrix for the poroelastic medium are obtained. As an example of the possible applications of the analytical formulation developed, formal solution is given for normal strip loading, normal line loading and shear line loading.