A method for use in seismic exploration comprises: accessing a set of seismic data representative of a subterranean geological formation and a subsurface attribute model of the subterranean geological formation; performing a wavefield extrapolation on the seismic data in the subsurface attribute model; applying the time-shift extended imaging condition to the extrapolated wavefields; forming shot-indexed, time shift gathers for each image pixel of the subsurface attribute model from the conditioned extrapolated wavefields; adaptively focusing the gathers; and stacking the adaptively focused gathers; and imaging the subterranean geological formation from the stacked, adaptively focused gathers. The method may, in some aspects, be realized by a computing apparatus programmed to perform the method or as a set of instructions encoded on a non-transitory program storage medium that, when executed by a computing apparatus, perform the method.