黃 飛，焦楊洋，米建宇，李樹清，趙志旗，劉 勇

(1.湖南科技大學 南方煤礦瓦斯與頂板災害預防控制安全生產(chǎn)重點實驗室，湖南 湘潭 411201；2.湖南科技大學 資源環(huán)境與安全工程學院，湖南 湘潭 411201；3.河南理工大學 河南省瓦斯地質與瓦斯治理重點實驗室——省部共建國家重點實驗室培育基地，河南 焦作 454003)

水射流割縫技術作為煤層增透抽采瓦斯的有效途徑之一，為預防煤與瓦斯突出、高效抽采煤層氣做出了卓越貢獻。其中，磨料水射流通過磨料的高速磨削作用，大幅提升了射流破碎煤巖的性能，在煤層割縫領域受到了廣泛的青睞。然而，通過大量工程實踐發(fā)現(xiàn)，傳統(tǒng)的高壓磨料水射流煤層割縫技術由于大量磨料顆粒(諸如石榴石、陶粒、石英等)的摻入帶來了一系列的問題：① 殘留磨料顆粒廣泛分布在割縫后的縫槽與裂隙內部，極容易堵塞瓦斯?jié)B流通道；② 沿壁面殘留分布的磨料顆粒增加了孔壁的摩擦阻力，容易導致卡鉆抱死等問題。因此，尋求一種既能增加射流沖擊破碎煤巖性能、又能有效避免上述問題的新型磨料射流成為該技術廣泛應用的瓶頸難題。文獻[4-5]研究發(fā)現(xiàn)，將冰塊破碎成細小顆粒后混入高壓水射流能有效進行物料清洗和表面除銹，這種新型射流被稱之為冰粒磨料射流。在相同工況下冰粒磨料射流的沖擊性能幾乎媲美常規(guī)磨料射流，且冰粒沖擊靶體后迅速融化可有效規(guī)避傳統(tǒng)磨料帶來的問題。由此可見，采用冰粒磨料水射流代替?zhèn)鹘y(tǒng)磨料水射流進行煤層割縫有望發(fā)揮磨料射流的磨削沖擊性能，同時避免傳統(tǒng)磨料射流帶來的堵塞滲流通道與卡鉆等問題。作為一種新型的磨料水射流技術，冰粒磨料射流沖擊性能的強弱是判定該技術能否有效破碎煤巖的關鍵，是該技術能否成功應用于煤層割縫的重要參數(shù)。針對高壓磨料射流沖擊下煤巖體的破碎規(guī)律問題，目前已有大量的文獻開展了研究。黃飛等通過大量的沖擊試驗研究了磨料水射流沖擊下具有層狀節(jié)理巖石的多尺度破碎形態(tài)。米建宇等采用SPH-FEM耦合算法模擬了后混合磨料水射流沖擊破碎巖石的規(guī)律，探討了沖擊壓力、磨料濃度與圍壓等因素對巖石破碎的影響規(guī)律。黃中偉等提出一種新型的液氮磨料射流，并開展了液氮磨料射流噴射高溫花崗巖室內實驗。劉勇等對比分析了磨料水射流與磨料氣體射流的破巖壓力，并基于統(tǒng)一強度理論構建了適用于磨料射流破巖的能量準則。劉送永等采用數(shù)值模擬手段研究了在有圍壓存在的情況下磨料射流的破巖規(guī)律，探討了不同大小圍壓對巖石破碎的影響。上述文獻表明，傳統(tǒng)磨料射流通過高速磨料的磨削作用與高速流體的沖刷作用，可以提升射流的沖擊破巖性能。然而，鑒于冰粒磨料的物性特征與常規(guī)磨料區(qū)別較大，其摻入射流后的動能大小與沖擊破巖規(guī)律都將與常規(guī)磨料射流有所差異。針對該方面的問題，目前鮮有文獻開展研究。因此，筆者擬開展冰粒磨料加速規(guī)律的理論推導，在此基礎上采用數(shù)值模擬的方法研究冰粒磨料射流的破碎煤巖規(guī)律。

1 冰粒磨料加速規(guī)律理論研究

根據(jù)冰粒磨料的摻入方式，可將冰粒磨料射流分為前混合與后混合冰粒磨料射流。其中，后混合磨料射流是將破碎的冰粒從噴腔直接引入，從而在高速水流的卷攜作用下獲得加速。在此過程中，細小冰粒受到溫度較高且速度極快水流的沖刷作用，極易融化。因此后混合冰粒磨料射流對于噴嘴內的低溫環(huán)境要求極為嚴格。目前，大多采用預混合的方式來制備冰粒磨料射流，即先將破碎冰粒與一部分冷卻水進行充分預混合，然后再摻入高速流體進行加速。冰粒獲得加速的過程幾乎都是在高壓管路中進行，最后高速冰粒磨料與水的混合流體在噴嘴中進行整形與束能后從噴嘴噴出。因此，本文后續(xù)所提到的冰粒磨料射流均是指預先混合的冰粒磨料射流。

1.1 冰粒磨料的受力分析

鑒于冰粒磨料的密度、硬度與水介質的親疏性以及冰粒消融等特性與常規(guī)磨料射流差異極大，因此冰粒磨料的加速規(guī)律對其最終的沖擊速度影響較大，并最終影響其沖擊動能的大小。

冰粒磨料的受力與加速過程極其復雜，為了方便推導其加速過程，筆者作如下假設：① 冰粒磨料加速過程主要發(fā)生在高壓管路內；② 冰粒磨料在高壓管路內的運動為一維流動；③ 冰粒之間不發(fā)生碰撞，且不影響水的流動；④ 高壓水在管路內的流動為定常流，且水溫保持恒定；⑤ 冰粒磨料為規(guī)則的球形顆粒，冰粒的溫度保持恒定，鑒于冰粒與高壓水環(huán)境之間的溫差恒定，本文假設冰粒的消融過程遵循線性消融規(guī)律，即球形冰粒的粒徑隨時間呈線性減小的趨勢。

=+

(1)

式中，為冰粒初始粒徑；為冰粒運動時間；為冰粒融化粒徑；對應為時刻時間節(jié)點；為冰粒隨時間的消融常數(shù)。

根據(jù)上述假設可知冰粒的加速過程主要發(fā)生在高壓管路內，依據(jù)文獻[18-19]可獲得冰粒在高壓管路內的主要受力情況如下：

(2)

其中，為慣性力；為壓差力；為Stokes阻力；為附加質量力；為冰粒密度；為冰粒速度；??為壓力梯度；為阻力系數(shù)；為射流密度；為射流速度；為冰粒初始面積；為冰粒融化面積。根據(jù)冰粒磨料在高壓管道內的運動特征可以推導出冰粒的平衡方程為

+++=0

(3)

將式(1)，(2)代入式(3)可以推導出冰粒磨料的加速度為

(4)

鑒于高速流體在高壓管路內的速度基本保持不變，流體的沿程動能損失可以忽略不計。依據(jù)伯努利方程可知=0。因此，式(4)可簡化為

(5)

1.2 冰粒磨料加速案例分析

依據(jù)微積分思想將冰粒磨料的運動時間切分成部分，冰粒磨料在第個微元時間段Δ內的加速度可視為常數(shù)。因此可以推導出冰粒磨料從時間節(jié)點到+1時的速度變化公式為

+1=+Δ

(6)

已知冰粒磨料的初速度為0，將式(5)代入式(6)便可計算出冰粒的速度表達式為

(7)

按照高壓管路內流速恒定的假設，射流的速度可表示為

(8)

其中，為管道內的流量，L/min；為管道的直徑；m。本文參照了常規(guī)前混合磨料射流發(fā)生裝置的相關參數(shù)，選取了流量、管徑阻力系數(shù)等射流參數(shù)，依據(jù)冰粒本身的物性選取了冰粒磨料參數(shù)，冰粒粒徑消融常數(shù)主要參考了文獻[16-17]，詳細參數(shù)見表1?；谑?6)～(8)編制計算程序，采用Matlab軟件進行迭代計算，獲得冰粒磨料的速度與加速度的變化規(guī)律如圖1所示。

表1 高壓管路中所用計算參數(shù)

冰粒磨料的加速度變化趨勢可以間接反映其在噴嘴內的受力情況。從圖1所示的冰粒磨料加速度在高壓管路內的變化曲線可知，冰粒磨料的加速大致可以分為3個階段：從冰粒磨料摻入高壓流體的高壓管路位置開始至0.14 m范圍段為冰粒磨料加速的第1階段，該階段的加速度迅速增加至峰值；從0.14～0.8 m范圍段為冰粒加速的第2階段，該階段的加速度從峰值呈指數(shù)形式急劇減?。粡母邏汗苈穬?.8 m位置至噴嘴出口的范圍段為冰粒加速的第3階段，該階段的加速度減小的趨勢變緩并逐漸趨于穩(wěn)定。根據(jù)冰粒磨料的加速度變化規(guī)律可以推測，冰粒磨料摻入高壓射流后的受力加速過程非常短暫，且主要集中在第1，2階段。

圖1 不同流量下冰粒磨料的加速度與速度變化曲線

通過對比常規(guī)磨料射流的加速度變化規(guī)律發(fā)現(xiàn)，常規(guī)磨料的加速度通常只有第2、第3階段，表明常規(guī)磨料一旦摻入高壓流體后其獲得的受力和加速度即達到峰值，隨后便逐漸減小。進一步分析表明，冰粒磨料的加速度之所以出現(xiàn)短暫的增加趨勢，主要是考慮了冰粒消融導致其粒徑逐漸減小，根據(jù)式(4)可知在較短時間內冰粒受力變化較小，而冰粒的粒徑有所減小，因此其加速度出現(xiàn)增加的趨勢；當達到峰值加速度后，冰粒受力急劇減小，此時雖然冰粒粒徑也按照式(1)的趨勢減小，但其減小幅度逐漸小于受力減小的幅度，加速度出現(xiàn)減小趨勢。

另一方面，根據(jù)冰粒磨料的速度變化曲線可知，冰粒磨料在混入高壓射流后，其速度急劇增加，約1 ms后速度增加的趨勢變緩并逐漸趨近于水的速度，最后與高壓水一起從噴嘴噴出形成高壓冰粒磨料射流。根據(jù)冰粒磨料的加速度和速度變化規(guī)律可知，鑒于冰粒的特殊性能獲得很好的加速過程，最終噴出速度幾乎可以接近水的速度，從而能夠極大地提升冰粒磨料射流的沖擊動能。

2 冰粒磨料射流沖擊破巖模型構建

2.1 SPH-FEM耦合方法

SPH-FEM的耦合主要采用罰函數(shù)法使質點的力作用在有限元上，從而實現(xiàn)SPH與有限元的耦合。在計算過程中，粒子與有限元表面保持接觸，避免穿透。如果粒子穿透了有限元表面，則粒子的位置和速度根據(jù)以下等式進行調整：

(9)

(10)

式中，為界面上有限元的壓力；為單元中心附近的粒子壓力；為元素中心附近的粒子數(shù)。

2.2 水的本構模型

筆者采用Null材料模型并賦予其狀態(tài)方程Mie-Grueisen，該狀態(tài)方程可表述為

(11)

其中，為流體所受的壓力；為水的初始密度；為水的黏滯系數(shù)；為單位體積內能；為質點速度與聲波速度(-)曲線截距；，與為-曲線斜率系數(shù)；為Grueisen常數(shù)；為一介體積修正量。筆者選取表2所示的材料參數(shù)來模擬水的本構關系。

表2 水的本構模型參數(shù)

2.3 冰的本構模型

筆者只考慮高速沖擊下冰的力學行為，因此采用CARNEY等所提出的本構模型來描述冰粒磨料。該本構模型不僅考慮了冰的拉伸屈服強度和壓縮屈服強度，而且還考慮了冰在拉伸狀態(tài)和壓縮狀態(tài)下的應變率效應。

(12)

(13)

式中，為冰的參考體積；為冰的體積；與分別為冰的壓縮與拉伸參數(shù)；為系數(shù)。

此外，該本構模型考慮了拉伸和壓縮過程中冰粒的階段壓力，計算公式為

(14)

表3 冰的材料參數(shù)

表4 冰的狀態(tài)方程參數(shù)

2.4 煤的本構模型

筆者采用HJC本構模型來模擬高速冰粒磨料射流沖擊下的煤巖，詳細材料參數(shù)見表5。其中，為煤巖密度；為單軸抗壓強度；為最大靜水抗拉強度；為剪切模量；為規(guī)范化最大強度；為允許塑性應變；為壓垮靜水壓力；為壓垮體積應力；為壓垮體積應變；為極限體積應變；，，，，，，，，為材料常數(shù)；為參考應變率。

表5 煤的本構參數(shù)

2.5 幾何模型構建與速度賦值

為了便于統(tǒng)計冰粒磨料濃度對射流沖擊破煤的影響規(guī)律，筆者將水與冰粒磨料采用SPH方法構建，將煤塊采用FEM方法構建。為了模擬真實的磨料在水中隨機分布的情況，采取如下建模方法：

第1步，建立一個磨料水射流圓柱幾何模型，在該幾何模型中按照體積分數(shù)比分別建立純水SPH模型和冰粒磨料SPH模型。

第2步，建立好2段SPH模型后，提取2段SPH粒子在空間中的坐標信息，然后采用隨機函數(shù)編程將2段SPH的坐標信息打亂，從而實現(xiàn)冰粒磨料在水射流中的隨機分布。構建的冰粒磨料射流沖擊破碎煤塊的幾何模型如圖2所示。

圖2 冰粒磨料水射流沖擊煤巖模型

根據(jù)2.4節(jié)的理論研究結果可知，冰粒磨料與純水介質從噴嘴噴出時，兩者的速度幾乎相同。本文忽略2個速度間的微小差異，將冰粒磨料與純水介質設置為相同的初速度進行模擬計算。

3 冰粒磨料射流破碎煤塊效果分析

按照上述方法建立數(shù)值模型，然后采用LS-DYNA進行模擬計算。圖3為冰粒磨料射流沖擊破碎煤塊的時間序列。從圖3可以看出，在冰粒磨料射流的沖擊下，煤塊中心形成V型的破碎坑，破碎坑外圍出現(xiàn)不同程度的損傷。隨著沖擊時間的增加，破碎坑持續(xù)向煤塊深部發(fā)展，損傷范圍也相應擴大。上述煤巖破碎現(xiàn)象與現(xiàn)有大量文獻獲得的破巖規(guī)律幾乎一致，表明本文所建立的數(shù)值計算模型具有一定可靠性。

圖3 煤巖損傷時序

另一方面，從圖3可以看出，當冰粒磨料高速沖擊煤巖后，部分冰粒迅速融化并轉換為液態(tài)的SPH水介質。這一現(xiàn)象與實際冰粒磨料射流沖擊煤巖后的消融規(guī)律相符，進一步表明本文所建立的冰粒磨料射流沖擊破碎煤巖的可行性。

為了定量獲取煤塊在冰粒磨料射流沖擊下的損傷規(guī)律，本文以煤巖受沖擊中心為基點，沿模型軸方向截取煤巖的1/2剖面，沿中心軸向下分別設置壓力監(jiān)測點1號與應力監(jiān)測點2號，記錄所選監(jiān)測點損傷隨時間的變化。此外，在距中心軸15 mm的徑向距離范圍內連續(xù)選擇38個損傷監(jiān)測點，標記為3～40號，并記錄這38個監(jiān)測點的最終損傷。同樣的，沿軸軸向距離處選取40個損傷監(jiān)測點標記為41～80號，以探究煤巖在3種水射流沖擊下的損傷范圍，詳細的監(jiān)測點布局如圖4所示。

圖4 監(jiān)測點布置示意

3.1 不同速度冰粒磨料射流破碎煤塊規(guī)律

為探究冰粒磨料射流速度對煤巖損傷規(guī)律的影響，筆者選取了速度分別為300，400，500，600，700，800，900 m/s的冰粒磨料射流進行破巖試驗，冰粒體積分數(shù)為12.5%，設置沖擊時間為2 s，不考慮圍壓作用。計算完成后提取了部分煤巖的損傷云圖，如圖5所示。由圖5可知，隨著射流速度的增大，煤巖V型破碎坑體積增大，破碎坑附近的損傷范圍也隨之增大。

圖5 不同速度沖擊下煤巖損傷云圖

為了定量獲取射流速度對煤巖宏觀破碎坑的影響規(guī)律，筆者采用軟件提取了煤巖的破碎深度與破碎深度數(shù)據(jù)，并繪制了兩者隨射流速度的變化關系，如圖6所示。從圖6可以看出，隨著射流速度的增加，煤巖的破碎深度呈現(xiàn)近似線性增長的趨勢。當冰粒磨料射流速度為300 m/s時，煤塊的最大破碎深度約為7.5 mm；當冰粒磨料射流速度增加至900 m/s時，煤塊的最大破碎深度達到了約30 mm。反觀煤塊的破碎寬度出現(xiàn)了不同的變化趨勢。具體來說，當冰粒磨料射流速度在300～500 m/s時，煤塊破碎深度隨速度增加呈現(xiàn)小幅度的增加趨勢；當冰粒磨料射流速度進一步增大，破碎寬度不再隨速度增加而有明顯的增加趨勢。分析研究表明，冰粒磨料射流速度增加，其沖擊動能=/2呈倍數(shù)關系增加，導致煤塊受到的沖擊能量增加，最終其破碎深度線性增加。然而，隨著冰粒磨料速度的增加，其沖擊磨削性能得到提升，從而能更好地磨削煤塊形成更加完整的切口，促使沖擊破碎向煤塊內部發(fā)展，而破碎寬度不再增加；另一方面，冰粒磨料射流的直徑不變，阻礙了破碎寬度的進一步擴展。

圖6 射流速度與破碎寬度和深度的關系

此外，為了掌握不同速度的冰粒磨料射流沖擊下煤塊的損傷發(fā)展規(guī)律，提取監(jiān)測點41～80號的損傷值，并繪制了不同速度下煤巖的損傷曲線，如圖7所示。定義損傷的取值范圍為0～1，其中，0代表沒有受到損傷，具體表現(xiàn)為網(wǎng)格單元未有任何形變；1代表煤塊被剝離，具體表現(xiàn)為單元網(wǎng)格出現(xiàn)大變形并達到損傷極限后被刪除。根據(jù)上述定義可知，煤塊縱向破碎坑的最大深度即損傷值為1時對應的縱向距離，而煤塊損傷的最大深度即損傷值為0時對應的縱向距離。

根據(jù)圖7所示，通過計算獲得如下關系式：

圖7 不同射流速度下煤巖損傷范圍

(15)

其中，Δ為縱向相對損傷增量，Δ=(-)/×100%，%。根據(jù)式(15)可以看出，冰粒磨料射流沖擊下，煤塊的縱向損傷增量保持在(0.129～0.133)內，且與冰粒磨料的速度沒有明顯關系。由此推測，冰粒磨料射流沖擊煤巖的最大損傷深度約為破碎坑深度的1.13倍。

3.2 不同體積分數(shù)冰粒磨料射流破碎煤巖規(guī)律

為研究冰粒磨料體積分數(shù)對煤巖損傷規(guī)律的影響，筆者選取了冰粒磨料的體積分數(shù)分別為0，12.5%，25.0%，37.5%，50.0%，進行破巖數(shù)值模擬試驗，射流速度為400 m/s，設置沖擊時間為2 s，不考慮圍壓作用。計算完成后提取了各冰粒體積分數(shù)下煤塊的破碎與損傷云圖，如圖8所示。當純水射流中加入體積分數(shù)為12.5%的冰粒磨料后，煤巖的破碎坑體積與損傷范圍明顯增大；然而，隨著冰粒磨料體積分數(shù)的繼續(xù)增加，煤巖破碎與損傷的范圍變化不明顯。根據(jù)上述結果可以推測，冰粒磨料的摻入能夠明顯提升射流沖擊破碎煤巖的性能。

圖8 不同磨料體積分數(shù)下煤巖損傷云圖

為了定量研究冰粒體積分數(shù)對煤巖宏觀破碎坑的影響規(guī)律，依次提取了不同冰粒體積分數(shù)下煤巖的破碎深度與破碎寬度數(shù)據(jù)，并繪制了2者隨射流速度的變化曲線，如圖9所示。

圖9 磨料濃度與破碎寬度和深度的關系

由圖9可知，在相同工況下的純水射流中摻入12.5%的冰粒磨料后，煤巖破碎坑的最大深度由約6 mm增加至約9.4 mm，增幅約為56.7%；破碎坑寬度由約9.6 mm增加至約13.6 mm，增幅約為41.7%。繼續(xù)增加冰粒體積分數(shù)至37.5%，煤巖破碎坑的最大深度約為11 mm，增幅約為17%；破碎坑寬度約為14.8 mm，增幅約為8.8%。繼續(xù)增加冰粒體積分數(shù)，煤巖破碎坑的最大深度和寬度幾乎不再增加。分析認為，出現(xiàn)上述變化趨勢的主要原因是純水射流與冰粒磨料射流沖擊破巖的作用不同。純水射流主要借助高速流體的水錘壓力和滯止壓力的沖刷作用破碎煤巖，而摻入冰粒磨料后，冰粒的磨削與切削作用主導了破巖作用，更加有利于煤巖的破碎。因此，在純水射流中摻入冰粒磨料后煤巖的破碎深度和寬度都大幅增加。由此可以推測，針對特定物性參數(shù)的煤巖體，存在一個破碎效果最佳的冰粒磨料體積分數(shù)。

為了進一步了解冰粒體積分數(shù)對煤巖損傷寬度的影響規(guī)律，提取監(jiān)測點3～40號的損傷值，并繪制了如圖10所示的損傷曲線。同理定義破碎坑的寬度為損傷值為1時對應的橫向距離，損傷的最大寬度為損傷值為0時對應的橫向距離。根據(jù)圖10進行計算獲得如下關系式：

圖10 不同體積分數(shù)冰粒下煤巖損傷范圍

(16)

其中，為冰粒磨料的體積分數(shù)，%；Δ為橫向相對損傷增量，Δ=(-)/×100%，%。根據(jù)式(16)可以看出，在磨料射流的沖擊下，煤巖橫向損傷直徑約為破碎坑直徑的1.3倍。

3.3 不同圍壓下冰粒磨料射流破碎煤巖規(guī)律

在具體的煤層割縫實踐中，煤巖均處于高圍壓限制約束狀態(tài)。為了使本數(shù)值模擬更加貼合工程實際，在數(shù)值模型中對煤巖4個側面均施加相同的壓力，模擬煤巖受到的圍壓狀態(tài)，分別對煤巖施加1，2，3，4，5，6，7，8，10 MPa的圍壓進行計算，選取射流速度為400 m/s，冰粒體積分數(shù)為12.5%，設置沖擊時間為2 s。計算完成后，同樣提取部分煤巖損傷云圖，如圖11所示。

從圖11可以看出，煤巖的破碎與損傷范圍隨著圍壓的增加出現(xiàn)明顯減小的趨勢。

圖11 不同圍壓狀態(tài)下煤巖損傷云圖

通過軟件提取破碎深度與寬度，并繪制2者隨圍壓的變化曲線如圖12所示。從圖12可以看出，在冰粒磨料射流沖擊下，煤巖的縱向破碎深度由無圍壓時的約15.2 mm降低至1 MPa圍壓時的13.1 mm；繼續(xù)增加圍壓，縱向破碎深度呈近似線性的趨勢減小。然而，煤巖的橫向破碎寬度隨著圍壓的變化始終在9.9 mm上下波動，未有表現(xiàn)出隨圍壓而變化的趨勢。通過上述現(xiàn)象可知，在冰粒磨料射流沖擊下，圍壓阻礙了煤巖的縱向破碎損傷，但是對于橫向的破碎損傷沒有影響。

圖12 圍壓與破碎寬度和深度的關系

提取監(jiān)測點41號在不同圍壓條件下的損傷隨時間的變化數(shù)據(jù)，并繪制該點的時序損傷，如圖13所示。從圖13可以看出，隨著圍壓的增加，煤巖在冰粒磨料射流沖擊下的損傷明顯受到了抑制作用，具體表現(xiàn)為圍壓越大，損傷需要的時間越長、損傷的程度越小。另外，從圖13還可看出，在冰粒磨料射流沖擊下，煤巖的損傷幾乎是階梯狀發(fā)展的。分析認為，這可能是由于煤巖被逐次剝離從而露出新的接觸面所導致的。

圖13 不同圍壓下煤巖的時序損傷

4 結 論

(1)冰粒磨料摻入高速流體后在高壓管路內獲得了先增大后減小的加速度，最終其速度幾乎接近純水速度，為后續(xù)數(shù)值模擬中冰粒磨料和純水的速度賦值提供了依據(jù)。

(2)冰粒磨料在沖擊煤巖前后出現(xiàn)消融現(xiàn)象，且煤巖出現(xiàn)V型破碎形態(tài)，表明本文所構建的冰粒磨料射流破碎煤巖的數(shù)值模型具可行性。

(3)煤巖縱向破碎深度隨冰粒磨料射流速度的增加而線性增加，而橫向破碎寬度隨速度呈現(xiàn)先增加后趨于穩(wěn)定的變化趨勢，煤巖縱向損傷深度約為縱向破碎深度的1.13倍；冰粒磨料摻入純水射流后，煤巖的縱向破碎深度和橫向破碎寬度均出現(xiàn)了大幅增加，隨著冰粒磨料體積分數(shù)的增加，破碎深度和破碎寬度增加的趨勢變緩，最后煤巖橫向損傷直徑約為橫向破碎寬度的1.3倍；圍壓越大縱向破碎深度越小，表明煤巖所受圍壓極大地抑制了其在冰粒磨料射流沖擊下的縱向破碎深度。