Numerical methods for hyperbolic SPDEs: a Wiener chaos approach