Orbits of pairs in abelian groups

We compute the number of orbits of pairs in a finitely generated torsion module (more generally, a module of bounded order) over a discrete valuation ring. The answer is found to be a polynomial in the cardinality of the residue field whose coefficients are integers which depend only on the elementary divisors of the module, and not on the ring in question. The coefficients of these polynomials are conjectured to be non-negative integers.