PREDICTION OF CO2 FREEZE-OUT CONDITIONS AND CO2 HYDRATE FORMATION CONDITIONS WITH A SINGLE EQUATION OF STATE

A recent study by Eggeman and Chaffin (2005 Eggeman , T. , and Chafin , S. ( 2005 ). Beware of pitfalls of CO₂ freezing predictions , Chem. Eng. Prog. , 101 ( 3 ), 39 – 44 . [Google Scholar]), which showed large discrepancies in CO₂ freeze-out conditions as predicted by several commercial simulators, prompted a reexamination of using the TBS equation of state for phase equilibrium calculations involving solids. Salim and Trebble (1994 Salim , P. , and Trebble , M. A. ( 1994 ). Modelling of solid phases in thermodynamic calculations via translation of a cubic equation of state at the triple point , Fluid Phase Equilib. , 93 , 75 – 99 .[Crossref] , [Google Scholar]) had previously presented a methodology for extending the Trebble-Bishnoi-Salim (TBS) equation of state (Salim, 1990 Salim , P. ( 1990 ). A modified Trebble-Bishnoi equation of state, M.Sc. thesis , University of Calgary . [Google Scholar]) to calculations involving a solid phase. In this study, the CO₂ freeze-out conditions in CO₂/CH₄ and CO₂/C₂H₆ mixtures are calculated from the TBS equation of state, and it is shown that they provide a better data fit than the traditional Poynting correction method. Furthermore, since the use of an equation of state in SLE/SVE calculations does not require the explicit assumption of a pure solid phase, the model was assessed for its ability to correlate CO₂ gas hydrate equilibrium conditions. Gas hydrates were simply treated as an impure solid phase, and it was seen that the predictions of gas hydrate equilibrium were in very good agreement with the experimental data. Computationally, the use of the TBS equation of state has the advantage, over the model of Yokozeki (2005 Yokozeki , A. ( 2005 ). Methane gas hydrates viewed through unified solid-liquid-vapour equation of state , Int. J. Thermophys. , 26 ( 3 ), 743 – 765 . [Google Scholar]), that it does not require a modifying factor (c_b) in the repulsive term to handle the presence of hydrates; they are instead handled using a unique binary interaction parameter for the hydrate phase.