未测试搜索对象风致漂移速率估计可接受误差的确定：基于风—流—漂移速率矩阵模拟的敏感性分析

《Journal of Marine Science and Engineering》：Determination of Acceptable Error for Leeway Rate Estimates of Untested Search Objects

　　

物体的风致漂移（leeway）直接影响漂移预测（drift prediction）与搜索救援（search and rescue，SAR）行动。计算物体leeway的标准方法是通过直接测量环境力（environmental forces）实现。然而，直接测量成本高昂且需要大规模野外测试，因此许多物体并不具备准确的leeway值。当出现涉及leeway未知物体的案例时，其漂移必须基于形状与尺寸相似的已知物体进行估计。leeway的不确定性——尤其是在使用替代物体（proxy object）时——可能导致施加错误的风强迫（wind forcing）量，并导致目标物体预测位置的差异。本研究针对leeway rate对漂移预测范围与位置的影响开展了敏感性分析（sensitivity analysis）。研究在风、流与leeway rate的矩阵化取值组合下完成了400次漂移模拟，以确定leeway rate的可接受误差限值（acceptable margin of error）。可接受误差定义为：在36小时（h）漂移过程中，真实物体与替代物体的粒子分布（particle distribution）之间至少保持50%重叠的任意leeway rate差异。敏感性分析的结果与真实环境数据的测试表明，leeway rate可存在相当于风速±1%的误差。

**研究背景**

漂流物体既是研究海洋环流的传统工具，也是搜索与救援行动的关键对象。自Theophrastus利用漂流瓶验证水体运动、Benjamin Franklin依据马尾藻（Sargassum）漂移绘制墨西哥湾流以来，漂流物体长期被用于测量洋流；Pingree率先开创物体leeway的直接测量；Allen与Plourde提出AP98方法，将leeway定义为风与波浪强迫引起的物体相对于海流的运动，其量值为leeway rate（LW，以风速百分比表示），方向由发散角（divergence angle）描述，并可分解为逆风漂移分量（downwind leeway，DWL）与侧风漂移分量（crosswind leeway，CWL）；Breivik等随后确立了leeway直接测量与间接测量的标准化方法。然而，直接测量需要走航式水流速剖面仪、气象站与定位设备，成本高且依赖大规模野外测试，因此已知leeway系数的物体十分有限；当案例涉及未知leeway的物体时，只能以形状与尺寸相似的已知物体作为替代物体（proxy）估计漂移。Cho等、Ni等的研究虽证明漂移预测对leeway高度敏感，但未回答"leeway rate需要多准确，预测搜索区域才不会在范围与位置上出现显著差异"这一作业层面的问题。

**研究内容与结论**

研究人员以一项为期五天的现代野外试验中19英尺panga（长尾小艇）的空载配置为基准搜索对象（leeway rate=5.98%，发散角=18.53°），在美国海岸警卫队使用的Search and Rescue Optimal Planning System（SAROPS，搜救最优规划系统）中开展了400次漂移模拟（风、流与leeway rate的矩阵化组合：风速5–25 kt、流速0.25–1 kt、leeway rate 0.48%–9.98%），系统量化leeway rate变化对预测搜索区域范围、位置与空间重叠的影响，并依据既定作业标准确立了leeway rate的可接受误差限值。

**关键方法**

研究人员运用以下关键技术方法：①SAROPS蒙特卡洛（Monte Carlo）粒子模拟（默认每次5000个粒子，以随机游走（random walk）表征风、流与leeway的不确定性，基准对象与测试对象各分配2500个粒子）；②以第一标准差椭圆（first standard deviation ellipse）定义搜索区域，采用Jaccard指数（Jaccard index）量化重叠率，将36小时漂移中重叠率首次低于50%的时间作为时间界限；③对10,900个椭圆面积拟合经验幂律函数，独立估计系数c_L、时间指数p_L与风速指数q_L；④建立两级线性回归模型预测质心分离距离；⑤以Grubbs检验（Grubbs test）识别异常值；⑥基于美国国家数字预报数据库（NDFD，风场）与缅因州Finite Volume Coastal Ocean Model（FVCOM，流场）等真实环境数据完成20次验证模拟。

**研究结果**

2.1 可接受leeway误差的定义（Defining Acceptable Leeway Error）：leeway rate变化改变施加于物体的风强迫量并改变预测位置；以第一标准差椭圆作为搜索区域的代理，用Jaccard指数J_n刻画重叠，当J_n＜0.5即重叠不足50%时判定为不可接受误差。

3.1 搜索区域增长与时间、风速及leeway的关系：流速不改变搜索区域大小但改变其位置；搜索区域总体随leeway rate增大而增大，5节风速下每1% leeway对应的面积增长差为16.7 km²，25节下达230 km²；仅20节风、1节流组合出现异常值（Grubbs统计量1.50＞临界值1.48），原因未知。幂律拟合显示时间指数p_L由0.48% leeway时的1.52升至约2%时的1.76后趋稳；风速指数q_L由0.495升至1.544，表明高leeway rate下面积增长对风速更敏感；36小时拟合面积由约545 km²（0.48% leeway）增至2470 km²（9.98%），基准值（5.98%）为1674 km²；5节风速下拟合系统性低估（NMAE=38.22%，归一化平均绝对误差），高于风速时为2.18%–5.28%。

3.2 可接受leeway rate界限（Acceptable Leeway Rate Bounds）：在所有风速下，相对基准值（5.98%）±1%的leeway rate变化在36小时漂移中始终保持至少50%重叠，因此±1%（占风速）为一般可接受误差。

3.3 粒子分布分离的预测（Predicting Separation of Particle Distribution）：质心逐小时分离距离与|LW_base?LW_test|呈强线性（R²=0.996–0.999），风速依赖斜率拟合R²=1.0；以真实环境数据验证，全感兴趣区域（AOI）时空平均风速预测最佳（MAE=0.21 km，平均绝对误差），表明该模型在真实时空变化风场下仍可合理预测分离。

3.4 预测模型测试：以无已知leeway系数的气胀式充气筏（inflatable dinghy）为测试对象，救生筏（+2.04%）与cuddy cabin（+1.43%）的leeway rate均超出±1%界限，搜索区域重叠在1.17小时内即降至50%以下；将充气筏发散角调整为与替代物一致后，cuddy cabin始终至少50%重叠，救生筏在7.67小时后才降至50%以下，表明发散角显著影响分离进程。

**讨论与结论**

讨论指出，此前leeway rate计算到两位小数，而本研究表明漂移预测在±1%风速的leeway rate误差范围内是稳健的，据此确立了leeway计算的最低精度要求，可辅助响应人员在缺乏物体实测数据时选择合适的替代物；研究仅限于leeway rate，未涉及发散角与DWL、CWL分量，充气筏案例表明还需建立发散角的可接受误差限值；±1%界限也为水动力模型（hydrodynamic model）等基于历史数据的计算方法提供了精度基准。研究结论为：leeway rate可允许相当于风速±1%的误差，只要其他变量相等，预测搜索区域仍可覆盖物体的实际位置。