We consider scalar hyperbolic conservation laws in several $(d \ge 1)$ spatial dimensions with stochastic initial data. We prove existence and uniqueness of a random-entropy solution and show existence of statistical moments of any order k of this random entropy solution. We present a class of numerical schemes of multi-level Monte Carlo Finite Volume (MLMC-FVM) type for the approximation of random entropy solutions as well as of their ...