Method for processing seismic data. In one implementation, the method includes converting a common midpoint (CMP) gather of seismograms into one or more interferogram common midpoint (ICMP) gathers, generating a semblance spectrum for each ICMP gather, stacking the semblance spectrum from each ICMP gather to generate a combined semblance spectrum and deriving a normal moveout (NMO) velocity profile from the combined semblance spectrum.