曾一鑫，白春華，王仲琦

(北京理工大學(xué) 爆炸科學(xué)與技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京 100081)

現(xiàn)階段國(guó)內(nèi)大型油氣田開(kāi)發(fā)多采用炸藥作為地震勘探的震源。隨著勘探精度的提高，為了獲得完整的地震資料，炸藥震源時(shí)常被埋放在離村莊、城鎮(zhèn)僅幾十米的位置。由此給房屋結(jié)構(gòu)帶來(lái)的振動(dòng)危害是不可忽視的。房屋結(jié)構(gòu)對(duì)爆炸地震的振動(dòng)響應(yīng)有著重要的研究?jī)r(jià)值。許多國(guó)家對(duì)于房屋結(jié)構(gòu)的振動(dòng)都建立了相應(yīng)的標(biāo)準(zhǔn)[1-3],這些標(biāo)準(zhǔn)有一個(gè)共同點(diǎn)：對(duì)不同頻率的振動(dòng)幅度有不同的允許閾值。Siskind[4-6]通過(guò)監(jiān)測(cè)采礦爆破周圍的房屋振動(dòng)發(fā)現(xiàn)，在接近房屋結(jié)構(gòu)固有頻率的頻段內(nèi)，房屋結(jié)構(gòu)上的振動(dòng)相對(duì)于臨近其地基的地面振動(dòng)有不同程度的放大。Pijush[7]指出結(jié)構(gòu)的振動(dòng)響應(yīng)取決于地表介質(zhì)的動(dòng)力學(xué)特性和房屋結(jié)構(gòu)的建造類型及自振特性。林鍵等[8]對(duì)3-4層房屋結(jié)構(gòu)的爆破振動(dòng)響應(yīng)進(jìn)行了測(cè)試，發(fā)現(xiàn)隨著樓層的增加，振動(dòng)速度峰值(PPV)會(huì)不同程度的增加。上述研究都是在采礦爆破的大前提下形成的，炸藥藥量通常在幾百公斤以上，且爆炸的位置距離房屋結(jié)構(gòu)在幾公里外。而作為地震探勘用的炸藥震源藥量通常小于10 kg，炸藥距離房屋較近，對(duì)于這類情況下房屋結(jié)構(gòu)的振動(dòng)響應(yīng)研究較少。此外，研究抗震性能好的房屋結(jié)構(gòu)和開(kāi)發(fā)減震材料是建筑學(xué)者始終關(guān)心的課題[9]。碳纖維板加固技術(shù)對(duì)于房屋結(jié)果的減震效果是十分顯著的[10-11],對(duì)爆炸地震效應(yīng)的減震效果的研究還較少。

基于上述原因，對(duì)一新建的二層鋼筋混凝土-水泥磚墻的房屋結(jié)構(gòu)進(jìn)行爆炸振動(dòng)測(cè)試試驗(yàn)，利用希爾伯特-黃變換(HHT)[12-13]對(duì)房屋結(jié)構(gòu)上不同位置出的振動(dòng)信號(hào)進(jìn)行頻譜分析。獲得了房屋結(jié)構(gòu)的對(duì)爆炸地表振動(dòng)的響應(yīng)特性，同時(shí)對(duì)比研究了碳纖維板加固對(duì)建筑結(jié)構(gòu)的減振效果。

1 振動(dòng)測(cè)試試驗(yàn)

1.1 試驗(yàn)現(xiàn)場(chǎng)條件

二層鋼筋混凝土-水泥磚墻房屋結(jié)構(gòu)建在平坦空曠的砂質(zhì)土上，如圖1。炸藥選用TNT，藥量為10 kg，炸藥埋覆于距離房屋結(jié)構(gòu)20 m，30 m與40 m的土中，深度為2 m。試驗(yàn)順序是從40m處到20m處，目的是為了防止爆炸后砂土成坑造成下一次試驗(yàn)產(chǎn)生較大的誤差。臨近房屋結(jié)構(gòu)地基的地面上振動(dòng)速度的測(cè)試系統(tǒng)選用美國(guó)勞雷公司的NZXP24型高精度地震儀，該儀器測(cè)試動(dòng)態(tài)范圍為144 dB。而在房屋結(jié)構(gòu)上的振動(dòng)測(cè)試選擇加速度傳感器(國(guó)內(nèi)702所生產(chǎn))，量程為10 g，靈敏度為24 PC/g，每個(gè)加速度傳感器外均附加鋼制的保護(hù)殼。傳感器的位置如圖2所示，將墻上的加速度傳感器命名為W1-W4，將鋼筋混凝土柱上的加速度傳感器命名為P1，將貼近房屋地基地面上的速度傳感器命名為G1。需要重點(diǎn)說(shuō)明的是：圖2中右側(cè)的墻W2和W4是通過(guò)碳纖維板加固的墻，而W1和W3是普通的磚墻。

圖1 二層水泥磚墻房屋結(jié)構(gòu)圖

圖2 傳感器位置及編號(hào)

1.2 振動(dòng)幅值特性

房屋結(jié)構(gòu)上振動(dòng)信號(hào)的峰值振動(dòng)速度與響應(yīng)的地面G1的峰值振動(dòng)速度如表1所示。可以看到，房屋結(jié)構(gòu)上的振動(dòng)速度峰值相對(duì)于G1來(lái)說(shuō)有不同程度的放大效果。

表1 峰值振動(dòng)速度(mm/s)

2 振動(dòng)信號(hào)分析

2.1 希爾伯特-黃變換(HHT)

傅里葉變換(FT)和小波分析(Wavelet)一直是運(yùn)用最廣泛的分析手段[14-15]。但二者有一個(gè)共同的缺陷：必須預(yù)先設(shè)定的先驗(yàn)基底函數(shù)。希爾伯特-黃變換(HHT)是近些年來(lái)信號(hào)分析研究領(lǐng)域的一個(gè)重大突破，相對(duì)比FT和Wavelet，HHT是一種自適應(yīng)的分析方法，不需要先驗(yàn)基底函數(shù)[12-13]；對(duì)信號(hào)的處理有更強(qiáng)的局部特性，能更好地表征信號(hào)的強(qiáng)間斷特性。尤其對(duì)于非線性系統(tǒng)產(chǎn)生的隨機(jī)非穩(wěn)態(tài)信號(hào)處理有很大的優(yōu)勢(shì)。HHT的本質(zhì)是對(duì)信號(hào)的平穩(wěn)化處理，主要包括了兩大模塊：一是對(duì)時(shí)程信號(hào)的經(jīng)驗(yàn)?zāi)B(tài)分解(Empirical Mode Decomposition, EMD)，模態(tài)分解后的數(shù)據(jù)序列被稱作固有模態(tài)函數(shù)(Intrinsic Mode Function, IMF)；二是對(duì)每個(gè)IMF進(jìn)行Hilbert變換，獲取其時(shí)間-頻率-能量三維特性。

2.2 測(cè)試信號(hào)的EMD分解

假定原始時(shí)程信號(hào)S(t)，找到所有極值點(diǎn)，用三次樣條曲線分別對(duì)極大值點(diǎn)、極小值點(diǎn)進(jìn)行差值獲取信號(hào)的上包絡(luò)線Smax和下包絡(luò)線Smin。對(duì)Smax和Smin求均值得到均值線M1(t)，再用原始信號(hào)S(t)減掉M(t)獲得第一個(gè)IMF分量C1(t)。

(1)

C1(t)=S(t)-M1(t)

(2)

判斷C1(t)是否滿足IMF函數(shù)的兩個(gè)要求[12]：①C1(t)中的數(shù)據(jù)極值點(diǎn)數(shù)量與曲線穿過(guò)零值點(diǎn)的數(shù)量相差必須≤1； ②C1(t)中任意數(shù)據(jù)點(diǎn)由Smax和Smin兩條包絡(luò)線確定的平均值為0。如果C1(t)滿足上述條件，則視C1(t)為第一個(gè)IMF函數(shù)。通常情況下，C1(t)還不能滿足上述條件，此時(shí)將C1(t)設(shè)為原始信號(hào)，重復(fù)上述算法：

C11(t)=C1(t)-M11(t)

(3)

這類算法有可能進(jìn)行k次，直到C1k(t)滿足IMF函數(shù)要求為止，才可以獲得信號(hào)的第一個(gè)IMF函數(shù)IMF1(t)：

C1k(t)=C1(k-1)(t)-M1(k-1)(t)

(4)

IMF1(t)=C1k(t)

(5)

這類重復(fù)算法有一個(gè)終止的判斷準(zhǔn)則，可以利用兩個(gè)連續(xù)處理結(jié)果之間的標(biāo)準(zhǔn)差R的值做判據(jù)：

(6)

通常R取值在0.2～0.3。第一個(gè)IMF函數(shù)產(chǎn)生之后，將從原始信號(hào)S(t)中減去IMF1(t)獲得的差值再重復(fù)上面的算法依次獲得余下的IMFi(t) (i=2,3,4…)。

利用上述方法對(duì)30 m處的試驗(yàn)時(shí)，G1、P1和W1的振動(dòng)信號(hào)進(jìn)行EMD分解，可以得到如圖3-圖5所示的結(jié)果。需要說(shuō)明的是P1和W1測(cè)到的原始信號(hào)為加速度信號(hào)，在進(jìn)行EMD分解之前利用數(shù)值積分程序?qū)⒓铀俣刃盘?hào)轉(zhuǎn)變?yōu)樗俣刃盘?hào)。每組信號(hào)的6個(gè)IMF分量IMF1-IMF6，包含了不同的頻域特征尺度，隨著分解的先后順序，頻率逐漸降低。分量的振幅越大，說(shuō)明該頻率段的信號(hào)能量越強(qiáng)。就地面振動(dòng)信號(hào)G1而言，IMF1-IMF3的幅值大，說(shuō)明在相對(duì)較高的頻率段內(nèi)振動(dòng)能量較高。而對(duì)于房屋結(jié)構(gòu)上振動(dòng)信號(hào)，IMF4-IMF6的幅值大，說(shuō)明在較低的頻率段內(nèi)振動(dòng)能量較高。此外P1和W1上的振動(dòng)幅度要高于G1，這個(gè)現(xiàn)象說(shuō)明房屋結(jié)構(gòu)對(duì)爆炸產(chǎn)生的地面振動(dòng)有放大效應(yīng)，定量的放大倍率和頻率特性會(huì)在下一節(jié)作詳細(xì)分析。

圖3 30 m處試驗(yàn)地面振動(dòng)信號(hào)G1及其IMF分量

2.3 振動(dòng)信號(hào)IMF分量的Hilbert頻譜分析

應(yīng)用EMD獲取的IMF分量再進(jìn)行Hilbert變換：

(7)

式中：PV代表柯西主值。運(yùn)用EMD算法獲取的IMF分量IMF(t) 和Hilbert變換獲取的變換譜H[S(t)]可以構(gòu)造解析信號(hào)z(t)，該解析信號(hào)的幅值a(t)，也稱包絡(luò)幅值，其表達(dá)式為：

z(t)=IMF(t) +iH[IMF(t)]=a(t)eiθ(t)

(8)

(9)

Hilbert時(shí)頻譜的表達(dá)式如下：

(10)

通過(guò)Hilbert變換還可獲得信號(hào)的瞬時(shí)頻率ωin：

(11)

(12)

將房屋結(jié)構(gòu)上振動(dòng)信號(hào)分解出的IMF分量進(jìn)行Hilbert變換，并求解頻率域內(nèi)與地面G1的比值，如圖6與圖7。

圖6 混凝土柱P1-P3振動(dòng)相對(duì)地面G1在頻率域內(nèi)的放大效果

不難發(fā)現(xiàn)，鋼筋混凝土柱和水泥磚墻，在靠近其固有頻率的一段頻率域內(nèi)，振動(dòng)信號(hào)的幅度相對(duì)于G1均有明顯的放大效果(Amplification Ratio > 1)，炸藥爆炸距離從40 m減小到20 m，能產(chǎn)生放大效果的頻率段變窄。這個(gè)現(xiàn)象產(chǎn)生的主要原因是離炸藥近的地面振動(dòng)頻率較高。對(duì)于混凝土柱，不同高度處振動(dòng)相對(duì)地面的放大倍率會(huì)有明顯不同，位置越高其放大倍率越大，P1>P2>P3。隨著炸藥更加靠近房屋結(jié)構(gòu)，P1，P2和P3在頻率域內(nèi)的振動(dòng)放大倍率都會(huì)降低，峰值放大倍率由11.6，10.1和6分別減小到4.8, 2.4和2.1，但始終維持在8-10 Hz；能夠產(chǎn)生放大效應(yīng)的頻率段從6-30 Hz變窄到6-18 Hz。對(duì)于墻體的振動(dòng)，處在第二層的W1和W2相對(duì)于G1的放大倍率要大于第一層的W3和W4，進(jìn)一步說(shuō)明了房屋結(jié)構(gòu)上位置越高振動(dòng)幅度越大。隨著炸藥更加靠近房屋結(jié)構(gòu)，W1,W2,W3和W4在頻率域內(nèi)的峰值放大倍率從12.3, 4.5, 8,2和4.4降低為3.8, 1.6,1.7和1.03。對(duì)比相同高度的墻體，可以發(fā)現(xiàn)由于W2和W4有碳纖維板的加固，相對(duì)于G1的放大倍率要分別遠(yuǎn)小于相同高度的W1和W3，其中W2約為W1的36.6%-49.8%，W4約為W3的34.1%-53.7%。事實(shí)上，在爆炸距離為30 m和20 m時(shí)，在第一層的墻體W4相對(duì)于G1幾乎沒(méi)有放大，而沒(méi)有碳纖維板加固的W3還能觀測(cè)到的放大倍率為4.2和1.7，這體現(xiàn)了碳纖維板有明顯的減震效果。值得注意的是，碳纖維板并未改變墻體的固有頻率，當(dāng)炸藥爆炸距離從40 m減小到20 m，墻體能產(chǎn)生放大效果的頻率段變窄，但出現(xiàn)峰值放大倍率的頻率維持在7-12 Hz。普通的水泥磚墻在受到地震作用時(shí)，其材料延性差，抗剪切能力低，尤其是構(gòu)造不規(guī)范的房屋在輕微振動(dòng)下便發(fā)生破壞；碳纖維板加固能極大的提高墻體的抗沖擊和抗剪切能力。

圖7 墻體W1-W4振動(dòng)相對(duì)地面G1在頻率域內(nèi)的放大效果

3 結(jié) 論

本文對(duì)新建的二層鋼筋混凝土-磚墻結(jié)構(gòu)進(jìn)行了爆炸振動(dòng)響應(yīng)測(cè)試試驗(yàn)，以房屋結(jié)構(gòu)上的振動(dòng)相對(duì)地基附近地面振動(dòng)的放大倍率為指標(biāo)，利用HHT分析了房屋結(jié)構(gòu)爆炸振動(dòng)的頻譜響應(yīng)特征。發(fā)現(xiàn)在靠近房屋結(jié)構(gòu)自振頻率的頻率域內(nèi)，房屋結(jié)構(gòu)呈現(xiàn)不同程度的振動(dòng)放大效應(yīng)。高度較高的位置處，包括鋼筋混凝土柱和墻體，它們?cè)谟騼?nèi)的峰值放大倍率要高于較低位置。而炸藥與房屋結(jié)構(gòu)的距離越小會(huì)造成峰值放大倍率的減小以及能夠產(chǎn)生放大效應(yīng)的頻帶變窄。對(duì)比分析了相同高度處經(jīng)碳纖維板加固后的墻體和普通墻體的振動(dòng)特性，結(jié)果表明碳纖維板加固對(duì)爆炸地震效應(yīng)有明顯的減振效果，有碳纖維板加固的墻體振動(dòng)幅度要低于相同高度普通墻體振動(dòng)的1/2。

[1]Siskind D E, Stagg M S, Kopp J W, et al. Structure reponse and damage produced by ground vibration from surface mine blasting[R]. US Bur Mines Rep Invest 8507, 1980.

[2]German Institute of Standards. Vibration of building-effects on structures[R]. DIN 4150, 1986,(3).

[3]Singh P K, Roy M P. Standardisation of blast vibration damage threshold for the safety of residential structures in mining areas[R]. Central Mining Research Institute (India), R&D report, 2006.

[4]Siskind D E, Stagg M S, Kopp J W, et al. Structure response and damage produced by ground vibration from surface mine blasting[R]. United States Bureau of Mines, RI 8507. 1980.

[5]Siskind D E, Stachura V J, Nutting M J. Low-Frequency vibration produced by surface mine blasting over abandoned underground mines[R]. United States Bureau of Mines,RI 9078. 1987.

[6]Dowding C H. Blast vibration monitoring and control[R]. Englewood Cliffs, NJ: Prentice-Hall, 1985.

[7]Pijush P R. Characteristics of ground vibrations and structural response to surface and underground blasting[J]. Geotechical and Geological Engineering, 1998, 16: 151-166.

[8]林鍵, 林從謀, 林麗群. 爆破振動(dòng)載荷作用下3～4層房屋結(jié)構(gòu)響應(yīng)測(cè)試研究[J]. 振動(dòng)與沖擊, 2010 29(3): 48-66.

LIN Jian, LIN Cong-mou, LIN Li-qun. Test research on structural response to blasting vibration loading of 3～4 layers building [J]. Journal of vibration and shock, 2010 29(3): 48-66.

[9]陳國(guó), 單波, 肖巖. 輕型竹結(jié)構(gòu)房屋抗震性能的試驗(yàn)研究[J]. 振動(dòng)與沖擊, 2011, 39(10): 136-142.

CHEN Guo, SHAN Bo, XIAO Yan. A seismic performance tests for a light glubam house [J]. Journal of vibration and shock, 2011, 39(10): 136-142.

[10]John J. Reinforcing masonry walls with composite materials [D]. Stanford University, 2000.

[11]Attard C K, Mario M. Modeling of unreinforced masonry walls under shear and compression [J]. Engineering structures, 2007, 29(9): 289-296.

[12]Huang N E. Introduction to the Hilbert-Huang Transform and its related mathematical problems[M]. Interdisciplinary mathematics, 2005.

[13]Huang N E, Long S R, Shen Z. The mechanism for frequency downshift in nonlinear wave evolution [J].Advances in applied mechanics. 1996, 32: 59-111.

[14]Grossmann P, Grossmann A, Morlet J. Cycle-octave and related transforms in seismic signal analysis[J]. Geoexploration,1984, 23: 85-102.

[15]林大超, 白春華. 爆炸地震效應(yīng)[M]. 北京：地質(zhì)出版社，2007.