抽油泵固定凡爾球運(yùn)動(dòng)規(guī)律的數(shù)值模擬

王小兵，呂雷綱，李 森，邊 斌，劉 陽，龔浩宇，王多琦

（1.常州大學(xué) 石油工程學(xué)院，江蘇常州 213164；2.Biological and Environmental Sciences and Engineering（BESE）Division，Water Desalination and Reuse Center（WDRC），King Abdullah University of Science and Technology，Thuwal 23955-6900，Kingdom of Saudi Arabia）

關(guān)鍵字：有桿抽油泵；固定凡爾球運(yùn)動(dòng)規(guī)律；多物理耦合；動(dòng)網(wǎng)格；接觸壓力方法

0 引言

國內(nèi)各油田大約80%的生產(chǎn)井使用有桿抽油技術(shù)，全國各大油田60%的產(chǎn)液量、75%的產(chǎn)油量是靠有桿抽油設(shè)備采出的［1］。抽油泵作為有桿抽油系統(tǒng)中核心的流體機(jī)械，工作原理是其固定凡爾球與游動(dòng)凡爾球的周期性啟閉控制泵內(nèi)流體的吸入與排出。

通過試驗(yàn)對抽油泵固定凡爾球井下運(yùn)動(dòng)規(guī)律的測定有很大的困難，因此需要借助于數(shù)值模擬等方法。試驗(yàn)研究［2］表明了固定凡爾球的實(shí)際運(yùn)動(dòng)狀態(tài)，但未能得出太多具有參考價(jià)值的結(jié)論；自1967年，ADOLPH［3］考慮固定凡爾球的受力以及運(yùn)動(dòng)特性，利用差分法構(gòu)建了固定凡爾球運(yùn)動(dòng)規(guī)律的微分方程后，孟英峰等［4-8］從數(shù)學(xué)模型方面解決微分方程中的奇點(diǎn)問題、魏氏效應(yīng)等。余海洋等［9］基于MATLAB，針對有桿抽油泵的上沖程階段進(jìn)行模擬研究，并探討了流體的動(dòng)力特性等因素對吸入特性的影響；吳建軍等［10］借鑒前者的理論成果，通過Simulink對抽油泵固定凡爾進(jìn)行仿真；李凌川等［11］基于FLUENT軟件，結(jié)合動(dòng)網(wǎng)格技術(shù)和UDF編程，對流體進(jìn)泵的過程進(jìn)行數(shù)值模擬得到了固定凡爾球的運(yùn)動(dòng)規(guī)律，在解決固定凡爾球與球罩的碰撞問題時(shí)采用了彈性恢復(fù)系數(shù)來描述固定凡爾球碰撞后的運(yùn)動(dòng)規(guī)律。以上文獻(xiàn)存在一些未能解決的問題：建立的數(shù)值模擬對柱塞位移引起泵內(nèi)體積增大的控制是依靠設(shè)置固定凡爾排出口為動(dòng)邊界，這對網(wǎng)格的影響極大，降低了計(jì)算的精度；對固定凡爾球與泵內(nèi)壁面接觸影響網(wǎng)格質(zhì)量的解決方法并未具體說明，且研究中運(yùn)用彈性回復(fù)系數(shù)描述凡爾球運(yùn)動(dòng)的方法忽略了流體阻力的影響。參考數(shù)值模擬在相關(guān)領(lǐng)域的研究［12-18］，為現(xiàn)有的問題提供解決方案。針對以上問題，應(yīng)用COMSOL Multiphysics 仿真軟件，對黏性流體中抽油泵固定凡爾球的運(yùn)動(dòng)規(guī)律進(jìn)行數(shù)值模擬；將CFD方法與固體力學(xué)耦合，并結(jié)合動(dòng)網(wǎng)格技術(shù)以及接觸壓力方法，建立固定凡爾球的流固耦合模型。解決了上沖程階段柱塞、固定凡爾球運(yùn)動(dòng)以及與其他固體壁面接觸引起的網(wǎng)格變形和奇點(diǎn)問題，進(jìn)而建立更加合理的固定凡爾球運(yùn)動(dòng)模型。研究結(jié)果得到了泵內(nèi)流場分布以及固定凡爾球的運(yùn)動(dòng)規(guī)律曲線。所得結(jié)果為實(shí)際生產(chǎn)中抽油泵工作制度優(yōu)化等提供了理論依據(jù)。

1 基本理論分析

1.1 假設(shè)條件

固定凡爾球做豎直向上的直線運(yùn)動(dòng)，而且伴隨兩種旋轉(zhuǎn)運(yùn)動(dòng)，若柱塞運(yùn)動(dòng)速度較小時(shí)，固定凡爾球繞水平軸自上而下轉(zhuǎn)動(dòng)［2，19］；當(dāng)柱塞的運(yùn)動(dòng)速度較大時(shí)，固定凡爾球繞垂向軸水平方向自轉(zhuǎn)，同時(shí)，沿閥座內(nèi)孔邊角繞球座孔軸線公轉(zhuǎn)；因?yàn)楸猛搀w積狹小，橫向的漂移很小。為簡化計(jì)算，提出以下幾點(diǎn)假設(shè)：（1）油井為直井，泵筒無傾斜；（2）閥球沿閥座中軸線運(yùn)動(dòng)，不考慮自轉(zhuǎn)及公轉(zhuǎn)；（3）柱塞與泵筒間不存在漏失；（4）泵內(nèi)流體為油水混合物，且不含溶解氣。

1.2 固定凡爾球的受力分析

在假設(shè)凡爾球不存在自轉(zhuǎn)和公轉(zhuǎn)的前提下，凡爾球的受力情況如圖1所示。凡爾球在垂向上受到自身產(chǎn)生的重力G、浮力FB以及流體對凡爾球的繞球作用力FD三者的共同作用。此外，凡爾球與壁面碰撞時(shí)還應(yīng)考慮接觸壓力。

圖1 凡爾球受力示意Fig.1 Schematic stress diagram of valve ball

2 基于COMSOL Multiphysics的多物理場耦合模型

單純的CFD方法不能解決流體域與固體之間的接觸壓力問題，還需要考慮固體力學(xué)解決流體域中兩固體的彈性碰撞問題。因此使用COMSOL Multiphysics建立流固耦合模型進(jìn)行瞬態(tài)研究，耦合類型選擇全耦合。完成模型選擇之后，在模型開發(fā)器界面中主要包含組件和研究兩個(gè)部分。組件中完成幾何模型、邊界條件及網(wǎng)格的設(shè)置，研究中選擇求解器和時(shí)間步長。

2.1 幾何模型

抽油泵固定凡爾是軸對稱結(jié)構(gòu)，所以選擇二維軸對稱的空間結(jié)構(gòu)進(jìn)行抽油泵固定凡爾幾何模型的建立能節(jié)約計(jì)算內(nèi)存且節(jié)省計(jì)算時(shí)間。研究外徑70 mm的固定凡爾，幾何模型如圖2所示。

圖2 抽油泵固定凡爾的幾何模型Fig.2 The geometric model of the standing valve of the oil well pump

2.2 抽油泵固定凡爾的物理模型

抽油泵固定凡爾的物理模型包括流體流動(dòng)模型和固體力學(xué)模型，流體流動(dòng)模型通過對雷諾數(shù)的計(jì)算選擇層流或是湍流；固體力學(xué)模型通過材料屬性選擇研究對象為線彈性材料、非線彈性材料或其他。抽油泵固定凡爾球的打開或關(guān)閉都會(huì)使泵內(nèi)形成狹窄流道，則流動(dòng)形態(tài)為湍流，湍流模型中的k-ε模型對于復(fù)雜幾何周圍外部流動(dòng)問題的求解效果良好，因此流體流動(dòng)模型選擇k-ε模型，考慮流體可壓縮性且包含重力。抽油泵及其配件的材料為鋼質(zhì)，則固體力學(xué)模型的研究對象確定為線彈性材料，對凡爾球設(shè)置重力，指定柱塞位移，其他部件設(shè)置為固定。

2.3 接觸壓力方法

柱塞未開始運(yùn)動(dòng)時(shí)，固定凡爾球由于重力原因落在球座上從而關(guān)閉吸入口，速度為零，此時(shí)固定凡爾球上下兩部分的流體不連通、流體連續(xù)性被破環(huán)，從而出現(xiàn)了奇點(diǎn)問題。根據(jù)機(jī)械密封理論［20］，固定凡爾球與球座的接觸面有一層液體薄膜，說明保證流體連續(xù)性即通過數(shù)學(xué)方法給固定凡爾球與泵內(nèi)會(huì)產(chǎn)生接觸的壁面設(shè)置偏移距離實(shí)現(xiàn)對奇點(diǎn)問題的解決。

首先，在模型開發(fā)器界面中組件的定義節(jié)點(diǎn)中添加接觸對，選擇凡爾球的壁面為源邊界、選擇球座及球罩的接觸壁面為目標(biāo)邊界如圖3所示。

圖3 定義壁面的接觸對Fig.3 Definition of the contact pair of the wall

其次，在組件中的固體力學(xué)模塊中添加接觸節(jié)點(diǎn)，選用定義的接觸對；接觸壓力方法選擇罰函數(shù)；偏移距離的設(shè)置需要重復(fù)測試，最后優(yōu)化為5 μm。該方法在不影響壁面接觸壓力的同時(shí)實(shí)現(xiàn)了流體的連續(xù)性，從而解決了凡爾球與泵內(nèi)壁面的接觸引起的網(wǎng)格問題。

2.4 動(dòng)網(wǎng)格及網(wǎng)格劃分

固定凡爾球運(yùn)動(dòng)時(shí)必須使用動(dòng)網(wǎng)格技術(shù)適應(yīng)邊界的運(yùn)動(dòng)或變形。在動(dòng)網(wǎng)格節(jié)點(diǎn)中，添加變形域節(jié)點(diǎn)并選擇所有流體域，網(wǎng)格平滑類型選擇Yeoh（更易于收斂）；添加固定邊界節(jié)點(diǎn)并選擇出入口邊界以及固體壁面；添加指定網(wǎng)格位移節(jié)點(diǎn)并選擇會(huì)發(fā)生移動(dòng)的所有邊界。

對柱塞段的網(wǎng)格劃分采用映射網(wǎng)格，選擇分布方法對各邊界上的網(wǎng)格數(shù)量進(jìn)行定量劃分，從而解決了大變形以及網(wǎng)格大量的重復(fù)劃分造成的內(nèi)存不足問題。對固定凡爾的網(wǎng)格劃分采用自由三角形網(wǎng)格，并添加邊界層和角細(xì)化對邊界進(jìn)行處理，對流體域進(jìn)行加密、固體域適當(dāng)粗化，網(wǎng)格如圖4所示。

圖4 網(wǎng)格劃分Fig.4 Mesh division

2.5 自動(dòng)重新劃分網(wǎng)格

固體域位置的不斷變化會(huì)導(dǎo)致網(wǎng)格質(zhì)量變差，所以對大位移的模型要進(jìn)行自動(dòng)重新劃分網(wǎng)格的設(shè)置。自動(dòng)重新劃分網(wǎng)格的操作步驟：在研究的步驟節(jié)點(diǎn)中，選用自動(dòng)重新劃分網(wǎng)格；然后獲取求解器配置的初始值，在求解器配置節(jié)點(diǎn)中的瞬態(tài)求解器中設(shè)置自動(dòng)重新劃分網(wǎng)格的條件，選擇條件類型為失真并設(shè)置超出失真范圍時(shí)停止為1。

2.6 固定凡爾球的初始位置

當(dāng)柱塞未開始工作時(shí)，固定凡爾球落于球座上并關(guān)閉吸入口。在幾何模型的繪制時(shí)將固定凡爾球與球座進(jìn)行配合，會(huì)破壞流體域連續(xù)性的網(wǎng)格拓?fù)浣Y(jié)構(gòu)，則凡爾球最初需要一段初始距離使固定凡爾球在重力作用下落在球座上，同時(shí)這一過程要對接觸壓力方法以及偏移進(jìn)行調(diào)整優(yōu)化。對柱塞位移施加時(shí)間條件，固定凡爾球因重力下落在球座上，時(shí)間控制函數(shù)為：（t［1/s］＞ 0），該函數(shù)表示與時(shí)間相關(guān)的變量從0 s開始施加，從而為模型計(jì)算提供了合理的初始值。

時(shí)間步設(shè)置為（-t0，Δt，t1），表示開始計(jì)算之前的t0時(shí)間段內(nèi)固定凡爾球落于球座上，完成凡爾球初始位置的設(shè)定。在開始計(jì)算時(shí)，突然施加入口壓力改變會(huì)導(dǎo)致計(jì)算錯(cuò)誤，因此在入口壓力邊界寫入階躍函數(shù)，從而解決壓力突變造成計(jì)算的不收斂。

3 結(jié)果分析

模型的計(jì)算參數(shù)為：固定凡爾外徑70 mm、腔體內(nèi)徑48 mm，凡爾球直徑40 mm，球座厚度10 mm、內(nèi)徑30 mm，吸入口直徑34 mm；泵筒上端流域長5 m、內(nèi)徑48 mm；固定凡爾、凡爾球、球座的材料均為鋼質(zhì)，密度為7 850 kg/m3、楊氏模量為210 GPa、泊松比為 0.3；泵內(nèi)流體黏度為 10 mPa·s、密度為900 kg/m3；上沖程階段游動(dòng)凡爾球關(guān)閉，以固體域近似為柱塞，長度為0.6 m。入口壓力為3 MPa，出口壓力為10 MPa，柱塞沖程為3 m、沖次為6 r/min，下泵深度為1 000 m。研究該工況下固定凡爾球的運(yùn)動(dòng)規(guī)律，分析柱塞沖次和泵內(nèi)流體粘度對固定凡爾球運(yùn)動(dòng)的影響。

3.1 固定凡爾運(yùn)動(dòng)規(guī)律分析

抽油泵固定凡爾球的位移曲線和打開程度反映了凡爾球打開的響應(yīng)速度，所得結(jié)果如圖5所示。

圖5 固定凡爾球運(yùn)動(dòng)規(guī)律曲線Fig.5 Motion curve of the standing valve ball

由圖5可知：泵內(nèi)流體為油水混合物且不含氣體時(shí)，泵內(nèi)較小的空間內(nèi)液體體積膨脹很微弱，凡爾球的打開過程幾乎不存在滯后現(xiàn)象，在柱塞上行的同時(shí)固定凡爾球打開；隨著柱塞的不斷上行，凡爾球上升高度越來越高，且運(yùn)動(dòng)過程中有一定的波動(dòng)并且在兩個(gè)時(shí)間點(diǎn)存在明顯的小幅度回落，最大升程為3.88 mm，所需開啟時(shí)間為2.43 s；在整個(gè)上沖程階段，固定凡爾球并未到達(dá)球罩所在高度而是在一定高度處懸浮，且柱塞上行的后半段（柱塞的減速上升階段），固定凡爾球逐漸回落且最終幾乎關(guān)閉吸入口，滯后關(guān)閉時(shí)間較短。

3.2 沖次對固定凡爾球位移的影響

研究的沖次范圍為 2，4，6，8，10 r/min，得到固定凡爾球的運(yùn)動(dòng)規(guī)律曲線如圖6所示。

圖6 沖次對固定凡爾球升程的影響曲線Fig.6 The influence curve of times of stroke on the motion law of standing valve ball

由圖6可知，沖次越高固定凡爾球運(yùn)動(dòng)過程中的波動(dòng)性越強(qiáng)且運(yùn)動(dòng)過程中出現(xiàn)的瞬間回落或者上升的幅度越大；在研究的沖次范圍內(nèi)，固定凡爾球都未到達(dá)球罩所在高度而是在一定高度處懸浮，沖次越高升程越高；在柱塞上行的后半段（柱塞的減速上升階段），固定凡爾球逐漸回落且最終幾乎關(guān)閉吸入口，沖次越高關(guān)閉程度越小。沖次越高升程越大，所需的開啟時(shí)間越短；當(dāng)沖次為2 r/min，最大升程為2.46 mm，所需開啟時(shí)間為7.38 s；當(dāng)沖次為4 r/min，最大升程為3.71 mm，所需開啟時(shí)間為3.8 s；沖次為6 r/min，最大升程為3.88 mm，所需開啟時(shí)間為2.43 s；當(dāng)沖次為8 r/min，最大升程為4.33 mm，所需開啟時(shí)間為1.7 s；當(dāng)沖次為10 r/min，最大升程為5.26 mm，所需開啟時(shí)間為1.27 s。

3.3 泵內(nèi)流體黏度對固定凡爾球位移的影響

研究的泵內(nèi)流體黏度為 10，30，50，70，90，110 mPa·s，并確定出凡爾球完全打開時(shí)的臨界黏度，得到固定凡爾球的運(yùn)動(dòng)規(guī)律曲線如圖7所示。

圖7 泵內(nèi)流體黏度對固定凡爾球位移的影響曲線Fig.7 The influence curve of viscosity of fluid in pump on the motion law of standing valve ball

如圖7所示，在研究的泵內(nèi)流體黏度范圍內(nèi)，黏度小于54 mPa·s（臨界黏度）時(shí)，凡爾球的打開程度隨泵內(nèi)流體黏度的增大而增加，但未完全打開；黏度大于54 mPa·s時(shí)，固定凡爾球完全打開，且黏度越大打開時(shí)間越短；當(dāng)凡爾球由于球罩限制，與其發(fā)生碰撞時(shí)，凡爾球經(jīng)過反彈后在流體作用力下懸浮在球罩處；在柱塞上行的后半段（柱塞的減速上升階段），固定凡爾球逐漸回落且最終幾乎關(guān)閉吸入口，黏度越高滯后關(guān)閉時(shí)間越長，且黏度大于54 mPa·s時(shí)，固定凡爾球在接近球座時(shí)波動(dòng)性更強(qiáng)。流體黏度越高固定凡爾球的升程越大，固定凡爾球達(dá)到最大打開程度所需開啟時(shí)間越短，黏度為10 mPa·s，最大升程為3.88 mm，所需開啟時(shí)間為2.43 s；當(dāng)黏度為30 mPa·s，最大升程為4.29 mm，所需開啟時(shí)間為2.16 s；黏度為50 mPa·s，最大升程為4.77 mm，所需開啟時(shí)間為2.17 s；黏度為54 mPa·s時(shí)，固定凡爾球完全打開，所需開啟時(shí)間為1.42s；流體黏度為70 mPa·s，固定凡爾球完全打開，所需開啟時(shí)間為1.28 s；流體黏度為90 mPa·s，固定凡爾球完全打開，所需開啟時(shí)間為1.17 s。

4 結(jié)論

（1）構(gòu)建抽油泵固定凡爾球的流固耦合模型運(yùn)用了動(dòng)網(wǎng)格技術(shù)及接觸壓力方法，對入口壓力邊界結(jié)合階躍函數(shù)解決了奇點(diǎn)問題和固定凡爾球與泵腔內(nèi)壁面的接觸問題，從而為柱塞泵球形凡爾的流體機(jī)械設(shè)計(jì)、優(yōu)化提供了數(shù)值模擬方法。

（2）沖次增加固定凡爾球的升程增加，但在研究范圍內(nèi)固定凡爾球完全打開；沖次增加對固定凡爾球上升過程的穩(wěn)定性有不利影響，沖次越高越易引起固定凡爾球在實(shí)際工作中與其他壁面碰撞的機(jī)會(huì)，造成球體損壞引起的流體漏失，降低泵效。

（3）泵內(nèi)流體黏度越高固定凡爾球的升程越高所需開啟時(shí)間越短；黏度54 mPs·s時(shí)為臨界黏度，小于54 mPa·s時(shí)固定凡爾球不能完全打開，大于54 mPa·s時(shí)固定凡爾球完全打開，與球罩產(chǎn)生碰撞并有較小幅度的反彈，最終穩(wěn)定在球罩出，泵內(nèi)流體黏度對固定凡爾球的開啟程度影響顯著。

（4）對沖次和泵內(nèi)流體黏度的分析為有桿抽油泵的工作制度優(yōu)化提供了理論參考。