天然氣水合物顆粒微觀受力及聚集特性研究進(jìn)展

劉海紅 ，李玉星 ，王武昌 ，陳 鵬 ，張慶東，樊新斌

（1 中國石油大學(xué)（華東）儲運(yùn)與建筑工程學(xué)院，山東 青島 266580；2 中國石油長慶油田分公司第八采油廠，陜西 西安 710015）

自1934年天然氣水合物首次被發(fā)現(xiàn)，氣體水合物逐漸被熟知并成為研究的熱點(diǎn)，尤其是近年來深海油氣田的開發(fā)，常用的熱力學(xué)抑制劑表現(xiàn)出了環(huán)境污染、用量大等問題，伴隨著動力學(xué)抑制劑和防聚劑的開發(fā)，水合物冷流技術(shù)及其對應(yīng)的風(fēng)險(xiǎn)管理策略已成為水合物防治的研究重心[1]，其核心在于不抑制水合物生成的前提下保證水合物顆粒分散在油相中安全輸送，這不僅保證了油氣管道的安全運(yùn)行，還能充分利用水合物的高儲氣性能達(dá)到了輸送天然氣的目的。

目前國內(nèi)外許多科研機(jī)構(gòu)如中國科學(xué)院廣州能源研究所、中國石油大學(xué)（華東）、中國石油大學(xué)（北京）、美國克羅拉多礦業(yè)學(xué)院、法國石油研究所等都致力于水合物冷流技術(shù)的研究，文獻(xiàn)調(diào)研發(fā)現(xiàn)當(dāng)今的研究大都是基于反應(yīng)釜和環(huán)道的水合物漿宏觀流動實(shí)驗(yàn)，觀測壓力、溫度、流速、含水量、顆粒粒徑等參數(shù)對水合物漿流變性的影響，對水合物微觀受力和聚集方面的研究較少，這主要是因?yàn)椋孩賹?shí)驗(yàn)條件局限，實(shí)驗(yàn)設(shè)備昂貴；②水合物聚集過程復(fù)雜，研究成果較少，沒有充足的理論基礎(chǔ)。但水合物顆粒受力及聚集特性是決定水合物漿宏觀流動的內(nèi)在原因，是保證安全流動的根本，下文就當(dāng)前國內(nèi)外關(guān)于水合物顆粒受力及聚集特性研究成果進(jìn)行總結(jié)。

1 水合物微觀測量裝置

為了研究水合物顆粒受力及聚集過程，美國克羅拉多礦業(yè)學(xué)院、法國石油研究所等研究機(jī)構(gòu)都引進(jìn)了一些微觀測量設(shè)備[2]，主要有聚焦光束反射測量儀（FBRM）、顆粒圖像顯微鏡（PVM）、高壓差示掃描量熱儀（HP-DCS）和微機(jī)械測力裝置（MMF），這些設(shè)備為水合物顆粒微觀受力及聚集過程的研究提供了有效的測試手段。

1.1 聚焦光束反射測量儀和顆粒圖像顯微鏡

聚焦光束反射測量儀和顆粒圖像顯微鏡是目前研究水合物顆粒和聚集體尺寸的常用設(shè)備，使用該設(shè)備研究水合物的生長、聚集過程，主要是通過測量水合物顆粒在反應(yīng)釜或環(huán)道實(shí)驗(yàn)中的粒徑變化，根據(jù)粒徑變化來判斷水合物聚集程度。圖1為研究反應(yīng)釜中水合物生成、聚集的一套常用設(shè)備，透明的窗口用于測試人員觀測和拍攝，F(xiàn)BRM和PVM 用于微觀粒徑測量。

圖1 高壓反應(yīng)釜

1.2 高壓差示掃描量熱儀

商用高壓差示掃描量熱儀操作壓力可達(dá)40 bar（1bar=105Pa），工作溫度范圍為-45～120℃，如圖2。根據(jù)水合物體系分解過程中的分解熱來判斷水合物的聚集程度，主要原理是水合物分解過程的分解熱與水合物/油的接觸面積有關(guān)，接觸面積越大分解熱就越多，說明水合物聚集程度越嚴(yán)重，基于這個(gè)特性判斷水合物的聚集程度。

1.3 微機(jī)械測力裝置

微機(jī)械測力裝置由數(shù)字視頻顯微鏡、纖維操縱器、冷卻單元、玻璃纖維懸臂和數(shù)據(jù)記錄處理系統(tǒng)組成，如圖3。該裝置可以測量水合物顆粒間黏附力，測量原理是彈簧胡克定律（圖4），主要操作步驟為：①滴狀液體置于玻璃纖維懸梁臂終端，將液滴置于液氮中生成球形水合物固體顆粒；②兩個(gè)水合物顆粒接觸-分離，用數(shù)字視頻顯微鏡追蹤分離過程中低彈簧常數(shù)懸梁臂終端顆粒的位移δ；③根據(jù)測得位移δ和懸梁臂彈性常數(shù)k，結(jié)合胡克定律F=kδ就能求得顆粒間黏附力；④每組實(shí)驗(yàn)進(jìn)行40次上面的接觸-分離實(shí)驗(yàn)，實(shí)驗(yàn)結(jié)果取平均值。

圖2 高壓差分掃描熱量儀

圖3 微機(jī)械實(shí)驗(yàn)裝置

圖4 微機(jī)械測力裝置原理圖

2 水合物顆粒微觀受力

水合物顆粒在流動體系中的微觀受力特性是研究水合物漿技術(shù)的基礎(chǔ)，水合物顆粒受力決定水合物顆粒的聚集過程，明確水合物顆粒間微觀相互作用才能從本質(zhì)上實(shí)現(xiàn)對水合物聚集過程的宏觀控制，保證水合物漿安全輸送。

2.1 水合物顆粒間的主要黏附力

許多學(xué)者都認(rèn)為水合物顆粒間的主要黏附力是毛細(xì)液橋力，為了證實(shí)這一點(diǎn)進(jìn)行了大量的實(shí)驗(yàn)和理論研究。

水合物顆粒間液橋如圖5，毛細(xì)液橋力計(jì)算公式如式（1）。

式中，γLL為油水界面張力；θ為固液接觸角，體系中無活性劑作用θ=0°；H為液橋長度；d為液橋?qū)λ衔镱w粒的浸沒高度；R*為兩個(gè)水合物顆粒的調(diào)和半徑。

圖5 水合物顆粒間液橋圖

Yang等[3]利用微機(jī)械測力裝置測量了四氫呋喃水合物在正癸烷中的黏附力，實(shí)驗(yàn)結(jié)果表明，顆粒間的黏附力測量值與毛細(xì)液橋力的計(jì)算數(shù)值吻合，可以用毛細(xì)液橋力來解釋四氫呋喃水合物顆粒間的黏附力。Taylor等[4]在Yang 實(shí)驗(yàn)的基礎(chǔ)上改進(jìn)實(shí)驗(yàn)裝置，測量了四氫呋喃水合物在正癸烷中的黏附力，向體系中加入防聚劑Span 20和動力學(xué)抑制劑PVCap，兩種物質(zhì)改變了顆粒表面特性，使得油水相間界面能降低，導(dǎo)致四氫呋喃水合物顆粒間黏附力降低。Dieker等[5]同樣利用微機(jī)械測力裝置測量了環(huán)戊烷水合物在添加少量原油情況下的黏附力，發(fā)現(xiàn)原油的存在大大降低了水合物顆粒間的黏附力，這是因?yàn)樵椭械氖退岷蜑r青質(zhì)降低了油水間界面能，改變了顆粒表面潤濕性，接觸角增大，從而使得水合物間黏附力降低。Aman等[6]對Taylor和Dieker的實(shí)驗(yàn)進(jìn)行了補(bǔ)充，利用微機(jī)械測力裝置和界面張力儀（IFT）測量了環(huán)戊烷水合物在多種添加劑下的黏附力和油水界面張力，指出不同添加劑的作用不同，可能同時(shí)改變油水界面能和固液接觸角，也可能是改變其中之一。Taylor、Dieker和Aman的實(shí)驗(yàn)結(jié)果都證明水合物顆粒間主要黏附力可以用毛細(xì)液橋力解釋，數(shù)值也吻合。

不僅實(shí)驗(yàn)如此，理論推導(dǎo)得到了相同的結(jié)論，Anklam等[7]為了探究水合物的聚集行為，對水合物顆粒的受力進(jìn)行了研究，詳盡地分析了水合物顆粒的立體勢力、范德華力、毛細(xì)液橋力和流體剪切力，給出了各力的計(jì)算公式，并結(jié)合實(shí)際運(yùn)行工況的物理?xiàng)l件進(jìn)行理論計(jì)算，定量分析了各力在水合物聚集過程中所起的作用，理論計(jì)算結(jié)果表明毛細(xì)液橋力是導(dǎo)致水合物顆粒聚集的主要黏附力。

McCulfor等[8]研究了毛細(xì)液橋力對分散在煤油中的玻璃顆粒流動特性的影響，實(shí)驗(yàn)結(jié)果表明毛細(xì)液橋力使得體系黏度急劇增大，但是表面活性Span 80和Arquad 2HT 可以降低顆粒間毛細(xì)液橋力，使得體系黏度降低，毛細(xì)液橋力對懸浮液流動性產(chǎn)生了重大影響，導(dǎo)致堵塞事故發(fā)生。

2.2 水合物受力平衡模型

Camargo等[9]基于水合物聚集體受黏附力和剪切力平衡，建立了聚集體臨界最大粒徑預(yù)測模型，該模型在已知初始水滴粒徑、連續(xù)相黏度、水合物體積分?jǐn)?shù)、剪切速率的條件下結(jié)合實(shí)驗(yàn)和理論確定的主要黏附力-毛細(xì)液橋力就可以預(yù)測流動中水合物的聚集尺寸。OLGA 軟件的CMSHYK 模塊[10-15]就是利用該模型計(jì)算水合物在油相中聚集體的粒徑，見式（2）。

式中，dA實(shí)際上是水合物顆粒最大臨界聚集粒徑dA,max，初始水合物顆粒聚集體較小，此時(shí)dA＜dA,max，水合物顆粒聚集體受到的聚集力大于分離力，聚集會繼續(xù)進(jìn)行，直到水合物顆粒聚集體dA接近或者等于dA,max，聚集過程才會停止，如果在這個(gè)過程中水合物漿能夠安全流動，那么在管道正常運(yùn)行的情況下將不會發(fā)生水合物堵塞事故。

3 水合物顆粒聚集特性

3.1 水合物顆粒聚集機(jī)理

輸油管道中防止水合物聚集是研究的重點(diǎn)，通常情況下管道中的自由水在流動剪切和天然或人工乳化劑的作用下以小液滴形式均勻地分散在油相中[16]。Turner等[17]利用聚焦光束反射測量（FBRM）原理的粒度儀測量了油包水乳狀液中的水合物顆粒，測量結(jié)果發(fā)現(xiàn)生成的水合物顆粒的粒徑與初始乳狀液水滴的粒徑一致，也就是說水滴直接作為反應(yīng)單元轉(zhuǎn)化為水合物顆粒，并且根據(jù)實(shí)驗(yàn)結(jié)果建立了以油相黏度、油水界面張力、剪切速率為函數(shù)的計(jì)算公式。在此基礎(chǔ)上，Turner等[18]提出水合物顆粒生長的“收縮核”模型（圖6），認(rèn)為水合物顆粒的生長過程大致為：在低溫、高壓作用下首先在油水界面成核，晶核在水滴表面不斷生長，形成包繞水滴的水合物殼，溶解在油相中的氣體穿過水合物殼使水合物殼內(nèi)未轉(zhuǎn)化水繼續(xù)生成水合物，水合物殼由外向內(nèi)生長，水滴內(nèi)核的轉(zhuǎn)化速率取決于氣體分子通過水合物殼的擴(kuò)散速率，最終水滴部分甚至完全轉(zhuǎn)化成固體水合物顆粒。初始乳狀液水滴粒徑是研究水合物聚集的基礎(chǔ)，Boxall[19]改進(jìn)了初始水滴粒徑計(jì)算公式，并應(yīng)用于 OLGA 軟件的CSMHYK 計(jì)算模塊，用于預(yù)測運(yùn)行管道中水合物的聚集情況。

圖6 油包水型乳狀液中水合物殼生長模型

圖7 油包水體系水合物生長聚集過程

水合物在管道中可以分為3個(gè)階段：成核、生長、聚集，如圖7所示，但至今對水合物顆粒的聚集機(jī)理認(rèn)識尚不充分，無統(tǒng)一定論。Austvik等[20]首先給出了一個(gè)水合物顆粒聚集機(jī)理的說法，他指出水合物顆粒在沒有自由水的體系很容易輸送而不聚集，聚集主要發(fā)生在水合物形成階段，一旦水合物生成完成，體系中的自由水大部分都轉(zhuǎn)化為水合物，這種情況下水合物可以被穩(wěn)定輸送而不會聚集，認(rèn)為是水合物形成階段有大量自由水存在，顆粒表面的強(qiáng)親水性使顆粒間產(chǎn)生液橋?qū)е旅?xì)液橋力，顆粒表現(xiàn)出“黏性”從而導(dǎo)致了水合物顆粒的聚集。Camargo和Palerm[9]認(rèn)為，水合物顆粒表面液橋是聚集的主要原因，毛細(xì)液橋力是聚集的主要黏附力，并指出顆粒間水橋也會固化形成水合物橋，導(dǎo)致不可逆聚集過程，水合物顆粒表面的潤濕性是產(chǎn)生聚集的主要因素，瀝青質(zhì)吸附在顆粒表面使得水合物由親水性轉(zhuǎn)化為親油性，可以防止聚集的發(fā)生。Fidel-Dufour等[21]認(rèn)為，水合物顆粒的聚集過程就像是一個(gè)準(zhǔn)化學(xué)反應(yīng)，聚集不是因?yàn)橐簶虻淖饔?，而是水合物顆粒（聚集體）與水滴接觸，水滴迅速成核結(jié)晶成為水合物顆粒并黏附在原有水合物顆粒（聚集體）上，或者水滴在兩個(gè)水合物顆粒（聚集體）之間成核結(jié)晶并起到黏結(jié)作用。Colombel等[22]在前人研究的基礎(chǔ)上，提出了接觸誘導(dǎo)聚集機(jī)理和剪切限制聚集機(jī)理統(tǒng)一于一體來解釋水合物聚集的理論，接觸誘導(dǎo)聚集機(jī)理指出水合物聚集是水滴與水合物顆粒接觸的結(jié)果；剪切限制聚集機(jī)理則主要從受力角度出發(fā)，考慮了水合物聚集體受到的流體動力與黏附力平衡對水合物聚集的限制作用。

3.2 水合物顆粒聚集模型

目前針對不同的聚集機(jī)理，眾學(xué)者提出了不同的聚集模型來預(yù)測水合物的聚集過程，其中Colombel等[22]提出的聚集模型認(rèn)可度最高。他以群體平衡模型為基礎(chǔ)，將接觸誘導(dǎo)聚集機(jī)理和剪切限制聚集機(jī)理兩個(gè)理論有機(jī)的結(jié)合，建立了聚集模型來說明水合物顆粒的聚集過程，該模型中接觸誘導(dǎo)聚集機(jī)理與聚集過程相關(guān)聯(lián)，剪切限制聚集機(jī)理與破碎過程相關(guān)聯(lián)，水合物顆粒的聚集過程可以用式（3）表示。

水合物聚集體的有效體積分?jǐn)?shù)見式（4）。

根據(jù)有效體積分?jǐn)?shù)φeff就可求得體系黏度變化情況，計(jì)算結(jié)果與流變儀測量黏度值吻合，該模型能夠在一定程度上準(zhǔn)確預(yù)測體系的黏度變化情況，但是該模型認(rèn)為水合物顆粒聚集過程是一個(gè)完全可逆的類化學(xué)反應(yīng)過程，這與液橋固化的不可逆性不符，模型仍待改進(jìn)。

4 結(jié)論和建議

水合物微觀特性是現(xiàn)今水合物研究的熱點(diǎn)和難點(diǎn)，水合物顆粒受力和聚集特性是研究水合物顆粒微觀特性的重點(diǎn)，目前研究表明毛細(xì)液橋力是水合物顆粒間的主要黏附力，是導(dǎo)致顆粒聚集的主要原因，實(shí)際生產(chǎn)中就可以從消除毛細(xì)液橋力入手，保證水合物漿安全流動。

水合物顆粒微觀受力和聚集特性是保證水合物漿穩(wěn)定輸送、管道安全流動的關(guān)鍵，具有重要的研究意義，作者就當(dāng)前的研究狀況為未來的研究提出幾點(diǎn)建議。

（1）繼續(xù)深入研究水合物顆粒微觀受力特性，顆粒受力是水合物顆粒聚集發(fā)生的決定因素，從受力角度入手提出水合物漿穩(wěn)定措施，保證安全流動。

（2）眾學(xué)者對水合物顆粒聚集機(jī)理和模型進(jìn)行了大量研究，但是至今仍無統(tǒng)一定論，實(shí)驗(yàn)觀測和理論研究結(jié)合是提出準(zhǔn)確水合物顆粒聚集機(jī)理和模型的關(guān)鍵。

（3）加強(qiáng)防聚集（AA）的研究，確定防聚集對水合物顆粒微觀受力和聚集過程的影響特點(diǎn)和本質(zhì)，為水合物漿穩(wěn)定流動提供指導(dǎo)。

符號說明

d——液橋浸沒高度，m

dA，dA，max，dp——分別為水合物聚集體、聚集體最大臨界、初始水合物顆粒的直徑，m

f——分形維數(shù)，量綱為1

FA——水合物顆粒（聚集體）間聚集力，N

FCBT——水合物顆粒（聚集體）間毛細(xì)液橋力，N

H——液橋長度，m

N0，N1，N2，Ni-1，Ni，Ni+1——水滴、含1、2、i-1、i、i+1個(gè)水合物顆粒聚集體的數(shù)濃度，個(gè)/m3

N0(t)，Ni(t)——t時(shí)刻水滴和含有i個(gè)水合物顆粒聚集體的數(shù)濃度，個(gè)/m3

R*——水合物顆粒（聚集體）調(diào)和半徑，m

t——時(shí)間，s

γ——剪切速率，s-1

γLL——油水界面張力，N/m

θ——接觸角，（°）

θk——結(jié)晶特征時(shí)間，s

μ——連續(xù)相流體黏度，Pa·s

φ——水合物體積分?jǐn)?shù)，%

φeff——水合物有效體積分?jǐn)?shù)，%

φmax——水合物聚集體最大的填充分?jǐn)?shù)，%

上角標(biāo)

A——聚集

F——破裂

*——調(diào)和平均值

下角標(biāo)

0，1，2 ，i-1，i，i+1——水滴、聚集體含有1、2、i-1、i、i+1個(gè)水合物顆粒

A——水合物聚集體或聚集過程

CBT——毛細(xì)液橋

eff——有效值

LL——液體與液體相互作用

max——最大值

p——水合物顆粒

[1]Lund A.Comments to some preliminary results from the exxon hydrate flow loop[J].Annals New York Academy of Sciences，1994，715：447-449.

[2]Johan Sjoblom，Bodhild ?vrevoll，Gunn Heidi Jentoft，et al.Investigation of the hydrate plugging and non-plugging properties of oils[J].Journal of Dispersion Science and Technology，2010，31：1100-1119.

[3]Yang Sung-oh，Kleehammer Derek M，et al.Temperature dependence of particle–particle adherence forces in ice and clathrate hydrates[J].Journal of Colloid and Interface Science，2004，277：335-341.

[4]Taylor Craig J，Dieker Laura E，Miller Kelly T，et al.Micromechanical adhesion force measurements between tetrahydrofuran hydrate particles[J].Journal of Colloid and Interface Science，2007，306：255-261.

[5]Dieker Laura E，Aman Zachary M，George Nathan C，et al.Micromechanical adhesion force measurements between hydrate particles in hydrocarbon oils and their modifications[J].Energy Fuels，2009，23：5966-5971.

[6]Aman Zachary M，Aspenes Guro，Dendy Sloan E，et al.The effect of chemistry and system conditions on hydrate interparticleadhesion forces toward aggregation and hydrate plug formation[C]//Woodlands，Texas：SPE International Symposium on Oilfield Chemistry，2011.

[7]Anklam Mark R，Dalton York J，Helmerich Luke，et al.Effects of antiagglomerants on the interactions between hydrate particles[J].American Institute of Chemical Engineers，2007，54（2）：565-574.

[8]McCulfor Jonathan，Himes Paul，Anklam Mark R.The effects of capillary forces on the flow properties of glass particle suspensions in mineral oil[J].American Institute of Chemical Engineers，2011，57（9）：2334-2340.

[9]Camargo Ricardo，Palerm Thierry.Rheological properties of hydrate suspensions in an asphaltenic crude oil[C]//Yokohama，Japan：Proceedings of the 6th International Conference on Gas Hydrates，2002.

[10]Keijo Kinnari，Catherine Labes-Carrier，Knud Lunde，et al.Hydrate plug formation prediction tool-an increasing need for flow assurance in the oil industry[C]//Vancouver，British Columbia，Canada：Proceedings of the 6th International Conference on Gas Hydrates，2008.

[11]Davies Simon R，Boxall John A，Koh Carolyn A，et al.Predicting hydrate-plug formation in a subsea tieback[J].SPE Production &Operations，2009，24（4）：573-578.

[12]Davies Simon R，Boxall John A，Dieker Laura E，et al.Predicting hydrate plug formation in oil-dominated flowlines[J].Journal of Petroleum Science and Engineering，2010，72：302-309.

[13]Davies Simon R，Boxall John A，Dieker Laura E，et al.Improved predictions of hydrate plug formation in oil-dominated flowlines[C]//Houston，Texas，USA ：Offshore Technology Conference，2009.

[14]Zerpa L E，Sloan E D，Koh C A，et al.Hydrate risk assessment and restart procedure optimization of an offshore well using a transient hydrate prediction model[C]//Rio de Janeiro，Brazil ：Offshore Technology Conference，2011.

[15]Zerpa Luis E，Dendy Sloan E，Sum Amadeu K，et al.Overview of CSMHYK ：A transient hydrate formation model[J].Journal of Petroleum Science and Engineering，2012，98-99：122-129.

[16]杜娟.低劑量抑制劑作用下水合物顆粒粘附作用及聚集機(jī)理研究[D].廣州：華南理工大學(xué)，2011.

[17]Turner D J，Miller K T，Sloan E D.Direct conversion of water droplets to methane hydrate in crude oil[J].Chemical Engineering Science，2009，64：5066-5072.

[18]Turner D J.Clathrate hydrate formation in water-in-oil dispersions[D].Golden，Colorado：Colorado School of Mines，2005.

[19]Boxall J A.Hydrate plug formation from ＜50% water content water-in-oil emulsions[D].Golden，Colorado：Colorado School of Mines，2009.

[20]Austvik T，Xiaoyun L，Gjertsen L H.Hydrate plug properties：formation and removal of plugs[C]//Salt Lake City，Utah，USA：Proceedings of the 3th International Conference on Gas Hydrates，1999.

[21]Fidel-Dufour A，Gruy F，Herri J M.Rheology of methane hydrate slurries during their crystallization in a water in dodecane emulsion under flowing[J].Chemical Engineering Science，2006，61（2）：505-515.

[22]Colombel Emilie，Gateau Patrick，Barre Lo?c，et al.Discussion of agglomeration mechanisms between hydrate particles in water in oil emulsions[J].Oil and Gas Science and Technology，2009，64（5）：629-636.