張辛 葛家旺 高嘉寶 孟翔

DOI：10.20031/j.cnki.0254?6094.202403017

摘 要? 針對海域天然氣水合物水射流破碎開采，對角形噴嘴、亥姆霍茲噴嘴、錐形噴嘴開展天然氣水合物儲層條件下空化效果對比，并深入研究進口壓力、出口壓力對亥姆霍茲噴嘴空化特性的影響。建立噴嘴內(nèi)流場的數(shù)學模型和物理模型，采用流體仿真軟件進行數(shù)值模擬，以氣含量作為判斷空化效果的主要依據(jù)開展研究。模擬結(jié)果表明：天然氣水合物儲層條件下，亥姆霍茲噴嘴空化性能優(yōu)于角形噴嘴和錐形噴嘴，噴嘴內(nèi)部最大氣含量達到0.95;保持其他條件不變，進口壓力越大，亥姆霍茲噴嘴產(chǎn)生的空化效果越好;天然氣水合物儲層的高壓環(huán)境會產(chǎn)生阻滯作用，降低亥姆霍茲噴嘴空化性能。

關(guān)鍵詞 天然氣水合物 空化射流 數(shù)值模擬 射流破碎

中圖分類號 TE95?? 文獻標志碼 A?? 文章編號 0254?6094（2024）03?0448?07

Numerical Simulation of the Flow Field in the Cavitation Nozzle for Marine Gas Hydrate

ZHANG Xin， GE Jia?wang， GAO Jia?bao， MENG Xiang

（College of Electromechanical Engineering， China University of Petroleum （East China））

Abstract?? Considering the water jet crushing exploitation of marine gas hydrate， the cavitation effects of angle nozzle， Helmholtz nozzle and conical nozzle under the conditions of natural gas hydrate reservoir were compared， and the influence of both inlet and outlet pressure on cavitation characteristics of Helmholtz nozzle was investigated， including having mathematical model and physical model of the flow field in the nozzle established， and the fluid simulation software adopted for the numerical simulation and the gas content based to judge cavitation effect. The simulation results show that， the cavitation performance of Helmholtz nozzle outperforms that of both angle nozzle and conical nozzle under the condition of natural gas hydrate reservoir. The maximum gas content inside the nozzle reaches 0.95. When other conditions keep unchanged， the greater inlet pressure brings about better cavitation effect produced by Helmholtz nozzle; the high?pressure gas hydrate reservoir produces blocking effect and reduces cavitation performance of Helmholtz nozzle.

Key words?? natural gas hydrate， cavitation jet， numerical simulation， jet crushing

作者簡介：張辛（1987-），副教授，從事石油機械工程、非常規(guī)能源裝備、機械電子工程的研究。

通訊作者：葛家旺（1998-），碩士研究生，從事空化射流技術(shù)的研究，1937731755@qq.com。

引用本文：張辛，葛家旺，高嘉寶，等.海域天然氣水合物空化噴嘴內(nèi)流場數(shù)值模擬[J].化工機械，2024，51（3）：448-453;476.

天然氣水合物是在低溫高壓條件下形成的籠型固態(tài)化合物，相對于常規(guī)能源，天然氣水合物分布范圍廣、儲存量大、利用率高，因此，實現(xiàn)天然氣水合物的高效安全開采至關(guān)重要[1，2]。水射流破碎開采法是近年來提出的天然氣水合物開采方法，利用噴嘴產(chǎn)生的射流打碎海底水合物儲層，不需要改變儲層的溫度和壓力條件，可以大幅提高開采儲層的安全性?？栈窃诰植康蛪簵l件下發(fā)生的液體轉(zhuǎn)變?yōu)闅怏w的現(xiàn)象，如果在水射流中引入空化作用，可以對天然氣水合物儲層產(chǎn)生射流沖擊和空泡沖蝕的雙重作用[3，4]。

為了提升空化射流的效果，前人對產(chǎn)生空化射流的關(guān)鍵部件空化噴嘴開展了大量研究。向美景等研究了亥姆霍茲上噴嘴出口結(jié)構(gòu)對空化射流自激振蕩特性的影響，結(jié)果表明，上噴嘴的出口凸出長度對自激振蕩特性產(chǎn)生顯著的影響[5]。萬春浩等研究了角形噴嘴對空化射流空化效果的影響，并基于FLUENT流體仿真軟件對角形噴嘴的內(nèi)外流場進行了數(shù)值模擬，結(jié)果表明，增大噴嘴進口壓力可以提高射流空化效果[6]。廖松等采用沖蝕靶材的方法對新型中心體空化噴嘴進行了實驗研究，結(jié)果表明，新型中心體空化噴嘴在壓力較低情況下性能較好，在相同泵壓下，沖蝕性能優(yōu)于角形噴嘴[7]。張逸群等研究了旋轉(zhuǎn)空化射流噴嘴對射流沖擊特性的影響，結(jié)果表明，旋轉(zhuǎn)空化射流的負壓分布范圍大，具有更強的空蝕效應(yīng)[8]。陳松平等研究了縮放型噴嘴結(jié)構(gòu)參數(shù)對射流流場的影響，結(jié)果表明，噴嘴擴張角對空化射流影響較大，射流剪切層內(nèi)產(chǎn)生許多漩渦結(jié)構(gòu)[9]。上述研究分析了不同結(jié)構(gòu)空化噴嘴的空化效果及其影響因素，豐富了空化射流相關(guān)理論，但現(xiàn)有針對空化噴嘴的研究大多局限于常壓條件，針對天然氣水合物儲層的高圍壓環(huán)境，需要進一步探究不同類型空化噴嘴的性能優(yōu)劣以及空化性能影響因素。

筆者利用計算流體力學對角形噴嘴、亥姆霍茲噴嘴、錐形噴嘴的內(nèi)流場進行數(shù)值模擬，以氣含量為判斷空化效果的主要依據(jù)，探究3種空化噴嘴在天然氣水合物儲層條件下的空化效果強弱，針對亥姆霍茲噴嘴，進一步研究了進口壓力、出口壓力對空化射流效果的影響，為空化射流在天然氣水合物開采領(lǐng)域的應(yīng)用提供參考。

1 物理模型

1.1 幾何模型

基于前人對空化噴嘴的研究，建立3種空化噴嘴內(nèi)腔結(jié)構(gòu)（圖1），對于角形噴嘴，入口段長7 mm，入口段收縮角13.5°，中間段長3 mm，中間段直徑2 mm，出口段長5 mm，出口段擴散角60°;對于亥姆霍茲噴嘴，入口段長7 mm，入口段直徑3.4 mm，中間段長3 mm，中間段直徑12 mm，中間段傾斜壁夾角120°，出口段長5 mm，出口段直徑4 mm;對于錐形噴嘴，入口段長7 mm，入口段收縮角25°，出口段長8 mm，出口段直徑2 mm。

1.2 邊界條件及網(wǎng)格劃分

基于噴嘴的對稱性，為了節(jié)約計算資源，采用二維軸對稱模型進行流場計算[10]。設(shè)置噴嘴入口邊界條件為壓力入口，設(shè)置噴嘴出口邊界條件為壓力出口，為保證流體在噴嘴中正常流動，需保證進口壓力高于出口壓力。湍流參數(shù)選擇湍流強度與水力直徑，進出口湍流參數(shù)均設(shè)置為5%。固體壁面條件滿足無滑移壁面條件，進壁面區(qū)域采用壁面函數(shù)處理[11，12]。流場迭代求解方法選擇Coulped，離散格式選擇QUICK方程。

用DesignModeler軟件進行建模后，利用Meshing對整個噴嘴進行網(wǎng)格劃分。精密的網(wǎng)格可以保證模擬流動的真實性，因此，在噴嘴壁邊界層區(qū)域進行局部加密，充分保證網(wǎng)格的質(zhì)量，最終網(wǎng)格的平均質(zhì)量大于0.9。

2 數(shù)學模型

空化射流涉及液體復雜的相變過程，因此，空化射流屬于湍流流動，筆者選取混合模型和RNG k湍流模型進行計算。

2.1 混合模型

混合模型允許不同相之間的滲透，計算精度高，不容易出現(xiàn)發(fā)散現(xiàn)象，可以同時和空化模型使用求解[13，14]。混合模型的連續(xù)方程為：

（ρ）+▽·（ρ）=0???? （1）

式中 ——質(zhì)量平均密度;

ρ——混合物的密度。

2.2 湍流模型

RNG k模型是一種兩方程湍流模型，不僅在漩渦捕捉方面有很大的優(yōu)勢，而且對雷諾數(shù)適應(yīng)性強，具有很高的湍流計算精度，適合復雜剪切流，在空化噴嘴的流場模擬中可以產(chǎn)生很好的模擬效果[15]。其k方程和ε方程分別為：

k方程+=（αμ）+G+G-ρε-Y （2）

ε方程+=（αμ）+C×（G+CG）-

Cρ-R? （3）

式中 C、C、C——經(jīng)驗常數(shù);

G——浮力影響產(chǎn)生的湍流動能;

G——平均速度梯度產(chǎn)生的湍流????????? 動能;

k——湍流動能;

R——空泡半徑;

Y——總耗散率;

α、α——普朗特常數(shù);

ε——湍流動能耗散率;

μ——有效黏度。

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

3.1 不同噴嘴空化效果對比

當流域局部壓力低于當前環(huán)境溫度下流體飽和蒸氣壓時，流體會發(fā)生空化現(xiàn)象。天然氣水合物在低溫高壓環(huán)境下生成，筆者以水為研究對象，設(shè)置水的飽和蒸氣壓（5 ℃時）為900 Pa，模擬天然氣水合物儲層低溫環(huán)境，設(shè)置噴嘴出口壓力5 MPa，模擬天然氣水合物儲層的高壓條件，設(shè)置噴嘴進口壓力40 MPa，保證流體的正常流動。對于工程應(yīng)用水射流，噴嘴出口速度理論計算公式如下：

v=44.77?? （4）

式中 p——射流壓力，MPa;

v——射流速度，m/s。

對3種不同結(jié)構(gòu)的空化噴嘴完成空化射流數(shù)值仿真，提取各噴嘴出口軸線處的速度值，表1為噴嘴出口軸線處速度仿真計算值與理論計算值的對比，由表可見：3種不同結(jié)構(gòu)的空化噴嘴出口軸線處速度仿真計算值與理論計算值基本吻合，仿真誤差均在1%以內(nèi)，證明采用的仿真計算方法合理[16]。

表1 空化噴嘴出口軸線處速度值

[噴嘴類型????? 出口軸線處速度/m·s-1????? 誤差/%?? 理論計算值??? 仿真計算值??? 角形噴嘴?????? 282.9????? 283.1????? 0.07 亥姆霍茲噴嘴??? 282.9????? 283.0????? 0.04 錐形噴嘴?????? 282.9????? 283.2????? 0.10 ]

圖2為3種不同結(jié)構(gòu)空化噴嘴的空泡分布云圖。由圖可見：3種空化噴嘴內(nèi)部均有空泡析出，其中，錐形噴嘴產(chǎn)生的空泡數(shù)量最少，僅在出口段壁面邊界層附近出現(xiàn)少量空泡，相比之下，角形噴嘴產(chǎn)生的空泡數(shù)量優(yōu)于錐形噴嘴，這是由于角形噴嘴出口段存在擴散角，可以為空泡提供足夠的發(fā)展空間。亥姆霍茲噴嘴產(chǎn)生的空泡數(shù)量最多，大量空泡集中在噴嘴中間段，且空泡對稱分布呈渦旋狀，造成這種現(xiàn)象的主要原因是亥姆霍茲噴嘴中間段存在傾斜壁結(jié)構(gòu)，在傾斜壁拐角處易形成漩渦低壓區(qū)，促進空化泡的產(chǎn)生[17，18]。

圖3為3種不同結(jié)構(gòu)空化噴嘴內(nèi)部最大氣含量，由圖可見：不同結(jié)構(gòu)空化噴嘴內(nèi)部最大氣含量有明顯差異，亥姆霍茲噴嘴內(nèi)部最大氣含量達到0.95，角形噴嘴內(nèi)部最大氣含量達到0.75，相比之下，錐形噴嘴內(nèi)部最大氣含量僅為0.55，氣含量越大，空化發(fā)生程度越劇烈，由此可見，空泡在亥姆霍茲噴嘴內(nèi)部發(fā)展狀態(tài)最佳，可以產(chǎn)生更強的空蝕作用，有利于對水合物沉積物的沖擊破碎[19，20]。綜合噴嘴空泡分布云圖可以看出，3種不同結(jié)構(gòu)空化噴嘴中亥姆霍茲噴嘴在天然氣水合物儲層條件下空化效果最好。

3.2 進口壓力對亥姆霍茲空化性能影響

空化噴嘴的工作性能受到多重因素的影響，當噴嘴結(jié)構(gòu)和出口壓力維持不變時，進口壓力對空化噴嘴性能產(chǎn)生很大的影響。保持出口壓力為5 MPa、使進口壓力分別為30、40、50、60 MPa，分析進口壓力對亥姆霍茲噴嘴氣含量的影響情況。由圖4可見：在不同進口壓力作用下，亥姆霍茲噴嘴中間段均產(chǎn)生聚集性空泡，亥姆霍茲噴嘴出口段在壁面邊界層附近產(chǎn)生空泡，空泡在噴嘴內(nèi)部呈對稱分布，當噴嘴進口壓力為30 MPa時，亥姆霍茲噴嘴出口段僅有少量空泡，空泡在噴嘴出口段長度占比不足二分之一，當噴嘴進口壓力為40 MPa時，亥姆霍茲噴嘴出口段空泡數(shù)量增多，空泡在噴嘴出口段長度占比超過二分之一，當進口壓力繼續(xù)上升，空泡發(fā)展至亥姆霍茲噴嘴出口端，空泡在噴嘴出口段長度占比達到百分之百。

圖5為不同進口壓力作用下亥姆霍茲噴嘴內(nèi)部最大氣含量，由圖可見：亥姆霍茲噴嘴內(nèi)部最大氣含量隨進口壓力的增加而上升，綜合不同進口壓力下亥姆霍茲噴嘴的空泡分布云圖，可以得出，進口壓力越大，亥姆霍茲噴嘴空化效果越好。這是由于進口壓力提高，射流剪切層作用變強，容易形成大尺度的漩渦。因此，為了保證天然氣水合物儲層條件下亥姆霍茲噴嘴的空化效果，可以適當提高噴嘴進口壓力，但過高的噴嘴進口壓力會造成對泵壓需求過大的情況，在實際應(yīng)用過程中，需要考慮噴嘴進口壓力與泵壓的平衡。

3.3 出口壓力對亥姆霍茲空化性能影響

天然氣水合物儲層的高壓環(huán)境會對噴嘴空化性能產(chǎn)生很大的影響。因此，保持進口壓力為50 MPa，使出口壓力分別為0、5、10、15 MPa，分析出口壓力對亥姆霍茲噴嘴氣含量的影響情況。由圖6可見：當出口壓力為0 MPa時，亥姆霍茲噴嘴內(nèi)部產(chǎn)生大量空泡，空泡幾乎充滿噴嘴中間段，隨著出口壓力的不斷提高，亥姆霍茲噴嘴內(nèi)部空泡數(shù)量逐漸減少，當出口壓力為15 MPa時，亥姆霍茲噴嘴內(nèi)部幾乎沒有空泡產(chǎn)生。

圖7所示為不同出口壓力作用下亥姆霍茲噴嘴內(nèi)部最大氣含量，由圖可見：隨著出口壓力的不斷增大，亥姆霍茲噴嘴內(nèi)部最大氣含量不斷降低。當出口壓力為0 MPa時，噴嘴內(nèi)部最大氣含量達到0.9以上，當出口壓力為5 MPa時，噴嘴內(nèi)部最大氣含量略有下降，但仍達到0.9以上，當出口壓力繼續(xù)提高，噴嘴內(nèi)部最大氣含量下降率增大，出口壓力為15 MPa時，噴嘴內(nèi)部最大氣含量已經(jīng)低于0.6。綜合不同出口壓力下亥姆霍茲噴嘴的空泡分布云圖，可以得出，較大的出口壓力會對亥姆霍茲噴嘴空化產(chǎn)生抑制作用。這主要是由于出口的高壓力會對射流速度和空泡發(fā)展產(chǎn)生阻滯作用。

4 結(jié)論

4.1 天然氣水合物儲層條件下，亥姆霍茲噴嘴空化性能優(yōu)于角形噴嘴和錐形噴嘴。亥姆霍茲噴嘴內(nèi)部產(chǎn)生聚集性空泡，空泡對稱分布呈渦旋狀，空泡潰滅產(chǎn)生壓力波，有利于水合物沉積物的沖蝕破壞。

4.2 當噴嘴結(jié)構(gòu)和出口壓力維持不變時，進口壓力對空化噴嘴性能產(chǎn)生很大的影響。進口壓力越大，亥姆霍茲噴嘴產(chǎn)生的空化效果越好。因此，在實際應(yīng)用中，根據(jù)高壓泵的工作范圍，可以適當提高噴嘴進口壓力。

4.3 天然氣水合物儲層的高壓環(huán)境會降低亥姆霍茲噴嘴空化性能。當出口壓力為0 MPa時，亥姆霍茲噴嘴內(nèi)部最大氣含量達0.9以上，當出口壓力為15 MPa時，亥姆霍茲噴嘴內(nèi)部最大氣含量低于0.6。出口壓力越小，亥姆霍茲噴嘴產(chǎn)生的空化效果越好。

參 考 文 獻

[1] 李根生，田守嶒，張逸群.空化射流鉆徑向井開采天然氣水合物關(guān)鍵技術(shù)研究進展[J].石油科學通報，2020，5（3）：349-365.

[2] 王旱祥，王姣姣，于長錄，等.天然氣水合物開采關(guān)鍵技術(shù)研究現(xiàn)狀[J].鉆采工藝，2022，45（5）：51-56.

[3] 楊林.海底水合物儲層在高壓水射流作用下的破碎過程及儲層改造增產(chǎn)研究[D].長春：吉林大學，2018.

[4] 王國榮，鐘林，周守為，等.天然氣水合物射流破碎工具及其配套工藝技術(shù)[J].天然氣工業(yè)，2017，37（12）：68-74.

[5] 向美景，王曉川，李登，等.亥姆霍茲上噴嘴出口結(jié)構(gòu)對射流振蕩特性的影響[J].振動與沖擊，2020，39（7）：74-80.

[6] 萬春浩，王瑞和，周衛(wèi)東，等.角形空化噴嘴內(nèi)外流場的數(shù)值模擬[J].斷塊油氣田，2019，26（2）：257-260.

[7] 廖松，鄧松圣，趙華忠，等.中心體空化噴嘴射流實驗研究[J].天然氣與石油，2021，39（2）：100-106.

[8] 張逸群，武曉亞，李根生，等.天然氣水合物旋轉(zhuǎn)空化射流沖蝕性能研究[J].中南大學學報（自然科學版），2022，53（3）：909-923.

[9] 陳松平，李偉，張文全，等.縮放型噴嘴空化流動特性的數(shù)值分析[J].排灌機械工程學報，2023，41（3）：1-7.

[10] 董景明，李帥，鐘蕭，等.兩段式噴嘴淹沒空化射流流場數(shù)值模擬[J].流體機械，2022，50（11）：87-95.

[11] 戚美，王立夫，陳慶光，等.塤型結(jié)構(gòu)噴嘴自激振蕩脈沖空化射流的數(shù)值模擬[J].山東科技大學學報（自然科學版），2018，37（6）：108-116.

[12] WU X Y，ZHANG Y Q，TAN Y W，et al.Flow?visua? lization and numerical investigation on the optimum design of cavitating jet nozzle[J].Petroleum Science，2022，19（5）：2284-2296.

[13] 張修占，侯文江，劉建成，等.用于船體表面清洗的風琴管噴嘴沖擊特性[J].船舶工程，2022，44（10）：142-148.

[14] YANG Y F，WANG G W，SHI W D，et al.Turbulence Characteristics in the Mixing Layer of a Submerged Cavitating Jet at High Reynolds Numbers[J].Sustainability，2022，14（19）：11963.

[15] 曹澤平，劉桓龍，謝遲新.新型環(huán)流式空化噴嘴的空化及清洗特性[J].液壓與氣動，2020（3）：136-144.

[16] 李鋒.空化水射流船舶清洗技術(shù)的研究[D].哈爾濱：哈爾濱工程大學，2009.

[17] 戚美，王立夫，趙見龍，等.雙腔室自激振蕩脈沖噴嘴空化射流外部流場的數(shù)值模擬[J].機械科學與技術(shù)，2019，38（4）：515-521.

[18] 方珍龍.亥姆赫茲式自激振蕩射流空化特性研究[D].武漢：武漢大學，2016.

[19] 王國榮，黃蓉，鐘林，等.固態(tài)流化采掘海洋天然氣水合物藏破碎參數(shù)的優(yōu)化設(shè)計[J].天然氣工業(yè)，2018，38（10）：84-89.

[20] 涂翊翔，樊辰星，錢亦楠，等.雙空化射流中心體結(jié)構(gòu)優(yōu)化及破煤巖特性[J].煤炭學報，2022，47（9）：3250-3259.

（收稿日期：2023-04-17，修回日期：2024-05-13）