CN115220026B - 一种基于冰面水平阵列的匹配场跨冰定位方法 - Google Patents
一种基于冰面水平阵列的匹配场跨冰定位方法 Download PDFInfo
- Publication number
- CN115220026B CN115220026B CN202210770497.XA CN202210770497A CN115220026B CN 115220026 B CN115220026 B CN 115220026B CN 202210770497 A CN202210770497 A CN 202210770497A CN 115220026 B CN115220026 B CN 115220026B
- Authority
- CN
- China
- Prior art keywords
- ice
- sound
- sound source
- field
- array
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 26
- 239000013598 vector Substances 0.000 claims abstract description 28
- 238000001228 spectrum Methods 0.000 claims abstract description 27
- 239000005457 ice water Substances 0.000 claims abstract description 25
- 230000005236 sound signal Effects 0.000 claims abstract description 8
- 230000005540 biological transmission Effects 0.000 claims abstract description 6
- 230000008878 coupling Effects 0.000 claims description 22
- 238000010168 coupling process Methods 0.000 claims description 22
- 238000005859 coupling reaction Methods 0.000 claims description 22
- 239000013535 sea water Substances 0.000 claims description 22
- 230000001360 synchronised effect Effects 0.000 claims description 11
- 230000003595 spectral effect Effects 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 6
- 230000004807 localization Effects 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 4
- 238000010606 normalization Methods 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 abstract description 3
- 230000000875 corresponding effect Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 238000003491 array Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000010792 warming Methods 0.000 description 1
Images
Classifications
-
- 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
- G01S11/00—Systems for determining distance or velocity not using reflection or reradiation
- G01S11/14—Systems for determining distance or velocity not using reflection or reradiation using ultrasonic, sonic, or infrasonic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/16—Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
- G01V1/20—Arrangements of receiving elements, e.g. geophone pattern
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Radar, Positioning & Navigation (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明的目的在于提供一种基于冰面水平阵列的匹配场跨冰定位方法，包括以下步骤：通过加速度计阵列采集冰下声源发出的跨冰声信号；提取跨冰声信号频谱向量；建立声源搜索距离和深度搜索网格序列；计算每个水平阵元与声源的实际距离序列；获取冰水环境冰声物理参数；求解每个水平阵元接收的跨冰声传播声场；提取不同期望位置的理论频谱向量；基于匹配场处理获得定位模糊曲面；输出最大模糊值对应的期望声源位置。本发明采用加速度参量作为定位参量，实现对冰下声源的跨冰定位，具有布放方便、坐标可精确校准以及经济等优势。
Description
技术领域
本发明涉及的是一种水声定位方法，具体地说是极地水声定位方法。
背景技术
随着全球变暖效应，北极航道的开发和利用越来越成为现实，探索北极、保障北极航行安全是北极水声学的重要研究方向，研究在北极海域的水声定位方法具有重大意义。
匹配场定位是水声定位的经典方法，原理是通过海水中布放垂直水听器阵或海底布放水平阵采集水中的声压场，并与基于水声传播模型计算的期望声源位置的拷贝声场进行匹配相关，获得不同声源深度和距离模糊度，模糊度最大值对应的深度和距离值即估计的声源位置。
水中匹配场定位关注的是水中声压的波动成分。由于北极常年被数米厚的冰层覆盖，常规的匹配场定位存在极地适应性问题。
一个是较厚冰层的阻隔，使得在海水中布放垂直水听器阵和海底布放水平阵工作难度大，冰下环境复杂，水听器等设备面临布放条件困难，阵元坐标难以校准等问题突出。
另一点是由于目前的匹配场定位方法仅考虑了海水环境中的声学参数，如海水声压和质点振速参量，没有涉及海冰冰层这种板状构型弹性介质的声学参数，如冰层位移和加速度分量，因此也不适应极地海域，特别是有数米厚冰层覆盖的极地环境。
因此亟待开发一种在冰上布放阵列用于匹配场定位的方法。
发明内容
本发明的目的在于提供能克服匹配场处理无法适应有冰层覆盖海域等难题的一种基于冰面水平阵列的匹配场跨冰定位方法。
本发明的目的是这样实现的：
本发明一种基于冰面水平阵列的匹配场跨冰定位方法，其特征是：
(1)在冰面布防加速度计阵列，通过加速度计阵列采集冰下声源发出的跨冰声信号；
(2)基于傅里叶变换提取跨冰声信号频谱向量；
(3)建立声源搜索距离ri和深度di搜索网格序列；
(4)依据声源搜索距离ri和方位角θ计算每个水平阵元与声源的实际距离序列
式中xn和yn分别为第n个水平阵元的x轴坐标和y轴坐标；
(5)获取冰水环境冰声物理参数；
(6)基于冰水耦合模型求解每个水平阵元接收的跨冰声传播声场：假设冰下声源发出的声波为简谐信号，海冰-海水耦合环境会改变声源信号的幅度和相位，改变程度取决于声源位置以及海冰-海水耦合环境参数，根据冰声简正波传播理论，接收信号的频域表达式为
式中G表示接收信号的频谱，z为加速度计的深度，zs为声源深度，r为声源相对于加速度计的距离，
表示声源频谱，kn为第n阶模态的水平波数，
为冰水耦合声场的第n阶简正波模态幅度；
(7)提取不同期望位置的理论频谱向量：将步骤(3)和(4)所得的搜索网格、以及步骤(5)的冰水环境物理参数带入步骤(6)求解每个阵元接收的跨冰声传播声场，采用冰水耦合简正波模型获得不同位置a＝(zs,r)跨冰声信号的理论频谱向并归一化为
式中，对于同步阵列，G(a)表示为频谱，包括频谱的幅度和相位，对于非同步阵列，G(a)表示为频谱的幅度；
(8)基于匹配场处理获得定位模糊曲面，常规匹配场处理(MFP)表示为
S(a)＝|wH(a)d(atrue)|2
式中S(a)表示位置匹配模糊函数，当拷贝场向量与真实场向量d(atrue)一致时，模糊函数为最大值，拷贝场向量对应的位置a即为估计的声源位置，该式被表述为
S(a)＝wH(a)K(atrue)w(a)
式中K(atrue)＝d(atrue)dH(atrue)，为真实场向量的协方差矩阵；
(9)根据步骤(8)结果，输出最大模糊值对应的期望声源位置：
2、步骤(5)中获取的冰水环境冰声物理参数包括：总深度、冰层纵波声速、横波声速、冰层密度、纵波衰减系数、横波衰减系数、海水声速、海水密度、沉积层纵波声速、沉积层密度、纵波衰减系数。
本发明的优势在于：
1、本发明提出在冰面上布放加速度计水平阵或者检波器阵列作为跨冰定位声信号接收设备，以替代在数米厚冰层覆盖的极地海域布放常规或矢量水听器等声呐设备，具有布放方便、坐标可精确校准以及经济等优势。
2、本发明提出了将匹配场处理引入有冰层覆盖的极地海域，并进行极地环境适应性的改进。
3、本发明在水平阵水中匹配场定位的基础上，开创性的提出了一种适用于极地冰区环境的跨冰定位方法，采用加速度参量作为定位参量，实现对冰下声源的跨冰定位。
附图说明
图1为本发明的跨冰定位示意图；
图2为本发明的流程图；
图3为本发明的水平阵仿真示意图；
图4为本发明中的非同步水平阵匹配场跨冰定位结果；
图5为本发明中的同步水平阵匹配场跨冰定位结果。
具体实施方式
下面结合附图举例对本发明做更详细地描述：
结合图1-5，本发明的应用场景为跨冰定位，通过在冰面布放加速度计等水平阵阵元，采集冰下声源发出的跨冰声信号，用于跨冰定位。
如图2所示，本发明包括以下步骤：
(1)在冰面布防加速度计阵列，通过加速度计阵列采集冰下声源发出的跨冰声信号；
(2)基于傅里叶变换提取跨冰声信号频谱向量；
(3)建立声源搜索距离ri和深度di搜索网格序列；
(4)依据声源搜索距离ri和方位角θ计算每个水平阵元与声源的实际距离序列
式中xn和yn分别为第n个水平阵元的x轴坐标和y轴坐标；
(5)获取冰水环境冰声物理参数；
(6)基于冰水耦合模型求解每个水平阵元接收的跨冰声传播声场：假设冰下声源发出的声波为简谐信号，海冰-海水耦合环境会改变声源信号的幅度和相位，改变程度取决于声源位置以及海冰-海水耦合环境参数，根据冰声简正波传播理论，接收信号的频域表达式为
式中G表示接收信号的频谱，z为加速度计的深度，zs为声源深度，r为声源相对于加速度计的距离，
表示声源频谱，kn为第n阶模态的水平波数，
为冰水耦合声场的第n阶简正波模态幅度；
(7)提取不同期望位置的理论频谱向量：将步骤(3)和(4)所得的搜索网格、以及步骤(5)的冰水环境物理参数带入步骤(6)求解每个阵元接收的跨冰声传播声场，采用冰水耦合简正波模型获得不同位置a＝(zs,r)跨冰声信号的理论频谱向并归一化为
式中，对于同步阵列，G(a)表示为频谱，包括频谱的幅度和相位，对于非同步阵列，G(a)表示为频谱的幅度；
(8)基于匹配场处理获得定位模糊曲面，常规匹配场处理(MFP)表示为
S(a)＝|wH(a)d(atrue)|2
式中S(a)表示位置匹配模糊函数，当拷贝场向量与真实场向量d(atrue)一致时，模糊函数为最大值，拷贝场向量对应的位置a即为估计的声源位置，该式被表述为
S(a)＝wH(a)K(atrue)w(a)
式中K(atrue)＝d(atrue)dH(atrue)，为真实场向量的协方差矩阵；
(9)根据步骤(8)结果，输出最大模糊值对应的期望声源位置：
本发明应用实例：
步骤1、通过布放在冰面的加速度计阵列采集冰下声源发出的跨冰声信号；相比于在冰下或海底布放基线，这一步骤具有布放简单、实施便捷、经济高效、基元坐标精确校准、作业安全性高的优势。如图3所示，加速度计水平阵的阵元数为50个，等间隔5m。声源距离原点阵元r＝700m处，深度d＝20m，方位角θ＝60°；
步骤2、基于傅里叶变换提取信号频谱向量；
步骤3、建立声源距离ri和深度di搜索网格序列；距离搜索范围1000，间隔1m，深度搜索范围50m，间隔1m。
步骤4、依据声源搜索距离ri和方位角θ计算每个水平阵元与声源的实际距离序列
式中xn和yn分别为第n个水平阵元的x轴坐标和y轴坐标。
步骤5、获取冰水环境冰声物理参数。环境参数设置：总深度50m，其中冰层厚度10m，水深35m，沉积层5m。冰层纵波声速3800m/s，横波声速1900m/s，密度0.9g/cm3，纵波衰减系数0.1dB/λ，横波衰减0.3dB//λ；海水声速1500m/s；密度1.0g/cm3，沉积层纵波声速1800m/s，密度1.2g/cm3，纵波衰减0.01dB/λ。
步骤6、基于冰水耦合模型求解每个水平阵元接收的跨冰声传播声场。假设冰下声源发出的声波为简谐信号，海冰-海水耦合环境会改变声源信号的幅度和相位，改变程度取决于声源位置以及海冰-海水耦合环境参数。根据海冰-海水耦合声传播理论，接收信号的频域表达式近似为
式中G表示接收信号的频谱，z为加速度计的深度，加速度计布放在冰面上，其实际值设置为0，zs为声源深度，r为声源相对于加速度计的距离，
表示声源频谱，kn为第n阶模态的水平波数，
为冰水耦合声场的第n阶模态幅度。
常规的水声简正波理论是建立在海水的波动方程上求解的，此处的一点改进是采用通过联合求解海冰中声传播的纳维(Navier)方程和海水中声传播的亥姆霍兹方程(Helmholtz)方程获得冰水耦合声场的水平波数kn和简正波模态
使得并不适合冰层覆盖海域进行水声定位的常规匹配场处理得到了极地环境适应性改进。
步骤7、提取不同期望位置的理论频谱向量。将步骤3和4所得的搜索网格、以及步骤5的冰水环境物理参数带入步骤6求解每个阵元接收的跨冰声传播声场。采用冰水耦合简正波模型获得不同位置a＝(zs,r)跨冰声信号的理论频谱向并归一化为
式中，对于同步阵列，G(a)表示为频谱，包括频谱的幅度和相位，对于非同步阵列，G(a)表示为频谱的幅度。
步骤8、基于匹配场处理获得定位模糊曲面。常规匹配场处理(MFP)表示为
S(a)＝|wH(a)d(atrue)|2
式中S(a)表示位置匹配模糊函数，当拷贝场向量与真实场向量d(atrue)一致时，模糊函数为最大值，拷贝场向量对应的位置a即为估计的声源位置。该式也可以被表述为
S(a)＝wH(a)K(atrue)w(a)
式中K(atrue)＝d(atrue)dH(atrue)，为真实场向量的协方差矩阵。如图4所示，为匹配场处理得到的归一化定位模糊曲面的分贝值。声源的实际位置为距离第一阵元700m，深度20m。
步骤9、根据步骤8结果，估计的声源位置为
如图4所示，为非同步水平阵列的跨冰定位结果图，如图中白圈标示，当期望声源距离700m，声源深度为20m时，定位模糊曲面的值最大，即估计的声源距离和深度分别为700m和20m，与实际声源位置一致。如图5所示为同步水平阵列的跨冰定位结果图，与图4相比，具有更低的背景噪声。同步阵列可以同时利用跨冰声响应的幅度和相位分量，而非同步阵列仅可以利用声场的幅度分量，导致同步水平阵列的定位模糊函数的背景噪声比非同步阵列的背景噪声更低，跨冰定位效果更好。
通过提出两种同步型水平阵以及非同步型冰面加速度计水平阵列匹配场定位方法，实现了针对极地冰区冰下声源的跨冰定位这一开创性应用。
Claims (2)
1.一种基于冰面水平阵列的匹配场跨冰定位方法，其特征是：
(1)在冰面布放加速度计阵列，通过加速度计阵列采集冰下声源发出的跨冰声信号；
(2)基于傅里叶变换提取跨冰声信号频谱向量；
(3)建立声源搜索距离ri和深度di搜索网格序列；
(4)依据声源搜索距离ri和方位角θ计算每个水平阵元与声源的实际距离序列
式中xn和yn分别为第n个水平阵元的x轴坐标和y轴坐标；
(5)获取冰水环境冰声物理参数；
(6)基于冰水耦合模型求解每个水平阵元接收的跨冰声传播声场：假设冰下声源发出的声波为简谐信号，海冰-海水耦合环境会改变声源信号的幅度和相位，改变程度取决于声源位置以及海冰-海水耦合环境参数，根据冰声简正波传播理论，接收信号的频域表达式为
式中G表示接收信号的频谱，z为加速度计的深度，zs为声源深度，r为声源相对于加速度计的距离，
表示声源频谱，kn为第n阶模态的水平波数，
为冰水耦合声场的第n阶简正波模态幅度；采用通过联合求解海冰中声传播的纳维方程和海水中声传播的亥姆霍兹方程获得冰水耦合声场的水平波数kn和简正波模态
使得并不适合冰层覆盖海域进行水声定位的常规匹配场处理得到了极地环境适应性改进；
(7)提取不同期望位置的理论频谱向量：将步骤(3)和(4)所得的搜索网格、以及步骤(5)的冰水环境物理参数带入步骤(6)求解每个阵元接收的跨冰声传播声场，采用冰水耦合简正波模型获得不同位置a＝(zs,r)跨冰声信号的理论频谱向量并归一化为
式中，对于同步阵列，G(a)表示为频谱，包括频谱的幅度和相位，对于非同步阵列，G(a)表示为频谱的幅度；
(8)基于匹配场处理获得定位模糊曲面，常规匹配场处理(MFP)表示为S(a)＝|wH(a)d(atrue)|2
式中S(a)表示位置匹配模糊函数，当拷贝场向量与真实场向量d(atrue)一致时，模糊函数为最大值，拷贝场向量对应的位置a即为估计的声源位置，该式被表述为
S(a)＝wH(a)K(atrue)w(a)
式中K(atrue)＝d(atrue)dH(atrue)，为真实场向量的协方差矩阵；
(9)根据步骤(8)结果，输出最大模糊值对应的期望声源位置：
2.根据权利要求1所述的一种基于冰面水平阵列的匹配场跨冰定位方法，其特征是：步骤(5)中获取的冰水环境冰声物理参数包括：总深度、冰层纵波声速、横波声速、冰层密度、纵波衰减系数、横波衰减系数、海水声速、海水密度、沉积层纵波声速、沉积层密度、纵波衰减系数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210770497.XA CN115220026B (zh) | 2022-06-30 | 2022-06-30 | 一种基于冰面水平阵列的匹配场跨冰定位方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210770497.XA CN115220026B (zh) | 2022-06-30 | 2022-06-30 | 一种基于冰面水平阵列的匹配场跨冰定位方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115220026A CN115220026A (zh) | 2022-10-21 |
| CN115220026B true CN115220026B (zh) | 2023-03-31 |
Family
ID=83610453
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210770497.XA Active CN115220026B (zh) | 2022-06-30 | 2022-06-30 | 一种基于冰面水平阵列的匹配场跨冰定位方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115220026B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116908914B (zh) * | 2023-07-19 | 2024-06-25 | 哈尔滨工程大学 | 一种高分辨偏振滤波弹性波极地方位估计方法 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8599649B1 (en) * | 2010-03-09 | 2013-12-03 | The United States Of America As Represented By The Secretary Of The Navy | Laser-based method of detecting underwater sound through an ice layer |
| CN112986902B (zh) * | 2021-02-23 | 2022-07-19 | 自然资源部第三海洋研究所 | 单检波器跨冰层对水中宽带声源分频段方位估计方法 |
| CN113359183B (zh) * | 2021-05-25 | 2023-09-29 | 哈尔滨工程大学 | 一种针对极地冰层的震源定位方法 |
-
2022
- 2022-06-30 CN CN202210770497.XA patent/CN115220026B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN115220026A (zh) | 2022-10-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111580048B (zh) | 一种利用单矢量水听器的宽带声源深度估计方法 | |
| CN102645265A (zh) | 一种基于虚拟时间反转镜的舰船辐射噪声级测量方法 | |
| CN112987004A (zh) | 一种浅海环境下基于水平阵列的水面水下目标分类方法 | |
| CN104091048A (zh) | 被动垂直双矢量水听器的水下目标深度测定方法 | |
| CN115220026B (zh) | 一种基于冰面水平阵列的匹配场跨冰定位方法 | |
| CN107016159A (zh) | 本征值确定方法及装置 | |
| CN112099019A (zh) | 一种水下声矢量探测装置 | |
| Hodgkiss et al. | Direct measurement and matched-field inversion approaches to array shape estimation | |
| CN108572349B (zh) | 一种深海环境下基于模型计算的声源定深方法 | |
| CN109632258A (zh) | 一种基于矢量传感器的收发分离的海洋内波声学检测方法 | |
| CN115236592B (zh) | 一种基于单阵元时频曲线匹配的冰声定位方法 | |
| CN114859420B (zh) | 一种浅海中目标分选及水下目标运动态势和深度判断方法 | |
| CN114909610B (zh) | 一种水下油气管道泄漏检测与定位的方法及控制系统 | |
| CN113126029B (zh) | 适用于深海可靠声路径环境的多传感器脉冲声源定位方法 | |
| Zhang et al. | Simulation of ship radiated noise field in deep sea based on statistical characteristics of sound source | |
| Zhao et al. | Improving Scholte-Wave vibration signal recognition based on polarization characteristics in coastal waters | |
| Kim et al. | Performance comparison of FFT-based and GCC-PHAT time delay estimation schemes for target azimuth angle estimation in a passive SONAR array | |
| Yayu et al. | Research on location of underwater sound source target in deep sea sound field based on bellhop model | |
| EP1089092A1 (en) | Method and apparatus for measuring physical parameters from the phase and amplitude histories of an acoustic signal | |
| RU2740334C1 (ru) | Способ приема сейсмоакустической и гидроакустической волн у дна водоема и устройство для его осуществления | |
| CN109100680B (zh) | 一种相切双弧线型的九元地声传感器阵列定位方法 | |
| CN110208745B (zh) | 一种基于自适应匹配滤波器的水声定位方法 | |
| Peng et al. | Modelling and Calculation of Seismic Wave Field Caused by Moving Ship | |
| Lee et al. | Study of Imaging of Submarine Bubble Plume with Reverse Time Migration | |
| Oh et al. | Prototype development of an underwater sound measurement module using four hydrophones |
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 |