This paper deals with the analytical description of single-phase flow caused by abstraction wells and governed by Darcy＇s law. Since we are mainly interested in the velocity field upon which an analysis of transport phenomena can be based, we may assume that the flow is quasi steady. A well may be composed of a number of point sinks, therefore the main problem is to determine the flow field caused by one point sink. For a perfectly layered subsurface the solution can be found by an approximation of the Hankel Transform using a Fast Fourier Transform. The final result is an algorithm that can be used to determine velocities and flow paths on any spatial scale, as will be shown by some examples.