為了探索新型轉子式油氣混輸泵出口球閥內流場(chǎng)規律,建立球閥流場(chǎng)的三維模型,利用Fluent軟件,將標準k-ε湍流模型與多相流技術(shù)相結合,采用SIMPLE 算法,對新型轉子式油氣混輸泵出口球閥內的三維氣液兩相流場(chǎng)進(jìn)行數值模擬。在容積含氣率為25% ,50% ,75% 的不同工況下,通過(guò)對球閥開(kāi)啟高度分別為3,5,7mm 時(shí)的速度場(chǎng)、壓力場(chǎng)與氣液相分布的分析,探討在氣液混輸過(guò)程中閥的開(kāi)啟高度及不同氣液比對閥內流場(chǎng)的影響規律。模擬結果表明:球閥開(kāi)啟高度越大,閥球上下壓差越??;閥隙流速隨著(zhù)開(kāi)啟高度的增大而減小。
球閥具有結構簡(jiǎn)單、互換性強、裝拆方便、便于清洗等優(yōu)點(diǎn)。為解決油田油氣混輸難題,將球閥與傳統外環(huán)流轉子泵結合,即在傳統外環(huán)流轉子泵出口增設了1組球閥,使其具有內壓縮功能,能更好地適應氣液兩相工況。目前,對于球閥的研究基本上是針對容積式往復泵球閥,主要建立球閥運動(dòng)規律的數學(xué)模型,研究球閥的開(kāi)啟特性等內容,且工況為純液態(tài)工況;對球閥閥口氣穴流場(chǎng)進(jìn)行的數值模擬與試驗研究也局限于液體介質(zhì)。 目前尚未見(jiàn)有關(guān)轉子式油氣混輸泵球閥運動(dòng)規律的研究報道。因此,對新型轉子式油氣混輸泵出口球閥的研究就顯得很有必要。
隨著(zhù)計算機技術(shù)和計算流體力學(xué)的發(fā)展,應用CFD方法對流場(chǎng)進(jìn)行分析已經(jīng)成為泵閥領(lǐng)域的研究熱點(diǎn)。因此,實(shí)現新型轉子式油氣混輸泵出口球閥三維流場(chǎng)的數值模擬,對于球閥的設計及優(yōu)化具有重要意義。
1 球閥結構及網(wǎng)格劃分
1.1 球閥結構
圖1為轉子式油氣混輸泵工作示意圖。新型轉子式油氣混輸泵在出口增設球閥以后,介質(zhì)要通過(guò)球閥才能輸送到出口管線(xiàn)中。當球閥關(guān)閉時(shí),閥球與兩轉子及端板形成封閉容積V。由于轉子不斷旋轉,封閉容積V不斷減小,容積中壓力不斷升高,直到封閉容積內的壓力達到開(kāi)啟壓力時(shí),閥球打開(kāi),介質(zhì)被排出。
圖1 轉子式油氣混輸泵工作示意圖
圖2為出口球閥結構示意圖。球閥由閥座和閥球組成,閥球開(kāi)啟后,介質(zhì)由閥座孔入口流入,通過(guò)閥隙進(jìn)入泵的排液腔。閥座孔直徑d=0.065m,閥座錐角α=45°,錐角長(cháng)度l=0.005m,閥球半徑R=0.045m。
圖2 出口球閥結構示意圖
1.2 建模與網(wǎng)格劃分
由于出口球閥尺寸相對整臺泵非常小,在整臺泵計算過(guò)程中,難以得到閥隙處的詳細流動(dòng)情況。因此,為了更全面地了解閥隙周?chē)c閥內的壓力和速度分布,選取閥座與閥球間隙及閥球兩側部分作為研究對象,進(jìn)行建模與分析。此外,球閥幾何形狀簡(jiǎn)單且為軸對稱(chēng)圖形,為了研究方便且減少計算量,采取三維軸對稱(chēng)模型,建立一半計算區域。利用Pro/E軟件建立開(kāi)啟高度為3mm時(shí)球閥的計算區域模型。將物理模型導入Fluent前處理軟件Gambit中進(jìn)行網(wǎng)格劃分。為了劃分質(zhì)量較好的網(wǎng)格,對模型進(jìn)行了分割并采用六面體/四面體混合單元,由于閥口的壓力梯度變化較大,因此對閥口加密了網(wǎng)格,使模擬結果更準確。三維模型及網(wǎng)格如圖3所示。同理可得到開(kāi)啟高度分別為5,7mm時(shí)的模型和網(wǎng)格。
圖3 開(kāi)啟高度為3mm時(shí)的模型與網(wǎng)格
2 模擬計算
2.1 邊界條件
介質(zhì)為原油和天然氣兩相混合物,原油的物理參數設置為ρoil=856kg/m3,動(dòng)力黏度ν=0.0072Pa·s,并假設原油不可壓縮;天然氣在Fluent自帶的材料里選擇。
1) 速度入口。新型轉子式油氣混輸泵的出口閥由3個(gè)球閥組成,已知泵的流量為100m3/h,假設通過(guò)每個(gè)球閥的流量相等且忽略泄漏,則由連續流條件可得通過(guò)每個(gè)閥座的速度為
(1)
式中:υ0為入口速度,m/s,方向與閥座入口邊垂直;Q為泵的流量,m3/h;d為閥座孔直徑,m。由入口速度和特征直徑計算得到入口雷諾數大于1.2×104,流動(dòng)為湍流,湍流強度設為10%,水力直徑為0.065m。
2) 壓力出口。已知出口絕對壓力為1.2MPa。
2.2 求解器與算法
模擬采用隱式壓力基求解器,流動(dòng)為穩態(tài)流動(dòng). 求解模型選擇兩相混合模型和標準k-ε湍流模型。壓力與速度耦合采用SIMPLE算法。
3 模擬結果與分析
3.1 壓力場(chǎng)分析
圖4-6為不同含氣率條件下,球閥在不同開(kāi)啟高度時(shí)對稱(chēng)面上的壓力分布云圖。
圖4 含氣率為25%時(shí)對稱(chēng)面上的靜壓分布
圖5 含氣率為50%時(shí)對稱(chēng)面上的靜壓分布
圖6 含氣率為75%時(shí)對稱(chēng)面上的靜壓分布
由壓力云圖可得,當開(kāi)啟高度為3mm,含氣率分別為25%,50% ,75% 時(shí),閥球上下壓差分別為0.06,0.04,0.02MPa;當開(kāi)啟高度為5mm,含氣率分別為25%,50% ,75% 時(shí),閥球上下壓差分別為0.04,0.02,0.01MPa;當開(kāi)啟高度為7mm,含氣率分別為25%,50% ,75% 時(shí),閥球上下壓差分別為0.02,0.01,0.01MPa. 以上分析表明:
1) 在同一含氣率的條件下,隨著(zhù)開(kāi)啟高度的增大,閥球上下壓差逐漸減小。
2) 在某一較小的固定開(kāi)啟高度時(shí),閥球上下壓差隨含氣率增大而減??; 開(kāi)啟高度較大時(shí),含氣率對閥球上下壓差影響較小。
3) 含氣率大時(shí),閥球上下壓差較小且受開(kāi)啟高度的影響較小。
4) 在球閥的整個(gè)流場(chǎng)中,閥隙處的壓強最小。
3.2 速度分析
圖7-9為不同含氣率和開(kāi)啟高度下流場(chǎng)Y-Z截面上的速度云圖和流線(xiàn)圖。
圖7 含氣率為25%時(shí)不同開(kāi)啟高度的速度云圖和流線(xiàn)圖
圖8 含氣率為50%時(shí)不同開(kāi)啟高度的速度云圖和流線(xiàn)圖
圖9 含氣率為75%時(shí)不同開(kāi)啟高度的速度云圖和流線(xiàn)圖
由速度云圖可知:氣液比一定時(shí),由于過(guò)流斷面突然減小,閥隙處的流速最大。隨著(zhù)開(kāi)啟高度的增大,閥隙流速不斷減小。圖7中,開(kāi)啟高度為5mm時(shí),閥隙流速為10m/s;開(kāi)啟高度為7mm時(shí)閥隙流速只有7m/s。
由流線(xiàn)圖可知,在閥隙附近有部分介質(zhì)由于壓差的作用回流,之后被閥隙的高速介質(zhì)帶出。例如圖7a中的流線(xiàn)所示,部分介質(zhì)從出口回流,但在閥隙附近流線(xiàn)方向突然改變,與從閥隙流出的介質(zhì)一起沿著(zhù)閥球壁附近流出。
另外,當開(kāi)啟高度為3mm時(shí),含氣率分別為25% ,50% ,75% 對應的閥隙流速均為15m/s。由上可知,同一開(kāi)啟高度下,含氣率對閥隙流速的影響不大。但同一開(kāi)啟高度下不同含氣率的流線(xiàn)不同,如開(kāi)啟高度為3mm時(shí),含氣率為75%的流線(xiàn)圖出現交叉流線(xiàn),不同于另外2種開(kāi)啟高度的流線(xiàn),說(shuō)明含氣率對介質(zhì)的流動(dòng)狀態(tài)有一定的影響。
3.3 相態(tài)分布分析
圖10為開(kāi)啟高度為3mm時(shí),不同含氣率的氣相體積分數分布云圖。
由圖10可知,氣相主要分布在閥球壁附近,遠離閥球的氣相介質(zhì)逐漸減少。通過(guò)模擬結果可知,氣相介質(zhì)密度較小,在閥球開(kāi)啟前,閥球底部分布的主要為氣體,球閥開(kāi)啟后,氣體介質(zhì)首先排出。這表明,氣液兩相分界較為明顯,有利于氣相介質(zhì)的單獨回收。
圖10 不同含氣率時(shí)氣相體積分數分布云圖
4 結論
1) 在含氣率一定的條件下,隨著(zhù)開(kāi)啟高度的增大,閥球上下壓差逐漸減??; 在球閥的整個(gè)流場(chǎng)中,閥球底部壓力最大; 閥隙處壓力梯度大,閥座倒角下端處較容易產(chǎn)生氣蝕。
2)在含氣率一定時(shí),閥隙流速隨著(zhù)開(kāi)啟高度的增大不斷減小。 部分介質(zhì)由于壓差的作用回流,之后被閥隙的高速介質(zhì)帶出。同一開(kāi)啟高度下,含氣率對閥隙速度的影響不大,但對流動(dòng)狀態(tài)有一定的影響。
3)新型轉子式油氣混輸泵在輸送氣液兩相介質(zhì)時(shí),氣相主要分布在閥球壁附近,遠離閥球氣相介質(zhì)逐漸減少。
經(jīng)典應用
產(chǎn)品中心