In this paper we perform a global analysis of the dynamics of the Chen systeṁwhere (x, y, z) ∈ R 3 and (a, b, c) ∈ R 3 . We give the complete description of its dynamics on the sphere at infinity. For six sets of the parameter values the system has invariant algebraic surfaces. In these cases we provide the global phase portrait of the Chen system and give a complete description of the α-and ω-limit sets of its orbits in the Poincaré ball, including its boundary S 2 , i.e. in the compactification of R 3 with the sphere S 2 of the infinity. Moreover, combining the analytical results obtained with an accurate numerical analysis, we prove the existence of a family with infinitely many heteroclinic orbits contained on invariant cylinders when the Chen system has a line of singularities and a first integral, which indicates the complicated dynamical behavior of the Chen system solutions even in the absence of chaotic dynamics.