A Chebyshev polynomial-based Galerkin method for the discretization of spatially varying random properties