CN112346038B - 一种基于浅海爆炸声传播信号的声信道响应估计方法 - Google Patents
一种基于浅海爆炸声传播信号的声信道响应估计方法 Download PDFInfo
- Publication number
- CN112346038B CN112346038B CN202011131870.4A CN202011131870A CN112346038B CN 112346038 B CN112346038 B CN 112346038B CN 202011131870 A CN202011131870 A CN 202011131870A CN 112346038 B CN112346038 B CN 112346038B
- Authority
- CN
- China
- Prior art keywords
- sound source
- signal
- shallow sea
- sound
- sea
- 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
- 238000004880 explosion Methods 0.000 title claims abstract description 56
- 238000000034 method Methods 0.000 title claims abstract description 26
- 230000004044 response Effects 0.000 title claims abstract description 25
- 238000002474 experimental method Methods 0.000 claims abstract description 23
- 239000002360 explosive Substances 0.000 claims abstract description 22
- 230000010349 pulsation Effects 0.000 claims abstract description 18
- 238000012545 processing Methods 0.000 claims abstract description 10
- 238000005259 measurement Methods 0.000 claims abstract description 5
- 238000001228 spectrum Methods 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 4
- 238000012360 testing method Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 claims description 3
- 230000035939 shock Effects 0.000 claims description 3
- 230000002159 abnormal effect Effects 0.000 claims description 2
- 239000006185 dispersion Substances 0.000 abstract description 6
- 230000005540 biological transmission Effects 0.000 abstract description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000007613 environmental effect Effects 0.000 description 3
- 239000013535 sea water Substances 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000005236 sound signal Effects 0.000 description 1
- 238000006467 substitution reaction 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Computer Networks & Wireless Communication (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
一种基于浅海爆炸声传播信号的声信道响应估计方法,涉及声纳定位。步骤:在深海海域利用同一批次的爆炸声源开展深海实验测量爆炸声源信号,获得爆炸声源的时域信号,用于浅海爆炸声传播实验的声源信号的估计;对深海采集的爆炸声源信号进行处理,获得实际声源信号;对浅海声传播接收信号进行自相关估计二次脉动时延,对实际声源信号的二次脉动时间进行调整,获得浅海声源信号的估计;利用估计得到的浅海声源信号,对浅海声传播接收信号进行解卷积处理,获得真实的浅海声信道响应。可获得真实的浅海声传播宽带时域信号;可获得真实的浅海声传播宽带频散结构。不需要实时采集每次使用的爆炸声源信号。
Description
技术领域
本发明涉及声纳定位,尤其是涉及一种基于浅海爆炸声传播信号的声信道响应估计方法。
背景技术
由于无线电波和光波在海水中传播时都要受到严重的衰减,不能有效地传递信息,因此水面舰艇和水下潜艇无法通过雷达等设备获得水下目标的信息,需要通过声学手段进行水下目标定位。声信号在浅海环境中传播受到海面、水体、海底等因素的影响,是一个复杂的声信道,研究浅海声信道响应对于水下目标定位具有重要意义。同时,由于浅海声信道响应受到水体声速剖面、海底声学参数的影响,因此可以用于浅海水体、海底声学反演。
根据海洋声传播的简正波理论,声信号在海洋中传播可以表达为一系列的简正波叠加在一起往前传播,每一号简正波具有自己的群速度,并且该群速度随着频率变化而变化。因此,在浅海远距离接收宽带声信号中,可以观察到明显的简正波频散结构,该结构与声源距离、海水声速剖面、海底声学参数息息相关,可用于声源和环境参数反演。
浅海声学实验经常采用爆炸声源。爆炸声源具有源级高、频带宽的优点,但爆炸声源存在明显的二次脉动,导致在接收信号的简正波频散结构混乱,难以深入分析利用。因此,开发浅海爆炸声传播信号的处理方法,获得真实的浅海声信道响应,对于浅海声学定位和环境参数反演研究都具有重要的应用价值。
发明内容
本发明的目的在于针对现有技术存在的上述问题,提供能获得真实的浅海声信道响应,用于浅海的声源定位和环境参数反演的一种基于浅海爆炸声传播信号的声信道响应估计方法。
本发明包括以下步骤:
1)在深海海域利用同一批次的爆炸声源开展深海实验测量爆炸声源信号,获得爆炸声源的时域信号,用于浅海爆炸声传播实验的声源信号的估计;
在步骤1)中,在深海海域利用同一批次的爆炸声源开展深海实验测量爆炸声源信号,获得爆炸声源的时域信号的具体方法可为:由实验船布放小艇,由小艇投放爆炸声源,实验船船舷布放水听器采集声源信号,小艇与母船保持一定距离,避免接收水听器限幅。
2)对深海采集的爆炸声源信号进行处理,采用复制直达部分、调整幅度和相位、抵消海面反射部分,对后续信号重复以上过程,即可去除海面反射的影响,获得实际声源信号;
3)对浅海声传播接收信号进行自相关估计二次脉动时延,对实际声源信号的二次脉动时间进行调整,即调整声源信号的冲击波与二次脉动之间无信号部分的时间长度,获得浅海声源信号的估计;
4)利用估计得到的浅海声源信号,对浅海声传播接收信号进行解卷积处理,即将接收信号经傅里叶变换后的频谱除去声源信号傅里叶变换后的频谱,其中异常的极大值频点做一定调整,之后经傅里叶逆变换,获得去除爆炸声源影响能体现浅海声传播效应的接收信号,即真实的浅海声信道响应。
本发明包括深海爆炸声源信号采集、声源信号估计、接收信号自相关估计二次脉动时延、调整声源信号的二次脉动时间、解卷积。与现有方法相比,本发明具有如下优点:
1)可以去除爆炸声源信号的影响,获得真实的声信道响应;
2)可以去除爆炸信号二次脉动的影响;
3)可以获得真实的浅海声传播宽带时域信号;
4)可以获得真实的浅海声传播宽带频散结构。
5)不需要实时采集每次使用的爆炸声源信号。
附图说明
图1是本发明中深海实验采集爆炸声源信号去除海面反射部分获得实际声源信号的示意图;其中,(a)为深海实验采集到的爆炸声源信号,(b)为去除海面反射部分后的爆炸声源信号,(c)为部分放大后的爆炸声源信号。
图2是本发明中对浅海声传播接收信号进行自相关估计二次脉动时延的示意图。
图3是本发明中浅海声传播接收信号解卷积前后的对比的示意图,包括时域信号、时频分析、频谱;其中,(a)为解卷积前的时域信号,(b)为解卷积前的时频分析,(c)为解卷积前的频谱和声源的频谱,(d)为解卷积后的时域信号,(e)为解卷积后的时频分析,(f)为解卷积后的频谱。
图4是本发明基于浅海爆炸声传播信号估计声信道响应的流程框图。
具体实施方式
以下实施例将结合附图对本发明作进一步的说明。
本发明基于浅海爆炸声传播信号的声信道响应估计方法,其基本原理是:
基于爆炸声源的浅海声传播实验,接收水听器采集到的信号是浅海信道响应和声源信号的卷积的结果,因此,为了获得实际的浅海声信道响应,需要去除声源信号的影响。但在实际实验中,特别是在浅海声传播测线实验中,较难直接采集每个爆炸声源的声源信号,因此采用同一批次的爆炸声源在深海实验中采集到的声源信号进行替代。但由于爆炸声源的二次脉动延迟时间不稳定,因此需要对声传播接收信号进行自相关估计二次脉动延迟时间,以进行调整。同时,深海实验采集的爆炸声源信号,较易与海面反射部分在时域上混叠在一起,因此需要去除海面反射部分的影响。
参见图4,本发明基于浅海爆炸声传播信号的声信道响应估计方法,基于上述的原理,具体按照以下步骤实施:
步骤一:在深海海域开展爆炸声源信号的测量实验,由小艇在母船附近投放爆炸声源,由布放于母船船舷的单水听器采集数据,获得爆炸声源的时域信号,如图1(a)所示,图中1指爆炸声源的冲击波,2指爆炸声源的二次脉动,3指冲击波的海面反射(经过一次海面反射,所以相位为负),4指二次脉动的海面反射(经过一次海面反射,所以相位为负);
步骤二:对深海采集的爆炸声源信号进行处理,采用复制直达部分、调整幅度和相位、抵消海面反射部分,去除海面反射的影响,获得实际声源信号,如图1(b)和图1(c)所示;
步骤三:开展浅海声传播实验,由自容式水听器组成潜标布放于接收点,实验船沿着声传播测线投放定深爆炸声源,获得不同距离的浅海爆炸声传播信号;
步骤四:由于每个爆炸声源的二次脉冲时间不固定,因此需要估计二次脉动时间,对浅海爆炸声传播信号进行自相关估计二次脉动时延,如图2所示,对实际声源信号的二次脉动时间进行调整,获得浅海声源信号的估计;
步骤五:利用估计得到的浅海声源信号,对浅海爆炸声传播接收信号进行解卷积处理,获得去除爆炸声源影响的浅海声信道响应,如图3所示。
本实施通过深海实验采集的爆炸声源信号,对浅海爆炸声传播信号的声源信号进行估计,去除爆炸声源信号的影响,获得真实的体现浅海声传播效应的声信道响应。和现有方法相比,本发明可以去除爆炸声源信号的影响,获得真实的声信道响应,可以去除爆炸信号二次脉动的影响,可以获得真实的浅海声传播宽带时域信号、宽带频散结构,不需要实时采集每次使用的爆炸声源信号。
本发明利用同一批次的爆炸声源开展深海实验测量声源信号,用于浅海爆炸声传播实验的声源信号的估计;对深海实验测量的爆炸声源信号,经处理后去除海面反射的影响,获得实际声源信号;对浅海爆炸声传播信号进行自相关估计二次脉动时延,对实际声源信号的二次脉动时间进行调整,获得浅海声源信号的估计;利用估计得到的浅海声源信号,对浅海爆炸声传播接收信号进行解卷积处理,从而得到真实的浅海声信道响应。和现有方法相比,本发明可以去除爆炸声源信号的影响,获得真实的声信道响应,可以去除爆炸信号二次脉动的影响,可以获得真实的浅海声传播宽带时域信号、宽带频散结构,不需要实时采集每次使用的爆炸声源信号。
以上所述仅为本发明的实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明专利保护范围内。
Claims (2)
1.一种基于浅海爆炸声传播信号的声信道响应估计方法,其特征在于包括以下步骤:
1)在深海海域利用同一批次的爆炸声源开展深海实验测量爆炸声源信号,获得爆炸声源的时域信号,用于浅海爆炸声传播实验的声源信号的估计;
2)对深海采集的爆炸声源信号进行处理,采用复制直达部分、调整幅度和相位、抵消海面反射部分,对后续信号重复以上过程,即可去除海面反射的影响,获得实际声源信号;
3)对浅海声传播接收信号进行自相关估计二次脉动时延,对实际声源信号的二次脉动时间进行调整,即调整声源信号的冲击波与二次脉动之间无信号部分的时间长度,获得浅海声源信号的估计;
4)利用估计得到的浅海声源信号,对浅海声传播接收信号进行解卷积处理,即将接收信号经傅里叶变换后的频谱除去声源信号傅里叶变换后的频谱,其中异常的极大值频点做一定调整,之后经傅里叶逆变换,获得去除爆炸声源影响能体现浅海声传播效应的接收信号,即真实的浅海声信道响应。
2.如权利要求1所述一种基于浅海爆炸声传播信号的声信道响应估计方法,其特征在于在步骤1)中,所述在深海海域利用同一批次的爆炸声源开展深海实验测量爆炸声源信号,获得爆炸声源的时域信号的具体方法为:由实验船布放小艇,由小艇投放爆炸声源,实验船船舷布放水听器采集声源信号,小艇与母船保持一定距离,避免接收水听器限幅。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011131870.4A CN112346038B (zh) | 2020-10-21 | 2020-10-21 | 一种基于浅海爆炸声传播信号的声信道响应估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011131870.4A CN112346038B (zh) | 2020-10-21 | 2020-10-21 | 一种基于浅海爆炸声传播信号的声信道响应估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112346038A CN112346038A (zh) | 2021-02-09 |
CN112346038B true CN112346038B (zh) | 2022-07-19 |
Family
ID=74359473
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011131870.4A Active CN112346038B (zh) | 2020-10-21 | 2020-10-21 | 一种基于浅海爆炸声传播信号的声信道响应估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112346038B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009068989A (ja) * | 2007-09-13 | 2009-04-02 | Nec Corp | アクティブソーナー装置およびアクティブソーナー装置による残響除去方法 |
CN105631194A (zh) * | 2015-12-20 | 2016-06-01 | 西北工业大学 | 一种利用模态频散曲线能量差反演海底衰减系数的方法 |
CN106610507A (zh) * | 2016-12-01 | 2017-05-03 | 中国石油天然气集团公司 | 一种去除虚反射干扰的方法和装置 |
CN110488253A (zh) * | 2019-09-05 | 2019-11-22 | 西北工业大学 | 基于互相关函数延迟求和的多途时延差估计方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2853417B1 (fr) * | 2003-04-01 | 2006-07-14 | Thales Sa | Procede de traitement de signaux, et sonar actif le mettant en oeuvre |
CA2829764A1 (en) * | 2011-03-15 | 2012-09-20 | Halliburton Energy Services, Inc. | Acoustic signal processing using model-based adaptive filtering |
-
2020
- 2020-10-21 CN CN202011131870.4A patent/CN112346038B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009068989A (ja) * | 2007-09-13 | 2009-04-02 | Nec Corp | アクティブソーナー装置およびアクティブソーナー装置による残響除去方法 |
CN105631194A (zh) * | 2015-12-20 | 2016-06-01 | 西北工业大学 | 一种利用模态频散曲线能量差反演海底衰减系数的方法 |
CN106610507A (zh) * | 2016-12-01 | 2017-05-03 | 中国石油天然气集团公司 | 一种去除虚反射干扰的方法和装置 |
CN110488253A (zh) * | 2019-09-05 | 2019-11-22 | 西北工业大学 | 基于互相关函数延迟求和的多途时延差估计方法 |
Non-Patent Citations (2)
Title |
---|
浅海声反射损失的实验研究;包青华等;《中国海洋大学学报(自然科学版)》;19870415(第02期);全文 * |
由爆炸声源及传播损失反演海底声学参数;赵梅等;《声学技术》;20090415(第02期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112346038A (zh) | 2021-02-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Demer et al. | 2016 USA–Norway EK80 Workshop Report: Evaluation of a wideband echosounder for fisheries and marine ecosystem science | |
CN101387701B (zh) | 一种基于前向预测的被动时间反转混响抑制方法 | |
CN111580048B (zh) | 一种利用单矢量水听器的宽带声源深度估计方法 | |
Park et al. | A study on propeller noise source localization in a cavitation tunnel | |
CN106249244A (zh) | 一种声矢量圆阵宽带相干信号源方位估计方法 | |
CN113109817A (zh) | 一种矢量水听器布放深度估计方法 | |
Renner | Packet-based ranging with a low-power, low-cost acoustic modem for micro AUVs | |
CN106872572A (zh) | 冰层粗糙表面垂直声反射系数测量方法 | |
Maguer et al. | SLITA: A new slim towed array for AUV applications | |
Cui et al. | Strong scattering targets separation based on fractional Fourier transformation in pulse-to-pulse coherent acoustical Doppler current profilers | |
CN112346038B (zh) | 一种基于浅海爆炸声传播信号的声信道响应估计方法 | |
CN117233737A (zh) | 主动目标深度估计方法、装置、计算机设备和存储介质 | |
CN115236592B (zh) | 一种基于单阵元时频曲线匹配的冰声定位方法 | |
CN100473959C (zh) | 海洋声场相速度的水听器测试方法 | |
CN113126029B (zh) | 适用于深海可靠声路径环境的多传感器脉冲声源定位方法 | |
CN103235312B (zh) | 目标回波强度测量方法 | |
CN112684437B (zh) | 一种基于时域warping变换的被动测距方法 | |
CN114325673A (zh) | 一种未知海洋环境基于双水听器的宽带脉冲声源被动测距方法 | |
CN112834623B (zh) | 一种基于压缩感知技术的声学覆盖层回声降低测量方法 | |
Lu et al. | Array shape estimation method based on extraction of single normal mode | |
CN110208745B (zh) | 一种基于自适应匹配滤波器的水声定位方法 | |
Felisberto et al. | Shallow water tomography with a sparse array during the INTIMATE'98 sea trial | |
Gao | Power cepstra measured in shallow water environments | |
Khyam et al. | OFDM based low-complexity time of arrival estimation in active sonar | |
Gao et al. | New acoustic positioning system for an underwater robot using multiple frequencies |
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 |