Denoising is an important preprocessing step to further analyze a hyperspectral image (HSI). The common denoising methods, such as 2-D (two dimensions) filter, actually rearrange data from 3-D to 2-D and ignore the spectral relationships among different bands of the image. Tensor decomposition method has been adapted to denoise HSIs, for instance Tucker3 (three-mode factor analysis) model. However, this model has problems in uniqueness of decomposition and in estimation of multiple ranks. In this paper, to overcome these problems, we exploit a powerful multilinear algebra model, named parallel factor analysis (PARAFAC), and the number of estimated rank is reduced to one. Assumed that HSI is disturbed by white Gaussian noise, the optimal rank of PARAFAC is estimated according to that the covariance matrix of the n-mode unfolding matrix of the removed noise should be approach to a scalar matrix. Then, the denoising results by PARAFAC decomposition are presented and compared with those obtained by Tucker3 model and 2-D filters. To further verify the denoising performance of PARAFAC decomposition, Cramer-Rao lower bound (CRLB) of denoising is deduced theoretically for the first time, and the experiment results show that the PARAFAC model is a preferable denoising method since the variance of the HSI denoised by it is closer to the CRLB than by other considered methods.Index Terms-Classification, Cramer-Rao lower bound (CRLB), denoising, hyperspectral image (HSI), parallel factor analysis (PARAFAC), Tucker3.