CO2置換CH4水合物的相平衡研究進(jìn)展

周琳淞，郭 平，杜建芬，霍建立

（1西南石油大學(xué)油氣藏地質(zhì)及開(kāi)發(fā)工程國(guó)家重點(diǎn)實(shí)驗(yàn)室，四川 成都 610500；2西南石油大學(xué)理學(xué)院，四川成都 610500）

天然氣水合物是在一定溫度、壓力下甲烷等氣體與水作用形成的籠形晶體化合物，廣泛分布于陸地永久凍土帶及水深大于300 m的海底，儲(chǔ)量豐富，被認(rèn)為是替代常規(guī)油氣資源的重要能源。目前開(kāi)采方法主要有熱激發(fā)法、降壓法、注入化學(xué)試劑法等，其中CO2置換開(kāi)采法具有實(shí)際可行、保護(hù)環(huán)境和資源等優(yōu)點(diǎn)，得到了各國(guó)科學(xué)家的廣泛重視。

由于CO2置換水合物中CH4的過(guò)程存在氣-液-固的相態(tài)轉(zhuǎn)換，而物質(zhì)相態(tài)轉(zhuǎn)換及達(dá)到平衡不僅受到溫度、壓力等熱力學(xué)條件的制約，而且受到氣體組成、溶液成分及鹽度、存在體系以及多孔介質(zhì)等的影響[1]，研究水合物相平衡點(diǎn)的意義在于分析影響置換過(guò)程中水合物相平衡條件的因素，得到提高CO2置換率的方法和措施，最大限度地增加天然氣等的產(chǎn)量。

1 置換機(jī)理

CO2置換CH4的過(guò)程涉及分子動(dòng)力學(xué)等微觀機(jī)理，進(jìn)而影響界面處的相平衡條件，因此不少學(xué)者對(duì)置換機(jī)理展開(kāi)了深入的研究。Schicks等[2]將合成的水合物置于CO2環(huán)境中進(jìn)行實(shí)驗(yàn)研究，微觀發(fā)現(xiàn)置換過(guò)程晶體結(jié)構(gòu)發(fā)生變化，水分子重排，水合物相與環(huán)境氣相之間的化學(xué)不平衡導(dǎo)致轉(zhuǎn)換發(fā)生。Ota Masaki等[3]研究認(rèn)為置換主要發(fā)生在水合物相表面，CH4水合物分解受水合物中水分子的重新排列控制，CO2水合物形成受水合物相擴(kuò)散控制，CH4水合物變得不穩(wěn)定并分解，擴(kuò)散入水合物相的CO2分子滲透進(jìn)入并保持水分子間氫鍵結(jié)構(gòu)不破壞，CH4可能開(kāi)始從水合物相中釋放傳遞到氣體 CO2相。近年來(lái)的研究[4-6]發(fā)現(xiàn)置換過(guò)程中水合物相不分解，置換主要在相界面進(jìn)行，依賴于氣體分子的擴(kuò)散、滲透作用，尤其是Baldwin和Jung等[7-8]在研究多孔介質(zhì)CO2(g)和CO2(l)注入條件下的置換反應(yīng)過(guò)程中，氣體分子在水合物表面的擴(kuò)散很慢，水合物保持固態(tài)，沉積物沒(méi)有硬度丟失。有學(xué)者通過(guò)分子模擬研究指出水分子間的氫鍵作用使得氣體分子從晶體籠子的進(jìn)入或溢出變得困難，進(jìn)而推測(cè)水合物相并不分解或者僅僅發(fā)生微小的溶解。Qi等[9]通過(guò)分子動(dòng)力學(xué)（MD）模擬研究 CO2兩種晶體占據(jù)（75%和 87.5%）的混合水合物的穩(wěn)定性，結(jié)果顯示平衡壓力、勢(shì)能等性質(zhì)在高占據(jù)體系中較高，說(shuō)明水合物相的分子結(jié)構(gòu)對(duì)平衡條件有著一定的影響。

2 相平衡實(shí)驗(yàn)

表1 實(shí)驗(yàn)步驟[17]

目前國(guó)內(nèi)外多采用可觀察的氣體水合物相平衡測(cè)試裝置（圖1）測(cè)定CO2置換水合物過(guò)程的相平衡條件，主要是采用等壓溫度搜索法或者恒溫壓力搜索法采集不同氣相組成的混合水合物對(duì)應(yīng)的平衡溫度和平衡壓力數(shù)據(jù)，該裝置由上海理工大學(xué)的楊光等[10]在總結(jié)水合物實(shí)驗(yàn)技術(shù)發(fā)展的基礎(chǔ)上，結(jié)合其它系統(tǒng)的優(yōu)點(diǎn)而設(shè)計(jì)，以恒溫壓力搜索法為例，實(shí)驗(yàn)步驟見(jiàn)表 1。不過(guò)該裝置無(wú)法準(zhǔn)確測(cè)試模擬多孔介質(zhì)條件下的混合氣體水合物相平衡條件，孫長(zhǎng)宇等[11]采用流量法改進(jìn)了這套裝置，將反應(yīng)釜設(shè)計(jì)為可填砂的管狀反應(yīng)器，通過(guò)控制氣體進(jìn)出口壓差來(lái)監(jiān)測(cè)流量的變化，判斷在實(shí)驗(yàn)溫度、壓力下是否有水合物生成，從而避免了多孔介質(zhì)難以觀察的問(wèn)題。近年來(lái)為了更直觀、更準(zhǔn)確地測(cè)試水合物生成的溫度、壓力、氣相組成、水合物相組成等的變化，有研究人員設(shè)計(jì)了可以連接光纖成像系統(tǒng)、色譜分析系統(tǒng)等的相平衡原位監(jiān)測(cè)實(shí)驗(yàn)裝置[12]，改善了實(shí)驗(yàn)測(cè)試的條件。在國(guó)外已經(jīng)開(kāi)始采用CT掃描等[34]成像技術(shù)觀察多孔介質(zhì)內(nèi)水合物的形成情況。

3 混合水合物相平衡的影響因素

CO2-CH4混合水合物的相平衡實(shí)驗(yàn)研究比較多，目的是獲取所研究條件的混合水合物相平衡實(shí)驗(yàn)數(shù)據(jù)，為水合物相平衡理論及模型研究提供實(shí)驗(yàn)依據(jù)，同時(shí)研究不同因素對(duì)水合物相平衡條件的影響。影響混合水合物相平衡條件的因素主要是水合物存在的體系（如電解質(zhì)溶液、有機(jī)物、表面活性劑等）、多孔介質(zhì)的孔徑分布、氣體組分、微波等介質(zhì)。

3.1 不同體系

電解質(zhì)溶液體系的水合物相形成條件研究是符合深海和凍土區(qū)沉積層實(shí)際的，表面活性劑對(duì)水合物相平衡的影響研究有利于研究促進(jìn)水合物相形成的方法。

3.1.1 純水體系

Anderson等[13]研究CH4、CO2、H2O 三元系統(tǒng)的相平衡（G-L-H）時(shí)發(fā)現(xiàn)，溫度低于283 K時(shí)CO2水合物的相平衡壓力比 CH4水合物低，說(shuō)明 CO2水合物的穩(wěn)定性更強(qiáng)。Adisasmito等[14]研究了溫度為 273～288 K、壓力為 1.2～11.2 MPa范圍內(nèi)的CO2-CH4混合水合物的H-V-Lw三相平衡條件。Seo等[15]不僅給出了混合氣體的三相平衡條件，還得到了混合氣體的較高四相點(diǎn)Q2（H-LW-LC-V），即水合物相液化點(diǎn)，四相點(diǎn)的存在必須滿足一定的氣相摩爾組成這一條件；隨著CO2氣相摩爾分?jǐn)?shù)的增加，三相平衡溫度將提高[16]；對(duì)于預(yù)測(cè)水合物的相平衡條件，可以采用van der Waals-Platteeuw模型，結(jié)合SRK-EOS和MHV2方程進(jìn)行計(jì)算。江承明等[17]在Seo等的基礎(chǔ)上測(cè)試了CO2含量較高時(shí)的混合氣體水合物的相平衡特性，指出四相共存狀態(tài)應(yīng)該為一條曲線。研究結(jié)果給出了隨CH4含量提高的水合物相平衡壓力特性，以及該混合氣體水合物體系較高四相點(diǎn)Q2穩(wěn)定區(qū)的邊界。結(jié)果表明，混合氣體中隨著CH4相對(duì)濃度的增加，Q2點(diǎn)隨之增加，四相共存狀態(tài)的范圍也隨著擴(kuò)大。畢崟等[18]研究了混合水合物四相區(qū)的分布特征，CH4組分超過(guò)臨界值時(shí)將只能形成三相狀態(tài)的混合水合物。研究四相共存狀態(tài)下的相平衡問(wèn)題對(duì)于液態(tài)CO2注入實(shí)現(xiàn)水合物置換有著重要的意義。

3.1.2 電解質(zhì)溶液體系

Dholabhai等[19]在273K附近溫度、壓力為1.5～10.6 MPa范圍內(nèi)測(cè)試了含有電解質(zhì)溶液條件下的CO2-CH4水合物的相平衡特征，在給定壓力下，平衡溫度隨著離子強(qiáng)度的增加而降低，NaCl和CaCl2具有抑制水合物形成的作用。張烈輝等[20]在?10～15 ℃溫度下混合氣體在純水及鹽類(lèi)體系下的水合物相平衡比較實(shí)驗(yàn)發(fā)現(xiàn)，鹽類(lèi)體系溶液使得水合物的生成溫度明顯下降，生成壓力明顯上升。王淑紅等[21]在研究中發(fā)現(xiàn)不同濃度的 NaCl溶液對(duì)天然氣水合物相平衡的影響也有所不同。隨著 NaCl濃度的增加，水合物生成的壓力增高，溫度降低，水合物生成區(qū)域逐漸變小。進(jìn)一步的研究認(rèn)為這是由于離子電荷的強(qiáng)電場(chǎng)作用對(duì)水分子產(chǎn)生強(qiáng)作用力，破壞了水的簇團(tuán)結(jié)構(gòu)，使得天然氣水合物分子形成晶格時(shí)需要付出一定的額外能量才能滿足水合物生成的相平衡條件，溶液濃度的增加在一定程度上加強(qiáng)了離子電荷對(duì)水分子的作用，促使水合物生成需要更多的能量補(bǔ)充，平衡條件不斷升高，但是在提到不同鹽類(lèi)溶液、不同離子對(duì)水合物形成條件的影響時(shí)，其影響程度和機(jī)理尚不清楚，目前的活度理論也不能完全解釋這種差異。

3.1.3 有機(jī)物體系

Servio等[22]在新己烷少量存在的情況下，測(cè)試了兩種CO2-CH4混合氣體在該體系中的起始水合物形成條件，新己烷的存在使得水合物的初始形成壓力有所降低，這可能與新己烷形成水合物的條件較低或者氣液界面的相互作用有關(guān)，但是需要得到進(jìn)一步的研究并證實(shí)。李遵照等[23]研究十二烷基硫酸鈉（SDS）對(duì)水合物置換過(guò)程的影響，這種表面活性劑能夠降低水的表面張力，加快水合物的形成速率。余匯軍等[24]繼續(xù)通過(guò)實(shí)驗(yàn)研究發(fā)現(xiàn)，十二烷基硫酸鈉（SDS）和十二烷基苯磺酸鈉（SDBS）的混合溶液（按1∶2的比例）能夠使CO2水合物的相平衡條件較之單一組分溶液進(jìn)一步降低，并且認(rèn)為表面活性劑降低了氣/液界面張力，提高了氣體的擴(kuò)散速率和傳質(zhì)率；陰離子活性劑吸附在帶有正電荷的金屬反應(yīng)器壁面使晶核更易形成，從而降低了水合物相平衡點(diǎn)。界面張力對(duì)水合物相平衡的影響得到了證實(shí)，曾志勇等[25]在研究多孔介質(zhì)條件下的CH4和CO2水合物相平衡預(yù)測(cè)模型時(shí)，將界面張力建立為與溫度和孔徑相關(guān)的經(jīng)驗(yàn)關(guān)系式更能準(zhǔn)確地?cái)M合實(shí)驗(yàn)所測(cè)得的平衡壓力值。多孔介質(zhì)內(nèi)氣-固-液相的接觸面非常大，因此界面張力的作用不可忽視。

3.2 多孔介質(zhì)

多孔介質(zhì)的毛管力作用不可忽視且對(duì)水合物形成的壓力、溫度條件有重要的影響，因此有必要對(duì)多孔介質(zhì)內(nèi)的混合水合物的相平衡條件展開(kāi)研究。多孔介質(zhì)對(duì)單一水合物的形成/分解過(guò)程影響的研究較為成熟。國(guó)內(nèi)外的理論模型和物理模擬研究[26-29]均表明，孔隙半徑越小，水合物形成的平衡壓力越高，平衡溫度越低。但是也有人認(rèn)為當(dāng)孔隙半徑較大或者介質(zhì)溫度較低時(shí)，孔徑分布對(duì)水合物形成的影響不大[30-32]。多孔介質(zhì)對(duì)混合氣體水合物形成的相平衡條件的影響是否同單一氣體水合物的規(guī)律一致，有學(xué)者進(jìn)行了深入的研究。Lee Huen[33]測(cè)試了模擬多孔介質(zhì)條件下的二元?dú)怏w（CH4-CO2）水合物相平衡，并同van der Waals-Platteeuw模型的計(jì)算結(jié)果進(jìn)行了比較，表明水合物形成的平衡壓力更高。多孔介質(zhì)對(duì)混合水合物相平衡的影響機(jī)理，目前主要從孔隙表面性質(zhì)、相接觸面、水合物形成等方面展開(kāi)研究。比如孔隙內(nèi)表面的潤(rùn)濕性導(dǎo)致了水的活度降低進(jìn)而使水合物形成的平衡條件提高；多孔介質(zhì)內(nèi)氣液界面的消失導(dǎo)致的有效孔隙半徑減小使得氣體需要在更高的壓力下形成水合物[34]；孔徑越小，對(duì)環(huán)境溫度的影響也越明顯，而水合物的聚集速率與溫度有關(guān)。這些研究從不同的角度解釋了多孔介質(zhì)對(duì)水合物相平衡的影響，至于這些因素的影響范圍和大小需要得到進(jìn)一步的驗(yàn)證，才能在實(shí)際水合物置換過(guò)程中更有效地避免不利因素的影響。

3.3 氣體組分

Akihiro[35]在測(cè)定 CH4與 CO2混合氣體水合物的平衡壓力和溫度時(shí)發(fā)現(xiàn)，在相同溫度下，隨著CH4含量的增加，平衡壓力顯著提高。王淑紅等[21]研究了外界條件對(duì)水合物相平衡曲線及穩(wěn)定帶厚度的影響，除氮?dú)馔獾钠渌鼩怏w都會(huì)使水合物的穩(wěn)定區(qū)域變大，且甲烷含量越少，水合物越容易形成；相同體積濃度C3H8、H2S的加入對(duì)水合物相平衡的影響程度要比C2H6、CO2明顯，表現(xiàn)為在相同溫度條件下，前者水合物相平衡時(shí)壓力降低更多，或相同壓力條件下前者水合物相平衡時(shí)溫度升高更多，使水合物穩(wěn)定帶增大更多。近期的研究[34]同樣指出，CO2的加入較之于純 CH4氣體促進(jìn)了 CH4水合物的生成，但是主要形成更穩(wěn)定的CO2水合物（實(shí)驗(yàn)過(guò)程中占水合物相的 65%以上）。氣體水合物的形成/分解條件有一定的差異，當(dāng)其它氣體的加入能夠形成更為穩(wěn)定的水合物時(shí)，相平衡條件降低，同時(shí)也說(shuō)明了氣體組分的相互作用對(duì)水合物相平衡的影響。

3.4 微波等介質(zhì)

國(guó)內(nèi)外在這方面的研究還較少。Fatykhov等[36]在研究中發(fā)現(xiàn)微波加強(qiáng)了丙烷水合物在管道的分解。梁德青等[37]實(shí)驗(yàn)研究表明，一定頻率的微波在適當(dāng)?shù)臈l件下可以降低水合物生成的過(guò)冷度和誘導(dǎo)時(shí)間促進(jìn)水合物的分解；水合物的分解產(chǎn)生的分解水與微波具有協(xié)同作用；相同壓力下，微波可提高水合物的相平衡溫度。在同樣的電場(chǎng)強(qiáng)度下，自由水的介電常數(shù)大于水合物/空籠的介電常數(shù)，使得自由水的化學(xué)勢(shì)增大更多，由固液相平衡理論吉布斯自由能的增大促使相平衡溫度提高。

4 混合水合物相平衡研究的應(yīng)用

研究氣體水合物相平衡點(diǎn)的意義在于尋找工程實(shí)際條件下能夠進(jìn)行CO2置換反應(yīng)的穩(wěn)定帶。根據(jù)深海區(qū)和永久凍土區(qū)不同的溫度、壓力變化情況，周錫堂等[38]結(jié)合現(xiàn)有研究成果繪制了兩種區(qū)域下的水合物相圖（圖2和圖3）以及水合物相存在的穩(wěn)定帶（HSZ），陰影區(qū)域Z為存在的穩(wěn)定帶。

從圖2和圖3可見(jiàn)，考慮到沉積層溫度、壓力是一定的，可能通過(guò)注熱或減壓改變沉積物中的溫度和壓力，使局部的溫度壓力條件進(jìn)入箭頭所指的區(qū)域A，該區(qū)域CH4水合物不能穩(wěn)定存在，CO2以液態(tài)形式注入實(shí)現(xiàn)置換。A區(qū)下方由CO2氣液平衡線、CO2水合物平衡線、CH4水合物平衡線所包絡(luò)的區(qū)域 CH4水合物也不能穩(wěn)定存在，但是該區(qū)域CO2呈氣態(tài)，對(duì)于工程實(shí)際注入條件有不便，并且也不利于CH4的采出分離。據(jù)相平衡實(shí)驗(yàn)研究，區(qū)域A甲烷水合物并不是完全不能存在，當(dāng)兩種氣體的比例及溫度、壓力值一定時(shí)可以形成穩(wěn)定的混合氣體水合物（圖4）[39]。Yuan等[40]在研究多孔介質(zhì)條件下的置換反應(yīng)過(guò)程中發(fā)現(xiàn)液態(tài)CO2可以在靠近A區(qū)左方由CO2氣液平衡線與CH4水合物平衡線包絡(luò)的區(qū)域?qū)崿F(xiàn)置換。

5 結(jié) 語(yǔ)

（1）CO2置換天然氣水合物中 CH4得到了廣泛的重視，已成為解決能源緊張和環(huán)境問(wèn)題的重要選擇之一。

（2）置換過(guò)程中的氣-液-固多相平衡問(wèn)題是目前國(guó)內(nèi)外研究的熱點(diǎn)，對(duì)于水合物的開(kāi)采有著重要的影響。研究階段經(jīng)歷了純水體系、電解質(zhì)溶液體系、有機(jī)物體系的混合水合物相平衡研究以及發(fā)展到目前的多孔介質(zhì)內(nèi)的綜合考慮實(shí)際儲(chǔ)層條件的水合物相平衡研究，而且在不斷探討各種因素對(duì)平衡條件的影響機(jī)理，建立了一些基于相態(tài)理論的水合物相平衡預(yù)測(cè)模型。

（3）研究多相平衡問(wèn)題的目的之一在于如何降低水合物形成的平衡條件，因此有不少學(xué)者在這方面做了大量的工作，如加入一定量的表面活性劑來(lái)降低界面張力達(dá)到降低水合物相平衡壓力的作用，研究了微波等介質(zhì)對(duì)水合物形成條件的影響等。近年來(lái)研究采用的CO2乳化液注入方式也是因其在界面性質(zhì)的優(yōu)勢(shì)得到了重視，置換效果也比較理想，但是乳化液置換過(guò)程的相平衡并不能直接用目前的氣-液-固相平衡理論，需要得到深入的研究。

（4）CO2置換法開(kāi)采水合物仍然需要結(jié)合注熱或者降壓等工程措施進(jìn)行，而且在開(kāi)采過(guò)程中水合物沉積層的結(jié)構(gòu)也在不斷變化，因此在今后的研究中必須加以考慮這些問(wèn)題，水合物多相平衡問(wèn)題不再是獨(dú)立的，應(yīng)該與滲流場(chǎng)、溫度場(chǎng)、應(yīng)力場(chǎng)等聯(lián)系在一起作為整體加以考慮。

[1]MajtnadarA，Mahmoodagh dam E，Boshinoi P R.Equilibria hydrate formation conditions for hydrogen，sulfide，carbon dioxide，and ethane in aqueous solutions of ethylene glycol and sodium chloride[J].J.Chem.Eng.，2000，45：20-22.

[2]Schicks J M，Spangenberg E，Giese R，et al.New approaches for the production of hydrocarbons from hydrate bearing sediments[J].Energies，2011，4（1）：151-172.

[3]Ota Masaki，Abe Yuki，Watanabe Masaru，et al.Methane recovery from methane hydrate using pressurized CO2[J].Fluid Phase Equilibria，2005，228-229：553-559.

[4]Bagherzadeh S Alireza，Englezos Peter，Alavi Saman，et al.Influence of hydrated silica surfaces on interfacial water in the presence of clathrate hydrate forming gases[J].J.Phys.Chem.C，2012，116：24907-24915.

[5]Jung Jong-Won，Carlos Santamarina J.Hydrate formation and growth in pores[J].Journal of Crystal Growth，2012，345（1）：61-68.

[6]Bai Dongsheng，Zhang Xianren，Chen Guangjin，et al.Replacement mechanism of methane hydrate with carbon dioxide from microsecond molecular dynamics simulations[J].Energy &Environmental Science，2012，5：7033-7041.

[7]Baldwin B A，Stevens J，Howard J J，et al.Using magnetic resonance imaging to monitor CH4hydrate formation and spontaneous conversion of CH4hydrate to CO2hydrate in porous media[J].Magn.Reson.Imaging.，2009，27（5）：720-726.

[8]Jung J W，Santamarina J C.CH4-CO2replacement in hydrate-bearing sediments：A pore-scale study[J].Geochem.Geophy.Geosy.，2010，11：Q0AA13.

[9]Qi Y X，Zhang H.MD simulation of CO2-CH4mixed hydrate on crystal structure and stability[J].Adv.Mater.Res.，2011，181-182：310-315.

[10]楊光，祁影霞，張華，等.CO2置換CH4水合物實(shí)驗(yàn)研究進(jìn)展[J].低溫與超導(dǎo)，2010，38（5）：70-75.

[11]孫長(zhǎng)宇，陳光進(jìn)，孫吉?jiǎng)?，?測(cè)定多孔介質(zhì)中氣體水合物相平衡條件的方法：中國(guó)，101477100A[P].2009-07-08.

[12]業(yè)渝光，劉昌嶺，賀行良，等.水合物相平衡原位監(jiān)測(cè)實(shí)驗(yàn)裝置：中國(guó)，202057608U[P].2011-11-30.

[13]Anderson R，Liamedo M，Tohidi B，et al.Experimental measurement of methane and carbon dioxide Clathrate hydrate equilibria in mesoporous silica[J].Phys.Chem.B，2003，107：3507-3514.

[14]Adisasmito S，F(xiàn)rank R J，Sloan E D.Hydrates of carbon dioxide and methane mixtures[J].J.Chem.Eng.Data，1991，36：68-71.

[15]Seo Y，Lee H.Hydrate phase equilibria of the carbon dioxide，methane，and water system[J].J.Chem.Eng.Data，2001，46：381-384.

[16]Seo Y，Lee H.Multiple-phase hydrate equilibria of the ternary carbon dioxide，methane，and water mixtures[J].J.Phys.Chem.B，2001，105：10084-10090.

[17]江承明，陽(yáng)濤，郭開(kāi)華.CO2-CH4混合氣體水合物相平衡實(shí)驗(yàn)研究[J].石油與天然氣工業(yè)，2010，39（6）：371-373.

[18]畢崟，陽(yáng)濤，郭開(kāi)華.二氧化碳-甲烷混合氣體水合物四相區(qū)實(shí)驗(yàn)研究[J].工程熱物理學(xué)報(bào)，2012，33（5）：739-742.

[19]Dholabhai P D，Bishnoi P R.Hydrate equilibria conditions in aqueou selectrolyte solutions：Mixtures of methane and carbon dioxide[J].J.Chem.Eng.Data，1994，39：191-194.

[20]張烈輝，李登偉，熊鈺，等.天然氣水合物相平衡的實(shí)驗(yàn)研究[J].天然氣工業(yè)，2006，26（3）：129-130.

[21]王淑紅，宋海斌，顏文.外界條件變化對(duì)天然氣水合物相平衡曲線及穩(wěn)定帶厚度的影響[J].地球物理學(xué)進(jìn)展，2005，20（3）：761-765.

[22]Servio P，Lagers F，Peters C，Englezos P.Gas hydrate phase equilibria in the system methane-carbon dioxide-neohexane and water[J].Fluid Phase Equilibria，1999，158-160：795-800.

[23]李遵照，郭緒強(qiáng)，陳光進(jìn)，等.CO2置換CH4水合物中CH4的實(shí)驗(yàn)和動(dòng)力學(xué)[J].化工學(xué)報(bào)，2007，58（5）：1197-1203.

[24]余匯軍，王樹(shù)立，江光世，等.表面活性劑對(duì) CO2水合物生成影響的實(shí)驗(yàn)研究[J].常州大學(xué)學(xué)報(bào)：自然科學(xué)版，2011，23（2）：55-59.

[25]曾志勇，李小森.基于PC-SAFT方程研究多孔介質(zhì)中水合物相平衡的預(yù)測(cè)模型[J].高等學(xué)?；瘜W(xué)學(xué)報(bào)，2011，32（4）：908-914.

[26]Smith D H，Seshadri K，Wilder J W.Assessing the thermodynamic feasibility of the conversion of methane hydrate into carbon dioxide hydrate in porous media[J].J.Energy Environ.Res.，2001，1（1）：101-117.

[27]Smith Duane H，Seshadri Kal，Uchida Tsutoma，et al.Thermodynamics of methane，propane and carbon dioxide hydrates in porous glass[J].American Institute of Chemical Engineers，2004，50（7）：1589-1598.

[28]Handa Y P，Stupin D.Thermodynamic properties and dissociation characteristics of methane and propane hydrates in 70 Aradius silica gel pores[J].J.Phys.Chem.，1992，96：8599-8603.

[29]Bower P G.Deep ocean field tests of methane hydrate formation from a remotely operated vehicle[J].Geology，1997，25：407-410.

[30]Turner D J，Cherry R S，Sloan E D.Sensitivity of methane hydrate phase equilibria the sediment pore size[J].Fluid Phase Equilibria，2005，228：505-510.

[31]Turner D J，Sloan E D.Hydrate phase equilibria measurements and predictions in sediments[C]//4th International Congress on Gas Hydrates，Yokohama，2002，5：19-23.

[32]李明川，樊栓獅，趙金洲.多孔介質(zhì)中天然氣水合物形成實(shí)驗(yàn)研究[J].天然氣工業(yè)，2006，26（5）：27-28.

[33]Lee Huen.Macro and micro approaches：Hydrate phase equilibria[C]//Proceedings of The Fifth Ocean Mining Symposium，Japan：2003.

[34]He Yuanyuan.Formation of methane and carbon dioxide hydrates in bulk and in porous media[D].Munich：Technische Universiteit Delft，2012.

[35]Akihiro Hachikubo.Phase equilibria and comparison of formation speeds of CH4and CO2hydrate below the ice point[C]//Proceedings of the Fifth International Conference on Gas Hydrates，Norway，2005.

[36]Fatykhov M A，Bagautdinov N Y.Experimental investigations of decomposition of gas hydrate in a pipe under the impact of a microwave electromagnetic field[J].High Temp.，2005，43（4）：614-619.

[37]梁德青，何松，李棟梁.微波對(duì)天然氣水合物形成/分解過(guò)程的影響[J].科學(xué)通報(bào)，2008，53（24）：3045-3050.

[38]周錫堂，俞志東.CO2置換開(kāi)采天然氣水合物相圖分析及數(shù)值模擬[J].Energy Development Frontier，2013，2（1）：40-47.

[39]Ohgaki K，Takano K，Sangawa H.Methane exploitation by carbon dioxide from gas hydrates-phase equilibria for CO2-CH4mixed hydrate system[J].J.Chem.Eng.Japan，1996，29（3）：478-483.

[40]Yuan Qing，Sun Changyu，Liu Bei，et al.Methane recovery from natural gas hydrate in porous sediment using pressurized liquid CO2[J].Energy Conversion and Management，2013，67：257-264.