A general mathematical model to solve the advection-dispersion transport equation for multiple solutes was developed, where the dual porosity mobile-immobile mass transfer, the two-site non-equilibrium model and firstorder transformation reactions were included. The two-site model was expressed with an equilibrium sorption term and a kinetic term. One of three kinetic models could be used: the non-linear, the bilinear and the pore diffusion model. The traditional Freundlich or Langmuir isotherms were employed to simulate no-interaction between the solutes, but with the extended Freundlich or extended Langmuir isotherms, a competitive sorption could be simulated. The transport equation was solved with the Moving Concentration Slope method. The mathematical model was tested and further simplified by using real data from soil column experiments, with 1,2-cis-dichloroethene and trichloroethene as model contaminants and silica gel and real soil samples as porous medium.