In this work, the bacterial consortium Bulab 5738 was used to simultaneously remove pyrene, phenanthrene and catechol from aqueous solutions. The bacterial population growth was estimated by means of optical density measurements while HPLC was used to quantify the pollutant concentration in the solution. The obtained data were used to model the systems, in term of biomass population and substrate concentration. The effects of pollutant concentration values were analysed using the outputs of a full factorial experimental design.