李小武，周　銘，胡　韜，代雨昊，何皇兵（昆明有色冶金設(shè)計研究院股份公司，云南　昆明　650051）

礦山臺階爆破震動速度預(yù)測方法的研究

李小武，周銘，胡韜，代雨昊，何皇兵
（昆明有色冶金設(shè)計研究院股份公司，云南昆明650051）

摘要：針對某露天礦山采場附近居民房屋出現(xiàn)裂縫的問題，進行礦山居民房屋現(xiàn)場測震，并現(xiàn)場實驗檢驗了傳統(tǒng)薩道夫斯基線性回歸方法與BP神經(jīng)網(wǎng)絡(luò)預(yù)測震動速度的兩種方法的可靠性，采用BP神經(jīng)網(wǎng)絡(luò)方法進行爆破震動質(zhì)點速度的預(yù)測，再通過調(diào)整爆破設(shè)計降低爆破震動，保障了房屋的安全。

關(guān)鍵詞：礦山爆破；爆破震動；神經(jīng)網(wǎng)絡(luò)預(yù)測；房屋安全

0　引　言

爆破技術(shù)在采礦、鐵路及國防建設(shè)等方面有著廣泛地應(yīng)用，爆破工作不可避免的會產(chǎn)生一些負面效應(yīng)，如爆破震動、沖擊波、飛石、噪聲等，而爆破震動是爆破施工作業(yè)中最嚴重的負面效應(yīng)。

減小爆破施工負面效應(yīng)的最重要方法就是預(yù)測爆破震動。實測結(jié)果表明：爆破震動強度與炸藥量、爆心距、巖土性質(zhì)以及場地條件等因素密切相關(guān)。國內(nèi)外采用預(yù)測爆破震動速度的公式不盡相同，見表1［1］。

表1　世界各國對爆破震動質(zhì)點振速峰值預(yù)測公式Tab.1　Worldwide prediction formula of the vibration velocity peak of blasting vibration particles

人工神經(jīng)網(wǎng)絡(luò)的出現(xiàn)，國內(nèi)外許多學(xué)者采用神經(jīng)網(wǎng)絡(luò)對爆破震動速度進行預(yù)測，在眾多神經(jīng)網(wǎng)絡(luò)的算法中，BP人工神經(jīng)網(wǎng)絡(luò)在爆破震動峰值振速預(yù)測中應(yīng)用最多，Mohamed、Amnieh、徐全軍、沈蔚等采用BP人工神經(jīng)網(wǎng)絡(luò)進行爆破震動峰值預(yù)報［3－8］。

1　工程背景

某礦山在開采過程中，與采場周邊的村民發(fā)生了比較突出的矛盾之一是采場臺階爆破作業(yè)對周邊的民用建筑房屋的震動損害與破壞，造成企業(yè)與村民關(guān)系緊張，給礦山的建設(shè)與生產(chǎn)帶來了一定的影響。

根據(jù)采場目前爆破設(shè)計的情況，計算得到的爆破震動最大影響距離為113 m，爆破飛石影響為115 m，而爆破確定的安全界線是爆破作業(yè)區(qū)200 m以外，因此從爆破震動的計算影響范圍來講爆破作業(yè)對房屋應(yīng)該是沒有影響的。針對村民多次反映的房屋墻面裂縫等問題，見圖1。礦山進行調(diào)查后認為主要的原因是：采礦場周圍的自建房屋結(jié)構(gòu)簡單、房齡偏長和存在隱性的裂縫，抗震等級難以保證；同時，采場爆破加劇了房屋的損壞。為了緩和企業(yè)與居民的關(guān)系，不至于影響礦山生產(chǎn)，在不考慮對房屋進行搬遷的情況下，尋求降震的方法。

圖1　房屋裂縫情況Fig.1　Cracks of house

2　降震技術(shù)措施

震動速度的影響因素也就是降震技術(shù)中需要處理的因素。根據(jù)國內(nèi)外許多爆破工作者及科研工作者的研究與總結(jié)，降震技術(shù)措施大致有以下幾種方法：

1）爆源控制降震

（1）減小單段最大起爆藥量

由薩氏公式分析，可知單段最大起爆藥量Q 是3個影響質(zhì)點速度V的因素中可以進行調(diào)控的因素，單段最大起爆藥量越大，質(zhì)點震動速度越大。在礦山實際生產(chǎn)中，可以進行孔內(nèi)微差起爆的方法，來減小單段最大起爆藥量。大量試驗研究表明，產(chǎn)生地震效應(yīng)跟爆破炮孔之間的微差時間有著密切的關(guān)系。微差時間的確定有3個原則：①使排間延發(fā)時間大于排內(nèi)延發(fā)時間；②選取微差時間應(yīng)使前后起爆的炸藥量產(chǎn)生的地震波互相干擾；③使前后起爆的炸藥量產(chǎn)生的地震波主震相不重疊。研究表明，微差時間一般選取30～50 ms為宜。需要根據(jù)不同的爆破場地的地質(zhì)條件和環(huán)境，通過試驗和長期觀察來確定，這種方案爆破技術(shù)人員在現(xiàn)場比較容易控制與調(diào)整。

[7]列寧：《列寧專題文集·論辯證唯物主義和歷史唯物主義》，北京：人民出版社，2009年，第161頁。

（2）合理的孔網(wǎng)參數(shù)

國內(nèi)外露天礦山臺階爆破孔網(wǎng)參數(shù)多采用“大孔距”爆破參數(shù)，大量試驗研究表明，在保持每個炮孔負擔(dān)面積基本不變的情況下，適當(dāng)?shù)販p小排距而加大孔距，不僅可以使爆破效果得到改善，而且能減少爆破地震波效應(yīng)。工程技術(shù)人員長期在現(xiàn)場負責(zé)，了解實際地質(zhì)情況，才能更好地調(diào)整好孔網(wǎng)參數(shù)。

（3）擁有良好的自由面

爆破試驗研究得知，爆破時越靠近自由面的炮孔，產(chǎn)生的爆破震動越小，因此，爆破作業(yè)中的每個炮孔必須有自由空間以配合微差技術(shù)，從而使后排炮孔爆破產(chǎn)生的壓縮波從這些自由面中反射，從而降低爆破震動。隨著雷管精度的不斷提高，逐孔起爆技術(shù)廣泛地應(yīng)用于礦山臺階爆破，為臺階爆破提供了良好的自由面。

（4）選擇合適的爆破作用指數(shù)

對露天各種大的爆破工程，爆破作用指數(shù)的大小很大程度上影響著爆破震動的強度，在一定的范圍內(nèi)成反比關(guān)系。因此礦山臺階爆破中應(yīng)盡可能獲得最大松動的爆破條件，以減少爆破震動強度。但該礦山由于生產(chǎn)情況，采用強制的壓碴拋擲爆破。

2）截斷傳播路徑進行隔震

采用減震溝減震可以干擾和阻隔爆破地震波的傳播，減震溝對爆破地震波起到反射和繞射，降低了爆破地震波的強度，從而加快了爆破地震波的衰減，達到了控制爆破震動的作用。但該方法操作起來復(fù)雜，礦山技術(shù)上的力量也不足，因此不予采納。

3　質(zhì)點震動速度預(yù)測

對爆破震動進行預(yù)測是采取有效爆破方案的前提。單段最大起爆藥量是控制爆破震動影響的關(guān)鍵。在幾所居民民房里安裝測震儀進行測震，包括爆心距、單段最大起爆藥量、臺階爆破處與測震點高程、爆心距方向線與自由面的夾角和測點處質(zhì)點震動速度，見圖2，共收集了16組數(shù)據(jù)，見表2。根據(jù)建筑物安全允許的震動速度，見表3，民房建筑物的允許震動速度為2.0 cm／s，可知其震動并未超過安全允許值，但由于房屋老化，參照相關(guān)經(jīng)驗適當(dāng)?shù)亟档推浒踩试S值為1.5 cm／s，再設(shè)計爆破方案，以達到降低爆破震動的原因。采用傳統(tǒng)的薩氏公式進行線性回歸與智能BP神經(jīng)網(wǎng)絡(luò)進行爆破震動的預(yù)測。

圖　2　爆心距與自由面的夾角Fig.2　Angle between blasting center distance and free surface

表2　爆破試驗震動速度Tab.2　Blasting vibration velocity

序號　爆心距（R／m）單段最大起爆藥量（Q／kg）臺階爆破處與測震點間的高程（H／m）爆心距與自由面的夾角　α／（°）震動速度V／（cm／s）12　286　580　16　335　1.11 13　295　600　14　231　1.09 14　278　500　7　256　1.05 15　204　680　24　67　1.95 16　289　640　25　117　1.31

表3　國內(nèi)建筑物允許震動速度Tab.3　Domestically permitted building vibration velocity

3.1線性回歸預(yù)測分析

已編入《爆破安全規(guī)程》的薩道夫斯基公式是國內(nèi)計算爆破震動的依據(jù)，對薩氏公式等號兩端取對數(shù)，得到：由于在同一地區(qū)，K與α都是常數(shù)可知，令y＝ lnV，x＝lnQ－lnR，b＝lnK，（1）式簡化為y＝b＋ ax，即lnV與－lnR具有線性關(guān)系。由線性回歸計算得到K＝57.2797，α＝1.1037。代入薩氏公式得到：

由表3的允許安全震動速度V＝1.5 cm／s和房屋離最近的爆破點的距離為204 m，代入公式（2）中得到單段最大起爆藥量為425.65 kg。

將K與α的值代入薩氏公式中得到預(yù)測爆破震動公式：

將公式（3）計算得到的線性回歸震速值與實際測量值進行誤差分析，以檢驗公式的準確性，如表4所示。誤差大多數(shù)均在10%以內(nèi)，只有第5、6、9、10及11次超過10%。近年來不少學(xué)者對誤差大的原因進行研究表明：這是由于測震點與爆破處的高程所引起的［9］。

表4　線性回歸震速計算值與實際測量值誤差分析Tab.4　Error analysis on the calculated vibration velocity by linear regression and the actual measured value

3.2神經(jīng)網(wǎng)絡(luò)預(yù)測分析

BP人工神經(jīng)網(wǎng)絡(luò)是采用反向傳播算法的一種神經(jīng)網(wǎng)絡(luò)，它是一種監(jiān)督訓(xùn)練多層網(wǎng)絡(luò)的算法，主要原理是通過學(xué)習(xí)，使網(wǎng)絡(luò)的輸出層的誤差平方和達到最小。以爆破震動主要影響的4個因素：單段最大起爆藥量Q、爆心距R、測震點與臺階爆破處的高程H以及爆心距與自由面的夾角α為作為神經(jīng)網(wǎng)絡(luò)輸入；采用3層BP神經(jīng)網(wǎng)絡(luò)，以爆破震動速度V作為神經(jīng)網(wǎng)絡(luò)的輸出。構(gòu)成1個多個輸入單個輸出的神經(jīng)網(wǎng)絡(luò)，隱含層神經(jīng)元個數(shù)采用6個，其拓撲結(jié)構(gòu)見圖3（圖中左側(cè)為輸入元素，右側(cè)為輸出元素）。網(wǎng)絡(luò)訓(xùn)練數(shù)據(jù)采用表2中的16組數(shù)據(jù)，神經(jīng)網(wǎng)絡(luò)訓(xùn)練采用的是trainlm函數(shù)，誤差傳遞采用的是tansig函數(shù)，誤差采用的是mse函數(shù)，訓(xùn)練步驟設(shè)置為4 000步，迭代訓(xùn)練精度達到0.001。用訓(xùn)練好的網(wǎng)絡(luò)就可以進行爆破震動速度的預(yù)測。

圖3　神經(jīng)網(wǎng)絡(luò)拓撲結(jié)構(gòu)Fig.3　Topology of neural network

4　現(xiàn)場檢驗

隨機進行了9次爆破測震試驗用于檢驗2種爆破震動預(yù)測方案的可靠性，試驗數(shù)據(jù)見表5。表中V測1由公式（4）計算得到，V測2由BP神經(jīng)網(wǎng)絡(luò)得到，V實由測震儀實際測量得到。從表5中可以看出線性回歸的最大相對誤差為20%，不超過10%的次數(shù)為6次，而神經(jīng)網(wǎng)絡(luò)的最大相對誤差為15.38%，不超過10%的次數(shù)為7次。神經(jīng)網(wǎng)絡(luò)預(yù)測值與實測值的對比及絕對誤差的對比分別見圖4、圖5。從圖4中可以看出神經(jīng)網(wǎng)絡(luò)預(yù)測值貼近與實際測量值附近，圖5可以神經(jīng)網(wǎng)絡(luò)預(yù)測值的絕對誤差曲線更接近X軸，因此相比較而言，實驗說明了神經(jīng)網(wǎng)絡(luò)在預(yù)測爆破震動強度時更準確一些。

表5　爆破質(zhì)點震動速度誤差對比Tab.5　Error comparison of vibration velocities of blasting particles

圖4　預(yù)測值與實測值的對比Fig.4　Comparison between predicted and measured values

圖5　預(yù)測值與實測值絕對誤差的對比Fig.5　Absolute error comparison between predicted and measured values

5　結(jié)　語

對設(shè)計方案進行爆破質(zhì)點的震動速度的神經(jīng)網(wǎng)絡(luò)方法預(yù)測后以調(diào)整單段最大起爆藥量，并采用孔間微差起爆技術(shù)，在房屋進行修補之后再未出現(xiàn)裂縫情況。

1）薩氏公式線性回歸分析與BP神經(jīng)網(wǎng)絡(luò)均可以用于爆破質(zhì)點速度的預(yù)測，但神經(jīng)網(wǎng)絡(luò)考慮對爆破震動強度的影響因素比薩氏公式更多，采用BP人工神經(jīng)網(wǎng)絡(luò)方法預(yù)測的爆破震動峰值振速更接近實測值，表明了BP人工神經(jīng)網(wǎng)絡(luò)方法預(yù)測爆破震動強度的優(yōu)越性。

2）采用在設(shè)計中提前預(yù)測爆破質(zhì)點的震動速度，再確定最終的爆破設(shè)計方案，為居民房屋的安全提供了保障，改善了企業(yè)與居民關(guān)系，加快了企業(yè)建設(shè)進程。

參考文獻：

［1］申旭鵬.研山鐵礦臺階炮孔爆破地表質(zhì)點振速峰值預(yù)測方法研究［D］.北京：北京科技大學(xué)，2012：2－3.

［2］劉殿中，楊仕春.工程爆破實用手冊［M］.北京：冶金工業(yè)出版社，2004.

［3］沈蔚，徐全軍，季茂榮，等.中深孔爆破震動參數(shù)的BP神經(jīng)網(wǎng)絡(luò)預(yù)報［J］.爆炸與沖擊，2002，22（4）：353 －357.

［4］徐全軍，張慶明，惲壽榕.爆破地震峰值的神經(jīng)網(wǎng)絡(luò)預(yù)報模型［J］.北京理工大學(xué)學(xué)報，1998，18（4）：472 －475.

［5］Mostafa Tantawy Mohamed.Artificial neural network for predic－tion and control of blasting vibrations in Assiut （Egypt）limestone quarry［J］.International Journal of RockMechanics＆MiningSciences，2009，46：426 －431.

［6］Manoj Kh，Singh TN.Prediction of blast induced ground vibrations and frequency in opencast mine：a neural network approach［J］.Journal of Sound and Vibration，2006 （289）：711－725.

［7］楊佑發(fā)，崔波.爆破震動速度峰值的預(yù)測［J］.震動與沖擊，2009，28（10）：195－198.

［8］陳德志，朱瑞賡.基于BP神經(jīng)網(wǎng)絡(luò)的路塹爆破對鄰近民房安全預(yù)測的研究［J］.巖石力學(xué)與工程學(xué)報，2002（2）：51－54.

［9］張?zhí)燔姡R銳，喬寶明，等.爆破震動中薩道夫斯基拓展式的回歸分析［J］.湖南科技大學(xué)學(xué)報：自然科學(xué)版，2012（3）：12－16.

中圖分類號：TD76

文獻標(biāo)識碼：B

文章編號：1004－2660（2015）01－0001－05

收稿日期：2015－01－20.

作者簡介：李小武（1980－），男，湖南人，工程師.主要研究方向：采礦工程.

Study on the Prediction Method of Vibration Velocity in Bench Blasting

LI Xiao－wu，ZHOU Ming，HU Tao，DAI Yu－h(huán)ao，HE Huang－bing
（Kunming Engineering＆Research Institute of Nonferrous Metallurgy Co.Ltd.，Kunming 650051，China）

Abstract：As cracks were observed in the wall of resident houses near the open stope，vibration measurement on site was carried out.The reliability of traditional Sadov’s linear regression and vibration velocity prediction by BP neural network was examined.The blasting vibration particle velocity prediction of by BP neural network was adopted and the blasting vibration was reduced by adjusting the blasting design，thus ensuring the housing security. Key words：mine blasting；blasting vibration；neural network prediction；housing security