SKI-210催化二甲苯異構(gòu)化反應(yīng)的動力學(xué)研究

馮培曦，周震寰，康承琳，蓋月庭

(中國石化石油化工科學(xué)研究院，北京 100083)

對二甲苯(PX)是重要的石油化工基礎(chǔ)原料，工業(yè)上為了滿足PX生產(chǎn)的需求，通常采用二甲苯異構(gòu)化工藝，將相對過剩的鄰二甲苯(OX)和間二甲苯(MX)轉(zhuǎn)化為PX[1]。催化劑是二甲苯異構(gòu)化工藝的核心，中國石化石油化工科學(xué)研究院通過在傳統(tǒng)二甲苯異構(gòu)化催化劑中引入改性組元，成功研發(fā)了SKI-210催化劑[2]。該催化劑具有活性高、選擇性好和穩(wěn)定周期長的優(yōu)點，產(chǎn)生了良好的經(jīng)濟效益[3]。反應(yīng)動力學(xué)是反應(yīng)器設(shè)計和工藝模擬優(yōu)化的重要依據(jù)，同時研究反應(yīng)動力學(xué)也有助于加深對整個反應(yīng)過程的理解[4]。

對于二甲苯異構(gòu)化過程，Okada等[5]使用固定床反應(yīng)器研究了在鎳離子交換的硅鋁催化劑上OX的異構(gòu)化反應(yīng)動力學(xué)模型，選擇基于Langmuir吸附的雙曲型動力學(xué)方程，利用初始速率法求解動力學(xué)參數(shù)。Corma等[6]使用類似的方法建立了二甲苯異構(gòu)化反應(yīng)整體的動力學(xué)模型。Iliyas等[7]使用提升管反應(yīng)器研究了二甲苯在USY分子篩上的異構(gòu)化反應(yīng)動力學(xué)，并加入了歧化副反應(yīng)，選用一級動力學(xué)方程，使用單一組分作為原料簡化反應(yīng)網(wǎng)絡(luò)，利用非線性擬合的方法求解動力學(xué)參數(shù)。陳金仙等[8]使用固定床反應(yīng)器研究了在HZSM-5分子篩上MX異構(gòu)化反應(yīng)的動力學(xué)，選用一級動力學(xué)模型，使用最優(yōu)化算法計算得到模型參數(shù)。

本課題利用固定床反應(yīng)器研究SKI-210催化二甲苯異構(gòu)化反應(yīng)動力學(xué)，選擇基于Langmuir吸附的動力學(xué)模型，利用初始速率法求解模型參數(shù)，并通過得到的模型研究反應(yīng)條件對二甲苯異構(gòu)化反應(yīng)的影響。

1 實 驗

1.1 原 料

試驗所用反應(yīng)原料為OX、MX，均為分析純，質(zhì)量分數(shù)大于99%，天津市大茂化學(xué)試劑廠提供。所用SKI-210催化劑為工業(yè)成品，中國石油撫順石化公司催化劑廠提供。

1.2 試驗裝置

試驗在固定床反應(yīng)器上進行。反應(yīng)器為長500 mm、內(nèi)徑6 mm的不銹鋼管，利用電熱爐加熱控制溫度。反應(yīng)器出口壓力由背壓閥控制。反應(yīng)氣相產(chǎn)物通過六通閥切入氣相色譜分析。氣相色譜為安捷倫7890A，HP-WAX毛細柱(60 m×0.25 mm×0.5 μm)，F(xiàn)ID檢測器。采用面積歸一化法定量。

1.3 試驗條件

整個動力學(xué)實驗均在催化劑活性穩(wěn)定期內(nèi)進行。反應(yīng)溫度選擇350～380 ℃，反應(yīng)壓力選擇0.3～0.6 MPa。催化劑的裝填量為0.2 g，反應(yīng)的質(zhì)量空速為26～78 h-1。本試驗為研究該催化反應(yīng)的宏觀動力學(xué)，因此將催化劑的內(nèi)擴散影響包含在動力學(xué)常數(shù)內(nèi)，試驗選擇的空速較大，因此忽略外擴散的影響。

2 反應(yīng)動力學(xué)模型

二甲苯異構(gòu)化反應(yīng)過程的主反應(yīng)為3種異構(gòu)體的相互轉(zhuǎn)化，對于SKI-210催化劑，二甲苯異構(gòu)體之間的轉(zhuǎn)化為順序型，即OX和PX之間不能直接轉(zhuǎn)化[9]。從試驗結(jié)果看，發(fā)生的副反應(yīng)為二甲苯的歧化。因此，整個反應(yīng)過程的反應(yīng)網(wǎng)絡(luò)為：

式中：k1，k-1，k2，k-2為主反應(yīng)速率常數(shù)，mol/(h·g)；kO，kM，kP為副反應(yīng)速率常數(shù)，mol/(h·g·MPa)；T為副產(chǎn)物，包括甲苯和三甲苯。

2.1 模型假設(shè)

在建立數(shù)學(xué)模型時，為了使模型簡化，便于計算，通常需要進行合理的假設(shè)。在本研究中作以下假設(shè)：①由于催化劑裝量較少，床層高度較低，因此假設(shè)催化劑床層內(nèi)的流體流動滿足固定床柱塞流模型，忽略沿反應(yīng)器方向的床層壓降；②二甲苯在催化劑上的吸附行為滿足Langmuir吸附，異構(gòu)化動力學(xué)模型選擇基于Langmuir吸附的雙曲型動力學(xué)模型；③反應(yīng)產(chǎn)物中甲苯和三甲苯的量很少，因此將副反應(yīng)視為一級不可逆反應(yīng)。

2.2 動力學(xué)模型的建立

根據(jù)以上假設(shè)，得到反應(yīng)體系內(nèi)各組分的反應(yīng)速率方程為：

-rM=[(k1+k2)KMPM-k-1KOPO-k-2KPPP]/
(1+∑KjPj)+kMPM

(1)

-rO=(k-1KOPO-k1KMPM)/
(1+∑KjPj)+kOPO

(2)

-rP=(k-2KPPP-k2KMPM)/
(1+∑KjPj)+kPPP

(3)

rT=kOPO+kMPM+kPPP

(4)

式中：r為反應(yīng)速率，mol/(h·g)；K為吸附平衡常數(shù)，MPa-1；P為各組分的分壓，MPa；下標O，M，P，T分別為OX、MX、PX和副反應(yīng)產(chǎn)物。可以通過繪制反應(yīng)物轉(zhuǎn)化率與停留時間的關(guān)系曲線，在曲線上求取某一點的斜率來確定該點對應(yīng)的反應(yīng)速率。根據(jù)柱塞流反應(yīng)器模型，停留時間可以通過W和F計算，其中：W為裝填催化劑的質(zhì)量，g；F為反應(yīng)物的流量，mol/h。

2.3 動力學(xué)參數(shù)的求解

由于反應(yīng)網(wǎng)絡(luò)中包含可逆反應(yīng)，導(dǎo)致動力學(xué)模型較復(fù)雜，因此采用初始速率法來簡化求解動力學(xué)參數(shù)。利用純OX或MX作為原料，在反應(yīng)初始時，可將反應(yīng)網(wǎng)絡(luò)簡化為：

反應(yīng)動力學(xué)方程可簡化為：

(5)

(6)

式中：ri為反應(yīng)原料初始反應(yīng)速率，mol/(h·g)；P0為反應(yīng)物的初始分壓，MPa。

2.3.1k1，k2，kM的求解

當進料為純MX時，由式(5)可知，產(chǎn)物OX和PX以及副產(chǎn)物的初始生成速率(r0)可表示為：

(7)

(8)

rT，0=kMPM，i

(9)

反應(yīng)的初始速率使用圖解法求得。通過繪制反應(yīng)物轉(zhuǎn)化率和W/F的關(guān)系曲線，在曲線上求取原點處的斜率即為初始反應(yīng)速率，其結(jié)果如表1所示。

對式(7)和式(8)進行變換可得：

(10)

(11)

繪制1/PM，i與1/rP，0和1/rO，0的關(guān)系曲線，擬合求得k1，k2，KM，如圖1和圖2所示。根據(jù)式(9)，繪制PM，i與rT，0的關(guān)系曲線，如圖3所示，擬合求得kM。

表1 以MX為原料時不同溫度和壓力下的產(chǎn)物初始生成速率

圖1 以MX為原料時PX初始生成速率與MX初始壓力關(guān)系擬合曲線溫度，℃：▲—350； ◆—360； ■—370； ●—380。圖2～圖5同

圖2 以MX為原料時OX初始生成速率與MX初始壓力關(guān)系擬合曲線

圖3 以MX為原料時副反應(yīng)產(chǎn)物初始生成速率與MX初始壓力關(guān)系擬合曲線

2.3.2k-1、kO的求解

以純OX為原料，根據(jù)試驗結(jié)果發(fā)現(xiàn)，k-1與kO相差很大，因此在計算k-1時忽略kO的影響。由式(6)可知，OX的初始反應(yīng)速率和副產(chǎn)物的初始生成速率分別為：

(12)

rT,0=kOPO,i

(13)

對式(12)變形，得：

(14)

OX的初始反應(yīng)速率通過繪制OX的轉(zhuǎn)化率與W/F的關(guān)系曲線，在曲線上求取原點處的斜率得到，其結(jié)果如表2所示。繪制1/PO，i與-1/rO，i的關(guān)系曲線，擬合求得k-1、KO，如圖4所示。繪制PO，i與rT，0的關(guān)系曲線，如圖5所示，擬合求得kO。

表2 以O(shè)X為原料時不同溫度和壓力下的反應(yīng)初始速率

圖4 以O(shè)X為原料時OX反應(yīng)初始速率與OX初始壓力關(guān)系擬合曲線

圖5 以O(shè)X為原料時副反應(yīng)產(chǎn)物初始生成速率與OX初始壓力關(guān)系擬合曲線

2.3.3k-2、kP求解

使用純PX為原料進行試驗時，由于PX能夠進入催化劑孔道，導(dǎo)致積炭嚴重阻塞孔道，催化劑失活過快。因此根據(jù)式(1)，用MX試驗數(shù)據(jù)通過非線性擬合的方法計算得到k-2。副反應(yīng)速率較小，對反應(yīng)整體影響不大，且反應(yīng)速率常數(shù)與OX和MX反應(yīng)速率常數(shù)相差不大，可直接用后二者的平均值計算。

3 結(jié)果與討論

3.1 動力學(xué)參數(shù)

根據(jù)上述方法求解得到各反應(yīng)溫度下的動力學(xué)參數(shù)及吸附平衡常數(shù)，如表3所示。動力學(xué)參數(shù)和溫度之間的關(guān)系滿足Arrhenius方程：

式中：k0為指前因子，mol/(h·g)；E為活化能，kJ/mol。利用不同溫度下的反應(yīng)速率常數(shù)，對1/T和lnkj擬合得到反應(yīng)的活化能和指前因子，結(jié)果如表4所示。

從表3和表4可以看出，二甲苯在催化劑上的吸附是放熱反應(yīng)，可自發(fā)進行。PX的吸附平衡常數(shù)比OX和MX的高，這是因為PX的特殊構(gòu)型使其能進入分子篩孔道內(nèi)吸附。PX到MX轉(zhuǎn)化反應(yīng)的活化能低于其他各主反應(yīng)，說明PX更易于轉(zhuǎn)化為MX。

表3 不同溫度下的反應(yīng)速率常數(shù)及吸附平衡常數(shù)

表4 各反應(yīng)速率常數(shù)對應(yīng)的活化能和指前因子對數(shù)值

3.2 模型驗證

為了驗證計算所得模型參數(shù)的可靠性，使用混合二甲苯為原料在壓力為0.5 MPa下進行試驗。原料組成如表5所示。將試驗所得產(chǎn)物的組成與模型計算值進行比較，結(jié)果如圖6所示。

表5 混合二甲苯原料組成 x，%

從圖6可以看出，模型計算值與試驗值吻合較好，誤差均在5%以內(nèi)，說明模型求解所得動力學(xué)參數(shù)是可靠的。

圖6 各溫度下試驗值與計算值比較試驗值：◆—MX； ▲—OX； ■—PX； ●—T計算值： —MX； —OX； —PX； —T

3.3 反應(yīng)條件的影響

對于二甲苯異構(gòu)化反應(yīng)，最重要的指標為PX收率。影響PX收率的因素主要有空速、溫度和壓力，以下以混合二甲苯(表5)為原料，利用求得的模型來考察這些因素對PX收率的影響。

3.3.1空速的影響

在反應(yīng)溫度分別為350，360，380 ℃、壓力為0.5 MPa時，考察空速對PX收率的影響，結(jié)果如圖7所示。從圖7可以看出：空速較高時，隨著空速越低，PX收率越高，這是由于空速越低，反應(yīng)物與催化劑接觸時間越長，反應(yīng)轉(zhuǎn)化率越高；在溫度為380 ℃、質(zhì)量空速約為12 h-1時，PX收率出現(xiàn)峰值，約為23%，這是由于在此空速下，主反應(yīng)達到了平衡；之后隨著空速降低，PX收率逐漸降低，副反應(yīng)產(chǎn)物較多。

圖7 進料空速對PX收率的影響溫度，℃： —350； —360； —380

3.3.2反應(yīng)溫度的影響

在反應(yīng)壓力為0.5 MPa、質(zhì)量空速分別為5，10，20 h-1時，考察反應(yīng)溫度對PX收率的影響，結(jié)果如圖8所示。從圖8可以看出：隨著反應(yīng)溫度升高，各空速下PX收率均升高，這是由于溫度升高，反應(yīng)速率常數(shù)增大，反應(yīng)加快，提高了反應(yīng)物的轉(zhuǎn)化率；在反應(yīng)溫度達到375 ℃后，隨著溫度進一步升高，各空速下PX收率增長均變緩，這是由于隨著反應(yīng)體系內(nèi)PX濃度增加，逆反應(yīng)加快，反應(yīng)趨于平衡；反應(yīng)達到平衡時，PX收率約為24%。

圖8 反應(yīng)溫度對PX收率的影響質(zhì)量空速，h-1： —5； —10； —20

3.3.3反應(yīng)壓力的影響

在反應(yīng)溫度為380 ℃、質(zhì)量空速分別為5，10，50，100 h-1時，考察反應(yīng)壓力對PX收率的影響，結(jié)果如圖9所示。從圖9可以看出：在質(zhì)量空速高于10 h-1時，隨著反應(yīng)壓力升高，PX收率均升高，這是由于壓力升高，反應(yīng)物濃度增大，反應(yīng)加快，反應(yīng)物轉(zhuǎn)化率升高；當空速較低時，如質(zhì)量空速為5 h-1時，隨著反應(yīng)壓力逐漸升高，PX收率出現(xiàn)下降趨勢，這是由于在低空速下，主反應(yīng)已達平衡，隨著壓力升高，副反應(yīng)加快，使得PX收率降低。為了使工業(yè)裝置有較高的處理量，同時能夠維持較高的PX收率，選擇質(zhì)量空速為10 h-1、反應(yīng)壓力為0.6 MPa較為合理，此時PX收率約為24%。

圖9 反應(yīng)壓力對PX收率的影響質(zhì)量空速，h-1： —5； —10； —50； —100

4 結(jié) 論

(1)根據(jù)二甲苯異構(gòu)化反應(yīng)產(chǎn)物分布，建立了二甲苯異構(gòu)化反應(yīng)網(wǎng)絡(luò)。在該反應(yīng)網(wǎng)絡(luò)中，主反應(yīng)為3種二甲苯異構(gòu)體之間的相互轉(zhuǎn)化，副反應(yīng)為二甲苯的歧化。

(2)利用基于Langmuir吸附的雙曲型模型建立了二甲苯異構(gòu)化反應(yīng)的動力學(xué)模型，使用初始速率法設(shè)計試驗并求解了模型的動力學(xué)參數(shù)。模型計算所得反應(yīng)體系組成與試驗結(jié)果吻合較好，相對誤差在5%以內(nèi)。

(3)利用所建模型考察操作條件對PX收率的影響，結(jié)果符合二甲苯異構(gòu)化反應(yīng)的試驗結(jié)果。在質(zhì)量空速為10 h-1、反應(yīng)溫度為380 ℃、反應(yīng)壓力為0.6 MPa時，可獲得較高的PX收率，約為24%，且該條件下裝置有較高的處理能力。