CN104215964B - 一种多列等差频率原波相互作用形成参量阵的声场获取方法 - Google Patents
一种多列等差频率原波相互作用形成参量阵的声场获取方法 Download PDFInfo
- Publication number
- CN104215964B CN104215964B CN201410409816.XA CN201410409816A CN104215964B CN 104215964 B CN104215964 B CN 104215964B CN 201410409816 A CN201410409816 A CN 201410409816A CN 104215964 B CN104215964 B CN 104215964B
- Authority
- CN
- China
- Prior art keywords
- parametric array
- frequency
- sound field
- sound
- parametric
- 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
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
- G01S7/523—Details of pulse systems
- G01S7/524—Transmitters
-
- 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
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/885—Meteorological systems
-
- 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
- G01S7/534—Details of non-pulse systems
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B13/00—Transmission systems characterised by the medium used for transmission, not provided for in groups H04B3/00 - H04B11/00
- H04B13/02—Transmission systems in which the medium consists of the earth or a large mass of water thereon, e.g. earth telegraphy
Abstract
本发明属于适用于参量阵工程应用领域,具体涉及一种多列等差频率原波相互作用形成参量阵的声场获取方法。本发明包括:读取多列高频原波组合而成声波波形;宽带高频换能器将预设的特定声波信号波形辐射出去,通过传播媒质的非线性解调作用,最终生成具有低频、高指向性的系列差频波束;设定参量阵声场的计算区域,并进行网格化离散;利用有限差分法获取多列声波相互作用条件下的参量阵声场。该方法提高了计算精度,更为准确的计算了参量阵辐射系统声场的分布特点,形象直观地显示出参量阵的声场分布情况,更加准确、全面地反映声场性能。
Description
技术领域
本发明属于适用于参量阵工程应用领域,具体涉及一种多列等差频率原波相互作用形成参量阵的声场获取方法。
背景技术
中国科学院东海研究站在2002年研制成功了一套“堤防隐患监测声纳”,该声纳是一种新型地质声呐,可以对江河湖底、海底及沉积层进行探测识别或对堤防的损毁程度进行探测和评估。原波频率为100kHz,差频为3kHz、6kHz、8kHz、10kHz、12kHz、15kHz、24kHz等七档,并且能根据工作需求进行档位切换。Lucilla Di Marcoberardino在文献“Nonlinearmulti-frequency transmitter for sea-floor characterization”认为频率多样性可以提供更为丰富的海底地层分布信息,并提出一种能同时生成多种谐波频率的声源模型,对沉没的邮轮、海底的地质特性、海底管道或电缆的铺设提供技术指导。John A.Birken在文献“Empirical results from frequency-scanning nonlinear sonar in deep water”中,同样基于海底地层的探测需求,对多列声波同时辐射形成系列差频声场进行了理论及实验研究。但其进一步的相关研究并未见诸于文献。
上述参考文献表明,改变参量阵的发射波形,充分利用由于非线性效应生成的二阶声波(如和频、倍频、差频等)能有效获取水下海底地层的声场信息。同时,理论研究发现,参量阵发射波形的改变可以对参量阵原波转换效率产生影响,参量阵辐射系统初始声波信号的改变可以有效提高参量阵的转换效率。比如,在换能器带宽允许的情况下,Merklinger通过理论推导得出如果峰值发射功率相同,改变辐射波形可在一定程度上提高参量阵声压级。基于此,本发明利用多列原波相互作用的激励模式,为参量阵的辐射系统提供一种新颖的思路,在有效提高参量阵的转换效率的同时,能生成系列低频、高指向性的差频波束,可为水声通讯、海底地质探测等领域的信号检测提高更加丰富的声场信息。
发明内容
本发明的目的在于提供一种提高参量阵的转换效率,还可利用生成的系列差频信号进行多频点信号检测和水声通讯的多列等差频率原波相互作用形成参量阵的声场获取方法。
本发明的目的是这样实现的:
(1)读取多列高频原波组合而成声波波形:其中,原波频率满足等差关系式f2-f1=f3-f2=…=fn-fn-1,通过信号发射系统向宽带高频换能器提供初始条件,任意信号发射器根据参量阵辐射系统的信号波形情况,生成指定的参量阵辐射信号;宽带功率放大器对参量阵辐射信号进行预处理,通过功率放大器的放大效应,使得有限振幅声波满足参量阵生成条件;宽带高频换能器将预设的特定声波信号波形辐射出去,通过传播媒质的非线性解调作用,最终生成具有低频、高指向性的系列差频波束;
(2)设定参量阵声场的计算区域,并进行网格化离散;利用有限差分法获取多列声波相互作用条件下的参量阵声场。
本发明的有益效果在于:
本发明将二阶龙格-库塔(Runge-Kutta)法与CNFD方法相结合,对多列等差频率原波相互作用形成参量阵声场进行推演。该方法一方面提高了计算精度,更为准确的计算了参量阵辐射系统声场的分布特点,形象直观地显示出参量阵的声场分布情况,更加准确、全面地反映声场性能。另一方面,该辐射系统能有效提高参量阵的转换效率,并且系列低频、高指向性的差频波束可为水声通讯信号检测提供更加丰富的声场信息。
附图说明
图1多频参量阵声场获取方法;
图2多频参量阵单通道辐射系统;
图3多频参量阵的原波激励及差频波生成图示;
图4五列原波相互作用的差频波声压幅值分布图示;(a)差频1;(b)差频2;(c)差频3;(d)差频4;(e)径向声压幅值比较;(f)轴向声压幅值比较;
图5四列原波相互作用的差频波声压幅值分布图示;(a)差频1;(b)差频2;(c)差频3;(d)径向声压幅值比较;(e)轴向声压幅值比较;
图6三列原波相互作用的差频波声压幅值分布图示;(a)差频1;(b)差频2;(c)径向声压幅值比较;(d)轴向声压幅值比较;
图7两列原波相互作用的差频波声压幅值分布图示;(a)差频1;(b)径向声压幅值比较;(c)轴向声压幅值比较;
图8多列原波相互作用形成差频波fd的声压幅值特性;(a)轴向;(b)径向;
图9多列原波相互作用形成差频波f2d的声压幅值特性;(a)轴向;(b)径向;
图10多列原波相互作用形成差频波f3d的声压幅值特性;(a)轴向;(b)径向。
具体实施方式
结合附图和实例对本发明进一步说明。
本发明提供了一种利用多列等差频率原波相互作用形成参量阵声场的获取方法。该发明将二阶龙格-库塔(Runge-Kutta)法引入计算参量阵声场分布的KZK理论模型,用于计算参量阵的近场声场。该法在有效提高参量阵的转换效率的同时,可以生成系列低频、高指向性的差频波束,能为水声通讯、海底地质探测等领域的信号检测提供更加丰富的声场信息。
本发明解决其技术问题所采用的方案包括以下步骤:
(a)根据参量阵发射阵元的形状确定参量阵近场声场的有限计算区域,采用和参量阵发射阵元边界拟合程度最佳的网格对计算区域分割;
(b)读取参量阵赖以传播媒介的声速、密度、非线性系数等物理参数;
读取参量阵发射系统的激励信号:多列高频原波组合而成声波波形。其中,所述原波频率满足等差关系式f2-f1=f3-f2=…=fn-fn-1,通过信号发射系统向宽带高频换能器提供初始条件,所述信号发射系统包括任意信号发生器、功率放大器、宽带高频换能器等部件;
(c)步骤(b)所述的多列原波相互作用形成的参量阵辐射系统为参量阵单通道辐射系统,所述系统包含三个系统:(c1)任意信号发射器:根据参量阵辐射系统的信号波形情况,生成指定的参量阵辐射信号;(c2)宽带功率放大器:对参量阵辐射信号进行预处理,通过功率放大器的放大效应,使得有限振幅声波满足参量阵生成条件;(c3)宽带高频换能器:将预设的特定声波信号波形辐射出去,通过传播媒质的非线性解调作用,最终生成具有低频、高指向性的系列差频波束。其中所述参量阵单通道辐射系统的任意信号发射系统、宽带功率放大系统及其宽带高频换能器系统的连接方式为线性连接;
(d)根据步骤(a)设定的离散网格模型以及适用于参量阵声场获取方法的有限振幅声波传播理论,利用二阶龙格库塔法和Crank-Ni col son法相结合的有限差分法获取多列声波相互作用条件下的参量阵声场。
实施例
(a)以圆形活塞声源形成的参量阵辐射系统为例,基于活塞声源的轴对称特性,化三维计算区域到二维roz平面,建立参量阵声场的有限计算区域模型,并进行离散化网格划分。沿z方向将0<z<zmax区间分成M段,轴向坐标下标m从1变化到M,把m=1称为第1层,且每层沿径向即r方向的径向坐标下标j从1变化到J,即每层都要计算J个声压值。
(b)读取参量阵赖以传播媒介的声速、密度、非线性系数等物理参数,读取参量阵发射系统的激励信号p(rj,z1,k)及谐波分量g(rj,z1,k)、h(rj,z1,k):
其中:其中:1≤j≤J-,闭区间[1,J-]表示活塞声源所处的径向位置,Kn表示n列原波相互作用,K2-K1=K3-K2=…=Kn-Kn-1;
(c)将初始条件p(rj,z1,k)及其谐波分量g(rj,z1,k)、h(rj,z1,k)输入反映参量阵声场传播非线性效应的守恒迎风格式计算模块,获得传播轴向步长dz后第二层网格z2处的声压pnon(rj,z2,k)及其谐波幅值分量gnon(rj,z2,k)、hnon(rj,z2,k);
(d)将z2层上的声压幅值pnon(rj,z2,k)及其分量gnon(rj,z2,k)、hnon(rj,z2,k)作为第一层网格的虚拟声源激励信号,输入描述参量阵声场传播衍射、吸收效应的DIRK+CNFD计算模块,获取轴向传播dz之后z2处的声压幅值:pabs,dif,non(rj,z2,k);
(e)重复步骤(c)-(d),按照z2→z3,z3→z4,…,zM-1→zM进行轴向逐步推进,依次循环计算获取参量阵声场分布的守恒迎风格式计算模块、DIRK+CNFD计算模块,即可获得多列原波相互作用激励条件下的参量阵声场。
实例参数设置如下:设换能器的半径为a=10cm,径向计算区域为(0,rmax),其中rmax=41a,换能器的轴向计算区域为(0,zmax)其中zmax=2.5d,d=πf0a2/c为换能器辐射原波中心频率对应的瑞利距离。为了减小边界反射对计算声场的干涉,设定区域[40a,41a]为PML区域,假定轴向每个瑞利距离划分为120等分,径向每个单位半径长度划分为30等分,可知计算区域的网格数为1230*300,径向间隔Δr=3.336mm,轴向步进长度dz=9.748mm,分别设置不同的原波初始条件:两列原波f1=53kHz、f2=57kHz,归一化声压幅值三列原波f1=51kHz、f2=55kHz、f3=59kHz;归一化声压幅值四列原波f1=49kHz、f2=53kHz、f3=57kHz、f4=61kHz,归一化声压幅值五列原波f1=47kHz、f2=51kHz、f3=55kHz、f4=59kHz、f5=63kHz,归一化声压幅值显然,上面不同列数原波形成的差频波是多样的:fd=4kHz、f2d=8kHz、f3d=12kHz、f4d=16kHz。下面组图给出不同原波列数形成的参量阵声场空间分布图示
由图4-7可知,n列原波的频率呈现等差形式时(如图3所示),由于传播介质的相互作用形成包含n-1(n≥2)个频率成分的参量阵声场,并且该多频参量阵的各个差频波束均可以实现低频、高指向性辐射,五列原波相互作用时差频波峰值声压幅值关系四列原波相互作用时差频波峰值声压幅值关系三列列原波相互作用时差频波峰值声压幅值关系如果对具有相同辐射功率的多列原波差频声压幅值轴向空间分布图进行分析(图8-10),可知,差频波fd的声压峰值关系为:pn=5>pn=4>pn=3>pn=2(图8);差频波f2d的声压峰值关系为:pn=5>pn=4>pn=3(图9);差频波f3d的声压峰值关系为:pn=5>pn=4(图10)。说明多列原波相互作用一方面形成频点信息更为丰富的差频波,另一方面随着原波列数的增加,参量阵功率转换效率亦同样得到提高。
Claims (1)
1.一种多列等差频率原波相互作用形成参量阵的声场获取方法,其特征在于:
(1)读取多列高频原波,而后组合而成声波波形:其中,原波频率满足等差关系式f2-f1=f3-f2=…=fn-fn-1,通过信号发射系统向宽带高频换能器提供初始条件,任意信号发射器根据参量阵辐射系统的信号波形情况,生成指定的参量阵辐射信号;宽带功率放大器对参量阵辐射信号进行预处理,通过功率放大器的放大效应,使得有限振幅声波满足参量阵生成条件;宽带高频换能器将预设的特定声波信号波形辐射出去,通过传播媒质的非线性解调作用,最终生成具有低频、高指向性的系列差频波束;
(2)设定参量阵声场的计算区域,并进行网格化离散;利用有限差分法获取多列声波相互作用条件下的参量阵声场。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410409816.XA CN104215964B (zh) | 2014-08-20 | 2014-08-20 | 一种多列等差频率原波相互作用形成参量阵的声场获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410409816.XA CN104215964B (zh) | 2014-08-20 | 2014-08-20 | 一种多列等差频率原波相互作用形成参量阵的声场获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104215964A CN104215964A (zh) | 2014-12-17 |
CN104215964B true CN104215964B (zh) | 2017-05-24 |
Family
ID=52097661
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410409816.XA Active CN104215964B (zh) | 2014-08-20 | 2014-08-20 | 一种多列等差频率原波相互作用形成参量阵的声场获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104215964B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106841382B (zh) * | 2017-01-23 | 2019-06-14 | 哈尔滨工程大学 | 基于三波耦合互作用非均匀混合介质非线性系数测量方法 |
CN113536554B (zh) * | 2021-07-01 | 2022-08-02 | 哈尔滨工程大学 | 一种采用压缩等效源法的封闭空间内声场预测方法 |
CN115665633B (zh) * | 2022-12-26 | 2023-03-31 | 中国人民解放军海军工程大学 | 一种参量阵扬声器基波调制的方法、记录媒体及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1300632A (zh) * | 1999-12-22 | 2001-06-27 | 许田 | 参量阵超声波医疗仪 |
CN102636786A (zh) * | 2012-05-17 | 2012-08-15 | 绵阳市浦发电子科技有限公司 | 一种基于参量阵的水下埋藏物探测声纳 |
CN103576574A (zh) * | 2013-10-10 | 2014-02-12 | 哈尔滨工程大学 | 一种三列声波非线性相互作用下声能量转换的相位控制方法 |
CN103575377A (zh) * | 2013-11-11 | 2014-02-12 | 哈尔滨工程大学 | 测量参量声场中差频波空间分布特性的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101081877B1 (ko) * | 2010-04-16 | 2011-11-09 | 국방과학연구소 | 파라메트릭 배열을 이용한 공기중 음향 송수신 장치 및 방법 |
-
2014
- 2014-08-20 CN CN201410409816.XA patent/CN104215964B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1300632A (zh) * | 1999-12-22 | 2001-06-27 | 许田 | 参量阵超声波医疗仪 |
CN102636786A (zh) * | 2012-05-17 | 2012-08-15 | 绵阳市浦发电子科技有限公司 | 一种基于参量阵的水下埋藏物探测声纳 |
CN103576574A (zh) * | 2013-10-10 | 2014-02-12 | 哈尔滨工程大学 | 一种三列声波非线性相互作用下声能量转换的相位控制方法 |
CN103575377A (zh) * | 2013-11-11 | 2014-02-12 | 哈尔滨工程大学 | 测量参量声场中差频波空间分布特性的方法 |
Non-Patent Citations (6)
Title |
---|
NUMERICAL METHODS FOR NONLINEAR WAVE PROPAGATION IN ULTRASOUND;GianmarcoF.Pinton;《Dissertations & Theses》;20071231;全文 * |
Underwater Multi-frequency Transmitter for Seabed Characterization;Lucilla Di Marcoberardino et al.;《OCEANS,2011 IEEE - Spain》;20111231;全文 * |
大振幅波非线性传播的频率特性;杨德森等;《哈尔滨工程大学学报》;20100731;第31卷(第07期);第928-931页 * |
有限声束非线性超声场的快速计算及其特性研究;先永利;《中国优秀硕士学位论文全文数据库 医药卫生科技辑》;20110915(第09期);第E080-45页 * |
水中声波非线性相互作用的声吸收声研究;杨德森等;《振动与冲击》;20121231;第31卷(第08期);第52-56页 * |
生物医学超声中若干非线性问题的研究;杜宏伟;《中国博士学位论文全文数据库 医药卫生科技辑》;20070915(第03期);第E059-3页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104215964A (zh) | 2014-12-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Hovem | Ray trace modeling of underwater sound propagation. Documentation and use of the PlaneRay model | |
CN103076594A (zh) | 一种基于互相关的水声脉冲信号双阵元定位的方法 | |
CN110399680B (zh) | 一种浅海弹性结构辐射声场计算方法 | |
CN104678384B (zh) | 一种波束域的声压差互相关谱分析水下目标速度估计方法 | |
CN104215964B (zh) | 一种多列等差频率原波相互作用形成参量阵的声场获取方法 | |
CN104765037A (zh) | 基于短垂直阵的水下目标定位稳健方法 | |
CN112987004A (zh) | 一种浅海环境下基于水平阵列的水面水下目标分类方法 | |
CN111830493A (zh) | 水中目标收发分置中高频声目标强度预报系统及方法 | |
CN107526105A (zh) | 一种波场模拟交错网格有限差分方法 | |
CN104793212A (zh) | 利用声波海底反射实现主动声纳远程探测的方法 | |
CN104597488A (zh) | 非等边长网格波动方程有限差分模板优化设计方法 | |
Yang et al. | Numerical simulation of acoustic fields in formation generated by linear phased array acoustic transmitters during logging while drilling | |
CN104035128A (zh) | 可控震源伪随机扫描信号生成方法 | |
CN104077479B (zh) | 一种基于守恒迎风格式获取参量阵声场空间分布的方法 | |
CN113866718A (zh) | 一种基于互质阵的匹配场被动定位方法 | |
CN103575377A (zh) | 测量参量声场中差频波空间分布特性的方法 | |
CN105676280A (zh) | 基于旋转交错网格的双相介质地质数据获取方法和装置 | |
CN102830424B (zh) | 一种检波器组合参数计算方法 | |
Qu et al. | Experimental study of a broadband parametric acoustic array for sub-bottom profiling in shallow water | |
CN110285876A (zh) | 一种海洋声场全波解的获取方法 | |
Rutenko et al. | Effect of a sloping bottom on sound propagation | |
CN103336302B (zh) | 基于高次余弦幅度加权的地震波束形成方法 | |
Rutenko | The influence of internal waves on losses during sound propagation on a shelf | |
CN104570131A (zh) | 一种估计大地电磁参数的方法和装置 | |
Petrov et al. | Problem of acoustic diagnostics of a damaged zone |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |