使用声学多普勒测速（ADV）的河口盐度反演：方法、敏感性与环境调制

《Journal of Marine Science and Engineering》：Estuarine Salinity Inversion Using Acoustic Doppler Velocimetry (ADV): Methodology, Sensitivity and Environmental Modulations

　　

从声学反向散射信号中获取水柱的水文物理特性对于实现河口环境中连续、非侵入式观测至关重要。本研究利用声学多普勒流速仪（ADV）中的脉冲相干技术，通过整合同步的温度和压力数据集，提出了一种估算河口水体平滑低频实用盐度的方法。实验室标定实验表明，在0至35 PSU的盐度范围内，梅德温（Medwin）公式的反演精度最高，但当盐度低于15 PSU时数据离散度增大。该方法进一步通过现场数据进行验证，其中声强剖面通过多探头空间平均和二次多项式拟合获得。现场应用结果显示，300分钟滑动平均滤波器可从湍流噪声中提取低频盐度入侵趋势，使框架能够追踪潮汐尺度盐度变化，相对于参考仪器，其低通趋势精度为±3.07 PSU，而原始未滤波反演的均方根误差（RMSE）为5.68 PSU。反演性能对环境动力学敏感：在流速为0.10–0.58 m/s、浊度为109.7–208.0 NTU的适中环境窗口内，误差偏差最小；而在较高流速（高达0.83 m/s）和严重浊度（高达278.1 NTU）期间，或流速低于0.10 m/s、声学反向散射降至90 dB以下的平潮期，不确定性增加。这些发现定量界定了声学盐度估算的环境约束条件，为恢复动态河口的低频盐度趋势提供了一种低成本、非侵入式的方法框架。

**论文解读文章：基于声学多普勒测速仪（ADV）的河口盐度反演：方法、敏感性与环境调制**

**研究背景与问题**

盐度是河口环境中控制水体物理化学性质的关键参数，其变化驱动密度梯度，影响环流、湍流混合、沉积物输运及溶解氧通量。传统盐度监测依赖电导率-温度-深度（CTD）仪器和光学反向散射传感器（OBS），但这些方法存在局限性：单点采样难以获取空间剖面，部署需高昂船时和人力成本，且物理探头侵入性会扰动微结构。卫星遥感虽能提供大范围表面盐度，但空间分辨率不足且无法解析底层垂向结构。声学技术（如声学多普勒流速仪ADV）通过反向散射信号可非侵入式监测流速、地形和悬浮沉积物浓度（SSC），但利用声信号估算盐度的研究仍较少。声速随温度、盐度和压力变化，为声学盐度反演提供了理论基础，但现有研究缺乏对不同反演公式的系统性评估，且复杂潮汐环境下的操作边界尚不明确。因此，本研究旨在构建基于ADV的盐度反演方法框架，通过实验室标定和河口现场观测，量化环境因素对反演精度的影响，并确定该方法的应用范围。

**研究内容与结论**

研究人员在《Journal of Marine Science and Engineering》发表论文，提出了一个整合实验室标定和现场观测的声学盐度反演方法。在实验室中，通过控制水槽实验评估了三种经典声速公式（Coppens、Chen–Millero–Li、Medwin）在0–35 PSU盐度梯度下的表现，并确定了Nortek Vectrino ADV的仪器常数。现场观测在珠江口黄茅海进行，持续115小时，利用6 MHz Nortek Vector ADV和OBS-3A传感器对比验证。结果表明，Medwin公式在整体盐度范围内反演精度最高，均方根误差（RMSE）为0.5–1.8 PSU，平均相对误差（MRE）最低（7.38%）。现场数据通过300分钟滑动平均滤波后，低通趋势精度为±3.07 PSU，而原始未滤波反演RMSE为5.68 PSU。反演精度受环境流场和浊度调控：最佳工况出现在流速0.10–0.58 m/s、浊度109.7–208.0 NTU的适中窗口；高流速（>0.58 m/s）和强浊度（>208.0 NTU）时误差增大；平潮期流速<0.10 m/s导致反向散射强度低于90 dB，误差达最大。该研究定量定义了声学盐度反演的环境约束，为动态河口提供了一种低成本的非侵入式低频盐度趋势监测方法。

**主要关键技术方法**

研究人员采用了以下关键技术方法：（1）基于ADV脉冲相干技术的声学反向散射强度剖面获取，通过多探头空间平均和二次多项式拟合确定峰值位置（D_max），进而反演声速。（2）利用Medwin声速公式结合同步温度和压力数据计算实用盐度，通过与OBS-3A参考传感器对比验证。（3）现场数据采集于珠江口黄茅海河口沙坝（113°05.127′E, 21°59.280′N），使用6 MHz Nortek Vector ADV和OBS-3A，以32 Hz采样率、300秒采样间隔、600秒爆发周期记录，持续115小时。实验室标定在80 cm × 60 cm × 60 cm玻璃水槽中进行，使用10 MHz Nortek Vectrino，以99.7% NaCl调节盐度梯度。

**研究结果**

**3.1 仪器常数标定**
通过纯净水标定实验，研究人员确定了Nortek Vectrino的仪器常数K_const为74 m²/s。不同时间窗口（2–30分钟）的统计表明该常数稳定，收敛于该值。

**3.2 实验室反演性能评估**
在7个盐度工况（0–35 PSU）下，三种声速公式均能反演盐度，但Medwin公式表现最佳，在整个盐度范围内RMSE为0.5–1.8 PSU，偏倚对称（-1.2至+1.5 PSU），平均相对误差MRE=7.38%。Coppens公式在>15 PSU时高估，Chen–Millero–Li公式在10–20 PSU时低估。低盐度区（<15 PSU）所有公式RMSE增大，因声速对盐度敏感度降低。

**3.3 现场反演性能评估**
现场数据覆盖115小时，从农历小潮到大潮过渡。对比ADV反演与OBS-3A参考盐度，整体相关性较高（R>0.8），但存在潮相依赖性。涨潮和落潮期间R>0.8，RMSE<5 PSU；平潮期R<0.7，RMSE达6 PSU。原始数据标准差为4.98 PSU，应用300分钟滑动平均滤波后降至3.07 PSU，有效提取了低频潮汐盐度趋势。

**讨论与结论**

**讨论部分总结**
反演精度对峰值位置D_max的定位敏感，Vector传感器由于几何设计具有更低的灵敏度系数（-8.4 PSU/mm），优于Vectrino（-27.2 PSU/mm），更适用于动态环境。温度变化增大灵敏度，压力在浅水区影响可忽略。误差传播模型预测总不确定度4.55 PSU，解释约91%的观察方差。300分钟滑动平均滤波在噪声抑制和信号保留间取得平衡，虽导致峰值振幅衰减约10–12%，但保持相位一致性。传感器空间偏移（水平1 m，垂直0.05 m）对盐度差异贡献可忽略。反演性能受流速和浊度协同调控：最佳工况（流速0.10–0.58 m/s，浊度109.7–208.0 NTU）下RMSE=4.31 PSU；高流速高浊度时RMSE=4.98 PSU；平潮期低流速低浊度时RMSE=5.69 PSU。方法局限性包括：实验室纯NaCl溶液与自然海水离子组成差异导致声速偏差<0.4 m/s，可忽略；需足够悬浮颗粒维持信噪比，不适用于清澈水体；适用于盐度变化大、背景盐度高的环境，不适用于低盐度稳定水体。

**研究结论部分翻译**
通过整合实验室标定与珠江口现场观测，本研究定量评估了从ADV声学反向散射剖面恢复平滑潮汐尺度实用盐度趋势的可行性。所提出的方法并非用于替代传统高精度或高频电导率传感器，而是为非侵入式提取长期水文物理宏观动力学提供了一种补充工具。主要结论如下：（1）在0–35 PSU盐度范围内，三种经典声速公式（Medwin、Coppens、Chen–Millero–Li）均展现出稳健反演能力，其中Medwin参数化一致性最高；但低盐度区（<15 PSU）跟踪不确定性增加。（2）未滤波原始反演RMSE=5.68 PSU，源于湍流扰动和电子噪声；通过空间滤波、二次曲线拟合和300分钟滑动平均滤波，误差降至±3.07 PSU，表明该方法适用于解析多日平滑水文波动和长期盐度入侵动态，而非高频瞬时变化。（3）局部水动力条件界定了反演性能的环境边界：静水压力在浅水区（<50 m）影响可忽略，而流速和悬浮沉积物通量强烈调制声学信噪比（SNR）。最佳跟踪精度出现在适中环境窗口（流速0.10–0.58 m/s，浊度109.7–208.0 NTU），误差在流速>0.58 m/s和浊度>208.0 NTU时增大，在平潮期流速<0.10 m/s时声学反向散射峰值降至90 dB以下，误差达最大。