水合物漿液在U形彎管中的流動特性

張 藝，李 聰，何騁遠，李 樂，周詩崠

（1. 常州大學 石油工程學院 油氣儲運技術省重點實驗室，江蘇 常州 213016；2. 江蘇城鄉建設職業學院 公用事業學院，江蘇 常州 213147）

天然氣水合物是水和天然氣分子在低溫高壓條件下形成的一類非化學計量型的籠形絡合物［1］，水合物顆粒極易在油氣管道運輸時生成，且在流動過程中會發生沉積，給管道安全運行帶來巨大隱患。因此，研究水合物在管道內的流動特性，可為預防水合物堵管提供基礎理論。水合物漿液的流動參數與顆粒直徑、顆粒體積分數、流動速度等均有一定關系［2］。但目前應用實驗方法對水合物漿液流動進行分析較難實現，因而應用數值模擬［3］的方法成為實驗方法的有效替代手段［4］。Balakin 等［5-6］研究了湍流條件下R11水合物漿液的流動特性，實驗結果表明漿液的流動特性幾乎不會被低濃度的水合物顆粒所影響。Wang等［7-9］用Fluent軟件模擬了水平管道中冰漿的流動，發現了壓降、顆粒粒徑、流速等參數在管道流動過程中相互之間的關系。趙鵬飛等［10］通過軟件分析了水平管路中漿液的流動特性，研究了固相顆粒粒徑、顆粒黏度和最大的內相堆積率等物性參數對漿液流動的影響。Camargo粒徑模型［11］比 Mulhe粒徑模型［12］更符合實際流動工況，與實驗數據吻合。Balakin等［13］通過數值模擬和實驗相結合的方法，還可準確預測管內沉積情況。上述模擬均未考慮水合物顆粒濃度的變化及彎管對流速和壓力的影響。

本工作以實驗室水合物高壓環路設備為基準，建立1∶1的U形管物理模型對直管和彎管進行對比分析，利用Fluent 15.0軟件對模型進行耦合求解，進而模擬流速、濃度及管型對管內水合物流動特性的影響。

1 數值模擬

1.1 幾何模型的建立

圖1為U形管的物理模型及網格劃分。

圖1 U形管物理模型（a）及網格劃分（b）Fig.1 Physical model(a) and grid division(b) of U-tube.

由圖1可知，物理模型是規則的U形管，模型由兩段相同的直管（長度為5.5 m，內徑為0.254 m）和一段半圓形的彎管（半徑為6.3 m）構成。網格采用多塊對接網格，為了準確地計算近壁面處的流動情況，從而提高模擬結果的精確性，對管道的出入口近壁面處采用加密網格，再進行網格無關性檢驗，確定網格數為1 200 000。

1.2 數值模型

為使模擬更加準確，在模擬計算過程中應用以下假設：水合物漿液在流動體系中達到平衡，流動為等溫過程且不計相間質量傳遞，即在流動過程中無水合物分解、生成、聚積和破碎，漿液只由水合物顆粒相和水相組成，忽略管壁對水合物顆粒的黏附作用。在上述條件下，為保證管內流動過程中水合物顆粒為連續介質，本工作選用歐拉-歐拉雙流體模型。

采用Fluent 15.0軟件進行求解計算。模擬入口條件設置為速度入口，根據國內外現有研究成果［14-15］，在入口處定義初始濃度分別為37.7%，27.5%，25.4%，20.7%，研究粒徑選取248，355，450 μm；出口條件設置為自由流出，壓力-速度耦合采用SIMPLE模型，壁面采用無滑移邊界條件［16］，因為在彎管處可能會出現二次流現象，所以湍流模型采用RNG k-ε湍流方程［17］，當非穩態數值計算各因子殘差收斂到10-5時，認為計算結束。所有的物理模型以及所有的流動過程都是以質量守恒定律為基礎。管輸水合物漿液流動體系的動量方程符合牛頓第二定律的條件。

1.3 相間耦合

水合物漿液流動屬于氣-液-固多相流動，不同于單相流體流動。在模擬過程中，假設氣體已消耗完成，即混合流動是液固兩相流。當固體顆粒載荷較高時，要考慮相間耦合及曳力作用，本工作采用模型為Wen-Yu曳力模型［18］。

1.4 有效性分析

管道漿液輸送時壓力梯度變化是衡量管道運輸的重要運行參數，根據本課題組高壓可視水合物漿液實驗環道的實驗結果與模擬結果進行對比［19-20］，發現模擬壓降與實際壓降的相對誤差均在20%以內，與王繼紅等［8，21-22］國內外研究報道結果相一致，說明通過模擬的方法進行管道水合物流動分析具有可行性。

2 結果與討論

2.1 流速分布

圖2為水合物顆粒粒徑及顆粒濃度對直管和彎管處截面Z方向軸線上的流速分布的影響。由圖2可知，隨著水合物顆粒粒徑或顆粒濃度的增大，水合物顆粒流速也會小幅度的增大。由圖2a，c可知，在直管橫截面處流速的分布表現為下拋物線狀，在管道中軸位置處流速最大；且流速在管道中分布不均勻，由管道中心向管壁處逐漸降低，這與Lahiri等［23］模擬結果相符；管壁處速降低的幅度較大，這是由于近壁面處流速低，對水合物攜帶能力差，所以水合物顆粒流動速度分布呈開口向下的拋物線；且呈現出很好的對稱性。由圖2b，d可知，由于離心力和阻力的影響，在彎管橫截面處漿液流速呈大M形分布，以管道中軸位置處為分界點，上下區域呈對稱分布，這與王曉婭［14］通過模擬水合物漿液在彎管處流速分布所得結果一致。

圖2 水合物顆粒粒徑及顆粒濃度對直管（a，c）和彎管（b，d）處截面Z方向軸線上的流速分布的影響Fig.2 Effects of hydrate particle diameter and particle concentration on velocity distribution of Z-axis of section at straight(a, c) and bend(b, d) pipe.

圖3為彎管入口和出口橫截面處流速分布。由圖3a可知，彎管入口處流速分布較均勻，以管軸線為中心，上下呈半圓形對稱分布。由圖3b可知，彎管出口部分的速度呈凹凸的橢圓形對稱分布，這是因為離心力的作用以及流體撞擊壁面受到阻力的影響，以管段中軸處向管壁處流速由大變小；最大速度偏向彎管截面內凹側的位置，且彎管截面部分外側速度急速下降。由于流體的動量動能不斷發生交換，導致漿液流速在管道內部分布不均勻，且彎管處流速大于直管處流速。這種彎管二次流現象，主要由于流體在管道中流動時，彎管處會受到離心作用，流速越大離心作用也就越明顯。受離心力作用的影響，水合物漿液首先會向彎管部分外凸側流動，然后由于內側出現真空，受壓力作用再流回內側。此流動現象與四氫呋喃在彎管處流動現象相似［14］。

2.2 顆粒濃度分布

圖4為水合物顆粒粒徑對管道中水合物顆粒濃度分布的影響。由圖4a可知，在直管截面Z方向上顆粒濃度分布可分為3部分：高濃度區域（管道的上部）、在管的中心部分幾乎與入口濃度一致的均勻分布的區域和在下管壁附近具有最低濃度的區域［24］；且隨著固相顆粒粒徑的增大，管道中水合物顆粒濃度變化的幅度也增加。因為中間流體的流速比其他區域大，剪切作用較為明顯，影響顆粒的分布，使顆粒具有離散分布的趨勢，且由于濃度的增加，導致流體的湍流特性增大，使液相與顆粒相混合更加充分，從而使管道中部區域處于均相懸浮流的狀態［25］。在近壁面處，漿液流速較低，顆粒具有聚集的趨勢，表現出不均勻流動狀態。同時因為管道中的水合物密度小于水的密度，會出現水合物顆粒更多地聚集在管道上部的現象。由圖4b可知，在彎管截面Z方向上管道中部區域處水合物顆粒濃度分布整體較為均勻，在靠近管壁處顆粒濃度變化較大，管道上、下部區域，顆粒濃度均減小。隨著顆粒粒徑的增大，顆粒濃度也隨之增大；顆粒粒徑為450 μm時，顆粒濃度波動的范圍變化較大，中部區域顆粒濃度范圍為27.5%～28.3%。彎管處顆粒濃度分布呈中部高于管壁處的趨勢，出現該現象主要是因為水合物漿液在流經彎管時，不僅受到慣性力和曳力作用，還受到離心力的作用，且離心力最大，起主導作用；與此同時由于水合物的密度小于水的密度，即水合物顆粒受到的離心力小于液相受到的離心力；因此出現管道中部顆粒濃度大于上下管壁處顆粒濃度的現象。

圖3 彎管入口（a）和出口（b）橫截面處流速分布Fig.3 Flow rate distribution at the cross section of the bend inlet(a) and outlet(b).

圖4 直管（a）和彎管（b）截面Z方向軸線上水合物顆粒濃度分布Fig.4 Hydrate particle concentration distribution on the Z-axis of straight(a) and bend(b) pipe section.

2.3 管道沿線漿液流場分布規律

圖5為顆粒粒徑和漿液流速對管軸線處漿液體積分數的影響。由圖5a可知，漿液流經管道入口處時，水合物漿液體積分數先降低后升高，然后趨于穩定狀態；在流經彎管處，漿液體積分數急速增大，最后減小并趨于穩定。由圖5b可知，在流體流動穩定狀態下，流速為1，2，3 m/s對應的漿液體積分數均在27.6%附近。整個流動過程中，粒徑對漿液體積分數在管道中的分布影響很小，彎管對漿液體積分數影響相對較大。而漿液流速對水合物漿液體積分數的分布影響很小，可忽略不計。流速越高，流體對水合物顆粒的攜帶和分散能力就越強，水合物顆粒分布越均勻。這與宋光春等［21］的模擬結果一致。

2.4 彎曲段壓力分布

圖6為流體流經彎管處沿X軸方向以及Z軸方向的壓力分布。

由圖6a可知，流體流經彎管處，X軸方向上壓力分布不均，從左（管道內凹側）到右（管道外凸側）壓力逐漸增大，在管壁外凸側0.012 7 m附近，壓力達到最大。由圖6b可知，在Z軸方向上，壓力呈倒U形分布，在管道中部區域，壓力存在波動，但波動范圍不大，在管壁處壓力降低。上述現象是因為流體流經彎管處，受離心作用的影響產生彎管二次流［26］，流體流動狀態發生改變；流體會向外凸側集中，導致內凹側流體流失，形成中空區域，為保持穩定流動，外凸側的流體會向內凹側流動填補中空區域，因此從左（管道內凹側）到右（管道外凸側）壓力逐漸增大。且由于流體在彎管處受到離心力的作用，彎管處壓力由內側到外側逐漸增大。這與紀宏超等［27］用Fluent模擬的彎管處流動特性一致。

圖5 顆粒粒徑（a）及漿液流速（b）對管軸線處顆粒濃度分布的影響Fig.5 Influence of particle diameter(a) and slurry flow rate(b) on particle concentration at pipe axis.a Slurry flow rate 3 m/s；b Particle dimeter 248 μm

圖6 彎管處沿X軸線方向上（a）及Z軸線方向上（b）壓力分布Fig.6 Pressure distribution along the X axis(a) and Z axis(b) at the bend pipe.

3 結論

1）在直管橫截面處流速呈拋物線分布，流速在管道中間區域最大，靠近管壁流速逐漸減小。在彎管橫截面處流速呈M形分布，流速從管道中心到管壁處，先增大后減小。相對于水合物顆粒濃度、顆粒粒徑對流速分布的影響，彎管對漿液流速分布影響更大。

2）水合物的密度小于水的密度且受顆粒自身重力等因素的影響，在管輸過程中，受流動剪切的影響，顆粒分布不均勻，靠近上管壁區域顆粒濃度較高，靠近下管壁區域顆粒濃度較低，管中部區域顆粒濃度均衡。

3）漿液流經彎管處時，產生二次流動現象，且由于離心力的影響，在管道X軸線處外凸側的壓力大于內凹側的壓力值；在管截面Z軸處壓力呈倒U形分布。