The recently developed discrete Boltzmann method (DBM), which is based on a uniform linear evolution equation and has high parallel efficiency, is employed to investigate the dynamic nonequilibrium process of Kelvin-Helmholtz instability (KHI). It is found that, the relaxation time always strengthens the global nonequilibrium (GNE), entropy of mixing, and free enthalpy of mixing. Specifically, as a combined effect of physical gradients and nonequilibrium area, the GNE intensity first increases but decreases during the whole life-cycle of KHI. The growth rate of entropy of mixing shows firstly reducing, then increasing, and finally decreasing trends during the KHI process. While the free enthalpy of mixing is opposite to the entropy of mixing. Detailed explanations are as below. (i) Initially, binary diffusion smooths quickly the sharp gradient in the mole fraction, which results in a steeply decreasing mixing speed. (ii) Afterwards, the mixing process is significantly promoted by the increasing length of material interface in the evolution of the KHI. (iii) As physical gradients are smoothed due to the binary diffusion and dissipation, the mixing speed reduces and approaches zero in the final stage. Moreover, with the increasing Atwood number, the global strength of viscous stresses on the heavy (light) medium reduces (increases), because the heavy (light) medium has a relatively small (large) velocity change. Furthermore, for a larger Atwood number, the peaks of nonequilibrium manifestations emerge earlier, the entropy of mixing and free enthalpy of mixing change faster, because the KHI initiates a higher growth rate.