A method for predicting a plurality of surface multiples for a plurality of target traces in a record of seismic data acquired in a survey area. The method includes selecting a target trace and identifying two or more desired traces for multiple prediction based on the target trace. After identifying the desired traces, the method identifies one or more recorded traces for each desired trace. Each identified recorded trace is described as being substantially close to one of the desired traces. The method then includes correcting the identified recorded traces for one or more geometry-related effects associated with the survey area and convolving the corrected recorded traces to generate a plurality of convolutions. After convolving the corrected recorded traces, the method then stacks the convolutions.