徐功慧,陳鴻志,王二慶,韓引海
(91388部隊95分隊,廣東 湛江 524022)
海洋環境噪聲觀測技術及數據處理方法
徐功慧,陳鴻志,王二慶,韓引海
(91388部隊95分隊,廣東 湛江 524022)
海洋環境噪聲是水聲信道中的干擾背景場,是影響主動聲納探測能力、限制被動聲納工作性能發揮的主要原因之一。文中詳細介紹了目前海洋環境噪聲觀測的技術手段、測量方法和數據處理方法,并根據實際測量的海洋環境噪聲數據,詳細分析了海洋環境噪聲1/3倍頻程帶寬內的噪聲頻譜級計算方法和步驟。
海洋環境噪聲;數據處理;噪聲功率譜級
實際海洋環境中,海底地震和火山爆發、海面風浪、海洋生物活動、海上航運和石油勘探等自然和人為活動產生的聲波,在傳播過程中與海面、海底、水體等發生相互作用形成一個復雜的背景噪聲場,這些背景噪聲就是通常所說的海洋環境噪聲。海洋環境噪聲是水聲信道中的干擾背景場,會對聲納裝備的探測形成干擾,是限制被動聲納工作性能發揮的主要因素之一。要提高聲納抗干擾的能力不僅需要掌握環境的噪聲譜,還要了解噪聲場的時空統計特性和差異,以及對海況的依賴關系,為此需要測量在不同季節(時間)和不同海區、不同海況條件下海洋環境噪聲頻帶聲壓級和噪聲譜級,深刻認識海洋環境噪聲的時域、頻域和空域等特性,才能正確評價并提高聲納裝備的戰術性能。
海洋環境噪聲測量是對指定的海洋環境在規定的頻率內進行噪聲測量。基本的方法是使用抗流感噪測試系統來測量海洋環境噪聲。目前采用三類測量系統:第一類是基于測量船的測量方法,與水聲傳播實驗相結合的方法;第二類是采用浮標/潛標的測量方法,用錨系浮標/潛標測量海洋環境噪聲具有公認的優越性,由于穩定的聲接收基陣以及無工業噪聲,為在很寬頻段范圍(包括很低的次聲頻段)測量海洋環境噪聲特性提供了最佳條件;第三類是利用海軍岸站聲納或者警戒系統(包括水下監測網絡)進行海洋環境噪聲特性的測量,也將是我國噪聲監測技術發展的一個主要趨勢,具有測量長期性和能夠直接為聲納系統服務的特點。

圖1 海洋環境噪聲采集和處理的基本流程
海洋環境噪聲測量所使用的水聽器,需要有較高的靈敏度。對于帶前置放大器的水聽器,要求其前置放大器的本底噪聲級別要比較低。單個水聽器陣元的技術要求如下:
(1) 自由場靈敏度不低于-185 dB,參考值 1 V/μPa;
(2)水聽器等效噪聲和前置放大器以及其它電路總噪聲量級低于Wenz噪聲譜曲線的低限;
(3)在測量頻率范圍內,各個水聽器自由場靈敏度的不均勻性在±2 dB之內;
(4)在頻率為50~20 kHz區間內,其水平指向性圖與理想無指向性圖的偏差在±1.5 dB之內;
(5)在頻率為50~20 kHz區間內,垂直指向性的-3 dB波束寬度大于60°。
為了獲得精確的噪聲測量數據,在進行A/D變換時,應采用12位以上的A/D變換器,采樣率要高于51.2 kHz。在采樣之前,應進行寬帶抗混疊濾波器處理。
獼猴桃鮮果→清洗稱量→破碎打漿→果膠酶酶解→調制→發酵→固液分離→二次發酵→低溫靜置→澄清過濾→調配→滅菌→成品[19-20]
海洋環境噪聲信號是一種寬帶的隨機信號,在短時間內(一般是數分鐘內)它可視為平穩過程,如圖2所示,海洋環境噪聲原始信號。海洋環境噪聲級的定義有3種,包括寬帶總聲級、窄帶聲級和譜級。寬帶總聲級是指從感興趣的低截止頻率到高頻截止頻率的這個寬頻段內的噪聲總功率,如圖3所示;環境噪聲窄帶聲級是指某一窄帶內的環境噪聲總功率;環境噪聲譜級是指在1 Hz帶寬內環境噪聲功率按頻率的分布。寬帶海洋環境噪聲級可以提供信號所在海域的噪聲總體量級,使研究者對該海域的噪聲量級建立一個初步的概念,而窄帶噪聲級可以提供聲納所感興趣頻段內噪聲的量級,為聲納的信號處理技術提供依據;環境噪聲譜級可以提供每個頻率處環境噪聲的具體量級,使研究者能對環境噪聲量級隨頻率的變化有詳細的了解,為聲納設計和使用者提供一個有效的“聲納作戰地圖”,有利于我方設計和選擇聲納使用頻段,進一步地隱蔽我方目標信號,提取敵方目標。
在測量系統回收后,要進行適當的數據處理,以獲得海洋環境噪聲譜級、指向性和空間相關特性。將采集到的噪聲數據以每小時為單位分段,分段后使用FFT方法作功率譜分析。譜分析時,應把數據再細分為多段,然后對多段數據分別作FFT,并保證FFT的頻率分辨率小于1 Hz。對每段數據應進行加窗處理。由于數據量較多,在數據分段時可以不用重疊。多段數據之間進行功率譜平均,獲得該時間段的噪聲譜級。具體步驟為:

圖2 海洋環境噪聲原始信號圖

圖3 海洋環境噪聲寬帶總聲級
⑴1/3倍頻程平均功率譜計算
設有效噪聲信號為x(n),長度為L,將其按覆蓋分段法分為I段,每段長度為N,N根據分辨率、方差要求選定。后按平滑平均周期圖法求出其平滑平均周期圖功率譜估計P(k)。設第 i段信號序列為 xi(n),則:將 xi(n)乘以窗函數 w(n),再求 DFT,即:

式中:i=1,2,……,I

則第i段修正(加窗)周期圖為:


長度為L的整個實平穩隨機序列 的功率譜估值Pn(k)是上述I個修正周期圖的指數平均,即

式中:1≤n≤I

對線性平均:

式中,k=0,1,……,N-l
將Pn(k)或P(k)校正到工程上使用 S(k),如100 mV正弦信號以 1 μV 為參考級,S(k)應為 100 dB,此時 Pn(k)或 P(k)與S(k)有一差值,應扣除。這個差值與不同窗及樣點數有關。
⑵求1/3倍頻程帶寬總功率

式中:fh,fl分別為1/3頻程中心頻率fi的上下限頻率對應的頻率位置;△f為頻率分辨率,等于fs/N。fs為采樣率。
⑶求帶寬歸一化的噪聲功率譜級SPL(fi)

⑷求噪聲譜源級

式中:M(fi)為水聽器靈敏度;K(fi)為系統頻響;r為水聽器到等效聲中心的距離;20lgr為距離r(m)處反平方規律衰減,必要時按本地水聲傳播特性進行修正。

圖4 海洋環境噪聲1/3倍頻程平均功率譜級
⑸根據PSD(fi)求所測頻帶內的總聲源級Lp

式中:n為帶寬內包含的1/3倍頻程平均功率譜數。
在最終給出噪聲譜級時,應注意要折算到1 Hz帶寬內。為了方便和國外著名的Wenz噪聲曲線比較,通過積分求出1/3倍頻程帶寬內的噪聲頻帶級,然后通過平均給出每個1/3倍頻程帶寬內的噪聲頻譜級,如圖4所示。
海洋環境噪聲作為一種背景場,具有相對小的強度級別,對實驗技術和測量設備性能要求比較高。環境噪聲同海區具體環境、航船分布、氣象、季節變化等因素密切相關,需要海洋水文氣象條件的同步觀測。從長期測量和記錄角度來看,可以考慮采用沉底布放的潛標接收形式,布放到感興趣、敏感的海域,進行長時間監測記錄,這樣的數據分析才更具有統計意義。
[1]閻福旺,凌青,等.海洋水聲實驗技術[M].北京:海洋出版社,1999.
[2]R J尤立 克.水聲原理[M].哈爾濱:哈爾濱船舶工程學院出版社,1990.
[3] 汪德昭,尚爾昌.水聲學[M].北京:科學出版社,1981.
Ocean Ambient Noise Observation Technology and Signal Processing Methods
XU Gong-hui,CHEN Hong-zhi,WANG Er-qing,HAN Yin-hai
(Unit 91388,Zhanjiang Guangdong 524001,China)
Marine environmental noise is the interference background field of acoustic channels,it’s also the main factor affecting the detecting ability of active sonar and limiting the working performance of passive sonar.The observation technological means,measuring methods and the signal processing methods for ocean ambient noise are introduced.Based on observation data,the calculation methods and procedures of noise spectrum level on 1/3 octave bandwidth are detailedly analyzed.
ocean ambient noise;signal processing;noise power spectrum level
P716+.41,P733.22
A
1003-2029(2011)02-0069-03
2010-06-23
軍內科研項目(HJ2009-04-05-01)
徐功慧(1968-),男,山東日照人,工程師,主要從事海洋環境與目標特性測量研究