Chun-Gang Xu
ab and
Xiao-Sen Li
*ab
aKey Laboratory of Renewable Energy and Gas Hydrate, Guangzhou Institute of Energy Conversion, Chinese Academy of Sciences, Guangzhou 510640, People's Republic of China
bGuangzhou Center for Gas Hydrate Research, Chinese Academy of Sciences, Guangzhou 510640, People's Republic of China. E-mail: lixs@ms.giec.ac.cn; Fax: +86-20-87034664
First published on 24th February 2014
Hydrate-based CO2 separation and capture from gas mixtures containing CO2 has gained growing attention as a new technology for gas separation, and it is of significance for reducing anthropogenic CO2 emissions. Previous studies of the technology include the thermodynamics and kinetics of hydrate formation/dissociation, hydrate formation additives, analytical methods, separation and capture progress, equipment and applications. Presently, the technology is still in the experimental research stages, and there are few reports of industrial application. This review examines research progress in the hydrate formation process and analytical methods with a special focus on laboratory studies, including the knowledge developed in analog computation, laboratory experiments, and industrial simulation. By comparing the various studies, we propose original comments and suggestions on further developing hydrate-based CO2 separation and capture technology.
Gas hydrates are non-stoichiometric compounds composed of water molecules and small gas molecules. Examples of gas hydrates include methane (CH4), carbon dioxide (CO2), nitrogen (N2), hydrogen (H2), sulfur dioxide (SO2), hydrogen sulfide (H2S), ethane (C2H6), propane (C3H8), iso-butane (i-C4H10), ethylene (C2H4), propylene (C3H6).5,6 Water molecules (hosts) connect with each other by hydrogen-bonds to form cavities, and small gas molecules (guests) are stabilized in the cavities by van der Waals interaction forces.7 Currently, the most common gas hydrate structures include structure I (sI), structure II (sII), structure H (sH) and semi-clathrate (sc).8,9 In general, the gas hydrate structure is mainly determined by size of the gas molecule if a single gas acts as a single guest. However, the structure is also affected by the composition and/or pressure when gas mixtures with multiple components act as guests.10 The structures of sI and sII hydrates were first determined by single crystal X-ray crystallography in the early 1950s. After the structures were discovered, van der Waals and Patteeuw were the first two people to describe the hydrates in terms of stability and composition.11
In 1810, Humphrey Davy first reported gas hydrates in the Bakerian lecture to the Royal Society.12 Hammerschmidt suggested that it was methane hydrate blocking gas pipelines rather than ice in 1934.13 However, researchers didn’t investigate the hydrates existing in nature until 50 years later.7,14 The effects of gas hydrates on energy and environmental applications have been intensively researched since the early 1990s. The applications include production and transportation of gas in subsea flow assurance, the potential energy recovery from naturally occurring gas hydrates deposits, the storage of new fuels (natural gas or hydrogen) in hydrate materials, the role of gas hydrates in environmental safety (stability of seafloor and climate change), gas separation and purification from gas mixtures, and desalination of sea water.8,15–17
Hydrate-based gas separation is based on the selective partitioning of the components in the hydrate phase and in the gas phase.3,18–20 At the same temperature, different gases have their individual equilibrium hydrate formation pressures. A gas with a relatively low equilibrium hydrate formation pressure at a certain temperature is expected to be preferentially entrapped in water cavities to form a gas hydrate with higher thermodynamic stability, and it results in a gas-rich hydrate phase while the residual gas phase is gas-poor.21–23 Relative to the other gas separation technologies, gas hydrate separation technology has the following advantages: (1) simple process, (2) low investment, (3) low material and energy loss, and (4) environmental friendliness.24 Previous studies focused on the separation of CO2, CH4,25–32 H2,33–35 N2,36 oil gas,37,38 and other greenhouse gases.39–44 The main research content includes kinetics and thermodynamics of gas hydrate formation/dissociation, separation and capture processes, analytical methods, and equipment and applications. For example, the studies of CO2 separation focus on the equilibrium conditions of hydrate formation, hydrate formation additive/promoter approaches, and separation equipment and processes. This review paper examines research progress in the hydrate formation process and analytical methods with a special focus on laboratory studies of hydrate-based CO2 separation, and its content includes the following parts: equilibrium hydrate formation conditions, additives for forming CO2 hydrates, molecular-level measurements of the hydrates containing CO2, process and apparatus, cost and comparison, and conclusions.
In fact, all prediction models are established based on parameters (e.g. gas fugacity, Langmuir constant) and correlations between the parameters. There are five major methods to determine the correlations. The first is the K-value method, which utilizes vapor–solid equilibrium constants to predict hydrate formation conditions.56,57 The second is a gas-gravity plot developed by Katz.58 The gas-gravity is defined as the apparent molecular weight of a gas mixture divided by the apparent molecular weight of air. The method is a simple graphical technique and it is useful for an initial estimate of hydrate formation conditions. The hydrate formation chart is made according to the limited experimental data and calculations based on the K-value method. But, the statistical accuracy analysis (Sloan59) showed that the method is inaccurate, and different gas mixtures can lead to about 50% deviation in predicted pressure for the same gas-gravity. However, copious amounts of experimental data have been collected in the past 60 years, and a more accurate gas-gravity plot may be developed based on the data. Therefore, Holder et al.60 and Makogon61 developed empirical correlations for selected pure gases, and this is the third method. Kobayashi et al.62 developed a new empirical equation for hydrate formation conditions of natural gases. However, the empirical equations have their limitation of temperature and pressure, and the method of the gas-gravity plot is seldom used for predicting the equilibrium conditions of gas hydrates containing CO2. The fourth method for determining correlations between parameters involves a chart of permissible expansion, and it is based on the range of the permissible expansion that a natural gas can undergo without the possibility of hydrate formation. The chart of permissible expansion is drawn based on the gas-gravity chart using the Joule–Thomson cooling curve.58 The method is suitable for rough design of throttling in valves and chokes where gas expansion normally occurs. However, it is not suitable for natural gases containing more CH4. An average error of 10% can be obtained in general. The fifth is a statistical thermodynamic approach, developed by van der Waals and Platteeuw,11 predicting the equilibrium hydrate formation temperatures and pressures on the criterion ΔμH = ΔμW at equilibrium.60 Thus, from a statistical mechanics point of view, those constants that can impact the chemical potentials between the hypothetical empty and fully filled hydrate lattice become key for the approach.11,63–65
The Langmuir constant is quite important, and different models have been developed to obtain more accurate values of the Langmuir constant, for example, the Kihara potential model.66–68 Parrish and Prausnitz developed a correlation using the Kihara potential and experimental data of the hydrate formation, and the accuracy of the correlation is approximately 0.2%.47 The parameters of the Kihara potential model are empirically regressed from the experimental data of the cavity occupancy and the phase equilibrium. Although the Kihara potential models can reproduce the experimental data, its capability for extrapolation is generally poor. For instance, the Parrish–Prausnitz model47 and CSMHYD program7 can only predict phase equilibria below 40–50 MPa. Thus, the models cannot predict the structure changes in mixture hydrates (e.g., the transition from sI to sII in CH4–C2H6 hydrate10). Furthermore, Tee et al.69 found that the Kihara potential surface calculated from the experimental data was inconsistent with that calculated from the second virial coefficient and viscosity data, and the inconsistency was due to the pre-treatment of the model overrating the occupancy fraction of non-spherical guests in a small cage,51,70 i.e., the Kihara potential model needs further optimization, especially in accurately describing the interaction between guest and water molecules. Accurate intermolecular potential between a guest molecule and a water molecule can be directly obtained by an ab initio quantum chemical method, and the intermolecular potential is considered to be strongly angle dependent.71–75 The ab initio potential model can be used for predicting the hydrate number and the cage occupancies. Duan et al.76 successfully used an atomic site–site Lennard-Jones formula plus an electrostatic term to fit the ab initio intermolecular potential energy surface of the CO2–H2O complex to account for the angle dependent molecular potential angle. They improved the model and predicted the equilibrium pressure of CO2-hydrate in a wide T–P range with absolute average deviation less than 3%.
The best method for determining the hydrate formation conditions is to experimentally measure the hydrate formation at the temperature, pressure and composition of interest. However, it is impossible to satisfy the infinite number of conditions for which measurement are needed. Thus, the hydrate formation prediction methods need to interpolate between the measurements. However, such experimental measurements are both time consuming and expensive. Therefore, a comprehensive artificial neural network model (ANN model) was developed to enable the user to accurately predict hydrate formation conditions for a given gas mixture without having to make experimental measurements.77
Previous studies on predicting CO2 hydrate phase equilibrium are summarized in Table 1. The most recent predictions are rooted in the vdWP theory for hydrates.11 Although the methods mentioned above have their individual characteristics for predicting the equilibrium conditions for forming or dissociating the different gas hydrates, they all have their limitations. To solve these problems, an integrated method combined with the above different methods must be developed.
Authors | Temperature (K) | Pressure (MPa) | Study | Na |
---|---|---|---|---|
a Number of measurements. | ||||
Deaton & Frost184 | 273–283 | 1.3–4.3 | K-charts, giving the vapor–solid equilibria for natural gases (including pure gases or gas mixtures) at lower than 273.15 K and higher than 273.15 K. | 19 |
Carson & Katz56 | 277–283 | 2.0–4.5 | Katz method, using vapor–solid equilibrium constants to predict the hydrate formation conditions. Katz correlation is not recommended above 100–150 MPa, depending on the composition of the gas mixtures. | 15 |
Katz58 | 273–322 | 0.2–42.0 | Method of gas-gravity plots, which relate the hydrate formation pressure and temperature to gas-gravity. The method was useful for an initial estimate of hydrate formation conditions and the prediction is rough. | 128 |
van der Waals & Platteeuw11 | van der Waals–Platteeuw model, which was based on a statistical thermodynamic approach, accounting for the interactions between gas molecules and water molecules forming gas hydrates. | |||
Larson186 | 257–283 | 0.5–4.5 | Predicted the equilibrium hydrate formation conditions of CO2 hydrates. | 45 |
Miller & Smythe192 | 151–193 | 0–0.000022 | Dissociation pressure of CO2 hydrate with equations for CO2 hydrate dissociation pressures and vapor pressures. | 8 |
Robinson & Mehta188 | 274–283 | 1.3–4.5 | The conditions for initial hydrate formation in the system of CO2–C3H8–H2O over a wide concentration range for the hydrate-water-rich liquid–gas phase region were measured and predicted in terms of solid–vapor K-factor. | 7 |
Falabella193 | 148.8–240.4 | 0.02–0.1 | At a low pressure range, hydrates of CH4, C2H6, C2H4, C2H2 and CO2 are involved. The van der Waals–Platteeuw model was employed to predict the equilibria associated with experimental measurements. | 5 |
Ng & Robinson78,190 | 279–284 | 2.7–14.5 | A modification of the Parrish and Prausnitz program, predicting hydrate forming conditions for pure gases in the presence of up to 20 wt% methanol solutions. | 9 |
Holder et al.60 | Empirical correlations developed in different forms and with various numbers of parameters. | |||
Adisasmito et al.85 | 273–288 | 1.2–11.0 | Verifying the work done by Unruh and Katz and by Berecz and Balla-Achs by experimental measurements. | 9 |
Englezos63 | 269–281 | 1.1–4.3 | Because the solubility of CO2 in salt solutions cannot be computed accurately using rigorous thermodynamic models, thus CSMHYD in conjunction with Trebble–Bishnoi equation is adopted to predict the incipient CO2 hydrate formation pressure in NaCl solutions and the average deviation is around 7.2%. | 57 |
Dholabhai et al.83 | 273–279 | 1.3–2.5 | Coupling model of statistical thermodynamic model of van der Waals and Platteeuw with coefficient models. Equilibrium conditions of CO2 hydrate in pure water and single and mixed electrolytes. | 4 |
Englezos & Hall194 | 275–283 | 1.5–4.2 | CSMHYD model predicting and measuring CO2 hydrate formation pressure in electrolyte, water-soluble polymers and montmorillonite. | 6 |
Breland & Englezos195 | 275–280 | 1.6–3.0 | An isothermal pressure search method is employed to measure the incipient equilibrium data for CO2 hydrate in glycerol solutions (10, 20, 30 mass%). | 2 |
Tohidi et al.81 | 268–284 | 1.0–5.0 | The model based on a thermodynamic approach, in which an equation of state is combined with a modified Debye–Huckel electrostatic term, with only one adjustable parameter for the water-rich phase. Predicting phase equilibrium conditions for CO2 hydrates in presence of saline water. | |
Nakano et al.82 | 289–294 | 100–500 | The high-pressure phase equilibrium for CO2 in pure water and saturated liquid CO2 and Raman spectrum of CO2 hydrate. | 13 |
Fan & Guo196 | 264–284 | 0.5–5.0 | Hydrate phase equilibria for CO2–CH4, CO2–C2H6, CO2–N2, CO2–CH4–C2H6–N2 in pure water and 10 mass% NaCl solution. | 13 |
Wendland et al.79 | 270–305 | 0.5–7.0 | Equilibrium conditions for CO2–H2O system focusing on three- and four-phase equilibria including fluid, hydrate and ice phases. The experimental data are correlated with the equations of Clausius–Clapeyron type. | 9 |
Seo & Lee80,197 | 272–284 | 1.5–5.0 | The three phase equilbria for aqueous phase containing CO2 and CH4 were predicted. The vapor and liquid phases were treated with SRK-EOS incorporated with the second-order modified Huron–Vidal (MHV2) mixed rule and hydrate phase with the van der Waals–Platteeuw model. | 26 |
Duan & Sun76,198 | 253–293 | 0.5–200 | Ab initio potential model predicting initial hydrate formation conditions for CH4 and CO2. Compared to the models employing the Kihara potential or Lennard-Jones potential, atomic site–site potentials was more accurate either under low pressure or under high pressure. | 20 |
Li & Englezos199 | 298–313 | 5.0–11.0 | SAFT equation of state was employed for the correlation and prediction of vapor–liquid equilibrium of eighteen binary mixtures. The predicted values agreed with the experimental data except for the H2O–CH3OH–CH4 at low CH3OH concentration in liquid phase of 60 wt%. | 4 |
Bahadori & Vuthaluru200 | 265–298 | 1.2–40.0 | A novel correlation based on the extracted data from Katz gas-gravity charts was proposed to predict the hydrate formation conditions for gases with weights of 16–29, the absolute deviation in average around 0.18%. | 44 |
Zeng & Li201 | 270–282 | 0.8–4.0 | PC-SAFT in conjunction with the van der Waals–Platteuw model and capillary Kelvin model was employed to predict CH4 and CO2 hydrates formation equilibrium conditions in porous media. The deviations for CH4 hydrate and CO2 hydrates were 1.66% and 2.76%, respectively. | 18 |
Sloan7 | 277–283 | Up to 400 MPa | Presenting an alternative set of K-values for the Katz method, which are dependent upon gas composition and hydrate structures, the valid pressure up to 400 MPa. | 20 |
Karamoddin & Varaminian97 | 260–330 | 0–5.0 | A method using PR equation of state and different mixing rules for predicting hydrate formation conditions for binary mixtures of CH4, C2H6, C3H8, i-C4H10, CO2 and H2S. The interaction parameters were optimized by using two phase equilibrium data (VLw), and then the optimized parameters were used for three phase equilibrium (VLwH) calculations. | 63 |
Elgibaly & Elkamel77 | 250–320 | 0.001–1000 | Firstly proposing ANN compositional models to predict hydrate formation conditions based on the K-value method and gas-gravity chart method. The ANN models consist of four models. The predicted results were more accurate than those obtained by the conventional models. The new models have to be updated by being retrained by using extra collected data. | 2387 |
Eslamimanesh et al. 202 | 279–295 | 0–120 | A thermodynamic model is proposed for representation/prediction of phase equilibria of semi-clathrate hydrates of the CO2, CH4, or N2 + tetra-n-butylammonium bromide (TBAB) aqueous solution. The van der Waals–Platteeuw (vdW–P) solid solution theory is used, revised with two modifications for evaluations of Langmuir constants and vapor pressure of water in the empty hydrate lattice, in which these values are supposed to be a function of TBAB concentration in aqueous solution. The Peng–Robinson (PR-EoS) equation of state along with re-tuned parameters of Mathias–Copeman alpha function is applied for calculation of the fugacity of gaseous hydrate former. For determination of the activity coefficient of the non-electrolyte species in the aqueous phase, the non-random two-liquid (NRTL) activity model is used. | 40 |
Eslamimanesh et al.203 | 276–294 | 2–500 | The model based on conventional Clapeyron model. Considering that the effect of pressure on molar volume of gas hydrate could not ignored, the “Clausius–Clapeyron” equation was improved from ![]() ![]() |
40 |
Shuker et al.95 | 270–295 | 0–2.5 | NN model was employed to predict hydrate formation conditions for pure gases and gas mixtures. The ANN model was relatively accurate for a given gas mixture and without having to do experimental measurements compared to the previous models of K-factor, HYDOFF, CSMHYD and HYSYS. | 20 |
Heydari et al.96 | 273–296.5 | 0.3–29.0 | ANN models were used for prediction of hydrate formation temperature. The results included the calculation relative error and R2 values between the experimental data and ANN predictions. The results showed that the ANN models could be improved based on the more collected data in a wider range of temperature and pressure. | 167 |
For the isothermal method, the system is first set at a pressure higher than the expected equilibrium pressure to form hydrates, and the system pressure must be maintained by an external reservoir for addition or withdrawal of gas. The pressure is reduced gradually after the hydrate formation. The equilibrium pressure can be obtained by the visual observation of the hydrate disappearance. Thus, the method requires an apparatus with windows mounted in both back and front.7,15,78–85
For the isobaric method, the system pressure is first kept constant by an external reservoir for the addition of gas, then the temperature is decreased until a significant addition of gas is noted from the external reservoir, which indicates the hydrate formation. Then, the temperature is slowly increased to dissociate the hydrates. During the process of hydrate dissociation, the system pressure must be kept constant. As the last of the hydrate disappears, the point is taken as the equilibrium temperature of the hydrate formation at the constant pressure. Similar to the isothermal method, the isobaric method must use an apparatus with windows. Fig. 1 shows a typical P–T diagram of the three-phase and four-phase equilibria in the system of CO2–H2O by the method.79
![]() | ||
Fig. 1 P–T diagram of the three-and four-phase equilibria in the carbon dioxide–water system: comparison of the results reported by Wendland et al.79 (·) with literature data by Kuenen and Robson183 (⊗); Deaton and Frost184 (×); Unruh and Katz185 (+); Larson186 (□); Takenouchi and Kennedy187 (○); Robinson and Mehta188 (Δ); Vlahakis et al.189 (▲); Ng and Robinson190 (◊); Nakayama et al.191 (□); Adisasmito et al.85 (▼). |
The isochoric method describes how the pressure in a closed vessel changes with the temperature and the phase transition. The temperature is first lowered from the vapor–liquid (V–L) region, and the isochoric cooling leads the pressure to slightly decrease. At a certain temperature, the hydrate forms, causing a remarkable pressure drop. Then the temperature is slowly increased to dissociate the hydrates, resulting in the pressure rising quickly. Then, we can find the point that is the intersection of the hydrate dissociation trace with the initial cooling trace, and this point is taken as the equilibrium point for the hydrate dissociation. The isochoric method is commonly used for hydrate formation at high pressure, as it does not need to be visually observed.
Based on the isochoric method, a method of isochoric step-heating (T-cycle method) was developed.19,86–89 The experimental results obtained by this method are quite reliable and repeatable. However, it takes more than 24 hours for a measurement.86 Herri and Kwaterski90 improved the isochoric step heating method by sampling small quantities of gas and liquid during the heating procedure, and analyzing them by gas chromatography, ionic chromatography and refractive index measurement. By the mass balance, the improved method (a modified step heating method) allows one to determine the gas and hydrate compositions.
Differential scanning calorimetry (DSC) has also been used to determine binary phase (solid/solid or solid/liquid) equilibrium data in various fields.91 Recently, DSC was applied to determine the equilibrium conditions of hydrate.92–94 The key is calibrating the temperature by measuring the melting temperatures of pure materials and then computing a correction function by taking into account the deviation between the measurements and the known melting temperatures at different heating rates within a special temperature range. Once the temperature calibration is complete, the DSC software can automatically correct the measured temperature according to the correction function. The method of DSC requires careful calibration, and the sample mass and heating rate must be chosen within reasonable limits.
Table 2 shows previous hydrate equilibrium measurements involving pure CO2 or gas mixtures containing CO2 in pure water, or electrolyte solution, or porous media systems. Presently, the hydrate equilibrium measurements mainly relate to the sequestration of pure CO2 in the seafloor, while CO2 separation and capture pertains to gas mixtures containing CO2. However, due to the equilibrium equipment limitations, the measurements are operated at a relatively narrow range of pressures. The equilibrium data at higher pressure (for example, >100 MPa) still needs computational models. As mentioned in Section 2.1, more accurate prediction needs more experimental data. The ANN method should be further developed because it integrates four different models and saves computational time across a wide range of temperatures and pressures.95–98
Author(s) | System | Study |
---|---|---|
Hashimoto et al.102 | H2–CO2–THF–water | Phase equilibria and Raman spectroscopic analysis for gas hydrate. |
Shin et al.103 | CO2–3M1B–water | Thermodynamic stability, spectroscopic identification and cage occupation of binary CO2 hydrates. |
CO2–THF–water | ||
CO2–DXN–water | ||
Kang & Lee21 | CO2–N2–water | Phase equilibrium measurements for CO2 recovery from flue gas using gas hydrate. |
CO2–N2–THF–water | ||
Dholabhai et al.84 | CO2–water | CO2 hydrate equilibrium conditions in aqueous solutions. |
CO2–electrolytes–water | ||
CO2–methanol–water | ||
Dholabhai & Bishnoi204 | CO2–CH4–electrolytes–water | Hydrate equilibrium conditions in aqueous electrolyte solutions. |
Matsui et al.205 | CO2–C2H6–water | Phase equilibrium for binary hydrate systems. |
CO2–CF4–water | ||
Makino et al.206 | CO2–CP–water | Phase equilibrium and structural transition in the CO2–CP mixed hydrates. |
Yang et al.207 | CO2–Doda glass–water | Characteristics of CO2 hydrate formation and dissociation in glass beads and silica gel. |
CO2–silica gel–water | ||
Li et al.116 | CO2–H2–TBAB–water | Phase equilibrium conditions for CO2–H2–TBAB–water mixed gas hydrate. |
Kumar et al.105 | CO2–H2–C3H8–water | Hydrate phase equilibrium conditions for CO2–H2–C3H8–water mixed gas hydrate. |
Belandria et al.208 | CO2–N2–water | Experimental and predicted hydrate phase equilibrium conditions for CO2–N2–water gas hydrate. |
Lee et al.209 | CO2–electrolyte–water | Phase equilibria and kinetic behavior of CO2 hydrate in electrolyte and porous media solutions. |
CO2–porous media–water | ||
Sabil et al.210 | CO2–THF–water | Phase equilibria in ternary system of CO2–THF–water. |
Maekawa211 | CO2–alcohols–water | Hydrate equilibrium conditions for CO2 hydrates in presence of alcohols, glycols, and glycerol. |
CO2–glycols–water | ||
CO2–glycerol–water | ||
Kim et al.155 | CO2–H2–TBAB–water | Hydrate-based CO2 capture for pre-combustion process in IGCC plant. |
Lin et al.212 | CO2–TBAB–water | Hydrate phase equilibrium and dissociation enthalpy for CO2–TBAB hydrate. |
Li et al.120 | CO2–TBAB (TBACl, TBAF) | Hydrate phase equilibrium for CO2 hydrate in presence of TBAB, TBACl, TBAF. |
Kang et al.100 | CO2–N2–THF–water | Hydrate phase equilibrium for CO2–N2–THF hydrate. |
Zhang et al.108 | CO2–H2–CP–water | Thermodynamic analysis of hydrate-based pre-combustion capture of CO2. |
Sugahara et al.213 | CO2–H2–water | Hydrate phase equilibria for CO2–H2 hydrate. |
Park et al.214 | CO2–N2–silica gel–water | Hydrate phase equilibrium and NMR analysis for CO2–N2–silica gel hydrate. |
Mayoufi et al.215 | CO2–TBPB–water | Hydrate phase equilibria for TBPB hydrate and CO2–TBPB hydrate. |
Mayoufi et al.216 | CO2–TBMAC–water | Hydrate phase equilibrium for CO2–TBMAC hydrate. |
Lee et al.22 | CO2–H2–THF–water | Hydrate phase equilibrium, gas consumed and CO2 separation efficiency in the process of pre-combustion capture of CO2. |
Lee et al.118 | CO2–CH4–N2–water | Thermodynamic stability, spectroscopic identification, and gas storage capacity of CO2–CH4–N2 hydrate. |
Nakano et al.82 | CO2–water–liquid CO2 | High-pressure phase equilibrium and Raman spectroscopic analysis for CO2 hydrate. |
Ohgaki et al.217 | CO2–CH4–water | Hydrate phase equilibrium for CO2–CH4 gas hydrate and replacement of CH4 by CO2. |
Bruusgaard et al.218 | CO2–CH4–water | Hydrate phase equilibrium for CO2–CH4 hydrate. |
Adisasmito et al.85 | CO2–CH4–water | PVT studies on dissociation conditions of CO2–CH4 hydrate. |
Belandria et al.219 | CO2–CH4–water | Phase equilibria in the CO2–CH4–H2O system has measured by the method of isochoric pressure-search method in the conditions of 233–373 K and up to 60 MPa. |
Beltran & Servio220 | CO2–CH4–water | PVT studies on dissociation conditions and composition measurement of gas phase. |
CO2–CH4–neohexane-emulsion | ||
Fan et al.221 | CO2–CH4–SDS–water | Experimental and modeling studies on the hydrate formation of CO2 and CO2-rich gas mixtures. |
Deschamps & Dalmazzone222 | CO2–N2–TBAB–water | Dissociation enthalpies and phase equilibrium for TBAB semi-clathrate hydrates of N2, CO2, N2–CO2 and CH4–CO2. |
CO2–CH4–TBAB–water | ||
Fan et al.223 | CO2–H2–TBAB–water | Efficient capture of CO2 from simulated flue gas by formation of TBAB or TBAF semi-clathrate hydrates. |
CO2–H2–THF–water | ||
Meysel et al.224 | CO2–N2–TBAB–water | Incipient equilibrium conditions for the formation of semi-clathrate hydrates from quaternary mixtures of (CO2–N2–TBAB–H2O). |
Mohammadi et al.225 | CO2–N2–TBAB–water | Experimental data for the hydrate dissociation conditions for the system comprising mixtures of CO2 (0.151/0.399 mole fraction) + N2 (0.849/0.601 mole fraction) + TBAB (0.05/0.15/0.30 mass fraction) in the conditions of 277.1–293.2 K and up to 16.21 MPa. |
Kamata et al.32 | CO2–H2S–TBAB–water | A high-pressure vessel of separation gas from gas mixture. |
Belandria et al.226 | CO2–H2–water | Molar compositions of carbon dioxide (and hydrogen) in the gas phase in equilibrium with gas hydrate and aqueous phases were measured for various (H2–CO2) gas mixtures t water systems in the temperature range of 273.6–281.2 K at pressures up to ∼9 MPa. The compositions of the gas phase were measured using an isochoric technique, in combination with the ROLSI capillary gas-phase sampling and a gas chromatography technique. |
The additives are classified into thermodynamic and kinetics types. The thermodynamic additives have the tendency to moderate the equilibrium conditions to higher temperature or lower pressure, and they often consist of organic compounds, including tetrahydrofuran (THF), propane (C3H8), cyclopentane (CP), and tetra-n-butyl-ammonium (bromide, fluoride or chloride; TBAB, TBAF and TBACl, respectively). Kinetic additives accelerate the hydrate formation and they typically consist of surfactants, including sodium dodecyl sulfate (SDS), and dodecyl trimethyl ammonium chloride (DTAC).7,99
Tetrahydrofuran (THF) solution has proved to be capable of significantly reducing the hydrate formation pressure at a given temperature.100–102 Kang et al.,100 Linga et al.,101 and Hashimoto et al.102 found that the equilibrium pressure of hydrates for both CO2–N2 and CO2–H2 mixtures in the presence of THF solution was considerably lower than that without the additive, and THF of 1.0 mol% proved to be the optimal concentration for CO2 separation from the mixture.22,100 However, THF can form sII THF-hydrates, competing with CO2 for occupying large cavities (51264).22,103 Therefore, THF cannot remarkably improve the gas consumption and CO2 separation efficiency although it can moderate the conditions of forming gas hydrates for either CO2–N2 or CO2–H2.104 Kumar et al.105 compared the hydrate equilibrium conditions for CO2–H2 and CO2–H2–C3H8, and found that 3.2 mol% C3H8 added to a CO2–H2 mixture reduced the pressure by approximately 50%. They verified that the addition of C3H8 into the CO2–H2 mixture reduced the hydrate phase equilibrium pressure without comprising the CO2 recovery and found that H2 existed in both sI hydrates formed by CO2–H2 mixture and sII hydrates formed by CO2–H2–C3H8 mixture (C3H8 of 2.6 mol%) by Raman spectroscopy.106 Cyclopentane (CP) can also reduce the hydrate phase equilibrium pressure. Zhang et al.107,108 found that the equilibrium of ternary CO2–H2–CP hydrates was significantly lower than that of ternary CO2–H2–THF at a vapor-phase CO2 mole fraction of 0.3204. However, because the CO2–H2–CP forms sII hydrate and CP occupies the large cavities of 51264 preferentially, CO2 can only compete with H2 to occupy the small cavities of 512, reducing the selectivity of CO2 over H2 in the hydrate phase.108 TBAB in water forms a semi-clathrate hydrate under moderate conditions. In the semi-clathrate hydrate system, the anions (Br−) are strongly incorporated with the host water lattice and a single TBA cation (TBA+) occupies four cavities, leaving dodecahedral cavities for small gaseous molecules.9,88,109–112 Arjmandi et al.,88 Oyama et al.113 and Duc et al.24 found that the equilibrium conditions for binary CO2–TBAB hydrates were considerably lower than those for pure CO2 hydrate. Furthermore, the hydrate phase equilibrium pressure shifts to lower with the increase of the TBAB concentration.32,114,115 TBAB of 0.29 mol% is considered as the optimum to recover CO2 from either flue gas or IGCC synthesis gas because TBAB of more than 0.29 mol% makes no more contribution to the CO2 recovery.19,23,88,116,117 However, adding TBAB reduces the gas consumption because the TBA+ occupies the big cavities in the sc hydrates. Therefore, there is some controversy over the use of TBAB as a gas hydrate promoter. Besides, both TBAF and TBACl can also reduce the hydrate equilibrium conditions.118–121 TBAF is rarely used because it is much more expensive than TBAB, although the effect of reducing pressure of TBAF is superior to that of TBAB.
SDS is widely used as a kinetic additive. Zhong and Rogers122 studied the SDS effects on gas hydrate formation and found the hydrate formation rate can be increased by multiple orders of magnitude in the presence of SDS solution or a related surfactant solution. Tajima et al.,123 Li and Chen,124 Rossi et al.,125 and Torre et al.126 further verified that the SDS of critical micelle concentration (CMC) has the best effect on enhancing the hydrate formation rate.
However, either using thermodynamic additives or kinetic additives cannot resolve all the problems (extreme conditions, low gas consumption, low hydrate formation rate and low CO2 recovery) of hydrate-based CO2 separation and capture from the gas mixtures. Therefore, researchers have focused on the synergistic effect of thermodynamic additives and kinetics additives. Li et al.117 added DTAC into the 0.29 mol% TBAB solution to investigate hydrate-based CO2 capture from the flue gas and considerably improved the formation conditions, the formation rate and CO2 recovery. Ricaurte et al.127 used THF in combination with SDS to investigate CO2 removal from a CH4–CO2 gas mixture by hydrate formation and found that the combinational additive decreased the hydrate formation pressure and improved the selectivity of CO2 capture. The combination of two thermodynamic additives or two kinetic additives has also been studied. Ding128 used SDS associated with an anionic fluorosurfactant (FS-62) (FS-62/SDS: 100/1000 ppm) as a joint additive to hydrate-based CO2 capture from a CO2–N2 mixture, but its effect on raising the gas consumption was limited. However, Li et al.19 found that the addition of CP into a 0.29 mol% TBAB solution could considerably improve hydrate-based CO2 separation from IGCC synthesis gas, and the result was further verified by powder X-ray diffraction (PXRD) spectroscopy.18
In fact, additives are still in the course of being screened, i.e., which kind of additive is the best for hydrate-based CO2 separation and capture from the gas mixtures is still undetermined. Currently, THF and TBAB are two of the most popular thermodynamic additives. Much work is needed to further screen additives, especially to screen one that can remarkably enhance the gas uptake.
Neutron diffraction studies are able to determine the positions of the guest and the host in a hydrate crystal, to trace the structural changes during the hydrate formation, and to measure the extent of guest–host interactions in a hydrate lattice.137–142 However, few studies on CO2 hydrate via neutron diffraction are reported. Henning et al.139,143 observed complete conversion from the hexagonal ice to the sI type CO2 hydrate as the temperature of the sample was slowly increased through the melting point of D2O.
Raman spectroscopy is a good way to identify the hydrate structure and composition because the Raman peak is determined by the inter-atomic vibration.144–146 Hydration number and relative cavity occupation can also be measured via Raman spectroscopy.103,147–152 Fig. 2 shows a typical Raman spectrum of CO2.153 The split peaks (located at 1276 cm−1 and 1384 cm−1, respectively) are caused by Fermi resonance effect corresponding to C–O symmetric stretching (ν1) and O–C–O bending (2ν2) modes of CO2 molecules. Due to the various conditions of measurements and hydrate structures, the Raman shifts of CO2 molecules may be changes in a small range.103,106,154 The Raman peaks of molecules of CH4, N2, H2 are around 2910 cm−1, 2325 cm−1 and 4120 cm−1, corresponding to their vibration modes of ν1 sym C–H stretching, ν1 sym N–N stretching and H–H oscillation (pure molecular vibration). Thus, it is quite simple to identify CO2 from CO2–CH4, CO2–N2, and CO2–H2 gas mixtures via Raman spectra.102,106,154–158
![]() | ||
Fig. 2 Raman spectra of solvated CO2 and CO2 hydrates. The black dotted lines located at the solvated CO2 peak positions show the peak shift associated with hydrate formation.153 |
NMR can be used to identify the hydrate structure and quantify the relative cavity occupancy. Presently, the study of NMR on hydrate-based CO2 separation focuses on the spectra of 1H, 129Xe and 13C.159–162 Among all the studies, 1H NMR has been used for ethane, propane, and isobutene hydrates;159 129Xe NMR has been used for identifying ratios of xenon atoms in small and large cavities;160–162 and 13C NMR has been applied to study hydrates of CO2, CH4 and C3H8.149,163–165 Seo et al.166 and Seo and Lee136 found, based on the 13C NMR spectra for CO2, that CO2 molecules occupied the large cavities of sI, and as the CO2 increased in the CO2–N2 vapor phase, the role of stabilizing both small and large cavities was transformed from N2 to CO2 molecules.
Duc et al.24 conducted experiments to separate CO2 from a CO2–N2 mixture in the presence of TBAB at under suitable operation conditions, and they proposed a continuous multi-stage separation process. However, hydrate-based CO2 separation becomes more difficult and the CO2 recovery becomes lower with the decrease of CO2 concentration in the gas mixtures, and the separation process with the single hydrate method cannot separate CO2 from gas mixtures such as flue gas or fuel gas completely and efficiently.168 Thus, various hybrid processes were proposed. In the hybrid processes, CO2 was first separated by the hydrate method, leaving lean-CO2 gas mixtures to be separated by other methods such as chemical adsorption, cryogenic separation, and membrane separation. Linga et al.167 proposed a hybrid process that combined the hydrate method and membrane separation method for separating CO2 from flue gas and fuel gas, respectively. Xu et al.169 conducted experiments to separate CO2 from IGCC synthesis gas by a hybrid process of two-stage hydrate separation in combination with a chemical absorption separation. Surovtseva et al.170 designed a process combined cryogenic and hydrate methods to capture CO2 from IGCC flue gases. Pure CO2 (>95%) could be obtained by the above hybrid processes, which were considered to be more efficient and economic compared to conventional CO2 separation methods.104,167,169–171 However, two problems retard even pilot-scale application of the hybrid processes: the gas hydrate capacity per volume of water (L) and the hydrate formation rate per unit of time (h).
The hydrate-based CO2 separation process does not use large absorber towers or steam-reboiled regenerators, so both the capital and operating costs of the process are primarily in refrigeration and compression systems. Before designing parameters for the process, researchers must focus on some parameters and key processes:
(1) Hydrate number of the hydrate: the number of water molecules required for removing a mole of CO2 has a great effect on the heat removal requirements in the reactor, and therefore the size of the refrigeration system.
(2) Slurry concentration: the amount of free water circulating in the system must be heated and cooled, and this also affects the size of the refrigeration system.
(3) Temperature of the reactor: the performance of the system is greatly affected by the temperature in either CO2 removal or heat transfer.
SIMTECHE designed a test apparatus for a plant of coal feed 5 kt per day. The mixed gases and conditioned water were metered to a reactor, mixing with the reactants, and the reactants flowed through a tail tube (ID: 4.8 mm). Under the conditions of around 10 MPa and 269 K and at the feed rates of CO2 (0.94 mol min−1) and H2O (28 mol min−1), the test results showed that the hydrate production rate was 0.271 mol min−1, which matched the design requirement.172 However, being subject to hydrate formation route, CO2 separation efficiency and gas–liquid-hydrate separation in a special pressurised vessel, there are as yet no mature hydrate-based CO2 separation process being applied in industry, even in a small-scale pilot application. Thus, we have a long way to go to achieve the final goal of hydrate-based CO2 separation and capture from gas mixtures in industry, especially for the aspects of the gas–liquid-hydrate separation processes, improvement of the CO2 separation efficiency and optimization of the process parameters.
![]() | ||
Fig. 3 Block flow diagram of the SIMTECHE CO2 capture process (for IGCC applications).177 |
![]() | ||
Fig. 4 Basic layout of flow rate reactors in the ETM system.177 |
Although it has been subject to gas–solid–liquid separation in special pressured vessel and low CO2 recovery, the continuous CO2 separation process is still immature. Furthermore, no complete set of equipment has been developed and utilized for CO2 hydrate-based separation and capture until now. Therefore, it is imperative to further develop the process associated with developing apparatus for hydrate-based CO2 separation and capture.
![]() | (1) |
The cost of disposing of CO2 consists of four factors, including separation (i.e. capture/separation of CO2 from combustion gases), compression, pipelining and injection (e.g. pumping and disposal wells). Capture/separation costs represent the largest financial impediment among the four factors, accounting for approximately three fourths of the total costs. Hence, it is necessary to develop efficient, cost-effective transportation and capture/separation technologies to allow large-scale use of geologic sinks. New CO2 separation technologies such as gas hydrate and membrane are developed under the driving force.
Spencer et al.181 first proposed a relatively comprehensive economic analysis on the basis of a hybrid process of on-stage hydrate in conjunction with chemical absorption proposed by SIMTCHE. The engineering basis of the analysis was divided into four parts: on-stream factor, major equipment costs estimated from ASPEN “Icarus”, installation factor, and contingencies. In the process, two-stage hydrate separation was carried out at 22 °F and 108 to 625 psia. Under the conditions for 90% removal of CO2, the cost for capturing one tonne CO2 from the flue gas was about 18 dollars. The cost was much lower than those by chemical absorption, membrane, cryogenic separation and solid physical adsorption. According to Wong et al.'s analysis, the range for the cost of capturing CO2 from flue gas using amine absorption is $30–$50 per tonne (t) of CO2, and the costs for the other technologies such as solid physical adsorption (including pressure swing adsorption (PSA) and temperature swing adsorption (TSA)), cryogenic separation, membranes, hybrid membrane/amine processes, electrical swing adsorption (ESA) and sorbent energy transfer system are higher.179,182 The gas hydrate technology has a certain advantage for the economy, however, there are some barriers to the technology, including the ability to release CO2 from the hydrate in an energy efficient manner, efficient capture of CO2, stale pre-hydrate, and trace contaminant interference with hydrate formation. The barriers limit the further development of the gas hydrate separation technology. Until now, there are quite few publications reporting on the process of CO2 separation by gas hydrate and the relevant cost analysis. Thus, in order to evaluate the CO2 separation cost by gas hydrate, the further development of gas hydrates, including kinetics and thermodynamics, is necessary, as well as CO2 hydrate separation processes and equipment.
The present computation models are mainly based on the van der Waals–Platteeuw model and are developed according to the various revisions of relevant parameters. Although the models are considered to be increasingly accurate, they still have individual limitations such as the limitation of pressure range and the dependency on the experimental data. ANN, a new computation model consisting of the four conventional models, has been developed to save computation time and to obtain more accurate predictions. The experimental data for equilibrium conditions for hydrates containing CO2 or gas mixtures have been obtained in different systems, including pure water, solutions with additives, and solutions with electrolytes in different measurement methods. The numerous data are the basis for hydrate-based CO2 separation and capture.
Neither thermodynamic additives nor kinetic additives can resolve all the problems of hydrate-based CO2 separation, thus, combination additives have developed. It was found that the combination additives of THF and SDS, TBAB and CP, THF and C3H8, etc., can considerably improve the CO2 separation compared to the single additives. However, until now, there is no conclusion to determine which kind of additive is the best for hydrate-based CO2 separation and capture from the gas mixtures.
The hydrate structures, compositions and the cage occupancies can be identified by diffraction methods and microscopic methods. However, due to the limitation of sampling, the molecular-level measurements are difficult to conduct. Thus, it is necessary to develop new apparatus or new ways to make the measurements simpler.
Because the equilibrium conditions become extreme as CO2 decreases in gas mixtures such as CO2–H2 or CO2–N2, it is quite difficult to thoroughly separate and capture CO2 from the gas mixtures by the hydrate method. Thus, integrated separation technologies such as hydrate/chemical absorption, hydrate/membrane and hydrate/cryogenic have been developed, and the experimental results indicate almost all the CO2 can be completely separated and captured. But, being subject to the hydrate formation method, CO2 separation efficiency and gas–liquid-hydrate separation in a special pressured vessel, there are no mature hydrate-based CO2 separation processes that are applied in industry, even in small scale–scale pilot applications.
New apparatuses that can promote the gas hydrate formation have been developed, including a new type of reactor, the continuous flow reactor. The objective of developing a new apparatus is expected to bring the hydrate-based CO2 separation from experiment to industrial practice. However, the present apparatus does not yet match the requirements.
This journal is © The Royal Society of Chemistry 2014 |