1. A method for use in seismic data processing, comprising: improving the continuity of the primary events in an original seismic dataset to create a continuous seismic dataset; predicting at least the first order multiples of the original seismic dataset from the continuous seismic dataset; deterministically matching the predicted multiples and one of the original seismic dataset or the continuous seismic dataset; and taking the difference between the matched datasets.