Seepage face boundary condition

Dear Users,
I am wondering is there is an option to define the seepage face boundary condition as described in this paper or in the following figure.
seepage_surface
I want to simulate the drawdown of the reservoir. Does anybody have tricks to simulate such a scenario?