梁海峰，朱耀劍，趙陽升，趙建忠

（1太原理工大學(xué)采煤工藝研究所，山西 太原 030024；2太原理工大學(xué)化學(xué)化工學(xué)院，山西 太原 030024）

?

水逸度模型預(yù)測THF添加劑體系下氣體水合物相平衡

梁海峰1，2，朱耀劍2，趙陽升1，趙建忠1

（1太原理工大學(xué)采煤工藝研究所，山西 太原 030024；2太原理工大學(xué)化學(xué)化工學(xué)院，山西 太原 030024）

摘要：加入添加劑降低水合物生成壓力是當(dāng)前水合物法分離混合氣體研究熱點(diǎn)。本研究以水的逸度模型為基礎(chǔ)，結(jié)合PRSV2狀態(tài)方程研究了CH4、O2、N2及其混合氣體水合物在純水體系下的相平衡條件；通過UNIFAC基團(tuán)貢獻(xiàn)法對添加劑四氫呋喃（THF）水溶液進(jìn)行基團(tuán)劃分，計(jì)算該體系下液相各組分的活度，理論研究了添加劑THF對氣體水合物相平衡條件的影響；結(jié)果表明在相應(yīng)的溫度范圍內(nèi)，與其他模型相比，在純水體系下該模型預(yù)測精度較高，在THF水溶液體系下該模型對單組分和雙組分氣體的預(yù)測精度平均相對誤差在7%左右，隨著THF濃度增加，氣體水合物相平衡壓力的降低幅度減?。划?dāng)THF摩爾分?jǐn)?shù)達(dá)到6%時(shí)，對氣體水合物相平衡影響達(dá)到最大。相關(guān)研究結(jié)果為混合氣體的大規(guī)模工業(yè)提純分離提供了理論基礎(chǔ)。

關(guān)鍵詞：模型；水合物；相平衡；活度系數(shù)；四氫呋喃；UNIFAC基團(tuán)貢獻(xiàn)法

項(xiàng)目。

第一作者及聯(lián)系人：梁海峰（1980—），男，博士，講師，主要從事多孔介質(zhì)內(nèi)傳熱傳質(zhì)及水合物開采研究。E-mail lianghaifeng@tyut.edu.cn。

氣體水合物是小分子氣體與水在一定的溫度、壓力條件下（通常為低溫、高壓條件）生成的非化學(xué)計(jì)量半籠型晶體，水分子構(gòu)成籠型結(jié)構(gòu)，氣體分子填充在籠型中，兩者之間通過范德華力保持穩(wěn)定。根據(jù)籠型結(jié)構(gòu)的不同，氣體水合物結(jié)構(gòu)分為Ⅰ、Ⅱ和H型。利用不同氣體之間生成水合物的溫度、壓力條件差異來實(shí)現(xiàn)混合氣體的分離是當(dāng)前氣體分離技術(shù)的研究熱點(diǎn)之一。

目前，國內(nèi)外學(xué)者對水合物法分離混合氣體進(jìn)行了大量的理論與實(shí)驗(yàn)研究，并取得了較大的進(jìn)展。HAPPEL等[1]首先在1994年利用新的分離裝置通過水合物法分離氮?dú)夂图淄?；馬昌峰等[2]等對甲烷和氫氣、二氧化碳和氫氣分別進(jìn)行了水合物法分離試驗(yàn)，結(jié)果表明采用水合物法提純氫氣效果良好；YAMAMOTO[3]、KANG[4]和孫強(qiáng)[5]等分別通過水合物法進(jìn)行了混合氣體的分離實(shí)驗(yàn)，都得到了較為理想的結(jié)果。而通過加入添加劑來降低水合物生成平衡條件也是許多學(xué)者研究的熱點(diǎn)[6-11]，ZHONG 等[12-13]分別研究了四丁基溴化銨（TBAB）和環(huán)戊烷/環(huán)己烷存在下低濃度煤層氣的水合物相平衡條件，結(jié)果表明同純水體系相比，混合氣體水合物相平衡條件向高溫低壓區(qū)域轉(zhuǎn)變，并且水合物相中包含的甲烷分?jǐn)?shù)增加為約0.5；孫登林等[14]研究了蒙脫石（MMT）對低濃度瓦斯水合分離效率的影響，加入MMT后甲烷提純率可提高30%左右；吳強(qiáng)等[15]研究發(fā)現(xiàn)，SDS-高嶺土條件下甲烷的提純率可提高至58.41%。在理論研究方面，更多學(xué)者研究通過理論模型來研究在純水體系下的水合物相平衡條件。1959年，VAN DER WAALS和PLATTEEUW[16]建立了基于統(tǒng)計(jì)熱力學(xué)的模型（vdW-P模型），成功預(yù)測了甲烷水合物的相平衡條件，隨后PARRISH[17]、ROBINSON[18]、JOHN[19]、DU[20]、高偉[21]、李小森[22]等對vdW-P模型進(jìn)行修正的基礎(chǔ)上，取得了更好的預(yù)測結(jié)果。我國學(xué)者CHEN等[23-26]提出了一個(gè)雙步驟水合物成核動力學(xué)機(jī)理模型（Chen-Guo模型），并成功預(yù)測了多種單組分及雙組分混合氣體的水合物相平衡條件。2000年，KLAUDA等[27-28]提出了基于各相中水的逸度相等的模型，對N2、CO2、CH4等單組分進(jìn)行了相平衡研究，取得了更好的預(yù)測效果。綜上，目前關(guān)于混合氣體在促進(jìn)劑條件下水合物相平衡的研究大多數(shù)僅限于實(shí)驗(yàn)，在理論方面的研究仍然側(cè)重在單組分及雙組分氣體在純水體系中的相平衡研究，尤其對于對水合物生成影響較大的添加劑四氫呋喃（THF）相平衡理論研究較少。本研究以水的逸度模型為基礎(chǔ)，結(jié)合PRSV2狀態(tài)方程研究了混合氣體水合物在純水體系下的相平衡條件；通過UNIFAC基團(tuán)貢獻(xiàn)法對添加劑THF水溶液進(jìn)行基團(tuán)劃分，計(jì)算該體系下液相各組分的活度，理論研究了添加劑THF對水合物相平衡條件的影響。

1 水逸度模型

水逸度模型，即在相平衡狀態(tài)下水合物相中水的逸度與液相或冰相中水的逸度相等，見式(1)。

式中，f為逸度，Pa；T為溫度，K；P為壓力，Pa；上角標(biāo)H、π分別代表水合物相和水或冰相，下角標(biāo)w代表水。

1.1 水合物相中水的逸度計(jì)算

水合物相中水的逸度計(jì)算公式如式(2)所示。

式中，vm為空水合物晶格中每個(gè)水分子所擁有的m型孔穴數(shù)，對于Ⅰ型水合物，大、小孔穴的值分別為2/46、6/46；對于Ⅱ型水合物，大、小孔穴的值分別為16/136、8/136；上角標(biāo)β為空水合物晶格；θml為客體分子l在m型孔穴的占有率，其表達(dá)式如式(3)所示。

式中，Cml(T)為Langmuir常數(shù)，將采用Du-Guo模型[20]來計(jì)算；下角標(biāo)l表示氣體，fwβ為空水合物的逸度，其表達(dá)式如式(4)所示。

氣相逸度采用修正的PR方程Peng-Robinson-Stryjek-Vera狀態(tài)方程[29]（PRSV2 EOS），如式(5)。

式中，R是氣體常數(shù)；T為溫度，K；V為摩爾體積，mol/m3；a和b分別表示為反應(yīng)分子間吸引力和體積的參數(shù)，見式(6)～式(9)。

式(6)～式(9)中，TC為臨界溫度，K；Pc為臨界壓力，MPa；ω為偏心因子；k1、k2和k3分別為經(jīng)驗(yàn)常數(shù)，相關(guān)數(shù)值如表1所示。

表1 PRSV2方程中參數(shù)值

對于混合氣體，采用Reid建議的混合規(guī)則，見式(10)、式(11)。

式(10)、式(11)中，下角標(biāo)i、j分別為氣體組分類別；m為混合氣體。

對于液相中添加劑存在時(shí)，本文對采用的逸度模型進(jìn)行稍微修正，即假定氣體小分子（CH4、N2、O2）同添加劑分子（四氫呋喃）形成Ⅱ型水合物，其中添加劑分子完全占據(jù)大孔，而氣體小分子只進(jìn)入小孔。相應(yīng)的水合物相水的逸度公式為式(12)～式(14)所示。

式(12)～(14)中，vsmall和vlarge分別為水合物晶格中每個(gè)水分子所占有的小孔和大孔數(shù)；θsl為氣體小分子在小孔穴的占有率；θlp為添加劑分子在大孔穴的占有率；Cl、Cp為Langmuir常數(shù)；下角標(biāo)p表示添加劑。添加劑逸度計(jì)算公式如式(15)所示。

式中，x為添加劑在液相中的摩爾分?jǐn)?shù)；γ為活度系數(shù)；下角標(biāo)p表示添加劑。

1.2 冰及液相中水的逸度

冰及液相中水的逸度表達(dá)式分別如式(16)、式(17)所示。

采用KRICHEVSKY等[30]提出的溶解度方程計(jì)算水合物存在下氣體的溶解度，簡稱K-K方程，其公式如式(18)所示。

式中，sl為氣體組分l的溶解度；HlW為氣體組分l的亨利系數(shù)；Vl∞為無限稀釋下氣體組分l的摩爾體積，mol/m3。其中無限稀釋下氣體組分l的摩爾體積和亨利系數(shù)的獲得是求解氣體溶解度的關(guān)鍵，采用Van’t Hoff方程[31]計(jì)算各種氣體的亨利系數(shù)。具體公式如式(19)所示。

表2 K-K方程中參數(shù)值

2 UNIFAC基團(tuán)貢獻(xiàn)法

UNIFAC基團(tuán)貢獻(xiàn)法是把溶液看成是由各個(gè)基團(tuán)混合而成的，即采用有限量的基團(tuán)數(shù)取代各種化合物及其組成的溶液，混合物各組分的活度系數(shù)可由組成混合物分子各基團(tuán)的基團(tuán)參數(shù)及相互作用參數(shù)得到。UNIFAC法是將活度系數(shù)lniγ表示為組合部分和剩余部分之和，見式(20)。式中，為考慮分子大小和形狀對活度系數(shù)的貢獻(xiàn)，其表達(dá)式如式(21)所示。

式中，Z為配位數(shù)，取值10；i、j代表組分；xi、xj分別為組分i、j的摩爾分?jǐn)?shù)；Φi為組分i的體積分?jǐn)?shù)；Θi為組分i的面積分?jǐn)?shù)；qi為組分i的面積參數(shù)。

剩余部分活度系數(shù)表達(dá)式如式(22)所示。

表3和表4分別給出了UNIFAC模型中部分基團(tuán)參數(shù)及基團(tuán)之間的交互作用參數(shù)[32-34]。

3 結(jié)果與討論

3.1 單組分氣體

為了驗(yàn)證模型的準(zhǔn)確性，對單組分氣體N2、O2和CH4水合物進(jìn)行了模型預(yù)測，并與實(shí)驗(yàn)值[35-36]及vdW-P模型、Chen-Guo模型預(yù)測結(jié)果[37]進(jìn)行對比，如圖1所示。從圖1中可以看出，CH4和N2水合物實(shí)驗(yàn)值與模型預(yù)測值非常接近，且吻合良好，對O2而言，溫度較低時(shí)誤差較大，而隨著溫度的升高，實(shí)驗(yàn)值與模型預(yù)測值較為吻合，這表明水逸度模型可以較好地預(yù)測單組分氣體的水合物相平衡條件。

表3 UNIFAC模型的相關(guān)參數(shù)

表4 基團(tuán)之間的交互參數(shù)

同時(shí)，圖1比較了3種模型對單組分氣體水合物相平衡條件的預(yù)測精度。從圖1中可以看出，對于CH4氣體水合物，3種模型預(yù)測精度較為接近；對于N2水合物，在溫度較低時(shí)逸度模型預(yù)測精度較差，隨著溫度升高，預(yù)測結(jié)果較為準(zhǔn)確；對于O2水合物，在相應(yīng)溫度條件下逸度模型預(yù)測精度較好。

圖1 不同模型條件下N2、O2和CH4氣體水合物生成相平衡曲線

圖2 CH4、N2和O2氣體在THF存在下水合物生成相平衡曲線

圖2 所示為單組分氣體CH4、N2和O2水合物在不同THF摩爾分?jǐn)?shù)下模型預(yù)測值與實(shí)驗(yàn)數(shù)據(jù)對比示意圖。從圖2中可以看出，模型預(yù)測值同實(shí)驗(yàn)參考數(shù)據(jù)較為吻合。同時(shí)，在添加THF條件下，各氣體水合物生成壓力大大下降，表明添加THF可以極大地降低氣體水合物生成壓力。

3.2 雙組分氣體

雙組分混合氣體水合物相平衡研究是目前許多學(xué)者研究的重點(diǎn)，本研究也就此進(jìn)行了相關(guān)研究。圖3所示為79%N2-21%O2在純水及THF摩爾分?jǐn)?shù)為5%條件下模型預(yù)測值與實(shí)驗(yàn)數(shù)據(jù)[36]對比圖。從圖3中可見，預(yù)測值與實(shí)驗(yàn)值吻合良好，且隨著溫度升高，壓力呈指數(shù)上升，而在THF條件下，在較高溫度下相平衡曲線趨緩，水合物相平衡壓力下降幅度更大。

3.3 相對誤差分析

為了更直觀地了解模型的預(yù)測精度，表5給出了各組分在不同摩爾分?jǐn)?shù)THF存在下的模型預(yù)測值同實(shí)驗(yàn)數(shù)值之間的平均相對誤差。從表5中可以看出，在相應(yīng)的溫度范圍內(nèi)，單組分和雙組分氣體水合物的預(yù)測精度平均相對誤差在7%左右，效果相對較好。

圖3 (79%N2-21%O2)+5%THF體系水合物相平衡曲線

表5 THF存在下不同組分的平均相對誤差

圖4 不同體系下低濃度煤層氣水合物生成相平衡曲線

3.4 煤層氣水合物相平衡預(yù)測

在上述研究基礎(chǔ)上，對不同摩爾分?jǐn)?shù)THF溶液條件下，低濃度煤層氣（26%CH4-60%N2-14%O2）水合物相平衡生成條件進(jìn)行預(yù)測并比較分析，以作理論參考。

圖4所示為不同THF摩爾分?jǐn)?shù)條件下低濃度煤層氣水合物相平衡曲線圖。通過圖4可以看到，無論溶液中THF存在與否，低濃度煤層氣水合物相平衡生成壓力都隨溫度升高而升高；THF的存在顯著降低了水合物的生成壓力；為了更形象直觀地看出THF的降低效果，選取溫度為284K時(shí)，THF摩爾分?jǐn)?shù)為5.4%、1.8%以及0.6%溶液同純水體系相比，相應(yīng)的煤層氣水合物相平衡壓力降低幅度分別為81.74%、76.89%、41.71%；隨著溶液中THF含量增加，低濃度煤層氣水合物相平衡壓力的降低幅度逐漸減小，這與相關(guān)學(xué)者研究[9-10]THF摩爾分?jǐn)?shù)為6%時(shí)對氣體水合物相平衡條件影響最大相一致。

4 結(jié) 論

在本研究中，以水的逸度模型為基礎(chǔ)，結(jié)合PRSV2狀態(tài)方程研究了單組分氣體N2、O2和CH4水合物在純水體系下的相平衡條件；通過UNIFAC基團(tuán)貢獻(xiàn)法對添加劑THF水溶液進(jìn)行基團(tuán)劃分，計(jì)算該體系下液相各組分的活度，理論研究了添加劑THF對水合物相平衡條件的影響，得出了以下結(jié)論。

（1）所建立的水逸度模型+PSRV2狀態(tài)方程能夠有效的預(yù)測純水體系下單組分及混合氣體的水合物相平衡條件。

（2）結(jié)合UNIFAC基團(tuán)貢獻(xiàn)法，采用上述模型能夠很好的預(yù)測促進(jìn)劑THF條件下單組分及混合氣體的水合物相平衡條件，為進(jìn)一步優(yōu)化添加劑選擇提供了理論參考。

（3）通過研究不同THF摩爾分?jǐn)?shù)對氣體水合物相平衡的影響，發(fā)現(xiàn)THF的存在可以顯著降低氣體水合物相平衡生成壓力；一定溫度下，隨著溶液中THF含量的不斷增加，低濃度煤層氣水合物相平衡生成壓力降低幅度逐漸減小。

參 考 文 獻(xiàn)

[1] HAPPEL J，HNATOW M A，MEYER H. The study of separation of nitrogen from methane by hydrate formation using a novel apparatus[J]. Annals of the New York Academy of Sciences，1994，715（1）：412-424.

[2] 馬昌峰，王峰，孫長宇，等.水合物氫氣分離技術(shù)及相關(guān)動力學(xué)研究[J].石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2002，26（2）：76-78.

[3] YAMAMOTO Y，KOMAI T，KAWAMURA T，et al. Studies on separation and purification of guest component[C]//Proc.4th International Conference on Gas Hydrates，Yokohama，Japan，2002：428-432. [4] KANG S P，LEE H. Recovery of CO2from flue gas using gas hydrate thermodynamic verification through phase equilibrium measurements[J]. Environmental & Science Technology，2000，34 （20）：4397-4400.

[5] 孫強(qiáng)，劉愛賢，郭緒強(qiáng). 水合物法分離合成氣實(shí)驗(yàn)研究[J].高?；瘜W(xué)工程學(xué)報(bào)，2010，24（5）：739-744.

[6] 呂秋楠，李小森，徐純剛，等.低濃度煤層氣分離提純的研究進(jìn)展[J].化工進(jìn)展，2013，32（6）：1267-1272.

[7] MAINUSCH S，PETERS C J，SWAAN A J，et al. Experimental determination and modeling of methane hydrate in mixture of acetone and water[J]. Journal of Chemical and Engineering Data，1997，42 （5）：948-950.

[8] HEUVEL M M，PETERS C J，ARONS J S. Influence of waterinsoluble organic components on the gas hydrate equilibrium conditions of methane[J]. Fluid Phase Equilibrium，2000，172（1）：73-91.

[9] SUN Q，GUO X Q，LIU A X，et al. Experiment on the separation of air-mixed coal bed methane in THF solution by hydrate formation[J]. Energy & Fuels，2012，26（7）：4507-4513.

[10] SUN Qiang，CHEN Guangyin，GUO Xuqiang，et al. Experiments on the continuous separation of gas mixtures via dissolution and hydrate formation in the presence of THF[J]. Fluid Phase Equilibria，2014，361（1）：250-256.

[11] FAKHARIAN H，GANJI H，NADERI F A，et al. Potato starch as methane hydrate promoter[J]. Fuel，2012，94（1）：356-360.

[12] ZHONG D L，ENGLEZOS P. Methane separation from coal mine methane gas by tetra-n-butyl ammonium bromide semiclathrate hydrate formation[J]. Energy & Fuels，2012，26（4）：2098-2106.

[13] ZHONG D L，DING K，YANG C，et al. Phase equilibria of clathrate hydrates formed with CH4+N2+O2in the presence of cyclopentane or cyclohexane[J]. Journal of Chemical and Engineering Data，2012，51 （12）：3751-3755.

[14] 孫登林，張強(qiáng)，張保勇，等.蒙脫石對低濃度瓦斯水合分離效率的影響[J].黑龍江科技學(xué)院學(xué)報(bào)，2012，22（4）：372-376.

[15] 吳強(qiáng)，朱玉梅，張保勇.低濃度瓦斯氣體水合分離過程中十二烷基硫酸鈉和高嶺土的影響[J].化工學(xué)報(bào)，2009，60（5）：1193-1198.

[16] VAN DER WAALS J A，PLATTEEUW J C. Clathrate solutions[J]. Advances in Chemical Physics，1959，2（1）：2-57.

[17] PARRISH W R，PRAUSNITZ J M. Dissociation pressures of gas hydrates formed by gas mixtures[J]. Industrial & Engineering Chemistry，1972，11（1）：27-35.

[18] NG H J，ROBINSON D B. The measurement and prediction of hydrate formation in liquid hydrocarbon-water system[J]. Industrial & Engineering Chemistry Fundamentals，1976，15（4）：293-297.

[19] 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.

[20] DU Y H，GUO T M. Prediction of hydrate formation for systems containing methanol[J]. Chemical Engineering Science，1990，45（4）：893-900.

[21] 高偉.天然氣水合物相平衡及其表面張力影響研究[D].南京：東南大學(xué)，2005.

[22] 李小森，吳慧杰，馮自平，等.利用統(tǒng)計(jì)締合流體理論狀態(tài)方程預(yù)測混合氣體水合物的平衡形成條件[J]. 化學(xué)學(xué)報(bào)，2007，65（1）：59-66.

[23] CHEN G J，GUO T M. A new approach to gas hydrate modeling[J]. Chemical Engineering Journal，1998，71（2）：145-151.

[24] 陳光進(jìn)，馬慶蘭，郭天明.氣體水合物生成機(jī)理和熱力學(xué)模型的建立[J]. 化工學(xué)報(bào)，2000，51（5）：626.

[25] SUN C Y，MA C F，CHEN G J，et al. Experimental and simulation of single equilibrium stage separation of （methane+hydrogen） mixtures via forming hydrate[J]. Fluid Phase Equilibria，2007，261（1/2）：85-91. [26] LIAO Zhixin，GUO Xuqiang，ZHAO Yanyan，et al. Experimental and modeling study on phase equilibria of semiclathrate hydrates of tetra-n-butyl ammonium bromide + CH4，CO2，N2，or gas mixtures[J]. Industrial & Engineering Chemistry Research，2013，52（51）：18440-18446.

[27] KLAUDA J B，SANDLER S I. A fugacity model for gas hydrate phase equilibria[J]. Industrial & Engineering Chemistry Research，2000，39（9）：3377-3386.

[28] MORADI G，KHOSRAVANI E. Modeling of hydrate formation conditions for CH4，C2H6，C3H8，N2，CO2and their mixtures using the PRSV2 equation of state and obtaining the Kihara potential parameters for these components[J]. Fluid Phase Equilibria，2013，338（1）：179-187.

[29] PENG D Y，ROBINSON D B. A new two-constant equation of state[J]. Industrial & Engineering Chemistry Fundamentals，1976，15（1）：58-64.

[30] KRICHEVSKY I R，KASARNOVSKY J S. Thermo dynamical calculations of solubilities of nitrogen and hydrogen in water at high pressures[J]. Journal of the American Chemical Society，1935，57 （11）：2168-2171.

[31] JOSEPH C M，RUBLN B，MOTHY R R，et al. Partial molar volumes of gases at infinite dilution in water at 298.15K[J]. Journal of Chemical and Engineering Data，1982，27（1）：22-24.

[32] 郭天民.多元汽-液平衡和精餾[M]. 北京：石油工業(yè)出版社，2002.

[33] FREDENSLUND A，JONES R L，PRAUSNITZ J M. Group-contribution estimation of activity coefficients in no ideal liquid mixtures[J]. AIChE Journal，1975，21（6）：1086-1099.

[34] POLING B E，PRAUSNITZ J M，O’CONNELL J. The properties of gases and liquid[M]. 5th ed. New York：McGraw-Hill，2001.

[35] SEO Y T，KANG S P，LEE H. Experimental determination and thermodynamic modeling of methane and nitrogen hydrates in the presence of THF，propylene oxide，1,4-dioxane and acetone[J]. Fluid Phase Equilibria，2001，189（1/2）：99-110.

[36] YANG H J，F(xiàn)AN S S，LANG X M，et al. Phase equilibria of mixed gas hydrates of oxygen+tetrahydrofuran，nitrogen+tetrahydrofuran，and air+tetrahydrofuran[J]. Journal of Chemical and Engineering Data，2011，56（11）：4152-4156.

[37] 朱耀劍. 水合物法分離低濃度煤層氣熱力學(xué)模型預(yù)測[D]. 太原：太原理工大學(xué)，2014.

研究開發(fā)

Phase equilibrium study of gas mixtures hydrate formation with additives THF based on water fugacity model

LIANG Haifeng1，2，ZHU Yaojian2，ZHAO Yangsheng1，ZHAO Jianzhong1

（1Mining Technology Institute，Taiyuan University of Technology，Taiyuan 030024，Shanxi，China；2College of Chemistry and Chemical Engineering，Taiyuan University of Technology，Taiyuan 030024，Shanxi，China）

Abstract：Gas hydrates have absorbed more and more attentions as a potential medium for gas separation. The additive promoters can reduce the hydrate equilibrium pressure at a given temperature. In this paper，we had mainly conducted on phase equilibrium of CH4、O2、N2and corresponding gas mixtures hydrate based on the water fugacity model integrated with the gas state equations of PRSV2. Furthermore，the groups of additive-aqueous solution were divided based on UNIFAC group contribution method to calculate activity of each composition in aqueous phase，and the effects of additives on phase equilibrium pressure of aforementioned gas mixtures hydrate were studied. The results showed that the water fugacity model is more accurate compared with other models in pure water system. Under the condition of tetrahydrofuran (THF) aqueous solution system，the water fugacity model prediction accuracy is about 7% for single or double-component gas hydrate. With the increase of THF concentration，the scale of the equilibrium pressure depression is reduced. When the THF aqueous solution concentration reaches 6% mole fraction，the phase equilibrium pressure of gas hydrate reaches the minimum.

Key words：model；hydrate；phase equilibria；activity coefficient；THF additives；UNIFAC group contribution method

基金項(xiàng)目：國家自然科學(xué)基金青年基金（51106104）、山西省自然科學(xué)基金青年基金（20120210022-5）及國家自然科學(xué)基金（51074111）

收稿日期：2015-09-02；修改稿日期：2015-10-23。

DOI：10.16085/j.issn.1000-6613.2016.03.007

中圖分類號：TQ 028.8

文獻(xiàn)標(biāo)志碼：A

文章編號：1000–6613（2016）03–0700–06