In this contribution we present an exploratory study of several novel methods for numerical stochastic perturbation theory. For the investigation we consider observables defined through the gradient flow in the simple {phi}^4 theory.