周德棕，毕明君,✉，章海东，张伟，张桂志

（1.南方海上风电联合开发有限公司, 广东 珠海 519080；
2.中国能源建设集团广东省电力设计研究院有限公司, 广东 广州 510663）

随着环境问题日益突出，优先发展新能源已引起各国的重点关注。综合分析国家能源转型要求、清洁能源消纳目标以及新能源成本快速下降等因素，预计“十四五”期间，我国年度新增风电装机有望达到25 GW，到2025年，全国新能源总装机规模在750 GW～800 GW，占全国电源总装机的26%～28%，发电量占比约为12%。根据江苏、广东、浙江、福建、上海等国家或地方政府已批复的海上风电发展规划进行测算，预计到2025年中国海上风电累计装机容量将达到30 GW左右，80%的装机集中在江苏、广东、福建等省份，江苏、广东有望建成千万千瓦级海上风电基地[1]。

单桩基础作为海上风电机组的支撑结构之一，具有结构形式简单、施工工艺成熟、建造成本较低等优势，广泛应用于国内外的海上风电工程。原本处于动态平衡的海床，由于单桩基础的安装，将会改变局部的水动力条件，使得维持动态平衡的外部条件被打破，其直接后果是导致结构物附近的泥沙发生局部冲刷。

桩基周围的局部冲刷，可能会致使基础结构频率改变而影响风电机组正常发电，甚至导致整体基础结构丧失稳定性而引发灾难。国内外已经对桩基周围的局部冲刷问题进行了广泛研究，研究方向主要集中在波浪作用或波流联合作用对桩基础周围海床的冲刷侵蚀[2-3]、冲刷深度和冲刷范围[2,4]开展物理模型试验或者数值计算等研究。高徐昌等[5]基于能量守恒原理，提出了1种潮流作用下的局部冲刷深度预测公式；
Liang等[6]根据796个现场实测数据，分析了中美规范的计算方法，发现两者的冲刷深度预测均偏于不安全；
韩海骞[7]通过总结杭州湾桥梁试验结果，建立了潮流作用下的大直径墩柱冲刷计算公式。Wang等[8]基于水槽试验，总结了适用于波浪引起冲刷计算的公式。张玮等[9]对比了几种不同的基础冲刷计算公式，改进了波浪作用的方式并且选择了一种推荐公式。杜硕等[10]结合江苏如东海上风电场的冲刷监测数据，应用4种统计学评价指标，提出适用于波流作用下的单桩基础局部冲刷深度预测公式，预测值与现场监测数据吻合较好。

局部冲刷会使得桩基基础的承载力下降，横向受力不均匀，最终会导致桩基基础发生在位失稳，严重威胁风机的安全运营。Mostafa[11]分别基于Plaxis软件和Lpile软件，分析了黏性土和砂土地基中单桩基础水平承载力的变化特点，并对比分析了两种结果；
马殿滨等[12]考虑了冲刷前后土体物理特性的变化特点；
杨少磊和马宏旺[13]建立了考虑冲刷情况下海上风电单桩基础的优化设计模型，结果表明，考虑冲刷的工况下适当增加单桩基础的径厚比值比增加单桩长度经济性更好。

减少桩基局部冲刷对桩基影响另一方面可以从加强对床面防护入手，和庆冬和戚建功[14]通过某海上风电项目运行1年后大直径单管桩基础出现严重冲刷现象引出问题，对桩基础冲刷原理、目前国内外常规采取的防冲刷方案进行了概述，对某海上风电项目采用淤泥固化方案进行桩基础防冲刷试验的原理、参数、方案、工程实施进行了详述。本文从风机基础砂被和固化土防护方案开展了物理模型试验，分析了不同防护方案的防护效果及存在不足，对未来桩基础防护方案给出建议。

本工程项目场址位于珠江口外的伶仃洋海域，场址属于广东省海上风电工程规划中场址，本工程规划容量为 198 MW，一期总容量为 120 MW，二期工程位于一期工程南侧，规模78 MW。一期建设34台3 MW及 3台6 MW风机，完成了 34台3 MW风机基础的土建及31台3 MW风机的安装调试工作。相应配套的集电海缆、三角岛升压站、集控中心、两回110 kV送出海缆中的首回均已建设完成。

本工程采用大直径单桩基础，考虑到工程所在区域的水文、地质条件及单桩基础周围土易冲刷的特点，本次针对直径8.0 m桩基础开展了波流局部冲刷试验研究，并对砂被等防护效果进行验证。主要研究内容和水动力条件见表1。

表1 试验内容及和水动力条件Tab.1 Test contents and hydrodynamic conditions

2.1 模型设计

本次物理试验在水槽中进行，水槽规格为：长×宽×深=44.0 m×3.0 m×1.2 m，水槽一端配有造波机系统，可根据拟定的波浪要素，自动产生所需的不规则波浪。造波机后侧和水槽的另一端设置消浪斜坡。水槽内安置了4台大流量潜水泵。在距离造波板 20.0 m 处布置 3.0 m×3.0 m 动床，具体布置见图1。

图1 模型试验布置图Fig.1 Layout of model test

波高测量采用电容式波高仪，采用DS30多功能自动采集系统采集，由计算机控制并进行数据处理。流速采用旋桨流速仪进行测量。地形测量为在水槽量测安装滑道及滑块，采用激光测距仪进行测量，见图2。

图2 试验动床及风电桩基础模型照片Fig.2 Picture of erodible beds and pile foundation model of wind turbine

2.2 模型相似及模型沙选取

1）模型相似

竖向圆柱形状的结构是海洋工程中的重要构筑物，其周围的局部冲刷直接关系到建筑物的安全与稳定。在圆柱周围冲刷坑形成过程中，其周围的水流、波浪运动使泥沙起动并输移，因此，在确定物理试验模型的相似比尺时，不但要考虑结构尺寸按比例缩尺，更应该考虑水流、波浪以及泥沙运动的相似比尺和冲淤形态的整体相似。本次试验将根据工程结构特点和工程海域的水动力条件，采用大比尺模型试验。模型设计比尺相似要求如下：

由水流平面二维运动方程：

主要比尺关系：

重力相似（流速比尺）：

阻力相似（糙率比尺）：

水流运动相似：

式中：

λh、λl−模型几何比尺；

λC、λn−速度比尺；

λt−时间比尺。

2）模型沙选取

由于本模型的主要任务是研究风机基础附近冲刷坑深度及范围，因此模型沙选择主要考虑泥沙起动相似λvc=λv。

由风电基础现场土层采样资料，海床表层1～7 m为淤泥，塑性指数（Ip）在 11.8～23.8；
其下 2～8 m 皆为淤泥质粉质黏土，塑性指数（Ip）在8.2~29.2。

工程区土层泥沙起动摩阻流速（u*c）可按下式计算：

式中：

Ip−土的塑性指数（考虑到本次试次风电基础的可能冲刷深度，取为8.2进行计算）；

τc−临界起动切应力（N/m2）；

ρ−水的密度。

结合以往的研究经验，根据比尺相似的原则，经过比较选择，本次物理模型试验的模型砂采用中值粒径d50为 0.30 mm，密度γs为 1.33 t/m3的煤粉，原来模拟工程海域淤泥土层的海床。

本次物理模型试验选用窦国仁公式计算纯水流作用下泥沙的起动流速：

式中：

d−泥沙中值粒径；

Δ−糙率，当d<0.50 mm 时，取Δ=0.50 mm；

εk−泥沙粘结力参数，原型沙取εk=2.56 cm3/s2，煤粉取εk=0；

δ−薄膜水厚度，取 0.21×10−4cm；

g −重力加速度，取 9.81 m/s2。

对于波浪作用下泥沙起动计算，原型沙起动条件仍然采用起动摩阻流速(u*c)。

Madsen和Grant[15]将泥沙起动标准曲线引入到波浪作用下。提出采用Jonsson[16]推荐的波浪摩擦系数fw来计算振荡流作用下动床床面的切应力，其中床面剪应力的瞬时最大值τcm为：

式中：

umc−微幅波理论确定的波浪水质点近底水平速度的最大值；

式中：

ν−粘滞系数；

am−波浪水质点近底水平振幅最大值，由微幅波理论确定。

本次物理模型试验采用Zhou等[17]提出的公式，进行波浪作用下的模型沙起动计算，即：

当d＜0.05 cm 时，为层流区：

当d≥0.05 cm 时，为紊流区：

式中：

近年来农村商业银行的不良贷款余额、不良贷款率均出现大幅攀升。2015年商业银行不良贷款率上升至1.59%，导致这种情况发生有两个层面的原因，第一个层面是农民受教育程度不高，文化程度和素质有限，这在一定程度上对农民的信用意识造成了影响，直接导致在农村地区躲债赖债成为一种普遍现象，这大大打击了金融机构向农户提供金融服务的积极性；
第二个层面就是邯郸市广大农村地区征信体系尚未健全。农村地区大部分中小企业和农户没有信用档案，缺乏对失信者的惩罚机制，金融机构也因整体信用环境较差而惜贷，因此严重影响了农村金融体系的发展和稳定。

s− 泥沙的比重；

d− 泥沙中值粒径；

v− 粘滞系数。

计算得到的水流、波浪作用下泥沙起动相似条件分别见表2和表3。

表2 水流作用下泥沙起动流速及其比尺Tab.2 Incipient velocity of sediment and model scale under the action of current

表3 波浪作用下泥沙起动流速及其比尺Tab.3 Incipient velocity of sediment and model scale under the action of wave

由表可见，水流、波浪作用下原型沙和模型沙起动流速相似比尺基本接近流速比尺值，本次试验模型比尺取为1:30，模型试验动力条件见表4。

表4 局部冲刷模型试验内容和动力条件Tab.4 Contents and hydrodynamic conditions of model test for local scour

2.3 试验过程

本次试验开始前，先在动床段周围布置波高测量仪，率定水槽内试验时的设计波浪要素，将率定系数记录到造波系统的计算机里，供冲刷试验时使用。

进行冲刷试验之前，先在动床段的中心位置安装缩尺后的桩基模型，然后将模型砂铺满动床段，最后开始放水至试验设计水深。

进行模拟波流共同作用冲刷试验过程中，首先，调试水槽内水流的流速，当实测的平均流速达到试验设计流速时；
然后，启动造波机，按前述设定的率定系数，施加波浪条件；
最后，连续监测单桩基础周围局部地形变化情况，当连续2次监测到地形数据（包括冲刷深度和冲刷范围）基本一致时，则认为本次试验的局部冲刷达到了冲淤平衡状态，停止试验，测量模型周围冲刷坑的范围和深度。

对于防护试验，铺平模型沙后铺设砂被或采用水泥砂浆构建固化土模型，以研究防护效果。

当流体流经圆柱结构时，由于结构的存在，导致其周围的流动结构发生深刻变化，具体包括圆柱前方的下降水流、前缘马蹄涡结构、后方尾涡脱落以及圆柱两侧的流线收缩，相关的流动结构如图3所示，流动结构的改变对圆柱的局部冲刷有着重要的影响作用。马蹄涡是诱发水下垂直圆柱发生局部冲刷的关键因素，众多学者对其开展大量的研究工作。由于圆柱的存在，其上游会产生负压梯度，从而使来流边界层发生流动分离，分离的边界层会在圆柱的前缘形成马蹄涡系结构[18]。图4分别为极端高水位叠加50 a一遇波浪和0.20 m/s水流条件变截面桩基受波流作用后海床照片和极端低水位叠加50 a一遇波浪和0.20 m/水流条件变截面桩基受波流作用后海床照片。

图3 单桩基础周围流动结构Fig.3 Fluxion structure around monopile foundation

图4 极端高、低水位时波流作用后底床情况Fig.4 Erodible beds after wave and current action at extreme high and low water levels

表5为本次试验所考虑组次最大冲刷深度，从结果可见，最大冲刷深度发生于极端高水位叠加50年一遇波流情况，最大冲刷深度为0.13 m，冲刷坑沿波流方向范围约1.0 m，垂直波流方向范围约1.34 m；
对于设计高水位，由于叠加波浪和水流条件较好，因此最大冲刷深度只有0.06 m，冲刷坑沿波流方向范围1.0 m，垂直波流方向范围约1.0 m。试验结果表明，对于本次所考虑水位和波流条件，最大冲刷深度不超过桩径0.5倍，总体来讲冲刷深度较小。

表5 变截面桩基模型冲刷深度结果Tab.5 Scour depth results of variable section monopile model

毕明君[2]等针对砂质海岸单桩基础的冲刷特征与防护措施开展了试验研究，结果表明砂被间的缝隙仍然会使得泥砂被淘出，本次试验中首先将冲刷坑采用3.7~7.4 g石子回填、整平，施工期间桩基冲刷坑深度取为5 cm，冲刷坑范围为直径0.5 m，在其上覆盖拼合后尺寸 1.0 m×1.0 m 厚度为 1 cm砂被，单块砂被尺寸为0.5 m×0.5 m，图5为砂被防护铺设过程照片。在极端高水位和极端低水位及相应波浪作用3 h后，砂被具有较好的稳定性，在砂被边缘处存在一定的冲刷，但是由于砂被为柔性结构，因此砂被边缘可以自由弯曲并与海床贴合，见图6。

图5 砂被铺设后的单桩模型照片Fig.5 Picture of monopile model with sand blanket

图6 试验后的砂被防护单桩模型照片Fig.6 Picture of monopile model with sand blanket after test

近年来，砂被、块石等原材料价格上涨较快，同时砂被铺设工艺复杂、施工难度较大，风电基础防护也逐渐尝试新的防护方法。固化土防护是国内外新推出的一种新型胶凝材料固化技术，将淤泥固化后产生防护层，其强度将稳定持续增长，是替代砂被、抛石防护的优良材料，可以用于基础防护工程[14]。

由于本次模型试验为缩尺试验，比尺的存在导致无法实现对固化土浇筑过程和凝固过程的模拟，本次试验中采用水泥砂浆对固化土进行模拟，主要关注防护范围对防护效果的影响。试验中面固化土厚度为 1.0 cm，防护范围直径分别为 5.3 cm、8.0 cm和 13.3 cm。

与圆柱结构存在会导致其周围的流动结构发生深刻变化的机理类似，固化土防护结构也会对水流产生影响，因此在固化土防护边缘处会形成一定冲刷。图7和图8分别为防护范围直径为5.3 cm和8.0 cm时固化土在波流作用前后情况。固化土防护后在其边缘处引起的冲刷和淘刷，导致固化土断裂和塌陷，从而降低对桩基周围底床的防护效果。针对该情况，对于防护范围直径13.3 cm情况，在固化土边缘铺设宽度为16.76 cm的1.7~3.4 g块石。由于块石较重，局部水流流速增加或者马蹄涡并不能使块石发生失稳，在波流作用下，边缘处块石可以形成稳定缓坡对固化土形成防护，图9为该防护方案的尺寸图及波流作用后情况。可见当固化土防护范围较大，桩基所引起的水流变化对固化土以外区域影响较小，同时由于块石的防护，杜绝了固化土下部形成淘刷，从而整个防护结构较为稳定。因此建议采用固化土进行防护时，防护范围不小于5.0倍桩基直径，在固化土边缘应铺设一定量块石以抑制固化土下部淘刷。

图7 固化土防护在波流作用前后的照片（防护直径0.53 m）Fig.7 Pictures of solidified soil protection before (a) and after (b)wave and flow action (the protection diameter is 0.53 m)

图8 固化土防护在波流作用前后照片（防护直径0.80 m）Fig.8 Pictures of solidified soil protection before (a) and after (b)wave and flow action (the protection diameter is 0.80 m)

图9 固化土防护在波流作用后照片（防护直径1.33 m）Fig.9 Picture of solidified soil protection after wave and flow action (the protection diameter is 1.33 m)

本文针对单桩风机基础开展了冲刷及防护研究工作，主要关注砂被和固化土防护的防护范围和防护效果，研究主要结论如下：

1）对于本次试验所考虑直径23.3 cm桩基基础结构，最大冲刷深度为13.3 cm，桩前冲刷坑沿波流方向范围约0.5 m，冲刷情况整体较好。

2）基于砂被防护开展了试验研究，结果表明拼合后尺寸 1.0 m ×1.0 m 厚度 0.01 m 砂被可以对本工程特定的海床底质条件形成较好的防护，可为具体工程提供借鉴。

3）基于原材料上涨和施工较为复杂的情况，本次也对固化土防护措施开展了试验研究，试验结果表明固化土为硬防护，防护后边缘处仍然存在冲刷和淘刷，可以导致固化土断裂从而大大降低防护效果，因此建议防护范围不应小于5.0倍桩直径，并应在固化土防护边缘铺设块石以抑制固化土边缘处发生淘刷而影响固化土结构稳定。

猜你喜欢波流模型试验冲刷波流耦合作用下渤海及北黄海风暴潮增水及海浪过程的数值模拟海洋通报(2021年5期)2021-12-21反推力装置模型试验台的研制及验证航空发动机(2020年3期)2020-07-24偶感东坡赤壁诗词(2018年6期)2018-12-22透射槽波探测技术对煤层冲刷带的研究与应用消费导刊(2017年24期)2018-01-31槽道内涡波流场展向涡的分布特征北京航空航天大学学报(2017年11期)2017-04-23气泡对点蚀状态管道冲刷腐蚀的数值模拟山东工业技术(2016年15期)2016-12-01台阶式短加筋土挡墙行为特征的离心模型试验西南交通大学学报(2016年6期)2016-05-04巨厚坚硬岩浆岩不同配比的模型试验研究采矿与岩层控制工程学报(2015年3期)2015-12-16水库坝区冲刷漏斗的形成机理浙江大学学报（工学版）(2015年2期)2015-05-30电渗—堆载联合气压劈烈的室内模型试验浙江大学学报（工学版）(2015年7期)2015-03-01