楊蕊 劉德俊 王鶴超

摘? ? ? 要： 以豎直傾斜彎管內(nèi)流動(dòng)的天然氣水合物漿液作為研究對(duì)象，運(yùn)用Fluent軟件進(jìn)行數(shù)值模擬，并通過(guò)FBRM實(shí)驗(yàn)進(jìn)行驗(yàn)證。主要研究了天然氣水合物漿液在管道內(nèi)流動(dòng)、水合物顆粒體積分?jǐn)?shù)的變化導(dǎo)致水合物顆粒不斷聚集情況，因此研究了不同流速、不同管徑對(duì)天然氣水合物顆粒體積分?jǐn)?shù)產(chǎn)生的影響，并對(duì)其影響穩(wěn)定流動(dòng)的情況進(jìn)行分析。研究結(jié)果表明：低流速下，受到浮力的影響，水合物顆粒會(huì)在進(jìn)入彎管前的直管段處在管道上方聚集，在彎管處體積分?jǐn)?shù)達(dá)到峰值;管徑不變的情況下，流速越低，彎管處水合物體積分?jǐn)?shù)越大，聚集范圍也越大;流速不變情況下，管徑越大，水合物體積分?jǐn)?shù)在整個(gè)流動(dòng)過(guò)程中變化越劇烈，彎管處水合物體積分?jǐn)?shù)越大，聚集的風(fēng)險(xiǎn)越高;模擬的結(jié)果與實(shí)驗(yàn)的誤差約為11%。模擬結(jié)果能夠?qū)λ衔餄{液輸送提供理論支持。

關(guān)? 鍵? 詞：水合物漿;數(shù)值模擬;顆粒體積分?jǐn)?shù);豎直彎管

中圖分類(lèi)號(hào)：TQ 022.4? ? ? ?文獻(xiàn)標(biāo)識(shí)碼： A? ? ? ?文章編號(hào)： 1671-0460（2020）11-2597-05

Simulation Study of Hydrate Slurry Flow in Vertical Bend

YANG Rui 1， LIU De-jun 1， WANG He-chao 2

（1. School of Petroleum and Natural Gas Engineering， Liaoning Shihua University， Fushun 113001， China;

2. China Petroleum Pipeline Engineering Co.， Ltd. Southeast Asia Project Department， Langfang 065000， China）

Abstract：? Taking natural gas hydrate slurry flowing in a vertically inclined elbow as a research object， numerical simulation was performed by using Fluent software， and simulation results were verified by FBRM experiments. The natural gas hydrate slurry flowing in the pipeline was mainly studied as well as the hydrate particles continuous accumulation caused by the change of the volume fraction of hydrate particles. Therefore， the effect of different flow rates and different pipe diameters on the volume fraction of natural gas hydrates was studied. The effect of the volume fraction of natural gas hydrates on the steady flow was also analyzed. The research results showed that： at low flow rates， under the effect of buoyancy， hydrate particles accumulated above the pipe at the straight pipe section before entering the bend， and the volume fraction peaked at the bend; When the pipe diameter was constant， the lower the flow rate， the larger the hydrate volume fraction at the bend， and the larger the range of aggregation; Under the same flow rate， the larger the pipe diameter， the more drastic the hydrate volume fraction change in the flow process， the larger the hydrate volume fraction at the bend， the higher the risk of aggregation; The error between the simulation result and the experiment result was about 11%. The simulation results can provide theoretical support for hydrate slurry transportation.

Key words： Hydrate slurry; Numerical simulation; Volume fraction of particles; Vertical elbow

天然氣水合物是由甲烷為主的烴類(lèi)氣體和水構(gòu)成的晶體化合物，由于其高度濃縮的特性非常適用于天然氣的運(yùn)輸和儲(chǔ)存[1]。我國(guó)的南海海域天然氣水合物資源豐富，具有較高的開(kāi)發(fā)價(jià)值，因此研究水合物漿液管道流動(dòng)特性是實(shí)現(xiàn)深海開(kāi)發(fā)利用天然氣水合物資源的重要環(huán)節(jié)，是實(shí)現(xiàn)管輸技術(shù)大規(guī)模工業(yè)應(yīng)用的重要基礎(chǔ)[2-5]。對(duì)于純水體系的水合物漿液在管道中的流動(dòng)特性國(guó)內(nèi)外專(zhuān)家學(xué)者進(jìn)行了廣泛的研究。王武昌[6]等進(jìn)行了水合物顆粒體積分?jǐn)?shù)為0～65.2%的四氫呋喃水合物漿液流動(dòng)實(shí)驗(yàn)，提出臨界水合物體積分?jǐn)?shù)，在臨界以下為牛頓流體，超過(guò)臨界時(shí)，壓降梯度隨體積分?jǐn)?shù)的增加急劇增加，漿體為Bingham流體。JOSHI[7]通過(guò)循環(huán)管路中水合物形成實(shí)驗(yàn)，提出純水系統(tǒng)管道中天然氣水合物堵塞的機(jī)理，將天然氣水合物形成分為3個(gè)區(qū)域，以1區(qū)到2區(qū)即均相到非均相時(shí)天然氣水合物濃度作為水合物堵塞的指標(biāo)，而這一指標(biāo)是由壓力、溫度、速度方向決定的。DING[8]等同樣根據(jù)水合物濃度并考慮了水合物顆粒的直徑，由此提出了定量估算水合物團(tuán)聚度和沉積度的幾種方法以保障水合物漿液流動(dòng)穩(wěn)定。

由此可見(jiàn)，水合物的體積分?jǐn)?shù)是決定水合物漿液安全流動(dòng)的重要參數(shù)。這在Ahmad和李玉星[9-10]等進(jìn)行的試驗(yàn)中也有所運(yùn)用?？紤]到在深海開(kāi)發(fā)天然氣水合物漿液輸送過(guò)程中，豎直彎管系統(tǒng)是傳輸中常見(jiàn)結(jié)構(gòu)，因而在本文中探討水合物漿液在管道內(nèi)的流動(dòng)時(shí)水合物體積分?jǐn)?shù)變化以及管徑、入口流速對(duì)豎直彎管內(nèi)水合物顆粒最大體積分?jǐn)?shù)的影響。

1? 模型建立

1.1? 幾何模型建立

如圖1所示， 選取直徑40 mm，彎徑比R/D=1.5，傾角135°的豎直彎管，前后各取500 mm以克服彎管端面效應(yīng)。

1.2? 數(shù)學(xué)模型建立

基于顆粒動(dòng)力學(xué)理論的歐拉雙流體模型，采用Fluent軟件建立模擬水合物漿液在豎直彎管內(nèi)流動(dòng)的數(shù)值模擬，假設(shè)流動(dòng)過(guò)程無(wú)相變并且相間無(wú)質(zhì)量和熱量傳遞，水合物顆粒為均勻球體，選用RNG k-ε湍流模型，選取gidaspow模型為相間耦合的曳力模型[11-14]。

1.2.1? 基本控制方程

在假設(shè)流體具有不可壓縮性的前提下，將連續(xù)性方程與N-S方程均化后，得到控制方程：

。? ? ? ? ? ? （1）

。（2）

式（1）與式（2）中： xi（i =1，2，3）與xj（j =1，2，3）代表坐標(biāo)分量;ui（i =1，2，3）與uj（j =1，2，3）代表時(shí)均速度分量;p代表時(shí)均壓力 ;μ代表動(dòng)力黏度 ;ρ代表流體密度 ; 為雷諾應(yīng)力分量。

1.2.2? 湍流基本方程

相比于標(biāo)準(zhǔn)的k-ε模型，重整化群RNG k-ε模型在ε方程中加了一個(gè)條件，提高了湍流漩渦方面的精度。這就使RNG k-ε模型可以更好地處理中等強(qiáng)度的旋流和低雷諾數(shù)流動(dòng)。其基本表達(dá)形式為：

。（3）

。（4）

在（3）式和（4）式中，k為湍動(dòng)能， m2·s-3;μ1為紊流黏度，Pa·s;G為湍動(dòng)能耗散相，G=μ1（▽?duì)?+（▽?duì)?）T） ▽?duì)?;σk為k方程的湍流普朗特?cái)?shù)，σk=1.0;σ?為?方程的湍流普朗特?cái)?shù)，σ_?=K^2/（[C_1-C_（2 ）]√（C_μ ））? ? ;C1=1.42;C2=1.68。

1.2.3? 相間耦合

CFD模擬中液固相間的耦合作用是通過(guò)相間動(dòng)量交換來(lái)實(shí)現(xiàn)的，選取gidaspow模型為相間耦合的曳力模型，符合水合物漿液的液固兩相間曳力作用，其基本表達(dá)式如下：

; （5）

;? （6）

。? ? ? ? ?（7）

在（5）、（6）和（7）中，β為固液傳熱系數(shù);Res為水合物顆粒雷諾數(shù);CD為曳力系數(shù)，與雷諾數(shù)相關(guān);ds代表水合物顆粒平均粒徑，m。

2? 模型驗(yàn)證

2.1? 實(shí)驗(yàn)驗(yàn)證

Joshi的實(shí)驗(yàn)中，水合物體積分?jǐn)?shù)在達(dá)到7%時(shí)由均相向非均相過(guò)渡，因此設(shè)置入口水合物體積分?jǐn)?shù)7%，超出這一數(shù)值有水合物沉積風(fēng)險(xiǎn)。為了驗(yàn)證模擬可行性，如圖2所示，利用自主搭建的實(shí)驗(yàn)平臺(tái)，對(duì)管道中水合物顆粒直徑與顆粒數(shù)進(jìn)行測(cè)定，實(shí)驗(yàn)平臺(tái)主要包括兩部分：第一部分由直徑為? ? 40 mm，長(zhǎng)度為3 000 mm的（4）PVC管路連接水槽與閥門(mén)，中間連接如圖2的實(shí)驗(yàn)管段;第二部分由梅特勒-托利多提供的（8）聚焦光束反射測(cè)量?jī)x（FBRM-D600L）測(cè)量系統(tǒng)以及型號(hào)為? ?QSF10-10-0.55的潛水泵（2）和型號(hào)為PS-61005的普斯變頻電源組成的環(huán)路系統(tǒng)。測(cè)量系統(tǒng)主要用于對(duì)管內(nèi)處于不同流動(dòng)速度下的水合物顆粒的直徑及數(shù)量分布進(jìn)行測(cè)量，測(cè)量位置分別為直管和傾管中點(diǎn)處（實(shí)驗(yàn)管段對(duì)應(yīng)處鉆取直徑為18 mm的圓形孔洞）。具體的實(shí)驗(yàn)流程是首先在水槽（3）中加入代替水合物晶體的冰粒子和去離子水，使冰粒子體積分?jǐn)?shù)達(dá)到7%，之后將探頭插入測(cè)量位置的孔洞中，之后啟動(dòng)潛水泵電源開(kāi)關(guān)，打開(kāi)實(shí)驗(yàn)程序開(kāi)始測(cè)量。測(cè)量時(shí)間為30 min，數(shù)據(jù)記錄間隔為10 s。

由于啟動(dòng)潛水泵后冰漿開(kāi)始流動(dòng)，測(cè)得粒子數(shù)不斷增加，之后粒子數(shù)不斷下降是因?yàn)檠h(huán)流動(dòng)的冰粒子被泵葉輪切割，粒徑不斷減小，10 min后趨于穩(wěn)定，因此數(shù)據(jù)從10 min之后開(kāi)始記錄，得到不同直徑粒子數(shù)分布圖（圖3）。實(shí)驗(yàn)結(jié)果發(fā)現(xiàn)A點(diǎn)較E點(diǎn)水合物體積分?jǐn)?shù)稍高，這與模擬所得數(shù)值情況相同。實(shí)驗(yàn)重復(fù)測(cè)量5次，由粒子數(shù)n與粒徑d計(jì)算該處水合物體積分?jǐn)?shù)φ，與相同條件下模擬所的數(shù)值進(jìn)行比較（圖4），實(shí)驗(yàn)所得φ值較模擬值高，是由于鏡頭探入管道的部分導(dǎo)致截面變小。

2.2? 管道截面處云圖分析

為了更好地研究水合物漿液在管道流動(dòng)時(shí)體積分?jǐn)?shù)的變化，在與粒徑分布實(shí)驗(yàn)相同的流速、入口體積分?jǐn)?shù)條件下進(jìn)行了數(shù)值模擬，首先得到了如圖5所示的直徑40 mm管道在流速0.55 m·s-1時(shí)A到E截面水合物體積分?jǐn)?shù)云圖。水合物漿在流入直管段后，A截面顯示水合物顆粒向管道上方聚集，到彎道入口B截面時(shí)聚集更明顯，這是由于水合物顆粒密度比水低的性質(zhì)導(dǎo)致在流動(dòng)的過(guò)程中受到浮力的影響不斷向管道上方聚集。在進(jìn)入彎管結(jié)構(gòu)后，渦流降低了聚集的效果，但是彎管結(jié)構(gòu)導(dǎo)致聚集范圍集中，這是由于水合物漿液進(jìn)入彎管后，彎管外側(cè)壓力減小，速度升高，顆粒體受到馬格努斯力和壓力梯度力的影響向高壓低速方向聚集[15]，在彎道出口D截面水合物顆粒小范圍聚集更加明顯。進(jìn)入傾管段后，管壁處沒(méi)有明顯聚集效果，水合物顆粒在管道軸線(xiàn)附近輕微聚集。

3? 水合物顆粒體積分?jǐn)?shù)影響因素分析

水合物漿液在流入豎直彎管后發(fā)生明顯的聚集現(xiàn)象，經(jīng)過(guò)一段時(shí)間后可能會(huì)有堵塞現(xiàn)象發(fā)生，為了水合物漿液的穩(wěn)定安全流動(dòng)，研究了0.5、1、? ?2 m·s-1不同流速，25、32、40 mm不同管徑下管道內(nèi)水合物顆粒體積分?jǐn)?shù)變化，水合物漿液在流入彎管和流出彎管進(jìn)入傾管段過(guò)后表現(xiàn)完全不同，在彎管處水合物體積分?jǐn)?shù)增長(zhǎng)最明顯，因此著重討論彎管處水合物顆粒的聚集情況。

3.1? 流速對(duì)水合物體積分?jǐn)?shù)的影響

水合物漿液流入彎管后沿彎管外側(cè)流體流速增大，壓強(qiáng)減小，導(dǎo)致內(nèi)測(cè)水合物聚集，水合物體積分?jǐn)?shù)增大，根據(jù)水合物漿液流動(dòng)云圖可知體積分?jǐn)?shù)在軸向變化更加明顯，選取B、D截面不同流速下軸向水合物體積分?jǐn)?shù)變化進(jìn)行比較。管徑不變條件下，流速越低，彎管入口處體積分?jǐn)?shù)變化越明顯，水合物聚集處體積分?jǐn)?shù)越高;彎管出口處水合物聚集處不僅體積分?jǐn)?shù)高，聚集的范圍也更小，即與高流速比較水合物聚集更集中。

3.2? 管徑對(duì)水合物體積分?jǐn)?shù)的影響

由于水合物聚集的方位都在管道頂部，因此對(duì)于不同管徑對(duì)水合物體積分?jǐn)?shù)的影響，根據(jù)模擬結(jié)果以管道頂部軸向水合物體積分?jǐn)?shù)變化作為參考依據(jù)，圖6顯示在流速1 m·s-1時(shí)不同管徑下水合物體積分?jǐn)?shù)發(fā)生明顯變化，主要分為以下幾點(diǎn)：

①相同流速下管徑越大水合物體積分?jǐn)?shù)變化越劇烈;②直管段位置，管徑越大水合物顆粒受浮力影響越明顯，水合物顆粒會(huì)更早地在管道上部聚集;③豎直彎管位置，管徑越大水合物體積分?jǐn)?shù)最大值越大，水合物聚集風(fēng)險(xiǎn)越高，同時(shí)進(jìn)出彎管水合物體積分?jǐn)?shù)變化劇烈，參考2.2分析在此處漿體發(fā)生的繞流更加劇烈。見(jiàn)圖7。

4? 結(jié) 論

根據(jù)Joshi的實(shí)驗(yàn)以7%為水合物顆粒向聚集發(fā)展的臨界值，超過(guò)這一數(shù)值的純水系水合物漿液會(huì)在水合物顆粒不斷聚集生長(zhǎng)的情況下有堵塞的風(fēng)險(xiǎn)，因此在實(shí)際輸送過(guò)程中應(yīng)避免天然氣水合物體積分?jǐn)?shù)超過(guò)臨界值，保障水合物漿液的穩(wěn)定流動(dòng)。經(jīng)過(guò)實(shí)驗(yàn)和模擬的研究，發(fā)現(xiàn)豎直彎管系統(tǒng)中，水合物顆粒直徑范圍為10～300 μm，根據(jù)粒徑分布發(fā)現(xiàn)大部分粒徑值為50 μm。水合物漿液流入管道后，水合物顆粒會(huì)在進(jìn)入彎管后開(kāi)始聚集，在流出彎管時(shí)達(dá)到最大值，在傾斜管道處只在中部輕微聚集。由于彎管處壓強(qiáng)差和流速差影響，水合物顆粒聚集明顯，因此整個(gè)系統(tǒng)中此處的堵塞風(fēng)險(xiǎn)更高。之后通過(guò)對(duì)3種不同流速、不同管徑的彎管體系進(jìn)行數(shù)值模擬，結(jié)果表明：

1）在管道直徑一定時(shí)，不同流速的水合物漿液流入彎管后在彎管上方聚集，流速越低，上方的水合物體積分?jǐn)?shù)越大。隨后的聚集范圍逐漸增大，流速越低，聚集的范圍也越大，因此在水合物漿液輸送過(guò)程中盡可能提高流速防止發(fā)生堵管。

2）在流速一定時(shí)，管道直徑越大進(jìn)入彎管前受到浮力影響越大，入口上方體積分?jǐn)?shù)初始值越大，小管徑幾乎不受浮力影響，這導(dǎo)致管徑越大彎管處聚集越強(qiáng)烈，在之后的傾管段也發(fā)生急劇的上升，發(fā)生堵塞的可能性更高。

3）為了降低管道內(nèi)水合物聚集、堵塞的風(fēng)險(xiǎn)，在實(shí)際情況中通常泵速固定的情況下盡量選擇小管徑的管道，同時(shí)也要考慮水合物顆粒對(duì)管道的沖蝕，在彎管處會(huì)產(chǎn)生較大的風(fēng)險(xiǎn)，隨著水合物體積分?jǐn)?shù)的增加，水合物漿液的黏度也會(huì)不斷增加向非牛頓流體轉(zhuǎn)化，因此會(huì)加劇堵塞的形成。

參考文獻(xiàn)：

[1]祝有海，吳必豪，劉玉山.天然氣水合物：21世紀(jì)的新能源[M].北京：海洋出版社，2017.

[2]閆柯樂(lè)，鄒兵，姜素霞，等. 水合物漿液流動(dòng)與流變特性研究進(jìn)展[J]. 化工進(jìn)展，2015（7）：1817-1825.

[3]朱超，王武昌，王琳.天然氣水合物漿流變性的研究現(xiàn)狀及發(fā)展趨勢(shì)[J].天然氣工業(yè)，2011，31（7）：77-81.

[4]李長(zhǎng)俊，黃婷，賈文龍. 深水天然氣水合物及其管道輸送技術(shù)[J]. 科學(xué)通報(bào)，2016，61（22）：2449-2462.

[5]王武昌，李玉星，樊栓獅，等. 四氫呋喃水合物漿流動(dòng)特性[J]. 化工進(jìn)展，2010，29（8）：1418-1422.

[6]張威，韓冊(cè)，吳迪，班久慶，等.天然氣水合物的生成機(jī)理與防治研究[J].當(dāng)代化工，2016，45（8）：1953-1954.

[7]JOSHI S V， GRASSO G A， LAFOND P G， et al. Experimental flowloop investigations of gas hydrate formation in high water cut systems [J]. Chemical Engineering Science， 2013（97）：198-209.

[8]DING L， SHI B H， LV X F， et al.Hydrate formation and plugging mechanisms in different gas-liquid flow patterns [J]. Industrial & Engineering Chemistry Research， 2017， 14 （56）： 4173-4184.

[9]MAJID A A A. An investigation on the viscosity and transportability of methane hydrate slurries using a high pressure rheometer and flowloop [D]. Colorado School of Mines，2015.

[10]李玉星，朱超，王武昌，等. 天然氣水合物漿液穩(wěn)定性影響因素的正交實(shí)驗(yàn)[J]. 油氣儲(chǔ)運(yùn)，2011，30（9）：685-689.

[11]魏丁，王武昌，李玉星，等. 管道CCl3水合物漿流動(dòng)特性的數(shù)值模擬[J]. 油氣儲(chǔ)運(yùn)，2016，35（8）：828-832.

[12]韋雪蕾，劉寶玉，潘振，等. 水平管道水合物漿液流動(dòng)特性的數(shù)值模擬[J]. 化學(xué)工程，2018，46（3）：41-46.

[13]江國(guó)業(yè)，王曉婭.水平彎管內(nèi)水合物漿的流動(dòng)特性[J].科技導(dǎo)報(bào)，2015，33（2）：64-68.

[14]孫賢，劉德俊，王文武，等.彎管體系內(nèi)漿液流動(dòng)分析[J].化學(xué)工程，2019，47（6）：58-63.

[15]郭烈錦.兩相與多相流動(dòng)力學(xué)[M].西安：西安交通大學(xué)出版社，2001.