Method for correcting OBC or deep-towed seismic streamer data for surface-related multiple reflections. The measured pressure data, preferably after conditioning (71), are simulated using a forward model that includes a water propagation operator between source locations and receiver locations and a term representing primary impulse responses (72). Other terms include direct arrivals and source wavelets. Iterative optimization of an objective function is used to minimize the difference between measured and simulated data, updating the primary impulse response term and optionally the source wavelets term each iteration cycle (73). The converged primary impulses (74) are used to construct simulated multiples and direct arrivals (75), which can be subtracted from the measured data. Optionally the measured data might be blended during the forward simulation (72), to save computational costs in the forward simulation (72) and in the inversion (73).