Ensemble modeling of the two-dimensional stochastic confined groundwater flow through the evolution of the hydraulic head's probability density function