丁超，邱家興，程健，程玉勝

船舶螺旋槳空化噪聲解調(diào)包絡(luò)形狀仿真研究

丁超，邱家興，程健，程玉勝

(海軍潛艇學(xué)院，山東青島 266000)

船舶螺旋槳空化噪聲解調(diào)包絡(luò)是一系列形狀基本一致、間隔基本相等、高度有一定隨機性、成組結(jié)構(gòu)的波包。而波包形狀會直接影響到解調(diào)譜的結(jié)構(gòu)形狀。通過分析比較實測信號、解調(diào)包絡(luò)形狀，將波包形狀分為兩大類：圓滑型和尖銳型，并分別設(shè)置了對應(yīng)的兩種仿真模型：半正弦模型和高斯模型。仿真信號和實測信號的檢驗結(jié)果表明，半正弦模型能夠更好地模擬圓滑型包絡(luò)的特點，高斯模型能夠更好地模擬尖銳型包絡(luò)的特點。

解調(diào)包絡(luò)；波包形狀；包絡(luò)模型

0 引言

當(dāng)船舶螺旋槳發(fā)生空化時，船舶輻射噪聲的主要成分就成為螺旋槳空化噪聲。在時域上，螺旋槳空化噪聲呈現(xiàn)出明顯的振幅調(diào)制。通過檢波處理，可以從時域中提取出振幅的起伏趨勢，即解調(diào)包絡(luò)。解調(diào)包絡(luò)受到船尾形狀、槳葉形狀、槳葉間差異等因素的影響會呈現(xiàn)出不同的形狀[1]，從而對解調(diào)譜的結(jié)構(gòu)形狀產(chǎn)生很大的影響。

本文觀察分析了大量實測船舶輻射噪聲的解調(diào)包絡(luò)，將波包形狀劃分為圓滑型和尖銳型兩類，并給出了與之對應(yīng)的半正弦模型和高斯模型兩種仿真模型。最后，通過仿真信號和實測信號驗證了兩種包絡(luò)模型的適用性。

1 解調(diào)方法對解調(diào)包絡(luò)形狀的影響

從螺旋槳空化噪聲中提取調(diào)制包絡(luò)，需要經(jīng)過檢波，流程如圖1所示。

圖1 檢波流程圖

該流程中，帶通濾波用來選擇合適的檢波頻段，低通濾波用來提取低頻包絡(luò)信號。常用的檢波方法有平方檢波、絕對值檢波和希爾波特檢波。絕對值檢波和希爾波特檢波效果基本一致，因此本文僅比較平方檢波與絕對值檢波的處理效果。此外，考慮到帶通濾波頻段可能對解調(diào)包絡(luò)形狀造成影響，本文在進行實測信號解調(diào)處理時，均采用1～5 kHz的帶通濾波頻段。

1.1 空化噪聲仿真

圖2 調(diào)制信號

獲得的仿真信號如圖3所示。

圖3 仿真信號

分別用平方解調(diào)方法和絕對值解調(diào)方法對圖3中的信號進行解調(diào)，得到的結(jié)果如圖4(a)、4(b)所示。

圖4 不同檢波方法處理效果比較

對比圖2與圖4可以發(fā)現(xiàn)，絕對值解調(diào)包絡(luò)與原始調(diào)制信號更為接近；平方解調(diào)方法不但使解調(diào)包絡(luò)形變更大，而且同組4個波包的幅值也由1:2:3:4變?yōu)榧s1:4:9:16。

1.2 蒙特卡洛實驗

為比較兩種檢波方法對解調(diào)包絡(luò)造成的影響，進行蒙特卡洛實驗。

首先定義平均歐氏距離：

分別計算平方解調(diào)和絕對值解調(diào)方法的平均歐氏距離。進行10 000次蒙特卡洛實驗，求10 000次平均歐氏距離的均值，結(jié)果如表1所示。

表1 蒙特卡洛實驗結(jié)果

由表1可知，絕對值解調(diào)方法得到的包絡(luò)形變更小，因此，本文在后續(xù)的信號解調(diào)過程中均采用絕對值解調(diào)方法。

2 解調(diào)包絡(luò)形狀分類

2.1 波包形狀分類

對船舶輻射噪聲進行絕對值解調(diào)，可以獲得該信號的解調(diào)包絡(luò)。本文對大量的商船輻射噪聲進行了絕對值解調(diào)處理，并根據(jù)解調(diào)包絡(luò)質(zhì)量進行了信號篩選，篩選標(biāo)準(zhǔn)如下：

(1) 解調(diào)包絡(luò)清晰可見、周期明顯；

(2) 允許存在毛刺，但毛刺不能掩蓋包絡(luò)的基本形狀。

通過比較篩選之后的船舶輻射噪聲解調(diào)包絡(luò)可以發(fā)現(xiàn)，包絡(luò)形狀大致可以分為兩類：一類如圖5(a)所示，包絡(luò)總體形狀趨勢起伏較平緩，稱其為“圓滑型包絡(luò)”；另一類如圖5(b)所示，包絡(luò)總體形狀趨勢起伏較劇烈，稱其為“尖銳型包絡(luò)”。

圖5 不同類型的螺旋槳空化噪聲解調(diào)包絡(luò)

篩選之后的船舶輻射噪聲共42組，其解調(diào)包絡(luò)形狀統(tǒng)計結(jié)果如表2所示。

表2 不同類型解調(diào)包絡(luò)的占比

基于統(tǒng)計結(jié)果，在信噪比較高、解調(diào)包絡(luò)形狀較規(guī)則的信號中，解調(diào)包絡(luò)是圓滑型的噪聲信號占比76.2%。

2.2 波包形狀模型

文獻[1]認(rèn)為，船舶輻射噪聲解調(diào)包絡(luò)形狀符合高斯型，如圖6(a)所示，即：

通過觀察實測船舶輻射噪聲解調(diào)包絡(luò)可以發(fā)現(xiàn)，高斯模型與尖銳型包絡(luò)形狀較為吻合；大多數(shù)情況下，包絡(luò)波形為圓滑型，形狀類似于半個正弦波，如圖6(b)所示，即：

其中，是幅值相關(guān)系數(shù)，f是信號頻率。

3 調(diào)制包絡(luò)功率譜

調(diào)制包絡(luò)對載波信號進行調(diào)制，通過信號解調(diào)可以獲得解調(diào)包絡(luò)。經(jīng)過解調(diào)之后，解調(diào)包絡(luò)與原始的調(diào)制包絡(luò)相比存在一定的失真，但信噪比較高時仍能保持較高的相似性。因此，對形狀規(guī)則的原始調(diào)制包絡(luò)進行理論推導(dǎo)，其結(jié)果仍然適用于解調(diào)包絡(luò)。

3.1 理論模型

為簡化表示，令：

則：

令：

則式(8)可以簡化為

根據(jù)隨機過程功率譜密度的基本定義，該隨機過程的功率譜密度為

將式(12)代入式(11)可得：

由式(13)可知，脈沖形狀將會影響到調(diào)制包絡(luò)功率譜的整體形狀。只有使用更符合真實情況的脈沖模型，才能得到更好的仿真結(jié)果。

3.2 理論模型仿真

如圖7所示，高斯模型包絡(luò)得到的功率譜軸頻高頻諧波強度較高；半正弦模型包絡(luò)得到的功率譜軸頻高階諧波強度微弱。

通過對大量實測船舶輻射噪聲的分析觀察可以發(fā)現(xiàn)：具有尖銳型包絡(luò)的噪聲信號，解調(diào)譜出現(xiàn)軸頻高階諧波的概率較高；具有圓滑型包絡(luò)的噪聲信號，解調(diào)譜出現(xiàn)軸頻高階諧波的概率較低。證明高斯模型更適用于尖銳型包絡(luò)，半正弦模型更適用于圓滑型包絡(luò)。

圖7 不同包絡(luò)模型調(diào)制信號的功率譜對比

4 仿真信號驗證

4.1 仿真模型建立

關(guān)于螺旋槳空化噪聲已經(jīng)有許多研究，目前有兩類典型的空化模型[2]：

第一類模型假定低頻調(diào)制信號對高頻載波進行調(diào)制，得到最終的空化噪聲為：

其中，為低頻調(diào)制信號，為高頻載波。

第二類模型將螺旋槳空化所形成的空泡看作是一個空泡集，每一葉片旋轉(zhuǎn)一次都是對這個空泡集的一次擾動，使得這個空泡集有一次輻射聲脈沖[3]。文獻[4]-[6]中所使用的方法也屬這一模型。

第一類模型中，通常將低頻調(diào)制信號假設(shè)為幾個成倍頻關(guān)系的正弦波，這種假設(shè)便于公式的推導(dǎo)并能獲得較為理想的解調(diào)譜，但其解調(diào)包絡(luò)與實際情況不符。本文將兩種模型結(jié)合，將低頻調(diào)制信號假設(shè)為圖2形式的成組結(jié)構(gòu)的脈沖，單個波包分別設(shè)置為高斯型和半正弦型。

4.2 仿真結(jié)果對比

圖8所示的高斯型仿真信號的解調(diào)譜中，軸頻高階諧波較為明顯；圖9所示的半正弦型仿真信號的解調(diào)譜中，軸頻高階諧波較為微弱。該結(jié)果與圖7結(jié)果基本一致，也與實際情況基本一致。

圖8 高斯型包絡(luò)調(diào)制信號的仿真

可以證明，高斯脈沖模型更符合具有尖銳型解調(diào)包絡(luò)的噪聲信號的特征規(guī)律；半正弦脈沖模型更符合具有圓滑型解調(diào)包絡(luò)的噪聲信號的特征規(guī)律。

5 實測信號驗證

5.1 信號處理流程

在船舶輻射噪聲解調(diào)包絡(luò)與兩個模型的比對過程中，為保證比對結(jié)果有效，需要信噪比足夠高，解調(diào)包絡(luò)毛刺較少。基于以上考慮，從前文篩選出的42組商船輻射噪聲信號中選取解調(diào)包絡(luò)形狀規(guī)則、毛刺較少的圓滑型包絡(luò)信號和尖銳型包絡(luò)信號各10個進行處理，解調(diào)包絡(luò)與包絡(luò)模型比對流程如圖10所示。

圖10 解調(diào)包絡(luò)與設(shè)置的包絡(luò)模型比對流程

其中，選取單個波包是選取信號包絡(luò)中形狀規(guī)則、干擾較少的信號。

5.2 處理結(jié)果分析

對每個信號各取10個形狀較規(guī)則的包絡(luò)，分別與高斯模型和半正弦模型對比，求出10個平均歐氏距離，并取均值，計算結(jié)果如圖11所示。

從圖11中可以發(fā)現(xiàn)，圓滑型包絡(luò)信號與半正弦模型的平均歐氏距離更小，尖銳型包絡(luò)信號與高斯模型的平均歐氏距離更小，說明半正弦模型更符合圓滑型包絡(luò)，高斯模型更符合尖銳型包絡(luò)。并且由表2可知，在統(tǒng)計中76.2%的信號具有圓滑型包絡(luò)。因此，半正弦模型更具有代表性。

圖11 兩種不同包絡(luò)調(diào)制信號的平均歐氏距離分布圖

6 結(jié)論

(1) 通過蒙特卡洛實驗發(fā)現(xiàn)，絕對值解調(diào)方法對調(diào)制信號的形變較?。?/p>

(2) 船舶輻射噪聲解調(diào)包絡(luò)可以分為圓滑型和尖銳型兩種，通過統(tǒng)計發(fā)現(xiàn)，大多數(shù)信號的解調(diào)包絡(luò)是圓滑型；

(3) 通過仿真包絡(luò)信號、實測信號包絡(luò)與半正弦模型、高斯模型平均歐氏距離分析發(fā)現(xiàn)，圓滑型包絡(luò)與半正弦模型更為吻合，尖銳型包絡(luò)與高斯模型更為吻合，且半正弦模型具有更廣泛的適用性。

[1] 陶篤純. 艦船噪聲節(jié)奏的研究(I)—數(shù)學(xué)模型及功率譜密度[J]. 聲學(xué)學(xué)報, 1983, 8(2): 65-69.

TAO Duchun. A study on ship radiated noise rhythms (I)– mathematical model and power spectrum density[J]. Acta Acust, 1983, 8(2): 65-69.

[2] 車永剛, 程玉勝. 螺旋槳空化噪聲譜模型及結(jié)構(gòu)特征[J]. 哈爾濱工程大學(xué)學(xué)報, 2010, 31(7): 837-841.

CHE Yonggang, CHENG Yusheng. Interpreting features of a propeller’s cavitation spectrum[J]. Journal of Harbin Engineering University, 2010, 31(7): 837-841.

[3] D. 羅斯. 水下噪聲原理[M]. 關(guān)定華, 譯. 北京: 海洋出版社, 1983: 371-373.

[4] 蔣國健, 林建恒, 馬杰, 等. 艦船螺旋槳空泡噪聲的數(shù)理模型[J]. 聲學(xué)學(xué)報, 1998, 23(5): 401-403.

JIANG Guojian, LIN Jianheng, MA Jie, et al. Theoretical model of noise caused by ship propeller cavitation[J]. Acta Acustica, 1998, 23(5): 401-408.

[5] 湯渭霖. 水下噪聲學(xué)[M]. 哈爾濱: 哈爾濱船舶工程學(xué)院, 1984: 100-113.

TANG Weilin. Underwater noise[M]. Harbin: Harbin Shipbuilding Engineering Institute, 1984: 100-113.

[6] 陶篤純. 螺旋槳空化噪聲譜[J]. 聲學(xué)學(xué)報, 1982, 7(6): 344-351.

TAO Duchun. Spectrum of propeller cavitation noise[J]. Acta Acustica, 1982, 7(6): 344-351.

Simulation research on the demodulation envelope shapes of vessel propeller cavitation noise

DING Chao, QIU Jia-xing, CHENG Jian, CHENG Yu-sheng

(Navy Submarine Academy, Qingdao 266000, Shandong, China)

The demodulation envelopes of vessel propeller cavitation noise is a series of wave packets with basically similar shapes, almost equal intervals and a certain randomness in amplitude, moreover, these wave packets are in groups. The shapes of wave packets can directly affect the structures and shapes of the demodulation spectrums. By analyzing and comparing the demodulation envelope shapes of measured signals, it is shown that the shapes of wave packets can be divided into two types: smooth type and sharp type, and correspondingly, two simulation models can be set up: semi-sine wave model and Gaussian model. The test results of simulated and measured signals show that the semi-sine wave model can better simulate the characteristics of smooth envelope, and the Gaussian model can better simulate the characteristics of sharp envelope.

demodulation envelope; the shape of wave packet; envelope model

TN911.7

A

1000-3630(2019)-06-0611-06

10.16300/j.cnki.1000-3630.2019.06.002

2018-09-18;

2018-10-20

丁超(1990－), 男, 山東青島人, 博士研究生, 研究方向為水聲目標(biāo)識別。

丁超, E-mail: 15964918058@163.com