催化学报  2015, Vol. 36 Issue (2): 153-159   PDF (25414 KB)    
扩展功能
加入收藏夹
复制引文信息
加入引用管理器
Email Alert
RSS
本文作者相关文章
樊铖
张强
魏飞
Controllable oxidation for oil recovery: Low temperature oxidative decomposition of heavy oil on a MnO2 catalyst
Cheng Fan, Qiang Zhang , Fei Wei    
Beijing Key Laboratory of Green Chemical Reaction Engineering and Technology, Department of Chemical Engineering, Tsinghua University, Beijing 100084, China
Abstract: Heavy oil is a readily available alternative energy resource with a reserve that is more than twice that of conventional light oil. In situ combustion is one of the most promising strategies for heavy oil exploitation, and the modulating of the oxidation behavior of heavy oil is an efficient way to expand the applicability of the in situ combustion method. MnO2 nanoparticles were employed to facilitate the cracking of heavy compounds, promote heat production, and improve recovery efficiency. The oxidative decomposition rate of heavy oil was doubled in the low temperature interval, and the heat release rate was accelerated in the high temperature interval. The increased weight loss at low temperature was attributed to the decomposition of heavy components. The detection of incomplete oxidation products by mass spectroscopy under excessive oxygen flow at high temperature indicated a diffusion controlled process of oil combustion. The same amount of CO2 from the combustion of less fuel demonstrated an increased oxidation degree of the products. The apparent activation energies of the oxidation reactions were decreased by 10-30 kJ/mol at low temperature and 20-40 kJ/mol at high temperature by the addition of MnO2. MnO2 can render in situ combustion more feasible for various oil reservoirs, and is also promising for other thermal recovery processes for improved oil recovery.
© 2015, Dalian Institute of Chemical Physics, Chinese Academy of Sciences.
Published by Elsevier B.V. All rights reserved.
Key words: Manganese dioxide     In situ combustion     Oxidative decomposition     Fuel deposition    
稠油的可控氧化热采:MnO2催化剂上的重油低温氧化分解行为
樊铖, 张强 , 魏飞    
清华大学化学工程系, 绿色反应工程与工艺北京市重点实验室, 北京100084
摘要:采用MnO2纳米催化剂促进稠油重组分的低温氧化裂解及高温氧化生热, 进而促进火烧油层技术在超稠油油藏开发中的应用, 提高火烧油层技术的采收率. 结果表明, MnO2促进了稠油的低温裂解, 高温生热量并不随燃料沉积量的减少而降低. MnO2的加入促进了低温下小分子的生成, 并使高温氧化反应更加温和, 提高了高温产物的氧化度. 相比无催化剂的氧化反应, MnO2使稠油的低温氧化和高温氧化的表观活化能分别降低了10-30 kJ/mol和20-40 kJ/mol.
关键词二氧化锰     火烧油层     氧化裂解     燃料沉积    
1. Introduction

Petroleum is the most important resource for fuels [1, 2] and chemicals [3, 4, 5]. The reliance on petroleum is still growing [6] because of the absence of a good alternative substance. Consequently, the exploitation of heavy oil, whose reserve is over twice that of conventional oil, has received much attention in recent years. Thermal recovery methods have been widely employed to reduce its viscosity and to simplify its transportation. Among the various strategies, in situ combustion, which is based on the exothermic oxidative reactions of heavy oil, is a promising method [6].

The control of an in situ combustion process is a big challenge because of the different geological characteristics as well as complex properties of heavy oil [7, 8, 9]. Neither the recovery efficiency nor product quality can be accurately predicted. As a result, the broad application of in situ combustion is severely hindered. Much effort has been devoted to the better control of in situ combustion. Moore et al. [10] and Weissman et al. [11] first proposed the concept of a combination of in situ combustion with down-hole catalytic upgrading in the 1990s. The addition of a Ni/Mo/Al catalyst by Moore et al. [10] in combustion tube experiments induced significant hydrodesulphurization and a gravity increase from 15.3° to 21.2° API. Catalytic upgrading using a conventional hydroprocessing catalyst combined with a hydrogen source under similar conditions was investigated by Weissman et al. [11]. A 50% sulfur removal and an 8-point API gravity increase were obtained. Together with recent advances in controllable horizontal drilling, the catalytic upgrading process in situ (CAPRI) integrated with toe-to-heel air injection (THAI), collectively called THAI-CAPRI [12, 13, 14, 15], was proposed, and a field trial has been implemented [16]. Nassar and co-workers [17] reported that water-in-vacuum gas oil microemulsions containing trimetallic (W, Ni, and Mo) ultra-dispersed colloidal nanoparticles penetrated into the porous medium, reacted with the bitumen and enhanced oil recovery. In addition to the investigations on improving product quality, the extension of its application to reservoirs with various APIs has also attracted much attention since the API of the oil affects the performance of this technology. In light reservoirs [18], it is difficult to achieve the desired high temperature region to ensure a self-sustaining process because not enough fuel is avai lable for combustion, while in heavy reservoirs, excessive fuel is deposited and combusted resulting in uneconomical exploitation [19]. Therefore, if heat production by combustion was enhanced or the cracking of heavy compounds at low temperature was promoted, fuel deposition for combustion will rely less on the property of the oil. Then, the in situ combustion can be controlled to be within an appropriate range, and is suitable for a wider variety of oil reservoirs.

Solid catalysts with extraordinary performance and multi-functions [20, 21, 22] have been employed to modulate the oxidation behavior of heavy oil [23, 24]. Castanier’s group [25, 26, 27, 28, 29] proposed water soluble metallic additives as catalyst to modify the fuel deposition reactions and to increase fuel laydown for light oils. This has been successfully applied in fields in Hungary [30] and Niemangu [31] in China. However, suitable additives for reducing fuel deposition have not been reported until very recently when Babadagli and coworkers [32] found that nickel ions had a significant effect on the low temperature reactions and increased the amount of mass loss by 10% at 300 °C, whereas this effect was not observed at 200 °C. A shift of the LTO products from oxygenated compounds to CO2 and H2O was also detected. Unfortunately, the catalytic effect on the kinetics and products and the heat-release performance of nickel ions was not fully illustrated in the work. Since heat is the most important factor for the thermal recovery method, the evaluation of heat release is critical, especially when the amount of fuel deposition is decreased. An additive is applicable only if it has no negative impact on the heat generation of oil combustion at high temperature. With a catalyst that promotes the cracking of heavy compounds at low temperature and enhances the heat production of the residue at high temperature, the in situ combustion process would give higher recovery efficiency and have wide applications in various oil reservoirs.

In this contribution, a metal oxide catalyst, MnO2, was used to facilitate the cracking of heavy compounds at low temperature and to enhance the heat-release rate at high temperature. In addition to the weight loss performance, special attention was also paid to the change in thermal effect and representative product concentrations due to the addition of the catalyst. The reason we selected MnO2 as the metal oxide catalyst was because of its excellent catalytic oxidation performance at low temperature [33, 34, 35, 36]. The MnO2 nanoparticle catalyst was characterized by X-ray diffraction (XRD), transmission electron microscopy (TEM), and X-ray photoelectron spectroscopy (XPS). The distributed activation energy model (DAEM) was employed to describe the kinetic behavior of the catalytic oxidation of heavy oil. This method is expected to modulate the oxidation behavior, regulate the oxidation characteristics of a crude oil, improve the quality of the produced oil, and enhance oil recovery.

2. Experimental
2.1. Feedstock

The heavy oil used was obtained from Xinjiang Oil Field. The properties and elemental analysis of the heavy oil are listed in Table 1. The crude oil was separated from the sand to eliminate catalytic effects arising from inorganic compounds. It was subsequently devolatilized under an inert atmosphere at 300 °C for 3.0 h to eliminate the influence of light component volatilization on weight loss and the heat flow profiles prior to oxidation experiments.

Table 1
Physical properties and elemental analysis of the heavy oil from Xinjiang Oil Field.

The MnO2 catalyst was obtained from Alfa Aesar. The mixture of MnO2 and heavy oil was obtained by the following steps. First, MnO2 nanoparticles were introduced into the heavy oil under vigorously mechanical stirring in an oil bath held at 80 °C to give an amount of 0.1 wt%. Second, the mixture was continuously stirred for 1 h, and then further treated in an ultrasonic bath at 50 °C for another 0.5 h. Finally, the mixture was transferred to a ceramic boat to devolatilize the light component under Ar at 300 °C for 3.0 h.

2.2. Catalytic oxidation of heavy oil

Thermal gravimetric analysis (TGA, Mettler Toledo, TGA/DSC 1) under an O2 atmosphere was carried out to investigate the oxidation performance of the heavy oil and the catalytic activity of the MnO2. The mass change and heat flow were measured in a temperature range of 30-550 °C with a ramp rate of 5 °C/min and a feed rate of 90 ml/min. The gas phase product evolved from oil oxidation were sampled into a mass spectrometer (MS, Hiden, QIC-20) by a capillary tube heated at 200 °C to measure the profiles of typical products and evaluate the catalytic effect of the MnO2 catalyst. Since the concentration of each product in the gas phase was very low, a secondary electron multiplier (SEM) detector was employed to improve the measuring accuracy. The products were identified by their molecular fragment ions at different mass-to-charge ratios (m/z). The m/z values of 15, 18, 43, and 44 were used for CH4, H2O, C3H8, and CO2, respectively. The same experiment was repeated at least twice to ensure repeatability and accuracy.

2.3. Catalyst characterization

The morphology of the catalyst was determined by a JEM 2010 TEM (JEOL Ltd., Japan) equipment operated at 120 kV. The crystal phase of the catalysts was determined by wide angle XRD on a D8 ADVANCE Rigaku D/max-RB diffractometer (Bruker, Germany) using Cu Kα radiation at a voltage and current of 40 kV and 120 mA. The chemical state of the elements on the surface of the catalyst was identified by XPS on an ESCALAB 250Xi spectrometer using Al Kα as excitation source. The binding energies (BEs) of the Au 4f, Mn 2p, and O 1s were calibrated against the C 1s signal (284.6 eV) of contaminant carbon.

3. Results and discussion
3.1. Morphology and electronic structure of the MnO2 catalyst

As observed in Fig. 1, the MnO2 catalyst displayed a high crystallinity of the β-MnO2 phase (JCPDS 44-0141). No diffraction peak from an impurity was detected. Irregular shapes and a polycrystalline structure with rough edges were observed in the high resolution TEM images (Fig. 2).

Fig. 1. XRD pattern of the MnO2 catalyst.

Fig. 2. TEM (a), HRTEM images (b), and selected area electron diffraction pattern (c) of the MnO2 catalyst. The inset in (a) exhibits the morphology of the MnO2 catalyst.

The chemical state of the O and Mn species on the surface of MnO2 was investigated by XPS to help understand its catalytic performance. The asymmetric peak located at 642.3 eV (Fig. 3(a)) indicated the existence of Mn(IV) [37], which was in good accordance with the XRD result. There were various surface O species since the O 1s spectrum could be deconvoluted into three peaks [37, 38]: at 529.8-530.1 eV (OI) assigned to lattice oxygen, 531.3 eV (OII) assigned to surface adsorbed oxygen, OH groups and oxygen vacancies, and 533.0 eV (OIII) associated with adsorbed molecular water (Fig. 3(b)). The OH groups and oxygen vacancies were the most important species for catalytic oxidation because of their high mobility and excellent activity of oxygen. From the data in Table 2, we observed that 15.7% of surface oxygen was OII species, indicating a high catalytic performance in the following experiments.

Fig. 3. Mn 2p (a) and O 1s (b) XPS spectra of the MnO2 catalyst.

Table 2
XPS results for the MnO2 catalyst.
3.2. Catalytic oxidation of heavy oil

The oxidation and catalytic oxidation performance of the heavy oil were evaluated by TGA. The TG curves in Fig. 4 exhibited that the oxidation of heavy oil involved three stages: low temperature oxidation (50-300 °C), plateau section (300-350 °C), and high temperature (350-500 °C) oxidation.

Fig. 4. TG (a) and DSC (b) of heavy oil oxidation with and without the MnO2 catalyst at low temperature.

At low temperature, partial oxidation of oil components was dominant and the weight loss rate was smooth. Oxidative addition reactions are believed to dominate the reactions because a large amount of data from both laboratory and field trials revealed a significant viscosity increase of the oil produced after the low temperature period. However, in this case, since the oil was devolatilized at 300 °C before the test, a significant weight loss was observed in the temperature interval from 200 to 350 °C (Fig. 4(a)) indicating the occurrence of oxidative decomposition reactions. Both oxidative addition and cracking reactions occurred at low temperature, and each was dominant in a different temperature interval. At temperatures less than 200 °C, oxidative addition was the dominant reaction leading to little weight loss, while at 200-350 °C, oxidative decomposition became the dominant reaction resulting in a large weight loss. This indicated the possibility for the decrease of fuel deposition at the end of low temperature period using oxidative decomposition promoted by MnO2 additives. The TG profile of catalytic oxidation by the MnO2 particles displayed a larger weight loss. Total weight loss of the oil sample with MnO2 at 350 °C was increased to 26.5%, almost twice that of the oil without the catalyst (14.3%). Despite its increase in weight loss, the amount of released heat at low temperature was also enhanced. In Fig. 4(b), the peak heat flow exhibited an increase of 17.6% from 1.7 W/g to 2.0 W/g. Figure 4 demonstrated that the MnO2 catalyst facilitated the oxidative decomposition of heavy compounds and gave a promotion in heat generation at low temperature (< 300 °C). However, further evaluation of its catalytic effect on the high temperature exothermic behavior was still necessary since the high temperature combustion contributes most of heat required for the in situ combustion.

Figure 5 illustrated the weight loss and heat release from oil oxidation in the high temperature range. The TG curves exhibited an extremely sharp weight loss as the main reactions shifted to combustion, and the end temperature was lowered because of the faster reaction rate and less residue from the low temperature oxidation. In comparison, the amount of heat released from high temperature combustion (Fig. 5 and Table 3) did not decrease despite less fuel deposition, indicating an enhancement of combustion heat by MnO2.

Fig. 5. TG (a) and DSC (b) of heavy oil oxidation with and without the MnO2 catalyst at high temperature.

The TG and DSC profiles illustrated that MnO2 decreased the amount of fuel deposition in the low temperature range, and gave the same heat from the combustion of less fuel in the high temperature range. Thus, the use of MnO2 is expected to reduce the deposition of excessive fuel in heavy reservoirs and will make in situ combustion more feasible.

3.3. Mass spectrometric (MS) analysis of oil oxidation product

The products of oil oxidation were sent to a mass spectrometer to track the catalytic oxidation of heavy oil. Inorganic products from complete oxidation (CO2, H2O) and organic small moleculars (CH4, C3H8) were selected as representative products. The comparison of the product flow between oil oxidation without and with MnO2 is presented in Fig. 6. The figures showed an increase in product both from complete oxidation and oxidative cracking at low temperature with the addition of MnO2, in accordance with the promotion effect on low temperature weight loss. Compared with the doubled weight loss in the low temperature range, the observed small increase of the four products indicated the formation of other gaseous products. This would probably lead to an improvement in recovery efficiency. The catalytic effects of MnO2 on the four products in the high temperature range were prominent (Fig. 6). Since the experiments were carried out with excessive O2 flow, the detection of products from incomplete oxidation at high temperature indicated that oil combustion was a diffusion controlled process and excess O2 does not always lead to complete combustion of the oil. The concentrations of CH4, C3H8, and H2O were observed to be significant decreased, while the production of CO2 averaged over the entire high temperature range gave a cumulative amount that revealed no obvious change. Considering that there was less fuel deposited for combustion, it was the nearly constant amount of CO2 that explained why there was no change in heat generation at high temperature compared with the experiment without the catalyst. These results confirmed the promotion in the oxidation degree of the product. On the basis of the MS data, it was speculated that MnO2 facilitated the oxidative decomposition of heavy components at low reaction temperatures, accelerated the transfer rate of oxygen in the oil phase, and boosted the complete oxidation of oil components at high temperatures, which consequently led to the promotion of low temperature cracking and high temperature oxidation heat.

Fig. 6. Mass spectra of molecular fragment ions at m/z = 44 (CO2) (a), m/z = 18 (H2O) (b), m/z = 15 (CH4) (c), m/z = 43 (C3H8) (d).

Table 3
TG and DSC results of heavy oil oxidation with and without the MnO2 catalyst.
3.4. The DAEM kinetic of catalytic oxidation

The distributed activation energy method (DAEM) was employed as the kinetic analysis to provide more insight into the catalytic mechanism. The DAEM method, which is widely used in the kinetic analysis of systems with complex reactions, assumes that many irreversible first order parallel reactions with different rate parameters are involved in the system and only one reaction occurs at a specified temperature. Therefore, the overall rate dV/dt is approximated by the rate of the j-th reaction at the temperature at which only the j-th reaction is occurring. The mathematical equation is given by

where V is the total volatile evolved by time t, V is total volatiles of the sample, f(E) is a distribution curve of the activation energies representing the differences in the activation energies of many first order irreversible reactions, and k0 is the frequency factor corresponding to the E value.

After mathematical transformation, Eq. (1) can be rewritten as,

where β is the heating rate. Thus, the E and k0 value corresponding to the V/V value is determined by the Arrhenius plot of ln(β/T2) vs. 1/T at the same V/V value. The activation energy of heavy oil oxidation with MnO2 catalyst was calculated by a method similar to that reported in previous publications [39, 40, 41, 42].

The activation energy versus temperature relationship in the whole temperature range is illustrated in Fig. 7. The symbols on the profiles of the activation energy stand for the percent of weight loss from 1% to 100% with an interval of 1%. Similarly to the TG profiles, the two profiles showed similar tendency and were divided into three regions: low temperature region, transition region, and high temperature region. The activation energy was less than 150 kJ/mol in the low temperature interval and 200 kJ/mol in the high temperature interval. With the addition of the β-MnO2 catalyst, the apparent activation energy decreased by 10-30 kJ/mol at low temperature and by 20-40 kJ/mol at high temperature. Since the apparent activation energy is defined as the minimum energy required for atoms or molecules to undergo chemical transformation, the decrease of activation energy was attributed to the reduction of the potential barrier and increase in oxygen transfer coefficient from the contribution of mobile oxygen vacancies on the surface of MnO2. The oil mixed with MnO2 exhibited an earlier weight loss and the limits of the transition region were also shifted to higher conversion. This indicated that more reactants were involved in the low temperature oxidation and more light products were available, which was consistent with the result of the TG experiment.

Fig. 7. Apparent activation energy of oxidation from the DAEM method.
4. Conclusions

A β-MnO2 nanoparticle catalyst with rough polycrystalline edges and abundant absorbed oxygen and oxygen vacancies was used for in situ combustion to enhance the low temperature decomposition of heavy oil components and to promote heat release at high temperature. TGA experiments demonstrated that the total weight loss at low temperature was almost doubled from 14.3% to 26.5% by the addition of the β-MnO2 catalyst, and the heat release from high temperature combustion did not decline despite less fuel deposition. MnO2 promoted the low temperature cracking and the high temperature oxidation heat by its catalytic effect on the low temperature oxidative decomposition of heavy components and the increase in the transfer rate of oxygen in the oil phase. The apparent activation energies obtained by the DAEM method were decreased, which confirmed the catalytic role of the MnO2 catalyst. The outstanding performance of oil oxidation with the addition of the MnO2 catalyst would expand the applications of thermal recovery strategies and afford improved produced oil with increased oil recovery efficiency.

Acknowledgments

The authors thank Dr. Cheng Zan at State Key Laboratory of Enhanced Oil Recovery, Research Institute of Petroleum Exploration & Development, China National Petroleum Corporation and Prof. Lin Shi at Department of Thermal Engineering, Tsinghua University for helpful discussion.

References
[1] Chu S, Majumdar A. Nature, 2012, 488: 294
[2] Al-Bahlani A M, Babadagli T. Chem Eng J, 2012, 181: 557
[3] Vogel G H. Chem Eng Technol, 2008, 31: 730
[4] Keim W. Petrol Chem, 2010, 50: 298
[5] Li Y F, Wang H F, Wang G, Gao J S. Chem Eng J, 2012, 211: 255
[6] Shah A, Fishwick R, Wood J, Leeke G, Rigby S, Greaves M. Energy Environ Sci, 2010, 3: 700
[7] Brüggemann P, Baitalow F, Seifert P, Meyer B, Schlichting H. Fuel Process Technol, 2010, 91: 211
[8] Zhao D W, Wang J, Gates I D. Fuel, 2014, 117: 431
[9] Chu Y, Fan C, Zhang Q, Zan C, Ma D S, Jiang H, Wang Y, Wei F. Chem Eng J, 2014, 248: 422
[10] Moore R G, Laureshen C J, Ursenbach M G, Mehta S A, Belgrave J D M. J Can Petrol Technol, 1999, 38(13): 96
[11] Weissman J G, Kessler R V, Sawicki R A, Belgrave J D M, Laureshen C J, Mehta S A, Moore R G, Ursenbach M G. Energy Fuels, 1996, 10: 883
[12] Shah A, Fishwick R P, Leeke G A, Wood J, Rigby S P, Greaves M. J Can Petrol Technol, 2011, 50(11-12): 33
[13] Xia T X, Greaves M. J Can Petrol Technol, 2002, 41(8): 58
[14] Xia T X, Greaves M. In: SPE International Thermal Operations and Heavy Oil Symposium. Margarita Island: The Society of Petroleum Engineers, 2001. 69693
[15] Greaves M, Xia T X. J Can Petrol Technol, 2004, 43(9): 25
[16] Hart A. Int J Petrol Sci Technol, 2012, 6(2): 79
[17] Hashemi R, Nassar N N, Almao P P. Energy Fuels, 2013, 27: 2194
[18] Reservoir Engineering Section on in situ Combustion. Department of Energy of the United States. California, 1998
[19] The Effect of Metallic Additives on the Kinetics of Oil Oxidation Reactions in in situ Combustion. Department of Energy of the United States. California, 1988
[20] Zhao M Q, Zhang Q, Huang J Q, Wei F. Adv Funct Mater, 2012, 22: 675
[21] Zaera F. ChemSusChem, 2013, 6: 1797
[22] Sun X Y, Wang R, Su D S. Chin J Catal (孙晓岩, 王锐, 苏党生. 催化学报), 2013, 34: 508
[23] Montes A R, Gutierrez D, Moore R G, Mehta S A, Ursenbach M G. J Can Petrol Technol, 2010, 49(2): 56
[24] Chao K, Chen Y L, Liu H C, Zhang X M, Li J. Energy Fuels, 2012, 26: 1152
[25] Castanier L M, Brigham W E. In Situ, 1997, 21: 27
[26] Shallcross D C, De los Rios C F, Castanier L M, Brigham W E. SPE Reservoir Engineering, 1991, 6: 287
[27] Castanier L M, Baena C J, Holt R J, Brigham W E, Tavares C. In: Proceedings of the 2nd Latin American Petroleum Conference. Caracas: The Society of Petroleum Engineers, 1992. 23708
[28] Nares H R, Schachat-Hernandez P, Ramirez-Garnica M A, Cabrera-Reyes M C, Noe-Valencia L, La Salle U. In: Latin American and Caribbean Petroleum Engineering Conference. Buenos Aires: The Society of Petroleum Engineers, 2007. 107837
[29] Castanier L M, Brigham W E. J Petrol Sci Eng, 2003, 39: 125
[30] Racz D. In: Proceedings of European Meeting on Improved Oil Recovery. Rome, 1985
[31] In situ Combustion Handbook Principles and Practices. Department of Energy of the United States. California, 1999
[32] Shokrlu Y H, Maham Y, Tan X, Babadagli T, Gray M. Fuel, 2013, 105: 397
[33] Ramesh K, Chen L W, Chen F X, Liu Y, Wang Z, Han Y F. Catal Today, 2008, 131: 477
[34] Shi F J, Wang F, Dai H X, Dai J X, Deng J G, Liu Y X, Bai G M, Ji K M, Au C T. Appl Catal A, 2012, 433: 206
[35] Jiang F, Zhu X W, Fu B S, Huang J J, Xiao G M. Chin J Catal (姜枫, 朱晓文, 符宝嵩, 黄金金, 肖国民. 催化学报), 2013, 34: 1683
[36] Wang M X, Zhang P Y, Li J G, Jiang C J. Chin J Catal (王鸣晓, 张彭义, 李金格, 姜传佳. 催化学报), 2014, 35: 335
[37] Liang S H, Teng F, Bulgan G, Zong R L, Zhu Y F. J Phys Chem C, 2008, 112: 5307
[38] Chen C M, Zhang Q, Yang M G, Huang C H, Yang Y G, Wang M Z. Carbon, 2012, 50: 3572
[39] Miura K. Energy Fuels, 1995, 9: 302
[40] Miura K, Maki T. Energy Fuels, 1998, 12: 864
[41] Fan C, Zan C, Zhang Q, Ma D S, Chu Y, Jiang H, Shi L, Wei F. Fuel Process Technol, 2014, 119: 146
[42] Liu X G, Li B Q, Miura K. Fuel Process Technol, 2001, 69: 1