三維前視聲吶信號處理方法

劉治宇，劉曉東，董飛，劉小剛

?

三維前視聲吶信號處理方法

劉治宇，劉曉東，董飛，劉小剛

(中國科學院聲學研究所海洋聲學技術實驗室，北京100190)

三維前視聲像聲吶是安裝在小型水下載體上的重要聲學探測設備。提出了基于波達方向估計技術的三維前視聲吶信號處理方法，接收陣水平方向采用波束形成技術，垂直方向采用波達方向估計技術。在波束內采用分裂孔徑相位法提高水平方向分辨率，形成聲吶陣前方的三維聲像。仿真結果和水池實驗結果表明：在接收陣面面積相同的條件下，能夠獲得優于常規技術的分辨率，可以實現點目標和連續目標的探測，適合于在水下載體上安裝使用。

波束形成；波達方向估計；相位法；三維前視聲吶

0 引言

隨著海洋開發和利用越來越被重視，應用于近海底觀測的各種水下載體發展迅速，如自治無人水下機器人、遙控水下機器人和載人水下機器人等。三維前視聲像聲吶用于探測載體前方三維空間內障礙物、目標和海底，并進行聲成像。由于水下載體速度的要求，以及空間和功耗的限制，安裝在小型水下載體上的三維前視聲像聲吶應該具有體積小、重量輕、功耗低和分辨率高的特點。

傳統的三維前視聲吶主要采用接收平面陣的二維波束形成技術，其分辨率受波束角寬限制，要求陣面尺度大，系統復雜[1]。如果采用二維波達方向估計技術雖然提高了分辨率，但是在空間中同時到達的回波個數較多的復雜情況下，不能有效地分辨出多個目標[2]。文獻[3]利用水平方向波束形成技術和垂直方向波達方向技術，垂直方向的波達方向(Direction of Arrival, DOA)估計技術可分辨多個目標，適合于近水面和近海底的應用，可以有效提高垂直方向的分辨率，但水平方向的分辨率較低。

為了克服以上幾種方法的缺點，本文采用了一種波束形成技術、相位法和波達方向估計技術相結合的信號處理方法。通過水平方向的波束形成將波束指向特定的方向，提高該方向的信噪比并減少同時到達聲吶陣的信源個數；在垂直方向采用波達方向估計技術，提高垂直方向的分辨率；并在水平方向的波束內，采用分裂孔徑相位法提高對連續目標的水平方向分辨能力。利用這樣的信號處理方法，可以在陣面尺寸較小的情況下，有效地探測聲吶陣前方目標和地形，并提高水平和垂直兩個方向的分辨率，使用該方法的聲吶陣更適合安裝在小型水下載體上。

1 信號處理方法

本文提出的信號處理方案是針對安裝在小型水下載體上的三維前視聲像聲吶，其信號處理流程中關鍵技術是結合波束形成技術、相位法和波達方向估計技術，信號處理總體流程如圖1所示。在水平平面內采用波束形成技術形成多個水平方向的波束，波束寬度約為7°；在波束內，采用相位法求出水平方向的入射角，提高探測連續目標的水平方向分辨率；同時在垂直平面內采用波達方向估計技術求出垂直方向的入射角，最后利用這兩個方向的入射角度和回波到達時間等信息計算出散射點在三維空間中的位置，得到聲像。

由于在波束內采用相位法，可以在較寬波束內提高水平方向分辨率，因此可減少水平基元個數，減小陣面尺寸。

下面分別對其中所使用的各種信號處理技術進行闡述。

1.1 水平方向波束形成

頻域波束形成采用快速傅里葉變換(Fast Fourier Transform, FFT)同時形成多個指向不同角度的波束，頻域波束形成的計算方法為[4]

與延時加權求和的時域波束形成相比，采用FFT波束形成方法能夠加快計算速度，提高計算效率，但是FFT中波束序號是整數，所以FFT只能形成有限個特定角度的波束。由于水下載體的姿態通常比較好，不需要做實時的姿態穩定，因此只要這些波束能夠覆蓋聲吶陣前方區域，固定的波束形成角度對聲成像不會產生影響。

波束形成技術的作用類似于空間濾波，可以保留波束指向角方向的增益并壓制偏離波束指向角的散射點，所以只需要考慮波束中心在海底投影附近的散射點的影響。波束主瓣寬度與波束指向角的大小有很大關系，波束指向正下方，即波束指向角為零時，波束主瓣最窄；而當波束指向角變大時，波束主瓣在海底的投影變寬，而且投影在波束中心的兩側不對稱。如果直接采用波束形成后的結果進行判別，則測深結果的分辨率較差。

1.2 分裂孔徑相位法

分裂孔徑相位法的基本原理[5-8]是：把接收陣分成兩個完全相同的子陣，子陣波束形成后，等效于將接收陣面轉向波束指向角，波束中心上的散射元正好位于兩個子陣的等效中心連線的垂直平分線上，也就是說波束中心點到達兩個子陣的相位差為零。

而偏離波束中心的水平入射角也可以利用它與兩子陣相位差之間的對應關系計算得到。令兩個子波束的同相分量和正交分量分別為1，2，1，2，可以利用共軛相乘的方法來提取相位差，計算公式為

(3)

其中：是兩個子陣間隔的基元間隔；是波數；是波束指向角；是實際角度與波束指向角的偏差。在遠場條件下，相位差為零時的方向對應于波束指向角，此時有。所以由相位法計算出水平方向的實際角度為

利用上述相位法估計回波入射角的優點是能夠有效地利用波束內幾乎所有相位差信息來估計所有散射點的回波入射角度，對于連續目標能夠得到比較精確的角度估計結果。但是它的缺點是只能分辨波束內同時到達的一個目標，無法分辨同時到達的多個目標。此外，對于聲吶陣正前方的目標，相位差變化的區間太小，相位差曲線任何一點細微的抖動都會給最終的估計結果帶來較大的誤差，所以相位法不適用于正前方的目標。

1.3 信源數估計

波達方向估計算法多采用基于特征子空間的算法，它利用信號子空間與噪聲子空間的正交性估計來波方向，需要預先知道信源個數。而當信號源數估計不正確時，也就是對信號子空間和噪聲子空間估計不準，即兩者不完全正交，會造成對估計信號源時的虛警或者漏警，也就會造成在估計信號方向時的偏差，所以波達方向估計技術中的首要問題是信源數估計問題。

文獻[9]提出了一種基于特征空間的信源數估計方法，它將陣列信號的協方差估計值分別投影到信號的特征子空間和噪聲子空間。由于信號子空間與噪聲子空間相互正交，易于由表征投影大小的判據值區分信號和噪聲的貢獻，得到了性能優良的信源數估計方法。該方法優于傳統的蓋氏圓法，而且不需要人工輸入信號閾值，可以實現自動的閾值設定，具有一定的實用價值。

1.4 波達方向估計

在垂直平面內對做完波束形成后的信號用波達方向估計技術求入射角，系統中采用旋轉不變子空間(Estimation of Signal Parameters via Rotational Invariance Techniques, ESPRIT)算法，利用數據協方差矩陣信號子空間的旋轉不變特性估計信號參數。與其他方法相比，ESPRIT算法不需要進行譜峰搜索，計算量比較小，計算速度比較快。

ESPRIT算法的原理[2]是：認為相鄰子陣間存在一個固定間距，這個固定間距能反映出各相鄰子陣間的旋轉不變特性。ESPRIT算法就是利用這個子陣間的旋轉不變性實現陣列的DOA的估計。

式中：是發射信號，是空間陣列的流型矩陣，旋轉不變關系，則子陣1的陣列流型，子陣2的陣列流型，是兩個子陣陣列流型的合并形式，噪聲一般是零均值高斯白噪聲，與信號不相關。

對接收信號的協方差矩陣進行特征分解可得：

(7)

本文采用實值空間總體最小二乘的ESPRIT算法[2](Unitary TLS-ESPRIT)，對常規ESPRIT算法的改進主要表現在以下幾個方面：

(2) 進行雙向平滑，減小估計誤差，顯著增強了對相關性較強的回波信號的處理。可以進行對估計結果的有效性判決。

2 仿真結果

為了驗證將波束形成技術、相位法和波達方向估計技術三者相結合的信號處理算法的性能，在Windows平臺上用Matlab 7.0進行仿真計算。仿真實驗中，平面接收聲吶陣選用16×6個陣元，相鄰陣元間隔為半波長，有效陣面大小為80 mm×30 mm。發射信號帶寬為10 kHz、中心頻率為150 kHz的chirp信號；實驗目標設定為聲吶陣前方三維空間中的一根桿，桿中心點P到參考點的距離為75 m；水平方向入射角為45°，垂直方向入射角為30°，信噪比為20 dB。產生仿真數據的方法見參考文獻[10]。

采用本文提出的信號處理方法仿真計算出的角度如圖2所示。與實際散射點的入射角度相比，水平方向入射角的均方根誤差為1.1037°，垂直方向入射角的均方根誤差為0.9141°。

對于有效陣面大小相同的平面聲吶陣，分別采用結合波束形成技術、相位法和波達方向估計技術的新方法與采用常規二維波束形成方法在入射角度估計上的誤差比較見表1。可以看出，結合波束形成技術、相位法和波達方向估計技術的新方法能夠獲得較高的水平分辨率和垂直分辨率，優于常規二維波束形成技術。

表1 兩種方法計算出角度的誤差

3 水池試驗結果

聲吶陣元布設圖見圖3，發射陣的中心基元用于發射，周圍8個基元為啞元，為中心發射基元提供必要的邊界條件，保證中心基元能夠實現寬覆蓋發射，水平、垂直方向波束開角可達到85°×85°。單基元發射能夠實現寬覆蓋，但發射靈敏度低，影響系統的作用距離。

接收陣為18×8基元的平面陣，中間16×6個基元接收，邊緣基元為啞元，每個陣元單獨引線。它可在水平平面內形成約7°的窄波束。聲吶陣尺寸較小，有效陣面尺寸為80 mm×30 mm，非常適合安裝在水下載體頭部。

為了驗證方法的有效性，課題組在某消聲水池開展了目標探測試驗和池壁探測試驗。該水池長約13 m，寬為7.5 m，深為5 m，其各側面及上下表面均覆蓋有消聲瓦。試驗時將聲吶艙固定于航車的吊架上，通過該吊架可以根據需要在水池內任意部位吊放聲吶系統。本實驗將聲吶艙吊放于水池中間且距離水池后壁約1.5 m處，入水深度為2.5 m。利用CPCI信號處理機采集并存儲原始AD數據，采用上述信號處理方法對原始AD數據處理得到聲吶陣前方圖像。

目標探測試驗選用兩個典型目標。目標1為空心圓柱，直徑為18 cm，高度為30 cm，吊放于聲吶陣右前方。目標2為空心球，直徑為33 cm，吊放于聲吶陣正前方。兩目標均懸于水中，距池底2.5 m，中心高度相同，水平距離為0.4 m，試驗具體布設情況見圖4。從圖5的圖像結果可以看出，該信號處理方法可以分辨出兩個目標，圖像中目標的位置與實際相吻合，水平方向采用的相位法技術可以恢復目標的多個散射點，具有高的水平分辨率。

池壁探測試驗的結果如圖6、7所示。由于水池前壁布有吸聲尖劈，部分散射點的回波較弱，圖像點呈不均勻分布，見圖6。但是從上方看，所有圖像點基本在同一個平面，見圖7。池壁探測試驗表明，該方法可以恢復出連續目標的圖像，具有高的水平分辨率和垂直分辨率，聲吶測距誤差為18.78 cm。利用池壁探測試驗采集的原始AD數據，采用二維波束形成技術得到的測距誤差為37.48 cm。如利用文獻[3]提出的方法，則測距誤差為31.63 cm。從以上比較可以看出，該方法在聲吶陣有效陣面比較小的條件下，能夠獲得優于常規技術的水平分辨率和垂直分辨率。

4 結論

應用在小型水下載體的三維前視聲像聲吶有體積小和分辨率高這兩個方面的技術要求，本文提出了一種將波束形成技術、分裂孔徑相位法和波達方向估計技術相結合的方法。水池試驗結果表明，該方法在聲吶陣有效陣面比較小的條件下，能夠獲得優于常規技術的水平分辨率和垂直分辨率，測距誤差小于20 cm，適合在小型水下載體上應用。

[1] 李啟虎. 聲吶信號處理引論[M]. 2版, 北京: 海洋出版社, 2000, 167-169.

[2] 王永良. 空間譜估計理論與算法[M]. 北京: 清華大學出版社, 2004, 26-28.

[3] Kraeutner P, Brumley B, Guo H, et al. Rethinking forward-looking sonar for AUV’s: combining horizontal beamforming with vertical angle-of-arrival estimation[C]// Oceans 2007 MTS/IEEE Conference proceedings, 2007: 1-7.

[4] Roger G Pridham, Ronald A Mucci. A novel approach to digital beamforming[J]. J Acoust Soc Am, 1978, 63(2): 425-435.

[5] Luren Yang, Torfinn Taxt. Multibeam sonar bottom detection using multiple subarrays[C]// Proc Oceans’ 97, 1997: 932-938.

[6] Yoshida Z, Asada A, Ikeda Y, et al. High precision survey by the multi-beam sonar in the dam site[C]// OCEANS 2004 MTS/IEEE Conference Proceedings, 2004: 1133-1138.

[7] Fairfiel N, Wettergreen D. Active localization on the ocean floor with multibeam sonar[C]// OCEANS 2008 MTS/IEEE Conference Proceedings, 2008: 1-10.

[8] 李海森, 周天, 徐超. 多波束測深聲吶技術研究新進展[J]. 聲學技術, 2013, 32(2): 73-80.

LI Haisen,ZHOU Tian, XU Chao. New developments on the technology of multi-beam bathymetric sonar[J]. Technical Acoustics, 2013, 32(2): 73-80.

[9] 朱維慶, 胡雋, 劉曉東, 等. 波達方向估計中的特征空間的信源數估計方法[J]. 聲學學報, 2009, 34(3): 97-102.

ZHU Weiqing, HU Juan, LIU Xiaodong, e al. Source number estimation using eigenspace in direction of arrival estimate[J]. Acta Acustica, 2009, 34(3): 97-102.

[10] Gueriot D, Sintes C. Forward looking sonar data simulation through tube tracing[C]// OCEANS 2010 IEEE Sydney, 2010: 1-6.

Signal processing method for 3D forward-looking sonar

LIU Zhi-yu, LIU Xiao-dong, DONG Fei, LIU Xiao-gang

(Laboratory of Acoustic Technology,Institute of Acoustics, Chinese Academy of Sciences,Beijng100190, China)

A signal processing method for 3D forward-looking Sonar based on Direction of Arrival Estimation is presented. Beamforming is used in horizontal direction of the receive array, and DOA estimation is employed in vertical direction. Phase estimation is also implemented to improve the resolution in horizontal direction for detecting continuous target. Finally the 3D image in the front of sonar array is formed. The results of tests show that the method can get higher resolution than the sonar using beamforming with the same area of receiver array. The sonar using the method is more suitable for underwater vehicles.

beamforming;direction of arrival estimation;phase estimation; 3D forward-looking sonar

TB565

A

1000-3630(2015)-04-0317-05

10.16300/j.cnki.1000-3630.2015.04.005

2014-07-16;

2014-10-12

國家863計劃項目(2006AA09Z119)

劉治宇(1978－), 男, 遼寧沈陽人, 博士, 副研究員, 研究方向為水聲信號處理、陣列信號處理。

劉治宇, E-mail: liuzhiyu@mail.ioa.ac.cn