魏永星，于金花，李　琦，周　瑩，?！≌埽Ｖ救A，高　超

（國家海洋技術(shù)中心，天津　300112）

實測海洋環(huán)境噪聲數(shù)據(jù)譜級特性研究

魏永星，于金花，李琦，周瑩，常哲，牛志華，高超

（國家海洋技術(shù)中心，天津300112）

通過對噪聲測量潛標系統(tǒng)測得的94 d海洋環(huán)境噪聲數(shù)據(jù)及相應的氣象數(shù)據(jù)進行聯(lián)合時頻、相關(guān)性、統(tǒng)計特性分析，得到試驗海區(qū)風速、降雨等自然現(xiàn)象與噪聲譜級之間的關(guān)系。分析結(jié)果表明，海洋環(huán)境噪聲譜級與風速變化在低頻段相關(guān)性較低，高頻段相關(guān)性較高，同時1 kHz以上頻段譜級值分布接近呈正態(tài)分布；降雨對環(huán)境噪聲的影響主要分布在500 Hz以上頻段。

海洋環(huán)境噪聲；風速；相關(guān)性；統(tǒng)計特性；降雨

海洋環(huán)境噪聲是永恒存在的聲場，包含了豐富的水體、海面、海底等環(huán)境信息，研究者可以通過充分研究海洋環(huán)境噪聲來獲得海洋環(huán)境的物理特性及豐富的海洋聲學信息，還可以利用海洋環(huán)境噪聲的寬頻帶和任意時間存在的特性提取相關(guān)的海洋環(huán)境參數(shù)。

根據(jù)Knudsen，Urick，Wenz與Carey等的海洋環(huán)境噪聲文獻，可以將不同環(huán)境噪聲源分為相互覆蓋的三段：（1）低頻段1～100 Hz，主要來源于潮、浪、涌的壓力脈動，大尺度湍流以及遠處的風暴、地震等；（2）10～500 Hz，總體來說較為平緩，來源于遠處的航船，強弱與航運的頻繁程度有關(guān)；（3）500 Hz～25 kHz頻段，主要來源于風動海面。其中，風關(guān)噪聲、降雨噪聲是海洋環(huán)境的主要噪聲源之一，也是本文的重點研究對象。

本文以海洋環(huán)境噪聲測量潛標系統(tǒng)在某海域獲取的94 d海洋環(huán)境噪聲數(shù)據(jù)及相應的氣象資料為依托，通過時頻變換、相關(guān)性、統(tǒng)計特性分析等研究海洋環(huán)境噪聲譜級的時間變化特性以及風速、降雨等因素對海洋環(huán)境噪聲的影響。

1　資料處理方法

1.1海洋環(huán)境噪聲譜級

海洋環(huán)境噪聲級（NL）是用來衡量環(huán)境噪聲強弱的一個量。本文將原始數(shù)據(jù)進行分段處理得到噪聲譜級值，分段時保證FFT的頻率分辨率不大于1 Hz，并對每段數(shù)據(jù)進行加窗處理。因為實測數(shù)據(jù)量充足，在數(shù)據(jù)分段時并未進行重疊處理。多段數(shù)據(jù)之間進行功率譜平均后進行帶寬歸一化，即可得到該時間段的1/3倍頻程噪聲譜級值［1］，具體公式如下：

設有效噪聲信號為x（n），將其分為I段，每段長度為N，設第i段信號序列為xi（n），將xi（n）乘以窗函數(shù)w（n），其傅里葉變換為：

式中：i=1，2，……，I；k=0，1，2……，N-1。

則第i段修正（加窗）周期圖為：

則x（n）的功率譜估值P（k）為上述I個修正周期圖的線性平均：

式中：k=0，1，2，……，N-1。

1/3倍頻程噪聲譜級值如式（4）所示：

式中：NL（fi）為1/3倍頻程中心頻率fi處的噪聲譜級值；fh，fl分別為fi的上下限頻率對應的頻率位置；Δf為頻率分辨率；M（fi）為水聽器靈敏度。

1.2相關(guān)理論

兩變量之間的相關(guān)現(xiàn)象一般以相關(guān)系數(shù)來表示，基本公式如式（5）：

1.3統(tǒng)計理論

許多統(tǒng)計理論是建立在正態(tài)分布的基礎上，一般來說，資料達到一定的數(shù)量，會呈現(xiàn)正態(tài)分布。正態(tài)分布的指數(shù)方程式如式（6）所示：

式中：x為實測噪聲譜級值；σ為標準差；μ為期望值。

2　資料來源及處理流程圖

本文中使用的海洋環(huán)境噪聲數(shù)據(jù)來源于海洋環(huán)境噪聲測量潛標系統(tǒng)在某海域進行的94 d的連續(xù)測量。風速、降雨資料來源于美國國家環(huán)境預報中心（NCEP）再分析資料，空間分辨率為1°×1°。

本文資料分析流程如圖1所示。

圖1　資料分析流程圖

3　數(shù)據(jù)處理結(jié)果分析

3.1風速對環(huán)境噪聲譜級的影響

海面粗糙度是高頻段自然噪聲的噪聲源，與風速有關(guān)［2］。風造成的氣泡、空化以及海面激勵所產(chǎn)生的浪花會產(chǎn)生噪聲影響高頻段。如果能準確給出海洋環(huán)境噪聲級與風速的關(guān)系，可以通過海洋環(huán)境噪聲的實測數(shù)據(jù)來反演海面風速。1948年Knudson等［3］整理了第二次世界大戰(zhàn)期間海洋環(huán)境噪聲的海上測量數(shù)據(jù)，獲得以海況或者風力為參數(shù)的著名的Knudson譜，此后發(fā)現(xiàn)，風速與海洋環(huán)境噪聲譜級的相關(guān)性比海況的相關(guān)性更好。1964年P(guān)iggott［4］指出，在淺海（水深40 m以內(nèi)）噪聲譜級與風速的對數(shù)呈簡單的線性關(guān)系，1972年Crouch和Bur［5］將此結(jié)果推廣到5 000 m深海。

3.1.1噪聲譜級與風速時頻及相關(guān)分析 本文將潛標系統(tǒng)得到的海洋環(huán)境噪聲數(shù)據(jù)進行濾波、剔除異常值后做時頻分析，得到風速與海洋環(huán)境噪聲譜級的時間變化圖2所示。

圖2　94 d風速與海洋環(huán)境噪聲譜級時間變化圖

由圖2可以看出，在頻率較高即500 Hz以上時，譜級值與風速的相關(guān)性較好，頻率較低即時500 Hz以下時相關(guān)性較差，說明風速與高頻相關(guān)程度高。

仔細觀察風速變化，發(fā)現(xiàn)在6月22日出現(xiàn)高峰值，著重分析6月17日00:00至6月30日23:00共14 d數(shù)據(jù)，得到其風速與海洋環(huán)境噪聲譜級的時間變化如圖3所示。

圖3　14 d風速與海洋環(huán)境噪聲譜級時間變化圖

由圖3可以看出，海洋環(huán)境噪聲譜級值的低頻段與風速的變化趨勢不太相似，主要是因為低頻段主要噪聲源為遠處航船噪聲，風速的影響較小，出現(xiàn)的峰值也是由于航船引起的；而高頻處在6月22日出現(xiàn)高峰，升高了近20 dB，此時對應的風速高達15 m/s。

為定量分析風速對噪聲譜級值的影響程度，圖4給出了譜級與風速的相關(guān)系數(shù)曲線：

圖4　海洋環(huán)境噪聲譜級與風速的相關(guān)系數(shù)

由圖4可以看出，500 Hz以下的頻段譜級值與風速的相關(guān)性較低，低于0.45，分析原因是該頻段遠處船舶噪聲占主要因素；500 Hz以上高頻段相關(guān)性較高，高于0.45，分析原因是該頻段風關(guān)噪聲占主要因素。

3.1.2噪聲譜級的統(tǒng)計特性分析由于海洋環(huán)境噪聲是隨機的，本文對94 d的海洋環(huán)境噪聲譜級值進行統(tǒng)計特性分析得到風關(guān)海洋噪聲級的統(tǒng)計特性。

50 Hz，100 Hz，1 000 Hz，8 000 Hz四個中心頻點的概率分布如圖5所示。

圖5　50 Hz，100 Hz，1 000 Hz，8 000 Hz頻譜值概率分布圖

由圖5可以明顯看出，海洋環(huán)境噪聲譜級在1 kHz以上大致符合正態(tài)分布，因為這一頻段的海洋環(huán)境噪聲主要是由海面風引起的，測量時間跨度長，所產(chǎn)生的噪聲接近于正態(tài)分布，此次海試結(jié)果與文獻［6］的結(jié)果類似。在1 kHz以下頻段為非正態(tài)分布，原因是航船噪聲是該頻段的主要噪聲源，其時空特性變化決定了非正態(tài)分布的特性。

3.2降雨對環(huán)境噪聲譜級的影響

當海面受到雨滴拍打，甚至產(chǎn)生氣泡時，就會產(chǎn)生降雨噪聲，使海洋環(huán)境噪聲譜級發(fā)生變化，其強度有時比風關(guān)噪聲大得多，是影響海洋環(huán)境背景噪聲級的一個重要干擾源，會極大降低聲納的檢測能力，影響水聲通信的應用頻段和降低水聲設備的性能。不同的降雨強度產(chǎn)生的雨噪聲具有獨特的譜形狀。軍事上，海面降雨產(chǎn)生的高強度寬頻帶的雨聲能夠降低聲納的檢測能力。

此次海試時間跨度較長，降雨信息十分豐富，故分析結(jié)果有較強的實用價值。本文統(tǒng)計該次海試不同降雨情況下的海洋環(huán)境噪聲譜級值，得到降雨與噪聲譜級之間的關(guān)系，如圖6所示。

圖6　降雨與噪聲譜級之間的關(guān)系

由圖6（a）可以看出，在500 Hz以下的低頻段降雨影響程度不大，原因為該低頻段船舶噪聲占主要因素；500 Hz以上高頻段不同的降雨影響程度變化較大；陣雨比晴朗時高10 dB左右，中雨比陣雨高5 dB左右，暴雨時噪聲譜近于“白噪聲”，比晴朗時高出約25 dB。此次分析結(jié)果與Heindsman等人在長島海峽觀測到的降雨自然噪聲譜結(jié)果（如圖6（b）所示）大致符合，譜線形狀較相似。

建立降雨強度與噪聲譜級之間的定量關(guān)系，不僅有助于掌握海洋環(huán)境噪聲源的相關(guān)知識，而且有利于實現(xiàn)利用海洋環(huán)境噪聲譜級的變化反演海上降雨的技術(shù)。

4　結(jié)論

本文通過對94 d連續(xù)海洋環(huán)境噪聲譜級值及風速、降雨等自然現(xiàn)象的時頻變換及相關(guān)性、統(tǒng)計特性分析得到以下結(jié)論：

（1）海洋環(huán)境噪聲譜級與風速變化在500 Hz以上高頻段相關(guān)性較高，500 Hz以下低頻段相關(guān)性較低，同時1 kHz以上譜級值分布大致呈正態(tài)分布，1 kHz以下非正態(tài)分布。原因分析是高頻段風生噪聲占主要因素，低頻段船舶噪聲占主要因素；

（2）降雨主要影響噪聲頻段為500 Hz以上。暴雨時噪聲譜近于“白噪聲”，比晴朗時高出約25 dB。

在本文的研究基礎之上可以進一步把噪聲信號當作有用的信號進行寬頻帶和時間特性加以分析利用，從而實現(xiàn)海況和降雨級別等海洋中氣象參數(shù)的反演，這將對海洋開發(fā)、天氣預報、海洋氣候?qū)W研究及國防軍事起到重要的作用。

［1］徐功慧，陳鴻志，王二慶，等.海洋環(huán)境噪聲觀測技術(shù)及數(shù)據(jù)處理方法［J］.海洋技術(shù)，2011，30（2）：70-71.

［2］劉伯勝，雷家煜.水聲學原理［M］.哈爾濱：哈爾濱工程大學出版社，2006.

［3］Knudsen VO，Alford R S，EmlingJ W.Underwater Ambient Noise［J］.Journal ofMarine Research，1948，7:410-429.

［4］PiggottCL.AmbietSeaNoiseatLowFrequenciesin ShallowWateroftheScotian Shelfs［J］.Journalofthe AcousticalSocietyofAmerica，1964，36（11）:2152-2163.

［5］Crouch W W，Burt P J.The Logarithmic Dependence of Surface-Generated Ambient-Sea-Noise Spectrum Level on Wind Speed［J］.Journal ofthe Acoustical SocietyofAmerica（Part 2），1972，51（3）:1066-1072.

［6］Hsiang-Chih C，Chi-Fang C，Ruey-Chang W.Analysis of the Significantly Statistical Duration to Clarify the Uncertainty of Ambient Noises Due toEnvironmental Changes［C］//OCEANS2008 MTS/IEEE Kobe Techno-Ocean，2008：1-5.

Analysis on the Characteristics of Measured Ocean Ambient Noise Spectrum

WEI Yong-xing，YU Jin-hua，LI Qi，ZHOU Ying，CHANG Zhe，NIU Zhi-hua，GAO Chao
National Ocean Technology Center，Tianjin 300112，China

Through 94 d ocean ambient noise data and corresponding meteorological data collected from a submersible buoy system，joint time-frequency，correlation and statistical analyses are conducted to derive the relations between the noise spectrum and natural phenomena such as wind speed and rainfall.The analysis results show a lower correlation between noise spectrum and wind speed change in the lower frequency band and a higher correlation in the upper frequency bands.The distribution of noise spectrum level is close to the normal distribution above 1 kHz frequency and the effects of rainfall on ambient noise spectrum mainly appear in the band above 500 Hz.

ocean ambient noise；wind speed；correlation；statistical analysis；rainfall

TB566

A

1003-2029（2016）03-0036-04

10.3969/j.issn.1003-2029.2016.03.007

2015-07-11

國家海洋公益性行業(yè)科研專項資助項目（201005004）

魏永星（1987-），女，助理工程師，主要研究方向為數(shù)字信號處理與研究。E-mail：weiyongxing1987@163.com