1.中海石油深海开发有限公司，广东深圳5180002.哈尔滨工程大学船舶与海洋工程学院，黑龙江哈尔滨150001
摘要：为提高张力腿平台和立管与筋腱之间的动态耦合计算结果的准确性，更好地模拟平台的筋腱系统，运用势流理论,依据南海油气田开发既有资料，使用SESAM软件对一座传统型TLP平台张力筋腱进行分析。考虑筋腱和立管的实际工程设计，采用分段构建系泊系统的方法建立分析模型，并得到了频域下的附加质量系数和系统刚度矩阵。通过静水偏移模拟得到了平台系统水平及垂向刚度曲线和张力筋顶部张力，校验了系统质量和刚度的准确性。通过自由衰减模拟，验证平台固有周期，并采用规则波时域分析的方法求解了张力筋腱张力幅值响应算子(RAO)。结果表明，分段建模能很好地模拟实际工程中的张力腿平台(TLP)平台刚度，系统水平刚度与垂向沉降校验和自由衰减模拟可以更好地保证计算模型的正确性。
关键词：张力腿平台；静水偏移；自由衰减；沉降；固有周期；筋腱张力；水平刚度；幅值响应算子
张力腿平台(TLP)是一种典型的顺应式平台，通过数条张力筋腱与海底相连[1]。TLP的六自由度运动可以简单分为两类，一类主要由系泊系统刚度控制，包括垂直面内的垂荡、横摇和纵摇；另外一类主要由结构的几何形式和系泊张力的水平分量决定，包括水平面内的纵荡、横荡和艏摇运动[2]。由于浮力和重力差筋腱产生了较大的刚度，使得TLP垂直面上运动的固有周期比较小，水平面上运动的固有周期一般比较大[3]，其垂荡运动的固有周期一般为2～4s内，水平方向的系统刚度只是预张力的小角度正选分量，纵横荡运动的固有周期一般在100～200s内，其特有的结构形式较大程度地避免了波浪激励载荷共振的发生[4]。
张力腿平台在环境力作用下有着较大的位移，TLP系统的非线性恢复刚度在其动态运动中起着非常重要的作用。纵荡、横荡、垂荡相对于自身位移是一个非线性函数[5]。通常在进行TLP数值模拟时把平台被视为刚体，张力筋腱模拟为准静态弹簧[6]，计算出的平台刚度矩阵便近似等效为时域分析中的系统刚度。但在实际结构中，系泊的张力筋腱是多成分的复合结构。平台在发生水平位移时导致平台垂向沉降发生，水平自由度和垂直向自由度耦合导致系统结构的刚度随空间位置产生非线性变化。而且，平台的水平运动导致吃水变化，系统环境回复刚度也发生了变化。其低频运动的相互耦合，结构在水平面上运动大变形导致了结构刚度矩阵和阻尼矩阵产生非线性变换[7-8]，从而导致TLP在波浪中产生非线性运动。因而确保模型线性刚度矩阵的准确性是进行张力腿平台动态分析的首要工作，对平台刚度进行系统验证也成为确保计算分析准确的工作之一。
文中使用SESAM软件，对一座传统结构形式TLP的张力筋腱响应进行了研究。考虑平台的结构细节，使用Genie模块建立平台下浮体湿表面模型以及莫里森杆元模型，使用HydroD建模并计算了系统的质量矩阵、刚度矩阵和附加质量系数。考虑实际张力筋腱及立管结构分段，在DeepC中建立全耦合模型，进行静水刚度验证模拟和自由衰减模拟，并通过时域分析方法求解张力筋腱张力响应RAO。
文中用于计算的传统张力腿平台由平台本体(平台上体、立柱和下体浮箱)、张力筋腱系泊系统以及立管系统三大部分组成，其主体结构由4根尺寸相同的圆柱形立柱与4个矩形浮箱连接构成，呈对称分布，主尺度见表1。
表1张力腿平台主参数
平台主尺度及组块有效载荷数值水深/m405排水体积/m350000主吃水/m31浮筒间距/m59浮箱尺寸/m19.5×56.5×25.5浮筒尺寸/m8.5×8.5总预张力/t16064总宽/m78.5干重/t32127
TLP共8根张力筋，12根立管(1根钻井立管，11根生产立管)，张力筋分段参数如表2所示。在莫里森力的计算中，张力筋的拖曳力系数Cd和附加质量系数Ca规定如下：Cd=1.1(上部94.3m长度部分)，Cd=1.0(张力筋余下部分)，Ca=1.0。立管拖曳力系数Cd取1.2，附加质量系数Ca取1。
表2张力筋分段参数
名称外径/m内径/m长度/m单位干重/(t·m-1)EA/NLAJ0.6600.4836.61.31273.19×1010TTS1.0160.93540.01.10352.49×1010MB1.0160.94282.01.03122.34×1010TBS1.0160.93541.51.10292.49×1010SHT0.6600.4985.11.24452.96×1010BTM0.6600.4980.188.66002.96×1010
在HyDrod计算中，张力腿平台在体坐标系下的运动方程如下：
在TLP平台的线性刚度矩阵中，主要包含两部分贡献，一部分为张力筋的贡献，另一部分为浮体共有的静水刚度(静水恢复力)[10]，张力筋的贡献在数量上起主导作用。具体刚度矩阵的各项系数中，起主要作用的仍然是对角阵上的各项，其他项主要是考虑到各自由度之间的相互耦合作用。实际上，刚度矩阵中的各项系数会随着平台位置的移动而发生变化，实质上是一个非线性项，只有平台运动幅度在一定范围内，可以近似认为是常数。本目标张力腿平台刚度矩阵中主要项的计算方法如下[11-12]：
对于纵荡自由度：
横荡自由度：
K32=0
垂荡自由度：
式中：Atendon为张力筋横截面面积，Etendon为张力筋轴向刚度，Apriser为生产立管的横截面积，Epriser为生产立管的轴向刚度，Adriser为生产立管的横截面积，Edriser为生产立管的轴向刚度，ρW为海水的密度，g为重力加速度，S为浮筒横截面积。
横摇自由度：
K34=0
式中：p1、p2为张力腿y方向上的2个位置坐标，q1为生产立管y方向坐标，q2为钻井立管y方向坐标。b为浮筒中心y方向距离。
纵摇自由度：
K35=0
式中D为浮筒直径。
艏摇自由度：
K36=0
在本节中，使用SESAM软件中GeniE模块建立了湿表面模型以及莫里森杆元模型，网格划分尺度1m，湿表面单元14517个，如图1所示。在HydroD模块中建立浮体-张力筋-TTR模型，如图2所示，计算频域下TLP刚度矩阵。
图1湿表面模型
图2频域计算模型
2.1附加质量及刚度矩阵
附加质量系数为一个6×6矩阵，文中分析对目标平台水动力性能起主要作用的主对角线上的4个值，横荡、纵摇、垂荡、艏摇附加质量如图3所示。
图3附加质量系数
由附加质量图可以看出，由于平台的对称性，纵荡和横荡的附加质量是相同的，附加质量取值范围为2.20×107～5.52×107，垂荡方向附加质量范围为2.48×107～2.71×107。
在频域刚度计算的模型中加入了筋腱和立管模型，将其平台运动的影响等效为线性的刚度，通过考虑筋及立管刚度的影响计算得到平台的运动幅值响应。静水回复力矩阵K1以及张力筋及TTR刚度矩阵K2如下所示：
2.2水平/垂向刚度模拟
使用DeepC模块建立全耦合模型，张力筋腱和TTR采用串联弹簧建模方式，每段有着不同的内径外径及刚度，如图4所示。
图4时域计算模型
为了评估此方法建模的精确程度，以及TLP平台的水平方向刚度及水平偏移对垂荡的影响，通过施加水平方向为225°的力来得到平台的准静态偏移曲线。图5为作用力施加在重心处时平台静态偏移曲线；图6为T1、T5(迎浪和顺浪)两根张力筋顶部张力随水平运动变化曲线。
图5平台静态偏移及垂荡曲线
图6张力筋顶部张力曲线
从图5中可以看出，随着作用在平台重心处水平力的增加，平台垂荡及水平刚度并非恒定值，平台整体刚度随着偏移/垂荡的增加而增大。当外力作用使平台产生较大偏移时，平台会产生绕重心的力矩。从图6可以看出，由于力矩的产生，迎浪和背浪的2根张力筋出现了明显不同的张力值。张力筋T1和T5在水平位移为0时，张力筋顶部张力为初始预张力。
这里给出了模型纵荡、横荡、纵摇及艏摇4个自由度上的进行的自由衰减模拟，如图7～10所示。通过后期数据处理，得到平台的运动的固有周期。
图7纵荡自由衰减
图8垂荡自由衰减
图9纵摇自由衰减
图10艏摇自由衰减
在随机海况下的二阶波浪力和频率ωi+ωj处振荡时，会与张力腿平台的垂荡、横摇纵摇产生共振。由于张力腿平台张力筋的刚度使张力腿平台的垂向运动(垂荡、横摇、纵摇)固有周期都相对较低，平台固有周期计算表达式为
式中：Kii是第i自由度上的系统刚度，Mii是i自由度上的结构质量/转动惯量，Aii是i自由度上的附加质量。结合前文求出的到的附加质量Aii及总刚度矩阵中的Kii，将结果与计算公式进行比较，如表3所示。
表3软件模拟与公式计算对比表
名称时域模拟/s公式估算/s误差/%纵荡84.684.410.23横荡84.684.410.23垂荡2.12.090.48横摇2.62.682.99纵摇2.62.724.41艏摇58.150.7114.57
从表3中可以看出，模型进行自由衰减模拟后，纵荡、横荡和垂荡软件拟合结果和数值计算理论结果基本一致。
本文采用时域分析方法计算规则波下所得的张力筋张力响应RAO。计算了在15个不同频率的规则波情况下，张力筋1、5两根张力筋张力响应。分别计算每种规则波(不同周期)情况下张力筋张力趋于稳定周期时的响应幅值，并通过单位化计算得到张力筋张力响应RAO。图11、12分别为T1、T5两根张力筋张力响应的RAO。
图11T1张力筋张力响应RAOs
图12T5张力筋张力响应RAOs
从图中可以看出，不同的来浪角度对不同的张力筋张力RAO有着不同的影响。张力筋张力的固有周期出现在6.5～7.5s处，在90°和270°来浪时极致出现最大值。迎浪和背浪的张力筋RAO也有着不同的响应，T1的极值出现在270°来浪时，极值为934.5kN/m；T5的极值出现在90°来浪时，极值为937.1kN/m。
张力腿平台是深海油气开发中广泛应用的平台型式之一。在张力筋腱的有效约束下张力腿平台可以避开波浪能量集中的频率，使得平台在纵、横摇及垂荡性能上明显改善，大大提高了平台的适用范围、舒适度和安全性。文中通过对一座传统TLP张力筋腱进行研究，结果表明：
1)筋腱及立管分段建模时，能更好地模拟实际工程中的TLP平台刚度，进行静水静力模拟测试，可以验证模型建模的准确性，为后续动力分析提供参考。
2)提出一种模拟验证方法，通过筋腱及立管分段建模，求解平台线性刚度矩阵和附加质量系数矩阵，通过理论公式计算TLP运动的固有周期并与自由衰减模拟结果进行对比，视其结果拟合是否准确，可验证模型计算的准确性。
参考文献：
[1]李军.关于南海海洋平台选型及研究方向探讨[J].广东造船,2010(1):33-35.
[2]ABOU-RAYANAM,EL-GAMALAR.Waveinducedmotionofatriangulartensionlegplatformsindeepwaters[J].Oceansystemsengineering,2013,3(2):149-165.
[3]任顺利,陈建民,梁艳霞.张力腿刚度对平台运动响应影响分析[J].石油矿场机械,2008,37(5):28-31.
[7]YANGCK,KIMMH.TransienteffectsoftendondisconnectionofaTLPbyhull-tendon-risercoupleddynamicanalysis[J].Oceanengineering,2010,37(8/9):667-677.
[8]EL-GAMALAR,ESSAA,ISMAILA.Effectoftetherstensionforceonthebehavioroftriangulartensionlegplatform[J].Americanjournalofcivilengineering&architecture,2014,2(3):107-114.
[9]阎岩,张崎,黄一.基于张力腿平台的顶张紧式立管碰撞分析[J].上海船舶运输科学研究所学报,2012,35(1):1-6.
[10]曾晓辉,沈晓鹏,刘洋,等.考虑多种非线性因素的张力腿平台动力响应[J].海洋工程,2006,24(2):82-88.
[11]杨慧.MOSESTLP在极端环境下的动力性能研究[D].哈尔滨:哈尔滨工程大学,2015.
[12]GUJiayang,YANGJianmin,LVHaining.ComprehensivestudyontendontensionofTLPconsideringwind,waveandcurrent[J].Journalofshipmechanics,2012,16(9):1041-1055.
本文引用格式：
冯丽梅，苏威，闫发锁.张力腿平台筋腱动力特性分析与校验[J].应用科技,2017,44(4):22-27.
FENGLimei，SUWei，YANfasuo.StudyonthedynamiccharacteristicsofTLPtendonswithverification[J].Appliedscienceandtechnology,2017,44(4):22-27.
StudyonthedynamiccharacteristicsofTLPtendonswithverification
FENGLimei1,SUWei2,YANFasuo2
1.CNOOCDeepwaterDevelopmentCo.Ltd.,Shenzhen518000,China2.CollegeofShipbuildingEngineering,HarbinEngineeringUniversity,Harbin150001,China
Keywords：TLPplatform;staticoffset;freedecay;set-down;naturalperiod;tendontension;horizontalstiffness;responseamplitudeoperator
DOI：10.11991/yykj.201703012
收稿日期：2017-03-28.网络出版日期：2017-04-28.
基金项目：工业和信息化部海洋工程装备科研项目——500米水深油田生产装备TLP自主研发.