Modelling the phreatic surface in regional flood defences