详情

原文传递一种船舶近壳体瞬态噪声源的水声和振动融合定位方法

专利名称：	一种船舶近壳体瞬态噪声源的水声和振动融合定位方法
摘要：	本发明属于异常声源定位技术领域，提供了一种船舶近壳体瞬态噪声源的水声和振动融合定位方法，首先在船体表面选择位置布放水听器和加速度计，当船舶舷外近壳体处发生瞬态噪声事件，根据水声信号判定瞬态噪声事件记录时刻t0，接着选取参考水听器以及观测水听器，并对数据进行预处理，本发明通过在船体表面选择位置布放水听器和加速度计，根据水听器可以确定噪声源可能存在的区域，进一步利用船体表面的加速度计，根据振动信号强度对比结果，缩小上述区域范围，具备定位精度高的优点同时通过利用时延数据计算绘制双曲面并求解与已知三维模型的交线，公式复杂程度较低，无需大量的计算，能够在线定位瞬态噪声源位置，具备计算速度快的优点。
专利类型：	发明专利
国家地区组织代码：	湖北;42
申请人：	中国人民解放军海军工程大学
发明人：	余文晶;武世纪;程果;黄秀峰
专利状态：	有效
申请日期：	2023-06-30T00:00:00+0800
发布日期：	2023-11-10T00:00:00+0800
申请号：	CN202310793772.4
公开号：	CN117031404A
代理机构：	北京精翰专利代理有限公司
代理人：	张丹
分类号：	G01S5/28;G01S5/18;G01H17/00;G;G01;G01S;G01H;G01S5;G01H17;G01S5/28;G01S5/18;G01H17/00
申请人地址：	430033 湖北省武汉市硚口区解放大道717号
主权项：	1.一种船舶近壳体瞬态噪声源的水声和振动融合定位方法，其特征在于：包括以下步骤： S1、首先在船体表面选择位置布放水听器和加速度计，用于采集舷外水声信号及壳体振动信号； S2、当船舶舷外近壳体处发生瞬态噪声事件，根据水声信号判定瞬态噪声事件记录时刻t0； S3、接着选取参考水听器以及观测水听器，并对数据进行预处理； S4、逐一计算步骤S3中的观测水听器与参考水听器之间的时延； S5、通过利用已知的船舶结构数据，建立船舶三维模型，然后根据双曲定位原理，且由步骤S4中的已知时延数据绘制双曲面，再结合船舶三维模型计算两者的交线，最后确定噪声源的可能区域； S6、在步骤S5的基础上，进一步利用船体表面的加速度计，按照振动信号强度对比，缩小步骤S5中的噪声源可能存在的区域范围。 2.如权利要求1所述一种船舶近壳体瞬态噪声源的水声和振动融合定位方法，其特征在于：在步骤S1中，所述水听器布放在船体外表面，加速度计布放在船体内表面，且水听器数量M≥3，加速度计数量L≥2。 3.如权利要求1所述一种船舶近壳体瞬态噪声源的水声和振动融合定位方法，其特征在于：在步骤S2中，所述的瞬态噪声事件判定的标准为：任一水听器监测的舷外水声信号总级超过报警阈值，或任一加速度计所监测的壳体振动信号总级超过报警阈值；单次超阈值的信号持续时间，在1s以内；且两次超阈值的信号之间，存在0.2s以上的间隔。 4.如权利要求1所述一种船舶近壳体瞬态噪声源的水声和振动融合定位方法，其特征在于：在步骤S3中，所述水听器为参考点水听器或者为观测点水听器，并对数据进行预处理的方法如下： a、选择包含t0时刻前、后各ΔT的数据作为瞬态噪声数据，其中，ΔT可为200ms； b、t0时刻向前1秒为背景噪声记录时刻，记作tb，选择包含tb时刻前、后各ΔT的数据作为背景噪声数据； c、对瞬态噪声数据和背景噪声数据分别进行傅里叶变换，并在频域上划分为N个频段，中心频率分别记作f1、f2...fN； d、分别对每个水听器，计算瞬态噪声、背景噪声各个频段的能量，分别记作E1m、E2m...ENm和E01m、E02m...E0Nm,其中，m表示水听器编号，有m＝1,2,…,M； e、计算每个水听器每个频段的能量超标量，记作ΔEnm＝Enm-E0nm,其中，n表示频段编号，有n＝1,2,…,N； f、对比得到ΔEnm中取最大值的频段，即ΔEGm＝max(ΔEnm),该频段编号fG为瞬态噪声事件的主要特征频段； g、对比选取ΔEGm中最大的三个值，记作ΔEGH、ΔEGI、ΔEGJ，即编号为H、I和J的3个水听器，其中任选1个作为参考点水听器，记作H号，其他2个作为观测点水听器，记作I号和J号。 5.如权利要求1和4所述一种船舶近壳体瞬态噪声源的水声和振动融合定位方法，其特征在于：在步骤S4中，逐一计算观测水听器与参考水听器之间时延的方法如下： Ⅰ.根据步骤g，选取参考点水听器和观测点水听器，其相应的瞬态噪声数据记作：PH、PI和PJ； Ⅱ.对PH、PI和PJ进行希尔伯特变换，并对变换的结果取绝对值，得到QH、QI和QJ； Ⅲ.分别计算QH、QI和QJ的最大值，记作QHM、QIM和QJM； Ⅳ.在0至1之间选取一个幅值系数a，在QH、QI和QJ曲线上，首次取值为QHM·a、QIM·a和QJM·a所对应的时刻，分别记作tH、tI和tJ； Ⅴ.H号水听器相对于I水听器的时延为：τHI＝tH-tI；H号水听器相对于J水听器的时延为：τHI＝tH-tI。 6.如权利要求1所述一种船舶近壳体瞬态噪声源的水声和振动融合定位方法，其特征在于：步骤S5通过以下过程具体进行： (1)对船舶壳体建立固定坐标系下的三维模型，标记参考水听器(xH yH zH)和观测水听器(xI yI zI)和(xJ yJ zJ)的坐标位置； (2)根据步骤(4)计算的时延差τHI和τHJ，建立定位方程组，并绘制定位双曲面；其中，v为声波传播速度； (3)结合步骤(2)得到的双曲面和步骤(1)的船舶壳体三维模型，计算两者交线，为瞬态噪声源可能存在的区域。 7.如权利要求1所述一种船舶近壳体瞬态噪声源的水声和振动融合定位方法，其特征在于：在步骤S6中，缩小步骤S5中噪声源可能存在的区域范围采用以下方法： (i)根据加速度计在船舶壳体分布，将其划分为若干网格区域； (ii)对所有加速度计，计算t0时刻前、后各ΔT时间段内，位于编号fG频段内的振动信号能量，记作εl，l表示频段编号，有l＝1,2,…,L； (iii)对所有加速度计，计算tb时刻前、后各ΔT时间段内，位于编号fG频段内的振动信号能量，记作ε0l，l表示频段编号，有l＝1,2,…,L； (iv)计算所有加速度计，特征频段内的振动信号能量增加量Δεl＝εl-ε0l； (v)对比选取不同加速度计特征频段内的振动信号能量增加量的最大值ΔεK，编号K加速度计所对应的网格区域为噪声源可能存在的区域。 8.如权利要求1所述一种船舶近壳体瞬态噪声源的水声和振动融合定位方法，其特征在于：步骤S6缩小噪声源可能存在的区域范围的方法如下： ①取权利要求6中步骤(3)和权利要求7中步骤(v)所得到的区域的相交部分，作为瞬态噪声源存在区域； ②如果权利要求6中步骤(3)和权利要求7中步骤(v)所得到的区域不存在相交部分，则对比选取不同加速度计特征频段内的振动信号能量增加量的第二大值ΔεK1，取编号K1加速度计所对应的网格区域和权利要求6中步骤(3)所得到的区域的相交部分； ③以此类推，直到选出两类区域相交部分作为瞬态噪声源存在区域。