A three-dimensional CFD-DEM model is proposed to investigate the aeolian sand movement. The results show that the mean particle horizontal velocity can be expressed by a power function of heights. The probability distribution of the impact and lift-off velocities of particles can be described by a log-normal function, and that of the impact and lift-off angles can be expressed by an exponential function. The probability distribution of particle horizontal velocity at different heights can be described as a lognormal function, while the probability distribution of longitudinal and vertical velocity can be described as a normal function. The comparison with previous two-dimensional calculations shows that the variations of mean particle horizontal velocity along the heights in two-dimensional and three-dimensional models are similar. However, the mean particle density of the two-dimensional model is larger than that in reality, which will result in the overestimation of sand transportation rate in the two-dimensional calculation. The study also shows that the predicted probability distributions of particle velocities are in good agreement with the experimental results.