关于使用双体船对浮动式风力涡轮机进行可达性概率评估的学术研究报告
本文旨在向各位研究人员介绍一项发表于2016年《Energy Procedia》第94卷，题为《Probabilistic assessment of floating wind turbine access by catamaran vessel》（使用双体船对浮动式风力涡轮机进行可达性概率评估）的研究工作。该研究由来自西班牙坎塔布里亚环境水力学研究所（Environmental Hydraulics Institute “IH Cantabria”）的Michele Martini, Alfonso Jurado, Raúl Guanche* 和 Iñigo Losada 共同完成。这项研究在2016年1月于挪威特隆赫姆举行的第13届深海海上风电研发会议（EERA DeepWind’2016）上被发表，并经过同行评审。
一、 学术背景与研究目的
本研究属于海洋可再生能源，特别是海上浮动式风力发电技术领域，聚焦于其运行与维护（O&M）中的关键挑战——人员可达性。随着海上风电向更深水域发展，浮动式风力涡轮机（Floating Wind Turbines）展现出巨大潜力，因其可在固定基础不经济的深水区利用更稳定、更强劲的风能。然而，浮动平台的动态运动特性为其运行维护带来了新的复杂性和风险，其中技术人员如何安全、可靠地抵达平台进行巡检和轻型维修是核心问题之一。通常，人员通过配备护舷（fender）的双体船船员转运船（Catamaran Crew Transfer Vessels, CTVs）以“步行上班”（walk-to-work）模式登临平台。目前，相关的安全指南和建模工具主要针对固定式基础风机，缺乏专门用于评估浮动平台人员可达性的分析模型。
因此，本研究的动机源于三点：1) 海上风电，尤其是浮动式风电的市场发展趋势与潜力；2) 海上风电场运行维护的高成本和高风险；3) 缺乏针对浮动平台“步行上班”可达性的建模工具。本研究的主要目标是开发一种数值方法，用于评估双体船在“步行上班”配置下对浮动式风力涡轮机的可达性（accessibility），即计算在特定海况下能够安全进行人员转移的时间百分比。
二、 详细研究流程与方法论
本研究采用了一种系统的频率域概率评估方法，主要流程可分为以下几个步骤：
1. 建立约束多体水动力系统模型： 研究将双体船（CTV）和浮动式风力涡轮机平台（本研究选用OC4半潜式平台）视为一个受约束的多体系统。每个物体具有六个自由度。研究在频率域内建立了该系统的运动方程。水动力系数（附加质量、势阻尼、静水恢复力、波浪激励力）使用势流求解器HydroD（DNV Sesam软件）进行计算。为了在频率域模型中考虑非线性效应，研究对两个关键非线性力进行了线性化处理：a) 通过谐波线性化方法计算了代表粘性阻力的等效线性阻尼矩阵；b) 使用有限差分法和准静态悬链线方程，在船舶系柱推力与系泊系统反力平衡的平衡点附近，对平台系泊力进行了线性化，得到了等效线性刚度矩阵。系统的核心约束是模拟双体船通过护舷与平台登陆结构接触的“球铰”连接，该约束不允许两个物体在接触点发生相对平移。通过引入约束矩阵，求解运动方程，可以得到系统各自由度运动以及接触点反作用力相对于波浪激励的传递函数。
2. 定义可达性判据与响应变量： 安全访问需要满足两个物理条件：第一，护舷与平台接触面不发生滑动；第二，两物体之间的相对转动角度在允许的安全限度内。基于此，研究定义了四个关键的响应变量来量化这些条件： * α 和 β： 这两个变量由接触点的水平和垂向反作用力（λ1 和 λ3）组合而成，分别用于判断护舷是否会向上滑动（α）或向下滑动（β）。无滑动条件要求 |α(t)| 和 |β(t)| 均小于（双体船系柱推力 × 静摩擦系数）。 * 相对横摇角 δρ 和相对艏摇角 δψ： 从船舶坐标系观察，两物体之间的相对横摇和艏摇角度。人员转移要求这些相对角度不超过设定的极限值（本研究设为5度）。 通过步骤1中获得的传递函数，可以计算出这四个响应变量对波浪的传递函数。
3. 短期海况极端最大值计算： 对于一个短期的平稳海况，假设波浪 elevation 为高斯过程，则线性系统的响应也是高斯过程，其波峰高度服从瑞利（Rayleigh）分布。对于每个响应变量x，首先通过其传递函数和波浪谱计算响应方差。然后，基于瑞利分布和响应波峰数量n，可以推导出在给定置信因子（confidence factor, fe，本研究取0.95）下，该海况持续时间内响应波峰极端最大值（即所有波峰都低于此值的概率为fe）的解析表达式。此步骤为每个海况下的每个响应变量计算出一个用于判断的极端值。
4. 可达性判定与长期统计分析： 对于每一个小时的海况，计算上述四个响应变量的极端最大值。将这些极端值与各自对应的阈值（基于摩擦系数、系柱推力和最大允许相对转角）进行比较。只有当四个条件同时满足时，才判定该小时海况下访问是可能的。研究将这一判定流程应用于苏格兰阿伯丁附近一个特定站点（水深90米）从1980年至2013年共34年的 hindcast（后报）波浪数据。数据来自IH Cantabria的全球海洋波浪（GOW）再分析数据库，包含每小时的有效波高（Hs）、谱峰周期（Tp）、平均波向（θm）和平均方向分布（σθ）等参数。通过统计所有可访问的小时数占总小时数的比例，得到长期的年平均可达性以及月度可达性变化等统计结果。
三、 主要研究结果
1. 系统响应特性分析： 通过求解约束多体系统方程，研究获得了双体船和OC4平台的运动幅值响应算子（RAOs）以及接触力、α、β等变量的RAOs。结果明确显示，由于两个物体在惯量和 submerged volume 上的巨大差异，它们的固有动力响应不同（CTV对短波敏感，OC4对长波敏感）。然而，在约束条件下，平台的运动显著影响了CTV的响应，反之则影响较小。接触力在平台运动较大的周期附近达到峰值，但值得注意的是，最大响应发生的波浪周期（16.5秒和24.0秒）并不完全对应于OC4的固有周期（18.0秒和27.5秒），这表明是两物体之间的相对运动，而非平台的绝对运动，主导着接触力和系统响应。分析还发现，对于短周期波浪，α的响应大于β，意味着此时发生向上滑动的风险更高。
2. 规则波下的极限波高： 研究绘制了在规则波假设下，能够实现安全访问的极限波高极坐标图。结果表明：对于短周期波浪，横浪（beam sea）是最佳访问方向，因为此时接触力较低；对于长周期波浪，顶浪（head sea）由于平台对波浪的遮蔽或吸收效应，比随浪（following sea）更有利于访问。在对应于OC4垂荡运动强响应的波浪周期附近，无论波向如何，可达性都很差。这些发现与针对固定式风机的研究结论有相似之处，但也揭示了浮动平台特有的动力耦合效应。
3. 长期可达性概率评估： 应用所提出的概率方法对34年的历史波浪数据进行计算后，得到该站点使用指定CTV访问OC4平台的年平均可达性为23.7%。可达性表现出强烈的季节性变化和显著的波动性： * 月度变化： 7月份可达性最高（平均约50%），1月份最低（平均约9%）。这与该海域波高的月度变化趋势一致（7月波高最小且最稳定，1月波高最大且变异性强）。 * 变异性： 可达性估计存在很大的不确定性（标准差较大）。有趣的是，夏季月份的可达性变异性比冬季月份更显著。这是因为冬季平均波高远高于CTV的承受极限，其波动对可达性概率影响不大；而夏季平均波高处于临界范围附近，微小的波动就会导致可达性概率发生剧烈变化。 * 与波高的关系： 整体上，可达性随有效波高Hs的增加而降低。当Hs在1.5米左右时，结果最为分散，此时波浪方向、分布和峰值周期等其他参数对船舶响应的影响变得显著。当波高非常低或非常高时，访问总是或永远不可能，其他波浪参数的影响很小。
四、 研究结论与价值
本研究成功提出并应用了一套完整的频率域概率方法，用于评估双体船对半潜式浮动风力涡轮机的“步行上班”可达性。研究结论指出，CTV的响应显著受到平台运动的影响，尤其是在平台自身运动较大的情况下。接触力的大小主要由两物体若自由运动时将产生的相对运动幅度决定。通过对苏格兰一个代表性站点34年数据的分析，量化了该地的年平均和月度可达性水平，并揭示了其高度的季节性和年际变异性。
这项研究的科学价值在于，它首次将约束多体水动力学系统分析与概率统计方法相结合，专门用于解决浮动式海上风电结构的人员访问评估问题，填补了该领域建模工具的空白。在应用价值上，该方法为海上风电场运营商和设计者提供了关键的决策支持工具，可用于： * 优化运维策略： 基于可达性概率的月度变化，规划最佳的维护时间窗口。 * 评估不同船舶/平台组合： 比较不同CTV或不同浮动平台设计的可操作性。 * 站点适应性评估： 在风电场选址和可行性研究阶段，评估特定海域的运维可达性。 * 降低风险和成本： 通过科学预测可工作时间，减少因天气窗口误判导致的计划延误和成本超支，并提升人员转移的安全性。
五、 研究亮点
六、 其他有价值的发现
研究在讨论中指出，当前的频率域线性模型存在局限性，未来的工作可以考虑纳入二阶波浪力并研究大位移情况下的响应。此外，研究发现的对于短波横浪更优、对于长波顶浪更优的访问方向规律，为实际操作中的船舶操控提供了有价值的参考。研究案例中选用的OC4平台和特定CTV的参数，也为同类研究提供了一个可比的基准。最后，作者在致谢中提到了该项目获得了欧盟第七框架计划等资助，体现了该研究属于大型国际合作项目的一部分，其成果具有更广泛的研发背景和应用前景。