張 辛 葛家旺 高嘉寶 孟 翔

(中國石油大學(xué)(華東)機(jī)電工程學(xué)院)

0 引 言

我國南海神狐海域天然氣水合物沉積層具有弱膠結(jié)、滲透性差及非成巖等特點(diǎn)，采用現(xiàn)有的開采方法很難實(shí)現(xiàn)效率與經(jīng)濟(jì)的均衡，需要進(jìn)行儲(chǔ)層改造以提升天然氣水合物低滲透儲(chǔ)層的生產(chǎn)性能[1-3]。空化射流是近年來發(fā)展起來的新型射流，目前已經(jīng)普遍應(yīng)用于石油領(lǐng)域。當(dāng)流場(chǎng)局部壓力低于當(dāng)前溫度下液體飽和蒸汽壓時(shí)，液體中產(chǎn)生空泡的過程被稱為空化?？栈瘹馀莅l(fā)生潰滅，會(huì)產(chǎn)生高頻的壓力波，天然氣水合物儲(chǔ)層吸收空化產(chǎn)生的壓力波能量，從而產(chǎn)生微裂紋[4-7]。因此，空化射流技術(shù)作為儲(chǔ)層改造方式，可以提高天然氣水合物儲(chǔ)層的局部滲透率和孔隙度，提高開采效率。

空化射流主要依靠空化噴嘴產(chǎn)生，國內(nèi)外常用的一種噴嘴結(jié)構(gòu)是風(fēng)琴管噴嘴，大量學(xué)者對(duì)這種結(jié)構(gòu)噴嘴進(jìn)行了仿真研究與試驗(yàn)驗(yàn)證。趙新澤等[8]對(duì)風(fēng)琴管噴嘴的結(jié)構(gòu)進(jìn)行改進(jìn)，設(shè)計(jì)了圓弧段擴(kuò)張口徑，并采用流體仿真軟件對(duì)噴嘴進(jìn)行流場(chǎng)分析，結(jié)果表明，風(fēng)琴管噴嘴空化效果顯著。韓健等[9]將打擊試驗(yàn)與噪聲分析方法相結(jié)合，分析了風(fēng)琴管噴嘴的空化特性，結(jié)果表明，風(fēng)琴管噴嘴能產(chǎn)生高頻壓力振蕩并促進(jìn)空化射流的產(chǎn)生。孫軍等[10]研究風(fēng)琴管噴嘴出口段內(nèi)表面結(jié)構(gòu)對(duì)空化特性的影響，對(duì)管道內(nèi)表面凸起結(jié)構(gòu)完成流場(chǎng)壓力分析，結(jié)果表明，凸起的存在對(duì)近壁面流場(chǎng)壓力有顯著影響，凸起達(dá)到一定數(shù)量時(shí)，壓力波動(dòng)呈現(xiàn)出類似交變變化。張修占等[11]采用數(shù)值模擬與試驗(yàn)研究相結(jié)合的方法，研究風(fēng)琴管噴嘴射流沖擊特性，結(jié)果表明，合理選擇射流靶距使空泡充分發(fā)展，可以有效增強(qiáng)風(fēng)琴管噴嘴的沖擊性能。于海濤等[12]利用Fluent軟件對(duì)風(fēng)琴管噴嘴空化特性展開研究，結(jié)果表明，適當(dāng)增加風(fēng)琴管噴嘴擴(kuò)散段長度可以增強(qiáng)空化作用。

以上研究表明，風(fēng)琴管噴嘴具有優(yōu)異的空化效果。然而目前對(duì)空化射流的研究主要應(yīng)用于高壓清洗和石油鉆采方面，在天然氣水合物開采領(lǐng)域還鮮有提及，針對(duì)空化射流在天然氣水合物儲(chǔ)層高圍壓條件下的應(yīng)用，需要進(jìn)一步研究風(fēng)琴管噴嘴的空化沖擊性能。為此，筆者建立了風(fēng)琴管噴嘴模型，通過Fluent流體仿真軟件開展風(fēng)琴管噴嘴空化的模擬研究，根據(jù)風(fēng)琴管噴嘴的速度分布、壓力分布及氣相體積分?jǐn)?shù)分布等流場(chǎng)特征，研究圍壓條件下風(fēng)琴管噴嘴的空化性能；搭建空化射流沖擊套管試驗(yàn)裝置，分析圍壓條件下噴嘴內(nèi)腔結(jié)構(gòu)參數(shù)對(duì)風(fēng)琴管噴嘴沖擊效果的影響規(guī)律，以期為天然氣水合物儲(chǔ)層條件下空化射流的應(yīng)用提供理論依據(jù)。

1 數(shù)值模型建立

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

空化射流涉及液體相變過程，因此空化射流屬于氣液兩相流，同時(shí)也屬于湍流流動(dòng)，本文選取混合模型和RNGk-ε湍流模型進(jìn)行計(jì)算?；旌夏Ｐ陀?jì)算精度高，收斂速度快，符合空化特性的計(jì)算要求[13-14]。RNGk-ε湍流模型不但在漩渦計(jì)算方面精度高，而且可以很好地模擬空化噴嘴近壁區(qū)域的空化現(xiàn)象[15-17]，所以選擇RNGk-ε湍流模型完成空化現(xiàn)象的仿真運(yùn)算?；旌夏Ｐ偷倪B續(xù)方程為：

(1)

1.2 噴嘴幾何參數(shù)

依據(jù)瞬態(tài)流理論對(duì)風(fēng)琴管噴嘴進(jìn)行結(jié)構(gòu)設(shè)計(jì)，噴嘴主要結(jié)構(gòu)尺寸如圖1所示，入口段直徑D1=26 mm，入口段長度L1=8 mm；諧振腔直徑D2=8 mm，諧振腔長度L2=21 mm；圓柱段直徑D3=4 mm，圓柱段長度L3=7 mm；擴(kuò)散段長度L4=14 mm，擴(kuò)散段角度α=60°。流體在入口段流入，經(jīng)過諧振腔時(shí)產(chǎn)生共振，促進(jìn)大結(jié)構(gòu)分離環(huán)狀渦流的產(chǎn)生，增強(qiáng)空化作用。

圖1 噴嘴主要結(jié)構(gòu)尺寸示意圖Fig.1 Main structure and dimension of nozzle

1.3 邊界條件

圖2 噴嘴計(jì)算域Fig.2 Calculation domain of nozzle

設(shè)置噴嘴出口為壓力出口，研究圍壓作用下流體的空化特性。為保證流體能夠從出口正常流出，需要保證進(jìn)出口存在一定的壓力差，設(shè)置噴嘴入口為壓力入口。進(jìn)出口的射流湍流強(qiáng)度均設(shè)為5%。選擇壓力-速度耦合求解算法為SIMPLEC，采用二階迎風(fēng)離散格式，壓力方程采用PRESTO！，近壁面區(qū)域采用壁面函數(shù)法處理[19-20]。

1.4 網(wǎng)格劃分

采用Meshing對(duì)計(jì)算域進(jìn)行網(wǎng)格劃分，生成面網(wǎng)格最小尺寸為0.1 mm。為了提高計(jì)算效率，對(duì)流體重要流動(dòng)區(qū)域的網(wǎng)格局部加密[21]。網(wǎng)格單元質(zhì)量數(shù)值范圍為0～1，數(shù)值越接近于1，網(wǎng)格單元質(zhì)量越高。通過網(wǎng)格無關(guān)性分析，計(jì)算域劃分網(wǎng)格單元質(zhì)量均大于0.9，網(wǎng)格單元質(zhì)量?jī)?yōu)異，滿足運(yùn)算要求。

2 數(shù)值模擬結(jié)果與分析

2.1 風(fēng)琴管噴嘴空化特性

圖3為入口壓力30 MPa、圍壓5 MPa時(shí)風(fēng)琴管噴嘴速度分布云圖。由圖3可見，流體速度在噴嘴的圓柱段達(dá)到較高水平，后經(jīng)噴嘴擴(kuò)散段形成高速射流。在流體噴出圓柱段的初始階段，存在噴嘴軸線速度保持不變的核心區(qū)，隨著高速射流的逐漸擴(kuò)張，射流能量逐漸衰減，等速核心區(qū)消失。等速核心區(qū)內(nèi)保持較大的軸向速度，有利于射流完成沖擊破碎[22]；同時(shí)根據(jù)伯努利方程，具有高流速的等速核心區(qū)內(nèi)壓力較低，有利于空化泡的產(chǎn)生。

圖3 風(fēng)琴管噴嘴速度分布云圖Fig.3 Nephogram for velocity distribution of organ-pipe nozzle

圖3中局部放大部分為外流場(chǎng)初始階段流體速度矢量。由圖3可見，在擴(kuò)張段出口處存在回流現(xiàn)象，這是由于高速射流與外流場(chǎng)靜止流體發(fā)生剪切作用，從而在流體作用區(qū)域形成漩渦。漩渦中心的壓力低于周圍壓力，在促進(jìn)空化泡生長的同時(shí)，有利于空化泡在外流場(chǎng)的運(yùn)輸，使空化泡在更遠(yuǎn)距離潰滅，增強(qiáng)空化泡的沖蝕作用。

空化現(xiàn)象產(chǎn)生的前提是存在局部低壓區(qū)，因此對(duì)圍壓條件下風(fēng)琴管噴嘴形成的壓力場(chǎng)進(jìn)行分析。圖4為入口壓力30 MPa、圍壓5 MPa時(shí)風(fēng)琴管噴嘴壓力分布云圖。由圖4可見，在風(fēng)琴管噴嘴內(nèi)部存在明顯的低壓區(qū)域，且均集中在圓柱段，圓柱段大部分區(qū)域的壓力都處于負(fù)壓狀態(tài)，形成“負(fù)壓區(qū)”。圓柱段內(nèi)“負(fù)壓區(qū)”對(duì)于空化非常關(guān)鍵，“負(fù)壓區(qū)”內(nèi)流體局部絕對(duì)壓力小于當(dāng)前溫度下飽和蒸汽壓，流體中的氣核在“負(fù)壓區(qū)”迅速膨脹，在流體內(nèi)部形成含有水蒸氣的空化泡[23]。產(chǎn)生“負(fù)壓區(qū)”是由于圓柱段的突縮結(jié)構(gòu)，使圓柱段流速迅速升高，從而產(chǎn)生壓降。

圖4 風(fēng)琴管噴嘴壓力云圖Fig.4 Nephogram for pressure of organ-pipe nozzle

2.2 不同圍壓對(duì)空化特性的影響

由壓力分布云圖可知，風(fēng)琴管噴嘴內(nèi)部負(fù)壓值是空化產(chǎn)生的關(guān)鍵因素，所以保持入口壓力30 MPa，使圍壓分別為0、5、10及15 MPa，分析圍壓對(duì)風(fēng)琴管噴嘴內(nèi)部負(fù)壓值的影響。圖5為不同圍壓下風(fēng)琴管噴嘴內(nèi)部產(chǎn)生的最小負(fù)壓值。由圖5可見：隨著圍壓的增大，風(fēng)琴管噴嘴內(nèi)部產(chǎn)生的最小負(fù)壓值逐漸上升，且圍壓達(dá)到的水平越高，最小負(fù)壓值隨圍壓變化的范圍越明顯；當(dāng)圍壓從0提高至5 MPa時(shí)最小負(fù)壓值僅上升0.04 MPa，而圍壓從10 MPa提高至15 MPa時(shí)最小負(fù)壓值上升0.38 MPa，“負(fù)壓區(qū)”中負(fù)壓值越小，流體中越容易形成大尺寸的空泡。因此，在一定的圍壓范圍內(nèi)，風(fēng)琴管噴嘴內(nèi)部能產(chǎn)生負(fù)壓，從而實(shí)現(xiàn)空化泡的初生；但是高圍壓條件在一定程度上會(huì)抑制風(fēng)琴管噴嘴空化“負(fù)壓區(qū)”的極值。

圖5 風(fēng)琴管噴嘴內(nèi)部最小負(fù)壓值Fig.5 Minimum negative pressure inside organ-pipe nozzle

空化現(xiàn)象一旦發(fā)生，流體就會(huì)變成氣液兩相流，所以通過氣相云圖可以觀察空化區(qū)域的大小。圖6為同一時(shí)刻風(fēng)琴管噴嘴在不同圍壓下的氣相體積分?jǐn)?shù)分布云圖。云圖顏色越深，表明這個(gè)區(qū)域的含氣體積分?jǐn)?shù)越大[24]。由圖6可見：在不同圍壓作用下，內(nèi)流場(chǎng)和外流場(chǎng)均產(chǎn)生大量氣體，產(chǎn)生的氣體沿軸心線呈對(duì)稱分布；隨著圍壓的提升，同一時(shí)刻外流場(chǎng)產(chǎn)生的氣體區(qū)域面積逐漸減小，這意味著隨水射流到達(dá)沖擊目標(biāo)的空泡數(shù)目減少，空泡潰滅產(chǎn)生的能量衰減。造成這種現(xiàn)象主要有以下2個(gè)原因：一是噴嘴出口高圍壓的擾動(dòng)，空化泡在風(fēng)琴管噴嘴內(nèi)部的初生受到限制，從而在風(fēng)管琴噴嘴內(nèi)部產(chǎn)生的空泡數(shù)量較少，隨水射流到達(dá)外流場(chǎng)的空泡數(shù)目便相對(duì)減少；二是外部環(huán)境較高的壓力，使空化泡不能充分流動(dòng)發(fā)展，大量空化泡在外流場(chǎng)的初始階段便因?yàn)楦邏鹤饔枚鴿纭Ｒ虼?，針?duì)高圍壓條件下的空化射流，可以合理控制射流靶距。

圖6 同一時(shí)刻不同圍壓下噴嘴氣相體積分?jǐn)?shù)分布云圖Fig.6 Nephogram for gas phase volume fraction distribution of nozzle under different confining pressures at the same time

2.3 不同諧振腔直徑對(duì)空化特性的影響

風(fēng)琴管噴嘴的諧振腔對(duì)于空化泡初生起關(guān)鍵性作用。保持入口壓力30 MPa、圍壓5 MPa，使諧振腔直徑分別為8、10、12和14 mm，分析諧振腔直徑對(duì)風(fēng)琴管噴嘴氣相分布的影響。圖7為不同諧振腔直徑噴嘴在外流場(chǎng)產(chǎn)生空泡的最大軸向長度和最大徑向長度。由圖7可見，在不同時(shí)刻，噴嘴在外流場(chǎng)產(chǎn)生的空泡數(shù)量明顯不同，T=2 ms時(shí)外流場(chǎng)中空泡的最大軸向長度和最大徑向長度均大于T=1 ms時(shí)。這主要是由于空化泡的形成、發(fā)展、破裂需要一定的時(shí)間。隨著諧振腔直徑不斷變大，同一時(shí)刻外流場(chǎng)產(chǎn)生的空泡數(shù)量呈現(xiàn)先增大后減小的趨勢(shì)，當(dāng)諧振腔直徑為10 mm時(shí)，外流場(chǎng)中產(chǎn)生的空泡數(shù)量最多。由此可見，諧振腔直徑在一定范圍內(nèi)存在最優(yōu)值，此時(shí)風(fēng)琴管噴嘴有更優(yōu)異的空化性能。

圖7 外流場(chǎng)空泡長度Fig.7 Cavitation bubble length in external flow field

3 空化射流沖擊套管試驗(yàn)

為了驗(yàn)證圍壓下風(fēng)琴管噴嘴的空化沖擊效果，以套管為沖擊靶物搭建試驗(yàn)裝置，對(duì)不同諧振腔直徑的風(fēng)琴管噴嘴開展沖擊套管試驗(yàn)。

3.1 試驗(yàn)裝置

圖8為圍壓條件下風(fēng)琴管噴嘴沖擊套管試驗(yàn)裝置流程圖。由圖8可見，混砂車完成磨料和清水的混合，泵車加壓將射流介質(zhì)送至高壓管線，最后進(jìn)入空化噴嘴完成射流切割內(nèi)套管。多組噴嘴均勻地布置在噴槍上，同時(shí)射流沖擊內(nèi)套管。試壓泵的主要作用是將液體注入內(nèi)外套管形成的環(huán)空間隙，為內(nèi)套管提供圍壓，其上安裝有壓力表，可以觀察并確保圍壓準(zhǔn)確；內(nèi)套管如果被射流沖穿，壓力表數(shù)值變動(dòng)；快速接頭與排氣閥的主要作用是排除內(nèi)外套管形成的環(huán)空間隙中的氣體。內(nèi)套管采用N80石油鋼管，壁厚14.27 mm。

圖8 試驗(yàn)裝置流程圖Fig.8 Flow chart of test device

3.2 試驗(yàn)結(jié)果與分析

保持試驗(yàn)條件不變，以相同時(shí)間內(nèi)射流對(duì)內(nèi)套管的沖擊坑深度為判斷依據(jù)，對(duì)不同諧振腔直徑的風(fēng)琴管噴嘴沖擊性能進(jìn)行評(píng)估，諧振腔直徑分別為8、10、12和14 mm。

圖9為相同時(shí)間內(nèi)不同諧振腔直徑風(fēng)琴管噴嘴對(duì)套管的沖擊形貌。

圖9 不同沖擊坑形貌Fig.9 Morphologies of different impact pits

由圖9可見，不同諧振腔直徑的風(fēng)琴管噴嘴產(chǎn)生的沖擊坑大小明顯不同，諧振腔直徑為10 mm的噴嘴形成的沖擊坑最大，諧振腔直徑為8 mm的噴嘴形成的沖擊坑最小。由此得出，諧振腔直徑會(huì)對(duì)風(fēng)琴管噴嘴的空化沖擊效果產(chǎn)生顯著影響。

對(duì)不同噴嘴產(chǎn)生沖擊坑的深度進(jìn)行測(cè)量，圖10為不同諧振腔直徑風(fēng)琴管噴嘴產(chǎn)生的沖擊坑深度變化曲線。

圖10 不同沖擊坑深度變化曲線Fig.10 Variation of different impact pit depths

由圖10可見：諧振腔直徑8 mm的風(fēng)琴管噴嘴產(chǎn)生的沖擊坑深度最小，沖擊坑的深度僅為7.4 mm；諧振腔直徑10 mm的風(fēng)琴管噴嘴產(chǎn)生的沖擊坑深度最大，沖擊坑的深度達(dá)到11.0 mm。相比之下，兩者形成的沖擊坑深度存在很大差值，相對(duì)于諧振腔直徑8 mm噴嘴，諧振腔直徑10 mm噴嘴產(chǎn)生的沖擊坑深度增加了48.6%。這與數(shù)值模擬結(jié)果一致，驗(yàn)證了仿真結(jié)果的準(zhǔn)確性。

4 結(jié) 論

(1)圍壓條件下風(fēng)琴管噴嘴在內(nèi)流場(chǎng)產(chǎn)生明顯的“負(fù)壓區(qū)”，有利于空化泡的初生，在外流場(chǎng)產(chǎn)生漩渦低壓區(qū)，可以促進(jìn)空化泡的生長擴(kuò)散，從而產(chǎn)生更強(qiáng)的射流沖擊力。

(2)高圍壓在一定程度上會(huì)抑制風(fēng)琴管噴嘴空化“負(fù)壓區(qū)”的極值，保持其他條件不變的情況下，圍壓值越大，外流場(chǎng)中產(chǎn)生的空泡數(shù)量越少。由于外部環(huán)境較高的壓力，存在空化泡在外流場(chǎng)的初始階段便因?yàn)楦邏鹤饔枚鴿?，所以針?duì)高圍壓條件下的空化射流，可以合理控制射流靶距，減少未到達(dá)沖擊目標(biāo)而破裂的空泡數(shù)目，形成沖蝕能力更強(qiáng)的空化射流。

(3)風(fēng)琴管噴嘴的諧振腔直徑明顯影響噴嘴的射流沖擊能力。在一定范圍內(nèi)，隨著諧振腔直徑變大，同一時(shí)刻外流場(chǎng)產(chǎn)生的空泡數(shù)量呈現(xiàn)先增大后減小的趨勢(shì)。與諧振腔直徑8 mm的噴嘴相比，相同時(shí)間內(nèi)諧振腔直徑為10 mm的噴嘴對(duì)套管產(chǎn)生的沖擊坑深度增加48.6%。