王火平，鞠霞，孙佳，于龙

（1.中海石油(中国)有限公司深圳分公司 深水工程建设中心,广东 深圳 518067；2.自然资源部第一海洋研究所,山东 青岛 266061；3.自然资源部 海洋环境科学与数值模拟重点实验室,山东 青岛 266061；4.山东省海洋环境科学与数值模拟重点实验室,山东 青岛 266061；5.青岛海洋科学与技术试点国家实验室 区域海洋动力学与数值模拟功能实验室,山东 青岛 266237）

南海是西北太平洋最大的边缘海，是内孤立波最频发的海域之一[1]。内孤立波是海洋内波的一种，是海洋密度跃层内产生的一种非线性内波，通常是由强流通过陡峭海底地形激发产生，在传播过程中波型近似保持不变[2]，振幅较大[3-4]。内孤立波是一种常见的海洋现象，在边缘海和陆坡海域出现得尤为频繁[1,5-6]。内孤立波具有大振幅和突发性强流等特征，其对海洋工程施工和水下航行安全等都具有重要影响[7-10]。

南海北部陆坡深水区油气资源潜力大，勘探前景好[11]，是我国深水油气开发的重点发展区域之一。由于南海台风于每年6 月份以后显著增多[12]，因此南海北部陆坡海区油气开采平台建设一般集中于春季进行。该海区频发的大振幅内孤立波，主要源于吕宋海峡，由潮地相互作用激发的内潮在西向传播过程中受非线性、非静力和旋转共同调制下激发形成。其中吕宋海峡处M2内潮影响更强，全日内潮发挥次要作用[13-16]。受复杂地形和水文环境影响，该海区内孤立波模态信号复杂、周期性显著[17-19]。南海北部陆坡海域频发的内孤立波过程严重威胁油气开采平台建设安全，是平台建设过程中最关注的海洋现象[20-21]。如2011 年11 月，一次强烈的内孤立波过程导致处于施工过程中的流花油田半潜式平台在短时间内产生了巨大的漂移，严重影响了工程施工安全[21]。因此，对南海北部陆坡区内孤立波特征进行分析具有重要的意义。

本文通过分析南海北部陆坡海域现场观测数据，对海域内春末第一模态内孤立波重现期以及流向和流速等动力结构特征进行统计分析，分析结果可为该海域工程水下施工提供参考和有力支撑。

1 数据

图1a 为南海北部陆坡海域布放的针对内孤立波过程开展观测所布放的一套浮标周边的水深地形特征，该浮标位于（116°32′32″E，21°26′54″N），水深305 m。为准确捕捉内波信号，浮标于距海表面0.8 m 处设置1 台SeaBird CTD（型号为SEB 37 SM）进行温盐连续观测，距海表面1.5 m 处设置1 台TRDI 75K ADCP（型号为WH LR 75）向下进行流速剖面观测，浮标设计图如图1b 所示。ADCP 流速剖面垂向分辨率为8 m，第一层流速数据位于水深18.3 m 处，共观测33 层，最底层位于274.3 m 处。由于内波过程持续时间短，为准确捕捉内波信号，浮标温度、盐度和流速观测的时间分辨率均为3 min。浮标观测时间由2021 年5 月5 日持续至6 月3 日，共30 d。图1 所用的地形数据为SRTM（Shuttle Radar Topography Mission，数据下载网址为https://portal.opentopography.org/datasetMetadata?otCollectionID=OT.122019.4326.1)[22]，由美国太空总署（National Aeronautics and Space Administration）和国防部国家测绘局（National Imagery and Mapping Agency）联合开发。本文采用的地形数据空间分辨率为15″×15″，包含了陆地高程和海洋深度数据。

图1 浮标位置及其周边地形和浮标结构设计Fig.1 Buoy location and its surrounding topography,schematic diagram of the buoy design drawing

图2 为观测时间内所获取的温度、盐度时间序列，以及ADCP 的仪器姿态和水平流速的剖面序列。由图2a 可知，整个观测期内，仪器纵摇和横摇基本在10°以内，数据良好率在95% 以上，ADCP 观测数据质量可靠。温度序列存在明显的日周期信号，海表温度在5 月15 日前逐渐升温，期间盐度无明显的线性变化趋势，5 月19 日至27 日期间温度升高，盐度同步降低(图2b)。图2c 和2d 分别为观测期内的纬向和经向水平流，存在较明显的日变化。由于ADCP 固定悬挂于水面浮体以下，随波浪起伏，垂向流观测误差偏大，因此本文只利用水平流观测数据进行内波信号提取。

图2 观测期ADCP 仪器姿态和温度、盐度及水平流速剖面序列分布Fig.2 ADCP attitude and the observed temperature,salinity,and meridional and zonal current velocities

2 内孤立波流场结构及判定

第一模态内孤立波根据流场结构可分为2 层，上层水体运动与传播方向一致，下层水体背向传播方向运动。以2021 年5 月24 日12:00 至16:00 的水平流观测结果为例（图3a 和3b），上层70 m 以浅为西北向流，100 m以深为东向流，经向水体运动与上层流向相反，70～100 m 深度范围存在一股较弱的西向水体流动。13:30左右，浮标位置潮流转向，随后伴随产生较为明显内孤立波波列信号，强度逐渐减弱，重现期约30 min。从海表温度、盐度时间序列数据可以发现，内孤立波过程期间温度和盐度曲线均有明显的下凹，结合流场结构可以判断为第一模态内孤立波信号。波列中内孤立波强度随时间降低，海表温度和盐度的波动信号逐渐消失（图3c）。对前4 个信号较强的内孤立波的流场剖面进行集合平均，可发现内波流速剖面为2 层结构，上层以西北向流为主，流向随深度逆时针变化，下层以NE—SE 向流为主（图3d）。以此结构为依据，本文利用2021 年5 月5 日至6 月3 日共30 d 的流场结合温盐场数据进行人工识别判定，将内孤立波波列中的子波作为1 次内孤立波过程，共选取179 次内孤立波过程，并对其特征进行统计分析。

图3 2021-05-24 内孤立波波列水平流速剖面序列分布Fig.3 Vertical structure of meridional and zonal current velocities of the internal solitary waves during 24 May,2021

3 内孤立波重现期

内孤立波在南海北部陆坡区域多以波列形式出现，重现期较短，强度较大（图3）。南海北部主要为不规则半日潮类型，内潮对内孤立波重现期有显著影响。对浮标捕捉到的内孤立波信号进行分析，发现南海北部的内孤立波主要存在3 个重现期信号（图4a）。其中，重现期为0～1 h 的内孤立波占总数的58.4 %，这说明波列形式是南海第一模态内孤立波出现的主要形式。由于南海背景环流对内孤立波流场的影响以及人工判别方式可能存在的误差，部分信号较弱的内孤立波过程可能未被记录，从而导致40～60 min 重现期占比较大。综合考虑上述因素的影响，波列中的内孤立波重现期应在20～30 min（图4b）。除波列形式的重现期信号外，该区域的内孤立波仍存在半日（11～12 h）和全日（24 h）两个重现期信号（图4c）。这2 个信号主要是由半日内潮和全日内潮在吕宋海峡处发生潮地相互作用，以及内潮非线性变陡机制激发效应所致[11-14]，显示了南海北部陆坡海区第一模态内孤立波的主要变化周期。激发生成的半日和全日周期的内孤立波由巴士海峡生成后以波列形式向西传播，影响南海北部陆坡区域海洋流场特征。

图4 内波重现期概率密度分布Fig.4 Histogram of the recurrence interval of internal solitary waves

4 内孤立波传播方向及流速特征

图5 为陆坡区域浮标位置处观测到的179 个内孤立波流向流速玫瑰图及表层（ADCP 有效观测结果第一层）流流速、最大流流速以及最大流速所在深度的概率密度分布图。由图5 可知，该区域内波传播方向角度多在281°～303°，以WNW 传播方向为主，NW 方向次之（图5a）；最大流速集中在上层区域，与内波传播方向一致，少部分发生于背向内孤立波传播方向的深层区域（图5b）。表层流流速多分布于20～80 cm/s范围内，占总数的81.5 % 以上（图5c）；最大流速多分布于40～120 cm/s，占总数的88.4 % 以上（图5d），较表层流速稍强，集中分布于100 m 以浅，以西北向流速结构为主（图5e）。此外，近底层存在一处最大流速汇集区（图5e），表明内孤立波可能引发底流的突变，并对海底附近的施工作业产生重要影响。

图5 内孤立波流向流速玫瑰图及表层流流速、最大流流速以及最大流速所在深度的概率密度分布Fig.5 Speed-direction rose maps of the internal solitary waves and histogram of surface current speed,maximum current speed,and depth of the maximum current

图6 为浮标位置处根据现场观测数据统计得到的内波传播方向示意图。浮标处水深近300 m，水深变化梯度大，等深线梯度由西逐渐向北偏转，浮标位置处地形梯度为北偏西方向。受地形变化影响，吕宋海峡传至此处的内孤立波传播方向发生改变，由西逐渐向北偏转，最终形成西偏北方向的传播特征。

图6 浮标周边海区地形、内孤立波波峰线及传播方向示意图Fig.6 Topography and the diagram of peak line and propagation direction of the internal solitary waves

由于不同内孤立波之间的结构差异性，因此本文对观测期内179 个内孤立波的流速进行集合分析。图7 为179 个内波期间的流速分量大小和流矢量剖面结构。由图7 可见，第一模态内孤立波普遍的双层结构特征，南海陆坡区域处上层水体为西北向流动，下层为东南向流动。受潮汐运动影响，流场存在明显半月周期信号。纬向流在5 月16 日和6 月1 日左右达到极大值，分别对应小潮和大潮。内孤立波流场由海表向下沿逆时针转动，流速极小值（流场转向处）主要集中在100～150 m 深度层（图7c）。

图7 内孤立波流速结构Fig.7 Vertical structure of the current associated with the internal solitary waves

5 结论

本文通过分析处理南海北部陆坡海区2021 年5 月5 日至6 月3 日期间1 套内孤立波浮标的流速和温盐观测数据，对该海区内孤立波过程进行了识别判定和结构分析。通过对该浮标捕获的179 次第一模态内孤立波信号的重现期、传播方向和流速结构等特征进行集合分析发现，南海北部陆坡海区春季第一模态内孤立波垂向流场为双层结构，上层为西偏北向流动，下层为东南向流动。受内潮的影响，南海北部陆坡区域内孤立波存在半日和全日两个重现周期，主要以20～30 min 间隔的波列形式向西偏北方向传播，传播方向受陆坡海区地形梯度制约。内孤立波期间，最大流速多发生于水深100 m 以浅的上表层，流速集中在60～120 cm/s，底层流增强。流向随深度逆时针变化，并于100～150 m 水深处发生层结流场转向。这些结果表明，南海北部陆坡海区内孤立波频发、规律较为明显。

内孤立波伴随的水下强流主要出现在上表层，容易引发作业船只定位漂移、对油气平台水下结构和施工人员形成猛烈冲击，造成重大安全隐患。因此，该海区油气平台建设施工时应主动进行内孤立波实时监测预警，从而有效保障生产安全。本文研究成果也提升了对该海区内孤立波时空变化规律的认识水平，丰富了内孤立波研究成果。