李興華，龍 源，紀(jì) 沖，鐘明壽，路 亮

(解放軍理工大學(xué) 野戰(zhàn)工程學(xué)院，南京 210007)

爆破方法是核電站建設(shè)、水利樞紐工程基礎(chǔ)、山區(qū)高速公路路塹開挖以及露天礦山開采等工程項目的重要施工手段。在爆破施工現(xiàn)場，地形一般起伏較大，爆破地點(diǎn)與周圍的建(構(gòu))筑物或者需要保護(hù)的目標(biāo)之間存在一定的高程差。研究表明在地形起伏較大的情況下，地形地貌對爆破地震波傳播具有較大影響［1］。

沈蔚等［2］采用灰關(guān)聯(lián)分析法對爆破震動參數(shù)的主要影響因素進(jìn)行分析，分析表明高程差為影響爆破震動參數(shù)的準(zhǔn)優(yōu)因素。目前對高程差對爆破地震波影響研究主要集中在高程差對爆破地震波峰值影響研究［3－6］，而對爆破震動信號能量分布影響研究較少。但是爆破震動對建(構(gòu))筑物的破壞是一個動態(tài)響應(yīng)的問題［7］。在動態(tài)響應(yīng)分析過程中不僅要考慮爆破地震波峰值，還應(yīng)該考慮爆破震動主頻、振動持續(xù)時間以及信號中的能量分布對建(構(gòu))筑物的破壞作用，因此有必要對高程差影響下的爆破震動信號能量分布做進(jìn)一步的研究。

本文依托某核電站核島基坑負(fù)挖工程采集了正高程差(測點(diǎn)位置高于爆源位置)影響下的爆破震動信號，采用小波包分析方法對采集的信號進(jìn)行了深入的分析。旨在深入研究高程差對爆破地震波的影響，為優(yōu)化爆破設(shè)計，減少爆破震動危害提供依據(jù)。

1 負(fù)挖爆破產(chǎn)生的震動信號采集

1.1 試驗背景

信號采集試驗結(jié)合核電站核島基坑負(fù)挖爆破項目進(jìn)行。核島基坑附近松散層主要由第四紀(jì)的沖積層、海積層、坡積層和殘積層組成，基巖主要是燕山期的侵入巖。巖石為微風(fēng)化斑狀、中細(xì)?；◢弾r，巖石硬度系數(shù)8～10。在基坑負(fù)挖爆破區(qū)域西側(cè)大約80 m處有一坡度約為70°的邊坡，邊坡坡腳標(biāo)高為+8 m，共有5個臺階。

圖1 爆破振動監(jiān)測測點(diǎn)布置Fig.1 Arrangement of blasting vibration monitoring points

爆破負(fù)挖工程分兩層進(jìn)行，第一層負(fù)挖深度為5.5 m，第二層負(fù)挖深度為4.5 m。爆破開挖方式為中間為掏槽爆破，四周為臺階爆破，基坑邊緣為預(yù)裂爆破。其中臺階爆破的孔徑為89 mm，單耗為0.52 kg/m3，最小抵抗線為1.2 m。

1.2 測點(diǎn)布置

在爆破負(fù)挖期間分別在邊坡32 m臺階和56 m臺階上設(shè)置了一個測點(diǎn)1和測點(diǎn)2，測點(diǎn)位于臺階的中央，兩測點(diǎn)上下在同一垂直平面內(nèi)，測點(diǎn)與爆區(qū)中心的連線為測量時的東西方向。試驗場地情況及傳感器布置示意圖如圖1所示。

1.3 爆破震動測試系統(tǒng)

本次爆破震動監(jiān)測所用儀器為美國WHITE公司生產(chǎn)Mini－SeisⅡ微型振動速度儀，該儀器的最大采樣頻率為2 048 Hz，最低響應(yīng)頻率可達(dá)1 Hz。Mini－SeisⅡ型爆破地震儀系統(tǒng)原理圖如圖2所示。

圖2 Mini-SeisⅡ型爆破地震儀系統(tǒng)原理圖Fig.2 The principle of Mini－SeisⅡtesting apparatus

表1 爆破地震波速度峰值回歸分析參數(shù)Tab.1 The regression analysis parameters of velocity peak of blasting seismic wave

圖3 第二層負(fù)挖爆破測點(diǎn)速度峰值對比圖Fig.3 Comparison of velocity peak between different measuring points at the second negative excavation blasting

2 高程差對爆破震動信號峰值的放大效應(yīng)

將表1中的K，α值代入經(jīng)驗公式，得到爆破震動速度峰值的計算公式。利用經(jīng)驗公式做出第二次負(fù)挖三個方向速度峰值曲線的比較圖。觀察圖3，發(fā)現(xiàn)第二層負(fù)挖時水平徑向速度峰值和垂直方向峰值速度，在一定的比例距離范圍內(nèi)，在相同的比例距離下，測點(diǎn)2的速度峰值大于測點(diǎn)1的速度峰值。這說明當(dāng)爆破地震波沿正的高程差地形傳播時，高程對爆破震動峰值具有放大效應(yīng)。因此在對具有高程差地形爆破震動進(jìn)行預(yù)報時，高度差是一個不可忽略的因素。

再觀察圖3水平切向方向速度峰值曲線，發(fā)現(xiàn)在相同的比例距離下，測點(diǎn)1的速度峰值大于測點(diǎn)2的速度峰值，即在高程差影響下震動峰值并沒有表現(xiàn)出放大現(xiàn)象，說明高程對爆破震動峰值的放大效應(yīng)具有方向性。

3 高程差對爆破震動信號能量分布的影響

3.1 爆破震動信號的小波包分析

小波包分析為信號提供一種更精細(xì)的分析方法，將頻帶進(jìn)行多層次劃分，對小波分析沒有細(xì)分的高頻部分進(jìn)一步分解，實(shí)現(xiàn)了頻帶的均勻劃分，把特征信號的能量集中到更小更均勻的頻帶中，時頻特性比小波變換更精確，使之在處理爆破地震波這種非平穩(wěn)信號方面具有極強(qiáng)的優(yōu)越性［8］。

一個給定信號x(t)若進(jìn)行i層小波包分解在該次分解中可以得到2i個子頻帶。若原始信號最低頻率成分為 0，最高頻率成分為ωm，每個子頻帶寬度為ωm/2i［9］。爆破震動信號的頻率一般在200 Hz以下，本次試驗的爆破震動記錄儀器的信號采樣頻率1 024 Hz，根據(jù)采樣定理，則其奈奎斯特(Nyquist)頻率為512 Hz。因此，對爆破震動信號進(jìn)行7層分解，對應(yīng)的最低頻帶為0～4 Hz。

3.2 基于小波包變換爆破震動信號能量表示方法

對一爆破震動信號x(t)進(jìn)行i層小波包分解，進(jìn)行小波包分解系數(shù)重構(gòu)，可以提取各頻帶范圍內(nèi)的信號，且總信號可以表示為:

式中:(i，k)表示第i層分解節(jié)點(diǎn)(i，k)上的重構(gòu)信號，其中，j=2i，k=0，1，2，…，j－1。

如果采用二次能量型時頻表示對應(yīng)于每個頻帶上的重構(gòu)信號，可以定義時頻譜［10］:

式中:fk為第k個頻帶的中心頻率，第k個頻帶信號的總能量為:

式中:vk，l為重構(gòu)信號 xi，k(t)離散點(diǎn)的幅值(l=1，2，3，…，n為爆破震動信號離散點(diǎn)個數(shù))。

第k個頻帶的能量占爆破震動信號總能量的比例為:

根據(jù)式(2)、式(3)可以近似得到在整個時域和頻域上爆破震動信號x(t)的連續(xù)分布時頻譜。由式(4)可以得到信號經(jīng)小波包分解后不同頻帶的能量百分比，從而可以找出爆破震動信號在傳播過程中能量的分布規(guī)律。

3.3 高程差對信號不同頻段能量分布影響的分析

(1)爆破震動信號的速度時程曲線

從第一層負(fù)挖爆破中選取用于小波包分析的爆破震動信號，相關(guān)的爆破參數(shù)如表2所示，相應(yīng)的速度時程曲線如圖4和圖5。

(2)爆破震動信號的能量分布

對圖3和圖4所示的爆破震動信號采用db8［11－12］小波基函數(shù)進(jìn)行深度為7層的小波包分解。根據(jù)式(1)～(4)編制相應(yīng)的MATLAB計算程序，從而獲得上述爆破地震信號的各頻帶能量分布情況圖，如圖6所示。

對比表2中爆破震動信號的主頻發(fā)現(xiàn)，測點(diǎn)2的主頻較測點(diǎn)1偏低。利用各頻段能量分布圖可以得到各主頻所在的頻段的能量在信號總能量中所占的比例較低，最高的為7.8%，最低的為0.7%。

由上述負(fù)挖爆破震動信號的頻帶能量分布圖可以看出爆破地震波的頻率分布較廣，但是能量絕大部分集中在20～100 Hz，占到總能量的80%以上。對比分析位于不同高程兩個測點(diǎn)的爆破震動信號不同頻帶能量分布的特征，可以得到在高程效應(yīng)的影響下，隨著高度的增大，爆破震動信號能量的分布越來越傾向低頻帶，主振頻帶有向低頻發(fā)展的趨勢。爆破震動信號中低頻帶的能量占整個信號能量的比例越來越大，而高頻帶的能量所占比例越來越小。說明高程差的存在有利于高頻成分的衰減。由于建(構(gòu))筑物的固有頻率一般都比較低，在工程爆破中，這顯然不利于周圍建(構(gòu))筑物的安全。

表2 爆破地震波參數(shù)Tab.2 Parameters of blasting seismic wave

圖4 測點(diǎn)1(高程32 m)爆破震動速度時程曲線Fig.4 Velocity curve of blasting vibration at measuring points 1

圖5 測點(diǎn)2(高程56 m)爆破震動速度時程曲線Fig.5 Velocity curve of blasting vibration at measuring points 2

圖6 爆破震動信號的頻帶能量分布對比Fig.6 Comparison of the energy distribution in different frequency bands of blasting vibration signal

圖7 第一層負(fù)挖爆破T方向0～100 Hz時頻分布Fig.7 Time-frequency distribution of 0 ～100 Hz of transverse velocity at the first negative excavation blasting

(3)爆破震動信號的三維時頻分布

根據(jù)式(2)、式(3)利用MATLAB可以得到爆破震動信號的三維時頻譜圖。根據(jù)爆破震動信號能量的分布圖顯示第一層負(fù)挖爆破炮次1的T方向的能量分布向低頻發(fā)展的趨勢較明顯。為了能比較清楚地顯示結(jié)果，第一層負(fù)挖爆破炮次1的T方向震動信號0～100 Hz頻帶范圍內(nèi)的時頻譜圖如圖7所示。爆破震動信號的大部分能量集中0～100 Hz，因此重點(diǎn)分析這一頻段的能量分布。

基于小波包分析的時頻圖上可以清楚的看出爆破震動信號的能量隨時間和頻率的分布。對比第一層爆破負(fù)挖炮次1兩個測點(diǎn)T方向0～100 Hz頻段范圍內(nèi)的能量分布，測點(diǎn)二整體能量與測點(diǎn)一相比較低，但是低頻成分中的子頻帶增多，相對來說低頻中的能量是增大了。測點(diǎn)二的能量分布與測點(diǎn)一相比有向低頻發(fā)展的趨勢。

4 結(jié)論

(1)依托工程實(shí)踐采集了高程差影響下的爆破震動數(shù)據(jù)，采用最小二乘法回歸分析得到不同高程下的經(jīng)驗公式。此包含高程因素的經(jīng)驗公式能夠為負(fù)挖爆破地點(diǎn)周圍核電設(shè)施的震動監(jiān)測提供理論依據(jù);

(2)對比分析結(jié)果表明，正的高程差對爆破震動峰值有放大效應(yīng)，并且這種高程的放大效應(yīng)具有一定的方向性和一定的范圍性;

(3)基于小波包分析得到了高程效應(yīng)下爆破震動信號的主頻和能量分布的規(guī)律，隨著高程的增大，爆破震動信號的能量有向低頻發(fā)展的趨勢，這對于爆點(diǎn)周圍的建(構(gòu))筑物的防護(hù)是不利的。

［1］唐 海，李俊如.凸形地貌對爆破震動波傳播影響的數(shù)值模擬［J］巖土力學(xué)，2010，31(4):1289 －1294.

TANG Hai，LI Jun-ru.Numerical simulation influence of protruding topography on blasting vibration wave propagation［J］.Rock and Soil Mechanics，2010，31(4):1289 －1294.

［2］沈 蔚，徐全軍，黃文華.用灰關(guān)聯(lián)分析法確定爆破震動參數(shù)的主要影響因素［J］.工程爆破，2000，6(4):6－8.

SHEN Wei，XU Quan-jun，HUANG Wen-hua.Determination of main influencing factors on blasting vibration parameters by grey correlation analysis［J］.Engineering blasting，2000，6(4):6－8.

［3］周同齡，楊秀甫，翁家杰.爆破地震高程效應(yīng)的實(shí)驗研究［J］.建井技術(shù)，1997，18(增刊):31 －35.

ZHOU Tong-lin， YANG Xiu-pu， WENG Jia-jie.Experimental research on height effect of blasting vibration［J］.Mine Construction Technology，1997，18:31－35.

［4］周同齡，李玉壽.反映高程的爆破震動公式及其應(yīng)用［J］.江蘇煤炭，1997，4:21 －22.

ZHOU Tong-ling， LIYu-shou.Application ofblasting vibration formula related to highness［J］.Jiangsu Coal，1997，(4):21－22.

［5］張 濤，郭學(xué)彬，蒲傳金.邊坡爆破振動高程效應(yīng)的實(shí)驗研究與分析［J］.江西有色金屬，2006，20(4):10－13.

ZHANG Tao，GUO Xue-bin，PU Chuan-jin.Experimental analysis and studies of elevate effect of slope blasting vibration［J］.Jiangxi Nonferrous Metals，2006，20(4):10 －13.

［6］唐 海，李海波，蔣鵬燦.地形地貌對爆破振動波傳播的影響實(shí)驗研究［J］.巖土力學(xué)與工程學(xué)報，2007，26(9):1817－1823.

TANG Hai，LI Hai-bo，JIANG Peng-can.Experimental study on effect of topography on propagation of blasting waves［J］.Rock Mechanics and Engineering，2007，26(9):1817－1823.

［7］婁建武，龍 源，方 向.基于反應(yīng)譜值分析的爆破震動破壞評估研究［J］.爆炸與沖擊，2003，23(1):41－46.

LOU Jian-wu，LONG Yuan，F(xiàn)ANG Xiang.Study on blasting vibration dam age based on response spectrum［J］.Explosion and shock waves，2003，23(1):41 －46.

［8］婁建武，龍 源，徐全軍.基于小波包技術(shù)的爆破地震波特征提取及預(yù)報［J］.爆炸與沖擊，2004，24(3):261－267.

LOU Jian-wu，LONG Yuan，XU Quan-jun.A study on the extraction and prediction of blasting seismic wave characteristics by wavelet packets technique［J］.Explosion and shock waves，2004，24(3):261 －267.

［9］中國生，傲麗萍，趙 奎.基于小波包能量譜爆炸參量對爆破振動信號能量分布的影響［J］.爆炸與沖擊，2009，29(3):300－305.

ZHONG Guo-sheng，AO Li-ping，ZHAO Kui.Influence of explosion parameters on energy distribution ofblasting vibration signal based on wavelet packet energy spectrum［J］.Explosion and shock waves，2009，29(3):300 －305.

［10］林大超，施惠基，白春華.爆破震動時頻分布的小波包分析［J］.工程爆破，2006，8(2):1－5.

LIN Da-chao， SHI Hui-jin， BAI Chun-hua. Ontimefrequency distribution of blasting vibration through wavelet packet transform［J］.Engineering blasting，2006，8(2):1－5.

［11］凌同華，李夕兵.單段爆破振動信號頻帶能量分布特征的小波包分析［J］.振動與沖擊，2007，26(5):41－43.

LING Tong-hua，LI Xin-bin.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.

［12］李興華.高程差對爆破震動信號特性影響分析［D］.南京:解放軍理工大學(xué)，2010.