CFD simulation with DES-SA turbulence model are conducted to investigate the wake velocity characteristics. The velocity distribution of axial, tangential and radial velocity component is analysed by comparing numerical results with previous experimental and theoretical results. The axial velocity component is continuously recovered along the axial direction. Based on the position of minimum axial velocity, turbine wake can be divided into two zones: zone of flow establishment and zone of established flow. In the zone of flow establishment, two-dipped velocity valleys of the axial velocity distribution appear along the radial direction. These two valleys are combined to one valley at the rotation axis in the zone of established flow. The tangential velocity component is the second largest and the radial velocity component accounts for only a small proportion. Both of tangential and radial velocity decreases along the axial direction. Two velocity peaks of both tangential and radial velocity appear along the radial direction. Several empirical equations are finally proposed to describe the velocity distribution of turbine wake.