天然氣水合物輸送管道氣液固三相流運動特性分析

吳官紀，潘 振，陳保東，楊 帆

（遼寧石油化工大學 石油天然氣工程學院, 遼寧 撫順 113001）

隨著經濟的不斷發展，能源問題日益凸顯，石油和煤等傳統的能源日益枯竭，這就要求我們尋找和探索新的可替代能源。在尋找過程中，人們發現，蘊藏在海底的天然氣水合物是一種高效、清潔的可代替能源，具有相當可觀的發展前景［1］。

雖然天然氣水合物的儲量很大，但是它通常都是儲藏在海底，而且覆蓋層也很淺，因此，對天然氣水合物的開采和運輸是天然氣水合物的一個應用難題。通過探索，人們發現固態開采法是目前開采天然氣水合物最有效、最合理的開采方法。在固態開采過程中，如何將天然氣水合物通過水力輸送技術運送至海面是最關重要的部分。因此本文利用Fluent軟件對天然氣水合物的豎直輸送管道進行了數值模擬，對管段內部的流場特性進行了分析［2-5］。

通過模擬過程得出天然氣水合物的顆粒濃度、顆粒密度、顆粒粒徑以及天然氣水合物輸送的流速對天然氣水合物的輸送影響。結果表明，天然氣水合物的顆粒粒徑偏小、適度的體積濃度以及較小的天然氣水合物顆粒密度有利于其水力輸送［6,7］。

1 計算流體動力學基本控制方程

1.1 基本方程

（1）質量守恒定律

質量守恒定律是指單位時間內流體微元體中質量的變化與同一時間間隔內流入該微元體的凈質量相等，方程式為：

式中：r —流體密度；

t —時間；

u、v、w —速度矢量在x、y、z方向上的分量。

（2）動量守恒定律

動量守恒定律是指微元體中流體的動量對時間的變化率與外界作用在該微元上的各種力之和相等，方程式為：

式中：u、v、w —速度矢量在x、y、z方向上的分量；

m —動力粘度；

sn、sm、s0是動量守恒定律的廣義源項。

（3）能量守恒定律

能量守恒定律的定義為：微元體中能量的增加率與進入微元體的凈熱流量加上體力與面力對微元體所做的功相等，方程式為：

式中：T —溫度；

cp—比熱容；

k —傳熱系數；

ST—耗散項。

1.2 基本假設

為了便于利用Fluent軟件模擬固、液、氣三相流的流動狀態，常需要做以下假設：

（1）固相天然氣水合物的固體顆粒看作球形，并且假設天然氣水合物的的顆粒粒徑均勻大小不變。

（2)三相流中的固相是連續運動介質，液相(水)是不可壓縮的流體，同時假設固液氣三相的物理性質均是常數。

（3)通常認為氣體為理想氣體狀態。

（4）流體與壁面間無滑移

1.3 三維有限元網格劃分

模型網格質量的優劣對軟件的計算結果的精確性及計算過程的速度起著決定性作用，因此在模擬過程中劃分網格形式有重要的意義。本論文以豎直上升管道中的某一段為模擬對象，利用Fluent 軟件中的 Gambit 軟件開始三維模型的建立工作，在圖1 中可以看到相關的三維模型。

圖1 為三維模型的網格劃分Fig.1 3 d model of meshing

以下是模擬過程中所需的主要幾何參數 D=0.3 m, H=5 m。非結構化網格和結構化網格是構成網格總體的兩部分。其中非結構化網格是指在網格區域內部相鄰兩點之間的毗鄰單元不盡相同；而對于結構化網格則恰好相反，在網格區域內部的毗鄰單元全部相同，這樣建立的網格質量很好，模型網格生成效率也很高，并且利用這種網格對于區域的邊界擬合過程的實現也變得相對輕松，通過使用參數化方法以及樣條插值的方法對空間或者曲面進行擬合，模擬的區域較為光滑，可以接近試劑模擬的模型。

非結構化網格和結構化網格的區別在于，在計算結果的收斂性上，結構化網格形式要比非結構化網格形式優越，但是在模擬復雜的幾何結構過程中，非結構化網格更具有實際的應用價值，可操作性強。在本篇論文中，筆者采用結構化網格形式建模。

1.4 邊界條件

（1）進口條件：將管段底部端進口定義為速度進口邊界，氣液固三相的速度大小根據流量確定，并且在入口方向速度分布均勻，方向垂直于進口面，另給出湍動能、耗散率的預估值及氣液固三相體積濃度及固體顆粒直徑。

（2）壁面條件：所有壁面上流體都滿足無滑移條件，計算時，近壁處采用標準壁面函數法處理。

（3）出口條件：將管道上端出口定義為自由出流出。

2 數值模擬結果與分析

設定天然氣固體顆粒粒徑為3 mm，流速為1 m/s，顆粒密度為2 000 kg/m3，體積分數為0.3，水的速度為1 m/s, 體積分數為0.4，CH4的速度為1 m/s，體積分數為0.3，考察顆粒直徑對管段壓力損失和阻力損失的影響，顆粒直徑變化范圍為3、4、5、6、7、8 mm。通過Fluent 模擬繪制出口附近水相流動速度分布云圖、出口附近固、液、氣流動速度分布云圖、出口附近固、液氣相壓力分布云圖、（圖中以顆粒直徑為3 mm為例）(圖2-7)。

圖2 出口附近氣相流動速度分布云圖Fig.2 Near the export gas flow velocity distribution nephogram

圖3 出口附近液相流動速度分布云圖Fig.3 Export liquid flow velocity distribution near the clouds

圖4 出口附近氣相流動速度分布云圖Fig.4 Near the export gas flow velocity distribution nephogram

通過比較三者的速度變化可得隨著管徑的變化，三相的速度均發生變化，而固相變化的最快，管道中心的速度均為最大，但固相的變化更快且最大速度為1 m/s,而氣相和液相的最大速度均超過開始速度1 m/s，通過細心地發現；氣相的速度在同一點的位置要大于或等于液相的速度，由于本論文管道選取的不夠長，氣相的速度變化與液相的速度變化近似相等，在壁面的速度均為零，此結論完全符合流體力學及多相管流學科。

圖5 出口附近固相流動壓力分布云圖Fig.5 Near the solid flow pressure distribution nephogram

圖6 出口附近液相流動壓力分布云圖Fig.6 Near the liquid flow pressure distribution nephogram

由圖中可以看到固液氣的壓力均是隨著管徑的變化而有規律的變化，在同一點固相壓力大于液相壓力，液相壓力大于氣相壓力，不難發現三者壓力呈對稱狀態，說明壓力變化的規律可以尋找到可計算的公式，為以后的開采天然氣水合物的壓力分析及開采的措施應對方法提供了現實的基礎。

圖7 為出口附近氣相流動壓力分布云圖Fig.7 Near the gas flow pressure distribution nephogram

3 壓降的計算公式及粒徑對壓降的影響

3.1 壓降的計算公式

在海洋中開采天然氣水合物的過程中，輸送管道內部的固液兩相流動的流動特性尚不明確，天然氣水合物的顆粒的密度、顆粒粒徑以及天然氣水合物的體積濃度等參數設置是否符合實際開采過程中的工作還尚未確定，天然氣水合物的顆粒粒徑、顆粒的密度以及體積濃度等參數變化對管道輸送造成的影響。

模型段管道壓降

式中：pin—表示輸送管道流體流進壓力；

pout—輸送管道流體的出口。

模型管道壓力損失

式中：l0—表示所選管道的長度，m；

r —表示液體的密度。

模型管道阻力損失

式中：rm—表示混合流體的密度。

3.2 顆粒粒徑對管段壓力損失和阻力損失的影響

預定顆粒密度為2 000 kg/m3，顆粒體積濃度為30%，流速為1.0 m/s考察顆粒粒徑對管段壓力損失和阻力損失的影響，顆粒粒徑的變化范圍為3、4、5、6、7、8 mm。

通過圖8可以發現壓力損失和阻力損失都是隨粒徑的增大而增大的，具有正相關關系。其中壓力損失與阻力損失的曲線趨近于平行。當天然氣水合物的粒徑在3～7 mm之間時，壓力損失和阻力損失呈現出有規律的變化，粒徑與壓力損失可以用一條正比例關系的函數來表示。但在顆粒粒徑為 8 mm時，圖像函數不在是一條光滑的曲線，壓力損失和阻力損失突然激增，說明當天然氣顆粒直徑為8 mm時，對二者損失的影響很大，是一個轉折點。

圖8 為天然氣固體顆粒直徑與壓力損失的關系Fig.8 The relationship between gas solid particle diameter and the pressure loss

4 結束語

本論文通過對天然氣水合物的形成及分解條件，建立了氣液固三相流的數學模型，運用質量守恒、動量守恒、能量守恒三大控制方程。在 fluent軟件中，找到適合的方法，模擬了天然氣水合物在豎直管道中壓力、密度、速度的變化，使人們在開采的過程中能夠發現更經濟、更合理的開采方法。其中在選取管道時只是選取了管道長度為5 m，發現若選擇的管道過長，會導致模擬結果的失敗，這是美中不足之處。

隨著天然氣水合物顆粒直徑的不斷增大，管段進出口壓力降也在不斷增大，而且壓力降增大的趨勢變得不在趨于平緩，這一點可以說明，大的天然氣水合物的顆粒直徑有利于天然氣水合物的有效輸送。在粒徑在3～7 mm之間的變化范圍內，管段內壓力損失和阻力損失的變化趨于平緩，而當粒徑為8 mm 時，管段內的壓力損失和阻力損失增大明顯。對于輸送天然氣水合物的流速，管段內的壓力損失隨著漿體流速的升高而呈現上升趨勢，管段內的阻力損失也隨著漿體流速的升高而不斷升高。

［1］ 周懷陽, 彭小彤, 葉瑛. 天然氣水合物[M]. 北京: 海洋出版社,2000.

［2］ 史斗, 孫成權, 朱岳年. 國外天然氣水合物研究進展[M]. 蘭州: 蘭州大學出社, 2002.

［3］ 孫志高, 樊栓獅, 郭開華, 等. 天然氣水合物分解熱的確定[J]. 分析測試學報, 2002, 21(3): 7-9.

［4］ 廖健, 梅東海, 楊繼濤, 等. 天然氣水合物相平衡研究的進展[J].天然氣工業, 1998, 18(3), 75-82.

［5］ 祝道平, 寧正福. 利用高能氣體壓裂技術開采天然氣水合物可行性分析[J]. 重慶科技學院學報(自然科學版), 2009, 11(3): 37-39.

［6］ 李淑霞, 陳月明, 杜慶軍. 天然氣水合物開采數值模擬的參數敏感性分析[J]. 現代地質, 2005, 19(1): 108-112.

［7］趙生才, 德國氣水合物研究計劃簡介[J].天然氣地球科學，2001,12(2), 63-67.

［8］ 馬文婧;南海天然氣水合物開發的風險因素分析[D].青島：中國海洋大學,2011.