An analytical solution for wave-induced seabed response in a multi-layered poro-elastic seabed

In this study, a new analytical solution for the wave-induced seabed response in a multi-layered poro-elastic seabed is developed. The seabed is treated as a multi-layered porous medium and characterized by Biot’s theory. The displacements of the solid skeleton and the pore pressure are expressed in terms of two scalar potentials and one vector. Then, the Biot’s dynamic equation can be solved using Fourier transformation and reducing to Helmholtz equations. To obtain the general solutions for the multi-layered poro-elastic seabed in the frequency-wave-number domain, the transmission and reflection matrices (TRM) method is used to form the equivalent stiffness. Using the boundary conditions and continuous conditions, the frequency-wave-number domain solutions are obtained. Finally, the time-space domain solutions for the multi-layered poro-elastic seabed are obtained by means of the inverse Fourier transformation with respect to the horizontal coordinate. Based on the new solution, a parametric study is carried out to examine the effects of soil characteristics (number of layers, permeability and shear modulus) and wave characteristics (water depth and wave steepness) on seabed responses. The results indicate that the seabed response is affected significantly by permeability, shear modulus and relative water depth.     

Analysis of seabed instability using element free Galerkin method

Wave-induced seabed instability, either momentary liquefaction or shear failure, is an important topic in ocean and coastal engineering. Many factors, such as seabed properties and wave parameters, affect the seabed instability. A non-dimensional parameter is proposed in this paper to evaluate the occurrence of momentary liquefaction. This parameter includes the properties of the soil and the wave. The determination of the wave-induced liquefaction depth is also suggested based on this non-dimensional parameter. As an example, a two-dimensional seabed with finite thickness is numerically treated with the EFGM meshless method developed early for wave-induced seabed responses. Parametric study is carried out to investigate the effect of wavelength, compressibility of pore fluid, permeability and stiffness of porous media, and variable stiffness with depth on the seabed response with three criteria for liquefaction. It is found that this non-dimensional parameter is a good index for identifying the momentary liquefaction qualitatively, and the criterion of liquefaction with seepage force can be used to predict the deepest liquefaction depth.     

Wave-induced response of seabed: Various formulations and their applicability

In this study, a set of generalized analytical solutions are developed for the wave-induced response of a saturated porous seabed under plane strain condition. When considering the water waves originating in deep water and travelling towards the shore, their velocities, lengths and heights vary. Depending on the characteristics of the wave and the properties of the seabed, different formulations (fully dynamic, partly dynamic, quasi-static) for the wave-induced response of the seabed are possible. The solutions for the response with these formulations are established in terms of non-dimensional parameters. The results are presented in terms of pore pressure, shear stress and vertical effective stress distributions within the seabed. For typical values of wave period and seabed permeability, the regions of applicability of the three formulations are identified and plotted in parametric spaces. With given wave and seabed characteristics, these regions provide quick identification of the appropriate formulation for an adequate evaluation of the wave-induced seabed response.     

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.     

Evaluation of storm wave-induced silty seabed instability and geo-hazards: A case study in the Yellow River delta

Models based on the theoretical framework of soil mechanics are presented to evaluate storm wave-induced silty seabed instability and geo-hazards through a case study in the Yellow River delta. First, the transient and residual mechanisms of wave-induced pore pressure are analyzed. Three typical models (i.e., elastic model, pore pressure development mode and elasto-plastic model) are proposed to calculate wave-induced stresses in the seabed. Next, mechanisms and calculation methods of wave-induced seabed instability modes such as scour, liquefaction, seepage instability and shear slide are proposed. Typical results of storm wave-induced excess pore pressure and seabed instability are given and relevant discussions are made. At last, the formation mechanism of geo-hazards in the Yellow River delta is analyzed based on the proposed mechanism and calculated results. Results and analysis indicate that both transient and residual mechanisms are important to storm wave-induced response of silty seabed and hence the elasto-plastic model is more appropriate. Complete liquefaction does not happen, while other types of instability occur mostly within 2–6 m under the seabed surface. Wave-induced scour, seepage instability and shear slide are all possible instability modes under the 1-year storm waves, and scour is predominant for the 50-year storm waves. The formation mechanism of geo-hazards such as shallow slide and storm wave reactivation, pockmarks, silt flow and gully, disturbed stratum and hard crust in the Yellow River are well explained based on the proposed mechanisms and calculated results of storm wave-induced silty seabed instability.     

波浪作用下有孔粘弹性海底的响应

本文提出了一种有孔粘弹性底面模型。根据Biot的三维固线理论,讨论了波浪作用下有孔粘弹性底面中孔隙水压力响应和固体颗粒承受的有效应力的响应。若忽略粘性,本文的模型则退化到yamamoto的有孔弹性底面模型。从物理意义来说,本文提出的模型更全面合理。     

A 3-D model for ocean waves over a Columb-damping poroelastic seabed

The evaluation of the wave-induced seabed instability in the vicinity of a breakwater is particularly important for coastal and geotechnical engineers involved in the design of coastal structures. In this paper, an analytical solution for three-dimensional short-crested wave-induced seabed instability in a Coulomb-damping porous seabed is derived. The partial wave reflection and self-weight of breakwater are also considered in the new solution. Based on the analytical solution, we examine (1) the wave-induced soil response at different location; (2) the maximum liquefaction and shear failure depth in coarse and fine sand; (3) the effects of reflection coefficients; and (4) the added stresses due to the self-weight of the breakwater.     

Transient soil response in a porous seabed with variable permeability

The soil permeability of many natural marine sediments decreases with depth because of consolidation under overburden pressure. This is accompanied by a decrease in porosity and void ratio that also affect the permeability. Conventional theories for wave-induced soil response have assumed a homogeneous porous seabed. This paper presents a new approach for the wave-induced response in a soil matrix, with variable permeability as a function of burial depth. The soil matrix considered is unsaturated and anisotropic, and is subject to a three-dimensional wave system. The pore pressure and effective stresses induced by such a system are obtained from a set of equations incorporating a variable permeability. Verification is available through reduction to the simple case of uniform permeability. The results indicate that the effect of variable soil permeability on pore pressure and vertical effective stress may be significant, especially in a gravelled seabed and for unsaturated sandy soils.     

Wave-induced seabed instability in front of a breakwater

The wave-induced soil response in a porous seabed has become an important factor for the stability of offshore facilities, because many marine structures may have failed due to seabed instability and concomitant subsidence. An analytical solution is presented for the wave-induced soil response under the action of a three-dimensional wave system. Based on this general solution, the mechanism of seabed instability is then investigated. The general solutions for pore pressure and effective stresses are readily reducible to two dimensions for progressive waves, and are compared to theoretical and experimental work available. Some dominant factors affecting the wave-induced seabed instability are discussed; including permeability, seabed thickness and degree of saturation.     

Numerical Simulation of Seabed Response and Liquefaction due to Nonlinear Waves

Based on Biot＇s consolidation theory, a two-dimensional model for computation of the seabed response to waves is presented with the finite element method. Numerical results for different wave conditions are obtained, and the effects of wave non-lineafity on the wave-induced seabed response are examined. Moreover, the wave-induced momentary liquefaction in uniform and inhomogeneous seabeds is investigated. It is shown that the wave non-linearity affects the distribution of the wave-induced pore pressure and effective stresses, while the influence of wave non-linearity on the seabed liquefaction potential is not so significant.     

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.     

Pore pressure response of seabed in standing waves and its mechanism

Experiments were carried out in a newly-developed rectangular flume to investigate the pore pressure response of silty seabed in standing waves, and to emphatically discuss the physical mechanism of wave-induced pore pressure. We first analyzed the features of wave-induced stress, laid out the experimental setup, and particularly monitored the temporal and spatial variation of pore pressure. Then, we summarized the physical essence of wave-induced pore pressure, discussed the applicability of the existing models, and proposed several key points about developing new models. Results indicate that regular standing waves were obtained in the rectangular flume. The oscillating pore pressure (OPP) and the residual pore pressure (RPP) were observed simultaneously in the seabed under wave nodes or antinodes, and were found with coupling effect. The OPP and RPP were related to the elastic and plastic volumetric strains, respectively. Both elastic and plastic volumetric strains may be caused by the wave-induced cyclic spherical stress or deviatoric stress. The elastic model could only simulate the OPP and the pore pressure development mode could only simulate the RPP. The yielding under isotropic compression and the shear expansion of geotechnical materials, as well as the wave-induced stress paths, should be considered in the elastoplastic model to simulate the wave-induced pore pressure accurately.     

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.     

Wave-induced seabed instability around a buried pipeline in a poro-elastic seabed

The subject of the wave–seabed–structure interaction is important for civil engineers regarding stability analysis of foundations for offshore installations. Most previous investigations have been concerned with such a problem in the vicinity of a simple structure such as a vertical wall. For more complicated structures such as a pipeline, the phenomenon of the wave–seabed–structure has not been fully understood. This paper proposes a finite-difference model in a curvilinear coordinate system to investigate the wave-induced seabed response in a porous seabed around a pipeline. Based on the present numerical model, mechanism of the wave-induced soil response is examined. Employing Mohr–Coulomb failure criterion, the wave-induced seabed instability is also estimated. The numerical results indicate the importance of the effect of pipeline on the seabed response.     

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.     

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.     

Principal stress rotation and cyclic shear strength for an elliptical rotation stress path

Based on the analytical solutions for wave-induced soil response in an unsaturated anisotropic seabed of finite thickness, an elliptical, i.e., noncircular rotation stress path is proven to be a more common state in a soil element with the cyclic shear stresses due to the traveling linear wave. The influence of an elliptical stress path on the characterization of failure strength is analyzed using three new parameters representing the shape, size, and the orientation of the ellipse. A series of cyclic rotational shear tests on the reconstituted specimens of Chinese Fujian Standard Sand have been performed to investigate the effect of elliptical stress path. A strength function in term of failure cycles is derived to quantify the failure strength of a given sand within a seabed subjected to regular wave loading. The results provide a basis for the evaluation of liquefaction potential of seabed but also point to a unique backbone cycle shear strength curve for soil under principal stress rotation.     

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.     

Numerical study on the interaction among a nonlinear wave,composite breakwater and sandy seabed

Most previous investigations related to composite breakwaters have focused on the wave forces acting on the structure itself from a hydrodynamic aspect. The foundational aspects of a composite breakwater under wave-induced cyclic loading are also important in studying the stability of a composite breakwater. In this study, numerical simulations were performed to investigate the wave-induced pore water pressure and flow changes inside the rubble mound of the composite breakwater and seabed foundation. The validity and applicability of the numerical model were demonstrated by comparing numerical results with existing experimental data. Moreover, the present model clearly has shown that the instantaneous directions of pore water flow motion inside the seabed induced by surface waves are in good agreement with the general wave-induced pore water flow inside the seabed. The model is further used to discuss the stability of a composite breakwater, i.e., the interaction among nonlinear waves, composite breakwater and seabed. Numerical results suggest that the stability of a composite breakwater is affected by not only downward shear flow generating on the seaward slope face of the rubble mound but, also, a high and dense pore water pressure gradient inside the rubble mound and seabed foundation.     

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.