A device, medium and method for deblending seismic data associated with a subsurface of the earth. The method includes receiving an input dataset generated by first and second sources S1 and S2 that are operating as simultaneous sources; arranging the input dataset based on the firing times of source S1; applying with a computing system an annihilation filter to the arranged input dataset to estimate cross-talk noise; convolving the cross-talk noise estimate with an operator to form a signal estimate using the firing times of S1 and S2; and generating an image of the subsurface based on the signal estimate.