The accurate prediction of the behaviour of geostructures is based on the strong coupling between the pore uid and the solid skeleton. If the relativ e acceleration of the uid phase to the skeleton is neglected, the equations describing the problem can be written in terms of skeleton displacements (or velocities) and pore pressures. This mixed problem is similar to others found in solid and uid dynamics. In the limit case of zero perme-abilit y and incompressibility of the uid phase, the restrictions on the shape functions used to approximate displacements and pressures imposed by Babuska-Brezzi conditions or the Zienkiewicz-Taylor patch test hold. As a consequence, it is not possible to use directly elements with the same order of interpolationfor the eld v ariables. This paper proposes tw o alternative methods allowing us to circumvent the BB restrictions in the incom- pressibilit y limit, making it possible to use elements with the same order of interpolation. The rst consists on in troducing the divergence of the momentum equation of the mixture as an stabilization term, the second is a generalization of the tw o step-projection method introduced by Chorin for uid dynamics problems.