The aim of this paper is to present some numerical simulations for different values of parameters and to compare their influence for a reaction-diffusion model which describes chemotaxis phenomena in a bioremediation process of a polluted environment (soil or water). We set a model of Keller-Segel type, and we consider for this model a small influence of the degradation rate of pollutant ( ) and its diffusion coefficient ( ). For the numerical simulations, we used the package COMSOL Multiphysics v3.5a, working with the finite element method.