基于PR立方型狀態(tài)方程普遍化溫度函數(shù)的研究與評(píng)價(jià)

趙文英，李文文，孫曉巖，曹曉榮，項(xiàng)曙光

（1 齊魯師范學(xué)院化學(xué)與化工學(xué)院，山東濟(jì)南250200； 2 青島科技大學(xué)過(guò)程系統(tǒng)工程研究所，山東青島266042）

引 言

立方型狀態(tài)方程是中等至高壓體系純組分蒸汽壓和流體混合物汽液平衡及焓、熵、熱容等熱力學(xué)性質(zhì)預(yù)測(cè)的重要方法。van der Waals[1]首次提出適用于真實(shí)氣體的立方型狀態(tài)方程，Redlich 和Kwong[2]改進(jìn)VDW EoS 引力項(xiàng)的形式并將引力項(xiàng)系數(shù)a 關(guān)聯(lián)為T(mén)r0.5的溫度函數(shù)得到RK EoS，但對(duì)純物質(zhì)蒸汽壓的預(yù)測(cè)精度不高。Soave[3]發(fā)現(xiàn)烴類物質(zhì)的α0.5與(1-Tr0.5)呈線性關(guān)系，并將系數(shù)m 普遍化為偏心因子的函數(shù)，得到著名的SRK EoS，提高了對(duì)非極性和弱極性物質(zhì)蒸汽壓的預(yù)測(cè)精度，但對(duì)液體體積的預(yù)測(cè)結(jié)果偏差較大。Peng 和Robinson[4]改進(jìn)RK EoS 引力項(xiàng)的形式并重新關(guān)聯(lián)Soave[3]溫度函數(shù)的系數(shù)m，得到PR EoS，能夠準(zhǔn)確預(yù)測(cè)非極性和弱極性物質(zhì)的蒸汽壓，且提高了液體體積的預(yù)測(cè)精度。SRK 和PR EoS 對(duì)非極性和弱極性物質(zhì)蒸汽壓的預(yù)測(cè)結(jié)果均比較準(zhǔn)確，但對(duì)極性物質(zhì)蒸汽壓預(yù)測(cè)結(jié)果的偏差較大。

引力項(xiàng)中溫度函數(shù)的形式影響狀態(tài)方程對(duì)蒸汽壓的預(yù)測(cè)精度。Mathias 和Copeman[5]在Soave[3]的基礎(chǔ)上增加(1-Tr0.5)的項(xiàng)數(shù)，提高了PR EoS 對(duì)極性物質(zhì)蒸汽壓的預(yù)測(cè)精度。Stryjek 和Vera[6]將Soave系數(shù)m 關(guān)聯(lián)為T(mén)r的函數(shù)，提高了低對(duì)比溫度區(qū)域蒸汽壓的預(yù)測(cè)精度。Twu等[7-8]、Gasem等[9]、Le Guennec等[10]、Mahmoodi 等[11]、Pina-Martinez 等[12]和Ghoderao等[13]不斷研究改進(jìn)溫度函數(shù)的形式，以提高溫度函數(shù)對(duì)蒸汽壓的預(yù)測(cè)精度。Zhao等[14]對(duì)適用于立方型狀態(tài)方程的多項(xiàng)式溫度函數(shù)進(jìn)行了概述。改進(jìn)的溫度函數(shù)能夠預(yù)測(cè)純物質(zhì)的蒸汽壓等熱力學(xué)性質(zhì)，但由于溫度函數(shù)的參數(shù)與物質(zhì)一一對(duì)應(yīng)，受參數(shù)的限制，適用范圍窄。Soave[3]最早提出普遍化溫度函數(shù)，此 后Twu 等[15]、Coquelet 等[16]、汪 萍 等[17]、Li 和Yang[18]、Forero 等[19]、Le Guennec 等[10]、Mahmoodi 等[11]不斷研究溫度函數(shù)普遍化的方法，用于非極性、弱極性和極性物質(zhì)蒸汽壓的預(yù)測(cè)。

本文擬在PR EoS 的基礎(chǔ)上，用70 種非極性、弱極性和極性物質(zhì)對(duì)17 種普遍化溫度函數(shù)預(yù)測(cè)蒸汽壓的能力進(jìn)行評(píng)價(jià)，總結(jié)普遍化溫度函數(shù)的特點(diǎn)，確定較好的普遍化溫度函數(shù)的形式和特點(diǎn)，為溫度函數(shù)普遍化方法的改進(jìn)提供借鑒。

1 普遍化溫度函數(shù)模型

自Soave[3]提出普遍化溫度函數(shù)以來(lái)，研究者不斷改進(jìn)溫度函數(shù)的形式，得到上百種溫度函數(shù)[14]，但普遍化的溫度函數(shù)形式較少。通過(guò)整理，得到17種適用于PR EoS 的普遍化溫度函數(shù)，其中5 種多項(xiàng)式溫度函數(shù)，12 種指數(shù)溫度函數(shù)。5 種多項(xiàng)式溫度函數(shù)中Peng-Robinson(1976)[4]和Robinson-Peng(1978)[20]是應(yīng)用較多的普遍化溫度函數(shù)；后來(lái)，Coquelet等[16]、汪 萍 等[17]、Valiollahi 等[21]分 別 對(duì)Mathias-Copeman(1983)[5]、Stryjek-Vera(1986)[6]和Valiollahi 等(2016)[21]溫度函數(shù)的參數(shù)普遍化。普遍化的指數(shù)溫度函數(shù)種 類 較 多，Gasem (2001)[9]、Joshipura (2009)[22]、Li-Yang (2011)[18]、Haghtalab (2011)[23]、Saffari-Zahedi(2013)[24]、Hou (2015)[25]、PR-PM (2017)[11]和PR-PM2(2017)[11]是新提出且直接普遍化的溫度函數(shù)，Le Guennec 等[10]、Forero 等[19]、Mahmoodi 等[11]則 分 別 對(duì)Twu (1988)[7]、Heyen (1980)[26]、Twu (1995)[15]和Coquelet (2004)[16]溫度函數(shù)重新普遍化。普遍化的溫度函數(shù)及參數(shù)關(guān)聯(lián)式見(jiàn)表1。

表1 用于評(píng)價(jià)的普遍化溫度函數(shù)Table 1 Generalized alpha functions for evaluation

續(xù)表

續(xù)表

2 評(píng)價(jià)方法

用11類共70種物質(zhì)，在0.42～1.0對(duì)比溫度范圍內(nèi)的4453 個(gè)實(shí)驗(yàn)飽和蒸汽壓數(shù)據(jù)對(duì)17 種普遍化溫度函數(shù)進(jìn)行評(píng)價(jià)。各物質(zhì)的實(shí)驗(yàn)蒸汽壓數(shù)據(jù)參考文獻(xiàn)見(jiàn)表2，臨界性質(zhì)和偏心因子等特性參數(shù)源自《Thermophysical Properties of Chemicals and Hydrocarbons》[27]。

將物質(zhì)的特性參數(shù)帶入普遍化的溫度函數(shù)關(guān)聯(lián)式，根據(jù)汽液相平衡時(shí)汽相逸度等于液相逸度的條件，計(jì)算該溫度下的飽和蒸汽壓。對(duì)蒸汽壓的計(jì)算借助于Matlab 軟件，用roots函數(shù)求解壓縮因子(Z)的三次方程，fsolve 函數(shù)求解汽液相平衡時(shí)的蒸汽壓。

飽和蒸汽壓的相對(duì)偏差(RD)和平均相對(duì)偏差(ARD)計(jì)算表達(dá)式見(jiàn)式(1)、式(2)。

表2 11類物質(zhì)基本信息Table 2 Basic information of 11 kinds of compounds

3 結(jié)果與討論

3.1 整體分析預(yù)測(cè)結(jié)果

17 種普遍化溫度函數(shù)對(duì)11 類物質(zhì)的預(yù)測(cè)結(jié)果見(jiàn)圖1。從圖1 可以看出，Hou (2015)溫度函數(shù)對(duì)各類物質(zhì)的預(yù)測(cè)結(jié)果偏差最大，平均相對(duì)偏差為13.90%，其次為Saffari-Zahedi(2013)溫度函數(shù)，平均相對(duì)偏差為6.56%,主要用于天然氣體系的汽液平衡計(jì)算，Hou等[25]在Saffari-Zahedi(2013)的基礎(chǔ)上重新關(guān)聯(lián)溫度函數(shù)的參數(shù)，用于富含CO2天然氣體系的汽液平衡預(yù)測(cè)。Saffari-Zahedi (2013)對(duì)C8以下的烷烴和芳烴、氣體、醚、酯等小分子物質(zhì)的預(yù)測(cè)結(jié)果較準(zhǔn)確，但對(duì)極性物質(zhì)和分子量較大的非極性物質(zhì)的預(yù)測(cè)結(jié)果偏差較大。單參數(shù)的Joshipura(2009)溫度函數(shù)對(duì)蒸汽壓預(yù)測(cè)結(jié)果的偏差也較大，平均相對(duì)偏差為6.18%，其中對(duì)正構(gòu)醇、正構(gòu)烷烴、氣體和酸類物質(zhì)預(yù)測(cè)結(jié)果的平均相對(duì)偏差超過(guò)9%。預(yù)測(cè)結(jié)果最準(zhǔn)確的是Forero (2016)溫度函數(shù)，對(duì)11 類物質(zhì)蒸汽壓預(yù)測(cè)結(jié)果的平均偏差為3.09%。 Peng-Robinson (1976)、 Robinson-Peng(1978)、Gasem (2001)、Li-Yang (2011)等13 種溫度函數(shù)對(duì)蒸汽壓預(yù)測(cè)結(jié)果的平均相對(duì)偏差接近，介于3.64%～4.50%。

3.2 正構(gòu)烷烴的預(yù)測(cè)結(jié)果

17 種普遍化溫度函數(shù)對(duì)C1～C20、C22和C24共22種正構(gòu)烷烴蒸汽壓的預(yù)測(cè)結(jié)果見(jiàn)圖2。從圖2 可以看出，Hou (2015)溫度函數(shù)對(duì)烷烴蒸汽壓的預(yù)測(cè)結(jié)果整體偏差，尤其對(duì)于C9以上烴類物質(zhì)蒸汽壓的偏差顯著增大，平均相對(duì)偏差超過(guò)20%，對(duì)C18以上重?zé)N類物質(zhì)蒸汽壓的預(yù)測(cè)結(jié)果的偏差超過(guò)45%，最高達(dá)77.9%。Saffari-Zahedi (2013)對(duì)烷烴蒸汽壓的預(yù)測(cè)結(jié)果略優(yōu)于Hou(2015)溫度函數(shù)，自正壬烷開(kāi)始，預(yù)測(cè)結(jié)果的偏差顯著增大，大于5%。開(kāi)發(fā)用于天然氣露點(diǎn)計(jì)算的Valiollahi(2016)溫度函數(shù)，對(duì)C18以上烷烴預(yù)測(cè)結(jié)果的最大偏差小于7.5%，但對(duì)C18以上重?zé)N預(yù)測(cè)結(jié)果的偏差顯著增大，高達(dá)29.4%；對(duì)22種烴類物質(zhì)蒸汽壓預(yù)測(cè)結(jié)果的平均相對(duì)偏差為6.28%，與Peng-Robinson (1976)溫度函數(shù)的預(yù)測(cè)結(jié)果接近(6.27%)。Joshipura (2009)溫度函數(shù)對(duì)22 種烴類物質(zhì)蒸汽壓預(yù)測(cè)結(jié)果的平均相對(duì)偏差為9.13%，高于Peng-Robinson (1976)溫度函數(shù)的預(yù)測(cè)結(jié)果。剩余12 種溫度函數(shù)對(duì)C8以下烴類物質(zhì)蒸汽壓預(yù)測(cè)結(jié)果的偏差接近，平均相對(duì)偏差均小于2.5%，但對(duì)C8以上烴類物質(zhì)的預(yù)測(cè)結(jié)果產(chǎn)生差異。PR-PM (2017)、PR-PM2 (2017)、Twu-Mahmoodi(2017)、Le Guennec (2016)、Haghtalab (2011),Gasem(2001)溫度函數(shù)對(duì)C8～C16之間的烷烴蒸汽壓的預(yù)測(cè)結(jié)果的平均相對(duì)偏差小于3%。Coquelet (2004)和Robinson-Peng (1978)溫度函數(shù)對(duì)C16以上烴類蒸汽壓的預(yù)測(cè)結(jié)果最準(zhǔn)確，平均相對(duì)偏差小于2%。由于汪萍等(2004)溫度函數(shù)的參數(shù)關(guān)聯(lián)式用C1～C16正構(gòu)烷烴的實(shí)驗(yàn)蒸汽壓數(shù)據(jù)擬合，對(duì)于22種正構(gòu)烷烴的預(yù)測(cè)結(jié)果最準(zhǔn)確，平均相對(duì)偏差為1.43%，其次為Robinson-Peng (1978) (1.86%)， Gasem (2001),Coquelet (2004)、Li-Yang (2011)、Haghtalab (2011)、Le Guennec (2016)、Forero (2016)、Twu-Mahmoodi(2017)、Coquelet-Mahmoodi (2017)、PR-PM (2017)和PR-PM2(2017)溫度函數(shù)對(duì)正構(gòu)烷烴蒸汽壓的預(yù)測(cè)結(jié)果介于2.05%～2.47%，預(yù)測(cè)結(jié)果接近。

圖1 參數(shù)普遍化溫度函數(shù)對(duì)不同類別物質(zhì)蒸汽壓預(yù)測(cè)結(jié)果的平均相對(duì)偏差Fig.1 ARDs of vapor pressures for different kinds of compounds with generalized alpha functions

以正丁烷為例，研究各溫度函數(shù)對(duì)蒸汽壓的預(yù)測(cè)偏差隨對(duì)比溫度的變化趨勢(shì)(圖3)。除Joshipura(2009)和Hou (2015)溫度函數(shù)外，其他溫度函數(shù)在0.6≤Tr＜1 溫度范圍內(nèi)對(duì)正丁烷的預(yù)測(cè)結(jié)果都比較準(zhǔn)確，相對(duì)偏差小于±1%；當(dāng)Tr＜0.6 時(shí)，隨著對(duì)比溫度的下降，溫度函數(shù)的預(yù)測(cè)結(jié)果相對(duì)偏差迅速增大，最高可達(dá)到15%。由此可見(jiàn)，對(duì)于烷烴等非極性物質(zhì)，普遍化溫度函數(shù)對(duì)較低溫度和較小壓力下蒸汽壓的預(yù)測(cè)偏差較大，在應(yīng)用時(shí)要注意溫度函數(shù)的應(yīng)用范圍。

總的來(lái)看，汪萍(2004)、Robinson-Peng (1978)、Gasem (2001)、Coquelet (2004)、Li-Yang (2011)、Haghtalab (2011)、Le Guennec (2016)、Forero (2016)、Twu-Mahmoodi (2017)、Coquelet-Mahmoodi (2017)、PR-PM (2017)和PR-PM2 (2017)溫度函數(shù)對(duì)正構(gòu)烷烴蒸汽壓的預(yù)測(cè)結(jié)果都比較準(zhǔn)確。

3.3 正構(gòu)醇的預(yù)測(cè)結(jié)果

17 種普遍化溫度函數(shù)對(duì)醇類物質(zhì)蒸汽壓預(yù)測(cè)結(jié)果的平均相對(duì)偏差見(jiàn)圖4。從圖4中可以看出，普遍化溫度函數(shù)對(duì)醇類物質(zhì)蒸汽壓的預(yù)測(cè)結(jié)果偏差明顯較大。Saffari-Zahedi (2013)和Hou (2015)溫度函數(shù)預(yù)測(cè)結(jié)果偏差的變化趨勢(shì)先略減小至平緩而后顯著增大，平均相對(duì)偏差分別為13.79% 和23.84%，其中Hou (2015)溫度函數(shù)預(yù)測(cè)結(jié)果的偏差最大。Forero (2016)溫度函數(shù)對(duì)醇類物質(zhì)蒸汽壓的預(yù)測(cè)結(jié)果最準(zhǔn)確，平均相對(duì)偏差為4.16%，顯著優(yōu)于其余溫度函數(shù)。Forero (2016)溫度函數(shù)對(duì)碳數(shù)小于C8醇類蒸汽壓預(yù)測(cè)結(jié)果的平均相對(duì)偏差小于4.2%，但對(duì)正壬醇和正十四醇的預(yù)測(cè)結(jié)果偏差較大，平均相對(duì)偏差分別為10.92%和10.78%，對(duì)正癸醇、正十二醇、正十六醇和正十八醇蒸汽壓預(yù)測(cè)結(jié)果的平均相對(duì)偏差介于2.18%～6.90%。除Saffari-Zahedi(2013)、Hou (2015)和Forero (2016) 3 種溫度函數(shù)外，其余溫度函數(shù)對(duì)醇類物質(zhì)蒸汽壓預(yù)測(cè)結(jié)果的偏差均呈現(xiàn)先增大后減小的趨勢(shì)，且預(yù)測(cè)結(jié)果的偏差接近，平均相對(duì)偏差介于9.07%～17.38%。

圖2 參數(shù)普遍化溫度函數(shù)對(duì)正構(gòu)烷烴蒸汽壓預(yù)測(cè)結(jié)果的平均相對(duì)偏差Fig.2 ARDs of vapor pressures for normal alkanes with generalized alpha functions

圖3 參數(shù)普遍化溫度函數(shù)對(duì)正丁烷蒸汽壓預(yù)測(cè)結(jié)果的相對(duì)偏差隨對(duì)比溫度變化趨勢(shì)Fig.3 RDs of generalized alpha functions on prediction of vapor pressures of n-butane with increasing reduced temperature

圖4 參數(shù)普遍化溫度函數(shù)對(duì)正構(gòu)醇蒸汽壓預(yù)測(cè)結(jié)果的平均相對(duì)偏差Fig.4 ARDs of vapor pressures for normal alcohols with generalized alpha functions

為了進(jìn)一步分析普遍化溫度函數(shù)對(duì)醇類物質(zhì)預(yù)測(cè)結(jié)果偏差較大的原因，以正丁醇為例，觀察蒸汽壓預(yù)測(cè)結(jié)果的相對(duì)偏差隨溫度變化的趨勢(shì)(圖5)。Forero (2016)溫度函數(shù)對(duì)正丁醇蒸汽壓預(yù)測(cè)結(jié)果在Tr=0.525～0.993溫度范圍內(nèi)偏差較小，相對(duì)偏差最大值為2.87%；Peng-Robinson (1976)等溫度函數(shù)的預(yù)測(cè)結(jié)果在Tr＜0.7 的低對(duì)比溫度范圍內(nèi)顯著增大，偏差最大值達(dá)60%。Saffari-Zahedi(2013)溫度函數(shù)對(duì)正丁醇蒸汽壓的預(yù)測(cè)結(jié)果比較準(zhǔn)確，平均相對(duì)偏差為6.52%。Saffari-Zahedi (2013)溫度函數(shù)對(duì)正丁醇至正十二醇的蒸汽壓預(yù)測(cè)結(jié)果優(yōu)于除Forero(2016)以外的其他溫度函數(shù)。

總的來(lái)看，普遍化溫度函數(shù)對(duì)醇類物質(zhì)蒸汽壓的預(yù)測(cè)結(jié)果均較大，只有Forero (2016)溫度函數(shù)的預(yù)測(cè)結(jié)果較準(zhǔn)確。通過(guò)分析溫度函數(shù)參數(shù)關(guān)聯(lián)式可以看出，F(xiàn)orero(2016)溫度函數(shù)是以偏心因子和極性因子為變量，通過(guò)極性因子校正Tr＜0.6 溫度范圍蒸汽壓的變化趨勢(shì)，提高了對(duì)低對(duì)比溫度蒸汽壓的預(yù)測(cè)精度。Valiollahi(2016)溫度函數(shù)以對(duì)比密度為變量普遍化溫度函數(shù)的參數(shù)，雖然普遍化方法改進(jìn)，但不適用于醇類等極性物質(zhì)蒸汽壓的預(yù)測(cè)。除上述兩種溫度函數(shù)外，PR-PM (2017)等溫度函數(shù)均采用偏心因子對(duì)參數(shù)普遍化。由此可見(jiàn)，僅用偏心因子普遍化溫度函數(shù)參數(shù)的方法不適用于醇類等極性物質(zhì)蒸汽壓的預(yù)測(cè)。

3.4 芳烴、鹵代烴、氣體和雜環(huán)類物質(zhì)的預(yù)測(cè)結(jié)果

普遍化溫度函數(shù)對(duì)芳烴的預(yù)測(cè)結(jié)果比較準(zhǔn)確，優(yōu)于正構(gòu)烷烴的預(yù)測(cè)結(jié)果，這是因?yàn)閰⑴c評(píng)價(jià)的芳烴類物質(zhì)的分子較小，分子最大的1-甲基萘也僅含有11 個(gè)碳。比較普遍化溫度函數(shù)對(duì)芳烴蒸汽壓的預(yù)測(cè)結(jié)果可以看出(圖1)，F(xiàn)orero (2016)溫度函數(shù)對(duì)芳烴蒸汽壓預(yù)測(cè)結(jié)果的偏差最小，平均相對(duì)偏差為1.40%，其次為PR-PR(2017)、PR-PM2(2017)和汪萍(2004)，分別為1.64%、1.69%和1.69%。Hou (2015)溫度函數(shù)的預(yù)測(cè)結(jié)果偏差最大，平均相對(duì)偏差為8.14%，其次為Joshipura(2009)和Coquelet(2004)，分別為3.82%和3.19%。其余溫度函數(shù)預(yù)測(cè)結(jié)果的平均相對(duì)偏差介于1.83%～2.93%之間。

圖5 參數(shù)普遍化溫度函數(shù)對(duì)正丁醇蒸汽壓預(yù)測(cè)結(jié)果的相對(duì)偏差隨對(duì)比溫度的變化趨勢(shì)Fig.5 RDs of generalized alpha functions on prediction of vapor pressures of n-butanol with increasing reduced temperature

比較圖1 普遍化溫度函數(shù)鹵代烴蒸汽壓的預(yù)測(cè)結(jié)果可以看出，只有Forero (2016)溫度函數(shù)預(yù)測(cè)結(jié)果的偏差最大，高達(dá)18.62%。Haghtalab (2011)溫度函數(shù)預(yù)測(cè)結(jié)果的偏差最小為1.83%，其余溫度函數(shù)的預(yù)測(cè)結(jié)果接近，平均相對(duì)偏差介于2.09%～4.67%。

普遍化溫度函數(shù)對(duì)小分子雜環(huán)類物質(zhì)的預(yù)測(cè)結(jié)果比較準(zhǔn)確，F(xiàn)orero(2016)、Twu-Mahmoodi(2017)、Coquelet-Mahmoodi (2017)、PR-PM (2017) 和PRPM2 (2017)預(yù)測(cè)結(jié)果的平均相對(duì)偏差均小于1%。比較普遍化溫度函數(shù)對(duì)氣體蒸汽壓預(yù)測(cè)結(jié)果比較準(zhǔn)確，除Joshipura (2009)和Hou (2015)溫度函數(shù)，其余溫度函數(shù)對(duì)氣體蒸汽壓預(yù)測(cè)結(jié)果的偏差介于1.4%～2.0%。

3.5 醚、酮、酯、酸類物質(zhì)的預(yù)測(cè)結(jié)果

17 種普遍化溫度函數(shù)對(duì)弱極性的醚、酮、酯類物質(zhì)蒸汽壓的預(yù)測(cè)結(jié)果比較準(zhǔn)確。Hou(2015)溫度函數(shù)對(duì)蒸汽壓預(yù)測(cè)結(jié)果的偏差最大，平均相對(duì)偏差分別為10.21%、10.86%、10.48%。其次Joshipura(2009)的偏差較大，平均相對(duì)偏差分別為4.47%、3.71%、3.35%。對(duì)醚、酮、酯類物質(zhì)預(yù)測(cè)結(jié)果最為準(zhǔn)確的是Forero (2016)溫度函數(shù)，平均相對(duì)偏差分別為1.31%、0.84%、1.22%。其次為Coquelet (2004)溫度函數(shù)，平均相對(duì)偏差分別為1.48%、1.31%、2.79%。其余溫度函數(shù)對(duì)醚、酮、酯類物質(zhì)蒸汽壓的預(yù)測(cè)結(jié)果介于0.7%～3.7%之間，預(yù)測(cè)結(jié)果接近。

普遍化溫度函數(shù)對(duì)酸類物質(zhì)的預(yù)測(cè)結(jié)果偏差較大，主要誤差來(lái)源于乙酸。普遍化溫度函數(shù)對(duì)乙酸蒸汽壓預(yù)測(cè)結(jié)果的最大偏差約為丙酸的10 倍。由于乙酸分子之間存在強(qiáng)締合作用，其性質(zhì)與丙酸、醇類等物質(zhì)性質(zhì)產(chǎn)生明顯的差異。難以用偏心因子普遍化的單參數(shù)或兩參數(shù)溫度函數(shù)準(zhǔn)確預(yù)測(cè)其蒸汽壓。大多數(shù)溫度函數(shù)對(duì)醇類物質(zhì)蒸汽壓預(yù)測(cè)結(jié)果的平均相對(duì)偏差在9%左右。其中以偏心因子和極性因子普遍化的Forero(2016)溫度函數(shù)對(duì)酸類物質(zhì)蒸汽壓的預(yù)測(cè)結(jié)果最準(zhǔn)確，平均相對(duì)偏差為1.07%，其次為Coquelet (2004)溫度函數(shù)，平均相對(duì)偏 差 為1.98%，Robinson-Peng (1978) 和Peng-Robinson (1976)預(yù)測(cè)結(jié)果的偏差分別為3.58%和6.41%。Saffari-Zahedi(2013)和Hou(2015)溫度函數(shù)預(yù)測(cè)結(jié)果的偏差最大，分別為19.31%和31.94%。

3.6 水的預(yù)測(cè)結(jié)果

普遍化溫度函數(shù)對(duì)水蒸氣壓預(yù)測(cè)結(jié)果的偏差較大。只有Forero(2016)溫度函數(shù)對(duì)水蒸氣壓的預(yù)測(cè)結(jié)果小于1%，為0.77%。其次為Robinson-Peng(1978)，平均相對(duì)偏差為2.15%，Li-Yang (2011)溫度函數(shù)的預(yù)測(cè)結(jié)果為3.55%。Le Guennec (2016)、Twu-Mahmoodi (2017)、Coquelet-Mahmoodi (2017)、PR-PM (2017)和PR-PM2 (2017) 5 種新普遍化溫度函數(shù)預(yù)測(cè)結(jié)果的平均相對(duì)偏差均在6%左右，偏差較大。由此可見(jiàn)，僅用偏心因子普遍化的溫度函數(shù)對(duì)水等極性物質(zhì)蒸汽壓預(yù)測(cè)結(jié)果的偏差較大。

4 結(jié) 論

11 類70 種物質(zhì)的4453 個(gè)實(shí)驗(yàn)蒸汽壓數(shù)據(jù)對(duì)適用于PR EoS 的17 種普遍化溫度函數(shù)的預(yù)測(cè)能力進(jìn)行了評(píng)價(jià)，結(jié)論如下：

(1) 以偏心因子和極性因子普遍化的Forero(2016)溫度函數(shù)對(duì)11類物質(zhì)預(yù)測(cè)結(jié)果的平均相對(duì)偏差最小，為3.09%；對(duì)醇、酸、醚、酯、水等弱極性和極性物質(zhì)蒸汽壓的預(yù)測(cè)結(jié)果的平均相對(duì)偏差分別為4.16%、1.07%、1.31%、1.22%和0.77%，明顯優(yōu)于以偏心因子普遍化溫度函數(shù)的預(yù)測(cè)結(jié)果。Forero(2016)溫度函數(shù)對(duì)烷烴、氣體和雜環(huán)類物質(zhì)蒸汽壓的預(yù)測(cè)精度與以偏心因子普遍化溫度函數(shù)的預(yù)測(cè)結(jié)果接近，但不適用于鹵代烴蒸汽壓的預(yù)測(cè)。

(2) 汪萍(2004)和Robinson-Peng (1978)溫度函數(shù)對(duì)正構(gòu)烷烴蒸汽壓預(yù)測(cè)結(jié)果的平均相對(duì)偏差分別為1.43%和1.86%，Li-Yang(2011)對(duì)氮?dú)狻⒁谎趸肌钡葰怏w預(yù)測(cè)結(jié)果的平均相對(duì)偏差分別為1.45%。PR-PM (2017)對(duì)芳烴預(yù)測(cè)結(jié)果的評(píng)價(jià)相對(duì)偏差為1.64%，僅次于Forero (2016)溫度函數(shù)。Haghtalab (2011)對(duì)鹵代烴的預(yù)測(cè)結(jié)果最準(zhǔn)確，平均相對(duì)偏差為1.83%。

(3)僅以偏心因子普遍化的溫度函數(shù)能夠準(zhǔn)確預(yù)測(cè)非極性和弱極性物質(zhì)的飽和蒸汽壓，但對(duì)極性物質(zhì)飽和蒸汽壓的預(yù)測(cè)結(jié)果偏差較大。對(duì)適用于極性物質(zhì)的溫度函數(shù)普遍化，需要表征極性性質(zhì)的極性因子、偶極矩等特性參數(shù)參與擬合，以提高極性物質(zhì)蒸汽壓預(yù)測(cè)結(jié)果的準(zhǔn)確度。

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

ARD——平均相對(duì)偏差

a——引力項(xiàng)參數(shù)

EoS——狀態(tài)方程

P——蒸汽壓，kPa

PR——Peng-Robinson

RD——相對(duì)偏差

RK——Redlich-Kwong

SRK——Soave-Redlich-Kwong

T——溫度，K

α——引力項(xiàng)中的溫度函數(shù)

ρ——密度，kg·m-3

ω——偏心因子

上下標(biāo)

c——臨界性質(zhì)

exp——實(shí)驗(yàn)值

sat——飽和狀態(tài)

r——臨界參數(shù)的對(duì)比態(tài)