Statistics for Multi-scale perturbation theory. Part I. Methodology and leading-order bispectrum corrections in the matter-dominated era