甲基環(huán)己烷催化裂化分子水平反應(yīng)動(dòng)力學(xué)模型的建立

張 旭，郭錦標(biāo)，周 祥，王鑫磊，于 博，付 翁

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

當(dāng)前，環(huán)保法規(guī)不但要求減少空氣污染物排放量，而且對(duì)燃油產(chǎn)品的各項(xiàng)質(zhì)量指標(biāo)提出了新的要求。例如，歐Ⅴ標(biāo)準(zhǔn)規(guī)定清潔汽油和柴油的硫質(zhì)量分?jǐn)?shù)上限值為10μg/g，對(duì)苯和整體芳烴含量也作出了明確的限制［1］，表明新的環(huán)保法規(guī)對(duì)車(chē)用燃油的質(zhì)量標(biāo)準(zhǔn)提出了更高的要求。另外，油品重質(zhì)化劣質(zhì)化加重的趨勢(shì)也是不爭(zhēng)的事實(shí)。目前主要是將重質(zhì)油加氫處理后再進(jìn)行催化裂化處理，生產(chǎn)清潔汽、柴油，這一加工過(guò)程首先是將芳烴變成環(huán)烷烴，再得到清潔汽、柴油。煉油工業(yè)經(jīng)歷了150多年的發(fā)展歷史，現(xiàn)在較難實(shí)現(xiàn)重大的技術(shù)創(chuàng)新。為了更高效地利用石油資源，應(yīng)將對(duì)石油及其轉(zhuǎn)化規(guī)律的認(rèn)識(shí)上升到分子水平。這就要求動(dòng)力學(xué)模型的發(fā)展要跟上這個(gè)新的要求，對(duì)復(fù)雜反應(yīng)體系中產(chǎn)物分布及其性質(zhì)作出更細(xì)致、準(zhǔn)確的預(yù)測(cè)，以反映出其與原料性質(zhì)、操作變量、催化劑之間的內(nèi)在關(guān)聯(lián)。

催化裂化反應(yīng)動(dòng)力學(xué)模型經(jīng)歷了關(guān)聯(lián)模型［2］、集總動(dòng)力學(xué)模型［3－5］到組分劃分更為精細(xì)的分子水平動(dòng)力學(xué)模型［6－9］的轉(zhuǎn)變。關(guān)聯(lián)模型在早期得到了廣泛的應(yīng)用，集總動(dòng)力學(xué)模型雖然仍是目前指導(dǎo)催化裂化過(guò)程的最主要?jiǎng)恿W(xué)方法，但是沒(méi)有賦予充分的化學(xué)意義，預(yù)測(cè)能力有限，并且無(wú)法提供詳細(xì)的石油分子組成及結(jié)構(gòu)信息。過(guò)去的二十年，分子水平動(dòng)力學(xué)模型開(kāi)發(fā)得到了快速的發(fā)展，該模型使人們能夠深入認(rèn)識(shí)反應(yīng)機(jī)理，定量預(yù)測(cè)前所未有的分子細(xì)節(jié)。

在20世紀(jì)80年代就有相關(guān)模型化合物的分子水平反應(yīng)動(dòng)力學(xué)模型研究報(bào)道，如烷烴［9－10］、烷基苯［11］的分子水平動(dòng)力學(xué)機(jī)理模型。由于烴類(lèi)化合物的催化裂化反應(yīng)涉及到大量化學(xué)反應(yīng)，構(gòu)成了一個(gè)復(fù)雜的反應(yīng)網(wǎng)絡(luò)，這也使得建立催化裂化動(dòng)力學(xué)模型變得十分困難。先后出現(xiàn)了結(jié)構(gòu)化模型（Structural model，SM）、單事件模型（Single event kinetic model，SEKM）、結(jié)構(gòu)導(dǎo)向型集總模型（Structure oriented lumping，SOL）、機(jī)理模型等。每一種模型都有其適用范圍，在適用范圍內(nèi)具有較好的預(yù)測(cè)能力。然而，部分模型外推性較差，有的模型需要較深厚的理論基礎(chǔ)，建立過(guò)程十分復(fù)雜，對(duì)建模人員的理論知識(shí)提出了更高的要求，限制了其運(yùn)用。因此，需要開(kāi)發(fā)一種操作性強(qiáng)、計(jì)算量小、適用范圍廣的催化裂化動(dòng)力學(xué)模型。

為此，筆者從單環(huán)環(huán)烷烴入手，以甲基環(huán)己烷作為模型化合物，探索建立基于反應(yīng)路徑層面的單環(huán)環(huán)烷烴催化裂化分子水平反應(yīng)動(dòng)力學(xué)模型，以期獲取一些關(guān)于單環(huán)環(huán)烷烴的定量動(dòng)力學(xué)參數(shù)，為建立較完善的單環(huán)環(huán)烷烴催化裂化分子水平反應(yīng)動(dòng)力學(xué)模型奠定基礎(chǔ)。

1 甲基環(huán)己烷催化裂化反應(yīng)原料和實(shí)驗(yàn)數(shù)據(jù)處理方法

1.1 原料和催化劑

甲基環(huán)己烷（MCHA），質(zhì)量分?jǐn)?shù)99%以上，Alfa Aesar公司產(chǎn)品。2種商業(yè)USY分子篩購(gòu)于Tosoh公司。

取質(zhì)量分?jǐn)?shù)為50%高嶺土、30%USY分子篩以及20%硅溶膠經(jīng)噴霧干燥、成型后，再通過(guò)焙燒、老化處理，得到實(shí)驗(yàn)用催化劑，按催化劑平均孔徑大小分別表示為 CAT-LC、CAT-SC，其中CAT-LC的孔徑較大。

1.2 實(shí)驗(yàn)數(shù)據(jù)處理方法

采用等溫固定床反應(yīng)器進(jìn)行催化劑活性評(píng)價(jià)實(shí)驗(yàn)［12］。根據(jù)建立分子水平動(dòng)力學(xué)模型的要求，在本研究中，采用摩爾流率（mol/s）來(lái)表示體系中各組分的濃度。（1）根據(jù)原料中各組分的質(zhì)量分?jǐn)?shù)和密度計(jì)算原料的平均密度；（2）根據(jù)原料的質(zhì)量流率計(jì)算產(chǎn)物的總質(zhì)量流率；（3）根據(jù)總質(zhì)量流率和質(zhì)量分?jǐn)?shù)計(jì)算各組分的質(zhì)量流率；（4）根據(jù)各組分的質(zhì)量流率計(jì)算各組分的摩爾流率。

2 建立甲基環(huán)己烷催化裂化反應(yīng)動(dòng)力學(xué)模型的基本步驟

2.1 產(chǎn)物分子分類(lèi)及命名

雖然MCHA催化裂化反應(yīng)產(chǎn)物的組成很復(fù)雜，但基本上可以將其分為烷烴、烯烴、環(huán)烷烴、芳烴4種基本化合物類(lèi)型。采用簡(jiǎn)單的烷烴、異構(gòu)烷烴、烯烴、環(huán)烷烴和芳烴（PIONA）分類(lèi)來(lái)表示數(shù)量龐大的反應(yīng)物、中間體和產(chǎn)物，不能滿足分子水平動(dòng)力學(xué)模型建立的需要。為了方便后期的反應(yīng)動(dòng)力學(xué)網(wǎng)絡(luò)書(shū)寫(xiě)，采取了系統(tǒng)命名法的簡(jiǎn)寫(xiě)形式，如Methane簡(jiǎn)寫(xiě)為MTA，具體寫(xiě)法見(jiàn)表1。

2.2 動(dòng)力學(xué)參數(shù)獲取

2.2.1 指前因子

根據(jù)過(guò)渡態(tài)理論、統(tǒng)計(jì)熱力學(xué)和阿倫尼烏斯方程，可得到反應(yīng)速率常數(shù)和指前因子的計(jì)算式 （1）和（2）。從過(guò)渡態(tài)理論出發(fā)，經(jīng)過(guò)一定的假設(shè)和引入矯正因子，用氣相分子的平動(dòng)熵近似代替難以獲得的反應(yīng)活化態(tài)的活化熵與反應(yīng)物熵差，大大簡(jiǎn)化了工作量，在一定程度上也保證了數(shù)據(jù)的準(zhǔn)確度［13］。

2.2.2 活化能

采用式（2）可以從理論的角度計(jì)算出所有反應(yīng)的指前因子值，而Ea可以通過(guò)量子化學(xué)計(jì)算軟件—Material Studio（MS）軟件的DMol3模塊計(jì)算獲得［14］。雖然可以通過(guò)模擬計(jì)算獲得所有反應(yīng)的Ea，但是催化裂化反應(yīng)數(shù)目眾多，其費(fèi)時(shí)較多而不具有現(xiàn)實(shí)可操作性。為此，運(yùn)用反應(yīng)族的概念來(lái)估算每一個(gè)基元反應(yīng)的速率常數(shù)。線性自由能關(guān)系能夠起到熱力學(xué)和動(dòng)力學(xué)的橋梁作用，Polanyi方程就是其中的一種，將反應(yīng)的焓變和活化能關(guān)聯(lián)，如式（3）所示。杜梅西克等［15］已經(jīng)證實(shí)Polanyi方程適合于酸催化的正碳離子化學(xué)反應(yīng)。因此每一個(gè)反應(yīng)就可以通過(guò)最多3個(gè)參數(shù)A、E0、α進(jìn)行描述。

每一個(gè)反應(yīng)的焓變就是產(chǎn)物和反應(yīng)物之間的焓差，如式（4）所示。對(duì)于催化裂化而言，由于反應(yīng)過(guò)程中存在許多正碳離子中間體，不能直接通過(guò)實(shí)驗(yàn)方式獲得，然而量化計(jì)算為反應(yīng)焓變計(jì)算提供了條件。目前，基于半經(jīng)驗(yàn)的VAMP量化計(jì)算可以獲得反應(yīng)的生成焓變。將Arrhenius方程和Polanyi方程關(guān)聯(lián)得到的反應(yīng)速率常數(shù)表達(dá)式如式（5）所示。

2.3 模型建立與求解

對(duì)前面的MCHA催化裂化反應(yīng)類(lèi)型進(jìn)行歸類(lèi)整理，進(jìn)一步細(xì)化反應(yīng)類(lèi)型?？偟膩?lái)說(shuō)，正碳離子主要經(jīng)歷了離子引發(fā)、離子轉(zhuǎn)移（變化）、離子消耗3大環(huán)節(jié)，此外反應(yīng)過(guò)程中有可能結(jié)焦，因此構(gòu)成了一個(gè)閉合的反應(yīng)網(wǎng)絡(luò)。

以A作為反應(yīng)物分子，B、C為產(chǎn)物分子進(jìn)行反應(yīng)網(wǎng)絡(luò)的編寫(xiě)。反應(yīng)網(wǎng)絡(luò)主要基于路徑層面，同時(shí)也兼顧機(jī)理層面。動(dòng)力學(xué)模型建立過(guò)程中主要涉及的方程為式（6）～式（9），它們分別是基元反應(yīng)速率、反應(yīng)速率常數(shù)、組分濃度及每一類(lèi)反應(yīng)的速率常數(shù)的表達(dá)式。

最后，使用對(duì)比圖的形式列出了組分摩爾流率的預(yù)測(cè)數(shù)據(jù)與實(shí)驗(yàn)數(shù)據(jù)對(duì)比，可以清晰地看出模型預(yù)測(cè)結(jié)果的準(zhǔn)確度，并可得到優(yōu)化的動(dòng)力學(xué)參數(shù)。如果擬合的結(jié)果不理想，可以進(jìn)一步調(diào)整算法參數(shù)。

3 甲基環(huán)己烷催化裂化反應(yīng)動(dòng)力學(xué)模型的建立

3.1 反應(yīng)數(shù)據(jù)處理

MCHA催化裂化反應(yīng)實(shí)驗(yàn)結(jié)果表明，MCHA的總轉(zhuǎn)化率隨著反應(yīng)溫度的升高和反應(yīng)時(shí)間的延長(zhǎng)而提高。在給定的反應(yīng)時(shí)間內(nèi)，MCHA的轉(zhuǎn)化率會(huì)隨著反應(yīng)溫度的升高而穩(wěn)步上升。此外，MCHA的總轉(zhuǎn)化率不受分子篩晶粒內(nèi)孔擴(kuò)散阻力影響。在所有反應(yīng)條件下，采用CAT-LC催化劑比采用CAT-SC催化劑得到的MCHA轉(zhuǎn)化率稍高。

采用GC-MS分析MCHA催化裂化反應(yīng)產(chǎn)物，檢測(cè)到了60多種產(chǎn)物。對(duì)產(chǎn)物數(shù)據(jù)進(jìn)行處理時(shí)遵循的原則是，烯烴組分不考慮順?lè)串悩?gòu)；烷烴組分只區(qū)分正構(gòu)、異構(gòu)，不考慮支鏈個(gè)數(shù)；環(huán)烷烴和芳烴只含1個(gè)側(cè)鏈；對(duì)于不確定結(jié)構(gòu)的烴類(lèi)按照異構(gòu)烴分子處理。根據(jù)這些原則，按照產(chǎn)物選擇性大小進(jìn)行排序，取選擇性較高的產(chǎn)物作為動(dòng)力學(xué)模型考慮的組分。將色譜分析數(shù)據(jù)歸并處理得到33種主要產(chǎn)物。表1列出了CAT-LC和CAT-SC 2種催化劑催化下，在反應(yīng)溫度550℃、停留時(shí)間7s時(shí)，MCHA催化裂化反應(yīng)主要產(chǎn)物的摩爾流率。

續(xù)表1

3.2 反應(yīng)機(jī)理分析

MCHA催化裂化反應(yīng)主要是由氫轉(zhuǎn)移、異構(gòu)化、質(zhì)子化開(kāi)環(huán)裂化、β－斷裂、脫烷基和烷基轉(zhuǎn)移6個(gè)主要反應(yīng)類(lèi)型所組成的復(fù)雜的串聯(lián)和平行反應(yīng)［16］。MCHA分子中包含5個(gè)仲碳原子、1個(gè)伯碳原子和1個(gè)叔碳原子，可以發(fā)生開(kāi)環(huán)反應(yīng)，生成烷烴和烯烴；以H＋進(jìn)攻叔碳原子形成叔正碳離子為主，進(jìn)一步發(fā)生β－斷裂和開(kāi)環(huán)反應(yīng)等。在酸性催化劑催化下，環(huán)烷烴的六元碳環(huán)經(jīng)過(guò)質(zhì)子化環(huán)丙烷異構(gòu)化，發(fā)生縮環(huán)反應(yīng)，生成帶有1個(gè)烷基取代的五元環(huán)烷烴。另外，六元環(huán)具有較強(qiáng)的脫氫趨勢(shì)，可形成苯、甲苯等芳烴分子，芳烴進(jìn)一步脫氫縮合形成焦炭［17－18］。

3.3 反應(yīng)規(guī)則的制定

反應(yīng)規(guī)則不僅是一個(gè)結(jié)構(gòu)能否發(fā)生某一反應(yīng)或者該結(jié)構(gòu)發(fā)生各種反應(yīng)順序的判據(jù)，也是判斷生成產(chǎn)物分子所有反應(yīng)的依據(jù)。它是對(duì)可能發(fā)生反應(yīng)的概括和歸納，不同的工藝過(guò)程需要針對(duì)性的制定不同反應(yīng)規(guī)則。

MCHA作為典型的單環(huán)環(huán)烷烴，根據(jù)催化裂化反應(yīng)化學(xué)知識(shí)，針對(duì)其催化裂化反應(yīng)制定了6大類(lèi)共計(jì)20條反應(yīng)規(guī)則。如，MCHA發(fā)生裂化反應(yīng)不生成甲基、乙基正碳離子，不會(huì)同時(shí)生成多個(gè)正碳離子，不生成二烯烴，取代基都為烷烴，不區(qū)分環(huán)烷環(huán)上取代基是直鏈還是支鏈；考慮正構(gòu)烯烴向異構(gòu)烯烴轉(zhuǎn)化，正構(gòu)烷烴向異構(gòu)烷烴轉(zhuǎn)化，六元環(huán)向五元環(huán)的縮環(huán)異構(gòu)的異構(gòu)化反應(yīng)；只考慮飽和環(huán)的開(kāi)環(huán)反應(yīng)；烯烴和環(huán)烷烴可以發(fā)生氫轉(zhuǎn)移反應(yīng)，生成芳烴和相應(yīng)的烷烴，環(huán)狀結(jié)構(gòu)的氫轉(zhuǎn)移反應(yīng)都以生成具有芳香性的穩(wěn)定結(jié)構(gòu)結(jié)束；環(huán)烷烴（有取代基的環(huán)烷烴）可以發(fā)生兩分子之間的烷基轉(zhuǎn)移和取代烷基的脫烷基反應(yīng)。

3.4 動(dòng)力學(xué)參數(shù)

過(guò)渡態(tài)理論結(jié)合統(tǒng)計(jì)熱力學(xué)可以求取MCHA催化裂化反應(yīng)的指前因子。筆者采取氣相分子平動(dòng)熵變值近似代替反應(yīng)過(guò)程中活化熵與反應(yīng)物熵的差值，引入一個(gè)校正因子θ，簡(jiǎn)化了求解過(guò)程，得到了MCHA催化裂化反應(yīng)中的指前因子數(shù)值［13］。在計(jì)算過(guò)程中，采用了一些歸并措施［19］，得到的指前因子數(shù)值與Watson等［20］的結(jié)果相吻合。不同類(lèi)別反應(yīng)的指前因子數(shù)值相差較大。

3.5 反應(yīng)網(wǎng)絡(luò)

以催化裂化反應(yīng)化學(xué)為指導(dǎo)，結(jié)合制定的反應(yīng)規(guī)則，對(duì)MCHA催化裂化反應(yīng)建立了74個(gè)基于反應(yīng)路徑層面分子水平化學(xué)反應(yīng)組成的反應(yīng)網(wǎng)絡(luò)。為了使反應(yīng)網(wǎng)絡(luò)具有封閉性，增加了7種重要的中間體分子，即反應(yīng)網(wǎng)絡(luò)中的分子數(shù)達(dá)到了40種。表2給出了各反應(yīng)族的反應(yīng)數(shù)目。與機(jī)理模型相比，反應(yīng)路徑層面反應(yīng)數(shù)目明顯減少，最多的反應(yīng)類(lèi)型是氫轉(zhuǎn)移反應(yīng)，這也與實(shí)際的環(huán)烷烴催化裂化反應(yīng)規(guī)律相一致。

表2 MCHA催化裂化反應(yīng)各反應(yīng)族反應(yīng)數(shù)目匯總Table 2 A summary of each reaction number of MCHA catalytic cracking

3.6 動(dòng)力學(xué)模型建立與求解

當(dāng)獲得完整的動(dòng)力學(xué)參數(shù)初值，搭建好反應(yīng)網(wǎng)絡(luò)后，運(yùn)用Levenberg-Marquardt（LM）算法對(duì)不同催化劑體系下實(shí)驗(yàn)得到的產(chǎn)物分布數(shù)據(jù)進(jìn)行回歸計(jì)算，旨在獲得一組具有普適性的動(dòng)力學(xué)參數(shù)值，使建立的模型能夠達(dá)到較好的預(yù)測(cè)效果。首先通過(guò)模型利用動(dòng)力學(xué)參數(shù)初值以一定時(shí)間步長(zhǎng)Δt對(duì)產(chǎn)品分布進(jìn)行模擬計(jì)算（Simulation），隨后對(duì)產(chǎn)物摩爾流率的預(yù)測(cè)值和實(shí)驗(yàn)值進(jìn)行比較，若方差達(dá)到最小要求，則停止計(jì)算，否則調(diào)整動(dòng)力學(xué)參數(shù)值，重復(fù)計(jì)算（Tuning）直到滿足要求為止。

4 MCHA催化裂化反應(yīng)動(dòng)力學(xué)模型預(yù)測(cè)結(jié)果與實(shí)驗(yàn)結(jié)果的比較

圖1為CAT-LC和CAT-SC 2種催化體系下MCHA催化裂化反應(yīng)實(shí)驗(yàn)結(jié)果與動(dòng)力學(xué)模型模擬計(jì)算的對(duì)比。從圖1可以看出，主要反應(yīng)產(chǎn)物如丙烯、異丁烷、甲苯等的摩爾流率的模擬數(shù)據(jù)與實(shí)驗(yàn)數(shù)據(jù)吻合良好，大部分結(jié)果分布于對(duì)角線附近，說(shuō)明建立的模型具有較高的準(zhǔn)確度。計(jì)算得到實(shí)驗(yàn)數(shù)據(jù)和預(yù)測(cè)數(shù)據(jù)的相關(guān)系數(shù)分別為0.92、0.82。雖然與Watson等［20］建立的辛基環(huán)己烷機(jī)理水平催化裂化反應(yīng)動(dòng)力學(xué)模型精度有一定的差距，他們建立的模型預(yù)測(cè)值與實(shí)驗(yàn)值誤差在5%范圍內(nèi)，但是筆者建立的動(dòng)力學(xué)模型不需要大量的理論計(jì)算，并且模型建立過(guò)程較為簡(jiǎn)單，擴(kuò)展性較好，動(dòng)力學(xué)參數(shù)可用于其他催化裂化反應(yīng)體系。MCHA催化裂化反應(yīng)過(guò)程中發(fā)生的化學(xué)反應(yīng)成百上千，筆者僅用40種石油分子建立的74個(gè)化學(xué)反應(yīng)來(lái)表示MCHA催化裂化反應(yīng)過(guò)程。

圖1中有個(gè)別數(shù)據(jù)相差較大，可能是由于在反應(yīng)過(guò)程中忽略了一部分的化學(xué)反應(yīng)，給模擬精度帶來(lái)了一定影響。以MCHA開(kāi)環(huán)裂化生成烯烴的反應(yīng)過(guò)程為例進(jìn)行詳細(xì)分析，該反應(yīng)過(guò)程示于圖2。從圖2可以看出，當(dāng)MCHA生成1個(gè)正碳離子后，一部分發(fā)生環(huán)烷環(huán)縮環(huán)反應(yīng)，一部分則發(fā)生開(kāi)環(huán)反應(yīng)。反應(yīng)生成的烷烴正碳離子以及烯烴型的正碳離子發(fā)生了質(zhì)子化、去質(zhì)子化、異構(gòu)化、甲基遷移、β－斷裂等一系列反應(yīng)，反應(yīng)過(guò)程十分復(fù)雜，并且中間發(fā)生的化學(xué)反應(yīng)大部分都是可逆反應(yīng)，增加了反應(yīng)的復(fù)雜程度，也給分子水平動(dòng)力學(xué)模型建立帶來(lái)了一定的誤差。

圖1 CAT-LC和CAT-SC催化劑作用下MCHA催化裂化反應(yīng)主要產(chǎn)物摩爾流率的預(yù)測(cè)值與實(shí)驗(yàn)值對(duì)比Fig.1 Comparison between predicted and experimental molar flow rates of main products from MCHA catalytic cracking over CAT-LC and CAT-SC catalysts

圖2 MCHA開(kāi)環(huán)裂化生成烯烴的反應(yīng)過(guò)程Fig.2 Reaction pathways of olefin formation from MCHA ring opening

表3是回歸計(jì)算得到的12類(lèi)MCHA催化裂化反應(yīng)動(dòng)力學(xué)參數(shù)值。將從優(yōu)化得到的動(dòng)力學(xué)參數(shù)與先前工作中得到的初值［13］進(jìn)行對(duì)比，發(fā)現(xiàn)兩者之間數(shù)量級(jí)無(wú)大的差別，僅個(gè)別數(shù)據(jù)的差別較大（如脫烷基反應(yīng)），說(shuō)明其初值具有較好的準(zhǔn)確性。對(duì)比2種催化劑E0數(shù)據(jù)可以看出，CAT-LC的催化裂化反應(yīng)活性要高于CAT-SC，β裂化、異構(gòu)化等反應(yīng)活化能差值不大，而對(duì)于氫轉(zhuǎn)移反應(yīng)，2種催化劑差值較大?；罨蹺0值可以反映催化劑的不同［12，21］。對(duì)氫轉(zhuǎn)移反應(yīng)而言，CAT-LC的E0值明顯低于CAT-SC，說(shuō)明CAT-LC更有利于氫轉(zhuǎn)移反應(yīng)的發(fā)生，這主要是因?yàn)镃AT-LC催化劑有較大的孔徑，有利于反應(yīng)物、產(chǎn)物分子的擴(kuò)散，提高雙分子氫轉(zhuǎn)移反應(yīng)速率［12，22］，表1中產(chǎn)物芳烴含量較多也證明了這一點(diǎn)。另外，經(jīng)過(guò)乙基環(huán)己烷、辛基環(huán)己烷等單環(huán)環(huán)烷烴的催化裂化實(shí)驗(yàn)數(shù)據(jù)驗(yàn)證，回歸得到的12類(lèi)動(dòng)力學(xué)參數(shù)值，也可以用于其動(dòng)力學(xué)模型的建立，這也為進(jìn)一步建立單環(huán)環(huán)烷烴催化裂化反應(yīng)分子水平動(dòng)力學(xué)模型提供數(shù)據(jù)支持。

表3 回歸得到的MCHA催化裂化反應(yīng)動(dòng)力學(xué)參數(shù)Table 3 Optimized kinetic parameters for MCHA catalytic cracking

5 結(jié) 論

采用MCHA作為單環(huán)環(huán)烷烴模型化合物和2種孔徑尺寸不同的催化劑進(jìn)行催化裂化反應(yīng)，對(duì)產(chǎn)物進(jìn)行了分析，選取了33種主要產(chǎn)物用于建立分子水平的催化裂化反應(yīng)動(dòng)力學(xué)模型。根據(jù)實(shí)驗(yàn)結(jié)果，簡(jiǎn)要分析了MCHA催化裂化6大反應(yīng)類(lèi)型，通過(guò)反應(yīng)原料和產(chǎn)物的分子簡(jiǎn)化命名管理、制定20條反應(yīng)規(guī)則、建立由74個(gè)化學(xué)反應(yīng)組成的反應(yīng)網(wǎng)絡(luò)、獲取動(dòng)力學(xué)參數(shù)以及模型建立與求解等步驟，建立了MCHA催化裂化基于反應(yīng)路徑層面的分子水平動(dòng)力學(xué)模型。建立的分子水平動(dòng)力學(xué)模型對(duì)主要的反應(yīng)產(chǎn)物摩爾流率的預(yù)測(cè)值與實(shí)驗(yàn)值吻合良好，相關(guān)系數(shù)分別為0.92、0.82?；貧w得到了12類(lèi)甲基環(huán)己烷催化裂化反應(yīng)動(dòng)力學(xué)參數(shù)值，可以用于其他單環(huán)環(huán)烷烴模型化合物動(dòng)力學(xué)模型的建立。

符號(hào)說(shuō)明：

A——指前因子，Pa1－n·s－1；

Ai——i類(lèi)反應(yīng)的指前因子，Pa1－n·s－1；

C——組分摩爾濃度，mol/g；

Cbulk，xi——i物質(zhì)在x類(lèi)反應(yīng)中的摩爾濃度，mol/g；

Cbulk，xj——j物質(zhì)在x類(lèi)反應(yīng)中的摩爾濃度，mol/g；

Ccations，xi——i物質(zhì)在x類(lèi)反應(yīng)中的離子摩爾濃度，mol/g；

Ccations，xj——j物質(zhì)在x類(lèi)反應(yīng)中的離子摩爾濃度，mol/g；

E0——活化能，J/mol；

E0，i——i類(lèi)反應(yīng)的活化能，J/mol；

Ea——反應(yīng)活化能，J/mol；

Ea，x——x類(lèi)反應(yīng)活化能，J/mol；

EA，forwardrxn——正向反應(yīng)活化能，J/mol；

EA，reverserxn——逆向反應(yīng)活化能，J/mol；

h——Planck常數(shù)，6.62×10－34（J·s）；

Hproducts——產(chǎn)物的摩爾生成焓，J/mol；

Hreactants——反應(yīng)物的摩爾生成焓，J/mol；

k——反應(yīng)速率常數(shù)，速率系數(shù)的單位取決于反應(yīng)的總級(jí)數(shù)；

kB——Boltzmann常數(shù)，1.38×10－23（J/K）；

kij——i類(lèi)反應(yīng)中j的反應(yīng)速率，速率系數(shù)的單位取決于反應(yīng)的總級(jí)數(shù)；

kcatalytic，x——x類(lèi)反應(yīng)的速率常數(shù)；

kcatalytic，y——y類(lèi)反應(yīng)的速率常數(shù)；

n——反應(yīng)分子數(shù)；

R——理想氣體常數(shù)，8.314J/（K·mol）；

t——反應(yīng)時(shí)間，s；

T——反應(yīng)溫度，K；

Δn——過(guò)渡態(tài)與反應(yīng)物之間的分子數(shù)變化，單分子反應(yīng)為0，雙分子反應(yīng)為－1；

ΔHrxn，j——j反應(yīng)的反應(yīng)焓變，J/mol；

ΔHrxn——反應(yīng)焓變，J/mol；

ΔS0＊——?dú)庀嗷罨夭睿琂/（mol·K）；

α——反應(yīng)熱與活化能之間的關(guān)聯(lián)系數(shù)，量綱為1；

αi——i類(lèi)反應(yīng)活化能與反應(yīng)熱之間的關(guān)聯(lián)系數(shù)，量綱為1。

［1］STANISLSUS A，MARAFI A，RAN M S.Recent advances in the science and technology of ultra low sulfur diesel（ULSD）production［J］.Catal Today，2010，153（1-2）：1-68.

［2］張結(jié)喜，齊艷華，邱建章.催化裂化關(guān)聯(lián)模型的研究［J］.計(jì)算機(jī)與應(yīng)用化學(xué)，2007，24（11）：1519-1522.（ZHANG Jiexi，QI Yanhua，QIU Jianzhang.Study on FCC correlation models［J］.Comput Appl Chem，2007，24（11）：1519-1522.）

［3］戴鑑，楊光福，王剛，等.催化裂化汽油改質(zhì)降烯烴并多產(chǎn)丙烯的反應(yīng)動(dòng)力學(xué)模型研究［J］.燃料化學(xué)學(xué)報(bào)，2008，36（4）：432-437.（DAI Jian，YANG Guangfu，WANG Gang，et al.Kinetic model of FCC naphtha reformulation for reducing olefin and enhancing propylene production［J］.J Fuel Chem Technol，2008，36（4）：432-437.）

［4］許友好，龔劍波，張久順，等.多產(chǎn)異構(gòu)烷烴的催化裂化工藝兩個(gè)反應(yīng)區(qū)概念實(shí)驗(yàn)研究［J］.石油學(xué)報(bào)（石油加工 ），2004，20 （4）：1-5. （XU Youhao， GONG Jianhong，ZHANG Jiushun，et al.Experimental study on“two reaction zone”conception connected with MIP process［J］.Acta Petrolei Sinica（Petroleum Processing Section），2004，20（4）：1-5.）

［5］王建平，許先焜，翁惠新，等.加氫渣油催化裂化14集總動(dòng)力學(xué)模型的建立［J］.化工學(xué)報(bào)，2007，58（1）：86-94.（WANG Jianping，XU Xiankun，WENG Huixin，et al.Establishment of 14lumps model for fluid catalytic cracking of hydrotreated residuum［J］.J Chem Ind Eng，2007，58（1）：86-94.）

［6］SOTELO-BOYáS R，F(xiàn)ROMENT G F.Fundamental kinetic modeling of catalytic reforming［J］.Ind Eng Chem Res，2009，48（3）：1107-1119.

［7］JAFFE S B，F(xiàn)REUND H，OLMSTEAD W N.Extension of structure-oriented lumping to vacuum residua［J］.Ind Eng Chem Res，2005，44（26）：9840-9852.

［8］WEI W，BENNETT C A， TANAKA R，et al.Computer aided kinetic modeling with KMT and KME［J］.Fuel Process Technol，2008，89（4）：350-363.

［9］WATSON B A，KLEIN M T， HARDING R H.Mechanistic modeling ofn-hexadecane cracking on rare earth Y［J］.Energ Fuels，1997，11（2）：354-363.

［10］WATSON B A，KLEIN M T， HARDING R H.Mechanistic modeling ofn-heptane cracking on HZSM-5［J］.Ind Eng Chem Res，1996，35（5）：1506-1516.

［11］WATSON B A，KLEIN M T， HARDING R H.Mechanistic modeling of a 1-phenyloctane/n-hexadecane mixture on rare earth Y zeolite［J］.Ind Eng Chem Res，1997，36（8）：2954-2963.

［12］AL-SABAWI M N.Heterogeneous kinetic modeling of the catalytic conversion of cycloparaffins［D］.Ontario：The University of Western Ontario，2009.

［13］張旭，郭錦標(biāo)，周祥，等.單環(huán)環(huán)烷烴催化裂化動(dòng)力學(xué)模型的建立——指前因子的計(jì)算［J］.石油學(xué)報(bào)（石油加工），2013，29 （2）：283-288.（ZHANG Xu，GUO Jinbiao，ZHOU Xiang，et al.Kinetic modeling of catalytic cracking of monocyclic cycloparaffins——Calculation of preexponential factors［J］.Acta Petrolei Sinica （Petroleum Processing Section），2013，29（2）：283-288.）

［14］沈本賢，田立達(dá)，劉紀(jì)昌.基于結(jié)構(gòu)導(dǎo)向集總的石腦油蒸汽裂解過(guò)程分子尺度動(dòng)力學(xué)模型［J］.石油學(xué)報(bào)（石油加工），2010，（S1）：218-225.（SHEN Benxian，TIAN Lida，LIU Jichang.A molecular kinetic model for naphtha steam cracking based on structure oriented lumping［J］.Acta Petrolei Sinica（Petroleum Processing Section），2010，（S1）：218-225.）

［15］杜梅西克，拉德，阿帕里西奧，等.多相催化微觀動(dòng)力學(xué)（沈儉一譯）［M］.北京：國(guó)防工業(yè)出版社，1998：25-40.

［16］張旭，周祥，郭錦標(biāo)，等.甲基環(huán)己烷催化裂化的研究進(jìn)展［J］.石油化工，2013，42（1）：104-110.（ZHANG Xu，ZHOU Xiang，GUO Jinbiao，et al.Advances in the development of methylcyclohexane catalytic cracking［J］.Petrochem Technol，2013，42（1）：104-110.）

［17］于珊，張久順，魏曉麗.環(huán)烷烴催化裂解生成乙烯和丙烯反應(yīng)探析［J］.石油學(xué)報(bào)（石油加工），2013，29（3）：475-481.（YU Shan，ZHANG Jiushun，WEI Xiaoli.Exploration and analysis on ethylene and propylene formation in naphthene catalytic cracking［J］.Acta Petrolei Sinica（Petroleum Processing Section），2013，29（3）：475-481.）

［18］YU S，ZHANG J S，WEI X L.Research on ethylene and propylene formation during catalytic pyrolysis of methylcyclohexane［J］.China Pet Process Pe，2012，14（4）：73-79.

［19］LEE J H，KANG S，KIM Y，et al.New approach for kinetic modeling of catalytic cracking of paraffinic naphtha［J］.Ind Eng Chem Res，2011，50（8）：4264-4279.

［20］WATSON B A，KLEIN M T， HARDING R H.Catalytic cracking of alkylcyclohexanes：Modeling the reaction pathways and mechanisms［J］.Int J Chem Kinet，1997，29（7）：545-560.

［21］Al-SABAWI M，DE LASA H.Kinetic modeling of catalytic conversion of methylcyclohexane over USY zeolites：Adsorption and reaction phenomena［J］.AIChE J，2009，55（6）：1538-1558.

［22］唐津蓮，許友好，汪燮卿，等.十氫萘在分子篩催化劑上的開(kāi)環(huán)反應(yīng)研究［J］.燃料化學(xué)學(xué)報(bào)，2012，40（12）：1422-1428. （TANG Jinlian， XU Youhao， WANG Xieqing，et al.Opening of naphthenic ring in decalin cracking over zeolite catalysts［J］.J Fuel Chem Technol，2012，40（12）：1422-1428.）