摘"" 要：針對大排量混砂離心泵，采用CFD-DEM耦合計(jì)算方法建立固液兩相流仿真模型，對離心泵內(nèi)兩相流流場進(jìn)行了研究。得到了不同葉輪轉(zhuǎn)速、流量、顆粒粒徑和混砂濃度等因素對離心泵外特性的影響。使用OKA、Archard磨損模型計(jì)算得到蝸殼及葉輪的磨損量，分析不同工況下顆粒運(yùn)動過程及顆粒對過流部件的定性定量的磨蝕規(guī)律，得到平均磨損率隨時間變化曲線，并對不同鍍層的離心泵進(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ī)律;鍍層

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

收稿日期： 2024-08-19

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

作者簡介： 彭俊威（1985-），男，高級工程師，碩士研究生，從事鉆井壓裂新備技術(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

文章編號：1001-3482（2025）01-0044-09

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

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

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

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

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

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

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

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

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

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

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

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

CFD-DEM耦合借助Fluent中的DDPM模型實(shí)現(xiàn)。其質(zhì)量守恒和動量守恒方程基本假設(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。

混砂離心泵工作時固相體積分?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為兩相相對速度，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分別為相間動量交換系數(shù)和顆粒動量變換系數(shù);SDPM，explicit為DPM顯式質(zhì)量源項(xiàng)，kg/（m3·s）。

顆粒的碰撞模型采用Hertz-Mindlin無滑移接觸模型，施加在顆粒上的法向力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為相對速度的法向分量，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為相對速度的切向分量;St為切向剛度。

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

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

τi=-μrFnRiωi（17）

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

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

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

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ù)，對于石英砂顆粒取0.12。sin（α）"""""" 計(jì)算由于重復(fù)變形引起的歸一化侵蝕;C（α）計(jì)算由于沖擊引起的材料移除過程。

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

Q=WFndt（22）

W=K/Hv（23）

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

1.2 模型建立

本文以寶石機(jī)械355.6 mm×304.8 mm×558.8 mm大排量混砂離心泵為研究對象。該離心泵葉片為開式葉片，蝸殼為雙流道結(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對離心泵幾何模型進(jìn)行簡化處理后進(jìn)行流體域抽取，并添加入口延長段和出口延長段，得到最終的離心泵流體域模型如圖1所示。

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

1.2" 網(wǎng)格劃分

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

1.3 邊界條件

流體域模型入口為速度入口，出口為壓力出口。湍流模型采用k-ω SST模型[16]。各部分流體域之間通過Interface邊界連接，設(shè)置葉輪繞-x軸以恒定角速度旋轉(zhuǎn)。設(shè)置時間步長1.0×10-4 s。同時EDEM中設(shè)置葉輪初始相位與流體域相同，并添加葉輪部件以相同角速度運(yùn)動。設(shè)置EDEM 時間步長為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)確性，通過對離心泵手冊中的特性曲線進(jìn)行處理，獲得清水工況下的離心泵外部特性曲線。在工作流量范圍內(nèi)等差選取7組工況進(jìn)行數(shù)值模擬，計(jì)算結(jié)果如圖4所示。在最大效率點(diǎn)附近（流量在15.14～22.71 m3/min），揚(yáng)程及效率計(jì)算誤差均在10%以內(nèi)，說明所建立得到計(jì)算模型滿足計(jì)算精度要求。

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

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

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

3 蝸殼及葉輪磨損研究

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

4 不同鍍層對離心泵磨損的影響

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

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

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

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

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

原離心泵材料為ASTM A532 Class Ⅲ Type A，采用高鉻合金Cr28。本身擁有良好的耐磨性能和耐腐蝕性能。為了增強(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工況為例，對比蝸殼及葉輪的磨損量以及磨損速率，如圖21～22所示。使用碳化鎢鍍層后，離心泵整體磨損量和磨損速率降低，葉輪磨損速率下降28.1%，蝸殼磨損速率下降39.95%。說明使用碳化鎢鍍層可以很好地降低離心泵磨損，提高使用壽命。

5 結(jié)論

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

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

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

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

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

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

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

參考文獻(xiàn)：

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

[2] 袁旭軍，吳漢川. 從我國壓裂市場現(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] 吳漢川，王峻喬，仇黎明. 混砂車吸入排出性能研究[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ù)測[J]. 蘭州理工大學(xué)學(xué)報(bào)，2015，41（2）： 55-60.

[8] 廖姣，賴喜德，廖功磊，等. 基于固液兩相流的離心泵內(nèi)部流場數(shù)值分析[J]. 熱能動力工程，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] 黃凱，劉棟，尤保健，等. 顆粒濃度對離心泵性能及磨損影響的研究[J]. 流體機(jī)械，2022，50（2）： 35-42.

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

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

37（2）： 181-189.

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

36（2）： 98-103.

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

44（11）： 53-60.

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

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

（編輯：韓睿超）