This paper investigated the problem of confined flow under dams and water retaining structures using stochastic modelling. The approach advocated in the study combined a finite elements method based on the equation governing the dynamics of incompressible fluid flow through aporous medium with a random field generator that generates random hydraulic conductivity based on log normal probability distribution. The resulting model was then used to analyse confined flow under a hydraulic structure. Cases for a structure provided with cutoff wall and when the wall did not exist were both tested. Various statistical parameters that reflected different degrees of heterogeneity were examined and the changes in the mean seepage flow, the mean uplift force and the mean exit gradient observed under the structure were analysed. Results reveal that under heterogeneous conditions, the reduction made by the sheet pile in the uplift force and exit hydraulic gradient may be underestimated when deterministic solutions are used.