李 越，白健華，蔣召平，于法浩，尚寶兵

(中海石油(中國(guó))有限公司 天津分公司，天津 300459)

機(jī)械采油技術(shù)是渤海油田井筒舉升工藝持續(xù)穩(wěn)產(chǎn)的基礎(chǔ)。截至2020年，渤海油田機(jī)械采油井?dāng)?shù)量占生產(chǎn)井總數(shù)的96%。

隨著渤海油田開發(fā)不斷進(jìn)行，油井均已進(jìn)入或即將進(jìn)入中高含水期。目前，渤海油田以高含水、高可采儲(chǔ)量采出程度并稱的“雙高”油田接近20個(gè)，儲(chǔ)量及年產(chǎn)量占比均達(dá)到50%以上。電潛泵技術(shù)具有占地面積小、排量大和舉升揚(yáng)程高等技術(shù)特點(diǎn)，成為了海上油田全生命周期穩(wěn)產(chǎn)和提液增產(chǎn)的主要手段[1]。在油田進(jìn)入中高含水期后，地層能量明顯下降，生產(chǎn)井提液穩(wěn)產(chǎn)的需求不斷增大。根據(jù)2015―2019年的統(tǒng)計(jì)結(jié)果，因生產(chǎn)井大泵提液需求導(dǎo)致的非故障檢泵比例已經(jīng)由2%提升至20%，增加了油井的生產(chǎn)成本。針對(duì)電潛泵結(jié)構(gòu)特點(diǎn)，提升葉導(dǎo)輪水力性能指標(biāo)，可直接拓寬電潛泵適用的排量范圍，適應(yīng)油井提液需求，降低換大泵頻次。葉導(dǎo)輪性能的好壞對(duì)離心泵性能影響很大，研究表明葉導(dǎo)輪的水力損失占離心泵水力總損失的40%～50%[2]。設(shè)計(jì)水力性能優(yōu)良的寬幅電潛泵葉導(dǎo)輪，就必須運(yùn)用流場(chǎng)分析軟件，并對(duì)設(shè)計(jì)結(jié)果進(jìn)行多次優(yōu)化[3-6]，比較通用且有效的方法為數(shù)值模擬結(jié)合樣件試驗(yàn)論證。本文應(yīng)用CFTurbo及有限元分析軟件，對(duì)寬幅電潛泵葉導(dǎo)輪開展了結(jié)構(gòu)設(shè)計(jì)和流場(chǎng)模擬分析，并采用3D打印技術(shù)和鑄造方法對(duì)樣件進(jìn)行試制，開展水力性能測(cè)試，以期為后續(xù)現(xiàn)場(chǎng)應(yīng)用提供參考。

1 實(shí)體建模

寬幅電潛泵舉升單元由1個(gè)葉導(dǎo)輪和1個(gè)葉導(dǎo)殼組成，流體舉升時(shí)，葉導(dǎo)殼靜止，軸承帶動(dòng)葉導(dǎo)輪單一方向旋轉(zhuǎn)，通過高速旋轉(zhuǎn)的離心力將動(dòng)能轉(zhuǎn)化為壓能，逐級(jí)傳遞將流體舉升至地面[7]。本文應(yīng)用CFTurbo軟件，選用三維模型建立寬幅電潛泵葉導(dǎo)輪及葉導(dǎo)殼模型，采用扭曲三角形法進(jìn)行繪圖[8]。在建模計(jì)算前，做如下假設(shè)[9-10]：①不考慮氣液兩相流及油水混相流，僅考慮以水為單一介質(zhì)的單相流；②流體在葉導(dǎo)輪內(nèi)做不可壓縮的穩(wěn)態(tài)流動(dòng)。

1.1 葉導(dǎo)輪參數(shù)設(shè)計(jì)

考慮渤海油田生產(chǎn)井產(chǎn)液量分布情況，產(chǎn)液量在100～400 m3/d的占比達(dá)50%以上，結(jié)合常規(guī)電潛泵的性能參數(shù)，設(shè)定寬幅電潛泵舉升單元額定排量288 m3/d，額定轉(zhuǎn)速2 850 r/min，額定點(diǎn)揚(yáng)程≥4 m，額定點(diǎn)泵效≥60%。導(dǎo)葉片數(shù)6片，導(dǎo)殼外徑103 mm，軸徑22.2 mm。圖1為寬幅電潛泵葉導(dǎo)輪和葉導(dǎo)殼的三維模型。

圖1 寬幅電潛泵葉導(dǎo)輪和葉導(dǎo)殼三維模型

在設(shè)計(jì)過程中，葉導(dǎo)輪和導(dǎo)殼都屬于三維復(fù)雜曲面結(jié)構(gòu)，設(shè)計(jì)材質(zhì)選用鑄鐵合金，砂模鑄造成型，并配合拋光和噴砂整形等工藝，保證葉導(dǎo)輪和導(dǎo)殼的精確性。

1.2 數(shù)學(xué)建模與網(wǎng)格劃分

應(yīng)用CFTurbo軟件針對(duì)葉導(dǎo)輪進(jìn)行數(shù)學(xué)模型建模。因葉導(dǎo)輪及葉導(dǎo)殼的葉片部分為三維扭曲結(jié)構(gòu)，需要對(duì)實(shí)體模型通過“二維→三維”的順序搭建實(shí)體模型，隨后設(shè)定葉導(dǎo)輪各項(xiàng)基本參數(shù)。研究表明，當(dāng)級(jí)數(shù)≥2時(shí)各級(jí)葉輪的進(jìn)口液體均是有旋運(yùn)動(dòng)[11-12]，隨著級(jí)數(shù)的增加(三級(jí)之后)，模型的單級(jí)揚(yáng)程、效率和軸功率的參數(shù)值與第2級(jí)葉輪基本一致[13]，考慮模擬的準(zhǔn)確性和時(shí)效性，本文采用三級(jí)模型進(jìn)行數(shù)值模擬。

利用Ansys對(duì)葉導(dǎo)輪和葉導(dǎo)殼的水體區(qū)域進(jìn)行網(wǎng)格劃分，相同模型的網(wǎng)格數(shù)量越多，模擬計(jì)算結(jié)果越準(zhǔn)確，同等數(shù)量的結(jié)構(gòu)化網(wǎng)格比非結(jié)構(gòu)化網(wǎng)格的計(jì)算精度更高，本文選用結(jié)構(gòu)化網(wǎng)格對(duì)模型進(jìn)行數(shù)值模擬，并對(duì)所有葉片的前端和尾部進(jìn)行網(wǎng)格加密處理，如圖2所示。

圖2 葉導(dǎo)輪和葉導(dǎo)殼水體三維網(wǎng)格劃分

1.3 邊界條件

選用標(biāo)準(zhǔn)k-ε湍流模型，假設(shè)流動(dòng)為完全湍流，忽略分子黏性的影響，這與預(yù)期設(shè)置流動(dòng)介質(zhì)為單相純水液流吻合，并且該模型對(duì)于三維旋轉(zhuǎn)流動(dòng)具有很好的效果[14]。在進(jìn)口端面給定進(jìn)口排量，進(jìn)口湍動(dòng)能k和湍動(dòng)能耗散率ε為：

(1)

(2)

式中：uin為進(jìn)口平均速度；lin為進(jìn)口混合長(zhǎng)度；Cμ為經(jīng)驗(yàn)常數(shù)，取值0.09。

收斂精度選用1×10-4，保證三維模擬的計(jì)算準(zhǔn)確率。邊界條件定義了進(jìn)口、出口和壁面3個(gè)類型。本文預(yù)置液體為單相純水液流，葉片的進(jìn)口處即為單一均相流動(dòng)[15-16]。導(dǎo)葉、葉片和輪轂等為壁面邊界，不考慮葉導(dǎo)殼本體和泵之間的空隙，因此模型為一個(gè)整體壁面模型。葉輪、葉片和流道內(nèi)均認(rèn)定為絕熱無滑移邊界[17]，壁面粗糙度為12.5 μm。

2 仿真結(jié)果與分析

對(duì)三維建模的寬幅電潛泵葉導(dǎo)輪和葉導(dǎo)殼內(nèi)部湍流，應(yīng)用有限元軟件進(jìn)行數(shù)模計(jì)算與結(jié)果分析，得到內(nèi)部的壓力場(chǎng)和速度場(chǎng)，并計(jì)算不同排量下的揚(yáng)程和效率，與樣件試驗(yàn)數(shù)據(jù)進(jìn)行比較。

本文模型針對(duì)3級(jí)葉導(dǎo)輪進(jìn)行模擬，為保證數(shù)據(jù)真實(shí)性和準(zhǔn)確性，參考電潛泵高效區(qū)范圍參數(shù)，對(duì)共計(jì)10個(gè)模型的不同工況點(diǎn)均進(jìn)行了模擬計(jì)算，以排量為變量(0.35Q額定～1.6Q額定)，其余參數(shù)保持不變，以最終設(shè)計(jì)泵(葉片包角60.8～60.5°均勻變化，葉片厚度1.2～1.5 mm均勻變化)進(jìn)行分析。

2.1 壓力場(chǎng)

葉導(dǎo)輪工作狀態(tài)下的壓力分布是判定葉導(dǎo)輪性能的重要指標(biāo)[18]。對(duì)葉導(dǎo)輪舉升單元流場(chǎng)進(jìn)行三維湍流數(shù)值模擬分析，得到葉導(dǎo)輪表面流體壓力場(chǎng)分布，如圖3所示。

從圖3中可以看出，在不同排量工況下，葉片對(duì)流體做功，將離心力轉(zhuǎn)化為壓力，葉導(dǎo)輪入口端至出口端壓力不斷增大，并在葉導(dǎo)輪出口端達(dá)到單級(jí)舉升單元的最大值，符合水力機(jī)械流動(dòng)規(guī)律。隨著排量的升高，出口壓力值逐漸降低，呈現(xiàn)出揚(yáng)程逐步下降的趨勢(shì)。葉片在小排量工況(0.35Q額定)時(shí)，葉片出口位置存在局部“壓力指進(jìn)”現(xiàn)象，表現(xiàn)為壓力梯度較大；在大排量工況(1.6Q額定)時(shí)，葉片背部存在明顯的低壓區(qū)域，這是由于在偏離額定工況較多時(shí)，液體流向和出口安放角不吻合，形成了沖擊現(xiàn)象，在小排量工況和大排量工況下均導(dǎo)致整體泵效降低。

圖3 寬幅葉導(dǎo)輪壓力場(chǎng)分析

2.2 速度場(chǎng)分析

圖4為葉導(dǎo)輪舉升單元內(nèi)流體速度矢量場(chǎng)分布圖，分別模擬了小排量工況、額定排量工況和大排量工況條件下的速度場(chǎng)。從圖4中可以看出，在不同排量工況下，模型葉片的入口處速度最低，從入口向出口方向逐漸增大，流道的外邊緣出口速度高于內(nèi)邊緣，這是離心力作用的結(jié)果。在小排量工況條件下，葉輪模型葉片的背面均出現(xiàn)了明顯的流動(dòng)分離現(xiàn)象，漩渦區(qū)域面積較大，這會(huì)導(dǎo)致泵效降低，水力效率下降；在大排量工況條件下，葉輪模型葉片背面的流動(dòng)分離現(xiàn)象減弱，漩渦區(qū)域面積減?。辉陬~定排量工況附近，流道內(nèi)均未出現(xiàn)漩渦，但仍有輕微的流動(dòng)分離現(xiàn)象，水力效率達(dá)到最大值。由于葉片的三維扭曲和流道內(nèi)外流動(dòng)速度不一致，葉片背面的流動(dòng)分離現(xiàn)象是不可消除的，但可通過結(jié)構(gòu)調(diào)整降低這種現(xiàn)象的影響，流動(dòng)分離現(xiàn)象是導(dǎo)致水力效率降低的主要因素。

圖4 寬幅葉導(dǎo)輪速度場(chǎng)分析

2.3 數(shù)值模擬與試驗(yàn)對(duì)比

為驗(yàn)證數(shù)值模擬計(jì)算結(jié)果是否準(zhǔn)確，應(yīng)用常規(guī)鑄鐵鑄造樣件對(duì)數(shù)模結(jié)果進(jìn)行1:1開模試制，同時(shí)應(yīng)用3D打印技術(shù)進(jìn)行樣件試制，測(cè)試該技術(shù)在葉導(dǎo)輪水力性能測(cè)試中的表現(xiàn)。圖5為寬幅葉導(dǎo)輪和導(dǎo)殼的3D打印樣件和鑄造樣件模型。

圖5 寬幅葉導(dǎo)輪和導(dǎo)殼3D打印樣件及鑄造樣件模型

試驗(yàn)選用水泵性能測(cè)試臺(tái)，測(cè)試泵的輸入轉(zhuǎn)速、轉(zhuǎn)矩、泵進(jìn)口壓力、泵出口壓力、泵出口排量、電機(jī)功率和揚(yáng)程等參數(shù)。測(cè)試臺(tái)主要由控制柜、電機(jī)、水箱和試驗(yàn)泵構(gòu)成，并配有轉(zhuǎn)矩傳感器、壓力傳感器、流量計(jì)、溫度計(jì)和出口流量調(diào)節(jié)閥等。

通過泵的性能測(cè)試，得到了3D打印樣件和鑄造樣件在不同排量下的揚(yáng)程、功率和水力效率，如表1～2所示。

表1 3D打印樣件水力性能測(cè)試參數(shù)

在3D打印樣件和鑄造樣件的試驗(yàn)過程中，葉輪采用單級(jí)浮動(dòng)安裝方式，并選用端面密封，試驗(yàn)密封性能良好。從試驗(yàn)結(jié)果可看出，在額定排量Q=288 m3/h條件下，3D打印樣件測(cè)試揚(yáng)程4.64 m，水力效率62.87 %；鑄造樣件測(cè)試揚(yáng)程4.87 m，水力效率65.83 %；數(shù)值模擬在相同排量工況下預(yù)測(cè)揚(yáng)程4.92 m，水力效率66.7 %。三者隨排量的變化趨勢(shì)一致，鑄造樣件和數(shù)值模擬結(jié)果誤差更小，證明本文數(shù)值模擬結(jié)果具有較高的精度和一致性。

表2 鑄造樣件水力性能測(cè)試參數(shù)

因試驗(yàn)臺(tái)試驗(yàn)?zāi)Ｐ捅眉?jí)數(shù)與數(shù)模泵級(jí)數(shù)不同，需將數(shù)值模擬與試驗(yàn)臺(tái)測(cè)試獲得的預(yù)測(cè)揚(yáng)程和功率等參數(shù)結(jié)果統(tǒng)一，再根據(jù)揚(yáng)程和功率計(jì)算出效率，并進(jìn)行參數(shù)對(duì)比分析。數(shù)值模擬計(jì)算結(jié)果與試驗(yàn)臺(tái)測(cè)試結(jié)果對(duì)比如圖6。

圖6 數(shù)值模擬與樣件水力性能測(cè)試結(jié)果對(duì)比

3種試驗(yàn)與數(shù)值模擬結(jié)果雖然趨勢(shì)一致，但仍存在誤差。鑄造樣件與數(shù)學(xué)模型結(jié)果誤差的主要原因?yàn)椋孩贁?shù)學(xué)模型對(duì)主要部件的建模精度較高，但對(duì)葉輪進(jìn)口和出口位置的參數(shù)進(jìn)行了簡(jiǎn)化，且未考慮流道壁面粗糙度對(duì)流體流動(dòng)的影響，造成了誤差；②數(shù)學(xué)模型忽略了泵的漏失影響，導(dǎo)致數(shù)學(xué)模型結(jié)果優(yōu)于實(shí)際測(cè)試結(jié)果。因此，在后續(xù)計(jì)算過程中，應(yīng)根據(jù)實(shí)際情況對(duì)建模參數(shù)再度優(yōu)化，提高模擬精度。

3D打印樣件與鑄造樣件存在誤差的主要原因?yàn)椋孩?D打印樣件成品的表面粗糙度和精度相比于鑄造樣件低，導(dǎo)致漏失增多，機(jī)械損失增大；②3D打印樣件材料為ABS，鑄造樣件材料為鑄鐵，減磨墊材質(zhì)為CEF酚醛布板，當(dāng)3D打印樣件裝配減磨墊進(jìn)行水力性能測(cè)試時(shí)，隨著測(cè)試時(shí)間延長(zhǎng)，3D打印樣件會(huì)出現(xiàn)一定程度的磨損，導(dǎo)致水力性能參數(shù)值偏差增大。

在泵額定排量工況條件下，數(shù)值模擬得到的效率與鑄造樣件的試驗(yàn)差約1.3% ，在誤差允許的范圍內(nèi)，這證明了用數(shù)值模擬方法來預(yù)測(cè)寬幅電潛泵水力特性的可行性，并驗(yàn)證了此設(shè)計(jì)方法準(zhǔn)確。同時(shí)，對(duì)比3D打印樣件和鑄造樣件的試驗(yàn)結(jié)果，誤差值約4.5%，數(shù)值模擬結(jié)果和3D打印樣件誤差值為6.1%。雖然結(jié)果存在一定誤差，但三者試驗(yàn)結(jié)果隨排量變化的趨勢(shì)一致，且3D打印樣件費(fèi)用低，周期短，適合進(jìn)行多輪次葉導(dǎo)輪優(yōu)化設(shè)計(jì)與性能測(cè)試，為最終設(shè)計(jì)方案節(jié)約時(shí)間，提供依據(jù)。

2.4 寬幅電潛泵樣件與常規(guī)電潛泵樣件對(duì)比

針對(duì)目前海上油田實(shí)際應(yīng)用的常規(guī)電潛泵高效區(qū)排量范圍窄、大排量點(diǎn)揚(yáng)程較低的現(xiàn)狀，結(jié)合提液增產(chǎn)實(shí)際需求并綜合上述研究，從影響葉導(dǎo)輪性能的結(jié)構(gòu)尺寸和三維性能參數(shù)等角度進(jìn)行葉導(dǎo)輪參數(shù)的優(yōu)化設(shè)計(jì)，研制出高效區(qū)流量100 ～450 m3/d的寬幅電潛泵葉導(dǎo)輪，該葉導(dǎo)輪與相近額定排量的常規(guī)電潛泵葉導(dǎo)輪性能對(duì)比參數(shù)如圖7。

由圖7可以看出：①相同排量下，寬幅電潛泵葉導(dǎo)輪水力效率均高于常規(guī)電潛泵葉導(dǎo)輪，額定點(diǎn)排量為288 m3/d，寬幅電潛泵水力效率比常規(guī)電潛泵水力效率高10.21%；②揚(yáng)程曲線中，在額定點(diǎn)排量之前，常規(guī)電潛泵單級(jí)葉導(dǎo)輪揚(yáng)程高于寬幅電潛泵，隨著排量增大，寬幅電潛泵單級(jí)葉導(dǎo)輪揚(yáng)程下降趨勢(shì)更緩，在排量超過380 m3/d后，寬幅電潛泵單級(jí)葉導(dǎo)輪揚(yáng)程高于常規(guī)電潛泵；③常規(guī)電潛泵單級(jí)葉導(dǎo)輪高效區(qū)流量150～418 m3/d，寬幅電潛泵單級(jí)葉導(dǎo)輪高效區(qū)流量100～458 m3/d，高效區(qū)拓寬15.7%，更有利于生產(chǎn)井提液穩(wěn)產(chǎn)時(shí)井液的舉升。

圖7 寬幅電潛泵和常規(guī)電潛泵單級(jí)葉導(dǎo)輪水力性能對(duì)比曲線

3 結(jié)論

1) 基于CFTurbo和有限元分析軟件給出寬幅葉導(dǎo)輪建模和邊界條件的界定方法，通過數(shù)值模擬得到了內(nèi)部壓力場(chǎng)和速度場(chǎng)分布規(guī)律，減少壓力局部“壓力指進(jìn)”現(xiàn)象和速度的渦流與流動(dòng)分離現(xiàn)象可提升葉導(dǎo)輪性能。

2) 數(shù)值模擬預(yù)測(cè)結(jié)果在小排量和額定排量工況與3D打印樣件及鑄造樣件試驗(yàn)結(jié)果吻合度較高，額定排量工況下，數(shù)值模擬預(yù)測(cè)結(jié)果與3D打印樣件及鑄造樣件差值分別達(dá)到6.1%和1.3%。在大排量工況下偏差增大，這是由于樣件在大排量工況下存在漏失，從而對(duì)泵的水力性能參數(shù)影響增大。在后續(xù)數(shù)值模擬過程中應(yīng)考慮泵的間隙漏失，提高模擬結(jié)果精度。

3) 3D打印樣件試驗(yàn)結(jié)果與鑄造樣件在額定排量工況點(diǎn)的誤差值約4.5%，誤差較小，且水力性能曲線隨排量變化的趨勢(shì)一致。 3D打印樣件費(fèi)用低，周期短，適合進(jìn)行多輪次葉導(dǎo)輪優(yōu)化設(shè)計(jì)與性能測(cè)試，為最終設(shè)計(jì)方案提供試驗(yàn)參數(shù)依據(jù)并縮短設(shè)計(jì)周期。

4) 對(duì)比寬幅電潛泵和常規(guī)電潛泵單級(jí)葉導(dǎo)輪性能測(cè)試結(jié)果，可以看出在額定點(diǎn)排量工況下，寬幅電潛泵比常規(guī)電潛泵水力效率高10.21%；寬幅電潛泵單級(jí)葉導(dǎo)輪高效區(qū)流量100～458 m3/d，相比于常規(guī)電潛泵單級(jí)葉導(dǎo)輪，高效區(qū)范圍拓寬15.7%，達(dá)到預(yù)期設(shè)計(jì)要求。寬幅電潛泵葉導(dǎo)輪舉升單元裝配形成寬幅電潛泵后，更有利于生產(chǎn)井提液穩(wěn)產(chǎn)時(shí)井液的舉升。