1. A method for predicting a plurality of surface multiples for a plurality of traces in a record of seismic data, comprising: providing a plurality of target traces at a nominal offset and a nominal azimuth; selecting a plurality of pairs of input traces, wherein the midpoints of the input traces are separated by half the nominal offset and the azimuth of a line connecting the midpoints of the input traces is equal to the nominal azimuth; convolving, using a computer, the selected pairs of input traces to generate a plurality of convolutions; and applying a three dimensional operator to the convolutions, wherein the three dimensional operator is a three dimensional demigration operator.