一种基于双随机基阵被动定位的大桥防撞方法
【专利摘要】本发明提供了一种基于双随机基阵被动定位的大桥防撞方法，该方法基于大桥防撞装置实现，所述装置包括两个随机换能器基阵，分别安装在大桥两侧的桥墩上，用于接收通航船舶的辐射噪声信号；所述方法包括：两个随机换能器基阵对接收到的通航船舶的辐射噪声信号进行预处理、频谱分析和带通滤波，得到船舶辐射噪声的有用信号；采用宽带频域非相干处理方法对船舶辐射噪声的有用信号进行角度估计，得到船舶相对于两个随机换能器基阵的角度θ1和θ2；从而计算出船舶的坐标；根据大桥两侧桥墩到船舶的距离判定船舶是否进入大桥防撞危险局域。本发明提出的方法测量精度高，装置安装简单方便，成本较低。
【专利说明】
-种基于双随机基阵被动定位的大桥防撞方法
技术领域
[0001] 本发明设及一种海洋领域的水声信号处理方法，具体设及一种基于双随机基阵被 动定位的大桥防撞方法。
【背景技术】
[0002] 跨河跨海大桥的建设，极大的缓解了路上交通的拥挤。但是大量跨河跨海桥梁的 修建，相当于给水上交通制造了许多障碍物，并且随着船舶向大型化发展和船舶通航密度 逐渐增大，导致通航船舶碰撞桥梁事故频繁发生，对海上交通造成了极大的威胁。现阶段应 用在防撞预警系统的方法主要有:激光感应、雷达定位、视频监控和水声定位等。
[0003] 水声目标定位主要分为被动目标定位和主动目标定位，与主动目标定位相比，被 动定位的隐蔽性较好。常用的被动定位方法有传统的被动定位方法、目标运动定位分析法、 匹配场被动定位方法、多基阵联合的纯方位被动跟踪方法等。但上述方法存在一定的不足， 如传统的被动定位方法，包括=子阵法和球面内插法，但其需要预先知道时延差，存在较大 的误差；目标运动定位分析法，在实时跟踪时应用困难;匹配场被动定位方法，在信道环境 复杂的条件下存在较大的误差;多基阵联合的纯方位被动跟踪方法，其设备较复杂，成本较 局等。
【发明内容】
[0004] 本发明的目的在于克服目前水声目标被动定位方法存在的上述缺陷，提出了一种 基于双随机基阵被动定位的大桥防撞方法，该方法将两条随机换能器基阵安装在大桥两侧 的桥墳上，被动接收靠近大桥的通航船舶的福射噪声，然后利用频域角度估计方法，对福射 噪声信号进行处理，进行船舶角度估计，最后根据双基阵间距W及与船舶的角度信息，融合 估计出船舶与桥墳的距离，并判断是否进入危险区，及时发出预警。该方法能实时检测出船 舶距离桥墳的距离，及时发出预警，检测精度高，并且简单方便，易于安装和维护，成本较 低。
[0005] 为了实现上述目的，本发明提出了一种基于双随机基阵被动定位的大桥防撞方 法，该方法基于大桥防撞装置实现，所述装置包括两个随机换能器基阵，分别安装在大桥两 侧的桥墳上，用于接收通航船舶的福射噪声信号;所述方法包括:两个随机换能器基阵对接 收到的通航船舶的福射噪声信号进行预处理、频谱分析和带通滤波，得到船舶福射噪声的 有用信号;采用宽带频域非相干处理方法对船舶福射噪声的有用信号进行角度估计，得到 船舶相对于两个随机换能器基阵的角度和02;从而计算出船舶的坐标(x〇，y〇);根据大桥 两侧桥墳到船舶的距离判定船舶是否进入大桥防撞危险局域。
[0006] 上述技术方案中，所述方法具体包括：
[0007] 步骤1)两个随机换能器基阵被动接收通航船舶的福射噪声信号，对船舶福射噪声 信号进行功率放大、带通滤波预处理；
[000引步骤2)对预处理之后的船舶福射噪声信号进行频谱分析，得到船舶福射噪声信号 的显著频率分布范围；
[0009] 步骤3)根据船舶福射噪声信号的显著频率分布范围，对船舶福射噪声信号进行带 通滤波，去除其它频率的无用噪声，筛选出船舶福射噪声的有用信号；
[0010] 步骤4)采用宽带频域非相干处理方法，对步骤3)得到的输出信号进行角度估计， 得到船舶相对于两个随机换能器基阵的角度和02;
[ocm]步骤5)根据两个角度值01和02,求解船舶的坐标(xo，yo)。
[0012] 步骤6)计算得到船舶距离两个桥墳的长度：
[0013]
[0014]
[001引其中，（財，71)，（料，72)分别为两个桥墳的坐标；
[0016] 步骤7)设定大桥避撞危险局域半径为D，若ri《D或n《D，则发出警报信息。
[0017] 上述技术方案中，所述步骤4)具体包括：
[0018] 步骤401)对步骤3)输出的数据进行时域数据分段，分为K个数据段，每段数据的采 样点数为J;对每段数据进行J点的FFT变化处理，取K个数据段中相同频率柜的数据组成频 域快拍数据，则宽带信号模型可W用频域表示为：
[0019] X(fj)=A(fj)S(fj)+N(fj) j = J
[0020] 式中，A(fj)是MXN的导向矢量矩阵，S(fj)为接收信号复包络，N(fj)为无用信号噪 声，M为阵元个数，N为信号源个数，船舶福射噪声频段共分为J个频率柜，f J为第j个频率柜；
[0021] 步骤402)根据最大似然估计定理，估计出每个频率柜上的采样数据协方差矩阵：
[0022] R(fj)=X(fj)X^(fj)/K
[0023] 步骤403)对协方差矩阵进行信号到达角度估计处理，估计得到每个频率柜所对应 的角度功率谱：
[0024]
[0025] 其中，a(0)为来波方向为0时，随机换能器基阵所对应的导向矢量；
[0026] 步骤404)对所有频率柜的角度功率谱进行累加求和，得出宽带船舶福射噪声的总 功率谱为：
[0027]
[0028] 步骤405)估计两个随机换能器基阵与船舶之间的角度；
[0029] 船舶福射噪声角度功率谱的尖峰对应着船舶相对于第一个随机换能器基阵的角 度：
[0030]
[0031] 通过上述步骤求解出船舶相对于第二个随机换能器基阵的角度02。
[0032] 本发明的优势在于：
[0033] 1、本发明提出的法可W对通航船舶进行定位，适合于对各类船舶；
[0034] 2、本发明提出的方法测量精度高，安装简单方便，成本较低。
【附图说明】
[0035] 图1为本发明中计算船舶位置的示意图；
[0036] 图2a为实例1的基阵1的功率谱图；
[0037] 图化为实例1的基阵2的功率谱图。
【具体实施方式】
[0038] 下面结合附图和具体实例对本发明做进一步详细的说明。
[0039] 如图1所示，一种基于双随机基阵被动定位的大桥防撞方法；该方法实现的条件 是:在大桥两侧的桥墳上分别安装一个随机换能器基阵;根据实际工程所需，合理设计随机 换能器基阵的参数，在两个阵列阵元个数相同的情况下，与传统的半波长布阵的均匀线阵 相比，随机换能器基阵的角度分辨率远大于均匀阵的角度分辨率;所述方法具体包括：
[0040] 步骤1)两个随机换能器基阵被动接收通航船舶的福射噪声信号，对船舶福射噪声 信号进行功率放大、带通滤波等预处理；
[0041] 步骤2)对预处理之后的船舶福射噪声信号进行频谱分析，得到船舶福射噪声信号 的显著频率分布范围；
[0042] 采用FFT变换，并采用频谱检测方法，检测出船舶福射噪声的频率分布范围；
[0043] 步骤3)根据船舶福射噪声信号的显著频率分布范围，对船舶福射噪声信号进行带 通滤波，去除其它频率的无用噪声，筛选出船舶福射噪声的有用信号；
[0044] 步骤4)采用宽带频域非相干处理方法，对噪声信号进行角度估计;具体包括：
[0045] 步骤401)对步骤3)输出的数据进行时域数据分段，分为K个数据段，每段数据的采 样点数为J;对每段数据进行J点的FFT变化处理，取K个数据段中相同频率柜的数据组成频 域快拍数据，则宽带信号模型可W用频域表示为：
[0046] X(fj)=A(fj)S(fj)+N(fj) j = J
[0047] 式中，A(fj)是MXN的导向矢量矩阵，S(fj)为接收信号复包络，N(fj)为无用信号噪 声，M为阵元个数，N为信号源个数，船舶福射噪声频段共分为J个频率柜，f J为第j个频率柜；
[0048] 步骤402)根据最大似然估计定理，估计出每个频率柜上的采样数据协方差矩阵：
[0049] R(fj)=X(fj)X^(fj)/K
[0050] 步骤403)对协方差矩阵进行信号到达角度估计处理，估计得到每个频率柜所对应 的角度功率谱：
[0化1 ]
[0052] 其中，a(0)为来波方向为0时，随机换能器基阵所对应的导向矢量；
[0053] 步骤404)对所有频率柜的角度功率谱进行累加求和，得出宽带船舶福射噪声的总 功率谱为：
[0化4]
[0055] 步骤405)估计两个随机换能器基阵与船舶之间的角度；
[0056] 船舶福射噪声角度功率谱的尖峰对应着船舶相对于第一个随机换能器基阵的角 度：
[0化7]
[0058] 通过上述步骤求解出船舶相对于第二个随机换能器基阵的角度02。
[0059] 步骤5)根据两个角度值01和02,求解船舶的坐标(xo，yo);
[0060] 如图1所示，实际桥墳的地理位置在位置规划时是已知的，分别为（XI，yi)，（X2， 72);根据；角形几何关系所得，1曰11目1=(7日-71)/(則-义1)，1:曰]1目2=(7日-72)/(則-如，结合；角 关系式，可W求得船舶的坐标(XO，yo);
[0061 ] 巧驢6 )结合桥谢的化理仿置，计貸得到船舶距离两个桥墳的长度分别为
[0062] 步骤7)设定大桥避撞危险局域半径为D，若ri《D或n《D，则会发出警报，使得船舶 远离大桥，避免船舶碰撞大桥事故的发生。
[0063] 实例 1:
[0064] 系统参数为:两个随机换能器基阵的阵元个数为10个;船舶福射噪声的频谱设为 25曲Z-35曲Z，采样率为400曲Z，基阵1的坐标为(2,5)米，基阵2的坐标为（10，1)米。
[0065] 由图2a和图化可知，两个基阵对接收的船舶福射噪声进行处理，得到的角度功率 谱可W检测出，船舶偏离基阵1的角度为40%偏离基阵2的角度为120%并根据基阵1和基阵 2的位置聲标可得：
[0066]
[0067]
[0068] 可计算得出船舶的位置坐标为(8.92,10.8)。
[0069] 进而求得船舶距离基阵1的距离
米，船舶距离基阵2 的距离天
。
[0070] 最后，将ri，n与实际的危险半径进行比较，得出是否需要发出警报。
【主权项】
1. 一种基于双随机基阵被动定位的大桥防撞方法，该方法基于大桥防撞装置实现，所 述装置包括两个随机换能器基阵，分别安装在大桥两侧的桥墩上，用于接收通航船舶的辐 射噪声信号;所述方法包括:两个随机换能器基阵对接收到的通航船舶的辐射噪声信号进 行预处理、频谱分析和带通滤波，得到船舶辐射噪声的有用信号;采用宽带频域非相干处理 方法对船舶辐射噪声的有用信号进行角度估计，得到船舶相对于两个随机换能器基阵的角 度Θ#ΡΘ 2;从而计算出船舶的坐标;根据大桥两侧桥墩到船舶的距离判定船舶是否进入大桥 防撞危险局域。2. 根据权利要求1所述的基于双随机基阵被动定位的大桥防撞方法，其特征在于，所述 方法具体包括： 步骤1)两个随机换能器基阵被动接收通航船舶的辐射噪声信号，对船舶辐射噪声信号 进行功率放大、带通滤波预处理； 步骤2)对预处理之后的船舶辐射噪声信号进行频谱分析，得到船舶辐射噪声信号的频 率分布范围； 步骤3)根据船舶辐射噪声信号的频率分布范围，对船舶辐射噪声信号进行带通滤波， 去除其它频率的无用噪声，筛选出船舶辐射噪声的有用信号； 步骤4)采用宽带频域非相干处理方法，对步骤3)得到的噪声信号进行角度估计，得到 船舶相对于两个随机换能器基阵的角度Θ#ΡΘ2; 步骤5)根据两个角度值01和02，求解船舶的坐标(XQ，yo); 步骤6)计算得到船舶距离两个桥墩的长度：其中，（XI，yi)，（ X2，y2)分别为两个桥墩的坐标； 步骤7)设定大桥避撞危险局域半径为D，若ri$DSr2<D，则发出警报信息。3. 根据权利要求2所述的基于双随机基阵被动定位的大桥防撞方法，其特征在于，所述 步骤4)具体包括： 步骤401)对步骤3)输出的数据进行时域数据分段，分为K个数据段，每段数据的采样点 数为J;对每段数据进行J点的FFT变化处理，取K个数据段中相同频率柜的数据组成频域快 拍数据，则宽带信号模型可以用频域表示为： X(fj)=A(fj)S(fj)+N(fj)j = l,2,··· ,J 式中，A(fj)是MXN的导向矢量矩阵，S(fj)为接收信号复包络，N(fj)为无用信号噪声，Μ 为阵元个数，Ν为信号源个数，船舶辐射噪声频段共分为J个频率柜，t为第j个频率柜； 步骤402)根据最大似然估计定理，估计出每个频率柜上的采样数据协方差矩阵： R(fj)=X(fj)XH(fj)/K 步骤403)对协方差矩阵进行信号到达角度估计处理，估计得到每个频率柜所对应的角 度功率谱：其中，a( Θ)为来波方向为Θ时，随机换能器基阵所对应的导向矢量； 步骤404)对所有频率柜的角度功率谱进行累加求和，得出宽带船舶辐射噪声的总功率 谱为：步骤405)估计两个随机换能器基阵与船舶之间的角度； 船舶辐射噪声角度功率谱的尖峰对应着船舶相对于第一个随机换能器基阵的角度：通过上述步骤求解出船舶相对于第二个随机换能器基阵的角度θ 2。
【文档编号】G01S15/18GK105954755SQ201610371527
【公开日】2016年9月21日
【申请日】2016年5月30日
【发明人】闫路, 许枫, 唐浩, 安旭东, 王梦宾
【申请人】中国科学院声学研究所