微差爆破試驗及爆破振動能量的小波包分析

何 理 鐘冬望 劉建程 朱 寬 黃小武

(1.武漢科技大學資源與環(huán)境工程學院，湖北 武漢430081;2.武漢科技大學理學院，湖北武漢430065;3.冶金工業(yè)過程系統(tǒng)科學湖北省重點實驗室，湖北武漢430065)

通過對爆破振動速度時程進行監(jiān)測與多角度分析，有利更深層次研究爆破地震波傳播理論、危害機制，優(yōu)化爆破參數(shù)，改善爆破效果，最終有效控制爆破地震危害。眾所周知，源于爆破地震波自身的復雜性、瞬時性，以及傳播介質(zhì)和周圍環(huán)境的多變性，導致爆破振動信號具有短時、突變的特點，是一種典型的非平穩(wěn)隨機信號［1-2］。針對其非平穩(wěn)的特點，許多學者近年將小波變換與小波包分析技術(shù)引入到爆破振動信號的分析中，其應(yīng)用價值已得到了充分的體現(xiàn)。徐國元等［3］通過小波變換技術(shù)從各頻帶能量分布的角度提出了基于小波變換的爆破地震安全能量分析法;趙明生等［4］結(jié)合AOK時頻分布與小波分析，研究得出，多段微差爆破中，隨著段數(shù)的增加爆破振動信號質(zhì)點振動速度峰值與主頻存在收斂性;凌同華等［5］基于小波分析的時－能密度法成功識別出了微差爆破中各段實際起爆間隔時間。

大量的工程事例揭示，在爆破地震強度超過甚至遠大于國家安全允許值的情況下，受控對象安全并未受到威脅;反而，在有些爆破地震強度未超過允許值時，受控對象卻受到了一定程度上的損害［6-7］。由此可見，在爆破地震安全評價時，采用單一強度因子作為爆破振動安全判據(jù)的合理性與適用性值得商榷。國內(nèi)外大量研究均已表明，爆破地震主要通過振動強度、振動主頻以及振動持續(xù)時間3個要素對結(jié)構(gòu)物造成危害。而基于小波包分析的各頻帶能量能夠同時涵蓋振動強度、振動主頻和振動持時3個要素。凌同華等［8］基于小波包分析對單段爆破振動信號的能量分布特征進行了研究，得出單段爆破中爆破振動信號成分主要以高頻(39－156 HZ)成分為主;徐國元、中國生等［3］通過小波包分析得到了地震信號各頻帶能量分布特征，研究得出基于爆破振動能量的爆破地震安全判據(jù)比用單一強度因子作為判據(jù)更加合理;在此基礎(chǔ)上劉敦文等［9］又利用小波包良好的時頻局部化特性，從爆破振動信號主振頻率所在頻帶能量和受控結(jié)構(gòu)自振頻率所在頻帶能量角度出發(fā)，提出了以爆破危害評定參數(shù)α作為爆破安全判據(jù)的新方法;趙明生、梁開水等［10］研究了段藥量對爆破振動信號不同頻帶能量分布特征的影響，為從能量的角度進行抗振、降振研究提供了分析基礎(chǔ)。綜上，針對爆破振動信號頻帶能量分布的研究，不僅為廣大爆破工作者開拓了一種爆破振動信號分析的新思路，而且為進一步改進和完善爆破地震安全判據(jù)積累了經(jīng)驗。因此，本研究基于前人研究成果，采用室內(nèi)混凝土邊坡相似微差爆破試驗，選取爆破振動監(jiān)測典型信號，結(jié)合小波包技術(shù)，利用MATLAB 7.0，編制相應(yīng)計算程序?qū)吰挛⒉畋普駝有?yīng)進行分析，研究振動信號能量和峰值速度隨距離的變化規(guī)律，分析減震溝對振動信號能量與峰值速度的削弱程度，以及對振動信號能量分布的影響，總結(jié)規(guī)律，為研究和控制邊坡工程爆破地震危害、實現(xiàn)爆破振動主動控制提供依據(jù)。

1 小波包分解與爆破振動信號頻帶能量表征

1.1 小波包分解特點

小波分析是把信號占有的頻帶一分為二，得到第一級低頻帶和高頻帶，下一級將低頻部分進一步分解成低頻帶和高頻帶，如此不停的將分解進行下去。從小波分解的結(jié)構(gòu)可以看出，不同級次具有不同的分辨率，隨著級次的升高，小波變換的頻率分辨率逐漸升高，同時，時間分辨率卻逐漸降低。針對小波分解時頻分辨率不能同時升高的不足，小波包分解將每一級次的所有濾波器子帶均一分為二，并如此傳遞至下一級。小波包分解能根據(jù)信號特征和分析要求自適應(yīng)地選擇相應(yīng)頻帶與信號頻譜相匹配。小波包分解是在小波分解的基礎(chǔ)上被提出的，具有嚴密的數(shù)學理論，是一種比小波分解更為精細的分解方法［11-12］。

1.2 爆破振動信號的小波包分解

將爆破振動信號進行小波包分解時，分解的層數(shù)應(yīng)視具體信號的爆破參數(shù)及振動測試儀的采樣率而定。通常工程爆破振動信號的頻率一般集中在5～300 Hz范圍內(nèi)，但都是在大藥量、爆破遠區(qū)情況下測得，考慮到試驗條件下的小尺度、單發(fā)雷管模擬單段藥量、混凝土材料堅硬完整等因素，都可能導致地震波呈現(xiàn)短持時、高頻率的特性。因此，基于采樣定理，將信號的采樣率設(shè)為8 kHz，則其奈奎斯特(Nyquist)頻率為4 kHz。將信號分解至10層，對應(yīng)的最低頻帶為0～3.91 Hz，根據(jù)小波包算法，采用二進尺度變換，進行n層分解后第j個頻帶重構(gòu)信號Sn，j對應(yīng)的頻率范圍為

式中，j=0，1，2….2n－ 1;f為分析頻率，f=4 kHz。

aa

1.3 信號各頻帶的能量表征

將被分析信號分解至10層，設(shè)S10，j對應(yīng)的能量為 E10，j，則有

式中，xj，i(j=0，1，2，…，210－ 1;i=1，2，…，k，k為信號的離散采樣點數(shù))為重構(gòu)信號S10，j各離散點對應(yīng)的幅值。

設(shè)分析信號的總能量為E0，則有

于是各頻帶能量占信號總能量的百分比為

式中，j=0，1，2….210－1。

通過式(1)～式(3)可計算得到信號各頻帶能量分布。

2 爆破振動測試

爆破振動測試試驗在實驗室混凝土邊坡上進行，根據(jù)爆破相似準則，取合適的幾何相似比CL，確定模型尺寸參數(shù);模型試件的制作采用425#硅酸鹽水泥和篩選后的細砂澆注而成，確定其質(zhì)量配比為0.5∶1∶2(水∶水泥∶砂子)，養(yǎng)護28 d;在邊坡模型底部加裝模板預制減震溝;為改善模型邊界條件的相似性，在模型兩邊加裝兩塊10 mm厚鋼板，鋼板與混凝土之間接觸通過安裝3 mm厚橡膠墊完成，兩塊鋼板通過4根鐵絲連接，通過扭緊鐵絲給模型施加預應(yīng)力。采用單發(fā)數(shù)碼電子雷管模擬單段炸藥爆炸，通過串聯(lián)相應(yīng)數(shù)量雷管并設(shè)置不同的起爆時間來模擬微差爆破，炮孔深度為15 cm，直徑5 mm，炮孔堵塞8 cm。振動測試儀器采用由成都中科測控研發(fā)的TC4850爆破振動測試儀，監(jiān)測點均分布在坡底同一高程，從實測數(shù)據(jù)中選取8條典型振動信號進行小波包分析，其中信號s7、s8為測點在減震溝內(nèi)側(cè)靠近坡腳處測得，減震溝寬度為20 mm，深度為15 cm，各信號的爆破條件及測點位置如表1所示。邊坡模型與尺寸如圖1所示。

表1 監(jiān)測點的爆破條件Table 1 Blasting condition of monitoring points

圖1 邊坡模型與尺寸(單位:cm)Fig.1 Slope model and its size chart

各振動信號對應(yīng)的豎直向振動速度時程曲線如圖2所示。

3 振動信號的小波包分析

3.1 小波基的選擇

圖2 各測點信號豎直向振動速度時程曲線Fig.2 Time curves about vertical vibration velocity of each measuring point

在對振動信號進行小波包分析時，所選小波基的匹配程度直接影響到最終的分析結(jié)果。Daubechies小波系列憑借其良好的緊支撐性、光滑性以及近似對稱性［13-14］，近年被廣泛的應(yīng)用于爆破振動信號分析中。Daubechies小波系列按正整數(shù)N具有不同系列(dbN)，目前db5和db8小波系列運用最為廣泛，本研究采用db8對振動信號進行小波包分析。

3.2 小波包分析

運用MATLAB 7.0對圖2中的爆破振動信號進行深度為10層的小波包分析。根據(jù)式(1)～式(3)編寫相應(yīng)計算程序，計算得到各信號總能量見表2。

表2 各信號總能量Table 2 The total energy of each signal

信號各頻帶能量分布如圖3所示。

各信號在能量集中頻帶的分布如表3所示。

4 數(shù)據(jù)分析

表1中，對比信號s4、s5、s6，在爆破參數(shù)相同的情況下，隨著爆心距的增加，各測點峰值振速逐漸降低，比較速度降低的幅度，可以得到隨著距離的增加，峰值速度衰減逐漸變緩，衰減趨勢呈明顯的非線性。

圖3 爆破振動信號各頻帶(F)能量(E)分布Fig.3 Energy distribution(E)of blasting vibration signals in each frequency band(F)

表3 信號在能量集中頻帶的分布Table 3 Signals distribution in energy concentrated bands

由圖2中s1、s2、s3、s5信號振動時程可以看出，當爆心距、單段最大藥量、總藥量相同的情況下，隨著段間間隔時間的增加，振動持續(xù)時間逐漸增長，分段振波互相干擾、疊加程度依次減弱，一定范圍內(nèi)適當增加延期時間可有效避免各段振波疊加，降低前一段振波波峰與后一段振波波峰相遇的概率，從而達到降低峰值振速的目的。

分別對比表1中信號s1與s7以及s4與s8，可以看到通過在測點與爆源之間設(shè)置一定尺寸的減震溝，可有效降低峰值振動速度，本研究試驗條件下，減震溝對速度的減震率達到了40%以上。

由表2中各信號總能量可以看出，各測點能量總體隨傳播距離的增加逐漸減小，且微差間隔時間對信號總能量影響較大;對比信號s4、s5、s6總能量，可以得到能量隨傳播距離增加同樣呈現(xiàn)出明顯的非線性趨勢衰減;對比信號s3與信號s7總能量可知，峰值速度大的信號對應(yīng)的能量不一定大，因為振動信號能量不僅僅是振動速度的表征，它是速度、頻率、振動持時綜合作用的結(jié)果。這也反映了地震動能量用于地震效應(yīng)安全評價的全面性、合理性。同時，受傳播路徑差異因素影響，爆破振動時程主振、次振周期均為時變的，導致相同距離處各測點各分段震波疊加效果復雜多變。這也解釋了其他條件相同情況下，信號s1、s2、s3、s5各自總能量并不隨間隔時間的增加而呈特定趨勢變化。分別比較信號s1與s7以及s4與s8的總能量，可以看出減震溝具有良好的減震效果，本試驗條件下，減震溝對能量的減震率達到了60%以上。

由圖3及表3可以看出，爆破振動信號的能量在頻域內(nèi)分布較為廣泛，大都集中在0～300 Hz頻帶范圍內(nèi)，在此基礎(chǔ)上，一般又可分為多個子振帶;對比圖2中信號s4、s5、s6可以看出，隨著傳播距離的增加，爆破振動信號高頻能量衰減較快，能量分布最終向低頻帶集中;分別對比表3中信號s1與s7以及s4與s8卓越頻帶，可以得到，減震溝對爆破地震波的高頻部分削弱較為顯著，在減震溝的作用下，爆破振動信號能量卓越頻帶向低頻集中。

相同爆破參數(shù)情況下，信號s4、s5、s6峰值速度與總能量隨傳播距離變化趨勢如圖4所示。

圖4 爆破振動信號能量與峰值速度衰減趨勢Fig.4 The decrease trend about energy and peak velocity of blasting vibration signal

5 結(jié)論

(1)邊坡微差爆破中，微差間隔時間對信號總能量影響較大，適當延長段間延期時間，可減小分段振波疊加、干擾程度，降低前后段波峰疊加的概率。

(2)減震溝對爆破地震波峰值振速與能量的削弱較為顯著，具有良好的減震效果;減震溝對爆破地震波的高頻部分削弱較低頻部分顯著，在減震溝的作用下，爆破振動信號能量卓越頻帶向低頻集中。

(3)一定范圍內(nèi)，爆破振動信號峰值振動速度與總能量隨傳播距離的增加衰減趨勢逐漸變緩;隨著傳播距離的增加，爆破振動信號高頻能量衰減較快，能量分布最終向低頻帶集中。

(4)爆破地震波能量主要集中在0～300 Hz頻帶范圍內(nèi)，在此基礎(chǔ)上，又可分為多個子振頻帶;峰值振動速度大的爆破振動信號攜帶能量不一定大，能量不僅僅是振動速度的表征。

［1］ 林大超，施惠基，白春華.基于小波變換的爆破震動時頻特征分析［J］.巖石力學與工程學報，2004，23(1):101-106.Lin Dachao，Shi Huiji，Bai Chunhua.Time-frequency characteristic analysis of blasting vibration based on wavelet transform［J］.Journal of Rock Mechanics and Engineering，2004，23(1):101-106.

［2］ 張雪亮，黃樹棠.爆破地震效應(yīng)［M］.北京:地震出版社，1981.Zhang Xueliang，Huang Shutang.The Blasting Seismic Effect［M］.Beijing:Seismological Press，1981.

［3］ 徐國元，中國生，熊正明.基于小波變換的爆破地震安全能量分析法的應(yīng)用研究［J］.巖土工程學報，2006，28(1):24-28.Xu Guoyuan，Zhong Guosheng，Xiong Zhengming.Study and application of energy analysis method for blasting seismic safety based on wavelet transform［J］.Journal of Geotechnical Engineering，2006，28(1):24-28.

［4］ 趙明生，梁開水，余德運，等.段數(shù)對爆破振動信號的時頻特性影響分析［J］.煤炭學報，2012，37(1):55-61.Zhao Mingsheng，Liang Kaishui，Yu Deyun，et al.Effect of segments on time frequency characteristics of blasting vibration signals［J］.Journal of China Coal Society，2012，37(1):55-61.

［5］ 凌同華，李夕兵，戴塔根，等.基于小波變換的微差爆破振動信號分離法［J］.地下空間與工程學報，2006，2(3):491-494.Ling Tonghua，Li Xibing，Dai Tagen，et al.The separating means of millisecond blasting vibration signal based on wavelet transform［J］.Journal of Underground Space and Engineering，2006，2(3):491-494.

［6］ 焦永斌.爆破地震安全評定標準初探［J］.爆破，1995(3):45-47.Jiao Yongbin.Study on blasting seismic safety assessment standard［J］.Blasting ，1995(3):45-47.

［7］ 陽生權(quán)，廖先葵，劉寶琛.爆破地震安全判據(jù)的缺陷與改進［J］.爆炸與沖擊，2001，21(3):223-228.Yang Shengquan，Liao Xiankui，Liu Baocheng.Default of the judging standard of blasting vibration safety abstract［J］.Explosion and Shock，2001，21(3):223-228.

［8］ 凌同華，李夕兵.單段爆破振動信號頻帶能量分布特征的小波包分析［J］.振動與沖擊，2007，26(5):41-43.Ling Tonghua，Li Xibing.Features of energy distribution of single deck blast vibration signals with wavelet packet analysis［J］.Journal of Vibration and Shock ，2007，26(5):41-43.

［9］ 劉敦文，粟 闖，龔運高.一種基于爆破振動信號小波分析的爆破危害評判新方法［J］.中南大學學報:自然科學版，2010，41(4):1574-1577.Liu Dunwen，Su Chuang，Gong Yungao.New method for blasting hazards evaluation based on wavelet analysis of blasting vibration signals［J］.Journal of Central South University:Science and Technology，2010，41(4):1574-1577.

［10］ 趙明生，梁開水，李本偉.段藥量對爆破振動信號時頻特性的影響研究［J］.振動與沖擊，2012，31(7):85-88.Zhao Mingsheng，Liang Kaishui，Li Benwei.Influence of deck charge on time-frequency characteristics of a blasting vibration signal［J］.Journal of Vibration and Shock，2012，31(7):85-88.

［11］ 趙明階，葉曉明，吳德倫.工程爆破振動信號分析中的小波方法［J］.重慶交通學院學報，1999，18(3):1-9.Zhao Mingjie，Ye Xiaoming，Wu Delun.The Wavelet-Fourier analysis of blasting seismic signal［J］.Journal of Chongqing Jiaotong U-niversity，1999，18(3):1-9.

［12］ 中國生，徐國元，熊正明.基于小波變換的爆破地震信號能量分析法的應(yīng)用研究［J］.爆炸與沖擊，2006，26(3):222-227.Zhong Guosheng，Xu Guoyuan，Xiong Zhengming.Application research of the energy analysis method for blasting seismic signals based on wavelet transform ［J］.Explosion and shock，2006，26(3):222-227.

［13］ 胡昌華，張軍波，夏 軍，等.基于MATLAB的系統(tǒng)分析與設(shè)計小波分析［M］.西安:西安電子科技大學出版社，2000.Hu Changhua，Zhang Junbo，Xia Jun，et al.Systems Analysis and the Design of Wavelet Analysis Based on MATLAB［M］.Xi'an:Xi'an University of Electronic Science and Technology Press，2000.

［14］ Daubechies I.The wavelet transform time-frequency localization and signal analysis［J］.IEEE Transactions on Information Theory 1990，36(5):961-1005.