費鴻祿，曾翔宇，楊智廣

（遼寧工程技術(shù)大學爆破技術(shù)研究院，遼寧阜新123000）

隧道掘進爆破振動對地表影響的小波包分析＊

費鴻祿，曾翔宇，楊智廣

（遼寧工程技術(shù)大學爆破技術(shù)研究院，遼寧阜新123000）

基于現(xiàn)場實測隧道掘進爆破振動信號數(shù)據(jù)，采用小波包分析技術(shù)對振動信號進行小波包能量譜分析，得到爆破振動信號能量在各頻帶上的分布。通過比較各點分析結(jié)果可以發(fā)現(xiàn)隧道掘進爆破地震波能量主要集中在切向與徑向。在沿隧道掘進方向隨著與掌子面之間距離的增大地震波的主頻下降，能量主要集中頻帶也越來越窄，并且向低頻帶發(fā)展；而在垂向上地震波主頻隨著埋深的減小而減小，同時地震波的能量主要集中頻帶也越來越窄并且向低頻帶發(fā)展。最后從能量的角度探討隧道掘進爆破地震波沿各個方向的衰減規(guī)律。

爆炸力學；小波包分析；能量分布；地震波；爆破振動

近年來，小波變換以其良好的高分辨率和自適應性以及在時域和頻域都具有表征信號局部特性的能力，被廣泛應用于對爆破振動信號的分析和處理中［4］，但是目前利用小波變換處理爆破振動信號的技術(shù)還處于起步階段［5－7］。凌同華等［8－9］利用小波理論確定了微差爆破的實際延遲時間并研究不同爆破參數(shù)對爆破主頻帶和能量分布的影響；婁建武等［10］研究了小波分析在結(jié)構(gòu)爆破振動響應能量分析法中的應用；中國生等［11］基于小波變換爆破地震信號能量分析法的應用進行了研究。但由于小波變換只能對信號的頻段進行指數(shù)等間距劃分，所以在高頻段其頻率分辨率較差［12］。小波包分析法對沒有細分的高頻部分進行進一步分解，能夠根據(jù)被分析信號的特征自適應地選擇相應頻段，使之與信號頻譜相匹配，從而提高了時頻分辨率，具有更廣泛的應用價值［13－14］。

本文中基于現(xiàn)場實測隧道掘進爆破振動數(shù)據(jù)，采用小波包分解和重構(gòu)的分析方法，從能量的角度探討沿隧道掘進方向及波的垂向（x）、切向（y）、徑向（z）地震波的傳播、衰減規(guī)律以及能量在不同頻帶下的分布情況。

1 小波包能量譜分析

1.1 小波包分析法

小波分析原理是把信號分解成低頻和高頻2個部分，并且在下一輪分解中，僅對分解出的低頻部分進行再次分解，以此類推至分解完成。從小波分析法中可以看出，在小波分解中低頻波的分辨率要高于高頻波，導致在波的高頻部分分辨率較低。然而在小波包分解過程中，每一輪分解不僅依次對波的低頻部分進行分解，對高頻部分同樣進行分解。由此可見小波包分析法提高了波在高頻段的分辨率，比小波分析更加精細［15－17］。

爆破振動數(shù)據(jù)采集采用TC－4850爆破測振儀，信號的記錄時長為2s，采集頻率為2kHz，則其奈奎斯特頻率為1kHz。采用小波包分解法可將數(shù)據(jù)分解8層，經(jīng)第1層分解后被分成低頻SA0（0，500Hz）和高頻SA1（500Hz，1 000Hz）等2個部分；第2層小波分解則是對第1層分解后的SA0段和SA1段進行進一步分解，對SA0分解得到低頻段SB0（0，250Hz）和高頻段SB1（250Hz，500Hz），對SA1分解得到低頻段SB2（500Hz，750Hz）和高頻段SB3（750Hz，1 000Hz）；以此類推，當分解到第8層即可得到SH0（0，3.906Hz）、SH1（3.906Hz，7.813Hz）、SH2（7.813Hz，11.718Hz）直至SH255（999.094Hz，1 000Hz）。

1.2 小波包能量譜分析原理

對爆破振動信號S（t）進行小波包分解后，第i層可以得到2i個子頻帶，則S（t）可以表示為［18－19］：

ST-SNE同樣使用Kullback-Leibler距離來最小化概率矩陣P與Q之間的差異，目標函數(shù)為：

式中fi，j（tj）為爆破振動信號運用小波包分解到第i層節(jié)點（i，j）上的重構(gòu)信號。若信號S（t）的頻率寬度為ω，則在第i層每個子頻帶的頻率寬度為ω／2i。

若Ei，j（tj）為fi，j（tj）所對應的能量，則有：

式中：xj，k（j＝0，1，2，…，2i－1；k＝1，2，…，m）為fi，j（tj）離散點的幅值，m為爆破振動信號采集點數(shù)。

由式（2）可得爆破振動信號S（t）的總能量E為：

爆破振動信號S（t）小波包分解到第i層時，各頻帶能量占信號總能量的比例為：

2 地震波在地表傳播衰減性質(zhì)

2.1 地震波信號現(xiàn)場監(jiān)測

鶴大高速小溝嶺隧道位于吉林省敦化市，該段圍巖主要為碎石，強風化杏仁狀玄武巖，中風化杏仁狀玄武巖，巖體飽和抗壓強度Rc＝6.2MPa，巖石完整系數(shù)Kv＝0.11，圍巖級別為Ⅴ級。該段巖性完全相同，未發(fā)現(xiàn)有斷層、溝坎的現(xiàn)象，有利于排除因地質(zhì)因素而對實驗結(jié)果產(chǎn)生的差異。監(jiān)測共進行5組，每組取點情況如圖1（a），每組監(jiān)測現(xiàn)場共設置5個點（A～E），監(jiān)測點的布置如圖1（b），其中B點位于掌子面正上方，相鄰兩點水平間隔為10m。現(xiàn)場共進行5次監(jiān)測，各次監(jiān)測的爆破參量如表1所示。表1中h為掌子面處地面高程，m；H為掌子面處設計高程，m；Q為總炸藥量，kg；q為單段最大藥量，kg；n表示雷管段數(shù)。

圖1 監(jiān)測點的布置Fig.1 Layout of survey points

由于第2組監(jiān)測地勢較為平緩（此時掌子面正上方坡率約為0.1），可以排除由于高差而產(chǎn)生的影響，所以抽取第2組監(jiān)測所得數(shù)據(jù)進行小波包分析，通過對比探討地震波沿隧道縱向傳播衰減規(guī)律。此外提取第4組中B點數(shù)據(jù)做同樣分析，通過與第2組中B點數(shù)據(jù)進行比較，探討地震波沿隧道垂向傳播衰減規(guī)律，為與第2組中的B點區(qū)分開，第4組中的B點在下文中稱為F點。

表1 現(xiàn)場測試相關(guān)參數(shù)Table 1 Site test related parameters

2.2 爆破振動信號的小波包能量譜分析

在利用小波包分析法對爆破振動進行分析的工程中，選擇恰當?shù)男〔ɑ陵P(guān)重要。目前db8以其良好的緊支撐性、光滑性以及近似對稱性，已成功地應用于分析爆破振動信號中［20］。本文中采用db8對爆破振動信號進行8層小波包分解。表2所示為各監(jiān)測點爆破振動信號通過小波包能量譜分析所得數(shù)據(jù)，圖2為經(jīng)過小波包分析法處理后所得的爆破振動信號能量譜分布圖。

表2 爆破振動信號小波包分析結(jié)果Table 2 Wavelet packet analysis results of blasting vibration signals

2.3 地震波在地表傳播衰減性質(zhì)

從表2中可以看出，15個信號的主頻帶都集中于0～200Hz，主頻帶能量在總能量中的比例普遍不低于60%，并且從圖4可知雖然爆破振動信號的能量在頻域上的分布較廣，但絕大多數(shù)能量主要集中在0～200Hz之間。結(jié)合表2及圖3中掌子面同側(cè)監(jiān)測點B～E各數(shù)據(jù)可以看出，隨著監(jiān)測點與掌子面之間距離的增大，爆破振動信號的主頻會逐漸降低，并且其衰減速度隨著距離的增大而逐漸減緩。

第2組的5個監(jiān)測點中，B點位于掌子面正上方，設其水平坐標為0，以隧道掘進方向為正向，5個點的水平坐標依次為10、0、－10、－20、－30。各點處總能量和主頻隨距離的變化關(guān)系如圖3～4所示。

圖2 各測點爆破振動信號能量譜分布Fig.2 Energy distribution of blasting vibration signals at measuring points

圖3 各監(jiān)測點總能量與距離的關(guān)系Fig.3 Relationship between energy and distance at measuring points

圖4 各監(jiān)測點主頻與距離關(guān)系Fig.4 Relationship between main frequency and distance at measuring points

比較圖4中近掌子面點B以及距離掌子面最遠點E的圖像可以看出隨著距離的增大，地震波的能量逐漸集中于低頻帶，由于建（構(gòu)）筑物的固有頻率普遍較低，所以當?shù)卣鸩ǖ念l率接近建（構(gòu)）筑物的固有頻率便產(chǎn)生危害。

比較表2、圖2及圖3中掌子面正上方的B點及掌子面前后對應點A與C振動波的各參數(shù)可以看出，雖然在掌子面前后對應位置爆破振動信號的主頻帶并沒有太大差別，但是在隧道未開挖一側(cè)信號的3個分量主頻均小于已開挖一側(cè)，而總能量均大于已開挖一側(cè)。由此可見在爆破近區(qū)隧道未開挖一側(cè)地表建筑物安全性比已開挖一側(cè)差。

通過表2、圖2及圖3中5個監(jiān)測點的數(shù)據(jù)不難看出，在同一個監(jiān)測點波的3個分量中，垂向（x方向）的主頻均大于切向（y方向）與徑向（z方向）的主頻，而波在垂向中所含的能量均小于切向與徑向波所含的能量。由此可見在爆破振動波的傳播過程中波在切向與徑向的頻率低、衰減慢、攜帶較多能量，是造成振動破壞的主要原因。

B點到F點，埋深從51.877m下降到15.074m，在2次炸藥用量相差不大的前提下，爆破地震波的主頻從50～90Hz下降到40～55Hz，并且隨著埋深的下降，爆破地震波的主振頻帶逐漸變窄，且能量有向低頻集中的趨勢。

3 結(jié) 論

基于小波包能量譜分析技術(shù)對爆破地震波進行技術(shù)分析，得到的地震波能量在傳播過程中一般衰減規(guī)律結(jié)論如下：

（1）爆破振動信號傳播到地表的能量在頻譜上分布較窄，主要分布在0～200Hz之間，其中切向與徑向相對于垂向主頻較低，能量較高，對地表建（構(gòu)）筑物影響相對較大。

（2）地震波總能量在爆破近區(qū)衰減較快、遠區(qū)衰減較慢；且隨著距離的增加高頻波被過濾，低頻波所含能量比例明顯上升，這對于處在爆破遠區(qū)、固有頻度較低的建筑及設施來說是不利的。

（3）地震波在隧道未開挖一側(cè)主頻小于已開挖一側(cè)，地震波所含能量大于已開挖一側(cè)，說明在掌子面對應位置上，隧道未開挖一側(cè)的建筑、設施相對于已開挖一側(cè)受破壞可能性更大。

（4）當隧道埋深較淺時，地震波的主頻也會相對較小，而且振動信號的主振頻帶越來越窄，能量主要集中在低頻帶。由此可見地震波在垂向傳播、衰減規(guī)律與沿隧道掘進方向有所不同。

［1］陳士海，魏海霞，張安康，等.基于小波包技術(shù)的爆破地震效應計算模型及安全判據(jù)研究［J］.爆炸與沖擊，2010，30（4）：377－382.Chen Shihai，Wei Haixia，Zhang Ankang，et al.Computational model and safety criterion of blasting vibration effect based on wavelet packet techniques［J］.Explosion and shock waves，2010，30（4）：377－382.

［2］管曉明，傅洪賢，王夢恕.隧道近距下穿山坡樓房爆破振動測試研究［J］.巖石力學，2014，35（7）：1995－2003.Guan Xiaoming，F(xiàn)u Hongxian，Wang Mengshu.Blasting vibration characteristics monitoring of tunnel under－passing hillside buildings in short－distance［J］.Rock and Soil Mechanics，2014，35（7）：1995－2003.

［3］閆鴻浩，劉智遠，李曉杰.煙囪爆破拆除及其振動測試分析［J］.爆破器材，2013，42（3）：21－24.Yan Honghao，Liu Zhiyuan，Li Xiaojie.Chimney blasting demolition and analysis of vibration measurement［J］.Explosive Materials，2013，42（3）：21－24.

［4］路亮，龍源，謝全民，等.爆破振動信號的提升小波包分解及能量分布特征［J］.爆炸與沖擊，2013，33（2）：140－147.Lu Liang，Long Yuan，Xie Quanmin，et al.Decomposition and energy distribution of blasting vibration signal based on second generation wavelet packet［J］.Explosion and shock waves，2013，33（2）：140－147.

［5］何軍，于亞倫，梁文基.爆破振動信號的小波分析［J］.巖土工程學報，1998，20（1）：47－50.He Jun，Yu Yalun，Liang Wenji.Wavelet analysis for blasting seismic signals［J］.Chinese Journal of Geotechnical Engineering，1998，20（1）：47－50.

［6］黃文華，徐全軍，沈蔚.小波變換在判斷爆破地震危害中的應用［J］.工程爆破，2001，7（1）；24－27.Huang Wenhua，Xu Quanjun，Shen Wei，et al.Application of wavelet transform in evaluation of blasting vibration damage［J］.Engineering Blasting，2001，7（1）：24－27.

［7］嚴鵬，盧文波，羅憶，等.基于小波變換時－能密度分析的爆破開挖過程中地應力動態(tài)卸載振動到達時刻識別［J］.巖石力學與工程學報，2009，28（S1）：2836－2845.Yan Peng，Lu Wenbo，Luo Yi，et al.Identification of arriving time of vibration induced by geostress dynamic unloading during blasting－excavation employing method of time－energy analysis based on wavelet transform［J］.Chinese Journal of Rock Mechanics and Engineering，2009，28（suppl 1）：2836－22844.

［8］凌同華，李夕兵.基于小波變換的時－能分布確定微差爆破的實際延遲時間［J］.巖石力學與工程學報，2004，23（13）：2266－2270.Ling Tonghua，Li Xibing.Time－energy analysis based on wavelet transform for identifying real delay time in millisecond blasting［J］.Chinese Journal of Rock Mechanics and Engineering，2004，23（13）：2260－2270.

［9］凌同華，李夕兵.地下工程爆破振動信號能量分布特征的小波包分析［J］.爆炸與沖擊，2004，24（1）：63－68.Ling Tonghua，Li Xibing.The features of energy distribution for blast vibration signals in underground engineering by wavelet packet analysis［J］.Explosion and Shock Waves，2004，24（1）：63－68.

［10］婁建武，龍源，徐全軍.小波分析在結(jié)構(gòu)爆破振動響應能量分析法中的應用［J］.世界地震工程，2001，17（1）：64－68.Lou Jianwu，Long Yuan，Xu Quanjun.Study on the wavelet analysis applied in structure response to blasting vibration［J］.World Information on Earthquake Engineering，2001，17（1）：64－68.

［11］中國生，徐國元，熊正明.基于小波變換的爆破地震信號能量分析法的應用研究［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 Waves，2006，26（3）：222－227.

［12］付曉強，張世平，張昌鎖.露天開采爆破振動信號小波包分析［J］.工程爆破，2013，19（1／2）：24－27.Fu Xiaoqiang，Zhang Shiping，Zhang Changsuo.Analysis of open－pit mining blasting vibration signals using wavelet packet method［J］.Engineering Blasting，2013，19（1／2）：24－27.

［13］Zhao M S，Zhang J H，Yi C P.Time－frequency characteristics of blasting vibration signals measured in milliseconds［J］.Mining Science and technology，2011，21（3）：349－352.

［14］Johansson D，Ouchterlony F.Shock wave interactions in rock blasting：The use of short delays to improve fragmentation in model－scale［J］.Rock mechanics and Rock Engineering，2013：46（1）：1－18.

［15］凌同華，李夕兵.多段微差爆破振動信號頻帶能量分布特征的小波包分析［J］.巖石力學與工程學報，2005，24（7）：1117－1122.Ling Tonghua，Li Xibing.Analysis of energy distributions of millisecond blast vibration signals using the wavelet packet method［J］.Chinese Journal of Rock Mechanics and Engineering，2005，24（7）：1117－1122.

［16］王利，王海亮.隧道爆破不同類型最大振速對建筑物的影響研究［J］.爆破，2012，29（4）：6－9.Wang Li，Wang Hailiang.Influence of different direction vibration peak velocity of tunnel blasting on building［J］.Blasting，2012，29（4）：6－9.

［17］謝全民，龍源，鐘明壽，等.小波與分形組合分析技術(shù)在爆破振動信號分析中的應用［J］.振動與沖擊，2011，30（12）：120－124.Xie Quanmin，Long Yuan，Zhong Mingshou，et al.Blasting vibration signal analysis with wavelet and fractal portfolio analysis technique［J］.Journal of Vibration and Shock，2011，30（12）：120－124.

［18］張德豐.MATLAB小波分析［M］.北京：機械工業(yè)出版社，2009.

［19］中國生，敖麗萍，趙奎.基于小波包能量譜爆炸參量對爆破振動信號能量分布的影響［J］.爆炸與沖擊，2009，29（3）：300－305.Zhong Guosheng，Ao Liping，Zhao Kui.Influence of explosion parameters on energy distribution of blasting vibration signal based on wavelet packet energy spectrum［J］.Explosion and Shock Waves，2009，29（3）：300－305.

［20］閆鴻浩，李曉杰，曲艷東，等.爆破振動速度測試精細分析［J］.巖土力學，2007，28（10）：2091－2094.Yan Honghao，Li Xiaojie，Qu Yandong，et al.Fine analysis of blasting vibration velocity testing［J］.Rock and Soil Mechanics，2007，28（10）：2091－2094.

Influence of tunnel excavation blasting vibration on earth’s surface based on wavelet packet analysis

Fei Honglu，Zeng Xiangyu，Yang zhiguang
（Institute of Blasting Technique，Liaoning Technical University，F(xiàn)uxin 123000，Liaoning，China）

In this work，based on the blasting vibration live data collected，the energy distributions at different frequency bands for the blasting vibration signals were obtained by the wavelet packet analysis technology.By comparing the results of each point analysis it was found that the tunnel excavation blasting seismic wave energy mainly concentrated in the tangential and radial directions；along the direction of the tunnelling，with the increase of the distance between the tunnel faces，the main frequency of the seismic wave decreased and the energy concentration frequency band became narrower and developed to the low frequency band；however，along the vertical direction of the tunnel，the main frequency decreases with the decrease of the depth and，at the same time，the wave energy concentration frequency band became narrower and developed to the low frequency band.The results show that wavelet packet analysis is an effective approach to estimating the structure safety under blasting vibration conditions.

mechanics of explosion；wavelet packet analysis；energy distribution；blasting seismic；blasting vibration

O383國標學科代碼：13035

A

10.11883／1001－1455（2017）01－0077－07

（責任編輯 王易難）

2015－04－30；

2015－08－25

費鴻祿（1963— ），男，博士，教授，feihonglu＠163.com。