A unified colored-noise approximation (UCNA) method has been widely used to solve the simple one-dimensional problem, and this paper attempts to extend this method to multi-dimensional systems. Firstly, a Fokker-Planck equation is obtained by UCNA method based on the Langevin equation for particle motion, then a two-order moment trajectory model is established on this basis. The new model can be successfully used to predict the backward-facing step two-phase flow, and the simulation results agree well with the measurements. This study shows that the UCNA method is still effective in dealing with multi-dimensional two-phase turbulent systems.