劉啟軍,邱家興,程玉勝
(海軍潛艇學(xué)院,山東 青島 266199)
船舶螺旋槳空化噪聲非均勻調(diào)制特性及其應(yīng)用
劉啟軍,邱家興,程玉勝
(海軍潛艇學(xué)院,山東 青島 266199)
從螺旋槳空化噪聲譜統(tǒng)計(jì)模型出發(fā),分析研究螺旋槳空化噪聲非均勻調(diào)制特性的機(jī)理和規(guī)律,并利用大量海上實(shí)錄船舶噪聲進(jìn)行驗(yàn)證。從形成機(jī)理看,認(rèn)為非均勻調(diào)制特性主要由船尾形狀、螺旋槳結(jié)構(gòu)以及船舶航行工況共同決定的尾流分布所導(dǎo)致,不同型號(hào)船舶的非均勻調(diào)制特性具有個(gè)體性;從艦船噪聲信號(hào)能量譜看,調(diào)制的非均勻特性,主要受能量譜峰值頻率的影響。最后基于船舶噪聲的非均勻調(diào)制特性,提出一種基于多子帶自適應(yīng)加權(quán)的DEMON增強(qiáng)算法,并利用實(shí)際數(shù)據(jù)進(jìn)行驗(yàn)證。
螺旋槳空化噪聲;非均勻調(diào)制;調(diào)制度;多子帶;線譜增強(qiáng)
船舶螺旋槳一旦空化,其輻射噪聲就會(huì)成為船舶航行時(shí)最主要的噪聲,螺旋槳在船尾非均勻尾流中旋轉(zhuǎn)并產(chǎn)生空化的過程中,空化噪聲能量又會(huì)被周期性調(diào)制,表現(xiàn)為在空化調(diào)制譜中有明顯的具有倍數(shù)關(guān)系的多根線譜,利用這些調(diào)制線譜可以判斷出物理意義明確的螺旋槳轉(zhuǎn)速和槳葉數(shù)等特征,是目前船舶噪聲識(shí)別中非常重要的識(shí)別特征[1–3]。
螺旋槳空化噪聲在不同頻帶上的調(diào)制程度不均勻,但這種非均勻調(diào)制特性的規(guī)律如何,調(diào)制強(qiáng)弱隨頻帶分布關(guān)系,乃至這種非均勻調(diào)制特性的利用等問題,目前已公開發(fā)表的文獻(xiàn)中未見有詳細(xì)研究和結(jié)論。本文將從螺旋槳空化噪聲譜統(tǒng)計(jì)模型出發(fā),將螺旋槳空化噪聲非均勻調(diào)制特性的機(jī)理和規(guī)律作為研究重點(diǎn)并進(jìn)一步考慮其運(yùn)用。
螺旋槳在非均勻尾流中轉(zhuǎn)動(dòng),每個(gè)時(shí)刻都有大量隨機(jī)產(chǎn)生和崩潰的空化氣泡。對(duì)于螺旋槳空化噪聲譜,本文綜合采用文獻(xiàn)[4]中的群空泡崩潰統(tǒng)計(jì)模型和文獻(xiàn)[3]中假設(shè)空泡群平均最大半徑服從正態(tài)分布的特點(diǎn),進(jìn)行推導(dǎo)[5–7],得到螺旋槳空化噪聲能量譜的表達(dá)式為[4]:
式中:C為一比例常數(shù);t為時(shí)間;ω為頻率;KR為一取決于氣體含量和環(huán)境壓力等因素的常數(shù)因子,海況和目標(biāo)航行深度等的變化會(huì)引起其變化;N為某一時(shí)刻發(fā)生崩潰的空泡個(gè)數(shù);Rmax為氣泡崩潰開始時(shí)最大半徑;ρ為水的密度;P∞為無窮遠(yuǎn)處水的壓力(靜壓力);為同一時(shí)刻大量空泡崩潰最大半徑Rmax的均值;σ為方差。
為了表征能量譜幅度被調(diào)制的程度,定義mp為能量譜在一個(gè)時(shí)間周期T內(nèi),能量譜幅值變化的深度,即幅度調(diào)制深度。螺旋槳空化噪聲可近似視為周期性局部平穩(wěn)隨機(jī)過程,短時(shí)間內(nèi)(幾個(gè)螺旋槳轉(zhuǎn)動(dòng)周期內(nèi))螺旋槳轉(zhuǎn)速、海洋環(huán)境、尾流等相對(duì)穩(wěn)定,mp可視為隨時(shí)間穩(wěn)定不變,因此mp可僅看作是頻率ω的函數(shù),記為mp(ω);在任一時(shí)間周期T內(nèi),中心時(shí)刻為t,能量譜的時(shí)間平均值可視為相對(duì)穩(wěn)定,記為G-(ω,t),能量譜最大值和最小值分別出現(xiàn)在t1和t2時(shí)刻,則
根據(jù)信號(hào)調(diào)制理論,能量譜被周期性幅度調(diào)制,可表示為:
推得能量譜幅度調(diào)制深度為:
從式(8)可看出,能量譜幅度調(diào)制深度完全可以由能量譜最大值與最小值之比(最大相對(duì)起伏)來表示,用分貝數(shù)表示并記為:
通過計(jì)算,結(jié)合圖1容易得出:
1)4 個(gè)參數(shù)中,僅N(t)變化時(shí),mG(ω)在不同頻率ω處為常數(shù),即空泡個(gè)數(shù)的變化引起了能量譜在不同頻帶間的均勻調(diào)制;而其他 3 個(gè)參數(shù)(t),KR,σ的變化均會(huì)引起能量譜在不同頻帶間的非均勻調(diào)制;
2)引起能量譜的調(diào)制深度mG(ω)較大主要因素為N(t),(t),而KR,σ的變化引起的調(diào)制深度相對(duì)較小,尤其是在能量譜峰值頻率以上的中高頻段上更是如此;
從以上分析可知,螺旋槳空化在不同頻帶的非均勻調(diào)制特性主要是受空泡平均最大半徑的影響,而影響空泡平均最大半徑的因素是螺旋槳槳葉和海水之間的相對(duì)來流速度[4,9-10],因此,最根本的原因是由船尾形狀、螺旋槳結(jié)構(gòu)以及船舶航行工況共同決定的尾流分布所導(dǎo)致的[3]。由于不同型號(hào)船舶的船尾形狀、螺旋槳結(jié)構(gòu)具有個(gè)體性,因此其也具有特有的螺旋槳空化噪聲非均勻調(diào)制特性。
此外,據(jù)圖1的仿真分析還可得出,從螺旋槳空化輻射噪聲能量譜看,調(diào)制的非均勻特性,主要受能量譜峰值頻率的影響。在峰值頻率的左側(cè)低頻段,調(diào)制深度最高,在峰值頻率的右側(cè)一小段頻帶內(nèi),調(diào)制深度最低,在高頻段調(diào)制深度基本保持不變,略高于最低調(diào)制深度。分析這種現(xiàn)象的原因:一是螺旋槳空化噪聲能量譜本身受空泡半徑的影響最大,其次是當(dāng)空泡半徑變大時(shí),不僅空化噪聲譜線整體上升,而且峰值頻率還要向低頻移動(dòng),勢(shì)必造成在空化噪聲能量譜在低頻段調(diào)制深度高于高頻段,并在緊靠峰值頻率而又高于峰值頻率一側(cè)具有最低的調(diào)制深度[1]。
下面通過海上實(shí)錄的船舶噪聲信號(hào)進(jìn)行調(diào)制譜分析,驗(yàn)證以上理論是否符合實(shí)際情況。為了定量描述噪聲調(diào)制度大小,采用文獻(xiàn)[4]提出的噪聲調(diào)制度σs計(jì)算方法。雖然螺旋槳空化噪聲不是船舶輻射噪聲的全部,但一旦螺旋槳發(fā)生空化,就將成為船舶輻射噪聲的最主要部分,因此可以用海上實(shí)錄的高工況航行船舶噪聲近似研究螺旋槳空化噪聲的調(diào)制特性。對(duì)船舶噪聲進(jìn)行DEMON分析,通過分析計(jì)算軸頻及其諧波線譜、葉頻及其諧波線譜、干擾線譜之間的大小關(guān)系,判定噪聲的調(diào)制度大小[11]。在分析噪聲不同頻帶的非均勻調(diào)制特性時(shí),首先分頻段解調(diào),然后對(duì)各頻段解調(diào)獲得的DEMON譜分別計(jì)算調(diào)制度大小,最后組成調(diào)制度矢量,得到非均勻調(diào)制特性,計(jì)算過程如圖2所示。
從實(shí)錄船舶噪聲庫中隨機(jī)選取工況較高的目標(biāo)513個(gè),每個(gè)目標(biāo)噪聲長(zhǎng)度40 s~3 min不等,對(duì)每個(gè)目標(biāo)按照?qǐng)D2所示步驟分別計(jì)算得到一個(gè)非均勻調(diào)制特性矢量圖,從統(tǒng)計(jì)結(jié)果看,本文認(rèn)為船舶噪聲的非均勻調(diào)制特性雖然看似無章可循,但其實(shí)是有一定規(guī)律的,大多數(shù)噪聲基本滿足本文上節(jié)所述的“從螺旋槳空化輻射噪聲能量譜看…在峰值頻率的左側(cè)低頻段,調(diào)制深度最高,在峰值頻率的右側(cè)一小段頻帶內(nèi),調(diào)制深度最低”的特點(diǎn)。
圖3為某大型船舶噪聲非均勻調(diào)制特性圖。圖中:①為目標(biāo)功率譜圖縱坐標(biāo)為信號(hào)頻率0~10 kHz;②為目標(biāo)全頻帶噪聲DEMON譜圖,橫坐標(biāo)為調(diào)制頻率,范圍為0~50 Hz;③為目標(biāo)DEMON譜歷程圖,橫坐標(biāo)為調(diào)制頻率,縱坐標(biāo)為信號(hào)頻率;④為目標(biāo)非均勻調(diào)制度矢量圖,縱坐標(biāo)為信號(hào)頻率,橫坐標(biāo)為各頻率子帶的調(diào)制度大小,dB。
該信號(hào)的功率譜圖較好滿足標(biāo)準(zhǔn)船舶噪聲功率譜特點(diǎn):峰值頻率很低,在峰值頻率以下的頻段,由于傳感器濾波的原因,看不到幅值的衰減,在峰值頻率以上的中高頻段大約以–6 dB/倍頻程的速率下降;其非均勻性調(diào)制度矢量圖基本符合上節(jié)分析的結(jié)論特點(diǎn),因信號(hào)頻率基本位于峰值頻率以上,其非均勻調(diào)制特點(diǎn)則只具有峰值以上頻率的調(diào)制特性,在整個(gè)信號(hào)頻率范圍內(nèi),噪聲調(diào)制度逐漸遞增,在高頻部分遞增的幅度較小,而在靠近且低于1 kHz的信號(hào)頻段,由于恰巧處于峰值頻率的后沿,調(diào)制度最弱。
傳統(tǒng)DEMON譜分析,一般是依據(jù)噪聲接收的工作頻段,設(shè)置一個(gè)固定的較為通用的解調(diào)頻段,所有噪聲信號(hào)均在該頻段內(nèi)濾波、檢波、DEMON譜分析,并據(jù)此進(jìn)行線譜提取和特征分析。而根據(jù)本文第1節(jié)結(jié)論,船舶噪聲的調(diào)制特性在不同頻帶上是非均勻的,那么在DEMON譜分析中,如果僅僅在調(diào)制度高的信號(hào)頻段進(jìn)行解調(diào),獲得的DEMON譜線譜清晰度是否會(huì)更高呢,這是本節(jié)要研究的問題。
根據(jù)第1節(jié)的分析,對(duì)艦船輻射噪聲的解調(diào),實(shí)際上為對(duì)某個(gè)頻段ω1噪聲能量Gi(t)的頻譜分析,即DEMON為:
則全頻段的DEMON譜可以近似表示為N個(gè)子頻帶DEMON的譜值積累:
考慮到將調(diào)制度大小差異較大的子帶DEMON譜值積累一起,與僅將調(diào)制度高的子帶DEMON譜值積累一起相比,調(diào)制特性應(yīng)該會(huì)有所減弱,因此本文在自動(dòng)計(jì)算各子帶調(diào)制度的基礎(chǔ)上,僅將調(diào)制度高的子帶DEMON譜積累作為最終的DEMON譜,如圖4所示。根據(jù)大量數(shù)據(jù)統(tǒng)計(jì)結(jié)果給定一個(gè)調(diào)制度可用門限?,若某個(gè)子帶的調(diào)制度σs(ωi)大于門限?,則判定該子帶權(quán)值為1,否則為0,最終DEMON譜D(f)為:
表1 多子帶自適應(yīng)加權(quán)的DEMON譜線譜增強(qiáng)效果Tab.1 The effect of DEMON line spectrum strength based onmulti sub-band adaptive weighted
仍然采用第2節(jié)使用的513個(gè)海上實(shí)錄船舶噪聲數(shù)據(jù)進(jìn)行分析計(jì)算,分別進(jìn)行全頻段解調(diào)和采用不同加權(quán)門限?時(shí)的多子帶自適應(yīng)加權(quán)解調(diào)。統(tǒng)計(jì)513個(gè)目標(biāo)的平均調(diào)制度和線譜平均增強(qiáng)幅度,如表1所示。2個(gè)典型目標(biāo)噪聲的DEMON譜線譜增強(qiáng)前后對(duì)比如圖5所示。
從實(shí)際噪聲數(shù)據(jù)的分析結(jié)果看,頻帶的非均勻調(diào)制特性對(duì)DEMON譜線譜幅值的計(jì)算結(jié)果影響還是比較明顯的,尤其是當(dāng)噪聲的調(diào)制特性僅分布在比較窄的頻帶內(nèi),影響更甚。通過多子帶自適應(yīng)加權(quán),可以獲得較優(yōu)的DEMON譜線譜清晰度。另外,權(quán)限?取較低值(5 dB)時(shí)就已經(jīng)取得的較好的增強(qiáng)效果,而繼續(xù)把權(quán)限值?升高,增強(qiáng)的幅度減緩,這是因?yàn)闄?quán)限?給定一個(gè)較低值時(shí),已經(jīng)將沒有調(diào)制的頻帶濾掉,基本上消除了無調(diào)制頻帶對(duì)調(diào)制特性的影響。
本文主要研究了螺旋槳空化非均勻調(diào)制特性的機(jī)理、規(guī)律和應(yīng)用。首先從船舶螺旋槳空化噪聲譜模型出發(fā),仿真分析了螺旋槳空化非均勻調(diào)制特性;然后利用大量海上實(shí)錄船舶噪聲,統(tǒng)計(jì)了實(shí)際船舶噪聲非均勻調(diào)制特性,驗(yàn)證了以上結(jié)論的合理性;最后基于船舶噪聲的非均勻調(diào)制特性,提出了一種基于多子帶自適應(yīng)加權(quán)的DEMON增強(qiáng)算法,試驗(yàn)表明其效果良好。主要形成以下結(jié)論:
1)從形成機(jī)理看,螺旋槳空化在不同頻帶的非均勻調(diào)制特性主要是由船尾形狀、螺旋槳結(jié)構(gòu)以及船舶航行工況共同決定的尾流分布所導(dǎo)致的,不同型號(hào)船舶的非均勻調(diào)制特性具有個(gè)體性。
2)從艦船噪聲信號(hào)能量譜看,調(diào)制的非均勻特性,主要受能量譜峰值頻率的影響。在峰值頻率的前沿低頻段,調(diào)制深度最高,在峰值頻率的后沿一小段頻帶內(nèi),調(diào)制深度最低,在高頻段調(diào)制深度基本保持不變,略高于最低調(diào)制深度。
3)船舶噪聲的非均勻調(diào)制特性對(duì)DEMON譜線譜幅值影響明顯,尤其是當(dāng)噪聲的調(diào)制特性僅分布在比較窄的頻帶內(nèi)時(shí),影響更甚。通過多子帶自適應(yīng)加權(quán),對(duì)DEMON譜線譜增強(qiáng)有明顯效果,可以獲得更優(yōu)的線譜清晰度與譜圖結(jié)構(gòu)特征。
[ 1 ]史廣智. 艦船噪聲調(diào)制譜諧波族結(jié)構(gòu)特性理論分析[J]. 聲學(xué)學(xué)報(bào), 2007, 32(1): 19–25.
SHI Guang-zhi. Theoretical analysis of the structure law of ship radiated-noise demodulation spectrum harmonic clan feature[J]. Acta Acustica, 2007, 32(1): 19–25.
[ 2 ]史廣智. 基于模板匹配的艦船螺旋槳噪聲特征精細(xì)分析[J].聲學(xué)學(xué)報(bào), 2009, 34(5): 401–407.
SHI Guang-zhi. A method for ship propeller blade-number recognition based on template matching[J]. Acta Acustica, 2009, 34(5): 401–407.
[ 3 ]陶篤純. 螺旋槳空化噪聲譜[J]. 聲學(xué)學(xué)報(bào), 1982, 7(6): 344–351.
TAO Du-chun. Spectrum of propeller cavitation noise[J]. Acta Acustica, 1982, 7(6): 344–351.
[ 4 ]邱家興, 程玉勝, 等. 船舶大側(cè)斜螺旋槳空化噪聲調(diào)制譜諧波結(jié)構(gòu)[J]. 聲學(xué)技術(shù)增刊-中國聲學(xué)學(xué)會(huì)水聲學(xué)分會(huì)2015年學(xué)術(shù)會(huì)議論文集, 2015, 34(2): 364–368.
QIU Jia-xing, CHENG Yu-sheng, et al. Demodulation spectrum harmonic line structure of high-skewed propeller cavitation noise[J]. Technical Acoustics, 2015, 34(2): 364–368.
[ 5 ]戚定滿, 魯傳敬, 何友聲. 空化噪聲譜特性研究[J]. 振動(dòng)與沖擊, 1999, 18(3): 32–36.
QI Ding-man, LU Chuang-jing, HE You-sheng. An investigation on the spectra of bubble noise[J]. Journal of Vibration and Shock, 1999, 18(3): 32–36.
[ 6 ]ROSS D. Mechanica of Underwater Noise[M]. 北京: 《水下噪聲原理翻譯組》譯. 1983: 277–394.
[ 7 ]陶篤純. 艦船噪聲節(jié)奏的研究(I)—數(shù)學(xué)模型及功率譜密度[J]. 聲學(xué)學(xué)報(bào), 1983, 8(2): 65–76.
TAO Du-chun. A study on ship radiated noise rhythms(i)-mathematical model and power spectrum densidy[J]. Acta Acustica, 1983, 8(2): 65–76.
[ 8 ]陶篤純. 艦船噪聲節(jié)奏的研究(II)—自相關(guān)函數(shù)及節(jié)奏信息的提取[J]. 聲學(xué)學(xué)報(bào), 1983, 8(5): 280–289.
TAO Du-chun. A study on ship radiated noise rhythms(ii)-autocorrelation function and extraction of rhythm information[J]. Acta Acustica, 1983, 8(5): 280–289.
[ 9 ]魏以邁. 均流與非均流中螺旋槳空泡噪聲的研究[J]. 聲學(xué)技術(shù), 1983, (4): 39–50.
WEI Yi-mai. Study on propeller cavitation noise under the steady and unsteady hydrodynamic[J]. Ship Science and Technology, 1983, (4): 39–50.
[10]何友聲, 王國強(qiáng). 螺旋槳激振力[M]. 上海: 上海交通大學(xué)出版社, 1987: 13–67.
HE You-sheng, Wang Guo-qiang. Propeller exciting force[M]. Shanghai: Shanghai Jiao Tong University Press, 1987: 13–67.
[11]邱家興, 程玉勝, 等. 船舶噪聲DEMON譜質(zhì)量評(píng)估方法研究[J]. 艦船科學(xué)技術(shù), 2014, 36(9): 46–49.
QIU Jia-xing, CHENG Yu-sheng, et al. Study of ship-radiated noise demon quality assessment methods[J]. Ship Science and Technology, 2014, 36(9): 46–49.
Non-uniform modulation feature of ship propeller cavitation noiseand its application
LIU Qi-jun, QIU Jia-xing, CHENG Yu-sheng
(Navy Submarine Academy, Qingdao 266199, China)
Based on the statistical model of ship propeller cavitation noise spectrum, the mechanism and laws of nonuniform modulation feature of propeller cavitation noise are analyzed and studied, and also be verified by a large number of ship radiated noise. From thepoint of formed mechanism, it shows that non-uniform modulation feature is mainly formed by the distribution of wake flow, which is formed by the combined action of stern shape, propeller structure and working condition. Non-uniform modulation feature vary in different class ship. Then from the point of ship radiated noise spectrum, nonuniform modulation feature is mainly decided by the energy spectrum peak frequency. Finally, an application method of nonuniformity envelope is proposed, that is the technology of DEMON line spectrum strength based on multi sub-band adaptive weighted, and it is also verified by real date.
propeller cavitation noise;non-uniform modulation;modulation degree;multi sub-band;line spectrum strength
TB533
A
1672 – 7619(2017)06 – 0018 – 05
10.3404/j.issn.1672 – 7619.2017.06.004
2017 – 03 – 14
劉啟軍(1975 – ),男,碩士,講師,研究方向?yàn)樗暷繕?biāo)識(shí)別。