A Model Reduction Method for Elliptic Pdes with Random Input Using the Heterogeneous Stochastic Fem Framework