Optimal base motion to compensate for the influence of sea currents during riser re-entry

One of the difficult operations, which consists in moving the riser and placing its end relatively close to a desired position, is the re-entry operation. Complex dynamic behavior of risers under different sea conditions requires efficient modelling methods. The model used in this paper applies a modification of the segment method using joint coordinates, in which it is possible to analyze only one selected deformation while neglecting the others. This enables a very high computational efficiency of the method to be achieved. The models developed take into account the impact of the environment in which the risers work. The model is validated by comparison of the authors' own results with those presented by other researchers and the simulations are concerned both with statics and dynamics of spatial risers. The numerical effectiveness of the method presented enables it to be applied in the solution of dynamic optimization problems, one of which is presented by the example of the re-entry process. The process of moving the riser is useful in emergency situations (evacuation) when it is necessary to disconnect the riser from the wellhead and move it together with the platform. This optimization task is a 3D problem due to the sea currents acting at different angles on the riser in relation to the direction defined by beginning and final positions of the bottom end of the riser. The calculations are carried out for a hang-off riser, and the optimal motion of the base for different conditions of the sea is defined. The influence of the LMRP (Lower Marine Riser Package) on this movement is also examined.     

深水钻井立管建模与角度安全控制

从水面钻井平台与水下立管联合作业的安全角度出发，提出一种将钻井立管的力学响应限制特性引入水面平台动力定位闭环控制中的位置保持方法，实现水面钻井平台（或船舶）基于立管角度响应的动态定位。利用有限元方法建立包括立管系统质量、系统刚度、结构阻尼和水动力载荷在内的立管运动控制模型。联合水面浮体和水下立管的低频运动特性建立水面浮体运动偏移与水下立管顶端角度及末端角度的相对运动关系模型。在此基础上，设计基于立管运行响应的动力定位控位方法，实现对立管顶端角度及末端角度的安全控制。仿真结果表明，所提出的方法可行，在外界突变的环境载荷瞬时作用于水面浮体时，能更快地跟踪新的期望最优位置，保证钻井立管运行在安全界限内。     

基于顶部张紧式立管动力分析的关键参数研究

文中提出了一种适用于工作于南中国海张力腿平台(TLP)的顶张紧式立管(TTR),并利用有限元软件ABAQUS/AQUA对该立管进行了模拟分析,计算用以预测立管波浪飞溅区和立管下部区域在波浪以及平台联合作用下的影响。计算时采取时域计算方法,计算过程中考虑了几何非线性的影响。并在计算完成TTR在位动力分析的基础之上,做了一定量的参数敏感性分析,重点考察立管侧向位移,弯矩以及等效应力水平对TTR一些重要参数的敏感性。结果表明,计算顶张力TTR动力分析时,必须考虑平台侧向位移跟轴向位移的耦合作用,立管响应受其顶部TLP平台以及张力比TTF的影响较大。随着立管顶部张力比TTF的增加,立管的变形逐渐变小,但顶部张力的增加同时增加了立管内部的等效应力,立管设计时应充分考虑到这两者之间的"平衡"关系。     

南海自由站立式立管设计与分析

自由站立式立管可作为南海油气开发的立管形式。根据设计参数建立自由站立式立管系统模型,根据南海某油田环境参数校核立管各工况下的强度结果并进行疲劳分析,探讨自由站立式立管动态强度与疲劳特性并验证是否满足规范要求。分析结果表明：该立管设计满足规范中的强度与疲劳要求,且疲劳性能较好;立管等效应力最大处为立管主体顶部应力节下端;最危险浪流方向为沿立管向平台方向;立管主体应力水平受平台运动影响较小,受顶部张紧力以及海流大小影响较大;疲劳损伤最大处位于立管主体顶部,且涡激疲劳损伤大于浪致疲劳损伤。     

平台运动对深水钻井隔水管非线性动力响应的影响

文章建立了波浪、海流和船体运动等载荷综合作用下的深水钻井隔水管非线性时域动力响应模型,模型考虑了顶部张力和隔水管内部泥浆对隔水管动力响应的影响,计人了张力沿隔水管长度的变化,并采用有限差分法求解了四阶偏微分运动方程.计算了随机波浪和规则波浪下船体运动对深水钻井隔水管非线性动力响应的影响,研究分析了船体平均偏移以及波频响应和低频响应的周期和幅值对深水钻井隔水管动力响应的影响.结果表明,船体运动对深水钻井隔水管动力响应有较大的影响,在随机波浪和海流作用下,考虑船体运动时,隔水管上部最大弯矩增大了59%.船体的平均偏移主要影响底部最大弯矩.船体波频运动对上部最大弯矩有较大影响,随着波频响应的周期和幅值的变化,上部最大弯矩变化较快,而低频响应的则主要影响隔水管上部弯矩.     

柔性立管涡激振动响应轨迹特性研究

为了深入研究细长柔性立管的涡激振动响应特性,进行了柔性立管的拖曳水池试验。由拖车拖动立管产生相对来流,根据应变测试得到的应变数据,基于模态叠加法得到位移响应。试验分析前,通过数值方法先针对刚性立管的涡激振动响应轨迹特性进行了分析。紧接着,通过试验方法对柔性立管的单模态以及多模态涡激振动响应轨迹特性进行了深入的分析和讨论。通过分析发现:柔性立管在低速下具有与刚性立管类似的轨迹响应特性,均呈现经典的8字形状;柔性立管在高速下,其轨迹开始变得混乱,这主要是由位移的多模态响应特性所产生。     

考虑平台作用的顶张式立管动力分析

上层浮式平台的运动幅度比固定式平台大,对立管的影响更明显.为了解平台对顶张式立管动力的影响,提出一种动力分析方法.在一定的条件下建立立管的数学模型,采用软件Tube2D对平台和立管进行动力响应和弯曲应力分析,得到考虑平台作用的立管扶正器布置间距参数的敏感性和立管弯曲应力的变化.结果表明:对于考虑平台作用的立管模型而言,内套管与油管在立管上部更容易发生碰撞,从避碰的角度考虑,应适当减小立管上部扶正器的布置间距;各层立管的弯曲应力按弯曲刚度之比分配的结果与实际计算结果不符,应根据耦合立管模型的计算结果对各层立管进行强度校核.     

Numerical study on the effect of current profiles on vortex-induced vibrations in a top-tension riser

In this paper, numerical simulations of vortex-induced vibrations in a vertical top-tension riser with a length-to-diameter ratio of 500 using our in-house code viv-FOAM-SJTU are presented. The time-dependent hydrodynamic forces on two-dimensional strips are obtained by solving the Navier-Stokes equations, which are, in turn, integrated into a finite-element structural model to obtain the riser deflections. The riser is discretized into 80 elements with its two ends set as pinned and 20 strips are located equidistant along the risers. Flow and structure are coupled by hydrodynamic forces and structural displacements. In order to study the effects of the shear rate, of the current profiles on the vortex-induced vibrations in the riser, vibrations, with varying shear rates, in both the in-line and cross-flow directions, are simulated. In addition to the time domain analysis, spectral analysis was conducted in both the temporal and spatial domains. Multi-mode vibration characteristics were observed in the riser. The relationship between dominant vibration mode number and the shear rate of current profiles is discussed. In general, the overall vibrations in the riser pipe include contributions from several modes and each mode persists over a range of shear rates. Moreover, the results suggest that with a larger shear rate the position of the maximum in-line time-averaged displacement will move closer to the end where the largest velocity is located.     

Experimental investigation on the dynamic responses of a free-hanging water intake riser under vessel motion

A large-scale model test of a free-hanging water intake riser (WIR) is performed in an ocean basin to investigate the riser responses under vessel motion. Top end of the WIR is forced to oscillate at given vessel motion trajectories. Fiber Brag Grating (FBG) strain sensors are used to measure the WIR dynamic responses. Experimental results firstly confirms that the free-hanging WIR would experience out-of-plane vortex-induced vibrations (VIVs) under pure vessel motion even for the case with a KC number as low as 5. Meanwhile, comparison between numerical results and experimental measurements suggests a significant drag amplification by out-of-plane vessel motion-induced VIV. What’s more, further study on WIR response frequencies and cross section trajectories reveals a strong correlation between vessel motion-induced VIV and local KC number distribution, owing to the small KC number effect. The presented work provides useful references for gaining a better understanding on VIV induced by vessel motion, and for the development of future prediction models.     

取消散货船横向舱口围板下横梁的结构强度分析

以某载重180 000t散货船为例,分析横向舱口围板下横梁对于结构强度所起的作用,探讨取消散货船此处横梁的可能性。在适当的情况下,取消此横梁后仍可保证船舶强度,同时可以减轻空船质量,增加实际舱容,为船舶设计优化提供新的思路。     

平台运动下悬链线立管涡激振动响应特性分析

悬链线立管会在半潜平台运动带动下出现涡激振动现象,对立管疲劳寿命造成影响.借助海工结构物动力分析软件Orcaflex建立立管模型,使用Iwan-Blevins尾流振子模型进行数值模拟.通过与试验结果对比验证了该方法的可行性.同时建立实尺度悬链线立管模型,并考虑背景海流的影响,分析平台垂荡、纵荡运动下的立管涡激振动响应和疲劳损伤.研究表明,平台运动产生的非定常流场将引起悬链线立管发生涡激振动现象,并在背景洋流激励的基础上增大疲劳损伤,在立管结构实际工程设计时应予以关注.     

超深水钻井隔水管动力分析理论研究与数值模拟

随着钻井作业向深水(500～1 500m)和超深水(1 500m以上)发展,在交变海洋环境载荷波浪力、海流力和浮式钻井平台运动的共同作用下,隔水管的动态响应更加显著.文中探讨了隔水管侧向振动的数学模型、动态特性分析中的结构与环境载荷建模技术及其非线性动力分析方法,研究并对比了不同分析方法在计算效率、计算精度和工程适用性等方面的差异.介绍了时域内应用ABAQUS软件进行超深水钻井隔水管非确定性动力分析的算法与详细流程,算例比较了不同边界条件对深水钻井隔水管动态特性的影响.研究表明,时域非确定性分析最为精确但需要时间最长,且只能采用线性AIRY波浪理论;理论上,海流主要引起隔水管动态响应的时不变部分,但该时不变部分不等同于海流引起的隔水管静态响应,一种简化方法只将海浪与钻井船运动作为动载荷而不考虑海流对动态响应的贡献;钻井船运动和波浪载荷是隔水管动态响应分析主要的动载荷,对于超深水隔水管来说,钻井船运动是首要的动载荷,其慢漂运动对隔水管性能有重要影响,而波浪仅对隔水管局部产生作用.     

深海系缆动张力的超谐动力响应

文章考虑深水系缆的拉伸和弯曲变形,取正交坐标和切向量描述任意结点的6个广义位移,考虑缆变形的几何非线性和水动力影响,基于细长杆理论推导缆索的单元刚度矩阵,引入单元内张力一致假设,推导出12×12的非线性刚度矩阵,并与有限元软件对接。针对1800m水深平台的系缆,考虑系缆端点垂荡运动,以及端点垂荡、横荡和海流作用两种状态下系泊系统的动力响应。计算结果表明,端点运动和海流作用,引起系缆系统的二倍超谐动力响应,并且随着流速增大,系缆超谐动力响应增大,系缆的超谐动张力分量对于缆索疲劳累积损伤具有重要影响,在系缆疲劳分析时应该予以考虑。     

Numerical analysis of vortex-induced vibration on a flexible cantilever riser for deep-sea mining system

The flexible cantilever riser, as a special form of the marine riser, can be encountered in a deep-sea mining system, where the bottom of the long vertical lifting pipeline is connected with the intermediate warehouse. The main objective of this paper is to investigate the effects of the bottom weight caused by the intermediate warehouse and the flow speed on the dynamic responses of the cantilever pipeline. A quasi-3D coupling algorithm based on the discrete vortex method and finite element method is employed to calculate the unsteady hydrodynamic forces and vortex-induced vibrations of this pipeline in the time domain, respectively. We first simulate the VIV of a long flexible riser with two fixed ends in a stepped flow to validate the feasibility of the present method. Then, systematic simulations of cross-flow VIV of the cantilever riser are carried out under a wide range of bottom weights and different current speeds. The number of the vibration mode shows the decreasing tendency with the increase of the bottom weight. In a certain range of the weight, the number of the dominant mode remains unchanged, while the vibration amplitude declines with increasing weight. An amplitude jump phenomenon can be observed when the transition of the dominant mode in two contiguous mode clusters occurs. Moreover, the higher-order modes are excited with the increase of the current speed.     

Diffraction of oblique water waves by small uneven channel-bed in a two-layer fluid

Obliquely incident water wave scattering by an uneven channel-bed in the form of a small bottom undulation in a two-layer fluid is investigated within the frame work of three-dimensional linear water wave theory. The upper fluid is assumed to be bounded above by a rigid lid, while the lower one is bounded below by a bottom surface having a small deformation and the channel is unbounded in the horizontal directions. Assuming irrotational motion, perturbation technique is employed to calculate the first-order corrections to the velocity potentials in the two fluids by using Fourier transform approximately, and also to calculate the reflection and transmission coefficients in terms of integrals involving the shape function representing the bottom deformation. Consideration of a patch of sinusoidal ripples shows that the reflection coefficient is an oscillatory function of the ratio of twice the component of the wave number along x-axis and the ripple wave number. When this ratio approaches one, the theory predicts a resonant interaction between the bed and interface, and the reflection coefficient becomes a multiple of the number of ripples. High reflection of incident wave energy occurs if this number is large.     

沿底定深拖曳系统运动仿真

主要研究分析沿底定深拖曳系统的拖体机动性,通过建立沿底定深拖曳系统的数学模型,采用运动仿真计算分析的方法,验证了沿底定深拖曳系统有足够的定深能力和避障能力,但在实际应用中还需要采取辅助手段来提高避障高度.     

海底阻抗对浅海波导中椭球体的声散射研究

基于Greenberg方法推导了含有阻抗边界条件的浅水Green函数,利用边界元法研究了含阻抗的三维海洋波导与散射体相互作用的问题.分析了海底阻抗对声波与结构体相互作用时的散射特性影响,数值分析表明:含有阻抗边界条件的浅海目标的散射特性与理想海底明显不同,海底阻抗对水平断面以及垂直断面的散射场都有一定的影响,目标距离海底越近,场点散射声压受海底边界条件的影响越大;阻抗越大,散射声压越大.指向性的分布对受海底阻抗的大小比较敏感.分析浅海声散射时,海底阻抗的影响不可忽视.     

Amplitude decay and energy dissipation due to the interaction of internal solitary waves with a triangular obstacle in a two-layer fluid system: the blockage parameter

Experiments were carried out in a wave flume to examine the propagation of depression-type internal solitary waves (ISWs) over a submerged triangular-shaped obstacle. We observed the distortion and breaking of the ISW depression in a two-layer stratified fluid system as induced by the obstacle on the bottom. Physical factors, such as the amplitude of the wave, the thickness of the two layers, and the height of the submerged obstacle were utilized in our analysis of the wave–obstacle interaction. Wave breaking was determined in relation to input parameters arranged before the experimental runs. A classification scheme based on the blockage parameter ζ was proposed. The various degrees of ISW-obstacle interaction were visually classified with three schematic forms (ζ < 0.55 for a weak encounter; 0.55 < ζ < 0.7 for a moderate encounter; and 0.7 < ζ for wave breaking). Furthermore, the reflection and transmission coefficients were utilized to explore the resultant influence on changes in amplitude and energy of an ISW. The wave characteristics (e.g., amplitude-based wave reflection, energy-based reflection, amplitude-based transmission, and energy-based transmission) during the wave–obstacle interaction showed an approximately linear relation with the blockage parameter. Influenced by the submarine obstacle, the transmitted waves were found to always consist of a leading pulse (solitary wave) followed by a dispersive wave train. Moreover, due to bottom friction the energy budget of the leading pulse was reduced by a factor of at least 5% which could approach 50% in cases with strong breaking.     

Wave scattering by undulating bed topography in a two-layer ocean

The problem of wave scattering by undulating bed topography in a two-layer ocean is investigated on the basis of linear theory. In a two-layer fluid with the upper layer having a free surface, there exist two modes of waves propagating at both the free surface of the upper layer and the interface between the two layers. Due to a wave train of a particular mode incident on an obstacle which is bottom-standing on the lower layer, reflected and transmitted waves of both modes are created by the obstacle. For small undulations on the bottom of the lower layer, a perturbation method is employed to obtain first-order reflection and transmission coefficients of both modes for incident wave trains of again both modes in terms of integrals involving the bed-shape function. For sinusoidal undulations, numerical results are presented graphically to illustrate the energy transfer between the waves of different modes by the undulating bed. U. BASU was born in 1949. She is a professor in the Department Applied Mathematics, Calcutta University, India. Her current research interests include water wave problems, continuum mechanics, etc.