2. 上海外高桥造船有限公司, 上海 200000;
3. 中国海洋大学 山东省海洋工程重点实验室, 山东 青岛 266100
2. Shanghai Waigaoqiao Shipbuilding CO., LTD., Shanghai 200000, China;
3. Shandong Provincial Key Laboratory of Ocean Engineering, Ocean University of China, Qingdao 266100, China
随着海洋油气开采逐步走向深海,为适应更加严峻的海洋条件,不同类型的海洋工程结构应运而生。钢悬链式立管 (steel catenary riser,SCR) 是20世纪90年代初期发展起来的一种新型深水立管系统,具有成本低和可承受高温高压等优点,是浮式平台输入输出的首选立管系统。钢悬链式立管集海底管线与立管于一身,一端连接井口,一端连接浮式结构[1],迄今为止,已有大量SCR安装在西非、墨西哥湾及巴西等海域。
钢悬链式立管由悬垂段和海底流线段组成,悬垂段和海底流线段是由钢管焊接而成的一条整体管线,悬垂段和海底流线段的分界点是触地点。当钢悬链式立管的悬垂段随浮式平台运动而发生运动时,或者在波浪等外荷载的作用下,海底流线段也随之产生运动,从而导致了钢悬链式立管与海床的相互作用[2-3]。研究表明,钢悬链式立管与海床相互作用对立管的动力响应和疲劳寿命有较大的影响[4-6]。然而,由于立管与海床相互作用的非线性及其对疲劳性能影响的不确定性,目前工程业普遍采用设计规范推荐的过度简化的方法来研究立管与海床土的相互作用[7],即采用线性刚度海床土来模拟管土相互作用。Brigde等对钢悬链式立管的研究结果表明疲劳寿命的评估与选用的海床土的刚度有关[8]。因此为了能够准确的评估立管触地点处的疲劳寿命,需要采用合理的海床模型。
根据相关的试验结果及数值分析结果,Aubeny提出了模拟立管与海床相互作用的非线性P-y曲线,其中P为土抗力,y为立管的贯入深度[9]。Nakhaee基于Aubeny提出的立管与海床土非线性相互作用模型对CABLE 3D程序进行了改编[10-11],以更真实地模拟海床土的刚度。在其研究中,立管被认为是经历周期性振动,即立管与海床土是完全的接触、分离及再接触。本文基于Aubeny提出的P-y曲线进行数学编程建模,程序考虑了钢悬链式立管与海床相互作用时可能位于弹性回弹、部分分离或完全分离的情况,是一个完整的钢悬链式立管与海床土非线性相互作用的力学模型。另外,对不同海床土壤刚度的计算结果进行了比较。
1 立管与海床相互作用模型 1.1 模型描述图 1为Aubeny提出的P-y曲线,立管与海床相互作用的P-y曲线可以分为两个阶段:1) 沿着路径0-1-1′的骨干曲线,即立管初次贯入土壤时遵循的路径;2) 大变形的回弹路径1-2-3-1,包括管土完全接触即弹性回弹路径1-2,管土的部分分离,即路径2-3以及沿着路径3-1的再接触与再加载。在这个过程中,立管向上运动时的管土完全分离并不会在每一种工况及触地点区域所有的点发生。在沿着滞回圈任何一点都可能发生反向的位移,如图 1中虚线所示。同时,周期荷载也可能完全在滞回圈内。其中从2点到P=0处的点以及部分分离阶段2-3,立管均受到吸力作用。Aubeny分别提出了骨干曲线,弹性回弹、部分分离及再接触阶段的经验公式[9, 12]。
1.2 程序编译本文用于计算锚泊或立管系统与浮式平台相互作用的非线性有限元分析程序基于大挠度柔性梁理论。柔性梁理论由Garrett首先提出,适用于不可伸长的杆[13]; 而后Paulling等[14-15]将其扩展为适用于小伸长的杆; Chen将其进一步拓展为适用于大伸长的锚泊线[16]。
根据Chen的论文,小伸长柔性杆的运动方程和约束条件分别为
$ \mathit{\boldsymbol{M\ddot r + }}{\left( {\mathit{\boldsymbol{Br''}}} \right)^{\prime \prime }} - {\left( {\bar \lambda \mathit{\boldsymbol{r'}}} \right)^\prime } = \mathit{\boldsymbol{q}} $ | (1) |
$ \mathit{\boldsymbol{r'}} \cdot \mathit{\boldsymbol{r'}} = {\left( {1 + \frac{{\bar \lambda - {P_f}{A_f} + {P_i}{A_i}}}{{{\rm{EA}}}}} \right)^2} $ | (2) |
式中:M为质量矩阵,r为杆的曲线坐标在笛卡尔坐标中的瞬时形态向量,B为弯曲刚度,
垂向海床土的支撑力可以表示为
$ {P^t} = P_i^t + k_i^t\delta y $ | (3) |
式中:Pt是在时间步t时的真实的土抗力,t是时间步,i是迭代次数。Pit和kit是在时间步t时,根据i-1步的迭代的位移结果计算得到的海床土抗力和海床土刚度。δy=yt-yit是真实的贯入深度和计算得到的贯入深度的差值。
在静态问题中,t恒等于0。计算时,首先将支撑力P00假设为等于立管在水中的重量,再根据骨干曲线得到土壤刚度k00,最终将其代入方程 (1) 中,即可得到新的贯入深度y10。然后再用新的贯入深度根据骨干曲线得到新的支撑力及土壤刚度,重复该迭代过程直到|yi+10-yi0|≤ε,其中ε为规定的精度要求。
在动态分析时,这个过程会更加复杂。具体分成以下几步:
1) t=1,i=1,立管的贯入深度y10根据静态分析的结果确定,支撑力P11和土刚度k11依然是根据骨干曲线确定。最终,新的贯入深度y11通过方程 (1) 得到。然后,i=i+1,移至步骤3;
2) t>1,i=1,首先根据Newmark-β方法预测贯入深度,然后利用其确定支撑力P1t和土壤刚度k1t。但是在这个过程中,要先判断立管正在经历的阶段,才能选择正确的经验公式计算支撑力及土刚度。最终通过方程 (1) 可以得到新的贯入深度y1t,然后i=i+1,移至步骤3;
3) t≥1,i>1,在上一次迭代过程中,方程 (1) 中代入Pi-1t和ki-1t计算得到的贯入深度yi-1t用于确定土抗力Pit和土刚度kit。这个过程同样需要判断立管正在经历的加载或卸载阶段。重复该迭代过程直到|yi+1t-yit|≤ε。然后t=t+1,并且i=1,移至步骤2。
2 模型计算结果及分析 2.1 模型验证为了验证所提出模型的正确性,采用上述模型编制了钢悬链式立管与海床土相互作用的分析程序CABLE PSI,该程序充分考虑了立管与海床土的非线性相互作用,考虑了海床土壤特征对于立管法向约束力的影响。将模型计算结果同现有商用软件OrcaFlex的计算结果进行比对。本文数值模拟选取的是一根典型的钢悬链式立管[17-18],立管的上端连接在船体上,在动态模拟过程中仅考虑船体的垂荡运动,垂荡幅值为2 m,垂荡周期为12 s。该立管长2 400 m, 外径0.355 6 m,内径0.304 8 m,密度7 850 kg/m3, 弹性模量2.07×108 kPa,模拟水深1 100 m,顶张力2 100 kN,悬挂角16°。数值模型的计算参数如下:时间步长为0.02 s,计算总时长为3 700 s,单元长度为6 m。
海床土的剪切强度是随着深度线性增加的,即满足
$ S = {S_0} + {S_g}y $ | (4) |
式中:S在深度y处的土壤剪切强度,S0是在海床处的土壤剪切强度,Sg是剪切强度梯度。本文中,S0为1.8 kPa,Sg为1.0 kPa/m。
分别用CABLE PSI和OrcaFlex对立管进行静态分析。图 2给出了他们的静态位置对比情况,x、y分别为立管节点的横坐标与纵坐标,通过比较发现,两者基本上呈现相同的静态位形,由于OrcaFlex中钢悬链式立管是通过没有质量的弹簧连接的质量块简化模拟的,而在CABLE PSI中钢悬链式立管是通过大挠度柔性梁模拟的,所以两者之间存在一定差异是合理的。
在静态分析结束后,分别用CABLE PSI和OrcaFlex对立管进行动态分析,由于两者所采用的管土相互作用模型不同,且OrcaFlex并不能考虑土壤的吸力作用,因此两者计算结果会有一定差异。图 3给出了两者立管动态分析中触地点区域的垂向位移时程图,通过比较发现,两者程序都趋于稳定之后,触地点区域的动态位移幅值虽然并不完全相同,但均为0.06 m左右。图 4给出了相应的频谱分析,从频谱图比较可以看出,图 4(a)和图 4(b)均有两个主要频率,且都在0.08 Hz和0.16 Hz附近。由此可见CABLE PSI程序的计算结果与商用软件OrcaFlex的计算结果能较好吻合。
2.2 载荷位移曲线
本文CABLE PSI程序所建立的数值模型采用载荷位移曲线,即P-y曲线来模拟立管与海床土的相互作用,将海床的土壤刚度看作是非线性的,而不是一个定值,那么利用该程序计算得到的海床土法向约束力与立管位移之间的关系曲线也应该与该模型的P-y曲线符合。
图 5给出了触地点区域某一节点的法向约束力与该节点位移之间的关系曲线,即实线部分,而虚线为利用骨干曲线经验公式绘制出的骨干曲线。土抗力大于0,表明此时立管受到土壤的支撑力作用;土抗力小于0,表明此刻立管发生上举运动,受到土壤的吸力作用,并且该节点处SCR与土壤接触后,发生完全分离,再与土壤再次接触。因此,由图 5可知,该节点的荷载位移曲线遵循采用的P-y曲线。从而验证了该程序在管土相互作用的过程中确实是采用P-y曲线控制海床土法向约束力与SCR位移之间的关系,即该程序模型能够对立管与海床土的非线性相互作用进行分析。
2.3 海床土非线性刚度与线性刚度对比线性刚度的海床其刚度只与SCR的单位长度湿重和直径有关,而非线性刚度的海床其刚度与真实的海床土地基有关,通过P-y曲线的形式表现出来。本节对海床土非线性刚度与线性刚度的计算结果进行了比较,其中非线性刚度海床土壤强度分为低强度、中强度和高强度三种情况,具体参数见表 1,其他计算参数与2.1节相同。
强度 | 参数 | 数值 |
低强度 | 未排水剪切强度/ kPa | 1.2 |
剪切强度梯度/ (kPa·m-1) | 0.8 | |
中强度 | 未排水剪切强度/ kPa | 1.8 |
剪切强度梯度/ (kPa·m-1) | 1.0 | |
高强度 | 未排水剪切强度/ kPa | 2.4 |
剪切强度梯度/ (kPa·m-1) | 1.2 |
图 6给出了线性刚度海床土、低强度海床土、中强度海床土和高强度海床土在SCR最大贯入深度时对应的触地点区域的立管位形图。可以看出海床土刚度的不同对SCR的位形有着较显著的影响。在垂荡幅值皆为2 m的情况下, 低强度海床土的立管最大贯入深度为0.113 m,中强度海床土的为0.058 m,高强度海床土的为0.038 m,线性刚度海床土的为0.189 m,可见低强度海床土是分析的三种非线性土壤中最大贯入深度与线性海床土结果最为相近的。因此通过线性海床土模型计算得到的结果其所能表征的是土壤刚度较低的海床,而非线性海床土模型能考虑不同海床土壤刚度对钢悬链式立管的影响。
图 7给出了分析时间为500~800 s的这300 s的时间历程中不同土壤刚度情况下,最大贯入深度节点处SCR垂向位移响应时程图和频谱图。从图中可以看出不同土壤刚度情况下,SCR该节点的垂向位移响应均呈现出明显的周期性,且幅值变化不大:
1) 线性刚度海床土情况下最大贯入深度节点处SCR垂向位移响应稳定在[-1 100.19, -1 100.145](以海平面为位移零点),幅值约为0.045 m;
2) 低强度海床土情况下该节点的垂向位移响应稳定在[-1 100.12, -1 099.97],幅值约为0.05 m;
3) 中强度海床土情况下该节点的垂向位移响应稳定在[-1 100.05, -1 099.92],幅值约为0.07 m;
4) 高强度海床土情况下该节点的垂向位移响应稳定在[-1 100.04, -1 099.99],幅值约为0.05 m。
可以看出不同土壤刚度情况下,SCR的垂向运动响应还是有一定差异的。
从图 7的频谱分析中可以看出不同土壤刚度情况下,SCR的振动频率成分皆为0.083 Hz和0.166 Hz,其中前者恰为立管上部的垂向运动频率,可见SCR的振动频率主要是受到外荷载的振动频率的影响,与海床土的性质关系不大。
表 2给出了线性刚度海床土与低强度海床土这两种情况下,钢悬链式立管的最大弯矩值、最大弯矩幅值及其位置,最大应力值、最大应力幅值及其位置和对应的最大年疲劳损伤率,其中位置是指从SCR上端开始沿立管轴向的距离。通过比较可知,两者的最大弯矩幅值和最大应力幅值差异非常显著,线性刚度土壤情况下,最大弯矩幅值为86.22 kN·m,而低强度土壤情况下,最大弯矩幅值为105.77 kN·m,后者增幅22.7%。线性刚度土壤情况下,最大应力幅值为43.06 MPa,而低强度土壤情况下,最大应力幅值为54.00 MPa,后者增幅25.4%。可见海床土刚度对于触地点区域弯矩幅值及应力幅值有较大影响。年疲劳损伤率的计算结果也验证了海床土刚度对于触地点区域的疲劳寿命影响较大。线性刚度土壤情况下年疲劳损伤率为0.023 8,而低强度土壤情况下年疲劳损伤率为0.051 1,后者增幅114.7%。另外,计算得到中强度土壤情况下立管的年疲劳损伤率为0.059 8,高强度土壤情况下为0.099 8。
计算结果 | 线性刚 度土壤 |
低强度 土壤 |
增幅/% |
最大弯矩值/(kN·m) | 194.2 | 197.95 | 1.9 |
最大弯矩位置/m | 1 392 | 1 392 | — |
最大弯矩幅值/(kN·m) | 86.22 | 105.77 | 22.7 |
最大弯矩幅值位置/m | 1 464 | 1 464 | — |
最大应力值/MPa | 101.80 | 103.03 | 1.2 |
最大应力值位置/m | 1 392 | 1 392 | — |
最大应力幅值/MPa | 43.06 | 54.00 | 25.4 |
最大应力幅值位置/m | 1 464 | 1 464 | — |
年疲劳损伤率 | 0.023 8 | 0.051 1 | 114.7 |
本文采用大挠度柔性梁模型模拟钢悬链式立管的运动,同时采用弹性基础梁模型模拟海床,而管土相互作用过程采用荷载位移曲线P-y进行模拟,建立了钢悬链式立管与非线性刚度海床的相互作用模型,具有一定先进性,开发了钢悬链式立管分析软件CABLE PSI。基于本文的研究,可以得出以下结论:
1) 通过与商用软件初步比较,CABLE PSI程序的计算结果与商用软件OrcaFlex的计算结果能较好吻合,验证了所建立模型的正确性。
2) 线性土壤与非线性土壤相比,线性土壤刚度SCR的贯入深度较大,触地点区域的SCR垂向位移幅值有一定差异。但是在上部结构的周期运动情况下,SCR的振动频率主要受外界荷载的周期影响,与海床土的性质关系不大。
3) 海床土壤的刚度变化对立管的运动响应、弯矩以及疲劳寿命有着显著影响,非线性土壤的立管疲劳损伤率是线性土壤的2~4倍多,这说明采用线性刚度土壤模拟管土相互作用与实际情况有一定差距。
钢悬链式立管与海床土之间的相互作用是一个比较复杂的过程,除了海床土的刚度非线性问题,还与上部结构的运动特征、海床土的物理力学性能、顶张力以及立管本身的结构形式等多种因素有关,因此,其模拟和计算方法有待深入研究。
[1] |
黄维平, 李华军. 深水开发的新型立管系统——钢悬链线立管 (SCR)[J].
中国海洋大学学报, 2006, 36(5): 775–780.
HUANG Weiping, LI Huajun. A new type of deepwater riser in offshore oil & gas production:the steel catenary riser, SCR[J]. Periodical of ocean university of China, 2006, 36(5): 775–780. |
[2] |
黄维平, 孟庆飞, 白兴兰. 钢悬链式立管与海床相互作用模拟方法研究[J].
工程力学, 2013, 30(2): 14–18.
HUANG Weiping, MENG Qingfei, BAI Xinglan. The simulation of the interaction between SCR and seabed[J]. Engineering mechanics, 2013, 30(2): 14–18. |
[3] |
刘娟, 黄维平. 钢悬链式立管涡激振动流固耦合非线性分析方法研究[J].
振动与冲击, 2014, 33(3): 41–45.
LIU Juan, HUANG Weiping. Fluid-structure interaction analysis for VIV of Steel catenary risers[J]. Journal of vibration and shock, 2014, 33(3): 41–45. |
[4] |
王坤鹏, 薛鸿祥, 唐文勇. 基于海床吸力和刚度衰减模型的深海钢悬链线立管动力响应分析[J].
上海交通大学学报, 2011, 45(4): 585–589.
WANG Kunpeng, XUE Hongxiang, TANG Wenyong. Dynamic response analysis of deepwater steel catenary riser based on the seabed-suction and stiffness-degradation model[J]. Journal of Shanghai Jiaotong university, 2011, 45(4): 585–589. |
[5] |
王坤鹏, 薛鸿祥, 唐文勇. 基于全耦合模型和管土作用模型的深海悬链线立管触地区域疲劳特性分析[J].
上海交通大学学报, 2014, 48(4): 576–582.
WANG Kunpeng, XUE Hongxiang, TANG Wenyong. Fatigue characteristics analysis of steel catenary riser near touchdown zone based on full coupled model and SCR-soil interaction model[J]. Journal of Shanghai Jiaotong university, 2014, 48(4): 576–582. |
[6] | WANG Kunpeng, XUE Hongxiang, TANG Wenyong, et al. Fatigue analysis of steel catenary riser at the touch-down point based on linear hysteretic riser-soil interaction model[J]. Ocean engineering, 2013, 68: 102–111. DOI:10.1016/j.oceaneng.2013.04.005 |
[7] | SHIRI H. Response of steel catenary risers on hysteretic non-linear seabed[J]. Applied ocean research, 2014, 44: 20–28. DOI:10.1016/j.apor.2013.10.006 |
[8] | BRIDGE C, LAVER K, CLUKEY E, et al. Steel catenary riser touchdown point vertical interaction models[C]//Proceedings of 2014 Offshore Technology Conference. Houston, Texas, USA:Offshore Technology Conference, 2014:16628-MS. |
[9] | AUBENY C P, BISCONTIN G, ZHANG J. Seafloor interaction with steel catenary risers (Finial Project Report)[R]. Texas:Texas A&M University, 2006. |
[10] | NAKHAEE A, ZHANG J. Dynamic interaction between SCR and the seabed[C]//Proceedings of ASME 2007 26th International Conference on Offshore Mechanics and Arctic Engineering. San Diego, CA, USA:ASME, 2007. |
[11] | NAKHAEE A. Study of the fatigue life of steel catenary risers in interaction with the seabed[D]. Texas:Texas A&M University, 2010. |
[12] | JIAO Yaguang. Non-linear load-deflection models for seafloor interaction with steel catenary risers[D]. Texas:Texas A&M University, 2007. |
[13] | GARRETT D L. Dynamic analysis of slender rods[J]. Journal of energy resources technology, 1982, 104(4): 302–306. DOI:10.1115/1.3230419 |
[14] | PAULLING J R, WEBSTER W C. A consistent, large-amplitude analysis of the coupled response of a TLP and tendon system[C]//Proceedings of the Fifth International Symposium on Offshore Mechanics and Arctic Engineering. Tokyo, Japan, 1986, 3:126-133. |
[15] | MA W, WEBSTER W C. Analytical approach to cable dynamics:theory and user manual[R]. Berkeley, CA:University of California, 1994. |
[16] | CHEN Xiaohong. Studies on dynamic interaction between deep-water floating structures and their mooring/tendon systems[D]. Texas:Texas A&M University, 2002. |
[17] |
刘娟. 考虑刚体摆动的钢悬链式立管出平面运动分析方法研究[D]. 青岛: 中国海洋大学, 2013: 58-60.
LIU Juan. Study of out-of-plane motion of SCRs with rigid swing[D]. Qingdao:Ocean University of China, 2013:58-60. |
[18] |
杨超凡. 钢悬链式立管与海床土非线性相互作用研究[D]. 青岛: 中国海洋大学, 2014: 35-40.
YANG Chaofan. Study on nonlinear interaction between SCR and seabed[D]. Qingdao:Ocean University of China, 2014:35-40. |