一种海底多跨管道涡激振动的预报方法
【技术领域】
[0001] 本发明涉及一种海底管道的研究方法，更具体地说，本发明涉及一种海底多跨管 道涡激振动预报方法。 技术背景
[0002] 随着我国对渤海、黄海、东海和南海等海洋油气资源的开采逐步进展，海底管道得 到了更广泛的应用。海底管道能够高效快速的将海底油井产出的石油和天然气等能源物质 运送到陆地上，其具有安装方便、造价低廉、高效运输等优点，已成为海底油气输送系统的 大动脉。
[0003] 由于海底地势起伏不平和海底洋流的冲刷作用，容易使海底管道出现悬跨段。当 外界来流流经悬跨管道所产生的漩涡脱落频率与管道结构的固有频率接近而达到锁定状 态时，悬跨管道就会发生涡激振动，从而出现较大的响应幅值使管道遭到严重的疲劳破坏。 管道一旦出现疲劳破坏，管内的油气等流体就会泄漏，以至造成严重的环境污染和生态破 坏。
[0004] 由于悬跨管道的两端由海底土壤支撑，而海底土壤的性质又决定着悬跨管道两端 端部的约束条件，进而决定着悬跨管道的固有频率，并影响着悬跨管道发生涡激振动。一般 的，海底土壤越坚硬，约束条件就越强烈。相应的悬跨管道固有频率越大，涡激振动就越难 发生。
[0005] 国内外很多专家学者对单跨管道涡激振动开展了大量的模型试验和理论研究，并 推出了VIVANA和Orcaflex等工程计算软件。而在实际工况中，当海底管道出现悬跨段时， 在其相邻的不远处往往还会出现其它悬跨段，这就涉及到海底管道出现两跨及两跨以上相 连的现象，即海底多跨管道。
[0006] 海底多跨管道在外界来流的作用下也会发生涡激振动，而且相邻悬跨段之间的悬 跨管道涡激振动存在着较强的相互影响。虽然现有技术能够对海底单跨管道的涡激振动进 行预报，但是对海底多跨管道的涡激振动还无法进行有效的预报。
【发明内容】
[0007] 本发明的目的就是克服现有技术的缺陷，并为此提供一种海底多跨管道涡激振动 的预报方法。该方法考虑了海底土壤对悬跨管道的一般边界约束条件，能够对海底多跨管 道的涡激振动进行有效的预报，为海底多跨管道涡激振动的研究、防范和治理提供了可靠 的理论依据和分析方法。
[0008] 本发明的技术方案是：
[0009] 一种海底多跨管道涡激振动的预报方法，包括以下步骤：
[0010] (1)建立如下管道结构和外界流场的流-固耦合作用模型：
[0013] 式中：x-管道轴向位移；
[0014]t-时间；
[0015] y(X，t) 一在X位置t时刻管道横向振动位移；
[0016] q(X，t) 一在X位置t时刻无量纲涡激升力系数；
[0017]Φi(X)-管道结构第i阶模态，在工程应用中可取前η阶模态；
[0018]Pl(t) -t时刻第i阶模态对应的管道广义坐标；
[0019]djt) -t时刻第i阶模态对应的无量纲涡激升力广义坐标；
[0020] m-质量项，包括管道结构质量，管内流体质量和管外流场附加质量ma，ma = CaπpD2/4,Ca为附加质量系数；
[0021] P-外界流体密度或海水密度；
[0022] D-管道直径；
[0023] c一阻尼项，包括结构阻尼cs和水动力阻尼cw，结构阻尼cs= 2πιωηζ，ζ为结构 阻尼比，水动力阻尼cw=CDPDV/2,CD为拖曳力系数；
[0024]ω;一管道第i阶固有圆频率；
[0025]V-外界来流速度；
[0026] -管道静止时的涡激升力系数；
[0027]ε-尾流振子模型参数；
[0028]Α-尾流振子模型参数；
[0029]ωs-漩涡脱落频率，ωs= 2πStV/D，St为斯特罗哈数；
[0030] 1 一海底多跨管道总长度；
[0031] (2)确定模态和尾流振子模型参数：
[0032]1)将海底管道简化为欧拉-伯努利梁模型，根据管道一般边界约束条件对结构固 有丰吴态进彳丁确定：
[0033]A.管道结构的固有模态：
[0036]EI-管道结构的弯曲刚度；
[0037]T一管道结构受到的轴向拉力；
[0038]B.海底多跨管道每一跨取一固有模态，其系数Cl]、c2]、〇3]和c4]需根据约束条件 确定，其中j表示第j跨，j= 1、2、3…；
[0040] 式中：kr一扭转弹簧系数；
[0041]kt 一拉伸弹簧系数；
[0043] 式中：y-连接点左侧位移；
[0044]y+-连接点右侧位移；
[0046] 根据管道的约束条件得到关于Cl]、c2]、C3，c4](j= 1、2、3···)的方程组，求解此 方程组系数的秩等于零的超越方程，可得到管道结构固有频率，将固有频率回代到方程组 中，求得模态系数(^、(^、(^和c4j(j= 1、2、3···)，进而可以得到海底多跨管道的各阶模态 和固有频率；
[0047] 2)根据试验数据对尾流振子模型中的参数ε和A进行标定：
[0048]当Vr>5,A= 10,当 0 <Vr〈5,A= 4;
[0053] 经1)和2)确定海底多跨管道涡激振动预报方法中流-固耦合作用模型的模态和 尾流振子参数，对流-固耦合作用微分方程组在时域上进行求解，得到海底多跨管道涡激 振动的时间历程响应，可以采用傅里叶变换将时域响应转换为频域响应。
[0054] 在以上步骤（2)中：
[0055] 模态的系数，根据管道的一般边界约束条件确定；
[0056] 尾流振子模型参数，根据试验数据确定。
[0057] 本发明采用扭转弹簧和拉伸弹簧模拟海底土壤对多跨管道的约束条件，并根据欧 拉-伯努利梁理论和经试验数据标定的尾流振子模型，提出一种能够对海底多跨管道涡激 振动进行有效预报的新方法。该方法克服了现有技术的诸多不足，其有益效果是：
[0058] (1)可以通过改变扭转弹簧或者拉伸弹簧的弹性系数模拟多种海底土壤对多跨管 道的约束条件，弹性系数取值越大，代表海底土壤越坚硬，对管道的约束条件越剧烈；
[0059] (2)利用试验数据对尾流振子模型中的参数进行标定，从而获得精确预报海底多 跨管道涡激振动的技术效果；
[0060] (3)可以得到海底多跨管道涡激振动响应的时间历程以及海底多跨管道涡激振动 响应的频率，为评估海底多跨管道涡激振动的疲劳寿命提供了可靠的理论依据；
[0061] (4)由于以上（1)、（2)、（3)的有益效果，可以为海底多跨管道涡激振动的研究、防 范和治理提供可靠的依据和有效的方法。
【附图说明】
[0062] 附图为本发明的方法流程图。
【具体实施方式】
[0063] 为了使本发明的优点和特征更容易被清楚理解，下面结合附图和实施例对其技术 方案作以详细说明。
[0064] 如前所述，海底地势高低起伏不平和海底洋流的冲刷作用，往往会使海底管道出 现两跨及两跨以上的悬跨现象，称为"多跨管道"。现有技术对海底悬跨管道涡激振动分析 的方法是仅截取其中的一个悬跨段，即以单跨管道为基础进行分析，却忽略了海底管道是 一个连续的整体以及相邻悬跨段之间的相互影响。而且现有技术尚缺乏对海底多跨管道涡 激振动进行有效预报的方法。
[0065] 本发明则考虑了海底土壤对海底多跨管道的一般边界约束条件，提出了一种基于 下述模型的海底多跨管道涡激振动预报方法。
[0066] 参照附图，该方法包括以下步骤：
[0067] (1)建立如下管道结构和外界流场的流-固耦合作用模型：
[0070] 式中：x-管道轴向位移；
[0071] t-时间；
[0072] y(X，t) 一在X位置t时刻管道横向振动位移；
[0073]q(X，t) 一在X位置t时刻无量纲涡激升力系数；
[0074]Φi(x) -管道结构第i阶模态，在工程应用中可取前η阶模态；
[0075]Pl(t)- t时刻第i阶模态对应的管道广义坐标；
[0076]山(t)一t时刻第i阶模态对应的无量纲涡激升力广义坐标；
[0077]m-质量项，包括管道结构质量，管内流体质量和管外流场附加质量ma，ma = CaπpD2/4,Ca为附加质量系数；
[0078] P-外界流体密度或海水密度；
[0079] D-管道直径；
[0080]c一阻尼项，包括结构阻尼cs和水动力阻尼cw，结构阻尼cs= 2πιωηζ，ζ为结构 阻尼比，水动力阻尼cw=CDPDV/2,CD为拖曳力系数；
[0081]ω;-管道结构第i阶固有圆频率；
[0082]V-外界来流速度；
[0083] CL。一管道静止时的涡激升力系数；
[0084]ε-尾流振子模型参数；
[0085]Α-尾流振子模型参数；
[0086]ωs-漩涡脱落频率，ωs= 2πStV/D，St为斯特