熊 凱， 盧春喜

(中國石油大學(xué) 重質(zhì)油國家重點實驗室， 北京 102249)

催化裂化(裂解)集總反應(yīng)動力學(xué)模型研究進(jìn)展

熊 凱， 盧春喜

(中國石油大學(xué) 重質(zhì)油國家重點實驗室， 北京 102249)

分別從催化裂化集總反應(yīng)動力學(xué)模型、催化裂解集總反應(yīng)動力學(xué)模型、集總反應(yīng)-流動耦合的催化裂化模型、集總反應(yīng)-催化劑性質(zhì)耦合的催化裂化模型、催化劑失活模型、建模數(shù)據(jù)來源、參數(shù)估算方法和集總動力學(xué)模型的關(guān)聯(lián)模型等方面，系統(tǒng)地介紹了催化裂化(裂解)集總反應(yīng)動力學(xué)模型的研究進(jìn)展。在分析總結(jié)的基礎(chǔ)上，針對催化裂化(裂解)集總反應(yīng)動力學(xué)模型研究的后續(xù)發(fā)展提出了幾點建議。

流化催化裂化(FCC)；催化裂解；集總模型；動力學(xué)；流體力學(xué)；催化劑性質(zhì)

在整個石油煉制的行業(yè)布局中，催化裂化一直以其良好的原料適應(yīng)性、高附加值產(chǎn)品的產(chǎn)出和較溫和的操作條件的優(yōu)點而成為原油加工的一個重要環(huán)節(jié)，是重質(zhì)進(jìn)料輕質(zhì)化的主要手段之一。據(jù)報道[1]，我國催化裂化裝置總加工能力已接近150 Mt/a，其中渣油占催化裂化進(jìn)料約40%，提供了國內(nèi)市場約70%的汽油、約30%的柴油和40%的丙烯。近年來，隨著原油重質(zhì)、劣質(zhì)化趨勢的不斷加劇，以及市場對更多清潔燃料和低碳烯烴的需求，催化裂化承受的壓力也在不斷加大。為了應(yīng)對這些新變化，催化裂化裝置急需優(yōu)化其原有的操作模式，改變操作條件，甚至進(jìn)行裝置和工藝改造。要想解決這些問題，歸根到底就是要弄清催化裂化反應(yīng)過程中數(shù)量龐大而復(fù)雜的化學(xué)反應(yīng)。目前，針對這一問題的研究，最成熟且效果較好的方法就是進(jìn)行集總動力學(xué)建模研究。

催化裂化反應(yīng)中涉及的反應(yīng)種類繁多，且多為平行-順序反應(yīng)，各反應(yīng)間偶聯(lián)性很強。20世紀(jì)60年代，Aris[2]與Wei等[3]提出了集總的方法，使催化裂化中組分繁多、多重體系的動力學(xué)研究有了新突破。所謂集總，即是將反應(yīng)系統(tǒng)中眾多的單一化合物，按動力學(xué)特性相似的原則，歸并為若干虛擬組分，稱為集總；而在動力學(xué)研究中，則把每個集總作為一個虛擬的單一組分，然后開發(fā)這些虛擬的集總組分的反應(yīng)網(wǎng)絡(luò)，建立簡化的集總反應(yīng)網(wǎng)絡(luò)的動力學(xué)模型。這種集總方法，首先在煉油催化加工過程的動力學(xué)研究中得到了成功的應(yīng)用，并在指導(dǎo)工業(yè)裝置的操作、設(shè)計、先進(jìn)控制和催化劑的選型方面成效顯著，為煉油廠創(chuàng)造了巨大的經(jīng)濟效益。

在市場對清潔燃料和低碳烯烴的需求逐年加大的背景下，集總動力學(xué)模型的研究顯得十分重要。筆者綜述了催化裂化(裂解)集總反應(yīng)動力學(xué)模型的發(fā)展歷程和研究現(xiàn)狀，并在此基礎(chǔ)上展望催化裂化(裂解)集總反應(yīng)動力學(xué)模型的后續(xù)發(fā)展，以期促進(jìn)該領(lǐng)域的研究工作。

1 催化裂化集總反應(yīng)動力學(xué)模型的發(fā)展

從催化裂化集總反應(yīng)動力學(xué)模型的發(fā)展歷程看，可以依次分為蠟油催化裂化集總動力學(xué)模型階段、重油催化裂化集總動力學(xué)模型階段和基于新工藝的催化裂化集總動力學(xué)模型階段。

1.1 蠟油催化裂化集總反應(yīng)動力學(xué)模型

20世紀(jì)60年代末，Mobil公司的Weekman等[4]基于Aris和韋潛光的集總理論，率先開發(fā)了蠟油催化裂化三集總模型。三集總模型將整個催化裂化反應(yīng)體系的所有組分劃分為原料油、汽油、氣體+焦炭3個集總，同時，還考慮了因催化劑表面生焦所引起的催化劑失活。利用三集總模型，可以較好地預(yù)測同一原料在不同反應(yīng)條件下的轉(zhuǎn)化率和產(chǎn)品分布，預(yù)測該原料生成汽油的選擇性，并根據(jù)生產(chǎn)需求選擇對應(yīng)的操作模式，實現(xiàn)優(yōu)化操作。但是，由于三集總模型在定義原料油集總時忽略了化學(xué)組成變化的影響，因而當(dāng)原料油組成有較大變化時，特別是采取回?zé)挷僮鲿r，模型的適用性較差。三集總模型是最早的針對蠟油為原料的催化裂化集總動力學(xué)模型，雖然有很大局限性，卻為處理催化裂化復(fù)雜反應(yīng)體系指明了方向，為催化裂化集總動力學(xué)模型的后續(xù)發(fā)展奠定了堅實的基礎(chǔ)。

Lee等[5]提出了蠟油催化裂化四集總模型。蠟油四集總模型是在三集總模型的基礎(chǔ)上，將氣體和焦炭分別集總，使煉油廠在生產(chǎn)過程中能夠分別監(jiān)測氣體和焦炭的生成，便于依據(jù)兩器熱平衡調(diào)節(jié)催化裂化裝置的操作，優(yōu)化產(chǎn)品分布。蠟油四集總模型因其簡潔的優(yōu)勢，后來成為在流動-集總反應(yīng)耦合催化裂化模型的開發(fā)中使用最多的集總反應(yīng)動力學(xué)模型。Ancheyta-Juarez等[6]提出了蠟油催化裂化五集總模型。蠟油五集總模型是在四集總模型的基礎(chǔ)上，將氣體集總細(xì)分為干氣和液化氣2個集總，以便在催化裂化反應(yīng)過程中能夠分別監(jiān)測干氣和液化氣的生成。Kraemer等[7]提出了蠟油催化裂化八集總模型，將原料油按餾程分為輕燃料油(LFO)(220～345℃)和重燃料油(HFO)(>345℃)2部分，然后將LFO和HFO分別細(xì)分為烷烴(P)、環(huán)烷烴(N)和芳烴(A)3個集總，再加上汽油和氣體+焦炭，共8個集總。任杰等[8]按照我國煉油廠對油品餾程的劃分習(xí)慣，對原料油和汽油集總餾程范圍重新進(jìn)行了劃分，在他們建立的蠟油八集總模型中，將原料油按餾程分為輕餾分油(LGO)(204～350℃)和重餾分油(HGO)(350～500℃)2部分。在五集總模型的基礎(chǔ)上，Hagelberg等[9]提出了一種新的蠟油催化裂化八集總模型。在此模型中，將汽油集總按照PONA組成細(xì)分為烷烴、烯烴、環(huán)烷烴和芳烴4個集總，而原料不細(xì)分。任杰等的蠟油八集總模型為原料集總細(xì)分模型，因而原料適應(yīng)性較強；Hagelberg等的蠟油八集總模型為產(chǎn)品細(xì)分模型，不僅能預(yù)測產(chǎn)品分布，還能預(yù)測汽油的質(zhì)量，但是原料適應(yīng)性不夠好。

Jacob等[10]基于Weekman催化裂化三集總模型開發(fā)了蠟油催化裂化十集總模型。將原料按結(jié)構(gòu)和餾程劃分為8個集總，再加上汽油、裂化氣+焦炭2個產(chǎn)物集總，使得整個反應(yīng)網(wǎng)絡(luò)劃分為10個集總。具體來說，是將原料油按餾程分為LFO(221～343℃)和HFO(>343℃)，LFO和HFO各自再分別細(xì)分為烷烴(P)、環(huán)烷烴(N)、單環(huán)芳烴(A)和多核芳烴(CA)4個集總，共8個集總。此外，還考慮了原料中堿性氮中毒、重芳烴吸附對催化劑活性的影響和催化劑的時變失活。與三集總模型相比，此模型最大的優(yōu)點在于不受原料組成變化的限制，原料適用性更強。針對我國催化裂化工藝的特點，朱炳辰等[11]在Jacob等開發(fā)的蠟油催化裂化十集總模型基礎(chǔ)上，提出了蠟油催化裂化十一集總模型。由于國外催化裂化裝置一般均采用小回?zé)挶然驘o回?zé)挷僮?，而我國催化裂化裝置普遍采用較大的回?zé)挶炔僮?，回?zé)捰椭泻写罅康碾y裂化的重芳烴，特別是短側(cè)鏈的多環(huán)芳烴，經(jīng)實驗驗證這些芳烴與新鮮原料油中所含芳烴的裂化性能具有很大差異。因此，在十一集總模型中，將HFO層中的芳核集總(CAh)按照環(huán)數(shù)細(xì)分為一、二環(huán)芳烴(CAh)與多環(huán)芳烴(PCAh)2個集總。催化裂化十一集總模型適應(yīng)我國催化裂化裝置較大回?zé)挶炔僮鞯膰?，成功地指?dǎo)了新工藝優(yōu)化、新裝置設(shè)計和最佳操作條件的選擇等。劉之定等[12]在蠟油催化裂化十一集總模型的基礎(chǔ)上，將氣體和焦炭分別集總，提出了蠟油催化裂化十二集總模型。

Araujo-Monroy等[13]提出了蠟油催化裂化十五集總模型。模型將原料油分為烷烴、環(huán)烷烴和芳烴3個集總，輕循環(huán)油(LCO)和汽油均按照其PONA組成各自分為4個集總，液化氣(LPG)分為烷烴和烯烴2個集總，再加上干氣和焦炭2個集總，總共15個集總。催化裂化十五集總模型的原料適應(yīng)性強，且能預(yù)測汽、柴油及液化氣的產(chǎn)品質(zhì)量。Pitault等[14]提出了蠟油催化裂化十九集總模型。模型將原料油分為烷烴、芳烴、1～3環(huán)環(huán)烷烴和3環(huán)以上環(huán)烷烴4個集總，LCO和汽油均分別按照烷烴、鏈烯烴、環(huán)烯烴、環(huán)烷烴和芳烴各自分為5個集總，LPG和干氣均按照烷烴和烯烴各自分為2個集總，再加上焦炭集總，共19個集總。通過對十九集總模型的研究，發(fā)現(xiàn)縮合反應(yīng)和氫轉(zhuǎn)移反應(yīng)與焦炭的生成及汽油的質(zhì)量有著重要的聯(lián)系。由于十九集總對于催化裂化反應(yīng)組分的描述較之于簡單集總模型細(xì)致了許多，Pitault等認(rèn)為，十九集總模型已介于簡單集總動力學(xué)模型與分子級別動力學(xué)模型之間。

以上綜述表明，集總劃分得越細(xì)致越能揭示催化裂化反應(yīng)的深層反應(yīng)規(guī)律。同時，越細(xì)致的集總劃分也將會給模型參數(shù)的求取帶來更大的難度。

1.2 重油催化裂化集總反應(yīng)動力學(xué)模型

20世紀(jì)80年代末，為了使催化裂化反應(yīng)動力學(xué)模型適用于渣油原料，Takatsuka等[15]在Weekman蠟油三集總模型的基礎(chǔ)上，開發(fā)了重油催化裂化六集總模型。與蠟油裂化相比，渣油裂化反應(yīng)過程的二次反應(yīng)和結(jié)焦反應(yīng)更劇烈。二次反應(yīng)對汽油的產(chǎn)率影響很大，而在提升管出口的返混現(xiàn)象對二次反應(yīng)的影響很顯著，因此，他們采用了多段全混流的提升管反應(yīng)器模型。此外，在反應(yīng)集總劃分方面，按照餾程將反應(yīng)體系劃分為渣油/油漿、餾分油/回?zé)捰?、柴油、汽油、裂化氣和焦?個集總，以適應(yīng)渣油催化裂化的情況。

牟盛靜等[16]考慮工業(yè)現(xiàn)場和在線計算的需要，基于提升管反應(yīng)器內(nèi)活塞流假設(shè)，提出了渣油催化裂化六集總模型。其集總劃分與Takatsuka等的渣油六集總模型相同。Xu等[17]進(jìn)一步將裂化氣集總細(xì)分為液化氣LPG和干氣2個集總，提出了渣油催化裂化七集總模型。Gao等[18]提出了重油催化裂化八集總模型，將反應(yīng)網(wǎng)絡(luò)分為原料油烷烴、環(huán)烷烴和芳烴3個集總，重循環(huán)油、輕質(zhì)油、液化氣、干氣和焦炭，共8個集總。

任杰[19]提出了重油催化裂化十一集總模型。將反應(yīng)網(wǎng)絡(luò)分為渣油反應(yīng)、重燃料油反應(yīng)、輕燃料油反應(yīng)和汽油反應(yīng)4層，含25個速率常數(shù)；渣油、重燃料油和輕燃料油均按照烷烴、環(huán)烷烴和芳烴分別細(xì)分為3個集總，再加上汽油和氣體+焦炭2個集總，共11個集總。鄧先梁等[20]和沙穎遜等[21]進(jìn)一步將氣體和焦炭分別集總，并將渣油層中的芳烴集總細(xì)分為芳烴和膠質(zhì)+瀝青質(zhì)2個集總，提出了重油催化裂化十三集總模型。孫鐵棟等[22]提出了渣油催化裂化十四集總模型。將回?zé)捰蛦为毤?，將原料分為渣油、蠟油和重循環(huán)油3層，并按照SARA組成將各層分別細(xì)分為飽和分、芳香分和膠質(zhì)+瀝青質(zhì)3個集總，再加上LCO、汽油、LPG、干氣和焦炭，共14個集總。

此外，國內(nèi)外許多服務(wù)于重油催化裂化裝置的軟件中都有各自的重油集總動力學(xué)模型。如中國石化洛陽石油化工工程公司的FCCM軟件、華東理工大學(xué)的RCCLK軟件、Aspen Tech的Aspen HYSYS軟件、UOP的UniSim軟件、KBC的FCC-SIM軟件、Shell的SHARC軟件、SimSci的PRO/II軟件等。Aspen HYSYS軟件的石油煉制軟件包中包含重油催化裂化二十一集總動力學(xué)模型[23]，21個集總分別為渣油(烷烴、環(huán)烷烴、芳烴側(cè)鏈、單環(huán)芳烴、雙環(huán)芳烴和三環(huán)芳烴)、蠟油(烷烴、環(huán)烷烴、芳烴側(cè)鏈、單環(huán)芳烴、雙環(huán)芳烴和三環(huán)芳烴)、柴油(烷烴、環(huán)烷烴、芳烴側(cè)鏈、單環(huán)芳烴和雙環(huán)芳烴)、汽油、裂化氣和焦炭(動力學(xué)焦和金屬焦)。很明顯，這些商業(yè)軟件的集總動力學(xué)模型能更細(xì)致地描述和調(diào)控催化裂化的生產(chǎn)過程。

1.3 基于新工藝的催化裂化集總反應(yīng)動力學(xué)模型

新時期，特別是2000年以后，催化裂化的加工工藝得到了豐富，產(chǎn)品分布也出現(xiàn)了一定程度的改變。因此，有必要將這部分內(nèi)容分離出來，單獨介紹。

王建平等[24]報道了加氫渣油催化裂化十四集總模型。他們直接以工業(yè)數(shù)據(jù)為建模數(shù)據(jù)，結(jié)合合適的分析手段和煉油廠對模型應(yīng)用方面的要求，提出了按照渣油SARA四組分劃分原料集總的方法。該模型的原料為加氫改質(zhì)后的渣油，按SARA四組分劃分為飽和分、芳香分、膠質(zhì)+瀝青質(zhì)3個集總；同時，考慮到清潔燃料生產(chǎn)對輕質(zhì)油品的要求以及煉油廠的實際需要，產(chǎn)品柴油和汽油分別劃分為烷烴、環(huán)烷烴、烯烴和芳烴4個集總，液化石油氣劃分為烷烴和烯烴2個集總，再加上氣體+焦炭，共14個集總；并通過分步法求取了動力學(xué)參數(shù)。通過工業(yè)數(shù)據(jù)驗證表明，該模型較好地反映了加氫渣油催化裂化的反應(yīng)規(guī)律，不僅能預(yù)測催化裂化產(chǎn)品分布，還能預(yù)測主要產(chǎn)品性質(zhì)，可為工業(yè)裝置操作的優(yōu)化提供指導(dǎo)。

吳飛躍等[25]基于靈活多效催化裂化(FDFCC)工藝，分別對汽油提升管和重油提升管中的反應(yīng)體系建立了汽油九集總模型和重油十二集總模型。汽油九集總模型充分考慮了汽油提升管中的二次反應(yīng)，重油十二集總模型的重油原料按照SARA族組成進(jìn)行了集總劃分。后續(xù)分析了兩提升管間的偶聯(lián)關(guān)系，并用數(shù)學(xué)方法將2個子反應(yīng)體系結(jié)合起來，最終建立了FDFCC工藝的集總反應(yīng)動力學(xué)模型。該模型除了能預(yù)測主產(chǎn)物的收率外，還能預(yù)測汽油中烯烴和芳烴以及液化氣中丙烯和丁烯的收率。該模型的成功建立證明了集總動力學(xué)模型的方法適于描述多反應(yīng)器偶聯(lián)的復(fù)雜反應(yīng)體系。

江洪波等[26]基于多產(chǎn)異構(gòu)烷烴(MIP)工藝，開發(fā)了重油催化裂化十二集總模型。模型中，原料和柴油均采用PONA結(jié)構(gòu)族組成進(jìn)行集總劃分，并對反應(yīng)路徑進(jìn)行了合理簡化。該模型除了能預(yù)測主產(chǎn)物的收率外，還能預(yù)測汽油中烯烴和芳烴以及丙烯的收率。

段良偉等[27]基于多產(chǎn)清潔汽油并聯(lián)產(chǎn)丙烯(MIP-CGP)工藝的特點，開發(fā)了重油催化裂化十集總模型。將原料按SARA四組分劃分，以工業(yè)數(shù)據(jù)為基礎(chǔ)，建立并驗證了模型。對主要產(chǎn)品汽油及汽油烯烴含量的預(yù)測值相對偏差在5%以內(nèi)，能滿足工業(yè)生產(chǎn)需要，并能為裝置的操作優(yōu)化提供指導(dǎo)。

劉熠斌等[28]基于兩段提升管催化裂化(TSRFCC)工藝，分別對一段提升管和二段提升管中的反應(yīng)體系建立了一段提升管重油六集總模型和二段提升管重油六集總模型；通過循環(huán)油物料衡算將2個子反應(yīng)體系結(jié)合起來，最終建立了TSRFCC工藝的集總反應(yīng)動力學(xué)模型。該模型能分別預(yù)測一段和二段提升管的轉(zhuǎn)化率和各產(chǎn)物的產(chǎn)率，提供操作優(yōu)化信息。

此外，還有一些比較特別的產(chǎn)物集總劃分。如Villafuerte-Macias等[29]報道了一個獨具特色的蠟油催化裂化七集總模型。為了研究催化裂化反應(yīng)中的硫平衡，以使原料中的硫元素盡可能多地轉(zhuǎn)化為酸性氣體H2S，而非有機硫化物，對H2S進(jìn)行了單獨集總。該模型能預(yù)測硫元素在產(chǎn)物中的分布情況，并能反過來促進(jìn)操作優(yōu)化，盡可能避免硫元素進(jìn)入汽油和焦炭產(chǎn)品中。

2 催化裂解集總反應(yīng)動力學(xué)模型的發(fā)展

從催化裂解集總反應(yīng)動力學(xué)模型的研究對象看，可以分為蠟油和重油催化裂解集總反應(yīng)動力學(xué)模型、FCC汽油改質(zhì)和輕油催化裂解集總反應(yīng)動力學(xué)模型2大類。

2.1 蠟油和重油催化裂解集總反應(yīng)動力學(xué)模型

隨著石化行業(yè)的蓬勃發(fā)展，市場對乙烯、丙烯、丁烯等低碳烯烴的需求逐年增加。為了多產(chǎn)低碳烯烴，中國石化石油化工科學(xué)研究院開發(fā)了深度催化裂解(DCC)工藝和催化熱裂解(CPP)系列催化裂解工藝，中國石化洛陽石油化工工程公司研究院開發(fā)了重油接觸裂解(HCC)工藝，中國石油大學(xué)開發(fā)了重油兩段提升管催化裂解多產(chǎn)丙烯(TMP)工藝。為了深入研究催化裂解的反應(yīng)機理，指導(dǎo)現(xiàn)有催化裂解工業(yè)裝置的操作和優(yōu)化催化裂解工藝設(shè)計，針對這些工藝的催化裂解反應(yīng)動力學(xué)模型多有報道。

許友好[30]在DCC深度裂解工藝基礎(chǔ)上，開發(fā)了DCC蠟油催化裂解四集總模型。以中型催化裂解試驗數(shù)據(jù)為建模數(shù)據(jù)，以汽油和輕烯烴為主要產(chǎn)物，將反應(yīng)體系劃分為VGO原料、汽油、LPG和干氣+焦炭4個集總，建立了反應(yīng)網(wǎng)絡(luò)，確定了動力學(xué)參數(shù)及催化劑失活方程。研究表明，模型能預(yù)測不同原料、不同操作條件下的催化裂解產(chǎn)品分布。張琪皓[31]使用下行床實驗裝置研究了丙烷脫瀝青油在DCC工藝催化劑CHP-1催化下的裂解性能，將丙烷脫瀝青油、汽油、丙烯、丁烯、丙烷+丁烷和干氣+焦炭分別集總，建立了催化裂解六集總模型。模型假定生成干氣和焦炭的反應(yīng)遵循自由基反應(yīng)機理，其他反應(yīng)則遵循正碳離子機理，并且根據(jù)不同反應(yīng)使用了不同的催化劑失活函數(shù)。研究表明，該模型的處理方法獲得了較大成功。

陳新國[32]基于CPP工藝，以大慶常壓渣油為原料，建立了常壓渣油催化裂解五集總動力學(xué)模型。模型的5個集總分別為重油、汽柴油、烯烴、烷烴+氫氣和焦炭。孟祥海等[33]在HCC工藝的基礎(chǔ)上，參照陳新國的CPP五集總模型，建立了大慶常壓渣油催化裂解HCC五集總動力學(xué)模型。以輕油和輕烯烴為主要目的產(chǎn)物，集總的劃分同CPP五集總模型。研究表明，HCC五集總模型對催化裂解輕油、氣體烷烴和氣體烯烴有較好的預(yù)測能力，短停留時間有利于多產(chǎn)低碳烯烴。在CPP五集總的基礎(chǔ)上，Meng等[34-35]又先后2次對集總數(shù)目進(jìn)行了擴大，分別建立了七集總和八集總模型。將烯烴集總細(xì)分為乙烯和丙烯+丁烯2個集總，將烷烴+氫氣集總細(xì)分為丙烷+丁烷和甲烷+乙烷+氫氣2個集總，成為七集總模型；再將重油原料細(xì)分為芳香碳和非芳香碳2個集總，成為八集總模型。這些嘗試都取得了較好的效果，對產(chǎn)品分布進(jìn)行了越來越細(xì)致的描述。由于大慶常壓渣油屬于石蠟基，為了研究不同質(zhì)量進(jìn)料的催化裂解反應(yīng)行為差異，Meng等[36]研究了環(huán)烷基屬的源自加拿大合成原油的重蠟油催化裂解行為。研究中使用催化裂解七集總模型，7個集總分別為重蠟油、汽油、丙烯+丁烯、丙烷+丁烷、乙烯、氫氣+甲烷+乙烷和焦炭。研究表明，原料集總比其他中間產(chǎn)物集總有更大的反應(yīng)速率常數(shù)和更低的表觀活化能。

王國良等[37]基于重油接觸裂解HCC工藝，建立了HCC重油催化裂解十六集總動力學(xué)模型。按餾程將原料劃分為重油層(H)、柴油層(L)(200～300℃)、汽油(GO)(C5+～200℃)、裂化氣(GAS)和焦炭(CK)，然后用PONA結(jié)構(gòu)族組成分別表示重油層和柴油層，裂化氣細(xì)分為氫氣、甲烷、乙烯、乙烷、丙烯、丙烷、丁烯和丁烷，共16個集總。模型體現(xiàn)了反應(yīng)溫度、壓力、劑/油質(zhì)量比、水/油質(zhì)量比、催化劑活性等因素的影響，各產(chǎn)物產(chǎn)率的計算值和實驗值絕對偏差一般在0.5%左右，模型精度和預(yù)測能力均較好。

Liu等[38]將大慶常壓渣油催化裂解反應(yīng)體系重新進(jìn)行了集總劃分，基于兩段提升管催化裂解多產(chǎn)丙烯的TMP工藝，提出了大慶常渣催化裂解八集總模型。模型的集總為重油(>350℃)、柴油(205～350℃)、汽油(IBP～205℃)、C4+丙烷、丙烯、乙烯、焦炭、氫氣+甲烷+乙烷，共8個集總。研究表明，產(chǎn)物產(chǎn)率的計算值和實驗值一致。杜峰等[39]在Liu等的八集總模型基礎(chǔ)上，將乙烷單獨集總，建立了大慶常壓渣油催化裂解九集總模型。研究表明，模型對柴油、汽油、乙烯、丙烯等主產(chǎn)物的產(chǎn)率預(yù)測效果較好，預(yù)測值與實驗值的相對誤差在4%以內(nèi)。郭菊花[40]在Liu等的八集總模型基礎(chǔ)上，將汽油細(xì)分為汽油烯烴、汽油芳烴和汽油飽和烴3個集總，并將丁烯單獨集總，建立了TMP催化裂解十一集總模型。研究表明，主要產(chǎn)品的計算值與實驗值的平均相對誤差在10%以內(nèi)。

2.2 催化裂化汽油改質(zhì)和輕油催化裂解集總反應(yīng)動力學(xué)模型

劉熠斌等[41]基于兩段提升管催化裂化TSRFCC工藝，開發(fā)了FCC汽油催化裂解多產(chǎn)低碳烯烴的六集總模型。模型的集總為柴油、汽油、液化氣(除丙烯外)、乙烯+丙烯、干氣(除乙烯外)和焦炭，共6個集總。研究表明，該模型具有較好的預(yù)測能力，計算值和實驗值吻合良好。

韓忠祥等[42]為了研究汽油二次裂解制低碳烯烴的反應(yīng)，建立了FCC汽油裂解八集總模型。他們認(rèn)為，F(xiàn)CC汽油二次反應(yīng)中只有少量柴油生成，且這部分柴油主要為芳烴，可將其歸到汽油芳烴集總中；此外，氫氣和焦炭的生成量非常少，建模中可以忽略。因而模型集總分為汽油烷烴、汽油烯烴、汽油環(huán)烷烴、汽油芳烴、丁烷+丁烯、丙烯、乙烯、甲烷+乙烷+丙烷，共8個集總。研究表明，模型的計算值和實驗值具有一致性。

吳飛躍等[43]基于靈活多效催化裂化(FDFCC)工藝，對汽油提升管中的反應(yīng)體系建立了汽油改質(zhì)九集總模型。模型的集總為柴油、汽油飽和烴、汽油烯烴、汽油芳烴、液化氣烷烴、丙烯、丁烯、干氣和焦炭，共9個集總。該模型充分考慮了汽油提升管中的二次反應(yīng)，能預(yù)測汽油的組成和丙烯產(chǎn)率，并能為汽油改質(zhì)提供操作優(yōu)化。

Wang等[44]為了研究FCC汽油二次轉(zhuǎn)化改質(zhì)的反應(yīng)規(guī)律，將汽油按照PONA組成劃分集總，開發(fā)了FCC汽油改質(zhì)八集總動力學(xué)模型。模型的集總為輕循環(huán)油、汽油鏈烷烴、汽油環(huán)烷烴、汽油烯烴、汽油芳烴、液化氣、干氣和焦炭，共8個集總。研究表明，產(chǎn)物產(chǎn)率的計算值和實驗值一致。

劉福安等[45]應(yīng)用經(jīng)典正碳離子反應(yīng)機理，提出FCC汽油改質(zhì)的十集總動力學(xué)模型。將汽油按照其PONA組成進(jìn)行了集總劃分，分為柴油、汽油正構(gòu)烷烴、汽油異構(gòu)烷烴、汽油環(huán)烷烴、汽油烯烴、汽油芳烴、液化氣烷烴、液化氣烯烴、干氣和焦炭，共10個集總。研究表明，模型對不同F(xiàn)CC汽油和操作條件具有較好的適應(yīng)性，產(chǎn)物產(chǎn)率的計算值和實驗值相對誤差為7%，能指導(dǎo)FCC汽油的改質(zhì)生產(chǎn)。

吳青等[46]按照PONA組成對汽油進(jìn)行了集總劃分，在固定流化床上研究了FCC汽油改質(zhì)過程中的反應(yīng)規(guī)律，提出了FCC汽油改質(zhì)八集總動力學(xué)模型。模型的集總分為柴油、汽油正構(gòu)烷烴、汽油異構(gòu)烷烴、汽油烯烴、汽油環(huán)烷烴、汽油芳烴、裂化氣和焦炭。采用該模型定量計算了汽油催化改質(zhì)過程中氫轉(zhuǎn)移、異構(gòu)化、環(huán)化、齊聚等二次反應(yīng)以及兩類氫轉(zhuǎn)移之間的比例，計算值和實驗值的相對誤差大多在2%以內(nèi)。吳青等[47]又以中型提升管數(shù)據(jù)為基礎(chǔ)，針對汽油組成中變化最大的烯烴和異構(gòu)烷烴的詳細(xì)碳數(shù)分布，構(gòu)建了更為詳細(xì)的反應(yīng)網(wǎng)絡(luò)，開發(fā)了FCC汽油改質(zhì)二十集總動力學(xué)模型。在八集總的基礎(chǔ)上，異構(gòu)烷烴和烯烴分別按照碳數(shù)細(xì)分為C3+C4、C5、C6、C7、C8、C9和C10+7個集總，共20個集總。模型的計算值與實驗值的相對誤差在6%以內(nèi)。

楊光福等[48]在催化裂化汽油輔助反應(yīng)器改質(zhì)降烯烴技術(shù)基礎(chǔ)上，建立了六集總反應(yīng)動力學(xué)模型。模型的集總分為柴油、汽油飽和烴、汽油烯烴、汽油芳烴、液化氣和干氣+焦炭，共6個集總。研究表明，該模型對原料和反應(yīng)條件變化具有較好的適應(yīng)性，能預(yù)測不同條件下產(chǎn)物的產(chǎn)率和汽油組成，有助于降低汽油烯烴含量的研究。為了在汽油降烯烴的同時聯(lián)產(chǎn)丙烯，戴鑑等[49]開發(fā)了相應(yīng)的六集總動力學(xué)模型。模型的集總分為汽油飽和烴、汽油烯烴、汽油芳烴、液化氣(除丙烯外)、丙烯和干氣+焦炭，共6個集總。結(jié)果表明，模型能預(yù)測不同條件下的汽油組成和丙烯產(chǎn)率。

鄒圣武等[50]在研究FCC汽油二次轉(zhuǎn)化過程中，將汽油和餾分油均按照PONA組成進(jìn)行集總劃分，提出了基于汽油烯烴轉(zhuǎn)化的十二集總模型。模型的集總分為餾分油(烷烴、烯烴、環(huán)烷烴和芳烴)、汽油(烷烴、烯烴、環(huán)烷烴和芳烴)、液化氣烷烴、液化氣烯烴、干氣和焦炭，共12個集總。研究表明，該模型能預(yù)測不同條件下的產(chǎn)物產(chǎn)率分布，包括汽油的詳細(xì)組成，有利于降低汽油烯烴含量的研究。

3 集總反應(yīng)-流動耦合的催化裂化模型的發(fā)展

集總反應(yīng)動力學(xué)模型一般只考察化學(xué)反應(yīng)，而忽略提升管反應(yīng)器內(nèi)流動、傳熱和傳質(zhì)的影響。但是，實際生產(chǎn)中提升管內(nèi)的傳熱、傳質(zhì)和氣固流動狀態(tài)必然對裂化反應(yīng)產(chǎn)生影響。催化裂化數(shù)學(xué)模型的關(guān)鍵和難點就在于精確地描述催化劑顆粒的流動狀態(tài)。顆粒流動模型的研究已有大量報道，特別是應(yīng)用計算流體力學(xué)CFD模擬軟件以后，這方面的研究發(fā)展迅速。

高金森等[51]借助重油催化裂化十三集總動力學(xué)模型和數(shù)值模擬的方法，系統(tǒng)地考察了提升管反應(yīng)器中的氣-固兩相流動、傳質(zhì)、傳熱及反應(yīng)等復(fù)雜因素及相互影響，建立了催化裂化提升管三維氣-固兩相流動-反應(yīng)模型。在后續(xù)研究中，高金森等[52]又將原料油液霧的流動、傳熱、氣化過程及其對反應(yīng)結(jié)果的影響全納入模型中，對催化裂化提升管反應(yīng)器進(jìn)行了三維微分模擬，考察了原料液霧流動氣化過程及其影響，建立了提升管氣-液-固三相流動-反應(yīng)模型。陳新國等[53]借助重油催化裂化十三集總動力學(xué)模型，結(jié)合催化劑顆粒流動方程、油氣湍流方程，耦合得到了催化裂化提升管反應(yīng)器數(shù)學(xué)模型。該耦合模型的建立有助于研究提升管內(nèi)的反應(yīng)情況。徐春明等[54]在重油催化裂化十三集總反應(yīng)動力學(xué)模型和流動模型的基礎(chǔ)上，綜合考慮流動、傳熱、反應(yīng)，建立了渣油催化裂化提升管反應(yīng)器的數(shù)學(xué)模型，通過對4套不同工況進(jìn)行數(shù)值模擬，考察了反應(yīng)器內(nèi)氣相組成分布和催化劑顆粒的溫度、速度分布，得到了與工業(yè)裝置基本一致的提升管出口參數(shù)，驗證了模型的可靠性。Wang等[55]借助蠟油催化裂化四集總模型和CFD模擬軟件，考察了FCC提升管中催化劑顆粒聚團現(xiàn)象對催化裂化反應(yīng)的影響。模擬結(jié)果表明，聚團顆粒與單個顆粒相比，由蠟油生成汽油、氣體和焦炭的反應(yīng)速率均降低，而由汽油生成氣體和焦炭的反應(yīng)速率則升高。催化劑顆粒聚團會導(dǎo)致反應(yīng)速率降低，氣體和汽油產(chǎn)率下降。

侯栓弟等[56]結(jié)合重油催化裂化十四集總動力學(xué)模型，考察了提升管內(nèi)反應(yīng)溫度、局部固體濃度變化以及流動對催化裂化反應(yīng)的影響，建立了重油催化裂化流動-反應(yīng)耦合模型?；谠摿鲃?反應(yīng)耦合模型，郭湘波等[57]進(jìn)行了模擬研究。結(jié)果表明，重油催化裂化反應(yīng)主要發(fā)生在噴嘴附近區(qū)域，在這一區(qū)域約45%的重油轉(zhuǎn)化為汽油和柴油。此外，針對FCC催化汽油改質(zhì)工藝，侯栓弟[58]報道了FCC汽油改質(zhì)反應(yīng)-流動耦合模型。針對不同結(jié)構(gòu)的反應(yīng)器對汽油改質(zhì)反應(yīng)和流動過程進(jìn)行了數(shù)值模擬，結(jié)果表明，對提升管反應(yīng)器，F(xiàn)CC汽油經(jīng)低溫改質(zhì)后，烯烴含量降低了約一半；對于提升管-床層反應(yīng)器，汽油中的烯烴質(zhì)量分?jǐn)?shù)能降至約5%。

Chang等[59]基于CFD軟件和FDFCC工藝的重油催化裂化十二集總動力學(xué)模型，綜合考察了FDFCC工藝重油提升管內(nèi)油劑低溫接觸、大劑/油質(zhì)量比操作條件下的進(jìn)料氣化、流體力學(xué)行為、傳熱和裂化反應(yīng)。模擬結(jié)果表明，提升管內(nèi)的流體力學(xué)行為、溫度和各物種的濃度均呈現(xiàn)巨大的不均一性，尤其在進(jìn)料霧化區(qū)。此后，他們基于提升管底部安裝氣升環(huán)流混合器的新結(jié)構(gòu)的FDFCC工藝重油提升管反應(yīng)器，并繼續(xù)借助十二集總動力學(xué)模型和CFD模擬軟件，考察了不同操作條件對進(jìn)料氣化、流動和反應(yīng)的影響，建立了重油提升管三相流動-反應(yīng)耦合催化裂化模型[60]。模擬結(jié)果表明，提升管底部安裝氣升環(huán)流混合器以后，冷、熱催化劑的混合更均勻，且使得提升管入口處的催化劑溫度更均一。此外還發(fā)現(xiàn)，噴嘴角度因素大于進(jìn)料液滴尺寸和反應(yīng)溫度對產(chǎn)品收率的影響。為了研究FDFCC工藝汽油提升管中的流體力學(xué)行為、傳熱和裂化反應(yīng)，Chang等[61]借助CFD模擬軟件和FDFCC工藝汽油催化裂化九集總動力學(xué)模型進(jìn)行了模擬研究，建立了汽油提升管流動-傳熱-反應(yīng)耦合催化裂化模型。借助該模型，得到了FCC汽油提升管最佳操作條件。

Li等[62]借助重油催化裂化十四集總動力學(xué)模型，考察了提升管反應(yīng)器的進(jìn)料情況對反應(yīng)器油-劑混合段的流體流動行為和催化裂化反應(yīng)的影響。模擬結(jié)果表明，進(jìn)料噴嘴出口速率決定了油-劑混合區(qū)的氣-固兩相流動情況，噴嘴的角度會對提升管內(nèi)油-劑混合區(qū)的流場和裂化反應(yīng)產(chǎn)生巨大影響。Chen等[63]借助蠟油催化裂化四集總動力學(xué)模型，重點考察了FCC提升管反應(yīng)器中催化劑的顆粒內(nèi)傳質(zhì)和傳熱行為，在綜合考慮質(zhì)量平衡、熱平衡、動量平衡、氣態(tài)方程、集總反應(yīng)動力學(xué)方程、多組分?jǐn)U散和對流傳熱的基礎(chǔ)上，提出了催化裂化的單顆粒模型。研究表明，沿著FCC提升管軸向存在3種典型的顆?，F(xiàn)象，并在傳質(zhì)、傳熱和反應(yīng)并存的情況下形成了傳遞控制區(qū)、反應(yīng)控制區(qū)和中間過渡區(qū)3種不同的催化和反應(yīng)操作區(qū)。該單顆粒模型能基于提升管反應(yīng)器中催化劑顆粒的位置來預(yù)測溫度、壓力、各物種質(zhì)量組成分布、反應(yīng)速率和顆粒內(nèi)有效擴散系數(shù)。Chen等[64]在后續(xù)研究中又借助多尺度CFD方法考察了單顆粒模型和宏觀氣、固兩相流間的關(guān)系，建立了催化裂化多尺度模型。模擬研究表明，多尺度模型能同時預(yù)測催化劑顆粒內(nèi)和提升管反應(yīng)器內(nèi)的流場分布，并能預(yù)測催化劑顆粒的粒徑分布。與傳統(tǒng)的CFD模型不同的是，多尺度模型的研究還發(fā)現(xiàn)，顆粒動力學(xué)(破碎和聚團等)對提升管反應(yīng)器內(nèi)的流場會產(chǎn)生明顯的影響。

Derouin等[65]結(jié)合蠟油催化裂化十九集總動力學(xué)模型，研究了提升管中的流體力學(xué)行為對催化裂化反應(yīng)的影響。研究表明，流動模型和反應(yīng)模型的結(jié)合能較好地預(yù)測產(chǎn)品分布，在研究流體力學(xué)行為對產(chǎn)品分布和產(chǎn)品質(zhì)量的影響方面意義較大。Das等[66]結(jié)合蠟油催化裂化十二集總動力學(xué)模型，考察了提升管中的三維流動狀態(tài)對催化裂化反應(yīng)的影響。結(jié)果表明，與反應(yīng)模型結(jié)合時，三維流動模型比一維流動模型能更準(zhǔn)確地預(yù)測產(chǎn)品分布。Lopes等[67]借助蠟油催化裂化四集總動力學(xué)模型，也對提升管中的流動和反應(yīng)行為進(jìn)行了三維模擬，并取得了良好效果。Zhu等[68]將蠟油催化裂化四集總動力學(xué)模型與提升管中的流體流動行為進(jìn)行耦合分析，考察了提升管中流體流動情況對催化裂化反應(yīng)的影響。由于催化劑的濃度從提升管底部到頂部逐漸降低，計算結(jié)果表明，提升管底部的裂化強度要遠(yuǎn)高于傳統(tǒng)的提升管模型中的情況。與傳統(tǒng)的提升管反應(yīng)模型比較，新開發(fā)的流動-反應(yīng)耦合模型能更好地預(yù)測工業(yè)生產(chǎn)中的產(chǎn)品分布情況。

以上綜述表明，集總反應(yīng)動力學(xué)模型與流動模型的耦合應(yīng)用能更好地反映催化裂化反應(yīng)器中的流動-反應(yīng)情況，并能更準(zhǔn)確地預(yù)測產(chǎn)品分布。

4 集總反應(yīng)-催化劑性質(zhì)耦合的催化裂化模型的發(fā)展

集總反應(yīng)動力學(xué)模型在建模過程中?；谀撤N特定的催化劑，而不考慮不同催化劑對模型動力學(xué)參數(shù)的影響，即將催化劑因素視為集總動力學(xué)模型內(nèi)部隱含的一個恒量。顯然，不同種類的催化劑常因為功能各異而對催化裂化的產(chǎn)品分布造成較大影響。在這種情況下，可以認(rèn)為已建立的集總動力學(xué)模型只適用于與之對應(yīng)的某類特定催化劑，因而針對不同的催化劑需要分別單獨地建立模型，即建模時需要做到“一劑一模型”。在涉及多種催化劑的情況下，如催化劑的對比篩選，常規(guī)的集總動力學(xué)建模方法將費時費力。如何使一個集總動力學(xué)模型同時適用于多種催化劑是目前面臨的一個難題。事實上，不同種類催化劑的根本區(qū)別在于其結(jié)構(gòu)組成和理化性質(zhì)上的不同。也就是說，如何開發(fā)一個包含催化劑結(jié)構(gòu)組成和理化性質(zhì)參數(shù)的集總動力學(xué)模型是問題的所在。目前，切中該問題實質(zhì)的研究報道較少，但亦不乏具有啟發(fā)性的工作。

De Lasa等[69]借助蠟油催化裂化八集總動力學(xué)模型在自制的提升管模擬反應(yīng)器上進(jìn)行了FCC催化劑的選型。該工作對裂化催化劑的構(gòu)效關(guān)系和典型的催化裂化動力學(xué)模型進(jìn)行了分析。歐陽福生等[70]基于重油催化裂化十一集總動力學(xué)模型，對3種典型的FCC催化劑分別建模并通過比較各催化劑的反應(yīng)速率常數(shù)確定了各自的最佳操作條件。Vieira等[71]借助蠟油催化裂化八集總動力學(xué)模型，研究了3種工業(yè)平衡FCC催化劑的動力學(xué)行為。研究發(fā)現(xiàn)，原料中重組分的裂化反應(yīng)的動力學(xué)常數(shù)與催化劑的表面積和微反活性呈現(xiàn)較好的一致性。Quintana-Solorzano等[72]結(jié)合蠟油催化裂化五集總動力學(xué)模型，對4種FCC催化劑分別建立了集總動力學(xué)模型。結(jié)合各集總間反應(yīng)的速率常數(shù)和集總濃度，確定了催化裂化反應(yīng)體系中重要的反應(yīng)路徑，并分析了催化劑的組成結(jié)構(gòu)與其裂化性能間的關(guān)系。江洪波[73]提出借助催化裂化集總反應(yīng)動力學(xué)模型構(gòu)建FCC催化劑結(jié)構(gòu)-反應(yīng)性能模型與催化劑性質(zhì)數(shù)據(jù)庫的概念，以期找出催化劑結(jié)構(gòu)與反應(yīng)性能間的定量數(shù)學(xué)關(guān)系和催化劑性能變化的內(nèi)在規(guī)律，提高新催化劑的開發(fā)效率。

以上這些研究工作和概念的提出都不同程度地推動著集總反應(yīng)-催化劑性質(zhì)耦合的催化裂化模型的發(fā)展，這方面的研究還不夠成熟，可以預(yù)見的是未來將有更多相關(guān)工作的報道。

5 催化劑失活模型的發(fā)展

在催化裂化反應(yīng)過程中，催化劑失活現(xiàn)象是嚴(yán)重影響反應(yīng)的因素之一。因此，在催化裂化集總動力學(xué)模型的研究中，也需要考慮催化劑失活的因素。廣義的催化劑失活主要包含(水)熱失活、積炭失活和中毒失活3個方面。在集總動力學(xué)模型的研究中，通常只考察積炭失活、芳烴吸附失活和堿性氮毒害失活3個因素，而(水)熱失活[74]和重金屬等毒害失活[75]則研究較少。其中，積炭失活是對反應(yīng)速率影響最大的催化劑失活方式。對催化劑失活的描述主要有催化劑停留時間理論(time-on-stream)和催化劑積炭理論(coke-on-catalyst)2種理論。相應(yīng)的集總動力學(xué)模型對于這2種催化劑失活理論典型的數(shù)學(xué)描述分別見式(1)、(2)。

φ(tc)=e-α·tc

(1)

φ(wc)=(1+βwc)-M

(2)

式(1)、(2)中，φ為催化劑積炭失活函數(shù)；tc為催化劑停留時間，s；wc為催化劑上積炭質(zhì)量分?jǐn)?shù)，%；α、β和M為參數(shù)。

為了更準(zhǔn)確地采用數(shù)學(xué)模型描述催化劑的失活行為，學(xué)者們進(jìn)行了詳細(xì)深入的研究。眾多研究表明，催化劑積炭失活理論能更準(zhǔn)確地描述催化劑的失活行為。此外，由于反應(yīng)體系中各原料間、原料與產(chǎn)物間以及各產(chǎn)物間對于催化劑失活的貢獻(xiàn)大小有較大差別，因此需要針對不同原料和不同產(chǎn)物，分別制定合適的失活表達(dá)式以分別描述，即催化劑選擇性失活。研究表明[76-78]，催化劑的選擇性失活描述方法應(yīng)用效果良好，能更好地反映不同集總對催化劑失活的貢獻(xiàn)。芳烴吸附和堿性氮毒害對催化劑失活的影響，在集總動力學(xué)模型中典型的描述分別見式(3)、(4)。

(3)

(4)

式(3)、(4)中，f(A)和f(N)分別為催化劑芳烴吸附效應(yīng)失活函數(shù)和堿氮毒害效應(yīng)失活函數(shù)；kA和kN分別為芳烴吸附失活常數(shù)和堿氮毒害失活常數(shù)；wA、wR、wASP和wN分別為進(jìn)料中的芳烴、膠質(zhì)、瀝青質(zhì)和堿氮質(zhì)量分?jǐn)?shù)，%；θ為催化劑的相對停留時間，s；φC/O為劑/油質(zhì)量比。

芳烴吸附的影響一般包括芳烴、膠質(zhì)和瀝青質(zhì)，堿氮毒害的影響一般還需要考慮催化劑的相對停留時間和劑/油質(zhì)量比。有學(xué)者認(rèn)為，如果劑/油質(zhì)量比較高或者原料堿氮質(zhì)量分?jǐn)?shù)低于400 μg/g時，可以忽略堿氮毒害的影響。

6 催化裂化集總動力學(xué)模型建模數(shù)據(jù)來源的發(fā)展

集總動力學(xué)模型建模的傳統(tǒng)方法是通過實驗室級別的實驗裝置獲取建模數(shù)據(jù)，如ACE、固定流化床、小型和中型提升管等反應(yīng)器。這種建模數(shù)據(jù)的獲取方法，能在較短時間內(nèi)獲取大量數(shù)據(jù)，并且能得到操作條件寬范圍變化下的數(shù)據(jù)。這類數(shù)據(jù)建立的集總動力學(xué)模型對實驗室級別的裝置有較好的適應(yīng)性，但是，當(dāng)使用這類模型描述工業(yè)提升管中的反應(yīng)情況時，往往出現(xiàn)較大的偏差。研究者們?yōu)榱藨?yīng)對這一問題，通常采用設(shè)置裝置因數(shù)的辦法進(jìn)行修正。隨著集總劃分得越來越細(xì)致，模型參數(shù)越來越多，實驗室數(shù)據(jù)用于工業(yè)裝置的矛盾更加凸顯出來。于是，研究者們開始直接使用工業(yè)實測數(shù)據(jù)建模，如王建平等[24]、吳飛躍等[25]和段良偉等[27]的研究。工業(yè)運行數(shù)據(jù)能更好地描述工業(yè)提升管反應(yīng)器中反應(yīng)情況，研究表明，其應(yīng)用效果良好。目前，已有不少研究者在集總動力學(xué)建模過程中全部或者部分采用工業(yè)數(shù)據(jù)。

在集總動力學(xué)模型建立的過程中，使用實驗室數(shù)據(jù)還是工業(yè)運行數(shù)據(jù)作為建模數(shù)據(jù)，要視情況來決定。如果集總模型只是作為實驗室研究用，實驗室數(shù)據(jù)更合適；如果集總模型要應(yīng)用于工業(yè)提升管上，則工業(yè)運行數(shù)據(jù)更合適。此外，還要綜合考慮獲取數(shù)據(jù)的難易度。

7 催化裂化集總動力學(xué)模型參數(shù)估算方法的發(fā)展

催化裂化集總動力學(xué)模型參數(shù)估算是基于建模數(shù)據(jù)和反應(yīng)速率方程組表達(dá)式，按照相關(guān)求解和優(yōu)化算法編制程序，求取最合適的模型動力學(xué)參數(shù)的過程。隨著集總動力學(xué)模型的發(fā)展，各種參數(shù)估算方法不斷產(chǎn)生，如Gauss-Newton法、Marquart法及其改進(jìn)型算法、單純形法、Powell法、共軛梯度法、變尺度法(B-F-G-S)、粒子群優(yōu)化算法PSO、遺傳算法等。這些方法各具特點，通常從收斂性、穩(wěn)定性、可靠性、有效性、適用性、易用性、可移植性等方面進(jìn)行比較后按需選用。此外，由于模型的集總劃分越來越細(xì)致，導(dǎo)致需要估算的動力學(xué)參數(shù)數(shù)量龐大。為了確保估算得到的參數(shù)的可靠性，往往一次估算的參數(shù)數(shù)量不能太多，因而需要分步或者分層估算。這樣得到的動力學(xué)參數(shù)既不失準(zhǔn)確性，又能高效快速地獲取。

8 催化裂化集總動力學(xué)模型關(guān)聯(lián)模型的發(fā)展

由于催化裂化集總動力學(xué)模型是基于多個理想假設(shè)前提而建成的模型，這就決定了其應(yīng)用會存在一些缺陷。為了彌補應(yīng)用方面的不足，拓展應(yīng)用范圍，開發(fā)配合集總動力學(xué)模型使用的關(guān)聯(lián)模型具有較大的現(xiàn)實意義。

江洪波等[79]為了滿足石油化工行業(yè)要求多產(chǎn)FCC輕端產(chǎn)品的需求，開發(fā)了配合重油催化裂化十一集總動力學(xué)模型使用的輕端產(chǎn)品關(guān)聯(lián)模型。為了滿足催化裂化產(chǎn)品結(jié)構(gòu)隨市場需求的變化，應(yīng)對汽油與柴油的餾程范圍經(jīng)常波動的情況，Liu等[38]建立了汽柴油切割點變化與柴油收率的關(guān)聯(lián)模型——餾程修正模型。沙穎遜等[80]為了解決渣油原料的結(jié)構(gòu)族組成數(shù)據(jù)表征需核磁共振波譜儀的問題，開發(fā)了一組可由重質(zhì)原料的常規(guī)物性數(shù)據(jù)求取其結(jié)構(gòu)族組成數(shù)據(jù)的關(guān)聯(lián)式，稱為原料結(jié)構(gòu)-性質(zhì)模型。

關(guān)聯(lián)模型是集總動力學(xué)模型研究中的重要部分之一，其靈活性、簡易性和實用性有效地彌補了集總動力學(xué)模型的應(yīng)用缺陷，后續(xù)研究應(yīng)繼續(xù)對其功能加以豐富。

9 總結(jié)與展望

在諸多學(xué)者的不斷努力下，針對催化裂化(裂解)在不同時期的實際生產(chǎn)需要，催化裂化(裂解)集總動力學(xué)模型一直進(jìn)行著自我改進(jìn)與完善，取得了良好的應(yīng)用效果，并為煉油廠創(chuàng)造了巨大的經(jīng)濟效益。經(jīng)過近半個世紀(jì)的發(fā)展，催化裂化(裂解)集總動力學(xué)模型的集總劃分更加細(xì)致合理，數(shù)據(jù)來源更加多元化，參數(shù)估算方法更加高效可靠，輔助模型更加豐富實用，與流動模型耦合后功能更強大，這些成就保證并不斷推動著集總模型在工藝操作優(yōu)化、新工藝研發(fā)、清潔燃料和低碳烯烴生產(chǎn)、催化劑優(yōu)選和在線控制等方面的應(yīng)用。

目前，由于環(huán)境法規(guī)的日益嚴(yán)格和石油化工行業(yè)精細(xì)化、多方向的快速發(fā)展，市場對燃料和低碳烯烴從量和質(zhì)2個方面都提出了要求。因此，集總動力學(xué)模型在集總劃分方面更細(xì)致并更具工藝特色是勢在必行。事實上，各種被稱為后集總動力學(xué)模型時代動力學(xué)模型的分子尺度動力學(xué)模型已經(jīng)受到了廣泛關(guān)注，并已有許多報道。但是，這些分子尺度動力學(xué)模型尚處研究階段，還不成熟。這種情況下，集總動力學(xué)模型更應(yīng)基于其成熟可靠的優(yōu)勢，不斷完善并拓展應(yīng)用，以期繼續(xù)占據(jù)動力學(xué)模型研究的核心位置。另外，集總模型與流動模型耦合后建立的模型，能更好地描述反應(yīng)器中流動-反應(yīng)耦合效應(yīng)對產(chǎn)品分布的綜合影響，其應(yīng)用效果已得到廣泛認(rèn)可。最后，集總反應(yīng)模型與催化劑性質(zhì)的耦合將使催化裂化模型同時適用于多種催化劑，從而極大地拓寬其應(yīng)用范圍。因此，集總反應(yīng)模型與流動、傳熱、傳質(zhì)等模型和催化劑性質(zhì)的耦合考察，也應(yīng)該繼續(xù)重點研究。

[1] 許友好. 我國催化裂化工藝技術(shù)進(jìn)展[J]. 中國科學(xué): 化學(xué), 2014, 44(1): 13-24. (XU Youhao. Advance in China fluid catalytic cracking (FCC) process[J]. Scientia Sinica Chimica, 2014, 44(1): 13-24.)

[2] ARIS R. Prolegomena to the rational analysis of system of chemical reactions[J]. Archive for Rational Mechanics and Analysis, 1965, 19(2): 81-99.

[3] WEI J, KUO J C W. Lumping analysis in monomolecular reaction systems. Analysis of the exactly lumpable system[J]. Industrial and Engineering Chemistry Fundamentals, 1969, 8(1): 114-123.

[4] WEEKMAN V W, NACE D M. Kinetics of catalytic cracking selectivity in fixed, moving, and fluid bed reactors[J]. AIChE Journal, 1970, 16(3): 397-404.

[5] LEE L S, CHEN Y W, HUANG T N, et al. Four-lump kinetic model for fluid catalytic cracking process[J]. Canadian Journal of Chemical Engineering, 1989, 67(4): 615-619.

[6] ANCHEYTA-JUAREZ J, LOPEZ-ISUNZA F, AGUILAR-RODRIGUEZ E. 5-Lump kinetic model for gas oil catalytic cracking[J]. Applied Catalysis A: General, 1999, 177(2): 227-235.

[7] KRAEMER D W, DE LASA H. Modeling catalytic cracking in a novel riser simulator[C]// Washington, D. C.: AIChE Annual Meeting, 1988.

[8] 任杰, 翁惠新, 劉馥英. 催化裂化反應(yīng)八集總動力學(xué)模型的初步研究[J]. 石油學(xué)報(石油加工), 1994, 10(1): 1-7. (REN Jie, WENG Huixin, LIU Fuying. Investigation of the lumped kinetic model for catalytic cracking reaction[J]. Acta Petrolei Sinica (Petroleum Processing Section), 1994, 10(1): 1-7.)

[9] HAGELBERG P, EILOS I, HILTUNEN J, et al. Kinetics of catalytic cracking with short contact times[J]. Applied Catalysis A: General, 2002, 223(1-2): 73-84.

[10] JACOB S M, GROSS B, WEEKMAN V W. A lumping and reaction scheme for catalytic cracking[J]. AIChE Journal, 1976, 22(4): 701-713.

[11] 朱炳辰, 翁惠新, 朱子郴. 催化反應(yīng)工程[M]. 北京: 中國石化出版社, 1999: 52-53.

[12] 劉之定, 胡九成. 催化裂化反應(yīng)動力學(xué)模型與優(yōu)化研究[J]. 化學(xué)反應(yīng)工程與工藝, 1994, 10(2): 118-123. (LIU Zhiding, HU Jiucheng. Research on the kinetic model of catalytic cracking and optimization[J]. Chemical Reaction Engineering and Technology, 1994, 10(2): 118-123.)

[13] ARAUJO-MONROY C, LOPEZ-ISUNZA F. Modeling and simulation of an industrial fluid catalytic cracking riser reactor using a lump-kinetic model for a distinct feedstock[J]. Industrial and Engineering Chemistry Research, 2006, 45(1): 120-128.

[14] PITAULT I, NEVICATO D, FORISSIER M. Kinetic model based on a molecular description for catalytic cracking of vacuum gas oil[J]. Chemical Engineering Science, 1994, 49: 4249-4262.

[15] TAKATSUKA T, SATO S, MORIMOTO Y, et al. A reaction model for fluidized-bed catalytic cracking of residual oil[J]. International Chemical Engineering, 1987, 27(1): 107-116.

[16] 牟盛靜, 蘇宏業(yè), 栗偉, 等. 工業(yè)渣油催化裂化過程六集總組合反應(yīng)器模型[J]. 高?；瘜W(xué)工程學(xué)報, 2005, 19(5): 630-635. (MU Shengjing, SU Hongye, LI Wei, et al. Reactor model for industrial residual fluid catalytic cracking based on six-lump kinetic model[J]. Journal of Chemical Engineering of Chinese Universities, 2005, 19(5): 630-635.)

[17] XU O G, SU H Y, MU S J, et al. 7-Lump kinetic model for residual oil catalytic cracking[J]. Journal of Zhejiang University: Science A, 2006, 7(11): 1932-1941.

[18] GAO H H, WANG G, LI R, et al. Study on the catalytic cracking of heavy oil by proper cut for higher conversion and desirable products[J]. Energy Fuels, 2012, 26(3): 1880-1891.

[19] 任杰. 重油催化裂化集總動力學(xué)模型研究[D]. 上海: 華東理工大學(xué), 1992.

[20] 鄧先梁, 沙穎遜, 王龍延, 等. 渣油催化裂化反應(yīng)動力學(xué)模型的研究[J]. 石油煉制與化工, 1994, 25(9): 35-39. (DENG Xianliang, SHA Yingxun, WANG Longyan, et al. Studies on a kinetic model of resid catalytic cracking[J]. Petroleum Processing and Petrochemicals, 1994, 25(9): 35-39.)

[21] 沙穎遜, 孟凡東, 閻遂寧, 等. 催化裂化反應(yīng)再生系統(tǒng)模擬優(yōu)化軟件的研究和應(yīng)用[J]. 石油煉制與化工, 1997, 28(9): 52-57. (SHA Yingxun, MENG Fandong, YAN Suining, et al. Development of simulation software FCCM and its application [J]. Petroleum Processing and Petrochemicals, 1997, 28(9): 52-57.)

[22] 孫鐵棟, 鐘孝湘, 陳昀, 等. 渣油催化裂化集總動力學(xué)模型的研究與應(yīng)用[J]. 石油煉制與化工, 2001, 32(4): 41-44. (SUN Tiedong, ZHONG Xiaoxiang, CHEN Yun, et al. Study and application of a lumping FCC kinetics model for vacuum residue[J]. Petroleum Processing and Petrochemicals, 2001, 32(4): 41-44.)

[23] PASHIKANTI K, LIU Y A. Predictive modeling of large-scale integrated refinery reaction and fractionation system from plant data. Part 2: Fluid catalytic cracking (FCC) process[J]. Energy Fuels, 2011, 25(11): 5298-5319.

[24] 王建平, 許先焜, 翁惠新, 等. 加氫渣油催化裂化14集總動力學(xué)模型的建立[J]. 化工學(xué)報, 2007, 58(1): 86-94. (WANG Jianping, XU Xiankun, WENG Huixin, et al. Establishment of 14 lumps model for fluid catalytic cracking of hydrotreated residuum[J]. Journal of Chemical Industry and Engineering (China), 2007, 58(1): 86-94.)

[25] 吳飛躍, 翁惠新. FDFCC集總反應(yīng)動力學(xué)組合模型的建立[J]. 華東理工大學(xué)學(xué)報(自然科學(xué)版), 2009, 35(3): 350-356. (WU Feiyue, WENG Huixin. Establishment of combined model for lumped reaction kinetics of FDFCC[J]. Journal of East China University of Science and Technology (Natural Science Edition), 2009, 35(3): 350-356.)

[26] 江洪波, 鐘貴江, 寧匯, 等. 重油催化裂化MIP工藝集總動力學(xué)模型[J]. 石油學(xué)報(石油加工), 2010, 26(6): 901-909. (JIANG Hongbo, ZHONG Guijiang, NING Hui, et al. Lumping kinetic model of heavy oil catalytic cracking for MIP technology[J]. Acta Petrolei Sinica (Petroleum Processing Section), 2010, 26(6): 901-909.)

[27] 段良偉, 孫鵬, 翁惠新. 新型催化裂化MIP集總反應(yīng)動力學(xué)模型[J]. 華東理工大學(xué)學(xué)報(自然科學(xué)版), 2012, 38(1): 17-22. (DUAN Liangwei, SUN Peng, WENG Huixin. Lumping kinetic model for new catalytic cracking of MIP process[J]. Journal of East China University of Science and Technology (Natural Science Edition), 2012, 38(1): 17-22.)

[28] 劉熠斌, 楊朝合, 山紅紅, 等. 提高汽、柴油收率的兩段提升管催化裂化動力學(xué)模型研究及應(yīng)用[J]. 石油學(xué)報(石油加工), 2007, 23(5): 7-14. (LIU Yibin, YANG Chaohe, SHAN Honghong， et al. Study and application of kinetic model of two-stage riser FCC-maximizing diesel and gasoline technology[J]. Acta Petrolei Sinica (Petroleum Processing Section), 2007, 23(5): 7-14.)

[29] VILLAFUERTE-MACIAS E F, AGUILAR R, MAYA-YESCAS R. Towards modelling production of clean fuels: Sour gas formation in catalytic cracking[J]. Journal of Chemical Technology and Biotechnology, 2004, 79(10): 1113-1118.

[30] 許友好. 催化裂解反應(yīng)動力學(xué)模型的建立及其應(yīng)用[J]. 石油煉制與化工, 2001, 32(11): 44-47. (XU Youhao. Establishment of kinetic model for DCC and its application[J]. Petroleum Processing and Petrochemicals, 2001, 32(11): 44-47.)

[31] 張琪皓. 下行床反應(yīng)器烴類催化裂解/裂化熱態(tài)試驗研究[D]. 北京: 清華大學(xué), 2003.

[32] 陳新國. 大慶渣油催化裂解反應(yīng)動力學(xué)及反應(yīng)系統(tǒng)模擬研究[D]. 北京: 中國石油大學(xué), 2002.

[33] 孟祥海, 高金森, 徐春明, 等. 大慶常壓渣油催化裂解反應(yīng)動力學(xué)模型研究[J]. 化學(xué)反應(yīng)工程與工藝, 2003, 19(4): 365-371. (MENG Xianghai, GAO Jinsen, XU Chunming, et al. Studies on the kinetic model of Daqing atmospheric residue catalytic pyrolysis[J]. Chemical Reaction Engineering and Technology, 2003, 19(4): 365-371.)

[34] MENG X H, XU C M, GAO J S, et al. Seven-lump kinetic model for catalytic pyrolysis of heavy oil[J]. Catalysis Communications, 2007, 8(8): 1197-1201.

[35] MENG X H, XU C M, GAO J S, et al. Catalytic pyrolysis of heavy oils: 8-Lump kinetic model[J]. Applied Catalysis A: General, 2006, 301(1): 32-38.

[36] MENG X H, XU C M, LI L, et al. Kinetics of catalytic pyrolysis of heavy gas oil derived from Canadian synthetic crude oil[J]. Energy Fuels, 2011, 25(8): 3400-3407.

[37] 王國良, 王龍延, 崔中強, 等. 重油接觸裂解制乙烯集總反應(yīng)動力學(xué)模型研究[J]. 石油化工, 2004, 33(2): 93-99. (WANG Guoliang, WANG Longyan, CUI Zhongqiang, et al. Study on lumped kinetic reaction network of heavy oil contact cracking process [J]. Petrochemical Technology, 2004, 33(2): 93-99.)

[38] LIU Y B, CHEN X B, ZHAO H, et al. Establishment of kinetic model for catalytic pyrolysis of Daqing atmospheric residue[J]. Chinese Journal of Chemical Engineering, 2009, 17(1): 78-82.

[39] 杜峰, 劉熠斌, 楊朝合, 等. 大慶常壓渣油催化裂解反應(yīng)集總動力學(xué)模型[J]. 中國石油大學(xué)學(xué)報(自然科學(xué)版), 2008, 32(5): 137-141. (DU Feng, LIU Yibin, YANG Chaohe, et al. Lumping kinetic model for catalytic pyrolysis of Daqing atmospheric residue[J]. Journal of China University of Petroleum, 2008, 32(5): 137-141.)

[40] 郭菊花. 重油兩段催化裂解多產(chǎn)丙烯集總動力學(xué)模型的初步研究[D]. 青島: 中國石油大學(xué), 2008.

[41] 劉熠斌, 陳小博, 楊朝合. FCC汽油催化裂解生產(chǎn)低碳烯烴的動力學(xué)模型研究[J]. 煉油技術(shù)與工程, 2007, 37(3): 5-9. (LIU Yibin, CHEN Xiaobo, YANG Chaohe. Study on kinetic model of light olefins production from FCC naphtha[J]. Petroleum Refinery Engineering, 2007, 37(3): 5-9.)

[42] 韓忠祥, 劉熠斌, 楊朝合. 催化裂化汽油二次裂解反應(yīng)集總動力學(xué)模型[J]. 化學(xué)反應(yīng)工程與工藝, 2011, 27(1): 73-78. (HAN Zhongxiang, LIU Yibin, YANG Chaohe. Study on lump kinetic model for secondary cracking of FCC gasoline[J]. Chemical Reaction Engineering and Technology, 2011, 27(1): 73-78.)

[43] 吳飛躍, 翁惠新, 羅世賢. FDFCC工藝中汽油提升管催化裂化反應(yīng)動力學(xué)模型研究[J]. 石油煉制與化工, 2007, 38(10): 59-63. (WU Feiyue, WENG Huixin, LUO Shixian. Study on the kinetic model for naphtha cracking in riser of FDFCC process [J]. Petroleum Processing and Petrochemicals, 2007, 38(10): 59-63.)

[44] WANG L Y, YANG B L, WANG Z W. Lumps and kinetics for the secondary reactions in catalytically cracked gasoline[J]. Chemical Engineering Journal, 2005, 109: 1-9.

[45] 劉福安, 侯栓弟, 武雪鋒, 等. 催化裂化汽油改質(zhì)反應(yīng)動力學(xué)模型研究[J]. 煉油技術(shù)與工程, 2005, 35(6): 25-30. (LIU Fuan, HOU Shuandi, WU Xuefeng, et al. Study on reaction kinetic model of FCC naphtha upgrading[J]. Petroleum Refinery Engineering, 2005, 35(6): 25-30.)

[46] 吳青, 何鳴元. 催化裂化汽油改質(zhì)過程的八集總動力學(xué)模型研究[J]. 煉油技術(shù)與工程, 2004, 34(11): 22-26. (WU Qing, HE Mingyuan. Study on eight lumping dynamic models for FCC gasoline upgrading[J]. Petroleum Refinery Engineering, 2004, 34(11): 22-26.)

[47] 吳青, 何鳴元. 催化裂化汽油改質(zhì)過程的20集總動力學(xué)模型研究 [J]. 石油化工, 2005, 34： 675-677. (WU Qing, HE Mingyuan. Study on twenty lumping dynamic models for FCC gasoline upgrading [J]. Petrochemical Technology, 2005, 34： 675-677.)

[48] 楊光福, 王剛, 田廣武, 等. 催化裂化汽油改質(zhì)反應(yīng)動力學(xué)模型研究[J]. 燃料化學(xué)學(xué)報, 2007, 35(3): 297-301.(YANG Guangfu, WANG Gang, TIAN Guangwu, et al. Kinetic model of FCC naphtha olefin reformulation[J]. Journal of Fuel Chemistry and Technology, 2007, 35(3): 297-301.)

[49] 戴鑑, 楊光福, 王剛, 等. 催化裂化汽油改質(zhì)降烯烴并多產(chǎn)丙烯的反應(yīng)動力學(xué)模型研究[J]. 燃料化學(xué)學(xué)報, 2008, 36(4): 431-436. (DAI Jian, YANG Guangfu, WANG Gang, et al. Kinetic model of FCC naphtha reformulation for reducing olefin and enhancing propene production[J]. Journal of Fuel Chemistry and Technology, 2008, 36(4): 431-436.)

[50] 鄒圣武, 侯栓弟, 龍軍, 等. 基于汽油烯烴轉(zhuǎn)化的催化裂化動力學(xué)模型[J]. 化工學(xué)報, 2004, 55(11): 1793-1798. (ZOU Shengwu, HOU Shuandi, LONG Jun, et al. Kinetic model of catalytic cracking based on olefin reaction mechanism[J]. Journal of Chemical Industry and Engineering (China), 2004, 55(11): 1793-1798.)

[51] 高金森, 徐春明, 楊光華, 等. 催化裂化提升管反應(yīng)器模型的建立[J]. 高?；瘜W(xué)工程學(xué)報, 1998, 12(3): 259-264. (GAO Jinsen, XU Chunming, YANG Guanghua, et al. Turbulent gas-solid two-phase flow-reaction model[J]. Journal of Chemical Engineering of Chinese Universities, 1998, 12(3): 259-264.)

[52] 高金森, 徐春明, 林世雄, 等. 催化裂化提升管反應(yīng)器氣液固3相流動反應(yīng)的數(shù)值模擬Ⅰ氣液固3相流動反應(yīng)模型的建立[J]. 石油學(xué)報(石油加工), 1999, 15(3): 16-22. (GAO Jinsen, XU Chunming, LIN Shixiong, et al. Numerical simulation on the gas-liquid-solid three-phase flow-reaction in FCC riser reactors Ⅰ Development of the gas-liquid-solid three-phase flow-reaction model [J]. Acta Petrolei Sinica (Petroleum Processing Section), 1999, 15(3): 16-22.)

[53] 陳新國, 徐春明, 郭印誠. 顆粒動力學(xué)模型的建立[J]. 計算機與應(yīng)用化學(xué), 2000, 17(6): 505-508. (CHEN Xinguo, XU Chunming, GUO Yincheng. Development of kinetic model of granular flow[J]. Computers and Applied Chemistry, 2000, 17(6): 505-508.)

[54] 徐春明, 王洪斌, 郭印誠, 等. 渣油催化裂化提升管反應(yīng)器組成分布的數(shù)值模擬[J]. 化工學(xué)報, 1999, 50(1): 57-64. (XU Chunming, WANG Hongbin, GUO Yincheng, et al. Numerical simulation of product yields in RFCC riser reactor[J]. Journal of Chemical Industry and Engineering (China), 1999, 50(1): 57-64.)

[55] WANG S Y, HE Y R, LU H L, et al. Simulation of performance of cracking reactions of particle clusters in FCC risers [J]. Industrial and Engineering Chemistry Research, 2008, 47(14): 4632-4640.

[56] 侯栓弟, 龍軍, 郭湘波, 等. 重油催化裂化反應(yīng)歷程數(shù)值模擬Ⅰ模型建立及流體力學(xué)性能模擬[J]. 石油學(xué)報(石油加工), 2005, 21(5): 19-27. (HOU Shuandi, LONG Jun, GUO Xiangbo, et al. Numerical simulation of reaction pathway in commercial riser of heavy oil catalytic cracking unit ⅠModel estimation and flow behavior simulation[J]. Acta Petrolei Sinica (Petroleum Processing Section), 2005, 21(5): 19-27.)

[57] 郭湘波, 侯栓弟, 龍軍, 等. 重油催化裂化反應(yīng)歷程數(shù)值模擬Ⅱ重油催化反應(yīng)歷程數(shù)值模擬[J]. 石油學(xué)報(石油加工), 2006, 22(1): 9-14. (GUO Xiangbo, HOU Shuandi, LONG Jun, et al. Numerical simulation of reaction pathway in commercial riser of heavy oil catalytic cracking unit Ⅱ Numerical simulation of heavy oil catalytic cracking process[J]. Acta Petrolei Sinica (Petroleum Processing Section), 2006, 22(1): 9-14.)

[58] 侯栓弟. 汽油催化改質(zhì)反應(yīng)過程數(shù)值模擬[J]. 化工學(xué)報, 2007, 58(3): 623-629. (HOU Shuandi. Numerical simulation of FCC gasoline catalytic conversion process[J]. Journal of Chemical Industry and Engineering (China), 2007, 58(3): 623-629.)

[59] CHANG J, ZHANG K, MENG F D, et al. Computational investigation of hydrodynamics and cracking reaction in a heavy oil riser reactor [J]. Particuology, 2012, 10(2): 184-195.

[60] CHANG J, MENG F D, WANG L Y, et al. CFD investigation of hydrodynamics, heat transfer and cracking reaction in a heavy oil riser with bottom airlift loop mixer[J]. Chemical Engineering Science, 2012, 78: 128-143.

[61] CHANG J, CAI W Z, ZHANG K, et al. Computational investigation of the hydrodynamics, heat transfer and kinetic reaction in a FCC gasoline riser [J]. Chemical Engineering Science, 2014, 111: 170-179.

[62] LI J, FAN Y P, LU C X, et al. Numerical simulation of influence of feed injection on hydrodynamic behavior and catalytic cracking reactions in a FCC riser under reactive conditions[J]. Industrial and Engineering Chemistry Research, 2013, 52(32): 11084-11098.

[63] CHEN G Q, LUO Z H, LAN X Y, et al. Evaluating the role of intraparticle mass and heat transfers in a commercial FCC riser: A meso-scale study[J]. Chemical Engineering Journal, 2013, 228: 352-365.

[64] CHEN G Q, LUO Z H. New insights into intraparticle transfer, particle kinetics, and gas-solid two-phase flow in polydisperse fluid catalytic cracking riser reactors under reaction conditions using multi-scale modeling[J]. Chemical Engineering Science, 2014, 109: 38-52.

[65] DEROUIN C, NEVICATO D, FORISSIER M, et al. Hydrodynamics of riser units and their impact on FCC operation[J]. Industrial and Engineering Chemistry Research, 1997, 36(11): 4504-4515.

[66] DAS A K, BAUDREZ E, MARIN G B, et al. Three-dimensional simulation of a fluid catalytic cracking riser reactor[J]. Industrial and Engineering Chemistry Research, 2003, 42(12): 2602-2617.

[67] LOPES G C, ROSA L M, MORI M, et al. Three-dimensional modeling of fluid catalytic cracking industrial riser flow and reactions[J]. Computers Chemical Engineering, 2011, 35(11): 2159-2168.

[68] ZHU C, JUN Y, PATEL R, et al. Interactions of flow and reaction in fluid catalytic cracking risers[J]. AIChE Journal, 2011, 57(11): 3122-3131.

[69] DE LASA H, DOGU G, RAVELLA A. Chemical Reactor Technology for Environmentally Safe Reactors and Products[M]. Amsterdam: Springer Netherlands, 1992: 71-131.

[70] 歐陽福生, 翁惠新. 重油催化裂化集總模型在指導(dǎo)裂化催化劑選用方面的研究[J]. 華東理工大學(xué)學(xué)報, 1996, 22(6): 666-674. (OUYANG Fusheng, WENG Huixin. Investigation of the lump model for heavy oil catalytic cracking on guiding the selection and usage of cracking catalysts[J]. Journal of East China University of Science and Technology, 1996, 22(6): 666-674.)

[71] VIEIRA R C, PINTO J C, BISCAIA E C, et al. Simulation of catalytic cracking in a fixed-fluidized-bed unit[J]. Industrial and Engineering Chemistry Research, 2004, 43(19): 6027-6034.

[72] QUINTANA-SOLORZANO R,RODRIGUEZ-HERNANDEZ A, GARCIA-DE-LEON R. Study of the performance of catalysts for catalytic cracking by applying a lump-based kinetic model[J]. Industrial and Engineering Chemistry Research, 2009, 48(3): 1163-1171.

[73] 江洪波. 重油催化裂化集總動力學(xué)模型的拓展與應(yīng)用研究[D]. 上海: 華東理工大學(xué), 1999.

[74] 任杰. 催化裂化催化劑水熱失活動力學(xué)模型[J]. 石油學(xué)報(石油加工), 2002, 18(5): 40-46. (REN Jie. Kinetic model of hydrothermal deactivation for catalytic cracking catalyst[J]. Acta Petrolei Sinica (Petroleum Processing Section), 2002, 18(5): 40-46.)

[75] LAROCCA M, FARAG H, NG S, et al. Cracking catalyst deactivation by nickel and vanadium contaminants[J]. Industrial and Engineering Chemistry Research, 1990, 29(11): 2181-2191.

[76] CORELLA J. On the modeling of the kinetics of the selective deactivation of catalysts. Application to the fluidized catalytic cracking process[J]. Industrial and Engineering Chemistry Research, 2004, 43(15): 4080-4086.

[77] BOLLAS G M, LAPPAS A A, IATRIDIS D K, et al. Five-lump kinetic model with selective catalyst deactivation for the prediction of the product selectivity in the fluid catalytic cracking process[J]. Catalysis Today, 2007, 127(1-4): 31-43.

[78] ZHOU X W, CHEN T, YANG B L, et al. Kinetic model considering reactant oriented selective deactivation for secondary reactions of fluid catalytic cracking gasoline[J]. Energy Fuels, 2011, 25(6): 2427-2437.

[79] 江洪波, 翁惠新. 催化裂化輕端產(chǎn)品關(guān)聯(lián)模型[J]. 化學(xué)反應(yīng)工程與工藝, 2001, 17(1): 39-43. (JIANG Hongbo, WENG Huixin. Common methods and principles of lumping kinetic model research[J]. Chemical Reaction Engineering and Technology, 2001, 17(1): 39-43.)

[80] 沙穎遜, 孟凡東, 閻遂寧. FCCU反再系統(tǒng)模擬優(yōu)化軟件的研究和應(yīng)用(I)[J]. 催化裂化, 1997, 16(2): 17-21. (SHA Yingxun, MENG Fandong, YAN Suining. Development of simulation software FCCM and its application (I)[J]. Catalytic Cracking, 1997, 16(2): 17-21.)

Research Progresses of Lump Kinetic Model of FCC and Catalytic Pyrolysis

XIONG Kai，LU Chunxi

(StateKeyLaboratoryofHeavyOilProcessing,ChinaUniversityofPetroleum,Beijing102249,China)

The research progresses of lump kinetic modelling of fluid catalytic cracking (FCC) and catalytic pyrolysis were reviewed from the following aspects of FCC lump kinetic model, catalytic pyrolysis lump kinetic model, lump reaction-hydrodynamics coupled FCC model, catalyst property-included FCC model, catalyst deactivation model, source of experimental data, estimation method of kinetic parameters and correlation model of lump kinetic model. Based on the analysis, several suggestions for further development of FCC and catalytic pyrolysis lump kinetic modeling were proposed.

fluid catalytic cracking (FCC); catalytic pyrolysis; lump model; kinetic; hydrodynamics; catalyst property

2014-08-14 第一作者： 熊凱，男，博士研究生，從事催化裂化反應(yīng)工程和工藝研究

盧春喜，男，教授，博士，從事多相流傳遞與反應(yīng)工程、催化裂化、流態(tài)化方面的研究；Tel：010-89733803；E-mail：lcx725@sina.com

1001-8719(2015)02-0293-14

TE624

A

10.3969/j.issn.1001-8719.2015.02.010