王文洋，羅志華，張宗國

(1.深圳市中金嶺南有色金屬股份有限公司 凡口鉛鋅礦， 廣東 韶關(guān)市 512000；2.中南大學(xué) 資源與安全工程學(xué)院， 湖南 長沙 410083)

0 引言

在開采過程中，由于孔底抵抗線過大或過小，使得爆炸能量過多的浪費(fèi)，容易出現(xiàn)塊度分布不均勻的爆破效果與造成堵孔沖孔等爆破危害（見圖1（c））。為了最大化合理利用爆炸能量以提高爆破效果，改善礦山經(jīng)濟(jì)效益，研究孔底抵抗線以控制爆炸效果對現(xiàn)場生產(chǎn)有很大的實踐意義。許多學(xué)者對爆破抵抗線進(jìn)行了分析研究，如：姜永恒等[2]針對某金礦大塊率高的問題，運(yùn)用數(shù)值模擬軟件對采場扇形孔進(jìn)行模擬，分析了巖體內(nèi)部質(zhì)點的有效應(yīng)力峰值，確定了中深孔爆破最優(yōu)抵抗線和孔底距；何闖等[3]研究抵抗線與自由面對爆破振動的影響，得出爆破振動受自由面與最小抵抗線的共同影響；沈慧明等[4]采用理論分析的方式，得出適當(dāng)調(diào)整抵抗線可以改善爆破效果；何士靄[5]通過合理的裝藥設(shè)計，使得爆破塊度、爆堆位移在合理范圍之內(nèi)；李育等[6]采用減小前排裝藥抵抗線來控制爆破振動；魏兆云等[7]通過理論計算的方式確定爆破中抵抗線參數(shù)來控制爆破效果。綜上可知，多數(shù)研究均采用數(shù)值模擬的手段研究抵抗線，為了研究大直徑深孔爆破中孔底抵抗線對爆破效果的影響，本文采用數(shù)值模擬方法，對不同孔底抵抗線下的等效漏斗半徑與等效漏斗體積進(jìn)行對比分析，優(yōu)選合理的孔底抵抗線以改善現(xiàn)場爆破效果。

圖1 采礦工藝

1 數(shù)值模擬

1.1 模型建立

為了研究孔底抵抗線對爆破效果的影響，采用數(shù)值模擬方法分析爆炸過程中試塊內(nèi)部應(yīng)力變化情況，數(shù)值模擬中的炸藥、模型材料等與現(xiàn)場實際基本一致，模擬抵抗線為變量。本文的主要研究對象為抵抗線，為了方便劃分網(wǎng)格，數(shù)值模擬中將原藥包長度為0.6 m通過幾何相似比例縮放為3 cm，相似系數(shù)為20，模擬共分為5組，孔底抵抗線分別為3 cm、4 cm、5 cm、6 cm、7 cm。

如圖2所示，模型長、寬、高分別為40 cm、40 cm、30 cm，炸藥直徑為8 mm，長度為3 cm；其中孔底抵抗線大小為變量，其中藥包上下層均為堵塞物，模型中ALE空間的直徑為5倍炮孔直徑，ALE與部分石塊單元重合。模型中忽略重力的影響。巖石模型單元大小為5 mm×5 mm，炸藥、堵塞、ALE采用漸變網(wǎng)格。

圖2 四分之一模型示意

1.2 材料參數(shù)選取

1.2.1 炸藥的材料模型與其狀態(tài)方程

在LS-DYNA程序之中，炸藥材料模型采用8號材料模型，炸藥的狀態(tài)方程為“*EOS_JWL”，此方程忽略炸藥的化學(xué)反應(yīng)方式而只考慮爆轟，經(jīng)過大量試驗數(shù)據(jù)修正，JWL狀態(tài)方程如下[8-10]：

式中，P為爆轟壓力，Pa；E0為單元最大面積，m2；V為爆轟產(chǎn)物相對體積，m3；ω為格林愛森參數(shù)；A、B、R1、R2均為表征炸藥物理特性的常數(shù)。

1.2.2 混凝土RHT材料模型

本文數(shù)值模擬采用混凝土 RHT本構(gòu)模型[12]，該模型中的彈性極限面方程、最大失效面方程和殘余強(qiáng)度方程可用于描述混凝土在爆炸沖擊荷載作用下的初始屈服強(qiáng)度、失效強(qiáng)度和殘余強(qiáng)度的變化規(guī)律。

1.2.3 炮泥材料模型

卷積神經(jīng)網(wǎng)絡(luò)可以用來處理并且分析多種多樣的圖像，卷積神經(jīng)網(wǎng)絡(luò)本身建立于深度學(xué)習(xí)技術(shù)的前提下。具體在交通運(yùn)用中，針對卷積神經(jīng)網(wǎng)絡(luò)有必要明晰其中的多層次網(wǎng)絡(luò)架構(gòu)，從而將其轉(zhuǎn)變成矩陣表示的完整神經(jīng)網(wǎng)絡(luò)。在此基礎(chǔ)上，對于相應(yīng)的原始圖像予以輸入處理，然后矩陣系統(tǒng)即可憑借當(dāng)前輸入的各項信息來判定整個矩陣形態(tài)，從而選擇不同顏色來表示各個相應(yīng)的圖片信息。針對輸出類別而言，對此設(shè)置了1000的默認(rèn)值。在現(xiàn)實運(yùn)用中，如果能靈活運(yùn)用上述圖像處理，則可以運(yùn)用圖片預(yù)測的方式來實現(xiàn)針對雨雪天氣的精確預(yù)測。

本文中炮孔采用炮泥進(jìn)行堵塞，選用土壤和泡沫材料模型（*MAT_SOLI_AND_FOAM）。相關(guān)參數(shù)見表1至表4。

表1 炸藥參數(shù)

表2 炸藥狀態(tài)方程參數(shù)

表3 試塊RHT材料模型的計算參數(shù)

表4 炮泥材料模型的計算參數(shù)

1.3 爆破過程

在后處理軟件（LS-PREPOST）中進(jìn)行分析，當(dāng)損傷值達(dá)到0.7時，巖石發(fā)生破壞，因此，模擬模型損傷值下線為 0.7，為了模擬爆炸整個過程，以孔底抵抗線為4 cm試驗為例，模型內(nèi)部損傷云圖演化過程如圖3所示。同時，為方便觀看漏斗形狀，將模型沿XZ面旋轉(zhuǎn)90°。

圖3 模型失效演化云圖

由圖3可知，在數(shù)值模擬過程中，藥包在10 μs時完全響應(yīng)，爆炸應(yīng)力波繼續(xù)往模型四周傳遞，藥包周圍發(fā)生破壞；10～40 μs時，損傷范圍逐漸沿著炮孔向自由面發(fā)展，模型內(nèi)部形成破碎圈；40～60 μs時，損傷范圍發(fā)展至模型下表面（自由面），爆破漏斗初步形成，漏斗半徑較小，有待進(jìn)一步擴(kuò)展；60～90 μs時，該時期為漏斗半徑擴(kuò)展時期，應(yīng)力波的反射拉伸導(dǎo)致爆生氣體膨脹快速排出孔外，致使爆破漏斗半徑進(jìn)一步擴(kuò)大并達(dá)到峰值。90～100 μs時，很明顯可以看出，該時期漏斗直徑基本沒有變化，模型表面的損傷區(qū)域增大是由于模擬過程中應(yīng)力波在模型中無法完全消耗，導(dǎo)致應(yīng)力波在模型面不斷發(fā)生反射，導(dǎo)致?lián)p傷區(qū)域較大。綜上，將100 μs看作試驗中爆炸過程的終止時間點。

1.4 不同抵抗線下爆破效果分析

通過對堵塞長度為2 cm，孔底抵抗線大小為3 cm、4 cm、5 cm、6 cm和7 cm時進(jìn)行數(shù)值模擬計算，并對其形成的爆破漏斗進(jìn)行分析。圖4中依次為3～7 cm堵塞長度下等效爆破漏斗模型。

由圖4可知，在孔底抵抗線為7 cm時，爆炸沒有形成爆破漏斗，模型表面基本沒有破壞，而孔底抵抗線為4 cm時，爆炸形成的爆破漏斗半徑最大，爆破漏斗的半徑隨孔底抵抗線的增大先增大后減小。以上結(jié)果說明，抵抗線能改變藥包爆炸能量的作用范圍，適當(dāng)?shù)牡挚咕€能有效減小能量在炮孔中的流失，最大化利用爆炸能量，進(jìn)而改善爆破效果。為定量描述爆破效果與孔底抵抗線之間的關(guān)系，將各方案形成的等效爆破漏斗半徑與等效漏斗體積進(jìn)行統(tǒng)計（D0＞0.7的單元），如圖4所示。

圖4 不同孔底抵抗線下爆破漏斗損傷模型

由圖5可知，孔底抵抗線較小時，大量的爆生氣體快速的從下部孔口沖出，導(dǎo)致炸藥能量作用于巖石的比例減小，造成了炸藥能量的損失浪費(fèi)。隨著抵抗線的增大，爆破漏斗的半徑與體積相應(yīng)的增大，但當(dāng)?shù)挚咕€大于4 cm時，爆破漏斗的半徑出現(xiàn)拐點，當(dāng)?shù)挚咕€大于5 cm時，爆破漏斗的體積也開始出現(xiàn)下滑。說明抵抗線的增大提高了藥包的作用深度，合理的抵抗線能夠有效改善爆破效果。

圖5 不同抵抗線下爆破漏斗半徑與損傷體積

1.5 不同測點應(yīng)變對比

在數(shù)值模擬中，各測點按照模型試驗中的應(yīng)變監(jiān)測點布置，提取各監(jiān)測點的應(yīng)變進(jìn)行分析，如圖6、圖7所示。

圖6 測點布置

由圖7可知，各監(jiān)測點在2.5 ms時，有效塑性應(yīng)變值趨向穩(wěn)定，監(jiān)測點 A、B、C的有效塑性應(yīng)變值與初始應(yīng)變率也較高，遠(yuǎn)遠(yuǎn)高于其他監(jiān)測點。監(jiān)測點D、E、F的應(yīng)變值變化較小，且時間短暫，其有效應(yīng)變數(shù)值大小排序為：A＞B＞C＞F＞E＞D，爆破出現(xiàn)漏斗，炸藥主體能量沿炮孔向下表面?zhèn)鬟f，模型下表面的應(yīng)變值遠(yuǎn)遠(yuǎn)高于其他表面質(zhì)點的應(yīng)變值。研究表明，應(yīng)變值與初始應(yīng)變率隨著炸藥的距離增大而減小。

圖7 測點有效塑性應(yīng)變

2 工程實例

由模擬結(jié)果可知，孔底抵抗線抵抗線為 5 cm時，其爆破效果較好，為了更好地驗證模擬成果，通過相似試驗，相似系數(shù)為20，現(xiàn)場試驗對應(yīng)的孔底抵抗線為1 m，其裝藥結(jié)構(gòu)如圖1（a）所示，采場爆破效果如圖8所示。圖8（a）中，采場上部硐室頂板錨網(wǎng)沒有出現(xiàn)明顯的破壞，圖8（b）中，下部硐室爆堆整體塊度較為均勻，達(dá)到生產(chǎn)過程中塊度破碎要求。試驗表明，該孔底抵抗線條件下能較好地改善現(xiàn)場爆破效果。

圖8 現(xiàn)場試驗爆破效果示意

3 結(jié)論

（1）通過數(shù)值模擬相似模擬試驗可知，隨著孔底抵抗線由3 cm增大至7 cm，模擬結(jié)果表明，其等效漏斗半徑與體積均呈先增大后減小，其中孔底抵抗線為 5 cm，等效漏斗體積最大，爆破效果最好。

（2）通過分析自由面上各監(jiān)測點有效塑性應(yīng)變可知，應(yīng)變值與初始應(yīng)變率隨著炸藥距離的增大而減小。

（3）將研究成果應(yīng)用于礦山采場爆破，現(xiàn)場爆破效果能較好地滿足生產(chǎn)需求，降低沖孔，改善爆破塊度，表明通過調(diào)整 ⅤCR爆破孔底抵抗線來優(yōu)化爆破效果具有一定現(xiàn)場指導(dǎo)意義。