王俊杰,潘艷秋,牛亞賓,俞路
(大連理工大學化工學院,遼寧 大連 116024)
《中國制造2025》明確指出:“新一代信息技術(shù)與制造業(yè)深度融合,正在引發(fā)影響深遠的產(chǎn)業(yè)變革,形成新的生產(chǎn)方式、產(chǎn)業(yè)形態(tài)、商業(yè)模式和經(jīng)濟增長點”[1]。在此背景下,石化產(chǎn)業(yè)的數(shù)字化、智能化轉(zhuǎn)型升級進程加快。
連續(xù)催化重整裝置是高辛烷值汽油、芳香烴的核心生產(chǎn)裝置[2],其高質(zhì)量運行對石化企業(yè)來說意義非凡且愈加重要。然而,由于催化重整反應體系復雜,組分數(shù)和反應眾多[3],導致通過常規(guī)方法構(gòu)建催化重整過程的動力學模型存在巨大困難。進入21 世紀后,隨著國家環(huán)保法規(guī)的不斷升級,以往僅控制石油產(chǎn)品整體平均理化性質(zhì)就能保障石油產(chǎn)品使用性能的傳統(tǒng)模式被打破,石油產(chǎn)品質(zhì)量標準已由餾分級上升到了分子級[4],目前常用的集總動力學模型[5]已無法滿足對催化重整產(chǎn)物進行分子水平預測的要求。因此,建立分子水平動力學模型勢在必行。隨著計算機性能的不斷提高和現(xiàn)代分析技術(shù)的顯著進步[6],石油分子工程[7]概念應運而生,使得分子水平催化重整動力學模型建立成為可能。
催化重整裝置分子水平模型構(gòu)建主要包括兩方面內(nèi)容:石腦油分子重構(gòu)模型和催化重整分子水平反應動力學模型。建立石腦油分子重構(gòu)模型的目的是通過一些石腦油常見物性作為輸入條件來計算石腦油餾分的詳細組成分布。對于石腦油這種組分種類復雜的混合物,構(gòu)建準確預測組分含量的分子重構(gòu)模型十分困難,目前有關(guān)研究主要集中在兩方面:一是調(diào)整預設分子庫中的分子數(shù)量;二是改進模型構(gòu)建方法。Van Geem 等[8]按質(zhì)量分數(shù)排序從173 種純組分中選出37種代表性組分,以熵值最大化為目標函數(shù),基于石腦油宏觀物性調(diào)整這些石腦油分子的摩爾分數(shù),建立石腦油分子重構(gòu)模型。Wang等[9]使用MTHS矩陣模型預測石腦油和汽油樣品的分子組成,用兩參數(shù)正態(tài)分布代替三參數(shù)伽馬分布,描述每個同系物內(nèi)部的統(tǒng)計分布,以更少的變量準確地預測了石腦油組成??傮w看來,目前石腦油分子重構(gòu)模型的發(fā)展相對成熟,但是作為保障模型計算準確性的最基礎(chǔ)工作,分子庫中分子選擇原則卻仍未明確提出,且雖然石腦油所含重組分較少,但幾十甚至上百個組分含量的優(yōu)化求解仍有很大難度,構(gòu)建高效的分子重構(gòu)模型仍是重中之重。
對于催化重整分子水平動力學模型構(gòu)建,研究者們嘗試了多種方法。周齊宏等[10]利用MTHS 矩陣模型描述重整進料組成,建立催化重整過程的分子水平模型,并以裝置經(jīng)濟效益最大化為目標來優(yōu)化裝置操作條件。Sotelo-Boyás 等[11]通過向量和布爾關(guān)系矩陣生成分子反應網(wǎng)絡,采用列文伯格-馬夸爾特法估計速率系數(shù)并將動力學模型應用于擬均相和非均相反應過程,模擬了具有向心徑向流動的工業(yè)三床絕熱催化重整爐,結(jié)果證明了模型的可靠性。王睿通等[12]基于結(jié)構(gòu)導向集總方法,建立裝置分子水平反應動力學模型,分析了反應溫度、壓力、空速的影響,實現(xiàn)了對于催化重整裝置中芳烴以及液體收率的雙目標優(yōu)化。Zhou等[6]基于KMT軟件建立了連續(xù)催化重整的分子動力學模型,對工業(yè)裝置的模擬結(jié)果證明了模型的可靠性。綜上所述,采用MTHS 矩陣、布爾關(guān)系矩陣、結(jié)構(gòu)導向集總、KMT 軟件等方法可較好地構(gòu)建催化重整的分子水平模型。但由于催化重整裝置的復雜性,上述文獻中仍存在僅對產(chǎn)物宏觀物性進行預測、僅對典型產(chǎn)物分子進行預測等問題,因此目前催化重整裝置分子水平模型的建立方法和模型本身均需要完善和提高。
本文基于國內(nèi)某石化企業(yè)智慧工廠建設需求,構(gòu)建基于工廠生產(chǎn)實際數(shù)據(jù)的催化重整分子水平反應動力學模型,探討催化重整裝置分子水平建模的思路和方法,實現(xiàn)通過石腦油宏觀物性預測催化重整產(chǎn)物分布的目的。
圖1為石腦油分子重構(gòu)模型的構(gòu)建思路。
圖1 石腦油分子重構(gòu)模型構(gòu)建思路
1.1.1 石腦油確定性分子庫構(gòu)建
構(gòu)建石腦油確定性分子庫的目的是確定構(gòu)成石腦油的純物質(zhì)種類及其物性。由于石腦油含有400余種組分,故目前石腦油分子重構(gòu)模型大多采用確定性方法,即選擇出的代表性分子在真實油品中實際存在且是在建模之初就確定好,不再隨模型的求解過程發(fā)生變化。
本文的背景裝置以石腦油為原料,采用連續(xù)重整工藝技術(shù)生產(chǎn)汽油、苯、混合二甲苯并副產(chǎn)氫氣?;谘b置特點,確定了分子選擇原則,即同時考慮到石腦油餾分的烴類組成及催化重整反應特性,選擇石腦油餾分中質(zhì)量分數(shù)高、參與催化重整反應以及用于評價催化重整原料反應性能的分子構(gòu)成石腦油確定性分子庫。詳細內(nèi)容如下。
(1)石腦油餾分烴類組成 石腦油餾分的主要成分是正構(gòu)烷烴和支鏈數(shù)少的異構(gòu)烷烴,其中只含有一個甲基支鏈的異構(gòu)烷烴含量占異構(gòu)烷烴總含量的一半以上;環(huán)烷烴通常只包含單環(huán)環(huán)烷烴,有環(huán)戊烷系化合物與環(huán)己烷系化合物兩類;芳香烴只含有單環(huán)芳香烴;烯烴的含量非常低但其同分異構(gòu)體眾多,僅選取少量常見的烯烴分子[13]。
(2)催化重整過程反應類型 包括脫氫反應、脫氫環(huán)化、異構(gòu)化反應以及加氫裂化[14]。需著重考慮參與這些反應的分子。
(3)重整原料反應性能 用“芳烴潛含量”來表征重整原料的反應性能,即當原料中的環(huán)烷烴全部轉(zhuǎn)化為芳烴時所獲得的芳烴量。需著重考慮計算重整原料的反應性能的分子。
從三個角度提出分子庫選擇原則,不僅保障分子重構(gòu)模型的可靠性,同時可延伸應用于后續(xù)動力學模型的建立,這種聯(lián)合構(gòu)建分子重構(gòu)模型-反應動力學模型的研究思路,可用于不同餾分油和不同生產(chǎn)裝置,提供出有效構(gòu)建分子水平動力學模型的思路。
1.1.2 集總組分物性計算
以上建立的分子庫中包含眾多分子,若直接求解各分子含量存在困難。為了減少變量數(shù),本文先采用MTHS 矩陣模型[15]描述石腦油餾分的組成及含量,即用矩陣表示石油餾分的組成,其中矩陣的行代表碳原子數(shù)目、列代表同系物(族組成),每個單元格代表對應碳數(shù)和分子類型的單個分子或同分異構(gòu)體集總的質(zhì)量分數(shù);再假設矩陣單元格集總中的同分異構(gòu)體處于熱力學平衡狀態(tài)(每個矩陣單元格的自由度都將降為1),借助Aspen Plus軟件求得每個組分的熱力學平衡分數(shù);最后利用混合規(guī)則[16],結(jié)合純物質(zhì)含量和物性,計算出矩陣單元格中每個集總組分的平均物性[17]。
1.1.3 分子重構(gòu)數(shù)學模型構(gòu)建
基于以上求取的集總組分物性,利用集總單元格的質(zhì)量分數(shù),即可求得石腦油混合物的物性。
上節(jié)將模型變量數(shù)從幾百降維到幾十,若同時調(diào)節(jié)這些變量存在耗時長等困難,故本節(jié)引入統(tǒng)計分布進一步減少變量數(shù)。假設集總單元格的質(zhì)量分數(shù)隨碳數(shù)遵循伽馬分布、同時可由分布函數(shù)得到分子重構(gòu)模型的約束條件,如式(1)、式(2)。
式(4)即為優(yōu)化模型的約束之一。另外,將完全滿足正構(gòu)烷烴、異構(gòu)烷烴、烯烴、環(huán)烷烴、芳香烴含量模擬值和實際值一致作為分子重構(gòu)模型的第二個約束條件,如式(5)所示。
為了保證石腦油分子重構(gòu)模型的準確性,本文以石腦油宏觀物性的模擬值與實際值的相對偏差的平方和為目標函數(shù),如式(6)所示。
分子水平反應動力學模型建立及求解思路如圖2所示。
圖2 分子水平反應動力學模型建立及求解思路
1.2.1 催化重整反應網(wǎng)絡生成
本文結(jié)合背景裝置情況、基于文獻[14, 18],提出了如下反應網(wǎng)絡簡化原則。
(1)烷烴環(huán)化 環(huán)化反應主要是烷烴直接環(huán)化,烷烴經(jīng)過五元環(huán)烷烴生成六元環(huán)烷烴,進而脫氫生成芳烴。
(2)芳烴反應 對于芳烴加氫裂化,只考慮芳烴側(cè)鏈的裂化;二甲苯和乙苯,三甲苯、甲乙苯和丙苯之間不相互轉(zhuǎn)化(轉(zhuǎn)化速率非常慢)。
(3)環(huán)烷烴反應 因環(huán)烷烴的異構(gòu)化和脫氫環(huán)化速率很快,忽略其裂化反應。
(4)烯烴反應 石腦油中烯烴含量極低,故不專門為其生成反應網(wǎng)絡。
(5)反應的可逆性 將異構(gòu)化、脫氫環(huán)化和芳構(gòu)化視為可逆反應,而將芳烴側(cè)鏈、烷烴加氫裂化視為不可逆反應。
對于由大量分子通過復雜的反應網(wǎng)絡相互關(guān)聯(lián)組成的復雜反應系統(tǒng),可利用自動反應網(wǎng)絡生成器來實現(xiàn)反應網(wǎng)絡構(gòu)建。目前已研發(fā)使用的自動反應網(wǎng)絡生成器有多種[19],其中RING 已被應用于生成單種物質(zhì)和簡單混合物的反應網(wǎng)絡[20-22],通過該軟件能夠在較短時間、人力消耗較少的情況下得到詳細的反應網(wǎng)絡。
1.2.2 分子水平催化重整反應動力學建模
建立分子水平動力學模型,包括動力學方程組構(gòu)建(將確定的反應網(wǎng)絡信息轉(zhuǎn)化為可求解的數(shù)學表達式)、反應速率方程確定兩部分內(nèi)容。
(1)動力學方程組構(gòu)建 以正辛烷為例,圖3和圖4分別為其局部反應網(wǎng)絡以及反應轉(zhuǎn)化為動力學常微分方程的流程。首先要獲得每一個反應的化學反應表達式;其次為提取反應中的分子與其化學計量數(shù)構(gòu)建轉(zhuǎn)化矩陣,矩陣的行數(shù)和列數(shù)分別表示分子個數(shù)和反應個數(shù);再次,將矩陣與每個反應的反應速率(r)相乘得到每個物質(zhì)的動力學方程。其他組分的動力學方程構(gòu)建與此類似。
圖3 正辛烷的局部反應網(wǎng)絡
圖4 反應網(wǎng)絡轉(zhuǎn)化為動力學方程流程圖
(2)反應速率方程確定 催化重整過程中各反應相對于烴組分呈簡單的一級反應,且與氫分壓呈指數(shù)關(guān)系[23]。但對于包含上萬個反應的催化重整反應網(wǎng)絡,求解前需要對速率常數(shù)進行降維處理。借鑒集總理論,認為動力學特性相似的一類組分的速率常數(shù)相同[5],將反應網(wǎng)絡中的分子集總到MTHS矩陣中,因此各反應反應速率常數(shù)與集總構(gòu)成及反應類型有關(guān)。在此基礎(chǔ)上,使用阿倫尼烏斯定律表示速率常數(shù)與溫度的關(guān)系,同時考慮其與壓力呈指數(shù)關(guān)系,計算如式(7)所示。
對反應頻率因子進行估計,反應活化能和壓力指數(shù)采用文獻[24]實驗結(jié)果,文獻中缺少的二甲苯轉(zhuǎn)化的18個參數(shù)采用文獻[25]的方法進行估算。
利用遺傳算法求解式(7)中的動力學參數(shù),目標函數(shù)為主要產(chǎn)物的質(zhì)量分數(shù)模擬值和實際值之差的平方和,如式(8)。
根據(jù)1.1 節(jié)中的石腦油確定性分子庫構(gòu)建原則,構(gòu)建包括C4~C12范圍內(nèi)的共270 個分子的石腦油分子數(shù)據(jù)庫。其中,正構(gòu)烷烴和異構(gòu)烷烴均包含C4~C12,烯 烴、環(huán) 烷 烴、芳 烴 分 別 包 含C4~C7、C5~C11、C6~C11。石腦油確定性分子庫中的分子種類、純物質(zhì)物性數(shù)據(jù)均可在NIST(https://webbook.nist.gov/)上獲得。
將上述270 個分子集總到MTHS 矩陣模型中,得到圖5(√和×分別表示在相應的碳數(shù)和同源物系列下有、沒有集總組分)的共計35個集總的矩陣。結(jié)合1.1 節(jié)中的方法計算出石腦油宏觀物性,目標函數(shù)式(6)中的權(quán)重因子設置為5%、10%、20%、40%蒸餾體積時為0.01,其他均為1。利用遺傳算法即可求得石腦油的詳細分子組成。
圖5 石腦油餾分MTHS矩陣
在 MATLAB 軟件平臺中求解石腦油分子重構(gòu)模型,模擬得到的蒸餾體積數(shù)據(jù)與實際值的比值如圖6所示??梢钥闯?,模擬蒸餾體積數(shù)據(jù)與實際數(shù)據(jù)基本吻合。表1中列出分子量和密度的模擬值與實驗值的比較,相對誤差均在0.5%以內(nèi)。
表1 物性實際值和模擬值的比較
圖6 模擬蒸餾體積數(shù)據(jù)與實際數(shù)據(jù)的比較
圖7顯示了該石腦油樣品分子組成的模擬值和實際值,可見正構(gòu)烷烴、異構(gòu)烷烴、環(huán)烷烴、芳香烴對應碳數(shù)的質(zhì)量分數(shù)與實際值總體上吻合較好,說明模型計算精度較高,可以在工廠實際中用來計算石腦油的詳細組成。
圖7 石腦油樣品分子組成的實際值和模擬值
基于1.2 節(jié)確定的催化重整反應類型以及反應網(wǎng)絡簡化規(guī)則,利用RING 軟件編寫代碼,RING將反應規(guī)則應用于分子及其生成的產(chǎn)物,從系統(tǒng)的初始反應物開始,結(jié)果是在反應規(guī)則的基礎(chǔ)上產(chǎn)生所有可能反應和物種的窮舉網(wǎng)絡。本文生成反應網(wǎng)絡的代碼編寫規(guī)則遵循文獻中所寫,得到包含18種反應類型、865個分子和6616個反應的石腦油催化重整反應網(wǎng)絡,反應類型參見表2。
表2 石腦油分子水平催化重整反應動力學模型中反應類型統(tǒng)計
依據(jù)獲得的反應網(wǎng)絡構(gòu)建分子水平催化重整反應動力學模型,按理想的平推流反應器建立物料衡算,采用工廠實際數(shù)據(jù)對反應動力學參數(shù)進行估計,結(jié)果見表3。該工廠生產(chǎn)工況較為穩(wěn)定,反應溫為532℃,反應壓力為0.35MPa,液時空速為1.58h-1,產(chǎn)物組成變化很小,故采用本文作者課題組[26]前期對相同背景裝置的模擬結(jié)果,來驗證模型在不同工況下模擬結(jié)果的可靠性,圖8為反應溫度分別為525℃、530℃、532℃、535℃、反應壓力為0.35MPa、液時空速為1.58h-1時27種產(chǎn)物組分的質(zhì)量分數(shù)模擬值與實際值的對比,模擬值與實際值的絕對誤差均在0.85%以內(nèi),模型對于芳香烴收率的模擬相對誤差小于0.35%,對于乙苯、對二甲苯、鄰二甲苯、間二甲苯的模擬絕對誤差小于0.5%,說明模型對于催化重整產(chǎn)物分布的預測可靠性較高。
表3 反應動力學參數(shù)估算值
圖8 產(chǎn)物組分的質(zhì)量分數(shù)實際值與模擬值比較
由建立的催化重整分子水平反應動力學模型進行重整反應過程分析,重整過程中反應溫度以及反應壓力對芳烴以及C5+液體收率的影響如圖9所示,C5+液體收率隨溫度升高而減小,原因是溫度升高會加劇重整過程中裂化反應的發(fā)生,從而導致C5+液體收率降低,此時芳香烴收率升高,主要是因為脫氫芳構(gòu)化反應是吸熱反應,溫度升高有利于反應的進行,因此其收率提高。提高反應壓力不利于烷烴環(huán)化反應與環(huán)烷烴芳構(gòu)化反應的進行,故提升反應壓力會導致芳烴收率降低,與芳烴相似,反應壓力的升高會使得C5+液體收率降低。
圖9 反應溫度、反應壓力對芳香烴、C5+液體收率的影響
(1)針對某催化重整裝置情況,構(gòu)建了包含270個分子的石腦油確定性分子數(shù)據(jù)庫,將分子集總到MTHS模型中,通過優(yōu)化算法建立石腦油分子重構(gòu)模型。模擬結(jié)果表明,石腦油混合物的分子量和密度模擬值與實際值的相對誤差均在0.5%以內(nèi),蒸餾體積數(shù)據(jù)及石腦油分子組成的模擬值與實際值吻合良好。
(2)基于以上石腦油分子重構(gòu)模型,運用RING 軟件生成反應網(wǎng)絡,構(gòu)建分子水平催化重整反應動力學模型并采用遺傳算法估算模型參數(shù)。模擬結(jié)果與工廠實際值絕對誤差較小,證明了模型的可靠性。
本文建立了基于工廠數(shù)據(jù)的催化重整裝置分子級反應動力學模型,實現(xiàn)了由已知原料宏觀物性預測催化重整產(chǎn)物分布的目的,為從分子水平上認識裝置反應機理進而支撐催化重整裝置的智能化建設奠定基礎(chǔ)。
符號說明
A6——6個碳的烷烴、正構(gòu)烷烴、異構(gòu)烷烴、五元環(huán)環(huán)烷烴、芳香烴,其余類似
D——石腦油的密度,kg/m3
E——反應活化能,kJ/mol
IP6——6個碳的異構(gòu)烷烴,其余類似
k——反應速率常數(shù),s-1
k0——反應頻率因子,s-1·MPa-b
M——同系物n中的一系列組分
Mw——石腦油的摩爾質(zhì)量,g/mol
N——同系物
6N6,5N6——6 個碳的六元環(huán)烷烴、五元環(huán)環(huán)烷烴,其余類似
NP6——6個碳的正構(gòu)烷烴,其余類似
O——石腦油分子重構(gòu)模型目標函數(shù)
P5——5個碳的烷烴,其余類似
PH——氫分壓,MPa
PR——分子水平催化重整動力學模型目標函數(shù)
pr——產(chǎn)物j的質(zhì)量分數(shù),%
pm,n——同系物n中單元格(m,n)的概率密度
q——催化重整產(chǎn)物組分數(shù)目
R——摩爾氣體常量,8.3145J/(mol·K)
T——反應溫度,K
t——積分變量
Vt——實沸點蒸餾數(shù)據(jù),%
w——權(quán)重因子
Xn——正構(gòu)烷烴、異構(gòu)烷烴、烯烴、環(huán)烷烴、芳香烴的含量值,%
ym,n——通過歸一化pm,n得到的對應于同系物n中組分m的分數(shù)
Z——碳數(shù)
αn——第n個同系物伽馬分布的形狀參數(shù)
βn——第n個同系物伽馬分布的尺寸參數(shù)
上角標
act——實際值
b——壓力指數(shù)
sim——模擬值
下角標
i——反應頻率因子數(shù)目
j——催化重整產(chǎn)物組分數(shù)
m——對應于同系物n的組分m
min,max ——最小、最大
n——同系物N的組成
t——t=1~12,表示5%、10%、20%、30%、40%、50%、60%、70%、80%、90%、95%、98%的蒸餾體積數(shù)據(jù)