王徽华

（江苏龙源振华海洋工程有限公司，江苏 南通 226014）

0 引 言

随着环保要求日益严苛，风能作为一种绿色能源越来越受到重视。由于海上的风况远远优于陆地，当前风力发电正逐步由陆地延伸到海上，海上风能的开发和利用已成为世界新能源发展的亮点。风电安装船作为建设海上风电场的关键装备，其开发利用也得到关注和重视。自升式风电安装船是一种全新的海洋工程船，主要用于运输和吊装海上风电设备。该船将运输船、海上作业平台、起重船及生活供给船的各项功能融为一体，可独立完成海上风电设备的运输和安装作业，因此在海洋风电安装领域得到广泛应用。该船通过将桩腿插入海底来支撑船体结构进行海上风机的吊装，桩腿入泥深度直接影响平台的吊装性能，因此开展自升式海上风电安装平台的入泥深度研究意义重大。

当前相关研究人员已针对海底土层承载力的计算开展较多工作。袁凡凡等[1]开展层状地基土的承载力计算，在迈耶霍夫和汉纳成层土地基极限承载力计算的基础上进行改进，提出多层土的极限承载力计算。杨军[2]采用数值模拟的方法开展自升式平台插拔桩土体数值模拟研究，得到入泥深度和拔桩力。张兆德等[3]采用有限元分析模型开展桩靴贯入土层的过程研究，得到桩土交界面处土体流动的动力反应特性。本文在经典土力学的基础上，采用层状土的承载能力计算理论开展海洋风电安装平台插桩入泥深度的计算，并将其与实际的施工记录及非线性有限元分析结果相对比，验证方法的准确性。

1 自升式风电平台桩靴入泥深度计算

1.1 桩靴入泥过程描述

预压加载时桩靴的贯入是一个连续的过程，贯入深度一般根据地基土承载力沿土层的变化曲线和桩靴承担的荷载来确定。地基土承载力沿土层的分布曲线是通过假想桩靴位于土层中不同的位置（见图1），由地基承载力公式计算得出的。

海上自升平台桩靴承载力计算的特点为：

1) 桩靴自身尺寸较大，在计算承载力时，需考虑桩靴承载影响范围内土层强度的变化，即按多层土（层内均质）地基或强度随土层深度线性变化的地基土进行承载力评价；

2) 桩靴贯入深度较大，需考虑地基土可能回流到桩靴上部对承载力的影响。

图1 桩靴入泥过程

SNAME[4]和ISO[5]推荐采用浅基础承载力计算公式评价桩靴承载力及贯入深度，并考虑层状地基对承载力的影响。

1.2 计算流程

1) 以基础位于层间分界面上为前提，计算各层土所能提供的承载力，计算结果用于下一步计算时判断土层的软硬和计算最终承载力；

2) 从海床泥面开始，针对不同的贯入深度，根据地基土层的情况（软硬的组合）计算承载力；

3) 根据不同贯入深度得到承载力，绘制承载力沿地基深度的变化曲线，根据桩腿预压荷载确定贯入深度，判断穿刺风险。

1.3 土层承载力计算

1.3.1 均质地基承载力计算

均质土地基极限承载力的一般计算式[6]为

式(1)中：c为地基土的黏聚力；q为基础底面以上荷载，一般为基础底面以上土重；Nc，Nq，Nγ为地基承载力系数，与地基土类型有关；sc，sq，sγ为基础形状修正系数，与基础形状及地基土类型有关；dc，dq，dγ为基础埋深修正系数，与基底深度及地基土类型有关。

1.3.2 层状地基承载力计算

1) 上硬下软的双层及多层地基的承载力采用MEYERHOF等[7]提出的层状地基承载力计算公式计算，即

式(2)中：Hi为第i层土的厚度；iφ为第i层土的内摩擦角；D为基础底面的埋深；iγ为第i层土的重度；ksi为第i层土的冲剪系数，根据上下土层的强度比和第i层土的内摩擦角，由MEYERHOF等[7]提供的计算图确定；cai为第i层土的黏结力，根据第i层土的黏聚力和上下土层的强度比，由MEYERHOF等[7]提供的计算图确定。

2) 采用式(1)和式(2)进行层状土的承载力计算。

2 穿刺破坏判断

对于存有软弱下卧层的海洋地基，上部坚实土层提供的承载力受软弱下卧层的影响。在平台吊装作业或风暴自存状态下，受环境荷载影响，桩靴对地压力可能暂时超过设定的对地压力值，或桩靴地基因处于偏心受荷状态而导致地基承载力有所下降。此时若坚实土层提供的极限承载力不足，桩靴地基可能会发生穿刺破坏。

图2为穿刺破坏示意，在桩靴地基极限承载力随桩靴贯入深度的变化曲线上出现承载力随贯入深度下降段，当桩靴承担的最大荷载超过上部土层的最大极限承载力时，桩靴地基会发生穿刺破坏。

由于环境荷载的不确定性，难以给出桩靴可能出现的最大对地压力及荷载偏心程度的确定值。一般采用在平台设计工况下得到的最大桩腿荷载乘以一定的安全系数K进行穿刺破坏的校核。

图2 穿刺破坏示意

式(3)中：Fmax为桩靴位于坚实土层中计算出的最大承载力；FDmax为设计指定（预计）的桩腿预压力。

根据SNAME[4]的推荐，当硬土层提供的极限承载力安全系数K≥1.5时，穿刺不会发生，该位置适合平台插桩。

3 计算结果

选取如东龙源试验风场进行计算分析，其中某机位的海底地质情况如下。

1) Q4-①层粉砂：主要由新近沉积的粉砂及粉土组成，含少量贝壳碎屑，层厚7.1m，饱和、松散，中等压缩性土。

2) Q4-②-2层粉砂：偶见贝壳碎屑，局部夹黏土薄层，层厚5.60m（本计算调整为3.6m），稍密，中等压缩性土。

3) Q4-③-夹1层状淤泥质粉质黏土：淤泥质粉质黏土与粉土呈互层状，其中淤泥质粉质黏土单层厚2～30mm，流塑；粉土单层厚1～20mm，淤泥质粉质黏土与粉土厚度比为2:1～5:1，见少量腐殖质及贝壳残片。层厚12.5m，属高压缩性土，工程性质差。

4) Q3-⑥-1层粉砂：见少量贝壳碎屑，局部夹粉土和薄层可塑状黏土，层厚19.80m，中密～密实，中等压缩性土。

5) Q3-⑥-3层粉细砂：含少量贝壳碎屑，局部夹粉土和薄层状黏土，本次揭露层厚14.6m，密实，属中等压缩性土。

风机采用“龙源振华2号”进行吊装，其设计桩靴压力为3.773×107N。各土层的参数见表1。

表1 各层土的参数

首先开展均匀地质下各层土的承载力计算，即不考虑各层土之间的相互影响，结果见表2。

表2 桩靴位于各层土表面时的承载力

考虑相邻土层的影响，进行层状土的承载力计算，结果见图3。

图3中实线为桩靴贯入土层后其土层极限承载力与入泥深度的变化曲线，点划线分别为3.773×107N和5.660×107N，可见该机位的入泥深度很小。考虑到表面浮泥的影响，入泥深度很小，仅为0.5m，实际的入泥深度为0.6m（根据施工记录），这与实际入泥深度相差很小，没有穿刺风险。

此外，采用非线性有限元数值模拟方法对自升式风电安装平台的桩腿入泥深度进行校核。地基土采用理想弹塑性模型描述，屈服准则采用ABAQUS内建的Mohr-Coulomb准则。有限元模型图见图4。在海底表面插桩，计算得到的载荷-位移曲线见图5。

图3 桩靴所在位置土壤承载力随桩靴贯入深度的变化

图4 有限元模型图

图5中，曲线在7.644×107N左右时急剧下降，可选用该值作为极限承载力，相应的桩靴沉降量为80mm。由图5可知，桩靴位于第一层土上时，该层土的承载力已经远大于桩腿的预压载荷，因此可以认为：由有限元计算可得的插深约为 0，考虑海底的浮泥影响，可认为插泥深度为0.5m。结果表明，本文提出的理论计算方法与非线性有限元分析结果及实际插桩入泥深度相吻合。

图5 载荷-位移曲线

4 结 语

针对自升式海洋风电安装平台的施工特点，提出一种用于评估其桩腿入泥深度的计算方法，并与实际的插桩入泥深度相比较，得出以下结论：

1) 利用该方法可较为快速地开展桩腿入泥深度评估，可为风电施工企业节省较多的成本；

2) 计算结果与实际的施工记录相比，误差较小，可为海上风电的施工提供较为准确的依据。

【 参 考 文 献 】

[1] 袁凡凡，闫澍旺，孙万禾.关于成层土地基极限承载力的计算方法[J].水利学报，2001, 32 (3): 41-45.

[2] 杨军.自升式平台插拔桩土体数值模拟及平台强度计算[D].北京：中国石油大学，2010.

[3] 张兆德，张心.土体流动对自升式平台桩靴极限承载力的影响[J].船舶与海洋工程，2015, 31 (3): 16-21.

[4] SNAME.Guidelines for site specific assessment of mobile jack-up units: Technical and research bulletin5-5A[S].1997.

[5] ISO.Petroleum and natural gas industries–offshore structures–part 4: Geotechnical and foundation design considerations: ISO 19901-4[S].2003.

[6] 钱家欢，殷宗泽.土工原理与计算[M].北京：中国水利水电出版社，1996.

[7] MEYERHOF G G, HANNA A M.Ultimate bearing capacity of foundations on layered soils under inclined load [J].Canadian Geotechnical Journal, 1978, 15: 565-572.