An alternative way for the derivation of the new KdV-type equation is presented. The equation contains terms depending on the bottom topography (there are six new terms in all, three of which are caused by the unevenness of the bottom). It is obtained in the second order perturbative approach in the weakly nonlinear, dispersive and long wavelength limit. Only treating all these terms in the second order perturbation theory made the derivation of this KdV-type equation possible. The motion of a wave, which starts as a KdV soliton, is studied according to the new equation in several cases by numerical simulations. The quantitative changes of a soliton's velocity and amplitude appear to be directly related to bottom variations. Changes of the soliton's velocity appear to be almost linearly anticorrelated with changes of water depth whereas correlation of variation of soliton's amplitude with changes of water depth looks less linear. When the bottom is flat, the new terms narrow down the family of exact solutions, but at least one single soliton survives. This is also checked by numerics.