2界面張力模擬近似計算

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

3結論

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

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

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

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