纖維爆炸索水下爆炸聲信號(hào)特征的小波分析

賈 虎，沈兆武

(1.中國(guó)科學(xué)技術(shù)大學(xué) 近代力學(xué)系，合肥 230027;2.南陽(yáng)師范學(xué)院 土木建筑工程學(xué)院，南陽(yáng) 473061)

含能材料水下爆炸產(chǎn)生的脈沖波是一個(gè)無(wú)指向性的高功率聲源，可以傳播很遠(yuǎn)的距離，同時(shí)會(huì)引起很強(qiáng)的混響效應(yīng)，因此含能材料有望成為一種新的噪聲干擾源［1]。

纖維爆炸索作為一種新型低能量導(dǎo)爆索，由于其側(cè)向約束力強(qiáng)，藥芯密度高，裝藥均勻，線裝藥密度低，爆速穩(wěn)定，質(zhì)量有保證，因此具有十分廣泛的發(fā)展前景［2]。為了探求合適的水聲干擾源，實(shí)現(xiàn)對(duì)水下武器的干擾，我們對(duì)纖維爆炸索進(jìn)行了水下爆炸實(shí)驗(yàn)研究，獲得了水下爆炸信號(hào)。

傳統(tǒng)水下爆炸壓力信號(hào)分析時(shí)往往采用時(shí)域分析方法，但是由于爆炸壓力信號(hào)持續(xù)時(shí)間短、突變快，屬于非平穩(wěn)信號(hào)，信號(hào)的結(jié)構(gòu)及其頻譜都是時(shí)變的，與Fourier變換建立的理論基礎(chǔ)有不一致［3，4]。為了更好地分析爆炸壓力信號(hào)的特點(diǎn)，一些研究者［5，6]采用小波變換處理水聲信號(hào)，并取得了比較理想的結(jié)果。本文采用小波變換方法對(duì)纖維爆炸索水下爆炸水聲特性進(jìn)行研究，并進(jìn)行小波重構(gòu)信號(hào)能量分布研究，以期有效分析爆炸聲信號(hào)特征，獲得有效水聲干擾源，為水下對(duì)抗提供新的手段。

1 小波分析原理

小波分析作為目前較為先進(jìn)的信號(hào)分析手段，已在包括爆破在內(nèi)的許多行業(yè)領(lǐng)域得到推廣應(yīng)用，其最大優(yōu)點(diǎn)是能對(duì)信號(hào)的“細(xì)則”部分像“顯微鏡”一樣同時(shí)在時(shí)域和頻域給予詳細(xì)的考察［7]。小波分析是一種窗口大小固定但形狀可改變，時(shí)間窗和頻率窗都可改變的時(shí)頻局部化分析方法。即在低頻部分具有較高的頻率分辨率和較低時(shí)間分辨率，在高頻部分具有較高時(shí)間分辨率和較小頻率分辨率［8]。由于小波分析在時(shí)域和頻域上都有良好的局部化性質(zhì)，并能對(duì)不同的頻率成分提供不同的分析分辨率，這使其性能大大優(yōu)于Fourier變換。在多分辨率條件下，采用二進(jìn)制小波時(shí)，信號(hào)被分解成低頻和高頻兩部分。在分解中，低頻部分失去的信號(hào)由高頻部分捕獲。下一層分解中再將分解出的低頻部分分解為低頻和高頻部分，低頻部分失去的信號(hào)仍有高頻部分捕獲，如此循環(huán)，最終完成更深層次的分解。

2 水下爆炸壓力試驗(yàn)

水下爆炸實(shí)驗(yàn)采用的水池由鋼板加工而成，壁厚30 mm，幾何尺寸為φ5 m ×5 m，水深4.4 m，水池立于地表并固定在減震設(shè)施上，如圖1所示。測(cè)試系統(tǒng)主要有美國(guó)PCB公司生產(chǎn)的W138A25長(zhǎng)型電氣石水下激波壓力傳感器(ICP)，PCB公司的482A22型恒流源和美國(guó)泰克DPO7054型示波器，信號(hào)采樣頻率為5 MHz。ICP傳感器為體積敏感型傳感器，安裝時(shí)不需將傳感器正對(duì)爆心。傳感器本身自帶的微型放大器，將高阻抗電荷轉(zhuǎn)換成＜100Ω的低阻抗電壓輸出，信號(hào)可用普通電纜傳送相當(dāng)遠(yuǎn)(幾百米)而無(wú)顯著信號(hào)衰減和噪聲干擾，得到的沖擊波p－t曲線光滑、平整。

試驗(yàn)樣品為纖維爆炸索，如圖2所示，外徑為1.5 mm，裝藥量為黑索金1.5 g/m，纖維爆炸索長(zhǎng)2.0 m。纖維爆炸索中點(diǎn)處固定在水池的中軸線水下2.9 m的位置，傳感器固定在直徑為1 mm的尼龍線上，尼龍線的下端掛配重垂直放入水下，傳感器入水深度2.9 m，與纖維爆炸索中點(diǎn)處放在相同的水深，距離纖維導(dǎo)爆炸索1.0 m。纖維爆炸索采用雷管一端起爆，為防止雷管爆轟影響水中流場(chǎng)，將雷管置于專門加工的抗爆器內(nèi)，有效削弱雷管影響。實(shí)驗(yàn)裝置平面布置如圖3所示。

圖4為測(cè)試得到水下爆炸沖擊波壓力和氣泡脈動(dòng)壓力曲線，圖5為沖擊波壓力時(shí)間曲線。

從圖5中可以看出，沖擊波壓力在3μs時(shí)間內(nèi)從零上升到峰值壓力Pm，數(shù)值為5.21 MPa，隨后近似成指數(shù)衰減。將壓力從Pm衰減到Pm/e(其中e為自然對(duì)數(shù)，取值為2.718 3)所需時(shí)間定義為特征時(shí)間θ，經(jīng)計(jì)算可得纖維爆炸索的特征時(shí)間θ為42μs，在θ之后，沖擊波的衰減要比θ時(shí)間內(nèi)緩慢得多。纖維爆炸索第一次氣泡脈動(dòng)周期約為15 ms，雖然氣泡脈動(dòng)壓力峰值不到?jīng)_擊波壓力峰值的15%，但是持續(xù)時(shí)間比沖擊波峰值壓力持續(xù)時(shí)間長(zhǎng)得多。

3 水下爆炸聲壓

纖維爆炸索水下爆炸沖擊波的瞬時(shí)壓力峰值非常高，持續(xù)時(shí)間非常短。圖6為水下爆炸實(shí)驗(yàn)獲得的纖維爆炸索水下爆炸聲壓級(jí)時(shí)域圖，從圖中可以看出，纖維爆炸索水下爆炸聲壓級(jí)完全可以達(dá)到200 dB以上，在持續(xù)15 ms以后時(shí)纖維爆炸索的爆炸聲壓級(jí)仍在200 dB上下擺動(dòng)，說(shuō)明纖維爆炸索水下爆炸具有很強(qiáng)的聲功率，這是其有望成為水下干擾聲源的顯著特點(diǎn)之一。

纖維爆炸索水下爆炸，會(huì)產(chǎn)生氣泡脈動(dòng)，并產(chǎn)生大量氣泡。由于水中介質(zhì)的不均勻性以及水下爆炸產(chǎn)生的氣泡脈動(dòng)現(xiàn)象，當(dāng)纖維爆炸索水下爆炸產(chǎn)生的爆炸聲源在傳播過(guò)程中與這些介質(zhì)相遇，會(huì)產(chǎn)生漫反射和散射，同時(shí)還要考慮到纖維爆炸索水中爆炸傳播的多路效應(yīng)，在傳播中產(chǎn)生迭加而出現(xiàn)混響效應(yīng)，提高水下爆炸脈沖壓力波的寬度和持續(xù)時(shí)間，從而達(dá)到水聲對(duì)抗的作用。從圖6中可以看出，纖維爆炸索在水下爆炸時(shí)會(huì)產(chǎn)生很強(qiáng)的混響效應(yīng)，且能保持在200 dB左右，持續(xù)一段時(shí)間。

4 纖維爆炸索水下爆炸聲信號(hào)小波分析

為了能夠獲得爆炸水聲信號(hào)的細(xì)節(jié)性質(zhì)，提供頻率成分的時(shí)間局部信息，采用一維離散小波變換對(duì)纖維爆炸索水下爆炸聲信號(hào)進(jìn)行變換，以得到信號(hào)的細(xì)部特征。

4.1 小波分析

小波變換具有時(shí)間分辨率，可以分析出不同頻率成分分量的時(shí)間衰減規(guī)律，這就為我們分析和提取纖維爆炸索水下爆炸聲信號(hào)的有效信息提供了幫助。

小波分析在實(shí)踐應(yīng)用中遇到的很重要的問(wèn)題就是小波基的選取，因?yàn)橛貌煌男〔ɑ治鐾粏?wèn)題會(huì)得到不同的結(jié)果。目前在非平穩(wěn)振動(dòng)信號(hào)分析中運(yùn)用比較多的是 Daubechies8 小波基函數(shù)［9，10]。

試驗(yàn)設(shè)置的信號(hào)采樣頻率為5 MHz，根據(jù)采樣定理，則其奈奎斯特(Nyquist)頻率為2.5 MHz。根據(jù)小波分析原理，采用圖5給出的水下爆炸沖擊波壓力試驗(yàn)曲線，利用Daubechies8小波函數(shù)對(duì)纖維爆炸索水下爆炸信號(hào)進(jìn)行離散小波變換，用db8小波函數(shù)將爆炸信號(hào)分解到第9層，對(duì)應(yīng)的最低頻帶為0～4.882 812 5。

根據(jù)圖5給出的爆炸壓力試驗(yàn)曲線，采用db8小波進(jìn)行分解后的重構(gòu)信號(hào)以及實(shí)測(cè)信號(hào)與重構(gòu)信號(hào)相對(duì)誤差分布如圖7所示。

從圖7可以看出實(shí)測(cè)信號(hào)與重構(gòu)信號(hào)之間相對(duì)誤差非常小，二者之間具有高度的一致性，說(shuō)明選取的小波基能真實(shí)反映信號(hào)情況，適合處理纖維爆炸索水下爆炸的短時(shí)非平穩(wěn)隨即信號(hào)問(wèn)題。

圖7 纖維爆炸索水下爆炸重構(gòu)信號(hào)及相對(duì)誤差分布圖Fig.7 Distribution map of underwater explosion signal reconstruction and relative error of fiber-based detonating cord

圖8 基于db8小波基的纖維爆炸索水下爆炸信號(hào)小波分層重構(gòu)信號(hào)圖Fig.8 Underwater explosion wavelet hierarchical reconstruction signal graph of fiber-based detonating cord based db8

表1 小波分解頻帶表Tab.1 Corresponding frequency band of wavelet decomposition

圖8為基于db8小波基的小波重構(gòu)信號(hào)，他們分別對(duì)應(yīng)10個(gè)頻率帶，圖中的a9及d1～d9為小波分量，其中a9為低頻分量，d1～d9為高頻分量。表1為采用Daubechies8小波函數(shù)進(jìn)行分解時(shí)對(duì)應(yīng)的頻帶表。

從圖8可以看出基于小波變換的分析方法可以獲得沖擊波壓力在各頻帶的分布和衰減信息。

圖8中，a9為小波分解重構(gòu)信號(hào)的第9層逼近信號(hào)，d1～d9為1到9層的細(xì)節(jié)信號(hào)。a9中24.6μs時(shí)幅值最大，且明顯高于各層細(xì)節(jié)信號(hào)的幅值，能流高度集中。d7～d9及a9的波形信號(hào)比較接近纖維爆炸索水下爆炸沖擊波信號(hào)，振動(dòng)幅值較大，能量主要集中在這部分，這說(shuō)明沖擊波壓力信號(hào)能量主要集中在低頻段。d1～d4的振動(dòng)幅值波動(dòng)較大，說(shuō)明纖維爆炸索水下爆炸沖擊波信號(hào)包含的頻率成分比較豐富。

4.2 小波變換后沖擊波信號(hào)各頻率帶能量統(tǒng)計(jì)

由小波分解得到了各頻帶范圍內(nèi)的小波分量，各頻帶內(nèi)的壓力分量仍為關(guān)于時(shí)間變化的曲線。圖8中將水下爆炸信號(hào)分析到了第九層，各層對(duì)應(yīng)的能量為:

式中:Ei為第i頻率帶信號(hào)對(duì)應(yīng)的能量;E0為分析信號(hào)的總能量;Si為第i頻率帶的小波分解信號(hào);xij為信號(hào)Si的離散點(diǎn)幅值;其中 i=1，2，…，10;j=1，2，…，m，m為信號(hào)的離散采樣點(diǎn)。

各頻率帶能量占被分析信號(hào)總能量的比例為:

為了進(jìn)一步分析水下爆炸沖擊波壓力信號(hào)各頻帶能量分布情況，根據(jù)小波變換分層重構(gòu)信號(hào)可以得到在不同頻帶上水下爆炸相對(duì)能量分布情況，如圖9所示。

從圖9中可以看出，在前四個(gè)頻帶內(nèi)，沖擊波壓力相對(duì)能量最大，占總能量的比例高達(dá)96.2%，說(shuō)明沖擊波壓力主要集中在40 kHz以下，尤以5 kHz以下頻帶能量最高。纖維爆炸索能量主要集中在40 kHz以下，如果作為水聲干擾源，完全可以覆蓋各類水聲器材的各工作頻段，或許能在水聲對(duì)抗中發(fā)揮重要作用。

圖9 各頻帶上相對(duì)能量分布情況Fig.9 Relative energy distribution of each band

5 結(jié)論

通過(guò)對(duì)纖維爆炸所水下爆炸壓力信號(hào)的分析，可得出以下結(jié)論:

(1)纖維爆炸索在水下爆炸時(shí)會(huì)產(chǎn)生很強(qiáng)的混響效應(yīng)，且能保持在200 dB左右，持續(xù)一段時(shí)間，說(shuō)明纖維爆炸索水下爆炸具有很強(qiáng)的聲功率，這一特點(diǎn)使其有望成為水下干擾聲源;

(2)采用小波變換的方法對(duì)纖維爆炸索水下爆炸壓力信號(hào)進(jìn)行分析，獲得了纖維爆炸索沖擊波壓力在各頻帶上的壓力－時(shí)間細(xì)節(jié)信號(hào)，發(fā)現(xiàn)纖維爆炸索水下爆炸沖擊波信號(hào)包含的頻率成分比較豐富，能量主要集中在低頻部分;

(3)為了進(jìn)一步分析水下爆炸沖擊波壓力信號(hào)各頻帶能量分布情況，根據(jù)小波變換分層重構(gòu)信號(hào)得到在不同頻帶上水下爆炸相對(duì)能量分布情況，發(fā)現(xiàn)聲信號(hào)能量主要集中在40 kHz以下，尤以5 kHz以下頻帶能量最高。纖維爆炸索能量主要集中在40 kHz以下，如果作為水聲干擾源，完全可以覆蓋各類水聲器材的各工作頻段，或許能在水聲對(duì)抗中發(fā)揮重要作用。

［1] 潘正偉，焦善武，顧曉輝.水下爆炸－高功率寬頻帶的水聲干擾源［J] .南京理工大學(xué)學(xué)報(bào)，1999，23(6):507－509.

［2] 梅 群.低能量導(dǎo)爆索關(guān)鍵技術(shù)及應(yīng)用研究［D] .合肥:中國(guó)科學(xué)技術(shù)大學(xué)，2007.

［3] 溫華兵，張 健，尹 群，等.水下爆炸船艙沖擊響應(yīng)時(shí)頻特征的小波包分析［J] .工程力學(xué)，2008，25(6):199－203.

［4] 李夕兵，張義平，劉志祥，等.爆破震動(dòng)信號(hào)的小波分析與HHT變換［J] .爆炸與沖擊，2005，25(6):528－535.

［5] 杜志鵬，汪 玉，楊 洋，等.艦艇水下爆炸沖擊信號(hào)擬合及應(yīng)用［J] .振動(dòng)與沖擊，2010，29(3):182 －184.

［6] 晏俊偉，龍 源，方 向，等.基于小波包變換的爆破地震波時(shí)頻特征提取及分析［J] .振動(dòng)與沖擊，2007，26(4):25－29.

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

［8] 李夕兵，凌同華，張義平.爆破震動(dòng)信號(hào)分析理論與技術(shù)［M] .北京:科學(xué)出版社，2009.

［9] 胡昌華，張軍波，夏 軍.基于MATLAB的系統(tǒng)分析與設(shè)計(jì)—小波分析［M] .西安:西安電子科技大學(xué)出版社，2004.

［10] 李夕兵，凌同華，張義平.爆破震動(dòng)信號(hào)分析理論與技術(shù)［M] .北京:科學(xué)出版社，2009.