无姿态测量条件下MEMS有效波高估算中的处理诱发不确定性

《Frontiers in Marine Science》：Processing-induced uncertainty in MEMS-derived significant wave height estimates without attitude measurements

　　

摘要专业翻译 引言： 由紧凑型微机电系统（MEMS）导出的有效波高（SWH）日益广泛地用于业务化波浪监测，但对同一惯性记录应用科学上可辩护的不同处理方案，所得估算结果可能有所差异。 方法： 本研究利用来自一个业务化浮标和一种MEMS配置的31天加速度记录，评估方法选择变异性，该记录中姿态观测不可用。对一条共同的倾斜轴加速度记录进行条件化处理，并通过受控的谱域和时域变体进行处理。敏感性评估以指定的参考工作流为基准，并对每个观测时段计算规定谱集合的经验方法间离散度。制造商的时间对齐SWH产品和区域NDBC浮标仅用于情境对比，不作为真值或正式验证。 结果： 针对所研究的记录，规定集合产生的平均处理离散度为0.023 m，中位数为0.018 m，95百分位为0.060 m；相应的95百分位相对离散度为6.42%。在所测试的配置中，时域稳定化修改产生的参考相对差异大于在共同加速度条件化后替换所评估的谱估计器所产生的差异。 讨论： 这些统计量仅量化规定集合的处理方法分量，不构成完整的测量不确定度预算。数值结果特定于所研究的数据集、传感器配置、处理方案和海况。

论文解读

一、研究背景与问题提出

有效波高（SWH）是海洋工程、海上作业、模型评估、卫星定标、灾害评估和波浪数据同化中广泛使用的重要参数。随波浮标通过加速度计、惯性测量单元（IMU）或全球导航卫星系统观测来估算海况，紧凑型微机电系统（MEMS）平台极大地拓展了空间分布式波浪测量的可行性。然而，从加速度恢复位移的过程对重力去除、传感器姿态、低频污染、滤波、去趋势、积分限和所选谱估计器均十分敏感。由于位移功率谱密度通过除以(2πf)⁴获得，微小的低频加速度分量可能被强烈放大。

以往研究已确立紧凑型波浪浮标的可行性，并比较了仪器、加速度恢复方法或谱估计器，但多数比较聚焦于与另一台仪器的吻合度或特定算法的性能，而对同一现场记录在保持观测输入固定条件下，应用一组可辩护的下游处理工作流所产生的逐时段离散度关注不足。本研究明确将所评估的量限定为处理选择分量，而非总测量不确定度。基于GUM的评估通过测量模型传播特征化输入不确定度，蒙特卡洛评估传播指定的输入概率分布，而重复性指特定条件下测量的一致性。

本研究利用一个无姿态观测的31天浮标部署，提出三个问题：处理选择对共同加速度记录导出的SWH影响有多强；所评估的条件化、重构和估计选择中哪些产生最大差异；规定谱集合的逐时段离散度是多少。

二、研究内容与主要结论

研究人员利用2023年10月在墨西哥湾部署的配备Murata SCC1300-D02惯性传感器的Viking波浪浮标采集的倾斜轴加速度测量数据。该传感器包含三轴加速度计和单轴陀螺仪，但分析数据集中仅三个加速度通道可用，滚转、俯仰、偏航和角速率观测均不可用，因此无法进行时变体坐标系到导航坐标系的姿态变换。加速度以4 Hz采样频率、10分钟观测窗口记录，窗口起始间隔15分钟。

研究将每个10分钟加速度窗口固定，同时变化SWH工作流的指定要素。程序包括：共同倾斜轴加速度信号的定义与条件化、谱域或时域重构、测试SWH估计器的应用、以及参考相对敏感性和逐时段集合离散度计算。采用单因子逐次变化设计，V0为参考工作流，V1–V5考察谱估计器或窗函数选择，V6–V7考察时域速度去趋势选择。

主要结论为：在所测试的配置中，时域稳定化修改产生的参考相对差异大于共同条件化后谱估计器替换所产生的差异。对于规定的V0–V3谱集合，平均和中位处理离散度分别为0.023 m和0.018 m，95百分位为0.060 m，相应95百分位相对离散度为6.42%。这些数值仅量化规定集合的处理方法分量，不构成完整测量不确定度预算。

三、关键技术方法

研究人员采用的主要技术方法包括：倾斜轴加速度重构，即在无姿态观测条件下通过去除标称重力获得工作加速度信号；信号条件化，包括带通滤波（0.102–0.695 Hz）、去均值和线性去趋势；谱域重构，分别采用Welch方法和直接快速傅里叶变换（FFT）配合不同窗函数（Hann、Blackman、Chebyshev、Hamming）估计加速度功率谱密度，再通过除以(2πf)⁴获得位移谱；时域重构，通过带通滤波后双重积分恢复位移，并对比速度去趋势启用与禁用两种配置；有效波高估计，谱域采用4√m₀，时域采用最高三分之一波高均值和四倍位移标准差。样本来源为2023年10月墨西哥湾Viking浮标31天部署记录，另有两个补充数据集（2023年8月和2026年3月）来自同一浮标和位置。

四、研究结果

处理诱发不确定性： 谱变体V1–V5与参考工作流V0的对比显示，谱变体紧密跟随V0，95百分位绝对相对差异低于约13%，均方根误差低于0.06 m。谱变体V0–V5相对集合均值的95百分位绝对集合相对差异不超过12.06%，相应绝对差异低于0.1 m。时域配置显示更大的参考相对偏离，V6（启用速度去趋势）的95百分位绝对相对差异为26.72%，而V7（禁用速度去趋势）为7.83%，表明启用速度去趋势增加了参考相对敏感性。规定的V0–V3集合产生平均处理离散度0.023 m、中位数0.018 m、95百分位0.060 m，相对离散度均值、中位数和95百分位分别为2.70%、2.21%和6.42%。

情境对比： 将规定谱集合、参考V0和代表性时域工作流与制造商同平台产品和NDBC 42040站对比。NDBC序列显示大致相似的区域演变，但99.26 km的空间间隔允许实质性空间差异。相对制造商同平台产品，集合、V0和V7的相关系数均大于0.99，平均绝对误差分别为0.082、0.079和0.046 m，95百分位绝对相对差异分别为28.28%、26.18%和17.21%。该一致性表明与制造商产品在所研究记录上的吻合，但不确立绝对精度。

五、讨论总结与研究结论

讨论部分指出，处理影响SWH估算，主要通过数据导出的有效频带、加速度条件化和所选配置。时域稳定化修改产生比谱估计器替换更大的参考相对差异，这与加速度到位移重构的物理机制一致：(2πf)⁴除法强烈放大低频加速度能量，因此条件化和频带限制在波能矩计算之前起作用。经验离散度仅隔离处理方法变异性，与外部验证误差有本质区别。该框架补充了传统验证，将内部处理变异性与外部一致性分离。对于所检查的三条记录，条件化工作流产生稳定的SWH时间序列，谱变体内部差异较小，支持该处理策略在该平台和传感器配置上的跨记录可重复性。局限性包括：数据导出的频带和主要定量分析基于单一浮标和单一MEMS传感器配置的31天记录；长期记录不可用；SCC1300系列已被Murata标记为不推荐用于新设计并计划停产。

研究结论为：对于所研究的31天部署和补充数据集，测试的谱估计器在共同加速度条件化后产生相对较小的差异，而时域稳定化配置产生较大的参考相对偏差。2023年10月数据集的V0–V3集合处理离散度均值、中位数和95百分位分别为0.023、0.018和0.060 m，95百分位相对离散度为6.42%。这些数值以规定集合为条件，不构成总测量不确定度预算。传感器误差、姿态效应、平台响应、环境变异性、时间和外部产品不确定度均未传播。在该数据集所代表的条件下，所选工作流产生与制造商同平台产品一致的稳定SWH时间序列，但因无共址独立参考，该一致性不确立绝对精度。该程序提供了一种可重复的方法来检验方法选择变异性，但数值结果和性能结论仍特定于所研究的部署、传感器配置、处理方案和海况。