卢钦先,王洪庆,刘晓建,郭辉群*
(1.中国能源建设集团广东省电力设计研究院有限公司,广东 广州,510663;2.珠江水利委员会珠江水利科学研究院,广东 广州,510610)
海上风电是新型可再生能源开发的前沿领域,最具有规模化发展潜力,受到沿海国家和区域的广泛关注。海上风电桩基础长期处于波浪、水流往复作用的海洋环境下,常伴有高速绕柱涡动,巨大水流挟沙力,形成局部冲刷坑。海上风电桩基础外形呈细长体,随着冲刷深度增加,强风强流作用下桩柱细长体的柔性自振频率会接近周边海洋动力的运动频率,极有可能形成共振,影响结构安全性。
许多学者对桩基冲刷进行了研究,主要集中在单纯水流和单纯波浪的情况[1-4]。波流作用下的桩基冲刷涉及到波浪、水流、桩柱与邻近底质之间复杂的相互作用,研究成果较少[5-6]。实验表明:波流联合动力并不是简单的叠加,不同动力组合情况下结构物周围涡流系存在较大的差异[7]。对于桩周最大冲深,部分学者们给出了相反观点。Qi等人[8]指出水流或波浪单独作用下的桩基周围冲刷坑深度明显小于波流共同作用下的。Sumer 等人[9]指出联合的冲刷深度小于纯水流的,但在大于纯波浪的情况下,该过程受KC 数(KC=UwmT/D,其中,Uwm是波浪最大水质点速度,T是波周期,D是桩直径)和速度比[Ucw(=Uc/(Uc+Uwm),其中,Uc是波流联合动力中的水流流速]的影响,并指出当Ucw>0.7 时,联合的冲刷深度接近单一水流的,因此,波流共同作用下单桩基础局部冲刷机理有待进一步研究。作者拟以珠海金湾风电基础受波流共同作用为例,通过物理模型实验,研究波浪、海流和波流共同作用下桩基周围的局部冲刷特征,以期为海上风电桩基设计相关冲刷参数的设定提供借鉴。
金湾风电场位于珠江口万山岛和高栏岛之间,海洋环境十分复杂,受波浪、海流、风暴潮等多种动力因素影响,其局部冲刷特性难以用现有的经验公式进行预测,潮汐属于不正规半日混合潮型。金湾风电场处潮汐特征选取三灶潮位站作为参证站,据自1965 年至今的验潮资料,历史最高潮位3.23 m,平均高潮位0.39 m,最低潮位-1.97 m,平均低潮位-0.70 m。在100 a 一遇高潮位和SE 向强浪作用下,有效波高在7.11~9.44 m 之间达到最大波高,W 向浪相对较弱。由于风电场范围水深差异较大,各点波况也明显不同。考虑工程区域波浪特征,5 a 和50 a 一遇海流设计流速分别为1.08、1.19 m/s,典型波浪要素见表1。
表1 试验原型波浪要素Table 1 Wave parameters of test prototype
根据工程结构和动力条件,采用大比尺物理模型试验,同时满足水流运动相似、波浪运动相似和泥沙运动相似。比尺相似要求为:重力相似(流速比尺),阻力相似(糙率比尺)λn=时间比尺流量比尺其中,λl为水平比尺,λh为垂直比尺。
模型中的波浪采用正态,根据模型和水槽造波能力,几何比尺确定为1:50,水流动力比尺(流速比尺λv、时间比尺λt)采用1:7.07,波要素比尺(波长比尺λL、波高比尺λH)采用1:50,(波速比尺λC、波周期比尺λT)采用1:7.07。
模型沙的选择以水流泥沙起动相似为主,根据风电场所处海域土层采样资料分析,区域主要以淤泥质土等软土层为主,中值粒径介于0.01~0.03 mm,淤泥层塑性指数(Ip)为18.9。通过泥沙起动流速计算公式为式(1)或式(2),求得底质起动流速情况。
式中:h为水深,m;d为泥沙中值粒径,取d=0.01~0.02 mm;Δ为糙率,d<0.5 mm 时,取Δ=0.5 mm;εk为泥沙黏结力参数,原型沙取εk=2.56 cm3/s2;δ为薄膜水厚度,取0.21×10-4cm;γ为泥沙,约2.65 t/m3;γs为水容重,约1 t/m3。
对于黏性土,临界起动切应力与土的塑性指数有关,可采用Smerdon公式计算:
式中:Ip为土的塑性指数;τc为临界起动切应力,N/m2。
起动摩阻流速可由式(4)进行计算:
式中:ρ为水的密度。
原型表层底质起动摩阻流速计算结果(u*c=0.043 8 m/s)可换算到平均流速为:
式中:u为垂线平均流速;u*c为摩阻流速;C为谢才系数,其中,R≈h(水深),n为糙率,g为重力加速度。
学者们将波浪条件和纯流Shields 泥沙起动标准曲线相结合,指出振荡流作用下床面切应力可利用波浪摩擦系数计算[10-11],床面剪应力τcm的最大瞬时值,定义为:
式中:umc为近底波浪水质点最大水平速度;fw为波浪摩擦系数,由Jonsson图表确定。
对于本研究选定的淤泥土层,由于粒径较细(d50=0.01~0.03 mm),而起动切应力又较大(u*c=0.043 8 m/s),处于光滑紊流区,fw可由式(7)计算。
fw= 0.09RE-0.2。 (7)
由式(1)~(7)计算,底质的起动流速介于1.03 ~1.81 m/s 之间,其中,最小为1.03 m/s(换算到模型中约为14.5 cm/s)。经过多次测试,本研究模型沙采用处理过的精煤屑,容重γs=1.35 t/m3,d50为0.4 mm,模型沙铺设简便,长期浸水不会板结,水下扰动不会密实,能够达到泥沙重复性试验的要求。该模型沙在试验水深时,起动流速为12~16 cm/s,满足原型沙起动准则,可用于后续冲刷试验研究。
本物理模型试验在珠江水利科学研究院里水科学试验基地的咸潮风浪流水槽中进行,水槽的长76 m×宽1.2 m×高1.5 m,实验布置如图1 所示。造波机位于水槽右端,水槽左端设有1:8 的消浪斜坡,斜坡上覆盖多孔材料减少波浪反射。同时,为了尽量避免反射波对试验区波浪的影响,试验中每造100个波后,造波机停止运行,待水面稳定后再重新造波,如此循环往复。水槽底的内部装有潜水泵,实验中通过控制潜水泵转机频率控制流量,水流稳定后使用流速仪测量水槽底部流速,不断调整转机频率,直到达到目标流速。
图1 物理模型试验布置(单位:m)Fig.1 Layout of physical model test(unit:m)
金湾风电场工程单桩基础柱径拟为7.5、8.0、8.5 m,因海洋环境中结构物的存在,改变了其周围的波浪和水流条件,桩基周围的海床常常出现淘刷现象。为了确定波流作用下桩周局部冲刷过程,首先针对原型7.5 m 桩基,在50 a 一遇波流条件下开展实验研究。根据现场经验,1 次台风过程约为3 h,则根据比尺缩放关系换算到试验中约为25 min,以此作为1 次冲刷过程标准。此后,每隔一段时间对桩前、后、两侧等4个测点高程进行量测,原型1~15 h 后,桩周测点和桩基前后冲刷深度测量结果如图2~3 所示,图中x为沿水槽方向,z为沿水深方向,θ为沿桩表面角度,其中,θ=0为迎水侧点。从图2~3 中可以看出,波流作用原型6 h后,冲刷深度基本达到稳定,原型第12 h和15 h各点冲刷深度仅有少量变化。因此,后续试验认为波流作用原型15 h 后冲刷坑达到稳定状态,其中,原型15 h 换到算模型试验约为2 h,以此作为海床稳定冲刷标准。
图2 桩周测点不同时刻冲刷深度Fig.2 Scour depth of measuring points around piles at different times
图3 桩基前后不同时刻剖面冲刷深度Fig.3 Scour depth of profile at different time before and after pile foundation
波浪、水流和波流耦合等多种动力下,桩基冲刷过程明显不同,为了掌握不同动力下水沙运动差异,首先采用高速相机对其进行捕捉(以7.5 m桩径为例)。
3.2.1 纯水流作用
在恒定流作用下,水流对桩基周围沙床的切应力受马蹄涡和两侧束水的影响显著增大。切应力大于泥沙的临界起动切应力时,发生冲刷,随着冲坑深度的增大,冲刷坑的深度不再变化,达到冲刷平衡,此时形成的桩基冲刷为清水冲刷,冲刷过程如图4所示。
从图4中可以看出,桩基的上游迎水面和两侧先发生冲刷,在纯流作用下,冲刷逐渐向下游发展,直至冲刷深度不再发生变化,冲刷坑的上游沙坑坡面较陡,下游沙坑坡面较缓,呈现不对称的“勺型”。桩周冲刷达到平衡后,测得柱前最大冲刷深度达到11 cm(即现场为5.5 m),约为0.73D;桩侧最大冲刷深度为10.5 cm(即现场为5.25 m),约为0.7D。
3.2.2 纯波浪作用
波浪作用下桩基周围的泥沙起悬过程如图5所示。波浪作用时,桩柱前产生反射,同时,在柱后掩护区形成衍射。从图5中可以看出,漩涡运动对大直径桩的冲刷过程影响较小,当波浪反射后,近底层流速增加,泥沙逐渐起动,反射和绕射导致水质点运动不对称性,引起泥沙的净输移,冲刷坑并没有发展至环绕全部基础,桩基周围形成明显沙纹,整个过程中冲刷坑并不显著。
图4 流作用下桩周泥沙起悬过程Fig.4 Sediment suspension process around pile under current
图5 波作用下桩周泥沙起悬过程Fig.5 Sediment suspension process around pile under wave
3.2.3 波浪水流共同作用
风机基础结构在波浪、水流共同影响下,桩周海床瞬时冲刷过程如图6所示。由于桩基的阻水作用,流体的行进速度变为零,转化为驻点压力。驻点压力导致水流向下流动,并且还导致桩两侧绕流的侧向流速明显增大,引起泥沙起悬。
图6 波流作用下桩周泥沙起悬过程Fig.6 Sediment suspension process around pile under wave and current
长时间的波流联合动力作用后,桩基根部出现明显的冲刷坑,如图7所示。桩前及两侧的冲刷程度大于背流/浪面的,冲刷形状呈现“勺”状,并伴有广泛分布宽度较窄的沙纹,桩背水面的中轴线附近为淤积区。桩基局部冲刷特征和水流、波浪动力特征相对应,兼具单一波/流两者引起的冲刷特性。桩周冲刷达到平衡后,测得柱前最大冲刷深度为9.50 cm(即现场为4.75 m),约为0.63D;桩测最大冲刷深度为9.20 cm(即现场为4.60 m),约为0.61D。从图7 中可以看出,波流作用下的冲刷值小于纯流作用下的,表明:波浪起到了一定回填作用。实际情况中,桩基所处的海洋环境波流一般共存,因此,后续研究中以波流动力为基础开展相关研究工作。
图7 波流作用下桩基近区冲刷形态Fig.7 Scour pattern around pile foundation under wave and current
基于超声波地形测量仪,对桩周海床冲淤形态进行测量,部分结果如图8 所示(已换算成原型,D=7.5m,h=14.0 m)。从图8 中可以看出,5 a一遇波流情况下的冲刷深度和冲刷范围均小于50 a一遇波流情况下的。其中,一次冲刷,仅桩周出现了较大冲刷,经过多次冲刷海床稳定后,桩周冲刷深度进一步增加,同时桩后两侧冲刷深度也明显增大。冲刷坑主要围绕在桩基周围,冲刷最深处在单桩迎水面和单桩结构后两侧,单桩背流侧冲坑略浅,尾流区有淤积。其他实验工况下波流引起的冲刷地形形态变化与该情况没有本质差异,仅冲刷幅度略有不同。表明:波浪在桩基局部冲刷过程中主要起掀沙作用,当冲刷深度较大时,波浪对冲刷坑的下切作用明显减弱。结合水体产生的周期性波动力,底部速度增大,迫使泥沙起动。引起局部冲刷的主要原因为:①桩柱阻水在柱前形成局部壅水,产生垂向水流下切床面;②波动力掀起泥沙输移;③桩柱两侧和后部地形的冲刷主要受马蹄形漩涡和尾流漩涡的影响。
图8 桩基近区冲刷形态Fig.8 Scour pattern around pile foundation
本试验结果表明:波流共同作用下桩基近区冲刷范围主要位于单桩迎水面和单桩结构后两侧,其中,迎水面冲刷最为显著。柱径为7.5 m、水深为13 m 的桩前、桩后床面地形剖面特征如图9 所示。基于比尺缩放模型试验成果,已换算到原型。从图9中可以看出,随着波流过程桩基冲深逐渐增大,而冲坑坡度未明显改变,柱基前后冲刷坡度约1:10~1:1,且桩前的冲刷坡度远大于桩后的。
图9 桩周冲刷剖面Fig.9 Scour profile around pile
对试验工况进行统计,可得不同工况下的桩周海床最大冲刷深度,见表2。
表2 50a一遇波况各工况最大冲刷深度Table 2 Maximum scour depth under each case in a once-in-50-year wave
由表2 可知,在现场柱径8.5 m 和水深25.0 m组合时,50 a一遇波流情况下桩基近区海床冲刷最为严重,在柱径8.5 m、水深25.0 m 组合的50 a 一遇波流情况下平面和三维冲刷形态如图10~11 所示。从图10~11 中可以看出,受尾涡脱落影响,桩基侧后方海床冲刷较明显,侧后方冲刷显著区与主流方向的夹角在30°~66°。沿波流方向上,冲刷深度达到-1.0 m 以上的范围最远距桩基中心约为35 m,冲刷深度达到-2.0 m 以上范围的约15 m,冲刷深度达到-3.0 m 以上范围的约9.0 m。垂直于波流方向上,冲刷深度达到-1.0 m 以上范围最远距桩基中心约为25 m,冲刷深度达到-2.0 m以上范围的约22.5 m,冲刷深度达到-3.0 m 以上范围的约7.5 m。
图10 桩周平面冲刷形态Fig.10 Plane scour pattern around pile foundation
图11 桩周三维冲刷形态Fig.11 Three-dimensional scour pattern around pile foundation
根据工程地质和动力条件,通过物理模型试验对珠海金湾风电场风机基础建设后的局部冲刷进行研究,得到结论为:
1)纯流作用下,受马蹄涡和两侧束水的影响,桩周出现明显冲坑,形成不对称的“勺型”。纯波作用下,桩基周围的冲刷坑不能发展至环绕全部基础,桩基周围出现明显沙纹,冲刷坑较为不明显。波流共同作用下,冲刷形态兼具这2种特征,冲刷坑形状似“勺”,并伴有广泛分布宽度较窄的沙纹,但本海域波流作用下的冲刷深度小于纯流的,表明波浪起到了回填作用。
2)冲刷坑主要围绕在桩基周围,冲刷最深处在单桩迎水面和结构后两侧,单桩背流侧冲坑略浅,尾流区有淤积,柱基前后冲刷坡度约为1:10~1:1,且桩前的冲刷坡度远大于桩后的。受尾涡脱落的影响,桩基侧后方海床冲刷较明显,侧后方冲刷显著区与主流方向的夹角在30°~66°。
3)考虑最不利情况,在沿波流方向上,冲刷深度-1.0、-2.0、-3.0 m以上范围最远距单桩桩基中心约为35、15、9.0 m,即4.12D、1.76D、1.06D。在垂直于波流方向上,冲刷深度-1.0、-2.0、-3.0 m 以上范围最远距单桩桩基中心约为25、22.5、7.5m,即2.94D、2.65D、0.88D。