王帷先，張瑞新，孫健東

(華北科技學(xué)院，北京 東燕郊 065201)

0 引言

目前，黑岱溝露天煤炭開(kāi)采主要采用拋擲爆破——拉斗鏟倒堆工藝綜合剝離高臺(tái)階巖層。其中，爆破拋擲形成的爆堆是關(guān)鍵分析指標(biāo)，它直接決定拉斗鏟站立平臺(tái)構(gòu)建的二次倒堆量，決定后續(xù)倒堆周期；而爆堆形態(tài)效果主要受爆破設(shè)計(jì)參數(shù)控制，因此也可用于評(píng)價(jià)并優(yōu)化設(shè)計(jì)方案。目前，有效拋擲率是露天拋擲爆破主要評(píng)價(jià)指標(biāo)，有關(guān)實(shí)驗(yàn)表明巖土爆破拋擲率在25%～50%之間，而近年來(lái)現(xiàn)場(chǎng)爆破有效拋擲率基本達(dá)到27%以上水平，實(shí)際平均穩(wěn)定在33%左右，還有較大改善空間。因此，研究爆堆形態(tài)對(duì)提升有效拋擲率、優(yōu)化爆破設(shè)計(jì)參數(shù)具有重要意義。

目前，大量研究圍繞拋擲爆破參數(shù)對(duì)爆堆形態(tài)的預(yù)測(cè)展開(kāi)，大部分研究通過(guò)Weibull模型描述爆堆形態(tài)，利用不同參數(shù)組合對(duì)形態(tài)進(jìn)行預(yù)測(cè)。如，通過(guò)拋擲爆破參數(shù)訓(xùn)練模型[1]，預(yù)測(cè)爆堆形態(tài)，并對(duì)爆堆加以分類(lèi)[2]，驗(yàn)證預(yù)測(cè)效果；通過(guò)開(kāi)發(fā)拋擲爆破設(shè)計(jì)系統(tǒng)[3,4]，實(shí)現(xiàn)自動(dòng)化爆破模擬仿真；通過(guò)拉斗鏟作業(yè)設(shè)計(jì)特性[5]，利用不同設(shè)計(jì)參數(shù)進(jìn)行形態(tài)預(yù)測(cè)。其他，通過(guò)重點(diǎn)研究網(wǎng)絡(luò)模型的拓?fù)浣Y(jié)構(gòu)參數(shù)[6]，改進(jìn)模型算法，優(yōu)化爆堆預(yù)測(cè)結(jié)果；通過(guò)利用極限學(xué)習(xí)機(jī)(ELM)[7]或引入遺傳算法[8]，擴(kuò)容小樣本數(shù)據(jù)實(shí)現(xiàn)爆堆形態(tài)預(yù)測(cè)；通過(guò)對(duì)比不同模型的爆堆預(yù)測(cè)結(jié)果[9]，論證模型適用性和準(zhǔn)確性。

上述研究基本實(shí)現(xiàn)以理論確定爆破參數(shù)并擬合預(yù)測(cè)爆堆形態(tài)，但爆破參數(shù)眾多，需要進(jìn)行大量爆堆形態(tài)預(yù)測(cè)模擬，而實(shí)際現(xiàn)場(chǎng)生產(chǎn)，更需要直接有效的參數(shù)設(shè)計(jì)方案。因此，本文提出以大量爆堆形態(tài)數(shù)據(jù)為基礎(chǔ)，建立從爆堆形態(tài)到爆破設(shè)計(jì)參數(shù)的關(guān)聯(lián)網(wǎng)絡(luò)預(yù)測(cè)模型，直接利用現(xiàn)場(chǎng)有效爆堆形態(tài)，反向預(yù)測(cè)爆破設(shè)計(jì)較優(yōu)參數(shù)值，優(yōu)化設(shè)計(jì)方案，形成現(xiàn)場(chǎng)設(shè)計(jì)的參考模板。

1 神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型

BP神經(jīng)網(wǎng)絡(luò)拓?fù)淠Ｐ褪墙鉀Q各類(lèi)數(shù)據(jù)之間合理關(guān)聯(lián)的分析模型，數(shù)據(jù)集通過(guò)各層不同節(jié)點(diǎn)進(jìn)行特征儲(chǔ)存和傳遞，利用誤差下降不斷修正權(quán)值，最后確定某類(lèi)現(xiàn)象與結(jié)果關(guān)系函數(shù)，從而實(shí)現(xiàn)預(yù)測(cè)。

1.1 模型的基礎(chǔ)單位

神經(jīng)元節(jié)點(diǎn)作為最基礎(chǔ)結(jié)構(gòu)體，具備了三個(gè)最基本也是最重要功能：加權(quán)、求和和轉(zhuǎn)移，各節(jié)點(diǎn)通過(guò)數(shù)量上的相互鏈接、擴(kuò)容組成功能完善的網(wǎng)絡(luò)結(jié)構(gòu)。圖1展示了第i個(gè)基本節(jié)點(diǎn)，其中，x1、x2…xj分別代表來(lái)自神經(jīng)元1、2 …j的輸入；w1、w2…wj則分別表示神經(jīng)元1、2 …j與第i個(gè)神經(jīng)元的連接強(qiáng)度，即權(quán)值；bi為閾值(若數(shù)據(jù)信號(hào)達(dá)到閾值將被激活)；最后數(shù)據(jù)信號(hào)通過(guò)激活函數(shù)f(·)得出y值，yk為第i個(gè)元節(jié)點(diǎn)的輸出。

圖1 單個(gè)神經(jīng)元功能

故對(duì)于該神經(jīng)元i的可推導(dǎo)出輸入后的值為：

(1)

式中，w為權(quán)重值；b為閾值。X= [x1，x2…xj]T，Wj= [wj1，wj2…wjn]。

輸入值Sj通過(guò)激活函數(shù)f(·)，得到第j個(gè)元節(jié)點(diǎn)輸出yj：

(2)

1.2 模型的傳播計(jì)算

通常在三層網(wǎng)絡(luò)構(gòu)造基礎(chǔ)上進(jìn)行模擬即可實(shí)現(xiàn)大部分復(fù)雜問(wèn)題的分析，而提高網(wǎng)絡(luò)訓(xùn)練精度的重點(diǎn)在于隱含層內(nèi)部神經(jīng)元節(jié)點(diǎn)的數(shù)量設(shè)置。模型主要通過(guò)正向和反向傳播交替循環(huán)，不斷分析誤差和更新權(quán)值，達(dá)到計(jì)算結(jié)果的逼近，如圖2所示。其中，反向傳播只在模型訓(xùn)練環(huán)節(jié)實(shí)現(xiàn)，它根據(jù)正向計(jì)算結(jié)果與預(yù)期值的對(duì)比差距作為反饋信號(hào)，反向逐層傳遞使訓(xùn)練沿著誤差減小的方向修正，不斷迭代各層連接權(quán)重值，最終完成數(shù)據(jù)信息的記憶。當(dāng)模型訓(xùn)練完保存后，將不再進(jìn)行反向傳播，訓(xùn)練和測(cè)試數(shù)據(jù)集的錯(cuò)誤都將降到最低。因此，誤差計(jì)算的反向傳播算法是調(diào)整權(quán)值提高精度的關(guān)鍵環(huán)節(jié)，常利用均方差計(jì)算：

(3)

(4)

圖2 網(wǎng)絡(luò)模型的兩種傳播過(guò)程

式中，ylp為網(wǎng)絡(luò)計(jì)算值；tlp為期望(理想)值；Ep為樣本均方差；E為網(wǎng)絡(luò)總誤差。

拋擲爆破分析涉及參數(shù)多樣、數(shù)據(jù)量大，要求模型具備足夠的分布儲(chǔ)存和并行處理能力，同時(shí)實(shí)現(xiàn)針對(duì)不同參數(shù)類(lèi)型進(jìn)行非線性映射處理，快速擬合并預(yù)測(cè)。BP神經(jīng)網(wǎng)絡(luò)模型滿足該分析要求，因此，搭建三層網(wǎng)絡(luò)結(jié)構(gòu)的爆堆形態(tài)預(yù)測(cè)模型。

2 爆堆形態(tài)及設(shè)計(jì)參數(shù)

為了建立爆堆和設(shè)計(jì)參數(shù)之間的關(guān)聯(lián)映射模型，提高預(yù)測(cè)效率和準(zhǔn)確度，統(tǒng)計(jì)可用于預(yù)測(cè)分析的關(guān)鍵參數(shù)有：臺(tái)階高度、煤層厚度、煤層傾角、炮孔角度、孔徑、孔距、排距、抵抗線、炸藥單耗、炸藥總藥量、爆堆形態(tài)數(shù)據(jù)、有效拋擲率、松散系數(shù)等。

其中，經(jīng)多年實(shí)踐論證，多項(xiàng)參數(shù)已基本達(dá)到最優(yōu)設(shè)計(jì)值，因此可不用于預(yù)測(cè)分析計(jì)算。如：現(xiàn)場(chǎng)采用的牙輪鉆機(jī)，炮孔直徑已統(tǒng)一為310 mm；設(shè)計(jì)的有效爆區(qū)寬度固定為85 m；統(tǒng)一鉆孔和爆區(qū)尺寸后，排距、孔距也可基本計(jì)算確定為7 m排距共計(jì)12排，每排孔距12 m若干個(gè)孔；各炮孔角度基本穩(wěn)定在65°。其他參數(shù)則受多因素影響，數(shù)值上呈現(xiàn)一定離散狀態(tài)，需通過(guò)建立輸入和輸出兩個(gè)部分的參數(shù)數(shù)據(jù)庫(kù)，用于模型分析。

2.1 輸入?yún)?shù)庫(kù)

模型以爆堆形態(tài)為核心，利用爆堆剖面數(shù)據(jù)描述曲線形態(tài)，因此輸入?yún)?shù)為：爆堆曲線形態(tài)數(shù)據(jù)(Hi)。

表1 輸入?yún)?shù)數(shù)據(jù)(部分)

2.2 預(yù)測(cè)庫(kù)建立

模型主要通過(guò)爆堆形態(tài)分析優(yōu)化爆破設(shè)計(jì)方案，因此確定輸出參數(shù)為：煤層厚度(h)、煤層傾角(φ)、裝藥量(Q)、單耗(q1,q2,q3)、有效拋擲率(η)、松散系數(shù)(ζ)。

表2 輸出參數(shù)數(shù)據(jù)(部分)

3 模型訓(xùn)練及預(yù)測(cè)

3.1 標(biāo)準(zhǔn)歸一

在模型分析中，由于數(shù)據(jù)維度多，數(shù)據(jù)值分布多樣，因此必須通過(guò)歸一化處理將輸入和輸出數(shù)值進(jìn)行無(wú)量綱化，使該分析模型中各數(shù)值轉(zhuǎn)換成相對(duì)值關(guān)系，以此縮小量值，使計(jì)算過(guò)程誤差減小。通常采用最大最小值法進(jìn)行標(biāo)準(zhǔn)歸一化處理，見(jiàn)式(5)。

(5)

3.2 網(wǎng)絡(luò)節(jié)點(diǎn)分析

針對(duì)爆堆預(yù)測(cè)模型，以輸入的爆堆曲線參數(shù)維度確定輸入節(jié)點(diǎn)數(shù)，以期望爆破參數(shù)數(shù)據(jù)維度確定輸出節(jié)點(diǎn)數(shù)。而隱含層節(jié)點(diǎn)數(shù)，常通過(guò)經(jīng)驗(yàn)公式計(jì)算參考范圍值，并通過(guò)訓(xùn)練實(shí)驗(yàn)進(jìn)行對(duì)比，最終定量確定最優(yōu)節(jié)點(diǎn)數(shù)。

(6)

式中，n為輸入節(jié)點(diǎn)數(shù)；m為輸出節(jié)點(diǎn)數(shù)；a為1～10之間常數(shù)。

經(jīng)計(jì)算，得出隱含層合理神經(jīng)元節(jié)點(diǎn)數(shù)為6～17個(gè)。在Matlab中直接分別進(jìn)行訓(xùn)練實(shí)驗(yàn)，通過(guò)對(duì)不同節(jié)點(diǎn)數(shù)的訓(xùn)練誤差判斷，見(jiàn)表3，確定最小誤差對(duì)應(yīng)的隱含層節(jié)點(diǎn)數(shù)為12個(gè)。

表3 不同隱含層神經(jīng)元節(jié)點(diǎn)數(shù)的網(wǎng)絡(luò)誤差分析

3.3 預(yù)測(cè)誤差分析

從樣本數(shù)據(jù)庫(kù)中隨機(jī)選取多組參數(shù)數(shù)據(jù)進(jìn)行實(shí)驗(yàn)，通過(guò)對(duì)比分析原始網(wǎng)絡(luò)和訓(xùn)練網(wǎng)絡(luò)的誤差可知，原始網(wǎng)絡(luò)模擬值中有12.5%的參數(shù)值誤差超過(guò)1，其中17N1-39剖面的后兩排炸藥單耗誤差數(shù)值接近7，數(shù)據(jù)差距極大；僅有31.2%的參數(shù)值誤差在0.1以下。訓(xùn)練后的網(wǎng)絡(luò)誤差明顯減小，參數(shù)值誤差水平小于0.1的結(jié)果高達(dá)71.8%，不再出現(xiàn)超過(guò)1的誤差值。對(duì)比說(shuō)明訓(xùn)練后的神經(jīng)網(wǎng)絡(luò)模型預(yù)測(cè)準(zhǔn)確程度提高。

表4 參數(shù)測(cè)試數(shù)據(jù)的真實(shí)值、初始值、仿真值比較

續(xù)表

如圖3所示，藍(lán)線和紅線分別為模型訓(xùn)練前和訓(xùn)練后的誤差分布，可明顯發(fā)現(xiàn)，訓(xùn)練后的誤差基本回落并穩(wěn)定在0.4以內(nèi)。

圖3 訓(xùn)練前、后測(cè)試誤差分布

從測(cè)試數(shù)據(jù)仿真值的單相誤差上分析，在獲得的所有數(shù)據(jù)誤差值中進(jìn)行分析，對(duì)每項(xiàng)指標(biāo)計(jì)算平均誤差值，煤層傾角和后兩排單耗誤差值仍然突出，整體樣本誤差均值在0.0798，而去掉突出誤差指標(biāo)后，整體樣本誤差均值下降到0.0453。通過(guò)結(jié)果分析，模型預(yù)測(cè)值基本處于合理范圍。

表4 測(cè)試樣本參數(shù)誤差

注：整體樣本均差值：0.0798

降噪樣本均差值：0.0453

4 案例應(yīng)用及優(yōu)化

4.1 案例應(yīng)用

根據(jù)露天礦高臺(tái)階的高度分布情況分組，現(xiàn)有三組滿足拉斗鏟倒堆作業(yè)平臺(tái)快速構(gòu)建的具有典型代表性意義的爆堆形態(tài)，見(jiàn)表5。將這三組爆堆形態(tài)用于爆堆形態(tài)預(yù)測(cè)模型，獲得相應(yīng)參數(shù)設(shè)計(jì)值，為現(xiàn)場(chǎng)爆破設(shè)計(jì)方案的優(yōu)化提供參考。

4.2 預(yù)測(cè)結(jié)果分析

經(jīng)模型預(yù)測(cè)計(jì)算，同時(shí)從爆堆樣本庫(kù)中挑選三組對(duì)應(yīng)相似的曲線形態(tài)樣本進(jìn)行對(duì)比，分析預(yù)測(cè)結(jié)果的有效性，見(jiàn)表5。

表5 預(yù)測(cè)數(shù)值分析

續(xù)表

其中，40 m曲線的預(yù)測(cè)效果比其他兩組要好，數(shù)值差異范圍基本處于可接受范圍，精度基本在90%以上。45 m曲線預(yù)測(cè)數(shù)值也相對(duì)匹配，而35 m曲線預(yù)測(cè)值存在部分偏差。三組對(duì)比案例中，后兩排炮孔單耗預(yù)測(cè)值都不理想，煤層傾角由于數(shù)值精確度的限制，利用相對(duì)誤差計(jì)算值進(jìn)行評(píng)價(jià)較為不準(zhǔn)。若排除以上兩類(lèi)參數(shù)影響，預(yù)測(cè)數(shù)值平均誤差降為0.0656水平，效果明顯變好。

表6 預(yù)測(cè)參數(shù)誤差

注：整體樣本均差：0.3014

降噪樣本均差：0.0656

從預(yù)測(cè)誤差的數(shù)值分布圖中也可確定，煤層傾角和后兩排炮孔的誤差較大，其余參數(shù)預(yù)測(cè)誤差分布平穩(wěn)，基本處于較低水平0.5以下。

圖4 預(yù)測(cè)誤差參考

4.3 方案優(yōu)化設(shè)計(jì)

通過(guò)預(yù)測(cè)分析，模型預(yù)測(cè)參數(shù)值基本處于合理范圍。將預(yù)測(cè)參數(shù)值用于方案預(yù)設(shè)計(jì)，獲得三組爆破方案。現(xiàn)場(chǎng)在相同高度范圍內(nèi)對(duì)臺(tái)階爆破作業(yè)設(shè)計(jì)時(shí)，則可參考該三組方案進(jìn)行整體設(shè)計(jì)、調(diào)整。

表7 預(yù)測(cè)設(shè)計(jì)參數(shù)

(1) 35 m臺(tái)階高度的拋擲爆破設(shè)計(jì)方案

該方案用于爆區(qū)臺(tái)階高度為33～37 m范圍內(nèi)的拋擲爆破，平均高度為35 m，爆區(qū)長(zhǎng)度為400 m以上，寬度85 m，孔深基本為36～41 m。

煤層厚度基本在32 m左右，煤層傾角在0°～2°范圍。拋擲爆破孔距12 m，排距6～7 m，共計(jì)12排孔，平均單耗為0.707 kg/m3，單排孔裝藥量在2127 t左右。

由于后兩排緩沖炮孔一般設(shè)計(jì)為銨油低密度炸藥—連續(xù)裝藥結(jié)構(gòu)，故平均單耗為0.520 kg/m3。

爆后爆堆有效拋擲率將處于31.5%的水平，松散系數(shù)為1.11。

(2) 40 m臺(tái)階高度的拋擲爆破設(shè)計(jì)方案

該方案用于爆區(qū)臺(tái)階高度為38 m～42 m范圍內(nèi)的拋擲爆破，平均高度為40 m，爆區(qū)長(zhǎng)度為400 m以上，寬度85 m，孔深基本為42 m～46 m。

煤層厚度基本在31 m左右，煤層傾角在0°～2°范圍。拋擲爆破孔距12 m，排距6～7 m，共計(jì)12排孔，平均單耗為0.767 kg/m3左右，單排孔裝藥量在2688 t左右。

由于后兩排緩沖炮孔一般設(shè)計(jì)為銨油低密度炸藥—連續(xù)裝藥結(jié)構(gòu)，故平均單耗為0.540 kg/m3。

爆后爆堆有效拋擲率應(yīng)處于33.5%的水平，松散系數(shù)為1.22。

(3) 45 m臺(tái)階高度的拋擲爆破設(shè)計(jì)方案

該方案用于爆區(qū)臺(tái)階高度為43～47 m范圍內(nèi)的拋擲爆破，平均高度為45 m，爆區(qū)長(zhǎng)度為400 m以上，寬度85 m，孔深基本為47～52 m。

煤層厚度基本在32 m左右，煤層傾角在0°～2°范圍。拋擲爆破孔距12 m，排距6～7 m，共計(jì)12排孔，平均單耗為0.790 kg/m3左右，單排孔裝藥量在2895 t左右。

由于后兩排緩沖炮孔一般設(shè)計(jì)為銨油低密度炸藥—連續(xù)裝藥結(jié)構(gòu)，故平均單耗為0.520 kg/m3。

爆后爆堆有效拋擲率應(yīng)處于32.2%的水平，松散系數(shù)為1.18。

5 結(jié)論

(1) BP神經(jīng)網(wǎng)絡(luò)模型具有較強(qiáng)的樣本數(shù)據(jù)記憶和分析能力，可實(shí)現(xiàn)大量參數(shù)樣本的計(jì)算，形成有效映射關(guān)系。在構(gòu)建拋擲爆破預(yù)測(cè)模型中，對(duì)爆堆數(shù)據(jù)的計(jì)算結(jié)果基本合理，能夠?qū)崿F(xiàn)爆破參數(shù)的預(yù)測(cè)。

(2) 由于爆堆形態(tài)直接決定拉斗鏟倒堆效率，以典型爆堆形態(tài)為直接研究對(duì)象，既能考慮到后續(xù)拉斗鏟倒堆工作量的控制，又能預(yù)測(cè)指導(dǎo)前項(xiàng)拋擲爆破基本參數(shù)的設(shè)計(jì)。

(3) 通過(guò)預(yù)測(cè)結(jié)果，可提前為新爆區(qū)的爆破設(shè)計(jì)及倒堆作業(yè)提供優(yōu)化數(shù)值參考。由此對(duì)原有設(shè)計(jì)方案進(jìn)行更新和優(yōu)化，在方案的不斷實(shí)踐檢驗(yàn)和迭代中，將更好提高有效拋擲率。