(范文)海洋锋区声传播数值仿真研究

1、海洋鋒区声传播数值仿真研究the numerical simulation of sound propagation in oceanfront zone学生姓名学号所在学院海洋与气象学院班级所在专业海洋技术申请学位理学学位指导教师职称讲师副指导教师职称答辩时间摘要abstractii1 引言12海洋锋的定义及其对海水声传播的影响23锋区声场仿真模型33.1 海洋锋特性33.2 声学计算模型简介43.3 海洋锋模型64 锋区声场仿真74.1 海洋锋区对不同声波频率的影响74.1.1 模型参数设置74.1.2 仿真结果及讨论85 总结10鸣谢11参考文献12摘要海洋锋区的海洋环境参数在水平方向和

2、竖直方向的梯度变化对海洋声传播有重要 的影响。环境参数在水平方向的梯度变化，特别是海水温度的水平梯度变化，会导致 海水声速的改变。声速的变化又会引起声线结构、传播损失、传播距离的改变。本文 通过建立海洋锋声场模型，应用ram模型來研究海洋锋区声传播特性。从仿真实验中, 观察到以下现象：声波穿越海洋锋后，海洋锋区会出现一个能量折射的现象，使得部 分能量在原有的传播方向上发生改变；海洋锋对不同频率的声传播影响不同，在经过 海洋锋区吋，高频声波发生的能量折射现象比低频声波的更加强烈。关键词：海洋锋；声传播；ram模型。abstractocean fronts have a great effecto

3、n acoustic wave propagation for it generates horizontal and vertical inhomogeneity atsound speed, which result in significant fluctuations of sound energy.firstly, according to the characteristics of ocean fronts, this paper used an ocean front model to simulate acoustic velocitywhere ocean fronts h

4、appen. secondly, the paper utilizes ram to produce 2d pictures of acoustic transmission loss in the area of ocean front.the following phenomena could be observed: the energy of sound wave will generate refraction when it goes through the area of ocean front:. this impact on the high-frequency sound

5、waves is more significant.keywords: ocean front; sound propagation; ram海洋锋区声传播数值仿真研究海洋技术，201212922229,曾志杰指导老师：屈科1引言声波是人类迄今为止已知唯一能够在水中远距离传播的能量形式，因其在在海水 中的良好传播特性使其成为人们进行海洋研究的主要信息载体，两次世界人战的爆发 促进了水声技术的发展，人们逐步认识了声波在海水屮传播的机理。海洋环境各种复 杂的因素都会对声信号的传播产生影响。例如，潘长明等人用声学调查实测传播损失 数据和umpe (the university of miami pa

6、rabolic)模型仿真研究东海近海温跃层对水声 场的影响，发现温跃层对低频声场较高频声场影响大，夏季正跃层海洋环境中传播损 失随深度增加而减小，春季逆跃层海洋环境屮传播损失随深度增加而增大；张旭等 人(2014)应用kraken简正波模型对实测中层冷水环境下的声场进行计算并分析， 发现声波在冷水层中传播时形成波导声场使其传播损失较常规环境明显减小。此外, 海洋环境对水声传播的影响还包括水声环境效应的不确定性，主要体现在：1、边界 条件变化引起的不确定性，包括海浪、潮汐运动、大洋环流潮流等引起的海面变化和 地质运动造成的海底地质、地形发生的改变；2、海洋水体中由于某种特性引起的不 确定性，如温

7、跃层的厚度变化；3、传播过程引起的不确定性，主要包括声波在海而、 海底和海水中的折射反射造成的多路径传播。海水中的声传播速度是最重要的声学参数，是影响声在海洋中传播的最基础的物 理量。在流体介质中声波是弹性纵波，它可以表示为c = m网。其中，目是绝热压缩 系数，p是海水密度，它们都是海水温度、盐度和静压力的函数。所以，海水声速是 温度t (c° )、盐度s (%。)和静压力p (kg/cm2)的函数，根据乌德公式(1)对上式进行温度、盐度和压强求导，可得qede= = 4.z110,0074?= = 1.14 a« = = 0.1751 qt，'卅5, b qp由

8、此可知，温度对声速的影响最大，每升高1c。，声速变化可达到4m/s；盐度每 增加1%0,声速增加1.14m/s；压强每增加1 atm,声速仅增加0.175m/s 4海洋锋区的诸多环境参数会发生急剧变化，这会直接影响海水的声学特性。海底、 海表的不均匀性以及海水温度和盐度的铅直分层特性是引起声场起伏的重要因素，此 外，气泡、冷暖水体、湍流和大洋屮的浮游生物和游泳动物群等对海洋声场的变化也 是实际中不可忽略的。例如，在声传播的过程中，气泡会对声波散射，其作用相当于 一层屏障；冷、暖水体对会使声传播吋发生折射；湍流使海水的温度、盐度产生局部 变化从而引起声速的局部变化；内波与远距离传播声信号的振幅和

9、相位有着密切关系。 但在研究海洋锋区声传播时，为了简化研究问题，需要将气泡、湍流等因素忽略，把 海洋锋区的温度、盐度和密度变化转化为声速的变化，建立海洋锋声场模型。通过观 察海洋锋的特性，获得海洋锋的重要参数如强度、尺度、位置，本文在此基础上进行 数值仿真，并使用声学模型进行计算、分析结果。2海洋锋的定义及其对海水声传播的影响在人类开发和利用海洋屮，各种以声波为主的海洋探测设备相继问世，例如，声 波测深仪和冋声鱼探仪、adcp地震剖面仪、侧扫声呐。然而作为声信息传输通道的 海水介质及其边界条件存在多种非均匀性，使得声波在海水介质中的传播规律复杂多 变。研究这样复杂多变的传播规律对水下探测和通信

10、系统的设计与使用是十分必耍的。 海洋锋、中尺度涡旋和内波是引起海洋水平非均匀性的典型现象2，从对声传播的 影响上看，中尺度涡旋可以看做是一种特殊类型的海洋锋。海洋锋一般是指性质明显不同的两种或几种水体之间的狭窄过渡带。根据空间 分布出现阶跃的海洋环境参数的不同，可将其分为温度锋、盐度锋、密度锋、声速锋、 水色锋等；从形成机制上可分为：沿岸流锋、河口羽状锋、浅水陆架锋、上升流锋和 强西边界流锋；从锋生尺度的大小上可分为：行星尺度锋（如南极锋、北极锋等）、屮 尺度锋（陆架坡折锋）和小尺度锋（局部海域的出现的锋）。海洋锋的持续时间较短时 仅数小时，较长时可达数月或更长。rtr丁锋带附近具有强烈的水平

11、辐合或辐散以及铅直方向的运动，因而是不稳定的, 其中存在着逐渐变性的过程和尺度的弯曲。通常情况下，海洋锋两侧会出现较大的 温度、盐度梯度，表现在如下声学异常叫生物的活动频繁增强了环境噪声和体积混 响，在墨西哥湾流中出现中等强锋时，其环境噪声可增大5-10db；声波在通过锋区 吋会产生不同程度上的反射或折射，增大能量损失从而导致距离减小和观测方位产生 误热声在海洋锋区声传播吋，由于水平方向上海水的温度、盐度和密度有着明显的水 平梯度变化，声速在水平方向上也会冇显著的梯度变化，使得声传播的方向和距离发 生改变。所以海洋锋区会对声呐的探测与反探测、水下通信和水下探测应用等产生重 要影响。研究海洋峰区

12、的声传播特性，冇助于海洋声传播理论更好地被完善。3锋区声场仿真模型3.1海洋锋特性为了弄清楚海洋环境对声传播的影响，众多国内外学者开展了许多声学实验，并 得出了重要的结论。weinberg和clark (1980)在实验屮观察到声传播在水平方向 上的折射角度可以超过1° 7o abbot等人(2001)曾在韩国沿岸区域开展了声学实验, 结杲表明不同的海底界面可以导致5-20db的声传播损失，不同的声源深度可以导致 5-15db声传播损失，不同的锋面结构可以导致5db的生传播损失。南明星，杨廷武，风雷(2003)用海洋锋模型仿真海洋锋，使用由美国国家与大 气管理局研制并开发的 harp

13、o (hamiltonian acoustic ray-tracing program for the ocean) 模型进行计算，得出了以下结论：海洋锋区会使声线发生水平偏转，水平偏转角度随 着距离的增加而变大，最大值达到1.0。左右叫他们还根据黑潮锋结构和分布特点模 拟黑潮海洋锋，认为台湾黑潮锋区会对声呐定位精度产生影响。lynch等人(2006)通过实验并得出以下结论：海洋小尺度运动(如内潮、内波) 会导致5-20db的传播损失冋。宋俊等人(2005)对我国东海浙闽沿岸锋和苏北沿岸潮 生锋进行声速剖面建模，采用motal (绝热简正波)为声场计算模型，仿真计算了低 频声信号(500hz)

14、穿越不同垂直结构、不同宽度和不同垂直结构的锋区传播损失情 况，仿真数据表明：浙闽沿岸锋对低频声传播影响很小，苏北沿岸潮生锋的影响却很 大；苏北沿岸潮生锋的宽度和强度不同引起的声传播损失差别不大，但它的层化水区 域的温度垂直梯度不同值引起的传播损失差别很大，他们对此现象的解释为，声线穿 越锋区后向下折射并且和有吸收衰减的海底发生多次反射引起的菅永军等人（2006）利用黑潮锋区断面的实测ctd数据获得锋面声速分布，通过 确定海洋锋面的宽度、位置、和陡度等特征参数，建立了与实测数据较符合的海洋锋 区声速计算理论模型，使用二维抛物方程模型研究海洋锋面的声传播特性。结论显示: 锋面对声传播损失的影响受到

15、声源位置、发射频率和接收深度等的影响，声源在50m 和发射频率/=150hz在不同接收深度,有锋面吋和无锋面吋传播损失最人相差20db,2o张旭等人（2013）运用bellhop水声数值模型分析东海黑潮锋环境下的声场特性, 结果显示：声线（声波频率1khz）穿越锋面发生干涉吋会使声能强弱空间分布出现偏 移，并且锋而引起的非均匀水文环境是引起声能空间异常变化的关键因素问。3.2声学计算模型简介影响声场的海洋环境因素有很多，为了研究不同海洋环境因素对声场的影响，人 们从实际的声传播问题中抽象出声传播模型，主要有：简正波模型、射线模型、抛物 方程模型、多途扩展模型、快速场模型、射线简止波模型、简止波

16、抛物方程模型和 射线简正波抛物方程模型，这些模型都符合最基本的波动方程，式（2）。这些模型 冇各自的优缺点和适用范围。射线理论适用丁处理高频声场的问题，口前主要的射线模型有：ray、bellhop. trimain和harpo等。它是建立在高频近似的基础上，通过对不同频率的衍射修正， 它也可以应用到较低频声场。射线理论既适用于海洋环境与距离无关的声传播问题， 也适用于海洋环境与距离有关的问题，但是对于求解非水平分层介质的问题上，它的 求解过程比较繁琐。简正波理论是海洋声学中比较完善的计算方法，得到了广泛的应用。目前很多已 经冇很多简正波计算程序已被开发，kraken、moatl、snap、co

17、uple等。简正波模 型是基于与距离无关的假设，通过"绝热近似和“模式耦合”方法可以将其扩展为与 距离冇关的模型。绝热简正波的优点是计算量小，计算速度快，但因为它的基本假设 是简止波保持绝热耦合，忽略了简止波z间的能量交换，所以它的适用范围局限在水 平方向变化缓慢的问题上。耦合简正波在处理二维、三维非均匀波动问题上，它的精 确程度比较高，它的结果可以作为其他检验其他模型的标准解。耦合简正波成像清晰、 精度高，而它的缺点是计算速度慢。抛物方程模型是针对水平变化海洋环境屮声场计算的一种快速近似算法，是波动 方程的窄角近似解，在水平较窄的角度范围内能计算的传播角度小于15。，后来发展 到了

18、可处理宽角、后向散射等问题的算法。由于模型假设上的限制，pe模型不能计算 水平变化比较剧烈的声场，但是它在处理水平变化和三维变化的声场的问题上具有优 越性。抛物方程主要使用了四种数值技术：1分裂步进傅里叶算法；2.隐式有限差分 (ifd)； 3常微分方程(ode)； 4有限元方法(fe)。多途扩展模型具冇与射线模型相同的部分性质，它通过积分无穷结合把波动方程 展开为声场积分表达式，每个声线路径可由一个积分表达式表示，它能正确估计焦散 区和声影区的声压场。主耍的程序模型冇fame和mule。快速场模型使用简止波近似方法来分离波动方程参数，并通过快速傅里叶变换算 法对变换结果进行数字估计。口前主要

19、的理论模型有：ffp、press、scooter、sparc 和 oases o射线简正波模型可以将简正波赋予射线的含义，以此得到相应的物理图像。主要 的射线简正波模型有广义相积分简正波模型(wkbz)和波束位移射线简正波模型 (bdrm)o简正波抛物方程方法(pmpe)是简正波和抛物方程方法的结合，包括绝 热简正波抛物方程方法(ampe),和耦合简正波抛物方程方法(cmpe)o射线简正 波抛物方程方法(rmpe)是建立在射线简止波模型和简止波抛物方程方法模型的基 础上建立的，在垂直方向上采用射线简止波模型进行简止波分析，在水平与方向角方 向采用抛物方程方法来求解简止波幅度系数方程。本文主要使

20、用抛物方程ram模型来计算仿真声场中的传播损失。nspe仃he navy standard parabolic equation)是研究水声传播算法的重耍成果，它是基丁抛物方程，用 于解决水声传播并且能够满足使用者对输入输出需求的算法。nspe声学模型的出现主 要是为了满足对深水域低频声传播的需要。在早期的几个版木里，ssf(spilt-step fourier solution method)结合了输入参数，能让操作者输入实测的海底数据。它可以模拟粗 糙的海表或者被冰覆盖的海表，能够仿真带宽和波束接收接通系统。nspe继承了 ssf 算法的假设：声传播媒介中没有垂直不连续性。但是在20世界9

21、0年代，海军研究的 重点放在了浅水，这就使得海底对声传播的作用更加明显，这使ssf的海底界面和液 态半无限界面近似不再起作用。幸运的是这个时候研究有了新的进展，产生了儿种能 在数值上处理垂直不连续性模型，提供了可靠地浅水预测。ram (range-dependent acoustic model)是认可程度比较高的模型,它是tl michael collins研究出来的,它使 用分步pade近似求解方程，这使得计算速度和计算精度比以往的pe模型冇了较人的 提高，得到广泛的应用。ram是发展的非常快的一个研究模型，因为它的求解允许步 长范围比声波波长大，并且不需要细微垂直网格。它的网格可以设置的

22、非常精细并且 计算速度能够达到军事环境的要求，它还允许并行计算从而使得计算速度犬犬提高。 所以，ram模型是漫长的科研探索z路上一个重要的里程碑z-o3.3海洋锋模型海洋锋三维模型表达式：(3)上式屮，g©代表求坐标位置；切4叭卩)指不存在海洋锋时各个位置点的声速；/指声扰动因子，控制着海洋锋的强度；、弘、如指海洋锋的中心位置，陀、叫、 指海洋锋的范围。将上式简化为二维模型悶：(4)式(4)屮，r, z是声场屮坐标位置点，wr是海洋锋在水平方向上的半径，wz 是海洋锋在铅直方向上的半径。co是在没有海洋锋时声场位置点的声速，本文把它设 置为munk深海声道声速剖而的声速，如图3-1黑

23、线。声速（m150015101520153015401550(a) munk深海声速剖面距离(km)15451540-1535-1530115251520151515101505图3-1(b)海洋锋声速场4锋区声场仿真4.1海洋锋区对不同声波频率的影响4.1.1模型参数设置在木次模拟实验中，木文把海洋锋中心设在水平距离60km、水深1000处，海洋 锋水平尺度为14km,垂直尺度为2000mo声扰动因子代表海洋锋强弱，本文将它设 为0.02o海底参数设置如下：海底声速为1600m/s,海底密度为卩二九*弓怦，压 缩波衰减系数为0.4db/xo声源深度为500m。为了观察不同频率声源在海洋锋的影

24、 响，本文将声源依次设为3oohz、400hz、500hz。在没有海洋锋吋，假设声速在水平方向没有梯度变化，垂直声速变化如图3-1 (a) 黑线所示。在深度1000米以内时，因为声速主要受到温度的影响，所以在深度为1000 米以上时，声速随着深度的增加而减小；在深度为1000米以下时，静压力的变化对 声速的影响较温度要大，所以在1000米以下的深度，声速随着静压力的增人而增人。 存在海洋锋时，声场如图(b)所示，芦速等值线向海洋锋中心位置弯曲，水平方向声速变化率上随着海洋锋中心位置的距离增人而减小。图3畀(a)蓝色部分绘制出 相同深度声速在水平方向上的变化范围，反映出存在海洋锋时，不同深度声传

25、播深度 的变化大小。1004060 80距离(km)200400e 600800 送10001200140020(a) f=300150200400-100 i 旦 600 姿800 送50 10001200-100150100120100200400e 600 弐800 炭100012001400204060 80(c) f=400150距离(km)80100140012020406080100120<b)仁300距离(km)m00500200400e 600 泾800 瑟100012001400150距离(km)(d)仁400图424.1.2仿真结果及讨论图42中(a)> (c

26、)、(e)是不同频率芦波300hz、400hz和500hz在没有海洋锋 时海水声传播的传播损失图，(b)、(d)、(e) 300hz、400hz> 500hz是在海洋锋区声 传播的传播损失图，蓝色虚线代表海洋锋的水平尺度。由图(a)、(c)、(d)可以看不同频率在海洋中的传播损失不同，高频声波在海洋中能量衰减的更快。rti (b)、(d)、 (f)可以看出不同频率声波经过海洋锋能量变化大小不同。但是他们都自一个共同点, 经过海洋锋后，能量发生折射。通过图42中有无海洋锋情况下的传播损失图对比， 可以看到在海洋锋区(蓝线区域)，高频声传播损失会比低频声传播损失减小的多，即 海洋锋区对高频声

27、波的能量折射作用比低频声波强烈。通过(c)、(d)图对比，可以明显的看到，能量会沿着海洋锋发生一个折射，在 海洋锋区能量从上层(800m以上)折射到下层(800米以下)，从而使上层能量减小, 下层能量增加。主耍体现在声波穿过海洋锋区后，即第二根蓝线z后，从颜色变化上 可以看出(d)图上层较(c)图上层传播损失明显增加，而下层传播损失减小。经过(a)、(b)仔细对比后，可以观察到经过海洋锋后，(b)图在深度为0500米范围内 的颜色较(a)图浅，而在500-1000米范围内的颜色(b)图较(a)图深，所以还是图4-3能看见能量分布上发生的变化。尽管图42可以看出不同声波在海洋锋传播时发生的能量折

28、射现象，为了更直观 的观测处海洋锋对不同频率声波的影响，木文给出了 500米处的传播损失图，这个深 度与声源深度相同。如图43所示，（a）、（b）、（c）是声波频率分别为300hz、400hz、 500hz在深度为500米处的传播损失图。两条黑色虚线代表海洋锋的水平尺度，绿色 实线和黄色虚线分别代表着不存在海洋锋时的传播损失和存在海洋锋时的传播损失。第一根黑色虚线前，黄色虚线和绿色实线重合，说明在没有经过海洋锋区域时两 种情况的传播损失是相同的，这也是模拟试验中必然的结果。但是在进入海洋锋以后, 黄色虚线偏离绿色实线，从数值上看说明了经过海洋锋以后传播损失增大。观察不同 频率（300hz、40

29、0hz和500hz）在有无海洋锋环境下的传播损失，可以清楚的看到， 在300hz声波经过海洋锋时，黄色线和绿色线偏离的程度比较低，它的传播损失增加 的相对较少，增加的幅度在020db。频率为400hz的声波经过海洋锋口寸，传播损失相 对300hz的声波增加，增加的幅度主耍在1020dbo 500hz的声波经过海洋锋时，传 播损失增加的主要幅度也是在10-20db,但是因为部分能量发生了重新分布，使得锋 区部分位置上的传播损失比无锋时更小。5总结本文通过建立不同尺度的海洋锋声场，使用ram模型对其求解，并画出了它们二 维声传播损失图和在深度为500米处（与声源深度相同）的声传播损失图，并得出以

30、下结论：声在海洋锋区传播时，会出现一个能量折射的现象，使得部分能量在原冇的传播 方向上发生改变；海洋锋对不同频率的声传播影响不同，在经过海洋锋区吋，高频声 波发牛的能量折射现象比低频声波的更加强烈。鸣谢本论文是在屈科老师的悉心指导下完成的。从论文选题、文献查找、算法模型的 使用和论文的修改屈科老师都给予我很多的帮助，他洞悉前沿科学的求新精神和一丝 不苟的研究精神使我受益颇深。rtr丁学生知识水平的限制和经验的不足，在完成论文 的过程中遇到了不少的疑问和困难，每次向屈科老师求助时，他总能够在繁忙的工作 屮挤出时间向我讲解，并给我提出很多建议，让我更好的完善整篇论文。在此向屈科 老师表示衷心的感谢

31、。在选题z前，朱凤芹老师也给了我许多宝贵意见，学生在她那获得不少的支持和 鼓励，在此致以诚挚的谢意。光阴往苒，大学时光转ii舜即逝。我要感谢海洋与气象学院每一个辛勤工作的老师, 谢谢你们给我的关心和帮助，谢谢你们在犬学这段吋间里对我的栽培和教育。我要感 谢曾在这段时间里给我帮助和理解的同学，谢谢你们陪伴我度过大学的美好时光。最后，我要感谢我的一直以来给默默给我支持的父母。正是他们的辛劳付岀才使 我完成了学业，是他们让我能够披荆斩棘勇往直前。参考文献1 潘长明，高飞，孙磊，王璐华，王本洪,李璨.浅海温跃层对水声传播损失场的影响.哈 尔滨工程大学,20142 张旭，程琛，刘艳.中层冷水环境卜的声场特性分析.海洋科学进展，20143