To obtain accurate migration velocities, we must estimate the velocity at migrated depth points. Wavefront focusing analysis with downward continuation yields the rms velocity at migrated depth points; however, the large amount of computation required for downward continuation limits use of this approach for routine processing. The purpose of this paper is to present an implementation of the Kirchhoff integral which makes the wavefront focusing analysis practical for time-migration velocity analysis. Downward continuation focuses the wavefront to the zero offset at the depth controlled by the velocity used for the continuation. The migration velocity is then determined from the depth where the focused wavefront attains the maximum amplitude. The flexibility of the Kirchhoff integral allows us to compute only the zero-offset trace at each depth point and lets us avoid most of the computation for the downward continuation of unstacked data. Furthermore, since the velocity is obtained from the location where the focused wavefront shows the maximum amplitude, prestack time migration with the velocity from this technique produces the maximum amplitude for the subsurface reflector.