漂浮式海洋牧场养殖装置系泊系统设计

海洋中蕴含着丰富的自然资源和能源，有用于发电的风能、波浪能和潮流能，还有用于养殖的渔业资源。随着人们对生活水平需求的不断提升，海洋渔业的市场也在不断扩展。为了满足不断增长的资源需求，越来越多的国家将养殖场地由海岸转移至深海，而随着海洋平台技术的不断发展，基于海洋能源平台建立的“海洋牧场”是推进养殖产业转型发展的重要方向。

近年来，越来越多的国家重视发展多功能海洋平台，以半潜式海洋平台结构为基础，综合风能-波浪能-潮流能海上互补发电海洋平台，深水钻井平台、海洋牧场都是国际上主要探讨和研究的话题。2012年，我国第六代深水半潜式钻井平台投入使用[1]，其最大作业水深达到2451m。2011年，日本开发制造了“Wind Lens”平台[2]，整个平台安装了2台风力机波浪能装置和太阳能。2018年，世界第一座半潜式智能装备海洋渔场“海洋渔场1号”于挪威海弗鲁湾海域中投放[3]，整个装置具备挪威先进的智能养殖技术同时结合了中国成熟的海洋平台建造技术，总高69m，直径110m，可抗12级台风。渔场安装有2万多个各类传感器以及100多个监控设备，在鱼苗投放、喂食、实时监控和渔网清洗等方面都实现了智能化和自动化。

在系泊方面，KIM等[4]针对一种FPSO进行系泊系统时域耦合，分析了不同风浪下的浮体运动响应和系泊动力分析，并与试验数据作对比。TANG等[5]通过建立网箱的时域数值模型，分析破损系泊系统下网箱的运动情况及系泊力的变化。LIN等[6]通过AQWA软件对半潜式浮式平台进行水动力性能及系泊系统分析，并研究了系泊对平台水动力的影响。结合前人所研究的结论，马勇等[7]考虑了水轮机、风力机与平台的相互作用；郭小天等[8]针对潮流能发电装置在各种外载荷下的运动性能，合理地设计了适用于潮流能发电站的弹性索-锚链组合系泊系统。周丙浩等[9]利用Fortran对AQWA进行二次开发，研究风力机、水轮机与平台的耦合运动效应。

本文以“海洋渔场1号”为母型设计了一型半潜式海洋牧场养殖装置，结合海洋能源利用和渔业养殖，在海洋牧场上安装风力机和潮流能水轮机。根据平台结构设计的2种系泊系统，利用间接时域分析法进行时域水动力性能的分析。研究成果可为海洋牧场系泊系统设计提供有效依据，为未来海洋牧场与能源综合利用平台的设计建造提供参考。

计算理论及模型

1.1 时域计算方法

对于研究浮体在波浪下的运动问题，首先需要求解流场速度势[10]。通过线型叠加入射势、辐射势及绕射势，以此来表示浮体周围流场的总速度势：

微信图片_20220630164549.png

式中：Φ₁(x,y,z,t)为入射波速度势；Φ_R(x,y,z,t)为辐射速度势；Φ_D(x,y,z,t)为绕射速度势；(x,y,z)为流场中的位置坐标；t为时间。

结合运动学边界条件、动力学条件、远场边界条件以及物体表面边界条件，最终得到总速度势定解条件，见式（2），

微信图片_20220630164552.png

其中：n为物面法向量；i为虚数单位；g为重力加速度；k为波数，满足k=ω²/g；v_n为浮体湿表面S₀的法向速度；ω为波浪频率。

式（2）为速度势的定解条件，入射波的速度势已知后，即可求得辐射速度势和绕射速度势。将浮体所受波浪力分成3种力，即入射、绕射和辐射，而入射波浪力与绕射波浪力构成波浪激励力。通过所得辐射势和绕射势，求解浮体的波浪力，进而建立浮体的频域运动方程，见式（3）。

微信图片_20220630164554.png

式（3）中：m为浮体质量；λ_a为额外阻尼系数；μij为附加质量；λ_ij为辐射阻尼；K为静水回复系数；K_a为额外刚度系数；x_i为浮体的位移；F为波浪激励力。

之后将浮体的频域运动方程转换为时域运动方程。当浮体进行六自由度微幅运动时，可以把运动看成一系列脉冲运动。将浮体周围的速度势进行叠加求解。浮体运动下，流体辐射势可以表示为

微信图片_20220630164556.png

式中：x_j(t) 为浮体在j方向运动的速度或者是角速度；ψ_j为浮体单位脉冲运动产生的速度势，且在j方向上；χ_j为浮体做单位脉冲运动后t时刻的速度势。

通过伯努利方程求得浮体表面动压力并进行湿表面积分，得到浮体辐射水动力和力矩：

微信图片_20220630164558.png

式中：m_ij为附加质量矩阵；K_ij为延迟函数；X_j(t)为浮体在j方向运动的加速度或者是角加速度。进而推导出时域运动方程表达式为

微信图片_20220630164559.png

式中：X_i(t)为第i个自由度上的波浪激励力；C_ij为浮体静水回复力矩阵；M_ij惯性矩阵。通过对式（6）换元得到：

微信图片_20220630164601.png

在求解时域运动方程前课通过面元法，获得浮体的频域水动力系数。在通过频域水动力系数求解时域运动方程的水动力系数。同时波浪激励力也需要进行频域到时域的转换，这种方法称为间接时域法。在建立浮体的时域运动方程后，利用数值方法4阶龙格-库塔求解方程，最终得到浮体各自由度时历运动数据。

1.2 计算模型

本文中的漂浮式海洋牧场养殖装置的母型是由中国船舶集团有限公司建造的“海洋渔场1号”养殖网箱，整个平台基础结构式为半潜式，主要分为四大构件：上层平台、立柱、网箱结构和下浮体。表1为海洋牧场主尺度；

微信图片_20220630164604.png

图1为漂浮式海洋牧场养殖装置的总体结构图；图2为海洋牧场水动力计算模型，其中，立柱和下浮体基于面元法计算，而横撑由于属于小尺度构建，采用莫里森单元进行计算。

系泊设计方案说明

2.1 环境参数

假定本文海洋牧场所处位置为黄海，该海域接近中国重要港口城市，具有运输方便的优点。综合水文环境资料，作业工况和生存工况环境参数见表2。

微信图片_20220630164612.png

2.2风力机载荷

本文的风力机使用沈阳工业大学开发的 200 kW 风力机，不同风速下风力机所受推力（见图 3）和侧向力（见图 4）通过查阅文献所得[11]。根据所得数据将风力机在相应风速下所得载荷按照定常力和力矩的形式加入数值模型中。

微信图片_20220630164616.png

微信图片_20220630164618.png

2.3 水轮机载荷

水轮机选择翼型为NACA0018、3叶片的 5kW 潮流能水轮机。为了确定水轮机载荷，通过STAR-CCM+ 数值模拟得到本文水轮机侧向力（见图 5）和推力曲线（见图 6）。根据所得数据将水轮机在相应流速下所得载荷按照定常力和力矩的形式加入数值模型中。

微信图片_20220630164621.png

微信图片_20220630164623.png

2.4 系泊设计方案

整个海洋牧场的系泊系统设计为悬链线形式。整个系泊系统的设计不仅要考虑浮体的运动，还需要考虑系泊材料的属性等性能。本文将对比纯锚链方案和锚链-钢缆-锚链组合式系泊方案，系泊材料参数见表3。

微信图片_20220630164626.png

根据海洋牧场形式，浮体的固定系泊点设置在结构物浮箱底部，每个系泊点连接两条锚链并与海底锚点连接。锚链对称分布，并且两两平行。海域水深为100m。设海洋牧场水平面中心为坐标原点，系泊线布置方式及系泊点数据具体坐标见图7和表4。

2.5 系泊设计要求

系泊系统的设计不仅要满足浮体所规定的位移要求，还需要满足系泊线的强度要求。在作业工况和生存工况下需要保证浮体不会发生倾覆，系泊系统不会断裂等情况。中国船社级的《海上移动平台入级规范》对锚链不同状态下张力的安全系数要求见表5。

微信图片_20220630164634.png

由于海洋牧场上并没有长期居住的作业人员，也没有对运动响应有严格要求的装置，因此运动响应只需要保证在作业工况下装置的正常工作，生存工况下系泊达到规定的安全系数即可。而对于合理的位移偏移量，也没有具体的设计规范，因此通过查数据以及经验，要求作业工况下，平台的纵荡和横荡运动限定在水深的10%左右。由于风力机具有自动偏航功能，横摇和艏摇对发电效率影响小，仅考虑纵摇影响。限定作业工况下纵摇角度小于5°。

算例分析

在作业工况和生存工况下，分别对海洋牧场的载体及其系泊系统进行时域耦合分析。在作业工况下，海洋牧场需要考虑风力机和水轮机作业时所受推力和力矩。而在生存工况下，海洋能装置停止工作，仅需要考虑海况对平台和系泊系统的影响。分析两种工况下海上牧场的六自由度运动和系泊张力。

3.1 作业工况结果分析

3.1.1 纯锚链方案

设计流向与风向为0°，设置不同的浪向角，研究海洋牧场在纯锚链系泊方案下3h内的六自由度运动情况以及系泊系统张力变化，统计结果见表6。

微信图片_20220630164638.png

观察表6数据可知：作业工况下，平台在浪向角为0°时的纵荡最大，为5.68m，符合浮体位移要求。而平台在其他自由度的运动数值较小，这是由于浪向角在0°时，浪向与风流同向，对平台纵荡方向的影响较大。

作业工况下，纯锚链系泊在各浪向角下的六自由度最大摇荡的时历曲线见图8，系泊系统的最大张力时历变化曲线见图9。由图8和图9可知：纯锚链系泊方案中，系泊锚链张力最大为665kN，最小安全系数为7.750，锚泊系统符合规范要求。在作业工况下，海洋牧场上的风力机与水轮机受到风速和流速的作用，在风速流速方向产生较大位移。