• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

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

    2015-04-07 02:52:57盧春喜
    石油學(xué)報(石油加工) 2015年2期
    關(guān)鍵詞:集總蠟油重油

    熊 凱, 盧春喜

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

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

    熊 凱, 盧春喜

    (中國石油大學(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é)模型的研究進展。在分析總結(jié)的基礎(chǔ)上,針對催化裂化(裂解)集總反應(yīng)動力學(xué)模型研究的后續(xù)發(fā)展提出了幾點建議。

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

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

    催化裂化反應(yīng)中涉及的反應(yīng)種類繁多,且多為平行-順序反應(yīng),各反應(yīng)間偶聯(lián)性很強。20世紀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è)計、先進控制和催化劑的選型方面成效顯著,為煉油廠創(chuàng)造了巨大的經(jīng)濟效益。

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

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

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

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

    20世紀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ǔ)上,將氣體集總細分為干氣和液化氣2個集總,以便在催化裂化反應(yīng)過程中能夠分別監(jiān)測干氣和液化氣的生成。Kraemer等[7]提出了蠟油催化裂化八集總模型,將原料油按餾程分為輕燃料油(LFO)(220~345℃)和重燃料油(HFO)(>345℃)2部分,然后將LFO和HFO分別細分為烷烴(P)、環(huán)烷烴(N)和芳烴(A)3個集總,再加上汽油和氣體+焦炭,共8個集總。任杰等[8]按照我國煉油廠對油品餾程的劃分習(xí)慣,對原料油和汽油集總餾程范圍重新進行了劃分,在他們建立的蠟油八集總模型中,將原料油按餾程分為輕餾分油(LGO)(204~350℃)和重餾分油(HGO)(350~500℃)2部分。在五集總模型的基礎(chǔ)上,Hagelberg等[9]提出了一種新的蠟油催化裂化八集總模型。在此模型中,將汽油集總按照PONA組成細分為烷烴、烯烴、環(huán)烷烴和芳烴4個集總,而原料不細分。任杰等的蠟油八集總模型為原料集總細分模型,因而原料適應(yīng)性較強;Hagelberg等的蠟油八集總模型為產(chǎn)品細分模型,不僅能預(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各自再分別細分為烷烴(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ù)細分為一、二環(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)組分的描述較之于簡單集總模型細致了許多,Pitault等認為,十九集總模型已介于簡單集總動力學(xué)模型與分子級別動力學(xué)模型之間。

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

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

    20世紀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]進一步將裂化氣集總細分為液化氣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)烷烴和芳烴分別細分為3個集總,再加上汽油和氣體+焦炭2個集總,共11個集總。鄧先梁等[20]和沙穎遜等[21]進一步將氣體和焦炭分別集總,并將渣油層中的芳烴集總細分為芳烴和膠質(zhì)+瀝青質(zhì)2個集總,提出了重油催化裂化十三集總模型。孫鐵棟等[22]提出了渣油催化裂化十四集總模型。將回?zé)捰蛦为毤偅瑢⒃戏譃樵?、蠟油和重循環(huán)油3層,并按照SARA組成將各層分別細分為飽和分、芳香分和膠質(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é)模型能更細致地描述和調(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族組成進行了集總劃分。后續(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)族組成進行集總劃分,并對反應(yīng)路徑進行了合理簡化。該模型除了能預(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進行了單獨集總。該模型能預(yù)測硫元素在產(chǎn)物中的分布情況,并能反過來促進操作優(yōu)化,盡可能避免硫元素進入汽油和焦炭產(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ù)目進行了擴大,分別建立了七集總和八集總模型。將烯烴集總細分為乙烯和丙烯+丁烯2個集總,將烷烴+氫氣集總細分為丙烷+丁烷和甲烷+乙烷+氫氣2個集總,成為七集總模型;再將重油原料細分為芳香碳和非芳香碳2個集總,成為八集總模型。這些嘗試都取得了較好的效果,對產(chǎn)品分布進行了越來越細致的描述。由于大慶常壓渣油屬于石蠟基,為了研究不同質(zhì)量進料的催化裂解反應(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)族組成分別表示重油層和柴油層,裂化氣細分為氫氣、甲烷、乙烯、乙烷、丙烯、丙烷、丁烯和丁烷,共16個集總。模型體現(xiàn)了反應(yīng)溫度、壓力、劑/油質(zhì)量比、水/油質(zhì)量比、催化劑活性等因素的影響,各產(chǎn)物產(chǎn)率的計算值和實驗值絕對偏差一般在0.5%左右,模型精度和預(yù)測能力均較好。

    Liu等[38]將大慶常壓渣油催化裂解反應(yīng)體系重新進行了集總劃分,基于兩段提升管催化裂解多產(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ǔ)上,將汽油細分為汽油烯烴、汽油芳烴和汽油飽和烴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汽油裂解八集總模型。他們認為,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組成進行了集總劃分,分為柴油、汽油正構(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組成對汽油進行了集總劃分,在固定流化床上研究了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)烷烴的詳細碳數(shù)分布,構(gòu)建了更為詳細的反應(yīng)網(wǎng)絡(luò),開發(fā)了FCC汽油改質(zhì)二十集總動力學(xué)模型。在八集總的基礎(chǔ)上,異構(gòu)烷烴和烯烴分別按照碳數(shù)細分為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組成進行集總劃分,提出了基于汽油烯烴轉(zhuǎn)化的十二集總模型。模型的集總分為餾分油(烷烴、烯烴、環(huán)烷烴和芳烴)、汽油(烷烴、烯烴、環(huán)烷烴和芳烴)、液化氣烷烴、液化氣烯烴、干氣和焦炭,共12個集總。研究表明,該模型能預(yù)測不同條件下的產(chǎn)物產(chǎn)率分布,包括汽油的詳細組成,有利于降低汽油烯烴含量的研究。

    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)器進行了三維微分模擬,考察了原料液霧流動氣化過程及其影響,建立了提升管氣-液-固三相流動-反應(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套不同工況進行數(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]進行了模擬研究。結(jié)果表明,重油催化裂化反應(yīng)主要發(fā)生在噴嘴附近區(qū)域,在這一區(qū)域約45%的重油轉(zhuǎn)化為汽油和柴油。此外,針對FCC催化汽油改質(zhì)工藝,侯栓弟[58]報道了FCC汽油改質(zhì)反應(yīng)-流動耦合模型。針對不同結(jié)構(gòu)的反應(yīng)器對汽油改質(zhì)反應(yīng)和流動過程進行了數(shù)值模擬,結(jié)果表明,對提升管反應(yīng)器,F(xiàn)CC汽油經(jīng)低溫改質(zhì)后,烯烴含量降低了約一半;對于提升管-床層反應(yīng)器,汽油中的烯烴質(zhì)量分數(shù)能降至約5%。

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

    Li等[62]借助重油催化裂化十四集總動力學(xué)模型,考察了提升管反應(yīng)器的進料情況對反應(yīng)器油-劑混合段的流體流動行為和催化裂化反應(yīng)的影響。模擬結(jié)果表明,進料噴嘴出口速率決定了油-劑混合區(qū)的氣-固兩相流動情況,噴嘴的角度會對提升管內(nèi)油-劑混合區(qū)的流場和裂化反應(yīng)產(chǎn)生巨大影響。Chen等[63]借助蠟油催化裂化四集總動力學(xué)模型,重點考察了FCC提升管反應(yīng)器中催化劑的顆粒內(nèi)傳質(zhì)和傳熱行為,在綜合考慮質(zhì)量平衡、熱平衡、動量平衡、氣態(tài)方程、集總反應(yīng)動力學(xué)方程、多組分擴散和對流傳熱的基礎(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é)合時,三維流動模型比一維流動模型能更準確地預(yù)測產(chǎn)品分布。Lopes等[67]借助蠟油催化裂化四集總動力學(xué)模型,也對提升管中的流動和反應(yīng)行為進行了三維模擬,并取得了良好效果。Zhu等[68]將蠟油催化裂化四集總動力學(xué)模型與提升管中的流體流動行為進行耦合分析,考察了提升管中流體流動情況對催化裂化反應(yīng)的影響。由于催化劑的濃度從提升管底部到頂部逐漸降低,計算結(jié)果表明,提升管底部的裂化強度要遠高于傳統(tǒng)的提升管模型中的情況。與傳統(tǒng)的提升管反應(yīng)模型比較,新開發(fā)的流動-反應(yīng)耦合模型能更好地預(yù)測工業(yè)生產(chǎn)中的產(chǎn)品分布情況。

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

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

    集總反應(yīng)動力學(xué)模型在建模過程中?;谀撤N特定的催化劑,而不考慮不同催化劑對模型動力學(xué)參數(shù)的影響,即將催化劑因素視為集總動力學(xué)模型內(nèi)部隱含的一個恒量。顯然,不同種類的催化劑常因為功能各異而對催化裂化的產(chǎ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)器上進行了FCC催化劑的選型。該工作對裂化催化劑的構(gòu)效關(guān)系和典型的催化裂化動力學(xué)模型進行了分析。歐陽福生等[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ī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ì)量分數(shù),%;α、β和M為參數(shù)。

    為了更準確地采用數(shù)學(xué)模型描述催化劑的失活行為,學(xué)者們進行了詳細深入的研究。眾多研究表明,催化劑積炭失活理論能更準確地描述催化劑的失活行為。此外,由于反應(yīng)體系中各原料間、原料與產(chǎn)物間以及各產(chǎn)物間對于催化劑失活的貢獻大小有較大差別,因此需要針對不同原料和不同產(chǎn)物,分別制定合適的失活表達式以分別描述,即催化劑選擇性失活。研究表明[76-78],催化劑的選擇性失活描述方法應(yīng)用效果良好,能更好地反映不同集總對催化劑失活的貢獻。芳烴吸附和堿性氮毒害對催化劑失活的影響,在集總動力學(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分別為進料中的芳烴、膠質(zhì)、瀝青質(zhì)和堿氮質(zhì)量分數(shù),%;θ為催化劑的相對停留時間,s;φC/O為劑/油質(zhì)量比。

    芳烴吸附的影響一般包括芳烴、膠質(zhì)和瀝青質(zhì),堿氮毒害的影響一般還需要考慮催化劑的相對停留時間和劑/油質(zhì)量比。有學(xué)者認為,如果劑/油質(zhì)量比較高或者原料堿氮質(zhì)量分數(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ù)的辦法進行修正。隨著集總劃分得越來越細致,模型參數(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)速率方程組表達式,按照相關(guān)求解和優(yōu)化算法編制程序,求取最合適的模型動力學(xué)參數(shù)的過程。隨著集總動力學(xué)模型的發(fā)展,各種參數(shù)估算方法不斷產(chǎn)生,如Gauss-Newton法、Marquart法及其改進型算法、單純形法、Powell法、共軛梯度法、變尺度法(B-F-G-S)、粒子群優(yōu)化算法PSO、遺傳算法等。這些方法各具特點,通常從收斂性、穩(wěn)定性、可靠性、有效性、適用性、易用性、可移植性等方面進行比較后按需選用。此外,由于模型的集總劃分越來越細致,導(dǎo)致需要估算的動力學(xué)參數(shù)數(shù)量龐大。為了確保估算得到的參數(shù)的可靠性,往往一次估算的參數(shù)數(shù)量不能太多,因而需要分步或者分層估算。這樣得到的動力學(xué)參數(shù)既不失準確性,又能高效快速地獲取。

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

    目前,由于環(huán)境法規(guī)的日益嚴格和石油化工行業(yè)精細化、多方向的快速發(fā)展,市場對燃料和低碳烯烴從量和質(zhì)2個方面都提出了要求。因此,集總動力學(xué)模型在集總劃分方面更細致并更具工藝特色是勢在必行。事實上,各種被稱為后集總動力學(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)用效果已得到廣泛認可。最后,集總反應(yīng)模型與催化劑性質(zhì)的耦合將使催化裂化模型同時適用于多種催化劑,從而極大地拓寬其應(yīng)用范圍。因此,集總反應(yīng)模型與流動、傳熱、傳質(zhì)等模型和催化劑性質(zhì)的耦合考察,也應(yīng)該繼續(xù)重點研究。

    [1] 許友好. 我國催化裂化工藝技術(shù)進展[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

    猜你喜歡
    集總蠟油重油
    基于撕裂法的變壓器繞組集總參數(shù)等效電路頻率響應(yīng)計算方法
    常減壓裝置蠟油熱供系統(tǒng)改造研究與應(yīng)用
    價值工程(2022年19期)2022-06-14 03:24:02
    重油加氫處理催化劑級配方法和重油加氫處理方法
    能源化工(2021年6期)2021-12-30 15:41:26
    船用調(diào)質(zhì)重油燃燒及減排技術(shù)的研究進展
    能源工程(2021年1期)2021-04-13 02:06:02
    船舶使用850cSt重油的設(shè)計研究
    蠟油加氫脫硫集總動力學(xué)模型研究
    化工進展(2020年4期)2020-05-08 10:23:50
    原料加氫預(yù)處理對催化裝置的影響
    化工管理(2015年14期)2015-12-21 01:26:32
    一種加載集總器件的可調(diào)三維周期結(jié)構(gòu)
    焦化蠟油絡(luò)合脫氮-催化裂化優(yōu)化加工工藝技術(shù)
    久久久精品94久久精品| 老熟妇乱子伦视频在线观看 | 97人妻天天添夜夜摸| 超色免费av| 中文字幕制服av| 黄色片一级片一级黄色片| 久久精品aⅴ一区二区三区四区| 亚洲av成人一区二区三| 国产成人影院久久av| 搡老岳熟女国产| 一级毛片女人18水好多| 亚洲欧美一区二区三区黑人| 午夜精品国产一区二区电影| 成人国产av品久久久| 免费在线观看影片大全网站| 一级片免费观看大全| 韩国精品一区二区三区| 人人澡人人妻人| 午夜福利免费观看在线| 亚洲精品中文字幕在线视频| 久久ye,这里只有精品| 大陆偷拍与自拍| 手机成人av网站| 大码成人一级视频| 性高湖久久久久久久久免费观看| 亚洲精品成人av观看孕妇| 欧美日韩亚洲国产一区二区在线观看 | 久久久久视频综合| 欧美黑人精品巨大| 久久影院123| 国产一区二区激情短视频 | 精品少妇一区二区三区视频日本电影| 色婷婷av一区二区三区视频| 男女午夜视频在线观看| 美国免费a级毛片| 人人妻人人澡人人看| h视频一区二区三区| 97人妻天天添夜夜摸| 黑人猛操日本美女一级片| 欧美精品一区二区大全| 在线观看一区二区三区激情| 99香蕉大伊视频| 搡老岳熟女国产| 我要看黄色一级片免费的| 少妇人妻久久综合中文| 欧美日韩亚洲综合一区二区三区_| 亚洲全国av大片| 丝袜美足系列| 国产97色在线日韩免费| 色94色欧美一区二区| 亚洲欧美激情在线| 亚洲情色 制服丝袜| 亚洲欧美成人综合另类久久久| 亚洲男人天堂网一区| 丰满饥渴人妻一区二区三| 久久人人爽av亚洲精品天堂| 成人av一区二区三区在线看 | 午夜福利乱码中文字幕| 夜夜夜夜夜久久久久| 精品一品国产午夜福利视频| 中文字幕精品免费在线观看视频| 日韩制服骚丝袜av| 久久人妻福利社区极品人妻图片| 老汉色∧v一级毛片| 日本撒尿小便嘘嘘汇集6| 狠狠婷婷综合久久久久久88av| 成年动漫av网址| 日韩 亚洲 欧美在线| 50天的宝宝边吃奶边哭怎么回事| 精品国产乱码久久久久久小说| 日韩三级视频一区二区三区| 久久久久久久国产电影| 2018国产大陆天天弄谢| 在线永久观看黄色视频| 国产精品免费视频内射| 两个人看的免费小视频| 日本vs欧美在线观看视频| 欧美精品一区二区免费开放| 成人亚洲精品一区在线观看| 国产精品一区二区精品视频观看| 欧美亚洲 丝袜 人妻 在线| 高清视频免费观看一区二区| 在线观看www视频免费| 精品一区二区三区四区五区乱码| 久久人妻福利社区极品人妻图片| h视频一区二区三区| 精品少妇一区二区三区视频日本电影| 欧美另类一区| 久久精品亚洲av国产电影网| 国产av一区二区精品久久| 一区二区三区激情视频| 看免费av毛片| 中文字幕色久视频| 国产成人影院久久av| 亚洲男人天堂网一区| 日韩欧美国产一区二区入口| 久久99热这里只频精品6学生| 69av精品久久久久久 | 黄色a级毛片大全视频| 欧美精品人与动牲交sv欧美| 十八禁网站网址无遮挡| 手机成人av网站| 91老司机精品| 三级毛片av免费| 电影成人av| 亚洲男人天堂网一区| 女人高潮潮喷娇喘18禁视频| av免费在线观看网站| 一级毛片女人18水好多| 国产精品亚洲av一区麻豆| 最近中文字幕2019免费版| 丝袜喷水一区| 国产免费视频播放在线视频| 国产成人欧美在线观看 | 不卡av一区二区三区| 午夜免费鲁丝| 欧美xxⅹ黑人| 不卡一级毛片| 女人爽到高潮嗷嗷叫在线视频| 高清av免费在线| 久久精品成人免费网站| 国产三级黄色录像| 97在线人人人人妻| 色婷婷久久久亚洲欧美| 亚洲激情五月婷婷啪啪| 一二三四社区在线视频社区8| 国产一级毛片在线| 一级毛片女人18水好多| 淫妇啪啪啪对白视频 | 啦啦啦在线免费观看视频4| 成年女人毛片免费观看观看9 | 亚洲午夜精品一区,二区,三区| 三级毛片av免费| 免费在线观看完整版高清| 亚洲av电影在线观看一区二区三区| 欧美日本中文国产一区发布| 亚洲全国av大片| 丰满迷人的少妇在线观看| 69精品国产乱码久久久| 亚洲国产精品999| 在线观看一区二区三区激情| 男人舔女人的私密视频| 国产男人的电影天堂91| 国产成人精品无人区| av视频免费观看在线观看| 一边摸一边做爽爽视频免费| 国产成人av激情在线播放| 免费高清在线观看视频在线观看| 国产精品久久久人人做人人爽| 精品卡一卡二卡四卡免费| 99香蕉大伊视频| 91av网站免费观看| 国产xxxxx性猛交| 下体分泌物呈黄色| 久久久久国产一级毛片高清牌| 性高湖久久久久久久久免费观看| 天天躁夜夜躁狠狠躁躁| 日日爽夜夜爽网站| 菩萨蛮人人尽说江南好唐韦庄| 国产欧美亚洲国产| 涩涩av久久男人的天堂| 亚洲va日本ⅴa欧美va伊人久久 | 在线观看人妻少妇| 色视频在线一区二区三区| 男人操女人黄网站| 女人被躁到高潮嗷嗷叫费观| 国产片内射在线| 人成视频在线观看免费观看| 男人添女人高潮全过程视频| 下体分泌物呈黄色| 人人澡人人妻人| 欧美日韩视频精品一区| videosex国产| 欧美在线黄色| 国产高清国产精品国产三级| 精品免费久久久久久久清纯 | 嫁个100分男人电影在线观看| 国产在线视频一区二区| 日韩中文字幕欧美一区二区| 久久国产亚洲av麻豆专区| 久久精品久久久久久噜噜老黄| 久久久精品国产亚洲av高清涩受| 久久久精品国产亚洲av高清涩受| 制服诱惑二区| 别揉我奶头~嗯~啊~动态视频 | 国产一区二区三区在线臀色熟女 | 999久久久国产精品视频| 在线十欧美十亚洲十日本专区| 美女主播在线视频| 在线观看免费日韩欧美大片| 欧美激情高清一区二区三区| 老汉色av国产亚洲站长工具| 1024香蕉在线观看| videos熟女内射| 国产成人欧美| 激情视频va一区二区三区| 下体分泌物呈黄色| 国产一区二区 视频在线| 大码成人一级视频| 亚洲精品国产一区二区精华液| 亚洲国产av新网站| 亚洲性夜色夜夜综合| 亚洲精品美女久久av网站| 免费在线观看完整版高清| 99国产综合亚洲精品| 午夜福利在线观看吧| 国产麻豆69| 中文字幕最新亚洲高清| 9191精品国产免费久久| 国产精品国产三级国产专区5o| 精品亚洲乱码少妇综合久久| 国产精品久久久人人做人人爽| 国产精品九九99| 男女之事视频高清在线观看| tocl精华| av天堂久久9| 亚洲国产欧美网| 国产无遮挡羞羞视频在线观看| 一区二区日韩欧美中文字幕| 亚洲欧美一区二区三区黑人| 99热国产这里只有精品6| 夫妻午夜视频| 亚洲成人手机| 亚洲全国av大片| 男人添女人高潮全过程视频| 女人久久www免费人成看片| 免费高清在线观看视频在线观看| 十分钟在线观看高清视频www| 真人做人爱边吃奶动态| 亚洲精品中文字幕在线视频| 亚洲国产看品久久| 51午夜福利影视在线观看| 国产成人免费观看mmmm| 啦啦啦视频在线资源免费观看| 国产精品久久久久久精品古装| 蜜桃在线观看..| 午夜福利在线观看吧| 久久久久久亚洲精品国产蜜桃av| 91九色精品人成在线观看| 亚洲精品中文字幕在线视频| 日韩精品免费视频一区二区三区| 亚洲七黄色美女视频| 各种免费的搞黄视频| 亚洲国产欧美网| 国产欧美日韩一区二区三区在线| 久久av网站| 丝袜人妻中文字幕| 亚洲精品国产色婷婷电影| 黑人巨大精品欧美一区二区mp4| 日韩视频在线欧美| 国产91精品成人一区二区三区 | 欧美97在线视频| 亚洲欧美一区二区三区黑人| 中文精品一卡2卡3卡4更新| 丝袜美足系列| 91成人精品电影| videosex国产| 极品人妻少妇av视频| 在线观看免费高清a一片| 亚洲中文字幕日韩| 久久亚洲精品不卡| 国产免费现黄频在线看| 成人黄色视频免费在线看| 日本av手机在线免费观看| 欧美激情久久久久久爽电影 | 中亚洲国语对白在线视频| 国产极品粉嫩免费观看在线| 多毛熟女@视频| 在线观看舔阴道视频| 日韩一区二区三区影片| 国产xxxxx性猛交| 美女高潮到喷水免费观看| 男女免费视频国产| 久久香蕉激情| 正在播放国产对白刺激| 午夜福利影视在线免费观看| 777久久人妻少妇嫩草av网站| 十八禁高潮呻吟视频| 99re6热这里在线精品视频| 国产成人精品久久二区二区91| 99久久国产精品久久久| 搡老乐熟女国产| 99国产精品一区二区三区| 老熟女久久久| 99国产精品免费福利视频| 亚洲七黄色美女视频| 亚洲五月婷婷丁香| 老司机午夜福利在线观看视频 | 亚洲人成电影观看| 两性夫妻黄色片| 国产亚洲欧美精品永久| 亚洲av电影在线观看一区二区三区| 超碰成人久久| 国产精品影院久久| 捣出白浆h1v1| 欧美亚洲 丝袜 人妻 在线| 亚洲,欧美精品.| 国产成人精品无人区| 性少妇av在线| 啦啦啦中文免费视频观看日本| 精品少妇一区二区三区视频日本电影| 在线观看一区二区三区激情| 性高湖久久久久久久久免费观看| 久久精品亚洲av国产电影网| 欧美一级毛片孕妇| 欧美乱码精品一区二区三区| 中文字幕制服av| 99精品欧美一区二区三区四区| 亚洲一码二码三码区别大吗| bbb黄色大片| 制服人妻中文乱码| 亚洲自偷自拍图片 自拍| 成年av动漫网址| 51午夜福利影视在线观看| 一本—道久久a久久精品蜜桃钙片| 国产精品一区二区免费欧美 | 欧美激情极品国产一区二区三区| 一区二区三区精品91| 国产男女内射视频| 夜夜骑夜夜射夜夜干| 十八禁高潮呻吟视频| av有码第一页| 中文字幕人妻熟女乱码| 国产99久久九九免费精品| 妹子高潮喷水视频| 人妻 亚洲 视频| 在线天堂中文资源库| 久久女婷五月综合色啪小说| 久9热在线精品视频| 亚洲av电影在线观看一区二区三区| 欧美日韩精品网址| 午夜福利一区二区在线看| 精品国产乱码久久久久久男人| 两人在一起打扑克的视频| 正在播放国产对白刺激| 手机成人av网站| 亚洲第一av免费看| 国产精品.久久久| 男人操女人黄网站| 黄片播放在线免费| 黄网站色视频无遮挡免费观看| 久久精品aⅴ一区二区三区四区| 51午夜福利影视在线观看| 首页视频小说图片口味搜索| 国产成+人综合+亚洲专区| 久久久国产欧美日韩av| 久久人妻熟女aⅴ| av又黄又爽大尺度在线免费看| 亚洲va日本ⅴa欧美va伊人久久 | 亚洲av日韩精品久久久久久密| 日韩有码中文字幕| 一区二区三区四区激情视频| 悠悠久久av| 在线精品无人区一区二区三| 国产黄频视频在线观看| 一级毛片女人18水好多| 国产日韩欧美视频二区| 久久毛片免费看一区二区三区| 又大又爽又粗| 久久综合国产亚洲精品| 亚洲三区欧美一区| 精品乱码久久久久久99久播| 黄片大片在线免费观看| av一本久久久久| 欧美成人午夜精品| 在线精品无人区一区二区三| 麻豆国产av国片精品| 啦啦啦视频在线资源免费观看| av线在线观看网站| 少妇裸体淫交视频免费看高清 | 欧美精品亚洲一区二区| 久久亚洲精品不卡| 啪啪无遮挡十八禁网站| 欧美精品av麻豆av| 在线观看www视频免费| 黄色视频不卡| 国产一区二区在线观看av| 欧美性长视频在线观看| 欧美黑人精品巨大| 亚洲一区二区三区欧美精品| 欧美久久黑人一区二区| 欧美 日韩 精品 国产| 美女高潮到喷水免费观看| 一区二区av电影网| 亚洲欧美清纯卡通| 精品国产一区二区三区四区第35| 国产精品久久久久久人妻精品电影 | 国产成人欧美| 亚洲中文av在线| 99热网站在线观看| 一本大道久久a久久精品| 男女无遮挡免费网站观看| 免费高清在线观看视频在线观看| 午夜免费鲁丝| 久久精品国产a三级三级三级| 中文字幕另类日韩欧美亚洲嫩草| av网站在线播放免费| 亚洲欧美日韩另类电影网站| 十八禁网站免费在线| 精品国产一区二区久久| 免费在线观看视频国产中文字幕亚洲 | 在线永久观看黄色视频| 18禁裸乳无遮挡动漫免费视频| 亚洲 国产 在线| 亚洲一区二区三区欧美精品| 一级毛片精品| 蜜桃国产av成人99| 91成年电影在线观看| 黄片播放在线免费| 午夜福利一区二区在线看| 亚洲免费av在线视频| 在线观看www视频免费| 亚洲av日韩在线播放| 免费一级毛片在线播放高清视频 | 久久精品亚洲熟妇少妇任你| 日本精品一区二区三区蜜桃| 自拍欧美九色日韩亚洲蝌蚪91| 97精品久久久久久久久久精品| 久久久久精品人妻al黑| 一级片'在线观看视频| 午夜精品国产一区二区电影| 丝袜人妻中文字幕| 免费日韩欧美在线观看| 国产欧美亚洲国产| 欧美另类一区| 丝袜喷水一区| 久久精品久久久久久噜噜老黄| 国产视频一区二区在线看| 国产一区二区三区在线臀色熟女 | 久久热在线av| 久久精品熟女亚洲av麻豆精品| 国产区一区二久久| 韩国高清视频一区二区三区| 男女之事视频高清在线观看| 高潮久久久久久久久久久不卡| 波多野结衣av一区二区av| 嫁个100分男人电影在线观看| 亚洲av成人一区二区三| 久久久水蜜桃国产精品网| 熟女少妇亚洲综合色aaa.| 老熟妇仑乱视频hdxx| 免费女性裸体啪啪无遮挡网站| 精品国内亚洲2022精品成人 | 一进一出抽搐动态| 亚洲欧美日韩另类电影网站| 少妇裸体淫交视频免费看高清 | 亚洲精品av麻豆狂野| 午夜福利影视在线免费观看| 国产欧美日韩精品亚洲av| 亚洲av男天堂| 97精品久久久久久久久久精品| 色94色欧美一区二区| 日韩一卡2卡3卡4卡2021年| 国产精品99久久99久久久不卡| 亚洲精品粉嫩美女一区| 午夜久久久在线观看| 亚洲专区中文字幕在线| 啪啪无遮挡十八禁网站| 精品少妇内射三级| 国产精品国产av在线观看| 夜夜夜夜夜久久久久| 亚洲熟女毛片儿| 免费观看av网站的网址| 50天的宝宝边吃奶边哭怎么回事| 国产亚洲av片在线观看秒播厂| 国产一区二区三区av在线| 女人爽到高潮嗷嗷叫在线视频| 午夜福利,免费看| 777米奇影视久久| 18禁观看日本| 亚洲激情五月婷婷啪啪| 国产亚洲精品一区二区www | 国产片内射在线| 亚洲国产毛片av蜜桃av| 一区在线观看完整版| 国产精品99久久99久久久不卡| 欧美日韩中文字幕国产精品一区二区三区 | 午夜久久久在线观看| 精品亚洲成国产av| 亚洲一区二区三区欧美精品| 亚洲欧美清纯卡通| 精品人妻熟女毛片av久久网站| 亚洲国产中文字幕在线视频| 日韩欧美国产一区二区入口| 日韩大码丰满熟妇| 在线观看舔阴道视频| 精品少妇黑人巨大在线播放| 法律面前人人平等表现在哪些方面 | 国产av一区二区精品久久| 国产成人免费观看mmmm| 久久久久久免费高清国产稀缺| 最近最新中文字幕大全免费视频| 精品高清国产在线一区| 国产一卡二卡三卡精品| 国产男人的电影天堂91| 欧美+亚洲+日韩+国产| 欧美精品一区二区免费开放| 手机成人av网站| 国产免费福利视频在线观看| 男男h啪啪无遮挡| 精品人妻熟女毛片av久久网站| 国产片内射在线| 国产又爽黄色视频| 久久久久久亚洲精品国产蜜桃av| 每晚都被弄得嗷嗷叫到高潮| 免费高清在线观看视频在线观看| 欧美精品啪啪一区二区三区 | 久久久国产成人免费| 一级a爱视频在线免费观看| 久久人妻熟女aⅴ| 亚洲中文日韩欧美视频| 欧美日本中文国产一区发布| 高清在线国产一区| 久久精品国产亚洲av香蕉五月 | 啦啦啦免费观看视频1| 日本wwww免费看| 黄色 视频免费看| 91国产中文字幕| 曰老女人黄片| 日本a在线网址| 国产成人精品久久二区二区免费| 国产精品二区激情视频| 夫妻午夜视频| 啦啦啦免费观看视频1| 午夜老司机福利片| 亚洲精品在线美女| 亚洲一区二区三区欧美精品| av不卡在线播放| 欧美激情高清一区二区三区| 国产成人免费无遮挡视频| 日韩三级视频一区二区三区| 久久午夜综合久久蜜桃| 80岁老熟妇乱子伦牲交| 大片电影免费在线观看免费| 巨乳人妻的诱惑在线观看| 国产男人的电影天堂91| 日韩 欧美 亚洲 中文字幕| 国产成人欧美在线观看 | 精品国产一区二区久久| 少妇粗大呻吟视频| 老司机靠b影院| 一区二区三区精品91| 国产精品免费视频内射| 狠狠狠狠99中文字幕| av电影中文网址| 国产精品一区二区在线观看99| 久久热在线av| 亚洲精品一二三| 91大片在线观看| 久久精品aⅴ一区二区三区四区| 亚洲综合色网址| 国产又爽黄色视频| 国产成人免费观看mmmm| 18禁观看日本| 精品一区二区三区四区五区乱码| 国产欧美日韩一区二区三区在线| 欧美日韩视频精品一区| 一本久久精品| 精品国产乱码久久久久久男人| 亚洲中文日韩欧美视频| 女性生殖器流出的白浆| 午夜福利视频精品| 日韩熟女老妇一区二区性免费视频| 一区二区三区四区激情视频| 十八禁网站免费在线| 国精品久久久久久国模美| 日韩熟女老妇一区二区性免费视频| 天堂8中文在线网| 在线观看免费视频网站a站| 少妇人妻久久综合中文| 丝瓜视频免费看黄片| 亚洲一卡2卡3卡4卡5卡精品中文| 国产无遮挡羞羞视频在线观看| 免费人妻精品一区二区三区视频| 欧美成人午夜精品| 青草久久国产| 欧美日韩黄片免| 成年美女黄网站色视频大全免费| 又黄又粗又硬又大视频| 天堂中文最新版在线下载| 国产精品久久久久成人av| 国产男女超爽视频在线观看| 国产成人免费无遮挡视频| 亚洲中文av在线| 精品高清国产在线一区| 99精国产麻豆久久婷婷| 女人爽到高潮嗷嗷叫在线视频| 免费观看人在逋| 免费日韩欧美在线观看| 老鸭窝网址在线观看| 久久久久国产精品人妻一区二区| 最新的欧美精品一区二区| 亚洲中文字幕日韩| 亚洲熟女精品中文字幕| 少妇粗大呻吟视频| 一区二区av电影网| 日本五十路高清| 一本色道久久久久久精品综合| 天天影视国产精品| 午夜久久久在线观看| 欧美激情 高清一区二区三区| 亚洲精品久久午夜乱码| 欧美黑人欧美精品刺激| 久久国产亚洲av麻豆专区| 91av网站免费观看| 久久久久国产精品人妻一区二区| 嫁个100分男人电影在线观看| 老汉色∧v一级毛片| bbb黄色大片| 男女午夜视频在线观看| 国产男人的电影天堂91| 久久亚洲精品不卡| 亚洲国产精品一区三区| 久久久水蜜桃国产精品网| 精品少妇内射三级|