Methods for processing seismic data are described. The method includes: obtaining seismic data; solving a series of partial differential wave equations, wherein a first partial differential wave equation describes propagation of a seismic wave going from a first reflector to a second reflector, wherein a second partial differential wave equation describes propagation of a seismic wave going from a second reflector to a third reflector, and wherein a third partial differential wave equation describes propagation of a seismic wave going from a third reflector to a seismic receiver, wherein outputting predicted internal multiples for further imaging or attenuation.