劉豫魯 劉清宇 劉雄厚，2 余 赟

（1.海軍研究院 北京 100161）（2.中科院聲學(xué)研究所 北京 100190）

淺海大孔徑線列陣子陣陣增益分析與評(píng)估?

劉豫魯1劉清宇1劉雄厚1，2余 赟1

（1.海軍研究院 北京 100161）（2.中科院聲學(xué)研究所 北京 100190）

聲信號(hào)在海洋中傳播時(shí)，由于受到復(fù)雜海洋環(huán)境的影響，波形會(huì)發(fā)生畸變。理想情況下常規(guī)波束形成的陣增益可以達(dá)到10logN，但是由于波形畸變，到達(dá)線列陣各陣元處的信號(hào)相關(guān)性下降，從而導(dǎo)致陣增益低于10logN。論文針對(duì)典型淺海環(huán)境，研究窄帶信號(hào)經(jīng)過淺海傳播，到達(dá)大孔徑線列陣各陣元處的信號(hào)水平縱向相關(guān)性，并分析此時(shí)大孔徑線列陣進(jìn)行子陣處理時(shí)的檢測(cè)性能。以常規(guī)波束形成為例，計(jì)算大孔徑線列陣的陣增益并分析影響陣增益的因素。

信號(hào)空間相關(guān)性；常規(guī)波束形成；陣增益

1 引言

在無界均勻介質(zhì)中，從一點(diǎn)發(fā)出的信號(hào)，將無畸變地傳播到各處，任何兩點(diǎn)接收的信號(hào)都是相關(guān)的。但是在實(shí)際海洋中由于海水介質(zhì)及其邊界在時(shí)間與空間上存在著有規(guī)律的與隨機(jī)起伏的變化，使聲場(chǎng)的空間相關(guān)性降低，進(jìn)而導(dǎo)致陣列處理增益下降［1］。因此研究信號(hào)的空間相關(guān)特性對(duì)提高信號(hào)處理的性能有重要意義。科研人員從20世紀(jì)70年代就開始對(duì)空間相關(guān)性進(jìn)行研究。研究初期得到一些實(shí)驗(yàn)和理論的結(jié)果，但沒有建立起與波導(dǎo)環(huán)境密切相關(guān)的能夠解釋客觀規(guī)律的較為完整的理論體系［2］。研究后期隨著聲傳播理論的逐漸完善，人們利用水聲波動(dòng)理論對(duì)聲場(chǎng)空間特性進(jìn)行分析［3～8］。陣增益定義為基陣輸出端的信噪比與單個(gè)陣元輸出信噪比的比值，是衡量傳感器陣列波束形成器性能的重要指標(biāo)之一［9～10］。本文以典型淺海環(huán)境為例，研究大孔徑線列陣陣元間的信號(hào)相關(guān)特性，分析影響信號(hào)空間相關(guān)特性的因素，通過仿真得出陣列處理增益，為大孔徑線列陣子陣處理能力的評(píng)估提供基礎(chǔ)。

2 信號(hào)空間相關(guān)性

2.1 信號(hào)空間相關(guān)性定義

從一點(diǎn)發(fā)出的信號(hào)，經(jīng)過傳播到達(dá)線列陣各陣元處，同一時(shí)刻兩個(gè)陣元接收到的信號(hào)的相似程度定義為信號(hào)的空間相關(guān)性。線陣列的信號(hào)空間相關(guān)性分為垂直相關(guān)性和水平相關(guān)性，水平相關(guān)性又分為縱向相關(guān)性和橫向相關(guān)性。本文討論的是水平縱向相關(guān)性，其大小由相關(guān)系數(shù)來衡量。如果是在無界均勻介質(zhì)中，那么到達(dá)任何兩點(diǎn)處的信號(hào)都是相關(guān)的，但是實(shí)際的海洋環(huán)境是不均勻且存在邊界的，所以信號(hào)波形會(huì)發(fā)生畸變，兩點(diǎn)間的信號(hào)相關(guān)性會(huì)下降，我們把相關(guān)系數(shù)下降到某一數(shù)值時(shí)（一般為0.707）對(duì)應(yīng)的縱向間隔定義為相關(guān)長(zhǎng)度［11］。位于同一深度且具有一定間隔的兩陣元(r,Zr)，(r+l,Zr+l)在同一時(shí)刻接收信號(hào)波形之間的互相關(guān)函數(shù)最大值與信號(hào)能量的比值為水平縱向相關(guān)系數(shù)，表達(dá)式為

其中Pr(t)和Pr+l(t)為接收信號(hào)，Pr+l(t+τ)為經(jīng)過延遲τ的信號(hào)Pr+l(t)，經(jīng)過傅里葉變換可得到頻域內(nèi)的表達(dá)式

其中，Pr(ω)和 P(ω)分別為兩陣元接受信號(hào)的頻譜，ω表示角頻率，ω1和ω2分別為窄帶信號(hào)的截止角頻率。

2.2 水平縱向相關(guān)性仿真

仿真環(huán)境采用一個(gè)典型的水平分層淺海波導(dǎo)環(huán)境，聲速剖面及海底底質(zhì)參數(shù)如圖1所示。

圖1 海洋波導(dǎo)環(huán)境

聲源信號(hào)為一帶寬10Hz的窄帶信號(hào)，中心頻率為380Hz。陣元間距為2m，陣元個(gè)數(shù)為100個(gè)。仿真中選取A，B，C，D四個(gè)點(diǎn)作為距離聲源最近的一號(hào)陣元的位置，他們和聲源的位置關(guān)系如圖2所示。

根據(jù)式（2）進(jìn)行仿真，得到線列陣一號(hào)陣元分別位于A，B，C，D四個(gè)位置時(shí)信號(hào)到達(dá)一號(hào)陣元和其他陣元的信號(hào)水平縱向相關(guān)性如下。

東北的冬儲(chǔ)在發(fā)生改變，從由大批發(fā)商、基層零售商、農(nóng)民組成的鏈條式儲(chǔ)備轉(zhuǎn)向了由大批發(fā)商、中小型復(fù)合肥企業(yè)組成的上游儲(chǔ)備。目前，東北地區(qū)基礎(chǔ)化肥缺口較大，但化肥價(jià)格繼續(xù)上漲的動(dòng)力不足，后市以窄幅波動(dòng)為主。

圖2 位置關(guān)系

圖3 A，B，C，D處的信號(hào)水平縱向相關(guān)性

可以看出：

1）距離聲源水平距離相同的A，C兩點(diǎn)，盡管所處深度不同，但相關(guān)半徑相差不大。

2）深度相同的A和B，C和D，A，C相對(duì)于B，D兩點(diǎn)相關(guān)半徑更大。

3）線列陣的一號(hào)陣元位于A，B，C，D四個(gè)位置時(shí)相關(guān)半徑分布在40m～52m的范圍內(nèi)。

3 陣增益

3.1 陣增益的定義及其計(jì)算方法

使用基陣的主要目的之一就是通過相干疊加信號(hào)和非相干疊加噪聲來改善基陣接收端的信噪比。改善的程度用陣增益來衡量，陣增益定義為基陣輸出信噪比和單個(gè)陣元信噪比的比值：

遠(yuǎn)場(chǎng)平面波情況下，由于存在信號(hào)到達(dá)方向角的問題，所以信號(hào)到達(dá)各個(gè)陣元的時(shí)間不同，信號(hào)產(chǎn)生延遲，存在相位差。常規(guī)波束形成就是通過設(shè)定合適的加權(quán)系數(shù)，對(duì)各陣元數(shù)據(jù)進(jìn)行延遲求和的過程。常規(guī)波束形成器加權(quán)向量w表示為

wi=exp(-jω(i-1)Δt)，i=1,2,…,N （4）其中，Δt=dcosθ/c，d為陣元間距，c為海水中聲速，θ為信號(hào)到達(dá)方向角，ω為信號(hào)角頻率，N為陣元個(gè)數(shù)。

考慮噪聲為各向同性白噪聲，則常規(guī)波束形成算法的陣增益可以表示為

式中ni為第i個(gè)陣元接收到的噪聲信號(hào)，則第i個(gè)陣元接收到的信號(hào)xi表示為xi=si+ni。

3.2 實(shí)驗(yàn)仿真

我們?cè)谇懊嫜芯康幕A(chǔ)上，根據(jù)式（5）對(duì)線列陣位于A，B兩點(diǎn)處的陣增益進(jìn)行仿真。對(duì)于A，B兩點(diǎn)處的信號(hào)到達(dá)角θ，需要進(jìn)行估計(jì)，而不是圖中簡(jiǎn)單的幾何關(guān)系。對(duì)于單個(gè)陣元的信噪比，我們用單個(gè)陣元信噪比的平均值來代替。最后得到線列陣一號(hào)陣元位于A，B兩點(diǎn)處的陣增益如下。

圖4 A，B兩點(diǎn)處的陣增益

從圖中可以看出

1）線列陣一號(hào)陣元位于A，B兩處時(shí)，隨著陣元數(shù)的增加，陣增益不斷增大，但增加到一定數(shù)值后再增加陣元數(shù)目，陣增益基本不變。線列陣位于兩處時(shí)的陣增益都小于10logN。

2）線列陣一號(hào)陣元位于A，B兩處時(shí)，位于A處時(shí)的陣增益稍大于位于B處時(shí)的陣增益。

根據(jù)文獻(xiàn)［12］中推導(dǎo)出的基陣信號(hào)增益為

式中M為水平接收陣的陣元個(gè)數(shù)，ρ是任意兩接收點(diǎn)之間聲場(chǎng)空間相關(guān)系數(shù)矩陣，h是一個(gè)由基陣各陣元上 p組成的向量，而 p的大小由各陣元上的ATL（平均傳播損失）決定?？梢钥闯?，信號(hào)增益由M，ρ，h三個(gè)量決定。A點(diǎn)和B相比，水平縱向相關(guān)系數(shù)更大，傳播損失更小，因此陣增益大于B點(diǎn)。

4 結(jié)語

本文以典型淺海環(huán)境為例，研究一窄帶信號(hào)經(jīng)過傳播到達(dá)大孔徑線列陣各陣元處的信號(hào)水平縱向相關(guān)性及線列陣的陣增益，發(fā)現(xiàn)與線列陣的深度相比，線列陣與聲源的距離對(duì)線列陣各陣元處的信號(hào)水平縱向相關(guān)性影響更大。仿真實(shí)驗(yàn)中A，C兩點(diǎn)相對(duì)于各自同深度的B，D兩點(diǎn)的相關(guān)半徑更大，線列陣位于四個(gè)點(diǎn)處的相關(guān)半徑分布在40m-52m的范圍內(nèi)。對(duì)于常規(guī)波束形成，影響陣增益的因素包括M，ρ，h三個(gè)量，仿真實(shí)驗(yàn)中，A點(diǎn)和B相比，水平縱向相關(guān)系數(shù)更大，傳播損失更小，陣增益大于B點(diǎn)，但都小于理想情況下的陣增益。仿真結(jié)果與理論相符合。

［1］張仁和，張雙榮，肖金泉，等.淺海遠(yuǎn)程聲場(chǎng)的空間相關(guān)性與時(shí)間穩(wěn)定性［J］.聲學(xué)學(xué)報(bào)，1981，6（1）：9-19.

［2］尚爾昌.隨機(jī)淺海信道中信號(hào)場(chǎng)橫向相干理論［J］.聲學(xué)學(xué)報(bào)，1979，4（3）：30-31.

［3］Rozenfeld I，Cable P，Carey W M.Spatial coherence in shallow water［C］∕Halifax：IEEE，1997：597-601.

［4］Zhu R.Spatial horizontal coherence of sound in shallow water［J］.J Acoust Soc Am，1992，92（2）：956-961.

［5］Ballard M S，Lin Y T，Lynch J F.Horizontal refraction of propagating sound due to seafloor scours over a range-dependent layered bottom on the New Jersey shelf［J］.J Acount Soc AM，2012，131（1）：2587-2598.

［6］Knobles D P，Stotts S A，Koch R A.Low frequency coupled mode sound propagation over a continental shelf［J］.J Acount Soc Am，2003，113（1）：781-787.

［7］Duda T F，Collis J M，Lin Y T.Horizontal coherence of low-frequency fixed-path sound in a continental shelf region with internal-wave activity［J］.J Acoust Soc Am，2012，131（1）：1782-1797.

［8］ Colosi JA，DudaT F，Lin Y T.Observationsof sound-speed fluctuations on the New Jersey continental shelf in the summer of 2006［J］.J Acount Soc Am，2012，131（1）：1733-1748.

［9］孫超.水下多傳感器陣列信號(hào)處理［M］.西安：西北工業(yè)大學(xué)出版社.2007.171-172.

［10］ Urick R J.PrinciplesofUnderwaterSound［M］.West-port：Peninsula Publishing.1983.33-34.

［11］邵炫，孫超.陸架斜坡海域信號(hào)水平縱向相關(guān)特性實(shí)驗(yàn)研究［J］.中國(guó)科學(xué)，2016，46（9）：1-8.

［12］謝磊，孫超，劉雄厚，等.陸架斜坡海域聲場(chǎng)特性對(duì)常規(guī)波束形成陣增益的影響［J］.物理學(xué)報(bào)，2016，65（14）：1-11.

Arraygain Analysis and Evaluation of Large-aperture Linear Array Subarray in Shallow Water

LIU Yulu1LIU Qingyu1LIU Xionghou1，2YU Yun1
（1.Naval Academy ，Beijing 100161）（2.Institute of Acoustics，Chinese Academy of Science，Beijing 100190）

When the acoustic signal is propagated in the ocean，the waveform will be distorted due to the complex marine environment.Ideally，the array gain of the conventional beamforming can reach 10 logN，but due to the distortion of the waveform，the signal waveform correlation is reduced at the linear array，resulting in a decrease for the array gain.In this paper，for a typical shallow sea environment，to study singal horizontal correlation when narrowband signal through the propagation arrive at a large aperture linear array，the signal relevant radius can be obtained.Taking the conventional beamforming as an example，the array gain of the linear array is calculated and the factors influencing the gain are analyzed.

signal spatial correlation，conventional beamforming，array gain

TB56

10.3969∕j.issn.1672-9730.2017.10.028

Class Number TB56

2017年5月15日，

2017年6月20日

劉豫魯，男，碩士研究生，研究方向：水聲工程。