成 芳，易先中，彭 灼，周元華，宋順平，盛治新，殷光品

(1.長江大學(xué) 機械工程學(xué)院，湖北 荊州 434000；2.中國石油集團川慶鉆探工程有限公司長慶鉆井總公司，陜西 西安 710021；3.青島金江源工業(yè)裝備有限公司，山東 青島 266111；4.湖北佳業(yè)石油機械股份有限公司，湖北 荊州434000)

0 引言

管道在城市給排水、儲運工程、化學(xué)工程及石油工業(yè)等各方面都發(fā)揮著重要的紐帶作用。頁巖氣在儲層上具有儲層致密、滲透率低等特征，此類特征會造成頁巖氣開發(fā)困難，水力壓裂技術(shù)是目前頁巖氣開發(fā)依靠的主要儲層改造手段[1-2]。高壓壓裂液和壓力管線是壓裂過程工藝中的重要組成部分，壓裂液中攜有的支撐劑固體顆粒會隨著高壓高速流體流動對管線造成不同程度的沖蝕磨損，一般認(rèn)為，沖蝕是流體介質(zhì)攜帶固體顆粒并以一定的速度或角度沖擊材料壁面所造成的磨損現(xiàn)象[3-5]。值得注意的是，沖蝕作為工業(yè)生產(chǎn)過程中引起材料磨損與設(shè)備失效的重要失效行為，具有突發(fā)性和危險性，磨損機理復(fù)雜，是一種極難預(yù)防的破壞形式[6]。有大量研究表明，彎管段的沖蝕磨損比直管段大約嚴(yán)重50倍[7]。同時，沖蝕會改變材料的性能，使得設(shè)備的剩余強度發(fā)生變化[8-9]，嚴(yán)重時將導(dǎo)致管匯的刺穿、破裂等重要經(jīng)濟損失和設(shè)備故障[10]。

由于計算流體動力學(xué)理論、CFD技術(shù)和實驗技術(shù)手段的不斷更新與發(fā)展，國內(nèi)外均對此有了大量的研究并取得了相應(yīng)成果。易先中等[5]對同一影響因素取平均值并比較其沖蝕速率的增長倍數(shù)，認(rèn)為流速是影響沖蝕磨損的主要因素；楊向前等[11]通過沖蝕磨損試驗機分析得到了不同直徑固體顆粒對35CrMo鋼沖蝕的影響，探討了在一定范圍內(nèi)的沖蝕角度與沖蝕速度時沖蝕的影響；孟杏等[12]對試樣進行了沖擊磨損試驗，分析了沖擊角度、沖蝕介質(zhì)和沖擊時間對沖蝕磨損量的影響；Lopez等[13]研究了顆粒速度和平均沖擊角對溶解的AISI 304和淬火回火AISI 420不銹鋼沖蝕的影響；Sugiyama等[14]對不同材料在鉆井液90°沖擊角下沖蝕后的表面輪廓進行了數(shù)值分析。總之，目前在水力壓裂對彎管的沖蝕研究方面，主要集中在流體的特性參數(shù)對管壁及不同材料的沖蝕影響分析，但在現(xiàn)場不可避免使用不同管徑的彎管，壓裂液的流速也會依據(jù)工程需要進行改變。而對于不同流速與不同管徑耦合時，沖蝕位置的變化及沖蝕規(guī)律尚不夠清楚，因此本文將對彎管進行數(shù)值模擬分析，探討沖蝕位置隨流速與管徑變化時的變化規(guī)律。

1 模型建立

1.1 幾何建模及參數(shù)

本文的研究對象為JY型彎管，其基本結(jié)構(gòu)尺寸是內(nèi)徑為D，曲率半徑為A=2D的彎曲度90°圓形彎管，材質(zhì)為PCrNi3Mo，建立模型進行沖蝕磨損分析。由于入口端面效應(yīng)的影響，為了確保流體在管道內(nèi)充分流動和流場穩(wěn)定，在彎管入口段和出口段延長距離為L=5D的直管段，如圖1所示為計算對象。網(wǎng)格采用六面體結(jié)構(gòu)性網(wǎng)格劃分方式，在彎管段進行網(wǎng)格加密處理，為保證更為真實準(zhǔn)確的流體流動狀態(tài)，需在管壁處設(shè)置合理的邊界層。壓裂液物理性能的基本參數(shù)如表1所示。

圖1 彎管幾何模型Fig.1 Geometric model of bend pipe

表1 壓裂液物理性能基本參數(shù)Table 1 Basic parameters for physical properties of fracturing fluid

1.2 計算模型的建立

水力壓裂過程中，壓裂液在快速流經(jīng)彎管時，管道內(nèi)部流體為復(fù)雜的液-固兩相三維紊流流場。因所計算的壓裂液支撐劑固體顆粒體積分?jǐn)?shù)小于10%，故對彎管的沖蝕數(shù)值模擬選用DPM模型。其中，連續(xù)相為壓裂液流體，離散相為壓裂液支撐劑顆粒，由于顆粒體積分?jǐn)?shù)小，可忽略顆粒間的相互作用，從而滿足顆粒負(fù)載流動模式。并在Eulerian坐標(biāo)系下求解連續(xù)相流體的流場，在Lagrangian坐標(biāo)系下計算離散相固體顆粒的運動軌跡方程。管內(nèi)部三維紊流流場的數(shù)值模擬采用標(biāo)準(zhǔn)k-ε方程模型。為了進行彎管沖蝕特性的研究，本文對分析做出如下假設(shè)：1)認(rèn)為管內(nèi)壁完整，不存在缺陷與凹坑；2)忽略溫度變化對流體黏度的影響，不考慮重力因素對流體流速的影響；3)流體不可壓縮并具有連續(xù)性，為穩(wěn)定定常流動，不存在氣泡，流體充滿管道內(nèi)部。

壓裂彎管內(nèi)壁面的沖蝕速率的大小與顆粒質(zhì)量流量、顆粒沖擊速度、顆粒粒徑大小、顆粒特性等因素有關(guān)，根據(jù)文獻定義[15]，即：

(1)

f(θ)=2.69θ+1.61θ2-8.84θ3+7.33θ4-1.85θ5

(2)

式中：Rerosion為沖蝕速率，kg·m-2·s-1；N為顆粒碰撞數(shù)目，無因次；mp為顆粒平均質(zhì)量流量，kg·s-1；C(dp)為顆粒粒徑的函數(shù)，一般情況下，C(dp)=1.8×10-9；θ為顆粒對管道壁面的沖擊角(入射角)，(°)；f(θ)為沖擊角(入射角)θ的函數(shù)；v為顆粒沖擊速度，m·s-1；b(v)為顆粒相對于速度v的函數(shù)，決定沖蝕快慢的重要參數(shù)，默認(rèn)為2.6；Af為沖蝕壁面計算單元的面積，m2。

在運用Fluent計算過程中，進口邊界條件定義為速度入口(velocity inlet)，出口邊界定義為自由流出(outflow)，設(shè)置入口處壓裂液流體與支撐劑固體顆粒的速度大小相等、方向相同，DPM模型中入口與出口均采用逃逸(Escape)類型，壁面條件為Wall壁面，在DPM模型中采用反彈(reflect)類型，反彈系數(shù)依據(jù)A.Forder和KR.Ahlert等[16-17]的研究成果取值，動能、能量、湍流耗散率的離散均選用二階迎風(fēng)格式。

2 計算結(jié)果分析

網(wǎng)格作為有限元計算的重要前提，網(wǎng)格數(shù)量的變化會導(dǎo)致計算結(jié)果的差異化。對所計算對象的網(wǎng)格進行無關(guān)性驗證，文獻[18]指出國際學(xué)術(shù)界認(rèn)為網(wǎng)格變化導(dǎo)致計算數(shù)據(jù)變化小于2%為可取值，因此本文以彎管的最大沖蝕率的變化為界定標(biāo)準(zhǔn)，當(dāng)彎管整體計算模型網(wǎng)格數(shù)量達(dá)到101 190時，最大沖蝕率數(shù)值變化已小于2%，計算結(jié)果滿足網(wǎng)格無關(guān)性的基本要求。且網(wǎng)格平均質(zhì)量高于0.8，等角斜率小于0.76，網(wǎng)格質(zhì)量較好。

2.1 彎管的流場分析

為了分析流體攜帶固體顆粒在流經(jīng)彎管時，管內(nèi)復(fù)雜的流場變化情況。依據(jù)已確定的初始數(shù)據(jù)、邊界條件及計算方法，本文選取壓裂液流速為5～25 m/s，管徑為50.8～152.4 mm的彎管設(shè)置合理梯度進行數(shù)值模擬分析，得到其壓力場分布與最大沖蝕速率。下面將選取在入口速度在5～25 m/s，管徑為127 mm的彎管進行分析說明。

從圖2可以得知，在管徑一定時，不論流速如何變化，在彎管段存在較為明顯的壓力梯度，具體分布為從彎管段的外側(cè)區(qū)域至內(nèi)側(cè)區(qū)域，壓力數(shù)值呈現(xiàn)逐漸減小的趨勢。同時，流體流經(jīng)彎管時的最大壓力均分布在彎管段內(nèi)壁面外側(cè)區(qū)域，最小壓力分布在彎管段內(nèi)壁面內(nèi)側(cè)區(qū)域。

圖2 彎管壓力云圖Fig.2 Nephogram of bend pipe pressure

為進一步闡述高壓壓裂液在流經(jīng)管道時的流場特性變化情況，將從2條流動路線進行考慮分析。壓裂液從“入口直管段→彎管段外側(cè)區(qū)域→出口直管段”流動時，壓力數(shù)值具體變化為“壓力減小→變大→減小”；壓裂液從“入口直管段→彎管段內(nèi)側(cè)區(qū)域→出口直管段”流動時，壓力數(shù)值具體變化為“壓力減小→減小→變大”?？梢哉J(rèn)為，流體流經(jīng)彎管段時，由于離心力的作用，使得彎管段的流體由內(nèi)側(cè)區(qū)域朝著外側(cè)區(qū)域運動，導(dǎo)致產(chǎn)生彎管段外側(cè)壓力大，內(nèi)側(cè)壓力小的現(xiàn)象。即彎管段壓力小的區(qū)域流體速度大，反之，壓力大的區(qū)域流體速度小。流體進入入口直管段至彎管段處壓力逐漸減小的原因為，在沿著流體流動方向有較大的沿程阻力，且在彎管段流體內(nèi)側(cè)區(qū)域向外側(cè)運動時產(chǎn)生了渦流。在彎管出口至出口直管處，且由于在彎管段渦流的產(chǎn)生和作用，使得有局部的壓力損失，即出口直管段的壓力變小。

2.2 彎管的沖蝕磨損分析

數(shù)值計算結(jié)果表明，彎管中沖蝕磨損發(fā)生區(qū)域有5處，分別為彎管的入口直管段、彎管段內(nèi)壁面外側(cè)(液體進入彎管后的正對區(qū)域)、彎管段靠近出口內(nèi)壁面的外側(cè)區(qū)域、直管段靠近彎管段的側(cè)方區(qū)域以及出口直管段。但是，在彎管的入口直管段僅有著零星狀的沖蝕磨損現(xiàn)象發(fā)生，在出口直管段也存在著一定程度的均勻點狀沖蝕磨損情況，入口直管段與出口直管段的沖蝕磨損>情況均不嚴(yán)重且不會隨著流速等特征參數(shù)變化而產(chǎn)生明顯的磨損程度加深和區(qū)域變大。

總體分析認(rèn)為，壓裂液固體顆粒對彎管的易沖蝕區(qū)域且沖蝕嚴(yán)重區(qū)域主要有3處，如圖3所示為彎管段靠近出口內(nèi)壁面的外側(cè)區(qū)域(A區(qū)域)、彎管段內(nèi)壁面外側(cè)(液體進入彎管后的正對區(qū)域)(B區(qū)域)以及直管段靠近彎管段的側(cè)方區(qū)域(C區(qū)域)。如圖4所示是管徑為127 mm的彎管在入口速度為5，15，25 m/s時的沖蝕云圖。

圖3 彎管易沖蝕區(qū)域示意Fig.3 Schematic diagram of erosion prone area in bend pipe

值得注意的是，無論入口速度如何變化，彎管段靠近出口內(nèi)壁面的外側(cè)區(qū)域(A區(qū)域)沖蝕磨損始終較為嚴(yán)重，而彎管段內(nèi)壁面外側(cè)(液體進入彎管后的正對區(qū)域)(B區(qū)域)和直管段靠近彎管段的側(cè)方區(qū)域(C區(qū)域)的沖蝕磨損情況隨著入口速度變化而改變?？梢赃@樣認(rèn)為，以上3處存在較大沖蝕磨損的原因在于，高壓壓裂液自直管段進入到彎管段后，由于彎管段存在壓力梯度、離心力等作用形式，使得壓裂液及大量的壓裂液固體顆粒均被甩向曲率半徑較大的外壁面，即顆粒由內(nèi)側(cè)朝著外側(cè)運動，且沖擊速度明顯變大，造成顆粒沖擊個數(shù)、壁面沖擊次數(shù)明顯變多，即彎管段外壁面(A區(qū)域)的沖蝕磨損較為嚴(yán)重。且隨著流速的變大，固體顆粒撞擊到壁面的速度也相應(yīng)變大，即使得顆粒的動能變大，在單位時間內(nèi)撞擊到壁面的顆粒個數(shù)與沖擊動能也有了較大的提升，使得彎管段內(nèi)壁面外側(cè)(液體進入彎管后的正對區(qū)域)(B區(qū)域)的沖蝕磨損程度加大。而在直管段靠近彎管段的側(cè)方區(qū)域(C區(qū)域)出現(xiàn)沖蝕磨損的原因在于二次流的影響使得顆粒沖擊到管壁側(cè)壁的速度變大，且壓裂液固體顆粒在沖擊到外側(cè)區(qū)域的壁面后，顆粒由于反彈效應(yīng)撞擊至彎管的側(cè)方，使得在側(cè)方也出現(xiàn)了一定的沖蝕磨損現(xiàn)象。

特別的，如圖4所示壓裂液流速從5 m/s至25 m/s時，經(jīng)分析得壓裂液流速的增大會導(dǎo)致直管段靠近彎管段側(cè)方區(qū)域沖蝕磨損情況減弱，而與直管段正對方向的彎管段區(qū)域的沖蝕磨損范圍不斷變大，即并非所有彎管段的側(cè)方區(qū)域(C區(qū)域)和直管段正對方向的彎管段區(qū)域(B區(qū)域)都會有較為嚴(yán)重的沖蝕磨損情況。為分析這一現(xiàn)象，本文將引入斯托克斯數(shù)(St)進行分析，斯托克斯數(shù)(St)是能用于度量顆粒的動量非平衡程度的無量綱量，其定義顆粒響應(yīng)時間與系統(tǒng)響應(yīng)時間之比為：

圖4 彎管沖蝕云圖Fig.4 Nephogram of bend pipe erosion

(3)

式中：ρd為顆粒密度，kg·m-3；dp為顆粒粒徑，mm；v為流速，m·s-1；μ為流體黏度，Pa·s；D為管徑，mm。

總體而言，流體攜帶固體顆粒從直管進入彎管區(qū)域后，固體顆粒會受到慣性力、曳力、二次流等作用力的共同作用，使得流體及固體顆粒的受力形式和運動軌跡更為復(fù)雜。所謂慣性力，是指在彎管中保持固體顆粒沿切向運動的作用力，即慣性力的作用會使得彎管在液固兩相流的沖蝕磨損發(fā)生在彎管內(nèi)壁面的外側(cè)區(qū)域，且為液體進入彎管后的正對區(qū)域；而曳力是指使得顆粒沿著流體流線運動的作用力，曳力的作用下將使得固體顆粒沿著液體流線流出彎管段；二次流是指流體或顆粒偏離主流運動方向上的平均流動現(xiàn)象，二次流的出現(xiàn)更增加了彎管中流體流動和沖蝕磨損的復(fù)雜性，在該作用力下顆粒將朝著彎管內(nèi)壁面的側(cè)方區(qū)域運動，使得在該處出現(xiàn)沖蝕磨損現(xiàn)象。

由式(3)斯托克斯數(shù)(St)的定義，斯托克斯數(shù)(St)也表征著顆粒的慣性力與曳力之比，分析將綜合考慮2種作用力的比值，當(dāng)St<1時，即慣性力的數(shù)值小于曳力時，顆粒響應(yīng)時間較短，顆粒的運動將緊跟連續(xù)相，顆粒的跟隨性較好。同樣的，當(dāng)顆粒在彎管段中慣性力的數(shù)值小于曳力時，此時的曳力將起主導(dǎo)作用，使得慣性力與曳力的方向發(fā)生偏移，合力方向指向彎管段靠近出口直管段內(nèi)壁面的外側(cè)區(qū)域(A區(qū)域)，大量固體顆粒向著該合力方向所指處沖擊，即沖蝕磨損區(qū)域主要發(fā)生在彎管段靠近出口內(nèi)壁面的外側(cè)區(qū)域；當(dāng)St>1時，即慣性力的數(shù)值大于曳力時，表示顆粒無足夠響應(yīng)時間來響應(yīng)流速的變化，顆粒的運動將獨立于連續(xù)相，即顆粒的運動受流體影響小，此時顆粒在彎管段中慣性力的數(shù)值大于曳力時，此時的慣性力將起主要作用，在慣性力與曳力的共同作用下，彎管段內(nèi)壁面外側(cè)(液體進入彎管后的正對區(qū)域)沖蝕磨損區(qū)域(B區(qū)域)的范圍將變大，即在該區(qū)域會有連續(xù)均勻的片狀沖蝕磨損。

而對于二次流對彎管側(cè)方區(qū)域的沖蝕影響，本文認(rèn)為應(yīng)以St的具體數(shù)值進行評價，在圖4中，管徑為127 mm的彎管在入口速度為5，10，15，20，25 m/s時，經(jīng)計算其St分別為0.334，0.688 9，1.032，1.376，1.72?？梢哉J(rèn)為，當(dāng)St不斷增長時，彎管側(cè)方區(qū)域的最大沖蝕速率也不斷下降，故認(rèn)為，由于St的增長，流體攜帶固體支撐劑顆粒的能力得到遏制，無法攜帶大量的固體支撐劑顆粒向著彎管內(nèi)側(cè)壁面運動，即沖擊彎管側(cè)方區(qū)域壁面的固體顆粒數(shù)量有限，沖蝕磨損程度也呈現(xiàn)下降趨勢。

同時，當(dāng)St=1左右時，側(cè)方區(qū)域的沖蝕磨損下降至幾乎與出口直管段的沖蝕磨損程度相近，但無論St數(shù)值怎么上升，依據(jù)圖4的彎管沖蝕云圖顯示，側(cè)方區(qū)域的沖蝕磨損情況已趨于穩(wěn)定。

值得注意的是，彎管段內(nèi)壁面外側(cè)(液體進入彎管后的正對區(qū)域)(B區(qū)域)與直管段靠近彎管段的側(cè)方區(qū)域(C區(qū)域)的最大沖蝕速率及磨損區(qū)域會隨著流動參數(shù)及管徑等結(jié)構(gòu)參數(shù)的變化而產(chǎn)生明顯的規(guī)律性差異，2處的沖蝕磨損情況呈現(xiàn)“此消彼長”的情況。

2.3 管徑與流速耦合沖蝕分析

由于在不同管徑與不同流速等特征參數(shù)下，與之對應(yīng)的最大沖蝕速率不斷變化，為綜合考慮管徑、流速與最大沖蝕速率的變化，得到如圖5所示的三維曲面關(guān)系示意圖。

圖5 不同管徑與不同流速下的彎管最大沖蝕速率關(guān)系Fig.5 Relationship between maximum erosion rate of bend pipe and different pipe diameters and flow velocities

由圖5可知，在任意管徑下，壓裂液流速的增長必然使得彎管的最大沖蝕速率呈指數(shù)關(guān)系增長。隨著流速的增長，依據(jù)式(3)斯托克斯數(shù)(St)將不斷變大，顆粒對彎管的最大沖蝕磨損區(qū)域?qū)⒂蒀區(qū)域向著A區(qū)域變化。

為比較分析不同管徑的彎管在不同壓裂液流速下最大沖蝕速率，當(dāng)管徑從50.8 mm至152.4 mm變化時，經(jīng)軟件Origin進行公式擬合得其指數(shù)函數(shù)關(guān)系分別為：

(4)

(5)

(6)

(7)

(8)

式中：y表示最大沖蝕速率，kg·m-2·s-1；x表示壓裂液流速，m/s；y1～y5分別是管徑為50.8，76.2，101.6，127，152.4 mm時與不同流速的函數(shù)關(guān)系，其擬合相關(guān)系數(shù)分別為R2=0.976 12，0.996 33，0.997 09，0.997 24，0.999 1，系數(shù)均接近于1，表明擬合關(guān)系較為合理，效果較好。

而隨著管徑從50.8 mm至152.4 mm變化時，在同一流速下，其斯托克斯數(shù)(St)也將不斷的減小，使得顆粒對彎管的最大沖蝕磨損區(qū)域發(fā)生在C區(qū)域與A區(qū)域。對于任一壓裂液流速，隨著管徑的增大，其最大沖蝕速率均呈現(xiàn)減小的趨勢。在管徑從50.8 mm至152.4 mm變化時，計算得知壓裂液流速為5，10，15，20，25 m/s下的彎管最大沖蝕速率分別減小了85.78%，87.53%，89.78%，91.53%，94.01%。在管徑為50.8 mm，壓裂液流速為25 m/s時，最大沖蝕速率達(dá)到了最大值，為0.001 35 kg·m-2·s-1?？梢哉J(rèn)為隨著彎管管徑的不斷變大，管徑的內(nèi)表面積變大，同等流速的壓裂液中固體支撐劑顆粒與管壁發(fā)生沖擊、碰撞的面積變大，則在單位面積上發(fā)生顆粒沖擊的次數(shù)減少，使得沖蝕磨損速率減小。即得知增大管徑，是減小沖蝕磨損的有效途徑。

3 結(jié)論

1)彎管的沖蝕磨損發(fā)生區(qū)域共有5處，易沖蝕且嚴(yán)重區(qū)域主要有3處，為彎管段靠近出口內(nèi)壁面的外側(cè)區(qū)域(A區(qū)域)、彎管段內(nèi)壁面外側(cè)(液體進入彎管后的正對區(qū)域)(B區(qū)域)以及直管段靠近彎管段的側(cè)方區(qū)域(C區(qū)域)。彎管中流場的作用，使得彎管段的流體攜帶固體顆粒由內(nèi)側(cè)區(qū)域朝著外側(cè)區(qū)域運動，而渦流的作用使得部分固體顆粒朝著彎管側(cè)方運動。

2)St<1時，沖蝕磨損主要發(fā)生在彎管段靠近出口內(nèi)壁面的外側(cè)區(qū)域(A區(qū)域)；St>1時，彎管段內(nèi)壁面外側(cè)(液體進入彎管后的正對區(qū)域)(B區(qū)域)的沖蝕磨損范圍將變大。St增大使得彎管段側(cè)方區(qū)域(C區(qū)域)的沖蝕磨損程度呈現(xiàn)下降趨勢，St=1左右時，側(cè)方區(qū)域的沖蝕磨損下降至幾乎與出口直管段的沖蝕磨損程度相近，St繼續(xù)上升，側(cè)方區(qū)域的沖蝕磨損情況已趨于穩(wěn)定。隨著St在0～1至St>1變化時，B區(qū)域與C區(qū)域的沖蝕磨損情況呈現(xiàn)出“此消彼長”的情況。

3)管徑越小，其最大沖蝕速率的增長幅度越明顯，管徑增大，最大沖蝕速率呈減小的趨勢，即增大管徑是減小沖蝕磨損的有效途徑之一。