舰船科学技术  2022, Vol. 44 Issue (18): 106-110    DOI: 10.3404/j.issn.1672-7649.2022.18.021   PDF    
半潜式海洋平台阻力预报数值模拟分析
易丛1, 宋春辉1, 李刚1, 陈勇军2, 来敏2     
1. 中海油研究总院有限责任公司,北京 100027;
2. 北京迪玛尔海洋技术有限公司,北京 100085
摘要: 基于STAR-CCM+,采用 $ k - \varepsilon $ 湍流模型对半潜式海洋平台进行三维数值模拟,计算海洋平台及压强分布,并通过与实验结果对比对本文算法进行验证。在此基础上,分析不同来流方向和不同来流速度下平台的阻力及表面压强的变化规律。结果表明:平台阻力受“遮蔽”效应和迎流面面积的影响较大;在相同来流速度下,不同方向的平台阻力随流向角度增大;相同来流方向下,平台阻力随着来流速度的增大而增大。
关键词: 半潜式海洋平台     阻力预报     湍流模型     遮蔽效应    
Numerical simulation analysis of resistance prediction of semi submersible offshore platform
YI Cong1, SONG Chun-hui1, LI Gang1, CHEN Yong-jun2, LAI Min2     
1. CNOOC Research Institute Co., Ltd., Beijing 100027, China;
2. Beijing Dimar Marine Technology Co., Ltd., Beijing 100085, China
Abstract: Based on STAR-CCM +, the three-dimensional numerical simulation of Semi-submerged platform is carried out by using turbulence model. The offshore platform and pressure distribution are calculated, and the algorithm is verified by comparing with the experimental results. On this basis, the variation laws of platform resistance and surface pressure under different inflow directions and velocities are analyzed. The results show that the platform resistance is greatly affected by the shielding effect and the area of the upstream surface. Under the same inflow velocity, the platform resistance in different directions increases with the inflow angle. Under the same inflow direction, the platform resistance increases with the increase of inflow velocity.
Key words: semi-submerged platform     resistance prediction     turbulence model     masking effect    
0 引 言

近年来,半潜式海洋平台因其良好的稳性、强大的抗风浪能力以及高速的作业效率等优点,在深海工程中发挥着不可或缺的作用。但海洋平台在海上作业时,受到风、浪、流等环境因素的作用,结构被破坏甚至发生倾覆等海洋事故时有发生,因此海洋平台阻力预报的准确性对海洋平台的整体安全性尤为重要。

为了实现海洋平台阻力的准确预报,国内外学者对平台阻力计算展开了大量研究。Mataymoto等[1]针对某大水线面立柱半潜平台提出计算其气隙响应与分布的完全非线性方法。Li等[2]基于一种新型深水平台建立平台主体与立管系统耦合模型,计算了风浪流联合作用下的平台的运动响应。范亚丽等[3]针对恶劣海洋环境下的半潜式支持平台,通过分析试验来掌握其运动以及阻力特性。张明霞[4]运用STAR-CCM+平台与OCIMF图谱插值方法对FPSO在多流向角度下的阻力数值进行计算。董斌等[5]在水池试验的基础上利用Ansys-AQWA软件数值仿真,准确模拟出平台在自存工况下的阻力。袁洪涛等[6]为了解系泊系统危险工况,对不同风浪流夹角对海洋平台系泊系统的影响进行研究。陈江华等[7]采用SWAN模型进行不同方向风流浪对平台稳定性影响的分析,确定平台危险工况。

现阶段,学者对半潜式平台环境阻力的研究内容主要集中在风速、浪高与平台本身结构对所受阻力的影响,对不同流向角度下的阻力预报的研究并不多见。因此本文以某半潜式钻井平台为研究对象,基于STAR-CCM+平台对不同流向角度及不同来流速度下的平台阻力和三维表面压力进行计算,并分析来流速度和来流方向对阻力及表面压力的影响,研究结果对平台设计具有一定指导作用。

1 数值方法

计算流体力学(CFD)是一种基于计算数学、近代流体力学和计算机科学3门学科,通过在计算机上模拟流体流动,实现数值计算与可视化处理,从而达到对所求物理现象能够较为准确分析研究的方法[8]。STAR-CCM+作为性能优良的CFD软件,在海洋结构物的水动力和流场模拟方面得到了广泛应用。

1.1 控制方程

本文数值模拟采用雷诺平均N-S方程,其笛卡尔坐标系下的张量表示形式为:

$ \frac{{\partial \rho }}{{\partial t}} + \frac{\partial }{{\partial {x_i}}}\left( {\rho {u_i}} \right) = 0,$ (1)
$ \begin{gathered} \frac{\partial }{{\partial t}}\left( {\rho {u_i}} \right) + \frac{\partial }{{\partial {x_i}}}\left( {\rho {u_i}{u_j}} \right) = - \frac{{\partial p}}{{\partial {x_i}}} + \frac{\partial }{{\partial {x_j}}} \\ \left[ {\mu \left( {\frac{{\partial {u_i}}}{{\partial {x_j}}} + \frac{{\partial {u_j}}}{{\partial {x_i}}} - \frac{2}{3}{\delta _{ij}}\frac{{\partial {u_l}}}{{\partial {x_l}}}} \right)} \right] + \frac{\partial }{{\partial {x_j}}}\left( { - \rho \overline {{{u'}_i}{{u'}_j}} } \right)。\\ \end{gathered}$ (2)

其中: $ {u_i} $ $ {u'_i} $ 分别为平均速度和脉动速度; $ - \rho \overline {{{u'}_i}{{u'}_j}} $ 为雷诺应力,表示湍流的影响。上述2个方程和瞬时N-S方程有相同的形式,只是速度或其他求解变量为时间平均量。

1.2 湍流模型

CFD软件中可选择的湍流模型有 $ k - \varepsilon $ 模型, $ k - \omega $ 模型,RNG $ k - \varepsilon $ 模型等,本文采用 $ k - \varepsilon $ 湍流模型。

湍动能 $ k $ 方程:

$ \rho \frac{{\partial k}}{{\partial t}} = \frac{\partial }{{\partial {x_i}}}\left[ {\left( {\mu + \frac{{{\mu _t}}}{{{\sigma _k}}}} \right)\frac{{\partial k}}{{\partial {x_i}}}} \right] + {G_k} + {G_b} - \rho \varepsilon - {Y_M}。$ (3)

湍流耗散率 $ \varepsilon $ 方程:

$ \rho \frac{{\partial \varepsilon }}{{\partial t}} = \frac{\partial }{{\partial {x_i}}}\left[ {\left( {\mu + \frac{{{\mu _t}}}{{{\sigma _k}}}} \right)\frac{{\partial \varepsilon }}{{\partial {x_i}}}} \right] + {c_{1\varepsilon }}\frac{\varepsilon }{k} + {G_k} + {G_{3\varepsilon }}{G_b} - {C_{2\varepsilon }}\rho \frac{{{\varepsilon ^2}}}{k}。$ (4)

湍流粘性系数方程:

$ {\mu _t} = \frac{{{c_\mu }\rho {k^2}}}{\varepsilon }。$ (5)

式中: $ k $ 为湍动能, $ \varepsilon $ 为湍流耗散率; $ {Y_M} $ 为可压速湍流脉动膨胀对总的耗散率的影响; $ {\mu _t} $ 为湍流粘性系数; $ {G_k} $ 为平均速度梯度所能引起的湍动能; $ {G_b} $ 为浮力影响所能引起的湍动能。

2 计算模型与边界条件 2.1 平台模型

以某半潜式海洋平台模型为研究对象,其水下部分结构如图1所示,平台由2个浮体、立柱以及横撑等模块组成,主要结构尺寸如表1所示。

图 1 半潜式平台水下部分示意图 Fig. 1 Schematic diagram of underwater part of semi submersible platform

表 1 平台模型主要结构尺寸 Tab.1 Main structural dimensions of platform model
2.2 计算域和边界条件

计算域和边界条件的设置如图2所示,计算域为正方体,长宽高分别为4 m×4 m×0.7 m。平台中心距离速度入口2 m,入口边界采用Velocity-inlet,出口边界采用完全发展出流边界条件Outflow[9],平台表面边界采用无滑移壁面边界条件,数值计算区域边界条件设置如图2所示。

图 2 边界条件 Fig. 2 Boundary conditions
2.3 网格划分

在划分网格时,将整个流域划分为若干区域进行网格划分,以此提高网格质量和控制网格数量,流场区域的网格划分如图3所示。为保证数值模拟精度,全域采用六面体结构化网格,平台壁面采用精细网格并保证第一层网格位于粘性底层内。半潜式海洋平台的轮廓特征线进行线控制,同时对平台表面设置面控制,并做好网格的光滑过渡。网格数约为196.8万个,半潜式海洋平台表面网格如图4所示。

图 3 网格区域网格划分 Fig. 3 Grid division of grid area

图 4 半潜式海洋平台网格划分 Fig. 4 Grid division of semi submersible offshore platform
3 计算结果及分析 3.1 模拟精度分析

为验证本文模拟的精度,对不同来流方向下的平台总阻力系数与相关试验结果进行对比,来流速度为1 m/s时对比如图5所示。可以看出,本文计算结果与试验值基本吻合,总阻力系数变化规律关于流向角度180°具有一定的对称性,总阻力系数在流向角度0°~60°时呈现上升趋势,总阻力系数在流向角度60°~75°时几乎不变,而75°~90°时总阻力系数略呈现减小趋势。

图 5 不同流向角度下阻力变化曲线对比图 Fig. 5 Comparison of resistance variation curves under different flow direction angles

计算值与试验值的误差范围如图6所示,误差范围均在3%以内,流向角度为15°,60°以及165°时的误差值最小,仅为0.06%左右;流向角度为30°及45°时误差较大为2.53%和2.35%。总体来看,模拟结果在工程允许的误差范围内,表明本文方法精确有效。

图 6 误差范围 Fig. 6 Error range
3.2 不同来流方向平台表面压力分析

图7为来流速度为1 m/s时不同来流方向下平台表面压力分布图,可见受水流撞击作用的影响,平台迎流面所受压力最大,平台以及立柱和横撑的背流面的压力较小,甚至表现为负压区域。当流向角度增大时,处于下游方向立柱所受阻力越来越小,这是由于随着流向角度增大,上游迎流方向立柱的尾流“屏蔽区”越来越大[10]。当流向角度为90°时,“屏蔽区”达到最大,下游方向结构完全处于右侧浮体的尾流中,其所受压力明显小于上游迎流方向。

图 7 不同流向角度压力分布图 Fig. 7 Pressure distribution at different flow directions and angles
3.3 不同来流方向下的流载荷分析

图8为来流速度为1 m/s时不同流向角度下不同方向阻力变化曲线图。可知,x方向阻力变化规律近似于余弦函数,阻力在流向角度0°~60°内减小速度较缓,在流向角度60°~90°内减小速度较快。

图 8 不同流向角度下阻力变化曲线 Fig. 8 Resistance variation curve under direction angles

y方向阻力变化规律近似正弦函数,阻力在流向角度0°~60°内增长速度较快,在流向角度60°~90°内增长速度较慢,当流向角度为60°时,y方向阻力最大为50.38 N。z方向阻力变化规律关于流向角度180°具有一定的对称性,且均为负压,这是由于平台在一定流向角度上迎流面积相似从而导致阻力变化曲线对称。当流向角度为0°~60°及90°~165°时,阻力随流向角度的增加而增加,当流向角度为60°~90°及165°~180°时,阻力随流向角度的增加而减小。

3.4 不同来流速度下的阻力分析

图9图11为不同来流速度下xyz方向的阻力,图12为不同来流速度下阻力最大值。由图9可知,不同流速下的x方向阻力变化规律均近似于余弦函数,x方向阻力的值域范围随着来流速度的增加而增加。当流向角度为90°时,x方向阻力值均为0。由图10图11可知,不同流速下的y方向、z方向阻力变化规律关于流向角度90°近似对称,流速越大,阻力越大。由图12可知,y方向和z方向的阻力最大值近似相等。

图 9 不同来流速度x方向阻力 Fig. 9 x direction at different inflow speeds

图 11 不同来流速度z方向阻力 Fig. 11 x direction at different inflow speeds

图 12 不同来流速度阻力最大值 Fig. 12 Max resistance at different velocities

图 10 不同来流速度y方向阻力 Fig. 10 ydirection at different inflow speeds
4 结 语

基于STAR-CCM+平台对不同流向角度及不同来流速度下的平台阻力和三维表面压力进行计算,并分析来流速度和来流方向对阻力及表面压力的影响。主要结论如下:

1)本文基于CFD数值模拟方法进行的数值分析有效,对半潜式海洋平台的阻力分析具有一定的参考价值;

2) $ k - \varepsilon $ 湍流模型对半潜式海洋平台阻力的预报偏差在2.6%以内;

3)半潜式平台所受阻力受“遮蔽”效应和迎流面面积的影响;

4)相同来流速度下,各方位阻力随流向角度增大规律各不相同;

5)相同流向角度下,阻力随着来流速度的增大而增大。

参考文献
[1]
MATAUMTO F T, WATAI R A, SIMOA A N. Wave run-up and air gap prediction for a large-volume semi-submersible platform[C]// ASME 2010 29th International Conference on Ocean, Offshore and Arctic Engineering, 2010.
[2]
LI B B, OU J P and TENG B. The full coupled effects of hull, mooring and risers model in time domain based on an innovative deep draft multi-spar[J]. China Ocean Engineering, 2010, 24(2).
[3]
范亚丽, 匡晓峰, 方田, 等. 半潜式支持平台恶劣海况下的运动性能试验研究[J]. 船舶工程, 2021, 43(2): 134-140.
Fan Yali, Kuang Xiaofeng, Fang Tian, et al. Experimental Study on the Motion Performance of Semi-submersible Support Platform in Severe Sea Conditions[J]. Ship engineering, 2021, 43(2): 134-140.
[4]
张明霞, 卢鹏程, 韩兵兵, 等. 基于STAR-CCM+平台和OCIMF方法的FPSO流载荷计算与分析[J]. 中国海洋平台, 2020, 35(4): 49-55.
Zhang Mingxia, Lu Pengcheng, Han Bingbing, et al.. Calculation and analysis of FPSO flow load based on STAR-CCM+ platform and OCIMF method[J]. China Ocean Platform, 2020, 35(4): 49-55.
[5]
董斌, 陈凯旋, 聂焱, 等. 半潜式平台横撑波浪抨击载荷下结构强度分析[J]. 舰船科学技术, 2019, 41(1): 80-84.
Dong Bin, Chen Kaixuan, Nie Yan, etc.. Structural strength analysis of cross brace of semi-submersible platform under wave attack load[J]. Ship Science and Technology, 2019, 41(1): 80-84. DOI:10.3404/j.issn.1672-7649.2019.01.015
[6]
袁洪涛, 陈正豪, 陈刚, 等. 风浪流联合作用下油气资源开发保障平台系泊系统响应的计算分析[J]. 船海工程, 2019, 48(2): 159-162+166.
Yuan Hongtao, Chen Zhenghao, Chen Gang, etc.. Calculation and analysis of mooring system response of oil and gas resources development support platform under wind and wave combined action[J]. ship & Ocean Engineering, 2019, 48(2): 159-162+166. DOI:10.3963/j.issn.1671-7953.2019.02.041
[7]
陈江华, 裴春, 尼鹏, 等. 风浪流联合作用下移动式水上作业平台的施工稳定性[J]. 船舶工程, 2021, 43(7): 147-152.
Chen Jianghua, Pei Chun, Ni Peng, etc.. Construction stability of mobile offshore platform under the combined action of wind and waves[J]. Ship Engineering, 2021, 43(7): 147-152.
[8]
安康. 半潜式海洋平台撑杆在拖航工况下的波浪砰击载荷预报[D]. 镇江: 江苏科技大学, 2020.
[9]
谷家扬, 杨琛, 刘为民, 等. 基于CFD技术的半潜式钻井服务支持平台风荷载预报[J]. 中国海洋平台, 2015, 30(6): 42-48.
Gu Jiayang, Yang Chen, Liu Weimin, et al.. Wind load forecast of semi-submersible drilling service platform based on CFD technology[J]. China Ocean Platform, 2015, 30(6): 42-48. DOI:10.3969/j.issn.1001-4500.2015.06.008
[10]
彭超, 冯光, 郑文涛, 等. 半潜式钻井平台风载荷CFD预报[J]. 船舶工程, 2020, 42(S1): 420-423.