正戊烷異構(gòu)化制異戊烷反應(yīng)熱力學(xué)分析

劉洪全, 于中偉, 馬愛(ài)增, 張秋平， 孫義蘭, 王子健

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

C5/C6異構(gòu)化是將C5和C6正構(gòu)烷烴轉(zhuǎn)化為相應(yīng)的異構(gòu)烷烴的過(guò)程。異構(gòu)化油具有硫含量低、不含烯烴和芳烴的特點(diǎn)，是理想的高辛烷值汽油調(diào)和組分，可提高發(fā)動(dòng)機(jī)的前端辛烷值和啟動(dòng)性能[1-3]。近年來(lái)，隨著我國(guó)汽油升級(jí)步伐的加快，C5/C6異構(gòu)化作為油品升級(jí)不可或缺的技術(shù)之一，越來(lái)越引起人們的重視。目前，已工業(yè)應(yīng)用的C5/C6異構(gòu)化技術(shù)，按催化劑體系可分為3類[4]，分別是中溫分子篩異構(gòu)化、固體超強(qiáng)酸異構(gòu)化和低溫異構(gòu)化。中溫分子篩異構(gòu)化技術(shù)通常以負(fù)載Pt或Pd的沸石分子篩為催化劑，反應(yīng)溫度為513.15～553.15 K、反應(yīng)壓力為1.6～3.0 MPa、氫/烴摩爾比為1～2。固體超強(qiáng)酸異構(gòu)化技術(shù)以負(fù)載Pt的硫酸化氧化鋯為催化劑，反應(yīng)溫度為413.15～473.15 K、反應(yīng)壓力為1.6～3.5 MPa、氫/烴摩爾比為1～2。低溫異構(gòu)化技術(shù)通常以負(fù)載Pt的氯化氧化鋁為催化劑，反應(yīng)溫度為393.15～443.15 K、反應(yīng)壓力為2.5～3.5 MPa、氫/烴摩爾比為0.05～0.5。

目前有關(guān)C5/C6異構(gòu)化反應(yīng)熱力學(xué)研究報(bào)道較少，僅有的幾篇文獻(xiàn)基本采用了假設(shè)反應(yīng)物為氣相，且為理想氣體的方式進(jìn)行計(jì)算和分析[5-6]。而從現(xiàn)有 C5/C6異構(gòu)化技術(shù)來(lái)看，反應(yīng)溫度范圍為393.15～553.15 K、壓力為1.6～3.5 MPa、氫/烴摩爾比為0.05～2，反應(yīng)溫度較低而壓力較高，在這樣的條件下，如果不考慮反應(yīng)物和產(chǎn)物的非理想性，可能會(huì)產(chǎn)生較大偏差，無(wú)法很好指導(dǎo)實(shí)驗(yàn)研究。筆者以正戊烷為研究對(duì)象，結(jié)合實(shí)際反應(yīng)體系，進(jìn)行正戊烷異構(gòu)化反應(yīng)熱力學(xué)分析，為C5/C6異構(gòu)化反應(yīng)工藝條件優(yōu)化提供參考。

1 熱力學(xué)分析方法

對(duì)于正戊烷異構(gòu)化反應(yīng)，主反應(yīng)為正戊烷異構(gòu)化生成異戊烷。采用現(xiàn)有的3種異構(gòu)化技術(shù)時(shí)，該反應(yīng)選擇性都很高，因此在熱力學(xué)分析中僅考察正戊烷異構(gòu)化生成異戊烷的反應(yīng)，不考慮其他副反應(yīng)。

1.1 反應(yīng)的和值的計(jì)算

假設(shè)反應(yīng)為氣相反應(yīng)，且正戊烷和異戊烷均為理想氣體，根據(jù)Kirchhoff公式，反應(yīng)體系的摩爾焓變?chǔ)Hm(T)的表達(dá)式如式(1)所示。

(1)

式(1)中，Cp,i為參加反應(yīng)各物質(zhì)的等壓摩爾熱容，其表達(dá)式如式(2)所示。

Cp,i=ai+biT+ciT2

(2)

式(2)中，ai、bi和ci是經(jīng)驗(yàn)常數(shù)，由各組分本身的特性決定。

通過(guò)查閱物理化學(xué)數(shù)據(jù)手冊(cè)[7-8]，可以得到正戊烷異構(gòu)化反應(yīng)中各組分的典型熱力學(xué)數(shù)據(jù)，將不同溫度下各組分的等壓摩爾熱容值代入式(2)，進(jìn)行線性回歸，可以得到常數(shù)ai、bi和ci的值，結(jié)果列于表1。

表1 由式(2)計(jì)算得到的正戊烷和異戊烷的等壓摩爾熱容的經(jīng)驗(yàn)常數(shù)Table 1 The empirical constants of constant pressure molar heat capacity of n-pentane and isopentane calculated by Eq.(2)

將式(2)代入式(1)，并進(jìn)行積分，可得式(3)。

(3)

Δa=∑νiai
Δb=∑νibi
Δc=∑νici

式(3)中，I1為積分常數(shù)，可將標(biāo)準(zhǔn)摩爾生成焓數(shù)據(jù)代入求得。將Gibss-Helmholtz公式(式(4))代入式(3),并進(jìn)行積分,可得式(5)。

(4)

(5)

式(5)中，I2為積分常數(shù)，可將標(biāo)準(zhǔn)摩爾生成吉布斯自由能數(shù)據(jù)代入求得。將理想氣體體系標(biāo)準(zhǔn)平衡常數(shù)的定義式(式(6))代入式(5),可得正戊烷異構(gòu)化反應(yīng)標(biāo)準(zhǔn)平衡常數(shù)與溫度的關(guān)系式(7)。

(6)

(7)

1.2 反應(yīng)體系的化學(xué)平衡和相平衡

(8)

(9)

(10)

由于正戊烷異構(gòu)化反應(yīng)壓力較高而反應(yīng)溫度較低，正戊烷和產(chǎn)物異戊烷有可能以氣-液兩相存在。在這種情況下，反應(yīng)過(guò)程不僅存在化學(xué)平衡，還存在相平衡。組分i相平衡的條件是其在各相中的化學(xué)勢(shì)相等，對(duì)于氣-液兩相體系，可用式(11)表示。將混合物中組分i逸度的定義式(12)[9]代入式(11)，得到另一個(gè)等效的以逸度表示的相平衡判據(jù)，如式(13)所示。將式(13)用逸度系數(shù)表示，則變?yōu)槭?14)。

μiV=μiL

(11)

(12)

(13)

(14)

汽化平衡比Ki的定義式如式(15)所示。從式(10)和式(15)可以看出，反應(yīng)平衡常數(shù)Kf和汽化平衡比Ki由各組元的平衡濃度和逸度系數(shù)決定。

(15)

1.3 組元逸度系數(shù)的計(jì)算方法

逸度系數(shù)可利用狀態(tài)方程來(lái)計(jì)算，對(duì)于烴類體系，選用Soave-Redlich-Kwong(SRK)狀態(tài)方程較為合適[10-11]。SRK狀態(tài)方程的標(biāo)準(zhǔn)形式見(jiàn)式(16)，其多項(xiàng)式形式見(jiàn)式(17)。

(16)

z3-z2+(A-B-B2)-AB=0

(17)

式(16)和(17)中，參數(shù)λ、γ、A和B的表達(dá)式分別如式(18)～(21)所示。

(18)

γ=0.08664RTc/pc

(19)

A=λαp/R2T2

(20)

B=γp/RT=0.08664pr/Tr

(21)

SRK狀態(tài)方程式(16)中的α稱為Alpha方程，Alpha方程為溫度的函數(shù)，在應(yīng)用時(shí)可根據(jù)物質(zhì)的種類選擇不同的Alpha方程用于計(jì)算。對(duì)于烴類，筆者采用Soave[12]給出的Alpha方程，如式(22)所示；對(duì)于氫氣，采用Twu等[13]提出的Alpha方程，如式(23)所示。

(22)

(23)

對(duì)于混合物，SRK狀態(tài)方程式(16)中λα和γ的表達(dá)式分別如式(24)和(25)所示。

λα=∑∑yiyj(λα)ij

(24)

γ=∑yiγi

(25)

式(24)中的(λα)ij可由式(26)表示，式(26)中的kij是與濃度無(wú)關(guān)的交互作用參數(shù)，對(duì)于烴類和氫氣，取kij值為0。

(26)

用SRK狀態(tài)方程求取組元逸度系數(shù)的通式見(jiàn)式(27)[14]。

(27)

式(27)中，Bi可表示為式(28)。

Bi=γip/RT

(28)

從上述公式可知，應(yīng)用SRK狀態(tài)方程求取逸度系數(shù)，需要各組元的物性參數(shù)，包括臨界溫度Tc、臨界壓力pc和偏心因子ω。本文中正戊烷和異戊烷的物性參數(shù)引自文獻(xiàn)[15]，氫氣的物性參數(shù)引自文獻(xiàn)[8]。

1.4 反應(yīng)平衡組成的計(jì)算

如前所述，正戊烷異構(gòu)化反應(yīng)可能是氣相、液相或氣-液兩相，要計(jì)算反應(yīng)平衡組成，首先要判斷反應(yīng)條件下的相態(tài)，而要判斷反應(yīng)相態(tài)，則要計(jì)算反應(yīng)體系的露點(diǎn)和泡點(diǎn)。

在特定的反應(yīng)條件下，由反應(yīng)進(jìn)度可得到如式(29)所示的各組元的物質(zhì)的量。

ni=ni0+ενi

(29)

式(29)中，νi對(duì)反應(yīng)物為負(fù)值，對(duì)產(chǎn)物為正值，由于表觀上氫氣不參與異構(gòu)化反應(yīng)，其賦值為0。

如果反應(yīng)體系為氣-液兩相，則各組元的總物質(zhì)的量為氣相物質(zhì)的量和液相物質(zhì)的量之和，由式(30)表示。

(30)

各組元的物料衡算如式(31)所示。

ni=nLxi+nVyi=[nL+Ki(nt-nL)]xi

(31)

由于摩爾分?jǐn)?shù)的總和應(yīng)等于1，由式(31)可得到式(32)。

(32)

分別將式(29)和式(30)代入式(32)，可得式(33)。

(33)

如果指定反應(yīng)溫度和起始組成，采用上述步驟可計(jì)算露點(diǎn)壓力。以nL值逼近nt為條件，采用上述步驟也可計(jì)算反應(yīng)體系的泡點(diǎn)。

如果反應(yīng)處于露點(diǎn)和泡點(diǎn)之間的兩相區(qū)，則可采用上述計(jì)算露點(diǎn)溫度方法中的(1)～(4)步來(lái)計(jì)算出平衡組成。

2 結(jié)果與討論

2.1 正戊烷異構(gòu)化反應(yīng)的熱效應(yīng)

表2 不同溫度下正戊烷異構(gòu)化反應(yīng)的標(biāo)準(zhǔn)摩爾焓變標(biāo)準(zhǔn)摩爾吉布斯自由能變和標(biāo)準(zhǔn)平衡常數(shù)值Table of n-pentane isomerization at different temperatures

2.2 反應(yīng)體系相態(tài)的判斷

按照1.4節(jié)的方法計(jì)算了各典型反應(yīng)壓力和氫/烴摩爾比條件下反應(yīng)體系的露點(diǎn)溫度，結(jié)果如表3所示。由于氫氣的存在，反應(yīng)體系的泡點(diǎn)溫度不易精確求出，但通過(guò)1.4節(jié)的計(jì)算可以判斷表3所列壓力和氫/烴摩爾比條件下的泡點(diǎn)溫度均低于273.15 K。

表3 不同壓力和氫/烴摩爾比(n(H2)/n(HC))下正戊烷異構(gòu)化反應(yīng)的露點(diǎn)溫度(Td)Table 3 Dew temperatures (Td) of n-pentane isomerization at different pressures and n(H2)/n(HC)

1), 2), 3), 4)—n(H2)/n(HC) value

從表3可知，反應(yīng)壓力和氫/烴摩爾比均對(duì)露點(diǎn)溫度有顯著的影響，相同壓力下氫/烴摩爾比越高露點(diǎn)溫度越低，而相同氫/烴摩爾比下壓力越高露點(diǎn)溫度越高。

利用表3數(shù)據(jù)還可以判斷現(xiàn)有異構(gòu)化技術(shù)的相態(tài)。在反應(yīng)壓力為1.6～3.5 MPa、氫/烴摩爾比為1～2的條件下，體系的露點(diǎn)溫度最高達(dá)408.05 K，低于現(xiàn)有中溫異構(gòu)化和固體超強(qiáng)酸異構(gòu)化的反應(yīng)溫度，因此可判斷在中溫異構(gòu)化和固體超強(qiáng)酸異構(gòu)化的反應(yīng)條件下，正戊烷異構(gòu)化為氣相反應(yīng)。從表3還可以看到，在反應(yīng)壓力為2.5～3.5 MPa、氫/烴摩爾比為0.05～0.5條件下，體系的露點(diǎn)溫度范圍為388.45～456.95 K，而低溫異構(gòu)化的反應(yīng)溫度為393.15～443.15 K，二者的溫度范圍有重疊，因此可以判斷低溫異構(gòu)化反應(yīng)體系可能為氣相或氣-液兩相。

2.3 反應(yīng)溫度對(duì)正戊烷異構(gòu)化反應(yīng)的影響

在反應(yīng)壓力為3.0 MPa、氫/烴摩爾比為0.5的條件下，計(jì)算了反應(yīng)溫度對(duì)正戊烷異構(gòu)化反應(yīng)的影響。

首先計(jì)算了不同溫度下反應(yīng)的相態(tài)、各相占比及組成，結(jié)果如表4所示。從表4可以看到：反應(yīng)溫度為323.15～413.15 K時(shí)，反應(yīng)體系呈氣-液兩相；隨著溫度升高，更多的液相組分被汽化，氣-液兩相的摩爾比逐漸增加。此外，在該溫度范圍內(nèi)，氣相中氫氣的含量隨著反應(yīng)溫度的升高而下降，而液相中的氫含量則先升高后下降，在393.15 K左右達(dá)到最高。

表4 不同反應(yīng)溫度下正戊烷異構(gòu)化反應(yīng)的相態(tài)、各相以及每相中各組分的平衡摩爾分?jǐn)?shù)Table 4 Phase states, equilibrium molar fractions of different phases and components in each phase for n-pentane isomerization at different reaction temperatures

p=3.0 MPa;n(H2)/n(HC)=0.5

圖1給出了不同溫度下正戊烷異構(gòu)化反應(yīng)的平衡異構(gòu)化率，曲線(1)是假設(shè)反應(yīng)體系為理想氣體的計(jì)算結(jié)果，曲線(2)是考慮體系的非理想性、用SRK方程進(jìn)行計(jì)算的結(jié)果。從圖1可以看出，無(wú)論采用哪種假設(shè)，正戊烷的平衡異構(gòu)化率都隨著反應(yīng)溫度的升高而降低，說(shuō)明正戊烷異構(gòu)化反應(yīng)在熱力學(xué)上是低溫有利的。反應(yīng)溫度高于423.15 K時(shí)，兩條曲線基本重合，說(shuō)明高于該溫度時(shí)反應(yīng)體系的非理想性并不明顯，可參照理想氣體來(lái)計(jì)算。而當(dāng)反應(yīng)溫度低于423.15 K時(shí)，采用SRK方程計(jì)算出的平衡異構(gòu)化率低于理想狀態(tài)的計(jì)算結(jié)果，且溫度越低，差異越大。在前述3種異構(gòu)化技術(shù)中，低溫異構(gòu)化的反應(yīng)溫度通常低于423.15 K，因此其反應(yīng)體系的非理想性不可忽略。

圖1 正戊烷異構(gòu)化反應(yīng)中平衡異構(gòu)化率隨反應(yīng)溫度(T)的變化Fig.1 Equilibrium isomerization ratio vs Tfor n-pentane isomerization(1) The reaction system is considered as ideal gas; (2) The nonidealityof the reaction system is considered using SRK equationp=3.0 MPa; n(H2)/n(HC)=0.5

2.4 反應(yīng)壓力對(duì)正戊烷異構(gòu)化反應(yīng)的影響

在反應(yīng)溫度為403.15 K、氫/烴摩爾比為0.5的條件下，計(jì)算了反應(yīng)壓力對(duì)正戊烷異構(gòu)化反應(yīng)的影響。

首先計(jì)算了不同溫度下反應(yīng)的相態(tài)、各相占比及組成，結(jié)果如表5所示。從表5可以看出，在低壓下，反應(yīng)體系呈氣相，隨著壓力升高，反應(yīng)體系由氣相轉(zhuǎn)變?yōu)闅?液兩相，且液相的比例隨著壓力的進(jìn)一步升高而增加，液相中氫的濃度也隨著壓力的升高而略有增加。

表5 不同反應(yīng)壓力下正戊烷異構(gòu)化反應(yīng)的相態(tài)、各相以及每相中各組分的平衡摩爾分?jǐn)?shù)Table 5 Phase states, equilibrium molar fractions of different phases and components in each phase for n-pentane isomerization at different reaction pressures

T=403.15 K;n(H2)/n(HC)=0.5

圖2是反應(yīng)溫度403.15 K、氫/烴摩爾比0.5條件下,正戊烷平衡異構(gòu)化率與反應(yīng)壓力的變化關(guān)系曲線。圖中曲線(1)是假設(shè)反應(yīng)體系為理想氣體的計(jì)算結(jié)果，曲線(2)是基于SRK方程的計(jì)算結(jié)果。從圖2可以看出，理想狀態(tài)下反應(yīng)壓力對(duì)正戊烷的平衡異構(gòu)化率沒(méi)有影響。當(dāng)采用SRK方程計(jì)算時(shí)，壓力低于1.6 MPa時(shí)的計(jì)算結(jié)果與理想氣體的結(jié)果基本一致，說(shuō)明反應(yīng)溫度為403.15 K、反應(yīng)壓力低于1.6 MPa時(shí)，反應(yīng)體系的非理想性也可以忽略。而當(dāng)反應(yīng)壓力高于1.6 MPa時(shí)，采用SRK方程計(jì)算出的正戊烷平衡異構(gòu)化率隨著反應(yīng)壓力的升高而下降，說(shuō)明高壓下反應(yīng)體系的非理想性逐漸顯現(xiàn)，且提高壓力對(duì)反應(yīng)平衡不利。

圖2 正戊烷異構(gòu)化反應(yīng)中平衡異構(gòu)化率隨反應(yīng)壓力(p)的變化Fig.2 Equilibrium isomerization ratio vs pfor n-pentane isomerization(1) The reaction system is considered as ideal gas; (2) The nonidealityof the reaction system is considered using SRK equationT=403.15 K; n(H2)/n(HC)=0.5

2.5 氫/烴摩爾比對(duì)正戊烷異構(gòu)化反應(yīng)的影響

在反應(yīng)溫度為403.15 K、反應(yīng)壓力為3.0 MPa的條件下，計(jì)算了氫/烴摩爾比對(duì)正戊烷異構(gòu)化反應(yīng)的影響。

首先計(jì)算了不同溫度下反應(yīng)的相態(tài)、各相占比及組成，結(jié)果如表6所示。從表6可以看出：在氫/烴摩爾比高于1.0時(shí)，反應(yīng)體系為氣相；隨著氫/烴摩爾比下降，反應(yīng)體系由氣相轉(zhuǎn)變?yōu)闅?液兩相，當(dāng)氫/烴摩爾比為0.5時(shí)，液相摩爾分?jǐn)?shù)為32.75%；進(jìn)一步降低氫/烴摩爾比，液相比例進(jìn)一步增加，當(dāng)氫/烴摩爾比為0.05時(shí)，液相摩爾分?jǐn)?shù)高達(dá)96.14%。此外，值得注意的是，在氫/烴摩爾比為0.05～0.5時(shí)，氣-液兩相中各組分的平衡摩爾分?jǐn)?shù)基本保持不變。

圖3給出了不同氫/烴摩爾比條件下正戊烷異構(gòu)化反應(yīng)的平衡異構(gòu)化率。曲線(1)是假設(shè)反應(yīng)體系為理想氣體的計(jì)算結(jié)果，曲線(2)是基于SRK方程的計(jì)算結(jié)果。從圖3可以看出，理想狀態(tài)下氫/烴摩爾比對(duì)正戊烷的平衡異構(gòu)化率沒(méi)有影響。采用SRK方程計(jì)算時(shí)，當(dāng)氫/烴摩爾比高于1.0時(shí)，其計(jì)算結(jié)果與理想狀態(tài)的基本一致；而當(dāng)氫/烴摩爾比低于1.0時(shí)，正戊烷平衡異構(gòu)化率隨氫/烴摩爾比的降低而下降，說(shuō)明低氫/烴摩爾比時(shí)反應(yīng)體系的非理想性比較明顯，且降低氫/烴摩爾比對(duì)反應(yīng)平衡不利。

表6 不同氫/烴摩爾比(n(H2)/n(HC))下正戊烷異構(gòu)化反應(yīng)的相態(tài)、各相以及每相中各組分的平衡摩爾分?jǐn)?shù)Table 6 Phase states, equilibrium molar fractions of different phases and components in each phase for n-pentane isomerization with different n(H2)/n(HC)

T=403.15 K;p=3.0 MPa

圖3 正戊烷異構(gòu)化反應(yīng)中平衡異構(gòu)化率隨氫/烴摩爾比(n(H2)/n(HC))的變化Fig.3 Equilibrium isomerization ratio vs n(H2)/n(HC)for n-pentane isomerization(1) The reaction system is considered as ideal gas; (2) The nonidealityof the reaction system is considered using SRK equationT=403.15 K; p=3.0 MPa

3 結(jié) 論

(1)正戊烷異構(gòu)化屬微放熱反應(yīng)，溫度是最敏感的反應(yīng)參數(shù)，降低反應(yīng)溫度可顯著提高正戊烷的平衡異構(gòu)化率。

(2)對(duì)于正戊烷中溫分子篩異構(gòu)化、固體超強(qiáng)酸異構(gòu)化和低溫異構(gòu)化3大反應(yīng)體系，前兩者為氣相反應(yīng)，氣體的非理想性不明顯，可假設(shè)其為理想狀態(tài)進(jìn)行熱力學(xué)分析；后者為氣相或氣-液兩相反應(yīng)，氣體的非理想性較明顯，采用理想狀態(tài)假設(shè)進(jìn)行熱力學(xué)分析偏差較大。

(3)當(dāng)體系非理想性不明顯時(shí)，反應(yīng)壓力和氫/烴摩爾比對(duì)正戊烷平衡異構(gòu)化率的影響可以忽略；當(dāng)體系非理想性較為明顯時(shí)，反應(yīng)壓力和氫/烴摩爾比對(duì)反應(yīng)平衡的影響不可忽略。適當(dāng)降低反應(yīng)壓力或提高氫/烴摩爾比均可提高正戊烷的平衡異構(gòu)化率。

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

ai——等壓摩爾熱容表達(dá)式中經(jīng)驗(yàn)常數(shù)，J/(K·mol)；

Δa——∑νiai的簡(jiǎn)略表達(dá)，J/(K·mol)；

A——SRK狀態(tài)方程參數(shù)；

Bi——SRK方程中i的參數(shù)B；

bi——等壓摩爾熱容表達(dá)式中經(jīng)驗(yàn)常數(shù)，J/(K2·mol)；

B——SRK狀態(tài)方程參數(shù)；

ci——等壓摩爾熱容表達(dá)式中經(jīng)驗(yàn)常數(shù)，J/(K3·mol)；

Cp,i——參加反應(yīng)各物質(zhì)的等壓摩爾熱容，J/(K·mol)；

ΔrCp,m——反應(yīng)的等壓摩爾熱容變，J/(K·mol);

ΔG——吉布斯自由能變，J/mol；

ΔH——焓變，J/mol；

ΔfHm,i(T0)——為溫度為T0時(shí)參加反應(yīng)各物質(zhì)的摩爾生成焓，J/mol；

ΔrHm(T)——反應(yīng)溫度為T時(shí)的摩爾反應(yīng)焓變，J/mol；

ΔrHm(T0)——反應(yīng)溫度為T0時(shí)的摩爾反應(yīng)焓變，J/mol；

I1——積分常數(shù)，J/mol；

I2——積分常數(shù)，J/(K·mol)；

Kf——用逸度表示的平衡常數(shù)；

kij——交互作用參數(shù)；

Ki——汽化平衡比；

ni——組分i物質(zhì)的量，mol；

ni0——組分i初始物質(zhì)的量，mol；

nL——液相物質(zhì)的量，mol；

nt——總物質(zhì)的量，mol；

nt0——初始總物質(zhì)的量，mol；

nV——?dú)庀辔镔|(zhì)的量，mol；

pc——臨界壓力，Pa；

pi——組分i的分壓，Pa；

pr——對(duì)比壓力；

p——壓力，Pa；

R——?dú)怏w常數(shù)，J/(K·mol)；

T0——初始溫度，K；

Tc——臨界溫度，K；

Td——露點(diǎn)溫度，K；

Tk——設(shè)定溫度，K；

Tr——對(duì)比溫度；

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

V——摩爾體積，m3/mol；

xG——?dú)庀嗟哪柗謹(jǐn)?shù)，%；

xH2——液相中氫氣的摩爾分?jǐn)?shù)，%；

xi——液相中組分i的摩爾分?jǐn)?shù)，%；

xi-C5——液相中異戊烷的摩爾分?jǐn)?shù)，%；

xL——液相的摩爾分?jǐn)?shù)，%；

xn-C5——液相中正戊烷的摩爾分?jǐn)?shù)，%；

yH2——?dú)庀嘀袣錃獾哪柗謹(jǐn)?shù)，%；

yi——?dú)庀嘀薪M分i的摩爾分?jǐn)?shù)，%；

yi-C5——?dú)庀嘀挟愇焱榈哪柗謹(jǐn)?shù)，%；

yj——?dú)庀嘀薪M分j的摩爾分?jǐn)?shù)，%；

yn-C5——?dú)庀嘀姓焱榈哪柗謹(jǐn)?shù)，%；

z——SRK狀態(tài)方程參數(shù)；

α——SRK狀態(tài)方程參數(shù)；

γ——SRK狀態(tài)方程參數(shù)，m3/mol；

γi——SRK方程中i的參數(shù)γ，m3/mol；

ε——反應(yīng)進(jìn)度，mol；

λ——SRK狀態(tài)方程參數(shù)，Pa·(m3/mol)2；

(λα)ij——SRK方程表達(dá)混合物時(shí)組分i與組分j的交叉參數(shù)，Pa·(m3/mol)2；

(λα)i——SRK方程中組分i的參數(shù)λ與α的乘積，Pa·(m3/mol)2；

(λα)j——SRK方程中組分j的參數(shù)λ與α的乘積，Pa·(m3/mol)2；

μiL——液相中組分i的化學(xué)勢(shì)；

μiV——?dú)庀嘀薪M分i的化學(xué)勢(shì)；

μi——組分i的化學(xué)勢(shì)；

νi——化學(xué)反應(yīng)式中各組分的計(jì)量系數(shù)，對(duì)反應(yīng)物為負(fù)值，對(duì)產(chǎn)物為正值，對(duì)不參與反應(yīng)的物質(zhì)為0；

ω——偏心因子。