Mapping prestack depth-migrated coherent signal and noise events back to the original time gathers using Fermat's principle

Because of its computational efficiency, prestack Kirchhoff depth migration is currently the method of choice in both 2-D and 3-D imaging of seismic d ata. The most algorithmically complex component of the Kirchhoff family of algorithms is the calculation and manipulation of accurate traveltime table s for each source and receiver point. Once calculated, we sum the seismic e nergy over all possible ray paths, allowing us to accurately image both spe cular and nonspecular scattered energy. Any seismic events that fall within the velocity passband, including reflected and diffracted signal, mode con versions, multiples, head waves, and aliases of surface waves, are imaged i n depth. The transformation of time gathers to depth gathers can be quite complicate d and nonintuitive to all but the seasoned imaging expert. In particular, e asily recognized head-wave events on common-shot gathers are often difficul t to differentiate from undermigrated coherent reflections on common-reflec tion-point depth gathers. In contrast, subsalt multiples that have propagat ed along complex ray paths are often easily recognized on common-offset dep th gathers but are indistinguishable from the distorted primaries on the in put common-shot or common-midpoint time gathers. In a related area, seismic reflection traveltime tomography is currently th e workhorse for 2-D and an active area of research and development for 3-D migration-driven velocity analysis. The objective function for this "veloci ty inversion" problem is to either minimize the temporal difference between picked and modeled time picks, or to maximize the similarity between, or f latness of, common-reflection-point depth picks. Once picked and associated with the correct reflector, time picks never need to be modified during th e velocity-model updating steps that ultimately lead to a feasible solution . In practice, such time picks are nearly impossible to make in those struc turally complex areas that justify the use of prestack depth migration. Ins tead, we almost always use the second objective function and pick reflector events in depth, where we can use our geologic insight to differentiate be tween signal and noise and where the difficulty of associating a picked eve nt with the velocity/depth model horizon completely disappears. The major d rawback of picking in depth is that these events need to be repicked each t ime any part of the overlying velocity/depth model has been updated. We show that by applying Fermat's principle, and by reusing the same travel time tables used in seismic prestack Kirchhoff depth imaging, we can map in terpreted events on the depth gathers to corresponding interpreted events o n the original time gathers. This technique, first introduced by J. van Tri er in 1990, is considerably more stable and, because we reuse the already c omputed migration traveltime tables, more economic than two-point ray-trace methods. In our first application of coherent noise suppression, we show how we can relate imaging artifacts seen on the depth image to their causative coheren t noise on the original time gathers. Once identified, these noise events c an be safely suppressed using conventional filtering techniques. In our sec ond application of reflection tomography, we show how we can pick partially focused reflectors in depth, and map them back to time, undoing the effect of the incorrect velocity/depth model used in prestack Kirchhoff depth mig ration such that the events never need to be repicked during subsequent vel ocity model updates.