A stochastic model reduction method for nonlinear unconfined flow with multiple random input fields