Integration of moment equations in a reduced-order modeling strategy for Monte Carlo simulations of groundwater flow