Non-linear analysis of nearly saturated porous media: theoretical and numerical formulation