浅海高频水声信道统计特性分析与仿真

丁 浩，李 伟，邓 鹏，崔沁青

(海军潜艇学院，山东 青岛 266042)

0 引 言

目前人类在海洋中的活动大都集中在浅海，而水声信道由于具备良好的传播特性使其成为海洋研究的主要信息载体，浅海水声信道也成为人们研究的一个热点问题[1]。浅海水声信道具有衰减严重、多径效应和频散效应等特性，会引起传输信号幅度的起伏和相位的波动[2]。

目前，国内外许多学者针对浅海水声信道开展了研究工作，取得了一系列的成果。如文献[3]利用AR(p)模型对浅水时变多途信道模型进行修正；文献[4]在建立系统等效模型的基础上采用用多项式拟合法得到Rice衰落信道一定范围内的信噪比近似解；文献[5 - 6]采用BELLHOP模型对水下信道进行仿真，并研究了对目标声源的探测方法和阵列最佳布放方法与声线的关系；文献[7]应用Rice衰落模型对浅海水声信道进行研究；文献[8 - 10]对水声信道的传播特性进行了分析研究。

掌握高频水声信道的传播特性在军用和民用领域都有十分重要的实用价值，本文基于虚源法和Xiao模型建立高频水声信道模型，并对其统计特性和传播特性进行仿真分析。

1 浅海高频水声信道模型

1.1 虚源法模型

虚源法就是把每一根声线等效为一个对应的虚源所射出的到达接收点的直达声线，到达接收点的信号即为各虚源射出的直达声线的总和。虚源法的信号传播示意图如图1所示，假设海深为H，海水为均匀介质层，上边界为自由界面，下边界为平整海底，发射器位于点O，接收接收器位于点S，R01为直达声线，R02为折射—海底反射声线，R03为折射—海面反射声线，R04为折射—海面反射—海底反射声线，θ为声线入射角。

图1 虚源法信号传播示意图Fig. 1 The signal propagation sketch map of image source method

以平面波为例进行研究，根据线性叠加原理，归一化的接收信号声压可表示为：

其中各本征声线的衰减系数为：

各本征声线相位为：

式中：k为波数，f0为 信号频率；c为传播速度；n为反射次数；Rni为本征声线的路径长度，可表示为：

式中：R0为收发点的水平距离。

海面反射系数V1可表示为：

对于均匀海底，其反射系数V2ni可用瑞利反射系数表示：

将发射器到接收器的几条不同信号传输路径称为本征路径。信号在本征路径上传输时，可看作一个稳定的主分量和许多随机的分散分量之和。其中主分量称为本征分量（本征声线），随机的分散分量称为多径分量。实验证明浅海水声信道传输近距离服从Rice模型，因此，每一条本征路径信道的传播特性可通过Rice衰落仿真体现。

1.2 Xiao模型

Xiao模型是一种用正弦波叠加法实现的Rice信道仿真模型，其优越之处在于它对所有正弦波的路径增益、多普勒频移、初始相位均引入随机变量，尤其是引入随机直射分量，改善了统计特性，能真实反映实际信道的物理特性。其模型如下：

1.3 浅海高频水声信道模型

为体现浅海高频水声信道的传播特性，这里用本征声线来代替Xiao模型中的直射分量，并考虑信道的衰减作用，则接收信号可以表示为：其中： ζm(t) 为第m条 本征路径信道复传播系数；x(t)为发射信号； τm为 传播时延；M为本征路径数量。则有：

式中：Cm， θm， φm分别为第m条本征声线的衰减系数、入射角和相位，其他参数含义同前。

2 水声信道模型统计特性分析

理想Rice衰落过程包络分布的概率密度函数为：

其中：I0(·) 为 第一类零阶修正贝塞尔函数；u为主分量的幅度；2δ2为随机分散分量的平均功率。令 Ω为ζm(t)的平均功率，则可以推导出：

任选一条本征路径，分析Rice因子K和谐波函数个数N对 ζm(t) 概率密度分布的影响。令N=10，当K取不同值时，利用式（11）仿真计算与式（14）计算所得概率分布曲线如图2所示。

图2K取不同值时 ζ m(t)的概率分布曲线Fig. 2 The probability distribution curve of ζm(t) under different K

对图2分析可知，当Rice因子K≤3时， ζm(t)的概率密度分布仿真曲线与理论计算曲线高度吻合，随K值的增大，仿真曲线与理论计算曲线的差距呈增大趋势。

令K=3 ， 当谐波函数个数N取不同值时，利用本文模型仿真与理论模型计算所得概率分布曲线如图3和图4所示。

图3N=5时， ζ m(t)的概率分布曲线Fig. 3 The probability distribution curve of ζm(t) whenN=5

图4N=10时， ζ m(t)的概率分布曲线Fig. 4 The probability distribution curve of ζm(t) when N=10

分析图3和图4可知，谐波函数个数N取不同值时，ζm(t)的概率密度分布仿真曲线与理论计算曲线都高度吻合。

3 浅海高频水声信道仿真

3.1 仿真参数

设发射器与接收器的水平距离R0=2500m，发射器深度z0=10m，接收器深度z=50m，海深H=100m，海水温度T=20 ℃，盐度S=35 ‰，静水压P=16，收发器相对速度v=10m/s，风速w=3kn，c=1.3，ρ=1.0×103kg/m3，c=1.6， ρ =1.8×103kg/m3。根据信道统计特性分析结果，令谐波函数个数N=10，Rice因子K=3。设发射信号为：

其中发射信号幅度A0=10，频率f0=20kHz，f1=19.98kHz，f2=20.02kHz，初始相位 θ0=π/3，则发射信号波形如图5所示。

3.2 传播损失

在理论上，本征声线的数量是无穷多的，随着反射次数n的增大，本征声线对声场的影响越来越小，因此本文只将与直达波声压之比大于0.1的本征声线纳入计算。

图5 发射信号波形图Fig. 5 The waveform image of transmitted signal

接收信号的传播损失可由射线理论计算得到，即

其中p0为距离声源1 m处的参考声压。接收信号的传播损失及其与球面扩展损失曲线比较如图6所示。

图6 传播损失示意图Fig. 6 The sketch map of transmission loss

由图6分析可知，由本文方法计算所得传播损失曲线沿球面扩展损失曲线上下震荡，这是由于本文设定的仿真条件与球面扩展的传播条件相吻合。

3.3 接收信号分析

由式（10）进行仿真计算可得接收信号的波形如图7所示。

由图7可知，与发射信号相比，接收信号的波形明显发生变化，这是由于经过信道传播后，信号产生衰减、相移以及频散的结果。此时，纳入计算的本征声线的数目为6，各本征声线与直达声线的幅度比值如表1所示。

图7 接收信号波形图Fig. 7 The waveform image of received signal

表1 各本征声线归一化幅度衰减值Tab. 1 The normalization amplitude loss value of each eigen acoustic ray

对表1进行分析可知，从3号本征声线开始，其幅度就出现明显下降，经计算可知，海面反射系数V1的 幅值仅为 0 .32，因此发生海面反射时，本征声线的幅度就会明显减小。

4 结 语

本文基于虚源法和Xiao模型建立高频水声信道模型，并分析其统计特性，当Rice因子K≤3时,与理想Rice衰落过程高度吻合，验证了模型的准确性。传播衰减符合球面扩展传播衰减规律，海面（海底）反射系数对本征声线幅度衰减有明显影响。由于信道的衰减、相移以及频散等作用，接收信号波形与发射信号相比发生明显变化。从而实现了对均匀浅海高频水声信道的仿真模拟，能够对高频水声信号的传播损失以及接收波形进行仿真预报。