3. 声速剖面数据获取与预处理

声速剖面建模的第一步不是“算声速”，而是拿到一条可信的声速—深度（或声速—压力）序列。海洋中温、盐、压的观测手段不同，数据格式、采样密度、误差结构和缺测模式也完全不同。若把原始观测直接代入声速经验公式，往往会把仪器启动瞬态、传感器响应滞后或定位误差放大成虚假的声速跃层。因此本节围绕三类主力数据源，讲清“拿到什么格式、怎么判断好坏、怎么剔除异常、怎么插值成可用剖面”。

三类数据源的格式与物理含义

CTD（温盐深仪）是最直接的剖面观测。下放式CTD通常以1～24 Hz采样，输出压力（dbar）、温度（℃）、电导率（S/m）或盐度、溶解氧等通道，常见交换格式为CNV（Sea-Bird）或NetCDF。CNV为文本格式，头部以*开头记录元数据，正文为多列数值，列名在# name行中给出。CTD的优点是垂向分辨率高、温盐同点同步测量；缺点是存在热滞后（温度传感器响应慢于压力）、电导池冲刷滞后，以及下放速度不均导致的非均匀采样。

XBT（抛弃式温深仪）只测温度—深度，不含盐度，常见格式为EDF或厂商自定义文本。其深度由下落时间经验公式换算，存在系统性深度偏差；探头类型（T-4、T-6、T-7、Deep Blue等）决定深度范围和误差量级。XBT最典型的问题是“温度偏移”和深度偏差，且没有盐度，必须借助气候态盐度或邻近CTD/Argo数据补盐后才能计算声速。

Argo浮标提供自动剖面，标准交换格式为NetCDF，核心变量包括PRES、TEMP、PSAL及质量控制标志PRES_QC、TEMP_QC、PSAL_QC。Argo剖面在上升过程中采样，垂向间隔约数米至十余米，且经过延迟模式（D-mode）再处理后可获得更高质量。其优势是全球覆盖、长期连续；劣势是垂向分辨率低于CTD、近表层数据受卫星定位与传感器稳定时间影响。

质量控制：先看标志，再做物理检查

对Argo数据，第一步永远是读QC标志，而不是直接读数值。约定俗成的标志含义是：1为“好”，2为“可疑”，3为“坏”，4为“校正后可用”。工程上通常保留1和4，对2做人工复核，对3直接置为缺测。对CTD的CNV文件，则要检查头部元数据中的# interval、# nvalues、# bad_flag，确认采样间隔与缺测标记，避免把-9.990e-29这类坏值当成真实温度。

物理检查至少包含四类：范围检查（温度约-2～40 ℃、盐度0～42、压力非负且单调增）；梯度检查（相邻层温度跳变超过阈值即判可疑，例如CTD常用0.5～1 ℃/m量级作粗筛）；尖峰检查（某点相对前后邻点的偏离超过3倍局部标准差）；密度倒置检查（由温盐压算出的位势密度应随深度基本不减，明显倒置往往意味着传感器滞后或盐度异常）。

异常值剔除：针对数据源“对症下药”

CTD最需要处理的是热滞后。温度传感器响应时间常数约0.1～0.5 s，而电导池冲刷时间常数更短，两者不同步会在强梯度层产生虚假的盐度尖峰。工程做法是对温度通道做一阶滞后校正：T_corr = T + τ·dT/dt，其中τ由传感器型号与下放速度标定，常见取值在0.1～0.5 s。校正后再算盐度，尖峰通常显著减弱。

XBT的核心问题是深度偏差与无盐度。深度偏差可按探头类型查表修正，例如T-4、T-6、T-7的深度误差量级不同，深水型需用Deep Blue校正。补盐时优先用同区域、同季节的Argo或CTD气候态盐度，若邻近有同步CTD则直接匹配；补盐后必须记录“盐度来源”，否则声速结果不可追溯。

Argo的异常多出现在近表层（传感器稳定时间不足）和深层（压力漂移）。处理原则是：QC=3的层直接剔除；QC=2的层结合邻层梯度判断，若与上下层趋势明显不符则剔除；对整条剖面做密度倒置检查，倒置层用上下有效层线性内插替换，而不是简单删除，以免破坏垂向连续性。

插值：把非均匀采样变成统一网格

三类数据的垂向采样都不均匀：CTD下放速度变化导致压力间隔忽大忽小，XBT按时间换算深度，Argo上升速度也不恒定。声速建模通常需要等压力间隔的剖面，例如每1 dbar或每5 m一层。插值方法的选择直接影响声速跃层的陡缓。

线性插值最简单、无过冲，适合采样较密（如CTD原始1 Hz数据）的情形；单调三次插值（PCHIP）在保持单调性的同时给出更平滑的曲线，适合Argo这类垂向间隔较大的剖面，能避免普通三次样条在跃层处产生的虚假振荡；按压力等间隔重采样则是先确定目标网格（如0～2000 dbar、步长1 dbar），再对每个目标层用相邻有效观测插值。对缺测段，应限制插值跨度，例如连续缺测超过20 dbar时标记为“不可插值”，而不是无限外推。

一个可复现的处理顺序是：读入原始文件 → 解析QC标志 → 范围/梯度/尖峰/密度倒置检查 → 热滞后或深度偏差校正 → 逐层剔除 → 按目标压力网格插值 → 输出带来源标记的剖面。以某次示例CTD剖面为例（以下为教学假设数据）：原始采样约1200层，经QC与尖峰检查剔除约3%的层，热滞后校正后盐度尖峰幅度由约0.15降至0.03，按1 dbar重采样后得到约1500层等间隔剖面，声速跃层位置与原始剖面偏差小于1 dbar。判断处理是否成功的依据不是“曲线好看”，而是跃层位置稳定、密度不再倒置、且被修改层均有记录。

完成上述步骤后，得到的才是一条“可用于声速计算”的剖面。后续无论用Mackenzie还是Chen-Millero公式，输入的温度、盐度、压力都应是经过质控、剔除与插值的一致网格数据，这样声速剖面的跃层结构才可信，建模结果才具备物理意义。