基于SPH法的套管沖蝕特性及結(jié)構(gòu)優(yōu)化研究

針對固井滑套壓裂中沖蝕特性及結(jié)構(gòu)優(yōu)化等問題，基于光滑粒子動力學(xué)方法（SPH）建立球形砂粒沖擊金屬表面的SPH數(shù)值模型，模擬多個砂粒沖擊套管材料42CrMo的微觀過程，分析砂粒的撞擊行為與沖蝕之間的關(guān)系；同時利用CFD模擬方法與Oka沖蝕模型結(jié)合現(xiàn)場實際工況，對不同端口特性的套管內(nèi)壁進(jìn)行沖蝕磨損的數(shù)值計算。研究結(jié)果表明：SPH法能較好地表現(xiàn)出砂粒在撞擊套管材料時的微觀變形，有助于研究微觀下材料的變形影響因素，且由SPH法模擬得到的無量綱沖蝕速率隨沖擊角度變化的模擬數(shù)據(jù)與Oka沖蝕公式中沖擊角函數(shù)趨勢相符。端口特性參數(shù)對套管內(nèi)壁的沖蝕有著較大的影響，具體表現(xiàn)在：隨著端口數(shù)量、端口面積的增加，平均沖蝕速率均有一定程度的減小，但存在著邊際遞減效應(yīng)；在研究的4種端口形狀中，膠囊形端口的平均沖蝕速率最小。研究結(jié)果可為套管端口優(yōu)化設(shè)計提供參考。

水力壓裂；光滑粒子流體動力學(xué)；微觀沖擊模型；套管沖蝕；端口特性參數(shù)

TE28

A

DOI： 10.12473/CPM.202401006

SPH-Based Research on Erosion Characteristics

and Structural Optimization of Casing

Zhang Lei" Wang Chenrui

（College of Mechanical and Transportation Engineering， China University of Petroleum （Beijing））

For purpose of erosion characteristic research and structural optimization in cementing sliding sleeve fracturing， based on smoothed particle hydrodynamics （SPH）， a SPH numerical model of spherical sand particle impacting metal surface was built. The micro process of multiple sand particles impacting casing material 42CrMo was simulated， and the relationship between sand particle impact behavior and erosion was analyzed. Meanwhile， combined with field conditions， the computational fluid dynamics （CFD） simulation method and the Oka erosion model were used for numerical calculation on the erosion wear of casing inwall with different port characteristics. The research results show that the SPH method can effectively reflect the microscopic deformation of sand particles when impacting casing materials， which is helpful for studying the factors affecting material deformation at the microscopic level. The simulation data of variation of dimensionless erosion rate with impact angle obtained by SPH method is consistent with the trend of impact angle function in the Oka erosion formula. The port characteristic parameters have a great influence on the erosion of casing inwall. As the number and area of ports increase， the average erosion rate decreases to some extent， but there is a marginal decline effect. Among the four port shapes studied， the capsule shaped port has the lowest average erosion rate. The research results provide references for the optimum design of casing port.

hydraulic fracturing;smooth particle hydrodynamics;micro impact model;casing erosion;port characteristic parameters

基金項目：中國石油天然氣股份有限公司科研項目“可開關(guān)固井滑套壓裂起裂機理及結(jié)構(gòu)優(yōu)化研究”（21GJ-FW-010）。

0" 引" 言

張磊，等：基于SPH法的套管沖蝕特性及結(jié)構(gòu)優(yōu)化研究

固井滑套壓裂技術(shù)是當(dāng)前的主流水平井壓裂技術(shù)，相較于其他壓裂技術(shù)有著高效率、低成本等經(jīng)濟技術(shù)優(yōu)勢［1］。在固井滑套壓裂過程中，含有砂?；蚱渌蝿┑膲毫岩涸诹鲃訒r會對套管端口及端口周圍內(nèi)壁產(chǎn)生不同程度的沖蝕損傷。合理的套管端口結(jié)構(gòu)設(shè)計可以改變端口附近流場，減輕含砂壓裂液對端口及其周圍的沖蝕磨損，達(dá)到安全壓裂以及預(yù)防相關(guān)部件失效的目的。沖蝕是指材料受到微小而松散的流動顆粒沖擊造成表面破壞的磨損現(xiàn)象，廣泛存在于石油、機械等各個領(lǐng)域［2］。利用光滑粒子動力學(xué)（Smoothed Particle Hydro-dynamics， SPH）方法可以對材料在微觀尺度上的撞擊特性進(jìn)行研究。

國內(nèi)外學(xué)者對沖蝕機理與壓裂工具優(yōu)化進(jìn)行大量研究。多數(shù)學(xué)者指出［3-5］，管壁沖蝕的重要因素是高速顆粒撞擊，并利用數(shù)值模擬驗證了當(dāng)部件幾何形狀、材料物性一定時，沖蝕速率取決于顆粒撞擊的速度、沖擊角分布等因素。王川等［6］總結(jié)了造成滑套沖蝕中流場與顆粒的運動特點。魏遼等［7］通過試驗與數(shù)值模擬得到了沖蝕量與砂比和排量之間的關(guān)系。WONG C.Y.等［8］對帶有中心孔的平板進(jìn)行沖蝕試驗，研究石英砂粒對孔邊緣的沖蝕現(xiàn)象，結(jié)果表明，復(fù)雜的幾何模型對沖蝕結(jié)果會有影響。XU Y.Q.、向正新等［9-10］對井下工具進(jìn)行研究，結(jié)果表明，可以利用CFD計算其內(nèi)部流場分布并對井下工具進(jìn)行結(jié)構(gòu)優(yōu)化。M.PARSI等［11］提出在沖蝕過程中存在擠壓、切削等多種機制。DONG X.W.等［12-13］與M.TAKAFFOLI等［14］驗證SPH方法在單顆粒撞擊目標(biāo)材料中的應(yīng)用，其模擬結(jié)果與試驗結(jié)果較為符合，驗證了SPH法在模擬大形變中的優(yōu)勢。M.AZIMIAN等［15］與V.HADAVI等［16］利用FEM法與SPH法分別對角粒子沖擊機理進(jìn)行研究，結(jié)果表明，沖擊粒子帶有尖角時，會強化沖蝕作用，且沖蝕作用與沖擊角與方位角均有關(guān)聯(lián)。廖翔云等［17］通過SPH方法建立數(shù)值模型，研究了組合射流對巖石的沖蝕效果。王方祥等［18］基于SPH與有限元的耦合方法，模擬了壓裂液返排期間放噴油嘴的磨損過程。

本文在上述研究的基礎(chǔ)上，將SPH方法與CFD計算工作相結(jié)合，考慮端口特性參數(shù)，圍繞固井滑套壓裂過程中的套管沖蝕現(xiàn)象開展數(shù)值計算，通過SPH法模擬顆粒以不同角度多次撞擊套管材料，驗證套管材料42CrMo鋼的沖擊角函數(shù)特性曲線，將對應(yīng)的結(jié)果與CFD數(shù)值計算結(jié)果相結(jié)合。同時基于現(xiàn)場給定的5組工況進(jìn)行數(shù)值計算，獲得不同工況下的沖蝕結(jié)果，并進(jìn)行端口特性參數(shù)對沖蝕作用影響規(guī)律的研究。研究結(jié)果可為套管端口優(yōu)化設(shè)計提供參考。

1" 計算模型

1.1" SPH模型

SPH法是一種拉格朗日無網(wǎng)格粒子方法。它利用的核函數(shù)對物理問題進(jìn)行近似處理，是眾多無網(wǎng)格方法中應(yīng)用最廣泛的一種方法。該方法利用攜帶有一定物理量（溫度、壓力、黏度等）的一組粒子集合描述固體與流體等連續(xù)相。由于粒子的幾何構(gòu)型直接刻畫材料的運動和變形，故SPH法可以直觀描述材料界面并追蹤材料變形，這些要素在沖擊力學(xué)數(shù)值模擬中極為重要［19–21］。砂粒撞擊材料在微觀層面上屬于沖擊，適用于SPH法。

利用SPH方法對目標(biāo)材料進(jìn)行多砂粒沖擊模擬，獲得砂粒對目標(biāo)材料的沖擊特性，驗證沖蝕模型中的沖擊角函數(shù)。在SPH法中，控制方程的表達(dá)式如下：

ρit=∑Nj=1mjvijWijxi（1）

vit=-∑Nj=1mjσiρ2i+σjρ2jWijxi（2）

uit=12∑Nj=1mjσiρ2i+σjρ2jvijWijxi（3）

式中：ρ為粒子的密度，kg/m3；t為時間，s；mj為粒子j的質(zhì)量，kg；vij=vi-vj，為粒子i與j之間的相對速度矢量的分量，m/s；Wij為核函數(shù)，m-3；σi、σj分別為粒子i與j的應(yīng)力張量的分量，Pa；vi為粒子 i的速度，m/s；ui為顆粒的單位質(zhì)量內(nèi)能，m2/s2。

本文數(shù)值計算中所使用的三次樣條核函數(shù)是一種常見的核函數(shù)［22］，其表達(dá)式如下所示：

WR，h=αd23-R2+12R3" 0≤R＜1

162-R3" 1≤R＜2

0" R＞2（4）

R=ri-rj/h（5）

式中：R為點i與j之間的相對距離， m；h為光滑長度，m；αd為核函數(shù)歸一化系數(shù)，在不同維度空間中取值不同，在三維空間中取值為1π32h3;ri和rj分別是粒子i和粒子j的位置向量，m。

圖1為SPH核函數(shù)示意圖。

1.2" Johnson-Cook模型

Johnson-Cook 模型適用于大變形、高應(yīng)變率下的材料力學(xué)響應(yīng)，被廣泛應(yīng)用于粒子沖擊的數(shù)值模擬中［23］。砂粒沖擊壁面在微觀結(jié)構(gòu)上屬于大形變沖擊問題，故使用該模型對目標(biāo)材料的力學(xué)響應(yīng)進(jìn)行描述。Johnson-Cook 模型［24］包括本構(gòu)模型與失效模型，本構(gòu)模型用于描述套管材料的彈塑性力學(xué)行為，失效模型用于計算材料的失效行為，其形式如下：

σeq=A+Bεneq1+Clnε*1-T*m（6）

T*=T-Tr/Tm-Tr（7）

εf=D1+D2expD3σ*1+D4lnε*

1+D5T*（8）

式中：σeq為材料的流動應(yīng)力，Pa；A、B為屈服強度和應(yīng)變硬化系數(shù)，Pa；ε*為等效塑性應(yīng)變率與參考應(yīng)變率的比值，無量綱；T為材料的當(dāng)前溫度，K；Tr為材料的室溫，K；Tm為材料的熔點，K；n為材料應(yīng)變硬化指數(shù)； C為應(yīng)變率敏感系數(shù)；m為溫度軟化指數(shù)；εeq為等效塑性應(yīng)變；εf為失效應(yīng)變；D1、D2、D3為應(yīng)力三軸度影響系數(shù)；D4為應(yīng)變率敏感系數(shù)；D5為溫度敏感系數(shù)；σ*為等效應(yīng)力與平均應(yīng)力的比值，無量綱。

套管材料為42CrMo，數(shù)據(jù)來源見文獻(xiàn)［25］，取值如表 1與表 2所示。

1.3" 顆粒追蹤模型

含砂壓裂液在套管中的流動屬于固-液兩相流動，其中壓裂液為連續(xù)相，砂粒為離散相。在湍流的作用下，顆粒與流體之間相互作用，其控制方程如下：

dvpdt=FDv-vp+gρp-ρsρp+Fz（9）

FD=3μCDRe4ρpd2p（10）

Re=ρdpvp-vμ（11）

式中：vp為顆粒流速，m/s；v為流場速度，m/s；ρp為顆粒密度，kg/m3;ρs為流體密度，kg/m3；g為重力加速度，m/s2；Fz為單位質(zhì)量附加外力，N/kg;FD（v-vp）為單位質(zhì)量曳力，N/kg；Re為相對雷諾數(shù)，無量綱；μ為流體動力黏度，Pa·s；dp為顆粒直徑，m；Cd為曳力系數(shù)，無量綱。

1.4" 沖蝕模型

沖蝕作用過程復(fù)雜，受到多種因素影響，在數(shù)值模擬中一般利用沖蝕公式進(jìn)行預(yù)測。選用的Oka沖蝕模型考慮了粒子沖擊速度、粒子粒徑、沖擊角函數(shù)與目標(biāo)材料物性參數(shù)。Oka沖蝕模型方程如下［26-28］：

ER=ρtkfαHvk1vpv′k2dpd′k3（12）

fα=sin αn11+Hv1-sin αn2（13）

式中：ER為沖蝕率，kg/（m2·s）；fα為沖擊角函數(shù)；k、k1、k2、k3、n1、n2為給定常數(shù)；Hv為材料維氏硬度，HV；v′為參考速度，vp為粒子實際速度，m/s；d′為參考粒徑，m。

沖擊角函數(shù)表征撞擊角度與沖蝕作用之間的關(guān)系。

2" SPH法模擬砂粒沖擊套管材料

構(gòu)建的三維微觀撞擊模型如圖 2所示。計算多個球形砂粒對同一位置的撞擊結(jié)果，砂粒之間互不影響；模擬中設(shè)置8個砂粒以不同角度對同一位置的多次撞擊，具體如圖 3所示。在t=0時刻對目標(biāo)材料進(jìn)行SPH粒子轉(zhuǎn)換。

撞擊角度與速度分布如表 3所示。沖擊速度取Oka沖蝕公式推導(dǎo)所用試驗數(shù)據(jù)，沖擊角度為速度方向與水平方向夾角。

對于高角度沖擊，垂直速度分量較大，材料受到擠壓發(fā)生變形形成兩邊微微凸起，中部下陷的沖擊坑，凸起在被其他粒子沖擊后，容易形成碎屑而脫落。對于低角度沖擊，水平方向速度分量更大，同時又具有一定程度的垂直速度分量，形成了剪切樣貌的沖擊坑，突起更高，會更容易引發(fā)沖蝕效果。在20°～45°的入射角度范圍內(nèi)，所形成的突起更高，在實際沖蝕情況下更容易被后續(xù)粒子沖擊發(fā)生剪切作用，導(dǎo)致材料損失。隨著沖擊角度逐漸降低，目標(biāo)材料中的最大應(yīng)力值也在增加。由于低角度沖擊材料后表面硬化，且砂粒與目標(biāo)材料之間發(fā)生滑移，砂粒在撞擊過程中僅部分動能轉(zhuǎn)化為目標(biāo)材料內(nèi)能，故圖4中低沖擊角度影響范圍較小，而高角度沖擊下，由于速度垂直分量較大，更多動能轉(zhuǎn)化為目標(biāo)材料內(nèi)能，影響范圍較大。

在沖擊過程中，材料發(fā)生塑性變形而導(dǎo)致材料失效，以塑性應(yīng)變作為失效參考依據(jù)，同時依據(jù)沖擊角函數(shù)處理方法處理模擬數(shù)據(jù)。無量綱化沖蝕量E依據(jù)90°下沖擊角對應(yīng)沖蝕量進(jìn)行標(biāo)準(zhǔn)化處理，計算公式如下：

E=EαE90（14）

式中：E為無量綱化沖蝕損失質(zhì)量；Eα為在沖擊角度為α?xí)r的沖蝕損失質(zhì)量，mm2/a；E90為沖擊角度在90°時的沖蝕損失質(zhì)量，mm2/a。

無量綱沖蝕速率與沖擊角度的關(guān)系如圖 5所示。由圖5可知，沖蝕作用隨著沖擊角度的增大而加強，達(dá)到最大沖蝕作用后則會出現(xiàn)減弱趨勢，該趨勢與Oka沖蝕模型中沖擊角函數(shù)趨勢一致。

erosion rate and impact angle

同時，模擬結(jié)果與Oka沖蝕模型曲線中均存在沖擊角度為20°～45°時沖蝕作用較為強烈，隨后逐漸減弱。雖然高角度沖擊比低角度下沖擊的影響范圍大，但低角度下的沖蝕量比高角度下沖擊大。在沖擊角函數(shù)中，也將目標(biāo)材料硬度作為一個影響因素。對于大多數(shù)金屬來說，重復(fù)變形部分為定值，故n1部分取原Oka模型中默認(rèn)值（n1=0.8）。沖擊角函數(shù)分為重復(fù)變形與切削，如圖 6所示。

3" 套管滑套流場與內(nèi)壁沖蝕計算

3.1" 幾何模型

套管滑套結(jié)構(gòu)如同7所示。根據(jù)現(xiàn)場套管滑套模型，套管最大外徑為177 mm，在端口處內(nèi)徑為154 mm，滑套內(nèi)徑為124 mm，總長度為1 136 mm，端口形狀為膠囊形，長為74 mm，寬為20 mm。提取的流體域利用六面體網(wǎng)格與多面體網(wǎng)格共節(jié)點，進(jìn)一步提升網(wǎng)格整體質(zhì)量。圖8為套管及滑套流場域計算網(wǎng)格。

3.2" 沖蝕計算

對端口進(jìn)行5組典型工況計算，如表4所示。計算中所采用壓裂液為EM30S，其密度為1 100 kg/m3，黏度為99 mPa瘙簚s，所用砂粒是粒徑為380、830 μm的石英砂，其體積密度為1 470 kg/m3，視密

度為2 620 kg/m3，含砂體積分?jǐn)?shù)均設(shè)置為參考值10%。以上各個參數(shù)均于離散相選項中進(jìn)行設(shè)置。

5組工況下套管沖蝕速率分布云圖如圖9所示。由圖9可知，在不同工況下套管損傷嚴(yán)重的區(qū)域幾乎一致，在端口附近其沖蝕作用更為明顯。從整體上看，沖蝕作用主要集中在套管的前半部分。由于流體幾乎從端口處流出，在后半段中因封閉導(dǎo)致流體幾乎無法流動，砂粒運動速度也較小，故其沖蝕作用不明顯。

5組工況下近端口附近的沖蝕情況如圖 10所示。由圖10可知：在低排量時，砂粒的速度相對較小，其對套管內(nèi)壁的沖蝕作用不強，導(dǎo)致靠近端口處沖蝕速率較小；隨著排量的增大，砂粒在套管中的速度也隨之增大，在第3、4、5組工況下最大沖蝕速率和平均沖蝕速率增長速度比第1、2組工況下更高，沖蝕作用也更為明顯。

port under 5 working conditions

3.3" 端口特性沖蝕影響因素

端口特性會影響套管端口出處壁面的沖蝕作用。這里選用端口數(shù)量、端口面積與端口形狀共3種因素進(jìn)行分析。由于篇幅所限，以下分析均基于第3組工況。

3.3.1" 端口數(shù)量

端口作為套管中的流體出口，其數(shù)量影響套管端口附近的流場?？紤]到現(xiàn)場生產(chǎn)，2端口套管通流面積過于狹小，不適合大排量的固井滑套壓裂；10端口對套管結(jié)構(gòu)破壞較大，容易引發(fā)套管破壞等安全問題。故在每個端口面積相同的情況下，設(shè)計了4端口、6端口、8端口模型，如圖11所示。

為了更好地分析端口數(shù)量對沖蝕的影響，將端口及其附近區(qū)域單獨截取，對該部分沖蝕速率做平均計算處理。端口數(shù)量對平均沖蝕速率的影響如圖 12所示。由圖12可知：隨著端口數(shù)量的增加，端口附近的平均沖蝕速率呈現(xiàn)下降的趨勢，平均沖蝕速率：由4端口的9.65×10-4 kg/（m2·s）下降至8端口的2.30×10-4 kg/（m2·s），下降幅度為73.8%；沖蝕速率隨著端口數(shù)目增多而呈現(xiàn)出較為明顯的下降，但其下降的趨勢逐漸減弱。雖然無法無限制增加端口數(shù)目，在允許范圍內(nèi)增加端口數(shù)目可以減輕端口沖蝕作用。

3.3.2" 端口面積

端口面積也是影響沖蝕的一個主要因素。為了探究端口面積對沖蝕速率的影響，保持其端口形狀為膠囊形，分別對4種面積的膠囊形端口進(jìn)行計算。端口尺寸如圖 13所示。

端口面積對平均沖蝕速率的影響如圖 14所示。由圖14可知，隨著端口面積的增大，平均沖蝕速率逐漸減小，面積為1 000 mm2端口的平均

沖蝕速率為3.64×10-4 kg/（m2·s），而面積為1 600 mm2端口的平均沖蝕速率為2.1×10-4 kg/（m2·s）。平均沖蝕速率隨端口面積增大而減小，故在允許范圍內(nèi)，應(yīng)當(dāng)盡量增大固井滑套端口處的面積。

3.3.3" 端口形狀

在研究中發(fā)現(xiàn)，端口形狀在一定程度上影響著端口及其附近的沖蝕情況。分別對矩形、菱形、膠囊形、橢圓形4種類型的端口（見圖15）進(jìn)行計算，結(jié)果如圖16所示。

如圖 16所示，在4種相同面積的端口中，平均沖蝕速率最高的端口形狀為矩形端口，橢圓形與菱形端口的平均沖蝕速率相近，這表明它們更加流線型的邊緣可能有助于減輕流體引起的沖蝕作用。膠囊形端口對應(yīng)的平均沖蝕速率最低，其流線型設(shè)計優(yōu)化了流體流過端口的路徑，最小化了流體對端口壁面的直接撞擊和湍流強度，在4種端口類型中顯示出最佳的抗沖蝕性能，相比于矩形端口其抗沖蝕性能提高了35.2%。說明端口形狀對端口周圍流場有著較大的影響，不同類型的端口形狀有著不同的抗沖蝕特性，在無其他特殊要求時，應(yīng)選擇膠囊形端口。

4" 結(jié)" 論

（1）通過建立連續(xù)砂粒沖蝕磨損仿真模型，利用光滑粒子法分析了砂粒連續(xù)沖擊材料時所發(fā)生的特征現(xiàn)象，結(jié)果表明，滑套材料42CrMo鋼在沖擊角度為20°～45°時有更大的沖蝕損失，其沖蝕規(guī)律與Oka沖蝕模型相符，也符合42CrMo鋼的金屬特性。

（2）通過CFD數(shù)值模擬5組不同工況下的套管沖蝕過程，結(jié)果表明，在整個固井滑套內(nèi)壁中，沖蝕最嚴(yán)重的區(qū)域出現(xiàn)在端口及其附近，主要原因為出口處出流速增加導(dǎo)致砂粒撞擊的速度增加，進(jìn)而使沖蝕速率明顯增大。

（3）端口的特性參數(shù)對套管及其附近的沖蝕有著明顯的影響，其影響因素包括端口數(shù)目、端口面積以及端口形狀。隨著端口數(shù)目的增加，平均沖蝕速率會在一定程度內(nèi)降低，但是降低的趨勢明顯減緩；隨著端口面積增加，平均沖蝕速率會降低，但是降低的趨勢也隨之減緩；端口形狀對沖蝕速率的影響也較大，膠囊形端口的抗沖蝕效果最好。合理選擇端口的特性參數(shù)能在一定程度上減輕固井滑套的沖蝕磨損。

［1］" "關(guān)皓綸，王兆會，劉斌輝．分段壓裂固井滑套的研制現(xiàn)狀及展望［J］．石油機械，2021，49（11）：84-92．

GUAN H G， WANG Z H， LIU B H. The technical states of cementing sliding sleeve for staged fracturing［J］. China Petroleum Machinery， 2021， 49（11）： 84-92.

［2］" 葉昊，陳林婭，周鴻翔，等．油氣管道彎管肘部二次流動侵蝕數(shù)值模擬［J］．石油機械，2023，51（4）：119-126．

YE H， CHEN L Y， ZHOU H X， et al. Numerical simulation on secondary flow erosion at elbow of oil and gas pipeline［J］. China Petroleum Machinery， 2023， 51（4）： 119-126.

［3］" 涂亞東，徐小兵，殷華鋒，等．高壓管匯沖蝕磨損規(guī)律研究［J］．石油機械，2018，46（2）：84-88．

TU Y D， XU X B， YIN H F， et al. Analysis on erosion wear laws of high pressure manifold［J］. China Petroleum Machinery， 2018， 46（2）： 84-88.

［4］" 張繼信，康健，樊建春，等．高壓彎管沖蝕失效分析及數(shù)值模擬［J］．潤滑與密封，2016，41（10）：40-46．

ZHANG J X， KANG J， FAN J C， et al. Failure analysis and numerical simulation of high pressure elbow erosion wear［J］. Lubrication Engineering， 2016， 41（10）： 40-46.

［5］" 莫麗，劉棲，許先雨．基于CFD數(shù)值模擬的異面三通管沖蝕磨損規(guī)律研究［J］．潤滑與密封，2022，47（8）：41-46．

MO L， LIU Q， XU X Y. Study on erosion wear law of tee in different planes based on CFD numerical simulation［J］. Lubrication Engineering， 2022， 47（8）： 41-46.

［6］" 王川，謝真強，王國榮，等．基于CFD的固井滑套沖蝕特性分析［J］．機械設(shè)計，2016，33（7）：66-70．

WANG C， XIE Z Q， WANG G R， et al. Erosion feature analysis of cementing sliding sleeve based on CFD simulation［J］. Journal of Machine Design， 2016， 33（7）： 66-70.

［7］" 魏遼，韓峰，陳濤，等．套管固井滑套沖蝕磨損模擬分析與試驗研究［J］．石油鉆探技術(shù)，2014，42（3）：108-111．

WEI L， HAN F， CHEN T， et al. Analysis and experimental research on erosion of cementing sliding sleeve［J］. Petroleum Drilling Techniques， 2014， 42（3）： 108-111.

［8］" WONG C Y， SOLNORDAL C， SWALLOW A， et al. Predicting the material loss around a hole due to sand erosion［J］. Wear， 2012， 276/277： 1-15.

［9］" XU Y Q， GUAN Z C， LIU Y W， et al. Structural optimization of downhole float valve via computational fluid dynamics［J］. Engineering Failure Analysis， 2014， 44： 85-94.

［10］" 向正新，李思行，錢利勤，等．壓裂球座沖蝕磨損規(guī)律研究和結(jié)構(gòu)優(yōu)化［J］．工程設(shè)計學(xué)報，2017，24（5）：555-562．

XIANG Z X， LI S H， QIAN L Q， et al. Study on erosion law and structure optimization of fracturing ball seat［J］. Chinese Journal of Engineering Design， 2017， 24（5）： 555-562.

［11］" PARSI M， NAJMI K， NAJAFIFARD F， et al. A comprehensive review of solid particle erosion modeling for oil and gas wells and pipelines applications［J］. Journal of Natural Gas Science and Engineering， 2014， 21： 850-873.

［12］" DONG X W， LIU G R， LI Z L， et al. A smoothed particle hydrodynamics （SPH） model for simulating surface erosion by impacts of foreign particles［J］. Tribology International， 2016， 95： 267-278.

［13］" DONG X W， LI Z L， FENG L， et al. Modeling， simulation， and analysis of the impact（s） of single angular-type particles on ductile surfaces using smoothed particle hydrodynamics［J］. Powder Technology， 2017， 318： 363-382.

［14］" TAKAFFOLI M， PAPINI M. Material deformation and removal due to single particle impacts on ductile materials using smoothed particle hydrodynamics［J］. Wear， 2012， 274/275： 50-59.

［15］" AZIMIAN M， SCHMITT P， BART H J. Numerical investigation of single and multi impacts of angular particles on ductile surfaces［J］. Wear， 2015， 342/343： 252-261.

［16］" HADAVI V， PAPINI M. Numerical modeling of particle embedment during solid particle erosion of ductile materials［J］. Wear， 2015， 342/343： 310-321.

［17］" 廖翔云，馬小晶，周新超，等．基于SPH方法的組合射流破巖模擬研究［J］．石油機械，2022，50（5）：1-8．

LIAO X Y， MA X J， ZHOU X C， et al. Numerical simulation of rock breaking by combined Jets based on SPH method［J］. China Petroleum Machinery， 2022， 50（5）： 1-8.

［18］" 王方祥，楊仲才，趙增權(quán)，等．壓裂液返排期間放噴油嘴的磨損規(guī)律研究［J］．石油機械，2020，48（1）：102-109．

WANG F X， YANG Z C， ZHAO Z Q， et al. Study on erosion wear of discharge nozzle during fracturing fluid backflow［J］. China Petroleum Machinery， 2020， 48（1）： 102-109.

［19］" 杜明超，李增亮，董祥偉，等．菱形顆粒沖擊材料表面沖蝕磨損特性分析［J］．摩擦學(xué)學(xué)報，2020，40（1）：1-11．

DU M C， LI Z L， DONG X W， et al. Analysis of material surface erosion characteristics due to rhomboid-shaped particle impact［J］. Tribology， 2020， 40（1）： 1-11.

［20］" 靳紹虎，劉科偉，黃進(jìn)，等．基于Lagrange及SPH算法的花崗巖侵徹數(shù)值模擬［J］．高壓物理學(xué)報，2021，35（5）：143-155．

JIN S H， LIU K W， HUANG J， et al. Numerical simulation of granite penetration based on Lagrange and SPH algorithm［J］. Chinese Journal of High Pressure Physics， 2021， 35（5）： 143-155.

［21］" HE Q G， CHEN X W， CHEN J F. Finite element-smoothed particle hydrodynamics adaptive method in simulating debris cloud［J］. Acta Astronautica， 2020， 175： 99-117.

［22］" 劉維平，張永祥．SPH方法在潰壩數(shù)值模擬中的應(yīng)用［J］．重慶大學(xué)學(xué)報，2017，40（11）：11-19．

LIU W P， ZHANG Y X. Application of SPH method to numerical simulation of dam-break［J］. Journal of Chongqing University（Natural Science Edition）， 2017， 40（11）： 11-19.

［23］" 舒暢，程禮，許煜．Johnson-Cook本構(gòu)模型參數(shù)估計研究［J］．中國有色金屬學(xué)報，2020，30（5）：1073-1083．

SHU C， CHENG L， XU Y. Research on parameter estimation of Johnson-Cook constitutive model［J］. The Chinese Journal of Nonferrous Metals， 2020， 30（5）： 1073-1083.

［24］" JOHNSON G R， COOK W H. Fracture characteristics of three metals subjected to various strains， strain rates， temperatures and pressures［J］. Engineering Fracture Mechanics， 1985， 21（1）： 31-48.

［25］" NING J Q， LIANG S Y. Model-driven determination of Johnson-Cook material constants using temperature and force measurements［J］. The International Journal of Advanced Manufacturing Technology， 2018， 97（1）： 1053-1060.

［26］" OKA Y I， OKAMURA K， YOSHIDA T. Practical estimation of erosion damage caused by solid particle impact： Part 1： Effects of impact parameters on a predictive equation［J］. Wear， 2005， 259（1/2/3/4/5/6）： 95-101.

［27］" OKA Y I， YOSHIDA T. Practical estimation of erosion damage caused by solid particle impact： part 2： mechanical properties of materials directly associated with erosion damage［J］. Wear， 2005， 259（1/2/3/4/5/6）： 102-109.

［28］" 卓仁燕， 馬新仿， 李建民， 等. 水平井限流壓裂對射孔孔眼沖蝕的影響［J］. 鉆采工藝， 2023， 46（2）： 77-82.

ZHUO R Y， MA X F， LI J M， et al. Effect of Limited Entry Fracturing for Horizontal Wells on Perforation Erosion［J］. Drilling amp; Production Technology， 2023， 46（2）： 77-82.

第一張磊，副教授，生于1983年，博士，研究方向為多相流數(shù)值模擬。地址：（102249）北京市昌平區(qū)。電話：（010）89733537。email： zlei@cup.edu.cn。

2024-01-28" 修改稿收到日期：2024-07-15

劉鋒