CN107651119B - 具有水下噪声源被动定位功能的海上浮标及定位方法 - Google Patents

具有水下噪声源被动定位功能的海上浮标及定位方法 Download PDF

Info

Publication number
CN107651119B
CN107651119B CN201710738570.4A CN201710738570A CN107651119B CN 107651119 B CN107651119 B CN 107651119B CN 201710738570 A CN201710738570 A CN 201710738570A CN 107651119 B CN107651119 B CN 107651119B
Authority
CN
China
Prior art keywords
array
noise source
positioning
array element
supporting rod
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201710738570.4A
Other languages
English (en)
Other versions
CN107651119A (zh
Inventor
卢涛
韩宝坤
魏国
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shandong University of Science and Technology
Original Assignee
Shandong University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shandong University of Science and Technology filed Critical Shandong University of Science and Technology
Priority to CN201710738570.4A priority Critical patent/CN107651119B/zh
Publication of CN107651119A publication Critical patent/CN107651119A/zh
Application granted granted Critical
Publication of CN107651119B publication Critical patent/CN107651119B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B63—SHIPS OR OTHER WATERBORNE VESSELS; RELATED EQUIPMENT
    • B63B—SHIPS OR OTHER WATERBORNE VESSELS; EQUIPMENT FOR SHIPPING 
    • B63B22/00—Buoys
    • B63B22/18—Buoys having means to control attitude or position, e.g. reaction surfaces or tether
    • B63B22/20—Ballast means
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B63—SHIPS OR OTHER WATERBORNE VESSELS; RELATED EQUIPMENT
    • B63B—SHIPS OR OTHER WATERBORNE VESSELS; EQUIPMENT FOR SHIPPING 
    • B63B22/00—Buoys
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/18—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using ultrasonic, sonic, or infrasonic waves
    • G01S5/20—Position of source determined by a plurality of spaced direction-finders
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J7/00—Circuit arrangements for charging or depolarising batteries or for supplying loads from batteries
    • H02J7/34—Parallel operation in networks using both storage and other dc sources, e.g. providing buffering
    • H02J7/35—Parallel operation in networks using both storage and other dc sources, e.g. providing buffering with light sensitive cells
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B63—SHIPS OR OTHER WATERBORNE VESSELS; RELATED EQUIPMENT
    • B63B—SHIPS OR OTHER WATERBORNE VESSELS; EQUIPMENT FOR SHIPPING 
    • B63B22/00—Buoys
    • B63B2022/006—Buoys specially adapted for measuring or watch purposes
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30—Assessment of water resources

Abstract

本发明公开了一种具有水下噪声源被动定位功能的海上浮标，包括浮标本体，在浮标本体的底部连接有锚链，在锚链的底部连接有重力锚，重力锚固定在海底，在浮标本体上还连接有尼龙绳，在尼龙绳上设置有噪声源定位阵元装置，在浮标本体上设置有摄像装置。采用上述浮标进行水下噪声源定位的方法包括以下步骤：步骤一、采用自适应门限检测方法对水下噪声源判别；步骤二、采用多阵元混合定位方法对判别出的水下噪声源进行被动声源定位；步骤三、启动摄像装置，对定位的噪声源进行拍照取证。本发明可用于对靠近浮标的水下噪声源进行定位，并可进行摄像等取证行为，而且功耗低。

Description

具有水下噪声源被动定位功能的海上浮标及定位方法
技术领域
本发明涉及一种海上浮标，具体地说是涉及一种具有水下噪声源被动定位功能的海上浮标，以及基于该海上浮标的定位方法。
背景技术
浮标，是指浮于水面上的一种航标，是指锚定在指定位置，用以标示航道范围，指示浅滩、碍航物或表示专门用途的水面助航标志。现在有的多功能浮标还装载有雷达应答器、无线电指向标、雾警信号和海洋调查仪器等设备，集成功能越来越多。但是随着当前航运量发展、领土争端等多方面原因，海上浮标容易遭到船舶撞击或者遭到蓄意破坏，尤其是搭载有精密调查仪器设备的海洋调查浮标，抗破坏能力更弱。浮标作为海上无人值守的平台，遭到破坏后不易搜集破坏证据。
发明内容
基于上述技术问题，本发明提供一种具有水下噪声源被动定位功能的海上浮标及定位方法。
本发明所采用的技术解决方案是：
一种具有水下噪声源被动定位功能的海上浮标，包括浮标本体，在浮标本体的底部连接有锚链，在锚链的底部连接有重力锚，重力锚固定在海底，在浮标本体上还连接有尼龙绳，在尼龙绳上设置有噪声源定位阵元装置，所述噪声源定位阵元装置包括第一支撑杆、第二支撑杆和第三支撑杆，所述第一支撑杆、第二支撑杆和第三支撑杆均相互交叉垂直布置，共同连接成三维坐标系结构，在第一支撑杆、第二支撑杆和第三支撑杆上分别设置有偶数个声压传感器，且第一支撑杆、第二支撑杆或第三支撑杆上的声压传感器均沿第一支撑杆、第二支撑杆和第三支撑杆的交点呈对称分布。
优选的，在第一支撑杆、第二支撑杆和第三支撑杆上均设置有4个声压传感器，沿交点呈两两对称分布，其中距离交点较远的声压传感器与交点之间的距离为距离交点较近的声压传感器与交点之间距离的两倍。
优选的，在尼龙绳上间隔设置有用于减少流噪声的布条流苏，在尼龙绳的底部连接有金属坠。
优选的，在浮标本体上设置有摄像装置。
一种水下噪声源被动定位方法，采用上述的海上浮标，包括以下步骤：
步骤一、采用自适应门限检测方法对水下噪声源判别
首先对各个声压传感器接收到的信号进行处理，对处理过的信号x(t)进行双门限能量检测，以λ1，λ2作为能量检测的上下门限，若累积能量Y＞λ1，则认为所选频谱被占用；若Y＜λ2，则认为所选频谱空闲；若λ2≤Y≤λ1，则认为利用能量检测法无法判决所选频谱的占用状态，定义该区域为不可靠区S，此时再对处理过的信号x(t)进行循环平稳特征检测；
步骤二、采用多阵元混合定位方法对判别出的水下噪声源进行被动声源定位
(1)利用三维坐标系X、Y、Z三个轴上的以原点为对称点的6m+n个阵元，m≥1,n≥0，且n为偶数，必须在一个轴上，且关于坐标系原点对称确立阵元坐标；
(2)以原点o为中心，对称分布于x,y,z轴上，其中M1-M6为6个阵元，组成XOY、YOZ、XOZ三个平面坐标系上的十字阵，此外有6(m-1)+n个阵元呈十字阵或线列阵布置在其中；
(3)设声源位置为S(x,y,z)，声源距坐标原点距离为r，声源S到阵元M1的距离为R1，声源S到阵元M2的距离为R2，d12,d13,d14为阵元之间的声程差，以d12为例进行说明，d12为目标到阵元1和阵元2的距离之差，则根据目标和声阵的几何位置关系有：
利用M1～M6 6个阵元组成空间坐标系中的三个四元十字阵
分别为xoy面，4个阵元组合M1,M2,M3,M4
分别为yoz面，4个阵元组合M2,M5,M4,M6
分别为xoz面，4个阵元组合M1,M5,M3,M6
A：xoy面十字阵，4个阵元M1,M2,M3,M4
Figure BDA0001388612480000021
得出
Figure BDA0001388612480000022
B：yoz面，4个阵元组合M2,M5,M4,M6
Figure BDA0001388612480000031
得出
Figure BDA0001388612480000032
C:xoz面，4个阵元组合M1,M5,M3,M6
Figure BDA0001388612480000033
得出
Figure BDA0001388612480000034
其中d12＝τ12c，d13＝τ13c，d14＝τ14c，d15＝τ15c，d16＝τ16c，d24＝τ24c，d25＝τ25c，d26＝τ26cτij为声源到阵元i和j的时延差，c为声速，目标声源坐标(x,y,z)；
(4)剩余6(m-1)+n个阵元可任意放在x轴，y轴，z轴上，并与x或y或z轴上的阵元组成间距为2L的线列阵或者间距为4L的十字阵；
以线列阵为例，该阵元位置为Mi(xi,yi,zi)，则有
Figure BDA0001388612480000035
c为声速；ti为第i个阵元所测时延；
以x轴三元线列阵为例
Figure BDA0001388612480000036
得到
Figure BDA0001388612480000041
得出
Figure BDA0001388612480000042
y存在二值：即线列声压阵存在的左右舷模糊问题；
(5)利用各十字阵和线列阵通过计算时延差得到值做求平均处理，得到一个较为准确的目标声源位置。
优选的，所述步骤一：在能量检测法中，判决门限λ是依据期望检测概率
Figure BDA0001388612480000043
的要求设置的，但由于信道环境的不确定性，当接收信号能量值在门限λ附近时很难做出可靠判断，为此，定义一个范围
Figure BDA0001388612480000044
为判决的不可靠区，即λ1＝kλ；λ2＝λ/k，式子中k≥1为不可靠区S的范围调整因子，且有Δ＝10logk；
通过检测信道保护带宽内的噪声能量可以近似估计该信道内的信噪比，从而得到信道状态，为了适应信道变化，依据估计的信道信噪比信息动态设置门限；如果实际检测概率比期望检测概率低，就适当增大λ1和λ2之间的距离，即增大k值，如果实际检测概率比期望检测概率高，就减小k值；当信道质量足够好时，令k＝1，此时λ1＝λ2＝λ，实行单门限检测。
优选的，所述步骤一：为了避免过于频繁的调整门限，引入期望检测概率范围因子a，当实际检测概率高于
Figure BDA0001388612480000045
低于
Figure BDA0001388612480000046
时才调整门限，若检测概率处于二者之间则不必调整，这里
Figure BDA0001388612480000047
为期望检测概率，a的取值应远小于1，以保证尽可能达到期望的检测概率。
优选的，所述阵元为声压阵元。
优选的，上述定位方法还包括：
步骤三、启动摄像装置，对定位的噪声源进行拍照取证。
本发明的有益技术效果是：
(1)浮标在海上使用时，一般通过蓄电池或者太阳能电池板提供电源，电源具有有限且不稳定的特性，而且浮标上搭载的各种设备需要长期工作，功耗较大，因此本发明在不改变现有浮标结构和功能的条件下，设计了一种可用于对靠近浮标的水下噪声源进行定位、且功耗低的浮标；其中噪声源定位采用声压被动定位方式(门限检测方式为自适应双门限频谱感知方式)，当水下噪声源累积功率超过自适应门限检测时，通过多声压阵元(声压传感器)对水下噪声源进行定位；定位后启动浮标上的摄像装置或者其它记录设备对靠近浮标的噪声源进行摄像等取证行为。
(2)本发明在浮标入水部分的下端通过软尼龙绳悬挂多阵元，阵元下方通过软尼龙绳挂着重物，尽可能的保证多阵元在波浪起伏或者恶劣海况条件下导致的浮标晃动的情况下多阵元稳定性。软尼龙绳上随机绑有布条流苏，其目的在于减少水流冲击尼龙绳造成的流噪声。
附图说明
下面结合附图与具体实施方式对本发明作进一步说明：
图1为本发明具有水下噪声源被动定位功能的海上浮标的结构示意图；
图2为在三维坐标系上确立阵元坐标的示意图。
具体实施方式
如图1所示，一种具有水下噪声源被动定位功能的海上浮标，包括浮标本体1，在浮标本体1的底部连接有锚链2，在锚链2的底部连接有重力锚3，重力锚3固定在海底。在浮标本体1上还连接有尼龙绳4，在尼龙绳4上设置有噪声源定位阵元装置5，所述噪声源定位阵元装置5包括第一支撑杆501、第二支撑杆502和第三支撑杆503，所述第一支撑杆501、第二支撑杆502和第三支撑杆503均相互交叉垂直布置，共同连接成XYZ三维直角坐标系结构。在第一支撑杆501、第二支撑杆502和第三支撑杆503上分别设置有偶数个声压传感器504，且第一支撑杆501、第二支撑杆502或第三支撑杆503上的声压传感器504均沿第一支撑杆、第二支撑杆和第三支撑杆的交点呈对称分布。
作为对本发明的进一步设计，在第一支撑杆501、第二支撑杆502和第三支撑杆503上均设置有4个声压传感器，沿交点呈两两对称分布，其中距离交点较远的声压传感器与交点之间的距离为距离交点较近的声压传感器与交点之间距离的两倍。声压传感器采用该布置方式能够较为精确的进行噪声源定位。
进一步的，在尼龙绳4上间隔设置有用于减少流噪声的布条流苏6，在尼龙绳的底部连接有金属坠7，以尽可能的保证多阵元在波浪起伏或者恶劣海况条件下导致的浮标晃动的情况下多阵元稳定性。
更进一步的，在浮标本体1上设置有摄像装置8，以在对噪声源进行定位后启动摄像功能进行取证。
利用上述海上浮标进行水下噪声源被动定位的方法，包括以下步骤：
步骤一、采用自适应门限检测方法对水下噪声源判别
因为海上浮标多采用电池或者太阳能电池板供电，电量有限且不稳定，因此不能采用实时探测的方法进行探测。此外，海面和海水中的噪声强度变化快，起伏大，因此在检测过程中设置检测门限时并未采用固定门限值来设置门限，而是采用自适应双门限频谱感知算法，具体介绍如下：
首先对各个阵元接收到的信号进行算数平均滤波得到x(t)，对处理过的信号x(t)进行双门限能量检测，以λ1，λ2作为能量检测的上下门限，若累积能量Y＞λ1，则认为所选频谱被占用；若Y＜λ2，则认为所选频谱空闲；若λ2≤Y≤λ1，则认为利用能量检测(固定门限方法)无法判决所选频谱的占用状态，定义该区域为不可靠区S，此时再对处理过的信号x(t)进行循环平稳特征检测。
自适应门限设置：在能量检测法(固定门限方法)中，判决门限λ是依据期望检测概率
Figure BDA0001388612480000061
的要求设置的，但由于信道环境的不确定性，当接收信号能量值在门限λ附近时很难做出可靠判断，为此，定义一个范围
Figure BDA0001388612480000062
为判决的不可靠区，即λ1＝kλ；λ2＝λ/k，式子中k≥1为不可靠区S的范围调整因子，且有Δ＝10logk。
通过检测信道保护带宽内的噪声能量可以近似估计该信道内的信噪比，从而得到信道状态，为了适应信道变化，依据估计的信道信噪比信息动态设置门限；如果实际检测概率比期望检测概率低，就适当增大λ1和λ2之间的距离，即增大k值，如果实际检测概率比期望检测概率高，就减小k值。当信道质量足够好时，令k＝1，此时λ1＝λ2＝λ，实行单门限检测。通过这种门限设置方法在保证检测概率的同时，又使得在信道质量好的时候降低算法复杂度，减少系统功耗。
为了避免过于频繁的调整门限，引入期望检测概率范围因子a，当实际检测概率高于
Figure BDA0001388612480000063
低于
Figure BDA0001388612480000064
时才调整门限，若检测概率处于二者之间则不必调整，这里
Figure BDA0001388612480000065
为期望检测概率，a的取值应远小于1，以保证尽可能达到期望的检测概率。
步骤二、采用多阵元混合定位方法对判别出的水下噪声源进行被动声源定位
传统的使用声压阵元进行声源定位的方法主要有线列阵、十字阵、圆阵等阵型，这些阵型在使用过程中都存在一些不可避免的问题，比如线列阵定位方法中存在声源目标左右舷不分的问题，十字阵在估计距离时有较大误差，其次在估计仰角时需要声源和多阵元之间的仰角较大时才有较高精度。因此本发明提出了一种利用直角坐标系轴上的多阵元进行被动声源定位的方法，利用直角坐标系X、Y、Z三个轴上的以原点为对称点的6m+n(m≥1,n≥0，且n为偶数，必须在一个轴上，且关于坐标系原点对称)个阵元，如图2所示。
则在空间坐标系中有阵元坐标为M1(L,0,0)，M2(0,L,0)，M3(-L,0,0)，M4(0,-L,0)，M5(0,0,L)，M6(0,0,-L)，M7(3L,0,0),M8(0,3L,0),M9(-3L,0,0),M10(0,-3L,0),M11(0,0,3L),M12(0,0,-3L)......以原点o为中心，对称分布于x,y,z轴上，其中M1-M6为6个阵元(最少为6个阵元)，组成XOY、YOZ、XOZ三个平面坐标系上的十字阵，此外有6(m-1)+n个阵元呈十字阵或线列阵布置在其中。
设声源位置为S(x,y,z)，声源距坐标原点距离为r。声源S到阵元M1的距离为R1，声源S到阵元M2的距离为R2，d12,d13,d14等为阵元之间的声程差，以d12为例进行说明，d12为目标到阵元1和阵元2的距离之差，则根据目标和声阵的几何位置关系有：
利用M1～M66个阵元组成空间坐标系中的三个四元十字阵
分别为xoy面，4个阵元组合M1,M2,M3,M4
分别为yoz面，4个阵元组合M2,M5,M4,M6
分别为xoz面，4个阵元组合M1,M5,M3,M6
A：xoy面十字阵。4个阵元M1,M2,M3,M4
Figure BDA0001388612480000071
得出
Figure BDA0001388612480000072
B：yoz面，4个阵元组合M2,M5,M4,M6
Figure BDA0001388612480000073
得出
Figure BDA0001388612480000074
C:xoz面，4个阵元组合M1,M5,M3,M6
Figure BDA0001388612480000081
得出
Figure BDA0001388612480000082
其中d12＝τ12c，d13＝τ13c，d14＝τ14c，d15＝τ15c，d16＝τ16c，d24＝τ24c，d25＝τ25c，d26＝τ26c
τij为声源到阵元i和j的时延差，c为声速，目标声源坐标(x,y,z)；
剩余6(m-1)+n个阵元可任意放在x轴，y轴，z轴上，并与x或y或z轴上的阵元组成间距为2L的线列阵(无法满足十字阵的情况下)或者间距为4L(随着阵元数目增加，可以更大，但是要小于半波长)的十字阵(线列阵方法如下，十字阵方法同上)。
以线列阵为例，该阵元位置为Mi(xi,yi,zi)，则有
Figure BDA0001388612480000083
c为声速；ti为第i个阵元所测时延；
以x轴三元线列阵为例
Figure BDA0001388612480000084
得到
Figure BDA0001388612480000085
得出
Figure BDA0001388612480000086
y存在二值：即线列声压阵存在的左右舷模糊问题。
最后利用各十字阵和线列阵通过计算时延差得到值做求平均处理，得到一个较为准确的目标声源位置。
上述方法中，最少需要六个阵元组成三个十字阵进行测量定位，以图中所示为例，12个阵元，组成共计有间距为2L的十字阵三个，间距为4L的十字阵三个，间距为2L的四元线列阵三条，每一个阵型都是进行单独时延定位，最后对每一个阵型的处理结果加权求平均，以有效解决传统水下声源被动定位的缺点。
以上阵元为声压阵元，具体为声压传感器，空气、水中均可使用，要求阵元间距小于半波长＝(声速/频率)/2。
上述方法中，还包括步骤三：启动摄像装置，对定位的噪声源进行拍照取证。
下面对本发明方法作如下补充说明：
1、声功率与声压之间的关系，用于门限检测，即根据能量进行上下门限检测
W＝(P^2*S)/ρc，其中P为声压，即声压传感器测得的值，单位为pa；S为声波垂直通过声压传感器的面积，即声压传感器的面积，单位为m2，这个值已知，属于对硬件设备的测量；ρc为声传播媒质的特性阻抗，单位为瑞利，即帕*秒/米(pa*s/m)，这个值也是已知，可以在需要使用的区域提前实际测量确定该值大小，或者使用默认的海水的特性阻抗值。
2、检测概率和虚警概率之间的关系(这里把涉及两者的所有因素给出了全面解释，包括判决概率，检测概率，虚警概率，漏报概率)；
噪声源实际目标有两种情况，有或者无；检测判决也有两种情况，有或者无；则共有四种情况，即有目标，可能检测为有，也可能检测为无，没有目标，可能检测为有，也可能检测为无。
(1)实际是H0假设为真，而判决检测为H0为真；
(2)实际是H0假设为真，而判决检测为H1为真；
(3)实际是H1假设为真，而判决检测为H0为真；
(4)实际是H1假设为真，而判决检测为H1为真；
则正确的判决检测为(1)(4)，错误的判决检测为(2)(3)。
对应于每一种判决检测结果，有相应的判决检测概率P(Dj|Hi)(i,j＝0,1)，假设Hi为真的条件下，判决检测Hj成立的概率；在假设Hi为真的条件下，判决检测量(x|Hi)的概率密度函数为f(x|Hi)，由于判决检测量(x|Hi)落在判决空间Di，则判决Hi成立，所以判决概率有
Figure BDA0001388612480000091
其中i,j＝0,1，对于判决概率，我们希望正确的判决概率尽可能大，而错误判决概率尽可能小。(判决概率是评价检测性能的重要因素之一)
代价函数Cij:表示实际是Hj假设为真，而判决为Hi假设为真所付出的代价。也称为风险函数。
正确的判决无代价：C00＝C11＝0；
检测概率：正确判决的概率P(D1|H1)和P(D0|H0)；实际有噪声源目标，检测为有噪声源目标，或者，实际没有噪声源目标，检测为没有噪声源目标的概率。
虚警：实际H0假设为真，而判决为H1假设为真(实际没有噪声源目标，检测为有噪声源目标)。又称为第一类错误。虚警引入的代价称为虚警代价C01。虚警发生的概率(虚警概率)为：P(D1|H0)称为虚警概率。
漏报：实际H1假设为真，而判决为H0假设为真。称为第二类错误。漏报引入的代价称为漏报代价C10。漏报发生的概率为：P(D0|H1)称为漏报概率。
3、信号异常检测算法
这里异常与否是通过比较概率大小得到；
传统的频谱感知算法，多利用接收信号的周期平稳特性实施信号的检测和辨识，虽然传统的方法如傅里叶变化等对时间平稳且规则的信号提供了完整的描述，但是对于自然界中大量存在的带有奇异性的不规则信号来讲，传统方法往往难以奏效，比如水声信道中，在许多情况下，目标声源和噪声频谱差异性很大，可以利用目标声源和噪声的频域特征检测异常。
异常检测感知算法：
其中M为小波分层的层数，m代表各个不同的小波层，N(m)表示第m层小波系数的总数。
Figure BDA0001388612480000101
分别表示噪声和异常(存在目标噪声源的小波系数的方差)，V0表示预先提取的噪声的方差参数，V1表示预先提取的异常(存在目标噪声源头)的方差参数，
Figure BDA0001388612480000102
为第m层的小波系数。
(1)对输入语音信号进行抽样，对离散化的语音信号进行分帧，记为Ri(0＜i＜D)，D为帧总数。
(2)对第i帧Ri进行M层的db N小波变换，N称为消失矩，N越大，小波分解得到的高频系数越小，就有更多的高频系数为0，去噪、压缩效果也更好，一般选用N值大的小波，N范围为0-9。这里所得的小波系数记为
Figure BDA0001388612480000103
信号为噪声(不存在目标噪声源)的概率：
Figure BDA0001388612480000111
信号为存在目标噪声源(异常)的概率
Figure BDA0001388612480000112
通过上述两式分别进行计算，如果
Figure BDA0001388612480000113
则为(存在目标噪声源)异常，否则标识此帧为噪声(无目标噪声源)。

Claims (8)

1.一种水下噪声源被动定位方法，采用具有水下噪声源被动定位功能的海上浮标，包括浮标本体，在浮标本体的底部连接有锚链，在锚链的底部连接有重力锚，重力锚固定在海底，其特征在于：在浮标本体上还连接有尼龙绳，在尼龙绳上设置有噪声源定位阵元装置，所述噪声源定位阵元装置包括第一支撑杆、第二支撑杆和第三支撑杆，所述第一支撑杆、第二支撑杆和第三支撑杆均相互交叉垂直布置，共同连接成三维坐标系结构，在第一支撑杆、第二支撑杆和第三支撑杆上分别设置有偶数个声压传感器，且第一支撑杆、第二支撑杆或第三支撑杆上的声压传感器均沿第一支撑杆、第二支撑杆和第三支撑杆的交点呈对称分布；
水下噪声源被动定位方法具体包括以下步骤：
步骤一、采用自适应门限检测方法对水下噪声源判别
首先对各个声压传感器接收到的信号进行处理，对处理过的信号x(t)进行双门限能量检测，以λ1，λ2作为能量检测的上下门限，若累积能量Y>λ1，则认为所选频谱被占用；若Y<λ2，则认为所选频谱空闲；若λ2≤Y≤λ1，则认为利用能量检测法无法判决所选频谱的占用状态，定义该区域为不可靠区S，此时再对处理过的信号x(t)进行循环平稳特征检测；
步骤二、采用多阵元混合定位方法对判别出的水下噪声源进行被动声源定位
(1)利用三维坐标系X、Y、Z三个轴上的以原点为对称点的6m+n个阵元，m≥1,n≥0，且n为偶数，必须在一个轴上，且关于坐标系原点对称确立阵元坐标；
(2)以原点o为中心，对称分布于x,y,z轴上，其中M1-M6为6个阵元，组成XOY、YOZ、XOZ三个平面坐标系上的十字阵，此外有6(m-1)+n个阵元呈十字阵或线列阵布置在其中；
(3)设声源位置为S(x,y,z)，声源距坐标原点距离为r，声源S到阵元M1的距离为R1，声源S到阵元M2的距离为R2，d12,d13,d14为阵元之间的声程差，以d12为例进行说明，d12为目标到阵元1和阵元2的距离之差，则根据目标和声阵的几何位置关系有：
利用M1～M6 6个阵元组成空间坐标系中的三个四元十字阵
分别为xoy面，4个阵元组合M1,M2,M3,M4
分别为yoz面，4个阵元组合M2,M5,M4,M6
分别为xoz面，4个阵元组合M1,M5,M3,M6
A：xoy面十字阵，4个阵元M1,M2,M3,M4
Figure FDA0004114252720000021
得出
Figure FDA0004114252720000022
B：yoz面，4个阵元组合M2,M5,M4,M6
Figure FDA0004114252720000023
得出
Figure FDA0004114252720000024
C:xoz面，4个阵元组合M1,M5,M3,M6
Figure FDA0004114252720000025
得出
Figure FDA0004114252720000026
其中d12＝τ12c，d13＝τ13c，d14＝τ14c，d15＝τ15c，d16＝τ16c，d24＝τ24c，d25＝τ25c，d26＝τ26cτij为声源到阵元i和j的时延差，c为声速，目标声源坐标(x,y,z)；
(4)剩余6(m-1)+n个阵元可任意放在x轴，y轴，z轴上，并与x或y或z轴上的阵元组成间距为2L的线列阵或者间距为4L的十字阵；
以线列阵为例，该阵元位置为Mi(xi,yi,zi)，则有
Figure FDA0004114252720000027
c为声速；ti为第i个阵元所测时延；
以x轴三元线列阵为例
Figure FDA0004114252720000031
得到
Figure FDA0004114252720000032
得出
Figure FDA0004114252720000033
y存在二值：即线列声压阵存在的左右舷模糊问题；
(5)利用各十字阵和线列阵通过计算时延差得到值做求平均处理，得到一个较为准确的目标声源位置。
2.根据权利要求1所述的一种水下噪声源被动定位方法，其特征在于：在第一支撑杆、第二支撑杆和第三支撑杆上均设置有4个声压传感器，沿交点呈两两对称分布，其中距离交点较远的声压传感器与交点之间的距离为距离交点较近的声压传感器与交点之间距离的两倍。
3.根据权利要求1所述的一种水下噪声源被动定位方法，其特征在于：在尼龙绳上间隔设置有用于减少流噪声的布条流苏，在尼龙绳的底部连接有金属坠。
4.根据权利要求1所述的一种水下噪声源被动定位方法，其特征在于：在浮标本体上设置有摄像装置。
5.根据权利要求1所述的一种水下噪声源被动定位方法，其特征在于，所述步骤一：在能量检测法中，判决门限λ是依据期望检测概率
Figure FDA0004114252720000034
的要求设置的，但由于信道环境的不确定性，当接收信号能量值在门限λ附近时很难做出可靠判断，为此，定义一个范围
Figure FDA0004114252720000035
为判决的不可靠区，即λ1＝kλ；λ2＝λ/k，式子中k≥1为不可靠区S的范围调整因子，且有Δ＝10logk；通过检测信道保护带宽内的噪声能量可以近似估计该信道内的信噪比，从而得到信道状态，为了适应信道变化，依据估计的信道信噪比信息动态设置门限；如果实际检测概率比期望检测概率低，就适当增大λ1和λ2之间的距离，即增大k值，如果实际检测概率比期望检测概率高，就减小k值；当信道质量足够好时，令k＝1，此时λ1＝λ2＝λ，实行单门限检测。
6.根据权利要求1所述的一种水下噪声源被动定位方法，其特征在于，所述步骤一：为了避免过于频繁的调整门限，引入期望检测概率范围因子a，当实际检测概率高于
Figure FDA0004114252720000041
低于
Figure FDA0004114252720000042
时才调整门限，若检测概率处于二者之间则不必调整，这里
Figure FDA0004114252720000043
为期望检测概率，a的取值应远小于1，以保证尽可能达到期望的检测概率。
7.根据权利要求1所述的一种水下噪声源被动定位方法，其特征在于：所述阵元为声压阵元。
8.根据权利要求1所述的一种水下噪声源被动定位方法，其特征在于还包括：
步骤三、启动摄像装置，对定位的噪声源进行拍照取证。
CN201710738570.4A 2017-08-25 2017-08-25 具有水下噪声源被动定位功能的海上浮标及定位方法 Active CN107651119B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710738570.4A CN107651119B (zh) 2017-08-25 2017-08-25 具有水下噪声源被动定位功能的海上浮标及定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710738570.4A CN107651119B (zh) 2017-08-25 2017-08-25 具有水下噪声源被动定位功能的海上浮标及定位方法

Publications (2)

Publication Number Publication Date
CN107651119A CN107651119A (zh) 2018-02-02
CN107651119B true CN107651119B (zh) 2023-04-11

Family

ID=61129034

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710738570.4A Active CN107651119B (zh) 2017-08-25 2017-08-25 具有水下噪声源被动定位功能的海上浮标及定位方法

Country Status (1)

Country Link
CN (1) CN107651119B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109747801A (zh) * 2019-01-25 2019-05-14 哈尔滨工程大学 一种用于水下滑翔机的四元十字阵声学检测系统
CN110116786B (zh) * 2019-04-30 2024-01-26 中国人民解放军海军工程大学 可定位海洋资源快速探测装置及其探测定位平台定位方法
CN110307896B (zh) * 2019-06-12 2021-08-06 中国舰船研究设计中心 一种船舶水下辐射噪声测量水听器悬挂系统及布放方法
CN117031402B (zh) * 2023-10-10 2023-12-12 南方海洋科学与工程广东省实验室(珠海) 一种基于小型浮标的声光融合监测设备及方法

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2688595B1 (fr) * 1992-03-13 1994-05-06 Thomson Csf Procede et dispositif de localisation de bruiteurs par une antenne constituee de bouees acoustiques passives.
JP2591490B2 (ja) * 1994-07-29 1997-03-19 日本電気株式会社 海域監視システム
USH1618H (en) * 1995-05-30 1996-12-03 United States Of America Coherent arrays of drifting sonobuoys
FR2821163B1 (fr) * 2001-02-16 2003-05-09 Thomson Marconi Sonar Sas Procede de detection d'objets mobiles au moyen de bouees sous-marines passives
US8279714B2 (en) * 2008-12-05 2012-10-02 Wood Hole Oceanographic Institution Compliant ocean wave mitigation device and method to allow underwater sound detection with oceanographic buoy moorings
CN104181505B (zh) * 2014-08-18 2017-09-01 吉林大学 一种基于近场源定位算法的多目标水声定位方法和系统
CN104237849B (zh) * 2014-09-26 2017-04-12 哈尔滨工业大学 一种双五元十字阵被动声定位融合方法
CN105691556B (zh) * 2016-01-21 2017-07-21 国家海洋局第三海洋研究所 海洋环境噪声源记录浮标
CN207129106U (zh) * 2017-08-25 2018-03-23 山东科技大学 一种具有水下噪声源被动定位功能的海上浮标

Also Published As

Publication number Publication date
CN107651119A (zh) 2018-02-02

Similar Documents

Publication Publication Date Title
CN107651119B (zh) 具有水下噪声源被动定位功能的海上浮标及定位方法
US10725149B1 (en) System and method for autonomous joint detection-classification and tracking of acoustic signals of interest
CN111239829A (zh) 一种拖曳式可控源电磁和水声复合水下目标探测系统及方法
US8195409B2 (en) Passive acoustic underwater intruder detection system
AU2019211746B2 (en) Devices, methods, and systems for underwater surveying
CN108254748A (zh) 基于激光测距和雷达图像的内河船舶超高报警系统及方法
CN110208806B (zh) 一种航海雷达图像降雨识别方法
CN110058245A (zh) 基于云模型的低频主动拖线阵声呐浅海探测效能评估方法
WO2018178913A1 (en) Dual accelerometer vector sensor
CN111983695B (zh) 水下磁探测阵列组建及其数据处理方法、装置和存储介质
Fedele et al. Space-time waves and spectra in the northern adriatic sea via a wave acquisition stereo system
CN108761470B (zh) 一种基于拖缆形态方程解析的目标定位方法
CN207129106U (zh) 一种具有水下噪声源被动定位功能的海上浮标
CN116523822A (zh) 一种基于侧扫声呐的海缆检测识别方法
CN116520247A (zh) 根据水下航行器上的水听器来估算目标声源位置的方法
CN107010181B (zh) 一种基于多波束仰扫的船舶吃水自动检测系统与方法
CN113075732B (zh) 一种消除高分辨率小多道地震地层异常起伏的方法
CN112083428B (zh) 基于声矢量场处理的海洋内波预警监测方法
CN211741599U (zh) 一种拖曳式可控源电磁和水声复合水下目标探测系统
RU2572085C1 (ru) Способ поиска, обнаружения и слежения за местоположением подвижных подводных объектов технической природы в акваториях морей и океанов
CN115291280A (zh) 基于海浪周期去除拖拽式浅剖地层异常起伏的方法
Li et al. Moving target detection and tracking interactive algorithm based on acoustic image
US20230213599A1 (en) Inverse estimation-based radius calculation method and system for ferromagnetic target detection
CN107621639B (zh) 一种基于连续线性和证据累加准则的水下障碍检测方法
Klauson et al. Ship source strength estimation in shallow water

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant