大排量離心泵固液兩相流磨損特性仿真分析

摘"" 要：針對(duì)大排量混砂離心泵，采用CFD-DEM耦合計(jì)算方法建立固液兩相流仿真模型，對(duì)離心泵內(nèi)兩相流流場(chǎng)進(jìn)行了研究。得到了不同葉輪轉(zhuǎn)速、流量、顆粒粒徑和混砂濃度等因素對(duì)離心泵外特性的影響。使用OKA、Archard磨損模型計(jì)算得到蝸殼及葉輪的磨損量，分析不同工況下顆粒運(yùn)動(dòng)過(guò)程及顆粒對(duì)過(guò)流部件的定性定量的磨蝕規(guī)律，得到平均磨損率隨時(shí)間變化曲線，并對(duì)不同鍍層的離心泵進(jìn)行數(shù)值計(jì)算。結(jié)果表明，葉輪主要磨損區(qū)在根部和壓力面中部，轉(zhuǎn)速每增加100 r/min，最大磨損量為葉輪沖擊磨損，增加28.5%;流量增加20%，最大磨損量為蝸殼的沖擊磨損，增加85%。使用碳化鎢鍍層，葉輪和蝸殼磨損速率分別降低28.1%和39.95%。研究結(jié)果可為混砂離心泵參數(shù)設(shè)計(jì)提供指導(dǎo)意義。

關(guān)鍵詞：混砂離心泵;CFD-DEM;磨蝕規(guī)律;鍍層

中圖分類(lèi)號(hào)：TE934.2"""""""" 文獻(xiàn)標(biāo)志碼：A"""""" doi：10.3969/j.issn.1001-3482.2025.01.008

收稿日期： 2024-08-19

基金項(xiàng)目： 中國(guó)石油天然氣集團(tuán)公司“十四五”前瞻性基礎(chǔ)性重大科技項(xiàng)目“大功率智能電驅(qū)壓裂裝備關(guān)鍵技術(shù)研究”（2022DJ4507）。

作者簡(jiǎn)介： 彭俊威（1985-），男，高級(jí)工程師，碩士研究生，從事鉆井壓裂新備技術(shù)研究工作，E-mail：2005pjw@163.com。

Simulation Analysis of The Wearing Characteristics of Solid-Liquid Two-Phase Flow in Centrifugal Pumps

PENG Junwei1，2，ZHOU Qing 1，2，GENG Peng3，DONG Fuqiang1，2，SUN Tianhui4，DENG Wei1，2，LI Sen1，2

（1.CNPC National Petroleum Drilling Equipment Engineering Technology Research Center Company Limited， Baoji 721002，China ; 2. CNPC Baoji Oifield Machinery Co.，Ltd.， Baoji 721002，China; 3.Xi’an Branch，China Petroleum Materials Company Limited ， Xi’an 710018，China;4.China University of Petroleum， Beijing 102249，China）

Abstract： In order to establish a solid-liquid two-phase flow simulation model for the" large-displacement sand mixing centrifugal pump， a CFD-DEM coupling calculation method was employed. As a result， the two-phase flow field within the centrifugal pump was subsequently studied. The impact of varying impeller speed， flow rate， particle size， and sand mixing concentration on the external characteristics of the centrifugal pump was determined. Additionally， the wearing of the volute and impeller was calculated using the OKA and Archard wear models. The qualitative and quantitative abrasion laws governing particle movement and particles on the overflow parts under different working conditions were analyzed. The average wearing rate variation curve with time was obtained， and the centrifugal pumps with different coatings were numerically calculated. The findings indicated that the predominant areas of impeller wear were located at the root and the midpoint of the pressure surface. For each 100 r/min increase in speed， there was a corresponding increase in wear due to impeller impact， which increased by 28.5%. Additionally， there was a 20% increase in flow rate， and the maximum wear amount was the impact wear of the volute， which increased by 85%. The application of a tungsten carbide coating resulted in a reduction in the wear rates of the impeller and volute by 28.1% and 39.95%， respectively. This finding provides valuable insight that can inform the design of sand-mixing centrifugal pumps.

Key words： sand-mixing centrifugal pump; CFD-DEM; abrasion laws; coating

文章編號(hào)：1001-3482（2025）01-0044-09

大排量混砂離心泵作為壓裂液輸送設(shè)備，在油田的壓裂作業(yè)中起到重要作用[1-2]。壓裂作業(yè)環(huán)境為固液兩相流工況，壓裂液攜帶大量支撐劑固體顆粒，運(yùn)行環(huán)境惡劣。內(nèi)部葉輪、蝸殼等長(zhǎng)期遭受堅(jiān)硬固體顆粒的沖刷，水力部件壁面磨損嚴(yán)重，從而導(dǎo)致泵水力性能明顯下降，輸送系統(tǒng)壽命遠(yuǎn)小于預(yù)期[3-4]。

為降低固液兩相流工況下的離心泵磨蝕，國(guó)內(nèi)外學(xué)者展開(kāi)了深入研究。其中數(shù)值模擬法具有邊界條件容易控制、耗費(fèi)少、節(jié)約成本，用時(shí)短等優(yōu)點(diǎn)，并且可反復(fù)模擬，優(yōu)化設(shè)計(jì)的過(guò)程也相對(duì)簡(jiǎn)單。因此在水力機(jī)械的研究中數(shù)值模擬有著重要的價(jià)值。目前對(duì)固液兩相流的數(shù)值模擬方法分為雙流體法和離散相法。

雙流體法是將固體顆粒視為與液體相同的連續(xù)介質(zhì)。其中，趙斌娟等[5]研究顆粒在離心泵內(nèi)的分布規(guī)律，研究表明隨著粒徑的減小，顆粒逐漸從葉片壓力面向吸力面流動(dòng)。Yifang等[6]研究了離心泵內(nèi)顆粒的分布規(guī)律，結(jié)果表明隨著顆粒體積分?jǐn)?shù)增大，顆粒逐漸向葉片吸力面聚集。程效銳等[7]對(duì)雙吸離心泵的磨損特性進(jìn)行了研究，結(jié)果表明增大粒徑會(huì)增加磨損，顆粒質(zhì)量濃度影響大于顆粒直徑影響，且質(zhì)量濃度變化會(huì)導(dǎo)致磨損位置的變化。廖姣[8]對(duì)粒徑0.05 mm，體積分?jǐn)?shù)為10%的顆粒下工況進(jìn)行了仿真分析，結(jié)果表明分流葉片可以減小葉片出口處的壓力脈動(dòng)，減少顆粒聚集現(xiàn)象的產(chǎn)生，減輕葉片出口處的摩擦。

離散相法是在拉格朗日坐標(biāo)系中采用獨(dú)立的方程來(lái)求解固相的運(yùn)動(dòng)，其中DEM法通過(guò)軟球模型計(jì)算顆粒的運(yùn)動(dòng)，能準(zhǔn)確模擬固液兩相流工況下離心泵內(nèi)顆粒的運(yùn)動(dòng)。Kafui[9]提出將DEM與CFD進(jìn)行耦合計(jì)算，以提高計(jì)算效率，擴(kuò)展了其在兩相流領(lǐng)域內(nèi)的應(yīng)用范圍。Tang等[10]采用CFD-DEM耦合方法研究了離心泵內(nèi)固體顆粒的運(yùn)動(dòng)特性，結(jié)果表明小顆粒具有較寬的速度范圍和速度峰值，而大顆粒具有較大的接觸力，餅狀顆粒的接觸次數(shù)最多，球形顆粒的接觸次數(shù)最少。黃凱[11]使用CFD-EDM耦合方法，監(jiān)測(cè)過(guò)流部件的磨損率和顆粒的運(yùn)動(dòng)，發(fā)現(xiàn)隨著質(zhì)量濃度的增加，揚(yáng)程和效率減小，葉輪的磨損率增大，葉片比前后蓋板更易被磨損，葉片的磨損區(qū)域由點(diǎn)狀衍變?yōu)槠瑺?。王田田等[12-13]對(duì)球形、類(lèi)球形、三角形三種顆粒在不同粒徑（1～3 mm）、體積分?jǐn)?shù)（2%～4%）下進(jìn)行了仿真計(jì)算，結(jié)果表明顆粒直徑增大，葉片磨損增加，磨損位置逐漸向葉片頭部轉(zhuǎn)移，顆粒體積分?jǐn)?shù)增加，磨損位置不變，磨損程度增大，顆粒尖銳程度增大也會(huì)加劇磨損。

現(xiàn)今中外學(xué)者借助兩相流理論和經(jīng)驗(yàn)?zāi)Ｐ?，開(kāi)展對(duì)于不同規(guī)格離心泵磨損特性的研究，通過(guò)固液兩相流場(chǎng)分析磨損機(jī)理[14-15]，但針對(duì)大排量高砂比泵的相關(guān)研究較少。本文使用CFD-DEM耦合計(jì)算的方法，針對(duì)石油行業(yè)的大排量混砂離心泵長(zhǎng)期遭受堅(jiān)硬固體顆粒沖刷的環(huán)境下，開(kāi)展對(duì)固液兩相流高砂比工況離心泵磨損機(jī)理的研究，研究離心泵性能和磨損趨勢(shì)變化，指導(dǎo)大排量混砂設(shè)備離心泵設(shè)計(jì)，提升壽命和經(jīng)濟(jì)性選型。

1 CFD-DEM兩相流計(jì)算模型

1.1 數(shù)學(xué)模型

壓裂液作為不可壓縮流體，在數(shù)值模擬過(guò)程中要保證質(zhì)量守恒和動(dòng)量守恒。連續(xù)性方程描述了流體質(zhì)量守恒的關(guān)系，即在單位時(shí)間內(nèi)流入單位體積的質(zhì)量等于流出單位體積的質(zhì)量，公式為：

+?·（ρ[→][v]）=Sm（1）

動(dòng)量守恒方程描述了兩相流體在流動(dòng)中的力學(xué)特性，包括流體的速度分布、流速和流動(dòng)阻力等，公式為：

（ρ[→][v]）+?·（ρ[→][v][→][v]）=-?p+?·（[=][τ]）+ρ[ɡ][-]+[→][F]（2）

式中：p為靜壓，[=][τ]為應(yīng)力張量，ρ[ɡ][-]為自重，[→][F]為相互作用力。

CFD-DEM耦合借助Fluent中的DDPM模型實(shí)現(xiàn)。其質(zhì)量守恒和動(dòng)量守恒方程基本假設(shè)是在顆粒離散相體積分?jǐn)?shù)足夠低情況，公式為：

"""""" +?·（ρ[→][v]）=SDPM + Sother（3）

+?·（ρ[→][v][→][v]）=-?p+?·τ+ρ[ɡ][→]+

"""" [→][FDPM] +[→][Fother]（4）

式中：ρ為連續(xù)相密度，kg/m3;[→][v]為連續(xù)相速度矢量，m/s;SDPM和Sother為DPM離散相及其他附加的質(zhì)量源項(xiàng)，kg/（m3·s）;ρ[ɡ][→]為重力，N/m3;[→][FDPM] 和[→][Fother]為DPM離散相及其他附加的體積力，N/m3。

混砂離心泵工作時(shí)固相體積分?jǐn)?shù)較高，為了克服上述方程限制，引入拓展體積分?jǐn)?shù)p解釋顆粒相：

（αpρp）+?（αpρp[→][vp]）=[q=1][nphas][∑]（[m][.]qp-[m][.]pq）（5）

（αpρp[→][vp]）+?（αpρp[→][vp][→][vp]）=

αpρp[ɡ][→]+Fvm，lift，user-αp?p+?·[αpμp（?[→][vp]+?[→][vp][T]）]+

KDP（[→][v]DPM - [→][vp]）+SDPM，explicit+

[q=1][nphas][∑]（[→][K]qp（[→][vq]-[→][vp]）-[m][.]qp[→][v]pq+[m][.]pq[→][v]pq）（6）

式中：αp為p相體積分?jǐn)?shù);ρp為p相密度，kg/m3;μp為p相黏度，kg/m2·s;[→][vp]和[→][vq]分別為p相和q相速度矢量，m/s;[→][v]pq為兩相相對(duì)速度，m/s;[m][.]qp和[m][.]pq為p相和q相之間的質(zhì)量傳遞，kg/（m3·s）;Fvm，lift，user為虛擬質(zhì)量力、升力及其他作用力的合力，N/m3;[→][v]DPM表示離散相顆粒的平均速度，m/s;[→][K]qp和KDPM分別為相間動(dòng)量交換系數(shù)和顆粒動(dòng)量變換系數(shù);SDPM，explicit為DPM顯式質(zhì)量源項(xiàng)，kg/（m3·s）。

顆粒的碰撞模型采用Hertz-Mindlin無(wú)滑移接觸模型，施加在顆粒上的法向力Fn模型由法向接觸力和阻尼力組成，為法向重疊量δn的函數(shù)：

"""""""""""""""" Fn=Fcn+Fdn=Eδn （7）

其中：

"""""""""""""""" =+（8）

"""""""""""""nbsp;"" =（9）

式中：Fcn為法向接觸力;Fdn為法向阻尼力;E為等效彈性模量;R為等效顆粒半徑;E1、E2分別為兩顆粒的彈性模量;R1、R2分別為兩顆粒的接觸半徑。

法向力中阻尼力Fdn表達(dá)式為：

Fdn=-2β（10）

其中：

Fdn==（11）

Fdn=β=（12）

Fdn=Sn=2E（13）

式中：m*為當(dāng)量質(zhì)量;m1、m2分別為兩顆粒的質(zhì)量，kg;vn為相對(duì)速度的法向分量，m/s;e為恢復(fù)系數(shù);Sn為法向剛度，N/m。

施加在顆粒上的切向力Ft由切向接觸力和阻尼力組成，為切向重疊量δt的函數(shù)：

Ft=Fct+Fdt=-St δt（14）

其中：

St=8G*（15）

式中：G*為當(dāng)量剪切模量，MPa。

切向力中阻尼力Fdt表達(dá)式為：

Fdt=-2β（16）

式中：vt為相對(duì)速度的切向分量;St為切向剛度。

顆粒之間的摩擦力μSFn會(huì)影響切向力，其中μS為靜摩擦因數(shù)。

顆粒的滾動(dòng)摩擦模型選擇標(biāo)準(zhǔn)滾動(dòng)摩擦模型，該模型通過(guò)在接觸表面添加力矩考慮滾動(dòng)摩擦情況：

τi=-μrFnRiωi（17）

式中：μr為滾動(dòng)摩擦因數(shù);Ri為顆粒質(zhì)心與接觸點(diǎn)的距離;ωi是接觸點(diǎn)處的角速度。

本研究中使用OKA磨損模型計(jì)算顆粒對(duì)過(guò)流部件的沖擊磨損，使用Archard模型計(jì)算顆粒對(duì)過(guò)流部件的摩擦磨損。

OKA磨損模型由Oka和Yoshida在2005年提出。主要用于計(jì)算顆粒沖擊設(shè)備表面時(shí)發(fā)生的侵蝕性磨損。該模型通過(guò)顆粒尺寸、顆粒沖擊速度和顆粒沖擊角度等參數(shù)的函數(shù)關(guān)系預(yù)測(cè)了顆粒沖擊引起的體積損失。

DωA=ɡ（α）E（α）mp（18）

E（α）=65W-" （） （）0.19（19）

""" ɡ（α）=sin（α）"""""" C（α）（20）

C（α）=1+Hv（1-sin（α））（21）

式中：Dω為沖蝕深度，m;A為沖擊面積，m2;ɡ（α）為歸一化的沖蝕角函數(shù);E（α）為單位質(zhì)量下的磨損體積，m3;mp為顆粒質(zhì)量，kg;W為材料磨損常數(shù);v為顆粒的沖擊速度，m/s;Hv為材料的維氏硬度，GPa;D為顆粒直徑，m;k1為試驗(yàn)推導(dǎo)的經(jīng)驗(yàn)系數(shù)，對(duì)于石英砂顆粒取0.12。sin（α）"""""" 計(jì)算由于重復(fù)變形引起的歸一化侵蝕;C（α）計(jì)算由于沖擊引起的材料移除過(guò)程。

Archard磨損模型由Archard于1940年提出。該公式通過(guò)考慮材料表面微觀突起和顆粒之間的接觸，描述了材料摩擦磨損程度與顆粒所作摩擦功之間的關(guān)系，方程定義為：

Q=WFndt（22）

W=K/Hv（23）

式中：Q為去除材料的體積，m3;Fn為顆粒沖擊力，N;dt為顆粒移動(dòng)的切向距離，m;W為磨損常數(shù);K為無(wú)量綱常數(shù);Hv為材料的維氏硬度，GPa。

1.2 模型建立

本文以寶石機(jī)械355.6 mm×304.8 mm×558.8 mm大排量混砂離心泵為研究對(duì)象。該離心泵葉片為開(kāi)式葉片，蝸殼為雙流道結(jié)構(gòu)。設(shè)計(jì)最大流量Q=29 m3/min，額定轉(zhuǎn)速1 200 ～1 400 r/min，清水條件下額定輸入功率≥370 kW，吸入口通直徑355.6 mm，排出口通徑304.8 mm。

使用SpaceClaim對(duì)離心泵幾何模型進(jìn)行簡(jiǎn)化處理后進(jìn)行流體域抽取，并添加入口延長(zhǎng)段和出口延長(zhǎng)段，得到最終的離心泵流體域模型如圖1所示。

EDEM中導(dǎo)入簡(jiǎn)化后的離心泵固體域模型，將顆粒工廠設(shè)置在入口延長(zhǎng)段靠近葉輪處，防止顆粒生成時(shí)流出流體區(qū)域?qū)е碌挠?jì)算發(fā)散。顆粒模型使用多個(gè)球形顆粒堆積模擬圓度為0.6的石英砂顆粒。EDEM中離心泵及顆粒模型如圖2所示。由于支撐劑顆粒體積較小，采用真實(shí)尺度對(duì)顆粒進(jìn)行建模會(huì)導(dǎo)致顆粒數(shù)量龐大，大大降低了計(jì)算效率。因此在仿真中使用了曳力放大的方法，將顆粒直徑成倍數(shù)放大后，在Fluent中的流體作用力也成倍進(jìn)行放大，使用大顆粒模擬多個(gè)小顆粒，降低顆粒個(gè)數(shù)提高計(jì)算效率。顆粒放大倍數(shù)根據(jù)不同工況進(jìn)行調(diào)整。

1.2" 網(wǎng)格劃分

流體計(jì)算模型由入口延長(zhǎng)段、葉輪流體域、蝸殼流體域、出口延長(zhǎng)段四部分組成。為保證網(wǎng)格質(zhì)量，縮短計(jì)算時(shí)間，本文采用六面體+多面體的方法進(jìn)行劃分。流道中心區(qū)域使用規(guī)則的六面體網(wǎng)格，在壁面附近使用多面體網(wǎng)格以適應(yīng)葉輪和蝸殼的曲面及倒角。模型網(wǎng)格劃分如圖3所示。

1.3 邊界條件

流體域模型入口為速度入口，出口為壓力出口。湍流模型采用k-ω SST模型[16]。各部分流體域之間通過(guò)Interface邊界連接，設(shè)置葉輪繞-x軸以恒定角速度旋轉(zhuǎn)。設(shè)置時(shí)間步長(zhǎng)1.0×10-4 s。同時(shí)EDEM中設(shè)置葉輪初始相位與流體域相同，并添加葉輪部件以相同角速度運(yùn)動(dòng)。設(shè)置EDEM 時(shí)間步長(zhǎng)為1.0×10-6 s。

仿真中壓裂液密度和黏度為1 000 kg/m3和40 mPa·s。支撐劑密度為2 650 kg/m3[17]。入口流速及顆粒生成速率根據(jù)工況條件進(jìn)行調(diào)整。

1.4 模型驗(yàn)證

為驗(yàn)證所建立的計(jì)算模型準(zhǔn)確性，通過(guò)對(duì)離心泵手冊(cè)中的特性曲線進(jìn)行處理，獲得清水工況下的離心泵外部特性曲線。在工作流量范圍內(nèi)等差選取7組工況進(jìn)行數(shù)值模擬，計(jì)算結(jié)果如圖4所示。在最大效率點(diǎn)附近（流量在15.14～22.71 m3/min），揚(yáng)程及效率計(jì)算誤差均在10%以?xún)?nèi)，說(shuō)明所建立得到計(jì)算模型滿(mǎn)足計(jì)算精度要求。

2 泵內(nèi)顆粒運(yùn)動(dòng)過(guò)程分析

以流量18.92 m3/min、顆粒尺寸30目、顆粒體積分?jǐn)?shù)10%、葉輪轉(zhuǎn)速1 200 r/min為例，提取顆粒進(jìn)入至流出葉輪區(qū)域過(guò)程的顆粒運(yùn)動(dòng)速度分布云圖，如圖5所示。顆粒進(jìn)入離心泵后0.03 s時(shí)，在離心力作用下快速向徑向擴(kuò)散，主要包裹了葉輪靠近中心的葉片根部邊緣區(qū)域。0.05 s時(shí)顆粒運(yùn)動(dòng)到葉輪流道中部位置，此時(shí)葉輪壓力面根部部分位置與顆粒碰撞。顆粒進(jìn)入離心泵后約0.07 s流出葉輪。

顆粒進(jìn)入至流出蝸殼區(qū)域過(guò)程的顆粒運(yùn)動(dòng)速度分布云圖，如圖6所示。顆粒進(jìn)入蝸殼后在慣性力作用下沿蝸殼流道內(nèi)壁向出口運(yùn)動(dòng)。分隔板將蝸殼分為上下兩個(gè)流道。0.15 s時(shí)分隔板內(nèi)側(cè)流道內(nèi)顆粒運(yùn)動(dòng)已經(jīng)基本達(dá)到穩(wěn)定，此時(shí)外側(cè)顆粒運(yùn)動(dòng)至蝸殼增壓管根部附近。0.20 s后，蝸殼內(nèi)部顆粒運(yùn)動(dòng)基本達(dá)到穩(wěn)定。顆粒進(jìn)入蝸殼后向外側(cè)偏移。

3 蝸殼及葉輪磨損研究

3.1 蝸殼及葉輪磨損區(qū)域

不同工況下的泵內(nèi)顆粒運(yùn)動(dòng)規(guī)律及分布基本相同，以流量為18.92 m3/min、顆粒尺寸30目、顆粒體積分?jǐn)?shù)10%、葉輪轉(zhuǎn)速1 200 r/min工況為例，提取了葉輪及蝸殼的摩擦磨損和沖擊磨損云圖如圖7所示。

蝸殼及葉輪表面磨損可分為三個(gè)區(qū)域。

蝸殼A區(qū)域?yàn)槲仛?nèi)流道底部。該區(qū)域不存在分隔板，顆粒運(yùn)動(dòng)速度較低，摩擦磨損和沖擊磨損量均比較少;B區(qū)域?yàn)榉指舭甯课恢?，從葉輪流出的高速顆粒在該位置向蝸殼內(nèi)外兩個(gè)流道分散，主要產(chǎn)生沖擊磨損;C區(qū)域?yàn)槲仛こ隹谠鰤汗芨?，分隔板上下兩個(gè)流道均出現(xiàn)磨損現(xiàn)象。由于分隔板下方流道的顆粒更靠近葉輪，運(yùn)動(dòng)速度更大，因此分隔板下方所受到的沖擊磨損和摩擦磨損相對(duì)于上方流道都更加嚴(yán)重。

葉輪A區(qū)域?yàn)槿~片根部邊緣，由于顆粒進(jìn)入離心泵時(shí)向軸向運(yùn)動(dòng)同時(shí)在離心力作用下向徑向擴(kuò)散，因此該位置同時(shí)發(fā)生了比較嚴(yán)重的摩擦磨損和沖擊磨損。B區(qū)域?yàn)槿~輪底面，顆粒進(jìn)入葉輪區(qū)域后，軸向分速度逐漸降低，運(yùn)動(dòng)至葉輪底部時(shí)主要為離心力作用下產(chǎn)生的徑向速度，因此在底部主要發(fā)生摩擦磨損，沖擊磨損量很低。葉輪C區(qū)域?yàn)槿~片壓力面中心區(qū)域，由于葉輪轉(zhuǎn)動(dòng)至蝸殼隔水角附近時(shí)顆粒速度低，在該位置葉片更容易與顆粒發(fā)生碰撞，因此主要產(chǎn)生沖擊磨損。

3.2 不同條件蝸殼及葉輪磨損速率

不同轉(zhuǎn)速下蝸殼的沖擊磨損及摩擦磨損如圖8～9所示。在其他條件相同時(shí)，葉輪轉(zhuǎn)速的提高增加了蝸殼分隔板下表面的沖擊磨損以及蝸殼下方流道的摩擦磨損，主要磨損區(qū)域未發(fā)生明顯變化。

圖10為不同轉(zhuǎn)速下的蝸殼沖擊磨損量及摩擦磨損量隨時(shí)間變化曲線。在其他條件相同時(shí)，轉(zhuǎn)速每增加100 r/min時(shí)，蝸殼的沖擊磨損量和摩擦磨損量分別平均增加13%和18.5%。

不同轉(zhuǎn)速下葉輪的沖擊磨損及摩擦磨損如圖11～12所示。在其他條件相同時(shí)，葉輪轉(zhuǎn)速的對(duì)葉輪壓力面的沖擊磨損影響較大，同時(shí)葉輪底面所受的摩擦磨損量也呈上升趨勢(shì)。

圖13為不同轉(zhuǎn)速下的葉輪沖擊磨損量及摩擦磨損量隨時(shí)間變化曲線。葉輪轉(zhuǎn)速的增加使顆粒獲得了更大的動(dòng)能，導(dǎo)致磨損更加嚴(yán)重。葉輪轉(zhuǎn)速的增加對(duì)葉輪的沖擊磨損影響更大。在其他條件相同時(shí)，轉(zhuǎn)速每增加100 r/min時(shí)，葉輪的沖擊磨損量和摩擦磨損量分別平均增加28.5%和8%。轉(zhuǎn)速主要影響了葉輪受到的沖擊磨損。

不同流量下蝸殼沖擊磨損及摩擦磨損如圖14～圖15所示。在其他條件相同時(shí)，流量的提高使得蝸殼所受摩擦磨損和沖擊磨損提高，主要磨損區(qū)域面

積增加，沖擊磨損和摩擦磨損范圍均向離心泵分隔板根尾部延伸。

圖16為不同流量下的蝸殼沖擊磨損量及摩擦磨損量隨時(shí)間變化曲線。在其他條件相同時(shí)，入口流量每增加20%，蝸殼的沖擊磨損量和摩擦磨損量分別增加85%和53%。

不同流量下葉輪沖擊磨損及摩擦磨損如圖17～18所示。在其他條件相同時(shí)，隨流量的提高葉輪主要沖擊磨損位置未發(fā)生明顯變化，葉片根部邊緣和底部所受的摩擦磨區(qū)域增加。

圖19為不同流量下的葉輪沖擊磨損量及摩擦磨損量隨時(shí)間變化曲線。在其他條件相同時(shí)，流量的增加使葉輪的沖擊磨損量上升幅度程增大趨勢(shì)、而摩擦磨損上升幅度減小。入口流量每增加3.785 m3/min，葉輪的沖擊磨損量和摩擦磨損量分別增加31.5%和38.5%。流量的增加對(duì)蝸殼的影響程度要大于對(duì)葉輪的影響程度。

由上述結(jié)果可知，在離心泵運(yùn)行0.3 s后，蝸殼及葉輪的磨損量基本呈線性增長(zhǎng)趨勢(shì)。因此對(duì)0.3 s后磨損數(shù)據(jù)進(jìn)行線性擬合處理，并將沖擊磨損量與摩擦磨損量進(jìn)行疊加，得到不同工況下蝸殼和葉輪的總體磨損速率，如圖20所示。由圖20可知，蝸殼及葉輪主要受到?jīng)_擊磨損影響，相同條件下摩擦磨損速率比沖擊磨損速率小得多。蝸殼磨損速率主要受流量影響，葉輪磨損速率主要受轉(zhuǎn)速影響。同時(shí)葉輪的磨損量要比蝸殼大4至5倍。

4 不同鍍層對(duì)離心泵磨損的影響

由上可知，蝸殼與葉輪所受的沖擊磨損要大于摩擦磨損。根據(jù)式（14）～（17），對(duì)于大排量離心泵的顆粒磨損受以下因素影響：

1） 顆粒本身物理屬性，如質(zhì)量、直徑。

2） 輸送參數(shù)引起的顆粒運(yùn)動(dòng)屬性，如顆粒運(yùn)動(dòng)速度、沖擊角度等。

3） 離心泵的材料屬性，如硬度及材料常數(shù)等。

顆粒本身物理屬性及輸送參數(shù)受離心泵作業(yè)指標(biāo)影響。為降低磨損，考慮以更換鍍層的方式提高離心泵過(guò)流部件表面硬度，提高蝸殼及葉輪的使用壽命。

原離心泵材料為ASTM A532 Class Ⅲ Type A，采用高鉻合金Cr28。本身?yè)碛辛己玫哪湍バ阅芎湍透g性能。為了增強(qiáng)其耐磨性能考慮在其表面增加碳化鎢涂層以提高耐磨性能。熱噴涂碳化鎢（WC）涂層硬度通常在1 200～1 600 HV之間，本研究中取1 400 HV。以流量為18.92 m3/min、顆粒尺寸30目、顆粒質(zhì)量分?jǐn)?shù)10%、葉輪轉(zhuǎn)速1 200 r/min工況為例，對(duì)比蝸殼及葉輪的磨損量以及磨損速率，如圖21～22所示。使用碳化鎢鍍層后，離心泵整體磨損量和磨損速率降低，葉輪磨損速率下降28.1%，蝸殼磨損速率下降39.95%。說(shuō)明使用碳化鎢鍍層可以很好地降低離心泵磨損，提高使用壽命。

5 結(jié)論

本文使用CFD-DEM耦合計(jì)算方法，計(jì)算了不同固液兩相流工況對(duì)大排量混砂離心泵顆粒運(yùn)動(dòng)過(guò)程及磨損特性進(jìn)行了研究。結(jié)論如下：

1） 蝸殼雙流道結(jié)構(gòu)會(huì)導(dǎo)致外流道顆粒運(yùn)動(dòng)滯后，且顆粒進(jìn)入蝸殼后向外側(cè)偏移。

2） 蝸殼主要磨損區(qū)域?yàn)殡p流道分隔板以及增壓管根部。葉輪主要磨損區(qū)域發(fā)生在葉片根部邊及葉片壓力面中部。

3） 在其他條件相同時(shí)，葉輪轉(zhuǎn)速每上升10 r/min，蝸殼的沖擊磨損量和摩擦磨損量分別平均增加13%和18.5%。葉輪的沖擊磨損量和摩擦磨損量分別平均增加28.5%和8%。

4） 在其他條件相同時(shí)，入口流量每增加20%，蝸殼的沖擊磨損量和摩擦磨損量分別增加85%和53%。葉輪的沖擊磨損量和摩擦磨損量分別增加31.5%和38.5%

5） 離心泵運(yùn)行0.3 s后，蝸殼及葉輪的磨損量基本呈線性增長(zhǎng)趨勢(shì)。相同條件下離心泵的沖擊磨損速率高于摩擦磨損速率比。蝸殼磨損速率主要受流量影響，葉輪磨損速率主要受轉(zhuǎn)速影響。同時(shí)葉輪的磨損量要比蝸殼大4～5倍。

6） 使用碳化鎢鍍層使葉輪及蝸殼的磨損速率分別降低28.1%和39.95%，可以很好地降低離心泵磨損，提高使用壽命。

參考文獻(xiàn)：

[1] 李文維. 小排量混砂增壓泵車(chē)的研制[J]. 科技與創(chuàng)新， 2020，（10）： 123-124.

[2] 袁旭軍，吳漢川. 從我國(guó)壓裂市場(chǎng)現(xiàn)狀談大型壓裂機(jī)組的研制[J]. 石油天然氣學(xué)報(bào)，2010， 32（3）： 383-385.

[3] 劉棟，倪子建，黃凱，等. 基于CFD-DEM的離心泵葉輪磨損特性研究和優(yōu)化設(shè)計(jì)[J]. 流體機(jī)械， 2022， 50（9）： 26-32.

[4] 吳漢川，王峻喬，仇黎明. 混砂車(chē)吸入排出性能研究[J]. 石油機(jī)械，2013，41（3）： 92-95.

[5] 趙斌娟，王澤. 離心泵葉輪內(nèi)流數(shù)值模擬的現(xiàn)狀和展望[J]. 農(nóng)機(jī)化研究，2002（3）： 49-52.

[6] SHI Yifang，Tang Lingfeng，TAN Yinwu，et al. Optimization of the Structural Parameters of a Plastic Centrifugal Pump[J]. Fluid Dynamics amp;amp;Materials Processing，2022，

18（3）：713-736.

[7] 程效銳，董富弟，楊從新，等. 基于顆粒摩擦和碰撞模型的離心泵葉片磨損預(yù)測(cè)[J]. 蘭州理工大學(xué)學(xué)報(bào)，2015，41（2）： 55-60.

[8] 廖姣，賴(lài)喜德，廖功磊，等. 基于固液兩相流的離心泵內(nèi)部流場(chǎng)數(shù)值分析[J]. 熱能動(dòng)力工程，2017，32（5）： 95-99，

139-140.

[9] Kafui D k，Johnson S，Thornton C，et al.Parallelization of a Lagrangian-Eulerian DEM/CFD code for application to fluidized beds[J].Powder Technology，2011，207（1/3）： 270-278.

[10] TANG cheng，YANG Youchao，LIU Pengzhan，et al. Prediction of Abrasive and Impact Wear Due to Multi-Shaped Particles in a Centrifugal Pump via CFD-DEMCoupling Method[J]. Energies，2021，14（9）：1-15.

[11] 黃凱，劉棟，尤保健，等. 顆粒濃度對(duì)離心泵性能及磨損影響的研究[J]. 流體機(jī)械，2022，50（2）： 35-42.

[12] 王田田，支嘉才，楊具瑞，等. 基于計(jì)算流體動(dòng)力學(xué)與離散元法的離心泵內(nèi)流場(chǎng)及磨損的數(shù)值模擬[J]. 湖南農(nóng)業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)版），2022，48（2）： 235-241.

[13] 王田田，池曉清，楊具瑞，等. 顆粒特性對(duì)離心泵葉輪葉片的磨損分析[J]. 水動(dòng)力學(xué)研究與進(jìn)展A輯，2022，

37（2）： 181-189.

[14] 趙偉國(guó)，鄭英杰，劉宜，等. 沙粒體積分?jǐn)?shù)對(duì)離心泵磨損特性影響的數(shù)值分析[J]. 排灌機(jī)械工程學(xué)報(bào)， 2018，

36（2）： 98-103.

[15] 汪家瓊，蔣萬(wàn)明，孔繁余，等. 固液兩相流離心泵內(nèi)部流場(chǎng)數(shù)值模擬與磨損特性[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào)，2013，

44（11）： 53-60.

[16] 李曉俊，袁壽其，潘中永，等. 離心泵邊界層網(wǎng)格的實(shí)現(xiàn)及應(yīng)用評(píng)價(jià)[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2012，28（20）： 67-72.

[17] 孫銳，楊釗.高速通道壓裂過(guò)程中支撐劑運(yùn)移規(guī)律研究[J]. 河南科學(xué)，2023，41（4）：516-523.

（編輯：韓睿超）