The vibroseis sources always produce harmonics, which are considered here as noise to be eliminated. An algorithm for removing such noise has been developed based on the previously obtained mathematical model of a vibroseis wavelet distorted by the harmonics. It is shown that the ideal inversion operator is recursive, i.e. is a filter with infinite impulse response. However, due to the peculiarity of the problem, this operator is successfully approximated by a filter with a short impulse response, allowing linearization of the optimization scheme. An objective being the energy of the noise attenuation result is formed. It allows application of the geometrical divergence correction despite the fact that it introduces distortions in the signal shape and the harmonics.