We solve the problem of constructing the mathematical model of a vibroseis signal distorted by harmonics. Basis function decomposition is used for the approximation. The basis decomposition coefficients are frequency-dependent, i. e. filters are applied instead of multipliers. The resulting model will be the starting point for further research aimed at developing processing algorithms for the suppression of harmonic distortion, or for utilizing the harmonics to extend the frequency range of the signal.