A reduced-order model for Monte Carlo simulations of stochastic groundwater flow