The S wave velocity structure within the depth of 0-100 km beneath digital seismic stations of Yunnan Province was obtained from teleseismic receiver function modeling. The results show that the crustal thickness changes greatly in Yunnan, reaching about 62 km in Zhongdian and Lijiang, and decreasing to 32-34 km in Jinghong, Simao and Changyuan in south. The thick crust extends from northwest to southeast, decreases in thickness and range, and is about 42 km thick around Tonghai. Its shape and range consist with Chuandian diamond block bounded by Xiaojiang and Yuanjiang faults. In the eastern and southern region, the crust is relatively thin, and the velocity contrast across Moho is obvious. In the areas with great variation of crustal thickness or thick crust area, the Moho is characterized by transition zone with high velocity gradient. In Yunnan region, S wave velocity structures show strong horizontal heterogeneity. Above the depth of 10 km, the S wave velocity in north area is obviously lower than that in south area, while within the depth of 10-20 km, the S wave velocity in north is higher than that in south. The velocity interface within crust is discontinuous, and the depth and range of lower velocity zones change with different seismic stations. There is no obvious crustal low velocity zone beneath nearly half of the stations. Influenced by the upper mantle in south area, within the depth of 40-50 km, S wave velocity in south area is higher than that in north, high velocity area extends to north, and the shape of low velocity zone tends to consist with Chuandian diamond block. The upper mantle velocity distribution within the depth of 70-80 km seems to be correlated to the distribution of strong earthquakes.