张海澜《理论声学》答案Theoretical Acoustics
海水吸收对远距离水下声场计算的影响
海水吸收对远距离水下声场计算的影响我国计划在“十三五”期间,向地球深处进军,全面实施深地探测、深海探测、深空对地观测战略,在“三深”领域跻身世界先进行列。
其中,深海探测就包括深海声传播的研究。
深海声道存在于全球各深海领域,具有良好的声传播效应,且十分稳定,不受季节影响。
深海声道的一个重要特点就是声速剖面存在一个极小值点,其所在深度被称为“深海声道轴”,声速在声道轴两侧分别呈现正梯度和负梯度,使声波能够沿着声道轴不断传播而不与海底或海面发生接触,从而传播数百甚至数千千米。
深海声道的另外一个特点就是深度足够大,使得海底声速达到或者超过水面声速,使声波可以在与海底接触之前发生反转。
当声波频率很低或传播距离较短时,海水吸收对水下声传播的影响几乎可以忽略。
但是当声波传播上千千米时,海水吸收的影响就不可以忽略了。
大多数常用的声场计算模型包含了海水吸收,例如基于简正波理论的COUPLE模型和基于波数积分理论的SCOOTER模型。
但是基于抛物方程理论的RAM模型忽略了海水吸收这个因素。
本文首先介绍了海水吸收系数的定义与经验公式,通过引入声速虚部来引入海水吸收,并改进了RAM模型,使之可以考虑海水吸收。
接着,考虑理想波导声传播问题,以解析解作为参考解,分别与COUPLE模型和改进后的RAM模型的结果进行对比,三种计算方式得到的结果几乎相同,验证了海水吸收对理想波导中声传播的影响。
然后计算Munk剖面下深海波导中的声场,分别使用SCOOTER模型、COUPLE模型和改进后的RAM模型进行计算,验证了改进后RAM模型的正确性,并展示了海水吸收对深海远距离传播的影响。
声波在海水中传播时,会有一部分声能量被海水吸收,转换成热能。
同时,当海水介质不均匀时,声波还会发生散射。
但是在实际的海洋环境中难以区分吸收和散射带来的影响,所以在计算中它们共同构成了海水吸收。
首先定义平面波吸收系数α[1]:其中A表示平面波幅值的均方根。
由式(1)可以得到A=A0exp(−αx),A0是x=0处的幅值的均方根,如果x的单位是m,α的单位就是Np·m-1(1Np·m-1=8.686dB·m-1)。
拉普拉斯方程从直角坐标系转换到圆坐标系下
2015/5/13
1
例 2.8 一个半径为 ρ0 的薄圆盘,上下两面绝热,内部无热源,圆周边 缘温度分布已知为 F (x, y) .求达到稳恒状态时圆盘内的温度分布。
设圆盘上坐标为 (x, y) 的点处的温度为 u(x, y) ,则上述问题归结为下 述定解问题:
ρ2 ( ∂2u + ∂2u )
=
∂2u
+
∂u ρ
+
ρ2
∂2u
∂x2 ∂y2
∂θ2 ∂θ
∂ρ2
(1.21)
另外,推导出极坐标系下的拉普拉斯方程,还可以用拉梅系数,采用一般正交坐标系下的哈密顿算符直接 运算(可能涉及《微分几何》)。 还可以参考《理论声学》(张海澜)第六章(6.1)对球坐标系梯度算子和散度算子 的推导过程。
∂u θ
∂ρ
+
ρ2
cos2
∂2u θ ∂ρ2
=
2
cos
θ
sin
θ
∂u ∂θ
+
sin2
θ
∂2u ∂θ2
+
∂u ρ
∂ρ
−
ρ
cos2
θ
∂u ∂ρ
−
ρ
cos
θ
sin
θ
∂2u ∂ρ2
−ρ sin θ
∂2u cos θ
∂ρ∂θ
+
ρ2
cos2
∂2u θ
∂ρ2
(1.17)
同样,对于 ∂u 有相同的过程,部分如下:
∂y
∂u ∂u
∂u
ρ = − sin θ + ρ cos x
∂x ∂θ
理论声学 Theoretical Acoustics
音频声学实验室
• 主体建筑建成 于1964、1985 年,国内外少有 • 大量的声场实 验和标准测量 •仪器设备改造 投资逾400万元
2018/10/10
理论声学(1) 绪论
30
Pneumatic loudspeaker
2018/10/10
理论声学(1) 绪论
31
Equipment Made in China
2018/10/10 理论声学(1) 绪论 1
主要参考书
马大猷,现代声学理论基础,科学出版社 (2004)。 杜功焕、朱哲民、龚秀芬著,声学基础(第2版), 南京大学出版社(2001年)。 Allan D. Pierce, Acoustics, McGraw-Hill Book Company P.M.莫尔斯、K.U.英格特著,理论声学(上、下 两册),科学出版社(1984、1986年)。 Dan Russell /~drussell/Demos.html,
2018/10/10
理论声学(1) 绪论
11
振动
• 频率 单位时间里波动的次数 单位 Hz次/秒 赫兹 • 周期 波动一次的时间 • 波长 波动一次传播的距离 • 声速 声速等于波长乘以频率 • 振幅
2018/10/10 理论声学(1) 绪论 12
声音的频谱
A4=440Hz
2018/10/10 理论声学(1) 绪论 13
2018/10/10
理论声学(1) 绪论
8
声波的定性描述
2018/10/10
理论声学(1) 绪论
9
声是弹性介质中传播的机械波
• 弹性介质: 空气 • 真空中没有声波 水 (流体) 固体
2018/10/10
水声学理论习题答案
《水声学》部分习题参考答案绪论1略2略3略4略5环境噪声和海洋混响都是主动声呐的干扰,在实际工作中如何确定哪种干扰是主要的?解:根据水文条件及声呐使用场合,画出回声信号级、混响掩蔽级和噪声掩蔽级随距离变化的曲线,如下图,然后由回声信号曲线与混响掩蔽级、噪声掩蔽级曲线的交点所对应的距离来确定混响是主要干扰,还是噪声为主要干扰,如下图,r R<r n,所以混响是主要干扰。
声信号级噪声掩蔽级R6工作中的主动声呐会受到哪些干扰?若工作频率为1000Hz,且探测沉底目标,则该声呐将会受到哪些干扰源的干扰。
解:工作中的主动声呐受到的干扰是:海洋环境噪声、海洋混响和自噪声,若工作频率为1000Hz,干扰来自:风成噪声、海底混响、螺旋桨引起的自噪声及水动力噪声。
7已知混响是某主动声呐的主要干扰,现将该声呐的声源级增加10dB,问声呐作用距离能提高多少?又,在其余条件不变的情况下,将该声呐发射功率增加一倍,问作用距离如何变化。
(海水吸收不计,声呐工作于开阔水域)解:对于受混响干扰的主动声呐,提高声源级并不能增加作用距离,因为此时信混比并不改变。
在声呐发射声功率增加一倍,其余条件不变的情况下,作用距离变为原距离的42倍,即R R 412 。
第一章 声学基础1 什么条件下发生海底全反射,此时反射系数有什么特点,说明其物理意义。
解:发生全反射的条件是:掠时角小于等于全反射临界角,界面下方介质的声速大于界面上方介质的声速。
发生全反射时,反射系数是复数,其模等于1,虚部和实部的比值给出相位跳变角的正切,即全反射时,会产生相位跳变。
2 略3 略第二章 海洋声学特性1 海水中的声速与哪些因素有关?画出三种常见的海水声速分布。
解:海水中的声速与海水温度、密度和静压力(深度)有关,它们之间的关系难以用解析式表达。
CCCz2 略3 略4 略5 略6 声波在海水中传播时其声强会逐渐减少。
(1)说明原因;(2)解释什么叫物理衰减?什么叫几何衰减?(3)写出海洋中声传播损失的常用TL 表达式,并指明哪项反映的主要是几何衰减,哪项反映的主要是物理衰减;(4)试给出三种不同海洋环境下的几何衰减的TL 表达式。
大连理工大学研究生培养方案
船舶与海洋工程(专业代码:0824授予工学硕士学位)一、培养目标为适应新世纪现代化造船技术的进步与发展、海洋工程装备研发和国防军工事业等对高层次专业人才的需要,船舶与海洋工程学科的硕士研究生培养,要求贯彻德智体全面发展的方针,以夯实理论基础,强化综合素质,提高专业基本技能为出发点,以科学技术研究能力培养为中心,突出硕士生综合素质和创新能力的培养,形成船舶与海洋工程专业的高层次人才培养特色。
本学科硕士学位获得者应掌握本学科的现状、发展方向和国内外学科的前沿发展动态,能较为熟练掌握一门外国语,阅读本专业的外文资料;具备一定的工程设计和实验能力,掌握基本测试技术、数据分析和计算机应用技术;具备本专业扎实的基础理论和系统的专门知识,有独立分析科学问题和解决工程实际问题的能力,达到国家学位条例对本学科硕士学位论文的要求。
硕士研究生毕业后能够从事船舶与海洋结构物设计制造、轮机工程和水声工程等相关领域的科研、设计、生产和管理等工作。
二、学科、专业及研究方向简介船舶与海洋工程学科为国家“211工程”和“985工程”重点建设学科,学科实力雄厚,现有船舶与海洋工程博士后科研流动站、船舶与海洋工程国家一级学科博士点和船舶与海洋结构物设计制造国家重点学科,拥有船舶与海洋结构物设计制造、水声工程和轮机工程3个二级学科的博士学位和硕士学位授予权。
本学科教学科研设施先进,拥有的船舶与海洋工程实验室是辽宁省高校重点实验室,设有船舶CAD工程中心、船模拖曳试验水池、造船工艺实验室、船舶结构振动实验室、声学实验室和结构环境损伤控制实验室,实验室的实验设施已达到国内先进水平。
船舶与海洋工程学科师资力量雄厚,已在人才培养和科学研究等方面做出了突出成绩。
近年来承担了国家863计划、国家科技攻关、国家自然科学基金、国防军工等方面的科研项目,获国家科技进步奖二等奖1项,省部级科技进步奖一等奖1项,二等奖2项,辽宁省科技转化一等奖2项,大连市政府科技进步奖一等奖4项;获国家专利技术15项。
理论声学习题10
理论声学 Theoretical Acoustics习题10.1 有一声波作用在赫姆霍兹共鸣器管口,求共鸣器内部的声压。
如果声波的频率等于共鸣器的共振频率,求共鸣器内的声压与口外声压的比。
管口声压0P ,0a a aQPi M Q R Q i C ωω=-++-共鸣器内声压P ,aQP i C ω=-01/1/aa a ai C P P i M R i C ωωω-=-+-,共振频率,ω=0a P P =习题10.2 画出左图所示带通滤波器的类比电路图,求出共振频率。
习题10.3 排气管截面积为S ,侧面接一长l ,截面积S 的短管,它能降低通过管道向外辐射的声波吗?入管()()()11121exp exp A P x P ikx P ikx =+-()()()111210exp exp A SQ x P ikx P ikx c ρ=+-⎡⎤⎣⎦ 出管()()131exp B P x P ikx =()()1310exp B SQ x P ikx c ρ=短管()()()24252exp exp C P x P ikx P ikx =+-()()()242520exp exp C SQ x P ikx P ikx cρ=--⎡⎤⎣⎦ 2x l =处0C Q =,()46exp P P ikl =-,()56exp P P ikl =120x x ==处,()()()000A B C P P P ==,()1234562cos P P P P P P kl +==+=()()()0000A B C Q Q Q -+=,123450P P P P P --+-=()()()123453632sin 1tan P P P P P P iP kl P i kl -=--=+=+⎡⎤⎣⎦ ()1322tan P P i kl =+⎡⎤⎣⎦,3P =习题11.1 在0z =界面两侧的两种介质的声速分别是1c 和2c 。
《水声学》习题作业-射线声学理论
一、简答题1.说明射线声学的基本方程、适用条件及其局限性,并说明球面波和柱面波传播时声线的传播方向。
2.水平分层介质中的“程函方程”如何表示?若海水中的声速分布如下图,试画出几条典型声线轨迹示意图。
3.声线弯曲满足的基本条件是什么?并定性说明它们之间的规律。
4.水平分层介质中已知某条声线的起始掠射角为0,起点处介质声速为0c ,海水介质的声速梯度为大于零的常数g 。
试证该声线轨迹为半径g c R 0 的一段圆弧。
5.海水中声速从海面的1500m/s 均匀地减小到100m 深处的1450m/s 。
求:(1)声速梯度;(2)从海面水平出射的声线传播到100m 深处时,声线传播的水平距离是多少?(3)上述声线到达100m 深处时的掠射角为多少?6.设海水中有负声速梯度,且其绝对值为常数g ,声源处的声速为0c 。
试证:水平发出的声线穿过的水层厚度为d 时,它在水平方向前进的距离()2102g d c r =。
7.设海水中有恒定负声速梯度,其绝对值为常数g ,海面声速为0c ,声源深度为d 。
试求恰巧在海面反转的声线的出射角(与水平面之夹角)?8.某浅海海域水深40m ,海面、海底都是平面。
声源深度10m ,声速梯度为常数,海面声速为1500m/s ,海底处为1480m/s 。
试计算并画出自声源沿水平方向发出的声线的轨迹,到第二次从海底反射为止。
9.驱逐舰要搜索一艘水中的敌潜艇,海水中声速梯度为-0.1/s ,海面声速为1500m/s 。
驱逐舰上声纳换能器的深度为10m ,当换能器的俯角为4.5o 时,发现水平距离1000m 处的潜艇,问潜艇的深度为多少?10.一艘潜艇位于180m 深处,该处声速为1500m/s ,海水中声速梯度为0.1/s 。
它的声纳换能器的仰角为10o 时探测到一艘水面船只。
问水面船只离潜艇的水平距离是多少?11.聚集因子F 是如何定义的,它有什么物理意义?举出二个F>1的场合。
利用声压场的浅海地声参数贝叶斯反演方法
第42卷第4期2021年4月哈㊀尔㊀滨㊀工㊀程㊀大㊀学㊀学㊀报Journal of Harbin Engineering UniversityVol.42ɴ.4Apr.2021利用声压场的浅海地声参数贝叶斯反演方法郑广学3,朴胜春1,2,3,祝捍皓4,5,6,张海刚1,2,3,李楠松1,2,3(1.哈尔滨工程大学水声技术重点实验室,哈尔滨150001;2.海洋信息获取与安全工信部重点实验室(哈尔滨工程大学),工业和信息化部,哈尔滨150001;3.哈尔滨工程大学水声工程学院,哈尔滨150001;4.浙江海洋大学海洋科学与技术学院,浙江舟山316022;5.浙江大学浙江省海洋观测-成像试验区重点实验室,浙江舟山316021;6.中国科学院声学研究所声场声信息国家重点实验室,北京100190)摘㊀要:针对浅海考虑弹性海底下的地声参数反演和对反演结果的不确定性分析问题,本文提出了一种基于贝叶斯理论的利用声压场信息的反演方法㊂首先通过快速场方法(FFM )计算声压场的理论预报值,其次根据贝叶斯理论基于似然函数建立实测声压场与理论声压场间的目标函数,最后依照Metropolis 准则采样给出待反演参数的后验概率密度(PPD )作为反演结果㊂结果表明:依据本方法获得的反演结果与经典遗传算法(GA )寻优结果基本吻合;反演声压场与水池实测声场的相关系数达0.89;反演结果的不确定性表明与横波声速衰减㊁纵波声速衰减相比,弹性海底中的横波声速㊁纵波声速以及密度的不确定性更小㊂关键词:地声反演;快速场方法;Metropolis 准则;贝叶斯理论;声压场;不确定性分析;遗传算法;水池实验DOI :10.11990/jheu.201906062网络出版地址:http :// /kcms /detail /23.1390.u.20210324.1209.008.html 中图分类号:O427.1㊀文献标志码:A㊀文章编号:1006-7043(2021)04-0497-08Bayesian inversion method of geo-acoustic parameter in shallow seausing acoustic pressure fieldZHENG Guangxue 3,PIAO Shengchun1,2,3,ZHU Hanhao4,5,6,ZHANG Haigang1,2,3,LI Nansong 1,2,3(1.Acoustic Science and Technology Laboratory,Harbin Engineering University,Harbin 150001,China;2.Key Laboratory of Marine Information Acquisition and Security (Harbin Engineering University ),Ministry of Industry and Information Technology;Harbin 150001,China;3.College of Underwater Acoustic Engineering,Harbin Engineering University,Harbin 150001,China;4.Institute of Marine Science and Technology,Zhejiang Ocean University,Zhoushan 316022,China;5.Key Laboratory of Ocean Observation-Ima-ging Testbed of Zhejiang Province,Zhejiang University,Zhoushan 316021,China;6.State Key Laboratory of Acoustics,Chinese Academy of Science,Beijing 100190,China)Abstract :The present study develops an approach for geoacoustic parameters inversion based on the Bayesian Theo-ry to solve the uncertainty analysis of geoacoustic parameters inversion in shallow sea.The fast field method (FFM)is used to calculate the theoretical prediction value of the sound pressure field.The objective function is then estab-lished based on the Bayesian theory to match the measured sound pressure field and the theoretical pressure field.Finally,the posterior probability density (PPD)is given as the inversion result according to the Metropolis rule.Inversion results obtained by this method are in good agreement with those obtained by the classical genetic algo-rithm (GA).The correlation coefficient between the inversion sound pressure field and the measured sound field in the tank experiments reaches 0.89.Besides,S-wave velocities,P-wave velocities,and seafloor density have fewer uncertainties than S-wave attenuation and P-wave attenuation.Keywords :Geo-acoustic inversion;the fast field method;Metropolis rule;Bayesian theory;sound pressure field;uncertainty analysis;genetic algorithm;tank experiment收稿日期:2019-06-18.网络出版日期:2021-03-25.基金项目:国家自然科学基金项目(11704337);浙江省海洋观察-成像试验区重点实验室项目(OOIT2018OF05);声场声信息国家重点实验室开放课题研究基金项目(SKLA201901).作者简介:郑广学,男,博士研究生;朴胜春,男,教授,博士生导师;祝捍皓,男,副教授.通信作者:祝捍皓,E-mail:zhuhanhao@.㊀㊀海洋水体环境参数与地声参数共同构成了水下声传播计算最重要的环境参数,相比水体环境参数,地声参数更难直接测量,因而对其的有效获取研究一直是水声学中的热点课题[1]㊂由于地声参数的变化会对声压场的分布产生重要影响,因此可以考虑以水听器接收到的复声压作为研究对象来反演地声参数㊂相比于直接测量的方法,利用声学方法反演可以快速㊁高效的获取大面积海区的底质参数,避哈㊀尔㊀滨㊀工㊀程㊀大㊀学㊀学㊀报第42卷免了人力㊁物力的浪费,因此受到了广泛的关注[2-5]㊂近年来,利用声学方法反演地声参数已取得了长足的发展㊂各种反演方法不断涌现,如:利用传播损失的地声参数反演方法[6];利用舰船辐射噪声等各类海洋噪声信息的地声参数反演方法[7-10];利用波导频散特性的反演方法等㊂但上述地声参数反演方法均着重研究反演问题中正演模型的选择,多仍采用传统寻优算法,如下降单纯形法(downhill simplex,DS)㊁遗传算法(genetic algorithm,GA)等对目标函数进行求解㊂但此类算法一般只能给出待反演参数的点估计,无法对参数的不确定性进行进一步描述,且在寻优过程中对算法初始模型依赖很强,容易陷入局部最优解㊂贝叶斯反演方法有效地估计最大后验概率(maximum a posteriori,MAP )模型参数,进而从统计角度定性和定量分析参数反演结果的不确定性㊂基于上述原因,本文以单水听器接收到各距离点上的声压场为研究对象,基于贝叶斯理论,建立了针对浅海考虑弹性海底下的地声参数反演方法,反演对象包括声速㊁声速衰减和密度3类海底声学参数㊂1㊀反演方法根据声场信息估计未知海底参数的声学反演方法一般由4部分组成:1)构建声传播模型;2)建立合适的目标函数;3)全局优化算法;4)反演结果的不确定性分析㊂在本文的反演方法研究中,同样针对上述4部分展开工作㊂1.1㊀浅海声场正演模型浅海波导中,海底横波声速对低频声信号的传播影响不可忽略[11],本文研究中选择了如图1所示的浅海参数化模型㊂其中,z =0表征海面,z 轴表深度,向下为正,r 正轴表示声场向外传播方向㊂设水深为H ㊁声源位于水深z s 处㊁密度ρ1㊁声速c 1;海底为半无限弹性介质,各向同性㊂纵波声速㊁横波声速㊁密度㊁纵波声速衰减和横波声速衰减分别用c p ㊁c s ㊁ρb ㊁αp 和αs 表示,该5项海底参数也是本文反演时的研究对象㊂图1㊀单层弹性海底的浅海参数化模型Fig.1㊀Parametric model in shallow sea with elastic bottom图1模型下,流体层中的各位置声压p (r ,z )为:p (r ,z )=ρ1ω2ʏɕ0Z (z ,ξ)J 0(ξr )ξd ξ(1)Z (z ,r )=A sin (βz ),0ɤz <z s B sin (βz )+C cos (βz ),z s <z <H 1{(2)式中:ξ为水平波数;r 为水平距离;J 0为贝塞尔函数;文献[12]给出了A ㊁B ㊁C ㊁β等参数的具体推导过程及表示式㊂对式(1)中声压场p (r ,z )的数值计算,通常采用简正波方法(normal mode method,NMM)及快速场方法(fast field method,FFM)[13]㊂考虑此模型为浅海环境,FFM 更适合用于声场的快速计算,故本文选择FFM 完成对上述参数化模型下声压场p (r ,z )的数值计算㊂在完成对声压场p (r ,z )数值计算的基础上,需建立恰当的目标函数来衡量实测声压p (r ,z )与理论预报声压pᶄ(r ,z )间的适配程度,并通过对其寻优求解便可实现该模型下各海底参数的获取㊂1.2㊀目标函数与优化算法贝叶斯方法用于反演问题时,观测数据向量d 和模型参数向量m 被视为随机的㊂根据贝叶斯理论有:P (m d )=P (d m )P (m )/P (d )(3)其中,P (m d )为后验概率密度(posterior probability density,PPD)㊂条件概率密度P (d m )定义为似然函数L (m ),为某一(固定)测量数据d 下m 的函数,代表数据提供的信息,用于衡量在模型参数m 下模型预报结果和数据的适配程度㊂P (m )为先验概率密度,表示模型参数先验信息,P (d )为测量数据的概率密度函数[14-15]㊂因此,后验概率密度既包含数据信息,也包含模型参数的先验信息㊂且P (d )与m 无关,可以看作1个常数,所以式(3)可化为:P (m d )ɖP (d m )P (m )(4)㊀㊀根据贝叶斯理论,目标函数是在高斯数据误差的假设下基于似然函数建立的,似然函数作为某一指定测量数据下模型参数的函数,代表了数据不确定性分布[14-15]㊂本文使用的目标函数为:E (m )=ln[B (m )|p |2](5)式中B (m )表示归一化的Bartlett 失配器:B (m )=1-|p ᶄ+(m )p |2|p ᶄ(m )|2|p |2(6)式中:p 表示单个水听器接收到的声压场(实测场),p ᶄ(m )为待反演参数向量m 的预报声压(拷贝场),上标 + 代表共轭转置㊂从式(6)可以看出,B (m )取得最优匹配时的条件为p (r ,z )=pᶄ(r ,z ),同时,E (m )取得极小值㊂目标函数的建立,将对m 的反演问题转化为在寻优范㊃894㊃第4期郑广学,等:利用声压场的浅海地声参数贝叶斯反演方法围内对式(5)的极小值求解问题㊂根据贝叶斯理论,上述待反演参数的PPD将依据Metropolis准则采样[16]给出,并从统计角度对其进行不确定性分析㊂2㊀声场正演模型仿真分析本文在完成了声场建模㊁确定了目标函数及寻优算法的基础上,利用数值模拟,对比图1模型下各海底参数与其反演结果,对研究方法进行验证㊂2.1㊀声传播特性分析图2给出了利用FFM计算所得图1所示浅海环境下,1~1500m100Hz声信号声压场的仿真结果,以传播损失(transmission loss,TL)曲线的形式给出㊂图2㊀5项参数对TL的影响程度分析Fig.2㊀The impact of five parameters on transmission loss㊀㊀仿真中,设声源深度为z s=20m,接收点z r位于水下10m处,表1给出了仿真中各海洋环境参数㊂图2(a)~(e)中的实线对应利用表1中各参数真值的计算结果,5项地声参数取讨论值时的计算结果分别用点线(较小值)和虚线(较大值)表示㊂图2 (f)中应用 距平 的定义[17],视设定的环境参数真值计算所得传播损失为均值,视各参数的讨论值计算所得传播损失为离散值来计算距平,用5项地声参数讨论值的距平大小来衡量各参数变化对传播损失的影响程度㊂表1中各参数的讨论值均设定为真值偏移ʃ10%,在此基础上讨论各参数变化对声压场的影响程度㊂从图2(a)~(e)的中可以看出,5项参数的变化对TL曲线影响程度各不相同,从图2(f)中5项参数讨论值的距平大小可以看出,在各参数均偏移ʃ10%的情况下,改变某一参数对TL曲线的影响程度从大到小依次为c p㊁c s㊁ρb㊁αp㊁αs㊂因此,5项海底参数对p(r,z)的影响程度可初步定性为:c p㊁c s>ρb>αp㊁αs㊂2.2㊀目标函数对各参数的敏感度分析在完成了正演建模㊁声压场数值计算的基础上,本文通过数值模拟仿真分析本文所给目标函数对5项地声参数的敏感度,并讨论该目标函数在反演中的可靠性㊂仿真中,海洋环境参数真值与各参数寻优区㊃994㊃哈㊀尔㊀滨㊀工㊀程㊀大㊀学㊀学㊀报第42卷间仍取表1中数值,取各参数设定真值下模拟计算得到p (r ,z )为参考真值㊂本文采用固定变量法来分析目标函数对单个地声参数的敏感度时,即其他参数固定,只在讨论区间内改变某一参数计算声压场理论预报值p ᶄ(r ,z ),由式(5)计算p (r ,z )与p ᶄ(r ,z )间的目标函数值E (m )作为敏感度评价标准㊂表1㊀海洋环境仿真参数Table 1㊀The simulation parameters of ocean environment 实验参数海深H /m 海水声速c 1/(m ㊃s -1)海水密度ρ1/(g ㊃cm -3)纵波声速c p /(m ㊃s -1)横波声速c s /(m ㊃s -1)海底密度ρb /(g ㊃cm -3)纵波衰减αp /(dB ㊃λ-1)横波衰减αs /(dB ㊃λ-1)设定真值10015001.025200010001.5000.5000.500讨论值1800㊁2200900㊁1100 1.350㊁1.6500.450㊁0.5500.450㊁0.550㊀㊀图3(a)~(e)分别给出了目标函数E (m )随表1中各单项参数变化时的数值变化规律㊂从图3可以看出,在5项地声参数讨论范围内,目标函数E (m )均只在其真值处取得极小值,且无伪峰或第2极小值,避免了局部最优解的陷阱,证明了该目标函数的可行性㊂但从图3(f)可知目标函数随各参数变化时的波动程度不一致,在5项参数的讨论范围内,目标函数值波动范围从大到小依次为c p ㊁c s ㊁ρb ㊁αp ㊁αs ㊂因此,5项海底参数对p (r ,z )的影响程度可定性为:c p >c s >ρb >αp >αs ,进一步证明了图2的讨论结果㊂图3㊀单参数敏感度Fig.3㊀The sensitivity of single parameter㊃005㊃第4期郑广学,等:利用声压场的浅海地声参数贝叶斯反演方法2.3㊀仿真反演结果分析在完成目标函数对5项海底声学参数敏感度的研究基础上,本节通过数值仿真,分析本文所研究的反演方法对目标函数寻优结果的可靠性㊂海洋环境参数如表1所设,当寻优结束时,各参数寻优结果与仿真真值的对比如图4给出㊂图中水平线段表示了各参数PPD 的均值位置及其方差大小,垂直于横坐标直线段标注仿真真值所在位置,各参数反演结果在表2中给出㊂图4㊀5项参数及其概率密度分布Fig.4㊀Five parameters and their probability density dis-tribution表2㊀寻优结果Table 2㊀The results of optimization实验参数设定真值寻优结果(均值ʃ方差)纵波声速c p /(m ㊃s -1)20001997.232ʃ48.878横波声速c s /(m ㊃s -1)10001002.761ʃ27.868海底密度ρb /(g ㊃cm -3)1.500 1.499ʃ0.045纵波衰减αp /(dB ㊃λ-1)0.5000.499ʃ0.014横波衰减αs /(dB ㊃λ-1)0.5000.498ʃ0.014㊀㊀从图4和表2中的结果可以看出,各参数的仿真真值与反演结果吻合程度较高,均处于其寻优结果的均值与方差的和或差范围内㊂且各参数的概率密度峰值的尖锐程度也反映了目标函数对各参数的敏感程度,如图4所示,目标函数对各参数的敏感程度为:c p ㊁c s >ρb >αp ㊁αs ,这也与2.1节㊁2.2节的讨论结果相吻合㊂为验证数值仿真中反演结果的可靠性,图5给出了仿真环境参数下TL 曲线与仿真反演结果计算所得TL 曲线的对比图㊂等间距取12处绘制误差棒,并标注反演结果计算所得TL 曲线的误差范围㊂从图中可以看出,仿真中TL 曲线均处于反推TL 曲线误差范围之内,证明了该反演方法在数值仿真应用中的可靠性㊂图5㊀TL 曲线对比Fig.5㊀The comparison of TL㊃105㊃哈㊀尔㊀滨㊀工㊀程㊀大㊀学㊀学㊀报第42卷3㊀消声水池缩比实验验证3.1㊀水池实验本文结合消声水池缩比实验[18],利用实验数据对本文所研究反演方法进行讨论㊂实验中选取质地均匀㊁高硬度的塑料板模拟 半无限弹性海底 ㊂表3给出了实验中各设备的布放参数,其中声源深度z s 为87mm,接收深度z r 为84mm,水深z 为182mm,水中声速c 1为1450.212m /s 由测量水中温度(11.15ħ)后通过经验公式求得㊂测量过程中,声源固定不动,发射信号为135kHz 的脉冲信号;接收水听器固定在可移动走架上,采集卡采样率fs 为20MHz,单个测量点分别测量10次以减少随机扰动带来的误差;完成一点测量后,走架带动水听器向远离声源方向移动2mm(误差不超过20μm)至下1个测量点,实验中测量点共720个㊂测量中TL 曲线如图6(b)所示,图6(a)给出了测量中所接收到的第1~100道接收信号时域图㊂图6㊀实验中接收到的信号Fig.6㊀The signals received in experiment3.2㊀反演结果验证与分析表3列出了各地声参数的寻优范围,以及利用本文所研究反演方法对上述实验数据反演所得到的 海底 地声参数值,同时还给出了利用GA 对该实验数据寻优的结果作为对比㊂表3㊀寻优结果Table 3㊀The results of optimization实验参数寻优区间GA 寻优结果寻优结果(均值ʃ方差)纵波声速c p /(m ㊃s -1)2200~25002379.0582379.058ʃ69.630横波声速c s /(m ㊃s -1)1100~13001191.1271191.127ʃ41.883海底密度ρb /(g ㊃cm -3)1.0~1.8 1.210 1.210ʃ0.178纵波衰减αp /(dB ㊃λ-1)0.1~1.10.5850.585ʃ0.212横波衰减αs /(dB ㊃λ-1)0.1~1.10.7730.773ʃ0.270㊀㊀图7给出了算法终止时,各参数在寻优区间上的分布及其概率密度分布,与图4相仿,图7中垂直于横坐标直线段标注GA 寻优结果所在位置,水平线段表示了各参数取值PPD 的均值位置及其方差大小㊂图8定量给出了参数间的相关系数矩阵图㊂从表4㊁图7和图8的结果中可以看到:首先,各待反演参数的GA 寻优结果均处于均值与方差的和或差范围内,证明了利用贝叶斯反演方法获得的结果与传统反演方法得到结果相吻合㊂另外,方差的相对大小从一方面也反映出该参数反演结果不确定性的大小㊂其次敏感度较小的参数对目标函数贡献也较小,在图8中表现为该参数与其他参数的相关性偏小,如αp 和αs ;换言之,较敏感参数间则表现出较明显的相关性,如c p ㊁c s ㊁ρb ,综上,可以看出5项待反演参数的敏感度大小分别为c p ㊁c s >ρb >αp ㊁αs ,这也与第2节结论一致㊂㊃205㊃第4期郑广学,等:利用声压场的浅海地声参数贝叶斯反演方法图7㊀算法终止时5项参数及其概率密度分布Fig.7㊀The five parameters and their probability densitydistribution at the termination of thealgorithm图8㊀参数间相关矩阵Fig.8㊀The correlation matrix of parameters㊀㊀为进一步验证表4反演结果的可信度,图9对比了利用反演结果计算得到的TL 曲线和实测TL 曲线;图10对比了同一测量点(第470号接收点),反演所得时域波形与该点实际接收信号的时域波形图(均进行归一化处理)㊂图9中,实线为实测TL 曲线,点线为应用表4反演结果计算所得的TL 曲线并等间距取12处绘制误差棒,从图中对比结果可以看到实验中TL 曲线基本处于反推TL 曲线误差范围之内,证明了该反演方法在实际应用中的可行性;图10中,实线为第470号接收点实测时域波形,点线为应用表4反演结果 预报 的第470号接收点时域波形,从图中对比结果可以看到时域波形计算预报结果与实测信号波形匹配较好,若不考虑信号 拖尾 ,相关系数可达0.89,证明了本反演方法在实际应用中的可靠性㊂图9㊀实测TL 与反演所得TL 对比Fig.9㊀The comparison of TL between measured and in-versed图10㊀时域波形对比Fig.10㊀The comparison of time domain waveform4㊀结论1)贝叶斯反演能够从统计角度给出反映地声参数向量综合信息的后验概率来表征反演结果,同时还能够分析各参数反演结果的不确定性,相比传统只能给出待反演参数点估计的反演方法,反演信息更加全面㊁科学㊂2)根据贝叶斯反演理论,利用本文给出的目标函数,运用Metropolis 准则采样,可以较为准确得到待反演参数的PPD㊂数值模拟和水池实验的结果均表明:PPD 分布特性与GA 寻优结果吻合程度较好,利用反演参数计算得到的声场特性与测量场基㊃305㊃哈㊀尔㊀滨㊀工㊀程㊀大㊀学㊀学㊀报第42卷本一致,时域波形相关系数达0.89㊂3)从5项地声参数反演结果的不确定性大小来看,在仿真模拟及水池实验的环境下,纵波声速c p㊁横波声速c s以及密度ρb的不确定性相对较小,对声压场的变化更敏感,5项参数依据敏感程度大小依次为c p>c s>ρb>αp>αs,声场研究结果与反演结果均证实了这一研究结果㊂参考文献:[1]臧涛,海伟.浅海海底声学参数反演方法研究[J].舰船电子工程,2015,35(11):154-156,172. ZANG Tao,HAI Wei.Inversion method of sea-bottom acoustics parameters in shallowwater[J].Ship electronic engineering,2015,35(11):154-156,172. [2]GERSTOFT P.Inversion of seismoacoustic data using ge-netic algorithms and a posteriori probability distributions [J].The journal of the acoustical society of America, 1994,95(2):770-782.[3]DALL'OSTO D R,DAHL P H.Geoacoustic inversion based on particle velocity[J].The journal of the acoustical society of America,2016,139(4):2125.[4]李倩倩,阳凡林,张凯,等.不确定海洋环境中基于贝叶斯理论的声源运动参数估计方法[J].物理学报, 2016,65(16):164304.LI Qianqian,YANG Fanlin,ZHANG Kai,et al.Moving source parameter estimation in an uncertain environment [J].Acta physica sinica,2016,65(16):164304. [5]李倩倩,李整林,张仁和.不确知海洋环境下的贝叶斯声源定位法[J].声学学报,2014,39(5):535-543. LI Qianqian,LI Zhenglin,ZHANG Renhe.Bayesian local-ization in an uncertain ocean environment[J].Acta acusti-ca,2014,39(5):535-543.[6]屈科,胡长青,赵梅.利用传播损失反演海底单参数[J].声学学报,2013,38(4):472-476.QU Ke,HU Changqing,ZHAO Mei.Single parameter in-version using transmission loss in shallow water[J].Acta acustica,2013,38(4):472-476.[7]CROCKER S E,NIELSEN P L,MILLER J H,et al.Geo-acoustic inversion of ship radiated noise in shallow water using data from a single hydrophone[J].The journal of the acoustical society of America,2014,136(5):EL362-EL368.[8]HARRISON C H,SIMONS D G.Geoacoustic inversion of ambient noise:a simple method[J].The journal of the a-coustical society of America,2002,112(4):1377-1389.[9]周建波,朴胜春,黄益旺,等.波浪起伏对噪声场空间特性的影响[J].哈尔滨工程大学学报,2017,38(7): 1056-1064.ZHOU Jianbo,PIAO Shengchun,HUANG Yiwang,et al. Effect of wave fluctuation on the spatial characteristics of noise field[J].Journal of Harbin Engineering University, 2017,38(7):1056-1064.[10]江鹏飞,林建恒,马力,等.一种海洋环境噪声分步反演地声参数方法[J].声学学报,2016,41(1):59-66.JIANG Pengfei,LIN Jianheng,MA Li,et al.A multi-step method of geo-acoustic inversion by ambient noise [J].Acta acustica,2016,41(1):59-66. [11]祝捍皓,郑广学,张海刚,等.浅海环境下低频声信号传播特性研究[J].上海交通大学学报,2017,51(12):1464-1472.ZHU Hanhao,ZHENG Guangxue,ZHANG Haigang,et al.Study on propagation characteristics of low frequency acoustic signal in shallow water environment[J].Journal of Shanghai Jiao Tong University,2017,51(12):1464-1472.[12]祝捍皓.基于声矢量场阻抗特征的地声参数反演研究[D].哈尔滨:哈尔滨工程大学,2014.ZHU Hanhao.Geoacoustic parameters inversion based on waveguide impedance in acoustic vector field[D].Har-bin:Harbin Engineering University,2014. [13]DINAPOLI F R.Fast field program for multi-layeredmedia[R].Middletown,Rhode Island:Naval Underwater Systems Center Newport Ri,1971:2-10. [14]BELCOURT J,DOSSO S E,HOLLAND C W,et al.Bayesian geoacoustic inversion of seabed reflection data at the New England mud patch[J].The journal of the acous-tical society of America,2017,142(4):2590.[15]于盛齐.基于反向散射强度的海底参数反演方法研究[D].哈尔滨:哈尔滨工程大学,2014.YU Shengqi.Geoacoustic parameter inversion based on backscattering strength[D].Harbin:Harbin Engineering University,2014.[16]高飞,潘长明,孙磊.Bayes匹配场地声参数反演:多步退火Gibbs采样算法[J].兵工学报,2017,38(7): 1385-1394.GAO Fei,PAN Changming,SUN Lei.Geoacoustic pa-rameters inversion of Bayes matched-field:a multi-annea-ling Gibbs sampling algorithm[J].Acta armamentarii, 2017,38(7):1385-1394.[17]黄嘉佑.气象统计分析与预报方法[M].北京:气象出版社,1990:5-6.HUANG Jiayou.Methods of meteorological statistical anal-ysis and forecast[M].Beijing:China Meteorological Press,1990:5-6.[18]祝捍皓,朱军,郑广学,等.一种浅海声简正波近场分离方法[J].声学学报,2019,44(1):39-48.ZHU Hanhao,ZHU Jun,ZHENG Guangxue,et al.A sep-eration method for normal modes in shallow water under near field[J].Acta acustica,2019,44(1):39-48.本文引用格式:郑广学,朴胜春,祝捍皓,等.利用声压场的浅海地声参数贝叶斯反演方法[J].哈尔滨工程大学学报,2021,42(4):497-504. ZHENG Guangxue,PIAO Shengchun,ZHU Hanhao,et al.Bayesian inversion method of geo-acoustic parameter in shallow sea using acoustic pressure field[J].Journal of Harbin Engineering University,2021,42(4):497-504.㊃405㊃。
海底声学探测技术应用考核试卷
A.海水温度降低
B.海水盐度增加
C.海水压力降低
D.声波频率增加
16.多波束测深系统与单波束测深仪相比,主要优势是?()
A.测量精度更高
B.测量范围更广
C.抗干扰能力更强
D.设备体积更小
17.以下哪个现象可能导致声波在海水中的传播路径发生弯曲?()
A.海水温度的垂直梯度
B.海水盐度的水平梯度
C.海水压力的垂直梯度
D.海水流速的水平梯度
18.在海底声学探测中,以下哪个参数可以用来描述声波传播过程中的能量损失?(")
A.声速
B.吸收系数
C.反射系数
D.透射系数
19.以下哪个设备主要用于探测海洋中的微小颗粒物质?()
A.单波束测深仪
B.侧扫声纳
C.多波束测深系统
D.回声探测仪
10.影响声波在海底沉积物中反射和吸收的主要因素是?()
A.沉积物的密度
B.沉积物的粒度
C.沉积物的湿度
D.沉积物的温度
11.在海底声学探测中,以下哪个因素可能导致声波传播距离变短?()
A.海水温度升高
B.海水盐度增加
C.海水压力增加
D.声波频率降低
12.以下哪个设备通常用于海洋油气资源勘探?()
C.沉积物粒度增大
D.声波频率降低
14.侧扫声纳在海底探测中的应用主要包括以下哪些方面?()
A.探测海底地形
B.探测海底物体
C.评估海底沉积物类型
D.测量海水温度
15.以下哪些因素会影响声波在海底沉积物中的反射系数?()
A.沉积物的粒度
B.沉积物的湿度
张海澜《理论声学》答案Theoretical Acoustics
1 Mg 1 Mg K x K xa cos 0t 2 K 2 K
2
2 Mg 1 K xa cos 0t Mgxa cos 0t 2 2K
2
重力势能: E3 Mgx Mgxa cos 0t
1 .2
一弹簧竖直悬挂,上端固定,下端系一质点组成简单质点振动系统。质点同时受到 向下的重力。分析质点的振动和能量的转换。
M
d 2x d 2x K Mg , , x 0 弹簧受力为零 Mg Kx x g ,特解 x0 2 2 dt K dt M
Mg Mg ,平衡位置 x ,简谐振动 K K
2 2 0 02 / 2 0 2 2 2 2 0 / 2 0 0 / 2 0 0 02 / 2 2 02
1
1 1 1 1 1,2 0 ,模式是 , ,模式是 0 1 2 / 2 2 0 2 2 1 1
张海澜的《理论声学》答案 1 .1 有一质点振动系统,固有频率 f 0 已知,质量和弹性系数待求。现在质点上增加已知 质量 m ,固有频率改变了 f 0' 。求原质点系统的质量和弹性系数。
2 f 0
K M
增加质量 m , 2 f 0 f 0
'
K M m
2 2
m f 0 f 0' 4 2mf 0 2 f 0 f 0' f0 m 可得 ,M ' ,K 1 f 0 f 0' M f 0 2 f 0 f 0' f 0' 2 f 0 f 0'
(整理)声学基础课后答案
习题11-1 有一动圈传声器的振膜可当作质点振动系统来对待,其固有频率为f ,质量为m ,求它的弹性系数。
解:由公式mmo M K f π21=得: m f K m 2)2(π=1-2 设有一质量m M 用长为l 的细绳铅直悬挂着,绳子一端固定构成一单摆,如图所示,假设绳子的质量和弹性均可忽略。
试问:(1) 当这一质点被拉离平衡位置ξ时,它所受到的恢复平衡的力由何产生?并应怎样表示? (2) 当外力去掉后,质点m M 在此力作用下在平衡位置附近产生振动,它的振动频率应如何表示? (答:lgf π210=,g 为重力加速度)图 习题1-2解:(1)如右图所示,对m M 作受力分析:它受重力m M g ,方向竖直向下;受沿绳方向的拉力T ,这两力的合力F 就是小球摆动时的恢复力,方向沿小球摆动轨迹的切线方向。
设绳子摆动后与竖直方向夹角为θ,则sin lξθ=受力分析可得:sin m m F M g M glξθ==(2)外力去掉后(上述拉力去掉后),小球在F 作用下在平衡位置附近产生摆动,加速度的方向与位移的方向相反。
由牛顿定律可知:22d d m F M tξ=-则 22d d m m M M g t l ξξ-= 即 22d 0,d gt lξξ+=∴ 20g l ω=即 01,2πgf l= 这就是小球产生的振动频率。
1-3 有一长为l 的细绳,以张力T 固定在两端,设在位置0x 处,挂着一质量m M ,如图所示,试问: (1) 当质量被垂直拉离平衡位置ξ时,它所受到的恢复平衡的力由何产生?并应怎样表示?(2) 当外力去掉后,质量m M 在此恢复力作用下产生振动,它的振动频率应如何表示?(3) 当质量置于哪一位置时,振动频率最低? 解:首先对m M 进行受力分析,见右图,0)(2202200=+-+--=εεx x Tx l x l TF x(0x 〈〈ε ,2022020220)()(,x l x l x x -≈+-≈+∴εε 。
水声学原理:第2章 海洋中声场的基本理论
哈尔滨工程大学
水声学原理
14
硕士学位课程
2.1 海洋中声场的射线理论
• 三维折射
– 当介质的折射率是三维坐标的函数时,声线将不能保 持在同一个平面内。
– 在海洋声学中,当分析内波对声场的影响时将会遇到 三维折射的情形。
– 在研究距离有关(Range-dependent)海洋中声波远 距离传播时也会遇到相同情况。
哈尔滨工程大学
水声学原理
26
硕士学位课程
2.1 海洋中声场的射线理论
• 海洋声层析(Ocean Acoustic Tomography)
– 海洋的变化特性对海洋气候、地球的天气都有显著的 影响,显著改变了声场的层状结构、导致了声信号的 起伏、扰动了声线路径。
– 从调查船和卫星获取的水团特征信息也是非常丰富的, 但尽管如此,这对实际生产来说还不够充分,因为目 前所获取的信息仍然只是海面和海面表层的信息。
(2A • W
A 2W
)
k
2 0
A[n 2
(W )2 ]
0
– 当声波频率足够高时,由上式可得程函方程和输运方 程(transport equation)
(W )2 n2
2A • W A2W 0
哈尔滨工程大学
水声学原理
16
硕士学位课程
2.1 海洋中声场的射线理论
• 三维折射
– 程函方程:
• 定义了声线几何坐标,声线垂直于等相位面
– 从某个距离开始,c 将超过海底的声速,在这种情况下,
声线将在海底发生反射。在远距离声传播中,由于这
种声线在海底的多次反射导致了较大的总声功率损失
哈尔滨工程或大者学 高吸收海底情况1 海洋中声场的射线理论
张海澜 理论声学 答案7-12
2Vcos f f 0 1 c
频移 300Hz 习题 7.5 给出由(7.25)和(7.26)到(7.30)的详细推导。
(7.25)
v
0 r
1
p 0
1 v p 0 0 r
0 p 0 r c2 r v
有效声压方向特性的指向性函数是sinsinkaka10148340kasin9014814807505sin0kakakaka20log07525db习题72点声源向空间辐射200hz的声波将一压差式传声器放在离声源001m或1mcosikrpsik200100037340kr37340kr103720log2600683db习题73有两个相同的压强式传声器相距d把它们的输出串联组成复合传声器分析复合传声器对不同方向入射的声波的响应
习题 8.4 矩形声波导的截面尺寸是 0.1m0.08m, 一端有声源发出 1000Hz 到 2000Hz 的声波, 另一端伸向无限远,分析其中的声波。
P ( x, y , z )
m ,n 0
AmnU mn ( x, y ) exp(ikmn z )
2 2 2
m x n y m n , kmn U m ( x, y ) mn cos cos a b c a b
11 2 ik r r r22 ikr1 1 r2 F1 1 1 2 F2 r1 ikr2 1 r 2 11 1 ik r r 1 1
k r k r
2 2 1
1
2 2 2
0.037 1 2600 1 10 3.7 1
波函数
铝熔体中超声空化区域分布及声压衰减系数
铝熔体中超声空化区域分布及声压衰减系数章武;李晓谦;蒋日鹏【摘要】利用激光测距仪测得铸造用超声辐射杆端面全局振幅及频率,基于 FFT 细化分析对所测振幅分析后得到辐射杆端面振幅整体分布特性,在铝熔体中,将钛板竖直置于超声辐射杆正下方,启动超声振动系统运行20 h,通过分析钛板空蚀区域探究超声在铝熔体中的空化区域,得出辐射杆的空化区域分布,空化区域为一椭球形且分界面明显,空化域最远可达变幅杆正下方68 mm,结合铝熔体中超声的空化阈值,推算出20 kHz 铸造用功率超声在铝熔体中传播的衰减系数。
%The global ultrasonic amplitude and frequency of an ultrasonic radiation rod's end face was measured using a laser rangefinder,and the data obtained were analyzed by utilizing FFT zoom analysis technique to obtain the global amplitude distribution of the rod's end face.A titanium plate was placed under the ultrasonic radiation rod in molten aluminum.An ultrasonic vibration system was activated for 20 h to find the ultrasonic cavitation area in molten aluminum by analyzing the cavitation corrosion area of the titanium plate,then the caviation area distribution of the ultrasonic radiation rod was obtained.It was shown that the cavitation area of the rod is an ellipsoid with obvious boundaries,and the longest distance between the boundaries and the ultrasonic radiation rod is 68mm.Moreover,the sound pressure attenuation coefficient of 20 kHz ultrasound was calculated with the ultrasonic cavitation threshold in molten aluminum.【期刊名称】《振动与冲击》【年(卷),期】2016(035)021【总页数】4页(P150-152,170)【关键词】超声振动系统;空化腐蚀;空化域;声压衰减系数【作者】章武;李晓谦;蒋日鹏【作者单位】中南大学轻合金研究院,高性能复杂制造国家重点实验室,长沙410083;中南大学轻合金研究院,高性能复杂制造国家重点实验室,长沙 410083;中南大学轻合金研究院,高性能复杂制造国家重点实验室,长沙 410083【正文语种】中文【中图分类】TG249.9;TB559超声铸造是一种无污染、高效率且极具发展潜力的新型铸造技术,目前的研究认为,超声波在熔体中产生的空化效应能显著细化晶粒,提高材料性能,在强超声场正负声压的高频交替作用下,形成空化泡,空化泡崩溃瞬间产生的高温高压以及高强度的微射流,促使铝合金熔体中的初生晶被打碎,异质结晶核数目增多,结晶核与固相间的润湿角被减小,实际结晶温度降低,过冷度增大,从而使凝固组织得到明显的细化[1-5]。
医学超声成像波束形成技术的演变
继续教育Development of beamforming for medical ultrasonic imaging医学超声成像波束形成技术的演变张海澜*,林伟军,王秀明(中国科学院声学研究所声场声信息国家重点实验室,北京 100190)[Key words] U ltrasonics;Beam for ming [关键词] 超声学;波束形成 [中图分类号] R445 1 [文献标识码] A [文章编号] 1003 3289(2010)07 1395 04[基金项目]国家自然科学基金(10874202)、国家高技术研究发展计划(2006A A060103)。
[作者简介]张海澜(1946 ),男,江苏苏州人,博士,研究员。
研究方向:声学,超声学,声成像。
[通讯作者]张海澜,中国科学院声学研究所声场声信息国家重点实验室,100190。
E mail:zhanghl@m [收稿日期]2010 04 06 [修回日期]2010 05 11在过去几十年中,随着电子和计算机技术的发展,医学诊断超声成像技术发生了巨大变化,其临床应用价值大为提高。
在组成超声成像系统的各个部分中,发射和接收超声波、采集检查部位数据的前端部分是制约整个系统性能的关键。
本文从原理上介绍这方面的一些发展。
1 波束形成目前的超声成像系统几乎均采用阵列探头及电子扫描和聚焦方式发射超声波束,其原理示意图见图1。
图1A 表示仪器产生的N 路激发电信号,它们的脉冲波形相同,但有不同的延迟;图1B 中的T 表示N 个单元组成的阵列探头,每个单元都可以激发和接收超声波。
图1A 中第1路到第N 路激发信号经过放大后分别激发阵列探头的第1到第N 个单元,向探头前方发出超声波。
由于不同单元的激发电信号有不同的延迟,各个单元发出的超声脉冲不同步。
在如图1所示的电信号激励下,两侧的单元发射时间比较早,中间的比较迟,各个单元发射的超声波向中间汇聚,先后到达图1B 中曲线1、2、3等所示的位置,并在F 点聚焦,F 点称为焦点。
理论声学习题5
习题 5.1 复习声波方程的推导。
基本假设基本方程,线性化0ρρ=-∇⋅v 0p ρ=-∇v20p c ρ=,02000,1s s S K dp c d ρρρβρρ'='⎛⎫=== ⎪'⎝⎭ 不可压缩流体s K →∞,0c →∞,ρ=0,∇⋅v =0,速度无散,声压调和20p ∇=一维的解(),p x t ,220p x∂=∂,()()00p p t f t x ρ=- ,x v f = ,()x v f t =。
声波()()()()()exp exp exp 1exp i kx t ikx i t ikx i t ωωω-=-→+-⎡⎤⎣⎦()()0011exp exp x v i kx t i t c cωωρρ=-→-⎡⎤⎣⎦ 习题5.2如果介质中有体力分布,作用在单位体积上的体力为()z y x F ,, ,求声波方程。
分析体力是重力的情况。
体力 (),,Vx y z dxdydz =⎰⎰⎰f F运动方程 dxdydz F dx x p dydz t v dxdydzx x +∂∂-=∂∂0ρ F v +-∇=∂∂p t 0ρ 波动方程 F ⋅∇=∂∂-∇22022tc p p 重力 zF ˆg ρ'= zc p g z g z g ∂∂=∂∂=∂'∂=⋅∇20ρρF z c p g t c p p ∂∂=∂∂-∇2022022估计重力的影响t ω∂∂ ,k z∂∂ ,因此2022220pgc z gk g p c c t ωω∂∂=∂∂ ,这个数一般很小,频率很低时才有作用 对于空气如果1g c ω=,相当于0.0052g f c π=≈⨯Hz习题5.3 设夏天的气温比冬天高40︒C ,如果大气密度近似认为不变,求这样的温差下同样声压的声波声强变化的百分比和声强级差。
理想气体绝热过程声速2000p c γρ=, 理想气体平衡状态方程0p V NkT =,N 分子数,T 绝对温度,231.3810k -=⨯J/K波尔兹曼常数0p Nk T mρ=,0c ∝声强2002e p I c ρ=7==% 合 2110log(/)0.3I I =dB习题5.4 空气中声学脉冲测距仪向目标发射一声脉冲,测量反射回波的到达时间,确定目标距离。
SEA理论
Physical Problem – Sources 系统的输入能量
气动声源 结构振动
hard drive fan
Copyright © ESI Group, 2007. All rights reserved.
Physical Problem – Path
传递路径
共振、连接、辐射……
hard drive fan structure-borne sound, air-borne sound
Copyright © ESI Group, 2007. All rights reserved.
Physical Problem – SourcePath-Receiver Model
能量如何在结构声学系统内部传递?
Copyright © ESI Group, 2007. All rights reserved.
Building
Marine
Copyright © ESI Group, 2007. All rights reserved.
Who uses VA One SEA?
International distribution More than 200 customers
Automotive DC Ford GM Honda Mitsubishi Rover-BMW Renault SAAB Isuzu Nissan Caterpillar Deere Aerospace Boeing Airbus NASA Aerospace Corp. Lockheed Hughes TRW Raytheon Northrop Grumman Daimler Bombardier Naval US Navy Newport News Electric Boat Lockheed Penn State BBN DRA - UK DCN - France Hitachi Mitsubishi More... Sandia IBM Seagate Johns Manville Owens Corning Carrier GE Otis Daikin Mitsubishi
声学所考试材料
数字信号处理
(美)奥本海姆(Oppenheim,A.V.),(美)谢弗(Schafer,R.W.)著;董士嘉,杨耀增译《数字信号处理》,科学出版社,1980.12
第二单元
HNC理论
黄曾旸《HNC理论》(概念层次网络)清华大学出版社,1998
第三单元
海洋声学
布列霍夫斯基《海洋声学基础》海洋出版社,1985
电子系统设计[专著] /俞承芳主编;宋万年等编著,复旦大学出版社, 2004
第三单元
半导体物理与技术器件
黄昆、韩汝琦《半导体物理基础》,科学出版社, 1979
夏海良、张安康等编等《半导体器件制造工艺》[专著],上海科学技术出版社, 1986
第三单元
微电子机械系统
格雷戈里.科瓦奇著,张文栋译《微传感器与微执行器全书》,科学出版社,2003
《通信原理》[专著]樊昌信,曹丽娜编著,6版,国防工业出版社, 2006
第三单元
语音信号处理
杨行峻、迟惠生等编著《语音信号数字处理》,电子工业出版社, 1995
第三单元
分层介质中的波
《分层介质中的波》,布列霍夫斯基赫著,杨训仁译,科学出版社
第三单元
声波导
《声波导》,M.R.雷特全特著,严仁博译,上海科学技术出版社
声学所2012年博士考试科目及参考用书
考试单元
考试科目
参考用书
备注
第一单元
英语
请参见中科院研究生院招生信息网/
科学院统考
第二单元
声学基础
杜功焕,朱哲民,龚秀芬著《声学基础》南京大学出版社,2001.3
P.M.莫尔斯,K.U.英格特,《理论声学》(上、下册),科学出版社,1984
马大猷《噪声控制学》科学出版社,1987年
海洋声场预报及战术应用
文章编号:1672-7649(2005)04-0066-04海洋声场预报及战术应用收稿日期:2004-03-03陈治国,张 星,章新华(海军大连舰艇学院,辽宁大连116018)摘 要: 随着科学技术的飞速发展,对海洋声场做出合乎要求的快速精确预报显得尤为重要。
射线模型、抛物线模型、简正波等模型适合于不同具体的海洋声场环境,这些传播模型可以实现声场数值预报的实时性。
基于上述模型建立的海洋声场预报系统及对其战术应用进行的研究,对声呐的设计、使用及舰艇作战有着重要的意义。
本文在上述模型基础上利用Matlab和F ortran P owerstation工具建立了海洋声场预报系统,使用声场预报系统对典型的海洋环境,即海洋锋、涡漩的声场能量传播规律进行了研究,并提出了在反潜作战中对海洋锋、涡漩、深海会聚区和影区等海洋环境的战术应用。
关键词: 传播模型;声场预报;战术应用中图分类号: P733.2;E843 文献标识码: AOcean acoustic field prediction and tactical applicationCHE N Zhi2guo,ZH ANG X ing,ZH ANG X in2hua(Dalian Naval Academy,Dalian116018,China)Abstract: T oday,with the fast development of technology,it is im portant to predict the acoustic field fast and correctly.Ray m odels,parabolic equation m odels,normal m ode m odel etc are adapt to different of concretely matter and can realize real time calculation.These m odels can realize the ocean acoustic field prediction in spots, on the basis of these m odels,equipment for prediction and tactical application,there are significance for devising or using of s onar and combating of naval ships.On the basis of these m odels,this paper establish the ocean a2 coustic field prediction system by Matlab and F ortran P owerstation.Numerical transmission losses in ocean front, v ortex are calculated by the prediction system.The paper offer the tactical application of ocean front,v ortex and deep ocean convergence zone in the antisubmarine cam paign.K ey w ords: propagation m odel;acoustic field prediction;tactical application0 引 言声波是目前惟一能够在海水介质中进行远距离信息传播的有效载体,因此,研究海洋中的声波特性具有极其重要的意义。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
张海澜的《理论声学》答案 1.1有一质点振动系统,固有频率0f 已知,质量和弹性系数待求。
现在质点上增加已知质量m ,固有频率改变了'0f 。
求原质点系统的质量和弹性系数。
02f π=增加质量m , ()'2ff π-=可得0'00f f f=-,()()2'00''02m f f M fff-=-,()()222'0''042m f ffK fffπ-=-1.2一弹簧竖直悬挂,上端固定,下端系一质点组成简单质点振动系统。
质点同时受到向下的重力。
分析质点的振动和能量的转换。
22d x M M g Kx dt=--,22d x K x g dtM+=-,特解0M g x K=-,0x =弹簧受力为零通解()0cos a M g x x t Kωφ=+-,平衡位置M g x K=-,简谐振动速度()00sin a v x t ωωφ=-+动能:()2210011sin 22a E M v M x t ωωφ==+⎡⎤⎣⎦弹黄势能:()222011cos 22a M g E K x K x t K ωφ⎡⎤==+-⎢⎥⎣⎦()()()22001cos cos 22a a M g K x t M gx t Kωφωφ=+-++⎡⎤⎣⎦重力势能:()30cos a M g E M gx M g x t K ωφ⎡⎤==+-⎢⎥⎣⎦()()20cos a M g M gx t Kωφ=+-。
0x '=取为平衡位置,M g x x K'=+通解()0cos a x x t ωφ'=+速度()00sin a v x t ωωφ=-+, 动能:()2210011sin 22a E M v M x t ωωφ==+⎡⎤⎣⎦弹黄势能:()22011cos 22a M g M g E K x K x t K K ωφ⎛⎫⎡⎤'=-=+- ⎪⎢⎥⎝⎭⎣⎦()()()22001cos cos 22a a M g K x t M gx t Kωφωφ=+-++⎡⎤⎣⎦重力势能:()30cos a E M gx M gx t ωφ==+1.3火车以速度0v 运动,车上有一简单质点振动系统,弹簧与火车运动方向平行,一端固定在火车上,另一端连接的质点沿火车运动方向振动。
分析质点的振动和能量转换的规律。
2112220X M K KX K M K ωω⎫⎛-⎫⎛=⎪⎪ -⎝⎭⎝⎭ ()4212120M M K M M ωω-+=()()210ω=, ()()221212M M K M M ω+=。
()()111200x x x v t ==+, ()()2211220M X M X +=()()()()()()()()11221001112200cos cos M x x v t A t M x x v t A t ωϕωϕ⎧=+++⎪⎨⎪=+-+⎩(0.1) 能量()()()()212212101sin 2T M M E v A t M ωωϕ⎡⎤=-+⎢⎥⎣⎦()()()()()()()()()()21221021222222201sin sin22M v M M v A t A t M ωωϕωωϕ=-+++()()()()()()()()()()2122021222222220sin sin22T M v M E M v A t A t ωωϕωωϕ=++++()()()22212211cos 22V K K M M E x x A t M ωϕ⎡⎤+=-=+⎢⎥⎣⎦1M 比2M 大很多,略去21MM , ()()222K M ω=,振型10X =,21X =()()()()()()2110212120sin 2T M v E M v A t ωωϕ==-+()()()()()()()()()()2122021222222220sin sin22T M v M E M v A t A t ωωϕωωϕ=++++()()222cos2V K E A t ωϕ=+()10v 不同,在不同的参考系中,2M 的动能的变化是不同的。
同时1M 的动能也有变化,1M 的速度实际上不变1.4有一动圈传声器,当作质点系统处理,测得振膜固有频率600Hz ,质量0.8克,求弹性系数和力顺。
()20211370K f M π = =,力顺518.810C K-==⨯(米/牛顿,秒2/千克)1.5质点系统受力()111sin sin ,f h t t ωωωω=+ ,求振动。
()()11sin cos cos 2h f t t t ωωωωω= +--+⎡⎤⎣⎦()()()()()()()()12112211113211exp()sin exp 2cos 2exp 2cos 2i i t t x M i R Kh i t h t x M i R Khi t ht x M i R K ωωωωωωωωωωωωωωωωωωωω- -> =--+--⎡⎤⎣⎦- -> =----+--+⎡⎤⎣⎦-+ -> =-+-++1.6有下列形式的作用力作用于质点振动系统上,求振动位移。
()120,1,2,112aa F kT t k T f k F k T t k T⎧⎛⎫ ≤≤+ ⎪⎪⎪⎝⎭= =⎨⎛⎫⎪- +≤≤+ ⎪⎪⎝⎭⎩把外力展开为傅利叶级数122cossinn n n n nt nt f a b TTππ∞∞===+∑∑0010Ta fdt T==⎰,022cos 0Tn nt a f dt TTπ⎫⎛==⎪⎝⎭⎰0 22sin 4 Tn a n nt b f dt F TTn n ππ⎧⎪⎫⎛==⎨⎪⎝⎭⎪⎩⎰偶奇1,3,5...1,3,5...441212sinexp aan n F iF nt i nt f nTn T ππππ∞∞==⎫⎛=→- ⎪⎝⎭∑∑221,3,5...212exp 2an TF i ntx n T n Z T πππ∞=⎫⎛=-- ⎪⎫⎛⎝⎭⎪⎝⎭∑()0221,3,5...2102an TF x x n n Z T ππ∞===-⎫⎛ ⎪⎝⎭∑()01,3,5...4102an i F v v n nZ T ππ∞===⎫⎛ ⎪⎝⎭∑瞬态部分,初始位移和初始速度为0x -和0v -1.7画出如右图所示弹簧并联系统的导纳型类比线路图,求出系统的等效弹性系数。
2mK C =并1.8如下图所示隔振系统,系统弹簧置于阻尼材料中,力阻1m R ,画出系统的导纳型类比线路图,分析1m R 的作用。
恒流源 1112m m mm V V V V F i C G V F i M V G ωω--⎧=+⎪-⎪⎨⎪=-+⎪⎩得到:()1121122111111m m m m m m m m mm m m V i R K i C G V M i R R K i M G i C G ωωωωωω+-+-==--++-+++-习题2.1 考虑如图3个相同的质点和4个相同的弹簧组成的一维的耦合共振系统。
求其共振频率和模式,验证模式的正交性。
0ω=()()()()1112221233332=0 00 M x K x K x x M x K x x K x x M x K x K x x ⎧++-⎪+-+-=⎨⎪++-=⎩22200222200022200/20/2/200/2ωωωωωωωωωω⎫⎛-⎪ -=⎪ ⎪-⎝⎭()10ωω=101⎫⎛⎪⎪⎪-⎭,()1,20ω=,模式是1121⎫⎛⎪ ⎝⎭正交:质量相同(111110021⎫⎛⎪=⎪⎪-⎭,(111110021⎫⎛⎪=⎪⎪-⎭,()111110221⎫⎛= ⎪ ⎝⎭习题2.2 考虑第一个例子中的两个质点三个弹簧组成的一维振动系统,两个质点的质量相同,三个弹簧的弹性系数相同。
两个质点分别受力()1f t 和()2f t ,试用简正方法写出运动方程,并和解非齐次方程组的方法的结果比较。
()111222120M x K x x Kx ⎨+-+=⎩22202M K Kk M K ωω⎫⎛-=⎪-⎝⎭()1ω=11⎫⎛⎪-⎝⎭,()2ω=,模式是11⎫⎛⎪ ⎝⎭()()1121221222M x Kx Kx f t M x Kx Kx f t +-=⎧⎨+-=⎩两式相加和相减()()()()()()112212a t x t x t a t x t x t =+⎧⎨=-⎩ ()()()11121M aKa f t f t q t +=+= , ()()()221223M aKa f t f t q t +=-= 傅里叶变换()()()1exp 2x t X i t d ωωωπ∞-∞=⎰,…112Q A M K ω=-+,2223Q A M Kω=-+()1211222123Q Q X A A M KM Kωω=+=+-+-+()()()()()()22121222323M KFF M KFF M KM K ωωωω-+++-+-=-+-+()()21122223M F KF KF M KM Kωωω-++=-+-+()1221222123Q Q X A A M K M Kωω=-=--+-+()()()()()()22121222323M KFF M KF F M KM K ωωωω-++--+-=-+-+()()22122223M F KF KF M KM Kωωω-++=-+-+用解非齐次方程的方法解()112122122M x Kx Kx f t ⎨+-=⎩()()2121212222M K X K X F K X M K X F ωω⎧-+-=⎪⎨-+-+=⎪⎩()()12221121222222322F K F M KM F KF KF X M KM KM KK KM Kωωωωωω--+-++==-+-+-+---+()()21222122222222322M KF KF M F KF KF X M KM KM KKKM Kωωωωωω-+--++==-+-+-+---+习题2.3 考虑如下图质点被两个相同的弹簧固定,弹簧中的张力是T ,质点可以在三维空间中运动,位移很小,分析其振动。
质点位置xyz ,平衡位置是原点,弹簧原长0l 平衡位置处弹簧长度0'T l l K=+,K 弹性系数固定点的坐标()',0,0l ±质点位移后弹簧长度1,2l =保留一阶小量1,2'l l x =≈± 质点在x 方向的受力()()102012''x l l x K l l K l l l l +---+-()12x K l l K x ≈--≈-质点在y 方向的受力()()()102001222'''r y y y T K l l K l l K l l y K yl l l l ----≈--≈-=-质点在z 方向的受力()()1020122'r z z T K l l K l l z K z l l l ----≈-=-2x K K=,2'r T K l =运动方程 000x r r M x K x M y K y M z K z +=⎧⎪+=⎨⎪+=⎩三个方程是解耦的,就是简正坐标 频率方程 ()()2220x r M K M K ωω--=()1ω==x 方向的振动,振动与张力T 无关()()23ωω===yz 平面内的椭圆振动,yz 平面内的任意方向都可以作为简正方向。