Leaky aquifers constitute complicated hydrological structures, whose inclusion in numerical models of hydrological systems is difficult, because of their three-dimensional nature. Methods for treating such systems can be classified as fully three-dimensional and quasi-three-dimensional. The latter have clear numerical advantages when applicable. In this paper a critical discussion of existing quasi-three-dimensional models is presented.