馬希金，張亞瓊，張 潮，崔生磊

(1.蘭州理工大學能源與動力工程學院， 甘肅 蘭州 730050；2.甘肅省流體機械及系統(tǒng)重點實驗室， 甘肅 蘭州 730050)

隨著陸上石油的日益減少，人們越來越重視海底石油的開發(fā)。海上原油產(chǎn)出物主要是油和氣的混合物，同時還含有少量的水和砂，是一種多相混合物[1]。在傳統(tǒng)的海上石油開采過程中，需要建立一套復(fù)雜的處理設(shè)施，使原油產(chǎn)出物在海上平臺或海底就分離，然后將油、氣分別輸送到陸上加工基地。這種開采工藝增加了油氣分離裝置，不僅在技術(shù)上帶來了困難，而且非常不經(jīng)濟[2]，由此催生了油氣混輸泵。油氣混輸泵作為軸流式多相泵的一種，它具備了液相泵與壓縮機的雙重性能，替代了一般原油輸送過程中的分離裝置與輸送管道，體現(xiàn)出了其顯著的經(jīng)濟效益。

早在20世紀初，國外的石油企業(yè)便開始研究油氣混輸技術(shù)，于 1984 年開始制定的海神計劃[3-6]，由法國石油研究院(IFP)、法國國家石油公司(Total)及挪威石油開發(fā)公司(Statoil)共同合作的一個計劃項目，所研究的主要課題之一便是多相混輸泵。海神泵選用的是NACA翼型，經(jīng)過好幾代的實驗與優(yōu)化，能夠有效地避免氣液分離，在性能上得到了很大的提高，并且在不斷的探索改進之后，成功地應(yīng)用于油田開采上。1994年螺旋軸流式多相混輸泵被 FRAMO公司應(yīng)用于挪威石油公司所屬的Draugen平臺，使此油田的油氣水總產(chǎn)量提高了 50%左右，其中原油凈產(chǎn)增加了40%[7]。

油氣混輸泵在油氣集輸方面有著重要的作用，而動葉輪的優(yōu)劣又直接影響油氣混輸泵的揚程與輸送效率。隨著多相混輸技術(shù)的發(fā)展，我國以中國石油大學、蘭州理工大學、清華大學為首的高等院校開始了對油氣混輸泵的研究。李清平等[8]、朱宏武等[9]對螺旋軸流式多相泵進行了現(xiàn)場試驗研究與優(yōu)化設(shè)計，成功地進行了該泵的增壓運轉(zhuǎn)試驗。馬希金等[10-11]研究了軸向和徑向間隙對油氣混輸泵性能的影響。馬希金等[12]研究了葉片重疊度數(shù)對油氣混輸泵的影響。曲鑫等[13]研究了葉片數(shù)及葉片相對包角在不同含氣率下的流場分布情況。羅先武等[14]研究了葉片傾斜對微型泵水力性能的影響，得出適當傾斜一定的角度可以有效抑制壓力面二次流損失，從而提高水力性能。吳泳敏等[15]分析了直葉片、前傾葉片、后傾葉片對風扇氣動性能的影響。而葉片傾斜角對油氣混輸泵性能影響的研究較少。

本文以自主研發(fā)的YQH-100型三級軸流式油氣混輸泵為例，改變動葉傾斜角度，對油氣混輸泵內(nèi)部流動特性進行數(shù)值模擬，分析葉片壓力及氣相分布情況，并分析其揚程、增壓及效率曲線，為油氣混輸泵的設(shè)計及水力優(yōu)化提供參考依據(jù)。

1 計算模型的建立及網(wǎng)格劃分

1.1 數(shù)值計算模型

油氣混輸泵的基本參數(shù)：設(shè)計流量Q=100 m3/h，單級揚程H=33 m，轉(zhuǎn)速n=2 950 r/min，效率η=33%，泵進口管徑Din=100 mm，泵出口管徑Dout=65 mm，葉輪輪轂進口直徑D1=170 mm，葉輪輪轂出口直徑D2=180 mm，葉輪輪緣直徑D3=230 mm，輪轂半錐角θ=4.08°，葉片數(shù)z=4，葉輪輪轂比dh=0.782 6，輪緣翼型安放角β=10°，輪轂翼型安放角ɑ=12.68°。

根據(jù)模型的基本參數(shù)，用Pro/E軟件建立三維模型，YQH-100型油氣混輸泵的三維模型如圖1所示，主要由吸入室、動葉輪、長短導(dǎo)葉及壓出室等部件組成。為了得到較好的網(wǎng)格質(zhì)量，對隔舌、隔板等倒圓角，對一些細小的倒圓角等進行簡化得到了泵的流體域的三維造型。

圖1 全流域計算模型圖

1.2 葉片傾斜角度設(shè)計方案

葉片與輪轂垂直時葉片的傾斜角度為0°，葉片傾斜方向與葉輪旋轉(zhuǎn)方向一致為正，反之為負，并規(guī)定葉片正傾斜為前傾葉片，反之為后傾葉片。下面主要針對以下5種方案進行討論。方案一(M1)：動葉傾斜角度為-10°；方案二(M2)：動葉傾斜角度為-8°；方案三(M3)：動葉傾斜角度為-4°；方案四(M4)：動葉傾斜角度為0°即為原模型；方案五(M5)：動葉傾斜角度為4°。

1.3 網(wǎng)格劃分及網(wǎng)格無關(guān)性檢驗

將如圖1所示的模型導(dǎo)入專業(yè)的網(wǎng)格劃分軟件ICEM中進行網(wǎng)格劃分，創(chuàng)建相應(yīng)的body、part以及interface面。

劃分網(wǎng)格是數(shù)值模擬計算的一個重要環(huán)節(jié)，網(wǎng)格質(zhì)量好壞及數(shù)目多少將影響數(shù)值計算的準確度。網(wǎng)格數(shù)目越多，則越接近準確值，但網(wǎng)格數(shù)目越多，對計算機的要求越高。下面對其進行網(wǎng)格無關(guān)性驗證，綜合考慮各種因素，取網(wǎng)格數(shù)206萬比較合適。其模型計算網(wǎng)格如圖2所示。

圖2 模型計算域網(wǎng)格

2 控制方程及邊界條件的設(shè)置

2.1 控制方程

用FLUENT 軟件對內(nèi)部流場進行模擬，采用Mixture模型作為氣液兩相流模型，主相定義為水，輔相為空氣，泵內(nèi)部流動為不可壓縮的湍流流動，并將不可壓縮流動的連續(xù)方程與動量方程作為控制方程[16]。

連續(xù)方程：

(1)

動量方程：

(2)

式中：αi為各相的體積分數(shù)，氣液兩相流中即為液相與氣相的體積分數(shù)，并滿足αl+αg=1；ρi為各相密度；wi為各相相對速度；τ為黏性應(yīng)力張量；Mi為各相所受的相間作用力；fi為質(zhì)量力，有離心力與科氏力。

2.2 邊界條件設(shè)置

進口設(shè)置為速度進口，其值由入口流量與進口截面面積計算得出，清水時為自由出流(outflow)，不同的含氣率時為壓力出口，出口壓力由清水條件下的壓力計算得出。模擬系數(shù)均采取默認值并設(shè)置旋轉(zhuǎn)坐標系，用封閉的N-S方程進行離散，均采用一階迎風格式進行離散。用SIMPLEC算法對速度和壓力進行耦合，殘差值的收斂標準皆為10-4，模型采用標準k-ε模型。湍動能k和湍動能耗散率ε根據(jù)公式(3)計算。

(3)

式中，lin為進口處的混合長度，為進口直徑的0.015倍，固壁邊界均為無滑移固壁邊界條件，近壁區(qū)的流體流動均采用標準壁面函數(shù)進行處理。

3 計算結(jié)果及分析

3.1 葉片傾斜角對混輸泵外特性的影響

圖3為含氣率GVF=0時泵的揚程及效率隨流量的變化曲線。

(a)不同方案下的揚程曲線

(b)不同方案下的效率曲線

由圖3(a)可知M1～M5曲線的走勢基本一致，皆隨著流量的增大，揚程降低。在設(shè)計流量下，M1～M4揚程逐漸增大，M4～M5揚程逐漸減小。葉片后傾的角度越大，揚程降低得越明顯，揚程的極差值為8m，葉片前傾同樣會使揚程降低。由圖可知在設(shè)計流量左右，M4的揚程最高，M3、M5次之，M1的揚程最低。M1、M2曲線的彎曲程度較大，在0.8Qd～0.6Qd之間隨著流量的減小，揚程急劇上升。這是由于泵內(nèi)部產(chǎn)生的“二次回流”所引起的，二次回流會造成流動損失，在外特性上表現(xiàn)為揚程降低。

由圖3(b)可知各方案曲線的整體趨勢一致，泵的效率隨著流量的增大而升高，在Qd～1.2Qd之間，泵的效率相差不大，而在小流量工況下，M4的效率較好，M1、M2的效率較差。分析得出在純水情況下，M4比其他幾種方案要好。

圖4為不同方案下泵的增壓及效率隨含氣率的變化曲線。由圖4(a)可知M1～M5曲線皆隨著含氣率的增大而降低。 M3、M4與M5增壓相差不大，M4的增壓稍高于M3、M5。這說明M4的增壓效果較其他幾種方案要好，M1的增壓最低，其增壓效果最差。

(a)增壓隨含氣率的變化情況

(b)效率隨含氣率的變化情況

圖4(b)是各方案的效率隨含氣率的變化情況。隨著含氣率的升高，效率逐漸降低，含氣率從0升高到10%過程中，混輸泵的效率急劇下降，這主要是氣液兩相流的不均勻流動造成的。含氣率大于10%后效率趨于穩(wěn)定，減小的趨勢變緩。由圖4(b)可以看出M3、M5方案的效率皆高于M4，其中M3的效率最好，M2的效率最差，從效率方面考慮，M3方案優(yōu)于其他幾種方案。

3.2 混輸泵內(nèi)部壓力分布

圖5為GVF=0.5時第三級葉輪葉片工作面壓力云圖。圖中可以看出5種方案下壓力的分布趨勢基本一致。動葉工作面的壓力從輪轂到輪緣逐漸增大，沿著葉高方向徑向分布。以下幾種方案在葉片尾部1/3處靠近輪轂側(cè)皆存在明顯的高壓區(qū)，這是流體在此處發(fā)生紊流流動所引起的， 故形成了高壓力梯度區(qū) 。圖中可以看出M1方案中整體壓力較小，M4方案中的整體壓力較大且較為均勻。這說明M1的增壓效果最差，M4的增壓效果最好，M3、M5次之，與圖4(a)的增壓曲線的結(jié)果一致。但M4方案的局部高壓區(qū)較M3、M5要大，則在此處發(fā)生紊流流動造成的損失也越大，表現(xiàn)在外特性上即為效率的降低。

圖5 GVF=0.5時第三級葉輪葉片工作面壓力云圖

圖6為GVF=0.5時第三級葉輪葉片背面壓力云圖?？煽闯鋈~片背面頭部壓力分布不均勻，這是來流與葉片頭部碰撞形成的，致使葉片頭部壓力分布不均勻。由圖還可看出葉片背面壓力從頭部到尾部沿軸向逐漸增大，從M1到M4葉片背面壓力分布逐漸變得均勻，方案M3、M4的壓力分布較好，在葉片尾部壓力過渡平緩，壓力驟增現(xiàn)象減弱。

圖6 GVF=0.5時第三級葉輪葉片背面壓力云圖

3.3 混輸泵內(nèi)部氣相分布

圖7為GVF=0.5時第三級葉輪葉片工作面的氣相分布情況。由圖可看出含氣率沿著徑向逐漸減小。這是由于慣性離心力的作用，液體沿著輪緣分布，而氣體則向輪轂聚集，從而形成了輪轂側(cè)的含氣率高于輪緣側(cè)的氣相分布?？拷~片尾緣輪轂側(cè)存在明顯的高含氣區(qū)，這是由于大量的氣體聚集在輪轂側(cè)形成氣團。幾種方案氣團分布趨勢一致，葉片后傾可有效減小氣團在輪轂側(cè)的聚集面積，致使氣團徑向移動。這是由于葉片后傾，向心力減小，葉輪對流體的控制能力減弱，致使氣團沿著徑向移動，但傾斜角度過大，氣團徑向分布阻塞流道不利于流體的流動。M1、M2、M5靠近頭部輪緣側(cè)皆存在明顯的低含氣區(qū)，此處液相聚集，氣液分布不均，造成的流動損失較大。M3、M4氣相分布較為均勻，變化梯度較小，效果較好。

圖7 GVF=0.5時第三級葉輪葉片工作面氣相分布圖

圖8為GVF=0.5時第三級葉輪葉片背面的氣相分布情況。由圖可直觀看出葉片頭部存在明顯的高含氣區(qū)。這說明此處氣體大量聚集，并且可看出M1—M4葉片頭部的高含氣區(qū)隨著傾斜角度的增大而減小。M1方案1/2弦長靠近輪緣側(cè)存在局部的高含氣區(qū)，氣液分離現(xiàn)象明顯，水力損失較大，說明過大的傾斜角不利于氣液的均勻混合。M3方案的氣相分布優(yōu)于其他幾種方案。

4 結(jié)論

1)GVF=0時，葉片后傾的角度越大，揚程下降得越明顯，原模型的揚程最高，葉片后傾4°次之。隨著流量的增大，泵效率逐漸升高，在大流量工況下葉片傾斜角對泵效率影響不大。

2)隨著含氣率的升高混輸泵的增壓降低，葉片傾斜，增壓效果減弱，原模型的增壓效果最好。

3)隨著含氣率升高，混輸泵效率降低，葉片適當傾斜一定的角度，可提高混輸泵效率，葉片后傾4°的效率較好。

4)隨著葉片后傾的角度增大，葉片工作面輪轂側(cè)氣團面積減小，氣團沿著徑向移動，葉片后傾4°方案工作面的氣相分布較好。

圖8 GVF=0.5時第三級葉輪葉片背面氣相分布圖

5)隨著葉片后傾角度的增大，葉片背面頭部氣團面積減小，傾斜角度過大葉片輪緣側(cè)會形成高含氣區(qū)，葉片后傾4°方案的背面氣相分布更加均勻。

6)綜合考慮增壓、混輸泵效率以及葉片工作面背面的氣相分布，葉片的最佳傾斜角范圍為-4°～0°。