To solve the vulnerability of singular value denoising to the effective number of singular value orders and influence of Hankel matrix structure, the paper has proposed a multiscale singular value denoising algorithm (M-SVD). Meanwhile, in concerning that when a fault of bearing occurs, the fault feature information contained in vibration signal tends to be weak and complex, the M-SVD method is blended with variational mode decomposition (VMD) algorithm in identification. Firstly, when bearing faults occur, raw vibration acceleration signals are decomposed using the superiority of VMD algorithm for non-stationary signal; in addition, considering that with the larger difference between singular values of signal, feature information will be more obvious, the method takes advantage of sensibility of variance to signal intensity and proposes to constitute Hankel matrix based on variance feature of signal. During the construction of Hankel matrix, the method to determine effective length of variance sequence (optimal scale) has been brought forward by sample entropy which is able to precisely measure the regularity of signal and detect the weak change of signal. Furthermore, component signals are denoised in line with reconstructed Hankel matrix and the spectrum of signals denoised are selected to describe fault characteristic information and judge the failure types. For the sake of verification of the accuracy and effectiveness of proposed method, the paper has analyzed and verified the data of rolling bearing in 3 different failure types and 2 different installation positions of accelerometers. The comparative analysis with classical singular value difference spectrum (SVDS) and singular value median decomposition (SVMD) denoising has further proved that presented VMD combined M-SVD method can effectively restrain noise while retaining fault characteristic information of signal. Therefore, the identification for rolling bearing fault types can be implemented, accurately.