氣體水合物分解熱的研究進展

柳秀婷+劉妮+張亞楠+由龍濤

摘要：氣體水合物作為新一代蓄冷介質，在空調蓄冷領域有良好的應用前景.其分解熱是重要的熱物性之一，對蓄冷系統的設計至關重要.介紹了國內外氣體水合物分解熱的研究進展，列出了氣體水合物分解熱的測量、計算方法以及影響氣體水合物分解熱的主要因素和影響機理，并指出氣體摩爾分數和添加劑摩爾分數均會影響水合物的分解熱，且氣體分子直徑是影響水合物分解熱的主要因素.研究為水合物的開采利用、穩定性研究以及蓄冷系統的設計提供了理論基礎.

關鍵詞：氣體水合物； 蓄冷； 分解熱

中圖分類號： TB 94 文獻標志碼： A

Abstract：As a new type of cool storage medium，gas hydrate has a good application prospect in air conditioning.The dissociation heat，one of the important thermal physical properties，is necessary to the designing of cool storage system.Recent advances in the dissociation heat of gas hydrate is presented in this paper.This work lists the calculated and measurement methods of the dissociation heat and elaborates the main factors affecting the dissociation heat and corresponding mechanism.This paper points out that the gas concentration and additive concentration have the effect on the dissociation heat.And the dominant factor that influences the dissociation heat of the hydrate is gas molecular diameter.The study on the dissociation heat provides the theoretical basis for the exploitation，the hydrate stability study，and the storage system design.

Keywords：gas hydrate； cool storage； dissociation heat

隨著社會和經濟的持續發展，空調使用越來越普及，工業和居民用電量增大，導致我國電網負荷和“峰谷”差越來越大.但是由于我國能源分布不均，能源利用效率普遍偏低，造成了能源供給的緊張.空調蓄冷技術能充分利用能量，節約費用.它利用白天與夜間的電價差，在夜間儲存冷量，白天釋放冷量制冷，有效平衡能源系統.

從外部而來的客體分子與水分子所形成的一種類似于包絡狀的結晶體被統稱為氣體水合物，它可以在冰點以上的條件

下進行結晶固化.制冷劑氣體水合物由常規制冷劑和水在5～12℃相互作用生成，相變溫度高，相變潛熱大，儲-釋冷時的傳遞效率高，克服了水、冰和共晶鹽等蓄冷介質的不足，被認為是一種比較有應用前景的空調蓄冷介質.

分解熱是選擇蓄冷材料時的一個重要依據.分解熱越高，蓄冷劑的用量就越小，相應地可以減小蓄冷槽的容積，從而節省初投資和占地面積.此外，水合物的形成會堵塞油氣輸送系統，需要加熱分解[1]；每150～180個體積的天然氣溶于每1個體積水中，可以儲存和運輸天然氣[2]；水合物的各種熱性質直接影響它的穩定性，存在于地下的各種氣體水合物的開采和一系列利用也需要利用水合物分解熱數據進行分析[3].因此，測定水合物的分解熱對于蓄冷材料的選擇、蓄冷系統的設計以及水合物的開采利用和穩定性研究都非常必要.

1 氣體水合物分解熱的研究現狀

氣體水合物作為蓄冷介質應用于空調領域中具有很多優點.美國、日本率先對其展開研究，我國自20世紀80年代開始了這一研究工作.Gawron等[4]研究了相變溫度在-70～0℃內的許多單元相變蓄冷材料的熱物性；Krichel[5-6]繪制了大量的相變材料圖表.20世紀70年代末，美國國家實驗室開始利用差示掃描量熱法（DSC）測試了多種相變材料的性能[7-8].Handa[9]和Lievois等[10] 都曾使用過量熱法對純組分和一些混合物的水合物分解熱進行一系列測量分析，如甲烷、乙烷和丙烷氣體水合物分解熱.Rueff等[11]在平均溫度為285 K時運用DSC對甲烷水合物的分解熱進行了測量，最終結果為54.49 kJ·mol-1.文獻[12-13]中也對CFC11水合物和HCFC141b水合物的分解熱進行了測試.利用DSC測量了低壓制冷劑HCFC141b氣體水合物的分解熱，結果表明，在相變溫度約為6.5℃時，其分解熱約為330 kJ·kg-1，和冰的潛熱基本相同.

2 氣體水合物分解熱的測量方法

水合物分解熱包含兩部分：水由水合物相轉變為液相所引起的焓差與氣體水合物相轉變為氣相引起的焓差.分解熱所包含的測量方法一般分為一般卡計法、差熱分析法和差示掃描量熱法.

2.1 一般卡計法

連續加熱絕對法和投下冷卻法為最常用的卡計法.當物質不發生過冷現象，并且融解熱和凝結熱保持一致時，可以使用投下冷卻法進行測量.該方法需要測量四個溫度，先測試融點Tm，然后分別對試樣在低于融點的溫度T1和高于融點的溫度T2、T3三點做三次測試，由每次的焓值計算出融解熱.該方法原理簡單但測試過程繁瑣，不適用于多相變點和有過冷現象的物質.在整個升溫過程中均使用定常功率加熱試樣的方法稱為連續加熱絕對法 [14].使用該方法時卡計的熱保護屏可以為試樣系統提供一個絕熱環境.忽略試樣的熱物性變化，在加熱一段時間后，試樣溫度在空間的分布相對不變，試樣各點的溫度隨時間的溫升速率均相等.在此過程中，布置在試樣任意點的熱電偶均能測到試樣的平均溫升.

2.2 差熱分析法

當溫度在一定程序控制下以恒定速度變化時，對作為參比物的某種穩定物質與等量試樣進行連續測量，得到溫度差與溫度之間的熱譜圖，并且通過分析從而得到所需要的數據被稱為差熱分析法（DTA）.在選定參比物時，要求參比物在一定試驗溫度下不會發生任何化學反應和物理變化.

DTA曲線只適用于定性或半定量分析，主要原因是在產生熱效應時試樣與參比物，以及周圍環境的溫度有很大差異，熱效應測量的靈敏度和精準度因三者之間發生了熱交換均有所降低，并且試樣升溫速率的非線性使得校正系數K不斷隨溫度變化，從而導致無法進行定量分析.

2.3 差示掃描量熱法

差示掃描量熱法是另外一種熱分析法.在程序控制溫度下，測量輸入至試樣與參比物內的熱功率差和時間或溫度之間的關系.DSC具有樣品使用量少、測試速度快和定量分析精確的優點.利用DSC測試水合物分解熱的研究在國內外都還很少，并且只進行了幾種水合物的測試，主要有氙、四氫呋喃、環氧己烷、甲烷和丙烷等.

與DTA曲線不同的是，K在DSC曲線中與溫度無關，不會隨著溫度變化而變化，因此DSC曲線可用于定量分析.然而，在使用DTA和DSC測量時，由于水合數均為假定，因此測量結果的準確性有待改進.

上述兩種方法是測量水合物分解熱最常用的方法.國內外已有學者利用這兩種方法對不同水合物的分解熱進行了測量分析.Fournaison等[15]利用DTA測得了CO2水合物的分解熱為507 kJ·kg-1，比冰的分解熱（333 kJ·kg-1）大.Lin等[16]在壓力為2 MPa下使用DTA測得質量濃度為9.01%的四丁基溴化銨（TBAB）水合物的分解熱為313.2 kJ·kg-1，比CO2水合物的分解熱小，但是與冰的分解熱接近.孫志高等[17]利用德國Netzsch的204HP/1/G Phoenix型DSC熱分析儀，采用熱循環法測量了TBAB水溶液的分解熱，測量結果為195.9 kJ·kg-1.Marinhas等[18]也通過熱循環法利用111type setaram 型高壓DSC裝置測得CO2水合物漿的分解熱為500 kJ·kg-1，大約是冰漿分解熱的1.5倍.此外，Sari等[19]通過實驗證實了在壓力為3 MPa時，水合物生成溫度為1～2℃，溶解溫度為8～10℃，并且利用高壓DSC設備測得在固相質量分數10.8%時，CO2水合物的分解熱為54 kJ·kg-1.此結果與Marinhas等[18]測得的結果一致.

3 氣體水合物分解熱的計算方法

水合物分解熱不僅不易直接測量，而且也不可能直接測量所有條件下水合物的分解熱，其主要原因是一些氣體水合物的形成壓力較高，以及形成的水合物純度不易確定等.因此，可以利用測得的水合物相平衡數據，運用克勞修斯克拉貝龍方程計算0℃以上水合物的分解熱.這為確定0℃以上氣體水合物的分解熱提供了一個簡單可行的方法.計算方程可描述為

式中：p、T分別為水合物形成/分解的壓力和溫度；ΔH為分解熱；Z為氣體壓縮因子；R為通用氣體常數.

當壓縮因子Z變化很小時，在很小的溫度范圍內都可以得到焓值.

然而，在使用克勞修斯克拉貝龍方程計算分解熱時，得出的理論計算值與實驗值往往會出現偏差.例如，理論上Ⅱ型水合物在形成時水分子與外來分子的分子數之比為17∶1，而在水合物的實際生成過程中，由于外來分子的晶格占有率可能稍有改變，從而導致理論計算值出現偏差.

為了得到更為準確的數據，在使用式（1）計算水合物的分解熱時應同時考慮相體積變化、非理想氣體、氣體溶解度和水飽和度等因素.然而，在一定的限制條件下，式（1）可以應用在單變系統中（水+氣體水合物）[1，20-21].

除了以上所介紹的方法外，還有很多測量水合物分解熱的方法，只是不同方法得到的結果存在差異.表1給出了CO2水合物的水合數和分解熱[22].但是目前還無法確定哪一種方法最精確.

4 影響水合物分解熱的因素

4.1 氣體分子摩爾分數對水合物漿分解熱的影響

形成水合物漿的氣體分子摩爾分數對漿體的相變潛熱有較大影響，最佳的氣體分子摩爾分數可以使水合物漿的表觀黏度最小，相變潛熱最大.CO2水合物漿的分解熱與CO2摩爾分數有關[22-23]，并且存在最佳摩爾分數.在四氫呋喃（THF）摩爾分數為5%，CO2摩爾分數分別為2%和3%時，發現CO2摩爾分數對CO2+THF水合物漿的分解熱影響不大；當CO2摩爾分數增加至9%時分解熱迅速增加，尤其是在低溫區域；然而，當CO2摩爾分數增加到19%以及29%時，水合物漿的分解熱明顯減小.圖1給出了不同THF摩爾分數下CO2+THF水合物漿的分解熱.

4.2 添加劑對水合物分解熱的影響

目前添加劑作為水合物生成促進劑已被廣泛應用在水合物生成技術中，最佳的添加劑濃度應既可促進水合物的形成，又能促使水合物的分解[23].在生成CO2水合物時加入THF不僅可以降低水合物的生成壓力，而且可以增加水合物的分解熱[24].當溫度在282.06～273.15 K內變化時，CO2水合物的分解熱為56.85～75.37 kJ·kg-1，CO2+THF水合物的分解熱為112.37～152.27 kJ·kg-1.水合物分解熱增加了兩倍，其主要原因是水合物的結構發生了改變（從Ⅰ型轉變為Ⅱ型）.并且，其分解熱與THF摩爾分數有關（如圖1所示），在摩爾分數為7%時得到最大值，大于7%時由于降低了水合物的穩定性反而對分解熱有抑制作用.

Delahaye等[25]也將THF作為添加劑加入CO2中形成CO2+THF水合物，發現CO2+THF水合物的分解熱約為CO2水合物分解熱的兩倍，此結果與Sabil等[23]的結果一致.

4.3 氣體分子直徑對水合物分解熱的影響

水合晶體的結構取決于氣體分子與晶穴內徑的相對尺寸大小，并且外來氣體分子的大小對水合物的形成起著至關重要的作用.劉永紅等[26]根據測定的制冷劑氣體水合物的形成條件和壓縮因子，運用克勞修斯克拉貝龍方程計算并分析比較了HFC134a和HCFC141b水合物在0℃以上的分解熱，結果表明水合物分解熱與其結構有關.雖然這兩種水合物的氣體組分不同，但是分子大小相差無幾，其分解熱也比較接近.孫志高等[27]使用克勞修斯克拉貝龍方程分別計算了三種不同水合物的分解熱，結果表明氣體分子對水合物的分解熱有很大影響.

5 結 論

本文指出了氣體水合物分解熱的研究進展、測量與計算方法及影響因素，指出了氣體水合物在空調蓄冷領域的廣闊前景.綜述了國內外的研究成果后發現，形成水合物的氣體摩爾分數和添加劑摩爾分數對水合物的分解熱均有一定的影響，但由于實驗過程中摩爾分數范圍寬，為得到最佳濃度，還需要進一步分析研究.影響水合物分解熱的主要因素并不是氣體水合物的氣體組分，而是形成水合物的氣體分子直徑.此外，水合物分解熱還受到其他因素的影響，例如壓縮因子、水合物的含水量等.在今后的研究工作中，應綜合考慮各方面的影響因素，為氣體水合物的穩定性研究以及蓄冷系統的設計等工作提供精確的基礎數據.

參考文獻：

[1] SLOAN E D，FLEYFEL F.Hydrate dissociation enthalpy and guest size[J].Fluid Phase Equilibria，1992，76：123-140.

[2] KHOKHAR A K.Storage properties of natural gas hydrates[D].Trondheim：Norwegian University of Science and Technology，1998.

[3] SLOAN E D J.Clathrate hydrates of natural gases[M].2nd ed.New York：Marcel Dekker Inc.，1998.

[4] GAWRON K，SCHROEDER J.Properties of some salt hydrates for latent heat storage[J].International Journal of Energy Research，1977，1（4）：351-363.

[5] KRICHEL K.Latentwarmespeicher：Teil 1，Eigenshaften und anwendungs moglichkeien von latentwarmespeichern：Report BMVGFBWT，79-12[R].Berghausen，FRG：Institute fuer Chemie der Treib Explosivstoffe，1979.

[6] KRICHEL K.Latentwarmespeicher：Teil 2，Speichermassen bayformen und werkstoffe fuer schmelzwarmespeich：Report BMVGFBWT，79-13[R].Berghausen，FRG：Institute fuer Chemie der TreibundExplosivstoffe，1979.

[7] CANTOR S.Application of DSC （differential scanning calorimetry） to the study of thermal energy storage[J].Thermochimica Acta，1978，26（1/2/3）：39-47.

[8] LANE G A.Solar heat storage：Latent heat materialVolume II：Technology[M].Florida：CRC Press Inc.，1986：50-65.

[9] HANDA Y P.Compositions，enthalpies of dissociation，and heat capacities in the range of 85 to 270 K for clathrate hydrates of CH4，C2H6，C3H8[J].The Journal of Chemical Thermodynamics，1986，18（10）：915-921.

[10] LIEVOIS J S，PERKINS R，MARTIN R J，et al.Development of an automated，high pressure heat flux calorimeter and its application to measure the heat of dissociation and hydrate numbers of methane hydrate[J].Fluid Phase Equilibria，1990，59（1）：73-97.

[11] RUEFF R M，SOLAN E D，YESAVAGE V F.Heat capacity and heat of dissociation of methane hydrates[J].AIChE Journal，1988，34（9）：1468-1476.

[12] 梁德青.新型蓄冷介質制冷劑氣體水合物相平衡熱力學研究[D].上海：上海交通大學，2001.

[13] 梁德青，郭開華，樊栓獅，等.HCHC141B氣體水合物融解熱的DSC測試[J].工程熱物理學報，2002，23（增刊）：47-49.

[14] 陳則韶，葛新石，顧沁.量熱技術和熱物性測定[M].合肥：中國科學技術大學出版社，1990.

[15] FOURNAISON L，DELAHAYE A，CHATTI I.CO2 hydrates in refrigeration processes[J].Industrial & Engineering Chemistry Research，2004，43（20）：6521-6526.

[16] LIN W，DELAHAYE A，FOURNAISON L，et al.Thermal properties of CO2 hydrate in the presence of TBAB[C].Beijing：International Congress of Refrigeration，2007.

[17] 孫志高，江承明，孫立.新型空調蓄冷材料四丁基溴化銨-水相變條件和相變熱實驗研究[J].制冷學報，2009，30（5）：24-26.

[18] MARINHAS S，DELAHAYE A，FOURNAISON L，et al.Modeling of the available latent heat of a CO2 hydrate slurry in an experimental loop applied to secondary refrigeration[J].Chemical Engineering and Processing，2006，45（3）：184-192.

[19] SARI O，HU J，BRUN F，et al.Insitu study of the thermal properties of hydrate slurry by high pressure DSC[C].Beijing：International Congress of Refrigeration，2007.

[20] GLEW D N.Some stoicheiometric gas hydrate[J].Nature，1959，184（4685）：545-546.

[21] BARRER R M.Validity of clapeyrons equation for phase equilibria involving clathrates[J].Nature，1959，183（4659）：462-463.

[22] HU J，SARI O，HOMSY P.Investigation of thermophysical and flow properties of CO2 hydrate slurry[J].International Journal of AirConditioning and Refrigeration，2011，19（3）：213-229.

[23] SABIL K M，WITKAMPA GJ，PETERSA C J.Estimations of enthalpies of dissociation of simple and mixed carbon dioxide hydrates from phase equilibrium data[J].Fluid Phase Equilibria，2010，290（1/2）：109-114.

[24] 劉莉，王世平，巫江虹，等.R142b氣體水合物分解放冷過程的實驗研究[J].制冷學報，1999，14（2）：1-5.

[25] DELAHAYE A，FOURNAISON L，MARINHAS S，et al.Effect of THF on equilibrium pressure and dissociation enthalpy of CO2 hydrates applied to secondary refrigeration[J].Industrial & Engineering Chemistry Research，2006，45（1）：391-397.

[26] 劉永紅，梁德青，郭開華，等.制冷劑氣體水合物分解熱的確定[J].西安交通大學學報，2003，37（3）：234-236.

[27] 孫志高，樊栓獅，郭開華，等.天然氣水合物分解熱的確定[J].分析測試學報，2002，21（3）：7-9.