CO hydrogenation combined with water-gas-shift reaction for synthetic natural gas production: a thermodynamic and experimental study

The hydrogenation of CO to synthetic natural gas (SNG) needs a high molar ratio of H2/CO (usually large than 3.0 in industry), which consumes a large abundant of hydrogen. The reverse dry reforming reaction (RDR, 2H2 + 2CO ↔ CH4 + CO2), combining CO methanation with water-gas-shift reaction, can significantly decrease the H2/CO molar ratio to 1 for SNG production. A detailed thermodynamic analysis of RDR reaction was carried out based on the Gibbs free energy minimization method. The effect of temperature, pressure, H2/CO ratio and the addition of H2O, CH4, CO2, O2 and C2H4 into the feed gas on CO conversion, CH4 and CO2 selectivity, as well as CH4 and carbon yield, are discussed. Experimental results obtained on homemade impregnated Ni/Al2O3 catalyst are compared with the calculations. The results demonstrate that low temperature (200–500 °C), high pressure (1–5 MPa) and high H2/CO ratio (at least 1) promote CO conversion and CH4 selectivity and decrease carbon yield. Steam and CO2 in the feed gas decrease the CH4 selectivity and carbon yield, and enhance the CO2 content. Extra CH4 elevates the CH4 content in the products, but leads to more carbon formation at high temperatures. O2 significantly decreases the CH4 selectivity and C2H4 results in the generation of carbon.


List of symbols
Greek symbols a ik Number of atoms of the k element present in each molecule of species i l i Chemical potential of species i (J/mol) / i Fugacity coefficient of species i k k Lagrange multiplier 1 Introduction Natural gas is a highly efficient and clean fossil fuel due to its high calorific value, low sooting tendency and slag free products, leading to its increasing consumption year by year (Gao et al. 2015;Meng et al. 2015a;Rönsch et al. 2016). In 2014, the consumption of natural gas in China increased to 197.3 billion cubic meters, with a growth rate of 30.9% every year in the last decade (BP 2016). Recently, the consumption of natural gas has raised a serious concern regarding its depletion because of its limited reserves (Kopyscinski et al. 2010;Huo et al. 2013), in comparison, coal is considered as a much more abundant energy resource in many countries. The production of synthetic natural gas (SNG) from coal has been developed to be a potential route to circumvent the limited supply of natural gas, especially in China (Li et al. 2014a, b;Lu et al. 2014).
Among the coal-to-SNG production processes, SNG is produced through the four major steps, i.e., coal gasification, water-gas-shift (WGS) reaction (CO ? H 2 O $ H 2 ? CO 2 ), gas cleaning and CO methanation (3H 2 ? CO $ CH 4 ? H 2 O) (Shinde and Madras 2014;Wang et al. 2015). The CO methanation reaction is a key process for increasing SNG production (Meng et al. 2015b;Götz et al. 2016;Gao et al. 2016). If one mole of CO is converted to methane, three moles of H 2 are stoichiometrically required. However, the content of carbon in coal is usually more than 60 wt% (up to more than 90 wt% in bituminous coal) whereas that of hydrogen is \5 wt% (Martelli et al. 2011;Shen et al. 2016). The high content of carbon in coal results in low H 2 /CO molar ratios, usually less than one, of produced gas from coal gasification (Messerle et al. 2016). For example, the produced gas of the British Gas-Lurgi (BGL) coal gasification process is composed of 60%-70% CO, 27%-30% H 2 , 0%-7% CH 4 , 1%-4% CO 2 , and trace amounts of O 2 and light hydrocarbons (Yu and Wang 2010). To increase the H 2 /CO ratio, the WGS reaction should be well controlled. It is worthwhile to mention, in order to enhance the CO conversion and CH 4 yield during industry processes, an even higher H 2 /CO ratio is usually used. For instance, the H 2 /CO ratio of the Lurgi process for methanation was optimized at about 3.2, and that of the Topsøe Recycle Energy Efficient Methanation (TREMP) process reached about 3.5 (Kopyscinski et al. 2010). More amount of CO needs to be converted to produce H 2 by WGS reaction in order to get a high H 2 /CO ratio, which results in the high operating cost and energy consumption.
The reverse dry reforming (RDR) reaction (2CO ? 2H 2 $ CH 4 ? CO 2 ), which is the combination of CO methanation with WGS reaction, can be used to produce SNG. Recently, many studies have focused on the RDR reaction with the H 2 /CO ratio of one. Yan et al. (2013) found that the catalyst preparation methodologies significantly affected the activity and stability of Ni/SiO 2 catalysts. Jiang et al. (2013Jiang et al. ( , 2014 investigated the stepwise sulfidation and sulfidation temperature on the catalytic activity of MoO 3 /CeO 2 -Al 2 O 3 . It is because there are many advantages of this reaction. First, the feed gas has a low H 2 / CO ratio of one, which needs less hydrogen; Secondly, the water-free products can diminish the damage of steam on methanation catalyst, and the by-product CO 2 can be easily removed by employing low-temperature methanol purification process. In other words, less H 2 is needed from the gas of coal gasification, which simplifies the SNG production process and reduces the cost. In literatures, data are available for the thermodynamic analysis of methanation reactions. Miguel et al. (2015) and Sahebdelfar et al. (2015) conducted a thermodynamic calculation of CO 2 methanation based on the method of Gibbs free energy minimization and compared with the experimental data. Gao et al. (2012) analyzed the thermodynamic properties of several reactions during the complete methanation of CO and CO 2 . However, these thermodynamic studies were carried out based on the complete methanation reactions. To our knowledge, there is little information on the thermodynamic analysis of the RDR reaction that occurs at low temperatures. Therefore, it needs to perform the calculations based on the Gibbs free energy minimization method and validate the data through experimental means.
It is well known that the produced gas from coal contains many impurities, such as steam, CO 2 , CH 4 , O 2 and light hydrocarbons of C 2 H 4 and C 2 H 6 . In order to increase the production of SNG and optimize the H 2 /CO ratio of the produced gas, effects of these substances on the catalytic performance of the RDR reaction have to be investigated. Moreover, the yield of solid carbon should be taken into account during the thermodynamic analysis.
The objective of this work is to elucidate, through a thermodynamic study supported by experimental data, the effects of temperature, pressure and the other factors affecting the RDR reaction, such as the H 2 /CO molar ratio and addition of H 2 O, CH 4 , CO 2 , O 2 , and C 2 H 4 in feed gas on the catalytic activity and selectivity and the yield. For this purpose, this study does not take into account of reaction kinetics, practical heat and mass transfer processes. It is expected to produce necessary thermochemical data to describe the effectiveness of the RDR reaction and to provide useful guidance to chemical engineers for optimizing the individual processes.

Thermodynamic analysis software
The HSC Chemistry software 6.0 allows simulating chemical reactions and processing on the thermochemical data basis. In this study, the modules of reaction equations and equilibrium compositions were utilized to calculate the effects of various substances in conversion, selectivity and yield. The calculations were performed based on an extensive thermochemical database, which contains enthalpy (H), entropy (S) and heat capacity (C p ) data of more than 17000 chemical compounds (Roine 2010; Kumar et al. 2016).

Thermodynamic analysis method
The equilibrium products at different temperatures and pressures were calculated using the Gibbs free energy minimization method, which has been widely applied for thermodynamic calculations (Adhikari et al. 2007;Nahar and Madhani 2010;López Ortiz et al. 2015). The detailed interpretation of this theory can be referenced by Wang et al. (Wang and Cao 2012;Wang et al. 2014).
The total Gibbs free energy is expressed as Eq.
(1). As to chemical reaction equilibrium state, f i H = p H , f i = y i u i p, and G i H = DG fi H are supposed. With the Lagrange multiplier method, Eqs. (2) and (3) indicate the minimum Gibbs free energy of each gas and total system without solid ones, respectively. Equation (4) is the constraining condition. The vapor-solid phase equilibrium is applied to the Gibbs free energy of carbon, as shown in Eq. (5). Thus, the minimization formation of Gibbs free energy Eq. (6) is obtained by substituting Eq. (1) with Eqs. (3) and (5).
In the HSC Chemistry software 6.0, the reaction system needs to be specified, in terms of its phases and species, and the amount of the reactants. The program calculates the amount of products at equilibrium in isothermal or isobaric condition for a heterogeneous system. At the equilibrium state, the free energy of the system is minimized.
It should be noted that these thermodynamic analyses do not include any reaction kinetic limitation or transport process in the real process. However, thermochemical calculations show a great importance in adjusting the feasibility of a reactive process under certain conditions. Here, possible reactions are summarized in Table 1 for the calculations which were carried out based on different types of gases including CO, H 2 , CO 2 and CH 4 , and the solid product of deposition carbon (graphite). Other substances such as alcohols, acids and high hydrocarbons are not taken into account due to their trace contents in the equilibrium  1 Equilibrium constants of the reactions as a function of temperature CO hydrogenation combined with water-gas-shift reaction for synthetic natural gas production… 441 gas mixture. The elemental mass balance is evaluated by carbon, hydrogen, and oxygen. The conversion of CO, selectivities of CH 4 and CO 2 , yields of CH 4 and solid carbon are calculated as follows: Here, i indicates all carbon containing species (CO, CO 2 , CH 4 and C 2 H 4 ) at inlet, and N i indicates the number of carbon atom of i-th species.

Experimental study
The alumina (191 m 2 /g, Shandong Aluminum Co., China) supported homemade Ni-based catalyst was prepared by the co-impregnation method, as described in Meng's works (Meng et al. 2017). The Ni-based catalyst, with the Ni loading of 20 wt% and La loading of 4 wt%, showed the specific surface area of 128 m 2 /g and pore size of 5.1 nm, and the catalyst was denoted as ExCat. To validate the thermodynamic calculations, the RDR reaction was carried out in a stainless steel, high-pressure fixed-bed tube reactor (10 mm 9 2 mm 9 500 mm) within the temperature  Fig. 3 Effect of temperature and pressure on catalytic performance. a CO conversion, b CH 4 selectivity, c CH 4 yield, and d carbon yield range of 300-550°C. 300 mg of Ni/Al 2 O 3 catalyst (20-40 mesh) was placed in the reactor. Prior to the RDR reaction, the catalyst was reduced at 550°C in a H 2 (99.99%, purchased from Taiyuan Iron & Steel (Group) Co., Ltd., China) flow diluted with 25% N 2 (99.995%, purchased from Taiyuan Iron & Steel (Group) Co., Ltd., China) for 6 h. A mixed feed gases of H 2 /CO = 1 (the gas of CO with a purity of 99.9% was purchased from Taiyuan Iron & Steel (Group) Co., Ltd., China) were introduced and controlled with the mass flow controller (MFC), preheat treatment was finished at 200°C in first oven at a space velocity of 20000 mL/(g h) -1 . In the second oven, two thermocouples are employed for the reaction. One is placed closely to the reactor, in the middle of the oven to control the oven temperature. The other one is placed inside of the catalyst bed for the measurement of reaction temperature of catalyst bed. The outlet gas steam was cooled by condenser (2°C) and quantitatively analyzed by an online gas chromatography (GC, Agilent 7890A) using helium (99.999%, purchased from Taiyuan Iron & Steel (Group) Co., Ltd., China) as the carrier gas. The GC equipped with a flame ionization detector (FID) with an HP-AL/S column was employed to analyze CH 4 , and a thermal conductivity detector (TCD) equipped with a Porapak-Q column, HP-PLOT/Q column, and HP-MOLESIEVE column was employed to analyze CO 2 , CO, and N 2.

Equilibrium analysis of the reactions
The equilibrium constants K of R1-R8 at various temperatures are shown in Fig. 1. The value of K is calculated using the Van't Hoff equation: It can be seen in Fig. 1, as the temperature increases, all the K values decrease except that of R6, which agrees with the Le Chatelier's principle. R1, R2, R3, R5, and R7 play important roles in the RDR reaction system. When the temperature is lower than 500°C, the equilibrium constant K reduces in the order of R1 Among all these reactions, R1 and R2 show elative high K values at low temperatures, which will lead to the high conversions of CO. CO 2 could be converted via reactions of R3 and R8; however, the CO 2 cannot be fully converted, which is due to that the reactions c Fig. 4 Effect of H 2 /CO molar ratio on catalytic performance. a CO conversion, b CH 4 selectivity, c CH 4 yield, and d carbon yield CO hydrogenation combined with water-gas-shift reaction for synthetic natural gas production… 443 of R1, R4, and R5 generate CO 2 . Moreover, the solid carbon generated from the reaction of R5 to R8, and the Boudouard reaction (R5) acts a dominant role due to its largest K value. Importantly, all these reactions may occur simultaneously in the system, resulting in a balanced composition of the products. Figure 2 shows the methanation products with their mole fractions at equilibrium temperatures at 0.1 MPa. The feed gas contains H 2 and CO with a H 2 /CO stoichiometric ratio of 1. The products mainly consist of CH 4 and CO 2 in the temperature range of 200-400°C. The mole fractions of CH 4 and CO 2 decrease as the temperature increases, whereas the mole fractions of H 2 and CO exhibit the opposite trend. This can be explained by that the methanegenerating reactions (R1-R3) are exothermic reactions and a higher temperature inhibits them. It is also found that the mole fraction of CH 4 is lower than that of CO 2 . At a low H 2 /CO ratio, the Boudouard reaction (R5) is more preferable, leading to a large amount of CO 2 and solid carbon.

Equilibrium compositions
On the other side, the amount of H 2 is higher than that of CO, as shown in Fig. 2. when the temperature over 550°C, the CO mole fraction increases linearly with temperature, which is due to the inhibition of the RDR reaction and the CO produced by the reverse Boudouard reaction (Nahar and Madhani 2010). The solid carbon is produced between 300 and 800°C, with a maximum amount around 575°C. Boudouard reaction will not occur when the temperature reaches 700°C or higher, as shown in Fig. 1, from which the solid carbon is mainly derived from the methane cracking reaction R6. Once the solid carbon produced, it may cover the active sites and results in catalyst deactivation (Takenaka et al. 2008). Thus to emphasize, a proper temperature below 600°C for RDR reaction is focused, in looking for a high yield of methane.

Effect of temperature and pressure
The effects of temperature and pressure on RDR reaction are shown in Fig. 3. Figure 3a shows that the CO can be fully converted when the temperature was below 400°C, further increasing the temperature resulted in the decrease of CO conversion, from the reason of the exothermic nature of RDR reaction. At a constant temperature, a higher pressure promotes the CO conversion due to the volume reducing nature of RDR reaction. The above results c Fig. 5 Effect of H 2 O added in the feed gas on catalytic performance. a CO conversion, b CH 4 selectivity, c CO 2 selectivity, and d carbon yield indicate that a lower temperature and a higher pressure are favorable for the RDR reaction, in terms of increasing CO conversion. On the other hand, the increasing range of CO conversion is not obvious when the pressure is higher than 1 MPa. In Fig. 3b, high CH 4 selectivity is obtained at low temperatures and high pressures. This is because all the methane producing reactions are volume reducing and exothermic. When the temperature is higher than 550°C, the K value of R5 is comparatively larger than that of R1 (Fig. 1), and the Boudouard reaction (R5) becomes dominantly in the reaction system, leading to a high CH 4 selectivity at relative low temperatures. In Fig. 3c, a CH 4 yield close to 50% is obtained in the temperature ranges of 200-300°C from 1 to 5 MPa. Thus, to get a comparatively high CO conversion and CH 4 yield, the conditions of high pressure and low temperature are recommended. Knowing that low temperature is not benefit to accelerate the reaction rate and high pressure is harmful for the equipment, a pressure range of 2-3 MPa and a temperature range of 300-500°C are favorable for the RDR reaction. The variation of carbon yield is presented in Fig. 3d. All these carbon yield curves exhibit a volcano characteristic, with less yield of carbon at high pressures. The solid carbon results from many reactions, including R5, R6, R7, and R8 (as shown in Table 1), from which have different K values. Since the K value of R6 is negative at 200-550°C and the value of R5 is higher than that of R7 and R8 at 200-800°C (Fig. 1), so R5 is the main reason for the deposition of carbon. At the point of 0.1 MPa and around 575°C, the carbon yield reaches the maximum (23%). Accordingly, at this condition, the occurrence of R6 triggered a higher production of carbon. However, further increase the temperature results in the decrease of carbon yield, possibly because the reverse reactions of R5, R7 and R8 consumes a comparable amount of solid carbon.

Effect of H 2 /CO ratio
Since the produced gas derived from coal gasification has a varying ratio of H 2 /CO (Zheng and Furinsky 2005), there is a need to investigate the effect of H 2 /CO ratio on performance of RDR reaction. Figure 4a exhibits the effect of H 2 /CO molar ratio on CO conversion. The CO conversion changes slightly as the H 2 /CO ratio increases. When the pressure is 3 MPa, the CO conversion reaches nearly 100% in the temperature range of 200-550°C. Figure 4b shows the changes of CH 4 selectivity. A high H 2 /CO ratio improves the selectivity of CH 4 , as it can be found when c Fig. 6 Effect of CH 4 imputed into the feed gas on catalytic performance. a CO conversion, b CH 4 selectivity, c CO 2 selectivity, and d carbon yield the H 2 /CO ratio is 0.8, the maximum CH 4 selectivity of 40% can be obtained. Considering the nearly 100% CO conversion below 450°C, there should be much carbon deposition (see Fig. 4d). When the H 2 /CO ratio increased to 1, the corresponding initial CH 4 selectivity increases to 50%. Further increase H 2 /CO ratio to 3, the CH 4 selectivity enhanced remarkably to 100%. In addition, high pressure is favorable to improve the CH 4 selectivity. So a high H 2 /CO ratio or pressure value is useful for SNG production. Figure 4c reveals the variations of CH 4 yield. Higher H 2 /CO ratio or pressure and lower temperature lead to a higher CH 4 yield. Figure 4d demonstrates the variation of carbon yield. When the H 2 /CO ratio is 0.8, a significant amount of carbon is produced, because a large amount of unreacted CO is converted to solid carbon via R5. The highest carbon yield is 28 mol% at 0.1 MPa. To avoid such carbon deposition, high H 2 /CO ratio and high pressure are recommended for the RDR reaction.

Effect of H 2 O content
Steam controls the H 2 /CO ratio via WGS reaction (R4), which is mostly used in methanation and ammonia synthesis industrial process. Moreover, it can be also used for eliminating the carbon deposition to some extent via reverse R7 and R8.
The effect of steam content in feed gas is shown in Fig. 5. As can be found in Fig. 5a, the introduction of steam slightly decreases the CO conversion at 0.1 and 3 MPa. Although the steam does not participate in the RDR reaction, it can inhibit the methanation reaction R2 and thus decrease the CO conversion. As shown in Fig. 5b, steam has a small effect on the selectivity of CH 4 . However, the selectivity of CO 2 significantly increases as the amount of steam increases at 0.1 MPa (Fig. 5c). From  Fig. 5d, the additional steam decreases the yield of carbon, especially at high pressures. When the steam ratio reaches 0.4 at 3 MPa, only trace amount of carbon is formed at the temperature higher than 650°C. This is because the added steam promotes the reverse R8 reaction. This is also verified by Fig. 5c, where the addition of steam results in the increase of CO 2 selectivity. In industrial methanation processes, such as High Combined Shift Methanation (HICOM) (Ensell and Stroud 1983) and Ralph M. Parsons (RMP) (G. A White et al. 1975), certain amount of steam is required to eliminate the solid carbon. c Fig. 7 Effect of CO 2 introduced into the feed gas on catalytic performance. a CO conversion, b CH 4 selectivity, c CO 2 selectivity, and d carbon yield

Effect of CH 4 content
The syngas, derived from gasification of coal, frequently contains a certain amount of CH 4 . In addition, the contained CH 4 in methanation products is generally cycled to dilute the feed gas, in order to avoid reaction temperature run away from the fixed-bed technology (Rönsch et al. 2016). Thus, to study the effect of CH 4 on the RDR reaction is of necessity. Figure 6a shows the effect of CH 4 content on CO conversion. Additional CH 4 results in slight decrease of CO conversion at both 0.1 and 3 MPa. The reason is that CH 4 is the product of reactions R1, R2 and R3. Introduction of CH 4 can inhibit these reactions (R1, R2 and R3), and thus decrease the CO conversion. In Fig. 6b and c, additional CH 4 greatly improves the CH 4 content in the product and reduces the CO 2 selectivity at both 0.1 and 3 MPa, especially at low temperatures. However, the introduction of CH 4 into the feed gas sharply exacerbates the deposition of solid carbon, especially at high temperatures (Fig. 6d). Increasing pressure is very effective to reduce the carbon deposition. Comparing to 0.1 MPa, the carbon yield has been greatly reduced at 3 MPa. Interestingly, increasing the CH 4 ratio from 1 to 2 increases slightly the carbon yield, when the temperature was lower than 600°C. Therefore, to enhance the CH 4 yield and avoid the deposition of carbon, a certain content of CH 4 in feed gas or product gas for recycling should not be ignored and the temperature should not exceed 600°C.

Effect of CO 2 content
Usually, the syngas needs to be purified to remove the acidic gas, such as CO 2 . In this section, the effect of CO 2 on the RDR reaction is elucidated. Figure 7 shows the effect of CO 2 content on the performance. In Fig. 7a, the increasing of CO 2 ratio results in a slightly decrease of CO conversion, because the addition of CO 2 inhibits the RDR reaction and probably accelerates the reaction rate of CO 2 methanation (R3), based on the Le Chatelier's principle. Figure 7b, c reveal the effect of CO 2 on the selectivities of CH 4 and CO 2 , respectively. As the CO 2 amount increases, the selectivity of CH 4 decreases gradually at 0.1 and 3 MPa. However, the selectivity of CO 2 exhibits the opposite trend, especially when the temperature is higher than 600°C at 0.1 MPa, the selectivity of CO 2 increases significantly under the high CO 2 amount. Meanwhile, the carbon yield (see Fig. 7d) drops gradually with an increasing of CO 2 ratio, due to the inhibition of Boudouard c Fig. 8 Effect of O 2 imported into the feed gas on catalytic performance. a CO conversion, b CH 4 selectivity, c CO 2 selectivity, and d carbon yield reaction (R5). To decrease the carbon yield, the addition of CO 2 could work at both pressures. In brief, high pressure inhibits carbon formation, however, it is better to remove the CO 2 in the syngas to get a high methane yield.

Effect of O 2 content
In addition to steam, air and oxygen are also used as the feed gases for the coal gasification to produce the syngas. Thus, the effect of O 2 amount on the RDR reaction is necessary to be studied. The results are shown in Fig. 8. Figure 8a depicts the effect of O 2 amount on CO conversion. The CO conversion exhibits almost no change at 200-500°C at a constant pressure. However, when the temperature is higher than 500°C, as the O 2 amount increases, the CO conversion decreases slightly at 3.0 MPa. It is probably due to that a small amount of O 2 reacts with the CH 4 (partial oxidation of methane, 1/2O 2 ? CH 4 $ CO ? 2H 2 ), which changes the CO conversion. Figure 8b shows the effect of O 2 amount on CH 4 selectivity. In the range of 200-500°C, as the O 2 content rises, the selectivity of CH 4 decreases especially at a high pressure. It is because that, when O 2 is introduced into the reaction system, it reacts with CO to generate CO 2 , resulting in the decrease of CH 4 selectivity. Another reason could be the reaction of partial oxidation of methane that consumes CH 4 , which also decreases the CH 4 selectivity. From Fig. 8c, the CO 2 selectivity increases as the O 2 content rises, especially when the temperature is higher than 500°C at 0.1 MPa. It is due to that the CO and the produced solid carbon are oxygenated by O 2 , especially at high temperatures. Both of the possibilities result in the increase of CO 2 selectivity. In Fig. 8d, the addition of O 2 reduces the carbon yield to a large extent at 0.1 and 3 MPa, due to the fact that solid carbon reacts with the O 2 . It is noticeable that no carbon is formed at 3 MPa when the O 2 ratio is 0.1 or 0.2. Therefore, to obtain a high CH 4 yield, from the thermodynamic point of view, the feed gas should contain none of O 2 .

Effect of C 2 H 4 content
In the process of coal gasification, some trace amounts of high hydrocarbons are produced, such as C 2 H 6 and C 2 H 4 . Here, C 2 H 4 is taken as a typical hydrocarbon to study the effect on performance. From Fig. 9a, b, C 2 H 4 slightly affects the CO conversion and CH 4 selectivity. This is mainly due to that, C 2 H 4 is not involved in the reactions c Fig. 9 Effect of C 2 H 4 brought into the feed gas on catalytic performance. a CO conversion, b CH 4 selectivity, c CO 2 selectivity, and d carbon yield (R1, R2, R4 and R5), and in which CO is a reactant. This observation agrees with Gao's results (Gao et al. 2012). However, as shown in Fig. 9c, the CO 2 selectivity decreases as the C 2 H 4 ratios increases. Figure 9d shows that the carbon yield rises rapidly with increasing the C 2 H 4 ratio, especially at 0.1 MPa. The temperature corresponding to the maximum carbon yield is around 600°C. It can be concluded that at low pressure (0.1 MPa) and high temperatures (about 600°C), the reaction of C 2 H 4 cracking (C 2 H 4 $ C ? 2H 2 ) should not be ignored. Whereas at high pressure (3 MPa), the carbon yield is remarkably reduced. Therefore, as an impurity, the C 2 H 4 content must be controlled to prevent the formation of solid carbon.

Comparison between thermodynamic calculations and experimental results
The RDR activity test was carried out on ExCat to compare the experimental results with thermodynamic calculations.
Both the results at various temperatures and pressures are shown in Fig. 10. Figure 10a shows the comparisons of CO conversion. Thermodynamically, the equilibrium conversion of CO is almost 100% between 200 and 400°C, further increasing the temperature decreases the CO conversion, especially at a low pressure (0.1 MPa). Experimental results show that the CO conversion firstly increases and then decreases as the temperature rises and the highest CO conversion was obtained around 400°C. The high pressure is benefit to improve the CO conversion. It should be noted that the experimental conversion of CO at 400°C or below is much lower than that of the calculated values, which may be due to the low reaction rate at low temperatures. When the temperature is 400°C or higher, the experimental results show good accordance with the calculated ones. Figure 10b shows the CH 4 selectivity at various temperatures and pressures. The calculated results shows that increasing temperature decreases the CH 4 selectivity, whereas the increasing pressure enhances the CH 4 selectivity. At low temperature and pressure, the experimental results are slightly higher than the calculation ones, while at high temperature, the experimental CH 4 selectivity is much higher than the calculated one. Figure 10c shows the comparison of CO 2 selectivity. The calculation results exhibit that the selectivity of CO 2 is constant at various temperatures and pressures, and the experimental results c Fig. 10 Comparison between experimental results and calculated ones at various temperatures and pressures. a CO conversion, b CH 4 selectivity, c CO 2 selectivity, d Carbon yield are slightly lower than the calculation ones. The above comparison show that more amount of CO 2 converted to CH 4 during the reaction. Figure 10d shows the comparison of carbon yields. The calculated carbon yields rose as the temperature increased, and decreased as the pressure increased, which means a large amount of CO is converted to carbon. However, the experimental results show that the yields of carbon were nearly zero, probably due to the catalyst inhibit the formation of carbon. The above discussion show that the experimental results are generally in accordance with the calculated ones at different temperatures and pressures. The result also indicates that the Gibbs free energy minimization method is an ideal tool for thermodynamic analysis of the RDR process.

Conclusions
A detailed thermodynamic equilibrium analysis of reverse dry reforming (RDR) reaction by minimizing the Gibbs free energy method in the range of 200-800°C and 0.1-3 MPa, and an experimental results in the range of 300-550°C and 0.1-3 MPa are studied. The calculation results demonstrate that low temperature and high pressure are beneficial for the CO conversion and CH 4 yield, and high H 2 /CO ratio (at least 1) promotes CH 4 yield and decreases carbon yield. In the range of 200-500°C and 1-5 MPa, the CO conversion and CH 4 yield reach 95%-100% and 43%-50%, respectively. Steam in the feed gas enhances the CO 2 selectivity and inhibits the generation of carbon, almost no carbon formed at the H 2 /CO/H 2 O ratio of 1/1/0.4, when the temperature is below 600°C at 3 MPa. CH 4 contained in the recycling product gas elevates the CH 4 content in the products, but also leads to more solid carbon at 500-800°C, especially at 0.1 MPa. CO 2 has a negative effect on CH 4 selectivity, but it could result in a slightly decrease of carbon yield at the temperature higher than 500°C. O 2 is not preferable for increasing CH 4 selectivity and decreasing the CO 2 selectivity although it decreases the carbon yield. C 2 H 4 is prone to crack, creating a high carbon yield. As impurities, O 2 and C 2 H 4 should be completely removed to get a high CH 4 yield. The experimental data are consistent with the calculation ones, indicating that minimizing the Gibbs free energy is effective to analyze the RDR reaction thermodynamically. This work is expected to provide a valuable suggestion in the process optimization for SNG production by combining CO methanation with WGS reaction.