The 3D conical Radon transform for seismic signal processing

Abstract

Recently, the attenuation of highly aliased broadband surface waves has received significant attention as their high-frequency and low-speed components render f- k x- k y filters ineffective. As a popular tool widely used in seismic signal processing, the 3D linear Radon ([Formula: see text]) transform does not match the surface wave events in the time domain. This leads to the dispersal of energy over ellipses on the zero-intercept slice, whereas the reflected events over the ellipses are reflected on each intercept slice of the 3D linear Radon transform (LRT) domain, decreasing the resolution in slowness. We have introduced a new type of Radon transform defined on circular cones called 3D conical Radon transform (CRT) for seismic signal processing. Unlike LRT, the CRT maps seismic data to surface integrals on circular cones in the 3D seismic records. Consequently, surface wave events are focused as points on the zero-intercept slice, whereas the reflected events are points on each intercept slice of the CRT domain, which significantly improves the resolution in slowness compared with the LRT. We determine the performance of the CRT for seismic signal processing for random noise attenuation, primary and multiple separations, and surface wave attenuation. Based on a comparison with the LRT, synthetic data and field data sets validate the effectiveness of the CRT method at improving the focusing properties.