It is deemed impossible to calculate the vacuum averages of gauge-invariant quantities in QCD in a closed form, e.g. using the path integrals. But standard perturbation theory techniques don't work at distances, where the running coupling constant reaches 1.
The stochastic vacuum model is based on the approximation of nonperturbative QCD as a Gaussian process. It allows to calculate Wilson loops.