Application of the cell potential method to predict phase equilibria of multicomponent gas hydrate systems

We present the application of a mathematical method reported earlier by which the van der Waals-Platteeuw statistical mechanical model with the Lennard-Jones and Devonshire approximation can be posed as an integral equation with the unknown function being the intermolecular potential between the guest molecules and the host molecules. This method allows us to solve for the potential directly for hydrates for which the Langmuir constants are computed, either from experimental data or from ab initio data. Given the assumptions made in the van der Waals-Platteeuw model with the spherical-cell approximation, there are an infinite number of solutions; however, the only solution without cusps is a unique central-well solution in which the potential is at a finite minimum at the center to the cage. From this central-well solution, we have found the potential well depths and volumes of negative energy for 16 single-component hydrate systems: ethane (C2H6), cyclopropane (C3H6), methane (CH4), argon (Ar), and chlorodifluoromethane (R-22) in structure I; and ethane (C2H6), cyclopropane (C3H6), propane (C3H8), isobutane (C4H10), methane (CH4), argon (Ar), trichlorofluoromethane (R-11), dichlorodifluoromethane (R-12), bromotrifluoromethane (R-13B1), chloroform (CHCl3), and 1,1,1,2-tetrafluoroethane (R-134a) in structure II. This method and the calculated cell potentials were validated by predicting existing mixed hydrate phase equilibrium data without any fitting parameters and calculating mixture phase diagrams for methane, ethane, isobutane, and cyclopropane mixtures. Several structural transitions that have been determined experimentally as well as some structural transitions that have not been examined experimentally were also predicted. In the methane-cyclopropane hydrate system, a structural transition from structure I to structure II and back to structure I is predicted to occur outside of the known structure II range for the cyclopropane hydrate. Quintuple (L(w)-sI-sII-L(hc)-V) points have been predicted for the ethane-propane-water (277.3 K, 12.28 bar, and x(eth,waterfree) = 0.676) and ethane-isobutane-water (274.7 K, 7.18 bar, and x(eth,waterfree) = 0.81) systems.