2021, Vol. 43
2. 天然气水合物国家重点实验室，北京 100028;
3. 国家能源深水油气工程技术研发中心，北京 100028
2. National Research and Development Center for Energy Engineering Technology in Deep Water, Beijing 100028, China;
3. National Energy Deepwater oil and Gas Engineering Technology Research and Development Center, Beijing 100028, China
FLNG是一种集海上天然气生产、储存、外输功能于一身的新型海工装备。目前FLNG采用的外输方式有旁靠外输和尾输2种，其中尾输方式由于作业时两船相对距离远，作业时容许的环境载荷更大，被认为更适合在环境条件相对恶劣的海区应用。
在过去的十年中，海洋工业界与学术界对FLNG尾输作业开展了大量研究工作并获得了许多有价值的成果[1-5]。然而在以往的研究中，尾输过程中FLNG于LNG运输船两船体的水动力响应被独立进行处理，但在实际过程中，FLNG和LNG运输船之间的距离在外输过程中一直变化，以technipFMC公司的ATOL外输系统为例，要求在外输过程中FLNG与LNG运输船之间距离保持在60～90 m范围，这一距离与FLNG船宽接近，可能发生两船体的水动力耦合现象，最不利条件下可能导致两船碰撞。
本文根据中场理论，对外输过程中FLNG和LNG运输船耦合水动力响应进行研究，并分析一些因素，如两船距离、装载工况组合以及不同环境条件的影响。结果表明，随着船体距离的降低，FLNG对LNG运输船的遮蔽效应增强，而非耦合的水动力分析结果将高估作用于LNG运输船上的低频漂移力与两船间距。风流载荷但有助于提升整体系统的稳定性，减小船体相对运动和张力的变化范围，从而降低两船碰撞的风险。
1 设计基础 1.1 FLNG与LNG运输船FLNG主尺度参数为船长340.0 m，船宽61.0 m，型深37.0 m，有10个GTT液舱，整体LNG容积为252000 m3。LNG运输船主尺度参数伟船长291.0 m，船宽47.0 m，有4个GTT液舱，LNG容积为173500 m3。FLNG与LNG运输船的总布置图如图1所示。FLNG与LNG运输船在外输过程中具有3种典型装载组合，具体参数如表1所示。
|
图 1 FLNG与LNG运输船总布置图 Fig. 1 Main arrangement plan of FLNG and LNG carrier |
|
表 1 外输过程中典型装载组合 Tab.1 Load condition combination during offloading process |
FLNG采用15根系泊缆绳呈3组进行定位，所有缆绳在分段结构上完全相同，主要分为3段，包括顶部锚链，中部钢缆与底部锚链。底部锚链长1392 m，直径为146.05 mm，破断拉力为18908 kN，中部由螺旋钢缆组成，长1859 m 直径139.7 mm，破断拉力为19186 kN，顶部锚链与底部锚链规格相同，长100 m。缆绳几部分之间采用卸扣连接，每根缆绳采用吸力锚定位，系泊水平跨距为2815 m，顶部预张力约为5150 kN。
1.3 尾输系泊系统尾输过程中，FLNG与LNG运输船之间采用4根连接缆控制相对运动，尼龙缆长90 m，直径139.00 mm，破断拉力为1819 kN，风暴刚度为6000 kN，空气重量为5.2 kg/m。
2 耦合系统数值模型 2.1 水动力分析模型采用Hydrostar软件建立带有液舱的FLNG与LNG运输船双船体耦合水动力模型，如图2所示。
|
图 2 FLNG和LNG运输船的水动力模型 Fig. 2 Hydrodynamic model of FLNG and LNG carrier |
|
图 3 两船体之间水面的划分 Fig. 3 Mesh model of water plane between FLNG and LNG carrier |
FLNG与LNG运输船船体上的2阶波浪漂移力根据中场理论计算，为了考虑两船体水动力的耦合作用，需要对船体之间的水面进行网格划分并给定能量耗散系数（通常根据水池模型试验获得），以抑制此区域内的波面共振。
2.2 耦合水动力分析分析采用法国船级社的Ariane软件进行FLNG尾输作业耦合运动分析，建立FLNG和LNG运输船、单点系泊系统与尾输临时系泊系统如图4所示，尾输系泊缆绳的布置如图5所示。
|
图 4 FLNG单点系泊系统布置 Fig. 4 Single point mooring system of FLNG |
|
图 5 尾输系泊系统 Fig. 5 Tandem mooring system |
在进行耦合运动分析之前，首先计算两船体的水动力性能。分别采用耦合分析法考虑两船体之间间距为40 m，60 m，80 m时计算水动力结果并与非耦合分析结果进行对比，如图6～图8所示。
|
图 6 两船垂荡运动RAOs对比 Fig. 6 Heave RAO comparison of two vessels |
|
图 8 两船纵荡漂移力RAOs对比 Fig. 8 Surge drift force RAO comparison of two vessels |
图6和图7表明：不同间隙距离对于FLNG船体垂荡与纵摇运动影响很小，且耦合分析与非耦合分析几乎没有差异；间隙距离对LNG运输船船体垂荡与纵摇运动影响相对较大，其中耦合分析的结果在0.4 rad/s～1.0 rad/s范围明显小于非耦合分析，特别是当波浪频率等于0.7 rad/s时。对于FLNG纵荡运动的波浪漂移力而言（见图8），仅在当波浪频率大于0.6 rad/s时，耦合分析结果略微小于非耦合分析。但对于LNGC而言，根据远场或近场理论的非耦合分析结果会明显大于根据中场理论得到的耦合分析结果。当两船相对距离减小时，RAO曲线的幅值也将明显降低，特别是当波浪频率大于1.0 rad/s后，这表明耦合效应在增强，同时FLNG对LNGC在波浪漂移力方面具有明显的遮蔽效应。
|
图 7 两船纵摇运动RAOs对比 Fig. 7 Pitch RAO comparison of two vessels |
|
图 9 两船垂荡运动RAOs对比 Fig. 9 Heave RAO comparison of two vessels |
|
图 11 两船纵荡漂移力RAOs对比 Fig. 11 Surge drift force RAO comparison of two vessels |
图9与图10表明在不同装载组合条件下，FLNG与LNG运输船的运动RAO变化较小。图11表明当LNG运输船增加排水量或吃水之后，作用于LNG运输船船体上的纵荡漂移力也将上升。
|
图 10 两船纵摇运动RAOs对比 Fig. 10 Pitch RAO comparison of two vessels |
|
图 12 两船垂荡运动RAOs对比 Fig. 12 Heave RAO comparison of two vessels |
|
图 13 两船纵摇运动RAOs对比 Fig. 13 Pitch RAO comparison of two vessels |
|
图 14 两船纵荡漂移力RAOs对比 Fig. 14 Surge drift force RAO comparison of two vessels |
由图12～图14可知，波浪入射角度对于两船体的运动与漂移力均具有很大的影响，当波浪入射角度增加时，两船体的垂荡与纵摇运动都将加剧。对于LNG运输船而言，当波浪呈180°入射时，其所受到的波浪漂移力幅值将明显小于波浪呈135°方向入射时，这也表明FLNG对LNG运输船在波浪漂移力方面影响显著，相反，LNG运输船对FLNG的影响却十分有限。
上述分析表明，在水动力影响因素中，波浪入射角度对水动力结果影响最大，其次为两船之间的间距，而装载工况组合的影响最小。此外非耦合的水动力分析将会过高估计LNG运输船上作用的平均波浪漂移力。
3.2 耦合分析结果由于外输工况下FLNG与LNG运输船之间的间距持续变化，因此分别采用不同间距下的水动力结果（分别用HD1，HD2，HD3，HD4以表示两船之间间距为40 m，60 m，80 m和非耦合的状态），并将计算结果输入至Ariane软件进行运动分析对比。
图15～图18为LC2装载工况时采用不同水动力结果下的耦合运动分析结果。环境条件假设为仅波浪载荷作用，有义波高Hs为5 m，谱峰周期Tp为10 s，波浪谱为Jonswap，谱峰因子为2.0。
|
图 15 两船之间距离时间历程 Fig. 15 Time history of distance between two vessels |
|
图 16 两船之间距离频谱 Fig. 16 Spectrum of distance between two vessels |
|
图 18 系泊缆张力距离频谱 Fig. 18 Spectrum of hawser tension |
图15表明，在仅波浪作用下，两船相对水平向运动有一个较大的范围，特别是当采用HD2水动力计算结果时，两船最小间距仅为30 m，极易发生碰撞。图17和图18显示了系泊张力基本有低频成分组成，采用非耦合水动力结果HD4所得到的系泊张力明显大于采用HD1～HD3的耦合水动力分析结果，这一现象也表明FLNG在外输作业过程中遮挡了作用在LNGC上的波浪载荷。
|
图 17 系泊缆张力时间历程 Fig. 17 Time history of hawser tension |
为了研究环境载荷入射角度的影响，根据参考文献[9]，选取5个典型的环境载荷方向组合，包括全平行工况，交叉流工况，交叉风工况，交叉风流工况以及对称风流工况，具体的入射角度如图19所示。风速Vwind为20.0 m/s，流速Vcurrent为1.0 m/s，并采用HD3的水动力计算结果，全部的计算工况与相应编号如表2所示，相对距离与系泊张力结果如表3所示。
|
图 19 5种典型外输作业环境载荷方向组合 Fig. 19 Five typical directional combination in tandem offloading operation |
|
|
表 2 全部计算工况 Tab.2 Total calculation cases |
|
|
表 3 相对距离与系泊张力结果 Tab.3 Results of distance and hawser tension |
由表3可知，在LC3装载工况下，系泊缆绳张力最小，而在LC1装载工况下，两船体之间相对运动范围最小。对于环境载荷组合而言，第4个工况（交叉风流工况）和第5个工况（对称风流工况）比其他工况产生了更加极端的响应，但根据缆绳的破断强度，最大系泊张力仍然可以满足强度要求。与图15和图17相比较发现，当考虑风与流载荷之后，船体相对运动和系泊缆绳张力的变化范围和最大值都有所降低，这说明风流载荷从某种程度上提升了整体系统的稳定性。
4 结 语本文采用数值计算方法研究FLNG尾输过程的动力特性。根据中场理论计算FLNG与LNG运输船的双船耦合水动力，并建立了包括FLNG，LNG运输船、单点系泊系统以及外输定位缆绳的耦合分析模型并进行运动求解，得到结论如下：
1）根据中场理论，在FLNG进行尾输过程中FLNG船体与LNG运输船将出现水动力耦合效应，非耦合水动力分析将高估作用于LNG运输船体上的纵荡漂移力。在对水动力分析的影响因素中，波浪入射角度对结果的影响最大，其次为两船之间的相对距离，而装载工况组合的影响最小。
2）在纯波浪工况下，两船之间的相对距离范围较大并可能导致两船碰撞，而风流载荷将在某种程度上降低两船体相对运动和系泊张力的变化范围，改善整体系统的稳定性。
| [1] |
JUN-BUMN Rho, ALEXANDER A. Korobkin, coupled analysis of deepwater floating system including VIV in time domain[C]// Proceedings of the ASME 32nd International Conference on Ocean, Offshore and Arctic Engineering, 2007, OMAE, 29523.
|
| [2] |
JAAP de Wilde, JOHAN Dekker, Review of hydrodynamic and nautical studies for offshore LNG operations[C]// Proceedings of the ASME 32nd International Conference on Ocean, Offshore and Arctic Engineering, 2013, OMAE, 11366.
|
| [3] |
ZHAO Wenhua, YANG Jianmin, HU Zhiqiang. Hydrodynamics of an FLNG system in tandem offloading operation[J]. Ocean Engineering, 2013, 57: 150-162. DOI:10.1016/j.oceaneng.2012.09.015 |
| [4] |
H CHEN, T MOAN, S HAVER. Prediction of relative motion and probability of contact between FPSO and shuttle tanker in tandem offloading operation[J]. Journal of Offshore Mechanics & Arctic Engineering, 2004, 126(3): 91-100. |
| [5] |
WANG Qiang, SUN Liping, MA Shan. Time-domain analysis of FPSO tanker responses in tandem offloading operation[J]. Journal of Marine Science and Application, 2010, 9(2): 200-207. DOI:10.1007/s11804-010-9070-4 |
| [6] |
CHEN Xiao-bo. Middle-field formulation for the computation of wave-drift loads[J]. Journal of Engineering Mathematics, 2007, 59(1): 61-82. DOI:10.1007/s10665-006-9074-x |
| [7] |
Research Department of Bureau Veritas[R]. Hydrostar for Experts User Manual, 2014.
|
| [8] |
Research Department of Bureau Veritas, Ariane7 User Guide[R]. 2011.
|
| [9] |
Bureau Veritas. Guidance note quasi-dynamic analysis of mooring systems using ARIANE[J]. Software, 1998.
|