Multiples suppression by adaptive subtraction methods based on second order statistic can obtain a perfect result in case the primaries and multiples are orthogonal, but the assumption is always not valid. In this paper, we combine independent component analysis (ICA) with inverse scattering series (ISS) method to suppress surface-related multiples. Multiples predicted by ISS and raw seismic records make up several observation components, then multiples are separated from the seismic records by ICA. Tests on synthetic data generated by convolution and finite-difference modeling show the advantages of ICA in multiples suppression even mismatching in amplitude between the predicted multiples and raw seismic data. But there is uncertain in multiples separation by ICA in case of low signal-to-noise ratio.