陳彬，辛峰，宋小飛，李鑫鋼，嚴(yán)紅超

（天津大學(xué)化工學(xué)院，天津 300072；天津化學(xué)化工協(xié)同創(chuàng)新中心，天津 300072）

相變漿液中甲烷水合物的生成過(guò)程強(qiáng)化

陳彬，辛峰，宋小飛，李鑫鋼，嚴(yán)紅超

（天津大學(xué)化工學(xué)院，天津 300072；天津化學(xué)化工協(xié)同創(chuàng)新中心，天津 300072）

利用相變材料（PCM）正十四烷的固液相變過(guò)程，吸收甲烷水合釋放的熱量，實(shí)現(xiàn)了直接換熱強(qiáng)化水合過(guò)程的目的。正十四烷與水混合制成相變?nèi)橐海≒CE），經(jīng)冷卻后形成漿液。在半間歇水合器中，測(cè)定并計(jì)算了甲烷水合物在此漿液中的收率和生成速率。為了提高計(jì)算的準(zhǔn)確性，設(shè)計(jì)了一套PVT裝置，通過(guò)減壓法實(shí)驗(yàn)測(cè)定了低溫條件下甲烷在正十四烷中的溶解度。實(shí)驗(yàn)結(jié)果表明：低溫條件下，甲烷在正十四烷中的溶解度與壓力基本呈線性關(guān)系；相比于間接傳熱方式下的水合過(guò)程，相變漿液中甲烷水合物收率及生成速率得到了有效提升。

甲烷水合物；相變漿液；溶解度；傳熱；水合物收率；水合速率

引　言

水合物，是由小分子氣體與水分子之間以非化學(xué)計(jì)量比形式形成的籠型晶體，常見的氣體有甲烷、乙烷和二氧化碳［1］。由于形成水合物的氣體分子尺寸大小不同，形成的水合物在晶體結(jié)構(gòu)上也存在一定的差異，如SⅠ型結(jié)構(gòu)（如CH4、CO2、H2等），SⅡ型結(jié)構(gòu)（如N2、乙烷等）和SH型結(jié)構(gòu)（如甲基環(huán)己烷）［2］。甲烷水合物，也稱作“可燃冰”，主要存在于深海和高緯度沉積巖層內(nèi)，根據(jù)目前已經(jīng)探明的儲(chǔ)量，以甲烷水合物形式存在的有機(jī)碳能源是當(dāng)前化石能源中儲(chǔ)存的有機(jī)碳能源的兩倍以上［3-4］。由于水合物在能源儲(chǔ)存及運(yùn)輸［5-6］、混合氣體分離［7-8］、海水脫鹽淡化［9］及二氧化碳捕獲［10］等方面的應(yīng)用前景，水合物的研究正越來(lái)越受到國(guó)內(nèi)外學(xué)者的廣泛關(guān)注。

甲烷水合物在生成過(guò)程中的水合熱效應(yīng)較為顯著，每摩爾甲烷水合物生成時(shí)釋放的熱量為（-54.44±1.45）kJ［11］，且容易發(fā)生團(tuán)聚，受限于甲烷水合物生成過(guò)程中極大的傳熱和傳質(zhì)阻力，甲烷水合物生成速率和水合物收率一般較低［12］。為了強(qiáng)化甲烷水合物的生成過(guò)程，實(shí)踐中已經(jīng)采用過(guò)多種方法，如添加表面活性劑［13］、添加活性炭和納米硅［14］、使用添加了硅砂和活性白土的多孔介質(zhì)作為水合發(fā)生場(chǎng)所［15］、采用高速攪拌［16］的水合體系等手段來(lái)強(qiáng)化水合過(guò)程中的傳質(zhì)。為了移除水合過(guò)程中產(chǎn)生的大量水合熱，Pang等［17］通過(guò)在水合發(fā)生器中內(nèi)置一組纏繞的導(dǎo)熱銅管以強(qiáng)化水合過(guò)程中的傳熱；Yang等［18］研究了在管式換熱器中水合過(guò)程的生成動(dòng)力學(xué)，然而在這樣的間接傳熱體系中，其傳熱效率依然較低。當(dāng)前已有一些利用體系內(nèi)組分相變，如冰顆粒［19-20］的融化來(lái)實(shí)現(xiàn)直接傳熱式水合過(guò)程的研究，其傳熱效率更高，但制冰過(guò)程能耗大，并且冰顆粒容易凝結(jié)，水合器易發(fā)生阻塞等問(wèn)題，不利于連續(xù)操作。近來(lái)，采用油水乳液［21］來(lái)研究氣體水合物生成過(guò)程的報(bào)道越來(lái)越多，在油水乳液中，氣體分子與水分子之間的接觸面積可以有效增大，同時(shí)甲烷等氣體在油相中溶解度遠(yuǎn)大于其在水中的溶解度，有利于氣體分子向水合物相的傳質(zhì)。在較為普遍的油包水［22］體系中，油相作為分散劑可以有效減緩水合物在形成過(guò)程中的團(tuán)聚，因而水合過(guò)程得以持續(xù)進(jìn)行，水合物收率也較純水體系有所提高。

目前采用以石蠟類材料作為相變材料［23］的水包油乳液來(lái)研究甲烷水合物生成的報(bào)道很少，本文采用了以正十四烷為相變材料的相變漿液［24］，利用正十四烷融化時(shí)的相變潛熱，可大為緩解水合過(guò)程中由水合熱造成的體系溫升，水合體系溫度可維持在正十四烷的熔點(diǎn)（277.6 K）附近；同時(shí)甲烷在正十四烷中的溶解度遠(yuǎn)大于其在水中的溶解度，利于水合過(guò)程中的甲烷傳質(zhì)；水合過(guò)程開始前，漿液中大量的固態(tài)正十四烷顆粒有利于誘導(dǎo)甲烷水合物的非均相成核，促進(jìn)水合物的生長(zhǎng)。同時(shí)，現(xiàn)有的甲烷在正十四烷中的溶解度數(shù)據(jù)中［25］，其在279～293 K條件下的溶解度數(shù)據(jù)尚不完備。為使水合過(guò)程計(jì)算更加準(zhǔn)確，本實(shí)驗(yàn)利用水合器設(shè)計(jì)了一套 PVT裝置，采用減壓法測(cè)定了低溫條件下，甲烷在正十四烷中的溶解度數(shù)據(jù)。本文還研究了相變漿液中的甲烷水合物生成過(guò)程，重點(diǎn)關(guān)注了不同壓力條件對(duì)甲烷水合物收率和生成速率的影響。

1　實(shí)驗(yàn)材料和方法

1.1材料

甲烷，純度99.99%，天津市六方工業(yè)氣體有限公司。正十四烷，純度98%，上海阿拉丁生化科技股份有限公司。脫水山梨醇單油酸酯（Span80）及脫水山梨醇單油酸酯聚氧乙烯醚（Tween60），化學(xué)純，天津市光復(fù)化學(xué)試劑有限公司。

1.2實(shí)驗(yàn)裝置

實(shí)驗(yàn)所用的裝置如圖1所示，水合器（內(nèi)徑35 mm，容積50 m l）是由316L不銹鋼加工制成。在水合器上方安置了一個(gè)轉(zhuǎn)速可達(dá)1200 r·min-1的電磁驅(qū)動(dòng)攪拌器。為使水合過(guò)程中水合器內(nèi)外溫差盡可能減小，在水合器外壁處纏繞了數(shù)圈導(dǎo)熱銅管，將低溫水浴槽中的冷卻液（乙二醇水溶液）與銅管閉合連通，通過(guò)控制冷卻液的溫度和流速可實(shí)現(xiàn)對(duì)水合器外壁溫度的控制。水合器及甲烷進(jìn)氣管路采用聚氨酯泡沫進(jìn)行封裝隔熱。采用兩個(gè)Pt100阻抗型溫度傳感器監(jiān)測(cè)水合器內(nèi)外的溫度，精度±0.1 K。通過(guò)CYB-20S型壓力傳感器（北京威斯特中航機(jī)電技術(shù)有限公司）和K9138F-1P-B型背壓閥（北京熊川閥門制造有限公司）對(duì)水合器內(nèi)的壓力進(jìn)行監(jiān)測(cè)及控制，壓力傳感器及背壓閥的精度為±10 kPa。溫度、壓力及質(zhì)量流量計(jì)（MFM）信息由數(shù)據(jù)采集系統(tǒng)每 5 s記錄一次；質(zhì)量流量控制器（MFC）信息通過(guò)D07-19BM流量積算儀由相應(yīng)軟件顯示并記錄。質(zhì)量流量計(jì)和質(zhì)量流量控制器（北京七星華創(chuàng)電子股份有限公司）精度為滿量程的2%。實(shí)驗(yàn)所用相變?nèi)橐菏怯筛咝阅躎25型分散機(jī)（德國(guó)IKA）制備得到。

1.3實(shí)驗(yàn)流程

1.3.1甲烷在正十四烷中溶解度測(cè)定使用恒流泵以5 m l·min-1的進(jìn)液速率泵入30 m l正十四烷，啟動(dòng)低溫恒溫水浴槽使裝置溫度穩(wěn)定在測(cè)量溫度，記錄此時(shí) PVT單元中的壓力（P1）和溫度（T1）。通過(guò)質(zhì)量流量控制器（MFC）通入甲烷氣體至一定壓力，關(guān)閉進(jìn)氣閥門，記錄MFC累積進(jìn)氣流量Qcul。啟動(dòng)磁力攪拌，設(shè)置轉(zhuǎn)速 250 r·min-1，使溶解達(dá)到平衡，當(dāng)PVT單元內(nèi)壓力穩(wěn)定時(shí)，關(guān)閉攪拌使裝置恒溫1 h，記錄此時(shí)的壓力（P2）和溫度（T2）。為減小誤差，在每個(gè)狀態(tài)點(diǎn)時(shí)的溶解度至少測(cè)量5次。

圖1　實(shí)驗(yàn)裝置示意圖Fig.1　Experimental apparatus

1.3.2甲烷水合物生成以質(zhì)量含量2.92%:1.08%:45%:51%的比例，精確稱取Tween60、Span 80、正十四烷和去離子水進(jìn)行混合，使用分散機(jī)在10000 r·min-1的轉(zhuǎn)速下將混合液分散5 m in使其形成均勻乳液。水合開始前，校準(zhǔn)圖1中所示的質(zhì)量流量控制器（MFC），質(zhì)量流量計(jì)（MFM）以及溫度、壓力傳感器等儀器。使用去離子水和制備的乳液相繼清洗水合器后，使用甲烷氣體多次換氣，使水合器內(nèi)處于甲烷氣氛。使用恒流泵泵入 16 m l待用乳液，啟動(dòng)低溫恒溫水浴槽，使冷卻液循環(huán)流通，待相變?nèi)橐褐械乃械恼耐橐旱文Y(jié)成固態(tài)小顆粒后，將漿液緩慢升溫至固態(tài)正十四烷顆粒的熔點(diǎn)（277.6 K）附近。將甲烷氣體快速通入水合器中，使水合器內(nèi)的壓力達(dá)到預(yù)設(shè)值，待水合器內(nèi)壓力值穩(wěn)定后，開啟攪拌，水合過(guò)程中不斷進(jìn)氣以使體系恒壓。當(dāng)MFC與MFM的瞬時(shí)流量差為零時(shí)，水合過(guò)程結(jié)束。

1.4實(shí)驗(yàn)原理

1.4.1甲烷在正十四烷中溶解度測(cè)定溶解前后，對(duì)PVT單元（水合器）中的甲烷氣體作物質(zhì)量衡算

其中nd為溶解在正十四烷中的甲烷氣體的量，mol；nin為甲烷進(jìn)氣的量，mol；ncs、ncf分別為溶解前后PVT單元?dú)庀嘀械募淄闅怏w的量，mol。

其中，20（m l）為體系內(nèi)甲烷氣體體積，22356.2 m l·mol-1是使用 P-T狀態(tài)方程［26］計(jì)算得到的甲烷氣體在標(biāo)況下的摩爾體積；Vmol（T1，P1），Vmol（T2，P2）分別為甲烷氣體在相應(yīng)狀態(tài)下的摩爾體積，可通過(guò)P-T狀態(tài)方程［26］計(jì)算。

甲烷在正十四烷中的溶解度x

其中，nn-tetra為泵入的正十四烷的量，mol。

1.4.2甲烷水合物生成甲烷水合物的生成過(guò)程是按非化學(xué)計(jì)量比進(jìn)行的相變過(guò)程［1］

其中，N為甲烷水合物水合數(shù)，理論最大值為5.75［2］。

水形成水合物后，體積約為初始體積的 1.25倍［27］。在0～t時(shí)間段內(nèi)，考慮到相變漿液的體積膨脹，進(jìn)入到相變漿液體系中的甲烷氣體的量為

其中，nMFC（t），nMFM（t）分別為0～t時(shí)間內(nèi)質(zhì)量流量控制器及質(zhì)量流量計(jì)檢測(cè)的甲烷氣體累積的量，mol；ΔnCH4（h）（t）為0～t時(shí)間內(nèi)用于形成水合物的甲烷氣體的量，mol；P、T為實(shí)驗(yàn)壓力和溫度；MW、ρW分別為水的摩爾質(zhì)量和密度，Z為實(shí)驗(yàn)條件下甲烷氣體的壓縮因子，可通過(guò)P-T狀態(tài)方程［26］計(jì)算得到。

甲烷在水中的溶解度相較于在正十四烷中的溶解度可忽略不計(jì)［28］，同時(shí)實(shí)驗(yàn)是在強(qiáng)攪拌條件下進(jìn)行，作如下假設(shè)：

（1）甲烷在正十四烷中的溶解可快速達(dá)到飽和；

（2）只考慮固態(tài)正十四烷顆粒融化和水合物形成過(guò)程中的熱效應(yīng)。

進(jìn)入到相變漿液中的甲烷氣體一部分用于參與水合，一部分溶解在正十四烷中，如式（5）所示

其中，ΔnCH4（h）（t）、ΔnCH4（d）（t）分別為形成水合物時(shí)消耗的甲烷氣體的量和溶解在液相正十四烷中的甲烷氣體的量，mol。

對(duì)水合過(guò)程進(jìn)行熱量衡算，如式（6）所示

其中，ΔHh、ΔHe分別為每摩爾甲烷水合物生成焓變和每摩爾正十四烷融化焓變。

將式（5）和式（6）代入到式（4）中，可計(jì)算得到在0～t時(shí)間內(nèi)，參與形成水合物消耗的甲烷氣體的量如式（7）所示

定義水合物收率?：在0～t時(shí)間內(nèi)形成甲烷水合物所消耗的甲烷氣體量與理論最大甲烷儲(chǔ)存量之間的比值。

其中，nW，i為相變漿液中水組分的初始量，mol。

定義水合速率rh（t）：?jiǎn)挝惑w積水中，在t～（t+Δt）時(shí)間內(nèi)參與形成甲烷水合物的甲烷氣體消耗速率。

其中，VW，i為相變漿液中水組分的初始體積。

2　實(shí)驗(yàn)結(jié)果與討論

2.1甲烷在正十四烷中溶解度測(cè)定

由于甲烷在正十四烷中溶解度較大，對(duì)相變漿液中甲烷水合物的生成過(guò)程影響較大。利用水合器，本實(shí)驗(yàn)設(shè)計(jì)了一套PVT裝置，采用減壓法測(cè)定了在臨近固態(tài)正十四烷熔點(diǎn)時(shí)的甲烷溶解度數(shù)據(jù)。

如圖2所示，在279.15 K和293.15 K條件下，溫度對(duì)于甲烷在正十四烷中的溶解度影響較小，采用279.15 K條件下甲烷在正十四烷中的溶解熱力學(xué)數(shù)據(jù)具有較強(qiáng)的可靠性。

圖2　甲烷在正十四烷中的溶解度測(cè)量Fig.2　Methane solubility measurements at 279.15 K and 293.15 K respectively

實(shí)驗(yàn)分別測(cè)定了在2～6 MPa壓力條件下的甲烷溶解度數(shù)據(jù)，壓力對(duì)甲烷溶解度的影響趨于線性，尤其是在壓力范圍較小的情況下，甲烷溶解度（x）與壓力基本呈線性相關(guān)，式（10）給出了溫度為279.15 K時(shí)甲烷溶解度與壓力的線性關(guān)系

Nourozieh等［25］測(cè)定了在溫度高于295 K條件下甲烷在正十四烷中的溶解度。與本實(shí)驗(yàn)在293.15 K溫度下測(cè)得的溶解度數(shù)據(jù)相比，壓力為2 MPa時(shí)，相對(duì)偏差為 8%，壓力為 5 MPa時(shí)，相對(duì)偏差為3.8%。因而本實(shí)驗(yàn)所測(cè)的甲烷溶解度數(shù)據(jù)可靠性較強(qiáng)。

表1　甲烷在正十四烷中的溶解度數(shù)據(jù)Table 1　Practical measurements for methane solubility in n-tetradecane

當(dāng)水合壓力為6.28 MPa時(shí)，水合過(guò)程結(jié)束后，溶解在正十四烷中的甲烷氣體量占進(jìn)入到相變漿液中的甲烷氣體總量的28.6%。結(jié)合所測(cè)溶解度數(shù)據(jù)和線性關(guān)系，實(shí)驗(yàn)采用了如表1所示的熱力學(xué)數(shù)據(jù)進(jìn)行后續(xù)水合過(guò)程的計(jì)算。

2.2甲烷水合物生成

在水合溫度277.6 K，攪拌速率800 r·min-1條件下，實(shí)驗(yàn)研究了在 5、5.25、5.5、6、6.28 MPa下的甲烷水合物生成過(guò)程。在相變漿液中，水合過(guò)程的誘導(dǎo)時(shí)間很短，主要包括恒速水合階段和快速水合階段兩個(gè)階段。

2.2.1甲烷水合物收率水合物的生長(zhǎng)是一種晶體生長(zhǎng)過(guò)程，在本實(shí)驗(yàn)使用的相變漿液中，固態(tài)正十四烷顆粒為水合物的生長(zhǎng)提供了大量成核點(diǎn)，因而水合過(guò)程可以快速開始。

如圖3所示，在各組實(shí)驗(yàn)條件下，水合過(guò)程的誘導(dǎo)時(shí)間一般在 3～10 min。當(dāng)經(jīng)過(guò)一定時(shí)間的恒速水合過(guò)程后，水/水合物/正十四烷（液相）三相接觸面積不斷增大，為進(jìn)一步水合提供了更多的生長(zhǎng)點(diǎn)，因而水合過(guò)程進(jìn)入快速水合階段。由于水合過(guò)程中，生成的水合物晶體容易團(tuán)聚，漿液黏度不斷增大，甲烷從氣相向水合物相的傳質(zhì)阻力增大，同時(shí)，固態(tài)正十四烷不斷融化，水合過(guò)程的傳熱效率也不斷降低。受限于傳質(zhì)和傳熱阻力的增大，水合過(guò)程繼而停止。

圖3　攪拌速率800 r·m in-1和初始溫度277.6 K條件下的甲烷水合物收率Fig.3　Methane hydrate yield ratio at agitation speed 800 r·m in-1and initial temperature 277.6 K

當(dāng)水合壓力在5～5.25 MPa下時(shí)，水合物收率隨壓力增大也相應(yīng)提高，最大水合物收率為92.2%（單位體積水合物可存儲(chǔ)標(biāo)況下甲烷 157體積）。當(dāng)水合壓力高于5.25 MPa時(shí)，隨著壓力的增大，水合過(guò)程的轉(zhuǎn)化率不斷降低，主要由于水合物生成速率隨著壓力的增大不斷加快，生成的水合物晶體顆粒容易團(tuán)聚，使得進(jìn)一步水合的傳質(zhì)阻力增大，因而當(dāng)水合壓力過(guò)大時(shí)，水合物收率卻會(huì)降低。

2.2.2甲烷水合物生成速率采用以正十四烷為相變材料（PCM）的相變漿液，甲烷水合物的生成速率大幅提升，水合過(guò)程耗時(shí)也相應(yīng)減少。

如圖4所示，當(dāng)甲烷水合物經(jīng)過(guò)誘導(dǎo)開始形成后，將以較為穩(wěn)定的生成速率進(jìn)行水合。隨著水合壓力的增大，水合過(guò)程的推動(dòng)力增大，恒速階段的水合速率也逐步提高。水合壓力為6.28 MPa時(shí)，恒速階段的水合速率為0.0267 mol CH4·（m l H2O·s）-1，較水合壓力 5 MPa時(shí)的恒速水合速率 0.014 mol CH4·（m l H2O·s）-1提升了90%。

圖4　攪拌速率800 r·m in-1和初始溫度277.6 K條件下的甲烷水合物生成速率Fig.4　Methane hydrate formation rate vs. time at agitation speed 800 r·m in-1and initial temperature 277.6 K

當(dāng)恒速水合進(jìn)行一段時(shí)間后，水合速率會(huì)出現(xiàn)一個(gè)陡升過(guò)程，水合轉(zhuǎn)入快速階段，隨著水合壓力的增大，水合過(guò)程將更快地進(jìn)入快速階段。當(dāng)相變漿液中的正十四烷固體顆粒完全融化后，水合過(guò)程趨于停止。

表2　快速水合階段的平均水合速率、水合過(guò)程時(shí)間及水合物收率Table 2　Average rate of methane hydrate formation during rapid formation stage and hydrate formation time as well as hydrate yield

如表2所示，水合壓力增大時(shí)，水合物生成推動(dòng)力也會(huì)相應(yīng)增大，但最終水合物收率不僅與水合推動(dòng)力有關(guān)，也與水合物生成過(guò)程時(shí)間有關(guān)。當(dāng)水合壓力為5.25 MPa時(shí)，水合過(guò)程持續(xù)時(shí)間可達(dá)50 m in，最終的水合物收率也最高。

表3　本研究中甲烷水合物生成過(guò)程與已有文獻(xiàn)數(shù)據(jù)之間的比較Table 3　Comparison of induction time， hydration time and hydrate yield between present work and previous works

在已有的關(guān)于甲烷水合物的研究中，多是采用水相溶液來(lái)研究其生成過(guò)程。如表3所示，在Ganji等［29］的實(shí)驗(yàn)中，純水體系下，甲烷水合物的生成誘導(dǎo)時(shí)間約為12 h，并且水合過(guò)程時(shí)間超過(guò)24 h。同時(shí)其研究了不同添加劑（SDS，LABS，CATB，ENP）對(duì)水合速率的影響，在最優(yōu)的實(shí)驗(yàn)條件下（8.3 MPa，500 mg·L-1SDS溶液中），生成的甲烷水合物最大收率為79.9%。Hao等［30］研究了在十二烷基硫酸鈉溶液中甲烷水合物的生成，最優(yōu)的水合時(shí)間約為370 min。從表3中的對(duì)比可以看出，采用正十四烷作為相變材料強(qiáng)化甲烷水合物生成過(guò)程中的傳熱，可大大減少水合時(shí)間，同時(shí)能得到較高的水合物收率。

3　結(jié)　論

在本實(shí)驗(yàn)工作中，通過(guò)測(cè)量低溫條件下甲烷在液態(tài)正十四烷中的溶解度，結(jié)合更為精確的計(jì)算過(guò)程，研究了在相變漿液中甲烷水合物生成過(guò)程。實(shí)驗(yàn)結(jié)果發(fā)現(xiàn)：（1）在279.15～293.15 K條件下，溫度對(duì)甲烷在正十四烷中的溶解度影響較小，在2～6 MPa壓力范圍內(nèi)，甲烷溶解度隨壓力基本呈線性變化；（2）在實(shí)驗(yàn)壓力5～6.28 MPa下，當(dāng)經(jīng)過(guò)3～10 min的誘導(dǎo)過(guò)程后，水合物開始生成，在各條件下，水合過(guò)程持續(xù)時(shí)間不超過(guò)50 m in；（3）水合壓力5.25 MPa時(shí)，可得到92.2%的最大水合物收率，較6.28 MPa下的甲烷水合物收率提高了49.6%。后期將更多地關(guān)注水合過(guò)程中的傳質(zhì)過(guò)程。

符號(hào)說(shuō)明

ΔHe——焓變，kJ·mol-1

M ——摩爾質(zhì)量，kg·mol-1

N ——水合數(shù)

n ——物質(zhì)的量，mol

P ——壓力，Pa

Q——體積流量，m l

R——Boltzmann常數(shù)

r——水合速率，mol CH4·（m l H2O·s）-1

T——溫度，K

t——時(shí)間，min

V——摩爾體積，m l·mol-1

x——摩爾分?jǐn)?shù)，%

Z——?dú)怏w壓縮因子

? ——水合物收率，%下角標(biāo)

cf——溶解后

cs——溶解前

cul——累積

d——溶解

h——水合

i——初始時(shí)

in——進(jìn)入

References

［1］ SLOAN E D， KOH C A. Clathrate Hydrates of Natural Gas［M］. 3rd ed. Boca Raton: CRC Press Inc.， 1998: 721.

［2］ STROBEL T A， HESTER K C， KOH C A， et al. Properties of the clathrates of hydrogen and developments in their applicability for hydrogen storage［J］. Chemical Physics Letters， 2009， 478（4）: 97-109.

［3］ LU S M. A global survey of gas hydrate development and reserves: specifically in the marine field［J］. Renewable and Sustainable Energy Reviews， 2015， 41: 884-900.

［4］ LEE S Y， HOLDER G D. Methane hydrates potential as a future energy source［J］. Fuel Processing Technology， 2001， 71（1）: 181-186.

［5］ SLOAN E D. Fundamental principles and applications of natural gas hydrates［J］. Nature， 2003， 426（6964）: 353-63.

［6］ CHA J H， SEOL Y. Increasing gas hydrate formation temperature for desalination of high salinity produced water w ith secondary guests［J］. ACS Sustainable Chemistry & Engineering， 2013， 1（10）: 1218-1224.

［7］ TOM ITA S， AKATSU S， OHMURA R. Experiments and thermodynamic simulations for continuous separation of CO2from CH4+ CO2gas mixture utilizing hydrate formation［J］. Applied Energy， 2015， 146: 104-110.

［8］ LEE H J， LEE J D， LINGA P， et al. Gas hydrate formation process for pre-combustion capture of carbon dioxide［J］. Energy， 2010， 35（6）: 2729-2733.

［9］ PARK K N， SANG Y H， JIN W L， et al. A new apparatus for seawater desalination by gas hydrate process and removalcharacteristics of dissolved minerals （Na+， Mg2+， Ca2+， K+， B3+）［J］. Desalination， 2011， 274（1）: 91-96.

［10］ DASHTI H， YEW L Z， LOU X. Recent advances in gas hydrate-based CO2capture［J］. Journal of Natural Gas Science and Engineering， 2015， 23: 195-207.

［11］ GUPTA A， LACHANCE J， SLOAN JR E， et al. Measurements of methane hydrate heat of dissociation using high pressure differential scanning calorimetry［J］. Chemical Engineering Science， 2008， 63（24）: 5848-5853.

［12］ RIBEIRO C P， LAGE P L C. Modelling of hydrate formation kinetics: state-of-the-art and future directions［J］. Chemical Engineering Science， 2008， 63（8）: 2007-2034.

［13］ ANDO N， KUWABARA Y， MORI Y H. Surfactant effects on hydrate formation in an unstirred gas/liquid system: an experimental study using methane and micelle-forming surfactants［J］. Chemical Engineering Science， 2012， 73（19）: 79-85.

［14］ GOVINDARAJ V， MECH D， PANDEY G， et al. Kinetics of methane hydrate formation in the presence of activated carbon and nano-silica suspensions in pure water［J］. Journal of Natural Gas Science & Engineering， 2015， 26: 810-818.

［15］ SAW V K， UDAYABHANU G， MANDAL A， et al. Methane hydrate formation and dissociation in the presence of silica sand and bentonite clay［J］. Oil & Gas Science & Technology， 2014， 70（6）: 1087-1099.

［16］ LINGA P， KUMAR R， LEE J D， et al. A new apparatus to enhance the rate of gas hydrate formation: Application to capture of carbon dioxide［J］. International Journal of Greenhouse Gas Control， 2010，4（4）: 630-637.

［17］ PANG W X， CHEN G J， DANDEKAR A， et al. Experimental study on the scale-up effect of gas storage in the form of hydrate in a quiescent reactor［J］. Chemical Engineering Science， 2007， 62（8）: 2198-2208.

［18］ YANG D， LE L A， MARTINEZ R J， et al. Kinetics of CO2hydrate formation in a continuous flow reactor［J］. Chemical Engineering Journal， 2011， 172（1）: 144-157.

［19］ STAYKOVA D K， KUHS W F， SALAMATIN A N， et al. Formation of porous gas hydrates from ice powders: diffraction experiments and multistage model［J］. The Journal of Physical Chemistry B， 2003，107（37）: 10299-10311.

［20］ HENNING R W， SCHULTZ A J， THIEU V， et al. Neutron diffraction studies of CO2clathrate hydrate: formation from deuterated ice［J］. The Journal of Physical Chemistry A， 2000， 104（21）: 5066-5071.

［21］ HUANG L， LIANG M， BO W， et al. Separation of ethylene from refinery dry gas via forming hydrate in w/o dispersion system［J］. Separation & Purification Technology， 2013， 116（37）: 342-350.

［22］ LI X G， CHEN C， CHEN Y， et al. Kinetics of methane clathrate hydrate formation in water-in-oil emulsion［J］. Energy & Fuels， 2015，29（4）: 2277-2288.

［23］ HUANG L， PETERMANN M， DOETSCH C. Evaluation of paraffin/water emulsion as a phase change slurry for cooling applications［J］. Energy， 2009， 34（9）: 1145-1155.

［24］ SONG X F， XIN F， YAN H C， et al. Intensification and kinetics of methane hydrate formation under heat removal by phase change of n-tetradecane［J］. AIChE Journal， 2015， 61（10）: 3441-3450.

［25］ NOUROZIEH H， KARIZNOVI M， ABEDI J. Vapor-liquid equilibrium measurement and thermodynamic modeling of binary systems （methane + n-tetradecane）［J］. Fluid Phase Equilibria， 2012，318: 96-101.

［26］ PATEL N C， TEJA A S. A new cubic equation of state for fluids and fluid mixtures［J］. Chemical Engineering Science， 1982， 37（3）: 463-473.

［27］ PARRISH W R， PRAUSNITZ J M. Dissociation pressures of gas hydrates formed by gas mixtures［J］. Industrial & Engineering Chemistry Process Design & Development， 1972， 11（1）: 26-35.

［28］ SERVIO P， ENGLEZOS P. Measurement of dissolved methane in water in equilibrium with its hydrate［J］. Journal of Chemical & Engineering Data， 2002， 47（1）: 87-90.

［29］ GANJI H， MANTEGHIAN M， OM IDKHAH M R， et al. Effect of different surfactants on methane hydrate formation rate， stability and storage capacity［J］. Fuel， 2007， 86（3）: 434-441.

［30］ HAO W， WANG J， FAN S， et al. Evaluation and analysis method for natural gas hydrate storage and transportation processes［J］. Energy Conversion and Management， 2008， 49（10）: 2546-2553.

Enhancement of methane hydrate formation process in phase change slurry

CHEN Bin， XIN Feng， SONG Xiaofei， LI Xingang， YAN Hongchao
（School of Chemical Engineering and Technology， Tianjin University， Tianjin 300072， China； Collaborative Innovation Center of Chemical Science and Engineering（Tianjin）， Tianjin 300072， China）

The phase change emulsion（PCE） containing phase change material（PCM） of n-tetradecane was used to directly absorb heat of the hydration by the solid-to-liquid phase change of n-tetradecane. Once n-tetradecane in PCE was solidified in the hydrator， the PCE would convert to slurry. The impacts of different pressures from 5.0 to 6.28 MPa on the yield of methane hydrate and hydration rate in the slurry at the temperature 277.6 K were researched at a sem i-batch hydrator. Considering the influence of methane dissolution in n-tetradecane on the hydration process， a PVT apparatus to measure the solubility of methane in liquid n-tetradecane was designed， and a more precise calculation process was established to describe the methane hydration. The results showed that：at low temperature conditions， the solubility of methane in liquid n-tetradecane distributed linearly versus pressure，and the methane hydrate formation process was intensified significantly compared to indirect heat removal hydration.

methane hydrate； phase change slurry； solubility； heat removal； hydrate yield； hydration rate

date: 2016-03-25.

Prof. XIN Feng， xinf@tju.edu.cn

supported by the National Basic Research Program of China （2012CB215005）.

TE 8

A

0438—1157（2016）08—3202—07

10.11949/j.issn.0438-1157.20160346

2016-03-25收到初稿，2016-06-07收到修改稿。

聯(lián)系人：辛峰。第一作者：陳彬（1991—），男，碩士研究生。

國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃項(xiàng)目（2012CB215005）。