正戊烷熱裂解自由基反應模型的研究

姬偉毅

（中國石油化工股份有限公司 化工事業部，北京 100728）

正戊烷熱裂解自由基反應模型的研究

姬偉毅

（中國石油化工股份有限公司 化工事業部，北京 100728）

在自由基反應機理的基礎上，構建了正戊烷熱裂解自由基反應網絡，根據量子化學計算結果和文獻數據，建立了正戊烷熱裂解自由基反應模型，該模型包括14個分子和17個自由基共計199個自由基反應。采用NM-Ⅱ型蒸汽熱裂解模擬實驗裝置對正戊烷進行不同溫度條件下的熱裂解實驗，以驗證模型模擬結果的精度。實驗結果表明，在爐管出口溫度低于750 ℃時，正戊烷熱裂解產物收率的模擬結果與實驗結果基本吻合；在爐管出口溫度高于750 ℃時，模擬結果與實驗結果存在偏差；建立的正戊烷熱裂解自由基反應模型能夠較好地模擬低溫工況下裂解產物的收率情況。

自由基反應模型；正戊烷；熱裂解；乙烯

乙烯工業是石油化工的基礎。目前，世界上乙烯主要是通過石油烴類在管式爐中以蒸汽熱裂解方式生產的。我國乙烯工業起步于20世紀60年代，隨著乙烯工業的快速發展，先后引進了30～100 kt/a規模的管式裂解爐。但在裂解反應動力學研究方面，國內比較滯后，目前還僅僅處于探索階段，成熟的機理模型非常少，主要是由實驗數據回歸得到的經驗模型［1-4］。為了滿足裂解爐設計開發的需要，經驗模型正逐步被分子模型所取代［5-6］。

中國石油化工股份有限公司是國內乙烯的主要生產商，2011年乙烯產量首次突破10 Mt。近年來，隨著乙烯原料結構優化、輕重原料分儲分裂、裂解爐改造以及先進過程控制推廣應用等要求，需要建立必要的裂解反應機理模型［7］。因此，構建較大分子烴類熱裂解的自由基反應網絡，完善反應參數的計算方法［8-10］，進而建立起一套完整的熱裂解反應機理模型，用于裂解爐的設計、改造以及操作優化，對于提高我國乙烯裝置的競爭力水平具有重要意義。

本研究在自由基反應機理的基礎上［11］，構建了正戊烷熱裂解自由基反應網絡，通過量子化學方法計算得到94個自由基反應動力學參數，從文獻［12-13］補充碳四及碳四以下烷烴熱裂解自由基反應及其動力學參數，建立了正戊烷熱裂解自由基反應模型，包括14個分子和17個自由基共計199個自由基反應，并通過實驗來驗證模型模擬結果的精度。

1 模型與實驗部分

1.1 正戊烷熱裂解自由基反應模型

烴類熱裂解的機理反應過程主要由鏈引發反應、氫奪取反應、自由基β分解反應、鏈終止反應構成，具體如下：

鏈引發反應：

氫奪取反應：

自由基β分解反應：

鏈終止反應：

根據上述反應類型，盡可能列舉出正戊烷熱裂解時各自由基與分子之間的基元反應，得到由94個自由基反應構成的正戊烷熱裂解反應網絡，反應產物主要為碳五及碳五以下組分。通過量子化學的方法，計算得到94個自由基反應動力學參數。以這94個自由基反應及其動力學參數為基礎，從文獻補充碳四及碳四以下烷烴熱裂解的自由基反應及其動力學參數，構成了共199個自由基反應組成的正戊烷熱裂解自由基反應網絡，反應過程包含14個分子和17個自由基。

采用Chemkin Pro軟件（版本4.1.1）對正戊烷熱裂解過程進行模擬。根據反應器的尺寸參數及溫度和壓力分布、物料流量（投油量和投水量）、化學反應（自由基反應），計算得到正戊烷熱裂解各產物的含量及其收率。

1.2 正戊烷熱裂解實驗

實驗以分析純正戊烷為原料，在中國石化北京化工研究院NM-Ⅱ型蒸汽熱裂解模擬實驗裝置上進行。該實驗裝置（見圖1）主要由1個隔離的水蒸汽發生單元和7個隔離的加熱裂解單元（分別對應T1～T7）組成，每個裂解單元可以分別控制不同的溫度，以模擬不同爐型產生的不同溫度分布。實驗中，前6個隔離加熱裂解單元設定為相同的溫度，將最后一個隔離加熱單元設定為不同的溫度，并通過控制投油量來保持物料在裂解模擬實驗裝置中的停留時間一致。

由于熱裂解反應的主要目的產物是乙烯和丙烯，因此乙烯和丙烯的收率是衡量裂解性能高低的標準。在熱裂解反應過程中，爐管出口溫度是裂解反應最重要的影響因素。本實驗通過研究不同爐管出口溫度對裂解反應的影響，從而驗證模型對正戊烷熱裂解反應模擬結果的精度。

圖1 NM-Ⅱ型蒸汽熱裂解模擬實驗裝置示意Fig. 1 Schematic diagram of the NM-Ⅱ type experimental equipment for the steam pyrolysis ofn-pentane.

正戊烷熱裂解實驗條件見表1。

表1 正戊烷熱裂解實驗的工藝條件Table 1 The experimental conditions for the steam pyrolysis ofn-pentane

按照表1給定的各項工藝參數設定模擬實驗裝置的控制方案，先開啟蒸餾水的控制閥，并開啟各段的加熱裝置，使模擬實驗裝置的各段溫度達到設定值后，開啟原料的控制閥開始投油，待各段溫度穩定在設定值后開始實驗。在實驗過程中取兩次氣相產物試樣，分析其成分和含量。待實驗結束后停止投油，并稱量冷凝的液相中水和焦油的含量，將分離出的焦油作為液相試樣，測定其成分和含量。通過反應前后的物料平衡，計算得到正戊烷熱裂解后的原料轉化率和各產物的產量。

2 結果與討論

2.1 正戊烷熱裂解的模擬結果

采用本研究建立的正戊烷熱裂解自由基反應模型，對應表1實驗條件對正戊烷熱裂解產物收率進行模擬，模擬結果見表2。

表2 不同實驗條件下正戊烷熱裂解各產物收率的模擬結果Table 2 The simulation results for the steam pyrolysis ofn-pentane under different experimental conditions obtained by the pyrolysis radical reaction model

2.2 正戊烷熱裂解的實驗結果

正戊烷熱裂解的實驗結果見圖2，顯示了正戊烷轉化率、主要產物收率隨爐管出口溫度（tcot）的變化趨勢。

圖2 正戊烷熱裂解實驗各產物收率隨tcot的變化Fig. 2 The changes of various product yields with coil outlet temperature(tcot) in the the steam pyrolysis ofn-pentane.

從圖2可看出，正戊烷熱裂解反應的主要產物為乙烯和丙烯，其次為甲烷；正戊烷的轉化率隨tcot的升高而增大，當tcot低于650 ℃時，正戊烷的轉化率低于5%；隨tcot的升高，乙烯收率不斷增大， 且極值出現在實驗設計的最高tcot處；丙烯收率呈現先增大后減小的趨勢，極值出現在850 ℃；甲烷收率和乙烯收率一樣，隨tcot的升高而增大，極值也出現在實驗中的最高tcot處；其他產物的收率均低于5%。

2.3 模擬結果與實驗結果的比較

將正戊烷熱裂解模擬結果和實驗結果對tcot做圖，結果見圖3。從圖3a可看出，正戊烷轉化率的模擬結果與實驗結果在750 ℃以下基本吻合，但在750 ℃以上模擬結果高于實驗結果。從圖3b可看出，在低于750 ℃時，甲烷、乙烯和丙烯收率的模擬結果與實驗結果差別不大（低于5%），但高于750℃后收率差別變大；乙烯和甲烷收率的模擬結果與實驗結果的變化規律基本一致，而丙烯收率的實驗結果出現最大值，但模擬結果卻未見最大值。出現這種情況是因為正戊烷在tcot高于750 ℃時，溫度越高，苯、甲苯、二甲苯等大分子產物生成量會更多，而本研究暫未考慮生成這些大分子的反應，導致模型在高溫時計算的產物收率出現偏差。

圖3 正戊烷熱裂解的模擬結果與實驗結果的比較Fig.3 Comparison between the calculation results and experimental data.

3 結論

（1）根據正戊烷的自由基反應網絡和量子化學計算得到的反應動力學參數，從文獻補充碳四及碳四以下烷烴的自由基反應網絡和參數，采用Chemkin Pro軟件模擬正戊烷的熱裂解反應過程，建立了正戊烷熱裂解自由基反應模型。

（2）建立的正戊烷熱裂解自由基反應模型主要考慮了碳五及碳五以下裂解產物，當tcot在750 ℃以下時，該模型模擬得到的正戊烷轉化率及甲烷、乙烯、丙烯收率與實驗結果差別不大，能較好地模擬低溫工況下裂解產物的收率情況；當tcot高于750 ℃時，由于生成苯等大分子的反應影響較大，使得模擬結果與實驗結果存在偏差。因此，在建立正戊烷熱裂解自由基反應模型的過程中，需要進一步考慮生成苯等大分子產物的反應，以提高更高tcot條件下裂解產物收率的預測精度。

［1］ 楊元一，曾清泉，許士興，等. 重質油裂解數學模型——裂解反應動力學模型［J］. 石油化工，1984，13（11）：703 -710.

［2］ 莊曉峰，俞金壽. 裂解深度的建模及應用［J］. 自動化儀表，2004，25（6）：31 - 35.

［3］ 王國清，杜志國，張利軍，等. 應用BP神經網絡預測石腦油熱裂解產物收率［J］. 石油化工，2007，36（7）：699 - 704.

［4］ 趙佩清，顏學峰. 基于新型蟻群算法優化的重油熱裂解模型［J］. 化工自動化及儀表，2007，34（6）：20 - 23.

［5］ Hirato M，Yoshioka S. Simulations of the Pyrolysis of Naphtha，Kerosene，and Gas Oil with a Tubular Reactor［J］.Int Chem Eng，1973，13（2）：347 - 354.

［6］ Kumar P，Kunzru D. Modeling of Naphtha Pyrolysis［J］.Ind Eng Chem Proc Des Dev，1985，24（3）：774 - 782.

［7］ 李蔚，張兆斌，周叢，等. 裂解爐管內自由基反應模型的研究進展［J］. 乙烯工業，2010，22（2）：1 - 6.

［8］ 李華，張兆斌. 熱裂解鏈引發-終止反應的計算——C—H鍵的斷裂-形成［J］. 石油化工，2006，35（7）：643 - 648.

［9］ Li Hua，Chen Bozhen，Huang Mingbao. CASPT2 Investigation of Ethane Dissociation and Methyl Recombination Using Canonical Variational Transi tion State Theory［J］.Int J Chem Kinet，2008，40（4）：161 - 173.

［10］ 李蔚，王國清，杜志國，等. 石油烴熱裂解反應速率常數估算方法的研究進展［J］. 乙烯工業，2009，21（2）：1 - 7.

［11］ 張兆斌，李華，張永剛，等. 丁烷熱裂解自由基反應模型的建立和驗證［J］. 石油化工，2007，36（1）：44 - 48.

［12］ Sundaram K M，Froment G F. Modeling of Thermal Cracking Kinetics：Ⅰ. Thermal Cracking of Ethane，Propane and Their Mixtures［J］.Chem Eng Sci，1977，32（6）：601 - 608.

［13］ Sundaram K M，Froment G F. Modeling of Thermal Cracking Kinetics：Ⅱ. Cracking ofiso-Butane，ofn-Butane，of Mixtures Ethane-Propane-n-Butane［J］.Chem Eng Sci，1977，32（6）：609 - 617.

（編輯 李治泉）

Radical Reaction Model for n-Pentane Pyrolysis

Ji Weiyi
（Chemical Department，China Petroleum & Chemical Corporation，Beijing 100728，China）

Based on radical reaction mechanism，a model for n-pentane steam pyrolysis was established，which included 14 molecules，17 free radicals and 199 free radical reactions in total，according to the result of quantum chemistry calculation and data from literatures. The experiments of the n-pentane steam pyrolysis were carried out in a NM-II type cracking simulation device under different conditions to validate the precision of the simulation results. The results indicated that the radical reaction model could simulate the n-pentane pyrolysis process at low temperature(below 750 ℃)，and there was some difference between the results from the experiments and from the simulation at higher temperature(above 750 ℃).

radical reaction model；n-pentane；pyrolysis；ethylene

book=27,ebook=27

1000-8144（2012）06 - 0633 - 04

TQ 203.8

A

2012 - 01 - 18；［修改稿日期］2012 - 03 - 24。

姬偉毅（1972—），男，陜西省富平縣人，大學，高級工程師，電話 010 - 59969552，電郵 jiweiyi@sinopec.com。