CN101183150A - 基于高频时反阵的时反发射聚焦反波束形成方法 - Google Patents
基于高频时反阵的时反发射聚焦反波束形成方法 Download PDFInfo
- Publication number
- CN101183150A CN101183150A CNA2007101644383A CN200710164438A CN101183150A CN 101183150 A CN101183150 A CN 101183150A CN A2007101644383 A CNA2007101644383 A CN A2007101644383A CN 200710164438 A CN200710164438 A CN 200710164438A CN 101183150 A CN101183150 A CN 101183150A
- Authority
- CN
- China
- Prior art keywords
- omega
- time
- psi
- formula
- emission
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 13
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 10
- 230000006835 compression Effects 0.000 claims abstract description 4
- 238000007906 compression Methods 0.000 claims abstract description 4
- 230000015572 biosynthetic process Effects 0.000 claims description 8
- 238000004613 tight binding model Methods 0.000 claims description 4
- 238000012544 monitoring process Methods 0.000 claims description 3
- 230000021615 conjugation Effects 0.000 claims description 2
- 230000000694 effects Effects 0.000 claims description 2
- 230000005284 excitation Effects 0.000 claims description 2
- 238000005070 sampling Methods 0.000 claims description 2
- 238000000926 separation method Methods 0.000 claims description 2
- 238000001228 spectrum Methods 0.000 claims description 2
- 238000006467 substitution reaction Methods 0.000 claims description 2
- 238000001514 detection method Methods 0.000 abstract description 7
- 238000010586 diagram Methods 0.000 description 8
- 238000005259 measurement Methods 0.000 description 4
- 238000013461 design Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 230000007613 environmental effect Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000003760 hair shine Effects 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 230000007306 turnover Effects 0.000 description 1
- 230000004304 visual acuity Effects 0.000 description 1
Images
Landscapes
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
一种基于高频时反阵的时反发射聚焦反波束形成方法,该波束形成方法采用高频时反阵,通过发射声波照射目标,再接收目标的回波信号,将接收到的回波信号在时间上进行翻转然后再发射出去,时反信号会在目标处聚焦,表现为时间波形压缩和空间聚焦,再次接收增强的目标回波信号,形成接收聚焦波束。本发明提供一种采用高频时反镜适应浅海声场、降低混响、提高对浅海的探测能力的基于高频时反阵的时反发射聚焦反波束形成方法。
Description
技术领域
本发明涉及海洋资源开发和利用中的声探测技术,尤其是一种浅海降低混响、提高信混比的波束形成方法。
背景技术
浅海声场是上有海面、下有海底的波导,声波在其传播存在直达和上下界面反射路径,而现有水声设备假设了海洋是个自由场、无边界的,采用平面波模型进行处理,因此,在浅海或近岸海域工作时,它们的性能都要下降。对于主动探测工作方式,存在的缺点是:1、多途导致时延扩展严重;2、混响大、回混比低;3、探测性能弱。
发明内容
为了抑制混响、适应浅海声场、并提高浅海小信号探测能力,本发明提供一种采用高频时反镜适应浅海声场、降低混响、提高对浅海小目标探测能力的基于高频时反阵的时反发射聚焦反波束形成方法。
本发明解决其技术问题所采用的技术方案是:
一种基于高频时反阵的时反发射聚焦反波束形成方法,该波束形成方法采用高频时反阵,通过发射声波照射目标,再接收目标的回波信号,将接收到的回波信号在时间上进行翻转然后再发射出去,时反信号会在目标处聚焦,表现为时间波形压缩和空间聚焦,再次接收增强的目标回波信号,形成接收聚焦波束。
作为优选的一种方案:在波导环境下,基于简正波模型,在频域的波动方程,即Helmholtz方程为(1):
[▽2-k2]G(r,z)=0. (1)
式中:G(r,z)为格林函数,r为声源到测量处的距离,k2=ω2/c2(z)为波数,ω为信号源角频率,c(z)为声速梯度,利用分离变量法,设G(r,z)=Φ(r)Ψ(z),将其代入式(1)经整理可得到深度方程(2)和距离方程(3):
式中:kr、kz分别为波数的水平分量和垂直分量,它们满足:
式(2)是经典的Strum-Liouville特征值问题;式(3)是一阶Bessel方程,其解为零阶Hankel函数;假设激励声源位于水下zs处,在远场条件下,忽略时间因子e-jωt,它产生的声场可近似为(5):
式中:ρ为介质密度,krm为第m号水平波数;
时反处理系统由垂直线阵(vertical line array,VLA)、时反阵(timereversal array,TRA)和探察源(probe source,PS)组成,如图1所示。它将接收的信号按先到后发、后到先发的顺序完成信号的时间反转发射,VLA用于监测时反发射信号在声场中的聚焦特性,PS位于水下zps处,与TRA的距离为R。
设PS发射信号为s(t),从PS到TRA中第j(j=1,2,…J)个阵元处的信道格林函数可表示为G(R,zj,ω),则第j个阵元接收的声压场在频域表示为(6):
Pj(ω)=S(ω)G(R,zj,ω). (6)
式中: 为信号的频谱,TRA将接收到的信号进行时间反转发射,对应于频域内进行共轭;因而时反发射信号表示为(7):
时反发射后在观测点r的声压场在频域表示为(8):
它是N个阵元产生的时反声场的综合效果,其对应的时域表示为(9):
参照式(5),整理式(8)得(10):
TRA垂直布满整个波导且充分采样,利用模深度函数的正交性质得到(11):
取n=m并在j上积分得到(12):
当r=R时,
式(13)对于能在波导中传播的有效简正模,krm近似为常数,同时,
最后,式(13)近似为(15):
PTR(R,z,ω)=Aδ(z-zs)S*(ω) (15)
式中A为一常量,求式(15)的傅立叶反变换得到时域的声压场如下(16):
pTR(R,z,t)=Aδ(z-zs)s(-t). (16)
当z=zs时,有下式(17)成立:
pTR(R,zs,t)=As(-t). (17)
由式(13)和式(17)可见,通过r=R和z=zs完成时反发射信号在声源位置处的空间聚焦,进而完成时反发射聚焦反波束形成。
从式(17)又可见,时反发射信号在声源位置处的信号波形表现为源信号在时间上的翻褶,克服了多途径传播产生的信号在时间上的扩展。因此,时反处理利用时反发射信号在声场中聚焦点所对应的距离和深度信息完成对目标的定位。
本发明的技术构思为:波导的多路径包含目标的位置信息,若能利用这种多路径结构知识则可以提高水声设备的探测性能。波动方程的时反不变性和声场的收发互易性确保了稳定声场具有时反聚焦特性。时反处理则利用时反聚焦特性完成对目标在距离和深度上的定位。
高频时反镜适应浅海声场分布特性,基于波动方程的时反不变性和稳定声场的收发互易性,利用时反发射聚焦反波束形成将更多的能量照射目标、实现混响降低。高频时反镜是收发合置的垂直线阵,可以单基元发射也可以全阵发射,通过发射声波照射目标,再接收目标的回波信号,将接收到的回波信号在时间上进行翻转然后再发射出去,时反信号会在目标处聚焦,表现为时间波形压缩和空间聚焦,再次接收增强的目标回波信号,通过接收聚焦波束形成结合匹配滤波完成对目标的检测与定位。
本发明的有益效果主要表现在:1、采用高频时反镜适应浅海声场;2、降低混响;3、提高对浅海的探测能力。
附图说明
图1是时反处理示意图。
图2是实验室波导声学参数示意图。
图3是基于波导环境模拟的时反聚焦结果示意图。
图4是实验室波导环境实测的时反聚焦结果示意图。
图5是抑制混响的湖上实验结果示意图。
图6是实验室波导探测弱目标实验结果示意图。
图7是BS发射情况下对目标进行估距的结果示意图。
图8是TR发射情况下对目标进行估距的结果示意图。
具体实施方式
下面结合附图对本发明作进一步描述。
参照图2,高频时反镜性能参数设计及实测结果:表1列出了高频时反镜的物理参数。图2为实验室波导实验环境的声学参数模型,采用时反阵、探查源以及监视声场时反聚焦的垂直线阵。在时反聚焦实验中各设备布置如图1所示。TRA和VLA从水下4cm布到水下128cm,两阵相距9.6m。PS与TRA相距9.6m,位于水下74cm处。PS源发射一个中心频率为20kHz,带宽为2.4kHz的窄带升余弦包络PCW信号。图3和图4分别给出了实验室波导环境模拟的和实测的时反处理空间(深度-距离)聚焦结果。表2为对应的垂直分辨力和水平分辨力的理论值和实际测量值。从空间分辨力上可以看出,仿真结果与实测结果比较吻合,说明设计的高频时反镜达到设计目标。
表1为高频时反镜性能参数。
环境 | 频率 | 阵长 | 阵元数 | 间距 |
实验室波导环境 | 20kHz | 1.24m | 32 | 0.04m |
表1
表2为中心频率为20kHz,带宽为2.4kHz时仿真和实测时反处理垂直、水平分辨力。
仿真垂直分辨力 | 实测垂直分辨力 | 仿真水平分辨力 | 实测水平分辨力 |
0.09m | 0.07m | 0.49m | 0.41m |
表2
基于时反发射聚焦反波束形成的混响抑制湖上实验:在湖上实验中,通过比较BS(broad side,BS)和TR(time reversal,TR)两种发射情况下的混响和到达目标处的声能量,说明时反发射聚焦反波束形成技术可以抑制混响、提高回混比。BS发射,即利用TRA阵在正横方向发射一个峰值归一化的信号;时反发射聚焦反波束形成:利用PS声源发射一个信号,用TRA阵接收,时反归一化后重新发射,要求其归一化后最大点的功率与BS发射时最大点的功率相同。
目标为三个直径为21cm、长度为51cm的圆柱桶并排组成,PS系在目标中间,位于水下9m,与TRA相距20m。发射信号形式为5ms、10-15kHz的LFM信号。图5为两种情况下的混响衰减曲线。从图上可以看出,TR发射产生的总混响级比BS发射产生的混响约低3-5dB,这是由于时反发射的功率比BS发射的功率小。图6为目标附近录制的时间波形,从图上可以看出,TR发射照射到目标的能量约比BS发射高出6dB。显然,若在目标处入射相同的能量,TR发射只需要较少的发射功率,因此,它所产生的混响也小。
基于高频时反镜的弱目标探测湖上实验:在BS发射和TR发射两种情况下,分别进行多路径补偿波束形成,将它们的输出与发射信号的拷贝做相关,对目标进行检测,再根据拷贝时延对目标距离进行估计。发射信号形式为5ms、10-15kHz的LFM信号,目标为前面三个圆筒的组合体,与TRA相距20m。图7和图8分别给出了BS发射和TR发射两种情况下对目标进行估距的结果。从图中可以看出,在两种情况下,对目标都能够进行有效的检测与估距。显然,在TR发射情况下,利用时反聚焦实现回混比增强后,其对目标的探测性能更具有优势。
Claims (2)
1.一种基于高频时反阵的时反发射聚焦反波束形成方法,其特征在于:该波束形成方法采用高频时反阵,通过发射声波照射目标,再接收目标的回波信号,将接收到的回波信号在时间上进行翻转然后再发射出去,时反信号会在目标处聚焦,表现为时间波形压缩和空间聚焦,再次接收增强的目标回波信号,形成接收聚焦波束。
2.如权利要求1所述的基于高频时反阵的时反发射聚焦反波束形成方法,其特征在于:在波导环境下,基于简正波模型,在频域的波动方程,即Helmholtz方程为(2):
[▽2-k2]G(r,z)=0. (2)
式中:G(r,z)为格林函数,r为声源到测量处的距离,k2=ω2/c2(z)为波数,ω为信号源角频率,c(z)为声速梯度,利用分离变量法,设G(r,z)=Φ(r)Ψ(z),将其代入式(2)经整理可得到深度方程(3)和距离方程(4):
式中:kr、kz分别为波数的水平分量和垂直分量,它们满足:
式(3)是经典的Strum-Liouville特征值问题;式(4)是一阶Bessel方程,其解为零阶Hankel函数;假设激励声源位于水下zs处,在远场条件下,忽略时间因子e-jωt,它产生的声场可近似为(6):
式中:ρ为介质密度,krm为第m号水平波数;
时反处理系统由垂直线阵VLA、时反阵TRA和探察源PS组成,TRA由J个阵元构成,它将接收的信号按先到后发、后到先发的顺序完成信号的时间反转发射,VLA用于监测时反发射信号在声场中的聚焦特性,PS位于水下zps处,与TRA的距离为R;
设PS发射信号为s(t),从PS到TRA中第j(j=1,2,…J)个阵元处的信道格林函数可表示为G(R,zj,ω),则第j个阵元接收的声压场在频域表示为(7):
Pj(ω)=S(ω)G(R,zj,ω). (7)
式中: 为信号的频谱,TRA将接收到的信号进行时间反转发射,对应于频域内进行共轭;因而时反发射信号表示为(8):
时反发射后在观测点r的声压场在频域表示为(9):
它是N个阵元产生的时反声场的综合效果,其对应的时域表示为(10):
参照式(6),整理式(9)得(11):
TRA垂直布满整个波导且充分采样,利用模深度函数的正交性质得到(12):
取n=m并在j上积分得到(13):
当r=R时,
式(14)对于能在波导中传播的有效简正模,krm近似为常数,同时,
最后,式(14)近似为(16):
PTR(R,z,ω)=Aδ(z-zs)S*(ω) (16)
式中:A为一常量,求式(16)的傅立叶反变换得到时域的声压场如下(17):
pTR(R,z,t)=Aδ(z-zs)s(-t). (17)
当z=zs时,有下式(18)成立:
pTR(R,zs,t)=As(-t). (18)
由式(14)和式(18)可见,通过r=R和z=zs实现时反发射信号在声源位置处的空间聚焦,进而完成时反发射聚焦反波束形成。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNA2007101644383A CN101183150A (zh) | 2007-11-30 | 2007-11-30 | 基于高频时反阵的时反发射聚焦反波束形成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNA2007101644383A CN101183150A (zh) | 2007-11-30 | 2007-11-30 | 基于高频时反阵的时反发射聚焦反波束形成方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101183150A true CN101183150A (zh) | 2008-05-21 |
Family
ID=39448489
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNA2007101644383A Pending CN101183150A (zh) | 2007-11-30 | 2007-11-30 | 基于高频时反阵的时反发射聚焦反波束形成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101183150A (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102132537A (zh) * | 2008-06-27 | 2011-07-20 | 法国电信公司 | 估计通过时间反演而预均衡的信号的聚集的质量的方法 |
CN101645264B (zh) * | 2009-07-17 | 2012-03-14 | 中国科学院声学研究所 | 一种基于迭代最小二乘-预滤波的最优主动时反聚焦方法 |
CN101592730B (zh) * | 2009-06-05 | 2012-06-20 | 浙江大学 | 基于参数调节随机共振及后处理的传感器阵波束域微弱信号处理方法 |
CN103164588A (zh) * | 2013-04-15 | 2013-06-19 | 西北工业大学 | 一种大功率聚能反射罩的设计方法 |
CN105301114A (zh) * | 2015-10-20 | 2016-02-03 | 浙江大学 | 一种基于多通道空时逆滤波技术的声学覆盖层插入损失测量方法 |
CN107306168A (zh) * | 2016-04-20 | 2017-10-31 | 中国科学院声学研究所 | 真实时反和虚拟时反相结合的信号发射及处理方法 |
CN108802683A (zh) * | 2018-05-30 | 2018-11-13 | 东南大学 | 一种基于稀疏贝叶斯学习的源定位方法 |
CN110530510A (zh) * | 2019-09-24 | 2019-12-03 | 西北工业大学 | 一种利用线性声阵列波束形成的声源辐射声功率测量方法 |
CN111273294A (zh) * | 2020-02-19 | 2020-06-12 | 西北工业大学 | 一种基于时间反转的水声脉冲回波测距方法 |
CN115079092A (zh) * | 2022-05-06 | 2022-09-20 | 哈尔滨理工大学 | 一种码头环境下水下声源定位方法 |
-
2007
- 2007-11-30 CN CNA2007101644383A patent/CN101183150A/zh active Pending
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102132537A (zh) * | 2008-06-27 | 2011-07-20 | 法国电信公司 | 估计通过时间反演而预均衡的信号的聚集的质量的方法 |
CN102132537B (zh) * | 2008-06-27 | 2014-10-22 | 法国电信公司 | 估计通过时间反演而预均衡的信号的聚集的质量的方法 |
CN101592730B (zh) * | 2009-06-05 | 2012-06-20 | 浙江大学 | 基于参数调节随机共振及后处理的传感器阵波束域微弱信号处理方法 |
CN101645264B (zh) * | 2009-07-17 | 2012-03-14 | 中国科学院声学研究所 | 一种基于迭代最小二乘-预滤波的最优主动时反聚焦方法 |
CN103164588A (zh) * | 2013-04-15 | 2013-06-19 | 西北工业大学 | 一种大功率聚能反射罩的设计方法 |
CN103164588B (zh) * | 2013-04-15 | 2015-10-14 | 西北工业大学 | 一种大功率聚能反射罩的设计方法 |
CN105301114A (zh) * | 2015-10-20 | 2016-02-03 | 浙江大学 | 一种基于多通道空时逆滤波技术的声学覆盖层插入损失测量方法 |
CN105301114B (zh) * | 2015-10-20 | 2018-01-12 | 浙江大学 | 一种基于多通道空时逆滤波技术的声学覆盖层插入损失测量方法 |
CN107306168A (zh) * | 2016-04-20 | 2017-10-31 | 中国科学院声学研究所 | 真实时反和虚拟时反相结合的信号发射及处理方法 |
CN107306168B (zh) * | 2016-04-20 | 2019-10-22 | 中国科学院声学研究所 | 真实时反和虚拟时反相结合的信号发射及处理方法 |
CN108802683A (zh) * | 2018-05-30 | 2018-11-13 | 东南大学 | 一种基于稀疏贝叶斯学习的源定位方法 |
CN110530510A (zh) * | 2019-09-24 | 2019-12-03 | 西北工业大学 | 一种利用线性声阵列波束形成的声源辐射声功率测量方法 |
CN111273294A (zh) * | 2020-02-19 | 2020-06-12 | 西北工业大学 | 一种基于时间反转的水声脉冲回波测距方法 |
CN115079092A (zh) * | 2022-05-06 | 2022-09-20 | 哈尔滨理工大学 | 一种码头环境下水下声源定位方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101183150A (zh) | 基于高频时反阵的时反发射聚焦反波束形成方法 | |
CN102539530B (zh) | 基于时间反转聚焦的水声无源材料回声降低/反射系数测量方法 | |
CN102183435B (zh) | 一种基于多路径反射理论的海底密度和声速测量方法 | |
CN103048642B (zh) | 基于频域最小二乘法的水声脉冲信号匹配场定位方法 | |
CN102590349B (zh) | 基于时间反转聚焦的水声无源材料插入损失/透射系数测量方法 | |
CN108828522A (zh) | 一种利用垂直阵lcmv波束形成的水下目标辐射噪声测量方法 | |
CN105301114B (zh) | 一种基于多通道空时逆滤波技术的声学覆盖层插入损失测量方法 | |
CN103076594B (zh) | 一种基于互相关的水声脉冲信号双阵元定位的方法 | |
CN101387701A (zh) | 一种基于前向预测的被动时间反转混响抑制方法 | |
CN103091677B (zh) | 一种基于时反的均匀线列阵波束形成方法 | |
CN102608609A (zh) | 一种基于mimo的高分辨水下目标探测装置及方法 | |
CN107656243A (zh) | 非均匀媒质中联合doa/toa海洋多径环境定位方法 | |
CN108107436A (zh) | 一种基于可靠声路径的水下目标主动分类与定位方法 | |
CN106443792B (zh) | 一种水下多缆声学网络定位方法 | |
CN108519602A (zh) | 海底冷泉快速成像系统及其阵列信号处理方法 | |
CN110058245A (zh) | 基于云模型的低频主动拖线阵声呐浅海探测效能评估方法 | |
CN107505596A (zh) | 基于双扩展水声信道环境下的mimo主动探测信号设计与检测系统和方法 | |
CN104180891A (zh) | 一种基于声阵的海洋声传播损失实时测量方法 | |
CN106324278A (zh) | 一种基于声参量阵的风速测量方法 | |
CN103926581B (zh) | 一种声纳目标回波亮点参数的测量方法 | |
CN109001297A (zh) | 基于单矢量水听器的大样本水声材料声反射系数测量方法 | |
CN104280737B (zh) | 加权宽带时反算子分解声成像方法 | |
CN103713289A (zh) | 基于分布式Phased-MIMO联合处理的目标探测方法 | |
CN103926586B (zh) | 一种使用发射子阵的mimo阵列测深方法 | |
CN103513249B (zh) | 一种宽带相干模基信号处理方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C12 | Rejection of a patent application after its publication | ||
RJ01 | Rejection of invention patent application after publication |
Open date: 20080521 |