砂土中吸力式筒形基础冲刷数值模拟分析_赵学亮.pdf

哈尔滨工程大学学报 Journal of Harbin Engineering University ISSN 1006-7043,CN 23-1390/U 哈尔滨工程大学学报网络首发论文哈尔滨工程大学学报网络首发论文 题目：砂土中吸力式筒形基础冲刷数值模拟分析 作者：赵学亮，隋淑环，陈馨睿，任彦忠，吴成恩，宋启明 收稿日期：2021-12-11 网络首发日期：2022-11-24 引用格式：赵学亮，隋淑环，陈馨睿，任彦忠，吴成恩，宋启明砂土中吸力式筒形基础冲刷数值模拟分析J/OL哈尔滨工程大学学报.https:/ 网络首发网络首发：在编辑部工作流程中，稿件从录用到出版要经历录用定稿、排版定稿、整期汇编定稿等阶段。录用定稿指内容已经确定，且通过同行评议、主编终审同意刊用的稿件。排版定稿指录用定稿按照期刊特定版式（包括网络呈现版式）排版后的稿件，可暂不确定出版年、卷、期和页码。整期汇编定稿指出版年、卷、期、页码均已确定的印刷或数字出版的整期汇编稿件。录用定稿网络首发稿件内容必须符合出版管理条例和期刊出版管理规定的有关规定；学术研究成果具有创新性、科学性和先进性，符合编辑部对刊文的录用要求，不存在学术不端行为及其他侵权行为；稿件内容应基本符合国家有关书刊编辑、出版的技术标准，正确使用和统一规范语言文字、符号、数字、外文字母、法定计量单位及地图标注等。为确保录用定稿网络首发的严肃性，录用定稿一经发布，不得修改论文题目、作者、机构名称和学术内容，只可基于编辑规范进行少量文字的修改。出版确认出版确认：纸质期刊编辑部通过与中国学术期刊（光盘版）电子杂志社有限公司签约，在中国学术期刊（网络版）出版传播平台上创办与纸质期刊内容一致的网络版，以单篇或整期出版形式，在印刷出版之前刊发论文的录用定稿、排版定稿、整期汇编定稿。因为中国学术期刊（网络版）是国家新闻出版广电总局批准的网络连续型出版物（ISSN 2096-4188，CN 11-6037/Z），所以签约期刊的网络版上网络首发论文视为正式出版。哈 尔 滨 工 程 大 学 学 报 Journal of Harbin Engineering University 收稿日期收稿日期：2021-12-11.基金项目基金项目：国家自然科学基金资助项目（52078128）.作者简介作者简介：赵学亮，男，副教授.通信作者：通信作者：赵学亮，E-mail：.砂土中吸力式筒形基础冲刷数值模拟分析 赵学亮1，隋淑环1，陈馨睿1，任彦忠2，3,吴成恩3,宋启明4（1.东南大学 混凝土及预应力混凝土结构教育部重点实验室（东南大学），江苏 南京 211189；2.华北电力大学 河北省低碳高效发电技术重点实验室，河北 保定 071003；3.中广核（福建）风力发电有限公司，福建 福州 350001；4.福建永福电力设计股份有限公司，福建 福州 350001）摘 要：为研究吸力式筒形基础在砂土中单向流作用下的冲刷机理和影响因素，该文通过建立三维数值模型，基于重正化群湍流模型和梅叶-彼德推移质输运公式对吸力筒基础的局部冲刷进行模拟计算。结果表明：当流速较低时，吸力筒局部冲刷经历快速冲刷-缓慢发展-冲刷平衡3 个阶段；当单向流流速大于 2 m/s 时，局部冲刷深度出现波动变化，且流速越大，流速对于最大冲深的影响越剧烈；吸力筒露出床面的部分越多，冲刷程度越大，达到冲刷平衡也更快；吸力式三筒基础筒间距影响水流下游方向筒基础内外侧局部冲刷深度的发展。本文研究结果可为工程中的吸力筒基础的设计运行和维护提供参考。关键词：吸力筒基础；三筒基础；局部冲刷；高径比；筒间距；单向流；数值模拟；砂土冲刷 中图分类号：U443.22 文献标识码：A Numerical Analysis of Scour of Suction Bucket Foundation in Sand Zhao Xueliang1,Sui Shuhuan1,Chen Xinrui1,Ren Yanzhong2,3,Wu Chengen3,Song Qiming4(1.Key Laboratory of Concrete and Prestressed Concrete Structure of Ministry of Education,Southeast University,Nanjing 211189,China;2.North China Electric Power University,Hebei Key Laboratory of Low Carbon and High Efficiency Power Generation Technology,Baoding 071003,China;3.China Guangdong Nuclear(Fujian)Wind Power Co.,Ltd.Fuzhou 350001,China;4.Fujian Yongfu Power Engineering Co.,Ltd.Fuzhou 350001,China)Abstract：To investigate the mechanism and influencing factors of the scour around suction bucket foundation,this paper established a three-dimensional numerical flume to simulate the scour of the suction bucket foundation under current in sand based on the Renormalized Group k-model and Meyer-Peter Muller formular.The results demonstrate that current velocity has great influence on maximum erosion depth.Under low velocity of current,the local scour process of the mono-bucket foundation shows three stages:rapid erosion stage,slow erosion stage,and equilibrium stage.The maximum scour depth increases with more bucket exposed above the seabed.When the current velocity is greater than 2m/s,the scour process shows different three stages：rapid erosion stage,fluctuating erosion stage,and equilibrium stage.The scour pit around the backward bucket of tripod-bucket foundation shows asymmetry shape,which is greatly influenced by the bucket spacing.Results from this paper have a certain practical guide for design,operation and maintenance of the suction bucket foundation.Key words：suction bucket foundation;tripod-bucket foundation;local scour;height-diameter ratio;bucket spacing;unidirectional flow;numerical simulation;scour in sand 吸力式筒形基础是近年来从深海基础中逐渐发展起来的一种新型海上风机基础，具有筒体入泥浅、海上施工周期短、施工成本低等优点，施工过程不需要复杂的嵌岩施工技术，安装过程无噪音，对环境友好。但筒网络首发时间：2022-11-24 10:46:51网络首发地址：https:/ 哈 尔 滨 工 程 大 学 学 报 Journal of Harbin Engineering University 形基础因阻水宽度大，在位工作期间受波浪和海流作用，基础周围土体极易受到冲刷，对海上风电机组的正常运行产生非常不利的影响。在单向水流作用下，最大冲刷深度的影响因素主要有希尔兹数1-2，相对水深 hw/D(hw为水深，D 为结构物直径)3，泥沙级配4，结构物尺寸5-7，海床坡度8等。但对于结构物冲刷的研究集中于单桩基础，筒形基础作为一种新型海上风电基础，对其冲刷的研究较少。对于筒型基础局部冲刷问题的研究中，赵雁飞9建立了通过建立数值模型模拟了单桩基础和筒型基础周围的冲刷情况。于通顺10对 3 种不同比尺的复合筒型基础模型在粉砂质海床产生的局部冲刷情况进行模型试验，得到包括冲刷平衡时间、模型周围冲刷坑的范围、最大冲刷深度等的冲刷规律，并依据系列比尺模型试验原理推算得出实际工程中筒基的冲刷深度。Stroescu11对比单桩和单筒基础的冲刷形态，发现筒型基础自带防冲刷性质，在各种环境条件下，其表现类似于冲刷防护单桩。刘茜茜12通过物理模型试验研究变截面圆柱尺寸对基础局部冲刷的影响，结果发现筒型基础的最大冲刷深度与宽度随上部结构直径增大而增加，且当上部结构与主筒直径比值在 0.20.5 时，结构的防冲刷能力较好。吸力式多筒基础在较深海域能表现出更好的安全和经济性，其筒与筒之间会因布置形式的不同导致基础周围水流特性的改变，从而影响局部冲刷。Sumer等13对波流共同作用下不同布置形式的群桩进行冲刷试验，分析了相对桩间距 G/D（G 为桩间距，D 为桩径）和 KC 数对群桩基础局部冲刷的影响。Qi 等14针对双柱基础研究其桩间距和水流攻角 对于冲刷深度的影响，并得出与 G/D 和 相关的双桩基础冲刷深度公式。刘明维等15对内河框架码头大直径桩柱串列及并列四桩柱绕流进行数值模拟，研究了桩柱周围的绕流流场形态和水动力特性，分析了流速、桩间距对遮流效应的影响。由于海洋环境下吸力筒基础周围流场相对复杂，且筒形基础直径较大，模型试验因缩尺效应不能准确反映结构物周围冲刷形态，因此数值模拟成为研究吸力筒冲刷问题最有效和相对准确的方式。目前研究中，对于吸力筒基础冲刷问题的研究尚不充分，且多为单筒吸力筒基础，对于新型三筒基础的冲刷问题有待研究。本文通过 CFD 软件建立三维数值水槽，模拟砂土中单向流条件下海上风电吸力式筒型基础周围的局部冲刷，改变入射水流要素和筒基尺寸，通过分析吸力筒周围冲刷坑深度和地形，研究吸力式筒型基础的局部冲刷机理和影响因素。1 局部冲刷数值理论 计算流体动力学（computational fluid dynamics，CFD）借助计算机对包含流体流动、热传导等相关物理现象的系统进行数值计算，并通过图像显示进行分析，利用流动基本方程对流动过程进行数值模拟。1.1 湍流模型 常用的湍流模型包括普朗特混合长度模型、一方程模型、标准 k-c 模型、重正化群（renormalized group,RNG）k-模型以及大涡（large eddy simulation,LES）模型。其中 RNGk-模型通过对湍动粘度的修正，考虑到了平均流动中的湍流涡旋及旋流流动情况，从而提高了模拟精度，因此本研究选用 RNGk模型。1.2 自由表面处理方法 基于界面重构技术的自由表面处理方法（volume of fluid method，VOF）基于流体体积的百分比，每个网格（控制体积）在本身具备的流体属性外又增加了一个流体体积百分比函数 F。计算自由表面随时间的变化，需要移动网格中的体积百分比，流体百分比的变化应满足的输运方程为：0FFFFuvwtxyz(1)式中：u、v、w分别为 x、y、z 方向上的速度分量。求解输运方程得到的每个单元体的 F 值可用以确定自由表面所在的单元。1.3 泥沙冲刷模型 本文依托实际工程海域条件和工程地质环境，冲刷模型针对无粘性沙，通过预测泥沙的侵蚀、对流及沉积来估计泥沙的运动。基于床面剪切应力的希尔兹数为：()sxfgd(2)式中：g 为重力加速度；ds为沉积物中值粒径；x为沉积物密度;f为水密度.床面泥沙启动可借助临界希尔兹数cr来确定。临界希尔兹数cr基于广泛使用的 Soulsby16提出的经验公式进行计算：*0.30.055 10.021 1.2crexpDD(3)考虑到沉积物坡度对泥沙起动的影响，公式（3）可以改写为：2222cossincostansinsintancrcr (4)式中：1 32*1/xfsDgvd；v是流体粘滞系数；是河床坡度，是水流方向和坡度方向之夹角；为泥沙静止角。悬移质沙通过与流体间的对流产生输运。如果不考虑 VOF 和 FAVOR 网格处理技术(fractional-area/volume obstacle representation，FAVOR)，悬移沙的输运方程可写为：()0ssCuCt(5)式中：sC为悬移质沙的浓度；u是悬浮泥沙流速。一般情况下，沙的密度大于水的密度，因此沙会下沉或飘移，漂移的几率与浮力、拖曳力间的平衡有关。沙子 哈 尔 滨 工 程 大 学 学 报 Journal of Harbin Engineering University 的动量平衡方程和水沙混合物的动量平衡方程可写为：1SssssruKuuPFutf (6)1SssuuuPFt (7)由方程（6）、（7）可以求出漂移速度。采用梅叶-彼德（Meyer-Peter Muller）公式对推移质输运进行计算。平面单宽体积输沙率由该模型预测可得：1.53()()sfqgdcrbfss(8)式中：s为推移质系数，用于预测床面剪应力超过临界剪应力时推移质沙输运的几率，其默认值为 8.0。2 吸力筒基础冲刷数值模型 本研究模拟砂土中恒定流条件下海上风电吸力式筒型基础周围的局部冲刷。图 1 为参考实际工程吸力筒尺寸建立的基础模型图，吸力筒单筒基础和多筒基础的各筒直径 D 均为 10 m，筒长 L 为 10 m，单筒基础上下结构直径比为 0.4，多筒基础的相对筒间距 G/D为 2.5（G 为三筒基础各筒中心点之间的距离）。(a)单筒 (b)三筒 图 1 吸力筒基础模型 Fig.1 Suction bucket foundation models 2.1 模型建立 模型尺寸考虑入射波流稳定以及消除边界效应的影响，选取 X 方向总长为 220 m（22D），Y 方向宽度为 80 m（8D），Z 方向高度为 38 m。设置结构物在模型中心，考虑我国近海岸实际情况，控制水深为 25 m。为了减小初始模拟时来流水流对入口和出口边界处河床的淘刷，在上游入水口处、下游出水口设置厚2 m、宽 70 m、高 8 m 的稳水挡板，模型底部铺设置8 m 厚的床沙（图 2）。网格划分 X 方向 220 个，Y 方向 80 个，Z 方向 90 个网格，为提高计算精度，采用嵌套网格法在筒周一倍筒径处和床附近进行局部网格加密17。图 2 数值模型与边界条件布置 Fig.2 Sketch of numerical model and boundary conditions 在边界条件的设置中，单向流工况的入口选择速度作为边界，给定水流边界为均匀入流 V，水流流速 Uc为 1.5 m/s，边界的入流水深设置为 25 m。设置出口边界条件为静水压力边界 P，流体高度以控制出流水深不变为条件，设置为水面高度。计算区域左右两侧为对称边界，上边界采用刚盖假定，各物理量的垂向梯度为0。底部边界条件设置为壁面无滑移条件。初始水体设定在挡板上下游之间，上游初始水位与下游保持一致，水面无坡降，采用静水压力。设置 Z方向为重力方向，重力加速度为-9.81 m/s2。冲刷模块中，泥沙中值粒径 ds为 0.135 mm，密度为 2 670 kg/m3，水下休止角为 30。模拟为单一不可压缩流体，双精度条件。2.2 数值模型验证 为验证数值模拟结果的合理性，选取 Zhao 等5的物理模型试验进行对比验证。以模型试验中 test C 工况为验证工况，结构物尺寸（长宽高）为 40cm20cm20cm，中值粒径0.135mmsd，泥沙比重2.67xf，临界希尔兹参数0.0625cr，水深0.5mh，平均流速0.325m sU。首先对模型流场进行验证，沿水深方向设置测点，得到的水槽断面流速沿水深分布曲线，各流速测点的采样数据可按照流速对数分布规律进行拟合，即：30()ln()fsuzu zk(9)式中：fu为摩阻流速，z 为流速测点距离理论底床零点的距离；sk为 Nikuradse 等效粗糙度。图 3 对桩前行进水流断面流速进行对比，模型试验和数值模拟的断面流速重合度较高，且都满足对数分布规律，故模型流场拟合较好。哈 尔 滨 工 程 大 学 学 报 Journal of Harbin Engineering University 0.10.20.30.40.50.00.10.20.30.40.5 Zhao试验5 本文模拟z/mUc/m s-1 图 3 水流速度沿水深分布对比 Fig.3 Current velocity contribution with water depth 本文分别从冲刷形态、最大冲刷深度 2 个方面对试验和模拟的局部冲刷进行验证。从图 4 的地形云图对比可以看出，不同阶段下模型试验和数值模拟的冲刷地形结果基本一致，前期冲刷先发生在迎流侧桩角处，桩后两侧有八字形淤积。t=10 min 时，模型试验的最大冲刷深度 S=5.8 cm，数值模型最大冲刷 S=5.5 cm，两者相差 5.5%，结果吻合较好。后期桩前侧冲刷坑逐步加深扩大，但最大冲刷深度仍发生在桩角处，桩后两侧淤积也逐步后移至正后方，桩周地形变化趋势与物理模型试验更为接近。t=1 h 时，冲刷进入稳定阶段，模型试验的最大冲刷深度 S=10.8 cm，数值模型最大冲刷深度 S=10.1 cm，两者相差 6.5%，结果基本吻合。由此，数值模型计算结果能够良好的反映物理模型试验的规律。(a)t=10 min(b)t=6 h (c)t=10 min (d)t=1 h 图 4 物理模型试验与数值模拟地形云图对比 Fig.4 Comparison of model test and numerical result of scour depth contour map 3 吸力筒局部冲刷机理与影响因素分析 3.1 吸力筒局部冲刷机理分析 根据实测潮流资料，选取水深为 25 m、水流流速为 1.5 m/s 的单向流作为代表性水流条件。实际工程中吸力筒安装时，由于土塞的发生，筒体不能完全下沉至床面处，因此，设置吸力筒露出床面高度为 0.1D，计算得到不同时刻的地形图见图 5。(a)t=600 s (b)t=1 800 s (c)t=3 600 s 图 5 单筒基础冲刷地形变化图 Fig.5 Process of scour around suction bucket foundation 通过对筒周冲刷坑形态变化的观察可以发现，侵蚀主要发生在吸力筒迎流侧两侧，最大冲刷位置发生在吸力筒迎流侧沿水流方向 4590。冲刷开始阶段，吸力筒两侧床面冲刷迅速发展，在 10 min 内最大冲深为 1.246 m，已达到最大冲深的 68.9%。随着床面高度的下降，冲刷深度随时间推移发展速度逐渐缓慢，且冲刷坑形态逐渐明显呈马蹄状，边坡与最大冲刷深度位置成一定角度保持不变，冲刷坑范围向外扩展。1 800 s 开始，筒前后侧淤积逐渐明显，前侧泥沙淤积高度约0.081 m，远小于筒后侧淤积高度 0.550 m。随着筒后淤积不断后移，冲刷坑的形态基本趋于稳定，3 600 s 时冲刷深度达到 1.809 m，约为 0.18 倍筒径，发生在筒周与水流方向夹角约为 67.5处，冲刷半径达到 0.97D，此时局部冲刷趋于平衡。露出床面部分同为0.1D的吸力式三筒基础在相同水流条件下，冲刷地形变化如图 6，水流方向自左向右，即 x 轴正方向。可以看出，吸力式三筒基础的各筒冲刷坑形态和单筒相似，都为两侧冲刷，后侧淤积。床面高度(m)床面高度(m)床面高度(m)床面高度(m)床面高度(m)床面高度床面高度（m）床面高度床面高度（m）哈 尔 滨 工 程 大 学 学 报 Journal of Harbin Engineering University (a)t=300 s (b)t=1 800 s (c)t=3 600 s 图 6 三筒基础冲刷地形变化 Fig.6 Process of scour hole around tripod-bucket foundation 各筒周围冲刷的发展趋势相似，都经历快速冲刷-缓慢冲刷-冲刷平衡的阶段，但平衡时达到的最大深度存在差异。3 600 s 内，三筒基础的最大冲深发生在后筒内侧，深度达到 2.441 m，比前筒最大冲刷深度 1.809 m 增加 34.8%；后筒外侧冲刷呈波动上升趋势，最终冲刷平衡深度和前筒大致相同。通过对比可明显发现三筒基础的后侧桩基周围冲刷并非对称，内侧比外侧冲刷严重 1.35 倍。图 79 分别为 t=1 800 s 时刻 XY、XZ、YZ 筒周平面流速分布图，从不同视角反映了冲刷过程中吸力筒基础周围的流场变化。可以看出，水流遇到吸力筒基础后发生绕流，桩前上部流速减小，经过基础两侧时流速增加，竖直方向上有向下的水流，靠近海床的底部有向上的水流，由于水流剪切层分离和桩前逆压梯度的影响，在桩侧形成马蹄形涡旋，导致床面切应力增大，从而在筒基两侧形成冲刷坑。水流经过桩基时，两侧边界层分离形成尾涡，尾涡在桩基下游脱落，桩后出现泥沙淤积。(a)单筒基础 (b)三筒基础 图 7 XY 平面流速分布(m/s)Fig.7 Current contribution of XY plane(m/s)(a)单筒基础 (b)三筒基础（前筒）(c)三筒基础（后筒）图 8 XZ 平面流速分布(m/s)Fig.8 Current contribution of XZ plane(m/s)(a)单筒基础 (b)三筒基础（前筒）(c)三筒基础（后筒）图 9 YZ 平面流速分布(m/s)Fig.9 Current contribution of YZ plane(m/s)吸力式三筒基础受塔筒及支撑桁架的影响，床面上部流态紊乱，但靠近海床底部的各筒周流场变化与单筒基础类似。由于筒间存在互相影响，水流方向下游各筒基周围流场不对称，由 XY 平面流速分布图可以看出，上游水流经过前筒后发生绕流，由于后筒的阻碍，流线在后筒内侧集中并形成强度更大的桩侧马蹄涡，导致后筒内侧床面附近的水流流速局部增大，这也是形成冲刷坑不对称的主要原因。3.2 吸力筒基础局部冲刷影响因素分析 3.2.1 高径比影响 露出海床部分基础尺寸对于结构物周围流场的形态有重要的影响。考虑实际工程中吸力筒下沉位置及筒顶结构高度，分别设置吸力筒露出床面高径比 H/D为 0、0.05、0.1、0.2（特别说明，这里的高径比是指吸力筒露海床面的高度 H 与筒径 D 的比值，与通常所说的整个吸力筒长度 L 与直径 D 比值的长径比 L/D 不同）。相同入射流条件下（水流流速 Uc=1.5 m/s）不同高径比的吸力筒最大冲深随时间发展见图 10。考虑到数值结果的稳定性和计算效率，各工况模拟冲刷时长统一取 1 h，通过对比 1 h 内吸力筒的冲刷模拟情况，分析各参数对吸力筒基础的局部冲刷的影响规律。因 1 h 时冲刷还未达到完全平衡，最终的平衡冲深采用 Sheppard18提出的单桩基础的冲刷深度与时间发展的指数函数关系来推算。()1exp1expS tabtcdt(10)式中：S(t)为时间 t 对应的冲刷深度函数；a，b，c，d为系数，可通过模拟所得数据采用最小二乘法进行拟合确定。当 t 取无穷大时，即可推算得平衡冲深eSac。图 10 为高径比分别为 0，0.05，0.1，0.2 的吸力筒基础的冲刷历程线，可以看出，不同工况下冲刷发展历史过程相似，分为 3 个阶段：快速冲刷、缓慢冲刷、冲刷平衡。冲坑形成初期，冲刷深度增加速率较快，600s后增速逐渐放缓，最终冲刷深度趋于平稳，冲刷达到平衡。各工况的冲刷历时线与公式(9)拟合的R2均大于0.9。床面高度(m)哈 尔 滨 工 程 大 学 学 报 Journal of Harbin Engineering University R2为拟合优度，代表公式与原始数据的拟合程度，数值越接近 1 拟合程度越好。由此可知，式(9)也同样适用于吸力筒基础的冲刷历程拟合。由拟合函数可以推算出高径比 H/D 为 0，0.05，0.1，0.2 的吸力筒平衡冲刷深度分别为 0.820 m、0.974 m、1.174 m、2.054 m。由图 11 的高径比和冲刷平衡深度的关系曲线，可以看出平衡冲深随高径比增大而增大，且增速也随高径比提高。分析原因是，吸力筒露出床面的部分越多，阻水面积越大，水流在筒周产生的马蹄涡强度越大，冲刷程度越大。090018002700360045000.00.51.01.52.02.5Se=0.820mSe=1.174mSe=2.054m冲刷深度S/m时间历程t/s H/D=0 H/D=0.05 H/D=0.1 H/D=0.2Se=0.974m 图 10 不同高径比吸力筒的冲刷历程及拟合曲线 Fig.10 Scour process around buckets with different high-diameter ratio and fitting curves 0.000.050.100.150.200.00.51.01.52.02.5冲刷平衡深度Se/m高径比H/D 图 11 高径比和平衡冲深的关系曲线 Fig.11 Curve of high-diameter ratio vs.equilibrium scour depth 3.2.2 水流因素影响 吸力筒基础工作区域水流运动复杂，故本模拟设定了不同流速的恒定流条件，研究流速对吸力筒周局部冲刷的影响。由上一节结论可得高径比为 0.2 时，平衡冲深能达到 0.21D，冲刷坑的形态也更明显，为了更清晰地研究水流因素对于冲刷的影响，后续模拟仍选择高径比为 0.2 的吸力筒模型。由图12可以看出，流速对于最大冲深有较大影响。当 Um1.5m/s 时，冲刷趋势遵循快速冲刷-缓慢冲刷-冲刷平衡的 3 个阶段。Uc=0.5 m/s、1 m/s、1.5 m/s 时，平衡冲刷深度 Se分别为 0.108 m、0.894 m、2.054 m。当 Uc=2 m/s 时，冲刷首先进入快速发展的阶段，但当450 s 后，冲刷深度开始波动上升，且冲刷速率减缓，拟合得冲刷平衡深度 Se为 4.970 m。当 Uc=2.5 m/s 时，冲刷深度波动上升的幅度更大，平衡冲刷深度可达7.979 m。以上分析可得，当单向流流速大于 2 m/s 时，冲刷发展过程变为快速冲刷-波动冲刷-冲刷平衡的 3个阶段，图13为水流流速和平衡冲刷深度的关系曲线，可直观看出，流速越大，其对于平衡冲深的影响越剧烈；通过图 12 各拟合曲线的对比可得，流速越大，吸力筒周达到冲刷平衡所需时间越长。012002400360048000123456789Se=0.108mSe=2.054mSe=0.894mSe=7.979mSe=4.970m冲刷深度S/m时间历程t/s Uc=0.5m/s Uc=1m/s Uc=1.5m/s Uc=2m/s Uc=2.5m/s 图 12 不同流速下冲刷时间历程线 Fig.12 Scour process around buckets under different current velocity and fitting curves 0.00.51.01.52.02.53.00246810平衡冲刷深度Se/m水流流速Uc/m s-1 图 13 流速与平衡冲深的关系曲线 Fig.13 Curve of current velocity vs.equilibrium scour depth 3.2.3 筒间距影响 吸力式三筒基础的筒间距 G 是影响局部冲刷的重要影响因素。筒间距小于 1.5D（D 为单筒直径）时，与单筒基础区别不明显，不具备组合式多筒基础分析意义；当筒间距大于 4D 时，整体结构不利于实际施工，没有现实研究意义。有研究认为多筒基础的筒间距对吸力筒结构整体用钢量和水平极限承载力有较大影响，在满足承载力的基础优化上推荐长径比为 1.25 的三筒基础相对筒间距 G/D 值为 319，故本研究建立四种吸力式三筒基础，筒间距 G 分别为 2D、2.5D、3D、3.5D。哈 尔 滨 工 程 大 学 学 报 Journal of Harbin Engineering University 0600120018002400300036000.00.51.01.52.0冲刷深度S/m时间历程t/s 单筒 G/D=2 G/D=2.5 G/D=3 G/D=3.5(a)前筒 0600120018002400300036000.00.51.01.52.02.53.0冲刷深度S/m时间历程t/s 单筒 G/D=2 G/D=2.5 G/D=3 G/D=3.5(b)后筒内侧 0600120018002400300036000.00.51.01.52.0冲刷深度S/m时间历程t/s 单筒 G/D=2 G/D=2.5 G/D=3 G/D=3.5(c)后筒外侧 图 14 不同筒间距的三筒基础各筒冲刷历程 Fig.14 Scour process around each bucket with different bucket spacing 图 14 对比了不同筒间距的三筒基础在单向流（Uc=1.5m/s）下，各筒的冲刷深度随时间变化曲线。图 18-(a)可以看出各曲线发展趋势相似，重合度较高，即筒间距对于水流方向上游筒基的局部冲刷影响不明显，且与单筒基础的冲刷历程比较相似。图 14-(b)中筒间距 G3D 时各曲线重合度较高，冲刷历程相似，且1h 时的冲刷深度明显高于单筒基础，筒间距 G=2D、2.5D、3D 时，后筒内侧冲刷深度对比单筒基础分别增加 38.7%、34.8%、31.5%，但当筒间距 G3.5D 时，后筒内侧冲刷深度只增加了 9.9%。这说明当筒间距G3D 时，其对于后筒内侧的局部冲刷影响较大，冲刷深度随筒间距的提高而减小；当筒间距 G3.5D 时，其对于后筒内侧冲刷坑深度的影响迅速下降。分析原因是，水流流经后侧两筒时，筒间压缩水流导致筒内侧局部水流流速增强，而该增强效应随着筒间距的增大而逐渐减弱。由图 14-(c)中可以看出不同筒间距的筒基外侧冲刷发展相似，都经历快速冲刷-波动发展-缓慢冲刷阶段，且冲刷 1h 时冲刷坑最大深度与单筒基础相差较小，可见筒间距对后筒外侧几乎没有影响。综上分析，吸力式三筒基础筒间距 G3D 时，水流方向下游的筒基内侧的局部冲刷受筒间距 G 影响较大，冲刷坑随 G 变大而缓慢变深，但当筒间距 G3.5D 时影响变小。4 结论 1）吸力式单筒基础的冲刷较相等直径的单桩基础程度较小，冲刷坑呈对称马蹄状，3600s 时冲刷深度 S达到 1.809m，约为 0.18D，发生在筒周与来流方向夹角约为 67.5处，冲刷半径达到 0.97D。三筒基础各筒周冲刷形态接近马蹄状，但形状不对称，后侧桩内侧冲深比外侧深 1.35 倍。2）不同高径比的筒形基础局部冲刷发展历史过程相似，分为 3 个阶段：快速冲刷-缓慢发展-冲刷平衡。但最大冲刷深度和平衡时间有所差异，吸力筒露出床面的部分越多，冲刷程度越大，达到冲刷平衡更慢。3）单向流流速较低时，冲刷趋势遵循快速冲刷-缓慢冲刷-冲刷平衡的 3 个阶段，而当单向流流速较高大于 2m/s 时，冲刷发展过程变为快速冲刷-波动冲刷-冲刷平衡的 3 个阶段，且流速越大，其对于最大冲深的影响越剧烈。4）吸力式三筒基础的筒间距 G 对水流方向下游的筒基的局部冲刷影响较大，会扩大筒基内侧冲刷坑的深度和宽度，且该影响随着 G 的扩大不断变弱，当 G大于 3.5D 时，其影响可忽略不计。参考文献：1 SUMER B M,CHRISTIANSEN N,FREDSE J.The horseshoe vortex and vortex shedding around a vertical wall-mounted cylinder exposed to waves J.Journal of Fluid Mechanics,1997,332:41-70.2 BAYRAM A,LARSON M.Analysis of scour around a group of vertical piles in the field J.J Waterw Port Coast Ocean Eng-ASCE,2000,126(4):215-220.3 杜升涛,梁丙臣,鲁有祥.淹没式箱形结构的几何参数对局部冲刷影响研究C/第十四届全国水动力学学术会议暨第二十八届全国水动力学研讨会文集(下册).中国吉林,2017,743-755.4 MELVILLE B W.The physics of local scour at bridge piersC/Proceedings 4th International Conference on Scour and Erosion(ICSE-4)November 5-7,2008,Tokyo,Japan,2008:28-40.5 ZHAO M,ZHU X,CHENG L,et al.Experimental study of local scour around subsea caissons in steady currents J.Coastal Engineering,2012,60:30-40.6 刘超,孙文,张积乐,等.海上风电场单桩基础结构冲刷实验研究 J.太阳能学报,2016,37(2):316-321.7 YAO W,AN H,DRAPER S,et al.Experimental 哈 尔 滨 工 程 大 学 学 报 Journal of Harbin Engineering University investigation of local scour around submerged piles in steady current J.Coastal Engineering,2018,142:27-41.8 程永舟,唐雯,李典麒,等.波浪作用下斜坡沙质海床上桩柱周围局部冲刷试验研究 J.水科学进展,2018,29(2):260-268.9 赵雁飞.海上风电支撑结构波浪力及基础冲刷的三维数值模拟研究 D.天津大学,2010.10 于通顺,练继建,齐越,等.复合筒型风电基础单向流局部冲刷试验研究 J.岩土力学,2015,36(4):1015-1020.11 STROESCU I-E,FRIGAARD P,FEJERSKOV M.Scour Development Around Bucket Foundations J.International Journal of Offshore and Polar Engineering,2016,26(1):57-64.12 刘茜茜,陈旭光,冯涛,等.波浪作用下桶形基础冲刷特性试验研究 J.海洋工程,2019,37(6):104-113.13 SUMER B M,FREDSE J.Wave scour around group of vertical piles J.Journal of Waterwa