海底管道整体式止屈器性能研究与穿越压力预测

整体式止屈器可有效抑制海底管道屈曲传播并避免管道的大规模屈曲失效。为探究整体式止屈器的止屈特性,预测其穿越压力,利用ABAQUS软件建立有限元模型,分析了整体式止屈器厚度、长度、径厚比等参数对穿越压力和穿越形式的影响机制,揭示了整体式止屈器的止屈特性。研究结果表明:在平行穿越阶段,增加止屈器厚度、有效长度可明显提升穿越压力;而在垂直穿越阶段,穿越压力增长速率明显放缓。结合有限元模拟结果,拟合得到了整体式止屈器在不同穿越形式下的穿越压力预测公式,其相较于已有经验公式准确性更高,对整体式止屈器的设计和建造具有重要的指导意义。     

整体式止屈器数值模拟方法与优化研究

海底管道是海洋石油主要的运输方式,但海底管道在外部静水压力作用下很可能发生局部屈曲压溃,造成灾难性后果。安装止屈器可有效阻止屈曲传播,保证管道安全,整体式止屈器在实际中应用最为广泛。整体式止屈器存在通过减小止屈器体积以节约成本进行优化的可能性,因此,在保证止屈器止屈效果的前提下,可通过减小止屈器的体积,以达到减轻重量、节约成本、减少安装难度等多种目的。研究过程中使用ABAQUS进行模型的仿真模拟,比较了广义弧长法的优缺点,决定选用更方便的静水流体单元法,计算得到不同止屈器的穿越压力以及相应的体积,对比各止屈器的穿越压力变化和体积变化,实现结构优化的目标。基于以上研究,可以得出减少止屈器体积的结构优化方案具有可行性的结论。     

Buckle propagation in pipelines with non-uniform thickness

Finite element solutions for the steady-state buckle propagation pressure in a pipeline with non-uniform thickness are given. The results are useful in finding buckle propagation pressures in corroded pipelines. It has been found that when corrosion is equal to or less than 10% of the original pipeline thickness, the pipeline collapses in an overall shell buckling mode; otherwise, the pipeline experiences local buckling. The propagation pressure decreases with both the thickness and angular extent of the reduced section, but the rate of decrease with thickness reduction is almost independent of the angular extent of corrosion when it is greater than 90°.     

Numerical analysis of cross-section ovalization in the deep-sea pipeline lateral buckling process

Abstract

The deep-water pipeline is the main means of transportation in offshore oil and gas development engineering. The deep-water pipeline may incur lateral global buckling due to the high temperature and pressure that are applied on the pipeline to ensure the contents’ liquidity. With the increasing operating water depth, a higher temperature and pressure are applied to the pipeline, causing large lateral deformation and a large bending moment. Due to the inhomogeneous distribution of the bending moment on the cross-section, different points on the cross-section will deform differently. This kind of deformation causes the cross-section to turn into an oval ring. The cross-section ovalization caused by global buckling was rarely analyzed in former engineering practice since the load is relatively low. With the increase in operation water depth and operation load, the ovality caused by global buckling is noticeable. This article analyzed cross-section ovalization caused by pipeline lateral global buckling with a numerical simulation method. The pipelines with different initial cross-section shapes were simulated, and the influence of several impact factors, including load, pipeline and soil factors on the ovality of the cross-section, were analyzed. The results show that the initial cross-section shape type has little effect on the pipeline ovalization pattern. The initial ovality of the pipeline with an oval ring cross-section shape has little influence on the residual ovality. Among all the factors analyzed in this paper, the pressure difference is the primary factor that should be considered in a pipeline ovalization check.     

A modified numerical calculation method of collapse pressure for thick-walled offshore pipelines

Considering the shear deformation and thickness stretching of large deformation, a modified numerical calculation method based on the thick shell theory is established to determine the collapse pressure of thick-walled pipes. Verification experiments are conducted on ten pipe specimens in hyperbaric chambers. The good agreement between experimental results and numerical predictions shows the validity and reliability of the new numerical calculation method. Combining DNV specification, the characteristic collapse pressure is also calculated for comparison. The difference between experimental results and DNV calculations illustrates the latter one is much conservative in predicting collapse pressure for thick-walled pipes. Sensitivity analysis on manufacturing imperfections and material properties is investigated for pipes with different D/t ratios. Thick-walled pipes are easier to be affected by initial ovality, residual stress and hardening factor. Based on the stress distribution at the moment of collapse, a novel discovery is found that the collapse pressure of thick-walled pipes is dominated by material plastic behavior.     

A numerical model for ultimate soil resistance to an untrenched pipeline under ocean currents

One of the main concerns for pipeline on-bottom stability design is to properly predict ultimate soil resistance in severe ocean environments.A plane-strain finite element model is proposed to investigate the ultimate soil resistance to the partially-embedded pipeline under the action of ocean currents.Two typical end-constraints of the submarine pipelines are examined,i.e.freely-laid pipes and anti-rolling pipes.The proposed numerical model is verified with the existing mechanical-actuator experiments.The magnitude of lateral-soil-resistance coefficient for the examined anti-rolling pipes is much larger than that for the freely-laid pipes,indicating that the end-constraint condition significantly affects the lateral stability of the untrenched pipeline under ocean currents.The parametric study indicates that,the variation of lateral-soil-resistance coefficient with the dimensionless submerged weight of pipe is affected greatly by the angle of internal friction of soil,the pipe-soil friction coefficient,etc.     

The structure of shelf and gulf stream motions in the Georgia bight

The structure of the current and temperature fields along 30°N over the mid-shelf and western Blake Plateau in the South Atlantic Bight has been investigated by combining two moored instrument experiments in the summer of 1981. The shelf moorings were part of the second Georgia Bight Experiment (GABEX-II) and the Gulf Stream mooring data on the Blake Plateau have been described by LEE and WADDELL (1983). Empirical Orthogonal Functions (EOF) in the frequency domain are used to extract shelf and Gulf Stream coherent current and temperature fluctuations in the two- to 14-day period band. Three modes are found, of which the first two are interpreted as Gulf Stream meander and frontal eddy circulations. The difference between them is chiefly in the shelf motions; the first mode is primarily restricted to the shelf edge, whereas the second mode penetrates to the 40m isobath. The third mode dominates at mid-shelf and is the only mode that shows strong coherence with the windstress and local sea-level fluctuations. The relationship of the modes to the occurrence of mid- and inner-shelf cold sub-surface intrusions, generated by shelf-edge Gulf Stream frontal eddies, is examined. All three modes are found to play a role in the initiation, growth and decay of these structures.     

刚性连接双层海底管道高温侧向屈曲分析方法研究

以惠州25-3/1双层海底管道为例,提出高温荷载下刚性连接双层海底管道的侧向屈曲分析方法,并研究该类管道的侧向屈曲特性.将解析分析与有限元技术相结合,在提出刚性连接双层海底管道屈曲分析方法的同时,发展了Hobbs公式的应用范围,使其成为屈曲分析有限元建模指导工具.以分析不稳定平衡问题的改进Riks方法被证明为高温海底管道屈曲分析的有效算法,得到的弧长-载荷比例因子曲线能够充分地描述双层管道的屈曲失稳过程.惠州管道的分析结果表明,刚性连接双层海底管道的整体热稳定性较高,在海床上一般不容易发生侧向屈曲,但热荷载作用下内管的后屈曲将有可能引发较高的应力集中,对此需要充分校核;另外环空间隙的减小有利于提高内管的临界屈曲载荷,从而提升管道的整体性能.     

Global buckling of pipelines in the vertical plane with a soft seabed

Offshore pipelines are usually buried to avoid damage from fishing activities and to provide thermal insulation. Provided that the buried pipelines are sufficiently confined in the lateral direction by the passive resistance of the trench walls, they may be subject to vertical buckling caused by a rise in temperature. Vertical buckling is usually called upheaval buckling because the heated pipeline is assumed to move upwards conventionally. However, the seabed may be very soft, especially where a pockmark or abyssal ooze appears. Consequently, under thermal compressive force, the pipeline may buckle downward and penetrate into the seabed because the downward soil resistance is small. In this study, we extended an analytical solution for vertical pipeline buckling on a rigid seabed to a soft seabed, and the effects of soil resistance on pipeline stability, buckling mode and amplitude are illustrated and analyzed.     

海洋纤维增强热塑性立管极限承载拉力预测方法

海洋新型纤维增强热塑性立管因其可盘卷、耐腐蚀、耐疲劳和轻质化等优点,在深水油气开发中应用前景十分广阔。热塑性立管具有复合材料的各向异性、受力耦合效应及复杂的本构关系,且承受浮体运动和复杂海洋环境载荷,其失效模式尚未明确。针对轴对称载荷作用下纤维增强热塑性立管极限承载力问题,进行热塑性管稳态热传导和热应力的理论推导,求解了稳态温度和应力分布,首次给出了在任意温度载荷作用下管体径向位移的解析解,并直接求解其径向、轴向、环向和剪切应力。采用各向同性层Von Mises和各向异性层最大应力（Max Stress）准则或Tsai-Hill准则判定热塑性管的失效,基于应力分布、失效准则和二分法计算了热塑性管的极限载荷。温度载荷、纤维铺设角度和径厚比对管道的应力分布影响显著。不同温度载荷会改变失效指数沿径向的变化趋势,增大轴向拉力将增大热塑性管的失效指数,选用不同的失效准则在管体失效判定上存在一定的差异。热塑性管温度越低、纤维铺设角越小及径厚比越大,管道对轴向拉伸载荷的承载能力越强。     

Instantaneous multi-mode identification and analysis of vortex-induced vibration via a mode decomposition method

The dynamic characteristics of marine risers/pipes often present serried modes with various frequencies due to high levels of structural flexibility and slenderness, especially when the flow velocity is non-uniformly distributed along the span. Therefore, the vortex-induced vibration (hence VIV) for slender risers/pipes is usually characterized by multi-mode motions. In this paper, by means of a newly developed empirical mode decomposition (EMD) method which contributes to more efficient instantaneous multi-mode identification and analysis, new characteristics of a multi-mode “lock-in” vibration process of a large-scale flexible pipe subject to shear flow were discussed. Because the two-degree vibration along the span can be analyzed simultaneously, the effects of multi-mode VIV were investigated systematically. From the given illustrative examples, it was found that the vibration energy diffusion between the fluid and the structure, and among the participating modes, may be repeatable and reversible, or even irreversible, which causes VIV to be highly intricate. The coexistence of multiple modes, energy transfer, and mode switching/jump is observed when the reduced velocity is relatively high. The multi-dominant mode phenomenon is also found in both cross-flow (CF) and in-line (IL) VIVs. Energy transfers between the CF and IL directions occasionally occur, and CF VIV is apt to dominate the vibration process, because it is superior to IL VIV with the increment of the reduced velocity.     

海底管线整体屈曲过程中塑性变形研究

海底管线是海洋油气工程中主要的输送手段。在工作状态下,受高温高压的影响,深海管线可能会发生水平向整体屈曲。随海洋油气作业水深的增大,施加于管线的温度和压强也逐渐增加,导致管线产生较大的屈曲位移和截面应力,使得截面产生塑性应变。本文采用数值模拟方法,对海底管线整体屈曲过程中塑性区的分布及其与整体屈曲影响因素的关联性、塑性变形对水平向变形的影响和塑性变形造成的截面椭圆度的变化规律进行分析,研究塑性变形对整体屈曲过程影响的规律。     

双拱初始缺陷海底管线水平向整体屈曲数值模拟分析

为了研究具有双拱反对称初始缺陷海底管线的整体屈曲特性,采用模态分析法将最可能出现的缺陷形态引入数值分析模型中。针对管线在高温高压作用下发生整体屈曲的动态变形特征,运用显式动力数值模拟方法研究了管线整体屈曲过程中水平向变形与轴向变形随温度和内压的变化规律,建立了在整体屈曲过程中屈曲管段与滑动管段轴力的变化过程与初始缺陷形态的关系。将数值模拟结果同经典解析解和室内模型实验结果进行对比,验证了本方法的可靠性。工程算例的分析结果表明,管线整体屈曲的发生是一个由低阶向高阶发展的过程,具有双拱缺陷的管线首先发生二阶模态的整体屈曲,而后过渡到四阶模态;管线整体屈曲的变形包括屈曲段的水平向变形和滑动管段的轴向缩进,其中水平变形释放了管壁内的轴力,轴力的释放量随初始缺陷尖锐程度的降低而增大;轴向缩进变形由于受到地基土的摩阻力使滑动管段内的轴力发生累积,轴力的累积量随初始缺陷的尖锐程度的降低而增加。以上研究成果对指导实际工程具有现实意义。     

Dynamic interaction between two parallel submarine pipelines considered vortex-induced vibration and local scour

Abstract

The mechanism of local scour under two vibrating pipelines is investigated numerically in this research. A sediment scour model is adopted to estimate the motion of sediment. The general moving objects model, which is dynamically coupled with fluid flow, is set up to simulate the vortex-induced vibration (VIV) of the pipeline. The sediment scour model and pipeline vibration model are verified with the previous experimental results and show good agreement. Then, the coupling effects between the pipeline vibration and the local scour are investigated numerically. The effects of G/D (the ratio of the distance between the two pipelines to the diameter of the pipelines) on the local scour and the VIV of the pipeline are examined. The results indicate that the maximum scour depth under the vibrating pipelines is much larger than the scour depth under the fixed pipelines. Due to the shadowing effect of the upstream pipeline, the maximum scour depth under the upstream pipeline is deeper than that under the downstream pipeline. The pipeline vibration magnitude is closely related to the strength of the vortex that sheds behind the pipeline. The effect of G/D on the shape and strength of the vortices that shed behind the pipelines is significant.     

On the external pressure capacity of deepwater sandwich pipes with inter-layer adhesion conditions

Sandwich pipes (SP), a composite structure consisting of two concentric steel tubes and a polymeric or cement-based core has been proposed to be well-insulated and withstand high installation and operational loads in deepwater oil and gas field. In the paper, the post-buckling responses and pressure capacity of sandwich pipes filled with the solid polypropylene annular were investigated. The degree of the inter-layer adhesion condition between the core layer and the surrounding pipes was modeled by the contact surfaces adopting different maximum shear strength values to allow the relative displacement between the layers. The effects such as inter-layer adhesion interactions, thickness-to-radius ratios, the core thickness, the material parameters, the relative initial ovality directions and the inelastic anisotropy on the collapse pressure of SPs were discussed. More than 2000 finite element (FE) models of the sandwich pipes with practical configurations were constructed and analyzed using the programming language Python within the FE software package ABAQUS. Based on the FE results, a simplified equation was developed and recommended for the collapse pressure calculation of sandwich pipes with the polypropylene annular.     

Water wave scattering by bottom undulations in an ice-covered two-layer fluid

The scattering of plane surface waves by bottom undulations in an ice-covered ocean modelled as a two-layer fluid consisting of a layer of fresh water of lesser density above a deep layer of salt water, is investigated here by using a simplified perturbation analysis. In such a two-layer fluid there exist waves of two different modes, one with higher mode propagates along the interface and the other with lower mode propagates along the ice-cover. An incident wave of a particular mode gets reflected and transmitted by the bottom undulations into waves of both the modes so that transfer of wave energy from one mode to another takes place. The first-order reflection and transmission coefficients of two different modes are obtained due to incident waves of again two different modes by employing Fourier transform technique in the mathematical analysis. For sinusoidal bottom topography these coefficients are depicted graphically against the wavenumber. These figures show how the transfer of energy from one mode to another takes place.     

Higher harmonics induced by a submerged horizontal plate and a submerged rectangular step in a wave flume

The decomposition of a monochromatic wave over a submerged plate is investigated experimentally in a wave flume. Bound and free higher harmonic modes propagating upstream and downstream the structure are discriminated by means of moving resistive probes. The first-order analysis shows a resonant behaviour linked to the ratio of the plate's width and the fundamental mode wavelength over the plate. The second-order analysis shows an energy transfer from the fundamental mode towards free harmonics propagating downstream the structure. This transfer is linked to the ratio of the width of the plate and the bound harmonic wavelength over the plate. We also performed experiments with a submerged step to compare the efficiency of both structures. The submerged plate is shown to be a more efficient breakwater than the step, at the first as well as the second-order.     

The dispersion relations for surface plasmon in a nonlinear—metal—nonlinear dielectric structure

In the asymmetric and symmetric nonlinear--metal--nonlinear dielectric structures, this paper studies the analytic dispersion relation for surface plasmon in a system consisting of a thin metallic film covered on two sides media of intensity-dependent refractive indexes by applying a generalised first integral approach. Especially in the symmetric waveguide structure, two possible modes can exist: the odd mode and the even mode. The dispersion relations of the two modes are obtained. Due to the nonlinear dielectric, the squared magnitude of the electric field at the interface appears and alters the dispersion relations. Numerical results are compared to those from a certain approximate treatment.     

管线钢在海泥中的腐蚀行为研究

利用线性极化电阻、极化曲线和电化学阻抗测量技术,研究了API Spec X60两种管线钢ERWC电阻焊管)和SML(无缝钢管)在取自中国东海长江口地区的不同类型海泥中的电化学腐蚀行为。结果表明,SML管线钢在海泥中的耐蚀性明显优于ERW管线钢,管线钢在海泥中的腐蚀行为主要受阴极去极化剂一氧的扩散控制。     

Analysis of hydrobiological pattern in the Bizerte lagoon (Tunisia)

This study deals with the use of a multivariate analysis method in order to understand the functioning of the Bizerte lagoon ecosystem. A set of hydrobiological parameters was collected during 10 monthly campaigns in 2004. A variant of the EOF (Empirical Orthogonal Function) technique is used to examine hydrobiological variability modes in this lagoon. The permanent features characterising the lagoon are obtained by applying the EOF technique to data after removing the annual cycle. Two major modes were obtained. The first one shows a contrast between the northern sector of the lagoon, influenced by the Mediterranean Sea, and the southern one influenced by continental waters. This mode is mostly prominent for temperature and salinity. The second mode shows a contrast between shallow waters along the coastal zones and the deep central area of the lagoon. This contrast is more pronounced towards the southwestern edge of the lagoon, a region of important freshwater input from the Tinja channel, which drifts towards the downstream area of Menzel Bourguiba. This second mode may also correspond to northward water circulation under southerly wind blowing frequently in summer. The spatial and temporal variability in the lagoon is highlighted by applying the EOF to the raw data without removing the annual cycle. A weight series is obtained for each variable by extracting the corresponding data and then projecting them on the corresponding eigenvector. Using this technique, only one mode was found dominant with more than 52% of the total variance. It was also found that the north–south contrast of temperature, salinity and phosphorous and the coast-centre contrast of Chlorophyll a are enhanced in summer. Other variables such as ammonia and nitrate also show significant annual cycles but with local patterns.