曹 琦，傅鶴林，王松周，劉運(yùn)思

(中南大學(xué)土木工程學(xué)院，湖南長沙 410075)

對于爆破振動速度衰減規(guī)律的研究，工程中主要采用蘇聯(lián)學(xué)者薩道夫提出的經(jīng)驗公式V=對實測數(shù)據(jù)進(jìn)行擬合，得到衰減公式。但此經(jīng)驗公式無法準(zhǔn)確地體現(xiàn)場地特征對衰減的影響，也無法分辨場地對不同頻率波的不同衰減作用，因而在應(yīng)用衰減公式預(yù)測非測點爆破振速時，得到的結(jié)果不夠準(zhǔn)確。國內(nèi)不少學(xué)者將小波包變換引入到爆破震動研究中，利用小波變換將爆破震動信號在時域、頻域上展開，根據(jù)不同的頻帶的地震波對結(jié)構(gòu)的危害程度取加權(quán)系數(shù)，將地震波信號在各個頻帶上的主震相的能量按加權(quán)系數(shù)合成，提出了用合成后的能量值作為爆破地震危害程度的判據(jù)。

數(shù)值分析可采用不同地層參數(shù)的模型，有效地體現(xiàn)圍巖特征對振速衰減作用的影響。小波包分析則從頻率這個角度出發(fā)，將復(fù)雜的爆破振動信號分解成不同頻段的振速波，進(jìn)而可以研究場地對不同頻段波的衰減影響。

本文擬運(yùn)用數(shù)值模擬和小波包分析方法，對一組實測振動數(shù)據(jù)分別進(jìn)行模擬和分析，并對非測點數(shù)據(jù)進(jìn)行預(yù)測，力求得到較為合理的預(yù)測值。

1 數(shù)值模擬

1.1 爆破荷載的確定

對于爆破振動的模擬，難點主要在于爆破荷載的計算，目前大部分研究都是針對單孔爆破，實際隧道開挖過程中，往往是多段多孔爆破，本文采用國際上慣用的爆破荷載模式［1－3］:

式中:P'max為多孔爆破峰值荷載;f(t)為時間滯后函數(shù)，取指數(shù)型表達(dá)式:

當(dāng) f(t)=1.0，t=tr時，可得 Ptr計算公式:

式中，tr為爆破沖擊波起始作用時間，計算公式為:

式中:n，m為阻尼系數(shù)，取決于爆破沖擊波起始作用位置和波形;ω為振動頻率，取決于巖體的縱波波速Cp和藥孔直徑Db，按下面公式計算:

由波動理論:

式中:μ為巖石的泊松比;ρr為巖石的密度;E為巖石的彈性模量。

爆破沖擊波在炸藥與巖體界面上給予巖體的最大荷載與巖體特性有關(guān):

式中:Pmax為單孔爆破在巖石表面產(chǎn)生的峰值荷載;V為炸藥的爆速;ρ0為炸藥的密度;P0為最大爆炸荷載，耦合情況下

炮孔內(nèi)的沖擊波作用在巖石圈上的峰值荷載隨著距離而衰減，第n個炮孔的單孔爆破峰值荷載隨著距離變化的衰減關(guān)系如下:

式中:rn為計算點離爆孔的距離;rb為藥孔半徑，是爆破峰值荷載衰減系數(shù)，按公式α=2+μ/(1－μ)計算。

爆轟波在巖石中以波的形式傳播，因此這里按波的疊加規(guī)律，將各個炮孔產(chǎn)生的爆破峰值荷載疊加，求得作用在爆心的爆破峰值荷載，即P'max=。波的疊加取決于波速、頻率、相位及空間傳播方向等，而炸藥類型、藥量、裝藥方式、鉆孔深度、鉆孔角度、鉆孔直徑、鉆孔布置相對位置、起爆時間、場地特征等都會影響波的疊加。本文中近似認(rèn)為同一斷面各個炮孔炸藥類型、藥量、裝藥方式、孔深、孔角、孔徑一致，并假定隧道處于同一種圍巖中。采用微差爆破時，只需選用最大段藥量產(chǎn)生的峰值荷載［4］。

圖1為某隧道III級圍巖上臺階(三臺階)開挖炮眼布置。使用2號巖石乳化炸藥，其爆速為3 000 m/s，炸藥的密度為 1 g/cm3，炮眼直徑為 0.04 m，孔深 2 m，根據(jù)文獻(xiàn)［4］，取阻尼系數(shù) n=0.055，m=0.035。此處圍巖為板巖，其巖石力學(xué)參數(shù)見表1。

圖1 III級圍巖上臺階(三臺階)開挖炮眼布置圖Fig.1 III class rock on the steps(three steps)excavation bore hole layout

表1 巖石物理力學(xué)性質(zhì)Table 1 Physical and mechanical properties of rock

由以上8個參數(shù)，求得:

Pmax=3.52 GPa，α =2.3，rb=0.02 m。再根據(jù)炮眼布置圖求得最大段MS5的峰值荷載:

爆破中心處的爆破荷載隨時間的表達(dá)式為:P(t)=10.95 × (－6.07)× (e－3718t－ e－2367t)GPa爆破荷載隨時間的歷程如圖2所示。

圖2 爆破荷載時間歷程Fig.2 Blast loading time history

1.2 模型建立

為了獲得地層振動特征值，先進(jìn)行特征值分析，在這個過程中使用曲面彈簧來定義彈性邊界。根據(jù)獲得的振動特征值，再重新定義模型的邊界條件，施加爆破荷載，進(jìn)行爆破振動時程分析。從結(jié)果中可以提取測點振速。每一次爆破都獨立進(jìn)行模擬，從而可以分別得到不同次爆破下測點的振速，模型如圖3。由數(shù)值模擬得到的爆破峰值振速模擬值見表2。

圖3 模型圖Fig.3 Model diagram

1.3 數(shù)值模擬結(jié)果

由數(shù)組模擬得到的爆破峰，振速模擬組見表2。

2 小波包分析

2.1 爆破振速實測數(shù)據(jù)

現(xiàn)場實測的爆破振速如表3所示。

2.2 小波包理論

小波包理論是在小波理論的基礎(chǔ)上發(fā)展而來的，能夠?qū)π〔ǚ治鲋械母哳l部分進(jìn)一步分解，并能夠根據(jù)被分析信號的特征，自適應(yīng)地選擇相應(yīng)頻段，使之與信號頻譜相匹配，從而提高時頻分辨率。

表2 爆破峰值振速模擬值Table 2 Blasting simulated peak velocity

表3 各測點隨距離變化的實測峰值振速及頻率Table 3 Measuring points changes with distance and frequency of the measured peak velocity

因而應(yīng)用小波包對地震波信號進(jìn)行頻譜分析，可以得到更為精準(zhǔn)的信息。若存在數(shù)列{hk}，使得函數(shù)族{μk}可以由以下方式構(gòu)成:

則稱函數(shù)族{μk}為一小波包。

特別的，當(dāng)n=0時，有正交小波函數(shù)μ1(t)及對應(yīng)尺度函數(shù)μ0(t)使得:

由此可知，函數(shù)族 {μk(t) }是由尺度函數(shù)μ0(t)確定的小波包。

利用式(11)和式(12)可以得到如下空間分解:

這樣的分解過程，信號的高頻部分也能像低頻部分一樣分解成2部分，從而得到比一般小波分析更精細(xì)的解。

2.3 小波包系數(shù)峰值衰減規(guī)律

爆破地震波在圍巖中傳播時，不同頻率范圍的地震波分量衰減程度是不一樣的，這種衰減反應(yīng)了爆源和場地特征。為了研究爆破地震不同頻率段波在介質(zhì)中的傳播規(guī)律，必須將不同頻率段地震波信號從總信號中提取出來。應(yīng)用小波包技術(shù)提取出不同頻段的地震波分量，然后再逐一分析各分量振動衰減規(guī)律與距離、藥量的直接關(guān)系，在此基礎(chǔ)，對所有頻段峰值速度的衰減規(guī)律進(jìn)行歸一化處理，就可以得到預(yù)測非測點處振動峰值速度的公式［5－7］。圖4為第1次爆破中測點1地震波采用‘db5’小波分解后，頻帶為16～64 Hz的小波包系數(shù)峰值與折算距離的關(guān)系。這表明各小波包信號的振動衰減規(guī)律仍然滿足薩道夫經(jīng)驗公式。

由表3可以看出，第1次爆破中，各測點爆破地震波頻率主要集中在45～50 Hz，70～80 Hz以及100～110 Hz，根據(jù)小波包分析的二分法原則，研究本次信號時，主要采用16 ～64 Hz，64 ～96 Hz，96 ～128 Hz和128～160 Hz 4個頻段。表4給出了以上各頻段振動效應(yīng)與比例距離的衰減系數(shù)以及之間的擬合相關(guān)系數(shù)。

圖4 16～64 Hz小波包系數(shù)峰值與折算距離的關(guān)系Fig.4 16 ～ 64 Hz wavelet packet coefficients of the relationship between the peak and converted from

2.4 建立在小波包系數(shù)峰值基礎(chǔ)上的爆破地震波預(yù)測

為了預(yù)測相同爆破條件下的爆破振動在某測點處的振動峰值波速，必須綜合各頻段的地震波衰減規(guī)律。這里采用如下公式進(jìn)行:

其中:j為頻段序號，這里分為4個頻段;i為原始信號點的序列號;sj(i)為第j頻段第i信號點的小波包系數(shù)。

在進(jìn)行振動峰值波速的預(yù)報時，先選取小波包各頻段的系數(shù)sj(i)，對其中的最高峰值歸一化處理。然后才能采用式(14)對距離為R處峰值波速進(jìn)行預(yù)報。

表4 各頻段振動效應(yīng)衰減參數(shù)Table 4 The effect of attenuation of the vibration frequency parameters

2.5 小波包分析結(jié)果

小波包分析結(jié)果匯總見表5。

3 對比分析

由圖5可知，小波包預(yù)測的測點爆破振速與實測值吻合較好，而數(shù)值模擬振速在離爆源60 m范圍內(nèi)均大于實測值，這是因為計算爆破振動峰值荷載時，是直接將各炮孔爆破產(chǎn)生的振動荷載峰值相加，從而導(dǎo)致計算的峰值荷載較大。并且在模型中加荷載時，是直接加在開挖后的隧道內(nèi)表面，而很多爆孔產(chǎn)生的爆破沖擊波到達(dá)隧道內(nèi)表面時已經(jīng)衰減成應(yīng)力波，應(yīng)力波的值要小于沖擊波，因此這也是導(dǎo)致模擬振速大于實測值的原因。隨著距離的增大，當(dāng)距離超過60 m時，模擬值與實測值趨于一致，這是因為爆破荷載衰減很快，峰值荷載的影響時間很短，所以當(dāng)距離較遠(yuǎn)時，地層對波的衰減作用很明顯，從而弱化了峰值荷載偏大對振速的影響。

表5 基于小波包分析的爆破振速預(yù)測值Table 5 Based on wavelet packet analysis of the predictive value of blasting vibration velocity

圖5 爆破振速對比圖Fig.5 Comparison chart blasting vibration velocity

4 結(jié)論

(1)目前，對于爆破振動荷載的計算大部分是建立在單孔或者單段爆破模型的基礎(chǔ)上，所以用波的疊加理論建立全斷面所有爆破段的爆破模型比較困難，從而導(dǎo)致數(shù)值模擬結(jié)果不夠準(zhǔn)確。

(2)當(dāng)距離超過60 m時，基于小波分析的爆破振速模擬結(jié)果與實測值基本吻合，表明地層對爆破振動的衰減作用明顯。

(3)基于小波包的爆破振速預(yù)測考慮了地層對不同頻段波的影響，從而得到了與實測結(jié)果較為吻合的結(jié)果。

［1］石洪超，丁 寧，張繼春，等.爆破動力作用下小凈距隧道圍巖振動效應(yīng)分析［J］.爆破，2008，5(1):74－78.SHI Hong-chao，DIN Ning，ZHANG Ji-chun，et al.Analysis of vibration effects on surrounding rock for small clear distance tunnel under the dynamic action of blasting［J］.Blasting，2008，5(1):74 －78.

［2］Dowding C H，Aimone C T.Multiple blast－ hole stresses and measured fragmentation ［J］.Rock Mechanic and Rock Engineering，1985(18):17 －36.

［3］ Rathore S S，Bhandari S.Controlled fracture growth by blasting while protecting damages to remaining rock［J］.Rock Mechanics and Rock Engineering，2007，40(3):317－326.

［4］榮 耀.公路隧道爆破荷載的計算分析［J］.公路交通技術(shù)，2005(1):91－94.RONG Yao.Calculation and analysis of explosion loading of highway tunnels［J］.Technology of Highway and Transport，2005(1):91 －94.

［5］徐學(xué)勇.爆破震動信號的小波分析方法研究［D］.武漢:武漢理工大學(xué)，2006.XU Xue-yong，Study on wavelet analysis method for blasting seismic signals［D］.Wuhan:Wuhan University of Technology，2006.

［5］劉 輝，李波，吳從師，等.巖溶隧道掘進(jìn)爆破震動效應(yīng)分析［J］.長安大學(xué)學(xué)報:自然科學(xué)版，2010，30(4):56－59.LIU Hui，LI Bo，WU Cong － shi，et al.Analysis on effects of vibration resulted from driving blasting in Karst tunnel［J］.Journal of Chang’an University:Natural Science Edition，2010，30(4):56 －59.

［6］何本國，朱永全，張志強(qiáng)，等.超小凈距隧道不同爆破方式現(xiàn)場試驗研究［J］.鐵道科學(xué)與工程學(xué)報，2010，7(5):65－ 65.HE Ben-guo，ZHU Yong-quan，ZHANG Zhi-qiang，et al.Study on in－situ test for different blasting methods of the super small－ distance tunnel［J］.Journal of Railway Science and Engineering，2010，7(5):65－ 65.

［7］林從謀，楊林德，崔積弘，等.淺埋隧道掘進(jìn)爆破振動特征研究［J］.地下空間與工程學(xué)報，2006，2(2):276－279.LIN Cong-mou，YANG Lin-de，CUI Ji-hong，et al.Study on the vibration characteristic of the shallow－tunneling blasting［J］.Chinese Journal of Underground Space and Engineering，2006，2(2):276 －279.

［8］周志強(qiáng)，易建政，王 波，等.控制爆破技術(shù)研究現(xiàn)狀及發(fā)展建議［J］.礦業(yè)研究與開發(fā)，2010，30(3):103－108.ZHOU Zhi-qiang，YI Jian-zheng，WANG Bo，et al.Present situation and suggestions of research and development of controlled blasting technology［J］.Mining Research and Development，2010，30(3):103 －108.