上海交通大学学报(自然版), 2021, 55(6): 631-637 doi: 10.16183/j.cnki.jsjtu.2020.108

孤立波和海流作用下单桩基础局部冲刷及保护的数值分析

黎蔚杰, 张琪, 廖晨聪, 周香莲,, 刘晨晨

土木工程系, 上海 200240

Numerical Analysis of Local Scour and Protection of a Single Pile Around a Seabed Under Solitary Wave and Current

LI Weijie, ZHANG Qi, LIAO Chencong, ZHOU Xianglian,, LIU Chenchen

Department of Civil Engineering, Shanghai Jiao Tong University, Shanghai 200240, China

通讯作者: 周香莲,女,副教授,博士生导师.电话(Tel.): 021-34204833; E-mail:zhouxl@sjtu.edu.cn.

责任编辑: 陈晓燕

收稿日期: 2020-04-11  

基金资助: 国家自然科学基金(41572243)
国家自然科学基金(41727802)

Received: 2020-04-11  

作者简介 About authors

黎蔚杰(1994-),男,江西省抚州市人,硕士生,从事海洋土力学的研究

摘要

采用数值模拟的方法研究了孤立波和海流共同作用下单桩周围海床的局部冲刷.通过求解RNG k-ε湍流模型封闭的N-S方程,准确地模拟了桩周流场的变化情况.泥沙输运模型考虑泥沙的夹带、悬移质输送、沉积以及推移质输送的泥沙冲刷机制.对于斜坡床面,考虑重力对泥沙颗粒的影响,并对临界希尔兹数进行修正.将数值结果与实验数据进行比较,验证了数值模型的准确性.在此基础上,详细分析在不同波高条件下的桩基冲刷发展过程,初步研究了桩基保护层对冲刷的影响.结果表明:波-流共同作用下桩基周围海床的冲刷深度比单纯流或孤立波作用下冲刷深度更大;倾斜海床下的局部冲刷深度更大;保护层的设置可以有效地降低桩周局部冲刷深度,保护层粒径的大小对防护效果的影响较大.

关键词: 单桩; 冲刷; 孤立波; 海流; 冲刷保护

Abstract

A numerical simulation is conducted to investigate the local scour of the seabed around the pile under the actions of both solitary waves and currents. By solving the N-S equations enclosed by the RNG k-ε turbulence model, the change of the flow field around the pile is simulated accurately. The sediment transport model considers all the sediment transport processes including entrainment, suspended load transport, deposition, and bed-load transport. For the slope seabed, the influence of gravity on the sediment particles is considered, and the shields parameter is modified. The accuracy of the model is verified by comparing the numerical results with the experimental data. Besides, the process of scour development at different wave heights is investigated in detail. The effect of pile foundation protection layer on scour is preliminarily studied. The results show that the combined solitary wave-current conditions have a significant effect on the scour depth around the pile than the current-only or wave-only. The scour depth around the pile on the inclined seabed is larger than that of the horizontal seabed. The setting of protection layer can effectively reduce the local scour depth around the pile. In addition, the particle size of the protective layer has a great influence on the protective effect.

Keywords: monopile; scour; solitary wave; current; scour protection

PDF (6015KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

黎蔚杰, 张琪, 廖晨聪, 周香莲, 刘晨晨. 孤立波和海流作用下单桩基础局部冲刷及保护的数值分析[J]. 上海交通大学学报(自然版), 2021, 55(6): 631-637 doi:10.16183/j.cnki.jsjtu.2020.108

LI Weijie, ZHANG Qi, LIAO Chencong, ZHOU Xianglian, LIU Chenchen. Numerical Analysis of Local Scour and Protection of a Single Pile Around a Seabed Under Solitary Wave and Current[J]. Journal of shanghai Jiaotong University, 2021, 55(6): 631-637 doi:10.16183/j.cnki.jsjtu.2020.108

随着海上风电场的发展和建设,近海地区建造了许多的海上桩基础.桩基的存在改变了周围的局部流场,导致桩基周围发生局部冲刷.随着冲刷深度的增加,桩的埋置深度减小,进一步降低了桩基的安全性甚至引起上部结构物的塌陷.因此,有必要研究海洋环境荷载作用下桩基周围冲刷机理.

近年来,国内外学者对桩基周围的局部冲刷问题进行了大量的研究.Ma等[1]通过实验的方法研究了潮流作用下群桩基础的冲刷规律,考虑了在单向水流作用和潮流作用下典型的哑铃型桩群的局部冲刷特性.Xu等[2]通过水槽实验研究了孤立波作用下桩式防波堤周围的冲刷规律,考虑了孤立波波高和局部水深对最大冲刷深度及最大沉积高度等的影响.刘红军等[3]通过实验的方法研究了局部冲刷对海上风机承载性能的影响.漆文刚等[4]建立桩基-土的数值模型研究了冲刷对海上风机单桩桩基础水平承载特性的影响.Baykal等[5]建立了三维数值模型,研究了波浪作用下圆柱桩周的冲刷和回填过程,结果表明,波浪作用下海床的冲刷和回填主要是由桩后的尾涡控制.段玮玮等[6]通过水槽实验研究了波浪作用下桥梁基础冲刷防护装置对桩周局部冲刷的影响,结果表明,桩周最大冲刷深度与水深和护圈直径成反比关系.Petersen等[7]通过一系列实验研究桩周围冲刷保护层自身的边缘冲刷过程,结果表明,平衡冲刷深度与桩的直径以及冲刷保护层的宽度与厚度之比成正比.

目前,已有一些关于波和水流联合作用下桩基周围冲刷的研究,但很少有人研究孤立波和海流作用下桩基周围的冲刷机理.本文建立了包含单桩-海床-波和流的三维数值模型,研究在孤立波和海流共同作用下位于倾斜海床中的桩周围的冲刷过程.数值结果表明该模型可以很好地模拟桩基局部冲刷.在此基础上,通过将不同波高的孤立波和不同的流速进行组合,研究波流共同作用下桩周围局部冲刷发展的规律,探讨了最大平衡冲刷深度以及冲刷坑形态变化的发展规律.另外,通过设置防冲刷层,研究冲刷保护层在冲刷过程中的作用.

1 计算模型和数值方法

1.1 水动力模型

流体运动的控制方程为N-S方程,包含描述流体运动的连续性方程和动量方程.其表达式为

(uAx)x+(vAy)y+(wAz)z=0
ut+1φFuAxux+vAyuy+wAzuz=-1ρpwx+gx+fxvt+1φFuAxvx+vAyvy+wAzvz=-1ρpwy+gy+fywt+1φFuAxwx+vAywy+wAzwz=-1ρpwz+gz+fz

式中:u、v及w分别为流体x、y及z方向的速度;φF 为流动流体的体积分数;Ai 分别为x、y、z方向上流动流体的面积分数(i=x,y,z);ρ为流体密度;pw 为平均动水压力;gi 分别为x、y、z方向的重力加速度;fi 分别为x、y、z方向上流体的黏滞力加速度.

1.2 湍流模型

为了准确地模拟流场和桩基的相互作用,采用 RNG k-ε湍流模型对桩周局部冲刷进行了研究[8],湍流模型的表达式为

ket+1φFuAxkex+vAykey+wAzkez=Te+Be+Dke-εt
εtt+1φFuAxεtx+vAyεty+wAzεtz=D1εtke(Te+D3Be)+De-D2εt2ke

式中: t为时间;ke 为与湍流速度波动相关的湍流动能;Te 为由于速度梯度引起的动能产生项;Be 为浮力引起的湍流动能产生项;εt 为湍流动能耗散率;Dke 和De 为与Ai 、φF 以及扩散系数等参数相关的扩散项;D1,D2 和D3 为无量纲参数,其中D1 和D3 默认值分别是1.42和0.2;D2 是通过ke 和Te计算而来,其表达式为

D2=η(1-η/η0)1+β'η3

式中: η=2Eijkeεt,Eij=12uixj+ujxi;η0=4.38;β'=0.012.

1.3 泥沙输运模型

泥沙输运模型通过对沉积物侵蚀、移送和沉积的预测,模型整个冲刷过程.泥沙的启动和沉积方程如下[9]:

ulift,i=αinsd*0.3(θi-θd,i)1.5gdiρi-ρρ
usetting,i=vfdi(10.362+1.049d*3)2-10.36

式中:ulift,i为泥沙的启动速度;usetting,i为泥沙的沉降速度;αi为与泥沙颗粒启动相关的启动系数,其参考值为0.018;ns为垂直于泥沙颗粒的外法向量;d*为泥沙颗粒的无量纲颗粒直径; θi 为计算希尔兹参数;θd,i 为临界希尔兹参数,对于斜坡床面,修正后的表达式θd,i1

θd,i1=θd,icosψsinβctanφi+cos2βctan2φi-sin2ψsin2βctanφi

式中:ψ为流场和床面的夹角;βc 为床面倾斜角;φi为泥沙休止角;ρ为水的密度;di为泥沙颗粒的直径;ρi为泥沙颗粒的密度;vf为水的运动黏度.

泥沙的悬移质运动的控制方程为

Cs,it+Δ·(us,iCs,i)=Δ·Δ(DfCs,i)

式中: Cs,i 为悬移质的质量分数;us,i 为悬移质的运输速度;Df为扩散系数.

泥沙的推移质输运描述泥沙颗粒在流体作用下沿海床表面的滚动、跳跃和滑动的运动过程,Flow-3D采用Meyer-Peter Muller公式[10]计算推移质运输速度ubed,i,控制方程如下:

ubed,i=qb,idiφb,ifb

式中:qb,i为单位宽度的推移质运输速率

qb,i=8gρi-ρfρfdi312

φb,i 为组成海床的不同泥沙种类所占的体积分数;fb 为泥沙临界填充分数;ρf 为对应泥沙的密度.

通过输沙平衡来计算床面高程h的变化,表达式如下:

ht=-11-p'qb,ixi

式中:p'为泥沙孔隙率.

1.4 波浪理论

波浪模型中,孤立波的波面高度Hw方程和波速c方程基于Munk的理论[11],表达式为

Hwd=NMsinM1+HwdcosM1+Hwd+coshMXd
c=u0+c0

式中:d为静水深度;N和M为波高和波速的修正值, M=3c(c为修正系数),N=2H/d(H为波高);u0 为未受扰动的水流速度; c0=g(d+H)为波浪在静水中的速度;X=x-ct.

2 边界条件和网格

数值模型的边界条件如图1所示.模型入流边界选择波浪边界,海床底部设置为壁面边界,模型的出口设置为出流边界,模型的两侧边界设置为对称边界,模型的顶部采用压力边界来模拟大气压环境.为了防止波浪的反射,在出口端设置消波层.

图1

图1   数值模型的边界条件

Fig.1   Boundary conditions of numerical model


采用流体体积(VOF)法捕捉流体表面的位置,采用浸没边界(FAVOR)法方法对三维数值模型进行网格划分,每一个单元都包含了所有变量的信息.在每一个时间步长范围内都计算了描述沉积物的面积分数Ai和体积分数φF,以此来反映沉积物冲刷模型几何形状的改变.与传统的有限差分法相比,避免了对结构物边界拟合程度不高的缺点,同时网格的数量大大减少,节约了模型的计算时间.

3 模型验证

在进行数值研究之前,对单桩周围局部冲刷进行验证.参考Roulund等[12]的物理模型实验建立了数值模型,光滑圆柱的直径为0.1m,水深0.4m,平均入流速度为0.46m/s.设置沙床厚度为0.2m,沙床中砂石平均粒径为0.26mm,密度为 2650kg/m3,休止角为32°.

图2是圆柱迎水侧冲刷深度计算值随时间变化的关系和实验结果的对比,图中S为冲刷深度.可以看出,冲刷深度的数值模拟结果和实验数据拟合较好.因此,本文建立的圆柱冲刷三维数值模型能较好反应桩周的实际冲刷情况.

图2

图2   圆柱迎水侧冲刷深度数模拟值与实验数据的对比图

Fig.2   Comparison of simulated values and experimental data of scour depth at the upstream side of the pile


4 数值结果和讨论

图3所示,数值模型的尺寸为长250m,宽16m,z方向高度为20m.笛卡尔坐标的原点位于倾斜海床的起点(x=0).海床的倾斜度m为1∶50,由直径为0.5mm和0.85mm两种不同直径的泥沙颗粒组成.泥沙颗粒的密度为ρs=2650kg/m3,临界希尔兹参数为θd,i=0.048,夹带系数设置为0.018.推移质系数根据默认值设置为8.0,泥沙的休止角为30°,桩的直径为2m.桩周附近网格尺寸在xy方向为 0.2m,其他区域是0.4m,海床的网格尺寸为0.15m.

图3

图3   数值模型整体示意图

Fig.3   Integrated sketch of numerical model


4.1 孤立波和流共同作用对冲刷的影响

图4所示为单纯水流作用、单纯孤立波作用及孤立波和水流共同作用下桩周最大冲刷深度随时间变化的曲线图.可以看出,孤立波和水流共同作用下的冲刷深度发展速度快于单独水流作用或单独波浪作用,最大冲刷深度大于波浪和水流单独作用下的最大冲刷深度的线性叠加.在孤立波还未抵达桩前的时候,冲刷深度的发展趋势和单纯流作用下一致,在孤立波抵达桩前,冲刷深度曲线在短时间迅速下降.

图4

图4   单纯孤立波、水流、波流共同作用下冲刷深度变化图

Fig.4   Changes of scour depth under the action of wave-only, current-only, and wave-current


4.2 海床倾斜、孤立波波高对冲刷深度的影响

图5为相同条件下(波高H=4m,流速v=2m/s)位于水平和倾斜海床上桩周冲刷深度随时间的变化过程,从图中可以看出,倾斜海床桩周最大冲刷深度大于水平海床桩周最大冲刷深度.

图5

图5   水平海床和倾斜海床下桩周最大冲刷深度变化图

Fig.5   Changes of maximum scour depth around pile located on flat seabed and inclined seabed


图6为不同孤立波波高下桩周海床冲刷高程图.可以看出,不同波高下,虽然局部冲刷的最大冲刷深度都出现在桩前侧45° 的方向,但是桩周冲刷坑的形态有很大的不同.同时,孤立波波高越大,冲刷深度越大,冲刷范围也更广泛.

图6

图6   桩周海床冲刷高程图

Fig.6   Bed elevation contours around pile


4.3 保护层对冲刷深度的影响

图7所示,添加保护层的数值模型的尺寸为 35m×16m×7m,桩的直径1.5m.在桩周设置冲刷保护层,保护层形状为直径6m,高度为1m的圆柱,为了精确地捕捉到保护层形状的变化,冲刷保护层区域附近网格尺寸xy方向为0.08m.海床的厚度为3m,最小尺寸为0.1m.流边界处水流速度为2m/s.其余冲刷参数设置保持不变.

图7

图7   数值模型整体示意图(m)

Fig.7   Integrated sketch of numerical model (m)


图8所示为添加桩基保护层后单桩桩周的三维冲刷坑形态图.从图8(a)可以看出,保护层的存在使得海床在一开始并没有发生冲刷,随着时间的发展,保护层先被侵蚀,此时海床开始发生冲刷,但是冲刷坑的发展非常缓慢(图8(c)、8(d)、8(e)),最终残余的保护层围绕在桩的周围,且冲刷坑的形态基本不变,而桩基后方的沙堆进一步向下游方向移动,形成较长的沙纹(图8(f)).

图8

图8   单桩桩周的三维冲刷坑形态随时间的冲刷过程

Fig.8   Scour process of three-dimensional scour pit shape around a single pile over time


图9所示为t=500s时三维桩周冲刷坑形态和桩周海床冲刷高程图.比较图9(b)、9(d)与图9(a)、9(c)可以看出,设置保护层后桩基周围的冲刷坑形态发生明显的改变,且海床的最大冲刷深度明显减少,为0.43m(图9(d),海床的实际冲刷深度要减去保护层的厚度1m), 这表明保护层对单桩局部冲刷的影响较大.

图9

图9   t=500s时三维桩周冲刷形态和桩周海床冲刷高程图

Fig.9   Three-dimensional scour morphology and elevation of seabed scour around pile at t=500s


为了进一步研究保护层对冲刷的影响,对3种不同粒径的保护层进行了计算.图10所示为3种情况下的桩周最大冲刷深度随时间变化的曲线图,图中百分数为碎石质量分数.可以看出,随着组成保护层碎石粒径的增大,桩周的最大冲刷深度也明显减少.这是因为保护层可以在一定程度上改变桩周附近流场,降低近海边界层的水流速度,使得泥沙颗粒难以被起动和冲刷,这种保护效果随着保护层粒径增大更加明显.

图10

图10   不同粒径的保护层下桩周局部冲刷的计算结果

Fig.10   Calculation results of local scour depth under scour protection at different grain sizes


5 结论

本文研究在孤立波和水流共同作用下,位于倾斜海床上的桩基周围的局部冲刷和冲刷保护措施.分析了孤立波波高和海床的倾斜对桩周局部冲刷最大冲刷深度和冲刷坑形态的影响.在此基础上,研究了设置冲刷防护措施对冲刷深度的影响.主要结论如下:

(1) 孤立波和水流共同作用下桩周的冲刷深度大于波浪和水流单独作用下的线性叠加.

(2) 在相同的波浪流条件下,倾斜海床桩周最大冲刷深度大于水平海床桩周最大冲刷深度.

(3) 为了减少冲刷深度,提高海床的抗冲刷能力,可以在桩周设置冲刷防护(砾石、石块等),保护层的粒径大小对冲刷深度影响较大.

参考文献

MA L L, WANG L Z, GUO Z, et al.

Time development of scour around pile groups in tidal currents

[J]. Ocean Engineering, 2018, 163(10):400-418.

DOI:10.1016/j.oceaneng.2018.06.035      URL     [本文引用: 1]

XU C H, HUANG Z H, YAO Y.

A wave-flume study of scour at a pile breakwater: Solitary waves

[J]. Applied Ocean Research, 2019, 82(3):89-108.

DOI:10.1016/j.apor.2018.10.026      URL     [本文引用: 1]

刘红军, 杨奇.

局部冲刷对风机支撑系统承载性能的影响

[J]. 岩土力学, 2018, 39(2):722-727.

[本文引用: 1]

LIU Hongjun, YANG Qi.

Influence of local scour on bearing performance of wind turbine supporting system

[J]. Rock and Soil Mechanics, 2018, 39(2):722-727.

[本文引用: 1]

漆文刚, 高福平.

冲刷对海上风力机单桩基础水平承载特性的影响

[J]. 中国科学: 物理学力学天文学, 2016, 46(12):83-92.

[本文引用: 1]

QI Wengang, GAO Fuping.

Effects of scour on horizontal bearing behavior of monopile foundations for offshore wind turbines

[J]. Scientia Sinica (Physica, Mechanica & Astronomica) , 2016, 46(12):83-92.

[本文引用: 1]

BAYKAL C, SUMER B M, FUHRMAN D R, et al.

Numerical simulation of scour and backfilling processes around a circular pile in waves

[J]. Coastal Engineering, 2017, 122:87-107.

DOI:10.1016/j.coastaleng.2017.01.004      URL     [本文引用: 1]

段玮玮, 孙佳, 朱海芬, .

波浪作用下跨海桥梁桩基冲刷防护装置特性研究

[J]. 浙江海洋大学学报(自然科学版), 2019, 38(5):447-452.

[本文引用: 1]

DUAN Weiwei, SUN Jia, ZHU Haifen, et al.

Pro-perties of novel local scour protective equipment around pile of bridge under wave action

[J]. Journal of Zhejiang Ocean University (Natural Science) , 2019, 38(5):447-452.

[本文引用: 1]

PETERSEN T U, MUTLU SUMER B, FREDSØE J, et al.

Edge scour at scour protections around piles in the marine environment: Laboratory and field investigation

[J]. Coastal Engineering, 2015, 106(9):42-72.

DOI:10.1016/j.coastaleng.2015.08.007      URL     [本文引用: 1]

YAKHOT V, SMITH L M.

The renormalization group, the ε-expansion and derivation of turbulence models

[J]. Journal of Scientific Computing, 1992, 7(1):35-61.

DOI:10.1007/BF01060210      URL     [本文引用: 1]

MASTBERGEN D R, VAN DEN BERG J H.

Breaching in fine sands and the generation of sustained turbidity currents in submarine canyons

[J]. Sedimentology, 2003, 50(4):625-637.

DOI:10.1046/j.1365-3091.2003.00554.x      URL     [本文引用: 1]

MEYER-PETER E.

Formula for bed-load transport

[C]// The 2nd meeting for the international association for hydraulic structures research. Delft, Netherlands, 1948: 39-64.

[本文引用: 1]

MUNK W H.

The solitary wave theory and its application to surf problems

[J]. Annals of the New York Academy of Sciences, 1949, 51(3):376-424.

DOI:10.1111/nyas.1949.51.issue-3      URL     [本文引用: 1]

ROULUND A, SUMER B M, FREDSØE J, et al.

Numerical and experimental investigation of flow and scour around a circular pile

[J]. Journal of Fluid Mechanics, 2005, 534(12):351-401.

DOI:10.1017/S0022112005004507      URL     [本文引用: 1]

/