文章信息
- 刘润, 郝心童, 李成凤, 彭碧瑶. 2020.
- LIU Run, HAO Xin-tong, LI Cheng-feng, PENG Bi-yao. 2020.
- 深海管道轴向定向移动消减方案研究
- Study on deep-sea pipeline walking mitigation
- 海洋科学, 44(9): 112-120
- Marine Sciences, 44(9): 112-120.
- http://dx.doi.org/10.11759/hykx20190320002
-
文章历史
- 收稿日期:2019-03-20
- 修回日期:2019-05-06
铺设于海床表面的深海管道, 在启闭操作和管内温压联合作用下, 极易发生轴向移动(pipeline walking)。海管轴向移动本身不会导致管道的破坏, 但由此造成的连接结构应力超标、钢悬链线立管(steel catenary riser, SCR)丧失预应力及侧向屈曲管段载荷增加是导致海底管道失效的重要原因。Bruton等[1]指出, 现已在很多管道中观察到轴向移动现象, 其中有一起案例, 由于海管的轴向移动, 导致端部连接处断裂, 进而使海管系统发生故障。由此可见, 采用必要的工程措施消减海底管道的轴向移动量至关重要。
关于海底管道轴向移动现象的研究最早见于2000年。Tornes等[2]首先提出了由于瞬态热梯度作用导致的轴向“棘轮”效应, 并且根据管道的膨胀行为, 将管道进行了“长短”分类。此后很多学者针对海管的轴向移动现象进行了相关研究, 早期海底管道轴向移动的研究主要集中在诱发因素方面。2003年, Carr等[3]正式将海管的轴向整体移动命名为“walking”(本文翻译为轴向移动), 并分析了在瞬态热梯度作用与SCR张力作用下管道发生轴向移动的机理; 2008年, Carr等[4]对walking现象的诱发因素(海床斜坡, SCR张力, 瞬态热梯度)进行了全面的描述与分析, 建立了各诱发因素下管道轴向移动的解析解; 2017年, 彭碧瑶等[5]阐述了高温高压海底管道发生轴向移动的运动机理, 建立了单因素与多因素诱发管道轴向定向移动的数值模拟方法, 运用多元非线性回归分析了多因素诱发管道的轴向移动规律。此外, 在管道整体屈曲与轴向移动的耦合[6]以及管土相互作用模型对轴向移动的影响[7]等方面也不断有新的研究成果。
但对消减海底管道轴向定向移动方法的研究相对较少。2003年Carr等[3]提出对管道进行锚固可以抑制管道的轴向移动量, 但提出的最大锚固力理论计算方法与实测值相差较大, 因此后续研究主要采用有限元模拟。Carneiro等[8]提出了一种简化的一维模型, 用来研究管道轴向移动现象, 结果显示一维模型计算得到的轴向移动量与复杂的三维模型相差在10%以内, 但锚固力结果与三维模型相差较大。为了探究不同模型锚固力结果差异的来源, Castelo等[9]在Carneiro等人的研究基础上, 针对同一工况建立了二维模型, 研究表明数据结果的差异源于模型的不同, 与采用的有限元软件无关, 2D模型结果与1D模型非常相近, 与3D模型相差较大, 同时Castelo等指出, 最常用消减管道轴向移动的方法是在管道终端(pipeline end termination, PLET)使用锚链锚固, 或是将PLET滑动基础与吸力桩刚性连接。其他关于消减海底管道轴向移动的措施多见于具体的工程项目, 例如Demel等[10]以西非海岸一个海底管道工程为例, 提出采用海底混合基础的双锚固系统限制管道轴向移动, 并介绍了此锚固系统的设计、制造和安装过程。Carneiro等[11]以Brazilian海底管道工程中两个6 km长的管道为例, 提出采用双枕木防止管道屈曲破坏, 采用端部锚链锚固限制管道轴向移动。综上所述, 对消减海底管道轴向移动措施的研究很少, 尚无锚固力的准确计算方法, 针对具体工程项目的研究并没有揭示规律, 因此不具普适性。
为解决不埋深海管道在温压荷载循环作用下, 轴向移动量不断累积的问题, 本文针对Castelo等[9]提出的消减管道轴向移动的两种方法开展有限元模拟分析, 并对锚链锚固的位置与形式进行优化, 针对一端为SCR, 一端为PLET的典型短管模型, 提出采用中部单向锚链锚固结合端部滑动基础的措施是消减管道产生移动的最有效方法。
1 短管轴向移动诱因的敏感性分析管道中有效轴力最大的点称为虚拟锚固点。升温时, 虚拟锚固点为管道膨胀的中心点, 即管道由虚拟锚固点向两端膨胀; 降温时, 虚拟锚固点为管道收缩的中心点, 即管道由两端向虚拟锚固点收缩。以往研究显示[4], 高温高压海底管道在受到SCR张力、海床倾角及内部瞬态热梯度作用时, 可能会发生轴向移动。其主要机理为, 在以上3种因素作用下, 升温时管道的膨胀中心点与降温时管道的收缩中心点不一致, 导致管道在冷却收缩后无法回到最初的位置, 因此产生轴向移动。管道轴向移动的产生依赖于管道的启闭循环, 若频繁启闭管道, 每一启闭循环产生的轴向移动将会累积, 最终管道会朝向受拉一端, 沿斜坡向下一端及升温滞后一端产生持续的严重的管道整体移动。
为比较SCR张力、海床倾角及内部瞬态热梯度作用对海管轴向移动量的影响, 开展了管道轴向移动诱因的敏感性分析。升温过程中, 流过管道的热流造成的管道温度在时间和空间上分布的不均匀性称为管道内部瞬态热梯度。实际工程中, 管道开启时热流由热端向冷端移动并逐渐加热管道, 管道热端总是先被加热至设计温度, 因此管道内部一定存在瞬态热梯度作用。故在敏感性分析时, 主要针对瞬态热梯度与海床倾角顺向耦合, 瞬态热梯度与SCR张力顺向耦合两种工况, 具体如图 1所示。
图 1两种工况瞬态热梯度方向相同, 将管道坐标x=0 m处设为管道热流入口, 即管道热端, 将管道坐标x=2 000 m处设为管道冷端。工况a中海床整体倾角假设为3°, 管道两端皆为PLET端; 工况b管道一端为PLET端, 一端为SCR端, 在风暴条件下, SCR底部张力可达1 600 kN, 不考虑风暴条件, 管道正常启动下产生的SCR张力一般为200 kN[12]。
应用ABAQUS软件, 采用通过与经典解析解结果对比验证合理的静力算法[5], 建立了管道轴向定向移动数值模型。管道平直无弯曲, 全长2.0 km, 管道外径为0.273 1 m, 管道壁厚为0.014 4 m, 单位管重为689.07 N/m, 膨胀系数取1.1×10–5℃–1, 弹性模量为2.06×1011 Pa, 泊松比取0.3, 建模时管道采用梁单元, 单元类型为PIPE31;土体全长2.1 km, 弹性模量为2.0×107 Pa, 内摩擦角为18°, 土体密度为1 556.6 kg/m3, 粘聚力为1 800 Pa, 泊松比取为0.3, 建模时土体采用实体单元, 单元类型为C3D8R。管土间接触选取罚接触和法向硬接触, 摩擦系数取为0.5。
SCR张力以管道终端集中力的方式施加, 取SCR底部张力为FSCR=200 kN。PLET端采用非线性弹簧模拟, 力-位移关系曲线见图 2。
图 2中以管道向冷端即x轴正向移动时,轴向移动量u为正;管道有效轴力P以拉为正。
管道开启时, 管道沿程逐渐升温, 将管道内部瞬态热梯度作用简化为线性曲线, 图 3为简化的升温曲线。
如图 3所示, 瞬态热梯度为20℃/km, 将此升温曲线写入inp文件。假定海底环境温度为20℃, 管道开启时, 为保证油气的正常输送, 需将管道逐渐加热至150℃; 管道关闭后, 在周围海水的冷却作用下管道均匀降温至环境温度20℃, 将此开启关闭的过程视为一个循环, 用n表示。
图 4为图 1中工况a的P-x图, A表示降温过程的虚拟锚固点, B表示升温过程的虚拟锚固点; 图 5为两种工况管道热端轴向移动量u与循环次数n的关系曲线。
由图 4可得, 降温过程的虚拟锚固点A和升温过程的虚拟锚固点B不重合, 锚固点两侧管道运动方向如图中箭头所示, 由于AB段升温降温过程管道运动方向均指向沿斜坡向下一端, 因此管道向下坡发生整体轴向移动。由图 5所示, 管道在瞬态热梯度与SCR张力顺向耦合情况产生的轴向移动量远大于瞬态热梯度与海床倾角顺向耦合情况。为进一步揭示SCR张力相比其他两种因素对轴向移动的影响规律, 采用图 6所示的三因素综合作用模型, 并分别取SCR张力为0, 100, 200 kN进行计算。
仍将海床倾角取为3°, 瞬态热梯度取为20℃/km, 且将管道坐标x=0 m处设为管道热端。海床倾角使管道产生朝向热端的轴向移动, SCR张力和瞬态热梯度使管道产生朝向冷端的轴向移动。图 7为FSCR= 200 kN时, 管道的P-x图。
由图 7所示, 由于AB段升温降温过程管道运动方向均指向SCR作用端, 因此管道向SCR作用端发生整体轴向移动。降温时, 最大的有效轴力为516 kN; 升温时, 最大的有效轴力为–301 kN。图 8为管道热端的u-n图。
由图 8可知, 当FSCR=0时, 管道会向热端移动, 这是由海床倾角作用产生的, 其每循环轴向移动增量Δu=–0.238 m。当FSCR=100 kN时, 管道会向冷端移动, 这是由于SCR张力和瞬态热梯度的联合作用大于海床倾角的作用, 每循环轴向移动增量Δu=+0.113 m。当FSCR=200 kN时, 每循环轴向移动增量Δu=+0.477 m, 10个温度循环管道产生的整体轴向移动增量达Δu=+4.77 m, 最大轴向移动量u达到+5.98 m, 且轴向移动量u会随着n不断累积, 产生了工程中难以控制的累积轴向移动量, 对海底管道系统安全造成了巨大的威胁。
综上所述, SCR张力对管道轴向移动的影响显著。因此以下针对如图 6所示模型(其中FSCR=200 kN)研究管道轴向移动的最佳消减措施。
2 最佳锚固方式的确定 2.1 锚固方案1) 常见锚链系统组成
减少管道轴向移动的原理大致分为增加轴向摩擦力及施加外力两种[13]。工程上常采用的方法有使用混凝土配重层, 开沟/埋放, 局部抛石/沉排等, 但是这些方法在深海中都存在局限性。工程经验表明, 轴向锚固是一种消除深海管道轴向移动的常用方法, 具有很高的成功率和稳定性。图 9为管道中部锚固示意图。
工程中常采用锚链连接管道与锚桩。升温时, 允许管道自由膨胀通过锚桩; 冷却收缩时, 管道向背离锚桩的方向移动, 当管道回到初始安装位置仍继续朝背离锚桩的方向运动时, 锚链会产生拉力以消减轴向移动量。
2) 锚固计算方案
实际工程中可在管道沿程任何位置进行锚固, 为了得到最有效的锚固方案, 可对锚的位置进行优化。本文提出六种锚固方案, 如表 1所示。其中单向锚仅限制管道发生朝向SCR端的运动; 双向锚则限制管道朝向PLET端和SCR端的运动。
6种锚固方案下分别计算了u与n的关系, 计算结果如图 10所示。
由图 10可得, 相比于图 8不对管道进行锚固的工况, 6种锚固方案每循环管道的轴向移动增量Δu均为0, 即6种锚固方案均有效限制了管道轴向整体移动趋势, 使管道轴向移动量不会随着循环而不断累积。冷端移动量由大至小分别为方案2, 方案1, 方案5、6, 方案3, 方案4;热端移动量由大至小分别为方案4, 方案3, 方案5、6, 方案1, 方案2, 即热端移动量与冷端移动量呈相反趋势, 其中方案5及方案6的热端与冷端u-n图相同。方案4冷端移动量为0, 这是由于锚固在SCR处的双向锚将管道锚固点的位移完全限制, 但方案4热端移动量最大, 这是由于锚固点距离冷端较远, 不能有效控制冷端移动量; 同理, 方案2热端移动量为0, 但冷端移动量最大。而方案5及方案6两端移动量均较小, 这是由于锚固点在中部时, 对两端的约束力较为均衡。虽然六种方案都有效控制了管道的轴向移动, 但每种方案中的管道有效轴力差异很大。
图 11为方案1及方案2的管道有效轴力分布图。图 12为方案1管道降温过程中的有效轴力变化图。
由图 11可知, 升温时, 管道由中点朝向两端膨胀, 方案1由于采用了单向锚固方式, 升温时锚链不受力, 因此有效轴力在管道中部达到最大, 最大值为–333 kN; 方案2由于采用了双向锚固方式, 管道在锚固位置移动量为0, 但有效轴力达到了–559 kN。降温时, 管道由两端向中点收缩, 方案1与方案2的有效轴力分布图相同。可见当热端收缩到初始安装位置时, 两种方案锚链均产生了拉力以限制管道继续收缩, 最大有效轴力均发生在锚固位置处, 约为814 kN。对比图 7可知, 锚固后的最大有效轴力显著增加, 增量为298 kN。
方案1升温时锚链未对管道运动产生限制, 因此锚固力产生于降温过程, 由力的平衡条件可知, 锚固力为814 kN。而方案2升温与降温过程锚链均限制了管道的移动, 升温时, 由1号锚桩提供559 kN拉力; 降温时, 由2号锚桩提供814 kN的拉力。锚固力越大, 对锚桩要求越高, 对工程安全越不利。由图 12可知, 方案1的卸载过程中锚链提供的阻力是逐渐增大的, 直至锚的作用被完全激发达到最大值, 此时管道P也达到了最大值, PLET端的收缩被完全限制。
图 13为方案3及方案4的P-x图。
图 13与图 11相似, 但是由于锚固位置的改变, 单向锚与双向锚在升温阶段均会发挥作用, 以限制管道朝向SCR端的移动。而升温阶段管道内的有效轴力为压力, 过大的压力累积将导致管道发生整体屈曲。对比图 13与图 11可得, 在SCR端锚固比在PLET端锚固积聚的压力更大, 因此其发生整体屈曲的风险更高。
方案3升温时锚链立即产生作用力限制管道移动, 由力的平衡可得, 管道在锚固处的锚固力应为SCR张力与P之和, 即964 kN。而方案4升温与降温锚链均提供力以限制管道移动, 升温时, 由1号锚桩提供765 kN拉力; 降温时, 由2号锚桩提供608 kN的拉力。
在端部进行锚固时, 有效轴力重新分布后, 往往由一端向锚固端逐渐累积, 直到在锚固处达到最大, 由于累积路径较长, 锚固处的有效轴力往往很大, 而在中部锚固, 由于累积路径缩短, 可以解决有效轴力过大的问题。图 14为方案5与方案6的P-x图。
在中部锚固时管道有效轴力分布图与无锚固时(图 7)管道有效轴力分布图相似, 都是由两端向中间逐渐累积, 在管道中部某一位置有效轴力达到最大, 因此相较于在一端锚固, 中部锚固能显著降低有效轴力。由图 14可知, 管道最大有效轴向压力仍位于锚固处, 其值为457 kN。中部双向锚固与单向锚固得到的u-n图及P-x图相同。
2.3 推荐的锚固方案综合考虑SCR端最大移动量、PLET端最大移动量、最大锚固力Fanc, max、最大有效轴向压力Feff, max、有效轴向压力影响长度leff, 对以上6种锚固方案进行比较, 见表 2。
消减方案 | 最大移动量/m | Fanc, max/kN | Feff, max/kN | leff/m | |
PLET端 | SCR端 | ||||
PLET端单向锚 | 0.610 | 2.151 | 814 | –333 | 1 476 |
PLET端双向锚 | 0 | 2.710 | 814 | –559 | 1 476 |
SCR端单向锚 | 2.479 | 0.610 | 964 | –763 | 2 000 |
SCR端双向锚 | 2.479 | 0 | 765 | –765 | 2 000 |
管道中部单向锚 | 1.304 | 1.434 | 278 | –457 | 1 475 |
管道中部双向锚 | 1.304 | 1.434 | 278 | –457 | 1 475 |
由表 2可知, 6种不同的锚固方式对两端最大轴向移动量的结果影响较小; 但不同锚固方式对锚固力的要求不同, 端部锚固产生的锚固力约为中部锚固的2.8~3.5倍。
综合考虑以上几方面因素, 在管道的中部锚固具有较大的优势。在中部锚固的单向锚与双向锚其锚固力与位移控制均具有同样的结果, 即使用双向锚是不必要的。因此可得出结论:对于不考虑整体屈曲的短管, 方案5为最佳锚固方案。
3 锚固与滑动基础的耦合 3.1 端部锚固耦合滑动基础滑动基础在设计之初主要是解决“长管”膨胀与收缩问题的。海底管道因为启闭循环, 常常会出现明显的轴向膨胀与收缩, 而管道的端部结构例如PLET, PLEM等由于具有较大的重量和面积, 导致与管道相比其受到的地基土体摩擦力更大。因此当管道运动受到端部结构的限制时, 会对管道系统施加相当高的负荷。滑动基础允许管端自由移动, 可以有效地解决这一问题。管道轴向移动作为近年来新发现的安全问题, 在设计中往往被忽略。如果滑动基础设计没有考虑管道轴向移动, 那么管道系统将承受设计中没有预料的载荷。Thompson等[14]对传统的滑动基础做出改动, 即将端部锚固与滑动基础联合以解决海管的整体轴向移动问题。滑动基础机构如图 15所示。
图 15中, L1用于管道受热膨胀时滑动基础的预留偏移, L2用于管道降温收缩时向中点的应急偏移。数值模拟中, 为模拟滑动基础与端部锚固耦合, L1=2.5 m、L2=1.0 m, 采用非线性边界条件。即当轴向移动在滑动基础设计的运动范围之内时, PLET可以在基础上自由移动, 当PLET移动到滑动基础端部并进一步发生位移时, 其位移将被完全约束。
图 16, 图 17是一端为SCR, 一端为滑动基础的海底管道模型在启闭循环中的P-x图及u-n图。
分析图 16, 图 17可得, 管道以每循环0.516 m的恒定速度向SCR端移动, 直至在第三循环管道收缩时达到滑动基础右端端部, 此后管道整体移动被限制, 每循环升温过程有效轴力不变, 但降温过程的有效轴力不断增加, 直至n=5时有效轴力趋于稳定, 最大值为808 kN。
3.2 中部锚固耦合滑动基础虽然锚链系统能较好地控制管道的整体移动, 但显著增加了管道的有效轴力; 滑动基础能够在一定程度上放开对位移的限制, 从而降低管道的有效轴力。将二者相结合并选用前述最佳锚固方式即在管道的中部锚固端部采用滑动基础。图 18为锚链系统与滑动基础联合方案与只在中部单向锚方案的管道两端轴向移动量-循环对比图, 图 19为有效轴力对比图。
如图 18, 两种方案位移差别很小。但由图 19所示两种方案有效轴力分布差别较大。采用滑动基础时锚固力约为127 kN, 管道中的最大有效轴向压力为307 kN; 而无滑动基础时, 锚固力约为278 kN, 管道中最大有效轴向压力为457 kN。因此, 相较于仅在中部锚固的情况, 锚链系统与滑动基础相结合可以有效降低锚固力和有效轴向压力, 是最佳的消减管道轴向移动量的方法。
4 结论本文研究了2 km长的海底管道一端连接SCR, 另一端连接PLET时轴向整体移动的消减方案, 确定了最佳锚固点及锚固方式, 分析了滑动基础对消减轴向移动的作用。具体结论如下: (1)相较于海床坡度和瞬态热梯度等诱发因素, SCR张力对海管轴向整体移动的影响更为显著。(2)对管道进行锚固能较好的消减管道轴向整体移动量。但不同的锚固方式所需的锚固力大小不同且显著影响管道内部的轴力分布, 研究提出最佳锚固点宜设置在管道的中部并采用单向锚固方式, 与其他锚固方式相比, 所需的锚固力可降低64%~71%。(3)采用滑动基础与锚桩刚性连接的方式能有效解决短管轴向移动问题。滑动基础能显著地降低管道有效轴力。在中部单向锚的基础上, 两端增加滑动基础能进一步减小锚固力与有效轴向压力, 所需锚固力降低了54%, 最大有效轴向压力减小了33%, 因此对于不易发生整体屈曲的短管, 端部采用滑动基础中部采用锚桩锚固的方式最为适合。
[1] |
Bruton D, Carr M, Sinclair F, et al. Lessons learned from observing walking of pipelines with lateral buckles, including new driving mechanisms and updated analysis models[J]. Iranian Journal of Pediatric Hematology Oncology, 2010, 2(5): 1-5. |
[2] |
Tornes K, Jury J, Ose B A, et al. Axial creeping of high temperature flowlines caused by soil ratcheting[C]//American Society of Mechanical Engineers. Proceeding Conference on Offshore Mechanics and Arctic Engineering. New Orleans: American Society of Mechanical Engineers, 2000: 1229-1240.
|
[3] |
Carr M, Bruton D, Leslie D. Lateral buckling and pipeline walking, a challenge for hot pipelines[C]//Offshore Pipeline Technology Conference. Amsterdam: [s.n.], 2003: 1-36.
|
[4] |
Carr M, Sinclair F, Bruton D. Pipeline walking-understanding the field layout challenges, and analytical solutions developed for the SAFEBUCK JIP[J]. Spe Projects Facilities & Construction, 2008, 3(3): 1-9. |
[5] |
彭碧瑶, 刘润, 王秀妍, 等. 多因素诱发海底管道轴向定向位移量计算方法[J]. 海洋工程, 2017, 35(6): 84-93. Peng Biyao, Liu Run, Wang Xiuyan, et al. Calculation method of multi-factor induced axial directional displacement of submarine pipelines[J]. The Ocean Engineering, 2017, 35(6): 84-93. |
[6] |
Cumming G, Druzynski A, Tørnes K, et al. Lateral walking and feed-in of buckled pipelines due to interactions of seabed features[C]//American Society of Mechanical Engineers. Proceeding Conference on Offshore Mechanics and Arctic Engineering. Hawaii: American Society of Mechanical Engineers, 2009: 739-751.
|
[7] |
Rong H, Inglis R, Bell G, et al. Evaluation and mitigation of axial walking with a focus on deep water flowlines[C]//Offshore Technology Conference. Houston: [s.n.], 2009: 19862.
|
[8] |
Carneiro D, Murphy D. Simple numerical models for pipeline walking accounting for mitigation and complex soil response[C]//American Society of Mechanical Engineers. Proceeding Conference on Offshore Mechanics and Arctic Engineering. Rotterdam: American Society of Mechanical Engineers, 2011: 651-656.
|
[9] |
Castelo A, Galgoul N S. Analyses of different FE models for HP/HT-pipeline walking with sliding end structures[C]//American Society of Mechanical Engineers. American Society of Mechanical Engineers Digital Collection. Rio de Janeiro: American Society of Mechanical Engineers, 2013: 513-519.
|
[10] |
Demel J, Wallerand R, Rebours N, et al. Pipeline walking mitigation by anchor in-line structures with the use of a hybrid subsea foundation[C]//Offshore Technology Conference. Houston: [s.n.], 2016: 1-22.
|
[11] |
Carneiro D, Gouveia J, Parrilha R, et al. Buckle initiation and walking mitigation for HP/HT pipelines[C]//Proceedings of the Deep Offshore Technology International Conference. Monte Carlo: [s.n.], 2009: 1-17.
|
[12] |
Brunner M, Qi X, Zheng J, et al. Combined effect of flowline walking and riser dynamic loads on HP/HT flowline design[C]//Offshore Technology Conference. Houston: [s.n.], 2006: 1-5.
|
[13] |
白勇, 白强. 海洋工程设计手册, 海底管道分册[M]. 上海: 上海交通大学出版社, 2014: 274-276. Bai Yong, Bai Qiang. Marine Engineering Design Manual, Subsea Pipeline[M]. Shanghai: Shanghai Jiao Tong University Press, 2014: 274-276. |
[14] |
Thompson H, Reiners J, Brunner M, et al. Tahiti flowline expansion control system[C]//Offshore Technology Conference. Houston: [s.n.], 2009: 1-22.
|