波浪作用下孔隙海床-管线动力相互作用分析

波浪作用下海床中的孔隙水压力与有效应力是影响海底管线稳定性的主要因素。然而,在目前的海床响应分析中一般将管线假定为刚性,并不能合理地考虑海床与管线之间的相互作用效应,同时也没有考虑土体和管线加速度对海床动力响应的惯性影响,从而无法确定由此所引起的管线内应力。为此考虑管线的柔性,分别采用饱和孔隙介质的Biot动力固结理论和弹性动力学理论列出了海床与管线的控制方程,进而采用摩擦接触理论考虑海床与管线之间的相互作用效应,基于有限元方法建立了海床-管线相互作用的计算模型及其数值算法。通过变动参数对比计算讨论了管线几何尺寸、海床土性参数对波浪所引起的管线周围海床孔隙水压力和管线内应力的影响。     

An integrated three-dimensional model of wave-induced pore pressure and effective stresses in a porous seabed: II. Breaking waves

The evaluation of the wave-induced liquefaction potential is particularly important for coastal engineers involved in the design of marine structures. Most previous investigations of the wave-induced liquefaction have been limited to two-dimensional non-breaking waves. In this paper, the integrated three-dimensional poro-elastic model for the wave-seabed interaction proposed by [Zhang, H., Jeng, D.-S., 2005. An integrated three-dimensional model of wave-induced pore pressure and effective stresses in a porous seabed: I. A sloping seabed. Ocean Engineering 32(5/6), 701–729.] is further extended to simulate the seabed liquefaction potential with breaking wave loading. Based on the parametric study, we conclude: (1) the liquefaction depth due to breaking waves is smaller than that of due to non-breaking waves; (2) the degree of saturation significantly affects the wave-induced liquefaction depth, and no liquefaction occurs in full saturated seabed, and (3) soil permeability does not only significantly affect the pore pressure, but also the shear stresses distribution.     

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.     

随机波作用下埋管海床动态响应及液化研究

基于广义Biot动力理论和Longuet-Higgins线性叠加模型,构建波浪-海床-管线动态响应的有限元计算模型,求解随机波作用下,多层砂质海床中管线周围土体孔隙水压力和竖向有效应力的分布。采用基于超静孔隙水压力的液化判断准则,得出液化区的最大深度及横向范围,从而判断海床土体液化情况。考虑海洋波浪的随机性,将海床视为多孔介质,海床动态响应计算模型采用u-p模式,孔隙水压力和位移视为场变量。并考虑孔隙水的可压缩性、海床弹性变形、土体速度、土体加速度以及流体速度的影响,忽略孔隙流体惯性作用。参数研究表明:土体渗透系数、饱和度以及有效波高等参数对海床土体孔隙水压力、竖向有效应力和液化区域分布有显著影响。     

Wave-induced pore pressure around a composite breakwater

The interaction between wave, seabed and marine structure is a vital issue in coastal engineering, as well as marine geotechnical engineering. However, most previous investigations have been focused on the wave forces acting on the structure from the aspect of hydrodynamics. In this study, we will examine the problem of wave-seabed-caisson interaction from the aspect of marine geotechnical engineering. Based on Biot's poro-elastic theory (Biot, M.A., 1941. General theory of three-dimensional consolidation. Journal of Applied Physics 12, 155–164), a two-dimensional finite element model is proposed to investigate the wave-induced soil response in the vicinity of a caisson. Based on the numerical model, the water wave driven pore pressure around a caisson will be examined through a parametric analysis.     

Stability and liquefaction analysis of porous seabed subjected to cnoidal wave

Cnoidal wave theory is appropriate to periodic wave progressing in water whose depth is less than 1/10 wavelength. However, the cnoidal wave theory has not been widely applied in practical engineering because the formula for wave profile involves Jacobian elliptic function. In this paper, a cnoidal wave-seabed system is modeled and discussed in detail. The seabed is treated as porous medium and characterized by Biot's partly dynamic equations (u–p model). A simple and useful calculating technique for Jacobian elliptic function is presented. Upon specification of water depth, wave height and wave period, Taylor's expression and precise integration method are used to estimate Jacobian elliptic function and cnoidal wave pressure. Based on the numerical results, the effects of cnoidal wave and seabed characteristics, such as water depth, wave height, wave period, permeability, elastic modulus, and degree of saturation, on the cnoidal wave-induced excess pore pressure and liquefaction phenomenon are studied.     

Numerical study of wave-induced soil response in a sloping seabed in the vicinity of a breakwater

In this study, a mathematical integrated model is developed to investigate the wave-induced sloping seabed response in the vicinity of breakwater. In the present model, the wave model is based on the Volume-Averaged/Reynolds Averaged Navier–Stokes (VARANS) equations, while Biot's consolidation equation is used to govern the soil model. The influence of turbulence fluctuations on the mean flow with respect to the complicated interaction between wave, sloping seabed and breakwater are obtained by solving the Volume-Averaged k − ϵ model. Unlike previous investigations, the phase-resolved absolute shear stress is used as the source of accumulation of residual pore pressure, which can link the oscillatory and residual mechanisms simultaneously. Based on the proposed model, parametric studies regarding the effects of wave and soil characteristics as well as bed slopes on the wave-induced soil response in the vicinity of breakwater are investigated. Numerical results indicate that wave-induced seabed instability is more likely to occur in a steep slope in the case of soil with low relative density and low permeability under large wave loadings. It is also found that, the permeability of breakwater significantly affect the potential for liquefaction, especially in the region below the breakwater.     

The numerical study of wave-induced pore water pressure response in highly permeable seabed

The coupling numerical model of wave interaction with porous medium is used to study waveinduced pore water pressure in high permeability seabed.In the model,the wave field solver is based on the two dimensional Reynolds-averaged Navier-Stokes(RANS) equations with a k-ε closure,and Forchheimer equations are adopted for flow within the porous media.By introducing a Velocity-Pressure Correction equation for the wave flow and porous flow,a highly efficient coupling between the two flows is implemented.The numerical tests are conducted to study the effects of seabed thickness,porosity,particle size and intrinsic permeability coefficient on regular wave and solitary wave-induced pore water pressure response.The results indicate that,as compared with regular wave-induced,solitary wave-induced pore water pressure has larger values and stronger action on seabed with different parameters.The results also clearly show the flow characteristics of pore water flow within seabed and water wave flow on seabed.The maximum pore water flow velocities within seabed under solitary wave action are higher than those under regular wave action.     

An experimental study on non-linear progressive wave-induced dynamic stresses in seabed

In this study, the dynamic stresses within the seabed induced by non-linear progressive waves were explored through a series of hydraulic model tests on a movable bed within a wave flume. By comparing Stokes’ 2nd-order wave theory with the theory of wave-induced dynamic stresses within the seabed as proposed by Yamamoto et al. [1978. On the response of a poro-elastic bed to water waves. Journal of Fluid Mechanics 87 (1), 193–206.] and Hsu and Jeng [1994. Wane-induced soil response in an unsaturated anisotropic seabed of finite thickness. International Journal for Numerical and Analytical Methods in Geomechanics 18, 785–807], the experimental results show that the pressure on the seabed surface, the pore water pressure within the seabed as well as the vertical and the horizontal stresses are all smaller than their theoretical values. If we were to obtain the characteristics of seabed soil, the analytical solution of Hsu and Jeng [1994. Wane-induced soil response in an unsaturated anisotropic seabed of finite thickness. International Journal for Numerical and Analytical Methods in Geomechanics 18, 785–807] might agree to the simulation of the wave-induced effective stresses and shear stress in the sandy seabed. A different phase shift exists among all the three soil stresses. Their influences on the three dynamic stresses within seabed soil are important for seabed stability, and can be used in the verification of numerical models. In the whole, the non-linear progressive waves and the naturally deposited seabed are found to have a strong interaction, and the behavior of the induced dynamic stresses within the seabed is very complicated, and should be investigated integrally.     

随机波作用下海床动态响应分析

真实的海洋波浪是随机的,而前人对海床的动态响应分析大都是选用线性波或者Stokes波理论,对海床的模拟大都采用Biot拟静力模型,忽略了流体速度及土体位移加速度的影响。联合使用Longuet-Higgins随机波模型(采用Jonswap谱)以及动力u-p形式的海床响应计算模型,使用COMSOL Multiphysics多场耦合软件的PDE模块输入方程进行有限元计算,得到随机波作用下整体海床动态响应结果。将随机波结果与一阶Stokes波和椭圆余弦波结果进行对比,并对渗透系数和饱和度进行参数分析,研究表明渗透系数和饱和度对于随机波作用下海床动态响应影响显著。     

Coupled analysis for response and instability of sloping seabed under wave action

Wave-induced instability of seabed may cause damage to coastal and offshore structures. This issue has been investigated mostly for mildly sloping ($<5^{°}$) seabed considering uncoupled or one-way coupled response of wave and seabed interaction. However, some of the marine structures are founded on seabed with steeper slopes. In this study, the wave-induced response and instability of sloping seabed are evaluated using a coupled finite element model. The interaction between fluid and porous seabed accounting for the effect of fluid motion on the seabed response, and conversely the effect of seabed response on the fluid motion (but not on the surface wave profile) is considered. The results indicate that the system response (fluid pressure, stresses, etc.) and the extent of instantaneously liquefied zone within the sloping seabed with significant steepness are lesser than those for horizontal seabed. Moreover, for typical sediment and wave characteristics, for the flat seabed, the response obtained from fully coupled analysis is not significantly different from those obtained by uncoupled analysis. For the sloping bed, such difference is slightly greater as compared to that for the flat bed.     

波浪引起的海床失稳机理及有关孔隙水压力的讨论

近年来在海洋工程中由于海床失稳和下沉而引起的工程事故时有发生,因此波浪、海床和建筑物的相互作用问题越来越引起工程界的关注。文章首先介绍了波浪引起的海床失稳的基本机理,同时,由于孔隙水压力在海床失稳中起着重要作用,因此又通过引用前人成果,讨论了波土参数对孔隙水压力的影响。     

Effects of dynamic soil behavior and wave non-linearity on the wave-induced pore pressure and effective stresses in porous seabed

Most previous investigations for the wave-induced soil response have only considered the quasi-static soil behavior under linear wave loading. However, it is expected that the dynamic soil behavior and wave non-linearity will play an important role in the evaluation of wave-induced seabed response. In this paper, we include dynamic soil behavior and wave non-linearity into new analytical models. Based on the analytical solution derived, the effects of wave non-linearity on the wave-induced seabed response with dynamic soil behavior are examined. Numerical results demonstrate the significant effects of wave non-linearity and dynamic soil behavior on the wave-induced effective stresses. The applicable range of dynamic and quasi-static approximations is also clarified for engineering practice.     

Effects of cross-correlated multiple spatially random soil properties on wave-induced oscillatory seabed response

The evaluation of seabed response under wave loading is important for prediction of stability of foundations of offshore structures. In this study, a stochastic finite element model which integrates the Karhunen-Loève expansion random field simulation and finite element modeling of wave-induced seabed response is established. The wave-induced oscillatory response in a spatially random heterogeneous porous seabed considering cross-correlated multiple soil properties is investigated. The effects of multiple spatial random soil properties, correlation length and the trend function (the relation of the mean value versus depth) on oscillatory pore water pressure and momentary liquefaction are discussed. The stochastic analyses show that the uncertainty bounds of oscillatory pore water pressure are wider for the case with multiple spatially random soil properties compared with those with the single random soil property. The mean pore water pressure of the stochastic analysis is greater than the one obtained by the deterministic analysis. Therefore, the average momentary liquefaction zone in the stochastic analysis is shallower than the deterministic one. The median of momentary liquefaction depth generally decreases with the increase of vertical correlation length. When the slope of the trend function increases, the uncertainty of pore water pressure is greatly reduced at deeper depth of the seabed. Without considering the trend of soil properties, the wave-induced momentary liquefaction potential may be underestimated.     

Wave-induced seabed response in shallow water

To simulate the wave-induced response of coupled pore fluids and a solid skeleton in shallow water, a set of solutions with different formulations (fully dynamic, partly dynamic, and quasi-static) corresponding to each soil behavior assumption is presented. To deal with Jacobian elliptic functions involved in the cnoidal theory, a Fourier series approximation is adopted for expanding the boundary conditions on the seabed surface. The parametric study indicates the significant effect of nonlinearity for shallow water wave, which also enhances the effect of soil characteristics. The investigation of the applicability of reduced formulations reveals the necessity of a partly or even fully dynamic formulation for the wave-induced seabed response problem in shallow water, especially for thickened seabed. The analysis of liquefaction in the seabed indicates that the maximum depth of liquefaction is shallower, and the width of liquefaction is broader under cnoidal wave loading. The present analytical model can provide more reasonable result for the wave-induced seabed response in the range of shallow water wave.     

Three-dimensional poro-elastic integrated model for wave and current-induced oscillatory soil liquefaction around an offshore pipeline

To obtain a better understanding of the oscillatory soil liquefaction around an offshore pipeline, a three-dimensional integrated model for the wave–seabed–pipeline interaction (WSPI) is proposed by combining the Reynolds-Averaged Navier–Stokes equations for flow simulations and the dynamic Biot’s equation (“u-p” approximation) for the poro-elastic seabed model. Compared with previous investigations, the wave–current interaction is included in the present WSPI system. At a given time step, the wave pressure extracted from the flow model is applied on the seabed surface to determine the corresponding oscillatory seabed response around an offshore pipeline. The integrated numerical model is first validated using previous laboratory experiments. Then, a parametric study is conducted to examine the effects of flow obliquity and pipeline burial depth on the soil response around an offshore pipeline. Numerical results indicate that the soil under the pipeline is more susceptible to liquefaction at a reduced flow obliquity and pipeline burial depth. Moreover, the liquefaction depth in the case where the wave travels along the current can increase by 10%–30% compared to that in the case where the wave travels against the current, when the magnitude of the current velocity is 1 m/s.     

Spatial and temporal variation process of seabed dynamic response induced by the internal solitary wave

Internal solitary wave(ISW) is often accompanied by huge energy transport, which will change the pore water pressure in the seabed. Based on the two-dimensional Biot consolidation theory, the excess pore water pressure in seabed was simulated, and the spatiotemporal distribution characteristics of excess pore water pressure was studied. As the parameters of both ISW and seabed can affect the excess pore water pressure, the distribution of pore water pressure showed both dissipation and phase lag...     

Numerical study for waves propagating over a porous seabed around a submerged permeable breakwater: PORO-WSSI II model

The phenomenon of the wave, seabed and structure interactions has attracted great attentions from coastal geotechnical engineers in recent years. Most previous investigations have based on individual approaches, which focused on either flow region or seabed domain. In this study, an integrated model (PORO-WSSI II), based on the Volume-Averaged/Reynolds-Averaged Navier-Stokes (VARANS) equations and Biot's poro-elastic theory, is developed to investigate the mechanism of the wave-permeable structure-porous seabed interactions. The new model is verified with the previous experimental data. Based on the present model, parametric studies have been carried out to investigate the influences of wave, soil and structure parameters on the wave-induced pore pressure. Numerical results indicated: (i) longer wave period and larger wave height will obviously induce a higher magnitude of pore pressure at the leading edge of a breakwater; (ii) after a full wave-structure interaction, the magnitude of pore pressure below the lee side of a breakwater decreases with an increasing structure porosity while it varies dramatically with a change of structure height; and (iii) the seabed thickness, soil permeability and the degree of saturation can also significantly affect the dynamic soil behaviour.     

波浪加载下海底土质特性变化的研究

通过不同的制样方法 ,在水槽中模拟了多种海底在波浪作用下的变化行为。试验发现 ,加压排水固结的砂质粉土海底在波浪作用下较稳定 ;加压不产生排水而固结的砂质粉土海底易受波浪的冲刷 ;自然条件和轻微振动的砂质粉土海底最易受到波浪的扰动破坏 ,形成塌陷凹坑。粘粒含量较高的粉质粘土对波浪的反映不敏感 ;下卧软土层土体在上覆压力下的变形量与含水量关系密切 ,含水量越高 ,变形越大。试验结果证明 ,波浪会引起沉积物性质发生改变 ,也是引起海底形态变化的主要原因。     

The effects of variable permeability on the wave-induced seabed response

The topic of wave-seabed interaction is important for civil engineers with regard to stability analysis of foundations for offshore structures. Most previous investigations of such problems have simply assumed a seabed with uniform permeability, even if the evidence of variable permeability has been reported in the literature. This paper presents a finite-element model for investigating the wave-induced seabed response in a porous seabed, with variable permeability as a function of burial depth. The present finite formulation is established by using a combination of semi-analytical techniques and the Galerkin method. Based on the present numerical model, together with the Mohr-Coulomb failure criterion, the wave-induced seabed instability is estimated. The numerical results indicate that variable permeability affects the wave-induced seabed instability significantly, especially for gravelled seabeds.