海底滑坡作用下滩海管道结构安全分析

海底滑坡作为常见的海洋地质灾害,对海洋油气工程安全产生巨大威胁。海床土体失稳引起滑坡体滑动,会对海底管道产生拖曳作用。基于计算流体动力学方法(CFD)建立海底滑坡体对管道作用的评估模型,采用H-B模型描述块状滑坡体并与试验比较验证,分析不同海床倾斜度滑坡对管道的作用并拟合表达式;研究了海底管道在滑坡作用下的力学响应,并采用极限状态方法开展海底滑坡作用下管道结构极限安全分析,探讨了管道埋地状态时的极限安全界限,建立滑坡作用下管道结构安全分析方法。研究表明:滑坡对管道作用力与海床倾角呈现正相关,而覆土层厚度对作用力影响较小;随着不排水抗剪强度的减小,允许的滑坡宽度和速度均增加,表明土体不排水抗剪强度与引起的拖曳力呈正相关;滑坡土体宽度对极限安全速度影响较大。     

Numerical analysis of the impact forces exerted by submarine landslides on pipelines

Submarine pipelines that transport crude oil and natural gas are often in a complex marine geological environment and may become unstable and fail upon impact by submarine landslides. Previous research has mostly focused on the impact forces exerted by submarine landslides on suspended pipelines, but the impact of submarine landslides on pipelines laid on the seafloor at various impact angles, θ, have been relatively infrequently discussed, and the effects of suspended height, H, on the impact forces exerted by submarine landslides on pipelines have not been thoroughly investigated. In this study, based on the Herschel–Bulkley model, the impact forces exerted by a submarine landslide on laid-on or suspended pipelines at various impact angles θ were simulated using the computational fluid dynamics (CFD) approach. Equations for calculating the axial and normal drag coefficients of a submarine pipeline were proposed. The CFD numerical simulation results were rearranged based on the soil mechanics approach. By comparing the parameters, an essentially corresponding relationship was found between the soil mechanics and CFD approaches when the equations were used to calculate the impact forces exerted by a submarine landslide on a pipeline. In addition, a semi-analytical expression for the failure envelope was provided. Furthermore, the effects of H on the forces on a pipeline were discussed, and an equation for calculating the acting forces on a pipeline along the flow direction of a submarine landslide that comprehensively accounts for the effects of θ and H was proposed. The lift force was discussed preliminarily and the results provide a basis for further investigation. The achievement of this study is applicable for selecting locations of submarine pipeline routes and for designing submarine pipelines.     

海底滑坡对海洋单桩冲击力试验研究

为了研究海底滑坡对海洋单桩的冲击力大小,首先通过调整高岭土、粉砂的不同含量,得到不同流变特性、不同密度的碎屑流,采用Herschel-Bulkley模型和幂率模型对流体流变性质进行描述;随后利用自制海底滑坡模型槽,模拟碎屑流在不同流速和黏度下对模型桩的冲击;并结合流体力学理论,建立阻力系数与非牛顿流体雷诺数之间关系表达式。试验数据表明:碎屑流黏度和流速是影响海底滑坡冲击力的主要因素,海底滑坡冲击力随着泥浆黏度和流速的增加而增大。同时,考虑碎屑流剪切稀释特性,得到管桩阻力系数随雷诺数变化的拟合公式,为海洋桩基础设计提供参考。     

荔湾3-1气田潜在海底滑坡滑动距离预测及动态特性分析

A large number of submarine landslides with different scales have been identified in the canyon area of the submarine pipeline route of Liwan 3-1 gas field. There is still much chance that submarine sl...     

Wave-induced uplift force on a submarine pipeline buried in a compressible seabed

A two-dimensional finite-element simulation of the wave-induced hydrodynamic uplift force acting on a submarine pipeline buried in sandy seabed sediments subject to continuous loading of sinusoidal surface waves is presented. Neglecting inertia forces, a linear-elastic stress-strain relationship for the soil and Darcy's law for the flow of pore fluid are assumed. The model takes into account the compressibility of both components (i.e., pore fluid and soil skeleton) of the two-phase medium.The results of numerical analysis are presented and discussed with respect to soil and pore fluid parameters where special attention is paid to the question of soil saturation conditions. The meaning of the results is also related to surface wave conditions. As a general conclusion, the practical, engineering recommendation is given in order to make a realistic, safe and economic estimation of the wave-induced uplift force acting on a buried submarine pipeline.     

海底管线稳定性影响因素分析

实现管线的稳定性是埋设海底管线并使其正常运营的前提。分析了波浪、海流、管线周围土体性质及液化发生等方面对海底管线稳定性的影响，并提出了相应的对策建议．     

Effect of Rate-Dependent Soil Strength on Cylinders Penetrating Into Soft Clay

This paper presents a predictive model for undrained penetration of cylinders into soft seafloor soils. The penetration depth will depend on the velocity of the cylinder as it touches down at the seafloor, and the net deceleration of the cylinder as it is acted on by forces of self-weight, soil buoyancy, and soil-shearing resistance. The soil-shearing resistance force increases as a function of penetration depth and, due to the dependence of undrained shear strength on strain rate, penetration velocity. This paper presents finite element (FE) simulations that quantify both effects and form the basis of a simplified soil-resisting force model. Strain-rate effects are modeled within a framework of rate-dependent plasticity, with shearing resistance increasing semilogarithmically with increasing strain rate above a certain threshold strain rate. With all forces acting on the cylinder, estimated penetration depths are predicted from simple equations of motion for a single particle. Comparisons to laboratory tests involving penetration of cylinders into soft reconstituted marine clay show reasonable agreement between model predictions and measurements.     

An improved analytical approach for simulating the lateral kinematic distress of deepwater offshore pipelines

Offshore pipelines are critical infrastructures and any possible damage may have devastating financial and environmental consequences. Earthquake-related geohazards (such as strong ground motion, active seismic faults, submarine landslides and debris flows) consist crucial threats that an offshore pipeline has to overcome. The main aim of the current study is to examine analytically a seabed-laid offshore pipeline subjected to a lateral kinematic distress due to a submarine landslide or a debris flow. Extra emphasis is given on the impact of pipe-soil interaction on the pipe response, by the realistic representation of the soil resistance via a tri-linear model. Firstly, the proposed analytical model is validated with a numerical model utilizing the finite-element method. Subsequently, various combinations of soil parameters and loading conditions that affect the examined problem are investigated with realistic input data taken from the offshore section of the high-pressure natural-gas pipeline TAP (Trans Adriatic Pipeline) in the Adriatic Sea. Finally, useful conclusions are drawn regarding the applicability and the efficiency of the proposed approach.     

Parametric study of the wave-induced residual liquefaction around an embedded pipeline

Wave-induced liquefaction in a porous seabed around submarine pipeline may cause catastrophic consequences such as large horizontal displacements of pipelines on the seabed, sinking or floatation of buried pipelines. Most previous studies in relation to the wave and seabed interactions with embedded pipeline dealt with the wave-induced instaneous seabed response and possible resulting momentary liquefaction (where the soil is liquefied instantaneously during the passage of a wave trough), using theory of poro-elasticity. Studies for the interactions between a buried pipeline and a soil undergoing build-up of pore pressure and residual liquefaction have been comparatively rare. In this paper, this complicated process was investigated by using a new developed integrated numerical model with RANS (Reynolds averaged Navier–Stokes) equations used for governing the incompressible flow in the wave field and Biot consolidation equations used for linking the solid–pore fluid interactions in a porous seabed with embedded pipeline. Regarding the wave-induced residual soil response, a two-dimensional poro-elastoplastic solution with the new definition of the source term was developed, where the pre-consolidation analysis of seabed foundation under gravitational forces including the body forces of a pipeline was incorporated. The proposed numerical model was verified with laboratory experiment to demonstrate its accuracy and effectiveness. The numerical results indicate that residual liquefaction is more likely to occur in the vicinity of the pipeline compared to that in the far-field. The inclusion of body forces of a pipeline in the pre-consolidation analysis of seabed foundation significantly affects the potential for residual liquefaction in the vicinity of the pipeline, especially for a shallow-embedded case. Parametric studies reveal that the gradients of maximum liquefaction depth with various wave and soil characteristics become steeper as pipeline burial depth decreases.     

海底水平管线回流区水流特性研究

针对管线绕流剪切层区的水流特征,利用射流分析的类似方法,通过简化控制方程,从理论上探讨了剪切层区的时均速度分布规律。物理试验中采用ADV测速仪对回流区流场进行了测量,对各流区流速和紊动强度变化趋势进行了讨论,并将试验结果与理论推导进行对比,二者基本吻合,为进一步研究海底管线防护措施提供了理论依据。     

高雷诺数下近壁圆柱绕流数值模拟与分析

近壁圆柱绕流问题在海底悬跨管道的研究中具有重要的意义。在绕流阻力、升力以及海底土壤的耦合作用下,海底管道所发生的移位、悬跨等现象对于海底管道的安全运行构成了很大的威胁。正确预测各种绕流条件下管流之间的作用力是保证油气管道安全的首要任务。海底管道在极端海洋环境条件下的管、流相互作用为高雷诺数绕流问题,处于高雷诺数下的绕流模拟比处于低雷诺数下的绕流模拟要复杂很多,它需要更精细的网格以及合适的湍流模型。此文对处于悬跨状态下的海底管道进行数值研究,给出不同间隙比下海流绕流海底管道的流场结构形态,分析了间隙比对绕流阻力和绕流升力的影响,为进一步研究海底悬跨管道的受力和变形提供载荷边界数据。     

Euler-Euler two-phase flow simulation of tunnel erosion beneath marine pipelines

In this study an Euler-Euler two-phase model was developed to investigate the tunnel erosion beneath a submarine pipeline exposed to unidirectional flow. Both of the fluid and sediment phases were described via the Navier-Stokes equations, i.e. the model was implemented using time-averaged continuity and momentum equations for the fluid and sediment phases and a modified k−ε turbulence closure for the fluid phase. The fluid and sediment phases were coupled by considering the drag and lift interaction forces. The model was employed to simulate the tunnel erosion around the pipeline laid on an erodible bed. Comparison between the numerical result and experimental measurement confirms that the numerical model successfully predicts the bed profile and velocity field during the tunnel erosion. It is evident that the sediments are transported as the sheet-flow mode in the tunnel erosion stage. Also the transport rate under the pipe increases rapidly at the early stage and then reduces gradually at the end of the tunnel erosion beneath pipelines.     

往复流作用下海底管道局部冲刷机制数值模拟

基于开源的计算流体力学模式REEF3D,建立了海底管道局部冲刷水槽数值模型,在验证单向流实验结果的基础上,进一步对往复流作用下的海底管道局部冲刷机制进行了研究,并作对比分析。研究表明,总体上往复流对管道所在海床的局部冲刷规模比单向流弱。当流向改变后,原先下游的堆积区转变成上游,优先受到冲刷,并填充到管道下方的冲刷坑,同时使水流在管道下方的作用减弱。这种回填过程使短周期下的往复流作用需要更长的冲淤平衡时间。在回填与冲刷的共同作用下,上下游的冲刷坑坡度会因流向变化而变化。     

Interaction between submarine landslides and suspended pipelines with a streamlined contour

Recently, the security and stability of submarine pipelines have attracted much attention in ocean engineering. In this paper, pipelines with a streamlined contour (wedge, airfoil, double-ellipse, and arc-angle hexagon) are designed in hopes of defending against the impact of submarine landslides, and the computational fluid dynamics (CFD) approach is used to investigate the interaction between submarine landslides and streamlined pipelines. The results show that the peak interactional force is more representative of the hazard level of pipelines imposed by submarine landslides. It is also found that the streamlined pipelines possess a significant advantage in reducing the drag force and lift force of landslide–pipeline interaction with a maximum lessening percentage of 66.32 and 40.17%, compared with a conventional circular pipeline. In addition, the influence of applying streamlined pipelines to engineering is briefly discussed, and the empirical equation for estimating the drag force and lift force of streamlined pipelines induced by landslides is recommended based on the numerical test results.     

Prediction of hydrodynamic forces on submarine pipelines using an improved Wake II Model

The hydrodynamic force model for prediction of forces on submarine pipelines as described includes flow history effect (wake effects) and time dependence in the force coefficients. The wake velocity correction is derived by using a closed-form solution to the linearized Navier–Stokes equations for oscillatory flow. This is achieved by assuming that the eddy viscosity in the wake is only time dependent and of a harmonic sinusoidal form. The forces predicted by the new Wake (Wake II) Model have been compared to Exxon Production Research Company Wake Model in terms of time histories (force shape) and magnitudes of peak forces. Overall, the model predictions by the Wake II Model are satisfactory and represent a substantial improvement over the predictions of the conventional models. The conventional force models representing adaptations of Morison's equation with ambient velocity and constant coefficients give predictions that are in poor agreement with the measurements especially for the lift force component. The Wake II Force Model can be used for submarine pipeline on-bottom stability design calculations for regular waves with various pipe diameters.     

Prediction of hydrodynamic forces on submarine pipelines using an improved Wake II Model

The hydrodynamic force model for prediction of forces on submarine pipelines as described includes flow history effect (wake effects) and time dependence in the force coefficients. The wake velocity correction is derived by using a closed-form solution to the linearized Navier–Stokes equations for oscillatory flow. This is achieved by assuming that the eddy viscosity in the wake is only time dependent and of a harmonic sinusoidal form. The forces predicted by the new Wake (Wake II) Model have been compared to Exxon Production Research Company Wake Model in terms of time histories (force shape) and magnitudes of peak forces. Overall, the model predictions by the Wake II Model are satisfactory and represent a substantial improvement over the predictions of the conventional models. The conventional force models representing adaptations of Morison's equation with ambient velocity and constant coefficients give predictions that are in poor agreement with the measurements especially for the lift force component. The Wake II Force Model can be used for submarine pipeline on-bottom stability design calculations for regular waves with various pipe diameters.     

The influence of scale,slope and channel geometry on the flow dynamics of submarine channels

Submarine channels are major morphological features of the sea floor and are important in the transport of sediment to the deep ocean. Although much is known concerning the large-scale distribution of sediment within and surrounding submarine channels, there is little understanding of the fluid dynamic processes that control this sedimentation. Direct measurement of flow velocities and concentrations has proved to be extremely difficult within submarine channels, with the resultant paucity of direct observations making physical laboratory modelling a critical technique for examining the processes that operate in, and control, submarine channel development.Recent experimental and numerical studies have proposed a new model of secondary circulation within submarine channel bends, characterised by a reversal in the orientation of the secondary circulation cell relative to that found in meandering rivers. This new paradigm for submarine channels thus predicts basal flow from the inside to the outside of the bend at a bend apex, with an upper return flow directed towards the inner bend. The reversal in orientation of the secondary flow cell has been linked to the vertical distribution of downstream velocity and associated changes in centrifugal and pressure gradient forces. However, previous work has additionally proposed that shearing of the within-channel flow by overbank flow may also generate secondary flow reversal.This study assesses the applicability of the proposed submarine bend flow model against a range of key channel parameters. We demonstrate that the sense of secondary circulation is the same for all experimental conditions, strongly supporting the new model of secondary flow in submarine channels. Furthermore, investigation of overbank shear induced secondary circulation confirms for the first time that this mechanism can occur, and identifies the channel styles most likely to exhibit this effect. Such shear-induced circulation is, however, shown to be a secondary mechanism, with the vertical distribution of downstream velocity the principal mechanism. In certain channel configurations, the two mechanisms may act to augment one another.     

内流对海底管线涡致振动与疲劳寿命的影响

在复杂的海洋环境条件下，管道的动力学特性受到内外部流体的作用而呈现出新的特点。本文通过对处在海底平稳流动作用下的悬跨管道的涡致振动进行分析，特别地考虑到管道内部流动的作用，给出了内流速度对管道响应幅度的影响，进而指出其对管道疲劳的意义。     

海底管道管跨段在内外流体作用下的竖向动力特性研究

根据在复杂的海洋环境条件下,管道的动力特性受到内外流体的综合作用影响,呈现与陆地管道不同的特点。研究了结构受到外界流体所产生的涡激作用,同时考虑管内恒定流的影响,利用有限元方法对海底管道管跨段竖向振动的微分方程进行求解。在此基础上得出了管跨段在不同外流流速情况下,内流流速与结构动力响应幅值的关系及其对管跨段振动频率的影响,并进一步探讨了管跨失稳时的极限管内流速与跨长的关系,这些结论对于海底管道设计具有一定实际指导意义。     

考虑流固耦合时的海底管道悬跨段非线性动力分析

通过对管道的涡激振动试验,提出了考虑流固耦合的非线性涡激升力表达式,并用该式进行了海底管线悬跨段非线性动力响应时程分析。对考虑流固耦合与未考虑流固耦合情况下得到的管道动力响应时程进行对比,算例表明:当管外流场流速与管道顺流向振动速度值较接近时,不考虑流固耦合时的计算结果明显小于考虑流固耦合时的计算结果。分析认为,在管外流场流速与管道顺流向振动速度值较接近的情况下,管道的涡激振动计算宜采用非线性涡激力模型。