Chinese Journal of Ship Research
铺管船运动响应及对铺设管线的耦合作用
1，李磊2，张兆德2，郭远志2王玉红
1 316021浙江大学 海洋学院，浙江 舟山2 316022浙江海洋大学 船舶与机电工程学院，浙江 舟山
摘 要：［目的］随着深海海洋资源的不断开发，铺管船在海底管道铺设任务中扮演着非常重要的角色，故有必101”号铺管船作业能力提升实际工程项目为要开展浅水铺管船提升作业能力的分析工作。［方法］以“中油海SESAM Orcaflex依托，运用 和 软件分析铺管船在增加作业水深时其铺设管道与船体及系泊系统的相互作用，包括管道和管道水下悬垂段对船体运动响应和锚链张力的影响、船体六自由度运动对管道张应力和弯曲应力的影响。［结果］结果表明：管道对船体运动和锚链张力的影响较大，管道水下悬垂段的长度对船体纵荡运动和锚链张力的影响较大，铺管船垂荡运动对管道受力的影响最大。［结论］研究成果可为浅水铺管船作业能力提升项目提供理论参考和评估依据。关键词：铺管船；系泊系统；管道；耦合作用
中图分类号：U674.34
文献标志码：A
DOI：10.19693/j.issn.1673-3185. 01096
0引言
在实施海洋强国战略的过程中，不仅需要研发先进的技术和装备，同时为了实现资源的合理配置，也需要对老旧设备进行升级改造以提升其作业能力。铺管船在海底管道铺设任务中扮演着非常重要的角色，目前我国具备深水铺管作业能201”，而“109”、力的铺管船仅有“海洋石油 滨海101”、“CPP601”等均为浅水铺管“蓝疆”、“中油海船［1-3］。随着深海海洋资源的不断开发，对大型深海起重铺管船的需求也在日益增加［4］ ，而浅水铺管船则面临着被淘汰的危机，因此有必要开展浅水铺管船提升作业能力的分析工作。101”铺管船的设计作业水深只由于“中油海40m有 ，为了保障西气东输台湾海峡项目，需进一步提升其铺管作业水深，为此很多研究人员对其ANSYS-AQWA进行了分析探讨。曾骥等［ 5 ］运用软件分析了不同航道水深条件下铺管船的运动响应和系泊系统的载荷特性，完成了香港支线铺设40m作业的安全性评估；付翔［6］对比分析了 作业水深条件下同时移动单锚和双锚的锚链受力情况，验证了移动双锚方案的可行性；李磊等［7］研究60m了在 作业水深条件下铺管船提升锚泊定位的能力，分析了最大作业水深和最优铺管方案；高101”铺等［8］针对“中油海兆鑫 管船的托管架优化设计开展了分析。101”铺管船的研究综上所述，目前“中油海大多集中于船体运动响应和锚链载荷特性等方面，未针对其铺设管道与船体及系泊系统的相互101”铺管作用开展研究。因此，本文以“中油海船作业能力提升工程项目为基础，拟采用水动力SESAM Orcaflex，基于势流理论，在时域范软件 和围内分析铺管船增加作业水深时其管道与船体及系泊系统的耦合作用，并探讨船体六自由度运动101”铺管对管道受力的影响，从而完成“中油海船作业能力提升项目的安全评估工作，用以为项目实施提供理论参考和依据。
1 理论基础
2本文的计算过程主要分为 个步骤：一是利SESAM用 软件计算铺管船系泊系统与铺设管道Orcaflex之间的耦合作用；二是利用 软件计算铺管船运动对铺设管道受力的影响。本文采用连续介质法对系泊缆受力进行动力分析，假设 r ，ṙ ，r̈ 分别为单元的位移、速度和加速度，则其离散有限元动力方程为［9］ RI(r R RS(r RE(r （1） r t) + D(r r t) + t) = r t)式中： RI(r r t)为惯性力矢量，其中 t为时间； R D(r r t) 为阻尼力矢量； RS(r t) 为内部单元之间的作用力；RE(r r t)为外部作用力，主要包括重力、浮力、系缆受到的浪流作用力和船体运动引起的相关作用力。船体在时域范围内的耦合运动方程为［10-11］ t
式中：M 和 μ分别为质量矩阵与附加质量矩阵； ẍ(t) 为船体加速度；ẋ(t - τ) 为船体速度，其中 τ 为
积分的时间变量；K (τ) 为延时函数矩阵；C 为回
复力矩阵； x(t) 为船体位移； F (1) ，F (2)分别为wave wave一阶波浪力和二阶波浪力；F为系泊力；F moor current
为流力；F 为风力。wind OrcaFlex在 软件中，管道采用悬链线模型并将其视为柔性缆，其中若干段模型之间采用质点1链接，其有限元模型如图 所示。图中：Sx1 ，Sy1 ， 3 Sz 分别为柔性缆分解的 个方向，其中 Sz 为柔性3个缆的轴向方向；Nx，Ny，Nz分别为质点分解的方向，其中Nz为质点的轴向方向；α为 Sz 与Nz方向的夹角。 可以参考式（2）来求解管道与船体的耦合运动方程，采用集中质量法计算管道受到的轴向拉Orcaflex压、横向弯曲和扭转等作用力［12］。根据 软件说明书，管道单元所受的有效张力为
（3） Te = T + (P A - P A i) w o o i式中：Te 为有效张力；T 为管壁张力；P 和 P 分w i o别为管道内、外部压力； A 和 A 分别为管道的i o
内、外受力截面面积。
2 仿真模型 2.1 模型参数
SESAM-GeinE首先，在 中建立铺管船的湿表HydroD面模型；然后，导入 中计算得到船体的水DeepC动力参数；最后，导入 中进行船体与系泊系统的时域耦合分析，并以此计算结果为基础开展铺管船与管道的耦合作用研究。101”的设计参数如表1 “中油海 所示，湿表面2模型如图 所示。
2.2 环境工况
2 13］，#1～#8铺管船的作业海况如表 所示［ 锚3 3链的布置方式如图 所示，锚链和管道参数如表所示，Orcaflex 4计算模型如图 所示。
3 计算结果
3.1 铺设管道对船体运动响应和锚链张力的影响
0°，45°和 90°，管道假设波浪入射方向分别为 4和对船体运动响应及锚链张力的影响分别如表5表 所示。4 0°时，船体纵荡由表 可知，当波浪入射角为幅值的变化较大，其他入射角条件下的幅值变化较小，这是因为在铺管过程中，管道从船艉延伸至海底，其自重分力会对铺管船产生较大的纵向拉0°时船体的纵荡幅值变化较力，故波浪入射角为大；在斜浪及横浪工况下，船体的横荡、艏摇运动幅值变化较大，而垂荡幅值的变化则较小；当波浪
90°时，横摇、纵摇幅值的变化较大。由入射角为此可见，管道对船体的运动响应有一定的影响且对纵荡运动的影响较大，该结论与参考文献［14-15］的计算结果相符。5 0°时，由于船艉由表 可知，当波浪入射角为的管道对船体存在拉力作用，故#1～#4锚链的最大张力值有所增加，而#5～#8锚链的最大张力值45°时，#1～#4则有所减小；当波浪入射角为 锚链的最大张力值存在不同程度的增加，而#5～#8锚链的最大张力值则均呈减小趋势；当波浪入射角90°时，#3锚链的张力值增幅最大，虽然#6为 锚链和#3锚链相对于船体布置，但由于船艉管道的存上#6锚链的张力变化值小于#3在，实际 锚链。综上可知，管道对锚链张力的影响较大，该结论与参考文献［14］的计算结果相符。
3.2 管道水下悬垂段对船体运动响应和锚链张力的影响
2铺管船在作业期间铺设的管线分为 段，即上弓段和悬垂段。从管线进入导管架到管线自身的反弯点位置即为上弓段，而悬垂段则是指反弯点到海底接触点的管线长度。假设管道水下悬垂103，105，107，109，111 m，本节将段的长度分别为分析水下管道不同悬垂段对船体运动和锚链张力5 6的影响，结果如图 和图 所示。
5由图 可知，随着悬垂段长度的增加，纵荡运动的增幅最为明显，垂荡、纵摇、艏摇运动的增幅较为平缓，而横摇运动的幅值则是减小后维持不变。因此，除了纵荡运动之外，悬垂段长度的变化5对其他 个自由度运动的影响均较小。6可知，随着悬垂段长度的增加，#5～#8由图锚链的张力增加，而#1～#4锚链的张力则呈减小#5 #8趋势，且变化幅值相对较大。由于 和 锚链、#6和#7锚链、#1和#3锚链、#2和#4锚链相对于船体对称布置，故其均随悬垂段长度的增加而呈现出相同的变化趋势。可见，管道水下悬垂段的长度对铺管船体纵荡运动和锚链张力的影响较大。
3.3 船体运动对管道受力的影响
本节将分析船体各个自由度运动对铺设管道受力的影响，以研究管道受力对船体运动的敏感程度。在管道铺设过程中，船艉处附近管道上端的张应力最大，而海底段附近的管道因弯曲变形6 Orcaflex导致的弯曲应力最大。表 所示为输入7软件的铺管船运动幅值，图 所示为管道上端的轴向张应力随船体各个自由度运动的变化曲线， 8图 所示为管道的最大弯曲应力随船体各个自由度运动的变化曲线。7由图 可知，船体六自由度运动均会对铺设管道上端的管道受力产生一定的影响，且管道受力随着各个自由度运动幅值的增加而增加。从计算结果分析得出，横摇、艏摇、纵摇运动对管道受力的影响较小；纵荡运动产生的管道张应力最大
值略小于横荡，而垂荡运动产生的管道张应力值 最大。综上所述，船体平移运动对管道张应力的
影响要远大于摇荡运动，船体六自由度运动对管道张应力的影响依次为：艏摇<横摇<纵摇<纵荡<横荡<垂荡，该结论与参考文献［14，16］的计算结果一致。8由图 可知，船体六自由度运动均会对管道底部的弯曲应力产生一定的影响，其中横摇、艏摇
运动对管道弯曲应力的影响较小；纵摇运动产生 的管道弯曲应力最大值小于纵荡和横荡运动；垂 荡运动产生的管道弯曲应力值最大。同时，平移 运动的幅值不同，管道的最大受力位置也有所不 同；随着平移运动幅值的增加，管道的最大受力位 置将会逐渐向上端转移。综上所述，平移运动对
管道弯曲应力的影响要远大于摇荡运动，船体六自由度运动对管道弯曲应力的影响依次为：艏摇<横摇<纵摇<横荡<纵荡<垂荡，该结论与与参考文献［14，16］的计算结果相符。
4结论
101”铺本文以“中油海 管船提升作业能力实SESAM际工程项目为依托，基于势流理论，运用Orcaflex和 软件分析了铺管船在增加作业水深时管道与船体及系泊系统的相互作用，得到如下结
论： 1）管道对船体运动和锚链张力的影响较大， 0°入射角时船体且受波浪入射方向的影响，其中纵荡运动的幅值较大，90°入射角时锚链张力的幅值较大。2 ）管道悬垂段的长度对船体纵荡运动的影响较大，对其他自由度运动的影响较小；悬垂段长
度对锚链张力的影响较大，当悬垂段长度增加时， #5～#8锚链的张力将随之增加，而#1～#4锚链的张力则会减小。3 ）船体平移运动对管道受力的影响要远大于摇荡运动，其中垂荡运动的影响最大。船体六自由度运动对管道张应力的影响依次为：艏摇<
横摇<纵摇<纵荡<横荡<垂荡；对管道弯曲应力的影响依次为：艏摇<横摇<纵摇<横荡<纵荡<垂荡。
参考文献：
1］ 黄维平，曹静，张恩勇. ［ 国外深水铺管方法与铺管J］. 2011，29船研究现状及发展趋势［ 海洋工程， （1）：135-142. HUANG W P，CAO J，ZHANG E Y. State of the art and developing trend of deepwater pipe-laying and pipe-laying vessel abroad［J］. The Ocean Engineer⁃ ing，2011，29（1）：135-142（in Chinese）. 2］ 支洪彬. ——CPP601 船［N］. ［ 海洋利器 铺管 石油管道报，2016-11-03（4）. 3］ 方晟. 铺管船现状概况及发展趋势分析［J］. ［ 无线互联科技，2017（8）：105-107. FANG S. Analysis on present situation and develop⁃ ment trend of pipe laying vessel［J］. Wireless Internet Technology，2017（8）：105-107（in Chinese）. 4］程维杰. 大型起重铺管船船型开发［J］. 船舶，2015 ［ （3）：28-34. CHENG W J. Development of hull line for huge crane and pipe-laying vessels［J］. Ship and Boat，2015（3）： 28-34（in Chinese）. 5］ 曾骥，李陈峰，任慧龙，等. ［ 铺管船铺管作业锚泊系
统分析［J］. 船舶工程，2013，35（5）：24-27. ZENG J， LI C F REN H L ，et al. Analysis on moor⁃
， ing system of pipe-laying vessel under pip-laying［J］. Ship Engineering，2013，35 （5）： 24-27 （in Chi⁃ nese）. 6］ 付翔.“中油海101”铺管船锚泊作业分析［］ J . ［ 天津航海，2014（3）：4-5.
7］ 李磊，吴迪，付春丽，等. ［ 浅水铺管船提升锚泊定位究［C］//2015年深海能源大会论文集.能力的研 海口：中国造船工程学会，2015.
8］ 高兆鑫，付春丽，唐继蔚，等.“中油海101”铺［ 管船J］. 2015，43（8）：托管架优化设计［ 石油机械， 93-97. GAOZX FUCL ，TANG J W，et al. Design optimiza⁃
， tion of the stinger of “ZHONG YOU HAI 101” pipe-laying vessel［J］. China Petroleum Machinery， 2015，43（8）：93-97（in Chinese）. 9］ 刘树晓，唐友刚，李焱，等. ［ 内转塔单点系泊刚度对
析［J］. 中国舰船研究，2015，不同参数的敏感性分10（4）：86-93. LIU S X ，TANG Y G， LIY ，et al. Sensitivity analysis of different parameters for the stiffness of internal tur⁃ ret mooring systems［J］. Chinese Journal of Ship Re⁃ search，2015，10（4）：86-93（in Chinese）. ［ 10］ 操戈，李旭，张咏鸥，等. FPSO液舱晃荡与船舶时域10（1）：88-96.耦合运动数值模拟［J］. 中国舰船研究， 2015， CAO G， LIX ，ZHANG Y O，et al. Numerical simula⁃ tion of ship motion coupled with tank sloshing for FP⁃ SO［J］. Chinese Journal of Ship Research，2015，10 （1）：88-96（in Chinese）. ［ 11］ 刘旭平，唐友刚，张少洋. 船艉开槽型安装船及其运动特性数值分析［ J］. 中国舰船研究， 2015，10 （1）：97-102，108. LIU X P ，TANG Y G，ZHANG S Y. Numerical study on the motion characteristics of a new slotted stern marine engineering installation vessels［J］. Chinese Journal of Ship Research，2015，10 （1）：97-102， 108（in Chinese）. ［ 12］ 徐普，龚顺风，代胜泽，等. 铺管船运动对深水S型铺设管道动力性能影响分析［C］//第十八届中国
海洋（岸）工程学术讨论会（上）. 舟山：中国海洋学
会海洋工程分会，2017. ［ 13］ 陈红霞，华锋，袁业立. 中国近海及临近海域海浪的季节特征及其时间变化［ J］. 海洋科学进展， 2006，24（4）：407-415. CHEN H X，HUA F，YUAN Y L. Seasonal character⁃ istics and temporal variations of ocean wave in the Chinese offshore waters and adjacent sea areas［J］. Advances in Marine Science， 2006， 24 （4 ）： 407-415（in Chinese）. ［ 14］ 庞兆坤. 考虑管道影响的铺管船及系泊系统动力响应研究［ D ］ . 镇江：江苏科技大学，2012. PANG Z K. Study on dynamic response of the pipe-lay and mooring system considering the pipeline ［D］. Zhenjiang：Jiangsu University of Science and Technology，2012（in Chinese）. ［］ 15 宋林峰，孙丽萍，王德军. 深水S型铺管托管架-船体-管线耦合分析［ J］. 哈尔滨工程大学学报， 2013，34（4）：415-420，433. SONG L F， SUN L P ，WANG D J. Coupling analysis of stinger-lay barge-pipeline of S-lay installation in deep wate［r J］. Journal of Harbin Engineering Univer⁃ sity，2013，34（4）：415-420，433（in Chinese）. ［ 16］ 李金玉.起重铺管船水动力性能研究［D］. 上海：上海交通大学，2010. LI J Y. Study on pipelay derrick vessel hydrodynamics ［D］. Shanghai： Shanghai Jiao Tong University， 2010（in Chinese）.