舰船科学技术  2016, Vol. 38 Issue (s1): 166-171   PDF    
基于Bellhop声场模型的多阵元信道匹配定位技术研究
赵宏伟, 凌青, 余赟     
海军装备研究院, 北京 100161
摘要: 基于Bellhop声场模型,围绕信道匹配定位的相关问题,对多阵元信道匹配的定位效果进行仿真。在此基础上,对几种失配情况下的定位效果进行讨论。结果表明,垂直阵信道匹配定位效果优于水平阵,海深失配、阵列深度失配、声速失配对信道匹配定位效果影响明显。
关键词: 多阵元     信道匹配     被动定位     失配分析    
Research on multi-sensor channel matched localization technology based on Bellhop acoustic pressure field model
ZHAO Hong-wei, LING Qing, YU Yun     
Naval Academy of Armament, Beijing 100161, China
Abstract: Based on the Bellhop acoustic pressure field model and the related question of channel matched localization, the localization effect of multi-sensor channel matched technology is studied. It has been seen that the localization effect of vertical array is better than horizontal array, it also showed that environmental mismatch such as sea deep, array deep and sound speed distrabution can effectively damage the effect of channel matched localization.
Key words: multi-sensor     channel match     passive localization     mismatch analyze    
0 引言

被动定位一直是水声定位研究领域致力于解决的难题,作为被动定位中的一项重要技术,信道匹配定位技术基于匹配场处理[1]的思想,在水声信道缓慢变化时[2],利用海洋环境参数和声信道传播特性,无需信道详细先验知识,即可消除介质不均匀和边界对探测的影响,实现声信号源位置点的时间、空间聚焦[3]。Bellhop克服了传统射线模型中声影区强度为0和焦散线截面积为0处声强度为无穷大的缺陷[4],通过高斯波束跟踪方法,计算水平非均匀环境中的声场,可以作为信道匹配定位的理想模型。以往的研究中[5-7],单阵元和两阵元信道匹配定位技术已经广泛研究并验证了其性能,本文以10元垂直线阵和10元水平线阵为研究对象,对其信道匹配定位性能进行研究分析。

1 垂直阵和水平阵信道匹配定位原理

从信道冲激响应函数出发,可以给出时域中信道匹配的聚焦特性[8]。信道具有多途时延扩展的特性,声信号沿不同路径的声线不同时刻到达接收点,总的接收信号是通过接收点的所有各声线所传送的信号的干涉叠加,产生复杂的空间滤波特性[9]。多途信道的冲激响应函数为:

$ {h_j}\left( t \right) = \sum\limits_{i = 1}^{{N_j}} {{A_{ji}}\delta \left( {t- {\tau _{ji}}} \right)} {\text {。}} $ (1)

式中:Nj 为通过j号阵元接收点的本征声线的数目;Aji τji 分别为第i途径到达接收点的信号幅度及信号延迟。

j号阵元接收到的信号为:

$ {{\rm{p}}_{rj}}\left( t \right) \!=\! p\left( t \right) \otimes {h_j}\left( t \right){\rm{ + Noise}} \!=\!\! \sum\limits_{i = 1}^{{N_j}} {{A_{ji}}\left( {p\left( {t- {\tau _{ji}}} \right)} \right)} {\rm{ + Noise}},\!\! $ (2)

式中:符号“ $ \otimes $”为卷积运算。

由互易原理可知,声源Sj号阵元间与j号阵元到声源S间的信道冲激响应相同。此处将声源Sj号阵元间的信道冲激响应作为测量场向量,将Bellhop模型得到的海区信道估计即j号阵元到声源S间的信道冲激响应作为匹配场向量。声源发出的信息通过海洋信道并迭加噪声后即为阵元接收信号,将阵元接收信号反转后与j号阵元到声源S间的信道冲激响应进行卷积,等效于将阵元接收信号反向经过模拟信道发回声源SS处接收的信号为:

$ {{\rm{r}}_j}\left( t \right) = {p_{rj}}\left( {- t} \right) \otimes {h_j}{\rm{'}}\left( t \right) = p\left( {- t} \right) \otimes {h_j}\left( {- t} \right) \otimes {h_j}{\rm{'}}\left( t \right), $ (3)

式中:rj t)为声源S处聚焦的信号。这一反转卷积的过程可以在计算机中实现,通过Bellhop模型得到海区的信道估计,并将prj (-t)通过虚拟信道hj ′(t),即实现了信道的匹配。

则在原声源处总的接收信号为:

$ {r}\left( t \right) = \sum\limits_{j = 1}^J {{r_j}\left( t \right)} = p\left( {- t} \right) \otimes \sum\limits_{j = 1}^J {{h_j}} \left( {- t} \right) \otimes {h_j}'\left( t \right){\text {,}} $ (4)

$ {P_T}\left( t \right) = \sum\limits_{j = 1}^J {{h_j}\left( {- t} \right) \otimes {h_j}{\rm{'}}\left( t \right)}{\text {。}} $ (5)

PT为声源S到各阵元之间信道冲激响应函数的自相关函数之和,当hj ′(t)建模误差可以忽略时可将PT近似为δ函数,具有相关峰(t=T时刻)和较低的旁瓣。

图 1展示了阵列信道匹配定位的原理框图,用垂直和水平布放的水听器a~j被动接收目标声源Sr0h0)辐射的噪声,目标声源和水听器之间的海洋声信道用ha~hj表示,a~j号水听器接收到的目标辐射噪声和本地干扰噪声迭加后,将各阵元信号分别与模拟信道haik,……,hjik 进行卷积,haik为水听器a与搜索的空间位置点(rihk )之间的信道,其他同理。

将各阵元接收信号与模拟信道进行卷积的这一过程,即实现了信道匹配,阵元接收信号通过与之相匹配的信道时信号能量同相叠加达到最大,通过不相匹配的信道时信号能量衰减或抵消。水听器信道匹配处理后对应空间位置点(rihk )的输出分别为zaik),……,zjik),将zaik)……zjik)进行叠加的结果为Rzik),Rzik)的峰值Eik 是能量关于空间坐标的函数。由信道匹配定位的原理可知,Eik 的最大值对应的空间位置为声源S位置的估计。

图 1 阵列信道匹配定位原理框图 Fig. 1 The flow chart of channel matched localization theory
2 仿真试验

仿真条件:海深200 m,阵元间距d=10 m,10元垂直阵布放深度20~110 m,10元水平阵布放深度80 m,目标深度80 m,目标距垂直阵水平距离5 700 m,目标距水平阵水平距离5 700~5 790 m,信噪比均为阵元处测量值。海底底质为粗砂、细砂等,底质相关参数为:粗砂密度ρ=2.034 g/cm3,声速c=1.201 ×c0,吸收系数α=0.864 2;细砂密度ρ=1.957 g/cm3,声速c=1.147 ×c0,吸收系数α=0.895 5;砂-粉砂-粘土密度ρ=1.583 g/cm3,声速c=1.033 ×c0,吸收系数α=0.174 0,c0代表海底声速。

通过对图 2图 3的观察分析可得,相同滤波频带,无论对于垂直阵还是水平阵,随着信噪比增大,目标处聚焦的信号能量峰值变大,旁瓣降低,目标点越来越明显,信道匹配定位效果会变好。同一信噪比下,垂直阵的信道匹配定位效果要远远优于水平阵。

图 2 滤波频带[800 1 000] Hz,信噪比-5、0 dB,垂直阵匹配定位效果 Fig. 2 Vertical array output with-5, 0 dB input SNR and [800 1 000] Hz smoothing

图 3 滤波频带[800 1 000] Hz,信噪比为-5、0 dB,水平阵匹配定位效果 Fig. 3 Horizontal array output with-5, 0 dB input SNR and [800 1 000] Hz smoothing

观察图 4图 6发现,垂直阵和水平阵在高信噪比时较高滤波频带[1 800 2 000] Hz下,信道匹配定位效果要好于较低频带[200 400] Hz,目标点更加突出,旁瓣更低。图 5图 7的对比可以看出,低信噪比时,较低滤波频段下匹配定位效果更好。

图 4 信噪比5 dB,滤波频带[200 400]、[1 800 2 000] Hz,垂直阵匹配定位效果 Fig. 4 Vertical array output with 5 dB input SNR and [200 400], [1 800 2 000] Hz smoothing

图 5 信噪比-5 dB,滤波频带[200 400]、[1 800 2 000] Hz,垂直阵匹配定位效果 Fig. 5 Vertical array output with-5 dB input SNR and [200 400], [1 800 2 000] Hz smoothing

图 6 信噪比5 dB,滤波频带[200 400],[1 800 2 000] Hz,水平阵匹配定位效果 Fig. 6 Horizontal array output with 5 dB input SNR and [200 400], [1 800 2 000] Hz smoothing

图 7 信噪比-5 dB,滤波频带[200 400],[1 800 2 000] Hz,水平阵匹配定位效果 Fig. 7 Horizontal array output with-5 dB input SNR and [200 400], [1 800 2 000] Hz smoothing

观察图 8可得,同一信噪比下,垂直阵不同滤波带宽的信道匹配定位效果差别不大,目标处聚焦的信号能量峰值会随着滤波带宽增加而提高。

图 8 信噪比-5 dB,滤波频带[800 1 000],[600 1 200],[400 1 400] Hz,垂直阵匹配定位效果 Fig. 8 Vertical array output with-5 dB input SNR and [800 1 000], [600 1 200], [400 1 400] Hz smoothing

图 9为声速失配时垂直阵匹配定位效果,分别为声速相差1 m/s,3 m/s,6 m/s,可以看出,信道匹配定位的效果对声速比较敏感,声速相差1 m/s,对信道匹配定位效果影响不大,但是随着声速差增加,匹配定位效果下降明显。

图 9 声速失配时,垂直阵信道匹配定位效果 Fig. 9 Vertical array output with sound speed mismatched

图 10为底质失配时,垂直阵信道匹配定位效果,图 10(a)为粗砂底质匹配细砂底质、图 10(b)为粗砂底质匹配砂-粉砂-粘土底质、图 10(c)为粗砂50 m厚底质匹配粗砂100 m厚底质。可以看出,不同的底质失配,对信道匹配定位效果影响不大。

图 10 底质失配时,垂直阵信道匹配定位效果 Fig. 10 Vertical array output with bottom characteristics mismatched

图 11为海深失配对信道匹配定位效果的影响,图 11(a)为实际200 m海深匹配210 m海深、图 11(b)为实际210 m海深匹配200 m海深,目标实际位于(5 700,80)m。可以看出,信道匹配定位对海深非常敏感,海深失配时,信道匹配定位的效果会显著下降。

图 11 海深失配时,垂直阵信道匹配定位效果 Fig. 11 Vertical array output with bathymetry mismatched

图 12为阵列深度失配时,信道匹配定位的效果。图 12(a)为阵列实际0~90 m深度匹配10~100 m深度和阵列实际20~110 m深度匹配10~100 m深度,可以看出信道匹配定位对阵列深度失配也非常敏感,定位效果下降明显。

图 12 阵列深度失配时,垂直阵信道匹配定位效果 Fig. 12 Vertical array output with array deep mismatched
3 结语

通过对仿真试验的结果进行分析,可以得出基于Bellhop声场模型的多阵元信道匹配定位的规律。垂直阵的信道匹配定位效果要远远优于水平阵;无论对于垂直阵还是水平阵,随着信噪比变大,目标处聚焦的信号能量峰值会变大,旁瓣会降低,信道匹配定位效果变好;高信噪比时,较高滤波频段下信道匹配定位效果要优于低频段,较低信噪比时,较低滤波频段下信道匹配定位效果要优于高频段;滤波带宽不同对信道匹配定位效果的影响不明显;信道匹配定位对环境失配等因素非常敏感,其中阵列深度失配、海深失配和声速失配对定位效果影响最显著。

选择垂直阵进行信道的匹配定位,能得到比水平阵更好的效果,选择合适的信噪比和滤波频段,能改善阵列进行信道匹配定位的能力,建模时减小环境失配误差,能获得更好更准确的定位效果。

参考文献
[1] 杨坤德.水声信号的匹配场处理技术研究[D].西安:西北工业大学, 2003.
YANG Kun-de. A study on matched field processing of underwater acoustic signals[D]. Xi'an:Northwestern Polytechnical University, 2003.
[2] 生雪莉.矢量反转镜时空滤波技术及其在水声中的应用[D].哈尔滨:哈尔滨工程大学, 2005. http://www.oalib.com/references/17028958
[3] KUPERMAN W A, HODGKISS W S, SONG H C, et al. Phase conjugation in the ocean:Experimental demonstration of an acoustic time-reversal mirror[J]. The Journal of the Acoustical Society of America , 1998, 103 (1) :25–40. DOI: 10.1121/1.423233
[4] 卢晓亭, 张林. 水声传播建模研究现状综述[J]. 海洋技术 , 2010, 29 (4):48–53.
LU Xiao-ting, ZHANG Lin. Review of sound propagation modeling in underwater warfare environment[J]. Ocean Technology , 2010, 29 (4) :48–53.
[5] 生雪莉, 鲍习中, 方华, 等.基于垂直二元阵时间反转被动定位技术研究[C]//2008年全国声学学术会议论文集.上海:中国声学学会, 2008:2.
SHENG Xue-li, BAO Xi-zhong, FANG Hua, et al. Study on passive orientation of two-sensor-vertical-array based on time-reversal-mirror[C]//Acoustic Academic Conference of 2008. Shanghai:Chinese Acoustical Society, 2008:2.
[6] 惠俊英, 马敬广, 李峰, 等. 二元阵被动时间反转镜定位技术研究[J]. 哈尔滨工程大学学报 , 2007, 28 (11):1247–1251.
HUI Jun-ying, MA Jing-guang, LI Feng, et al. Passive localization technology using time reversal with a two-sensor array[J]. Journal of Harbin Engineering University , 2007, 28 (11) :1247–1251.
[7] 刘家亮, 王海燕, 姜喆, 等. 垂直线列阵结构对PTRM阵处理空间增益的影响[J]. 鱼雷技术 , 2010, 18 (4):263–267.
LIU Jia-liang, WANG Hai-yan, JIANG Zhe, et al. Influence of ULA structure on PTRM array spatial gain[J]. Torpedo Technology , 2010, 18 (4) :263–267.
[8] 殷敬伟.多途信道中Pattern时延差编码水声通信研究[D].哈尔滨:哈尔滨工程大学, 2007. http://www.oalib.com/references/18942081
[9] 余赟.浅海多途信道中声聚焦与声屏蔽技术研究[D].哈尔滨:哈尔滨工程大学, 2009.
YU Yun. Researches on acoustic focusing and acoustic covering in multi-path channel of shallow water[D]. Harbin:Harbin Engineering University, 2009. http://www.cqvip.com/QK/92273A/200903/29899409.html