催化裂化回?zé)捰图託浞磻?yīng)動(dòng)力學(xué)模型

張 奎， 戴立順， 劉 濤， 邵志才， 聶 紅

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

催化裂化工藝是石油煉制過程中重油輕質(zhì)化的主要技術(shù)之一。為進(jìn)一步提高催化裂化處理能力和催化裂化裝置輕油收率，研究人員已開展關(guān)于催化裂化回?zé)捰?即重循環(huán)油)的高效利用研究。研究表明，催化裂化回?zé)捰途哂忻芏却蟆?碳比低、芳烴含量高的特點(diǎn)。如果將催化裂化回?zé)捰椭苯臃祷靥嵘茉俅瘟鸦?，含量較高的芳烴在發(fā)生裂化的同時(shí)也發(fā)生縮合反應(yīng)，導(dǎo)致大量焦炭的生成[1-3]，降低石油資源的利用效率。

鑒于此，許友好等[1]提出將催化裂化回?zé)捰拖冗M(jìn)行加氫預(yù)處理，再進(jìn)入提升管進(jìn)行催化裂化。具體工藝流程是：在一定加氫工藝條件下，將催化裂化回?zé)捰团c氫氣的混合物通過加氫催化劑床層，對(duì)催化裂化回?zé)捰瓦M(jìn)行加氫處理；加氫處理后，加氫油中大分子芳烴的含量降低，可作為催化裂化裝置的優(yōu)質(zhì)原料進(jìn)一步進(jìn)行催化裂化反應(yīng)；催化裂化產(chǎn)物分餾后，所得催化裂化回?zé)捰宛s分再循環(huán)進(jìn)入加氫預(yù)處理裝置。在催化裂化回?zé)捰图託溥^程中，通過加氫飽和反應(yīng)(HDA)后，芳烴特別是多環(huán)芳烴含量下降，而環(huán)烷烴等飽和烴含量增加；部分S、N元素通過加氫脫硫反應(yīng)(HDS)和加氫脫氮反應(yīng)(HDN)被脫除。因此，對(duì)催化裂化回?zé)捰瓦M(jìn)行循環(huán)加氫預(yù)處理后，催化裂化輕質(zhì)油總收率提高，生焦量降低，SOx、NOx的綜合排放量下降。

為進(jìn)一步提高對(duì)催化裂化加氫反應(yīng)的認(rèn)識(shí)，科研人員愈發(fā)重視重餾分油加氫反應(yīng)動(dòng)力學(xué)模型的相關(guān)研究。Ancheyta等[4-6]先對(duì)常壓蠟油加氫處理過程中的HDS、HDN、HDA和烯烴加氫等反應(yīng)動(dòng)力學(xué)進(jìn)行了研究，后對(duì)減壓蠟油HDS、HDN、加氫脫瀝青質(zhì)(HDAsph)和加氫脫殘?zhí)?HDCCR)等反應(yīng)進(jìn)行了動(dòng)力學(xué)反應(yīng)研究[7]。Dalai等[8]利用滴流床反應(yīng)器研究了H2純度對(duì)蠟油加氫反應(yīng)(HDS、HDN和HDA)動(dòng)力學(xué)的影響，結(jié)果表明各反應(yīng)級(jí)數(shù)均為擬一級(jí)。Verstraete等[9]采用8集總化學(xué)反應(yīng)動(dòng)力學(xué)揭示了固定床渣油加氫處理過程中C、H、S、N、O、Ni、V元素的反應(yīng)遷移規(guī)律。Ferreira等[10]將渣油中瀝青質(zhì)、膠質(zhì)和芳烴劃分為19個(gè)集總，研究了渣油加氫動(dòng)力學(xué)反應(yīng)規(guī)律；Mujtaba等[11]采用滴流床反應(yīng)器研究了全餾分油的HDAsph反應(yīng)動(dòng)力學(xué)。然而，烴類相關(guān)反應(yīng)(HDA、HDAsph和HDCCR等)的反應(yīng)動(dòng)力學(xué)模型較為復(fù)雜，且對(duì)催化裂化回?zé)捰偷募託浞磻?yīng)動(dòng)力學(xué)相關(guān)研究較少，因此亟需建立簡(jiǎn)單實(shí)用的催化裂化回?zé)捰图託浞磻?yīng)模型，以便深入研究該加氫反應(yīng)的動(dòng)力學(xué)。

鑒于此，筆者開展了催化裂化回?zé)捰图託涮幚碓囼?yàn)，研究了催化裂化回?zé)捰椭蠸、N及烴類隨加氫處理工藝條件變化而遷移的規(guī)律。將加氫反應(yīng)重新分為3類主要反應(yīng)：烴類反應(yīng)、HDS反應(yīng)和HDN反應(yīng)，構(gòu)建了催化裂化回?zé)捰图託浞磻?yīng)宏觀動(dòng)力學(xué)模型，預(yù)測(cè)了催化裂化回?zé)捰徒?jīng)加氫預(yù)處理后所得加氫油的C、H、S、N元素含量。研究為加氫預(yù)處理反應(yīng)器設(shè)計(jì)、催化劑選擇、裝置能耗計(jì)算，及優(yōu)化反應(yīng)條件提供數(shù)據(jù)支持與參考。

1 試驗(yàn)部分

1.1 原料油

試驗(yàn)采用催化裂化回?zé)捰虯作原料油，其芳烴、S、N及膠質(zhì)的質(zhì)量分?jǐn)?shù)均較高，主要性質(zhì)如表1所示。

表1 催化裂化回?zé)捰虯的主要性質(zhì)Table 1 Main properties of FCC recycle oil A

1.2 裝置、催化劑及試驗(yàn)條件

試驗(yàn)采用滴流床加氫反應(yīng)器裝置，如圖1所示。催化劑采用中國(guó)石化石油化工科學(xué)研究院(RIPP)研發(fā)的加氫保護(hù)劑和重油加氫精制催化劑RDA-1，其主要活性組分為Ni-W，催化劑裝填總體積為100 mL。試驗(yàn)操作條件：氫/油體積比800，氫分壓8～14 MPa，反應(yīng)溫度633.15～663.15 K，液時(shí)空速0.6～5.0 h-1。

1.3 分析方法及儀器

油品中碳、氫含量的測(cè)定采用元素分析法，分析儀器為元素分析儀(型號(hào)為vario EL Ⅲ，德國(guó)elementar公司)。測(cè)定過程中，樣品高溫燃燒，樣品中的C、H分別轉(zhuǎn)化為CO2、H2O，然后載氣攜帶反應(yīng)氣進(jìn)入色譜柱進(jìn)行分離，最后經(jīng)熱導(dǎo)池進(jìn)行檢測(cè)。

油品中硫含量的測(cè)定采用能量色散X射線熒光光譜法，分析儀器為能量色散X射線熒光光譜儀(型號(hào)為Epsilon3XLE，荷蘭PANalytical公司)。分析過程中，將樣品置于X射線束中，通過測(cè)定2.3 keV的硫Kα特征譜線強(qiáng)度，并與校準(zhǔn)樣品的強(qiáng)度進(jìn)行對(duì)比，進(jìn)而計(jì)算樣品硫含量。

油品中氮含量的測(cè)定采用化學(xué)發(fā)光法，分析儀器為硫氮分析儀(型號(hào)為MultiTek，德國(guó)ANTEK公司)。分析過程中，將樣品在室溫下放入樣品舟中，由進(jìn)樣器將樣品舟送至高溫燃燒管，在O2氣氛下，N被氧化成NO，然后NO與O3接觸轉(zhuǎn)化為激發(fā)態(tài)的NO2，激發(fā)態(tài)NO2衰減時(shí)的發(fā)射光由光電倍增管檢測(cè)，從而計(jì)算樣品中的氮含量。

圖1 滴流床加氫反應(yīng)器裝置Fig.1 Trickle bed reactor device for testing hydrogenation of FCC recycle oil1—Feed tank; 2—Feed pump; 3—H2 mass flow meter; 4—Buffer tank; 5—Fixed-bed reactor;6—High-pressure separator; 7—Preheater; 8—Stabilizer; 9—Product oil tank; 10—Water tank

2 催化裂化回?zé)捰图託浞磻?yīng)動(dòng)力學(xué)模型

2.1 油品組分的定義及計(jì)算

根據(jù)石油碳?xì)浠衔镏谢鶊F(tuán)飽和度的差異，可以將石油及其產(chǎn)品中的碳?xì)浠衔锝M成進(jìn)行重新定義。在本研究中，石油及其產(chǎn)品中的碳?xì)浠衔锉怀橄蠖x為由“虛擬組分”CH基團(tuán)和CH2基團(tuán)兩部分組成。其中，CH基團(tuán)代表不飽和烴，其碳/氫原子比為1，用符號(hào)CH表示；CH2基團(tuán)代表飽和烴，其碳/氫原子比為1/2，用符號(hào)CH2表示。根據(jù)此定義，油品中碳?xì)浠衔镏蠧H基團(tuán)和CH2基團(tuán)兩組分的質(zhì)量分?jǐn)?shù)可由其C、H元素的質(zhì)量分?jǐn)?shù)計(jì)算得到，計(jì)算式分別為：

(1)

wCH2=(1-wCH)×100%

(2)

2.2 動(dòng)力學(xué)模型的構(gòu)建

根據(jù)石油原料及其產(chǎn)品中碳?xì)浠衔锝M成定義，不飽和烴加氫反應(yīng)即為CH基團(tuán)加氫生成CH2基團(tuán)的反應(yīng)。因此，在催化裂化回?zé)捰图託溥^程中主要發(fā)生3類反應(yīng)：CH與CH2間的加氫-脫氫可逆反應(yīng)、HDS反應(yīng)和HDN反應(yīng)。各反應(yīng)的化學(xué)反應(yīng)方程如下：

(3)

(4)

(5)

2.2.1 模型方程

文獻(xiàn)[12-16]研究表明，一級(jí)反應(yīng)動(dòng)力學(xué)模型可以較好地描述芳烴的加氫-脫氫可逆反應(yīng)、HDS反應(yīng)和HDN反應(yīng)。催化裂化回?zé)捰椭蟹紵N約占78%，因而回?zé)捰椭蠧H基團(tuán)以芳烴為主。因此，采用一級(jí)反應(yīng)動(dòng)力學(xué)模型適用于描述催化裂化回?zé)捰偷募託溥^程，并可以較好的模擬預(yù)測(cè)加氫過程中C、H、S、N元素的遷移規(guī)律。另外，由于吸附在催化劑活性中心上的氣態(tài)氫在加氫過程中可以轉(zhuǎn)化為活性氫[17]，參與加氫反應(yīng)，因而氫分壓對(duì)該加氫反應(yīng)有一定影響[18]：氫分壓高可以促進(jìn)氣態(tài)氫轉(zhuǎn)化為活性氫，進(jìn)而提高催化劑表面活性氫的含量[19-20]。因此，在各加氫反應(yīng)模型方程中引入了氫分壓項(xiàng)來描述氫分壓對(duì)催化裂化回?zé)捰图託浞磻?yīng)的影響。

根據(jù)各反應(yīng)化學(xué)反應(yīng)方程式，可以得出CH與CH2組分間的加氫-脫氫可逆反應(yīng)、HDS反應(yīng)和HDN反應(yīng)的反應(yīng)動(dòng)力學(xué)模型微分方程，如式(6)～(8)所示。

(6)

(7)

(8)

對(duì)模型微分方程進(jìn)行積分，可以得到模型積分形式，如式(9)～(11)所示。

(9)

wS=wS,0×exp(-k3pα3t)

(10)

wN=wN,0×exp(-k4pα4t)

(11)

其中，ki為反應(yīng)速率常數(shù)，可以用阿倫尼烏斯方程表示：

(12)

2.2.2 模型參數(shù)求解流程

根據(jù)模型方程積分式和試驗(yàn)數(shù)據(jù)，圖2列出了求解模型參數(shù)的計(jì)算流程。由圖2可知，首先將原料油中CH組分和S、N元素含量數(shù)值，反應(yīng)操作條件(不同氫分壓下系列條件)數(shù)值分別代入各模型方程，得到模型數(shù)據(jù)計(jì)算值，并將其與試驗(yàn)值進(jìn)行比較。以模型數(shù)據(jù)值與試驗(yàn)值間的總方差∑R2為目標(biāo)值，求得氫分壓指數(shù)αi值，并將其代入模型方程。然后，利用相同方法可得不同反應(yīng)溫度下的模型數(shù)據(jù)值。根據(jù)阿倫尼烏斯方程，可得各反應(yīng)指前

圖2 模型參數(shù)計(jì)算流程Fig.2 Calculation process chart of model parameter

因子A0,i和反應(yīng)活化能Ei值，最后整理得到模型方程參數(shù)αi、A0,i和Ei值。

2.3 模型計(jì)算與討論

2.3.1 CH基團(tuán)隨工藝條件的變化規(guī)律

按工藝試驗(yàn)數(shù)據(jù)，CH基團(tuán)隨工藝條件的變化規(guī)律如圖3所示。其中，圖3(a)為氫分壓的變化對(duì)CH組分質(zhì)量分?jǐn)?shù)的影響；圖3(b)～(e)分別為反應(yīng)溫度為633.15 K、643.15 K、653.15 K和658.15 K時(shí)CH組分質(zhì)量分?jǐn)?shù)隨停留時(shí)間t的變化。由圖3(a)可以看出，隨著氫分壓增加，CH組分質(zhì)量分?jǐn)?shù)下降，說明氫分壓升高可以促進(jìn)氣態(tài)氫轉(zhuǎn)化為活性氫，使催化劑表面活性氫濃度提高，進(jìn)而促進(jìn)CH組分加氫轉(zhuǎn)化為CH2組分。由圖3(b)～(e)可以發(fā)現(xiàn)，當(dāng)反應(yīng)溫度在633.15～658.15 K時(shí)，CH組分質(zhì)量分?jǐn)?shù)隨停留時(shí)間t增加而降低，說明無論從反應(yīng)熱力學(xué)還是反應(yīng)動(dòng)力學(xué)講，延長(zhǎng)停留時(shí)間均有利于CH組分加氫轉(zhuǎn)化為CH2組分。同時(shí)，由圖2模型參數(shù)計(jì)算流程得到CH與CH2間加氫-脫氫可逆反應(yīng)的氫分壓指數(shù)α1=0.50，而該反應(yīng)化學(xué)方程式中H2計(jì)量系數(shù)也為0.5，說明CH與CH2間加氫-脫氫可逆反應(yīng)過程更能反映和接近不飽和烴加氫的本征反應(yīng)過程。反應(yīng)k2/k1值隨反應(yīng)溫度升高而增大，說明CH與CH2間加氫-脫氫可逆反應(yīng)的逆反應(yīng)受反應(yīng)溫度影響比正反應(yīng)大，即隨溫度升高該反應(yīng)受熱力學(xué)影響限制也增大。

2.3.2 S含量隨工藝條件的變化規(guī)律

圖4為S質(zhì)量分?jǐn)?shù)隨工藝條件的變化規(guī)律。圖4(a)為S質(zhì)量分?jǐn)?shù)隨氫分壓的變化；圖4(b)～(e)分別為反應(yīng)溫度633.15 K、643.15 K、653.15 K和658.15 K 時(shí)S質(zhì)量分?jǐn)?shù)隨停留時(shí)間t的變化。由圖4(a)可以看出：油品中S質(zhì)量分?jǐn)?shù)隨氫分壓增加出現(xiàn)輕微下降，說明在HDS反應(yīng)過程中，通過加氫飽和-氫解脫硫反應(yīng)脫硫比例很低，而通過直接氫解脫硫反應(yīng)脫硫比例很高，因此氫分壓對(duì)HDS反應(yīng)的的影響不明顯；同時(shí)，計(jì)算得到α3=0.02，其數(shù)值很小也說明氫分壓對(duì)HDS反應(yīng)的影響很弱。由圖4(b)～(e) 可以發(fā)現(xiàn)：當(dāng)反應(yīng)溫度在 633.15～658.15 K 時(shí)，S質(zhì)量分?jǐn)?shù)隨停留時(shí)間t增加而下降，說明無論從反應(yīng)熱力學(xué)還是反應(yīng)動(dòng)力學(xué)講，延長(zhǎng)停留時(shí)間均利于HDS反應(yīng)進(jìn)行；同時(shí)，k3隨反應(yīng)溫度升高而增大，說明升高反應(yīng)溫度，提高了HDS反應(yīng)速率，促進(jìn)了HDS反應(yīng)進(jìn)行。

圖3 CH組分質(zhì)量分?jǐn)?shù)隨工藝條件的變化規(guī)律Fig.3 The change of CH mass fraction with operation conditions(a) Reaction conditions： T=633.15 K， LHSV=1.0 h-1， V(H2)/V(Oil)=800;(b) T=633.15 K; (c) T=643.15 K; (d) T=653.15 K; (e) T=658.15 K；(b)-(e) Reaction conditions： LHSV is 0.6-1.6 h-1， V(H2)/V(Oil)=800， p=10 MPa

2.3.3 N含量隨工藝條件的變化規(guī)律

圖5為N質(zhì)量分?jǐn)?shù)隨工藝條件的變化規(guī)律。其中，圖5(a)為N質(zhì)量分?jǐn)?shù)隨氫分壓的變化；圖5(b)～(e)為當(dāng)反應(yīng)溫度在633.15 K、643.15 K、653.15 K和658.15 K時(shí)N質(zhì)量分?jǐn)?shù)隨停留時(shí)間t的變化。由圖5(a)可以看出，油品中N質(zhì)量分?jǐn)?shù)隨著氫分壓升高而下降，說明氮化物加氫飽和過程受氫分壓影響較大，提高氫分壓可以促進(jìn)HDN反應(yīng)。這也說明回?zé)捰虷DN反應(yīng)途徑主要為芳環(huán)加氫-氫解脫氮，即氮化物的芳環(huán)結(jié)構(gòu)部分優(yōu)先發(fā)生加氫飽和反應(yīng)，然后發(fā)生氫解脫氮。同時(shí)，計(jì)算得到α4=1.16，其數(shù)值較大也反映了氫分壓對(duì)HDN反應(yīng)的影響較大。由圖5(b)～(e)可以發(fā)現(xiàn)：當(dāng)反應(yīng)溫度為633.15～658.15 K時(shí)，油品中N質(zhì)量分?jǐn)?shù)隨停留時(shí)間t增加而下降，說明無論從反應(yīng)熱力學(xué)還是反應(yīng)動(dòng)力學(xué)講，延長(zhǎng)停留時(shí)間均利于HDN反應(yīng)進(jìn)行；同時(shí)，k4值隨反應(yīng)溫度增加而增大，說明提高反應(yīng)溫度，提高了HDN反應(yīng)速率，促進(jìn)了HDN反應(yīng)進(jìn)行。

圖4 S質(zhì)量分?jǐn)?shù)隨工藝條件的變化規(guī)律Fig.4 The change of S mass fraction with operation conditions(a) Reaction conditions： T=633.15 K， LHSV=1.0 h-1， V(H2)/V(Oil)=800;(b) T=633.15 K; (c) T=643.15 K; (d) T=653.15 K; (e) T=658.15 K；(b)-(e) Reaction conditions： LHSV is 0.6-1.6 h-1， V(H2)/V(Oil)=800， p=10 MPa

利用阿倫尼烏斯方程可以求得指前因子A0和活化能E值。圖6為lnk隨1/T變化曲線，根據(jù)圖中線性擬合方程參數(shù)，可得指前因子A0,i和活化能Ei值，其中i為1～4。

根據(jù)圖2所示模型計(jì)算流程和圖3～圖6所示模擬計(jì)算數(shù)據(jù)，最終整理可得模型參數(shù)αi、A0,i和Ei的值，見表2。表2列出了催化裂化回?zé)捰图託浞磻?yīng)動(dòng)力學(xué)模型參數(shù)，根據(jù)表2模型參數(shù)數(shù)據(jù)可以進(jìn)行模型驗(yàn)證和預(yù)測(cè)。

2.4 動(dòng)力學(xué)模型驗(yàn)證

為了驗(yàn)證所建模型可靠性和增強(qiáng)模型外推性，在以上操作范圍外另取2組數(shù)據(jù)。此2組數(shù)據(jù)的操作條件分別為：反應(yīng)溫度663.15 K，氫分壓10 MPa，氫/油體積比800，液時(shí)空速4.0 h-1；反應(yīng)溫度663.15 K，氫分壓10 MPa，氫/油體積比800，液時(shí)空速5.0 h-1。然后根據(jù)所建模型可以計(jì)算得到此2組操作條件下加氫油中C、H、S、N和CH組分含量的預(yù)測(cè)值。表3為加氫油中C、H、S、

圖5 N質(zhì)量分?jǐn)?shù)隨工藝條件的變化規(guī)律Fig.5 The change of N mass fraction with operation conditions(a) Reaction conditions： T=633.15 K， LHSV=1.0 h-1， V(H2)/V(Oil)=800;(b) T=633.15 K; (c) T=643.15 K; (d) T=653.15 K; (e) T=658.15 K；(b)-(e) Reaction conditions： LHSV is 0.6-1.6 h-1， V(H2)/V(Oil)=800， p=10 MPa

圖6 ln k隨1/T變化曲線Fig.6 The curve of ln k to 1/T

表2 催化裂化回?zé)捰图託浞磻?yīng)動(dòng)力學(xué)參數(shù)Table 2 Kinetic parameters of hydrogenationreaction of FCC recycle oil

Here the symbol ‘i’ refers specially to a calculated integer.

N和CH組分試驗(yàn)值與預(yù)測(cè)值間的對(duì)比。從表3可以看出，加氫油中C、H、S、N和CH組分含量的試驗(yàn)值與預(yù)測(cè)值比較吻合，說明所構(gòu)建的催化裂化回?zé)捰图託浞磻?yīng)動(dòng)力學(xué)模型可以較好地預(yù)測(cè)加氫油中C、H、S、N和CH組分含量。

表3 催化裂化回?zé)捰图託溆椭蠧、H、S、N和CH組分試驗(yàn)值與預(yù)測(cè)值的對(duì)比Table 3 Comparison between measured and predictedC, H, S, N and CH content in hydrogenation FCC recycle oil

1#—The test under the conditions ofT=633.15 K， LHSV=4.0 h-1，V(H2)/V(Oil)=800,p=10 MPa;

2#—The test under the conditions ofT=633.15 K， LHSV=5.0 h-1，V(H2)/V(Oil)=800,p=10 MPa

3 結(jié) 論

(1)構(gòu)建了描述催化裂化回?zé)捰图託溥^程模型，將加氫反應(yīng)過程簡(jiǎn)化為3類反應(yīng)：CH組分與CH2組分間的加氫-脫氫可逆反應(yīng)、HDS和HDN反應(yīng)。各反應(yīng)反應(yīng)級(jí)數(shù)設(shè)定為1級(jí)，可較好的反映催化裂化回?zé)捰偷募託浞磻?yīng)規(guī)律，并簡(jiǎn)化了模型構(gòu)建流程。

(2)“虛擬組分”CH基團(tuán)和CH2基團(tuán)間的加氫-脫氫可逆反應(yīng)動(dòng)力學(xué)模型很好地描述了回?zé)捰图託溥^程中C、H變化規(guī)律；HDS和HDN反應(yīng)動(dòng)力學(xué)模型也較好的描述了回?zé)捰图託溥^程中S、N含量變化。通過構(gòu)建的動(dòng)力學(xué)模型可以預(yù)測(cè)回?zé)捰图託溥^程中的C、H、S、N含量。

符號(hào)說明：

A0,i(i為1～4)——各反應(yīng)反應(yīng)指前因子；

Ei——各反應(yīng)表觀反應(yīng)活化能，kJ/mol；

ki(i為1～4)——各反應(yīng)反應(yīng)表觀速率常數(shù)，其中ki=Aie-Ei/RT；

LHSV——液時(shí)體積空速，h-1；

p——?dú)浞謮?，MPa；

R——阿伏伽德羅常數(shù)，R=8.314 J/(mol·K)；

R2——模型值與試驗(yàn)值間的方差；

∑R2——模型值與試驗(yàn)值間的總方差；

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

t——停留時(shí)間(t=1/LHSV)，h；

wC——原料油或加氫油中C元素質(zhì)量分?jǐn)?shù)，%；

wH——原料油或加氫油中H元素質(zhì)量分?jǐn)?shù)，%；

wCH——加氫油中CH組分質(zhì)量分?jǐn)?shù)，%；

wCH,0——原料油中CH組分質(zhì)量分?jǐn)?shù)，%；

wCH2——加氫油中CH2組分質(zhì)量分?jǐn)?shù)，%；

wS,0——原料油中S元素質(zhì)量分?jǐn)?shù)，%；

wS——加氫油中S元素質(zhì)量分?jǐn)?shù)，%；

wN,0——原料油中N元素質(zhì)量分?jǐn)?shù)，%；

wN——加氫油中N元素質(zhì)量分?jǐn)?shù)，%；

αi(i=1,3,4)——各反應(yīng)中氫分壓指數(shù)項(xiàng)。