A finite volume solution of wave forces on submarine pipelines

In this paper, a numerical model is established for simulating the wave forces on a submarine pipeline. A set of two-dimensional Navier–Stokes equations is discretized numerically with a finite volume method in a moving mesh system. After each time step, the mesh is modified according to the changed wave surface boundary. The deffered correction second-order upwind scheme (SUDC) is adopted here to discretize the convective fluxes. The effects of the clearance between the pipeline and the seabed, water depth and wave height on wave forces are studied, respectively. The results by the numerical simulation agree well with the experimental data and theory value.     

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.     

泥沙突变理论在海底管线冲刷中应用

基于泥沙突变理论,针对海底管线冲刷下泥沙的运动特征,建立恒定流作用下的泥沙起动模式,确定希尔兹数、无量纲参数、冲刷坑深度之间相互作用的非线性方程,推导了恒定流作用下海底管线冲刷坑深度的预测公式。将相同条件下该公式的计算结果与前人的试验资料进行了对比,可发现尽管计算结果存在一定的误差,但也基本能满足对冲刷坑深度的预测和判断,从而证明了泥沙突变模型在预测海底管线冲刷坑深度中的适用性。     

Exchange flow of oil and sea-water in a ruptured submarine pipeline

The rupture of a submarine oil pipeline starts various mechanisms leading to an oil spill. Among these mechanisms the leakage of oil driven by the difference in specific gravities of oil and sea-water is difficult to estimate. A simple mathematical model has been developed and laboratory experiments have been carried out to obtain an insight into the density-driven exchange flow and to determine the leak rate. The mathematical model is predictive and takes account of the effects of friction, inclination of the pipeline, and inertia of the fluid. The experiments were done in a horizontal model pipeline. Theoretical and experimental results are in satisfactory agreement.     

海流作用下海底管道溢油扩散实验研究

采用物理模型试验研究了海流作用下海底输油管道溢油扩散特性。水槽中的原油从管道的泄漏点溢出,并从海底扩散到海面。采用非接触式光学测量技术,设计了图像采集与数据分析系统,研究原油泄漏机理和特性。通过该系统获得了溢油的轨迹、速度和海面上升时间。从定性和定量两个角度分析了水流和溢油量对溢油行为的影响。之后给出了特征物理量对各种因素的敏感性研究。溢油在水下主要以不同粒径的散状颗粒的形式分布。原油的上升过程可分为三个阶段：饱满、分散和聚集阶段。溢油量是影响油粒子上升时间的主要因素。在油粒子的上升过程中,原油的垂直速度主要受溢油量的影响,横向速度则更依赖于水流速度。不同水流与溢油量下,横向溢油速度的偏差远大于海面上升时间和垂直溢油速度的偏差。该研究可为溢油应急预测系统提供理论参考。     

海底管线地震应力影响分析

基于Biot动力固结理论建立海床-管线相互作用的计算模型。利用粘弹性人工边界,以大型有限元软件ADINA为平台对El Centro地震波作用下的海底管线的动力响应进行分析,重点讨论管线覆盖层几何形状以及海床土各向异性对海床土体的孔隙水压力和管线内应力的影响,更为实际地反映了海底管线的地震应力作用。     

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

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

波流共同作用下海底子母管线水动力的物理模型试验研究

通过物理模型试验研究海底子母管线分别在规则波加流和不规则波加流作用下的水动力特性。基于Morison方程,采用"等效直径法"分析得到子母管线拖曳力系数C_D,惯性力系数C_M和升力系数C_L(C_D+,C_L-)。试验分别考察了流速比U_c/U_w,母管与海床间隙比e/D及子母管间的相对缝隙G/D对海底子母管线水动力系数的影响。结果表明水动力系数随U_c/U_w的增大而减小;当e/D<0.5时,海床对子母管线受力的影响比较明显,C_D,C_M及C_L+均随e/D增大而减小,|C_L-|随e/D增大而增大;对子母管间的相互影响也不可忽略,C_D,C_M和|C_L-|均随G/D增大而减小,C_L+值随G/D增大而增大。     

波流共同作用下海底子母管线水动力的物理模型试验研究

通过物理模型试验研究海底子母管线分别在规则波加流和不规则波加流作用下的水动力特性。基于Morison方程,采用"等效直径法"分析得到子母管线拖曳力系数C_D,惯性力系数C_M和升力系数C_L(C_D⁺,C_L^-)。试验分别考察了流速比U_c/U_w,母管与海床间隙比e/D及子母管间的相对缝隙G/D对海底子母管线水动力系数的影响。结果表明水动力系数随U_c/U_w的增大而减小;当e/D<0.5时,海床对子母管线受力的影响比较明显,C_D,C_M及C_L⁺均随e/D增大而减小,|C_L^-|随e/D增大而增大;对子母管间的相互影响也不可忽略,C_D,C_M和|C_L^-|均随G/D增大而减小,C_L⁺值随G/D增大而增大。     

多波束测深技术在海底管道检测中的应用

为查明海底管道的位置和掩埋状态,在海底管道检测中使用多波束测深技术。介绍了多波束的系统组成、调查方法和数据处理等内容。重点对数据处理进行了论述,提出了确定海底管道位置和掩埋状态的分析方法,并对影响多波束检测效果的频率、波束角和水深等因素进行了探讨。通过分析认为,与单波束相比,多波束在海底管道检测中具有直观、高效和高分辨率的特点,具有较好的工程应用前景。     

海底管线局部冲刷机理研究综述

局部冲刷造成的悬空是海底管线失效的最大诱发因素。对均流及波浪作用下铺设在海床表面管线的局部冲刷起动机理、临界条件以及2D局部冲刷机理、流—管—土耦合机理和3D局部冲刷机理等做了分析和总结。     

A numerical study on dynamic properties of the gravity cage in combined wave-current flow

Recent work in the area of open-ocean aquaculture-system dynamics has focused either on the response of fish cages in waves or the steady drag response from ocean currents, not on them combined. In reality, however, the forces bearing on these open-ocean structures are a nonlinear, multidirectional combination of both waves and current profiles. In this paper, a numerical model has been developed to simulate the dynamic response of the gravity cage to waves combined with currents. When current flows are combined with regular waves, gravity-cage motion response (including heave, surge, and pitch) and mooring-line forces have been calculated. To examine the validity of simulated results, a series of physical model tests have been carried out. The results of our numerical simulation are all in close agreement with the experimental data.     

Experimental study of the transformation of bound long waves over a mild slope with ambient currents

The effect of currents on the variation of cross-shore bound long waves forced by bichromatic waves over a plane slope was investigated in the laboratory. In still water the growth rate of the shoaling bound long waves over the slope is proportional to h^– 5/2 (h is still-water depth). It was found that the opposing current makes the amplitudes of the bound long waves greater than those of still water for all cases. However, the amplitudes of bound long waves in a following current are reduced in the weakly modulated cases but are enhanced in the fully modulated case.     

海底泥流冲击悬跨管道拖曳力系数分析

海底滑坡是海洋油气工程最危险的地质灾害之一,直接影响海底管线运营安全。滑坡体失稳滑动过程中,由于海水掺和作用逐渐加速转变成快速滑动泥流,冲击海底悬跨管道。基于当前国际通用的计算流体动力学(CFD)方法,采用赫巴模型描述快速滑动泥流,计算分析海底滑坡冲击悬跨管道的受力特性,重点分析悬跨高度对管道法向拖曳力系数的影响。研究发现,管道法向拖曳力系数随悬跨高度的增长而增大,当达到某一悬跨高度时,管道拖曳力系数保持稳定。     

Experimental determination of wavelengths in the presence of a mean current

This paper describes a simple method for determining the wavelength of small amplitude waves under laboratory conditions where reflected wave components are present both with and without a mean current flow superimposed. It assumes a locally horizontal bed but requires no a priori assumption concerning the form of the dispersion relation with a coexisting current. Synchronous measurements of the water surface recorded along any straight line are analysed to yield Fourier coefficients at each location. It is then shown that for all practical conditions excluding a perfect standing wave, the average rate of change of wave phase in the chosen direction can be related directly to the component of incident wave number in that direction, irrespective of reflection coefficient or relative current strength. The technique has been applied to regular and bichromatic waves in a flume with an absorbing wave generator, and can also be applied in 3-D wave basins where waves and currents intersect at arbitrary angles. In combined wave–current experiments, by assuming the linear dispersion relation, it is also possible to estimate the effective current velocity.     

A numerical study of three-dimensional wave interaction with a square cylinder

In this study, a three-dimensional numerical model is used to study the wave interaction with a vertical rectangular pile. The model employs the large eddy simulation (LES) method to model the effect of small-scale turbulence. The velocity and vorticity fields around the pile are presented and discussed. The drag and inertial coefficients are calculated based on the numerical computation. The calculated coefficients are found to be in a reasonable range compared with the experimental data. Additional analyses are performed to assess the relative importance of drag and initial effects, which could be quantified by the force-related Keulegan and Carpenter (KC) number: KC_f=UT/(4πL). Here U is the maximum fluid particle velocity, T the wave period and L the length of structure aligned with the wave propagation direction. For small KC_f, the effective drag coefficient is proportional to 1/KC_f, provided the wavelength is much longer than the structural length. When wavelength is comparable to the structure dimension, the effective drag coefficient would be reduced significantly due the cancellation of forces, which has been demonstrated by numerical results.     

Separation of incident and reflected waves in wave–current flumes

A technique is developed to separate the incident and reflected waves propagating on a known current in a laboratory wave–current flume by analyzing wave records measured at two or more locations using a least squares method. It can be applied to both regular and irregular waves. To examine its performance, numerical tests are made for waves propagating on quiescent or flowing water. In some cases, to represent the signal noise and measurement error, white noise is superimposed on the numerically generated wave signal. For all the cases, good agreement is observed between target and estimation.     

Flow resistance in the coastal zone

The apparent bed roughness, the roughness value experienced by a mean flow outside the wave-boundary layer, is deduced from the physical bed roughness and the wave–current interaction mechanism. Both the physical bed roughness and the wave–current interaction are described by a (combination of) model(s). Modelling of the apparent bed roughness leads to realistic results, however, the final results are rather sensitive to the particular choice of these models. Four bed form models and two wave–current interaction models were implemented in a 1-DV flow model to calculate near-bed velocities. A comparison between measured and predicted velocities shows that reasonable results can be obtained in this way. A constant bed roughness of 0.1 m, however, leads to even better results at this site during all conditions. This can be explained by the reversed influence of the form roughness and the wave–current interaction on the apparent bed roughness value for varying wave conditions.     

后挖沟深度对深海海底管道屈曲影响数值分析

鉴于海底管道的服役水深越来越深,主要采用犁式挖沟机对预铺设于海床之上的海底管道采取后挖沟的方式将海底管道埋设于海床之下,以保护其免受不必要的损伤。针对后挖沟深度H是海底挖沟机的重要设计参数,也是影响管道悬跨的重要因素的问题,对SMD(UK)犁式挖沟机展开参数优化,确保作业过程中悬跨段管道在外部静水压力作用下,海底管道不会发生屈曲破坏。采用ABAQUS软件,分别建立了作业前和作业中两种工况下的悬跨模型,分析机械手对接触部分管道的损伤,结果显示,作业中的机械手对悬跨管道的损伤更大;同时,建立了作业中不同管径下,后挖沟深度对管道损伤的安全裕量关系曲线。进一步,结合作业中不同挖沟深度下的管跨段屈曲数值模型,对处于外部静水压力作用下的悬跨管的屈曲失效展开分析,结果显示,随着后挖沟深度的加大,不同管径下的悬跨段管道局部出现塑性压溃的临界压力值不断降低;管道外径的增大,降低了同一后挖沟深度下发生屈曲失效的压力值。最后,在后挖沟深度与外部静水压力组成的区域内,建立屈曲失效临界关系曲线,并划分出工作区和压溃区,为深海管道后挖沟埋管的施工提供工程参考。     

Deformation of monochromatic water wave trains propagating over a submerged obstacle in the presence of uniform currents

A nonlinear numerical model based on depth averaged equations and a relevant physical model have been investigated for the deformation of the water wave propagating over a submerged parabolic obstacle in the presence of uniform current. Physical and numerical modeling for wave with both following and opposing currents are done to explore the wave evolutions during passage over the submerged obstacle. A third-order Stokes dispersion relation is utilized in some cases in the computation. Separated flow zone is taken into consideration by two empirical equations obtained from the physical model testing done by the authors. Verification and validation of the numerical model by other published theoretical and experimental data are presented.