Based on the features of diurnal variation of ice surface albedo in the Lake Ulansuhai, the relationship between solar elevation angle and geographic coordinates and Julian days, together with normalized time, to express the diurnal variations of the ice albedo. Linear combinations of four probability density distribution functions with exponential forms, including Laplace, Gauss, Gumbel, and Cauchy equations, were used to simulate the diurnal variations in observed ice albedo from sunrise to sunset with solar elevation angle more than 5°. The results reveals that the Laplace combined statistical model is the best fit to the observations. It can not only clearly show the double-peak distribution in the diurnal variation in ice surface albedo when the elevation angle is greater than 5°, but also express the U-type distribution between the two peaks as the elevation angle is more than 15°. This model has advantages of simple form in expression and clear physical meaning. The length parameter is close to half day, and the peak location is associated with the time of sunrise, and the asymmetry of the two peaks can be also shown. It provides a solid foundation for the parameterization schemes on the diurnal variations in lake ice surface albedo in different regions.