馬 磊，張新波，高 振

（1.西南化工研究設(shè)計院有限公司，四川 成都 610225；2.中海石油氣電集團有限責任公司技術(shù)研發(fā)中心，北京 100028）

CNJ甲烷化催化劑反應(yīng)本征動力學(xué)模型

馬 磊1，張新波1，高 振2

（1.西南化工研究設(shè)計院有限公司，四川 成都 610225；2.中海石油氣電集團有限責任公司技術(shù)研發(fā)中心，北京 100028）

選取280～580℃，1.0～2.5MPa和汽氣比0.1～0.7的反應(yīng)條件，研究了CNJ型合成氣甲烷化催化劑的本征動力學(xué)。實驗所用原料氣體積組成為：H256%～58%，CO 6%～7%，CO2≈3%，CH424%～25%，C2H61.5%～2%，其余為N2。利用Matlab對實驗數(shù)據(jù)采用冪函數(shù)模型回歸得到本征動力學(xué)方程。將數(shù)據(jù)擬合得到的本征動力學(xué)方程利用ASPEN PLUS軟件進行模擬，并與工業(yè)裝置結(jié)果進行對比驗證。結(jié)果表明，出口關(guān)鍵組分模擬值與實驗值接近，證明模型的合理性。

甲烷化；本征動力學(xué)；催化劑

目前，國內(nèi)天然氣供應(yīng)的缺口逐年增大，對外依存度近兩年呈快速上升之勢。據(jù)發(fā)改委能源研究所預(yù)測，到2020年我國天然氣的消費量將達到2500億m3，國內(nèi)天然氣缺口達1000億m3以上[1]。在我國富煤貧油少氣的資源環(huán)境下，利用低階褐煤分層利用生產(chǎn)天然氣具有重要意義。而且在煤炭深加工的幾個發(fā)展方向中，煤制天然氣的能量轉(zhuǎn)化效率高，單位能量水耗低、投資省，工程化建設(shè)成熟度高，是有效的煤炭利用方式和煤制能源產(chǎn)品較好途徑[2]。

甲烷化是煤制天然氣及焦爐煤氣制天然氣項目中的核心工段，而甲烷化催化劑又是其中的關(guān)鍵。本文的目的在于研究甲烷化催化劑的動力學(xué)特征[3]，為針對催化劑特征設(shè)計最佳反應(yīng)器提供基礎(chǔ)條件。

1 實驗部分

1.1 實驗裝置

實驗裝置流程[4]如圖1所示，試驗用原料氣采用配制的鋼瓶氣，經(jīng)過減壓計量后分別進入預(yù)熱器及蒸發(fā)器，進料水經(jīng)平流泵加壓進入蒸發(fā)器蒸發(fā)，同時與原料氣充分混合后進入反應(yīng)器發(fā)生反應(yīng)，反應(yīng)器采用外加熱形式控制出口溫度，出口氣體經(jīng)冷卻、分離、計量后一部分去色譜進行組分分析，另一部分直接排空。實驗采用PID控溫，以保證反應(yīng)器熱點溫度的一致性。

圖1 甲烷化實驗裝置示意圖Fig.1 Schematic diagram of methanation experimental installation

1.2 實驗方法

(1)催化劑還原

以催化劑床層出口處熱電偶為控制熱電偶，關(guān)閉放空閥，用氫氣緩慢升壓至 0.5MPa，時間為10min，系統(tǒng)定壓，緩慢打開放空閥，以4℃/min的速度升溫至400℃（見結(jié)果討論），氫氣空速2000h-1，溫度升至400℃時開始計時，還原時間為2h。

(2)甲烷化反應(yīng)

打開氮氣鋼瓶頂閥及控制面板開關(guān)通入氮氣，關(guān)閉穩(wěn)壓閥2-2，系統(tǒng)壓力升至試驗壓力時定壓，打開穩(wěn)壓閥2-2；打開總電源，進行程序升溫設(shè)置，以4℃/min的速率升至目標溫度：第一段，室溫～280℃，時間60min；第二段，280℃保溫，時間10min，同時投入原料氣；第三段，280℃至目標溫度。

1.3 分析方法

采用安捷倫7820A氣相色譜儀及微量CO2紅外分析儀（0～50×10-6）為分析設(shè)備，并對分析設(shè)備進行了標定。

常量組分分析采用氣相色譜儀進行，工作條件見表1。

表1 色譜工作條件Table 1 Working conditions of chromatography

微量CO2分析采用微量CO2紅外分析儀，操作條件為氣體流速為0.5mL/min。

2 結(jié)果與討論

2.1 實驗數(shù)據(jù)

選用CNJ型工業(yè)化甲烷化催化劑進行甲烷化實驗，條件實驗采用的原料氣為典型的合成氣組成，見表2。

表2 合成氣條件實驗原料氣組成Table 2 Composition of feed gas

在催化劑性能穩(wěn)定的基礎(chǔ)上，進行了動力學(xué)測試。測試選用的反應(yīng)溫度280～580℃，壓力1.0～2.5MPa，汽氣比（物質(zhì)的量比）0.1～0.7。反應(yīng)器催化劑裝填量為10.98g。實驗結(jié)果如表3所示。實驗期間，根據(jù)催化劑特點，選擇入口溫度為280℃進行測試。

表3 實驗結(jié)果表Table 3 Experimental data

2.2 數(shù)據(jù)回歸

甲烷化反應(yīng)是指合成氣中CO、CO2和H2在一定溫度、壓力和催化劑作用下生成CH4的反應(yīng)[5]。主反應(yīng)如下：

主反應(yīng)式（1）和（3）可以代替主反應(yīng)式（2），在數(shù)據(jù)回歸的過程中考慮（1）和（3）即可。通過出口入口分析數(shù)據(jù)，可以得知CNJ型催化劑同時具有促進變換反應(yīng)的功能，利用冪函數(shù)模型進行回歸。

由于甲烷化反應(yīng)為偏離平衡的反應(yīng)，可以得到動力學(xué)方程式如下[2]：

利用MATLAB軟件[6]采用4階Runge-Kutta法解微分方程組（4）和（6）[7-8]，對進出口條件進行擬合，阻尼最小二乘法進行參數(shù)估計，得到：

將 a3、a4、b3、b4圓整為零后得到的甲烷化反應(yīng)動力學(xué)模型如下：

2.3 模擬驗證

表4 甲烷化入口氣組成Table 4 Inlet gas composition of methanation reactor

為驗證甲烷化反應(yīng)動力學(xué)方程的可靠性，將使用動力學(xué)方程（7）和（8）的模擬軟件結(jié)果與采用CNJ催化劑的工業(yè)裝置結(jié)果進行比較[1]。模擬過程采用工業(yè)裝置運行實際數(shù)據(jù)，并盡可能排除其他因素對結(jié)果的影響。為避免偶然性，采取同一裝置不同時間的3組數(shù)據(jù)（表4）進行模擬[6]。

氣體流量Qn≈57000m3/h，壓力2.83MPa（G），入口溫度280℃，汽氣比0.30。參照以上條件利用ASPEN PLUS模擬軟件搭建反應(yīng)器模型（圖2）。模擬過程采用了RKS-BM熱力學(xué)方程[5]。該方程是為了提高SPK方程所進行的修正，其方程式為：

式中，ω為偏心因子。

此方程是一個廣泛被工程界接受和使用的立方形方程，它能準確的計算極性物質(zhì)及含氫物質(zhì)的p-V-T數(shù)據(jù)。

圖2 反應(yīng)模型示意圖Fig.2 Schematic diagram of reactor

與工業(yè)裝置一致，選擇絕熱反應(yīng)器為模擬研究對象進行模擬實驗。模擬結(jié)果與實際結(jié)果見表5。

表5 模擬結(jié)果與實驗數(shù)據(jù)對比Table 5 Comparison of simulated results and experimental data

由表5可以看出，模擬結(jié)果與工業(yè)裝置實際值差別極小，甲烷化反應(yīng)遠離平衡。

通過分析動力學(xué)方程的特性，可以得知甲烷化反應(yīng)中氣體組分的分壓對于反應(yīng)存在較大的影響，反應(yīng)過程中通過增加原料氣中H2O、CH4等含量會抑制甲烷化反應(yīng)的進行。由于甲烷化反應(yīng)為強放熱反應(yīng)，通過動力學(xué)方程的研究可以得知，可通過改變進入反應(yīng)器的H2O、CH4等含量來控制反應(yīng)器熱點溫度不超過設(shè)備的耐受極限；將動力學(xué)方程運用在工藝設(shè)計的過程中，通過計算不同空速條件下的反應(yīng)情況，可以得到甲烷化反應(yīng)過程中的最佳空速條件；根據(jù)動力學(xué)方程，可以得到壓力與反應(yīng)程度存在正相關(guān)，但過高的壓力條件會造成較大的能耗。甲烷化反應(yīng)的動力學(xué)方程為甲烷化工藝的設(shè)計、設(shè)備的選型、控制方案的確定提供了參考。

3 結(jié)論

(1)在280～580℃，1.0～2.5MPa及汽氣比0.1～0.7的反應(yīng)條件下，研究得出了CNJ型合成氣甲烷化催化劑的本征動力學(xué)方程；

(2)甲烷化催化劑本征動力學(xué)方程與實驗結(jié)果和工業(yè)數(shù)據(jù)結(jié)果吻合度高；

(3)動力學(xué)方程的研究能夠為工業(yè)裝置的工藝開發(fā)、設(shè)備設(shè)計、控制方案提供可靠的分析。

[1]侯建國,高振,王秀林,等.中國甲烷化工藝技術(shù)專利現(xiàn)狀及分析 [J].天然氣化工·C1化學(xué)與化工,2015,40(2): 66-70.

[2]Bell D A,Towler B F,Fan M H.Coal Gasification and its applications[M].2011.

[3]于建國,于遵宏,沈才大.半水煤氣耐硫甲烷化催化劑的研制及其本征動力學(xué)的測試 [J].化工學(xué)報,1991,(1):96-103.

[4]華新雷,王立剛,徐勇華.利用高通量技術(shù)研究鈷基費托催化劑的宏觀動力學(xué)[J].催化學(xué)報,2009,30(8):740-747.

[5]福格勒.化學(xué)反應(yīng)工程原理[M].第四版.北京:化學(xué)工業(yè)出版社,2006.

[6]Hahn B H,Valentine D T.Essential MATLAB for Engineers and Scientists[M].5th Ed.Academic Press, 2013.

[7]王健紅,馮樹波,杜增智.化工系統(tǒng)工程理論與實踐[M].北京:化學(xué)工業(yè)出版社,2009.

[8]包宗宏,武文良.化工計算與軟件應(yīng)用[M].北京:化學(xué)工業(yè)出版社,2011.

Intrinsic kinetics model of CNJ methanation catalyst

MA Lei1,ZHANG Xin-bo1,GAO Zhen2
(1.The Southwest Research and Design Institute of the Chemical Industry Co.,Ltd.,Chengdu 610225,China;2.Research and Development Center of CNOOC Gas and Power Group,Beijing 100028,China)

The intrinsic kinetics of methanation reaction over the CNJ catalyst was studied at the reaction conditions of 280-580℃,1.0-2.5MPa and the steam/syngas molar ratio of 0.1-0.7.A raw gas consisting of 56%-58%H2,6%-7%CO,3%CO2,24%-25%CH4,1.5%～2%C2H6,and N2for balance by mole was used for the experiments.Matlab was used to fit the experimental data by the power function model to get the intrinsic dynamic equation.The intrinsic kinetic equation obtained by fitting the data was simulated by ASPEN PLUS,and the results were compared with the data of the industrial plant.The results show that the simulated values of the key components at outlet are close to the experimental values,which proves the rationality of the model.

methanation;intrinsic kinetics;catalyst

TQ013.2；TQ221.11；TQ426.94

：A

：1001-9219(2016)06-87-04

2016-06-18；

：馬磊（1986-）,男，工程師，從事煤制天然氣甲烷化、低階煤分級利用等方面的研究，電話：028-85962275，E-mail：victor.mau@outlook.com。