In time-lapse seismic applications, the signal produced by changes in the properties of subsurface rocks is generally obscured by noise associated with imperfect repeatability between surveys. A particularly important obstacle in the formation of time-lapse difference images is variation in the effective source wavelet between baseline and monitoring data sets. However, the partially separable influence of the wavelet within Green’s function model of seismic data permits two frequency-domain matching filters to be designed, which act to reduce source wavelet nonrepeatability. One is based on the spectral ratio of the baseline and monitoring wavelets and can be applied when prior estimates of the wavelets are available; the other is the average spectral ratio of the baseline and monitoring traces and can be applied when prior estimates are unavailable. After balancing the data sets with either of these filters, we further prepare for the imaging step with time-shift corrections, using a published fast local crosscorrelations algorithm, preparing the difference data for use in an imaging algorithm. Reverse time migration is engaged for the imaging task, but we observe that residual repeatability errors tend to be magnified at this stage when source-normalized crosscorrelation imaging conditions are used. Testing indicates that replacing this with a Poynting vector imaging condition strongly suppresses remaining errors in a robust manner. This includes stability within simulated data environments to noise-free data and data with random noise, up to signal-to-noise ratios of roughly one. Furthermore, our method illustrates better performance when compared with the conventional least-squares matching filter and common-depth-point-domain warping. At present, there is no common workflow for seismic imaging directly using time-lapse shot gathers. Our research contribution lies not only in the two matching filters but also in a novel workflow for time-lapse seismic imaging.