Three equations of state with a group contribution model for binary interaction parameters were employed to calculate the vapor-liquid equilibria of synthetic and real natural gas mixtures with heavy fractions. In order to estimate the binary interaction parameters, critical temperatures, critical pressures and acentric factors of binary constituents of the mixture are required. The binary interaction parameter model also accounts for temperature. To perform phase equilibrium calculations, the heavy fractions were first discretized into 12 Single Carbon Numbers (SCN) using generalized molecular weights. Then, using the generalized molecular weights and specific gravities, the SCN were characterized. Afterwards, phase equilibrium calculations were performed employing a set of (nc + 1) equations where nc stands for the number of known components plus 12 SCN. The equations were solved iteratively using Newton's method. Predictions indicate that the use of binary interaction parameters for highly sour natural gas mixtures is quite important and must not be avoided. For sweet natural gas mixtures, the use of binary interaction parameters is less remarkable, however.