The efficiency of naphthalene degradation by Pseudomonas putida G7 in soil was assessed using a mathematical model. The number of microorganisms and the concentration of naphthalene in soil samples were monitored. The feasibility of a spectrofluorometric method for naphthalene assay in soil samples was compared with high pressure liquid chromatography. A proposed mathematical model described the growth of the naphthalene-degrading strains and the consumption of substrates (naphthalene, naphthalene degradation intermediates and soil organic substances) in soil. To describe the growth kinetics of microorganisms having high affinity to substrates with low solubility, two differential equations with substrate exponent 2/3 were proposed. These equations were used to describe utilization of soil organic matter. The model parameters characterize the growth rates for different substrates and respective yield coefficients, specific bacterial death and adaptation rates, and also the rates of PAHs degradation and evaporation. These characteristics can be used in choosing the bacterial strains for biopreparations and efficient clean-up biotechnology of polluted soils.