徐建寧 趙義鵬 李萬鐘 朱端銀
(西安石油大學機械工程學院)
徐建寧,趙義鵬,李萬鐘,等.密封式防積液短節(jié)設(shè)計及其攜液效果分析.石油機械,2022,50(11):119-125,132.
在雙碳目標下,天然氣成為替代煤炭的重要能源,其開發(fā)需求的提升也促進了輸氣管線的建設(shè)[1-2]。由于天然氣為含濕氣體,再加之輸氣管線沿地形起伏,便容易在低洼拐角處及上傾段形成積液[3-5],導(dǎo)致氣體的流通面積減小,進而造成管線內(nèi)壓力分布不均,降低了管線的輸運能力。渦流工具[6-14]最初用于氣井排水采氣領(lǐng)域,近些年也被投入到天然氣的地面輸運中,依靠其發(fā)生螺旋流而得以改變管線起伏段的液態(tài)分布。
但由于渦流工具的最大外徑需小于管線內(nèi)徑才能使其裝配在管線中,導(dǎo)致本體外側(cè)與管線內(nèi)壁之間存在間隙,部分氣體并沒有經(jīng)過螺旋流道便從間隙流走;再加之端面處流體碰撞引起的能量損失,使形成的螺旋流強度較低,沒有達到最佳的攜液效果;且現(xiàn)有關(guān)于渦流工具的研究也都是建立在非密封的基礎(chǔ)上[15-21]。鑒于此,筆者設(shè)計了一種密封式防積液短節(jié),分別建立了光滑管線、有渦流工具管線和有防積液短節(jié)管線的有限元模型,并使用計算流體力學的方法模擬研究3種模型流場速度、靜壓和含液體積分數(shù)的分布及變化規(guī)律。研究結(jié)果可為密封式防積液短節(jié)的開發(fā)及應(yīng)用提供一定的理論依據(jù)。
密封式防積液短節(jié)結(jié)構(gòu)如圖1所示,由導(dǎo)流錐、旋流體和套筒裝配而成。導(dǎo)流錐為圓錐狀實心體,在其頭部的柱狀體上有外螺紋,與旋流體尾部槽內(nèi)的內(nèi)螺紋連接,其功能是對流道入口處的流體形成導(dǎo)流,避免碰撞引起的能量損失。旋流體和套筒自帶連接法蘭,與輸氣管線上的法蘭連接;套筒的內(nèi)徑與輸氣管線的內(nèi)徑相同;旋流體外側(cè)有4條旋流翼,旋流翼側(cè)表面與套筒內(nèi)表面接觸時形成密封;未被密封的區(qū)域則形成了螺旋流道,流道入口和出口分別連通兩側(cè)的輸氣管線。密封式防積液短節(jié)的主要結(jié)構(gòu)參數(shù)如下:連接法蘭直徑為128 mm,法蘭孔內(nèi)徑為10 mm,套筒內(nèi)徑為88 mm,旋流體半徑為34 mm,旋流翼翼高為10 mm,旋流翼螺距為110 mm,導(dǎo)流錐長度為75 mm,導(dǎo)流錐角度為30°。
圖1 密封式防積液短節(jié)結(jié)構(gòu)圖Fig.1 Schematic structure of sealed anti-effusion short joint
密封式防積液短節(jié)的工藝原理如圖2所示。在管線起伏前的位置安裝短節(jié),通過其自身的連接法蘭與兩側(cè)的管線法蘭連接。含濕氣體輸運至短節(jié)處,在導(dǎo)流錐的引流下進入流道入口,隨后流經(jīng)螺旋流道,并逐漸引發(fā)螺旋流,兩相流體受到離心力的作用開始加速旋轉(zhuǎn)并流出出口。由于氣、液兩相密度的差異,含濕氣體分離形成了環(huán)狀液膜和柱狀氣柱,并以螺旋狀的形式向管線上傾段運動。這種流動方式使得氣液兩相的摩擦阻力減小,氣體攜液能力增強,可帶動液體輸運更遠的距離。
圖2 密封式防積液短節(jié)工藝原理圖Fig.2 Technical principle of sealed anti-effusion short joint
分別對光滑管線、有渦流工具的管線和有密封式防積液短節(jié)的管線建立流體域模型。模型管長為1 100 mm,直徑為88 mm,管線起伏角度為20°,水平高差為168 mm。在模型網(wǎng)格劃分時,對過流面積較小的區(qū)域進行網(wǎng)格加密,模型網(wǎng)格劃分結(jié)果如圖3所示。
圖3 流體域模型及網(wǎng)格劃分Fig.3 Fluid domain model and mesh division
開展無關(guān)性驗證,因管線拐角處的流場相比其他流域具有較為明顯的不同,故以此處前后的壓降大小作為網(wǎng)格劃分質(zhì)量的評定指標。分別以單個網(wǎng)格尺寸為10~4 mm的7種不同水平對3種模型進行網(wǎng)格劃分,驗證分析結(jié)果如圖4所示。
圖4 網(wǎng)格無關(guān)性驗證Fig.4 Grid independence verification
由圖4可以看出,當模型網(wǎng)格尺寸小于等于6 mm時,3種模型拐角處的前后壓差都能夠保持穩(wěn)定,證明此時的網(wǎng)格水平已經(jīng)滿足計算精度的要求,故3種模型均以網(wǎng)格尺寸6 mm的水平進行網(wǎng)格劃分。
多相流模型選用Euler模型,其對每一相求解動量方程和連續(xù)性方程,通過壓力和時間交換系數(shù)實現(xiàn)耦合,具有較高的精度。湍流模型選用Realizable k-ε模型,其對旋流、流動分離和強流線彎曲等特征有較強的適應(yīng)性,更加符合防積液短節(jié)的工作特點。Realizable k-ε模型相關(guān)的理論計算公式如下:
式中:Gk為層流速度梯度產(chǎn)生的湍動能項,Pa/s;Gb為浮力影響產(chǎn)生的湍動能項,Pa/s;YM為可壓縮湍流脈動膨脹對總的耗散率的影響,Pa/s;C1、C2、C1ε、G3ε為常量;σk、σε分別為k方程和ε方程相對應(yīng)的普朗特數(shù);u為速度,m/s;μ為動力黏度,Pa·s;μt為湍流黏度,Pa·s;ρ為流體密度,kg/m3;ν為運動黏度,m2/s;Sk、Sε為自定義的湍動能項(Pa/s)和湍流耗散源項(Pa/s2);k為湍流動能,m2/s2;ε為湍流耗散率,m2/s3。
設(shè)定模型左端為速度入口,入口速度為2 m/s;右端為壓力出口,出口壓力為5 MPa。設(shè)定甲烷氣體為第1相,體積分數(shù)為0.85;液態(tài)水為第2相,體積分數(shù)為0.15。甲烷的密度為0.668 kg/m3,黏度為0.013 4 mPa·s;水的密度為998.2 kg/m3,黏度為1.003 0 mPa·s。
另外,設(shè)定重力加速度為9.8 m/s2,其他流域表面為無滑移邊界。
氣液在經(jīng)過渦流工具和密封式防積液短節(jié)后的流動跡線如圖5所示。由圖5可以看出:流體流型都由直流轉(zhuǎn)變?yōu)榱谁h(huán)狀螺旋流,但在渦流工具管線中,流體與工具端面發(fā)生碰撞,能量遭到損失,且一部分流體沒有進入螺旋流道便直接從間隙中流走;而在有密封式防積液短節(jié)管線中,導(dǎo)流錐的導(dǎo)流作用使流體流動平滑,減少了碰撞,螺旋流道的密封作用使流動的能量被充分利用,減少了間隙流失。因此防積液短節(jié)在管線中引發(fā)的螺旋流強度明顯高于渦流工具引發(fā)的螺旋流強度。
圖5 流動跡線圖Fig.5 Flow trace diagram
3種模型的軸向速度分布如圖6所示。由圖6可以看出:受到重力及管線起伏的影響,流場分布不均,光滑管線拐角處下方位的軸向速度遠小于上方位,導(dǎo)致氣液在此處動力不足,容易形成積液;而經(jīng)過渦流工具和防積液短節(jié)后流場呈現(xiàn)均勻的環(huán)狀分布,且軸向速度得到了提升,后者產(chǎn)生的軸向速度大于前者。
圖6 軸向速度分布云圖Fig.6 Cloud chart of axial velocity distribution
在管線模型的不同軸向位置截取平面,觀測截面各徑向位置的軸向速度,得到3種模型的截面軸向速度變化,如圖7所示。由圖7可以看出:光滑管線水平段(z<400 mm)速度維持不變,拐角處(z=400 mm)和上傾段(z>400 mm)上方位速度大于下方位速度,不同截面的軸向速度大小基本一致;而在經(jīng)過渦流工具和防積液短節(jié)后,盡管在上傾段初始有反向的軸向速度,但隨著軸向位置的遞進,反向速度逐漸消失,形成正向的軸向速度,并呈現(xiàn)遞增趨勢;在相同軸向位置,后者的軸向速度大于前者,證明密封式防積液短節(jié)比渦流工具產(chǎn)生了更高的軸向速度,有利于增加攜液距離。3種模型的切向速度分布云圖如圖8所示。由圖8可以看出:在光滑管線水平段基本不具備切向速度,從拐角處至上傾段逐漸具備一定的切向速度,但不足以使流型發(fā)生轉(zhuǎn)變;而經(jīng)過渦流工具和防積液短節(jié)后,切向速度得到了提升,后者產(chǎn)生的切向速度大于前者。
圖7 不同軸向位置的軸向速度對比Fig.7 Comparison of axial velocity at different axial positions
圖8 切向速度分布云圖Fig.8 Cloud chart of tangential velocity distribution
3種模型的截面切向速度變化如圖9所示。由圖9可以看出:光滑管線水平段(z<400 mm)基本不具備切向速度,拐角處(z=400 mm)切向速度出現(xiàn)了極差,上傾段(z>400 mm)不同軸向位置截面的切向速度大小基本一致;而在經(jīng)過渦流工具和防積液短節(jié)之后,切向速度明顯提高,并且在拐角處(z=400 mm)的切向速度為最大值,證明此時氣液的流型已經(jīng)具備發(fā)生轉(zhuǎn)變的條件,隨著軸向位置的遞進,截面的切向速度呈現(xiàn)遞減趨勢;在相同軸向位置,后者的切向速度大于前者,證明防積液短節(jié)比渦流工具產(chǎn)生了更高的切向速度,因此旋流強度更高,有利于轉(zhuǎn)變流體的流型。
圖9 不同軸向位置的切向速度對比Fig.9 Comparison of tangential velocity at different axial positions
3種模型的靜壓分布云圖如圖10所示。
由圖10可以看出:光滑管線整體靜壓分布不均,上方位的靜壓小于下方位靜壓,拐角處和上傾段有明顯的波動;而經(jīng)過渦流工具和防積液短節(jié)后,由于過流面積的驟減,導(dǎo)致過流后的靜壓減小,盡管防積液短節(jié)所在管線的壓降值相比光滑管線的壓降值和渦流工具管線的壓降值增大,但靜壓沿線分布均勻,總體無較大波動,因此不會影響管線的正常輸運。
圖10 靜壓分布云圖Fig.10 Cloud chart of static pressure distribution
3種模型的截面靜壓變化如圖11所示。由圖11可以看出,光滑管線不同截面的靜壓沿著徑向位置發(fā)生變化,上方位的靜壓小于下方位,拐角處(z=400 mm)靜壓變化梯度最大,上傾段(z>400 mm)變化梯度較小,隨著軸向位置的遞進,靜壓呈現(xiàn)緩慢減小的趨勢;而在經(jīng)過渦流工具和防積液短節(jié)之后管線雖急劇產(chǎn)生了壓降,但不同截面靜壓分布都較為均勻,隨著軸向位置的遞進,截面靜壓無較大波動,因此流體在上傾段流動平穩(wěn),證明防積液短節(jié)能起到穩(wěn)流降壓的作用。
圖11 不同軸向位置的靜壓對比Fig.11 Comparison of static pressure at different axial positions
3種模型的含液體積分數(shù)軸向分布云圖如圖12所示。
由圖12可以看出:光滑管線水平段積液并不明顯,但在拐角處和上傾段的下方位由于受到重力及流動速度較弱的影響,形成了積液,氣體則主要分布在上方位;而氣液經(jīng)過渦流工具和防積液短節(jié)后,氣液兩相的流型發(fā)生了轉(zhuǎn)變,液相流動為環(huán)狀流,氣相則為柱狀分布在管線中心,但可以看出前者的氣柱周圍仍有部分霧狀流液體,而后者的環(huán)狀流更加凝聚,氣柱體積也明顯大于前者,形成了更加明顯的氣液分界面。
圖12 含液體積分數(shù)軸向分布云圖Fig.12 Cloud chart of liquid volume fraction distribution
3種模型的截面含液體積分數(shù)變化如圖13所示。由圖13可以看出:光滑管線拐角處(z=400 mm)和上傾段(z>400 mm)不同軸向位置截面上方位的含液體積分數(shù)小于下方位;而在經(jīng)過渦流工具和防積液短節(jié)之后,截面液相的分布情況轉(zhuǎn)變?yōu)橹行牡?、周向?在有渦流工具管線中,不同截面的含液體積分數(shù)隨著軸向位置的遞進呈現(xiàn)先增后減的趨勢,說明部分液相產(chǎn)生了回流;而在有防積液短節(jié)的管線中,截面含液體積分數(shù)隨著軸向位置的遞進呈現(xiàn)遞增趨勢,說明液體持續(xù)向管線上游流動,防積液短節(jié)的攜液距離更遠;在周向位置,后者的含液體積分數(shù)高于前者,而在中心位置,后者的含液體積分數(shù)低于前者,證明防積液短節(jié)對流型的轉(zhuǎn)變比渦流工具更徹底,總體的攜液效果更佳。
圖13 不同軸向位置含液體積分數(shù)對比Fig.13 Comparison of liquid volume fraction at different axial positions
(1)密封式防積液短節(jié)的導(dǎo)流錐和螺旋流道起到了導(dǎo)流作用及密封作用,有效減少了端面流體碰撞和間隙流失,其引發(fā)的螺旋流強度高于渦流工具引發(fā)的螺旋流強度。
(2)密封式防積液短節(jié)在管線中產(chǎn)生的流體軸向速度和切向速度相比渦流工具得到了提升。
(3)密封式防積液短節(jié)能夠有效改善管線靜壓分布不均的現(xiàn)象,起到了穩(wěn)流降壓的作用。
(4)經(jīng)過渦流工具后液相仍有回流,而經(jīng)過密封式防積液短節(jié)后,管線氣液分界面更加明顯,攜液距離增加,流型轉(zhuǎn)變更徹底。