3. 海洋环境数据获取与预处理

水下滑翔机不靠螺旋桨推进，而是通过改变自身净浮力做锯齿状滑翔，因此它的航迹对海流极为敏感。规划一条从A点到B点的航线，如果只按直线距离计算，滑翔机很可能被横向流推离预定断面，甚至耗尽电池仍无法到达目标点。所以，航线规划的第一步不是画线，而是把“海洋环境”变成一张可计算、可查询的场。本节围绕数据源、剖面数据、流场格式与分辨率、插值网格化，以及最终如何构建面向航线规划的环境场展开。

一、常用海洋环境数据源及其角色

面向滑翔机任务的环境数据大致分四类，它们的时间、空间特性差异很大，决定了各自在规划中的用途。

• CTD（温盐深）：提供温度、盐度随深度的剖面，是计算海水密度、声速以及滑翔机净浮力与俯仰姿态的基础。滑翔机的浮力驱动效率直接取决于密度剖面，CTD数据不可替代。
• ADCP（声学多普勒流速剖面仪）：直接测量分层流速流向，是局部流场最可信的来源。船载ADCP走航观测可给出沿航线的剖面，锚系ADCP则提供时间序列。
• 卫星高度计：通过海面高度异常反演地转流，覆盖范围广、时间连续，适合获取大尺度背景流场，但只能给出表层附近的地转分量，且近岸精度下降。
• 再分析数据：如HYCOM、CMEMS等产品，融合卫星、浮标和模式结果，提供规则网格上的三维温盐流场，时间跨度长、时空连续，是航线规划最常用的背景场，但分辨率通常较粗（水平约1/12°量级，垂向数十层）。

二、温盐深剖面与流场数据的格式和分辨率

CTD剖面通常以“站号—压力（深度）—温度—盐度”的列式记录保存，NetCDF是主流格式，一个文件可含多个剖面。处理时先做质量控制：剔除压力倒转、盐度尖峰和温盐不匹配点，再按1 m或2 m间隔垂向平均，形成等间距剖面。滑翔机自身也携带CTD，其观测剖面可与历史数据比对，用于判断背景场是否仍然有效。

流场数据的组织方式更需注意。再分析产品一般给出u（东向）、v（北向）分量，网格为规则经纬度格点，垂向为sigma层或固定深度层。使用前必须明确两点：一是流向约定，u正方向为东、v正方向为北，流向角从正北顺时针度量；二是深度基准，是相对海面还是相对海底。分辨率方面，水平1/12°约合8～9 km，对滑翔机（时速约0.5～1 kn）而言，一天行程约20～40 km，跨越数个网格，因此背景场足以描述大尺度输运，但无法刻画中尺度涡的精细结构，这正是需要实测校正的原因。

三、插值与网格化方法

原始观测是散点，规划算法需要规则网格上的场，因此要做插值。常用方法有三类：

• 反距离加权（IDW）：实现简单，适合观测点较密的局部修正，但会产生“牛眼”伪影，不宜外推。
• 最优插值（OI）：引入背景场误差协方差和观测误差协方差，能在稀疏观测下给出较合理的分析场，是海洋数据同化的经典方法。
• 克里金（Kriging）：基于变差函数，可给出插值误差估计，适合评估哪些区域数据可信度低。

实际操作中，更常见的做法是“背景场+增量校正”：以再分析网格为背景，把ADCP实测流速减去背景值得到残差，对残差做IDW或OI插值后加回背景场。这样既保留了背景场的大尺度结构，又吸收了实测信息。垂向上，若实测剖面层数与背景场不一致，需先线性插值到统一深度层，再做水平插值。

四、面向航线规划的环境场构建

规划所需的环境场不是单一变量，而是一组按深度分层的二维网格：每一层包含u、v分量，另附密度剖面用于估算滑翔机的垂向运动与能耗。构建流程可概括为：确定任务时空范围→选取再分析产品并裁剪→统一深度层与坐标约定→融合实测数据→生成规则网格→质量检查。

质量检查的关键是看流场是否“物理合理”：检查是否有孤立异常矢量、流速是否超过区域合理上限、相邻网格梯度是否突变。对滑翔机而言，还要特别关注目标航线上的横向流分量——它直接决定航迹偏移量。一个实用的判断是：若某段航线的横向流长期接近滑翔机的水平滑翔速度，则该段航线需要预留修正角或改线。

上图概括了从四类数据源到预处理、再到融合生成规划环境场的主线，下方示意了“背景场加实测残差”的融合思路：背景场提供大尺度结构，实测残差修正局部细节，二者叠加得到既连续又贴近观测的流场。

五、一个可操作的构建示例

假设任务区为某近海断面，规划深度取0～200 m，时间窗口为未来5天。可按以下参数组织：

需要强调，表中数值为教学示例，实际任务应根据区域动力特征和数据可得性调整。构建完成后，环境场应能回答规划算法最关心的三个问题：某点某深度的流速流向是多少、沿航线累积的横向偏移有多大、密度剖面是否支持预期的滑翔姿态。只有这三个问题都能被稳定查询，环境场才算真正可用于航线规划。