Double compression of images occurs when one compresses twice, possibly with different quality factors, a digital image. Estimation of the first compression parameter of such a double compression is of a crucial interest for image forensics since it may help revealing, for instance, the software or the source camera. This paper proposes an accurate method for estimating the primary quantization steps in double-compressed JPEG images. This original methodology is based on an accurate statistical model of discrete cosine transform (DCT) coefficients that has been proposed in our previous works. We also present a thorough analysis of the double compression properties, taking into account carefully the effect of round-off noise. This analysis is used to derive an accurate range of possible value for quantization of primary DCT coefficients with respect to the secondary quantization step. Using both the statistical model of quantized DCT coefficients and the range of possible values of first quantization step, a model of the twice quantized DCT coefficients is established. Eventually, it is proposed to estimate the primary quantization value by finding, among a set of possible candidates, the one that best match the proposed statistical model in terms of minimal symmetrized Kullback-Leibler (KL) divergence. The numerical experiments on large databases of real images and comparisons with state-of-the-art approaches emphasize the relevance of the proposed method.