制冷劑氣體水合物相平衡分解條件預(yù)測

楊行 李璞 張龍明 李娜（西安交通大學(xué)化工學(xué)院熱流科學(xué)與工程教育部重點實驗室 西安 710049）

制冷劑氣體水合物相平衡分解條件預(yù)測

楊行 李璞 張龍明 李娜
（西安交通大學(xué)化工學(xué)院熱流科學(xué)與工程教育部重點實驗室 西安 710049）

本文從理論方面探究制冷劑氣體水合物的相平衡分解條件，應(yīng)用經(jīng)典的van derWaals-Platteeuw水合物熱力學(xué)模型預(yù)測了R22，R23，R125和R143a水合物的分解條件。在模型預(yù)測過程中，應(yīng)用SRK狀態(tài)方程對氣相和液相進(jìn)行了模擬計算。該模型預(yù)測結(jié)果與實驗數(shù)據(jù)誤差分別為1.21%，2.84%，2.23%和1.02%，并得到了制冷劑水合物相平衡圖及四相平衡點。同時對制冷劑水合物分解熱進(jìn)行了計算，對比發(fā)現(xiàn)Ⅱ型制冷劑水合物的分解熱大于Ⅰ型。

氣體水合物；制冷劑；相平衡；熱力學(xué)模型

隨著經(jīng)濟(jì)的快速發(fā)展，電力系統(tǒng)的供需矛盾日趨嚴(yán)重，白天用電高峰期電力資源緊缺，夜晚用電低谷期電力資源過剩，電力資源得不到有效利用。電力供需矛盾的問題為蓄冷空調(diào)技術(shù)提供了廣闊的發(fā)展前景，蓄冷工質(zhì)的研究對蓄冷空調(diào)技術(shù)的發(fā)展至關(guān)重要，常用的蓄冷工質(zhì)如冰、共晶鹽和水等，它們的缺點分別是蓄冷效率低、蓄冷密度低和熱交換效率低［1-2］。研究發(fā)現(xiàn)，制冷劑氣體水合物在冰點以上發(fā)生相變（5～12℃），具有較大的蓄冷能力（蓄冷密度與冰相近）。自從被提出應(yīng)用于蓄冷空調(diào)技術(shù)以來，制冷劑氣體水合物作為一種理想的蓄冷工質(zhì)而受到廣泛的關(guān)注和研究［3］。氣體水合物是一種籠型包合物，由水分子（主體分子）與其他氣體分子（客體分子）組成。水分子由于氫鍵作用在高壓、低溫條件下形成大小不同的多面體空穴，氣體分子如甲烷、乙烷及其他碳?xì)浠衔镒鳛榭腕w分子填充在空穴之中，形成一種穩(wěn)定的水合物結(jié)構(gòu)［4-5］。根據(jù)水分子的空間分布特征以及客體水合物分子的大小可以把目前已發(fā)現(xiàn)的氣體水合物晶體結(jié)構(gòu)分為三類，即Ⅰ型、Ⅱ型和H型［6］。

相平衡熱力學(xué)是制冷劑氣體水合物研究的一個重要領(lǐng)域，熱力學(xué)研究的主要目標(biāo)是獲得水合物的相平衡數(shù)據(jù)及其穩(wěn)定存在的條件范圍，故相平衡熱力學(xué)的研究有助于制冷劑氣體水合物應(yīng)用于基于水合物蓄冷的蓄冷空調(diào)技術(shù)。圖1所示為制冷劑氣體水合物相平衡示意圖［7］，圖中R，W和H分別表示制冷劑、水、水合物三種物質(zhì)，G，L和S分別表示氣、液、固三種相態(tài)，三平衡線的交點Q1和Q2分別為 R（G）-W（S）-H-W（L）和 R（G）-R （L）-H-W（L）四相共存點。CQ2和Q2D線近似可以認(rèn)為是制冷劑的氣液平衡線（水合物相和水相中水的蒸氣壓很低，可以忽略）。制冷劑氣體水合物平衡相圖的R（G）-H-W（S）和R（G）-H-W（L）為兩條三相平衡線，即圖中的AQ1和Q1Q2線，對這兩條平衡線進(jìn)行模型預(yù)測可以得到制冷劑水合物的分解條件。

圖1　制冷劑氣體水合物相平衡示意圖Fig.1 Phase equilibrium diagram of refrigerant gas hydrates

Javanmardi J等［8］通過實驗獲得了 R22的相平衡數(shù)據(jù)，并與理論模型進(jìn)行了對比驗證，結(jié)果表明高濃度的丙酮對R22水合物形成有抑制作用。Hashimoto S等［9］對制冷劑R134a，R125和R143a水合物進(jìn)行了相平衡實驗，獲得了這三種水合物三相平衡數(shù)據(jù)以及四相平衡點Q2。Liang Deqing等［10］利用恒溫操作法對制冷劑R134a，R141b和R152a氣體水合物的熱力學(xué)生成條件進(jìn)行較完整的實驗測量，首次獲得R141b和R152a水合物在冰點以下的三相平衡數(shù)據(jù)以及四相平衡點Q2。Eslamimanesh A等［11］提出了一種簡化的熱力學(xué)模型計算R152a，R141b和R32水合物的分解條件，該模型在vdW-P模型的基礎(chǔ)上假設(shè)氣相為理想氣體，并忽略氣相中的水分和水相中的氣體組分，結(jié)果表明該模型可以較合理地預(yù)測制冷劑氣體水合物的分解條件。

本文利用經(jīng)典的vdW-P氣體水合物熱力學(xué)模型結(jié)合 SRK狀態(tài)方程模擬預(yù)測 R22，R23，R125和R143a等制冷劑氣體水合物的分解條件，并繪制了平衡相圖，根據(jù)相圖計算得到了氣體水合物的四相平衡點Q2，并利用Clapeyron方程計算得到了制冷劑氣體水合物的分解熱值。結(jié)合本文的研究工作，可以對制冷劑氣體水合物進(jìn)行模型預(yù)測，有助于基于水合物蓄冷空調(diào)系統(tǒng)的設(shè)計。

1　熱力學(xué)模型

1.1vdW-P模型

van der Waals和Platteeuw利用巨正則配分函數(shù)描述氣體分子在水合物晶體孔穴中的分布情況，建立了具有統(tǒng)計熱力學(xué)基礎(chǔ)的相平衡理論模型（vdW-P模型），他們因此成為水合物熱力學(xué)理論的創(chuàng)始人，vdW-P模型也被看作是將統(tǒng)計熱力學(xué)用于描述實際體系的成功范例，到目前為止幾乎所有的氣體水合物相平衡計算模型都是在vdW-P模型的基礎(chǔ)上發(fā)展起來的［12］。

根據(jù)相平衡理論，處于相平衡體系中的每個組分在各相中的化學(xué)勢是相等的。當(dāng)制冷劑氣體水合物處于R（G）-H-W（L/S）三相平衡時，水在水合物相中的化學(xué)勢等于其在富水相（液相水或冰）中的化學(xué)勢，即:

以水在空水合物晶格（水合物晶格孔穴未被占據(jù)的假定狀態(tài)）中的化學(xué)勢μβ為參考態(tài)，相平衡約束條件可以寫為:

vdW-P模型利用統(tǒng)計熱力學(xué)中的巨正則配分函數(shù)理論推導(dǎo)得到水合物相中水的參考化學(xué)勢的計算公式。在模型的推導(dǎo)過程中，van der Waals和Platteeuw對氣體水合物的形成過程做了如下假定［12］:

1）氣體分子占據(jù)孔穴后不會引起水合物晶格的變形，也就是說水分子對水合物自由能的貢獻(xiàn)與氣體分子的填充狀況、尺寸大小和種類無關(guān)；

2）每個孔穴最多容納一個氣體分子，并且氣體分子之間不交換位置；

3）不同孔穴中的氣體分子不存在相互作用，氣體分子只與包絡(luò)在其周圍的水分子存在相互作用；

4）模型推導(dǎo)過程不考慮量子效應(yīng)；

5）氣體分子內(nèi)運動配分函數(shù)和理想氣體分子相同；

6）用球形引力勢描述氣體分子在孔穴中的位能，即把包絡(luò)在氣體分子周圍的水分子均勻分散在圓球面上。

基于這些假設(shè)，van der Waals和Platteeuw推導(dǎo)出水在水合物相中的參考化學(xué)勢計算公式為:

式中:R為氣體常數(shù)；T為水合物體系溫度，K；νm為水合物晶胞中m型孔穴的數(shù)目與晶胞中水分子個數(shù)的比值，例如，對于Ⅱ型氣體水合物（8M·16N· 136H2O）而言，大孔穴的νm為8/136；N為生成水合物的氣體組分?jǐn)?shù)；θjm為氣體分子j在m型孔穴中的占有率；fj為氣體分子j在水合物相中的逸度，可以利用一定的狀態(tài)方程進(jìn)行計算，Pa；Cjm為Langmuir常數(shù)，這一常數(shù)反映了水合物晶體中孔穴對于氣體分子吸附能力的大小，Cjm僅是溫度的函數(shù)。而富水相（液態(tài)水或冰）中水的參考化學(xué)勢可以用以下公式計算得到［12］:

式中:ΔhW，ΔvW和Δcp分別表示空水合物晶格與純水的摩爾焓差、摩爾體積差、摩爾比熱容差；上標(biāo)0是指在溫度T0=273.15 K時所對應(yīng)的值；λ為比熱容的溫度系數(shù)。對于不同結(jié)構(gòu)的氣體水合物而言，和λ等參數(shù)數(shù)據(jù)有所不同，如表1所示。

表1?、裥?、Ⅱ型氣體水合物富水相化學(xué)勢計算參數(shù)［13］Tab.1 Calculation properties of structure-I and structure-II hydrates［13］

Parrish W R等［14］考慮到Cjm只是溫度的函數(shù)，根據(jù)方阱勢能函數(shù)模型提出一個計算Cjm的經(jīng)驗表達(dá)式［14］:

式（10）極大地簡化了Cjm的計算過程，A、B可以通過實驗進(jìn)行擬合，Parrish W R等［14］利用Kihara勢能函數(shù)模型擬合出15種氣體的計算參數(shù)，Munck等擬合了8種氣體的計算參數(shù)。

1.2氣相組分逸度計算

SRK方程是1972年Soave對R-K方程提出的修正式，他認(rèn)為R-K方程未能準(zhǔn)確反映溫度對分子作用力的影響，并用一般化的a（T）代替了R-K方程中的a/T0.5項，并在a（T）中引入了物質(zhì)的偏心因子參數(shù)，其表達(dá)式為:

本文將運用SRK方程計算水合物相中組分逸度。上式相對溫度Tr，臨界壓力pc，臨界溫度Tc以及偏心因子等參數(shù)如表2所示。

表2　不同制冷劑的臨界參數(shù)以及偏心因子Tab.2 The critical parameters and acentric factor of different refrigerants

2　預(yù)測結(jié)果與分析

利用vdW-P模型計算制冷劑氣體水合物相平衡數(shù)據(jù)的基本流程為:

1）輸入溫度值計算該溫度下水合物的生成或分解壓力，輸入壓力初值p；

2）根據(jù)狀態(tài)方程計算氣相組分逸度fj；

3）計算Langmuir常數(shù)Cjm；

4）計算制冷劑分子在水合物晶格孔穴中的占有率θjm；

5）計算水合物相中水的參考化學(xué)勢；

6）計算富水相中（液態(tài)水或冰）水的參考化學(xué)勢；

7）對比通過步驟5和步驟6得出的計算值，判斷兩者之差的絕對值是否小于設(shè)定精度。若小于，則輸出壓力值作為該溫度下水合物的生成或分解壓力，若不滿足條件，調(diào)整壓力值返回步驟2）重新計算，直至符合條件。

對制冷劑R22，R23，R125和R143a氣體水合物的三相平衡區(qū)分別進(jìn)行模擬計算，計算結(jié)果如圖2～圖5所示，并將計算結(jié)果與文獻(xiàn)報道的實驗數(shù)據(jù)進(jìn)行對比。表3所示為文獻(xiàn)中這四種制冷劑水合物相平衡的實驗數(shù)據(jù)。

表3　文獻(xiàn)中有關(guān)水合物相平衡的實驗數(shù)據(jù)Tab.3 Experimental hydrate dissociation data from literature used in this study

正如前述模型，Langmuir常數(shù)是通過方程（10）計算所得，方程中所用到的參數(shù)是根據(jù)文獻(xiàn)實驗數(shù)據(jù)進(jìn)行優(yōu)化擬合得到，具體擬合參數(shù)如表4所示。

表4　方程（10）的擬合參數(shù)Tab.4 Fitting parameters of equation（10）

圖2所示為R22水合物相平衡圖。本文對三相平衡區(qū)R（G）-H-W（S）和R（G）-H-W（L）分別進(jìn)行模擬計算，其模擬結(jié)果與文獻(xiàn)數(shù)據(jù)相關(guān)性較高，說明該模型能夠較好地對R22三相平衡條件進(jìn)行預(yù)測。圖3所示為R23水合物相平衡圖，由圖中可以看出，模型計算值與文獻(xiàn)［7］和文獻(xiàn)［17］的實驗數(shù)據(jù)擬合程度較高，而與文獻(xiàn)［18］的數(shù)據(jù)誤差稍大，可能是由于實驗裝置不同以及裝置密封性的差異使實驗結(jié)果有偏差。

圖4和圖5所示分別為R125和R143a相平衡圖，其中R125模擬計算值與文獻(xiàn)［9］誤差較小，而與文獻(xiàn)［17］稍大，但總體變化趨勢是一致的。

圖2　R22水合物平衡相圖Fig.2 Equlibrium phase diagram of R22 hydrate

圖3　R23水合物平衡相圖Fig.3 Equlibrium phase diagram of R23 hydrate

圖4　R125水合物平衡相圖Fig.4 Equlibrium phase diagram of R125 hydrate

從圖2～圖5中可以看出，本文所給熱力學(xué)模型在計算水合物相平衡條件時具有良好的精度。模型計算值與前面所給文獻(xiàn)實驗值的平均絕對誤差如表5所示，其中最大平均誤差為2.84%，最小為1.02%。

將計算模型所得到的相平衡曲線與制冷劑氣體水合物的氣液平衡線結(jié)合，得到制冷劑氣體水合物的四相平衡點Q2，如圖2～圖5中所示，即水合反應(yīng)過程中只有當(dāng)溫度低于Q2點時水合物才有可能生成，Q2具體預(yù)測值見表6。由表6可知模型計算結(jié)果與實驗結(jié)果誤差較小，Q2點的精確預(yù)測有助于確定制冷劑氣體水合物的臨界生成條件。

圖5　R143a水合物平衡相圖Fig.5 Equlibrium phase diagram of R143a hydrate

表5　實驗值與計算值誤差Tab.5 Absolute average deviations between the calculated and experimental values

表6　制冷劑水合物四相平衡點Q2Tab.6 The quadrup le point，Q2，for four refrigerants used in this study

在基于制冷劑氣體水合物蓄冷系統(tǒng)的設(shè)計中，制冷劑氣體水合物的分解熱值也是一個重要參數(shù)。制冷劑氣體水合物的分解熱值可以通過Clausius-Clapeyron方程進(jìn)行計算:

式中:Z為氣體壓縮因子，可以通過SRK方程計算得到；R為通用氣體常數(shù)，dln p/d（1/T）可以通過三相平衡數(shù)據(jù)獲得。

圖6所示為制冷劑氣體水合物的ln p-1/T對比圖。從圖中可以看出，ln p和1/T呈現(xiàn)良好的線性關(guān)系。

圖6　制冷劑氣體水合物的ln p-（1/T）對比圖Fig.6 ln p-（1/T）contrast figure of refrigerant gas hydrates

圖7所示為制冷劑氣體水合物分解熱隨溫度變化的趨勢圖。從圖中可知，在低溫區(qū)域制冷劑氣體水合物分解熱值趨近一個恒定值。在接近275 K時，R22，R23，R125和R143a的分解熱值分別約為85 kJ/mol，83 kJ/mol，144 kJ/mol和145 kJ/mol。通過對比劉永紅等［19］計算所得制冷劑氣體水合物的分解熱值，發(fā)現(xiàn)Ⅱ型制冷劑氣體水合物的分解熱值高于Ⅰ型制冷劑氣體水合物。

圖7　制冷劑氣體水合物的分解熱Fig.7 Decomposition heat of refrigerant gas hydrates

3　結(jié)論

1）應(yīng)用經(jīng)典的vdW-P氣體水合物熱力學(xué)模型預(yù)測了R22，R23，R125和R143a水合物的相平衡分解條件，模型預(yù)測值與實驗數(shù)據(jù)平均絕對誤差分別為1.21%，2.84%，2.23%和1.02%，誤差較小，預(yù)測精度高。

2）利用計算得到的三相平衡數(shù)據(jù)繪制了R22，R23，R125和R143a水合物三相平衡相圖，通過結(jié)合相平衡圖與制冷劑氣液相平衡數(shù)據(jù)，獲得制冷劑水合物的四相平衡點Q2值，并與文獻(xiàn)中的實驗數(shù)據(jù)進(jìn)行對比，誤差較小，預(yù)測精度高。

3）計算得到制冷劑氣體水合物分解熱值，并與文獻(xiàn)結(jié)合，對比發(fā)現(xiàn)II型制冷劑氣體水合物的分解熱值高于I型。結(jié)合本文的工作，可以得到制冷劑水合物穩(wěn)定存在的條件范圍，以及通過比較分解熱的大小來選擇合適的制冷劑氣體水合物，為制冷劑氣體水合物蓄冷空調(diào)系統(tǒng)的設(shè)計提供了相關(guān)參考。

［1］ 李剛，謝應(yīng)明，劉道平.蓄冷空調(diào)新工質(zhì)-氣體水合物的替代選擇概述［J］.制冷學(xué)報，2008，29（3）:18-23. （LIGang，XIE Yingming，LIUDaoping.New type gas hydrate cool-storage media［J］.Journal of Refrigeration，2008，29（3）:18-23.）

［2］ 樊拴獅，謝應(yīng)明，郭開華，等.蓄冷空調(diào)及氣體水合物蓄冷技術(shù)［J］.化工學(xué)報，2003，54（Suppl.）:131-135. （FAN Shuanshi，XIE Yingming，GUO Kaihua，et al.Advance of gas hydrate cool storage technologies used in air conditioning system［J］.CIESC Journal，2003，54（Suppl.）:131-135.）

［3］ 方貴銀.蓄冷空調(diào)工程實用新技術(shù)［M］.北京:人民郵電出版社，2000.

［4］ Sloan E D J.Physical/chemical properties of gas hydrates and application to world margin stability and climatic change［J］.Geological Society，London，Special Publications，1998，137（1）:31-50.

［5］ Ballard A L，Sloan Jr E D.Hydrate phase diagrams for methane+ethane+propanemixtures［J］.Chemical Engineering Science，2001，56（24）:6883-6895.

［6］ Englezos P.Clathrate hydrates［J］.Industrial&Engineering Chemistry Research，1993，32（7）:1251-1274.

［7］ Kubota H，Shimizu K，Tanaka Y，et al.Thermodynamic properties of R13（CClF3），R23（CHF3），R152a （C2H4F2），and propane hydrates for desalination of sea water［J］.Journal of Chemical Engineering of Japan，1984，17（4）:423-429.

［8］ Javanmardi J，Ayatollahi S，Motealleh R，et al.Experimentalmeasurement and modeling of R22（CHClF2）hydrates in mixtures of acetone+water［J］.Journal of Chemical&Engineering Data，2004，49（4）:886-889.

［9］ Hashimoto S，Makino T，Inoue Y，et al.Three-phase equilibrium relations and hydrate dissociation enthalpies for hydrofluorocarbon hydrate systems:HFC-134a，-125，and -143a hydrates［J］.Journal of Chemical and Engineering Data，2010，55（11）:4951-4955.

［10］Liang Deqing，Guo Kaihua，Wang Ruzhu，et al.Hydrate equilibrium data of 1，1，1，2-tetrafluoroethane（HFC-134a），1，1-dichloro-1-fluoroethane（HCFC-141b）and 1，1-difluoroethane（HFC-152a）［J］.Fluid Phase Equilibria，2001，187:61-70.

［11］Eslamimanesh A，Mohammadi A H，Richon D.Thermodynamic model for predicting phase equilibria of simple clathrate hydrates of refrigerants［J］.Chemical Engineering Science，2011，66（21）:5439-5445.

［12］陳光進(jìn)，孫長宇，馬慶蘭.氣體水合物科學(xué)與技術(shù)［M］.北京:化學(xué)工業(yè)出版社，2007:41-49.

［13］John V T，Papadopoulos K D，Holder G D.A generalized model for predicting equilibrium conditions for gas hydrates ［J］.AIChE Journal，1985，31（2）:252-259.

［14］Parrish W R，Prausnitz JM.Dissociation pressures of gas hydrates formed by gas mixtures［J］.Industrial&Engineering Chemistry Process Design and Development，1972，11（1）:26-35.

［15］Wittstruck T A，Brey Jr W S，Buswell A M，et al.Solid hydrates of some halomethanes［J］.Journal of Chemical and Engineering Data，1961，6（3）:343-346.

［16］Chun M K，Yoon JH，Lee H.Clathrate phase equilibria for the water+deuterium oxide+carbon dioxide and water +deuterium oxide+chlorodifluoromethane（R22）systems ［J］.Journal of Chemical&Engineering Data，1996，41 （5）:1114-1116.

［17］Hashemi H，Babaee S，Mohammadi A H，et al.Experimentalmeasurements and thermodynamic modeling of refrigerant hydrates dissociation conditions［J］.The Journal of Chemical Thermodynamics，2015，80（11）:30-40.

［18］Heuvel M V D，Sawirjo NM，Peters C J.Influence of fluoroalkanes on the phase behaviour ofmethane gas hydrate systems［J］.Fluid Phase Equilibria，2006，241（1/2）: 124-137.

［19］劉永紅，梁得青，郭開華，等.制冷劑氣體水合物分解熱的確定［J］.西安交通大學(xué)學(xué)報，2003，37（3）:234-237.（LIU Yonghong，LIANG Deqing，GUO Kaihua，et al.The determination of decomposition heat for refrigerant gas hydrate［J］.Journal of Xi′an Jiaotong University，2003，37（3）:234-237.）

Li Na，female，Ph.D.，professor.School of Chemical Engineering and Technology，Xi′an Jiaotong University，+86 29-82665836，E-mail:lina@mail.xjtu.edu.cn.Research fields:gas hydrate of refrigerants，chemical separation processes.

Thermodynam ic M odel for Predicting Phase Equilibrium of Refrigerant Gas Hydrates

Yang Hang Li Pu Zhang Longming Li Na
（School of Chemical Engineering and Technology，Key Laboratory of Thermo-fluid Science and Engineering，Ministry of Education，Xi′an Jiaotong University，Xi′an，710049，China）

This study aims to investigate the phase equilibrium of refrigerant gas hydrates based on thermodynamic theory.A thermodynamicmodel based on the van der Waals-Platteeuw model is used to predict the hydrate dissociation conditions.Refrigerantsmodeled in this study include R22，R23，R125 and R143a.The SRK equation of state is employed formodeling the vapor and fluid phases.The deviation values betweenmodel predictions and the experimental data are1.21%，2.84%，2.23%and 1.02%，respectively.The phase equilibrium diagram of refrigeranthydrates and quadruple points are obtained.Decomposition heatof refrigeranthydrates is calculated.It is found that the decomposition heat of typeⅡrefrigerant hydrates is higher than that of typeⅠ.

gas hydrate；refrigerant；phase equilibrium；thermodynamic model

About the

TB61+1；TB61+2

A

0253-4339（2016）05-0033-06

10.3969/j.issn.0253-4339.2016.05.033

國家自然科學(xué)基金（51176154）資助項目。（The projectwas supported by the National Natural Science Foundation of China（No. 51176154）.）

2016年3月15日

簡介

李娜，女，博士，教授，西安交通大學(xué)化工學(xué)院，（029）82665836，E-mail:lina@mail.xjtu.edu.cn。研究方向:制冷劑氣體水合物，化工分離過程。