Paper: A Thermodynamic Analysis of Methanation Reactions of Carbon Oxides For The Production of Synthetic Natural Gas
Paper: A Thermodynamic Analysis of Methanation Reactions of Carbon Oxides For The Production of Synthetic Natural Gas
www.rsc.org/advances PAPER
A thermodynamic analysis of methanation reactions of carbon oxides for the
production of synthetic natural gas
Jiajian Gao,ab Yingli Wang,a Yuan Ping,a Dacheng Hu,a Guangwen Xu,a Fangna Gu*a and Fabing Su*a
Received 24th August 2011, Accepted 7th December 2011
DOI: 10.1039/c2ra00632d
Published on 30 January 2012. Downloaded on 9/17/2023 6:23:17 AM.
Synthetic natural gas (SNG) can be obtained via methanation of synthesis gas (syngas). Many
thermodynamic reaction details involved in this process are not yet fully understood. In this paper, a
comprehensive thermodynamic analysis of reactions occurring in the methanation of carbon oxides
(CO and CO2) is conducted using the Gibbs free energy minimization method. The equilibrium
constants of eight reactions involved in the methanation reactions were calculated at different
temperatures. The effects of temperature, pressure, ratio of H2/CO (and H2/CO2), and the addition of
other compounds (H2O, O2, CH4, and C2H4) in the feed gas (syngas) on the conversion of CO and
CO2, CH4 selectivity and yield, as well as carbon deposition, were carefully investigated. In addition,
experimental data obtained on commercial Ni-based catalysts for CO methanation and three cases
adopted from the literature were compared with the thermodynamic calculations. It is found that low
temperature, high pressure, and a large H2/CO (and H2/CO2) ratio are favourable for the
methanation reactions. Adding steam into the feed gas could alleviate the carbon deposition to a large
extent. Trace amounts of O2 in syngas is unfavourable for SNG generation although it can lower
carbon deposition. Additional CH4 in the feed gas almost has no influence on the CO conversion and
CH4 yield, but it leads to the increase of carbon formed. Introduction of a small amount of C2H4, a
representative of hydrocarbons in syngas, results in low CH4 yield and serious carbon deposition
although it does not affect CO conversion. CO is relatively easy to hydrogenated compared to CO2 at
the same reaction conditions. The comparison of thermodynamic calculations with experimental
results demonstrated that the Gibbs free energy minimization method is significantly effective for
understanding the reactions occurring in methanation and helpful for the development of catalysts
and processes for the production of SNG.
1. Introduction fuel cell and ammonia synthesis industry, methanation can also
be used for the removal of trace CO from H2-rich stream to
Recently, because of the increasing demand for natural gas and prevent catalyst poisoning.6,7 The catalysts and processes for
the reduction of greenhouse gases,1,2 interests have focused on methanation have been intensively explored and summarized in
the production of synthetic natural gas (SNG), which is the reviews.1,8–11
suggested as an important future energy carrier. The conven- Experimental results showed that the methanation reactions of
tional route for SNG production is based on the gasification of carbon oxides were sensitive to the many factors such as
coal or biomass to synthetic gas (syngas). After gas cleaning and temperature, pressure, composition of reactants (syngas).2,12,13
conditioning, the syngas containing H2, CO, CO2, H2O, CH4, Simultaneous methanation of CO and CO2 is frequently
and hydrocarbons3,4 can be converted to SNG by catalytic encountered in a coexistent system14,15 and there is a competitive
hydrogenation of carbon oxides (CO and/or CO2), otherwise reaction between CO and CO2 methanation.6,16 As we know, the
known as methanation. At the moment, there are at least 15 syngas obtained through gasification of coal and biomass
coal-to-SNG plants proposed in the U.S., which are in different inevitably includes high hydrocarbons such as C2H6, C2H4,
stages of development,5 and more than 20 coal-to-SNG plants and C2H2, and methanation of syngas involves many side-
are in construction and proposed in China because of China’s reactions that could result in complex by-products, which bring
abundant coal and poor natural gas. On the other hand, in the about a negative effect on the methanation process. The carbon
deposition always occurs on the surface of catalysts during
a
State Key Laboratory of Multiphase Complex Systems, Institute of methanation,17 resulting in the quick deactivation of catalysts,
Process Engineering, Chinese Academy of Sciences, Beijing, China. and this is one of the most serious problems in industry. It is
E-mail: [email protected] (F.Gu); [email protected] (F.Su);
Fax: +86-10-82544851; Tel: +86-10-82544850 reported that adding steam to the reactants could suppress
b
Graduate University of Chinese Academy of Sciences, Beijing, China carbon deposition.1 In the methanation industry, the recycling of
2358 | RSC Adv., 2012, 2, 2358–2368 This journal is ß The Royal Society of Chemistry 2012
View Article Online
occurring in such a complex reaction system. Thermodynamic products are given, the distribution of the products under a
equilibrium calculations of complex chemical systems based on minimum free energy can be established by using a general
the Gibbs free energy minimization can give answers to mathematical technique without knowing the chemistry of the
numerous questions, including the type of thermodynamically reactions. This method also can be used for condensed species,27
stable reaction products produced, together with their selectivity non-reacting species as well as reacting species.28 Therefore, the
and yield, whether a chemical reaction proceeds endothermically Gibbs free energy minimization method is widely used to
or exothermically, the impact of reaction parameters like perform thermodynamic analysis29–31 and should be more
temperature, pressure, or reactant ratios, and whether observed suitable for the present work because the methanation process
reaction temperatures are thermodynamically determined. By involves complex species and reactions. Here, the software
comparing the calculations with experimental results, the CHEMCAD was used for carrying out the calculation, and
identification of kinetic hindrances, i.e. thermodynamically meanwhile, the Soave–Redlich–Kwong method was employed as
allowed but somehow suppressed chemical reactions, should be the equation of state to determine fugacity in the calculation.
helpful in understanding the methanation process. It is noted The total Gibbs free energy GT of the system is given as:
that some thermodynamic analysis of methanation reactions has X
m X
m X fi
been done recently. Greyson and co-workers calculated the GT ~ ni mi ~ ni mH
i zRT ni ln (1)
i~1 i~1
PH
thermodynamics equilibrium compositions for CO methanation
system only for three reactions (CO + 3H2 « CH4 + H2O, CO +
H2O « CO2 + H2, 2CO « C + CO2).22 Anderson and co-
Introducing Lagrange multipliers for species i, subject to mass
workers studied the thermodynamics of the hydrogenation of
balance constraints:
carbon oxides for the production of hydrocarbons and
oxygenated organic molecules and found that it involves X
k
This journal is ß The Royal Society of Chemistry 2012 RSC Adv., 2012, 2, 2358–2368 | 2359
View Article Online
identical. At this point, the free energy of the system is at a the measurement within the temperature range 200–600 uC.
minimum. Firstly, 200 mg of each catalyst (20–40 mesh) was mixed with
Published on 30 January 2012. Downloaded on 9/17/2023 6:23:17 AM.
As a theoretical study, catalyst, reaction kinetics, and the 1.0 g quartz sand (20–40 mesh) and placed in the reactor. The
transport process are not considered here. However, the results thermocouple was directly inserted into the catalyst bed to
from the thermodynamic analysis could be used as a benchmark measure the actual reaction temperature. Prior to the reaction,
to compare with the experiment results. Here, possible reactions catalysts were reduced by H2 with a flow rate of 50 mL min21 at
involved in the methanation of carbon oxides are given in Table 1 450 uC for 2 h. Feed gas with a molar ratio of H2/CO/N2 = x/1/1
for the calculation.8,28 Our calculation was conducted based on (x = 1, 3, and 5) controlled with the mass flowmeters was
the gaseous compounds containing H2, O2, N2, CO, CO2, CH4, introduced into the reactor at a weight hourly space velocity
H2O, C2H4 (as a typical representative of high hydrocarbons), and (WHSV) of 30 000 mL h21 g21. The reaction temperature
solid carbon. Other oxygen-containing compounds such as was controlled by the heat controller with a heating rate of
methanol, methanoic acid, acetic acid and other high hydro- 5 uC min21. After one hour of steady-state operation at a set
carbons are not considered here because of their trace amounts in temperature, the outlet product gas stream was cooled down and
the equilibrium mixture gas. Carbon (C), hydrogen (H) and analysed online using a MicroGC (3000A, Agilent).
oxygen (O) atoms are used for the elemental mass balance. CO
and CO2 conversion, CH4 selectivity and yield, and carbon yield 3. Results and discussion
are defined in Table 2. It should be noted that the molar flow rate
of solid carbon used in the definition can be obtained through 3.1. The equilibrium constants of reactions involved in methanation
calculation by the Gibbs free energy minimization method in Fig. 1 shows the equilibrium constants (K) of the first eight
accordance with other gaseous compounds.27 possible reactions (see Table 1) occurring in the methanation
process as a function of temperature. K is calculated below by
2.2. Experimental procedure the van’t Hoff equation:
For a comparison with the thermodynamic calculation, CO d ln K H Dr Hm
H
methanation was carried out on two commercial nickel-based ~ 2
(5)
dT RT
catalysts marked as C1 (N112, NiO 67 wt%, JGC Catalyst and
Chemicals Ltd., Japan) and C2 (HT-1, NiO 57 wt%, LiaoNing
HaiTai SCI-TECH Development CO., LTD, China). Two fixed It can be seen that all exothermic reactions are suppressed as
bed reactors made of a quartz tube (at 1 atm) and stainless steel temperature is increased except for the endothermic methane
(at high pressures) both with a diameter of 10 mm were used in cracking reaction (R6). R1, R2, R3, R4, and R7 play a key role
Table 2 Definition of CO and CO2 conversion, CH4 selectivity and yield, and carbon yielda
Entry Definition
2360 | RSC Adv., 2012, 2, 2358–2368 This journal is ß The Royal Society of Chemistry 2012
View Article Online
This journal is ß The Royal Society of Chemistry 2012 RSC Adv., 2012, 2, 2358–2368 | 2361
View Article Online
3.2.3. Effect of H2/CO ratio. Since syngas has a variable ratio 3.2.4. Effect of H2O. Industrial experience showed that adding
of H2/CO, according to the stoichiometric ratio of R1, H2/CO is water steam into reactants could avoid carbon formation to a
required to be at around 3, which is normally controlled through large extent on the methanation catalysts.1 Fig. 5 shows the
water-gas shift reaction (R5).1 However, it is very hard to tune effect of adding steam into the feed gas at 1 and 30 atm. Fig. 5a
this value to be exactly 3 in the methanation industry. Thus, it is indicates that additional steam with different ratios leads to a
necessary to know the influence of this ratio on the methanation slight decrease in CO conversion at the high temperature (1 atm
process. Fig. 4a exhibits the effect of different H2/CO ratios at and 30 atm) since steam is one of the reaction products inhibiting
different pressures on CO conversion. It can be seen that the CO CO methanation (R1) to some extent according to Le Chatelier’s
conversion is not obviously changed at different H2/CO ratios (1, principle.
3, and 5). At 30 atm, a higher ratio seems to be beneficial Similarly, Fig. 5b shows that the introduction of steam results
towards CO conversion. Fig. 4b reveals the CH4 selectivity at in a small difference in CH4 selectivity, in which, small decreases
different H2/CO ratios. It can be seen that a higher H2/CO ratio with an increase of steam ratio can be found at 30 atm and at a
leads to a higher CH4 selectivity at both 1 and 30 atm. Reaction temperature more than 500 uC. Accordingly, Fig. 5c shows the
at high pressure also results in high CH4 selectivity at different effect of adding steam on the CH4 yield. The trend is similar to
2362 | RSC Adv., 2012, 2, 2358–2368 This journal is ß The Royal Society of Chemistry 2012
Published on 30 January 2012. Downloaded on 9/17/2023 6:23:17 AM. View Article Online
Fig. 5 Effect of water steam added into the feed gas on CO conversion Fig. 6 Effect of O2 contained in the feed gas on CO conversion (a), CH4
(a), CH4 selectivity (b), CH4 yield (c), and carbon yield (d). selectivity (b), CH4 yield (c), and carbon yield (d).
that in Fig. 5a and 5b. Fig. 5d reveals the effect of steam on the necessary to completely remove the O2 in feed gas for the
carbon yield at 1 atm. It can be seen that additional steam methanation process.
remarkably reduces the carbon deposition at 1 atm especially at
high temperatures although the temperature at which maximum 3.2.6. Effect of CH4. Syngas often contains a certain amount
carbon yield is obtained is not changed. When steam ratio of CH4.1,3,4 On the other hand, recycling the methanation gas
reaches to 0.5, there is no carbon formed. This is because adding product containing CH4 and H2O is usually employed to dilute
water would inhibit R7 and R8 significantly. Therefore, adding a the feed gas to control the reaction temperature in the industrial
small amount of water into the CO methanation system could methanation process such as the Lurgi process,33 TREMPTM
reduce carbon deposition. This is in agreement with the technology,18 and HICOM process.21 Thus, it is necessary to
HICOM21 and RMP32 methanation processes, in which excess understand the effect of CH4 on the methanation reaction as
steam is required for the first methanation reactor in order to shown in Fig. 7. Fig. 7a, 7b, and 7c show the effect of CH4 when
avoid carbon deposition. In addition, we did not do the it is introduced into the feed gas on the CO conversion, CH4
calculation at 30 atm because at this pressure, there is no carbon selectivity, and CH4 yield, respectively. It can be seen that
deposition even without adding steam regarding above Fig. 3d additional CH4 with different ratios does not lead to a large
and 4d. difference in conversion, selectivity and yield at both 1 and
30 atm. Accordingly, a high yield of CH4 still can be achieved at
3.2.5. Effect of O2. It is known that trace amount of O2 in the low temperature and high pressure even with high CH4 ratios
syngas is inevitable in the gasification process. At the present, (H2/CO/CH4 = 3/1/5) as shown in Fig. 7c. However, Fig. 7d
there is no report on the effect of O2 on methanation. Fig. 6 shows that deposited carbon yield sharply increases with a high
depicts the effect of O2 introduced into the feed gas on CO ratio of CH4 above 400 uC at both 1 and 30 atm mainly due to
methanation. Fig. 6a shows the O2 effect on CO conversion, the methane cracking reaction (R6), which occurs at high
from which, it can be seen that CO conversion decreases slightly temperature (Fig. 1). In brief, to get a high CH4 yield and a low
with an O2 increase at both 1 and 30 atm, owing to the reactions carbon yield, as well as optimized the energy efficiency, a large
of H2 and CO with O2, especially at high temperatures, recycling of the product gas is not suggested, and CH4 in the feed
respectively. Fig. 6b shows the effect of O2 on CH4 selectivity, gas should be controlled at a low level. The recycle ratio of
suggesting that O2 could lead to the obvious decline of CH4 product gas would be controlled in an appropriate range, for
selectivity in the whole temperature range because CO is partially example, 0.5–3.0 for the methanation process. This may explain
oxygenated to CO2. As a result, CH4 yield decreases with an O2 why the TREMPTM technology is dedicated to decreasing the
increase in reactants as shown in Fig. 6c. However, as shown in recycle ratio.18
Fig. 6d, the introduction of O2 can lead to a large reduction of
deposited carbon at 1 atm because of the reaction of O2 with 3.2.7. Effect of trace C2H4. It is also known that the presence
carbon to form CO2, especially at the ratio of H2/CO/O2 = 3/1/ of gaseous hydrocarbon compounds in syngas, such as C2H6,
0.5. Similarly, again at 30 atm, no carbon deposition occurs at C2H4, and C2H2, drastically raises the risk for carbon deposi-
high temperature. Therefore, considering the CH4 yield, it is tion.1 Fig. 8 shows the effect of C2H4 introduced into the feed
This journal is ß The Royal Society of Chemistry 2012 RSC Adv., 2012, 2, 2358–2368 | 2363
Published on 30 January 2012. Downloaded on 9/17/2023 6:23:17 AM. View Article Online
2364 | RSC Adv., 2012, 2, 2358–2368 This journal is ß The Royal Society of Chemistry 2012
Published on 30 January 2012. Downloaded on 9/17/2023 6:23:17 AM. View Article Online
Fig. 10 Effects of pressure and temperature on CO2 methanation: (a) Fig. 11 Effect of different H2/CO2 ratios on CO2 methanation: (a) CO2
CO2 conversion, (b) CH4 selectivity, and (c) CH4 yield. conversion, (b) CH4 selectivity, (c) CH4 yield, and (d) carbon yield.
which may explain why no significant carbon deposition is found in promote CO2 methanation, which is similar to that predicted by
CO2 methanation. Similarly, compared with CO methanation, CO2 the calculation for CO methanation (see Fig. 3c). Furthermore, it
methanation is also favourable in thermodynamics. In order to should be mentioned that no carbon generation is found under
obtain a high CH4 yield at 1 atm, the reaction temperature should our simulation conditions (not shown here). This is because
not exceed 300 uC. However, the eight-electron reduction of CO2 to during the CO2 methanation, water is formed that can suppress
CH4 by hydrogen is difficult to achieve because of significant the carbon deposition as mentioned above.
kinetic barriers.37 Thus, a much more active catalyst is needed to
hydrogenate CO2 to CH4 especially at the low temperature. 3.3.3. Effect of H2/CO2 ratio. Fig. 11 describes the effect of
different H2/CO2 ratios on CO2 methanation. It can be seen that
3.3.2. Effects of pressure and temperature. Fig. 10 shows the CO2 conversion in Fig. 11a and CH4 selectivity in Fig. 11b are
effects of pressure and temperature on CO2 methanation. remarkably affected by H2/CO2 ratio. High H2/CO2 ratio
Fig. 10a shows the effect on CO2 conversion. It can be seen generally leads to high CO2 conversion and CH4 selectivity at
that the CO2 conversion decreases with increasing temperature both 1 and 30 atm. When H2/CO2 is equal to 2, CO2 conversion of
and increases with the pressure, at temperatures below 600 uC, only about 50–70% can be obtained at either 1 atm or 30 atm, and
similar to the trend of CO methanation as shown in Fig. 3a. This maximum CH4 selectivity of 73% and 88% is achieved at 1 atm
is because CO2 methanation is also a volume reducing and and 30 atm, respectively. Fig. 11c shows the effect on CH4 yield.
exothermic reaction. When CO2 methanation takes place at When the H2/CO2 ratio is 2, CH4 yield is as low as about 40% at
1 atm and beyond 600 uC, the CO2 conversion gradually 1 atm and 45% at 30 atm, respectively, and CH4 yield can be
increases, different from what is seen in CO methanation as enhanced markedly by increasing H2/CO2 ratio at the same
shown in Fig. 3a. This is mainly because the reversed reaction of conditions. Fig. 11d shows the carbon yield. When H2/CO2 ratio is
R5 dominates above 600 uC and consumes CO2. Fig. 10b shows 2, abundant carbon deposition (up to 50%) is found below 500 uC
the effect on CH4 selectivity. It is observed that the trend with at both 1 and 30 atm. When H2/CO2 ratio is equal to or more than
variations in the pressure and temperature is similar to CO 4 (the curve at the ratio of 6 is not shown here), carbon deposition
methanation (see Fig. 3b). It should be noted that our is not found. Thus, in order to obtain a high CH4 yield and avoid
calculation result at 1atm is in accordance with Beuls’s25 and carbon deposition, H2/CO2 ratio should not be lower than 4 even
Ocampo’s26 simulations for CO2 methanation, and also in good at 30 atm. Conclusively, low temperature, high pressure, and
agreement with the experimental results of Du and co-workers proper H2/CO2 ratio is required for optimized CO2 methanation.
based on Ni/MCM-41 catalysts for CO2 methanation.38
Comparing Fig. 10a with Fig. 3a, we find that CO2 is more 3.3.4. Effect of H2O. Fig. 12 gives the effect of adding steam into
difficult than CO to be methanized at the same temperature and the feed gas (H2/CO2/H2O = 3/1/x, x = 0, 0.2, and 0.5, here x = 0
pressure although the CH4 selectivity below 500 uC in the CO2 means that no additional H2O is added). It is seen that at 1 and
methanation system is relatively higher than that in CO 30 atm, the addition of steam leads to a little decrease in CO2
methanation. Fig. 10c shows effect of temperature and pressure conversion (Fig. 12a), and no significant difference in CH4
on the CH4 yield. Low temperature and high pressure also selectivity (Fig. 12b) and CH4 yield (Fig. 12c). This is mainly
This journal is ß The Royal Society of Chemistry 2012 RSC Adv., 2012, 2, 2358–2368 | 2365
View Article Online
2366 | RSC Adv., 2012, 2, 2358–2368 This journal is ß The Royal Society of Chemistry 2012
View Article Online
possibly because the catalyst used here is not active enough, which
is a big challenge for methanation at low temperatures. However,
above 350 uC, the two catalysts show high activity for CO
conversion, which is highly consistent with the calculation values
at different pressures. Fig. 15b gives the CH4 selectivity on the two
catalysts. The experimental selectivity below 450 uC at 1 atm is
reasonably smaller than that of the calculation, but above 450 uC,
the selectivity is much higher. This may be because the actual
pressure in the reactor may have inevitably increased at 1 atm due
to the inherent system pressure drop as mentioned above. In
addition, the temperature gradient distribution also causes an
error in the catalyst bed temperature recorded, resulting in a
further gap between experimentation and calculation. Moreover,
other factors such as the reactants transport process may also have
Published on 30 January 2012. Downloaded on 9/17/2023 6:23:17 AM.
Table 3 Calculative results for syngas methanation together with the experiment data reported in the literature
Case 1 Case 2 Case 3
Entry
In Out Outcal In Out Outcal In Out Outcal
T/uC 300 604 604 398 729 729 325 590 590
P (atm) 27.2 27.1 27.1 30.0 30.0 30.0 30.0 30.0 30.0
H2 36.88 20.96 19.90 25.64 20.60 19.92 20.6 10.42 9.72
CO 4.28 1.17 1.10 18.61 8.40 8.03 8.40 2.21 2.02
CO2 6.13 4.46 4.28 14.74 23.30 23.58 23.30 27.74 27.84
CH4 28.12 37.44 38.02 0.048 4.94 5.29 4.94 10.20 10.48
H2O 19.19 29.82 30.59 40.23 41.96 42.38 41.96 48.56 49.06
N2(Ar) 5.41 6.15 6.10 0.72 0.79 0.80 0.79 0.86 0.87
This journal is ß The Royal Society of Chemistry 2012 RSC Adv., 2012, 2, 2358–2368 | 2367
View Article Online
(CO and/or CO2) using minimization of the Gibbs free energy 5 M. Chandel and E. Williams, Synthetic Natural Gas (SNG):
Technology, Environmental Implications, and Economics, Climate
method. The effects of temperature (200–800 uC), pressure (1– Change Policy Partnership, Duke University, 2009.
100 atm), ratio of H2/CO (and H2/CO2), and addition of water 6 A. H. Chen, T. Miyao, K. Higashiyama, H. Yamashita and M.
steam (and O2, CH4, and C2H4) on the methanation reactions were Watanabe, Angew. Chem., Int. Ed., 2010, 49, 9895–9898.
comprehensively investigated in terms of their effects on the 7 C. Galletti, S. Specchia, G. Saracco and V. Specchia, Chem. Eng. Sci.,
2010, 65, 590–596.
conversion of CO and CO2, CH4 selectivity and yield, as well as on 8 G. A. Mills and F. W. Steffgen, Catal. Rev. Sci. Eng., 1974, 8,
the deposition of carbon. The results indicate that a high CH4 159–210.
yield (> 99%) can be obtained from CO methanation at low 9 V. Ponec, Catal. Rev. Sci. Eng., 1978, 18, 151–171.
temperatures, high pressures, and high H2/CO ratio without 10 W. Wei and G. Jinlong, Front. Chem. Sci. Eng., 2010, 5, 2–10.
11 R. W. Dorner, D. R. Hardy, F. W. Williams and H. D. Willauer,
significant carbon deposition. Steam in the feed gas can remarkably Energy Environ. Sci., 2010, 3, 884–890.
suppress carbon formation. Impurity such as O2 and C2H4 should 12 J. Kopyscinski, T. J. Schildhauer and S. M. A. Biollaz, Chem. Eng.
be completely removed in the feed gas to maximise CH4 yield and Sci., 2011, 66, 924–934.
minimise carbon formation. A low recycle ratio of product gas 13 J. Kopyscinski, T. J. Schildhauer and S. M. A. Biollaz, Ind. Eng.
Chem. Res., 2011, 50, 2781–2790.
containing CH4 and H2O is needed to reduce carbon deposition. 14 V. Jimenez, P. Sanchez, P. Panagiotopoulou, J. L. Valverde and A.
Published on 30 January 2012. Downloaded on 9/17/2023 6:23:17 AM.
Compared to CO methanation, CO2 is relatively difficult to be Romero, Appl. Catal., A, 2010, 390, 35–44.
hydrogenated at the same reaction condition. The comparison of 15 M. M. Zyryanova, P. V. Snytnikov, Y. I. Amosov, S. A.
experimental data in the case studies with the calculations Ven’yaminov, E. Z. Golosman and V. A. Sobyanin, Kinet. Catal.,
2010, 51, 907–913.
demonstrated that the Gibbs free energy minimization method is 16 S. Eckle, H. G. Anfang and R. J. Behm, J. Phys. Chem. C, 2011, 115,
significantly effective for thermodynamic analysis of the methana- 1361–1367.
tion process. Our work can principally give a useful guidance in the 17 C. Bartholomew, Catal. Rev. Sci. Eng., 1982, 24, 67–112.
catalyst development and process control of methanation. 18 HaldorTopsoe, Topsøe Recycle Energy-efficient Methanation
Process, http://www.topsoe.com/business_areas/gasification_based/
Processes/Substitute_Natural_Gas.aspx, Accessed 08.02, 2011.
Nomenclature 19 GPGP, Practical experience gained during the first twenty years of
operation of the great plains gasification plant and implications for
aij the number of atoms of element j in species i future projects, Office of Fossil Energy, 2006.
20 Davy, Methanation Technology for Coal-to-SNG Project in China,
fi fugacity coefficient of species i http://www.davyprotech.com/pdfs/080110.pdf, Accessed 08.22.2011.
Fi,in molar flow rate for products i at inlet 21 R. L. Ensell and H. J. F. Stroud, The British gas HICOM
Fi,out molar flow rate for products i at outlet methanation process for SNG production, 1983.
GT total Gibbs free energy 22 M. Greyson, J. J. Demeter2, M. D. Schlesinger, G. E. Johnson, J.
Jonakin and J. W. Myers, Synthesis of Methane, D. O. Interior
ni mole of species i Report 5137, Bureau of Mines, U.S., 1955.
Ni the number of carbon atoms of species i 23 R. B. Anderson, C. B. Lee and J. C. Machiels, Can. J. Chem. Eng.,
PH standard-state pressure of 101.3 kPa 1976, 54, 590–594.
R molar gas constant 24 R. B. Anderson, J. Phys. Chem., 1986, 90, 4806–4810.
25 A. Beuls, C. Swalus, M. Jacquemin, G. Heyen, A. Karelovic and P.
mi chemical potential of species i Ruiz, Appl. Catal., B, 2011, DOI: 10.1016/j.apcatb.2011.02.033.
miH standard chemical potential of species i 26 F. Ocampo, B. Louis and A. C. Roger, Appl. Catal., A, 2009, 369, 90–96.
l Lagrange multiplier 27 M. K. Nikoo and N. A. S. Amin, Fuel Process. Technol., 2011, 92,
678–691.
28 G. A. Nahar and S. S. Madhani, Int. J. Hydrogen Energy, 2010, 35,
Acknowledgements 98–109.
29 S. Adhikari, S. Fernando and A. Haryanto, Energy Fuels, 2007, 21,
The authors gratefully acknowledge the supports from the 2306–2310.
Hundred Talents Program of the Chinese Academy of Sciences 30 S. Ozkara-Aydinoglu, Int. J. Hydrogen Energy, 2010, 35,
12821–12828.
(CAS), National Key Basic Research Program of China (No.
31 C. Rossi, C. G. Alonso, O. A. C. Antunes, R. Guirardello and L.
2011CB200906), National Key Technology R&D Program of Cardozo, Int. J. Hydrogen Energy, 2009, 34, 323–332.
China (No. 2010BAC66B01), State Key Laboratory of Multiphase 32 A. White G, R. Roszkowski T and W. Stanbridge D, in Methanation
Complex Systems of China (No. MPCS-2009-C-01), Knowledge of Synthesis Gas, AMERICAN CHEMICAL SOCIETY, 1975, vol.
146, ch. 8, 138–148.
Innovation Program of the CAS (No. KGCX2-YW-396), and
33 U. S. D. o. Energy, Practical experience gained during the first twenty
National Natural Science Foundation of China (No. 21031005). F. years of operation of the great plains gasification plant and implications
Gu is grateful to the support of K.C. Wang Post-doctoral for future projects, Office of Fossil Energy dakota gasification
Fellowships of the CAS and China Postdoctoral Science company prepared for US Department of Energy, 2006.
34 J. L. Li, N. Fu and G. X. Lu, Chin. J. Inorg. Chem., 2010, 26,
Foundation (No. 20100480026). 2175–2181.
35 T. Abe, M. Tanizawa, K. Watanabe and A. Taguchi, Energy Environ.
References Sci., 2009, 2, 315–321.
36 J. N. Park and E. W. McFarland, J. Catal., 2009, 266, 92–97.
1 J. Kopyscinski, T. J. Schildhauer and S. M. A. Biollaz, Fuel, 2010, 89, 37 K. R. Thampi, J. Kiwi and M. Gratzel, Nature, 1987, 327, 506–508.
1763–1783. 38 G. A. Du, S. Lim, Y. H. Yang, C. Wang, L. Pfefferle and G. L.
2 M. C. Seemann, T. J. Schildhauer and S. M. A. Biollaz, Ind. Eng. Haller, J. Catal., 2007, 249, 370–379.
Chem. Res., 2010, 49, 7034–7038. 39 H. Harms, B. Hohlein and A. Skov, Chem. Ing. Tech., 1980, 52, 504–515.
3 A. Haryanto, S. D. Fernando, L. O. Pordesimo and S. Adhikari, 40 C. WOODWARD, High-temperature methanation catalyst for SNG
Biomass Bioenergy, 2009, 33, 882–889. applicationshttp://www.anl.gov/PCS/acsfuel/preprint%20archive/
4 M. Jurascik, A. Sues and K. J. Ptasinski, Energy Environ. Sci., 2009, Files/21_4_SAN%20FRANCISCO_08-76_0022.pdf, Accessed 08.02,
2, 791–801. 2011.
2368 | RSC Adv., 2012, 2, 2358–2368 This journal is ß The Royal Society of Chemistry 2012