A method and apparatus for noise attenuation. The method includes receiving seismic data associated with at least two vintages (di, dj) collected for a same subsurface, wherein the first and second vintages (di, dj) are taken at different times; calculating a set of filters (fi, fj) that minimizes an energy function (E), wherein the energy function (E) includes a term representing a 4D difference between the first and second vintages (di, dj); calculating primaries (pi, pj) corresponding to the first and second vintages (di, dj) based on the set of (fi, fj); and calculating a 4D difference (Δij) based on the primaries (pi, pj). The 4D difference (Δij) is minimized.