We modeled the firing rate of populations of olfactory receptor neurons (ORNs) responding to an odorant at different concentrations. Two cases were considered: a population of ORNs that all express the same olfactory receptor (OR), and a population that expresses many different ORs. To take into account ORN variability, we replaced single parameter values in a biophysical ORN model with values drawn from statistical distributions, chosen to correspond to experimental data. For ORNs expressing the same OR, we found that the distributions of firing frequencies are Gaussian at all concentrations, with larger mean and standard deviation at higher concentrations. For a population expressing different ORs, the distribution of firing frequencies can be described as the superposition of a Gaussian distribution and a lognormal distribution. Distributions of maximum value and dynamic range of spiking frequencies in the simulated ORN population were similar to experimental results.