2界面張力模擬近似計算

在懸浮床加氫體系中，原料油的界面張力的大小直接影響體系中氣泡的大小和分布范圍及尺寸大小，這對于體系中氣體如H2從體相向催化劑表面、油品表面的擴(kuò)散有很大的影響，隨著溫度的變化，其界面張力也呈現(xiàn)出不同的變化規(guī)律。因此，計算工況條件下渣油在H2氣氛下的界面張力變化對于建立渣油懸浮床加氫氣泡尺寸模型尤為重要。

但是，目前關(guān)于氣體/渣油界面張力的研究集中在CO2/N2等氣體驅(qū)油的領(lǐng)域中，關(guān)于加氫領(lǐng)域中尤其是H2氣氛下渣油的界面張力研究還處于空白狀態(tài)。因此，本文主要以與H2分子結(jié)構(gòu)較為相似的N2為主要研究考察對象，考察N2氣氛下渣油的界面張力。另外，對于不同氣氛下的渣油的界面張力變化來說，氣體的沸點(diǎn)越低，其在高壓下溶解入渣油之中導(dǎo)致界面張力變小的幅度變化越小。

查表可知，在標(biāo)準(zhǔn)狀態(tài)下，CO2沸點(diǎn)-56.5℃，N2沸點(diǎn)-195.6℃，H2沸點(diǎn)-252.7℃，因此，相同壓力下，H2氣氛下降低的幅度應(yīng)該比N2氣氛相對應(yīng)的界面張力下降幅度小，因此，可以得到在高壓條件下的H2/渣油的界面張力應(yīng)當(dāng)處于該渣油在N2氣氛下常壓以及高壓條件下界面張力之間。

根據(jù)上述分析，并考慮到在較高溫度下渣油組成的劇烈變化導(dǎo)致的計算偏差情況，可以計算得到N2在常壓條件下以及高壓條件下的界面張力，從而對H2在高壓條件下的界面張力進(jìn)行限定。

2.1 N2氣氛常壓條件下的界面張力計算

首先，界面張力的模擬計算與油品的餾程等有很重要的關(guān)系，需要通過油品的平均沸點(diǎn)等數(shù)據(jù)，可以根據(jù)油品的蒸餾曲線得到的平均沸點(diǎn)數(shù)據(jù)。體積平均沸點(diǎn)，是平均沸點(diǎn)的表示方法之一，是恩氏蒸餾10%、30%、50%、70%、90%五個餾出溫度的算術(shù)平均值。用于求定其他物理常數(shù)。

由于減壓渣油過重，超過質(zhì)量分?jǐn)?shù)80%的油品在實(shí)際操作過程中是無法餾出，這給計算油品的平均沸點(diǎn)造成了困難。考慮到目前得到的油品蒸餾曲線均為典型的S型曲線，因此將現(xiàn)有曲線采用Boltzmann計算模型來進(jìn)行模擬擬合。

通過Boltzmann計算模型得到的紅色擬合曲線的相關(guān)性Adj R2=0.999 76≈1，說明擬合曲線擬合優(yōu)度非常好，基本可以用擬合曲線作為真實(shí)蒸餾曲線得到平均沸點(diǎn)數(shù)據(jù)。

有研究者提供了一種利用原料油各項(xiàng)可測理化參數(shù)進(jìn)行空氣氣氛、常壓條件下界面張力估測的方法，在實(shí)際應(yīng)用中被證明具有較高的準(zhǔn)確度，本研究將以此為基礎(chǔ)得到最初的界面張力數(shù)據(jù)，隨后進(jìn)行接下來的計算近似。

根據(jù)表一提供的Toledo減渣的物性數(shù)據(jù)，以及得到的減壓渣油的模擬蒸餾曲線，可得到Toledo減渣特性因數(shù)為11.6，真臨界溫度為760℃，據(jù)此查表計算可得在常壓、較高溫度條件下常壓條件、N2氣氛下Toledo減渣的氣液界面的界面張力，逐點(diǎn)計算后將結(jié)果匯總。

可以看到，根據(jù)文獻(xiàn)方法計算結(jié)果，對于Toledo減壓渣油來說，其在常壓條件下N2氣氛條件下的氣液界面的界面張力，當(dāng)不在真臨界溫度（760℃）較近的溫度范圍內(nèi)，都可以認(rèn)為隨溫度呈現(xiàn)線性下降的趨勢。

對于接近真臨界溫度范圍后的界面張力的變化，多篇文獻(xiàn)報道了較較低溫度范圍斜率更低的線性變化關(guān)系，直至在真臨界溫度條件下，氣液兩相體系最終演變?yōu)闅庖翰糠煮w系，此時體系的不存在界面，則顯然界面張力數(shù)據(jù)為零，與計算結(jié)果趨勢一致。

2.2 N2氣氛高壓條件下的界面張力估算

在較高壓力條件下，氣液界面中氣相一邊的壓力越大，氣體分子從氣相向液相中擴(kuò)散的動力也就越大，最終導(dǎo)致一部分氣體會溶解在渣油之中，由此導(dǎo)致氣液界面的界面張力下降，一般來說，壓力越高，溶解氣體越多，氣液界面的界面張力就越低。對于氮?dú)鈦碚f，文獻(xiàn)中報道了其在不同溫度下的壓力-界面張力關(guān)系，根據(jù)結(jié)果發(fā)現(xiàn)，在一定溫度下，常壓-20 MPa條件下的渣油界面張力與壓力為線性下降關(guān)系。

可以發(fā)現(xiàn)，不同溫度條件下渣油界面張力線性下降的斜率與溫度存在較為明顯的線性關(guān)系，即，溫度越高，則氣液兩相界面張力隨氣相壓力升高而降低的趨勢則越小，且呈線性變化趨勢。因此，根據(jù)此結(jié)果，可以以較低溫度條件下的氣液界面界面張力關(guān)系，推導(dǎo)在較高溫度條件下的氣液界面界面張力隨外界條件的變化關(guān)系。

不同溫度條件下的線性下降的斜率與溫度的變化也基本呈現(xiàn)出一定的的線性。另外，對于不同壓力下的界面張力-溫度關(guān)系，理論上應(yīng)當(dāng)都是線性變化，因此，斜率變化的線性變化規(guī)律也是符合一般規(guī)律的。由此可以推測得到在相應(yīng)溫度條件下的線性下降斜率，從而由N2常壓條件下Toledo減渣的界面張力推測得到N2在20 MPa氣氛條件下Toledo減渣的界面張力。

綜合上述分析，根據(jù)計算結(jié)果，黑色、紅色分別為常壓、20 MPa下N2與減壓渣油的界面張力隨溫度變化的變化趨勢，與文獻(xiàn)規(guī)律吻合，均呈現(xiàn)線性變化規(guī)律。

對于根據(jù)前文對于不同氣體在高壓條件下氣液界面的界面張力隨溫度變化的規(guī)律的推導(dǎo)，可知，H2氣氛條件下其與減壓渣油的界面張力隨溫度變化的規(guī)律曲線應(yīng)處在N2常壓以及20 MPa條件下的界面張力-溫度兩條曲線之間，即圖中的黑色直線以及紅色直線之間的陰影區(qū)域。

可以看到，對于本文給定的界面張力計算模擬的范圍來說，同時考慮到數(shù)量級，誤差基本可以控制在±3 mN/m范圍內(nèi)，另外考慮到計算過程的偏差，具有較大的準(zhǔn)確性，這一結(jié)果具有一定的可參考價值。

3結(jié)論

本文建立了一種估算在較為苛刻的條件下，減壓渣油在氫氣氣氛條件下的黏度、界面張力等隨溫度變化的關(guān)系的方法，在缺少相應(yīng)的測量方法的條件下，可以對油品的相應(yīng)參數(shù)進(jìn)行簡單的估算。

總體來說，油品的黏度、界面張力隨著溫度的變化都呈現(xiàn)出一定的變化規(guī)律，根據(jù)相關(guān)工作的經(jīng)驗(yàn)性公式以及相關(guān)模型可以對于油品在工況條件下的黏溫曲線以及界面張力進(jìn)行一定程度的預(yù)測。對于界面張力計算模擬的范圍來說，同時考慮到數(shù)量級，基本在±3 mN/m范圍內(nèi)，另外考慮到計算結(jié)果的誤差性，這一結(jié)果基本可以參考;對于黏溫關(guān)系來說，則存在數(shù)倍的差距，這一工作有實(shí)測工作的限制因素，以及工況條件下的減壓渣油組成變化等因素的影響同時存在。

懸浮床加氫工藝條件下界面張力、油品黏度模擬近似計算（一）

懸浮床加氫工藝條件下界面張力、油品黏度模擬近似計算（二）