楊蕊 劉德俊 王鶴超

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

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

中圖分類號：TQ 022.4? ? ? ?文獻標識碼： A? ? ? ?文章編號： 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

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

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

1? 模型建立

1.1? 幾何模型建立

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

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

基于顆粒動力學(xué)理論的歐拉雙流體模型，采用Fluent軟件建立模擬水合物漿液在豎直彎管內(nèi)流動的數(shù)值模擬，假設(shè)流動過程無相變并且相間無質(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）代表坐標分量;ui（i =1，2，3）與uj（j =1，2，3）代表時均速度分量;p代表時均壓力 ;μ代表動力黏度 ;ρ代表流體密度 ; 為雷諾應(yīng)力分量。

1.2.2? 湍流基本方程

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

。（3）

。（4）

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

1.2.3? 相間耦合

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

; （5）

;? （6）

。? ? ? ? ?（7）

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

2? 模型驗證

2.1? 實驗驗證

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

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

2.2? 管道截面處云圖分析

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

3? 水合物顆粒體積分數(shù)影響因素分析

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

3.1? 流速對水合物體積分數(shù)的影響

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

3.2? 管徑對水合物體積分數(shù)的影響

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

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

4? 結(jié) 論

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

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

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

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

參考文獻：

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

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

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

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

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

[6]張威，韓冊，吳迪，班久慶，等.天然氣水合物的生成機理與防治研究[J].當代化工，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)定性影響因素的正交實驗[J]. 油氣儲運，2011，30（9）：685-689.

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

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

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

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

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