王小兵 王多琦 李森 龔浩宇 劉陽(yáng) 呂雷綱 杜明智
(常州大學(xué)石油工程學(xué)院 江蘇常州 213000)
抽油泵是有桿抽油系統(tǒng)重要的組成部分,提高抽油泵泵效對(duì)于提高采油效率至關(guān)重要,而抽油泵柱塞與泵筒之間的間隙漏失是影響抽油泵效率的主要因素之一。國(guó)內(nèi)外學(xué)者對(duì)于抽油泵柱塞和泵筒之間的間隙漏失率的計(jì)算及影響因素的研究做了大量工作。研究初期國(guó)外學(xué)者通過(guò)理論的方法提出不同的漏失模型[1-3],對(duì)于漏失率的計(jì)算主要是公式的推導(dǎo),并對(duì)各項(xiàng)系數(shù)進(jìn)行修訂[4-7]。這種方法不能全面分析漏失率影響因素[8],對(duì)反映真實(shí)復(fù)雜井底條件對(duì)漏失率的影響有限。后期隨著CFD軟件的飛速發(fā)展,數(shù)值計(jì)算方法越來(lái)越多的被運(yùn)用到復(fù)雜流體研究領(lǐng)域[9],學(xué)者們利用有限元分析的方法對(duì)抽油泵間隙液體的流體流動(dòng)情況有了越來(lái)越詳細(xì)的了解[10-12]。
本文通過(guò)參考抽油泵標(biāo)準(zhǔn)基本參數(shù),根據(jù)井底有桿泵工作特點(diǎn),建立起在多種邊界條件下的多物理場(chǎng)耦合模型,通過(guò)數(shù)值計(jì)算結(jié)果分析得出在不同沖次、沖程、進(jìn)出口壓差、間隙寬度下間隙流體的速度場(chǎng)分布及研究其在不同工作制度下的漏失率影響因素。同時(shí)通過(guò)改變開(kāi)采流體物性參數(shù)來(lái)研究不同含水率及粘度流體在相同工況下的漏失率的變化。這對(duì)優(yōu)化抽油泵結(jié)構(gòu)、減少抽油泵漏失率、提高抽油泵效率有重要的指導(dǎo)和參考意義[13]。
為使得計(jì)算結(jié)果盡量與實(shí)際結(jié)果的誤差最小,本文根據(jù)井下抽油系統(tǒng)結(jié)構(gòu)特點(diǎn),建立起與現(xiàn)場(chǎng)抽油結(jié)構(gòu)相同的1∶1簡(jiǎn)化幾何模型。由于幾何模型呈軸對(duì)稱(chēng),所以可以利用模擬軟件特有的軸對(duì)稱(chēng)建模方式進(jìn)行幾何建模,并充分考慮采出液的流體物性參數(shù),利用模擬軟件對(duì)模型進(jìn)行流體和固體耦合。通過(guò)對(duì)軟件中所構(gòu)建的幾何模型,進(jìn)行網(wǎng)格劃分,邊界條件的設(shè)定,對(duì)模型進(jìn)行相對(duì)應(yīng)的數(shù)值計(jì)算,最后對(duì)模擬計(jì)算結(jié)果進(jìn)行可視化后處理,對(duì)相應(yīng)的結(jié)果進(jìn)行圖表和曲線(xiàn)的繪制與描述,并具體分析其對(duì)漏失率影響機(jī)理。
本文為簡(jiǎn)化計(jì)算的復(fù)雜性,便于構(gòu)建計(jì)算模型,做出如下假設(shè):
(1)將泵內(nèi)流體看作為只有油水兩相均勻混合的均相流體,不含氣體,不考慮原油中其他沙石等雜質(zhì)對(duì)流體流動(dòng)產(chǎn)生的影響。
(2)由于泵的柱塞與泵筒間間隙漏失只發(fā)生在上沖程階段,此時(shí)游動(dòng)閥球視作全程處于關(guān)閉狀態(tài),而固定閥球處于打開(kāi)狀態(tài),不考慮閥球滯后所帶來(lái)的進(jìn)出口壓差的變化。
(3)在柱塞運(yùn)動(dòng)過(guò)程中,僅考慮柱塞在沿泵筒方向的運(yùn)動(dòng),不考慮由于柱塞在徑向方向上受力不均勻所造成的柱塞在徑向方向上的偏移。
(4)由于柱塞的運(yùn)動(dòng)形式幾乎按照正弦運(yùn)動(dòng)規(guī)律,在柱塞運(yùn)動(dòng)時(shí),由于粘度、壓力等因素的影響,大部分流體的流動(dòng)狀態(tài)處于湍流。由于流體在間隙內(nèi)呈剪切流形式,因此在近壁面處流動(dòng)速度較低,流動(dòng)處于非湍流,因此,在進(jìn)行邊界層設(shè)定時(shí)將流體流動(dòng)設(shè)定為以湍流的形式進(jìn)行流動(dòng)。
(5)不考慮溫度對(duì)流體產(chǎn)生的影響,在泵工作時(shí),溫度的變化很小,幾乎對(duì)泵內(nèi)流體產(chǎn)生的影響很小,因此所涉及到的邊界均為絕熱邊界。
為保障計(jì)算結(jié)果的可靠性,建立起與實(shí)際抽油泵結(jié)構(gòu)參數(shù)相同的簡(jiǎn)化幾何模型。由于抽油泵柱塞和泵筒間間隙漏失只發(fā)生在柱塞上沖程過(guò)程中,此時(shí)游動(dòng)凡爾在壓差的作用下處于關(guān)閉狀態(tài),因此在建立模擬模型時(shí)簡(jiǎn)化并省略此部分以此來(lái)提高計(jì)算效率,抽油泵泵筒與柱塞間隙簡(jiǎn)化模型示意圖如圖1所示。
柱塞運(yùn)動(dòng)方程:
(1)
式中,v為柱塞運(yùn)動(dòng)速度,m/s;S為沖程,m;T為周期,s。
流體的漏失率為:
(2)
A1=π(R1/2)2
(3)
A2=π(R2/2)2
(4)
式中,Q為漏失率,m3/min;n為沖次;T1為上死點(diǎn)時(shí)刻,s;T2為下死點(diǎn)時(shí)刻,s;R1為泵筒內(nèi)徑,mm;R2為柱塞外徑,mm。
湍流是流體在時(shí)間和空間上的一種非線(xiàn)性的流體運(yùn)行形勢(shì),而泵內(nèi)流體除了在近壁面處由于粘度和壓力的影響處于層流,其他大多流體處于湍流狀態(tài),因此本文流體的流動(dòng)用湍流以及湍流能耗散率運(yùn)輸方程k-ε來(lái)描述流體流動(dòng)。
(5)
ρ▽·(u)=0
(6)
K=(μ+μT)(▽u+(▽u)T)
(7)
(8)
(9)
ε=ep
(10)
(11)
pk=μT[▽u:(▽u+(▽u)2)]
(12)
式中,ρ為流體密度,kg/m3;u為流速,m/s;μ為粘度,mPa·s;g為重力,N/s2;F為質(zhì)量力,N;l為管長(zhǎng),m;▽為拉普拉斯算子;p為壓力,Pa;ε為湍流耗散率,%;k為湍流動(dòng)能,J;Cμ為常數(shù),常取0.09。
流體在流動(dòng)過(guò)程中都滿(mǎn)足連續(xù)性方程與動(dòng)量守恒方程:
連續(xù)性方程:
(13)
式中,xi為在i方向的分量;uj為在j方向上的速度分量。
動(dòng)量守恒方程:
(14)
式中,τij為剪切應(yīng)力張量;ρgi為i方向上的體積力;Fi為源項(xiàng)。
由于實(shí)際抽油泵幾何的長(zhǎng)度與直徑比較大,通過(guò)軟件自動(dòng)劃分的網(wǎng)格的質(zhì)量較差,且占用計(jì)算機(jī)內(nèi)存較大,使得計(jì)算的時(shí)間較大,所以通過(guò)手動(dòng)劃分網(wǎng)格可以得到較好的網(wǎng)格質(zhì)量,而抽油泵柱塞與間隙的縫隙幾何結(jié)構(gòu)較小,所以對(duì)此部分進(jìn)行加密處理,可以使得運(yùn)算較高的收斂性,流體域整體的網(wǎng)格密度大于固體域網(wǎng)格密度,這樣在保證精確解的同時(shí),也減小了計(jì)算量,提高了計(jì)算效率,部分網(wǎng)格劃分示意圖如圖2所示。
圖2 部分網(wǎng)格劃分示意
通過(guò)對(duì)幾何模型最復(fù)雜處的網(wǎng)格質(zhì)量信息進(jìn)行統(tǒng)計(jì)可以得出表1。
表1 部分模型網(wǎng)格質(zhì)量統(tǒng)計(jì)信息
通過(guò)表1可以發(fā)現(xiàn),最小網(wǎng)格質(zhì)量和平均網(wǎng)格質(zhì)量均達(dá)到了1.0,這說(shuō)明網(wǎng)格的劃分已達(dá)到完美網(wǎng)格級(jí)別,這使得計(jì)算結(jié)果具有較高的收斂性。
本文所研究的是在瞬態(tài)情況下抽油泵工作時(shí)的流固耦合問(wèn)題,幾何參數(shù)根據(jù)φ70 mm抽油泵進(jìn)行設(shè)定,計(jì)算所需的其他參數(shù),依據(jù)實(shí)際邊界參數(shù)設(shè)定進(jìn)行相應(yīng)的調(diào)整。由于影響抽油泵泵筒和柱塞間間隙漏失率的影響因素較多,為分析各個(gè)影響因素對(duì)漏失率的影響,通過(guò)控制其他變量,觀(guān)察單一變量的影響。
2.1不同間隙寬度對(duì)抽油泵泵筒和間隙間漏失率的影響
由于柱塞泵間隙大小的變化對(duì)其間的漏失率影響較為直接,所以根據(jù)表2中5個(gè)不同的間隙等級(jí)[14],來(lái)各取5個(gè)間隙值大小(0.03,0.07,0.11,0.15,0.18 mm),對(duì)比其各自的速度場(chǎng)及漏失率的大小的變化情況,以此來(lái)分析間隙大小對(duì)柱塞泵漏失機(jī)理的影響。
表2 柱塞泵間隙標(biāo)準(zhǔn)等級(jí) mm
為分析柱塞泵間隙大小這一影響因素對(duì)抽油泵漏失率的影響,設(shè)定在相同的抽油泵工作制度下(沖次為6次/min,沖程為1.8 m)及相同的流體物性參數(shù)下,抽油泵柱塞與泵筒間隙內(nèi)流體的流動(dòng)情況,在上沖程過(guò)程中,截取3個(gè)時(shí)刻(0.5,2.5,7.5 s),柱塞附近的液體流動(dòng)速度變化圖即可得出圖3。
(a)0.5 s
如圖3所示,可以發(fā)現(xiàn)當(dāng)柱塞運(yùn)動(dòng)到不同位置時(shí),流體的速度分布有明顯的不同,且變化較大,但同時(shí)由于間隙較小,從整體上無(wú)法觀(guān)察到里面的流體流場(chǎng)的變化,因此將間隙做了局部放大處理,得出圖4。
(a)0.5 s
如圖4所示,是5個(gè)不同時(shí)刻間隙內(nèi)的流速分布圖,圖4(a)是指在柱塞運(yùn)動(dòng)的一個(gè)沖程過(guò)程中,處于柱塞位于下死點(diǎn)位置的時(shí)刻,此時(shí)流體流速在垂直軸向方向上是分布均勻的,而在徑向方向上存在明顯的分層現(xiàn)象,這是由于間隙內(nèi)液體具有粘性,會(huì)對(duì)泵筒壁及柱塞表面產(chǎn)生粘性力的作用,而在這種粘性力的作用下,使得間隙內(nèi)的流體以剪切流的形式進(jìn)行流動(dòng)。而隨著間隙流道從外側(cè)向中間發(fā)展時(shí),由于此時(shí)越往流道中間的液體受到泵筒壁及柱塞表面對(duì)其產(chǎn)生的粘性力影響越小,所以間隙中間的流速明顯是最大的。隨著柱塞從下死點(diǎn)到上死點(diǎn)運(yùn)動(dòng)過(guò)程中,柱塞間隙內(nèi)流體與柱塞的相對(duì)運(yùn)動(dòng)速度逐漸增大,而粘性力對(duì)間隙內(nèi)流體流速影響逐漸減小。如圖4(b)、圖4(c)、圖4(d)所示,1.5,2.5,3.5 s這3個(gè)時(shí)刻間隙內(nèi)流體剪切流的形式[15],并不像0 s和5 s時(shí)刻處這么明顯,因?yàn)檫@兩個(gè)時(shí)刻柱塞處于下死點(diǎn)和上死點(diǎn)位置,流體與柱塞的相對(duì)運(yùn)動(dòng)速度最小。
通過(guò)測(cè)得5個(gè)不同的間隙大小(0.03,0.07,0.11,0.15,0.18 mm)在4個(gè)時(shí)刻點(diǎn)(1.5,2.5,3.5,4.5 s)時(shí)的流速分布可以得到,如圖5所示不同時(shí)刻泵筒與柱塞間隙內(nèi)流體流速分布曲線(xiàn)。
(a) 1.5 s時(shí)不同間隙內(nèi)流體流速分布曲線(xiàn)
由圖5(a)、圖5(b)、圖5(c)、圖5(d)可知,不同間隙寬度對(duì)其內(nèi)流體的流速影響較為明顯,4個(gè)不同時(shí)刻的流速分布圖比較可以發(fā)現(xiàn),隨著柱塞運(yùn)動(dòng)速度的增大,不同寬度的間隙內(nèi)流體的流速差異逐漸減小,這也同時(shí)說(shuō)明液體粘度對(duì)速度影響較大。在間隙入口處曲線(xiàn)產(chǎn)生不規(guī)則上下波動(dòng),說(shuō)明入口處的速度較大波動(dòng),進(jìn)入間隙內(nèi)的流體經(jīng)過(guò)充分發(fā)展流體流速在整個(gè)間隙速度達(dá)到穩(wěn)定,其在整個(gè)間隙內(nèi)分布均勻。
為直觀(guān)的了解不同間隙大小對(duì)漏失率的影響,通過(guò)計(jì)算在上沖程過(guò)程中,5個(gè)不同間隙值下的漏失率的大小,并可以得出如圖6所示的漏失率隨間隙變化的關(guān)系曲線(xiàn)。
圖6 不同間隙漏失率曲線(xiàn)
圖6所示是在抽油泵沖次為6次/min,沖程為2.4 m,進(jìn)出口壓差為12.5 MPa,流體物性參數(shù)一致的情況下根據(jù)抽油泵柱塞與泵筒不同間隙寬度下的漏失率變化曲線(xiàn),由曲線(xiàn)的變化可以看出,間隙的漏失率幾乎隨間隙寬度呈線(xiàn)性變化趨勢(shì),而且這說(shuō)明間隙寬度是影響其漏失率的重要影響因素之一。
2.2不同沖程對(duì)抽油泵泵筒和間隙間漏失率的影響
為研究不同沖程對(duì)抽油泵間隙漏失率影響,設(shè)定在其他外界條件一定的情況下,通過(guò)計(jì)算在相同的沖次(6次/min),分布建立了8個(gè)不同沖程(1,1.5,2,2.5,3,3.5,4,4.5 m)的抽油泵柱塞與泵筒間隙間的漏失幾何模型,并得出其漏失每分鐘內(nèi)漏失率,如圖7所示。
圖7 不同沖程漏失率曲線(xiàn)
由圖7可知,漏失率隨沖程的變化趨勢(shì)幾乎呈線(xiàn)性變化,而且曲線(xiàn)斜率較大,這說(shuō)明較小的沖程變化都會(huì)引起較大的漏失率的變化,這是由于在相同時(shí)間內(nèi),在沖次不變的情況下,增加沖程會(huì)使得柱塞運(yùn)動(dòng)速度加快。增加柱塞沖程一方面增大了漏失的距離,另一方面加快了柱塞的運(yùn)動(dòng)速度,兩者同時(shí)作用下,都會(huì)使得抽油泵泵筒和柱塞間的漏失率增大。
2.3不同沖次對(duì)抽油泵泵筒和柱塞間隙間漏失率的影響
與以上研究方法相同,在研究抽油泵沖次對(duì)漏失率的影響時(shí),通過(guò)不改變其他外界條件和流體物性參數(shù),分別計(jì)算在不同沖次(4次、5次、6次、7次、8次)下漏失率的大小,并繪制出如圖8所示,漏失率隨沖次變化的曲線(xiàn)。
由圖8可知,在一個(gè)周期過(guò)程中,漏失率隨著沖次的增加而增加,但是兩者不呈線(xiàn)性關(guān)系,并且隨著沖次的增加,抽油泵漏失率增加的趨勢(shì)逐漸放緩。這是由于抽油泵沖次的變化會(huì)使得柱塞運(yùn)動(dòng)速度發(fā)生變化,而柱塞的運(yùn)動(dòng)速度會(huì)使得對(duì)間隙內(nèi)流體的剪切運(yùn)動(dòng)產(chǎn)生影響,說(shuō)明由于抽油泵柱塞運(yùn)動(dòng)速度的變化會(huì)引起的間隙流體的剪切速度相應(yīng)的變化,當(dāng)柱塞的運(yùn)動(dòng)速度增加一定值后,間隙內(nèi)流體的剪切速度不再隨柱塞的運(yùn)動(dòng)速度的變化發(fā)生較大的變化。
2.4不同含水率對(duì)抽油泵泵筒和間隙間漏失率的影響
在油田開(kāi)采過(guò)程中,隨著采收率的減小,會(huì)使得油田的產(chǎn)量減小,所以對(duì)地層進(jìn)行注水會(huì)提高井底能量,是增加采收率的一種有效方式。同時(shí)這也增加了采收液的含水率,所以研究不同含水率對(duì)抽油泵泵筒和柱塞間隙的液體流動(dòng)情況及漏失率的大小,對(duì)提高抽油泵效率降低漏失率有重要的研究意義。
為分析不同含水率對(duì)漏失率的影響,設(shè)定抽油泵的工作制度一定的情況下(沖程1.8 m,沖次6次,進(jìn)出口壓差11.5 MPa)不同含水率(30%,40%,50%,60%、70%,80%,90%,100%)對(duì)應(yīng)的抽油泵泵筒與柱塞間隙的漏失率,并得出含水率對(duì)間隙漏失率的影響曲線(xiàn)如圖9所示。
由圖9可知,隨著含水率的增加,抽油泵柱塞和間隙間的漏失率會(huì)逐漸增大且增大趨勢(shì)逐漸減緩,當(dāng)含水率到達(dá)60%時(shí)漏失率又會(huì)減??;隨著含水量達(dá)到80%左右時(shí),抽油泵間隙的流體漏失率又會(huì)反彈上升,且上升的速率很快,直到含水率達(dá)到100%。這是由于含水率的變化主要影響流體的密度和粘度,而這兩個(gè)參數(shù)都是流體在間隙內(nèi)流動(dòng)的主要影響因素,隨著含水率的增加液體粘度會(huì)減小。密度的變化對(duì)間隙中間的流體影響較大,粘度的變化對(duì)靠近柱塞及泵筒壁面的流體速度的影響較大,因此這兩個(gè)因素共同影響間隙漏失率。具體粘度這一流體物性參數(shù)對(duì)抽油泵泵筒和柱塞間隙內(nèi)流體的影響,會(huì)在下面詳細(xì)給出。
2.5不同粘度對(duì)抽油泵泵筒和間隙間漏失率的影響
粘度是影響抽油泵泵筒與柱塞間隙的主要影響因素之一,當(dāng)柱塞在泵筒中做簡(jiǎn)諧運(yùn)動(dòng)時(shí),流體在間隙流道內(nèi)呈剪切流,而不同粘度的流體對(duì)剪切流的流速影響較大,進(jìn)而會(huì)對(duì)間隙的漏失率產(chǎn)生影響,所以分析不同粘度的流體對(duì)漏失率的影響極其重要。
設(shè)定在一定的抽油泵工作制度下(沖程1.8 m、沖次6次、進(jìn)出口壓差11.5 MPa),通入不同黏度的流體(0.27,0.22,0.18,0.14,0.09 Pa·s)分別計(jì)當(dāng)柱塞運(yùn)動(dòng)到不同位置(1.5,2.5,3,3.5,4.5 m)時(shí)的間隙內(nèi)液體在徑向方向上的流速分布,并可得出不同黏度對(duì)抽油泵柱塞和泵筒間間隙流體徑向流速分布曲線(xiàn)如圖10所示。
(a)距泵筒底1.5 m處間隙徑向流速分布
由圖10可以得出,柱塞從下死點(diǎn)往上死點(diǎn)的運(yùn)動(dòng)過(guò)程中,在1s這個(gè)時(shí)刻時(shí)柱塞中點(diǎn)大概位于距泵筒底1.5 m處,可以從圖10(a)、圖10(b)、圖10(e)觀(guān)測(cè)到,流體的粘度越大,流體流速整體趨勢(shì)越高,且越靠近間隙中間處的流速差異越明顯,而在越靠近柱塞及泵筒壁面處的流速差異越小,在壁面上幾乎無(wú)差異。這是由于粘性越高的流體相互之間及與柱塞和泵筒壁面之間的依附性越好,所以會(huì)使粘度越高的流體在徑向流速分布的差異越小。這充分解釋了流體粘度對(duì)間隙間的流速分布影響較大。隨著柱塞繼續(xù)向上運(yùn)動(dòng)的過(guò)程中,由于柱塞在做簡(jiǎn)諧運(yùn)動(dòng),所以在這個(gè)過(guò)程中速度逐漸增大,如圖10(b)所示這時(shí)不同粘度的流體之間的流速總體差異越明顯。從圖10(c)知當(dāng)2.5 s時(shí)柱塞運(yùn)動(dòng)到3 m時(shí)差異達(dá)到最大,同時(shí)也可以發(fā)現(xiàn)當(dāng)柱塞的運(yùn)動(dòng)速度不同時(shí),不同粘度的流體在靠近柱塞壁處的分布也不會(huì)趨于相同,且差異越來(lái)越大,同樣是因?yàn)橹乃俣仍娇?,流體在壁面的依附性越差,其粘度越小的流體在壁面的依附性越弱,而在靠近泵筒壁面這一側(cè)不同粘度的流體流速依然趨于一致。
通過(guò)計(jì)算5個(gè)不同粘度的流體的漏失率,可以繪制出如圖11所示的不同粘度對(duì)漏失率的影響曲線(xiàn),可以直觀(guān)的觀(guān)察到液體粘度對(duì)漏失率的影響。
圖11 不同粘度的流體對(duì)間隙漏失率的影響曲線(xiàn)
從圖11可以發(fā)現(xiàn),隨著流體粘度的增加,抽油泵泵筒與柱塞間間隙的漏失率逐漸增加,而且增長(zhǎng)趨勢(shì)逐漸放緩,且趨于穩(wěn)定。這說(shuō)明在一定程度上,粘度對(duì)漏失率的影響不再明顯。
(1)抽油泵柱塞與泵筒間間隙大小對(duì)間隙漏失率影響較大,即隨間隙越大漏失率增加,在一定范圍內(nèi),兩者近似呈線(xiàn)性變化。
(2)通過(guò)分析不同沖次和沖程對(duì)漏失率大小的影響時(shí)發(fā)現(xiàn),由于抽油泵柱塞和泵筒間隙內(nèi)流體主要是以剪切流的形式存在,沖程或沖次的變化均會(huì)引起間隙內(nèi)流體的剪切流形式發(fā)生變化,進(jìn)而影響間隙漏失率發(fā)生變化,在一定范圍內(nèi)抽油泵沖次或沖程的的增加都會(huì)使得漏失率增加,但隨著沖次的增加漏失率增加的趨勢(shì)會(huì)減緩,因此在保證產(chǎn)量的情況下,盡量選擇合適的抽油泵沖次與沖程可以有效的降低抽油泵間隙漏失。
(3)本文分析了含水率和粘度這兩個(gè)流體的物性參數(shù),含水率在不斷增加的過(guò)程中會(huì)呈現(xiàn)先增加,再減小,然后再增加的變化趨勢(shì),而隨著粘度的增加,間隙漏失率逐漸增加,但增加的趨勢(shì)逐漸減緩,因此本文的計(jì)算結(jié)果會(huì)對(duì)開(kāi)采不同物性的流體,在降低間隙漏失方面,提供理論參考依據(jù)。