We review the mathematical models that have been proposed for bioremediation processes, confining ourselves to the case of homogeneous soils with only one species of bacteria and one species of pollutant. We focus mainly on the case of unsaturated soil, trying to clarify some aspects, since most of existing studies refer to saturated soils. The unsaturated case is of great interest, because it is precisely the situation of the vadose zone. Therefore we will try to present a sufficiently complete model for the vadose zone, including some original feature. Finally, a dimensional analysis will lead to a simplified formulation.