CN103575377A - 测量参量声场中差频波空间分布特性的方法 - Google Patents

测量参量声场中差频波空间分布特性的方法 Download PDF

Info

Publication number
CN103575377A
CN103575377A CN201310553700.9A CN201310553700A CN103575377A CN 103575377 A CN103575377 A CN 103575377A CN 201310553700 A CN201310553700 A CN 201310553700A CN 103575377 A CN103575377 A CN 103575377A
Authority
CN
China
Prior art keywords
frequency
sound
wave
parameter
difference
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.)
Granted
Application number
CN201310553700.9A
Other languages
English (en)
Other versions
CN103575377B (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.)
Harbin Engineering University
Original Assignee
Harbin Engineering University
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 Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN201310553700.9A priority Critical patent/CN103575377B/zh
Publication of CN103575377A publication Critical patent/CN103575377A/zh
Application granted granted Critical
Publication of CN103575377B publication Critical patent/CN103575377B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明提供的是一种测量参量声场中差频波空间分布特性的方法。首先通过两列高频声波f1和f2组成参量声场,测得参量声场的作用区域并将参量声场均匀划分为n个测量网格,在任一网格点A处测得差频波幅值
Figure DDA0000411498520000011
再单独发射低频声波f3,使其在A点的幅值与差频波幅值相等;最后同时发射三列声波f1、f2和f3,得到功率谱值
Figure DDA0000411498520000012
根据功率谱最小值
Figure DDA0000411498520000013
对应的f1初始相位可以得到该点位置处差频波的相位;重复以上步骤测得其它点的差频波幅度以及相位信息。本发明同时描绘出了差频波幅度以及相位在参量声场的分布特性,解决了现有的实验方法以及工程技术难以描述水下参量声场中差频波相位变化规律的问题,为分析声参量场的特性奠定了基础。

Description

测量参量声场中差频波空间分布特性的方法
技术领域
本发明涉及的是非线性声学领域的一种利用三波声互作用测量参量声场中差频波空间分布特性的实验方法。
背景技术
两列大振幅声波在介质中传播时会产生其它频率的声波,包括二倍频、和频以及差频等等,其中以差频波的应用最为广泛。众所周知,声参量阵是非线性声学领域中的重要应用之一,近年来在水声工程和超声医学领域取得了较大的进展,并被广泛应用到水声通讯、鱼群探测、隐蔽水雷探测以及超声成像等方面。它是指利用声波互作用生成差频波以实现在小的基阵尺寸条件下获得低频波束,由于指向性尖锐、频带宽,故能提高空间分辨率、抗混响并可以获得较高的信号处理增益,因此研究参量声场中差频波的空间分布特性至关重要。
文献“不同波束形成的参量阵近场比较,声学技术,2000,vol.2,p81-p83”公开了一种计算声参量场中差频波幅度分布规律的方法。该方法是通过对KZK方程进行数值求解,得到了均匀平面阵、聚焦阵及非衍射阵形成的参量阵波束近场声场,然后比较了不同波束声轴上的差频声压级的分布规律。文献“利用反射的差频声波测量媒质非线性参量B/A,科学通报,2003,vol.21,p27-p28”给出了一种两波非线性互作用产成的差频波幅度的计算方法,其中考虑了媒质对声波的吸收作用。但以上文献只给出了计算差频波幅度的方法,并未涉及到差频波相位的提取方法。
发明内容
本发明的目的在于提供一种可实现对差频波相位特征参数提取的利用三波声互作用测量参量声场中差频波空间分布特性的方法。
本发明的利用三波声互作用测量参量声场中差频波空间分布特性的方法为:
(1)利用高频发射换能器同时发射两个频率分别为f1和f2的声波以组成参量声场,通过两波生成的差频波声压级随距离的变化曲线得到参量声场的作用区域,并均匀划分测量网格;
(2)利用低频发射换能器发射频率为f3的声波,频率关系应满足f3=f2-f1,并保证低频换能器到高频换能器的距离满足远场条件L>>d2/λ,其中,L为低频发射换能器与高频发射换能器的距离,λ为低频声波的波长,d为低频发射换能器辐射面的直径;
(3)在任意一个网格点A处利用声压水听器测得同时发射f1和f2时的差频波幅值
Figure BDA0000411498500000011
然后单独发射低频声波f3,使其在A点的幅值与差频波幅值相等,并记录对应的发射工况;
(4)同时发射三列声波f1、f2和f3,并保证f1、f2、f3的幅度以及f1、f2初始相位与步骤(3)中的发射工况一致,并以
Figure BDA0000411498500000021
为间隔改变f3的初始相位一周并记录其在网格点A处的功率谱值
Figure BDA0000411498500000022
根据功率谱最小值
Figure BDA0000411498500000023
所对应的f1初始相位得到该点位置处差频波的相位;
(5)测量其它点的差频波幅度以及相位信息,重复步骤(3)~(4),最终得出差频波幅度以及相位在参量声场的分布特性。
本发明的利用三波声互作用测量参量声场中差频波空间分布特性的方法,同时描绘出了差频波幅度以及相位在参量声场的分布特性,解决了现有的实验方法以及工程技术难以描述水下参量声场中差频波相位变化规律的问题,为分析声参量场的特性奠定了基础。
附图说明
图1是本发明的差频波声压级随距离的变化规律图;
图2是本发明的水下测量差频波声场的俯视图;
图3是本发明的系统装置连接示意图;
图4是本发明的同时发射中频和高频声波时在A点处的声压级图;
图5是本发明的单独发射低频声波时在A点处的声压级图;
图6是本发明的改变Δ(x=0)对低频声波能量的影响图;
图7是本发明的低频初始相位为350°时三波合发在A点的声压级图。
具体实施方式
下面以测量高频40kHz和45kHz生成的差频波5kHz空间分布特性为例,结合附图对本发明做更详细地描述。
(1)本次实验在消声水池中进行,使高频发射换能器同时发射两个频率分别为f1和f2的声波以组成参量声场。其中f1=40kHz,f2=45kHz,则生成的差频波fd=5kHz。
通过行车控制水听器以2m/min的速度匀速运动,采集5min的数据并处理得出差频波声压级随距离变化的曲线如图1所示。由于水池的背景噪声在100dB附近,因此判断出参量声场截止的位置为6m处,并根据已知高频发射换能器的辐射角可以得出参量声场的作用区域。按照具体的工程要求选择一定的步长划分测量网格,设网格的个数为n。图2为水下差频波声场的俯视图。A为任意一个网格点。在水介质中,差频波声压随距离的衰减规律与线性声学不同:在S1区域内,差频波不断生成,与此同时吸收也逐渐增加,但是差频波随距离增加的效应超过介质对它的吸收效应,此时认为非线性作用大于声吸收作用,参量声场的起始阶段差频波幅值随距离线性增加。在S2区域内,由于中频声波和高频声波能量被吸收,生成的差频波越来越小,此时非线性作用小于声吸收作用,差频波幅值随距离增加而减小。在S3区域内,差频波幅值已经接近背景噪声,视为声参量场的截止区。
(2)利用低频发射换能器发射频率为f3的声波,使频率关系满足f3=f2-f1=fd=5kHz。
系统装置连接示意图如图3所示。在布放收发系统的过程中,应保证低频与高频发射换能器的声中心处于同一水平直线上;将水听器与接收支架刚性连接,使连接的实验仪器共地。低频发射换能器与高频发射换能器的距离L应满足下列关系:L>>d2/λ,其中,λ为低频声波的波长,d为低频发射换能器辐射面的直径。
(3)同时发射f1和f2,在任意一个网格点A处利用声压水听器测得此位置的差频波幅值A∈n。然后单独发射低频声波f3,通过调节发射端的信号源和功放使其在网格A点的幅值与测得的差频波幅值
Figure BDA0000411498500000036
相等。图4为同时发射40kHz和45kHz的声波,生成声压级为128.5dB的差频波。图5为单发5kHz低频声波时,在A点处产生的128.6dB声压级。并记录对应的发射工况。
(4)同时发射三列声波f1、f2和f3,保证f1、f2、f3的幅度以及f1、f2初始相位与步骤(3)中的发射工况一致。由于三波的相位关系对实验测量过程十分重要,因此必须使用同步的信号源发射三列声波。然后以为间隔改变f3的初始相位一周,并通过数据采集仪实时观察功率谱读值,记录其在网格点A的功率谱值m=1,2…表示初始相位的序号。
表1改变初始相位时5kHz声波的功率谱值
初始相位(°) 0 20 40 60 80
功率谱值(dB) -66.4 -59.7 -56.1 -53.7 -52.1
初始相位(°) 100 120 140 160 180
功率谱值(dB) -51 -50.3 -49.7 -49.5 -49.5
初始相位(°) 200 220 240 260 280
功率谱值(dB) -49.8 -50.5 -51.3 -52.9 -54.7
初始相位(°) 300 320 340
功率谱值(dB) -57.4 -62.2 -73.5
已知中频40kHz和高频45kHz合发时产生的差频波5kHz功率谱读值为-55.6dB。表1中的功率谱值表明在某些初始相位条件下,三波互作用后5kHz声波的功率谱值会大于或者小于差频波值。这是由于低频声波叠加到了差频声波上,若低频声波相位与差频波同向则会实现差频波的放大;若低频声波相位与差频波相反则会降低差频波的幅度。因此可以通过功率谱最大值
Figure BDA0000411498500000033
或最小值
Figure BDA0000411498500000034
对应的低频声波相位得出差频波的相位信息。下面通过对三波互作用方程进行数值计算得出更容易获得准确差频波相位信息的方式。
dB 1 ( x ) dx + α 1 B 1 ( x ) = - ω 1 β c 0 2 B 3 ( x ) B 2 ( x ) sin Δ ( x ) dB 2 ( x ) dx + α 2 B 2 ( x ) = ω 2 β c 0 2 B 3 ( x ) B 1 ( x ) sin Δ ( x ) dB 3 ( x ) dx + α 3 B 3 ( x ) = - ω 3 β c 0 2 B 1 ( x ) B 2 ( x ) sin Δ ( x ) dΔ ( x ) dx + β c 0 2 ( ω 3 B 1 ( x ) B 2 ( x ) B 3 + ω 1 B 3 ( x ) B 2 ( x ) B 1 - ω 2 B 3 ( x ) B 1 ( x ) B 2 ) cos Δ ( x ) = 0
其中,B1(x)、B2(x)、B3(x)分别为中频声波、高频声波以及低频声波随距离x变化的振速幅值;α1、α2以及α3分别为对应各列声波在水中的声吸收系数;ω1、ω2和ω3为对应声波的角频率;β为纯水介质的非线性参量;c0为介质中的声速;Δ(x)=s2(x)-s1(x)-s3(x)为三波的相位差,s1(x)、s2(x)、s3(x)分别为对应声波的相位。
根据实验测量选择数值计算的声源参数为:中频声波频率f1=40kHz,声源级spl1=190dB;高频声波频率f2=45kHz,声源级spl2=190dB;低频声波频率f3=5kHz,其声源级为spl3=120dB;初始相位差Δ(x=0)变化范围为0~π。图6为改变Δ(x=0)时低频声波能量ΔSPL1随距离s的变化量。可以得出
Figure BDA0000411498500000042
对应的相位区域较窄且尖锐,因此更容易找到准确的相位值点;而对应的相位区域较宽,若想找到较为准确的相位点需要更大的工作量。
根据以上规律,进而寻找5kHz能量下降最大的点。在表1中找出下降最大的点对应的角度
Figure BDA0000411498500000044
在这个相位点附近进行以10°为间隔的扫描(根据工程需要可以适当的调整扫描间隔),得到330°时的功率谱值为-65.5dB,350°的功率谱值为-88dB。由于350°的功率谱值与背景噪声相近,因此可以认为低频5kHz的声波相位与在A点处的差频波相位完全相反。图7为低频声波的初始相位为350°时的三波声压级图,分析其结果与实时观测一致。
根据球面波扩展规律,计算得到的f3的初始相位传播到网格点A时的相位值:
Figure BDA0000411498500000045
其中
Figure BDA0000411498500000046
为低频初始相位值,k3=2πf3/c0为低频的波数,xA为等效声中心到A点的距离,
Figure BDA0000411498500000047
为传播至A点低频声波的相位。
根据以上的测量结果得到
Figure BDA0000411498500000048
xA=0.7m,k3=20.94,代入式
Figure BDA0000411498500000049
中计算得出
Figure BDA00004114985000000410
取反向即为声参量场中该点差频波的相位信息,得到A点差频波的相位信息为-364°。

Claims (1)

1.一种测量参量声场中差频波空间分布特性的方法,其特征是:
(1)利用高频发射换能器同时发射两个频率分别为f1和f2的声波以组成参量声场,通过两波生成的差频波声压级随距离的变化曲线得到参量声场的作用区域,并均匀划分测量网格;
(2)利用低频发射换能器发射频率为f3的声波,频率关系应满足f3=f2-f1,并保证低频换能器到高频换能器的距离满足远场条件L>>d2/λ,其中,L为低频发射换能器与高频发射换能器的距离,λ为低频声波的波长,d为低频发射换能器辐射面的直径;
(3)在任意一个网格点A处利用声压水听器测得同时发射f1和f2时的差频波幅值
Figure FDA0000411498490000014
然后单独发射低频声波f3,使其在A点的幅值与差频波幅值相等,并记录对应的发射工况;
(4)同时发射三列声波f1、f2和f3,并保证f1、f2、f3的幅度以及f1、f2初始相位与步骤(3)中的发射工况一致,并以
Figure FDA0000411498490000011
为间隔改变f3的初始相位一周并记录其在网格点A处的功率谱值
Figure FDA0000411498490000012
根据功率谱最小值
Figure FDA0000411498490000013
所对应的f1初始相位得到该点位置处差频波的相位;
(5)测量其它点的差频波幅度以及相位信息,重复步骤(3)~(4),最终得出差频波幅度以及相位在参量声场的分布特性。
CN201310553700.9A 2013-11-11 2013-11-11 测量参量声场中差频波空间分布特性的方法 Active CN103575377B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310553700.9A CN103575377B (zh) 2013-11-11 2013-11-11 测量参量声场中差频波空间分布特性的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310553700.9A CN103575377B (zh) 2013-11-11 2013-11-11 测量参量声场中差频波空间分布特性的方法

Publications (2)

Publication Number Publication Date
CN103575377A true CN103575377A (zh) 2014-02-12
CN103575377B CN103575377B (zh) 2015-09-30

Family

ID=50047652

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310553700.9A Active CN103575377B (zh) 2013-11-11 2013-11-11 测量参量声场中差频波空间分布特性的方法

Country Status (1)

Country Link
CN (1) CN103575377B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104215964A (zh) * 2014-08-20 2014-12-17 哈尔滨工程大学 一种多列等差频率原波相互作用形成参量阵的声场获取方法
CN104501938A (zh) * 2014-11-19 2015-04-08 哈尔滨工程大学 一种在置于空气中的矩形混响水槽内测量水下声源低频辐射声功率的方法
CN104501939A (zh) * 2014-11-19 2015-04-08 哈尔滨工程大学 一种利用单水听器测量非消声水池混响时间的反演方法
CN106768303A (zh) * 2017-02-08 2017-05-31 哈尔滨理工大学 两列平面声波非线相互作用实验测试系统及方法
CN112198515A (zh) * 2020-10-13 2021-01-08 湖南国天电子科技有限公司 一种参量阵浅剖差频转换性能优化方法

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104215964A (zh) * 2014-08-20 2014-12-17 哈尔滨工程大学 一种多列等差频率原波相互作用形成参量阵的声场获取方法
CN104215964B (zh) * 2014-08-20 2017-05-24 哈尔滨工程大学 一种多列等差频率原波相互作用形成参量阵的声场获取方法
CN104501938A (zh) * 2014-11-19 2015-04-08 哈尔滨工程大学 一种在置于空气中的矩形混响水槽内测量水下声源低频辐射声功率的方法
CN104501939A (zh) * 2014-11-19 2015-04-08 哈尔滨工程大学 一种利用单水听器测量非消声水池混响时间的反演方法
CN104501939B (zh) * 2014-11-19 2017-05-24 哈尔滨工程大学 一种利用单水听器测量非消声水池混响时间的反演方法
CN106768303A (zh) * 2017-02-08 2017-05-31 哈尔滨理工大学 两列平面声波非线相互作用实验测试系统及方法
CN106768303B (zh) * 2017-02-08 2019-10-01 哈尔滨理工大学 两列平面声波非线性相互作用实验测试系统及方法
CN112198515A (zh) * 2020-10-13 2021-01-08 湖南国天电子科技有限公司 一种参量阵浅剖差频转换性能优化方法
CN112198515B (zh) * 2020-10-13 2021-06-29 湖南国天电子科技有限公司 一种参量阵浅剖差频转换性能优化方法
US11237258B1 (en) 2020-10-13 2022-02-01 Hunan Guotian Electronic Technology Co., Ltd. Method for optimization of a parametric array shallow profile difference frequency conversion performance

Also Published As

Publication number Publication date
CN103575377B (zh) 2015-09-30

Similar Documents

Publication Publication Date Title
CN106501795B (zh) 一种利用混响水池进行水声换能器互易校准的方法
CN103575377B (zh) 测量参量声场中差频波空间分布特性的方法
US11391863B2 (en) Method of free-field broadband calibration of hydrophone sensitivity based on pink noise
CN102590349B (zh) 基于时间反转聚焦的水声无源材料插入损失/透射系数测量方法
CN107272005B (zh) 基于可靠声路径下目标回波到达时延和到达角度的主动定位方法
CN104501938B (zh) 一种在置于空气中的矩形混响水槽内测量水下声源低频辐射声功率的方法
CN103776524A (zh) 一种测量强聚焦超声非线性声场分布的系统及其方法
CN113868583B (zh) 一种子阵波束聚焦的声源距离计算方法及系统
CN102879785A (zh) 基于频谱认知和分段跳频调频的水下物体探测方法及系统
Brown et al. A point-based scattering model for the incoherent component of the scattered field
CN110988850B (zh) 一种基于目标散射的换能器指向性测量方法及装置
Huang et al. A Gaussian finite-element method for description of sound diffraction
CN110109125B (zh) 球面聚焦相控参量阵声探测装置
Voloshchenko et al. Experimental study of the transmission of low-frequency acoustic waves through a water–air interface
CN114446277A (zh) 一种全向宽带声聚能装置
CN106768264B (zh) 基于焦域轴向最大振速的聚焦超声功率测量系统及方法
CN202995039U (zh) 基于频谱认知和分段跳频调频的水下物体探测系统
CN109916497B (zh) 一种在混响水槽测量水下声源甚低频辐射特性的方法
CN101615392B (zh) 一种利用随机阵列实现混响声场聚焦的方法
CN110244305A (zh) 一种水下目标信号散射的仿真方法
CN104215964A (zh) 一种多列等差频率原波相互作用形成参量阵的声场获取方法
Irarrazabal et al. Initial infrasound source characterization using the phase and amplitude gradient estimator method
CN113126029A (zh) 适用于深海可靠声路径环境的多传感器脉冲声源定位方法
CN113624330A (zh) 一种水下目标辐射噪声测量组合体积阵及测量方法
CN113253284A (zh) 基于目标散射特性的主动声呐干涉条纹生成方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant