基于统计方法的数值模拟再审视MV Estonia海难

《Journal of Marine Science and Engineering》：Revisiting the MV Estonia Accident Using Numerical Simulations with a Statistical Approach Shinwoong Kim and Petri Valanto

　　

1994年9月爱沙尼亚号渡轮（MV Estonia）失事以来，多家机构对该事故进行了调查。事故根因被认为已知， consequent 沉没过程也已建立。然而，2020年9月，水下遥控潜水器（ROV）拍摄的新视频记录显示，沉于海底的MV Estonia残骸右舷存在一处此前未知的穿透性损伤。基于这一新证据，爱沙尼亚安全调查局（ESIB）于2020年10月启动了对MV Estonia事故新信息的初步评估。该右舷新侧向损伤（NSD）是在MV Estonia尚处于水面漂浮状态时即已存在，还是由随后下沉过程中船体与海底碰撞所致，成为需要澄清的重要问题：若为前者，此前研究中提出的事故根因结论的有效性可能被证明为过早。汉堡船模试验池（HSVA）开展的本项调查目标之一即在于阐明此问题：研究呈现了针对多种损伤构型在波浪中进行沉没过程数值模拟的结果，不仅采用单一模拟，还运用了统计方法。

研究人员采用HSVA Rolls数值代码对受损船舶在波浪中的运动及舱室进出流进行时间域模拟。该程序将船舶视为六自由度系统，在定常海况中以给定平均航向角航行。车辆甲板（vehicle deck）上flooding水的动态运动通过浅水方程（shallow water equations）求解。为应对实际遭遇波序列及中心围壁（center casing）舱门状态未知带来的不确定性，研究人员开发了统计方法：对每种损伤构型，采用100组随机生成的波序列进行模拟，计算不同船舶行为类型的发生概率。通过调整中心围壁门排放系数（discharge coefficients），使多数模拟结果呈现与MV Estonia事故相符的行为特征（Category 3型行为，即船舶完成转向后在尾斜浪中因 hull 多处开口progressive进水而逐渐横倾最终倾覆）。概率超越函数（probability of exceedance, POE）用于评估倾覆时间分布，包括无条件POE及基于Jeffreys方法的95%置信区间。

研究发现，在"参考"中心围壁门条件下（2/3门在横倾角>40°时开启），无NSD时100次模拟中50次倾覆，其中14次（14%）落入实际事故时间范围（950–1500秒）；而加入NSD后，69次倾覆中无一次落入该时间范围，绝大多数在转向阶段（235–665秒）即告倾覆。仅NSD损伤的模拟显示，船舶在直航或新轨迹下均保持较小横倾，与实际事故明显不符。对比研究表明，减速比转向更能有效提高船舶 survivability。综合统计结果与个别案例，NSD在MV Estonia尚漂浮时即存在的概率极低，该损伤更可能产生于沉船与海底的碰撞。研究结论支持传统认知的仅 bow 损伤场景为最可能事故情景，并揭示了船员应急决策对船舶 survivability 的关键影响。

**一、研究背景与问题提出**

1994年9月28日，爱沙尼亚号渡轮（MV Estonia）在从爱沙尼亚塔林 freight 往瑞典斯德哥尔摩途中沉没，造成989人中仅137人幸存。联合事故调查委员会（JAIC）于1997年发布最终报告，认为事故根因为船首 visor 脱落导致车辆甲板（vehicle deck）进水。此后，SSPA联合体与HSVA联合体在2006–2008年间分别进行EDA独立调查，总体确认了JAIC的事故情景重构。

然而，2020年9月发现的新证据改变了调查格局：水下ROV视频显示，MV Estonia残骸右舷水线附近存在两处此前未知的穿透性损伤，其中位于82–83号肋骨之间、连通车辆甲板、舱室区及右舷压载舱的新侧向损伤（New Side Damage, NSD）尤为关键。该NSD总面积约1.239 m2（车辆甲板开口0.54 m2、舱室区0.60 m2、压载舱0.099 m2），其成因与时间成为争议焦点——若NSD在船舶尚漂浮时即已存在，将根本挑战JAIC及后续研究确立的事故认知体系。

ESIB因此委托HSVA与基尔应用科学大学开展数值模拟研究，核心科学问题为：NSD是船舶漂浮时已存在，还是沉船与海底碰撞后产生？此问题涉及对已有事故结论有效性的根本性质疑，具有重大技术与法律意义。研究人员于《Journal of Marine Science and Engineering》发表本项研究，运用升级的HSVA Rolls代码，结合统计方法系统评估了多种损伤场景。

**二、关键技术方法**

研究采用HSVA Rolls时域数值模拟程序，该程序基于准静力学/水动力学建模，求解船舶六自由度运动与舱室flooding的耦合过程，车辆甲板flooding水运动通过浅水方程描述。关键技术创新在于统计方法的应用：针对JAIC（1997）定义的海况（有义波高H_s 4.2 m、峰值周期T_p 8.3 s、JONSWAP谱），每种损伤配置均采用100组独立随机波序列进行蒙特卡洛式模拟，通过概率超越函数（POE）量化倾覆时间分布特征。

输入数据更新包括：基于2023年高分辨率声纳调查构建的HSVA新轨迹（6003 m，visor脱落至最终沉没时间42分钟，考虑风漂流与0.5节ENE表面流影响）；2024年孤寡实验获得的船桥上层建筑窗户破坏载荷（小窗23.1 m水柱、大窗7.5 m，修正2008年低估值）；ROV调查确认的中心围壁门状态（左舷7门中2门在位、右舷11门中2门在位，其余状态未知）；以及基于ROV观察的1100吨货物1.5 m横向位移（较2008年增加1.1 m）。中心围壁门排放系数经统计微调，使多数模拟呈现与事故相符的Category 3型行为（完成转向后渐进倾覆）。

**三、研究结果**

**3.1 仅 bow 损伤的统计模拟**

"参考"中心围壁门条件下，100次模拟中50次倾覆，分布特征为：33%于转向期间（235–665秒）、3%于转向后初期（665–950秒）、14%落入事故时间范围（950–1500秒）、50% survive。POE分布显示Weibull拟合（形状参数1.74、尺度参数674.10）与实际事故时间高度吻合，表明该配置物理合理。

"全关闭"条件下倾覆概率降至22%，几乎全发生于转向期，无事故时间范围案例，表明该中心围壁门配置不可信。"全开"条件下43次倾覆，5次落入事故时间范围但位于区间边缘，较"参考"条件可信度低。

**3.2 加入NSD的统计模拟**

"参考"条件加NSD后，倾覆次数激增至69次，其中64%发生于转向期，0%落入950–1500秒区间，仅2例晚于1800秒。POE尾部偏离主趋势，Weibull拟合（形状参数0.54、尺度参数66.33）显示倾覆极可能于900秒内完成。"全关闭"与"全开"条件加NSD后分别达70次倾覆，同样无事故时间范围案例。

**3.3 仅NSD损伤模拟**

直航定常速条件下，100次模拟均无倾覆，横倾约10°后稳定；直航减速条件下类似，峰值约25°后恢复。新轨迹减速条件下5次倾覆，但发生于4100秒左右，远晚于实际事故。

**3.4 航速航向对 survivability 的影响**

模拟显示，bow损伤伴NSD直航定常速时100%倾覆（250–900秒），减速后降至65%；bow损伤无NSD直航定常速时100%倾覆，减速后仅20%。转向引入附加动态横倾力矩，降低 survivability；减速显著优于转向作为应急措施。

**3.5 代表性沉船情景分析**

选取与事故最匹配的波序列进行个案比较：新轨迹更新输入与2008年旧轨迹结果高度相似，验证新轨迹可靠性。NSD存在时，无论中心围壁门如何配置，均于4–7分钟内快速倾覆，与幸存者证词重构的渐进过程根本不符。

**四、讨论与结论总结**

研究人员基于统计方法得出以下核心结论：NSD在MV Estonia漂浮期间存在的概率极低，该损伤最可能产生于沉船与海底的碰撞；仅bow损伤的新HSVA轨迹场景为最可能事故情景；中心围壁门"参考"条件（微调后2/3开启概率）为最可信配置；减速优于转向作为应急决策。个别模拟无法改变总体统计趋势，100次随机波序列的蒙特卡洛框架有效克服了单一样本的局限性。

具体结论翻译如下（ adopted 自原文10.2–10.3节）：

"基于新轨迹参考条件无NSD的模拟与实际事故时间匹配最佳，无条件区间POE为14%，条件区间POE达28%。相比之下，其他新轨迹场景（无论有无NSD）的无条件区间POE均为0%。加入NSD后，无论中心围壁门为'参考'、'全关闭'或'全开'条件，无一次模拟倾覆时间落入950–1500秒区间。因此，NSD在基于海底残骸位置构建的轨迹上漂浮期间即已存在的假设极不可信。若船舶保持原直航航线，无论是否减速，该总体结论不变。仅bow损伤的新HSVA轨迹场景为最可能损伤情景。新轨迹较旧轨迹有所改进，但横倾角时间历程保持高度相似。bow损伤仅直航场景虽不能完全排除，但基于海底残骸分布可予排除。仅NSD场景（直航或新轨迹）导致与实际事故明显不同的横倾过程：直航时船舶 survive，新轨迹时5/100次倾覆但远晚于实际事故，故该场景极不可能。"

"在 bow 损伤伴或不伴NSD两种条件下，直航减速均较定常速显著提高船舶 survivability。转向引入附加动态横倾力矩，不利于船舶生存。在经典bow损伤场景中，转向显著降低 survivability。总体而言，保持直航并减速（不转向）似乎为船舶提供最佳生存机会。"