Based upon the nonlinear model of random sea waves, the statistical distribution of wave surface elevation exact to the fourth order is derived as the truncated Gram-Charlier series, by calculating directly each order moment. The phenomenon found by Huang et al. that the agreement between observed data and investigated series deteriorates much more when the series is kept to λs is explained, The effect of the appmximation order on the truncation of series and the determination of coefficients is investigated. For the mth order approximation, the derived series is truncated at H3m-3 with the absence of H3m-4, and the coefficients of H3m-3 and H3m-6 are connected by a simple algebraic relation.