中圖分類號：TE937文獻(xiàn)標(biāo)識碼：A DOI：10.12473/CPM.202406039

Zhong Yin，Chen Shu，Zhao Yuhong，et al.Numerical simulation studyon the descaling mechanismof high-pres-sureswirling jet[J].ChinaPetroleumMachinery，2025，53（7）：112-121，148.

Numerical SimulationStudy on theDescalingMechanism of High-Pressure Swirling Jet

Zhong Yin1 Chen Shu2Zhao Yuhong'Xue Yongzhi3Yao Lin3 （1.SantanghuOilProductionManagementAreaofPetroChinaTuhaOilfieldCompany；2.SharedSeruiceCenterofPetroChina TuhaOilfield Company；3.SchoolofMechanical Engineering，Southwest Petroleum University）

Abstract： To further explore the active mechanism of swirling jet，a numerical model of swirling jet descaling was built，and its effctivenesswasverified bycomparing with experimental results.Then，the difference in flow field characteristics between swirling jetand conventional straight jetwas compared，and the mechanismofoverflow and erosionin the process ofswirling jet descaling was revealed.Finally，the influenceof pump pressure，target distance and jetangleon the descaling performance was analyzed.The study results show that the special nozle structuredetermines thethre-dimensional velocitydistribution characteristicsoftheswirling jet.Theswirling jethas a similarimpact crushing capacityand alarger impact area as theconventional straight jet.An increase in pump pressure can significantly improve the descaling performance of the swirling jet，the target distance has anoptimal value for affecting the descaling performance，and moderately increasing the jet angle can further promote the descaling effect. The study results provide reference forthe optimal application ofswirling jet indescaling operations.

Keywords： swirling jet descaling；nozzle structure； flow field characteristics；working parameters；evaluation index

0 引言

在油氣開采作業(yè)中，油管的結(jié)垢現(xiàn)象已成為制約油氣資源高效開發(fā)利用的重要因素。以吐哈油田為例，其油田區(qū)塊具有多斷塊、低滲透和弱揮發(fā)性油藏特征，油藏埋藏深且地層水富含成垢離子，為結(jié)垢提供了化學(xué)條件。油管結(jié)垢會導(dǎo)致其內(nèi)徑縮小，管道阻力增加，流體輸運(yùn)效率降低，進(jìn)而影響油氣井的正常作業(yè)[1-2]。因此，研究油氣井管道除垢機(jī)理，開發(fā)更為高效、實(shí)用的除垢技術(shù)，對提升油氣開采作業(yè)的效率和安全性具有重要意義。

當(dāng)前，常用的除垢方法包括機(jī)械除垢、化學(xué)除垢和物理除垢。其中，機(jī)械除垢方法包括管道清潔器除垢和機(jī)械刮除等。該類方法工作周期短，操作簡單且成本較低，但除垢效果欠佳，且除垢過程需要停機(jī)操作，影響經(jīng)濟(jì)效益?；瘜W(xué)除垢方法是指向管道內(nèi)注入酸性物質(zhì)實(shí)現(xiàn)除垢。該方法操作簡單，但酸性物質(zhì)極易腐蝕管壁，因此化學(xué)除垢方法通常和緩蝕劑聯(lián)合使用。物理除垢方法中最常見的是高壓水射流除垢，該技術(shù)利用射流的沖擊作用將水垢從管壁剝離3。與其他幾類除垢方法相比，射流除垢具有成本低且除垢效果好的優(yōu)勢，無需添加任何化學(xué)試劑，不會對地質(zhì)環(huán)境和生產(chǎn)過程產(chǎn)生影響4。大量的試驗(yàn)研究和實(shí)際應(yīng)用都表明，高壓水射流技術(shù)在工業(yè)除垢領(lǐng)域具有廣泛的應(yīng)用前景。劉錯沅1采用試驗(yàn)和數(shù)值模擬相結(jié)合的方式開展了連續(xù)管空化射流除垢作業(yè)參數(shù)的優(yōu)化研究。WANGA.M.等4通過單軸壓縮試驗(yàn)和巴西劈裂試驗(yàn)測量了煤層開采現(xiàn)場排水管道中水垢的抗拉強(qiáng)度，這是極少數(shù)關(guān)注水垢力學(xué)性能的研究。Y.TAMURA等5通過試驗(yàn)研究了除垢過程中射流結(jié)構(gòu)的變化情況和液滴速度沿噴射方向的衰減規(guī)律。趙金等通過試驗(yàn)研究了射流速度、噴嘴參數(shù)、移動速度和除垢劑對除垢效果的影響，研發(fā)了新型多孔射流工具，提高了油管的除垢效率。在現(xiàn)場應(yīng)用方面，Z.GUIMARAES等介紹了巴西海上油田使用新型連續(xù)管工具組合去除硫酸鋇結(jié)垢的案例。值得注意的是，現(xiàn)場引入的脈沖旋轉(zhuǎn)噴射工具顯著提升了除垢效率。斯倫貝謝公司針對地面管線的清潔問題，創(chuàng)新性地推出了連續(xù)管清洗技術(shù)及相應(yīng)的評估體系，目的是精確評估施工排量與液體攜帶雜質(zhì)能力的契合度[8。然而，利用高壓水射流進(jìn)行除垢作業(yè)時，操作人員往往僅憑借經(jīng)驗(yàn)來增加工作壓力和流量，這不僅使得射流除垢作業(yè)的可靠性降低，還導(dǎo)致了不必要的能耗和用水量的增加。因此，針對關(guān)鍵作業(yè)參數(shù)的研究對優(yōu)化射流除垢技術(shù)的工業(yè)應(yīng)用具有重要的指導(dǎo)意義。

噴嘴結(jié)構(gòu)作為射流除垢的關(guān)鍵部件之一，對除垢效果具有重要影響。一般來說，直射流噴嘴的沖擊能量集中，適用于切割作業(yè)或鉆小孔。與直射流相比，旋轉(zhuǎn)射流具有短距離沖擊面積大、壓力分布均勻等特點(diǎn)，在除垢解堵作業(yè)中展現(xiàn)出良好的適應(yīng)性[9-10]。針對旋轉(zhuǎn)射流的獨(dú)特性質(zhì)，王瑞和等[將旋轉(zhuǎn)射流應(yīng)用于破巖鉆井領(lǐng)域，率先通過試驗(yàn)和數(shù)值模擬相結(jié)合的方法研究了旋轉(zhuǎn)射流的結(jié)構(gòu)特性和流場分布。ZHANGD.等基于射流沖擊的水錘效應(yīng)建立了旋轉(zhuǎn)射流的數(shù)值模型。高軍紅等13采用數(shù)值模擬的方法比較了直射流、旋轉(zhuǎn)射流和直旋混合射流的鉆進(jìn)性能。LIH.等[14-15]對旋轉(zhuǎn)磨料水射流進(jìn)行了研究，提出了一種創(chuàng)新上端面結(jié)構(gòu)的改進(jìn)葉輪，改進(jìn)后的葉輪流動阻力小，擴(kuò)散角增大了 37% ，磨料的沖擊性能提高了 637.5% 。綜上，旋轉(zhuǎn)射流在鉆孔破碎領(lǐng)域的應(yīng)用已經(jīng)得到了廣泛且深入的研究；然而，其并未廣泛應(yīng)用于工業(yè)除垢領(lǐng)域，作用機(jī)理也尚不明確。

針對新疆吐哈油田油管的結(jié)垢問題，為了深入理解旋轉(zhuǎn)射流除垢過程中的流場特性，揭示不同工作參數(shù)對射流除垢性能的影響，筆者進(jìn)行了仿真研究。首先，建立了旋轉(zhuǎn)射流除垢的數(shù)值模型，并通過旋轉(zhuǎn)射流的流場擴(kuò)散性試驗(yàn)對該模型進(jìn)行了驗(yàn)證。然后，通過比較旋轉(zhuǎn)射流和常規(guī)直射流的流場差異，揭示了旋轉(zhuǎn)射流除垢過程中的漫流和沖蝕機(jī)理。在此基礎(chǔ)上，分析了泵壓、靶距和噴射角度等實(shí)際應(yīng)用場景中的工作參數(shù)對除垢性能的影響。該項(xiàng)研究旨在全面揭示旋轉(zhuǎn)射流除垢機(jī)理，為旋轉(zhuǎn)射流的應(yīng)用和高效除垢技術(shù)的開發(fā)提供理論指導(dǎo)。

數(shù)值模擬計算方法

1.1旋轉(zhuǎn)射流噴嘴

使用的旋轉(zhuǎn)射流噴嘴由噴嘴體和旋轉(zhuǎn)葉輪組成，如圖1所示。噴嘴體包括收縮段、出口段和擴(kuò)散段3部分。旋轉(zhuǎn)葉輪有5個葉片，葉輪螺旋角為 60^° 。噴嘴體和旋轉(zhuǎn)葉輪組合形成旋轉(zhuǎn)射流噴嘴，共有5條螺旋通道。噴嘴詳細(xì)結(jié)構(gòu)參數(shù)如表1所示。

1.2物理模型和邊界條件

為了模擬旋轉(zhuǎn)射流除垢的流場特性，建立了如圖2所示的三維射流沖擊幾何模型。該模型將水作為工作流體，包含旋轉(zhuǎn)射流噴嘴內(nèi)的流體域和外部流場。模型的邊界被設(shè)置為壓力入口和壓力出口，其詳細(xì)參數(shù)如表2所示。模型底部設(shè)為沖擊面，且被定義為絕熱無滑移壁面。射流沖擊底部面后沿徑向偏轉(zhuǎn)，并最終從側(cè)面流出，因此將外部流場區(qū)域的側(cè)面設(shè)置為出口。模型的其余邊界被設(shè)置為絕熱壁面。

網(wǎng)格質(zhì)量對計算結(jié)果有重要影響，良好的網(wǎng)格結(jié)構(gòu)可提高收斂速度并獲得更準(zhǔn)確的結(jié)果。采用FLUENTMeshing軟件對模型進(jìn)行離散化處理。對于幾何形狀規(guī)整且適合六面體網(wǎng)格劃分的區(qū)域，優(yōu)先采用六面體網(wǎng)格；其余部分則采用多面體網(wǎng)格劃分，且沿流體域壁面設(shè)置了邊界層網(wǎng)格。模型的網(wǎng)格策略如圖3所示。

1.3 基本假設(shè)

為了分析旋轉(zhuǎn)射流除垢過程中的流場特性，得到不同工作參數(shù)對射流除垢性能的影響規(guī)律，做出如下基本假設(shè)：

（1）采用三維模型，認(rèn)為計算過程射流是不

可壓縮流體，并以定常流動為主。（2）流場計算過程中，忽略溫度變化對流體的影響，并認(rèn)為流體和噴嘴壁面間不存在熱交換。（3）不考慮射流沖擊時的能量損失。

1.4控制方程

采用RNG k-ε 模型進(jìn)行湍流計算，控制方程包括連續(xù)性方程、動量方程及 k 和 ε 的輸運(yùn)方程。

連續(xù)性方程為：

式中： 、 ?_w 分別為流體在 、 z 方向上的速度分量， m/s 。

動量方程為：

式中： 為流體在讠 j 方向的空間坐標(biāo)變量， m δu_i?j u_j 為湍流時均速度沿讠 j 方向分量， m/s u_i^′ 、 u_j^′ 為湍流脈動速度沿讠 j 方向分量， m/s . ρ 為流體密度， kg/m³ p 為流體壓力，Pa; μ 為流體的動力（分子）黏度，Pa·s; 項(xiàng)為雷諾應(yīng)力項(xiàng)， Pa 。

為獲得準(zhǔn)確的模擬結(jié)果，采用RNG k-ε 模型進(jìn)行計算。該模型通過引入額外的修正項(xiàng)，能夠更好地處理湍流流動中的非均勻性。與Standard k-ε 模型相比，RNG k-ε 模型在捕捉湍流中的旋轉(zhuǎn)效應(yīng)方面表現(xiàn)更佳，因此廣泛應(yīng)用于射流沖擊和渦旋流動的模擬[16]。k和 ε 的輸運(yùn)方程為：

G_k+G_b-ρε-Y_M+S_k

式中： k 和 ε 分別為湍流動能 （m²/s²） ）和它的耗散率（m²/s³） ； μ_eff 為有效黏度， Pa?s ，是層流黏度和湍流黏度的綜合體現(xiàn)； G_k 為由平均速度梯度產(chǎn)生的湍流動能， Pa/s G_b 為由浮力產(chǎn)生的湍流動能， Pa/s ; Y_M 為可壓縮湍流中波動膨脹對總耗散率的貢獻(xiàn)， Pa/s C_1ε 、 C_2ε 和 C_3ε 均為常數(shù)，其中 C_1ε 為1.42， C_2ε 為1.68; α_k 為 k 的逆有效普朗特數(shù)； α_s 為 ε 的逆有效普朗特數(shù)； R_s 為附加湍流耗散項(xiàng)， Pa/s² ; S_k 和 S_ε 為用戶自定義源項(xiàng)， Pa/s ， Pa/s² ，在本模型中被忽略。

逆有效普朗特數(shù)通過下式計算：

式中： α_a 為無量綱變量，表示當(dāng)前逆有效普朗特數(shù)；

a₀ 為常數(shù)，取值為1.0； μ_mol 為層流黏度，Pa·s。

在高雷諾數(shù)極限下，此時湍流占主導(dǎo)， μ_eff 遠(yuǎn)大于 μ_mol ，式（5）右邊項(xiàng)趨近0，此時逆有效普朗特數(shù)約為1.393。

RNG k-ε 模型和Standard k-ε 模型之間的主要區(qū)別在于 ε 的輸運(yùn)方程中的附加項(xiàng) R_ε°R_ε 的計算式為：

其中：

η=S_k/ε

式中： C_Lμ=0. 089 5;η_L0=4. 38;β=0.012

將該附加項(xiàng)代入方程（4）中，原輸運(yùn)方程變?yōu)椋?/p>

其中：

式中： C_2ε 、 C_2ε^* 為修正前、后的經(jīng)驗(yàn)常數(shù)。湍流黏度的微分方程為：

其中：

式中： 為有效黏度之比，即有效黏度和動力黏度之比； C_v 為常數(shù)，取值趨近 100 。

在高雷諾數(shù)極限下，湍流黏度 μ_t 計算式為：

湍流易受旋流或者渦流的影響，RNG k-ε 模型可以通過適當(dāng)修改湍流黏度來考慮旋流或者渦流的影響，即：

式中： μ_?10 為不考慮渦流修正情況下的湍流黏度值； 為特征渦流數(shù)； a_s 為一個渦流常數(shù)，其值取決于流動是以渦流為主導(dǎo)還是僅有輕微渦流，默認(rèn)情況下取0.07。

2 模擬方法驗(yàn)證

2.1 網(wǎng)格和時間步長獨(dú)立性驗(yàn)證

為最大限度保證模擬結(jié)果與網(wǎng)格數(shù)量和時間步長無關(guān)，選擇了距離噴嘴 10mm 處的流體速度作為評價指標(biāo)，并對其進(jìn)行了獨(dú)立性驗(yàn)證（見圖4）。

使用了不同尺寸的網(wǎng)格對模型進(jìn)行劃分，得到網(wǎng)格數(shù)量分別為57516、89439、123849、178543和265999等5種方案。計算結(jié)果如圖4a所示。從圖4a可以看出，當(dāng)網(wǎng)格數(shù)量從178543增加至265999時，流體速度變化的程度極小。綜合考慮模擬精度和效率，后續(xù)計算選用數(shù)量為178543的網(wǎng)格劃分方案。針對時間步長，選擇了0.001、0.01、0.1、0.5和1s等5種方案進(jìn)行模擬。計算結(jié)果如圖4b所示。從圖4b可以看出，當(dāng)時間步長小于0.01s時，流體速度變化的程度極小。因此，后續(xù)計算案例中的時間步長設(shè)置為0.01s。

2.2計算方法驗(yàn)證

在開展旋轉(zhuǎn)射流模擬研究之前，將這里提出的數(shù)值模擬方法與周秋成等的試驗(yàn)結(jié)果進(jìn)行了對比，以驗(yàn)證計算方法的可靠性。鑒于旋轉(zhuǎn)射流的切向速度與射流擴(kuò)散形態(tài)密切相關(guān)，通過圖5所示的試驗(yàn)裝置開展了旋轉(zhuǎn)射流流場的擴(kuò)散性試驗(yàn)。

在驗(yàn)證案例中，旋轉(zhuǎn)射流噴嘴葉輪旋轉(zhuǎn)角度設(shè)定為 360^° 。模型入口采用 10m/s 的速度入口條件，模型出口采用 0.1MPa 的壓力出口條件。上述試驗(yàn)中拍攝的淹沒條件下流場擴(kuò)散照片如圖6a所示。從圖6a可以看出，射流擴(kuò)散角約為 18^° 。圖6b為驗(yàn)證案例中淹沒條件下旋轉(zhuǎn)射流的速度云圖。在保證試驗(yàn)結(jié)果照片中的比例尺與驗(yàn)證案例結(jié)果圖像尺寸基本一致的條件下，其射流擴(kuò)散角約為 16.9^° ，與試驗(yàn)所得流場的擴(kuò)散情況基本相符，說明這里提出的計算方法能夠有效預(yù)測旋轉(zhuǎn)射流流場的基本特性。

3 結(jié)果與討論

3.1旋轉(zhuǎn)射流的漫流作用

為了探究旋轉(zhuǎn)射流除垢的作用機(jī)制，分別模擬了旋轉(zhuǎn)射流和直射流的沖擊過程。雖然直射流與旋轉(zhuǎn)射流的噴嘴結(jié)構(gòu)不同，但兩者具有相同的壓力條件。圖7為2種射流的速度云圖。由圖7可知，流體通過噴嘴加速，根據(jù)伯努利原理，此時流體的壓力能轉(zhuǎn)化為動能，并在噴嘴出口段速度達(dá)到最大。圖7中的旋轉(zhuǎn)射流最大速度為257.7m/s ，直射流的最大速度為 247.3m/s 。高速射流從噴嘴噴出后充分發(fā)展，對水垢產(chǎn)生沖擊作用。隨后，射流方向沿徑向偏轉(zhuǎn)，形成徑向射流區(qū)。流體的徑向流動產(chǎn)生漫流，能夠?qū)λ副砻娈a(chǎn)生侵蝕效應(yīng)，有利于增大除垢面積。

圖8為射流速度沿軸向的分布規(guī)律。其中 L₁ L₂ 、 L₃ 、 L₄ 和 L₅ 分別代表混合段、收縮段、出口段、擴(kuò)散段和自由射流區(qū)域。從圖8可見：隨著軸向距離的增大，射流速度先增大后減小，但在 L₁ 和 L₂ 區(qū)域內(nèi)，旋轉(zhuǎn)射流的速度明顯高于直射流；在 L₃ 區(qū)域內(nèi)，由于噴嘴流道的截面積不斷減小，旋轉(zhuǎn)射流和直射流的速度迅速增長；高壓射流逐步轉(zhuǎn)化為高速射流，并在 L₄ 區(qū)域內(nèi)達(dá)到峰值，此時旋轉(zhuǎn)射流的速度仍高于常規(guī)直射流，在評估射流的沖擊能力時，沿軸向的速度起著決定性作用，這意味著旋轉(zhuǎn)射流的沖擊強(qiáng)度得以保證；在 L₅ 區(qū)域，射流充分發(fā)展，該階段旋轉(zhuǎn)射流的速度持續(xù)衰減，而直射流能夠保持穩(wěn)定流動的狀態(tài)，這是由于旋轉(zhuǎn)射流的夾帶效應(yīng)驅(qū)動射流與周圍流體相互擾動。與直射流相比，旋轉(zhuǎn)射流較大的切向速度增強(qiáng)了其膨脹率和湍流脈動強(qiáng)度，導(dǎo)致其速度迅速下降[18]。在近壁面區(qū)域，射流速度迅速下降，此時射流的動能逐漸向壓力能轉(zhuǎn)化；在滯止區(qū)，射流沿軸向方向的速度變?yōu)?。

旋轉(zhuǎn)射流的速度分布特性由噴嘴結(jié)構(gòu)決定，且直接影響除垢效率。為了確認(rèn)射流內(nèi)部流場的速度變化規(guī)律，選擇了噴嘴內(nèi)部多個不同位置的截面，如圖9所示。其中，平面1為旋轉(zhuǎn)葉輪部分的截面，平面2為噴嘴混合段的截面，平面3為噴嘴收縮段的截面，平面4為噴嘴出口段上部的截面，平面5為噴嘴出口段下部的截面，平面6為噴嘴擴(kuò)散段的截面。

從圖9可知：平面1的速度矢量顯示了流體通過5個葉輪槽時的狀態(tài)，總體來看，5條流道呈現(xiàn)出相似的分布特征；射流沿順時針方向旋轉(zhuǎn)并出現(xiàn)了明顯的速度梯度分布，此時，旋轉(zhuǎn)射流保持低速狀態(tài)；平面2的速度矢量顯示了流體離開旋轉(zhuǎn)葉輪進(jìn)人混合段后的狀態(tài)，在該流場中，中心速度出現(xiàn)并保持旋轉(zhuǎn)狀態(tài)，表明不同流道內(nèi)的射流開始匯聚，在混合段射流發(fā)展的初始階段，仍可清晰地辨識出由每個葉輪槽產(chǎn)生的獨(dú)立旋轉(zhuǎn)射流，它們之間的邊界明顯可辨；平面3的速度矢量顯示了流體進(jìn)入收縮段后的狀態(tài)，該區(qū)域不同流道內(nèi)的射流匯聚完畢，邊界也不再明顯，但仍保持旋流狀態(tài)，在射流發(fā)展的過程中，其速度分布呈現(xiàn)出從中心向外逐漸減小的趨勢，這歸因于旋轉(zhuǎn)射流通過收縮段錐形壁面時逐漸受到壓縮，使得能量向中心區(qū)域集中；平面4和平面5顯示了射流在出口段的狀態(tài)，該區(qū)域射流的速度達(dá)到最大，在壓差為 29.9MPa 的條件下，旋轉(zhuǎn)射流的最大速度達(dá)到 255.5m/s ，表明高壓射流已經(jīng)轉(zhuǎn)化為了高速射流，且射流高速區(qū)主要集中于射流核心；平面6的速度矢量顯示了射流在噴嘴擴(kuò)散段的狀態(tài)，由于旋轉(zhuǎn)射流具有更強(qiáng)的擴(kuò)散系數(shù)，因而易與周圍流體發(fā)生相互作用，此時，射流擴(kuò)散角不斷增大，其核心區(qū)域也會逐漸膨脹。

漫流速度是指射流在沖擊壁面中心沿徑向流動時表現(xiàn)出的速度特性，反映了流體從中心向四周擴(kuò)散的動態(tài)過程。圖10為旋轉(zhuǎn)射流和直射流在沖擊壁面處的速度矢量分布。從圖10可見：對于直射流，其速度分布呈規(guī)則的環(huán)形向外擴(kuò)展，但高速區(qū)域的面積較為有限；相比之下，旋轉(zhuǎn)射流的速度分布呈不規(guī)則的橢球形，這與其射流特殊的旋流狀態(tài)有關(guān)。旋轉(zhuǎn)射流的漫流作用顯著擴(kuò)大了其高速區(qū)的漫流作用范圍，使其在相同條件下可形成大于直射流的除垢覆蓋面積。

3.2旋轉(zhuǎn)射流的沖蝕作用

圖11a為2種射流在沖擊壁面處的切向速度分布。從圖11a可以看出，切向速度沿徑向先增大后減小，且旋轉(zhuǎn)射流的切向速度值遠(yuǎn)大于直射流。切向速度表征了流體的離心擴(kuò)散能力。這種離心作用有助于將垢體從管道表面剝離。峰值處即為水垢的最大去除區(qū)域，其切向速度的峰值出現(xiàn)在距離壁面中心 2.2mm 處。圖11b為旋轉(zhuǎn)射流和直射流的徑向速度分布。從圖11b可見，2種射流具有相似的速度分布形式。綜上，在旋轉(zhuǎn)葉輪的導(dǎo)流作用下，旋轉(zhuǎn)射流的切向速度和軸向速度共同作用，形成三維速度分布，使得旋轉(zhuǎn)射流具有更強(qiáng)的除垢效果。此外，盡管徑向速度對除垢面積的影響較大，但其未能突出旋轉(zhuǎn)射流的特性，因此在后續(xù)的研究中不予考慮。

圖12為旋轉(zhuǎn)射流和直射流的壓力分布情況。綜合圖7和圖12可以觀察到，射流壓力的變化是一個動態(tài)過程，即當(dāng)高壓射流從噴嘴入口流入時，噴嘴內(nèi)的流體在初始階段大部分維持了高壓狀態(tài)；隨著射流進(jìn)入收縮段，壓力能逐漸轉(zhuǎn)化為動能；當(dāng)射流到達(dá)近壁面區(qū)域后，其動能再次向壓力能轉(zhuǎn)化，形成了滯止區(qū)，滯止區(qū)的壓力顯著高于環(huán)境壓力，且沿沖擊壁面的方向產(chǎn)生了正壓力梯度。射流的發(fā)展是動能和壓力能相互轉(zhuǎn)化的過程。水垢表面不均勻的壓力分布能夠誘使其產(chǎn)生拉伸應(yīng)力并借此實(shí)現(xiàn)除垢。

沖擊壓力是射流除垢過程中的關(guān)鍵因素。提取沖擊壁面上流體的壓力分布，結(jié)果如圖13所示。從圖13可見：旋轉(zhuǎn)射流的最大壓力達(dá)到11.11MPa ，與直射流相差約 1MPa ，這說明旋轉(zhuǎn)射流具有與直射流相接近的沖擊破碎能力；隨著與壁面中心距離的增加，2種射流的壓力迅速減?。慌c直射流相比，旋轉(zhuǎn)射流在近壁面處的壓力下降趨勢較為平緩，這意味著其能夠在更廣泛的區(qū)域內(nèi)維持較高的壓力分布，表明旋轉(zhuǎn)射流沖擊的作用范圍更大。

在射流作用下，水垢會產(chǎn)生剪切應(yīng)力和拉伸應(yīng)力，當(dāng)這些應(yīng)力超過水垢的極限強(qiáng)度時，水垢將發(fā)生破壞，具體表現(xiàn)為剪切裂紋和拉伸裂紋的形成。射流的持續(xù)作用會導(dǎo)致裂紋進(jìn)一步擴(kuò)展和連通，最終在宏觀上表現(xiàn)為水垢的體積破碎，即沖蝕破壞。圖14為旋轉(zhuǎn)射流和直射流在沖擊壁面處的剪切應(yīng)力云圖。從圖14可見：在直射流作用下，壁面中心的應(yīng)力水平較低，其最大剪切應(yīng)力位于滯止區(qū)外圍的環(huán)形區(qū)域，這種應(yīng)力分布與射流的流動特性密切相關(guān)。由于近壁面處的徑向流動導(dǎo)致壓力梯度和流體速度發(fā)生劇烈變化，高水平的剪切應(yīng)力集中在切向速度較大的位置。不同于直射流，旋轉(zhuǎn)射流的漫流特征使其剪切應(yīng)力在壁面上的分布呈現(xiàn)出不規(guī)則的橢球形，并且不存在應(yīng)力水平突出的環(huán)形區(qū)域，這是由于旋轉(zhuǎn)射流到達(dá)壁面時速度耗散更為顯著；此外，旋轉(zhuǎn)射流同時具有較大的切向速度和徑向速度，這2種速度的共同作用導(dǎo)致剪切應(yīng)力在壁面上分布不均勻。

提取了沿壁面徑向的應(yīng)力分布，結(jié)果如圖15所示。從圖15可見：隨著徑向距離的增加，剪切應(yīng)力在局部呈現(xiàn)出增大趨勢，并分別在距離壁面中心2.2和 4.1mm 處達(dá)到峰值，這與不同射流沿壁面的速度分布特性有關(guān)，隨后剪切應(yīng)力逐漸減??；與直射流相比，旋轉(zhuǎn)射流在沖擊壁面處的剪切應(yīng)力下降趨勢表現(xiàn)得更為平緩，說明高水平應(yīng)力作用的范圍更廣，有利于增大射流除垢面積。

3.3工作參數(shù)對除垢性能的影響

在優(yōu)化旋轉(zhuǎn)射流噴嘴的性能方面，泵壓、靶距以及噴射角度等工作參數(shù)的合理選用是提升旋轉(zhuǎn)射流除垢效率的關(guān)鍵。在此，將上述3類參數(shù)作為研究旋轉(zhuǎn)射流除垢性能的特征參數(shù)，選定射流在近壁面處的切向速度、壁面沖擊壓力和剪切應(yīng)力作為評價指標(biāo)，基于相同的仿真條件開展研究。

3.3.1 泵壓

圖16為在靶距為 15mm ，泵壓分別為30、35、40和 45MPa 條件下射流的切向速度、壁面沖擊壓力和剪切應(yīng)力的變化。從圖16可知，泵壓的增加顯著提高了旋轉(zhuǎn)射流在近壁面處的切向速度，當(dāng)泵壓從 30MPa 增加至 45MPa 時，切向速度的峰值提升了 59.1% 。泵壓的升高也導(dǎo)致了壁面沖擊壓力的增大。在射流的發(fā)展過程中，更高的泵壓條件能夠賦予流體更大的速度。此外，泵壓對旋轉(zhuǎn)射流的湍流強(qiáng)度也有顯著影響，隨著泵壓的增加，湍流強(qiáng)度增強(qiáng)，流體混合和動量交換更加劇烈，從而使流體在不同方向上的速度分量增大。泵壓對壁面剪切應(yīng)力也有促進(jìn)作用，增大泵壓會影響流體的黏性和剪切效應(yīng)。在高壓條件下，流體的動能更大，流體顆粒之間的剪切力更大，對垢體的沖蝕效果更強(qiáng)。總的來說，在射流除垢作業(yè)中，增大泵壓對提升旋轉(zhuǎn)射流的除垢性能具有積極影響，但過高的泵壓可能會影響泵體、管柱和噴嘴的機(jī)械強(qiáng)度?，F(xiàn)場應(yīng)用應(yīng)基于除垢效果的具體需求，合理選擇并調(diào)整泵壓條件。

3.3.2 靶距

圖17為在泵壓為 30MPa ，靶距分別為15、20、25和 30mm 條件下射流的切向速度、壁面沖擊壓力和剪切應(yīng)力的變化。從圖17可知，當(dāng)靶距從 15mm 增加至 25mm 時，近壁面處旋轉(zhuǎn)射流的切向速度增大，此時，若繼續(xù)增加靶距，射流的切向速度開始降低。這是由于旋轉(zhuǎn)射流的切向速度依賴于其初始旋轉(zhuǎn)動量。在射流發(fā)展過程中，旋轉(zhuǎn)動量會因流體內(nèi)部的黏性和外部空氣阻力而逐漸衰減。隨著靶距增加，這種衰減效應(yīng)愈加顯著，

Fig.17Influence of diferent target distances ontangential velocity，impactpressure and shear stress致使射流在到達(dá)壁面時的切向速度降低，這表明靶距對旋轉(zhuǎn)射流切向速度的影響存在一個最優(yōu)值。靶距的增大也會導(dǎo)致壁面沖擊壓力的降低，這是由于射流在傳播過程中會逐漸擴(kuò)散，能量逐步衰減，從而降低了沖擊壓力。旋轉(zhuǎn)射流的擴(kuò)散效應(yīng)也會降低壁面的剪切應(yīng)力。總體而言，在射流除垢作業(yè)中，過大的靶距不利于除垢性能的提升。

3.3.3噴射角度

圖18為在泵壓 30MPa ，靶距 15mm ，噴射角度分別為 0^° 、 5^° 、 10^° 和 15^° 條件下射流的切向速度、壁面沖擊壓力和剪切應(yīng)力的變化。從圖18可知：噴射角度過大時，旋轉(zhuǎn)射流在近壁面處的切向速度會逐漸減小，而壁面沖擊壓力基本不受噴射角度的影響；噴射角度和壁面剪切應(yīng)力之間存在顯著的正相關(guān)關(guān)系，隨著噴射角度的增大，與壁面平行的切向動量分量減小，導(dǎo)致近壁面處射流的切向速度減小；噴射角度的變化也會引起靶面射流的分流現(xiàn)象和水墊效應(yīng)，使壁面剪切應(yīng)力分布的不均勻性更加顯著[19]。較大的噴射角度會增強(qiáng)壁面射流的湍流強(qiáng)度，從而增加壁面處的局部剪切應(yīng)力。在實(shí)際應(yīng)用中，適當(dāng)增大噴射角度，可以增強(qiáng)旋轉(zhuǎn)射流的除垢性能。

4結(jié)論

探討了旋轉(zhuǎn)射流除垢的作用機(jī)理，通過數(shù)值模擬手段，系統(tǒng)分析了旋轉(zhuǎn)射流的漫流效應(yīng)和沖蝕作用。針對實(shí)際應(yīng)用中泵壓、靶距和噴射角度等關(guān)鍵工作參數(shù)，采用切向速度、沖擊壓力以及剪切應(yīng)力作為評價指標(biāo)，對射流除垢性能進(jìn)行了系統(tǒng)評估。主要結(jié)論如下：

（1）射流的速度分布特性由噴嘴結(jié)構(gòu)決定，并會對除垢效率有顯著影響。在相同的泵壓條件下，旋轉(zhuǎn)射流的高速區(qū)更易耗散；在近壁面區(qū)域，旋轉(zhuǎn)射流的速度分布呈不規(guī)則的橢球形，但其高速區(qū)的影響范圍更大，因而具有更大的除垢面積。（2）旋轉(zhuǎn)射流除垢過程中具有三維速度分布。在近壁面區(qū)域，其同時具有較大的切向速度和徑向速度，這2種速度的共同作用導(dǎo)致剪切應(yīng)力在壁面上分布不均勻。旋轉(zhuǎn)射流具有和直射流相近的沖擊破碎能力。（3）以切向速度、壁面沖擊壓力和剪切應(yīng)力作為評價指標(biāo)時，提高泵壓能夠有效增強(qiáng)旋轉(zhuǎn)射流的除垢效果。靶距對除垢性能的影響存在一個最優(yōu)值，適當(dāng)增大噴射角度，可以有效增強(qiáng)旋轉(zhuǎn)射流的除垢能力。

參考文獻(xiàn)

[1]劉鍇沅．連續(xù)管套管除垢空化射流噴頭優(yōu)化設(shè)計與應(yīng) 用[J].石油機(jī)械，2024，52（2）：44-49. LIUKY.Optimized design and application of cavitation jetnozzle for CTcasingdescaling[J].China Petroleum Machinery，2024，52（2）：44-49.

［2]李善建，何浩軒，王澤坤，等．氣井井筒堵塞原因分 析及解堵工藝研究進(jìn)展［J]．西安石油大學(xué)學(xué)報（自 然科學(xué)版），2024，39（1）：56-65. LISJ，HEHX，WANG ZK，et al.Analysis of reasons for wellbore blockage in gas wellsand research progressin blockage removal technology[J].Journal ofXi’an Shiyou University（NaturalScience），2024，39（1）：56-65.

[3]ZHAO HJ，HUANG YH，DENG S，et al.Research progress on scaling mechanism and anti-scaling technologyof geothermal well system[J]. Journal of Dispersion Science and Technology，2022，44 （9）：1657-1670.

[4] WANGAM，SUNPF，ZHANGZW，etal.Numerical simulation and experimental study on high-pressure water jetdescalingin coal mine drainage pipeline[J]. Iranian Journal ofScience and Technology， Transactionsof Mechanical Engineering，2022，46（1）：113-130.

[5]TAMURA Y，UEOKA S，KIMURAY，et al. Influence of injection distance on water droplet behavior in high pressure descaling［J].ISIJ International，2O20，60 （1）： 128-135.

[6］趙金，趙星，蔡鵬，等．連續(xù)管多孔噴嘴射流工具研 制與清洗參數(shù)優(yōu)化[J]．石油鉆探技術(shù)，2023，51 （3）：83-89. ZHAO J，ZHAO X，CAI P，et al. Development of coiled tubing jet tools with multi-hole nozzesand cleaning parameter optimization[J].Petroleum Drilling Techniques，2023，51（3）：83-89.

[7]GUIMARAES Z，F(xiàn)RANCA AB，DUQUE LH，et al. Case histories of barium sulfate scale removal in offshore wells， Brazil，using anew enginered combination of coiled-tubing tools[C]//SPE/ICoTA Coiled Tubing and Well Intervention Conference and Exhibition.The Woodlands，Texas，U.S. A：SPE，2007： SPE 107063-MS.

[8]HASSIG FONSECA S，SERRANO D，VILLACRESC， et al.Methodology to remediate and evaluate surface flowline capacity with coiled tubing cleanouts[J]. SPE Productionamp;Operations，2022，37（3）：520-532.

[9]WU DY，ZHANG HD，ZHANG SH，et al. Determining whether the swirling jet is applicable in botom-hole cleaning：a CFD study[J]. Journal of Petroleum Science and Engineering，2022，208，PartA：109317.

[10］竺彪，吳德松，郭宏峰，等．大尺寸旋轉(zhuǎn)水射流解 堵工具的研制及特性測試[J]．石油機(jī)械，2015，43 （4）：105-108. ZHU B，WU D S，GUO HF，et al. Design of alarge size rotating water jet tool for wellbore plug removal and its performance test[J]. China Petroleum Machinery， 2015，43（4）： 105-108.

[11］王瑞和，周衛(wèi)東，沈忠厚，等．旋轉(zhuǎn)射流破巖鉆孔 機(jī)理研究[J]．中國安全科學(xué)學(xué)報，1999，9（增刊 1）：1-5. WANGRH，ZHOU WD，SHEN ZH，et al.Study on mechanism of Rock-Breaking by swirling water jet[J]. China Safety Science Journal，1999，9（Sup.1）：1-5.

[12]ZHANG D，GE Z L，ZHOU Z，et al. Water hammer characteristics of self-rotating cascade jets impacting on rock ［J].Rock Mechanicsand Rock Engineering， 2024，57 （4）：2681-2693.

[13]高軍紅，付必偉，董宗正．自進(jìn)式直旋混合噴嘴的 鉆進(jìn)性能研究[J]．石油機(jī)械，2022，50（9）： 72-78. GAO JH，F(xiàn)U B W，DONG Z Z. Drillng performance of self-propelled straight-rotary mixing nozzle[J]. China Petroleum Machinery，2022，50（9）：72-78.