参激与涡激耦合作用下海洋立管振动响应模型建立及预报方法

1.本发明涉及一种海洋立管振动响应模型建立方法和预报方法，属于海洋装备技术领域。


背景技术：

2.海洋立管作为海洋油气资源开发工程中的关键性基础装备，其上端与平台相连，下端与井口或水下生产系统相连，主要承担着传递能量、信息以及输送流体等重要作用。当海洋油气开采工作水深超过300米时，传统的固定式平台(如导管架平台)已不再适用，取而代之的是各式各样的浮式平台(如张力腿平台、半潜式平台、以及spar平台等)。与浅水固定式平台相连的立管相比，与深水浮式平台相连的立管将会面临新的挑战。深水立管除了承受外流绕结构流动引起的漩涡激励作用外，还会承受来自上端浮式平台的边界激励。浮式平台在周期性环境载荷(如波浪)作用下会产生周期性运动(比如垂荡)，这种运动会通过上端悬挂点作为边界激励作用在立管顶端，从而进一步导致立管内部张力发生周期性地变化，这种激励一般称为参数激励。参数激励引起的结构振动称为参激振动，属于非线性振动领域，其共振条件是一系列连续的区间，一旦系统的固有频率与参数激励的频率接近或相等时，极有可能引发剧烈的共振现象。此外，参激振动也将会与外流激励引起的涡激振动互为耦合，从而激发出更为复杂的振动响应特性。为了适应深海油气开采工程的发展需求，有必要建立一种准确可靠地参激与涡激耦合作用下海洋立管的振动响应预测方法。
3.目前绝大多数针对海洋立管振动响应特性的研究，均是将海洋立管轴线上的张力看作是恒定张力(即张力与时间变量以及空间变量均无关)。但是海洋立管在实际工作时，轴向张力不再为恒定张力，而是时变张力，变化部分的张力可由两部分组成：第一部分变化的张力(只与空间无关，而与时间无关)由湿重引起；第二部分变化的张力(只与时间有关，而与空间无关)由上端浮式平台的周期性运动产生。与轴向恒张力作用相比，轴向时变张力作用下，立管结构固有频率不再是恒定值，而是时变的。与此同时，对于真实的海洋来流，从海底到海面，流速变化会随着海流位置越来越大。换言之，实际海洋工程应用中的流体阻尼也是沿立管轴向呈动态变化的。针对这种真实流剖面下参激和涡激耦合作用下深海立管振动响应的研究，需要建立新模型、提出新方法，才能准确预报其振动响应特性，从而为深水立管早期合理设计以及服役期安全工作提供理论基础和技术保障。


技术实现要素：

4.本发明是为了解决目前的海洋立管振动响应特性的研究均是将海洋立管轴线上的张力看作是恒定张力导致海洋立管振动响应存在非常大的误差的问题。
5.参激与涡激耦合作用下海洋立管振动响应模型建立方法，包括以下步骤：
6.将柔性的海洋立管记为柔性立管，柔性立管圆柱体两端采用铰接边界条件，取坐标原点o为圆柱体的下端，x方向为顺流方向，z方向为铅直方向，y方向则是横向振动方向，
x,y以及z三个方向形成右手直角坐标系；
7.建立柔性立管的振动方程：
[0008][0009]
其中，y表示横向振动位移，t表示时间，z表示轴向坐标；ei为柔性立管的弯曲刚度；rs和rf分别表示结构阻尼系数以及流体阻尼系数；m
total
为单位长度的振动系统质量；θ(z,t)为柔性立管上受到的张力；p(z,t)为单位长度立管上受到的升力；
[0010]
所述柔性立管上受到的张力如下：
[0011]
θ(z,t)＝t
top-lwr zwr kacos(ωet)
[0012]
其中，t
top
为立管顶张力；l为立管长度；k为刚度补偿系数；wr为浸没重量，即湿重；
[0013]
所述单位长度立管上受到的升力如下：
[0014][0015]cl
(z,t)＝c
l0
·
q(z,t)/2；
[0016]
其中，c
l
(z,t)为升力系数；u(z)表示为任意位置海洋来流速度；d为立管外径；c
l0
为静止圆柱的升力系数；q(z,t)表示圆柱体z处在y方向的尾流振子的运动；
[0017]
采用改进的van der pol方程来满足尾流振子的非线性特性，尾流振子方程表达式如下：
[0018][0019]
其中，a、ε为经验参数；ωf为局部漩涡脱落频率。
[0020]
进一步地，所述的流体阻尼系数rf＝γωfρd2；
[0021]
其中，ωf为局部漩涡脱落频率，γ为粘滞力系数，ρ立管外部流体密度。
[0022]
进一步地，所述局部漩涡脱落频率ωf＝2πstu(z)/d；
[0023]
其中，u(z)是外部流体的速度沿立管轴向的函数，st为斯脱哈尔数。
[0024]
进一步地，所述粘滞力系数
[0025]
其中，为平均拖曳力系数。
[0026]
进一步地，所述单位长度的振动系统质量如下：
[0027][0028]
其中，ρs、ρf、ρw分别为柔性立管材料密度、立管内部流体密度、立管外部流体密度；d为立管外径，d为立管内径；cm为附加质量系数。
[0029]
进一步地，所述立管顶张力t
top
＝f
top w
r l；
[0030]
其中，f
top
为顶张力系数。
[0031]
进一步地，所述的刚度补偿系数k＝l wr/ac；
[0032]
其中，ac为系统补偿系数所定义临界赋值。
[0033]
进一步地，所述浸没重量
[0034]
其中，g为重力加速度。
[0035]
一种参激与涡激耦合作用下海洋立管振动响应预报方法，包括以下步骤：
[0036]
步骤1、利用所述的参激与涡激耦合作用下海洋立管振动响应模型建立方法建立参激与涡激耦合作用下海洋立管振动响应模型；
[0037]
步骤2、基于4阶runge-kutta方法对参激与涡激耦合作用下海洋立管振动响应模型中的振动方程与尾流振子方程组合进行联立求解，进而得到输出参激与涡激耦合作用下海洋立管振动响应特性。
[0038]
进一步地，步骤2所述的求解过程具体包括以下步骤：
[0039]
将步骤1所得振动方程与尾流振子方程组合，并将张力项展开，将时间二阶导项移至等号左边，其余项均移至等号右边，得到式(8)和式(9)：
[0040][0041][0042]
先对式(8)进行展开：假设柔性圆柱体总长度l可划分为m段，则计算空间步长为δz＝l/m；被离散后的m 1个空间点可表示为：z＝zi，i＝0,1,2,
……
,m；
[0043]
空间偏导数项的二阶精度差分格式表达式如下：
[0044][0045]
将其用矩阵形式表示：
[0046]
[0047][0048][0049]
式(11)至式(13)中出现的系数矩阵a、b、c，均进行了边界条件的处理；
[0050]
令结合四阶rk法展开如下：
[0051][0052][0053][0054]
[0055][0056]
至此，将立管参数、初始条件和边界条件带入式(14)和式(15)中，完成对式(8)的离散；
[0057]
按照相同的方式对式(9)进行离散；
[0058]
将离散后方程耦合，进而得到参激与涡激耦合作用下海洋立管振动响应特性。
[0059]
有益效果：
[0060]
本发明针对不同流剖面下参激与涡激耦合作用下深海立管振动响应特性进行了研究，建立了一套完成的海洋参激-涡激耦合模型，先是基于张力变化因素，给出参激与涡激联合作用下海洋立管的结构振动控制方程。随后基于4阶runge-kutta方法对结构振动控制方程以及尾流振子方程构成的耦合振动方程组进行联立求解。最终实现参激与涡激耦合作用下海洋立管振动响应分析和预报。该模型实时考虑了立管轴向张力变化引起的结构固有频率以及流体阻尼特性，可更为精确地去预报实时变张力作用下海洋立管振动响应特性。
[0061]
本发明针对一真实海洋立管模型，对其在均匀流、线性剪流、以及真实阶梯流作用下结构振动响应特性进行了分析，可以确定本发明能够实现快速预测参激与涡激耦合作用下海洋立管振动响应特性，通过不同工况的对比研究，为工程上进一步探究复杂海洋环境下结构体动力特性提供参考。
附图说明
[0062]
图1为线性剪切来流下两端铰接的细长柔性圆柱体模型示意图；
[0063]
图2为各类流剖面下不同工况下结构体rms值，其中图2(a)为均匀流剖面的结果，图2(b)为线形剪切流剖面的结果，图2(c)为阶梯剪切流剖面的结果。
具体实施方式
[0064]
具体实施方式一：
[0065]
本实施方式为参激与涡激耦合作用下海洋立管振动响应模型建立方法，包括以下步骤：
[0066]
如图1所示，针对深海立管，同时考虑海洋平台的沉浮运动以及海洋来流对立管的冲击作用，建立参激-涡激耦合振动模型，圆柱体两端采用铰接边界条件，取坐标原点o为圆柱体的下端，x方向为顺流方向，z方向为铅直方向，y方向则是横向振动方向，x,y以及z三个方向形成右手直角坐标系。柔性立管上受到的张力为θ(z,t)，由两部分组成：一是由顶张力和考虑湿重引起立管轴向随空间位置变化的张力；二是由平台沉浮引起的随时间周期性变化的张力，柔性立管的弯曲刚度为ei。本模型提供了一种更贴近实际工况的快速预测手段，因此海洋来流剖面可设置为均匀流、线性剪切流以及阶梯剪切流，图1中仅用线性剪切流作为展示。
[0067]
将图1中的细长柔性圆柱体看作细长张力梁模型，建立如下振动方程：
[0068][0069]
式(1)中，y表示横向振动位移，t表示时间，z表示轴向坐标；rs和rf分别表示结构阻尼系数以及流体阻尼系数。这里假设结构阻尼系数rs取为0，流体阻尼系数表示为：rf＝γωfρd2，其中γ为粘滞力系数，ρ立管外部流体密度，ωf为根据斯脱哈尔关系式计算得到的局部漩涡脱落频率，d为立管外径；为平均拖曳力系数，这里取为1.2；ωf＝2πstu(z)/d，u(z)是外部流体的速度沿立管轴向的函数，st为斯脱哈尔数；m
total
为单位长度的振动系统质量，包括三部分，分别为柔性立管结构质量、立管内部流体质量、立管外部流体附加质量；m
total
可表示如下：
[0070][0071]
式(2)中，ρs、ρf、ρw分别为柔性立管材料密度、立管内部流体密度、立管外部流体密度；d为立管内径；cm为附加质量系数，对于圆柱体，cm＝1.0。
[0072]
θ(z,t)表示轴向张力，由两部分组成：一是由顶张力和考虑湿重引起立管轴向随位置变化的张力；二是由平台沉浮引起的随时间变化的张力：
[0073]
θ(z,t)＝t
top-lwr zwr kacos(ωet)(3)
[0074]
其中，t
top
为所考虑的立管顶张力，t
top
＝f
top w
r l，f
top
为顶张力系数取为1.3；l为立管长度；k为刚度补偿系数，k＝l wr/ac，ac为系统补偿系数所定义临界赋值，取为10；wr为浸没重量(湿重)可表示如下：
[0075][0076]
其中，g为重力加速度；
[0077]
式(1)中，p(z,t)为单位长度立管上受到的升力，可表示如下：
[0078][0079]
式(5)中c
l
(z,t)为升力系数，可表示为：c
l
(z,t)＝c
l0
·
q(z,t)/2，其中c
l0
为静止圆柱的升力系数，这里取为0.3；u(z)表示为任意位置海洋来流速度；q(z,t)表示圆柱体z处在y方向的尾流振子的运动。
[0080]
采用改进的van der pol方程来满足尾流振子的非线性特性，尾流振子方程表达式如下：
[0081][0082]
式(6)中a以及ε为经验参数，a＝12；ε＝0.3。
[0083]
具体实施方式二：
[0084]
本实施方式为一种参激与涡激耦合作用下海洋立管振动响应预报方法，包括以下步骤：
[0085]
步骤1、建立参激与涡激耦合作用下海洋立管振动响应模型；
[0086]
步骤2、基于4阶runge-kutta方法对参激与涡激耦合作用下海洋立管振动响应模
型中的振动方程和尾流振子方程进行联立求解，进而得到输出参激与涡激耦合作用下海洋立管振动响应特性。
[0087]
步骤2所述的联立求解的具体过程包括以下步骤：
[0088]
runge-kutta方法是目前最为有效的一种求解微分方程初值问题的数值方法，其中因三阶方法常常不能满足求解需求，五阶方法对精度提升不大，所以常用局部截断误差为o(h5)的四阶方法，简称为四阶rk法。通过运用多个点前一步的值yn来计算下一时间步后的值y
n 1
，因其计算量较大，一般结合计算机进行求解。四阶rk法的常用基本格式如下：
[0089][0090]
将步骤1所得结构振动方程与尾流振子方程组合，并将张力项展开，为适应数值方法，将时间二阶导项移至等号左边，其余项均移至等号右边，得到式(8)和式(9)：
[0091][0092][0093]
因式(8)和式(9)处理思想同根同源，现仅对式(8)进行展开展示。对位置导数部分运用中心差分思想展开，假设柔性圆柱体总长度l可划分为m段，则计算空间步长为δz＝l/m。被离散后的m 1个空间点可表示为：z＝zi(i＝0，1，2，～，m)。空间偏导数项的二阶精度差分格式表达式如下：
[0094][0095]
将其用矩阵形式表示成式(11)，需要注意因首尾两个端点位置所用数据点存在虚拟点，故在程序运行时需通过边界条件进行覆盖，通过不同的边界条件运行出的结果会明显不同，式(11)至式(13)中出现的系数矩阵a、b、c，均进行了边界条件的处理。
[0096][0097][0098][0099]
令结合四阶rk法展开如下，需要注意式(14)中表述为向量形式，即整个立管上m 1个空间点：
[0100][0101]
[0102][0103][0104][0105]
至此，可将立管参数、初始条件和边界条件带入式(14)和式(15)中，完成对式(8)的离散。依照上述思路对式(9)进行离散，将离散后方程耦合，编写成程序语言，通过计算机输出参激与涡激耦合作用下海洋立管振动响应特性。
[0106]
实施例：利用上述方法得到参激与涡激耦合作用下海洋立管振动响应特性，基于分析数据，对实例进行计算分析。
[0107]
(1)基本参数(国际单位)
[0108]
本发明针对一真实海洋立管，具体参数如下，模拟其在不同流剖面(均匀流、线性剪切流以及阶梯剪切流)下参激与涡激耦合作用时海洋立管振动响应。
[0109]
表1立管计算参数
[0110][0111]
(2)计算结果
[0112]
计算结果如图2所示，图2为各类流剖面下不同工况下结构体rms值，其中图2(a)为均匀流剖面的结果，图2(b)为线形剪切流剖面的结果，图2(c)为阶梯剪切流剖面的结果。需要注意，图2(a)至图2(c)对应的各类流剖面中(b)～(e)工况顺序分别为f
top-a-ωe(顶张力系数-变张力幅值-变张力频率)，图2(a)为均匀流剖面，即整个轴向上流速均相等；图2(b)为剪切流剖面，顶部流速最大，底部流速为零，轴向上流速呈线性分布；图2(c)为阶梯流剖面，流剖面系数见表2。
[0113]
表2真实阶梯流剖面归一化系数
[0114][0115][0116]
从图中可以看出，均匀流剖面下，顶张力系数增大则立管振动位移均方根值有收
缩趋势，模态阶数明显降低，随着变张力幅值和频率增加，立管均方根值略微增加，总体变化不大。线性剪切流下，立管中上部也会出现振动位移峰值，同时可以发现立管均方根值对顶张力系数更为敏感。阶梯流剖面下，由于整个轴向上流速变化斜率不尽相同，可以发现立管对顶张力系数、变张力幅值以及变张力频率都比较敏感，轻微的增加都会引起均方根值的改变。对比三种流剖面，相同组合工况下均方根分布差异较大，说明立管振动响应对流速较为敏感，在预测真实海洋立管时应首先考虑模拟流速和真实流速的相似度。
[0117]
通过应用上述方法，可以实现快速预测参激与涡激耦合作用下海洋立管振动响应特性，通过不同工况的对比研究，为工程上进一步探究复杂海洋环境下结构体动力特性提供参考，提高海洋资源开发能力，降低海洋生态污染。
[0118]
本发明还可有其它多种实施例，在不背离本发明精神及其实质的情况下，本领域技术人员当可根据本发明作出各种相应的改变和变形，但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。