基于FIO-COM的海洋声学预报系统的构建与应用.pdf

基于 FIO-COM 的海洋声学预报系统的构建与应用杨春梅1,2,3,4，刘宗伟1,2,3,4，姜莹1,2,3,4，吕连港1,2,3,4*，肖斌1,2,3,4（1.自然资源部 第一海洋研究所，山东 青岛 266061；2.自然资源部 海洋环境科学与数值模拟重点实验室，山东 青岛 266061；3.山东省海洋环境科学与数值模拟重点实验室，山东 青岛 266061；4.崂山实验室 区域海洋动力学与数值模拟功能实验室，山东 青岛 266237）摘要：将海洋模式与声传播模型结合在一起，设计开发了一种适用于高性能计算机的全球海洋声学预报系统 FIO-GOAFS，该系统以自然资源部第一海洋研究所全球 0.1分辨率海浪-潮流-环流耦合模式（FIO-COM）为基础，利用海洋模式预报的温、盐、深参数计算声速剖面，并对声速剖面进行水声环境特征诊断，之后将海洋模式与水下声场传播模型协同连接，结合地声模型（海底地形和底质参数），实现了全球海域的水声环境特征诊断及水下声场及相关结果的预报。海洋模型提供水下声学预报所需的水体声速、海浪波高等参数，地声模型提供海底地形、底质声速、密度以及衰减等参数，通过调用海洋-声学连接模块提取声传播路径的地形及海洋环境参数剖面，实现海洋模型和声学模型的有效连接。全球海洋声学预报系统在高性能计算机上并行实现，主要包括声场计算中的频点、方位角并行以及声学预报时针对地理空间区域的并行。最后，利用该系统预报并分析了全球海域的水声环境特性及声呐作用距离的季节变化和空间分布特征，为现代声呐的设计、操作和水下应用提供参考。关键词：声学预报；海洋模型；声传播特性；声道；水声环境特性中图分类号：P733.2 文献标志码：A 文章编号：1671-6647（2023）03-0510-10doi：10.12362/j.issn.1671-6647.20220428001引用格式：杨春梅,刘宗伟,姜莹,等.基于 FIO-COM 的海洋声学预报系统的构建与应用J.海洋科学进展,2023,41(3):510-519.YANG C M,LIU Z W,JIANG Y,et al.Construction and application of ocean acoustic forecastingsystem based on FIO-COMJ.Advances in Marine Science,2023,41(3):510-519.近几十年来，各国对水声设备、水声理论、海洋环境等方面的研究提出了更高的要求。声呐技术取得跨越式发展的重要途径在于深入挖掘海洋环境水声特性。因此，研究海洋声速场及声场的分布、了解声速剖面及声传播的时空分布变化规律、掌握声道的分类特征，以及探索深浅海声道特征量的地理分布状况等，对提高反潜、探测能力和建立深海预警体系有着非常重要的意义。海洋声学的研究范畴主要集中在：探索以海面波浪、海水非均匀性，以及海底结构为代表的海洋环境在时空变化方面对声场作用的规律；研究声波在海洋探测和反演等方面的应用，即如何 收稿日期：2022-04-28资助项目：崂山实验室科技创新项目（LSKJ202204800）；国家自然科学基金项目（42206195）；山东省自然科学基金项目（ZR2020MA090）作者简介：杨春梅（1984）,女,助理研究员,博士,主要从事海洋声传播特性及声场建模方面研究.E-mail:*通信作者：吕连港（1974）,男,研究员,博士,主要从事物理海洋方面研究.E-mail:（胡筱敏编辑）第 41 卷第 3 期海洋科学进展Vol.41No.32023 年 7 月ADVANCESINMARINESCIENCEJuly,2023利用声波来探测海洋结构及海洋中物体的位置与特性1。在以往进行的海洋环境对声场影响研究中2-4，大多对海洋环境进行相对简单的假设，这必定会对分析结果造成一定偏差。随着时代与科学的进步，很多优秀成熟的海洋模式被开发，其研究预报对象也能分辨中小尺度海洋运动，可以与声场理论结合进行研究。近 20 a 来国外在这方面得到了极大的发展5-10。Botseas 等将隐式-有限差分 PE 模型与 HOOM 模型产生的海洋预报联系起来，实验结果成功地证明了 这 种 耦 合 方 式 的 有 效 性。Mellberg 等11于 1990 年 提 出 使 用 海 洋 描 述 预 报 系 统（the OceanicDescriptive Predictive System，ODPS）与有限差分模型耦合建立中尺度涡条件下的声场参数，对湾流区域冷、暖涡旋影响下的声传播特点以及产生的声传播时间的影响进行了深入的分析与总结。2016年 Heaney 和 Campbell12使用海洋预报数据库 ECCO2（Estimating the Circulation and Climate of the Ocean,Phase II）结合全三维抛物模型系统研究了海洋中尺度现象对全球尺度上的低频声传播影响，发现中尺度现象能导致 1.8的水平偏转，这可能对远距离低频声源的定位产生重大影响。国内在声学和海洋学的耦合研究方面也有了一定发展。笪良龙等13 利用海洋-声学耦合模式预报声速场时空变化，获取声速垂直结构不确定性分布规律，提出了经验正交函数-随机多项式展开方法以降低不确定参数维度，得到声场不确定性分布。肖瑶14利用遥感数据与全球再分析数据，结合涡旋探测和追踪方法分析我国南海北部和墨西哥湾中尺度涡旋的参数特征和分布情况，基于简正波和射线理论对中尺度涡引起声场时空变化机理进行了解释。自然资源部第一海洋研究所基于多年的积累，开发了全球 0.1分辨率海浪-潮流-环流耦合模式（Wave-Tide-Circulation Coupled Ocean Model by the First Institute of Oceanography，FIO-COM）。本文在该模式的基础上，构建了一种新型的适用于高性能计算机的全球海洋声学预报系统，该系统集合了一种水声环境特征诊断方法和几种声传播模型，使其具有海洋环境水声特征诊断和海洋声场相关预报功能；并将该系统应用于全球海域，分析了全球海域的水声环境和声呐作用距离的空间特征及季节变化。1海洋声学预报系统的构建和功能基于自然资源部第一海洋研究所自主研发的全球 0.1分辨率海浪-潮流-环流耦合模式（FIO-COM），本文提出了一种水声环境特征诊断方法，将其与海洋动力模式、水下声传播模型协同运作，构建了一种适用于全球海域的水下声学预报系统（Global Ocean Acoustic Forecasting System，FIO-GOAFS），该系统能够针对全球海域的水声环境进行声道特征诊断，并提供水下声场及相关结果预报。水声环境声道特征诊断结果包括表面声道深度（Sonic Layer Depth，SLD）及表面声道截止频率（Surface Duct Cutoff Frequency，SFD）、半声道条件（Half Channel Conditions，HAF）、深海声道轴深度（Deep Sound Channel Axis Depth，DSC）和浅海声道轴深度（Shallow Sound Channel Axis Depth，SSX）；水下声场及相关预报结果包括：传播损失断面分布、传播损失水平分布、汇聚区距离和声呐作用距离等。海洋声学预报系统功能模块的基本结构见图 1。本系统海底地形水深文件数据采用通用大洋水深制图（The General Bathymetric Chart of the Oceans，GEBCO）提供的 30分辨率全球大洋水深网格数据，该数据基于多波束数据插值而成，其中包含了源标识数据（Source ID-entified code，SID），用来标识真实原始声呐点水深数据和插值数据。海底底3 期杨春梅，等：基于 FIO-COM 的海洋声学预报系统的构建与应用511BOTSEAS G,LEE D,SIEGMANN W L.IFD:interfaced with Harvard open ocean model forecasts.United States:NavalUnderwater Systems Center Technical Report 8367,1989.质参数一部分利用自然资源部第一海洋研究所自主调查的海底环境参数资料，一部分采用反演获得的底质参数数据。海洋环境数值预报系统使用自然资源部第一海洋研究所自主研发的全球 0.1分辨率海浪-潮流-环流耦合数值预报系统（FIO-COM），该预报系统能够提供未来 7 d 全球海浪波高、波向、海表面高度、谱峰周期、海流流向和流速、海水温度、盐度、密度、水位和海冰密集度等海洋环境要素的预报产品，预报结果每隔 3 h 输出 1 次，其中海浪波高、温度、盐度、密度以及海冰密集度等均可转化为特定格式的海洋环境参数输入到声学模型中。海洋-声学连接模块通过将海洋模型的温度、盐度、密度、海浪波高等输出转化为声学模型的参数输入，将海洋模型和声学模型结合在一起。海洋声学传播解决模块集合了一种海洋环境水声特征诊断方法以及几种声传播模型，可用于诊断预测海洋环境水声特性以及声学的相关结果预报。本系统在高性能计算机上并行实现，在进行声场计算时，针对频点和方位角按进程个数进行均衡分配；在进行水声环境特征诊断及声呐作用距离预报时，将地理空间区域按进程个数均衡划分。本文针对全球海洋区域实现了水声环境特征诊断及声呐作用距离预报，并对预报结果进行了分析。参数初始化地形提取产品制作模型选择并行计算海洋环境参数提取地形数据海洋环境数值预报系统海洋环境预报模块海洋-声学连接模块海洋声学传播解决模块传播损失断面分布传播损失水平分布汇聚区距离声呐作用距离表面声道深度表面声道截止频率半声道条件深海声道轴深度浅海深度轴深度图 1 海洋声学预报系统的功能模块图Fig.1 Diagram of the functional modules of the Ocean Acoustic Forecasting System 2水声环境特征诊断及预报结果分析海洋中声速剖面的某些特定的垂直结构，能够限制声能量在一定深度范围内很少向外泄漏，声能量在该深度范围内远距离传播的这个水层结构称为声道。除高纬度、赤道等特殊区域外，典型深海声速剖面可分为三层结构：表面等温层、温跃层和深海等温层。极地区域深海一般不存在温跃层，声速从海面到海底呈现正梯度分布。水声学的浅海深度一般为几十米到数百米，其间没有固定的深海声道。海洋这些特定的分层现象及其产生的不同声传播模式与声呐的工作原理密切相关。水声环境特征诊断主要针对水体中的声速剖面进行。海洋模型提供全球海域的温度、盐度及对应的深度等环境参数预报，根据这些海洋环境参数求取对应网格点的声速剖面，本系统中声速剖面利用 Mackenzie 公式求解：c=1 448.96+4.591t5.304102t2+2.374104t3+1.340(S 35)+1.630102D+1.675107D21.025102t(S 35)7.1391013tD3，（1）tSD式中：为温度（）；为盐度；为深度（m）。本系统采用垂直梯度法对声速剖面进行特征诊断。512海洋科学进展41 卷水声特征诊断流程图如图2 所示。本系统利用声速垂直梯度确认声速剖面所有极小值点和极大值点的深度（图 3）。从海面开始，若初始梯度为负时，则判定该网格点不存在表面声道；若初始梯度为正时，查找确认自海面起第 1 个声速极大值点的深度，判定该深度为表面声道深度（SLD），计算表面声道截止频率（SFD）；若正梯度从海面一直延伸到海底，则确认该网格点满足半声道条件（HAF）。针对声速垂直梯度法确认的声速剖面极小值点深度（从负梯度转正梯度过零点的深度），查找确认声道轴深度。本系统以 350 m 深度作为浅海声道轴深度（SSX）和深海声道轴深度（DSC）分界线。2.1表面声道及表面声道截止频率海-气热通量变化过程和风浪搅拌作用使海洋近表层产生一定厚度的等温层。层内温度均匀，压力随深度增加，声速呈正梯度分布。近表层声源发出的声线在该层中传播时几乎被完全限制在表面层中传播，该层被称为表面声道。表面声道是海洋近表层较有效的声信道，能够实现声能的远距离传播。声能在表面声道中传播存在低频截止现象，只有高于某频率的声能才能传播。表面声道的截止频率（f0）为f0=9cs8H8aH，（2）csaH式中：为海面处的声速；为相对声速梯度；为表面声道深度。表面声道深度和截止频率是表面声道的重要参数。表面声道越深，捕获的声能越强，声传播距 05001 0001 5002 0002 5003 0003 5004 000深度/m1 480 1 4851 4951 5051 5151 4901 5001 5101 520 1 525声速/(ms1)(a)声速声速极大值点声速极小值点05001 0001 5002 0002 5003 0003 5004 000深度/m0.0100.00500.0050.0100.0150.020声速梯度/s1(b)声速梯度梯度从正变负过零点梯度从负变正过零点图 3 声速及声速梯度Fig.3 Sound velocity and its corresponding gradient 水深及对应声速参数判断表面声速梯度诊断声速极大值点深度并计算表面声道截止频率，判断是否为半声道诊断声速极小值点深度根据极小值点深度，确定深浅海声道轴输出水声环境诊断参数正梯度负梯度图 2 水声环境特征诊断流程图Fig.2 Flow chart of the diagnosis of underwater acousticenvironment characteristic 3 期杨春梅，等：基于 FIO-COM 的海洋声学预报系统的构建与应用513Haf0H3/2离越远。本文构建的海洋声学预报系统（FIO-GOAFS）能够给出全球海域的表面声道深度和表面声道截止频率。本系统中，表面声道深度定义为近表层声速极大值所在的深度；之后，由声速数据计算层中的相对声速梯度，并根据式（2）进一步计算截止频率。接下来利用系统给出了 2 月和 8 月的预报结果（图 4 和图 5）。图 4 给出了全球海域的表面声道深度（2 月和 8 月），水平网格为 0.10.1。从预报结果可以看出，低纬度海域常年表面声道很弱，厚度约几十米，赤道附近基本不存在表面声道。中、高纬度海域在夏季表面声道很弱，即使存在，表面声道深度也只有几十米。在冬季，混合层厚度加深，低纬度海域表面声道厚度可达 150200 m；中纬度海域表面声道深度可伸展至大洋主温跃层；高纬度大部分海域存在表面声道，表面声道深度在 100 300 m。表面声道截止频率与表面声道深度成反比，表面声道深度越浅，截止频率越高；表面声道深度越深，截止频率越低。冬季的中、高纬度海域表面声道截止频率很低，约几十赫兹，夏季海域几十米表面声道厚度的表面声道截止频率在 6001 000 Hz（图 5）。由于表面声道位于海洋表层，受到太阳辐射和风雨搅拌作用，因此，表面声道深度和截止频率具有明显的空间分布特征和季节性变化。(a)2 月5004003002001000深度/m8040NS04080180120WE600601201808040NS04080180120WE60060120180(b)8 月5004003002001000深度/m注：白色区域代表不存在表面声道的区域。图 4 全球海域表面声道深度预报结果Fig.4 Forecasting results of the sonic layer depth in the global ocean (a)2 月5004003002001000截止频率/Hz8040NS04080180120WE600601201808040NS04080180120WE60060120180(b)8 月5004003002001000截止频率/Hz注：白色区域代表不存在表面声道截止频率的区域。图 5 全球海域表面声道截止频率预报结果Fig.5 Forecasting results of the surface duct cutoff frequency in the global ocean 2.2半声道条件如果声速正梯度从海面一直延伸到海底，则称为半声道。声速随深度加深逐渐增大，最大声速位于海底。近表层声源发出的声线在传播过程中不断向上折射，在海面向下反射，然后再次向上折射，循环往复，向前传播。半声道效果相当于加强版的表面声道，声能量在半声道条件下能传播更远的距离。图 6 为 FIO-GOAFS 系统给出的全球海域 2 月和 8 月的半声道分布状况，由图 6 可见，半声道主要分布于冬季的极地区域以及部分浅海区域。514海洋科学进展41 卷(a)2 月8040NS04080180120WE600601201808040NS04080180120WE60060120180(b)8 月注：红色区域代表半声道区域。图 6 全球海域半声道分布预报结果Fig.6 Forecasting results of half channel conditions in the global ocean 2.3声道轴声速极小值所在的深度称声道轴。声道轴分为深海声道轴和浅海声道轴两类：深海声道轴是深海的一个稳定的特征，位于主跃层之下；浅海声道轴发生在主跃层，通常跟锋面和涡旋有关。2.3.1深海声道轴深海声道轴深海声道轴位于主跃层之下的最小声速深度。声源位于深海声道轴附近时，其在一定角度范围内出射的声线被限制于深海声道内传播。这部分声线不受海面散射和海底反射的影响，能够传播很远的距离。深海声道又被称为 SOFAR（Sound Fixing and Ranging Channel,SOFAR）声道，利用深海声道可以有效地定位和测距。图 7 为预报系统给出的 2 月和 8 月的深海声道轴深度预报结果。由图 7可见，2 月和 8 月的深海声道轴深度基本一致，说明深海声道轴受季节变化的影响较小，深海声道效应比较稳定，但是深海声道轴深度具有明显的空间分布特征。在纬向上，赤道海域主温跃层强而薄，深度约 300 m，因而赤道海域深海声道轴深度较浅，大约在 800 m；副热带海域主温跃层深度加深，厚度加大，深海声道轴变深，在 1 0001 200 m；高纬度海域主温跃层强度增大，厚度减小，水层变浅，声道轴深度在 350500 m；极地海域不出现深海声道轴。声道轴深度的分布随纬度的变化大体呈“W”形状分布。但是在一些高温、高盐附近的特殊海域，声道轴深度的分布不符合这个规律。例如，在北大西洋东部，由于高温高盐的地中海水溢出直布罗陀海峡下沉，该海域在 1 200 m的水层上出现了大片高温区；同样，在印度洋北部，由于红海和波斯湾的高温、高盐水下沉，使得印度洋北部 1 000 m 水层出现相应的高温区。1 000 m 水层高温区的出现使得北大西洋东部和印度洋北部原有的深海声道变得不明显，声速梯度变得缓和，最小声速值深度出现在 1 5002 000 m。(a)2 月2 0001 5001 0000深度/m8040NS04080180120WE600601201808040NS04080180120WE60060120180(b)8 月2 0001 5001 0000深度/m注：白色区域代表不存在深海声道轴的区域。图 7 全球海域深海声道轴深度预报结果Fig.7 Forecasting results of the depth of the deep sound channel axis in the global ocean 3 期杨春梅，等：基于 FIO-COM 的海洋声学预报系统的构建与应用5152.3.2浅海声道轴浅海声道轴浅海声道轴出现在混合层下界，主温跃层的上层区域，通常与锋面和涡旋有关。地中海和红海声道轴深度常年位于 350 m 以浅，本系统定义其为浅海声道轴，2 个海域不存在深海声道轴。图 8为预报系统给出的 2 月和 8 月的浅海声道轴深度预报结果。由图 8 可见，2 月和 8 月，日本海、阿拉斯加湾的浅海声道轴也一直存在，位于 200300 m。日本海从 200300 m 深度延伸至海底，声速线性增大，声速最小值出现在 200300 m。8 月，浅海声道轴出现在温带太平洋海域，声道轴深度位于 100300 m。某些海域常年存在 2 个声道轴。例如，直布罗陀海峡以西的北大西洋海域，2 月和 8 月均存在 2 个声道轴：浅海声道轴深度约 300 m，深海声道轴深度约 1 800 m，此声速分布的形成是由于地中海高盐暖水团入侵至大西洋1 200 m 深水层的缘故。印度洋 40S 部分海域存在双声道轴，浅海声道轴深度位于 120130 m，深海声道轴深度位于1 362 m 左右。(a)2 月300250200150100500深度/m8040NS04080180120WE600601201808040NS04080180120WE60060120180(b)8 月300250200150100500深度/m注：白色区域代表不存在浅海声道轴的区域。图 8 全球海域浅海声道轴深度预报结果Fig.8 Forecasting results of the depth of the shallow sound channel axis in the global ocean 3声场及相关预报结果海洋声学预报系统除具有上述水声环境特征诊断功能，还能提供声场及相关结果预报。本系统结合了几种模型用于声场计算；之后，在声场计算的基础上，设定声呐优质因数，对声场结果进一步处理，使得系统具备了声呐作用距离预报的功能。该系统最终可以提供 3 种形式的计算结果：距离-深度垂直断面的声场分布；距离-方位角水平平面的声场分布；全球或区域任意深度的声呐作用距离估计结果。相较于垂直剖面和水平剖面的声场预报，全球或区域的声呐作用距离预报的计算量较为庞大，计算过程中使用高性能计算机以提高计算效率。在本节中介绍针对全球海域的声呐作用距离预报。本次预报中设置声呐作用距离水平分辨率为 1。将海域中的计算任务按照海域面积均衡分配给每个进程，并且每个进程读取自己计算范围内的地形、温度、盐度、海面粗糙度等参数。计算前判定声源位置是否位于海洋区域，声源位于陆地直接跳过。对于单个网格点对应的声呐作用距离的求解：首先，假定目标声源在网格点的存在；其次，围绕目标声源计算各个方位角不同深度的水平距离传播损失，将计算结果与声呐优质因数进行比较；最后，对所有方位角的作用距离进行求和平均，得到该网格点的声呐作用距离。声呐系统的作用距离与工作频率息息相关。如果选定声呐优质因数，单从传播损失来考虑，相同海洋环境下，声信号频率越高，传播过程中衰减越大，传播损失也就越大。本文中的计算，只针对中心频率 100 Hz、1/3 倍频程带宽的目标声源；计算时设定声呐优质因数为 80 dB，对于被动声呐516海洋科学进展41 卷来说，该优质因数规定了最大允许传播损失。图 9 为本系统给出的 2021 年 8 月 5 日全球海域的声呐作用距离预报图。由图 9a 可见，当声呐阵列布设于深度 50 m、目标深度为 50 m 时，对于南北极区域，该深度正好处于表面声道深度，声呐在表面声道的作用距离很远，可达 100 km；收、发均在50 m 深度时浅海的声呐作用距离大于深海的声呐作用距离，这是因为声线在浅海中柱面传播，能量分布比较均匀，而在深海中，声线在近场区球面传播，声能衰减很快，声源近距离区域存在声影区。还可看出，8 月南半球（冬季）的声呐作用距离优于北半球（夏季）的声呐作用距离，这是因为在冬季，混合层的存在使得表面声呐作用距离变远。由图 9b 可见，当声呐阵列位于海底、目标深度为 50 m 时，深海中（水深大于临界深度）目标和声呐阵列之间存在着一种重要的传播信道可靠声路径，该路径传播距离远（海深的 57 倍）、传播损失低。因此，收、发分别位于海底、海面时深海区域（4 0005 000 m）声呐的作用距离要远大于表面声呐的作用距离；而且深海区域（4 0005 000 m）的声呐作用距离大于次深海海域（2 0004 000 m）的声呐作用距离。浅海海域的声呐作用距离与海底底质参数密切相关，硬质海底对声能吸收小，声呐作用距离较远；软质海底容易被声线穿透，声呐作用距离近。(a)声呐阵列布设在 50 m 深度302520151050作用距离/km8040NS04080180120WE600601201808040NS04080180120WE60060120180(b)声呐阵列布设在海底302520151050作用距离/km注：目标深度均为 50 m，目标中心频率均为 100 Hz，预报结果水平分辨率为 1，每个网格点计算方位角数为 8 个。图 9 全球海域声呐作用距离预报结果Fig.9 Forecasting results of sonar action rang in the global ocean 4结语基于自然资源部第一海洋研究所自主研发的全球 0.1分辨率海浪-潮流-环流耦合模式(FIO-COM)，本文提出了一种适用于全球海域的水下声学预报系统 FIO-GOAFS。该系统结合了一种水声环境特征诊断方法及几种水下声传播模型，能够提供全球海域的水声环境特征诊断结果及水下声场及相关结果预报。本文给出了 2021 年 2 月和 8 月的水声环境特征诊断结果，包括表面声道深度及表面声道截止频率、半声道条件、深海声道轴深度、浅海声道轴深度。诊断结果表明：表面声道深度和截止频率、半声道均具有明显的空间分布特征和季节性变化；半声道主要分布在冬季极地区域及部分浅海区域；深海声道轴深度常年比较稳定，不受季节变化的影响；浅海声道轴主要跟锋面和涡旋有关，北大西洋东部和印度洋 40S 海域部分存在双声道轴。本文最后给出了 2021 年 8 月声呐阵列分别布放在水深 50 m 深度和海底时全球海域的声呐作用距离预报结果。声呐阵列布放在 50 m 深度时，声呐作用距离主要与表面声道的分布有关；声呐阵列布放在海底时，声呐作用距离主要与可靠声路径有关。需要指出，声呐系统的作用距离与工作频率密切相关，本文所示声呐作用距离只针对声呐优质因数 80 dB、目标声源中心频率 100 Hz 的情况。在同等的声呐优质因数、只考虑传播损失的前提下，声源频率越高，声呐作用距离越近；声源频率越低，声呐作用距离越远。3 期杨春梅，等：基于 FIO-COM 的海洋声学预报系统的构建与应用517参考文献(References)：列布列霍夫斯基赫,扬雷桑诺夫.海洋声学基础M.北京:海洋出版社,1982:14-38.BREKHOVSKIKH L,LYSANOV Y.Fundaments of ocean acousticsM.Beijing:Ocean Press,1982:14-38.1 HALL M V,IRVING M A.Application of adiabatic mode theory to the calculation of horizontal refraction through amesoscale ocean eddyJ.Journal of the Acoustical Society of America,1989,86(4):1465-1477.2 HALL M V.Depth variation of acoustic horizontal refraction through a cold-core ocean eddyJ.Full Field Inver-sion Methods in Ocean and Seismo-Acoustics,1995,12:267-272.3 JIAN Y J,ZHANG J,LIU Q S,et al.Effect of mesoscale eddies on underwater sound propagationJ.Applied Acous-tics,2009,70(3):432-440.4 ROBINSON A R,WALSTAD L J.The Harvard open ocean model:calibration and application to dynamical process,forecasting,and data assimilation studiesJ.Applied Numerical Mathematics,1987,3(1):89-131.5 WARN-VARNAS A,CHIN-BING S A,KING D B,et al.Effects on acoustics caused by ocean solitons,Part A:OceanographyJ.Nonlinear Analysis:Theory,Methods&Applications,2009,71(12):e1807-e1817.6 CHIN-BING S A,WARN-VARNAS A,KING D B,et al.Effects on acoustics caused by ocean solitons,Part B:AcousticsJ.Nonlinear Analysis:Theory,Methods&Applications,2009,71(12):e2194-e2204.7 CHIU C S,EHRET L L.Three-dimensional acoustic mode propagation in the Gulf StreamJ.Journal of the Acousti-cal Society of America,1988,84(Supp.l):S92.8 ROBINSON A R,LEE D.Oceanography and acoustics:prediction and propagationM.New York:American Insti-tute of Physics,1994.9 WARN-VARNAS A C,CHIN-BING S A,KING D B,et al.Ocean-acoustic solitary wave studies and predictionsJ.Surveys in Geophysics,2003,24(1):39-79.10 MELLBERG L E,ROBINSON A R,BOTSEAS G.Modeled time variability of acoustic propagation through a GulfStream meander and eddiesJ.Journal of the Acoustical Society of America,1990,87(3):1044-1054.11 HEANEY K D,CAMPBELL R L.Three-dimensional parabolic equation modeling of mesoscale eddy deflectionJ.Journal of the Acoustical Society of America,2016,139(2):918-926.12 笪良龙,崔宝龙,过武宏,等.海洋-声学耦合过程不确定性快速预报J.声学学报,2017,42(5):535-542.DA L L,CUIB L,GUO W H,et al.The rapid uncertainty prediction of the ocean-acoustic coupled modelJ.Acta Acustica,2017,42(5):535-542.13 肖瑶.中尺度涡旋下的声传播效应及其反演研究D.北京:中国科学院大学,2019.XIAO Y.Effects of mesoscale ed-dies on deep-water sound propagation and tomographyD.Beijing:University of Chinese Academy of Sciences,2019.14518海洋科学进展41 卷Construction and Application of Ocean Acoustic ForecastingSystem Based on FIO-COMYANG Chunmei1,2,3,4，LIU Zongwei1,2,3,4，JIANG Ying1,2,3,4，L Liangang1,2,3,4，XIAO Bin1,2,3,4（1.First Institute of Oceanography,MNR,Qingdao 266061,China；2.Key Laboratory of Marine Science and Numerical Modeling,MNR,Qingdao 266061,China；3.Shandong Key Laboratory of Marine Science and Numerical Modeling,Qingdao 266061,China；4.Laboratory for Regional Oceanography and Numerical Modeling,Laoshan Laboratory,Qingdao 266237,China）Abstract:A global ocean acoustic forecasting system(FIO-GOAFS)for high performance computer is designed and established bycombining acoustic propagation models with the surface wave-tide-circulation coupled ocean model with 0.1 resolution developedby the First Institute of Oceanography,MNR.Sound velocity profiles are calculated with temperature,salinity,and pressure for thediagnosis of underwater acoustic environment characteristic,then the ocean model and underwater acoustic propagation model arecoupled.Combined with geo-acoustic model(seafloor topography and substrate parameters),the diagnosis of underwater acousticenvironment characteristics and the prediction of underwater acoustic field are realized.Ocean model simulated the sound velocityand wave height parameters,and geo-acoustic model provides the parameters including seafloor topography,sound velocity,density,and attenuation of substrate.By calling the ocean-acoustic coupled module,the seafloor topography and ocean environmentalparameters as well as sound propagation path are thus ext