Abstract: In some embodiments, a seismic processing method comprises assembling a specularity gather by determining a specularity value at each of a plurality of subsurface locations, and summing trace amplitudes into a plurality of bins, each bin characterized by a range of specularity values. The specularity value at a subsurface location is computed according to an angle between a normal to a local reflector and a direction of a total (source+receiver) traveltime gradient. For example, the specularity may be proportional to (e.g. equal to) a magnitude of the cosine of the angle. A diffraction image may be generated by summing specularity gather data over specularity, with specular event amplitudes attenuated relative to diffractive event amplitudes.

Abstract: In some embodiments, a time migration diffraction imaging method includes computing a pseudo-depth characterizing a subsurface seismic event by scaling a vertical traveltime using a scaling velocity. A specularity value for a subsurface seismic event is determined according to the pseudo-depth, and a contribution weight for a corresponding seismic trace amplitude is determined according to the specularity value. The specularity value may be determined according to an angle between a traveltime gradient and a normal to a local reflector surface. A diffraction image is generated according to a weighted sum of seismic trace amplitudes. The weighted sum attenuates the contribution of specular events relative to diffraction events.