Experimental and thermodynamic modelling of systems containing water and ethylene glycol: Application to flow assurance and gas processing

Accurate knowledge of hydrate phase equilibrium in the presence of inhibitors is crucial to avoid gas hydrate formation problems and to design/optimize production, transportation and processing facilities. In this communication, we report new experimental dissociation data for various systems consisting of methane/water/ethylene glycol and natural gas/water/ethylene glycol. A statistical thermodynamic approach, with the Cubic-Plus-Association equation of state, is employed to model the phase equilibria. The hydrate-forming conditions are modelled by the solid solution theory of van der Waals and Platteeuw. The thermodynamic model was used to predict the hydrate dissociation conditions of methane and natural gases in the presence of distilled water or ethylene glycol aqueous solutions. Predictions of the developed model are validated against independent experimental data and the data generated in this work. A good agreement between predictions and experimental data is observed, supporting the reliability of the developed model.     

Calculating the separation coefficients in argon, krypton and xenon gas mixture separation by gas hydrate crystallization

The separation of gas mixtures based on noble gases by means of gas hydrate crystallization is considered for various schemes of the process. The separation coefficient and the separation factor of the basic mixture components at 273 K and at different pressures and compositions of the gas mixture are calculated. It is shown that xenon can be almost completely separated from a mixture in two stages of gas hydrate crystallization even without the use of booster gases. Directional crystallization allows us to separate xenon more completely and crystallization at constant pressure enables us to obtain a more highly purified product.     

Sorption equilibria of CO2/CH4 mixture on activated carbon in presence of water

The sorption isotherms of CO2 + CH4 mixtures on an activated carbon were collected in the presence of water at a temperature suitable for hydrate formation. The equilibrium composition of both phases was determined. The initial concentration of CO2 in mixtures was set at 33, 38 and 42%, and the total pressure was up to 10 MPa. CO2 hydrates were firstly formed following the increase of total pressure, and CO2 dominates the sorbed phase composition. CO2 concentration in the sorbed phase begins to decrease when the partial pressure of methane allows for the formation of methane hydrates. Competition for hydrate cavities was observed between CO2 and CH4 as reflected in the isotherm shape and phase composition at equilibrium. The formation pressure of hydrates is lower for mixtures than for pure gases, and the highest sorption capacity of each gas decreased in the mixture sorption either.     

A new gas chromatographic methodology for the estimation of the composition of binary gas mixtures

The relatively new technique of reversed-flow gas chromatography (RFGC) is used to determine the diffusion coefficients of pure gases into gas mixtures (D(mix)(exp)). The pure gases are CO and CO(2), and the mixtures consist of H(2) and He in various volume percentage compositions. A linear regression analysis of D(mix)(exp) of CO and CO(2) in various mixtures of H(2) and He against the percentage composition (X(H2) or X(He)) of the mixtures at different temperatures results in an empirical equation relating D(mix)(exp) to the corresponding theoretical values of the diffusion coefficients of CO and CO(2) in the pure gases H(2) and He, as they are calculated from the Fuller-Schettler-Giddings equation. The empirical equation shows that the diffusion coefficient of an analyte gas in a gas mixture is the partial sum of its diffusion coefficients in the component gases, therefore making possible the determination of the mole fractions of the components of the mixture. The found percentage volume compositions are very close to those determined independently by routine gas chromatography, indicating that the proposed RFGC methodology could be successfully applied to the accurate determination of the volume composition of binary gas mixtures.     

Predicting hydrate forming pressure of pure alkanes in the presence of inhibitors

An inherent problem with natural gas production or transmission is the formation of gas hydrates, which can lead to safety hazards for production/transportation systems, and substantial economic risks. Hydrate inhibition with different inhibitors such as, methanol, ethylene glycol （EG）, triethylene glycol （TEG）, and sodium chloride solution continues to play a critical role in many operations. An understanding of when the hydrates form in the presence of these hydrate inhibitors, is therefore necessary to overcome hydrate problems. Several thermodynamic models have been proposed for predicting the hydrate formation conditions in aqueous solutions containing methanol/glycols and electrolytes. However, available models have limitations that include the types of liquid, compositions of fluids, and inhibitors used. The aim of this study is to develop a simple-to-use correlation for accurate prediction of hydrate-forming pressures of pure alkanes in the presence of different hydrate inhibitors, where the obtained results illustrate good agreement with the reported experimental data.     

Structures of hydrocarbon hydrates during formation with and without inhibitors

The formation of hydrates from a methane-ethane-propane mixture is more complex than with single gases. Using nuclear magnetic resonance (NMR) and high-pressure powder X-ray diffraction (PXRD), we have investigated the structural properties of natural gas hydrates crystallized in the presence of kinetic hydrate inhibitors (KHIs), two commercial inhibitors and two biological ice inhibitors, or antifreeze proteins (AFPs). NMR analyses indicated that hydrate cage occupancy was at near saturation for controls and most inhibitor types. Some exceptions were found in systems containing a new commercial KHI (HIW85281) and a recombinant plant AFP, suggesting that these two inhibitors could impact the kinetics of cavity formation. NMR analysis confirmed that the hydrate composition varies during crystal growth by kinetic effects. Strikingly, the coexistence of both structures I (sI) and II (sII) were observed in NMR spectra and PXRD profiles. It is suggested that sI phases may form more readily from liquid water. Real time PXRD monitoring showed that sI hydrates were less stable than sII crystals, and there was a conversion to the stable phase over time. Both commercial KHIs and AFPs had an impact on hydrate metastability, but transient sI PXRD intensity profiles indicated significantly different modes of interaction with the various inhibitors and the natural gas hydrate system.     

A model for prediction of gas hydrate formation conditions in aqueous solutions containing electrolytes and/or alcohol

In this contribution, a generalized method for predicting gas hydrate formation conditions in the presence of aqueous solutions, HL₁V calculations, is developed. Each phase is characterized separately. In this respect, the equation of Nasrifar et al. for calculating the activity of water in the presence of electrolytes and an alcohol is extended to water in various mixtures of electrolytes, alcohol and dissolved gas. An equation for calculating the activity of water in the presence of mixtures of electrolytes is also developed. These equations are then used to characterize the aqueous phase in HL₁V calculations. The modified Patel–Teja equation of state is used to characterize the vapour phase and the statistical model of van der Waals and Platteeuw for the hydrate phase. The proposed model is then compared with experimental results and other available models. No adjustable or curve-fitting parameters are used. The agreement with experimental results is very good. The comparison with other models also indicates that the proposed model predicts incipient hydrate formation conditions as good as the other models and in most cases even better.     

Hydrate-based carbon dioxide capture from simulated integrated gasification combined cycle gas

The equilibrium hydrate formation conditions for CO2/H2 gas mixtures with different CO2 concentrations in 0.29 mol% TBAB aqueous solution are firstly measured.The results illustrate that the equilibrium hydrate formation pressure increases remarkably with the decrease of CO2 concentration in the gas mixture.Based on the phase equilibrium data,a three stages hydrate CO2 separation from integrated gasification combined cycle (IGCC) synthesis gas is investigated.Because the separation efficiency is quite low for the third hydrate separation,a hybrid CO2 separation process of two hydrate stages in conjunction with one chemical absorption process (absorption with MEA) is proposed and studied.The experimental results show H2 concentration in the final residual gas released from the three stages hydrate CO2 separation process was approximately 95.0 mol% while that released from the hybrid CO2 separation process was approximately 99.4 mol%.Thus,the hybrid process is possible to be a promising technology for the industrial application in the future.     

Effect of various types of equations of state for prediction of simple gas hydrate formation with or without the presence of kinetic inhibitors in a flow mini-loop apparatus

This paper compares the effects of using various types of equations of state (PR,¹ SRK,² ER,³ PT⁴ and VPT⁵) on the calculated driving force and rate of gas consumption based on the Kashchiev and Firoozabadi model for simple gas hydrate formation for methane, carbon dioxide, propane and iso-butane with experimental data points obtained in a flow mini-loop apparatus with or without the presence of kinetic inhibitors at various pressures and specified temperatures. For this purpose, a laboratory flow mini-loop apparatus was set up to measure gas consumption rate when a hydrate forming substance (such as C₁, C₃, CO₂ and i-C₄) is contacted with water in the presence or absence of dissolved inhibitor under suitable temperature and pressure conditions. In each experiment, a water blend saturated with pure gas is circulated up to a required pressure. Pressure is maintained at a constant value during experimental runs by means of the required gas make-up. The total average absolute deviation was found to be 15.4%, 16.3%, 15.8%, 17.8% and 17.4% for the PR, ER, SRK, PT and VPT equations of state for calculating gas consumption in simple gas hydrate formation with or without the presence of kinetic hydrate inhibitors, respectively. Comparison results between the calculated and experimental data points of gas consumption were obtained in flow loop indicate that the PR and ER equations of state have lower errors than the SRK, VPT and PT equations of state for this model.     

Intensification of the performance of kinetic inhibitors in the presence of polyethylene oxide and polypropylene oxide for simple gas hydrate formation in a flow mini-loop apparatus

The main objective of the present work is enhancement of the performance of gas hydrate kinetic inhibitors in the presence of polyethylene oxide (PEO) and polypropylene oxide (PPO) for simple gas hydrate formation in a flow mini-loop apparatus. PEO and PPO are high molecular weight polymers that are not kinetic inhibitors by their self. For this investigation, a laboratory flow mini-loop apparatus was set up to measure the induction time and rate of gas hydrate formation when a hydrate-forming substance (such as C1, C3, CO₂ and i-C4) is contacted with water containing dissolved inhibitor in presence or absence of PEO or PPO under suitable temperature and pressure conditions. In each experiment, water containing inhibitors blend saturated with pure gas is circulated up to a required pressure. Pressure is maintained at a constant value during experimental runs by means of required gas make-up. The effect of PEO and PPO on induction time and gas consumption during hydrate formation is investigated in the presence or absence of PVP (polyvinylpyrrolidone) and l-tyrosine as kinetic inhibitors. Results were shown that the induction time is prolonged in the presence of PEO or PPO compared to the inhibitor only. Inclusion of PPO into a kinetic hydrate inhibitor solution shows a higher enhancement in its inhibiting performance compare to PEO. Thus, the induction time for simple gas hydrate formation in presence of kinetic hydrate inhibitor with PPO is higher, compare to kinetic hydrate inhibitor with PEO.     

A novel correlation for estimation of hydrate forming condition of natural gases

An inherent problem with natural gas production or transmission is the formation of gas hydrates, which can lead to safety hazards to production/ transportation systems and to substantial economic risks. Therefore, an understanding of conditions where hydrates form is necessary to overcome hydrate related issues. Over the years, several models requiring more complicated and longer computations have been proposed for the prediction of hydrate formation conditions of natural gases. For these reasons, it is essential to develop a reliable and simple-to-use method for oil and gas practitioners. The purpose of this study is to formulate a novel empirical correlation for rapid estimation of hydrate formation condition of sweet natural gases. The developed correlation holds for wide range of temperatures (265–298 K), pressures (1200 to 40000 kPa) and molecular weights (16−29). New proposed correlation shows consistently accurate results across proposed pressure, temperature and molecular weight ranges. This consistency could not be matched by any of the widely accepted existing correlations within the investigated range. For all conditions, new correlation showed average absolute deviation to be less than 0.2% and provided much better results than the widely accepted existing correlations.     

Initial estimation of hydrate formation temperature of sweet natural gases based on new empirical correlation

Production,processing and transportation of natural gases can be significantly affected by clathrate hydrates.Knowing the gas analysis is crucial to predict the right conditions for hydrate formation.Nevertheless,Katz gas gravity method can be used for initial estimation of hydrate formation temperature (HFT) under the circumstances of indeterminate gas composition.So far several correlations have been proposed for gas gravity method,in which the most accurate and reliable one has belonged to Bahadori and Vuthaluru.The main objective of this study is to present a simple and yet accurate correlation for fast prediction of sweet natural gases HFT based on the fit to Katz gravity chart.By reviewing the error analysis results,one can discover that the new proposed correlation has the best estimation capability among the widely accepted existing correlations within the investigated range.     

Solution and Diffusion Behavior of Pure Gases and Gas Mixtures in Glassy Polymer Membranes

A general model for the solution and diffusion behavior in pure gas-polymer membrane systems and gas mixture-polymer membrane systems has been developed. Proved by experiments on different glassy and rubbery polymer membranes at various temperatures and pressures, this model can achieve the prediction of permeation behavior of pure gases and gas mixtures in polymer membranes only using the model parameters obtained from experiments on pure gases. The calculated results are in good agreement with experimental.     

Considering multiple occupancy of cavities in clathrate hydrate phase equilibrium calculations

One of the major assumptions of the original van der Waals–Platteeuw (vdWP) model is the single occupancy of hydrate cavities. In this work, the vdWP model is modified to also account for multiple occupancies of hydrate cavities by small molecules. The developed model is evaluated by calculating the hydrate equilibrium conditions with either oxygen or nitrogen as guest molecules in pure form, as well as mixtures of nitrogen and propane (molecules of these pure gases and those in (nitrogen + propane) have double occupancy in large cavities of structure II up to a certain concentration of propane). The results of this modified model show good agreement with the experimental data reported in the literature.     

On the thermodynamic stability and structural transition of clathrate hydrates

Gas mixtures of methane and ethane form structure II clathrate hydrates despite the fact that each of pure methane and pure ethane gases forms the structure I hydrate. Optimization of the interaction potential parameters for methane and ethane is attempted so as to reproduce the dissociation pressures of each simple hydrate containing either methane or ethane alone. An account for the structural transitions between type I and type II hydrates upon changing the mole fraction of the gas mixture is given on the basis of the van der Waals and Platteeuw theory with these optimized potentials. Cage occupancies of the two kinds of hydrates are also calculated as functions of the mole fraction at the dissociation pressure and at a fixed pressure well above the dissociation pressure.     

Excess Gibbs potential model for multicomponent hydrogen clathrates

A new thermodynamic calculation procedure is introduced to predict the equilibrium conditions of multicomponent gas hydrates containing hydrogen. This new approach utilizes an excess Gibbs potential term to account for second- or higher-order water-cavity distortions due to the presence of multiple guest species. The excess Gibbs potential describes changes in reference chemical potentials according to different compositions of guest mixtures in the hydrate phase. To determine the equilibrium conditions of multicomponent gas hydrates, the excess Gibbs potential term is incorporated to the Lee-Holder model along with the Zele-Lee-Holder cell distortion model. For binary gas hydrates between hydrogen and the other gas molecule, the predicted equilibrium pressure deviates within 10-20% from the experimental value. For the ternary and quaternary mixture hydrates, the model prediction is reasonably good but its error increases with increasing pressure and temperature under the presence of THF.     

Influence of Gas Supply Changes on the Formation Process of Complex Mixed Gas Hydrates

Natural gas hydrate occurrences contain predominantly methane; however, there are increasing reports of complex mixed gas hydrates and coexisting hydrate phases. Changes in the feed gas composition due to the preferred incorporation of certain components into the hydrate phase and an inadequate gas supply is often assumed to be the cause of coexisting hydrate phases. This could also be the case for the gas hydrate system in Qilian Mountain permafrost (QMP), which is mainly controlled by pores and fractures with complex gas compositions. This study is dedicated to the experimental investigations on the formation process of mixed gas hydrates based on the reservoir conditions in QMP. Hydrates were synthesized from water and a gas mixture under different gas supply conditions to study the effects on the hydrate formation process. In situ Raman spectroscopic measurements and microscopic observations were applied to record changes in both gas and hydrate phase over the whole formation process. The results demonstrated the effects of gas flow on the composition of the resulting hydrate phase, indicating a competitive enclathration of guest molecules into the hydrate lattice depending on their properties. Another observation was that despite significant changes in the gas composition, no coexisting hydrate phases were formed.     

An efficient model for the prediction of CO2 hydrate phase stability conditions in the presence of inhibitors and their mixtures

A thermodynamic model for the prediction of CO₂ hydrate phase stability conditions in the presence of pure and mixed salts solutions and various ionic liquids (ILs) is developed. In the proposed model van der Waals and Platteeuw model is used to compute the hydrate phase, Peng–Robinson equation of state (PR-EoS) for the gas phase and the Pitzer–Mayorga–Zavitsas-Hydration model is employed to calculate the water activity in the liquid water phase. This model is an extension of the model developed by Tumba et al. (2011) for the prediction of methane and CO₂ hydrate phase stability conditions in the presence of tributylmethylphosphonium methylsulfate IL solution. Shabani et al. (2011) mixing rule is modified by incorporating the water–inhibitor (salt/IL) interaction parameter to calculate the water activity in mixed salt solutions. The model predictions are also calculated using the Pitzer–Mayorga model separately and compared with predictions of the developed model. The model predictions are compared with experimental results on the phase stability of CO₂ hydrate in the presence of ILs, pure and mixed salts as reported in literatures. The ILs are chosen from imidazolium cationic family with various anion groups such as bromide (Br), tetrafluoroborate (BF₄), trifluoromethanesulfonate (TfO), and nitrate (NO₃) and the common salts such as NaCl, KCl and CaCl₂. Good agreement between the developed model predictions and the literature data is observed. The overall average absolute deviation (AARD%) with Pitzer–Mayorga–Zavitsas-Hydration model is observed to be within ±1.36% while Pitzer–Mayorga model accuracy were about ±1.44 %. Further, the model is extended to calculate the inhibition effect of selected inhibitors (ILs and salts) on CO₂ hydrate formation.     

Comparison of three models for permeation of CO2/CH4 mixtures in poly(phenylene oxide)

Two models for the permeability of pure gases have been extended to include binary gas mixtures. The first is an extension of a pure gas permeability model, proposed by Petropoulos, which is based on gradients of chemical potential. This model predicts the permeability of components in a gas mixture solely on the basis of competition for sorption sites within the polymer matrix. The second mixed gas model follows an earlier analysis by Barrer for pure gases which includes the effects of saturation of Langmuir sites on the diffusion as well as the sorption processes responsible for permeation. This generalized “competitive sorption/diffusion” model includes the effect of each gas component on the sorption and diffusion of the other component in the mixture. The flux equations from these two models have been solved numerically to predict the permeability of gas mixtures on the basis of pure gas sorption and transport parameters. Both the mixed gas Petropoulos and competitive sorption/diffusion model predictions are compared with predictions from the earlier simple competitive sorption model based on gradients of concentration. An analysis of all three models is presented for the case of CO₂/CH₄ permeability in poly(phenylene oxide) (PPO). As expected, the competitive sorption/diffusion model predicts lower permeability than either of the models which consider only competitive sorption effects. The permeability depression of both CO₂ and CH₄ predicted by the competitive sorption/diffusion model is roughly twice that predicted by the competitive sorption model, whereas the mixed gas Petropoulos model predictions for both gases lie between the other two model predictions. For the PPO/CO₂/CH₄ system, the methane permeability data lie above the predictions of all three models, whereas CO₂ data lie below the predictions of all models. Consequently, the competitive sorption/diffusion model gives the most accurate prediction for CO₂, while the simple competitive sorption model is best for methane. The effects of mixed gas sorption, fugacity, and CO₂-induced dilation were considered and do not explain the inaccuracies of any of the models. The relatively small errors in mixed gas permeability predictions using either of the three models are likely to be related to “transport plasticization” of PPO owing to high levels of CO₂ sorption and its effect on polymer segmental motions and gas diffusivity.     

Kinetics and stability of CH4-CO2 mixed gas hydrates during formation and long-term storage.

The formation of CH4-CO2 mixed gas hydrates was observed by measuring the change of vapor-phase composition using gas chromatography and Raman spectroscopy. Preferential consumption of carbon dioxide molecules was found during hydrate formation, which agreed well with thermodynamic calculations. Both Raman spectroscopic analysis and the thermodynamic calculation indicated that the kinetics of this mixed gas hydrate system was controlled by the competition of both molecules to be enclathrated into the hydrate cages. However, the methane molecules were preferentially crystallized in the early stages of hydrate formation when the initial methane concentration was much less than that of carbon dioxide. According to the Roman spectra, pure methane hydrates first formed under this condition. This unique phenomenon suggested that methane molecules play important roles in the hydrate formation process. These mixed gas hydrates were stored at atmospheric pressure and 190 K for over two months to examine the stability of the encaged gases. During storage, CO2 was preferentially released. According to our thermodynamic analysis, this CO2 release was due to the instability of CO2 in the hydrate structure under the storage conditions.