懸浮床加氫裂化集總動力學模型的研究

張萍萍，杜 峰

(中國石油大學(華東)化學工程學院，山東 青島 266580)

懸浮床加氫裂化集總動力學模型的研究

張萍萍，杜 峰

(中國石油大學(華東)化學工程學院，山東 青島 266580)

在不同反應溫度、氫初壓條件下，通過高壓反應釜對克拉瑪依常壓渣油進行加氫裂化實驗，以此模擬懸浮床加氫裂化過程，并根據實驗數據及實際工藝中對各種輕油產品收率預測的需求建立了懸浮床加氫裂化六集總(氣體、汽油、柴油、蠟油、減壓渣油、焦)動力學模型，用Matlab軟件進行編程，采用最小二乘法對動力學參數進行估算，并進行誤差分析。結果表明，建立的六集總動力學模型能很好地對各集總產品收率進行預測，計算結果與實驗值基本吻合，大部分誤差在5%以內。

懸浮床 加氫裂化 集總動力學 參數估算 誤差分析

近年來，世界各國對輕質油品的需求日益增加，然而常規原油資源日益枯竭，世界原油不斷向重質化、劣質化的趨勢發展，再加上環保要求的不斷提高，因此，應大力發展重油加氫工藝來提高重油的加工深度、提高輕油收率、增加經濟效益[1-2]。懸浮床加氫裂化因其原料適應性強、輕油收率高、催化劑分散效果佳、經濟效益好等優點[3]，具有很好的發展前景，但也遇到了易生焦等瓶頸，加大對懸浮床加氫裂化動力學的研究對于反應器設計、操作條件優化等有重要的意義，可以為懸浮床加氫裂化工藝的工業化推廣奠定理論基礎。

懸浮床加氫裂化反應體系復雜，組分數目較多且各組分之間強偶聯，各種可逆反應、順序反應及平行反應同時存在[4]，所以一般采用集總的方法。所謂集總，就是按各類分子的動力學特性，將反應體系劃分成若干個集總組分，在動力學研究中把每個集總作為虛擬的單一組分來考察，建立集總動力學模型。Mosby等[5]建立了渣油加氫處理七集總動力學模型，楊朝合等[6]建立了渣油加氫裂化窄餾分集總動力學模型。劉晨光等[7]建立了孤島渣油在分散型催化劑存在下的加氫裂化六集總動力學模型。本課題以劉晨光等的部分研究結果為基礎，為滿足實際工藝中對輕油各產品收率預測的要求，對產品進行更加細致的集總劃分，建立懸浮床加氫裂化六集總動力學模型。

1 實 驗

1.1 原料油性質

對克拉瑪依常壓渣油(KAR)的基本性質進行分析，結果如表1所示。

表1 KAR的組成與性質

1.2 實驗流程

通過對KAR進行高壓釜反應來模擬懸浮床加氫裂化反應過程，其中原料油的添加量為(200±1) g，實驗流程[8]如圖1所示。

圖1 高壓釜操作及產物分離流程示意

1.3 催化劑、助劑與硫化劑

催化劑采用實驗室自制的油溶性鉬鎳復配催化劑(鉬鎳質量比為3∶2)，加入量為150 μg/g；助劑為十二烷基苯磺酸鈉(SDBS)，添加量為200 μg/g；硫化劑為硫粉，添加量為400 μg/g。

2 集總動力學模型的建立

根據沸程劃分集總，將反應體系劃分為氣體、汽油、柴油、蠟油、減壓渣油、焦6個集總，并作如下假設：①生焦反應為表觀二級反應，其它反應均為表觀一級反應；②忽略催化劑的失活反應。

根據上述集總劃分及假設，得到如圖2所示的反應網絡。

圖2 六集總動力學模型反應網絡

根據上述反應網絡推導出如下的動力學方程組：

dx1/dt=k1x2+k2x3+k3x4+k4x5

dx2/dt=-k1x2+k5x3+k6x4+k7x5

dx3/dt=-k2x3-k5x3+k8x4+k9x5

dx4/dt=-k3x4-k6x4-k8x4+k10x5

dx5/dt=-k4x5-k7x5-k9x5-k10x5-k11x52

dx6/dt=k11x52

式中：xi(i=1～6)代表各集總質量分數，t代表反應時間，kj(j=1～11)代表集總間的速率常數。

3 模型參數估算

參數估算實際上是個最優化的問題，目標函數為實驗值與計算值之間的殘差，即：

采用最小二乘法，用Matlab軟件進行編程，根據實驗數據對參數進行估算。其中，動力學常數的普遍關系式可以表示成溫度和氫初壓的函數，即：

ki=Aie-Ei/RTPai

式中：Ai為指前因子；Ei為反應活化能；ai為氫初壓指數；T為反應溫度；P為氫初壓。

由此可得如下關系式：

lnki與lnp、lnki與-1/T成一定的線性關系，根據估算得到的動力學參數k可以線性回歸得出氫初壓指數、活化能及指前因子。

4 結果與討論

對氫初壓為7 MPa、不同反應溫度(420，430，440，450 ℃)條件下的KAR加氫裂化反應情況進行研究，反應時間分別為0，20，40，60 min，各產物收率如表2所示。

表2 不同反應溫度下的產物分布

反應時間控制步驟如下：在其它實驗準備工作完善之后，通過高壓反應釜溫度控制器對釜內的原料油進行程序升溫，當反應體系達到要求的反應溫度(如430 ℃)時，開始計時，這時的反應時間被認為是0，其它反應時間都是以此為基準進行相應的延長。

由表2可以看出，在不同的反應溫度下，反應時間為0時，都有一定量的氣體生成，在430，440，450 ℃條件下，反應時間為0時還有焦生成，說明反應時間為0時已經進行了一定程度的加氫裂化反應。在同一反應條件下，隨著反應溫度的升高，氣體、汽油餾分、柴油餾分收率不斷增加，蠟油餾分收率則略有減少，生焦量增加。因為懸浮床加氫裂化反應主要按自由基熱裂化機理進行[9-10]，自由基熱裂化隨著反應溫度的升高而被促進，大分子自由基裂化成小分子化合物，氣體、輕油收率增加，對于強放熱的加氫裂化反應[11]，溫度升高不利于反應的進行，溫度升高促進了裂化反應，大分子化合物，一方面裂化為小分子化合物一方面又發生了縮聚反應[9]，生焦量增加。

對反應溫度為430 ℃、不同氫初壓(7，8，9 MPa)條件下的KAR加氫裂化反應情況也進行了研究，反應時間分別為0，20，40，60 min，各產物收率如表3所示。

表3 不同氫初壓條件下的產物分布

由表3可以看出：在同一反應條件下，隨著氫初壓的增加，氣體、汽油餾分、柴油餾分收率表現出減少的趨勢；蠟油餾分收率增加，但變化較小；生焦量減少。

由于加氫裂化總體來說是分子數減少的反應，氫初壓增加，對加氫裂化反應是有利的[11]，但是氫初壓增加導致氫自由基增多，會大大抑制烴自由基的裂化速率，后者的影響要大于前者，所以隨著氫初壓增加，氣體、輕油收率減少，生焦量減少。

采用最小二乘法，通過Matlab編程對不同溫度及氫初壓條件下的反應速率常數進行估算，同時根據線性回歸得到氫初壓指數、反應活化能及指前因子，如表4、表5所示。

表4 不同氫初壓條件下的動力學參數

由表4、表5可以看出，反應速率常數k隨著溫度增加而出現不同程度的增加趨勢，隨著氫初壓的增加，k值略有減少，但幅度不大。從反應活化能來看，在本實驗所提供的反應條件下，柴油→汽油、尾油→汽油以及柴油→氣體的反應活化能較小，比較容易發生，而尾油→柴油、尾油→氣體等反應的活化能較大，比較難發生反應。

對氫初壓為7 MPa、不同溫度及反應時間條件下的估算誤差進行了分析，結果見圖3，其中實線表示計算值，點表示實驗值。

以柴油集總為例，對其在各個反應條件下的估算誤差[估算誤差=(計算值-實驗值)/實驗值×100%]進行計算，結果如表6所示。由表6可知，除了420 ℃、20 min和430 ℃、20 min時的估算誤差大于5%外，其余誤差均小于5%，取其所有誤差絕對值的平均值，得到的平均估算誤差為3.03%。

經計算，其它集總的誤差情況與柴油集總相似，各集總在各個條件下的估算誤差大部分在5%以內，汽油的個別估算誤差大于5%，但小于10%，各集總在所有反應條件下的平均估算誤差分別為氣體3.92%、汽油4.79%、蠟油0.79%、減壓渣油4.25%、焦3.07%，可見估算結果較理想。

對435 ℃、7 MPa、反應時間0，20，40，60 min條件下的汽油、柴油、蠟油收率進行預測，結果如圖4所示。由圖4可知，數據點均勻分布在直線的兩側，汽油、柴油、蠟油在不同反應時間下的ycal/yexp

表5 不同反應溫度下的產物分布動力學參數

圖3 各集總收率的計算值與實驗值對比

表6 柴油集總估算誤差 %

平均值分別為0.980 4，1.004 5，1.005 5，說明預測值和實驗值較吻合，所建立的六集總動力學模型能滿足實際應用中對輕油產品收率的預測要求。

圖4 435 ℃、7 MPa條件下汽油、柴油、蠟油收率的預測值與實驗值對比

5 結 論

(1) 根據油溶性鉬鎳復配催化劑下KAR懸浮床加氫裂化反應規律，按照沸程劃分集總，建立了六集總(氣體、汽油、柴油、蠟油、減壓渣油、焦)動力學模型。

(2) 通過Matlab編程，采用最小二乘法對模型參數進行估算，估算誤差基本小于5%，所建動力學模型的合理性得到驗證。

(3) 運用該模型對其它反應條件下的輕油收率進行預測，預測值與實驗值吻合，說明該模型能滿足產品分布預測需求，對懸浮床加氫裂化工藝開發具有一定的指導意義。

[1] 劉海燕，于建寧，鮑曉軍.世界石油煉制技術現狀及未來發展趨勢[J].過程工程學報，2007，7(1)：176-185

[2] 張德義.含硫含酸原油加工技術進展[J].煉油技術與工程，2012，42(1)：1-2

[3] 方磊，郭金濤，吳顯軍，等.渣油懸浮床加氫研究現狀及發展趨勢[J].化工中間體，2008，4(9)：4-8

[4] 翁惠新，毛信軍.石油煉制過程反應動力學[M].北京：烴加工出版社，1987：4-23

[5] Mosby J F，Buttke R D，Cox J A，et al.Process characterization of expended bed reactor in series[J].Chemical Engineering Science，1986，41(4)：989-995

[6] 楊朝合，徐春明，杜峰.渣油加氫裂化反應的窄餾分集總動力學模型：Ⅰ.加氫裂化反應的六集總動力學模型[J].石油學報(石油加工)，1999，15(5)：44-49

[7] 劉晨光，周家順，闕國和，等.孤島渣油在分散型催化劑存在下加氫裂化反應動力學的研究[J].燃料化學學報，1994，22(1)：28-36

[8] 尚猛.油溶性催化劑和助劑在渣油懸浮床加氫裂化中的研究[D].青島：中國石油大學(華東)，2010

[9] 徐春明，林世雄.渣油加氫裂化反應機理的研究[J].石油煉制與化工，1995，26(9)：52-57

[10]張數義，鄧文安，羅輝，等.渣油懸浮床加氫裂化反應機理[J].石油學報(石油加工)，2009，25(2)：147-148

[11]張數義，羅輝，鄧文安，等.遼河渣油懸浮床加氫裂化反應條件的考察[J].石油化工高等學校學報，2008，21(3)：57-59

STUDY ON LUMPED KINETIC MODEL OF HYDROCRACKING IN SLURRY-BED

Zhang Pingping， Du Feng

(Chemical Engineering School， China University of Petroleum(East China)， Qingdao， Shandong 266580)

Under the conditions of different reaction temperatures and initial hydrogen pressures， the hydrocracking of Karamay atmospheric residue is done in an autoclave in order to simulate the hydrocracking process in slurry-bed. Based on the experimental data and the demand of predicting the yields of light oil products in actual process， a six-lumped kinetic model (gas， gasoline， diesel， gas oil， vacuum residue， coke) of hydrocracking process in slurry-bed is established. The kinetic parameters are determined by the least square method with matlab software and the error is also analyzed. The results show that the model can well predict the yield of each lumped product and the calculation results agree with the experimental value and most of the errors are within 5%.

slurry-bed； hydrocracking； lumped kinetics； parameter estimation； error analysis

2014-12-29； 修改稿收到日期： 2015-03-10。

張萍萍，碩士研究生，研究方向為石油與天然氣加工。

杜峰，E-mail：dufeng@upc.edu.cn。

中央高校基本科研業務費專項資金項目。