This paper derives the arbitrary order globally hyperbolic moment system for a non-linear kinetic description of the Vicsek swarming model by using the operator projection. It is built on our careful study of a family of the complicate Grad type orthogonal functions depending on a parameter (angle of macroscopic velocity). We calculate their derivatives with respect to the independent variable, and projection of those derivatives, the product of velocity and basis, and collision term. The moment system is also proved to be hyperbolic, rotational invariant, and mass-conservative. The relationship between Grad type expansions in different parameter is also established. A semi-implicit numerical scheme is presented to solve a Cauchy problem of our hyperbolic moment system in order to verify the convergence behavior of the moment method. It is also compared to the spectral method for the kinetic equation. The results show that the solutions of our hyperbolic moment system converge to the solutions of the kinetic equation for the Vicsek model as the order of the moment system increases, and the moment method can capture key features such as vortex formation and traveling waves.the notation Id is the identity operator, Ω(t, x) is the mean velocity, and J (t, x) denotes the mean flux at x and is defined by(2.3)Here K(y − x) is the characteristic function of the ball B(0, R) = {x : |x| ≤ R}, i.e. K(x) = 1 |x|