新型射流振蕩減摩阻工具設(shè)計(jì)及內(nèi)部流場(chǎng)特性數(shù)值模擬分析與實(shí)驗(yàn)驗(yàn)證*

田家林 何虹志 楊 琳 楊應(yīng)林 宮學(xué)成 胡志超 李居瑞

（1.西南石油大學(xué)機(jī)電工程學(xué)院 四川成都 610500； 2.中海油海洋石油工程（青島）有限公司 山東青島 266520）

在油氣井實(shí)際鉆井過程中，由于地層環(huán)境差異導(dǎo)致井下工況條件非常復(fù)雜，隨著新型鉆井技術(shù)的發(fā)展，各種新型鉆井工藝不斷被開發(fā)應(yīng)用，導(dǎo)致鉆柱與井壁之間的鉆井摩阻影響因素呈現(xiàn)多樣性，使井下工具安全評(píng)價(jià)、鉆井過程減摩降阻與提速增效面臨新的挑戰(zhàn)[1-4]。越來越多的學(xué)者與研究機(jī)構(gòu)也認(rèn)識(shí)到該問題的重要性，經(jīng)過不斷探索研究，研制出不同類型的井下減摩降阻工具，應(yīng)用于石油天然氣勘探開發(fā)領(lǐng)域，以實(shí)現(xiàn)鉆井過程中減摩降阻與提速增效[5-8]。

現(xiàn)有的減摩阻工具通常利用軸向沖擊與周向沖擊作用，將鉆井液的液壓能轉(zhuǎn)變?yōu)楣ぞ叩臋C(jī)械能，使工具產(chǎn)生振動(dòng)，通過在鉆具組合中加入減摩阻工具，能夠有效降低鉆井摩阻，提高鉆進(jìn)效率[9-10]。軸向減摩阻工具通過工具工作產(chǎn)生軸向振動(dòng)，形成軸向沖擊力，改變鉆柱軸向受力條件，實(shí)現(xiàn)軸向振動(dòng)減摩作用，相關(guān)工具有水力振蕩器、防滯動(dòng)工具、連續(xù)管減摩器等[11-14]；周向減摩阻工具產(chǎn)生的周期性扭矩傳遞到鉆頭，為鉆頭破巖提供一個(gè)額外的周期性扭矩，可提高鉆頭破巖效率，常見的周向工具為扭力沖擊器[15]。但是常用的減摩阻工具存在結(jié)構(gòu)復(fù)雜、耐沖蝕性能差、易受地層環(huán)境影響等缺點(diǎn)，限制了減摩阻工具的推廣應(yīng)用。

因此，為簡(jiǎn)化工具結(jié)構(gòu)，提高工具使用壽命，筆者設(shè)計(jì)了一種射流振蕩減摩阻工具，并結(jié)合該工具的工作原理，進(jìn)行了射流短節(jié)內(nèi)部流體流動(dòng)特性研究，通過CFD-Post處理得到數(shù)值模擬結(jié)果，最后通過實(shí)驗(yàn)測(cè)試，驗(yàn)證了數(shù)值模擬結(jié)果的正確性。研究結(jié)果可為射流振蕩減摩阻工具結(jié)構(gòu)優(yōu)化設(shè)計(jì)提供幫助，也可為進(jìn)一步展開射流振蕩減摩阻工具相關(guān)理論研究提供參考。

1 新型射流振蕩減摩阻工具結(jié)構(gòu)設(shè)計(jì)

新型射流振蕩減摩阻工具由減震器和射流振蕩器兩部分組成，其具體結(jié)構(gòu)如圖1所示。射流振蕩器安裝在減震器前端，在射流振蕩器內(nèi)部設(shè)置射流短節(jié)。由于射流短節(jié)存在特定形狀的射流腔，當(dāng)鉆井液進(jìn)入射流振蕩器后，通過在射流腔內(nèi)循環(huán)，按照特定的循環(huán)路線完成循環(huán)后，從流體出口流出，形成鉆井液工作壓降。

圖1 射流振蕩減摩阻工具結(jié)構(gòu)示意圖Fig.1 Schematic diagram of jet oscillation friction reducing tool

當(dāng)流體在射流短節(jié)內(nèi)循環(huán)時(shí)，若首先形成向上的入射流，經(jīng)過射流短節(jié)上部輸入通道進(jìn)入振蕩室，產(chǎn)生順時(shí)針渦流；當(dāng)完成依附壁面轉(zhuǎn)換后，流體經(jīng)過下部輸入通道進(jìn)入振蕩室，產(chǎn)生逆時(shí)針渦流，形成周期性的附壁與切換，形成流體的康達(dá)效應(yīng)。在工作過程中，射流振蕩器將產(chǎn)生周期性的脈沖壓力波，并傳遞到射流振蕩器后端的減震器上，在脈沖壓力的作用下，減震器內(nèi)部串聯(lián)的碟簧組出現(xiàn)周期性地壓縮與復(fù)位，使工具產(chǎn)生周期性的軸向振動(dòng)，從而實(shí)現(xiàn)減摩阻作用。相對(duì)現(xiàn)有常用的水力振蕩器、射流工具等減摩工具，本文提出的射流振蕩減摩阻工具的射流振蕩器內(nèi)部無運(yùn)動(dòng)件，依靠流體在循環(huán)流道內(nèi)自動(dòng)切換產(chǎn)生射流振蕩效果，從而實(shí)現(xiàn)減摩阻功能，而不是依靠運(yùn)動(dòng)部件的轉(zhuǎn)動(dòng)來改變流體的流動(dòng)特性，因此結(jié)構(gòu)更簡(jiǎn)單。

2 流場(chǎng)特性數(shù)值模擬分析

2.1 模型建立

2.1.1 幾何模型建立

根據(jù)射流振蕩減摩阻工具的結(jié)構(gòu)設(shè)計(jì)，按照表1中的主要結(jié)構(gòu)參數(shù)，建立射流短節(jié)內(nèi)部流體計(jì)算域結(jié)構(gòu)模型，結(jié)構(gòu)模型包括入水管流域、上下控制端口流域、上下輸入流道流域、上下直線流道流域、彎曲流道流域、振蕩室流域與出入口流域等，如圖2所示。

表1 射流短節(jié)主要結(jié)構(gòu)參數(shù)Table1 Main structural parameters of jet joint

圖2 射流短節(jié)內(nèi)部流體計(jì)算域結(jié)構(gòu)模型Fig.2 Computational domain structure model of fluid in jet joint

2.1.2 控制方程

射流振蕩減摩阻工具以鉆井液為輸送介質(zhì)，鉆井液進(jìn)入射流短節(jié)后，經(jīng)過內(nèi)部流道，射流振蕩器產(chǎn)生振蕩作用。忽略鉆井液所受到的體積力，并且鉆井液不與外界發(fā)生熱交換。結(jié)合射流短節(jié)內(nèi)流體的流動(dòng)狀態(tài)，選擇RNGk-ε湍流模型，該模型由瞬時(shí)納維-斯托克斯方程（N-S）演化而來，針對(duì)N-S方程進(jìn)行時(shí)間平均后即可得到雷諾平均N-S方程[16]：

式（1）中：ρ為流體密度，kg／m3；t為時(shí)間，s；x i、x j、x k分別表示分量符號(hào)（i=1、2、3，j=1、2、3，k=1、2、3）分別表示雷諾平均速度分量，m／s；μ為動(dòng)力黏度，N·s／m2ˉp為平均壓強(qiáng)，Pa′分別表示平均脈動(dòng)流速分量，m／s；δij為克羅內(nèi)克張量分量。

通過布西涅斯克假設(shè)，可以將雷諾應(yīng)力與平均速度梯度相關(guān)聯(lián)起來，其雷諾應(yīng)力模型為[16]

其中

式（2）、（3）中：μt為湍流黏度，Pa·s；k為湍流動(dòng)能，m2／s2；ε為湍流動(dòng)能耗散率，m2／s3；C為經(jīng)驗(yàn)常數(shù)，取C=0.09。

在RNGk-ε模型中，湍動(dòng)能k和湍流動(dòng)能耗散率ε分別由式（4）、（5）表示[17]

式（4）、（5）中：C1、C2為經(jīng)驗(yàn)常數(shù)，取C1=1.44、C2=1.92；σk為湍流動(dòng)能的普朗特?cái)?shù)，取σk=1；σε為湍流動(dòng)能耗散率的普朗特?cái)?shù)，取σε=1.3。

2.1.3 定解條件

在同一網(wǎng)格模型中，改變?nèi)肟诹髁看笮。?、4、6、8和10 L／s）進(jìn)行計(jì)算，得到射流振蕩減摩阻工具工作時(shí)射流短節(jié)內(nèi)流體部分參數(shù)的變化規(guī)律。在計(jì)算過程中，計(jì)算模型的入口定義為速度邊界條件，出口定義為壓力出口，出口壓力定義為一個(gè)大氣壓。計(jì)算模型入口、出口湍流強(qiáng)度和水力直徑可以通過式（6）計(jì)算得到[18]。

式（6）中：I為湍流強(qiáng)度，%；Re為雷諾數(shù)；Dh為水力直徑，m；υ為流體流速，m／s；A為過流面積，m2；Cl為過流截面周長(zhǎng)，m。

不同入口流量條件下的邊界條件參數(shù)見表2。

表2 不同入口流量條件下的邊界條件參數(shù)Table2 Boundary condition parameters under different inlet flow conditions

2.2 模擬計(jì)算結(jié)果分析

通過數(shù)值模擬計(jì)算，得到不同入口流量條件下射流振蕩減摩阻工具射流短節(jié)內(nèi)部流體工作參數(shù)，包括壓力場(chǎng)、速度場(chǎng)等計(jì)算結(jié)果，并對(duì)不同位置的流體工作參數(shù)進(jìn)行監(jiān)測(cè)和分析，主要對(duì)流體入口、流體出口、分流中心點(diǎn)以及流體中心截面的參數(shù)變化進(jìn)行監(jiān)測(cè)，具體位置如圖3所示。

2.2.1 速度場(chǎng)分析

結(jié)合數(shù)值計(jì)算過程，流體在射流短節(jié)內(nèi)的流動(dòng)存在附壁與切換現(xiàn)象，以入口流速為3.056 m／s的初始條件為例，整個(gè)附壁與切換過程中的流速變化如圖4所示。由圖4可以看出，在初始階段，流體經(jīng)過入水管直線流道進(jìn)入擴(kuò)散段后，首先完成均勻射流，隨著時(shí)間的增加，流體開始沿著上部輸入通道進(jìn)行切換，直到大部分流體依附于上部輸入通道；然后進(jìn)行中間切換，再逐漸依附于下部輸入通道，完成流體在射流短節(jié)內(nèi)的附壁切換過程。從模擬過程可以看出，射流振蕩器工作正常，流體在射流短節(jié)內(nèi)能夠順利完成附壁與切換過程。

圖3 模擬監(jiān)測(cè)界面位置示意圖Fig.3 Schematic diagram of simulation monitoring interface location

圖4 射流振蕩減摩阻工具模擬過程（入口流速為3.056 m／s）Fig.4 Simulation process of jet oscillation friction reducing tool（inlet velocity is 3.056 m／s）

通過建立監(jiān)測(cè)截面，分別得到不同初始條件下射流短節(jié)內(nèi)部流體中心截面上的速度分布云圖。通過CFD-Post得到不同入口流速條件下的x-y平面內(nèi)的速度計(jì)算結(jié)果，如圖5所示。從圖5可以看出，當(dāng)流速較低時(shí)，流體進(jìn)入射流短節(jié)后在振蕩室內(nèi)做無規(guī)律運(yùn)動(dòng)，大部分流體通過流體出口流出，較少部分流體通過上下直線流道回流至彎曲流道（圖5a）；當(dāng)流體入口流速增加到一定數(shù)值大小后，流體在振蕩室內(nèi)形成規(guī)則的旋流，在上下直線流道內(nèi)不存在回流（圖5d）；當(dāng)流體入口流速繼續(xù)增大，振蕩室內(nèi)流體的旋流狀態(tài)變紊亂，上下直線流道內(nèi)再次出現(xiàn)回流現(xiàn)象（圖5f）。

圖5 不同入口流速下的速度計(jì)算結(jié)果Fig.5 Results of velocity calculation at different inlet velocities

2.2.2 壓力場(chǎng)分析

通過速度場(chǎng)分析，完成了射流短節(jié)內(nèi)部流體的流動(dòng)狀態(tài)研究，使流體在射流短節(jié)內(nèi)的循環(huán)路線得到驗(yàn)證，為進(jìn)一步開展射流振蕩減摩阻工具內(nèi)部流場(chǎng)研究，進(jìn)行射流短節(jié)內(nèi)流體的壓力場(chǎng)分析，得到不同入口流速下的壓力場(chǎng)計(jì)算結(jié)果，如圖6所示。

圖6 不同入口流速下的壓力計(jì)算結(jié)果Fig.6 Pressure calculation results at different inlet flow rates

由圖6可以看出，流體入口位置壓力最大，隨著流體在射流短節(jié)內(nèi)循環(huán)，流體壓力逐漸降低，在流體出口位置，流體壓力最低；在上下輸入通道位置壓力分布差異比較明顯，在流體流速更大位置壓力反而相對(duì)更低，這是由于流體在射流短節(jié)內(nèi)存在附壁切換現(xiàn)象；在上下直線流道和控制端口位置，流體壓力也與流體流速分布規(guī)律相反；流體經(jīng)過入水管流道和擴(kuò)散段后，由于與分流劈之間存在水力沖擊作用，在分流中心點(diǎn)處壓力大小明顯高于附近位置的壓力。

2.2.3 分流中心點(diǎn)壓力分析

結(jié)合壓力場(chǎng)分析結(jié)果，分析分流中心點(diǎn)的壓力變化，得到不同入口流量條件下分流中心點(diǎn)的壓力變化曲線，結(jié)果如圖7所示，其中在數(shù)值計(jì)算過程中，計(jì)算時(shí)間步均為5 000步，每一步時(shí)長(zhǎng)為0.01 s。

圖7 不同入口流速條件下，分流中心點(diǎn)壓力變化曲線Fig.7 Pressure change curve of the shunt center point under different inlet velocity conditions

由圖7可以看出，隨著入口流量增大，即入口流速增加，分流中心點(diǎn)處壓力逐漸增大；隨著計(jì)算時(shí)間的增加，一段時(shí)間后模擬的6種入口流量條件下的分流中心點(diǎn)壓力均出現(xiàn)周期性壓力波動(dòng)，與模擬得到的附壁切換現(xiàn)象保持一致；不同入口流速條件下，當(dāng)流體在射流短節(jié)內(nèi)均勻射流時(shí)分流中心點(diǎn)壓力波動(dòng)范圍與壓力值相對(duì)較小，當(dāng)流體在射流短節(jié)內(nèi)產(chǎn)生附壁切換現(xiàn)象后，分流中心點(diǎn)位置的壓力開始出現(xiàn)周期性波動(dòng)，但周期性壓力波動(dòng)平穩(wěn)后，其周期均保持在10 s左右。

3 實(shí)驗(yàn)驗(yàn)證

為了分析研究仿真分析結(jié)果的準(zhǔn)確性，按照設(shè)計(jì)參數(shù)加工了一套射流振蕩減摩阻工具進(jìn)行實(shí)驗(yàn)測(cè)試，其中射流短節(jié)如圖8所示。實(shí)驗(yàn)方案如圖9所示，射流振蕩減摩阻工具以水平方式放置，采用一端固定一端軸向游動(dòng)，將減震器端使用機(jī)械機(jī)構(gòu)進(jìn)行夾持，保證左端不產(chǎn)生軸向移動(dòng)，射流振蕩器端游動(dòng)。采用雙泵系統(tǒng)對(duì)測(cè)試工具進(jìn)行供液，泵1和泵2從水箱中吸水，經(jīng)過出口總閥、入水軟管進(jìn)入射流振蕩減摩阻工具，由右端流出，經(jīng)過出水管流回水箱。實(shí)驗(yàn)過程中，所用的流體介質(zhì)均為清水，通過流量壓力監(jiān)測(cè)控制臺(tái)實(shí)現(xiàn)流量與壓力的監(jiān)測(cè)與調(diào)節(jié)。

結(jié)合實(shí)驗(yàn)測(cè)試方案進(jìn)行實(shí)驗(yàn)測(cè)試，實(shí)驗(yàn)現(xiàn)場(chǎng)如圖10所示。每次調(diào)節(jié)泵輸入流量后，待射流振蕩減摩阻工具工作平穩(wěn)后，通過數(shù)據(jù)采集儀器進(jìn)行壓力數(shù)據(jù)采集。

圖8 新型射流振蕩減摩阻工具射流短節(jié)Fig.8 Jet joint of the new jet oscillation friction reducing tool

圖9 測(cè)試設(shè)備安裝方案示意圖Fig.9 Schematic diagram of test equipment installation scheme

圖10 射流振蕩減摩阻工具實(shí)驗(yàn)現(xiàn)場(chǎng)Fig.10 Jet oscillation friction reducing tool test site

將壓力測(cè)試結(jié)果與數(shù)值模擬計(jì)算結(jié)果進(jìn)行對(duì)比，得到實(shí)驗(yàn)測(cè)試結(jié)果與仿真計(jì)算結(jié)果對(duì)比圖，如圖11所示。由圖11可以看出，實(shí)驗(yàn)測(cè)試結(jié)果與仿真計(jì)算結(jié)果得到的壓力參數(shù)與入口流速之間的變化規(guī)律基本一致。在射流振蕩減摩阻工具工作過程中，隨著入口流量的增加，射流振蕩減摩阻工具的入口壓力與出口壓力均逐漸增加，出口壓力大小增加緩慢，而入口壓力增加速度較快，測(cè)試壓降變化規(guī)律與入口壓力變化規(guī)律保持一致。但實(shí)驗(yàn)測(cè)試得到的壓力參數(shù)大小低于數(shù)值模擬得到的壓力參數(shù)，其原因在于實(shí)驗(yàn)測(cè)試過程中，壓力測(cè)試點(diǎn)為測(cè)試管線進(jìn)出口位置，而數(shù)值模擬計(jì)算結(jié)果的壓力參數(shù)對(duì)應(yīng)于射流短節(jié)流體出入口位置，并且在測(cè)試過程中，測(cè)試管線存在一定的漏失，都將導(dǎo)致測(cè)試結(jié)果低于數(shù)值模擬計(jì)算結(jié)果。

圖11 實(shí)驗(yàn)測(cè)試結(jié)果與仿真結(jié)果對(duì)比Fig.11 Comparison of experimental test results and simulation results

在其他條件不變的情況下，依次改變泵的輸入流量，測(cè)試不同流量下工具的頻振動(dòng)率與加速度，測(cè)試現(xiàn)場(chǎng)如圖12所示，測(cè)試結(jié)果見表3。從表3測(cè)試結(jié)果可以看出，隨著入口流量增大，工具的振動(dòng)頻率變快，振動(dòng)幅度變大。

圖12 射流振蕩減摩阻工具振動(dòng)測(cè)試Fig.12 Vibration test of jet oscillation friction reducing tool

表3 不同流量條件下的振動(dòng)測(cè)試結(jié)果Table3 Vibration test results under different flow conditions

4 結(jié)論

1）針對(duì)現(xiàn)有鉆井工藝中井下減摩阻工具結(jié)構(gòu)復(fù)雜、工作過程易受地層環(huán)境影響等問題，設(shè)計(jì)了一種新型射流振蕩減摩阻工具，該工具通過射流短節(jié)內(nèi)產(chǎn)生的工作壓降，形成脈沖壓力波，作用在減震器上產(chǎn)生軸向振動(dòng)，以改變鉆柱受力狀態(tài)，降低鉆井摩阻。新型射流振蕩減摩阻工具的射流振蕩器內(nèi)部無運(yùn)動(dòng)件，結(jié)構(gòu)更加簡(jiǎn)單，依靠流體在循環(huán)流道內(nèi)自動(dòng)切換產(chǎn)生射流振蕩效果，實(shí)現(xiàn)減摩阻功能。

2）數(shù)值模擬分析和實(shí)驗(yàn)研究結(jié)果表明，本文設(shè)計(jì)的新型射流振蕩減摩阻工具內(nèi)部流場(chǎng)隨著入口流量的增加，入口壓力與出口壓力均逐漸增加，出口壓力增加緩慢，入口壓力增加速度則較快；實(shí)驗(yàn)測(cè)試壓降變化規(guī)律與入口壓力變化規(guī)律保持一致，驗(yàn)證了射流振蕩減摩阻工具結(jié)構(gòu)設(shè)計(jì)的合理性和可靠性，可為實(shí)現(xiàn)射流振蕩減摩阻工具的現(xiàn)場(chǎng)實(shí)際應(yīng)用提供參考。