We present a finite difference-(Multi-level) Monte Carlo algorithm to efficiently compute statistical solutions of the two dimensional Navier-Stokes equations, with periodic bound- ary conditions and for arbitrarily high Reynolds number. We propose a reformulation of statistical solutions in the vorticity-stream function form. The vorticity-stream function for- mulation is discretized with a finite difference scheme. We obtain a convergence ...