The transient flow behavior of ground-water in aquifer-aquitard system with viscoelastic properties is studied. On the basis of previous works (Hantush, Neuman, Brutsaert, Corapcioglu), the new partial differential-integral equations are derived. The well-known equations (Hantush, Brutsaert) are the special cases of the new equations. The new equations describe the flow of a slightly compressible groundwater in layers with viscoelastic properties.Analytical solutions of the partial differential-integral equations are obtained by using the method of Laplace transform. The viscoelastic properties enhance the heterogeneities of elastic aquifer system which have delay and feed qualities. The agreements between the numericl inversion results of Laplace transform and the analytical solutions are good. The formulae predict the transient flow behavior of groundwater in this heterogeneous layers.