船舶力学

名称	数值	名称	数值
质量/kg	2.12205E7	甲板宽度/m	54.72
浮筒长度/m	80.6	立柱直径/m	12.9
浮筒宽度/m	16	吃水/m	13
甲板长度/m	54.72	排水量/m³	20786.65

名称	数值	名称	数值
质量/kg	2.12205E7	甲板宽度/m	54.72
浮筒长度/m	80.6	立柱直径/m	12.9
浮筒宽度/m	16	吃水/m	13
甲板长度/m	54.72	排水量/m³	20786.65

名称	数值	名称	数值
锚链（M1-M6）的长度/m	851.5	立管（R4-R7）的长度/m	490
锚链（M7-M12）的长度/m	1176.5	立管（R18）的长度/m	469
立管（R1）的长度/m	504	立管（R19）的长度/m	469
立管（R2）的长度/m	490	立管（R20）的长度/m	553.6
立管（R3、R8-R17）的长度/m	490	立管（R21-R22）的长度/m	480

名称	数值	名称	数值
锚链（M1-M6）的长度/m	851.5	立管（R4-R7）的长度/m	490
锚链（M7-M12）的长度/m	1176.5	立管（R18）的长度/m	469
立管（R1）的长度/m	504	立管（R19）的长度/m	469
立管（R2）的长度/m	490	立管（R20）的长度/m	553.6
立管（R3、R8-R17）的长度/m	490	立管（R21-R22）的长度/m	480

参数	数值	参数	数值
弹性模量E/GPa	1.0	海冰与结构间的摩擦系数μ_s	0.25
海水密度ρ_w/（kg·m^-3）	1035.0	颗粒法向粘结强度σ_bⁿ/MPa	0.6
海冰密度ρ_i/（kg·m^-3）	920.0	颗粒切向粘结强度σ_b^s/MPa	0.6
颗粒间的摩擦系数μ_b	0.25	颗粒抗压强度σ_c/MPa	2.53

参数	数值	参数	数值
弹性模量E/GPa	1.0	海冰与结构间的摩擦系数μ_s	0.25
海水密度ρ_w/（kg·m^-3）	1035.0	颗粒法向粘结强度σ_bⁿ/MPa	0.6
海冰密度ρ_i/（kg·m^-3）	920.0	颗粒切向粘结强度σ_b^s/MPa	0.6
颗粒间的摩擦系数μ_b	0.25	颗粒抗压强度σ_c/MPa	2.53

类型	网格总数
模型1	139 92
模型2	244 16
模型3	385 80

类型	网格总数
模型1	139 92
模型2	244 16
模型3	385 80

参数	年平均值
空气温度	最高	6 °C
最低	-35 °C
波浪	有义波高	7 m
谱峰周期	11.9 s
水流	近水面最大速度	1 m/s
水温	水面年最大值	7.5 °C
海冰（平整冰）	冰厚	1.5~2.0 m
冰速	0.2~0.4 m/s

参数	年平均值
空气温度	最高	6 °C
最低	-35 °C
波浪	有义波高	7 m
谱峰周期	11.9 s
水流	近水面最大速度	1 m/s
水温	水面年最大值	7.5 °C
海冰（平整冰）	冰厚	1.5~2.0 m
冰速	0.2~0.4 m/s

类型	波谱类型	有义波高/m	谱峰周期/s	类型	波谱类型	有义波高/m	谱峰周期/s
工况1	Jonswap	7	11.9	工况3	Jonswap	11	14.7
工况2	Jonswap	9	13.6	工况4	Jonswap	13	15.8

类型	波谱类型	有义波高/m	谱峰周期/s	类型	波谱类型	有义波高/m	谱峰周期/s
工况1	Jonswap	7	11.9	工况3	Jonswap	11	14.7
工况2	Jonswap	9	13.6	工况4	Jonswap	13	15.8

类型	海冰厚度/m	海冰速度/（m·s^-1）	类型	海冰厚度/m	海冰速度/（m·s^-1）
工况5	1.5	0.2	工况9	1.6	0.1
工况6	1.6	0.2	工况10	1.6	0.15
工况7	1.8	0.2	工况11	1.6	0.25
工况8	2.0	0.2	工况12	1.6	0.3

类型	海冰厚度/m	海冰速度/（m·s^-1）	类型	海冰厚度/m	海冰速度/（m·s^-1）
工况5	1.5	0.2	工况9	1.6	0.1
工况6	1.6	0.2	工况10	1.6	0.15
工况7	1.8	0.2	工况11	1.6	0.25
工况8	2.0	0.2	工况12	1.6	0.3

类型	海冰速度/（m·s^-1）	海冰厚度/m	海冰密集度	冰块平均面积/m²
工况13	1	2	80%	20
工况14	1	2	80%	50
工况15	1	2	80%	100
工况16	1	1.8	80%	50
工况17	1	1.8	60%	50
工况18	1	1.8	40%	50
工况19	1	1.8	20%	50
工况20	1	1.5	60%	50
工况21	1	2	60%	50

类型	海冰速度/（m·s^-1）	海冰厚度/m	海冰密集度	冰块平均面积/m²
工况13	1	2	80%	20
工况14	1	2	80%	50
工况15	1	2	80%	100
工况16	1	1.8	80%	50
工况17	1	1.8	60%	50
工况18	1	1.8	40%	50
工况19	1	1.8	20%	50
工况20	1	1.5	60%	50
工况21	1	2	60%	50

极地半潜式海洋平台在不同环境载荷下的动态响应分析

PDF下载

张奥博 ¹ , 闯振菊 ¹ , 刘社文 ¹ , 李春郑 ² , 季顺迎 ³ , 屈衍 ⁴

船舶力学 | 结构力学 2024,28(3): 422-434

收起

船舶力学 | 结构力学 2024, 28(3): 422-434

极地半潜式海洋平台在不同环境载荷下的动态响应分析

全屏

张奥博¹, 闯振菊¹, 刘社文¹, 李春郑², 季顺迎³, 屈衍⁴

作者信息

^1.大连海事大学　船舶与海洋工程学院，辽宁　大连　116026

^2.上海勘探设计研究院有限公司，上海　200335

^3.大连理工大学　工业装备结构分析国家重点实验室，辽宁　大连　116024

^4.华南理工大学　海洋科学与工程学院，广州　511442

张奥博（1996-），男，硕士研究生

闯振菊（1983-），女，博士，副教授，通讯作者，E-mail：zhenjuchuang@dlmu.edu.cn。

通讯作者:

通讯作者，E-mail：zhenjuchuang@dlmu.edu.cn

Dynamic response analysis of arctic semi-submersible offshore platform under different environmental loads

Ao-bo ZHANG¹, Zhen-ju CHUANG¹, She-wen LIU¹, Chun-zheng LI², Shun-ying JI³, Yan QU⁴

Affiliations

^1.Naval Architecture and Ocean Engineering College, Dalian Maritime University, Dalian 116026, China

^2.Shanghai Investigation, Design & Research Institute Co., Ltd., Shanghai 200335, China

^3.State Key Laboratory of Structure Analysis of Industrial Equipment, Dalian University of Technology, Dalian 116024, China

^4.School of Marine Science and Engineering, South China University of Technology, Guangzhou 511442, China

出版时间: 2024-03-20 doi: 10.3969/j.issn.1007-7294.2024.03.010

文章导航

摘要

收起

北极地区拥有丰富的油气资源，在该海域进行油气资源的开发和利用具有巨大的发展前景，但也极具挑战性。当前海洋平台的设计规范并未充分考虑海冰的影响。北极海域存在无冰期和有冰期，为了实现平台全年在海上进行作业，需要联合考虑波浪载荷和冰载荷的作用。本文研究某半潜式海洋平台在波浪、平整冰以及碎冰不同工况下的动态响应，采用势流理论和莫里森方程计算波浪载荷，利用离散元方法求解冰载荷，使用细长体有限元的方法分析锚泊系统，通过非线性时域计算求解得出平台的运动以及锚泊系统受力。结果显示：相较于其他环境载荷，平整冰作用下锚链的张力最大；在较小的海冰厚度、海冰密集度和冰块平均面积下的碎冰区，平台的运动符合现行安全规范。本文的研究可为北极地区油气资源的安全开发提供参考价值。

关键词

北极油气资源 / 半潜式海洋平台 / 波浪载荷 / 冰载荷 / 冰与结构相互作用

Abstract

收起

The Arctic region is rich in oil and gas resources. The exploitation and utilization of oil and gas resources in this area have great prospects, but are also very challenging. Arctic areas have ice-free and ice-covered seasons. As current design specifications for offshore platforms do not fully consider the impact of sea ice, the effect of wave loads and ice loads should be considered together to enable the platform to operate all year round. In this paper, the dynamic response of a semi-submersible offshore platform under different conditions was studied, such as wave, level ice, and ice floes. Potential flow theory and Morrison equation were used to calculate wave loads while the discrete element method was used to solve the ice loads. The slender body finite element method was used for analysis of the slender system. The motion of the platform and the force of the slender system were solved by nonlinear time-domain calculation. The results show that compared with other environmental loads, the tension of the mooring line is the highest under level ice. In a broken ice field with smaller ice thickness, ice concentration, and ice floe size, the platform motion conforms to current safety specifications. The research of this paper will provide a reference for the safe development of oil and gas resources in the Arctic region.

Key words

Arctic oil and gas resource / semi-submersible offshore platform / wave load / ice load / ice-structure interaction

引用本文

张奥博, 闯振菊, 刘社文, 李春郑, 季顺迎, 屈衍. 极地半潜式海洋平台在不同环境载荷下的动态响应分析. 船舶力学, 2024 , 28 (3) : 422 -434 . DOI: 10.3969/j.issn.1007-7294.2024.03.010

Ao-bo ZHANG, Zhen-ju CHUANG, She-wen LIU, Chun-zheng LI, Shun-ying JI, Yan QU. Dynamic response analysis of arctic semi-submersible offshore platform under different environmental loads[J]. Journal of Ship Mechanics, 2024 , 28 (3) : 422 -434 . DOI: 10.3969/j.issn.1007-7294.2024.03.010

正文

收起

0　引言

收起

全球变暖的现象加剧了北极地区的冰雪融化，也为北极油气资源的开采创造了条件。2008年，美国地质调查局完成了对北极圈以北所有地区未发现的常规油气资源的评估。报告推断出北极未发现的常规油气资源平均储量约为900亿桶石油、1669亿万立方英尺天然气和440亿桶液化天然气。其中，大约84%未被发现的油气资源位于海上^[1]。然而，在北极地区进行石油勘探是非常困难的，难点在于大量未发现的石油储存在水深大于100米^[2]且小于500米^[3]的水域，再加上海冰对结构物的影响，极地油气资源的开采并未成为常态。在极地进行海上作业的平台需要考虑作业水深、可变载荷、抗冰能力、恶劣环境的适应能力等多种因素^[4]。半潜式海洋平台因为其较好的运动性能和能够适应深水作业的特点被看作是在极地开采油气资源的最佳方案。

当前学术界对波浪载荷作用下海洋平台的研究取得了很多成果。刘海霞和肖熙^[5]得出了对于小结构物，波浪的拖曳力和惯性力是主要分量，对于大结构物，波浪的惯性力和绕射力是主要分量。大型半潜式平台的波浪载荷计算采用Morison公式和势流理论相结合的方法；张威等^[6]分别采用频域分析法与时域分析法对南海某深水半潜式平台的运动响应进行了计算，对比分析了两种方法在深水半潜式平台运动响应预报结果上的异同；宋安科^[7]对第五代半潜式钻井平台进行了水动力分析，并对系泊系统和平台进行了静力分析和时域耦合动力分析，得到了平台运动位移和系泊线张力的时历响应曲线。

国内外众多学者对海冰与平台的相互作用已开展了深入的分析。挪威科技大学Hansen和Loset^[8]采用二维圆盘单元对破碎冰场在外力作用下的运动学和动力学进行了模拟，并对定位在碎冰区的海洋结构物进行了分析；Zhou等^[9]建立了一种二维数值模型来模拟平整冰和锚泊结构相互作用的过程，应用船舶冰阻力的半经验公式模拟了冰的旋转和滑动过程，数值模拟的结果通过与Wright^[10]的实测数据进行对比得到了验证；刘璐等^[11]基于闵可夫斯基原理的扩展多面体离散元方法模拟船舶与海洋平台结构的相互作用过程，采用离散元方法分析了船舶结构绕行过程中的船舶和海洋平台结构冰载荷。

一些学者对海洋平台的抗冰性能也进行了研究。张帅^[12]通过数值模拟计算了圆柱及两种抗冰形状结构的波浪载荷，讨论了这两种结构在无冰期对海洋平台桩腿的影响；龙雪等^[13]采用具有粘结破碎功能的离散元方法计算海冰与锥体作用的破坏过程，将离散元计算冰载荷及海冰破坏模式与渤海现场实测数据进行对比验证，得出倒锥体结构可有效降低冰载荷从而具有较好的抗冰性能；张大勇等^[14]对典型半潜式平台与海冰相互作用机理和冰载荷模型进行分析，提出了该类结构在海冰作用下的主要失效模式及评价方法。

海洋平台在极地进行海上作业的过程中需要进行海冰管理。本研究的主要目的是研究半潜式生产平台在北极有冰和无冰时的动态响应，包括海冰管理前的平整冰工况以及海冰管理后的碎冰工况。通过时域模拟，研究平台在波浪荷载和冰荷载作用下的运动、锚链张力的变化规律，并对不同工况的影响加以分析。

1　数值模型介绍

收起

本文选用的半潜式海洋平台如图1所示，主要由包括钻井设备和生活区等大多数功能模块的上部结构、提供浮力的浮筒、四根连接浮筒和上层甲板的立柱以及立柱间的横撑组成。平台的更多信息可参照表1。

数值模型中的锚泊系统包括12根锚链及22根立管，与半潜式海洋平台构成一个耦合的整体。平台的运动会影响锚链的受力。朝着不同方向呈放射型布置的锚链可以提供给平台回复力和力矩，抵抗外部载荷以保证平台可以正常作业。锚泊系统的具体布置及外部载荷的入射角度如图2所示（M表示锚链，R表示立管），锚链和生产立管的具体情况见表2。

2　理论方法介绍

收起

2.1　波浪荷载分析

本文采用势流理论^[15]和Morison理论^[16]来计算分析波浪荷载。平台的水动力特性分析以及结构运动响应计算是基于SESAM软件^[17]来完成的。WADAM模块^[18]可以用来进行波浪载荷下的频域响应分析，通过计算得到一阶波浪力和二阶波浪力的传递函数。SIMA模块^[19-20]可以用来进行平台与锚泊系统耦合下的时域响应分析，通过计算不同工况下平台的运动以及锚泊系统的受力，得到整体的运动响应特性。

在研究半潜式海洋平台的运动时，忽略它的柔性变形，将其看作刚体，主要分析平台六自由度的运动情况。OXY平面与静止的水面重合，OZ轴垂直于水面，指向上方。

根据势流理论，可以得到半潜式海洋平台在频域下的一阶运动方程为

(1)

式中，M_ij为质量矩阵，μ_ij为附加质量矩阵，λ_ij为阻尼系数矩阵，C_ij为回复力系数矩阵，X_j为平台的运动，F_i为平台所受到的一阶波浪力。

Morison方程主要分为两个部分，分别是拖曳力计算和惯性力计算。Morison方程的表达式为

(2)

式中，F_D表示拖曳力，F_I表示惯性力。

线性的一阶波浪力是通过将速度势分解为入射势、绕射势和辐射势，利用伯努利方程求得。二阶波浪力对系泊系统的运动有着重要的影响，因为平台在波浪的作用下会产生漂移。平均漂移力是一个非线性的二阶力，计算时采用远场法^[21]和压力积分法^[22]。文中选择的半潜式海洋平台是双对称结构，波浪方向选为0°~90°，间隔为15°。传递函数计算周期为4~38 s，步长为2 s。

2.2　冰荷载分析

本文采用离散元^[23]的方法来分析冰荷载。该方法可以模拟冰从完整状态到破裂状态的整个过程^[24]，从而分析海冰的破坏模式以及冰载荷的大小。在离散元的方法中，海冰颗粒间的相互作用采用了基于莫尔-库仑剪切摩擦定律的弹性-粘性接触模型^[25]来确定。该方法采用了平行粘结模型将球形颗粒结合成平整冰。平行粘结模型中的两个粘结颗粒单元间设定一个弹性粘结圆盘，这样的模型既能传递力又能传递力矩^[26]。

平行粘结模型的基本情况如图3所示。在平行粘结模型中，

和

分别为颗粒间的法向力和切向力，

和

分别为颗粒间的法向力矩和切向力矩，作用在粘结圆盘上的最大正应力和最大剪应力^[24]可分别表示为

(3)

(4)

式中，R为粘结圆盘的半径、A为粘结圆盘的横截面积、I和J分别为粘结圆盘的惯性矩与极惯性矩。具体计算可表示为

(5)

判断平行粘结圆盘的破坏失效采用了断裂准则。在模拟海冰的断裂过程中，可以通过颗粒单元的法向和切向受力状态将单元之间的粘结失效方式分为拉伸破坏和剪切破坏^[27]，如图4所示。当粘结圆盘的最大正应力σ_max超过颗粒单元之间的拉伸破坏强度σ^t，或者最大剪应力τ_max超过颗粒单元间的剪切破坏强度τ^s时，粘结圆盘会失效。

碎冰区的构建是基于二维Voronoi图的切割算法来完成的。采用二维Voronoi图可以将碎冰区域快速划分为不规则的多边形集合。在此基础上，按照一定比例调整各个多边形的尺寸以生成碎冰区的平面几何模型，由此通过切割球体单元构造的平整冰即可得到碎冰块的离散元模型^[28]。此方法可以对不同海冰厚度、海冰密集度和冰块平均面积下的碎冰区进行模拟，从而进行离散元计算。

用离散元方法模拟海冰所需要的主要参数^[29]如表3所示。

本文将综合比较半潜式海洋平台在波浪载荷和冰载荷（包括平整冰和碎冰）下六个自由度的运动和锚链的受力，图5给出了数值研究的流程。

2.3　模型的网格收敛性分析

网格的收敛性分析是为了验证不同的网格数目对计算结果产生的影响。本文选用了三组不同的模型进行频域下水动力性能的计算和分析。这三组不同的模型是基于水下平均湿表面在纵向和周向显示不同的网格来创建的，三组网格模型的具体信息如表4所示。

本文通过对比三组网格模型的平台运动幅值响应算子RAO、平台受到的一阶波浪力以及二阶平均漂移力来验证网格的收敛性。图6~7为4~38 s的计算周期内平台纵荡、纵摇的一阶运动传递函数；图8~9为平台受到的一阶波浪力；图10~12为在远场法和压力积分法的计算下得到的平台二阶平均漂移力。

从图中可以看出，三组模型的一阶波浪力和二阶平均漂移力几乎一致，平台在纵荡、垂荡以及纵摇方向上的一阶运动幅值响应算子RAO略有偏差。模型3在存在偏差处的数据更接近三组模型数据的平均值，而且此模型的网格最精细。考虑到计算精度，本文最终选定模型3（如图13）作为后续的网格模型。

3　数值计算与结果分析

收起

3.1　计算工况介绍

本文选定了坐标在北纬60°~70°、西经45°~60°的格陵兰海域作为研究对象。格陵兰海域的气象条件、海洋情况以及海冰条件都可以在ISO 19906（2010）规范^[30]中找到。具体数据如表5所示。

本文需要分别对比波浪载荷和冰载荷作用下锚链的拉力变化以及海洋平台的运动响应。不同环境载荷下平台和锚泊系统的相关数据都能在SIMA中计算得到。根据格陵兰海域的环境情况，最终确定4组波浪作用下的工况、8组平整冰作用下的工况以及9组碎冰作用下的工况，具体情况如表6~8所示。波浪和海冰的入射方向均为180°。碎冰的速度与该海域水的流速一致，确定为1 m/s，计算碎冰时不考虑波浪的影响。波浪载荷的不同工况是通过有义波高和谱峰周期加以区分的；平整冰的影响因素包括冰厚和冰速，计算区域为2000 m×200 m；碎冰的变量分别为海冰厚度、海冰密集度和冰块平均面积，计算区域为4000 m×200 m。图14~19展示了不同密集度和平均面积下的碎冰场。

3.2　环境载荷对平台系统的的影响分析

图20展示了在不同波浪载荷的作用下，5号锚链受到拉力的最大值、最小值和平均值。选取5号锚链是因为在180°外部载荷的作用下，5号锚链始终处于张紧状态，比较具有代表性。从图中可知，随着有义波高和谱峰周期的增加，锚链受到的最大拉力也随之变大。工况4下锚链受到的最大拉力比工况1、工况2、工况3分别增大了55.6%、39.0%、13.2%。

讨论平整冰的影响因素可以参考图21，从图中可以看出不同冰厚和冰速下产生的差异，海冰厚度的影响相较于海冰速度更明显。冰的厚度越大，锚链受到的拉力就越大。而冰速变化对锚链的受力影响不大。这也与Wright^[31]与Zhou等^[32]的研究成果相符，即参照Kulluk的经验，冰速与冰载荷大小之间并不存在相关性；冰速较低时，锚链拉力的最大值受冰速的影响很明显，而拉力的平均值与之关系不大。

碎冰不同参数的作用效果如图22~24所示。图22中冰块的平均面积作为变量，在冰厚和密集度相同的条件下，平均面积越大，锚链受到的拉力就越大。从图中可以看出，在平均面积为100 m²的碎冰工况下锚链受到的拉力远大于在平均面积为20 m²的碎冰工况下的受力。图23中碎冰场的密集度作为变量，在冰厚和冰块平均面积相同的条件下，密集度越大，锚链受到的拉力就越大。从图中可以看出不同密集度下锚链拉力的差异非常明显。碎冰场的密集度越来越大，冰块可运动的空间就越小。冰块会在结构附近大规模堆积，会引发冰载荷的增大。图24中海冰厚度作为变量，在碎冰区海冰密集度和冰块平均面积相同的条件下，海冰越厚，锚链受到的拉力就越大。因为冰块的质量会随着冰厚的增大而线性增大，质量会影响动能进而影响在冰与结构物之间的作用力^[33]。从图中可以看出，本文选取的三组不同厚度海冰的作用效果差别不大但符合规律。

3.3　波浪载荷和冰载荷下的动态响应对比

本文的重点是为了验证半潜式海洋平台在结冰和开放水域两种不同的环境下是否都能进行正常的海上作业，需要对比在波浪载荷和冰载荷作用下的锚链受力情况以及平台的运动情况，分析时采用了时域比较法和统计学两种方法。

通过折线图来观察时域分析下锚链的张力。时域分析的优点在于可以直观了解到一定时间范围内力的变化情况（本文选用的是3600 s内的数据），样本值更多，作对比分析时就更具有说服力。本文对比了波浪载荷作用下的工况1、工况4，平整冰作用下的工况6、工况12以及碎冰作用下的工况15、工况18。具体情况如图25所示。从图中可以看出，平整冰作用下锚链受到的拉力远大于波浪载荷下受到的拉力。碎冰工况下的整体趋势会略高于波浪下的工况，从中体现出碎冰作用下锚链受力的平均水平略大于波浪作用下的受力。

图26为利用统计学进行分析的波浪、平整冰和碎冰工况下锚链最大拉力的对比图。三类不同的工况在柱状图中分别用三种不同的颜色加以区分。从图中可以看出，波浪和碎冰对锚链产生的最大拉力在同一水平线上，都远远小于平整冰的影响。

沿180°方向的载荷对平台的纵荡和纵摇影响最大，本文着重对比不同载荷的纵荡运动和纵摇运动。图27~30展示了同一入射方向上波浪、平整冰、碎冰三种不同环境下平台的运动情况。其中，图27为时域下纵荡运动的对比情况，从图中可以很明显地看出，在平整冰作用下平台的纵荡运动远大于波浪载荷下的运动，甚至相差数倍；碎冰作用下的纵荡运动介于平整冰与波浪之间。规范对深水区作业的半潜式平台的运动响应有具体要求。在正常的工况下，平台的纵荡运动一般不超过工作海域水深的10%，在恶劣环境的工况下，对平台运动响应的要求没那么严格，只需确保平台极限水平偏移量不超过水深的15%即可。本文计算的水深为330 m。图28为统计学散点图，红线下方浅蓝色的点为符合规范的工况。从图中可知，波浪载荷作用下平台最大的纵荡运动可以控制在10%水深的范围内，而在平整冰的作用下，纵荡运动已经远远超过了水深的15%，甚至超过了水深的35%。对于碎冰来说，合适的海冰厚度、海冰密集度以及冰块平均面积会使得平台的运动能够满足现有规范的要求。在本文中，工况13、工况18、工况19以及工况20就是能保证正常海上作业的碎冰区。

图29为时域下纵摇运动的对比情况。与纵荡的结论类似，在平整冰作用下平台的纵摇运动要大于波浪载荷下的运动，碎冰作用下的纵摇运动介于二者之间。对于深水区的半潜式平台，摇角一般要控制在10°以内。如图30所示，波浪载荷下平台的纵摇运动满足此规定，但是平整冰作用下的纵摇角度在某些工况下达到了25°，严重影响正常的工程作业。文中选择的碎冰工况，其作用下的纵摇运动均符合当前规范的要求。

4　结论

收起

本文研究了无冰与有冰两种不同情况下半潜式海洋平台在极地海上作业时的动态响应。无冰期时，平台在波浪载荷作用下进行生产作业，在波浪方向相同的情况下，有义波高和谱峰周期是需要考虑的影响因素。结冰期时，平台会受到海冰的撞击。对于平整冰来说，冰厚相较于冰速对平台运动和锚链张力的影响更大。为了保证海洋平台能够进行正常的海上作业，海冰管理是必不可少的一环，因此对碎冰的研究非常重要。研究碎冰时选择的变量包括海冰厚度、碎冰区的密集度以及冰块的平均面积。通过分析波浪载荷和冰载荷作用下平台的运动和锚链的张力，可以得到以下结论：

（1）对碎冰来说，在冰厚和密集度相同的条件下，冰块平均面积越大，锚链的张力就越大，平台的纵荡和纵摇也越明显；在冰厚和冰块平均面积相同的前提下，随着碎冰区密集度的增大，锚链张力和平台运动也会相应增大。在碎冰区的海冰密集度和冰块平均面积都相同的情况下，海冰越厚，对平台和锚泊系统产生的作用效果就更剧烈。归根结底，这三个变量都与海冰的质量息息相关。海冰的质量会影响海冰的动量，海冰动量的增大会最大程度影响海冰与结构物之间的作用力。

（2）结合时域分析和统计学两种方法可知，锚链在平整冰作用下的张力远大于波浪载荷和碎冰下的张力。在平整冰的作用下平台的纵荡运动和纵摇运动也远大于波浪载荷和碎冰下平台的运动。在平整冰的工况下进行海上作业危险性极大，而海冰管理后产生的碎冰区也无法完全保证平台的安全。碎冰对平台纵荡运动的影响大于纵摇运动，因为在有些碎冰工况下，平台的纵摇运动满足规范的要求，而纵荡运动则超过规范规定的最大偏移。故在极地进行油气开采时除了海冰管理的要求外，还要尽可能地减小冰载荷对平台纵荡运动的影响。

基金

收起

国家自然科学基金项目(52301311)
中央高校基本科研业务费专项资金资助(3132023516)
辽宁省教育厅基本科研项目(LJKMZ20220365)

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

[1]

Bird

K J

, Charpentier

R R

, Gautier

D L

, et al. Circum-Arctic resource appraisal: Estimates of undiscovered oil and gas north of the Arctic circle[J]. Geological Survey (US), 2008.

[2]

Hamilton

J M

. The challenges of deep water Arctic development[C]//The Twenty-first International Offshore and Polar Engineering Conference, International Society of Offshore and Polar Engineers, 2011.

[3]

Donald

L G

, Kenneth

J B

, Ronnald

R C

, et al. Assessment of undiscovered oil and gas in the Arctic[J]. Science, 2009, 324 (5931): 1175-1179.

[4]

刘大辉. 半潜式钻井平台抗冰性能多维优化设计及评估方法研究[D]. 杭州：浙江大学，2020.

Liu

Dahui

. Research on multi-dimensional optimization design and evaluation method of anti-ice performance of semi-submersible drilling platform[D]. Hangzhou: Zhejiang University, 2020.

[5]

刘海霞，肖熙.半潜式平台结构强度分析中的波浪载荷计算[J]. 中国海洋平台，2003(02)：3-6.

Liu

Haixia

, Xiao

. The wave load calculation in the analysis of structural strength of a semi-submersible[J]. China Offshore Platform, 2003(02): 3-6.

[6]

张威，冯玮，谢彬.深水半潜式平台运动响应预报方法的对比分析[J]. 中国海上油气，2009，21(01): 61-64+68.

Zhang

Wei

, Feng

Wei

, Xie

Bin

. Comparative analysis on forecasting procedures of deepwater semi-submersible platform motion-response[J]. China Offshore Oil and Gas, 2009, 21(01): 61-64+68.

[7]

宋安科.深水半潜式钻井平台系泊系统方案设计与分析[D]. 哈尔滨：哈尔滨工程大学，2008.

Song

Anke

. Design and analysis of mooring system of deepwater drilling semi-submersible[D]. Harbin: Harbin Engineering University, 2008.

[8]

Hansen

E H

, Loset

. Modelling floating offshore units moored in broken ice: Model description[J]. Cold Regions Science and Technology, 1999, 29(2): 97-106.

[9]

Zhou

, Su

Biao

, Kaj

Riska

, Torgeir

Moan

. Numerical simulation of moored structure station keeping in level ice[J]. Cold Regions Science and Technology, 2012, 71: 54-66.

[10]

Wright

. Full scale experience with Kulluk station-keeping operations in pack ice (with reference to Grand Banks Developments)[C]//Working Report 25-44 for the Program of Energy, Research and Development (PERD). National Research Council, Ottawa, Canada, 2000.

[11]

刘璐，尹振宇，季顺迎.船舶与海洋平台结构冰载荷的高性能扩展多面体离散元方法[J]. 力学学报，2019，51(06): 1720-1739.

Liu

, Yin

Zhenyu

, Ji

Shunying

. High performance extended polyhedron discrete element method for ice loads on ships and offshore platforms[J]. Chinese Journal of Theoretical and Applied Mechanics, 2019, 51(06): 1720-1739.

[12]

张帅.冰与直立型及锥型结构的相互作用力分析[D]. 哈尔滨：哈尔滨工程大学，2013.

Zhang

Shuai

. Analysis of interaction force between ice and vertical and cone structures[D]. Harbin: Harbin Engineering University, 2013.

[13]

龙雪，刘社文，季顺迎.水位变化对正倒锥体冰载荷影响的离散元分析[J]. 力学学报，2019，51(01): 74-84.

Long

Xue

, Liu

Shewen

, Ji

Shunying

. Influence of water level on ice load on upward-downward conical structure based on DEM analyss[J]. Chinese Journal of Theoretical and Applied Mechanics, 2019, 51(01): 74-84.

[14]

张大勇，于东玮，王国军，等. 半潜式海洋平台抗冰性能分析[J]. 船舶力学，2020，24(2): 208-220.

Zhang

Dayong

, Yu

Dongwei

, Wang

Guojun

, et al. Analysis of ice resistance of semi-submersible offshore platform[J]. Journal of Ship Mechanics, 2020, 24(2): 208-220.

[15]

Faltinsen

O M

, Minsaas

K J

, Liapis

, Skjordal

S O

. Prediction of resistance and propulsion of a ship in a seaway[C]//Proceedings of the 13th Symposium on Naval Hydrodynamics, Tokyo, 1980: 505-523.

[16]

Morison

J R

, Johnson

J W

, Schaff

S A

, et al. The force exerted by surface waves on piles[J]. Journal of Petroleum Technology, 1950, 2(5): 149-154.

[17]

Maruo

. The drift of a body floating on waves[J]. Journal of Ship Research, 1960, 4(3): 1-10.

[18]

DNV. SESAM User manual Sesam manager[R]. 2018.

[19]

DNV. SESAM User manual WADAM[R]. DNV Software Report No.: 92-7052, 2017.

[20]

SINTEF Ocean. RIFLEX 4.16.0 User guide[R]. 2019.

[21]

SINTEF Ocean. SIMO 4.16.0 User guide[R]. 2019.

[22]

Pinkster

J A

, Van Oortmerssen

. Computation of the first and second-order wave forces on oscillating bodies in regular waves[C]//Proceedings of the 2nd International Conference on Numerical Ship Hydrodynamics, University of California, Berkeley, 1977: 136-159.

[23]

Shen

H H

, Hibler

W D

, Lepparanta

. On applying granular flow theory to a deforming broken ice field[J]. Acta Mechanica, 1986, 63(1-4): 143-160.

[24]

Lau

, Lawrence

K P

, Rothenburg

. Discrete element analysis of ice load on ships and structures[J]. Ships and Offshore Structure, 2011, 6(3): 211-221.

[25]

Potyondy

D O

, Cundall

P A

. A bonded-particle model for rock[J]. International Journal of Rock Mechanics and Mining Sciences, 2004, 41: 1329-1364.

[26]

龙雪. 海洋结构物作用下海冰破坏模式及冰载荷的离散元分析[D]. 大连：大连理工大学，2019.

Long

Xue

. Discrete element analysis of sea ice failure mode and ice load under ocean structures[D]. Dalian: Dalian University of Technology, 2019.

[27]

Campbell

C S

. Stress-controlled elastic granular shear flows[J]. Journal of Fluid Mechanics, 2005, 539: 273-297.

[28]

朱红日，季顺迎，刘璐. 基于Voronoi切割算法的碎冰区构造及离散元分析[J]. 计算力学学报，2019，36(04): 454-463.

Zhu

Hongri

, Ji

Shunying

, Liu

. Construction of broken ice field with Voronoi tessellation algorithm and its DEM simulations[J]. Chinese Journal of Computational Mechanics, 2019, 36(04): 454-463.

[29]

孔帅，季顺迎，季少鹏，等. 基于高性能离散元方法的极区浮式平台冰载荷数值分析[J]. 中国舰船研究，2021，16(05): 64-70.

Kong

Shuai

, Ji

Shunying

, Ji

Shaopeng

, et al. Numerical analysis of ice load on floating platform in polar region based on high-performance discrete element method[J]. Chinese Journal of Ship Research, 2021, 16(05): 64-70.

[30]

ISO 19906:2010, Petroleum and natural gas industries-Arctic offshore structures[S]. 2010.

[31]

Wright

. Evaluation of full scale data for moored vessel station-keeping in pack ice[R]. PERD/CHC Report, 1999.

[32]

Zhou

, Su

, Riska

, et al. Numerical simulation of moored ship in level ice[C]//Proceedings of the International Conference on Offshore Mechanics and Arctic Engineering (OMAE). 2011.

[33]

刘璐. 扩展多面体离散元方法及其在海洋结构冰载荷分析中的应用[D]. 大连：大连理工大学，2019.

Liu

. Extended polyhedral discrete element method and its application to ice load analysis of marine structures[D]. Dalian: Dalian University of Technology, 2019.

2024年第28卷第3期

PDF下载

引用本文

BibTeX

文章信息

doi: 10.3969/j.issn.1007-7294.2024.03.010

接收时间：2023-09-27
首发时间：2026-03-21
出版时间：2024-03-20

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2023-09-27

基金

国家自然科学基金项目(52301311)

中央高校基本科研业务费专项资金资助(3132023516)

辽宁省教育厅基本科研项目(LJKMZ20220365)

作者信息

^1.大连海事大学　船舶与海洋工程学院，辽宁　大连　116026

^2.上海勘探设计研究院有限公司，上海　200335

^3.大连理工大学　工业装备结构分析国家重点实验室，辽宁　大连　116024

^4.华南理工大学　海洋科学与工程学院，广州　511442

通讯作者:

通讯作者，E-mail：zhenjuchuang@dlmu.edu.cn

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/cblx/CN/10.3969/j.issn.1007-7294.2024.03.010

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT

名称	数值	名称	数值
质量/kg	2.12205E7	甲板宽度/m	54.72
浮筒长度/m	80.6	立柱直径/m	12.9
浮筒宽度/m	16	吃水/m	13
甲板长度/m	54.72	排水量/m³	20786.65

名称

数值

名称

数值

质量/kg

2.12205E7

甲板宽度/m

54.72

浮筒长度/m

80.6

立柱直径/m

12.9

浮筒宽度/m

吃水/m

甲板长度/m

54.72

排水量/m³

20786.65

名称	数值	名称	数值
锚链（M1-M6）的长度/m	851.5	立管（R4-R7）的长度/m	490
锚链（M7-M12）的长度/m	1176.5	立管（R18）的长度/m	469
立管（R1）的长度/m	504	立管（R19）的长度/m	469
立管（R2）的长度/m	490	立管（R20）的长度/m	553.6
立管（R3、R8-R17）的长度/m	490	立管（R21-R22）的长度/m	480

名称

数值

名称

数值

锚链（M1-M6）的长度/m

851.5

立管（R4-R7）的长度/m

490

锚链（M7-M12）的长度/m

1176.5

立管（R18）的长度/m

469

立管（R1）的长度/m

504

立管（R19）的长度/m

469

立管（R2）的长度/m

490

立管（R20）的长度/m

553.6

立管（R3、R8-R17）的长度/m

490

立管（R21-R22）的长度/m

480

参数	数值	参数	数值
弹性模量E/GPa	1.0	海冰与结构间的摩擦系数μ_s	0.25
海水密度ρ_w/（kg·m^-3）	1035.0	颗粒法向粘结强度σ_bⁿ/MPa	0.6
海冰密度ρ_i/（kg·m^-3）	920.0	颗粒切向粘结强度σ_b^s/MPa	0.6
颗粒间的摩擦系数μ_b	0.25	颗粒抗压强度σ_c/MPa	2.53

参数

数值

参数

数值

弹性模量E/GPa

1.0

海冰与结构间的摩擦系数μ_s

0.25

海水密度ρ_w/（kg·m^-3）

1035.0

颗粒法向粘结强度σ_bⁿ/MPa

0.6

海冰密度ρ_i/（kg·m^-3）

920.0

颗粒切向粘结强度σ_b^s/MPa

0.6

颗粒间的摩擦系数μ_b

0.25

颗粒抗压强度σ_c/MPa

2.53

类型	网格总数
模型1	139 92
模型2	244 16
模型3	385 80

类型

网格总数

模型1

139 92

模型2

244 16

模型3

385 80

参数	年平均值
空气温度	最高	6 °C
最低	-35 °C
波浪	有义波高	7 m
谱峰周期	11.9 s
水流	近水面最大速度	1 m/s
水温	水面年最大值	7.5 °C
海冰（平整冰）	冰厚	1.5~2.0 m
冰速	0.2~0.4 m/s

参数

年平均值

空气温度

最高

6 °C

最低

-35 °C

波浪

有义波高

7 m

谱峰周期

11.9 s

水流

近水面最大速度

1 m/s

水温

水面年最大值

7.5 °C

海冰（平整冰）

冰厚

1.5~2.0 m

冰速

0.2~0.4 m/s

类型	波谱类型	有义波高/m	谱峰周期/s	类型	波谱类型	有义波高/m	谱峰周期/s
工况1	Jonswap	7	11.9	工况3	Jonswap	11	14.7
工况2	Jonswap	9	13.6	工况4	Jonswap	13	15.8

类型

波谱类型

有义波高/m

谱峰周期/s

类型

波谱类型

有义波高/m

谱峰周期/s

工况1

Jonswap

11.9

工况3

Jonswap

14.7

工况2

Jonswap

13.6

工况4

Jonswap

15.8

类型	海冰厚度/m	海冰速度/（m·s^-1）	类型	海冰厚度/m	海冰速度/（m·s^-1）
工况5	1.5	0.2	工况9	1.6	0.1
工况6	1.6	0.2	工况10	1.6	0.15
工况7	1.8	0.2	工况11	1.6	0.25
工况8	2.0	0.2	工况12	1.6	0.3

类型

海冰厚度/m

海冰速度/（m·s^-1）

类型

海冰厚度/m

海冰速度/（m·s^-1）

工况5

1.5

0.2

工况9

1.6

0.1

工况6

1.6

0.2

工况10

1.6

0.15

工况7

1.8

0.2

工况11

1.6

0.25

工况8

2.0

0.2

工况12

1.6

0.3

类型	海冰速度/（m·s^-1）	海冰厚度/m	海冰密集度	冰块平均面积/m²
工况13	1	2	80%	20
工况14	1	2	80%	50
工况15	1	2	80%	100
工况16	1	1.8	80%	50
工况17	1	1.8	60%	50
工况18	1	1.8	40%	50
工况19	1	1.8	20%	50
工况20	1	1.5	60%	50
工况21	1	2	60%	50

类型

海冰速度/（m·s^-1）

海冰厚度/m

海冰密集度

冰块平均面积/m²

工况13

80%

工况14

80%

工况15

80%

100

工况16

1.8

80%

工况17

1.8

60%

工况18

1.8

40%

工况19

1.8

20%

工况20

1.5

60%

工况21

60%