The paper is devoted to a new computational method for problems of transport in highly non-uniform media. In particular, the method is applied to the problem of anomalous contaminant transport in a field with a randomly distributed permeability, which was modelled as a stochastic advection process governed by a stochastic advection model. The stochastic advection model is used to generate different realisations of micro-dispersion parameters required for direct numerical simulations. The new numerical method combines the merits of finite-volume and finite-difference approaches and is demonstrated to be efficient and robust in several benchmark advection tests. For the stochastic advection problem considered the results of the new computational method are in a good agreement with analytical predictions available for different stochastic advection regimes.