Variance reduction in calculations of multiple integrals by Monte Carlo method is studied. Invariance reasons are applied to the method suggested by Ermakov and Zolotukhin for the variance reduction. Random invariant cubature formulas, which appear to combine high efficiency with simplification of the modelling problem, are studied and constructed.