Based on ECMWF ERA-Interim reanalysis data, 8 factors (T2 m, T1000, T975, T950, T925, T850, H700-850 and H850-1000) for identifying precipitation phases were obtained through analyzing the temperature and geopotential thickness of precipitation phases (rain, snow, sleet) in winter half year from 2008 to 2017 in Shandong Province, and the threshold indicators of the 8 factors were provided. The discriminant equation for precipitation phase identification was established and the deep learning DNN model was trained using the 8 factors and their threshold values, and the forecast accuracy of rain, snow and sleet increased by 1.9%, 0.2% and 21.6% using DNN method through randomization test, respectively. The inspection using ECMWF fine grid model products indicated that among a total of 106 stations of rain, snow and sleet, the discriminant equation and DNN method carried out wrong identifications for 29 and 14 stations, respectively. The results show that the DNN method performed better than the discriminant equation, and in particular, it significantly improved the identification ability of sleet.