張亞賓 甘德清 盧宏建 侯永康

(1.河北聯(lián)合大學礦業(yè)工程學院，河北 唐山063009;2.河北省礦業(yè)開發(fā)與安全工程實驗室，河北 唐山063009)

上世紀60 年代我國從瑞典引入了無底柱分段崩落法，由于該采礦方法高效、高度機械化、結(jié)構(gòu)簡單、生產(chǎn)安全、成本較低等優(yōu)點，在我國地下金屬礦山得到了廣泛應(yīng)用，尤其在地下鐵礦開采中，用該方法所采礦石占地下礦總產(chǎn)量的一半以上［1］。無底柱分段崩落法采用扇形中深孔布置落礦，其爆破質(zhì)量的好壞直接影響到回采指標的高低，爆破質(zhì)量的高低與巖石類別、炸藥類別、炸藥單耗以及爆破參數(shù)密切相關(guān)［2-6］。目前，大部分礦山根據(jù)本礦實際情況對爆破參數(shù)進行優(yōu)化研究，來提高爆破質(zhì)量，對爆破參數(shù)優(yōu)化研究的方法主要包括數(shù)值模擬及進行爆破漏斗物理實驗和建立爆破參數(shù)優(yōu)化經(jīng)濟模型，如人工神經(jīng)網(wǎng)絡(luò)BP 模型等［7-9］。

廟溝鐵礦是河北鋼鐵集團的主采礦山之一，原為露天開采礦山，經(jīng)過20 余a 的開采，已進入凹陷開采階段。礦山目前生產(chǎn)規(guī)模245 萬t/a，露天境界內(nèi)尚有資源儲量1 300 多萬t，露天開采還可以維持5 a。其深部還賦存礦石儲量為4 960 余萬t，需轉(zhuǎn)入地下開采。地下開采設(shè)計年產(chǎn)量為300 萬t，采用中深孔無底柱分段崩落法開采，采礦方法結(jié)構(gòu)參數(shù):進路斷面尺寸寬×高=4.5 m×3.8 m;進路間距為18 m，分段高度為20 m。設(shè)計中為了確定和優(yōu)化炮孔的爆破參數(shù)，采用數(shù)值模擬計算和模糊綜合評價的方法對爆破參數(shù)進行了研究。

1 爆破參數(shù)數(shù)值模擬分析

針對廟溝鐵礦設(shè)計地下開采結(jié)構(gòu)參數(shù)，參考類似礦山經(jīng)驗，選取抵抗線長度為1.8 m。抵抗線長度確定后，扇形炮孔設(shè)計時，主要考慮孔底距的大小對爆破效果的影響。設(shè)計了炮孔孔底距為1.98 m、2.16 m、2.34 m、2.52 m 這4 個中深孔爆破參數(shù)方案，并對其進行了數(shù)值模擬分析。圖1 為孔底距為1.98 m 爆破炮孔設(shè)計圖，表1 為4 個不同孔底距的設(shè)計參數(shù)。

圖1 孔底距為1.98 m 的爆破炮孔設(shè)計Fig.1 Blasting design when hole bottom distance at 1.98 m

表1 4 個方案的設(shè)計參數(shù)Table 1 Design parameters of the four schemes

根據(jù)上述4 種炮孔設(shè)計的炮孔參數(shù)進行數(shù)值模擬。在有限元程序ANSYS 中加載命令流的方式建立實體模型，劃分網(wǎng)格(如圖2 所示)，設(shè)置邊界約束條件和無反射邊界，定義運算時間、步長等參數(shù)，最后輸出計算文件。計算選擇取的炸藥為巖石型粉狀乳化炸藥，其主要計算參數(shù):密度1.26 g/cm3，爆速0.55 cm/μs，爆壓3.43 GPa。礦石參數(shù):密度3.4 g/cm3，彈性模量40 GPa，泊松比0.27，抗拉強度12.3 MPa。

為了減少運算時間，建立二分之一模型。在模型平面上設(shè)置5 個記錄點(5 個單元)，記錄點在2 條進路中垂線上每隔5 m 取1 個單元，作為時間歷程記錄單元，用以記錄該單元有效應(yīng)力場的變化歷程。

以孔底距a=1.98 m 為例，以條形藥包的形式起爆，起爆后不同時刻爆炸應(yīng)力云圖如圖3 所示。

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

圖3 不同時刻爆炸應(yīng)力Fig.3 Explosion stress at different time

通過LS－DYNA 軟件的LSPOST 程序?qū)Y(jié)果文件進行后處理，用有效應(yīng)力來表征介質(zhì)的應(yīng)力特征。從圖4 可得監(jiān)測線上的5 個單元最大有效應(yīng)力，具體數(shù)值見表2 所示。

圖4 監(jiān)測線上有效應(yīng)力隨時間變化曲線Fig.4 The effective stress variation curves with time on monitoring line

表2 記錄單元最大有效應(yīng)力Table 2 The maximum effective stress of the record element

爆炸荷載作用下介質(zhì)所受的應(yīng)力狀態(tài)是非常復雜的，從表2 中記錄的最大有效應(yīng)力值可以看出，各單元最大有效應(yīng)力均遠大于礦石的最大抗拉強度12.3 MPa，炸藥爆炸時能有效破壞炮孔之間的礦巖?？椎拙?.98 m 時，單元最大有效應(yīng)力值為108.2 MPa，孔底距為2.52 m 時，單元最大有效應(yīng)力值為107.9 MPa。

通過對4 種方案的相同單元最大有效應(yīng)力值進行比較可以看出:增加孔底距長度對有效應(yīng)力值的大小影響較小，當孔底距增加時，質(zhì)點最大有效應(yīng)力值會降低，可能會增加大塊率，但都大于巖石屈服強度，都可以將炮孔間的巖石破碎，當孔底距減小時，大塊率可能減小，但增加了炮孔數(shù)，使總裝藥量增加，增加了炸藥成本，所以對于上述4 種爆破參數(shù)設(shè)計需進行綜合評價。

2 爆破參數(shù)綜合評價

采用模糊數(shù)學的方法對上述4 個爆破參數(shù)方案的爆破效果進行綜合評價，根據(jù)礦山的實際情況，用于評價爆破效果的(U)主要有3 個:炸藥成本(U1)、采幅(U2)、適宜塊度率(U3)，可歸納為

爆破效果分為3 個等級，分別為“好，一般，不好”，可歸納為

評價爆破效果的3 個因素都屬于好到不好之間，構(gòu)成一個模糊集合，可用0 ～1 區(qū)間中的某個數(shù)值來表示，這個數(shù)值稱為評價因素的隸屬度，一般用μ 來表示隸屬度的大小，作如下規(guī)定:

先確定各評價因素的隸屬函數(shù)，選擇合適函數(shù)形式，再根據(jù)現(xiàn)場統(tǒng)計資料用待定系數(shù)法求出各函數(shù)的隸屬度，設(shè)待定系數(shù)為b，各因素的隸屬函數(shù)如下。

(1)炸藥成本:

(2)采幅:

(3)適宜塊度率:

表3 所示為4 個炮孔設(shè)計的經(jīng)濟技術(shù)指標。

表3 4 個爆破參數(shù)設(shè)計的經(jīng)濟技術(shù)指標Table 3 Eco-Tech indicators of the four blasting parameters

以2.52 m 為例，根據(jù)爆破設(shè)計的統(tǒng)計資料，分別將其代入上式，得出隸屬度μ 的值，根據(jù)隸屬度的大小對其進行分級。得出爆破效果的評判矩陣可表示為

根據(jù)專家評分，可得權(quán)值A(chǔ) 為

A =[ 0.396 0.193 0. 133] ，

由此可得試驗采場爆破技術(shù)經(jīng)濟指標模糊綜合評判值

B1= A·R = (0.396，0.193，0).

同理，可得孔底距為1.98 m 時，其他3 種炮孔設(shè)計的綜合評判值如下:

B2= (0.133，0，0.396);

當孔底距為2.12 m 時，

B3= (0.133，0.396，0.193);

當孔底距為2.34 m 時，

B4= (0，0.396，0.193).

根據(jù)最大隸屬度原則，可見綜合評判結(jié)果為2.52 m 好，2.12 m 和2.34 m 一般，1.98 m 不好。

綜上所述，經(jīng)過對4 種爆破參數(shù)的綜合評價得出:當孔底距為2.52 m 時，爆破效果最好，所以建議廟溝鐵礦在實際生產(chǎn)中中深孔爆破采用2.52 m 的孔底距。

3 結(jié) 論

(1)通過數(shù)值模擬分析得出，4 個不同孔底距的爆破參數(shù)，炸藥爆炸時均能有效破壞炮孔之間的礦巖，并且增加孔底距對有效應(yīng)力值影響不大，隨著孔底距的增加有小幅度的降低。

(2)針對4 種不同孔底距的爆破參數(shù)，根據(jù)它們的經(jīng)濟技術(shù)參數(shù)，用模糊數(shù)學的方法進行了綜合評價，根據(jù)最大隸屬度原則，經(jīng)對比得出:孔底距2.52 m 時為最優(yōu)爆破參數(shù)。

［1］ 王文杰，任鳳玉.無底柱分段崩落法礦石貧化原因分析［J］. 中國礦業(yè)，2008(3):69-72.

Wang Wenjie，Ren Fengyu. The ore depletion reasons analysis of sub-level caving method［J］. China Mining Magzine，2008(3):69-72.

［2］ 綦曉磊，蘇永定.無底柱分段崩落法在謙比西銅礦的應(yīng)用［J］.金屬礦山，2013(5):31-33.

Qi Xiaolei，Su Yongding. Application of none-pillar sublevel caving method in Chambisi Copper Mine［J］.Metal Mine，2013(5):31-33.

［3］ 任鳳玉，何榮興，Thiemo Amadou Mouctar Sow，等. 北洺河鐵礦改進裝藥結(jié)構(gòu)試驗研究［J］.金屬礦山，2011(2):9-12.

Ren Fengyu，He Rongxing，Thiemo Amadou Mouctar Sow，et al.Experimental study on improving charge structure of blasting hole in Beiminghe Iron Mine［J］.Metal Mine，2011(2):9-12.

［4］ 任鳳玉，周宗紅，穆太升，等. 夏甸金礦中深孔爆破參數(shù)優(yōu)化研究［J］.金屬礦山，2005(11):4-6.

Ren Fenyu，Zhou Zonghong，Mu Taisheng，et al.Parameter optimization for blasting with medium-to-deep holes［J］. Metal Mine，2005(11):4-6.

［5］ 楚立申.中深孔爆破在大冶鐵礦的應(yīng)用［J］. 礦業(yè)工程，2008，6(2):39-40.

Chu Lishen.Application of medium-deep blasting holes in Daye Iron Mine［J］.Mining Engineering，2008，6(2):39-40.

［6］ 吳 亮，張發(fā)勇. 淺析礦山開采中爆破參數(shù)優(yōu)化與災害控制［J］.現(xiàn)代礦業(yè)，2009(8):19-21.

Wu Liang，Zhang Fayong.Preliminary analysis of blasting parameter optimization and hazard control in mining［J］.Modern Mining，2009(8):19-21.

［7］ 周 馭，周文海，樓曉明.青海德爾尼銅礦爆破效果評估及優(yōu)化［J］.金屬礦山，2015(1):43-46.

Zhou Yu，Zhou Wenhai，Lou Xiaoming.Blasting effect evaluation and optimization of Deerni Copper Mine in Qinghai［J］. Metal Mine，2015(1):43-46.

［8］ 葉海旺，朱瑞賡.基于模糊綜合評判的巖石分級系統(tǒng)研究［J］.武漢理工大學學報，2003(5):47-49.

Ye Haiwang，Zhu Ruigeng. Study on classification system based Fuzzy integration judge［J］.Journal of Wuhan University of Technology，2003(5):47-49.

［9］ 祝文化，朱瑞賡，夏元友.堆石壩料開采爆破參數(shù)優(yōu)化的神經(jīng)網(wǎng)絡(luò)法研究［J］.工程爆破，2001(1):20-23.

Zhu Wenhua，Zhu Ruigeng，Xia Yuanyou. Research on optimization of blasting parameters for mining of dam rock-fill material by neural network method［J］.Engineering Blasting，2001(1):20-23.

［10］ 程 峰，王杰光，靳麗輝.模糊數(shù)學理論在金屬礦山安全評估中的應(yīng)用［J］.金屬礦山，2007(3):77-80.

Cheng Feng，Wang Jieguang，Jin Lihui.Application of fuzzy mathematics theory in mine safety assessment［J］.Metal Mine，2007(3):77-80.