带通海减动结构的圆筒型FDPSO 阻尼系数计算

2023-07-22黄孟丽陈林烽

舰船科学技术 2023年11期

张博，李清，黄孟丽，陈林烽

(1. 江苏科技大学船舶与海洋工程学院，江苏镇江 212003；2. 上海外高桥造船海洋工程有限公司，上海 200131)

0 引言

浮式钻井生产储卸油装置(floating drilling production storage and offloading system，FDPSO)是21 世纪海上油气开采的重要装备之一[1]，具有建造成本低、储油能力强、能适应多种海况、转移方便可重复使用等优点[2]。我国南海存在许多储量不高、分布散乱的边际油田，FDPSO 能够极大改善我国对边际油田开发不足的状况[3]。

从结构类型上可以将FDPSO 分为船型和圆筒形，但是船型FDPSO 存在垂荡和横摇运动性能差、船体疲劳问题严重、经常需要维修保养影响生产效率等缺点，而圆通型FDPSO 能够改善横摇和纵摇运动性能，在运动性能上圆筒型FDPSO 要优于船型FDPSO，但两者垂荡性能都较差[4]。在抑制圆筒型FDPSO 垂荡运动方面，曲志森等[5–6]提出了一种带通海延伸筒体和减动结构的新型FDPSO，能够有效抑制其垂荡运动，但是在利用水动力计算软件分析其阻尼系数时，由于通海孔尺寸远远小于FDPSO 基本尺寸，在数值仿真计算时存在局部网格加密区域多、网格数量大大增加导致计算效率低、计算结果误差增大等问题。

本文提出用无通海孔 FDPSO 代替有通海孔FDPSO计算垂荡阻尼系数的方法，并在理论上验证了其可行性。通过数值计算验证有通海孔 FDPSO和无通海孔FDPSO的等效垂荡质量相等，因而可以用无通海孔FDPSO代替有通海孔FDPSO进行后续阻尼系数分析。

1 分析理论

1.1 FDPSO 运动理论

在只考虑浮体单自由度运动的情况下，基于达朗贝尔原理的直接平衡法建立FDPSO 运动方程[7]：

式中：m为浮体的质量；c为阻尼系数；k为刚度系数；F为浮体的干扰力。

当FDPSO 做垂荡自由衰减运动时，干扰力F由粘性流体水动力Fv和垂荡附加质量引起的附加质量惯性力Fa构成。

式中，ma为 FDPSO 的垂荡附加质量。

将式（2）和式（3）代入式（1），可得到FDPSO垂荡衰减运动方程为：

有通海孔FDPSO 垂荡衰减运动方程为：

无通海孔FDPSO 垂荡衰减运动方程为：

式中：m1为有通海孔FDPSO 的排水质量，m2为无通海孔FDPSO 的排水质量，m1a为有通海孔FDPSO 的垂荡附加质量，m2a为无通海孔FDPSO 的垂荡附加质量。

式中，mδ为FDPSO 等效垂荡质量。

1.2 固有周期

FDPSO 固有周期的表达式可表示为[8]：

在计算FDPSO 各自由度的固有周期时，都可以按照式（8）计算。其中，j取1～6 分别为了垂荡、纵荡、纵摇、横荡、横摇、首摇；Tj为了某自由度方向的固有周期；kj为了某自由度方向的刚度；Mj为了某自由度方向的等效质量，由于本文只分析垂荡运动，Mj即为垂荡等效质量mδ。

2 FDPSO 有无通海孔对比

2.1 分析模型

以一种带通海结构的圆筒型FDPSO 为分析对象，如图1 所示。在传统圆筒型FDPSO 主筒体下部增加延伸筒体和减动结构。有通海孔FDPSO 如图1(a) 和图1(c)所示，通海孔直径为2 m，通海孔位于减动结构顶部、底部以及延伸筒体底部，通气孔贯穿主筒体连通延伸筒体，减动结构上下通海孔以及延伸筒体上通气孔和下通海孔在垂直方向上不在同一条直线上。在进行水动力数值计算时，有通海孔FDPSO 下部延伸筒体和减动结构视为中空透水结构，排水量可以忽略不记。无通海孔FDPSO 如图1(b)和图1(d)所示，在进行数值计算时可将整体构型视为实体结构，排水量较大，使浮体质量大大增加。2 种模型的基础参数如表1 所示。

表1 物理参数Tab. 1 Physical parameters

图1 FDPSO 模型示意图Fig. 1 Schematic diagram of FDPSO model

2.2 计算结果

在水动力仿真软件Aqwa 中对有通海孔和无通海孔的新型圆筒型FDPSO 进行频域分析。频域计算时，取波浪周期5～64 s，间隔1 s，共60 个波浪周期。计算结构垂荡运动的附加质量和运动响应传递函数。结果如图2 和图3 所示。

图2 运动响应传递函数（RAOs）Fig. 2 Motion response transfer function (RAOs)

图3 垂荡附加质量Fig. 3 Heave additional mass

运动响应传递函数（RAOs）和垂荡附加质量计算结果，如表2 所示。

表2 频域计算结果分析Tab. 2 Analysis of frequency domain calculation results

可知，无通海孔的FDPSO 垂荡附加质量m2a稳定在2.17×108kg，有通海孔的FDPSO 垂荡附加质量m1a稳定在3.11×108kg。而无通海孔的FDPSO 排水量m2为2.25×108kg，有通海孔的FDPSO 排水量m1为1.22×108kg。将FDPSO 的垂荡附加质量和排水量相加，得到垂荡等效质量，无通海孔FDPSO 的垂荡等效质量m2δ为4.42×108kg，有通海孔FDPSO 的垂荡等效质量m1δ为4.33×108kg，即2 种模型的垂荡等效质量几乎相等，m1δ≈m2δ。另外，计算结果显示有通海孔和无通海孔垂荡固有周期T 一样，均为19 s。结合固有周期的定义可知，2 种模型在运动过程中的刚度系数是一样的。

当有通海孔FDPSO 和无通海孔FDPSO 的基础尺寸一致时，由于通海孔处流体流动对整体结构水动力Fv的求解影响很小，2 种模型在运动中形成的粘性流体水动力几乎相同，即

结合以上2 种模型垂荡等效质量和固有周期的分析结果，可以得出2 种模型的垂荡衰减运动控制方程几乎是等效的。因此，在后续采用计算流体力学方法分析有通海孔FDPSO 的阻尼系数时，可选择分析更简单、准确度更高的无通海孔FDPSO 代替有通海孔FDPSO。

3 FDPSO 垂荡阻尼分析

利用水动力仿真软件STAR-CCM+对新型圆筒型FDPSO 的垂荡方向上运动衰减过程进行仿真计算，分析其运动情况和受力情况，通过自由衰减时历图分析其阻尼系数[9]。在分析有通海孔FDPSO 的阻尼系数时，选择更简单、准确度更高的无通海孔FDPSO代替有通海孔FDPSO。

3.1 分析模型

为了与文献[5]曲志森的试验测量结果进行对比，以无通海孔FDPSO 实际尺寸按照1∶77.8 的缩尺比建立模型，模型尺寸如表3 所示。基于SolidWorks 建立的平台三维几何模型如图4 所示。

表3 物理参数Tab. 3 Physical parameters

图4 FDPSO 几何建模Fig. 4 FDPSO geometric modeling

3.2 数值仿真计算

3.2.1 计算水域的建立

使用计算流体力学软件STAR-CCM+对FDPSO 模型进行研究，计算水域是以模型质心为中心建立，其长10.0 m、宽4.0 m、高5 m，如图5 所示。

图5 计算域尺寸Fig. 5 Calculation domain size

3.2.2 物理模型

在FDPSO 垂荡阻尼系数的计算中，采用粘性流体的基本理论计算FDPSO 的水动力，即求解不可压缩流体的连续性和Navier-Stokes 方程[10]。计算域内的水和空气两相流体采用欧拉-欧拉模型(Euler-Euler Model)进行描述，并采用求解基于VOF 方法的流体体积输运方程来捕捉水和空气的界面位置，即

式中：U为流体速度，Ur是用于压缩界面的速度场；α为流体体积分数。

3.2.3 网格划分

网格划分采用动网格技术，主要涉及计算区域和运动区域2 个区域的网格划分。计算区域的网格划分采用切割体棱柱层体网格技术，同时对流场复杂区域的网格进行细化加密。对自由液面进行3 层加密，其网格尺寸x方向和y方向均为基础尺寸的20%，z方向分别为基础尺寸的5%，10%，15%，如图6(a)所示。对减动结构间隙处的网格单独细化加密，网格尺寸x，y，z方向均为基础尺寸的2%，如图6(b)所示。对月池部分的网格单独细化加密，网格尺寸x，y，z方向均为基础尺寸的5%，如图6(c)所示。

图6 网格划分Fig. 6 Meshing

3.3 数值仿真结果及阻尼分析

计算中，使浮体初始位置低于其质心0.1 m，即初始液面比浮体标定水线高0.1 m，让浮体能够在不受外力的条件下做垂荡自由衰减运动。

3.3.1 数值计算结果

FDPSO 上升过程和下降过程中速度最大时的局部流场如图7 所示。

图7 FDPSO 上升和下降过程速度最大时局部流场图Fig. 7 Local flow field diagram of FDPSO at maximum speed during rising and falling process

可以看出：上升时，减动结构外侧和通海间隙下方存在涡；下降时，减动结构外侧和通海间隙上方存在涡。运动过程中，流体从减动结构间隙通过，与间隙外流场产生速度差，形成涡。当FDPSO 在流体中运动时，由于流体存在粘性，带动周围流体一起运动，此时会受到流体的粘性剪切力，粘性剪切力方向与其运动方向相反，并随着其运动速度增大而增大，所以流体的粘性剪切力有抑制运动的效果[11]，导致FDPSO 垂向运动持续衰减。

FDPSO 垂荡自由衰减运动位移和受力如图8 所示。