何 琪，宮在曉，李整林，李風(fēng)華，郭良浩

（1.中國(guó)科學(xué)院聲學(xué)研究所聲場(chǎng)聲信息國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京 100190；2.中國(guó)科學(xué)院大學(xué)，北京 100049；3.中山大學(xué)海洋工程與技術(shù)學(xué)院，廣東珠海 519000）

0 引言

陣形校準(zhǔn)精度對(duì)目標(biāo)測(cè)向性能有著顯著影響[1]。布放于海底的柔性水平陣，由于受到布放、內(nèi)波、浪涌等因素的影響，實(shí)際陣形往往與預(yù)設(shè)陣形有較大差異。進(jìn)行陣列形狀校準(zhǔn)工作，可以提高后續(xù)信號(hào)處理算法的探測(cè)性能，尤其是各種自適應(yīng)算法的性能。

已有的陣形校準(zhǔn)方法大致分為兩類(lèi)：一類(lèi)為無(wú)源校準(zhǔn)[2-3]，無(wú)需額外的校準(zhǔn)聲源，利用噪聲源就可以實(shí)現(xiàn)，但是反演得到的是陣列形狀，只知道兩兩陣元之間的距離，需要額外信息對(duì)初步結(jié)果進(jìn)行平移、旋轉(zhuǎn)、翻轉(zhuǎn)等操作才能獲取各陣元的準(zhǔn)確位置；另一類(lèi)為有源校準(zhǔn)[4]，通過(guò)利用校準(zhǔn)聲源，結(jié)合聲源與參考陣元的位置信息，測(cè)量各相鄰陣元間接收信號(hào)的時(shí)間延遲，從而獲得各陣元相對(duì)參考陣元的準(zhǔn)確位置；但由于淺海信道存在縱向相關(guān)振蕩現(xiàn)象[5]，以及受到海洋環(huán)境噪聲等因素影響，傳統(tǒng)方法[4]會(huì)因?yàn)闀r(shí)延估計(jì)誤差累積造成校準(zhǔn)陣形的偏差。

本文提出基于時(shí)延冗余測(cè)量的陣形校準(zhǔn)方法，通過(guò)估計(jì)全部?jī)蓛山M合的陣元間的信號(hào)時(shí)延，充分利用數(shù)據(jù)，可以有效提高校準(zhǔn)精度，并通過(guò)數(shù)值仿真與實(shí)驗(yàn)對(duì)方法進(jìn)行了驗(yàn)證。

1 理論分析

1.1 信號(hào)模型

根據(jù)簡(jiǎn)正波理論[6]，在水平不變分層的海洋波導(dǎo)環(huán)境中，點(diǎn)聲源激發(fā)的單頻聲場(chǎng)可以表示為M階簡(jiǎn)正波的線性疊加

其中：krm為第m階模態(tài)的水平波數(shù),Ψm為第m階模態(tài)的深度函數(shù),ρ(zs)為聲源處水介質(zhì)的密度,zs和z分別為聲源深度和接收深度，r為陣元與聲源間的水平距離。

假定在頻帶[fL,fH]內(nèi)，聲源信號(hào)的頻譜為S(f)，則第i號(hào)陣元接收該頻段信號(hào)的時(shí)域波形為

其中：ri為第i號(hào)陣元與聲源間的水平距離。

1.2 陣形校準(zhǔn)方法

校準(zhǔn)方法示意圖如圖1 所示，S1、S2分別為校準(zhǔn)聲源的位置，O為海底水平陣參考陣元的位置，Ai為某待校準(zhǔn)陣元的位置。

圖1 校準(zhǔn)方法示意圖 Fig.1 Schematic diagram of calibration method

已知聲 源S1(X1,Y1)、S2(X2,Y2)與參考陣元O(x1,y1)的GPS 位置信息，通過(guò)估計(jì)各陣元較參考陣元的信號(hào)到達(dá)時(shí)延，以求解Ai點(diǎn)坐標(biāo)(xi,yi)：

選取目標(biāo)函數(shù)如下：

通過(guò)優(yōu)化搜索，可得Ai點(diǎn)坐標(biāo)：

由于淺海信道存在縱向相關(guān)振蕩現(xiàn)象[5]，同時(shí)受海洋環(huán)境噪聲等因素影響，簡(jiǎn)單地以參考陣元為基準(zhǔn)估計(jì)時(shí)延存在誤差。

1.3 時(shí)延估計(jì)

任意兩個(gè)陣元間的信號(hào)時(shí)延可以根據(jù)相關(guān)峰的位置來(lái)估計(jì)：

其中：上標(biāo)*表示取共軛，最大相關(guān)系數(shù)對(duì)應(yīng)的δ即為陣元Au相對(duì)于陣元Aw接收到聲源Sk的信號(hào)時(shí)延估計(jì)值。

測(cè)量全部陣元組合后，可得兩兩陣元間的時(shí)延估計(jì)矩陣Gk：

1.3.1 時(shí)延累計(jì)測(cè)量

在高信噪比和信號(hào)空間相關(guān)較好的理想條件下，文獻(xiàn)[4]近似為通過(guò)兩兩相鄰陣元的信號(hào)相關(guān)來(lái)估計(jì)時(shí)延，累計(jì)可得第i個(gè)陣元相對(duì)參考陣元的時(shí)延量為

進(jìn)而估計(jì)陣形參數(shù)。

1.3.2 時(shí)延冗余測(cè)量

當(dāng)信噪比偏低、或者聲場(chǎng)相關(guān)存在起伏的情況下，僅僅測(cè)量(I? 1)個(gè)相鄰陣元的時(shí)延并沒(méi)有充分利用信號(hào)的空間信息。通過(guò)全部?jī)蓛申囋M合進(jìn)行次時(shí)延估計(jì)，然后由最小二乘方法計(jì)算第i個(gè)陣元相對(duì)參考陣元的時(shí)延量，可更好地對(duì)時(shí)延進(jìn)行估計(jì)：

1.4 陣形校準(zhǔn)誤差

定義陣形校準(zhǔn)誤差為

2 數(shù)值仿真

2.1 波導(dǎo)環(huán)境參數(shù)

通過(guò)數(shù)值仿真來(lái)對(duì)比上述兩種時(shí)延估計(jì)方法用于陣形校準(zhǔn)的差異。根據(jù)實(shí)驗(yàn)中測(cè)量的數(shù)據(jù)來(lái)設(shè)置聲場(chǎng)仿真所用的信道參數(shù)，接收陣位置處實(shí)測(cè)聲速剖面如圖2 所示，水層介質(zhì)密度設(shè)為1 g·cm-3，海深設(shè)為 96 m，海水的聲吸收衰減系數(shù)設(shè)為0 dB·λ-1。海底假設(shè)為平坦半無(wú)限基底，聲速1 630 m·s-1，密度為 1.76 g·cm-3，衰減系數(shù)0.33 dB·λ-1。參考聲速選取為1 568 m·s-1，聲源深度設(shè)為50 m，接收陣列布放于海底。仿真信號(hào)頻段選為20～30 Hz，分辨率為0.2 Hz，采樣率設(shè)置為16 000 Hz，使用Kraken 程序[7]計(jì)算聲場(chǎng)。

圖2 實(shí)測(cè)的聲速剖面 Fig.2 A measured sound velocity profile

2.2 陣形校準(zhǔn)

根據(jù)文獻(xiàn)[4]校準(zhǔn)聲源的最佳位置關(guān)系，布設(shè)校準(zhǔn)聲源位置如下：以x軸正半軸為0°方向，逆時(shí)針?lè)较蛐D(zhuǎn)角度遞增，S1在參考陣元O的0°方向，距離為20 km；S2在參考陣元O的270°方向，距離為20 km。

信噪比為10 dB 時(shí)陣形校準(zhǔn)結(jié)果如圖3 所示，校準(zhǔn)陣列是陣元間距10 m 的11 元直線陣。

圖3 信噪比為10 dB 時(shí)兩種方法的陣形校準(zhǔn)結(jié)果 Fig.3 Array shape calibration results of the two methods when the SNR is 10 dB

由圖3 中的對(duì)比可以看出，基于時(shí)延冗余測(cè)量的陣形校準(zhǔn)結(jié)果明顯優(yōu)于基于時(shí)延累積測(cè)量的校準(zhǔn)結(jié)果。進(jìn)一步考察不同信噪比下，兩種時(shí)延估計(jì)方法的陣形校準(zhǔn)精度。給接收信號(hào)添加不同大小的高斯白噪聲，并進(jìn)行100 次蒙特卡洛仿真實(shí)驗(yàn)。圖4 給出了信噪比在5～30 dB 范圍內(nèi)變化時(shí)，兩種方法陣形校準(zhǔn)誤差的均值與標(biāo)準(zhǔn)差。

圖4 不同信噪比下兩種方法的陣形校準(zhǔn)誤差的均值與標(biāo)準(zhǔn)差 Fig.4 Means and standard deviations of the array shape calibration errors of the two methods under different signal to noise ratios

從圖4 可以看出，在信噪比大于30 dB 時(shí)，時(shí)延冗余測(cè)量方法和時(shí)延累積測(cè)量方法都可以準(zhǔn)確地實(shí)現(xiàn)陣形校準(zhǔn)，而信噪比降低時(shí)，時(shí)延冗余測(cè)量的優(yōu)勢(shì)更為明顯，不僅誤差均值小，而且標(biāo)準(zhǔn)差也更小。如以誤差校準(zhǔn)精度0.1 m 為門(mén)限，冗余測(cè)量方法對(duì)信號(hào)的信噪比要求比累積測(cè)量方法低約10 dB。

3 實(shí)驗(yàn)數(shù)據(jù)處理

3.1 實(shí)驗(yàn)描述

2020 年9 月某近海實(shí)驗(yàn)中，8 元水平接收陣布放于海底，陣列的信號(hào)采樣率為16 kHz，靈敏度為?170 dB，實(shí)驗(yàn)船沿著測(cè)線投擲爆炸聲源，標(biāo)稱(chēng)深度為50 m。實(shí)驗(yàn)期間的聲速剖面、波導(dǎo)環(huán)境等參數(shù)與上文數(shù)值仿真所用參數(shù)一致。兩輔助聲源的布設(shè)位置與數(shù)值仿真接近，S1在參考陣元O的358.6°方向、距離19.67 km 處；S2在參考陣元O的269.4°方向、距離20.78 km 處，聲源的位置信息由船載差分GPS 獲得。由文獻(xiàn)[4]可知，當(dāng)校準(zhǔn)聲源按照最佳位置關(guān)系布放時(shí)，聲源位置的偏差對(duì)陣形校準(zhǔn)結(jié)果的影響較小。

3.2 處理結(jié)果

對(duì)接收信號(hào)進(jìn)行20～30 Hz 帶通濾波處理，各陣元接收到的帶內(nèi)信號(hào)波形如圖5 所示，圖5（a）為對(duì)S1的接收信號(hào)，圖5（b）為對(duì)S2的接收信號(hào)。

圖5 陣列接收的來(lái)自?xún)蓚€(gè)聲源的信號(hào)波形 Fig.5 Signal waveforms received by the array from two sources

對(duì)濾波后信號(hào)分別用兩種方法進(jìn)行時(shí)延估計(jì)，得到陣形校準(zhǔn)結(jié)果如圖6 所示，可以看出校準(zhǔn)后陣形與實(shí)驗(yàn)預(yù)設(shè)的直線陣形差別較大。

圖6 兩種方法進(jìn)行接收信號(hào)時(shí)延估計(jì)得到的的陣形校準(zhǔn)結(jié)果 Fig.6 Array shape calibration results obtained by using the two methods to estimate the time delay of array received signals

為了驗(yàn)證校準(zhǔn)陣形的有效性，對(duì)位于[參考陣元O269.4°方向、距離21.98 km]的測(cè)試聲源Stest的接 收信號(hào)進(jìn)行波束掃描，并與實(shí)際結(jié)果進(jìn)行比較。

圖7 給出了三種陣形數(shù)據(jù)對(duì)測(cè)試聲源信號(hào)的常規(guī)波束形成（Conventional Beamformer,CBF）結(jié)果，為了更好體現(xiàn)陣形準(zhǔn)確性對(duì)目標(biāo)測(cè)向的影響，選擇頻段為200～300 Hz。圖7（a）為根據(jù)實(shí)驗(yàn)預(yù)設(shè)陣形CBF 波束掃描結(jié)果，其中目標(biāo)的方位角估計(jì)為88.6°與271.4°，不具備左右舷分辨能力。圖7（b）為根據(jù)冗余測(cè)量時(shí)延估計(jì)陣形校準(zhǔn)后CBF 波束掃描結(jié)果，其中目標(biāo)的方位角估計(jì)結(jié)果為269.4°。圖7（c）為根據(jù)相鄰累積時(shí)延估計(jì)陣形校準(zhǔn)后CBF 波束掃描結(jié)果，目標(biāo)的方位角估計(jì)結(jié)果為269.5°。由GPS 測(cè)量得到的方位角為269.4°，因此可以看出校準(zhǔn)后陣形CBF 測(cè)向更準(zhǔn)確，驗(yàn)證了所提方法的有效性。

圖7 對(duì)圖6 中三種不同陣形的常規(guī)波束形成結(jié)果 Fig.7 Conventional beamforming results for the three different array shapes shown in Fig.6

3.3 進(jìn)一步處理

從圖5 可以看出，由于各陣元接收信號(hào)的信噪比很高，兩種陣形校準(zhǔn)方法的差別并不明顯。對(duì)接收信號(hào)添加10 dB 的高斯白噪聲，加噪聲后陣列接收信號(hào)的時(shí)域波形如圖8 所示。進(jìn)行50 次蒙特卡洛仿真實(shí)驗(yàn)，以分析兩種方法的陣形校準(zhǔn)結(jié)果的穩(wěn)定性。陣形校準(zhǔn)結(jié)果如圖9 所示，當(dāng)接收信號(hào)信噪比降低后，圖9（a）中充分利用了信號(hào)的空間信息進(jìn)行冗余時(shí)延估計(jì)，相較圖9（b）中時(shí)延累積方法的陣形校準(zhǔn)結(jié)果，校準(zhǔn)得到的陣元位置更加穩(wěn)定。

圖8 加10 dB 噪聲后陣列接收的來(lái)自?xún)蓚€(gè)聲源的信號(hào)波形 Fig.8 Signal waveforms received by the array from two sources after adding 10 dB noise

圖9 加10 dB 噪聲后的陣形校準(zhǔn)結(jié)果 Fig.9 Array shape calibration results obtained by using the two methods to estimate the time delay of array received signals after adding 10 dB noise

4 結(jié)論

本文提出了一種基于時(shí)延冗余測(cè)量的陣形校準(zhǔn)方法。數(shù)值仿真結(jié)果顯示，在信噪比偏低時(shí)，相較相鄰陣元時(shí)延累積的陣形校準(zhǔn)方法，所提的基于時(shí)延冗余測(cè)量的陣形校準(zhǔn)方法，更充分地利用了陣元的空間信息，可以提高陣形校準(zhǔn)精度與穩(wěn)定性。通過(guò)2020 年某近海實(shí)驗(yàn)的數(shù)據(jù)分析，對(duì)所提基于時(shí)延冗余測(cè)量的陣形校準(zhǔn)方法進(jìn)行了驗(yàn)證，用校準(zhǔn)后的陣形進(jìn)行波束掃描得到的聲源方位角與實(shí)際值更接近，表明了所提方法可以更準(zhǔn)確地給出各陣元相對(duì)參考陣元的位置，驗(yàn)證了該方法的有效性。

致謝：感謝參加實(shí)驗(yàn)的全體工作人員，他們的辛苦工作為本文提供了寶貴的實(shí)驗(yàn)數(shù)據(jù)。