The effect of near surface interfaces with irregular topography on the target reflections from deeper layers is a well known problem in digital seismic data processing. Due to strong contrast of velocities imposed by such interfaces, theses reflections lose their hyperbolic characteristic. As a consequence, the CMP stack procedure, which main assumption is the hyperbolic shaped reflections, fails in its attempt to improve the signal to noise ratio of the seismic section. To overcome this difficulty, a preprocessing technique which recovers fully or partially the hyperbolic pattern of these reflections is necessary. In this context, this work is divided in two parts. First, three known different algorithms used for the solution of this problem, with distinct levels of precision and computational cost, are discussed. Following that, a new algorithm which integrates accurate wave field extrapolation mechanisms with ray tracing techniques, is proposed. The performance of the new method is demonstrated on synthetic common mid point gathers and common shot gathers. Finally, these results are then compared with those obtained from one of the algorithms previously described, in order to assess the limitations and advantages of the proposed technique.