齊慶杰，趙尤信，周新華，沈志遠(yuǎn)，賈新雷

(1.遼寧工程技術(shù)大學(xué) 安全科學(xué)與工程學(xué)院，遼寧 阜新 123000；2.礦山熱動力災(zāi)害與防治教育部重點實驗室，遼寧 阜新 123000)

0 引言

我國煤炭賦存的地質(zhì)條件比較復(fù)雜，其開采難度較大，我國重點煤礦中，有三分之二以上的煤礦屬于高瓦斯突出礦井，大部分煤層是低透氣性煤層，其瓦斯抽放難度較大是我國煤礦安全開采的一大難題[1]。為了提高煤層的滲透率，國內(nèi)外學(xué)者采取了多種措施進(jìn)行瓦斯抽放[2-7]，其中深孔聚能爆破增透技術(shù)[8-15]提高了煤層的滲透率以及裂隙發(fā)育，在平煤六礦、十礦和焦煤九里山礦、鄭州大平煤礦進(jìn)行了應(yīng)用，取得了良好的效果。

為了提高深孔聚能爆破瓦斯抽放效率及其炸藥的利用率，很多學(xué)者對爆破參數(shù)[15-16]進(jìn)行研究，其中在優(yōu)化爆破孔孔徑、軸向裝藥長度、封孔參數(shù)、炸藥性質(zhì)以提高煤層的透氣性方面取得了成功，并在許多煤礦進(jìn)行應(yīng)用，但深孔聚能爆破增透技術(shù)參數(shù)優(yōu)化中沒有對裝藥位置進(jìn)行研究。裝藥位置的選擇直接影響深孔爆破增透的增透半徑、裂隙發(fā)育程度和爆破安全。

筆者基于對深孔聚能爆破理論研究和現(xiàn)場實踐，以平煤十礦24130工作面為例，對深孔聚能爆破裝藥位置進(jìn)行了數(shù)值模擬研究，結(jié)合現(xiàn)場試驗，對比不同裝藥位置爆破后周圍應(yīng)力分布和裂隙變化，以期為優(yōu)化裝藥位置提供理論參考。

1 數(shù)值模擬理論研究

1.1 深孔聚能爆破作用

普通裝藥模式的結(jié)構(gòu)圖如1(a)所示，其爆破產(chǎn)生的爆生氣體壓力分布均勻，爆生氣體不定向地向四周擴(kuò)散開，這種現(xiàn)象導(dǎo)致爆炸所產(chǎn)生的裂隙范圍較小，效果不明顯。聚能爆破的結(jié)構(gòu)圖如1(b)所示，其爆生氣體壓力更加集中于聚能方向，爆生氣體集中向聚能槽移動，因此，聚能爆破更好地避免了能量的分散，在聚能流方向上的動能和速度更高，致裂效果更加明顯。

圖1 爆破裝藥結(jié)構(gòu)Fig.1 Explosive charge structure

相比于普通裝藥模式，聚能裝藥模式下，爆炸后在聚能流方向上的動能和速度更高，聚能流沖擊煤體產(chǎn)生初始裂隙，爆破產(chǎn)生的裂隙方向由聚能方向決定，裂隙的深度不僅受聚能流的能量控制，而且受煤體強(qiáng)度的影響，由于煤體的強(qiáng)度是定值，則其裂隙深度可以用以下公式表示[17-18]：

(1)

式中：L為裂隙深度，m；L′為有效射流長度，m；ρj為聚能流密度，kg/m3；ρt為巖體密度，kg/m3。由式(1)可見，聚能流對煤體裂隙的發(fā)育起到定向作用，而后在爆炸能量作用下有利于裂隙在聚能方向充分?jǐn)U展，聚能爆破在非聚能方向上煤巖體致裂范圍變大。

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

聚能爆破數(shù)值模擬是聚能流、煤體和空氣之間相互作用的流固耦合過程，本文采用ALE算法，此算法能夠效結(jié)合Lagrange方法和Euler方法的優(yōu)點，解決大變形計算的難題。

1.3 材料和參數(shù)設(shè)置

模擬的炸藥選取的是LS-DYNA中的高能爆炸燃燒材料，此種材料爆轟后滿足JWL狀態(tài)方程，JWL狀態(tài)方程對爆轟產(chǎn)物的膨脹做功進(jìn)行準(zhǔn)確的描述，得到任意時刻被爆轟物體的爆轟壓力P為：

(2)

式中：P為爆轟壓力，MPa；A,B,ω,R1,R2均為炸藥性質(zhì)參數(shù)，ω=0.3；R1=4.12；R2=0.95；V為氣體產(chǎn)物的質(zhì)量體積，m3;E0為氣體產(chǎn)物爆炸時的能量，MJ。

炸藥相關(guān)參數(shù)如下：密度為0.95 ×103kg/m3；爆速為2 800 m/s；A=347 GPa；B=0.433 GPa。

上邊界頂壓如式(3)：

q=γgH

(3)

式中：q為頂壓，kN/m2；γ為巖土層的平均體積力，取γ=25 kN/m3；H為煤層埋深，480 m。計算可得q=117 600 kN/m2。

煤體的物理力學(xué)參數(shù)通過現(xiàn)場數(shù)據(jù)得到，具體參數(shù)如表1所示。

表1 煤體物理力學(xué)參數(shù)Table 1 Coal mechanical parameters

2 數(shù)值模擬

2.1 模型建立與網(wǎng)格劃分

針對平煤十礦己15，16-24130工作面，應(yīng)用ANSYS/LS-DYNA軟件，共建立3個幾何模型，如圖2所示。模型尺寸為40 m×40 m×60 m，孔深45 m，炸藥長度13 m，封孔段長度12 m，爆破孔孔徑為89 mm，藥包直徑為42 mm。其中，裝藥位置分別在爆破孔里端(距離封孔段20 m)、中間(距離封孔段10 m)和外端(距離封孔段0 m)，圖2(a)為裝藥位置在里端的幾何模型(另外2個模型只改變裝藥位置)，聚能方向為X平面方向，其網(wǎng)格劃分如圖2(b)所示，單元總數(shù)為361 210個，節(jié)點總數(shù)為378 970個。

圖2 模型及網(wǎng)格劃分Fig.2 Model and grid partition

2.2 模擬結(jié)果與分析

不同裝藥位置的爆破應(yīng)力YZ平面應(yīng)力云圖如圖3所示。由圖3可得，聚能爆破后，應(yīng)力波以球面波的形式向四周擴(kuò)散，炸藥中心附近應(yīng)力值最大，能量最集中。裝藥位置不同，應(yīng)力作用范圍有明顯差異。采取里端裝藥時，煤體里端所受應(yīng)力范圍較大，一定范圍之外的上部煤體并沒有受到明顯應(yīng)力作用；采取外端裝藥時，煤體外端所受應(yīng)力范圍較大，一定范圍之外的下部煤體并沒有受到明顯應(yīng)力作用；采取中端裝藥時，除了中部煤體所受應(yīng)力之外，上下部分煤體均受到應(yīng)力作用，范圍較大。分析認(rèn)為，炸藥爆炸后產(chǎn)生的爆生氣體除了作用于周圍煤體外，還會與炮孔內(nèi)空氣層相互作用，在空氣內(nèi)形成反射波持續(xù)作用于煤體。中端裝藥形式下，爆破孔上下部分均存在空氣層，爆轟波與上下部分空氣發(fā)生反射，作用時間較長，能量傳遞范圍廣；外端裝藥時藥柱緊鄰封孔段，容易產(chǎn)生能量損失；里端裝藥時藥包緊鄰下方煤體，不存在空氣層，爆生氣體作用于下方煤體后不存在應(yīng)力波的反射，作用范圍較小。

圖3 t=150 ms時不同裝藥位置爆破應(yīng)力云圖Fig.3 Blasting stress cloud map at different loading locations at t=150 ms

為了更準(zhǔn)確地分析出3種裝藥形式下爆破孔一定范圍內(nèi)應(yīng)力的變化，在爆破孔附近選取A，B，C，D，E5個觀察點，C點距離藥包中心O水平距離3 m，其余4點與C點在同一條豎直線上，各觀察點間距6 m。應(yīng)用LS-PrePost進(jìn)行后處理分析，提取出不同裝藥位置時5個觀察點的應(yīng)力時程變化。由提取結(jié)果可知，不同裝藥位置時，爆破有效區(qū)域內(nèi)的觀察點應(yīng)力線變化趨勢大體相同，都存在2個應(yīng)力峰值，且第1個應(yīng)力峰值要明顯大于第2個應(yīng)力峰值。這是因為炸藥起爆后聚能射流首先作用于煤體，產(chǎn)生較大應(yīng)力，而后爆轟波與空氣層產(chǎn)生反射波第2次作用于煤體，產(chǎn)生第2個應(yīng)力峰值。里端裝藥時，距離炸藥中心最近的C點最先在100 ms時到達(dá)第1個應(yīng)力峰值，約為106 MPa，隨后在空氣層的作用下，約195 ms后出現(xiàn)第2個應(yīng)力峰值；162 ms時B，D這2個點到達(dá)第1個應(yīng)力峰值，大小約為29 MPa，隨后在空氣層的作用下53 ms后第2次出現(xiàn)應(yīng)力峰值，大小約為34 MPa；距離炸藥較遠(yuǎn)處的A，E兩點應(yīng)力值并未發(fā)生明顯變化，這表明A，E兩點不在爆破影響范圍內(nèi)。外端裝藥時，距離炸藥中心最近的C點第1次達(dá)到應(yīng)力峰值的時間為105 ms，大小約為86.2 MPa，較里端裝藥時應(yīng)力值小，這是因為藥柱緊鄰封孔段，產(chǎn)生了一定的能量損失，160 ms時，D點第1次出現(xiàn)應(yīng)力峰值，大小約為40 MPa，隨后在空氣層的作用下50 ms后第2次出現(xiàn)應(yīng)力峰值，A，B，E點均未發(fā)生明顯應(yīng)力變化，表明A，B，E點不在爆破影響范圍內(nèi)。中端裝藥時，C點第1次達(dá)到應(yīng)力峰值的時間為95 ms，應(yīng)力峰值約為154 MPa，A，B，D，E點處均存在應(yīng)力的變化，應(yīng)力平均值約為49 MPa，表明所有觀察點均在爆破影響范圍內(nèi)，距離爆破點較近的B，D點約在87 ms時出現(xiàn)第1次應(yīng)力峰值，大小約為98 MPa；距離較遠(yuǎn)處的A，E點在110 ms時出現(xiàn)第1次應(yīng)力峰值，大小約為48 MPa。由以上分析可以得出，中端裝藥時，爆破有效影響范圍更大，應(yīng)力作用時間更持久，里端裝藥較小，外端裝藥時最小。

應(yīng)用LS-PrePost提取出不同時刻不同裝藥位置下煤體爆破產(chǎn)生的軸向和徑向的裂隙深度，繪制裂隙深度變化曲線圖，如圖4所示。

圖4 不同裝藥位置裂隙變化Fig.4 Curve diagram of fracture change in different loading position

由圖4可以看出：采取里端裝藥形式時，爆破產(chǎn)生的軸向裂隙(爆破孔方向)約在380 ms時趨于穩(wěn)定，最大值可以達(dá)到22 m，徑向裂隙在350 ms時趨于穩(wěn)定，最大值可以達(dá)到6.4 m；采取中端裝藥時，爆破產(chǎn)生的軸向裂隙在450 ms時趨于穩(wěn)定，應(yīng)力波作用于煤體的時間比較長，有效促進(jìn)了裂隙發(fā)育，最大值達(dá)到27.5 m，相比于里端裝藥和外端裝藥分別增加了5.5和7.4 m，平均增加6.5 m。徑向裂隙在450 ms時趨于穩(wěn)定，最大深度達(dá)到6.8 m，相比于里端裝藥和外端裝藥分別增加了0.4和1 m，平均增加0.7 m；采取外端裝藥時，爆破產(chǎn)生的軸向裂隙在370 ms時趨于穩(wěn)定，最大值達(dá)到20.1 m，徑向裂隙在400 ms時趨于穩(wěn)定，最大值達(dá)到5.8 m。

3 工程應(yīng)用

3.1 工作面概況

24130工作面所采煤層為己15.16煤層，該工作面己15與己16煤層合層，厚度為2.5～3.5 m，平均煤厚3.2 m，煤層傾角10°～12°，最大瓦斯原始壓力值為1.0 MPa，己15.16煤層瓦斯原始含量值為12.5 m3/t。

3.2 工程試驗

深孔聚能爆破的目的就是增大媒體內(nèi)裂隙，為瓦斯流動提供通道，提高煤層透氣性。因此，瓦斯?jié)舛群屯咚沽髁康脑黾臃仁强疾炀勰鼙菩Ч闹匾笜?biāo)。

具體實施聚能爆破鉆孔和抽采鉆孔布置圖如5-6所示。

圖5 鉆孔布置平面圖Fig.5 Implementation of borehole plan

圖6 鉆孔布置剖面圖Fig.6 Implementation of borehole layout profile

在24130工作面風(fēng)巷無斷層影響區(qū)域，進(jìn)行3組深孔聚能爆破，其中1#裝藥位置中端裝藥、2#在里端裝藥和3#在外端裝藥，孔深45 m，炸藥長度13 m，封孔段長度12 m，爆破孔孔徑為89 mm，藥包直徑為42 mm。聚能爆破鉆孔周圍4個鉆孔為抽采鉆孔，抽采負(fù)壓為20 kPa，鉆孔半徑為94 mm，采用“囊袋式兩堵一注”封孔器，1#爆破孔周圍4個抽采鉆孔為第1組，2#爆破孔周圍4個抽采鉆孔為第2組，3#爆破孔周圍4個抽采鉆孔為第3組。

得到深孔聚能爆破不同裝藥位置抽采鉆孔的平均瓦斯?jié)舛群推骄咚辜兞孔兓鐖D7-8所示。

圖7 不同裝藥位置聚能爆破孔周圍抽采鉆孔瓦斯?jié)舛茸兓疐ig.7 Gas concentration variation curves around boreholes in shaped charge holes at different charging locations

圖8 不同裝藥位置聚能爆破孔周圍抽采鉆孔瓦斯純量變化Fig.8 Variation curves of gas purity in boreholes around shaped charge holes at different charging locations

由圖7-8可知，聚能爆破前，爆破孔周圍抽采鉆孔的瓦斯?jié)舛仍?0%～37.5%之間，說明爆破孔周圍的所選區(qū)域媒體裂隙發(fā)育程度基本相同。爆破后裝藥位置在中間時，爆破后瓦斯?jié)舛茸畲竽軌蜻_(dá)到83.2%，瓦斯純量0.009 7 m3/min,17 d內(nèi)平均瓦斯?jié)舛葹?2.3%，平均瓦斯純量為0.008 4 m3/min，該區(qū)域抽采鉆孔平均瓦斯?jié)舛仁潜魄暗?.99倍，瓦斯純量為爆破前的4.08倍。爆破后裝藥位置在里端時，爆破后瓦斯?jié)舛饶軌蜻_(dá)到75.1%，瓦斯純量0.007 8 m3/min，17 d內(nèi)平均瓦斯?jié)舛葹?3.4%，平均瓦斯純量為0.006 7 m3/min，該區(qū)域抽采鉆孔平均瓦斯?jié)舛仁潜魄暗?.61倍，瓦斯純量為爆破前的2.96倍。爆破后裝藥位置在外端時，爆破后瓦斯?jié)舛茸畲竽軌蜻_(dá)到66.2%，瓦斯純量0.005 7 m3/min，17 d內(nèi)平均瓦斯?jié)舛葹?4.5%，平均瓦斯純量為0.004 6 m3/min，該區(qū)域抽采鉆孔平均瓦斯?jié)舛仁潜魄暗?.32倍，瓦斯純量為爆破前的2.23倍。依據(jù)抽采瓦斯?jié)舛群屯咚辜兞?，不同的裝藥位置中間(距離封堵段10 m)部分裝藥效果最好，具有良好的增透效果。

4 結(jié)論

1)在24130工作面風(fēng)巷進(jìn)行的不同裝藥位置深孔聚能爆破實驗中，中部裝藥17 d內(nèi)的平均瓦斯抽采濃度和瓦斯純量分別是爆破前的1.99倍和4.08倍，瓦斯抽采濃度和瓦斯純量分別是里端裝藥的1.23倍和1.37倍，是外端裝藥的1.5被和1.83倍。平煤十礦24130工作面中間裝藥在深孔聚能爆破的應(yīng)用中增透效果最好，工程實驗與數(shù)值模擬的結(jié)果基本一致。

2)不同裝藥位置相比，中端裝藥時爆破有效影響范圍較大、應(yīng)力作用時間較長，里端裝藥時較小，外端裝藥時最小。

3)聚能爆破中端裝藥，其爆轟波和應(yīng)力波作用于煤體的時間較長且應(yīng)力值更大，爆破產(chǎn)生的裂隙深度主要在軸向上加深，軸向裂隙和徑向裂隙比其他裝藥方式均有增加。爆炸產(chǎn)生的爆生氣體與空氣波二次作用產(chǎn)生反射波，延長了作用時間，有助于裂隙擴(kuò)展，加速了賦存瓦斯解吸和流動，爆破效果更好。