Both a theoretical and an empirical model were developed for predicting the formation of soluble microbial products (SMP) during drinking water biofiltration. Four pilot-scale biofilters with ceramsite as the medium were fed with different acetate loadings for the determination of SMP formation. Using numerically simulated and measured parameters, the theoretical model was developed according to the substrate and biomass balance. The results of this model matched the measured data better for higher SMP formation but did not fit well when SMP formation was lower. In order to better simulate the reality and overcome the difficulties of measuring the kinetic parameters, a simpler empirical model was also developed. In this model, SMP formation was expressed as a function of fed organic loadings and the depth of the medium, and a much better fit was obtained.