海上吸力基础式风机结构动力响应研究.pdf

海洋工程石油工程建设D0I:10.3969/j.issn.1001-2206.2023.05.002海上吸力基础式风机结构动力响应研究陈实，李书兆，邳帅，顾雨晨，王懿31.中海油能源发展股份有限公司清洁能源分公司，天津116 0 2 62.中海油研究总院，北京10 0 0 2 83.中国石油大学（北京）安全与海洋工程学院，北京10 2 2 4 9摘要：由于海上风机受到风浪流荷载的长期作用，因而其动力响应问题成为风机结构设计的关键问题。针对上述情况，以5MW吸力基础固定式海上风机为研究对象，充分考虑风机基础与土的相互作用等非线性荷载，利用有限元方法建立了综合考虑空气动力、水动力和土壤约束力作用的海上风机整体动力耦合分析模型，在多荷载工况组合下进行了海上风机结构的动力响应分析。研究结果表明，海上风机整体结构位移最大区域在塔筒顶部，应力最大区域在吸力基础与风机塔架的连接处。通过对线性叠加法和Turkstra 准则的对比分析验证可以看出，采用Turkstra 准则能更准确地获得对风机结构动力响应最不利的荷载组合方式。关键词：风浪流载荷；动力响应；吸力基础；有限元法Research on dynamic response of offshore wind turbine structure with suctionfoundationCHEN Shi,LI Shuzhao,PI Shuai,GU Yuchen,WANG Yi31.Clean Energy Branch,CNO0C Energy Development Co.,Ltd.,Tianjin 116026,China2.CNO0C Research Institute Co.,Ltd.,Beijing 100028,China3.College of Safety and Ocean Engineering,China University of Petroleum(Beijing),Beijing 102249,ChinaAbstract:Because of the long-term effect of wind and wave loads on offshore wind turbines,the dynamic response of offshore windturbines has become a key problem in the structural design of the wind turbine.In view of the above situation,the stationary offshorewind turbine with 5 MW suction foundation is taken as the research object,and the nonlinear load such as the interaction between thewind turbine foundation and the soil is fully considered.The overall dynamic coupling analysis model of the offshore wind turbine isestablished by using the finite element method,which comprehensively considers the aerodynamic,hydrodynamic,and soil bindingeffects.The dynamic response analysis of offshore wind turbine structure is carried out under multi-load conditions.The results showthat the maximum displacement area of the whole structure of the offshore wind turbine is the top of the tower,and the maximum stressarea is the connection between the suction foundation and the wind turbine tower.Through the comparison and analysis of the linearsuperposition method and Turkstra criterion,it can be seen that Turkstra criterion can more accurately obtain the most unfavorable loadcombination mode to the dynamic response of the wind turbine structure.Keywords:wind and wave loads;dynamic response;suction foundation;finite element method海上风机系统是一个受气动荷载、水动力荷载和控制荷载等多荷载作用的强耦合系统，相对于仅受风载的陆上风机而言，海洋环境和地质条件更为复杂。在服役期间，风、波浪及海流等环境荷载对海上风机系统的作用，对其结构的影响极大。因此利用风机结构动力响应分析来确保复杂环境条件下海上风机结构安全是海上风电机组稳定运行的关键环节。对海上风机结构进行动力响应研究首先要解决动态环境荷载的时域仿真问题，风荷载、波浪荷载和海流荷载是风机结构相关研究最多也是最主要的部分。在对风机结构进行动力响应分析时，大多针对多种工况和荷载组合情况进行动力特性的分析研究。SCHWARTZS等通过海洋环境荷载的疲劳荷载谱对海上风机塔架进行了频域分析和时域分析，将风机整体视为一个由叶片和塔架组成的多体动力系统，在风荷载作用下，风轮叶片挥舞振动与塔架的前后振动之间存在动力耦合效应。王明超等2 对比了p-y曲线法、m法和接触单元法等多种桩土相互作用模型下的单桩式海上风机，在风、浪荷载单独作用和联合作用下正常工况的运动响应情况。张小玲等13 根据Turkstra准则082023年10 月海洋工程陈实等：海上吸力基础式风机结构动力响应研究将荷载进行叠加组合，研究对单桩风机结构最不利的荷载组合方式。沈玉光4 对3MW筒型基础风机结构进行动力响应分析时，通过线性叠加法和Turkstra准则构建了多个荷载组合，以此来分析比较结构的最大响应位移和应力。Turkstra准则是美国国家标准A58推荐的多荷载组合方法，该方法基于不同类型的荷载不会同时达到最大值的假设，目的是通过考虑不同荷载的作用方式和方向，来确定一个适当的荷载组合，以保证结构在不同荷载情况下的安全性和可靠性。该方法组合形式较为充分，同时是ABS、D NV G L 等机构在海工设计计算中推荐的荷载组合方法5-6 。由上述分析可知，现有对于海上风机的动力响应分析研究多采用近似或解耦模拟，该方法对于常规海上结构物具有一定的适用性，研究对象通常为海上单桩式风机基础，但对于动态耦合较复杂的海上风机结构显然是不够准确的7 。在海上风机建模方面，目前对于固定式海上风表15MW的NREL风机设计参数类别设计参数转子质量/110机舱质量/240塔筒质量/347.46轮毂中心高度/m90部。塔筒参数按照5MW的NREL风机参数取值，单桩和吸力基础结构的详细参数见表2。塔筒、单桩和吸力基础全部选用Q345钢，材料参数见表3。基于Mohr-Coulomb土体本构模型建立吸力基础与土相互作用模拟模型。整体模型的网格划分单元类型均使用C3D8R实体单元。Y机结构的动力响应分析多采用分离式建模的分析方法8-1。虽然该分析方法仿真时间短，但在一定程度上忽略了多荷载耦合作用下风机上部结构和基础结构之间的耦合效应，且由于风机属于结构自身负有激励的海上结构物，所以对其进行动力响应分析不仅要考虑复杂海洋环境荷载的影响因素，还要合理地考虑风机运行工况下所产生的运行荷载12-14 。综上所述，本文在系统研究吸力基础式海上风机所受主要环境荷载的模拟及其计算方法基础上，提出一种综合考虑风机塔筒、单桩和基础结构的海上风机动力耦合分析方法。运用此方法针对吸力基础固定式海上风机结构进行了动力响应分析，研究其在多种荷载组合方式下的动力特性，并探究了对结构响应的最不利荷载组合情况，对固定式风机一体化设计具有重要的参考价值。1吸力基础风机模型建立本文选取美国可再生能源国家实验室（Na-tional Renewable Energy Laboratory，NR EL）开发的5MW基准风机模型，主要基本参数如表1所示。采用ABAQUS软件建立了考虑风机塔筒、单桩和吸力基础结构的有限元分析模型，如图1所示，风机叶片和机舱简化为质量块加载在风机顶ZX（a）吸力桶基础网格划分风荷载海平面浪荷载一流荷载（b）吸力桶基础荷载施加图1吸力基础模型及加载方式表2 单桩一吸力基础结构参数类别设计参数/m单桩外径6单桩壁厚0.05吸力基础外径20吸力基础高度20塔筒高度87.6塔筒底端直径6塔筒底端厚度0.027表3Q345钢材材料参数密度p/弹性模量E(kg/m)/GPa7850210类别吸力基础高度吸力基础壁厚肋板内侧高度肋板外侧高度塔筒质心位置塔筒顶端直径塔筒顶端厚度泊松比屈服强度/MPa0.3295第4 9 卷第5期09设计参数/m200.0550.538.2343.870.019海洋工程石由工程建设在动力分析步中，通过设置2 0 0 s的时间步长对结构进行动力响应分析；定义相应的边界条件和荷载设置，将塔筒结构进行分段处理以施加风荷载，水动力荷载作用于单桩结构上，其中对于波浪荷载的施加应考虑波高，将荷载作用区域进行相应的扩展。根据材料参数对风机结构和土体赋予相应的材料属性。2海上风机整体耦合分析方法2.1多荷载作用下海上风机结构耦合分析方法在海洋环境中，固定式风机主要受气动荷载和水动力荷载作用，且通过海底土层的反力支撑整个结构物。在这个过程中空气动力与水动力同时作用在结构上，使得结构产生加速度和速度，引起位移与应力的变化，结构的运动反过来又影响外部环境荷载的计算，这是一个空气动力、水动力和结构动力之间的强耦合过程。因此，海上风机结构的动力响应模型是一个空气动力-水动力-非线性土体-弹性体耦合的模型，其气动荷载与水动力荷载联合作用下风机结构的耦合运动方程如式（1）所示。M(i)+c(u)+K(u)=Fa(t)+(Filyla(t)+(Fual+(Fe)式中：【M为海上风机结构质量矩阵，kg；(i)为结构的加速度矢量，m/s；C 为海上风机结构阻尼矩阵，N/（m/s)；(i)为速度矢量，m/s；K 为海上风机结构刚度矩阵，N/m；u)为位移矢量，m；(F.（)为作用在风机结构上的随机气动荷载时程，N；(Fi y d（t)为作用在基础上的随机水动力荷载时程，N；(F)为地基土对基础结构的作用力，N；(Fc)为重力荷载，N。上述动力方程包含的环境荷载具有很强的随机性，且基础与土相互作用对于风机结构动力响应具有显著影响，因此本文将基于ABAQUS的动力隐式分析方法，采用直接积分法进行动力方程的求解，建立具有可行性的海上风机结构动力响应耦合分析方法。2.2海上风机载荷计算及模型建立本文对水深为50 m的阳江海域进行研究，取表4 所示的风机运行工况与极端工况参数作为后续模拟的基础参数。表4海况环境参数平均风速风机有义波高波浪周期海流流速海况/(m/s)运行11.4极端51.32.2.1随机风荷载计算风荷载主要指由气流运动作用在结构上的动压力，现实中风对于风机结构的作用是一个随机变量，主要包含风机顶部叶片风荷载和塔架上风荷载。1）风机顶部气动荷载计算。风机顶部叶片结构复杂，相关研究中通常采用经典叶素-动量理论进行风机顶部气动荷载计算。通过查阅5MW风机设备相关设计手册Definitionofa5-MWRefer-ence Wind Turbine for Offshore System Develop-ment0 可知，在平均风速超过风机额定风速时，机舱中的变桨距系统会调节叶片桨距角，对输出转矩和功率进行控制，使顶部结构减小受力，在平均风速达到极端风况时，风机进入停机空转阶段，对结构起到保护作用。因此在不同风况下，风机顶部受载可总结为：T逐行项部=F风荷(领定)+F选行荷较（2）(1)(3）式中：T运行顾部为风机运行工况下的顶部所受荷载，N；F风有载（M完）为风机受到的额定风载荷，N；F进行格为风机运行时产生的荷载，N；w 油顾部为风机在极端工况下的顶部所受荷载，N：F及格数（股）为风机在极端工况下受到的风荷载，N。本文的相关研究以5MW的NREL基准风机为例，基于稳态风和脉动风的模拟方法，结合叶素-动量理论，对其进行气动荷载计算，计算结果如图2 所示。2）风机塔筒风荷载计算。根据NB/T10105一2018海上风电场工程风电机组基础设计规范I中的设计要求，在对风机塔筒进行风荷载计算时，考虑到塔筒结构的高特性，需将其沿高度方向分段进行风荷载计算，在对结构进行力学分析时，塔筒风荷载的施加为将每段的风荷载进行叠加。作用在风机塔筒上的风荷载其计算公式如式（4)所示。状态/m运行4.2停机9.6均值/s7.89.2/(m/s)0.81.0102023年10 月海洋工程陈实等：海上吸力基础式风机结构动力响应研究式（5)所示。750TD2700dF=p4650式中：dF为波浪载荷，N；p 为流体密度，kg/m；N/薄桌60055505004504003500dzD为结构直径，m；CM 为惯性力系数；u,为波浪速度，m/s；z 为塔筒高度，m。将单位长度为n的波浪力沿桩长度方向积分，计算总的波浪力，见式（6)。50100时间s(5)150200F=dF0(6)（a）运行风况下气动荷载时程曲线180150120N/90603000（b）极端风况下气动荷载时程曲线图2风荷载计算2式中：F,为作用在风机塔筒上的风荷载，N；p 为空气密度，kg/m，取1.2 2 6；C,为拖曳力系数，塔筒为圆柱截面，取0.4 7；A为迎风面积，m；(t)为风速，m/s；t 为风荷载作用时间，s。将高度为8 7.6 m的塔筒结构由下至上分为9 段进行处理，其中第9 段塔筒为8 0 8 7.6 m，取每段中点处的高度为风荷载时程计算的水平高度，分别为5、15、2 5、35、4 5、55、6 5、7 5、8 3.8 m。2.2.2非线性波浪荷载计算波浪的周期、波高等因素关系着结构的动力响应，本文采用谐波叠加法描述波浪的非线性特征，再结合Morison方程对其进行波浪荷载的时域性计算。对于桩结构的波浪荷载计算采用工程上典型的Morison方程，结合上述研究所得的波浪运动时程曲线，即可得到作用在塔筒上的波浪荷载时程曲线。作用在单位长度塔筒上的波浪力表达式如式中：F为总的波浪力，N。基于Morison公式和模拟的随机波高时程曲线，以5MW海上风机为例，得到波浪荷载时程曲线，以极端工况为例，如图3所示。32NN/薄鼎050100时间/s150200（4）-20（a）运行海况下波浪荷载时程曲线64NN/薄20-20（b）极端海况下波浪荷载时程曲线图3波浪荷载计算2.2.3海流荷载计算海流荷载计算结果见图4。海流力的大小与流速与水深的变化直接相关，工程上针对无实测数据的情况下，通常以ABS与DNV推荐使用的公式计算海流流速15，即：5050100时间s100时间/s150150200200第4 9 卷第5期11海洋工程石油工程建设+Umh一Hh三HUmo=:V式中：h为距海底面h处海流的流速，m/s；U t o 为海面表层的潮流速度，m/s；U m o 为海面表层的风海流速度，m/s；h 为流质点距海底面的计算深度，m；H 为目标海域水深，m；为风流系数，0.0 2 4 0.05；V为海面10 min内平均风速，m/s。806040200-20400（a）运行海况下流荷载时程曲线150100500-50-1000（b）极端海况下流荷载时程曲线图4 流荷载计算3风浪流荷载耦合作用下结构动力响应分析3.1线性叠加法本文选取的第一种多荷载作用组合方式为线性叠加法，其在时域上将风浪流荷载进行线性组合，将叠加后的环境荷载直接作用到结构上，不考虑其他特殊情况。线性叠加法组合形式见式(11)。(7)(8)（9）(10)50100时间/s150F线(t)=F岚(t)+F%(t)+F流(t)式中：F线(t)、F凤(t)、F浪(t)和F流(t)分别表示线性叠加荷载、风荷载、浪荷载和流荷载，N。风浪流荷载线性组合作用下风机结构响应情况如图5所示。0.70.60.50.40.30.20.10.0160.0140.012150200100150时间/s(11)一运行工况.一极端工况50100时间/s（a）塔筒顶部位移时程曲线0.010三0.008位0.0060.0040.0020.000-0.002图5基于线性叠加法的风浪流荷载耦合作用下风机动力响应情况综上所述，从以风浪流荷载线性叠加进行动力响应分析所输出的塔筒顶部位移和吸力基200础顶部位移时程曲线来看，可以判断出风机结构在极端工况下的运动响应更为剧烈。为进一步对该结论进行验证，同时分析多荷载耦合作用下风机结构的受力情况，结合吸力基础结构受力主要集中在桩体的顶部与肋板的连接处，并且吸力基础顶部也有明显的承力特性，针对这两个关键结构的最大受力部位进行分析研究。桩体和吸力基础最大应力部位的应力时程曲线如图6 所示。通过基于传统理论的线性叠加法在风浪流荷载联合作用下风机结构的动力响应分析可知，极端工况下环境荷载对风机结构影响极大，结构的150-一极端工况一运行工况50100时间s（b）吸力基础顶部位移时程曲线200150200122023年10 月海洋工程陈实等：海上吸力基础式风机结构动力响应研究运行工况极端工况10080604020120110100908070605040图6 风机结构关键部位最大应力时程曲线动力特性也更具有研究代表性，因此在后续采用Turkstra准则进行荷载组合时，选取最有代表性的响应情况更剧烈的极端工况进行进一步的分析研究。3.2Turkstra准则本文选取的第二种多荷载作用组合方式为国际上普遍采用的Turkstra组合方法，Turkstra建议轮流地以一个荷载效应在O，T 的极值与其余荷载效应的瞬时值组合。故本节所计算的荷载组合为表5的组合工况。表5组合工况组合风工况组合1极值组合2运行组合3运行组合4极值组合5极值组合6运行分别通过上述6 种常见的不同极端荷载组合方式，对海上风机结构进行动力响应分析，同上文所述，本节依旧以塔筒顶部位移、吸力桶顶部位移时程曲线以及桩体和吸力基础最大应力部位的应力时程曲线作为主要分析因素，通过有限元模拟，计算所得结果如图7 所示。各组合下风机结构动力响应分析结果表明，以组合5对风机结构在分析时间内的运动幅度影响最显著。因此选取组合5作为Turkstra准则组合下的50100时间/s(a）桩体50100时间/s(b）吸力桶浪工况运行极值运行极值运行极值150一运行工况极端工况150200流工况运行运行极值运行极端极值200代表组合，将其与上节中通过线性叠加法计算的极端工况结果进行对比。极端工况下的线性叠加法与Turkstra准则的对比分析如图8 所示，可以观察到，采用线性叠加法所获得的结构响应特性与组合5作用下的计算结果相似度较高，从结构的位移响应和应力情况可以判断，组合5的荷载作用对风机结构动力响应的影响更大，并且Turkstra准则结果更加保守，因此通过上述分析，以外部激励对风机结构的位移和应力响应的极值、幅度的影响，以及以结构所表现出的动力特性作为依据，可以确定组合5（风、流荷载为极限工况）为具有研究代表性的最不利荷载组合方式。在风机设计和选址时，应尽可能避免这些荷载的组合情况。1.4组合1组合2 一组合3组合4 组合5组合61.21.00.60.40.20.000.0300.0250.0200.015位0.0100.0050.000-0.005M50100时间/s（a）塔筒顶部位移时程曲线一组合1组合2 组合3一组合4 组合5一组合650100时间s（b）吸力基础顶部位移时程曲线150150200200第4 9 卷第5期13海洋工程石由工程建设250T20015010050180160140120100806040200图7 基基于Turkstra准则的风浪流荷载耦合作用下风机动力响应情况0.9一线性叠加-Turkstra组合50.80.70.60.50.40.30.20.100.020m0.0150.0100.0050.000-0.0050（b）吸力桶顶部位移时程曲线140一组合1一组合2 一组合3一组合4 组合5一组合650100时间s（c）桩体最大应力时程曲线一组合1组合2 一组合3组合4 组合5一组合650100时间/s（d）吸力桶基础最大应力时程曲线50100时间/s（a）塔筒顶部位移时程曲线50100时间s线性叠加-Turkstra组合512010080604020150200150200150200线性叠加-Turkstra组合51502000120110100/9080706050400图：线性叠加法与Turkstra准则对比分析4结论本文以5MW吸力基础固定式海上风机为研究对象，构建了多荷载作用下的海上风机结构耦合分析方法，采用不同荷载叠加方法，通过有限元软件开展了风浪流联合作用下海上风机结构动力响应分析，获得以下主要结论。1）海上风机受到空气动力、水动力、基础与土非线性作用力的耦合作用，应采用建立海上风机基础和塔筒耦合模型方法开展动力响应分析。从本文选择的吸力基础海上风机结构的动力响应分析结果可以看出，海上风机整体结构位移最大区域在塔筒顶部，应力最大区域在吸力基础与风机塔架的连接处。2）通过对线性叠加法和Turkstra准则的对比分析验证，可以看出，采用Turkstra准则能更准确地获得对风机结构动力响应最不利的荷载组合方式，可以为实际工程中海上风机结构设计提供可靠依据。50（）桩体最大应力时程曲线一线性叠加.-.Turkstra组合550100时间/s（d）吸力桶最大应力时程曲线100时间/s150150200200142023年10 月海洋工程陈实等：海上吸力基础式风机结构动力响应研究参考文献1 SCHWARTZ S,ARGYRIADIS K.Analysis of the fatigue load-ing of an offshore wind turbine using time and frequency do-main methodsC)/Proceedings of the EWEC:EWEC,2001.2王明超，赵永生，杜炜康，等.海上风力机桩-土相互作用模型分析J.中国海洋平台，2 0 14，2 9（6)：4 1-4 7.3张小玲，刘建秀，杜修力，等.风浪流共同作用下海上风电基础与海床的动力响应分析J.防灾减灾工程学报，2 0 18，38(4):658-668.4沈玉光.海上风电筒型基础风机结构体系动力响应分析D.天津：天津大学，2 0 12.5刘建秀.风浪流共同作用下海上风电基础动力响应与承载性能分析D.北京：北京工业大学，2 0 18.6金伟良.工程荷载组合理论与应用M.北京：机械工业出版社，2006.7 VAN DER VALK P L C,VOORMEEREN S N,DE VALK P C,et al.Dynamic models for load calculation procedures of off-shore wind turbine support structures:Overview,assess-ment,and outlookJ.Journal of Computational and Nonlin-earDynamics,2015,10(4):041013.8张伟，冯景，姚明辉，等.水平轴风机旋转叶片非线性动力学模型的建立及分析J.科技导报，2 0 10，2 8（5）：4 0-4 4.9 ZHU J,ZHANG W.Numerical simulation of wind and wavefields for coastal slender bridgesJJ.Journal of Bridge Engi-neering,2017,22(3):04016125.10 JONKMAN J,BUTTERFIELD S,MUSIAL W,et al.Defi-nition of a 5-MW reference wind turbine for offshore sys-（上接第7 页）3 SNAME.Example(GO-BY)calculations using recom-mended practice for site specific assessment of mobile jack-up unitsS.Society of Naval Architects and Marine Engineers,Technical and Research Bulletin5-5A,New Jersey:2002.4】王建华，兰斐.钻井船插桩对邻近桩影响的耦合欧拉-拉格朗日有限元方法研究J.岩土力学，2 0 16（4）：112 7-1136.5 THO K K,LEUNG C F,CHOW Y K,et al.Eulerian finite el-ement simulation of spudcan-pile interactionJ.CanadianGeotechnical Journal,2013(6):595-608.6 LI Y,ZHOU Y,TAO J,et al.Numerical study of spudcanpenetration induced loading on an adjacent jacket Pile in softclaysC/The Twenty-fifth international offshore and polarengineering conference.Beijing:International society ofoffshore and polar engineers,2015:66-68.7】刘占阁，王栋，栾茂田.钻井船基础压入对邻近平台桩基影响的有限元分析CI/中国土木工程学会第九届土力学及岩土工程学术会议论文集（下册），北京：中国土木工程学会，2 0 0 3：1 439-1 442.8吴永韧，鲁晓兵.桩靴贯入对固定平台基础扰动的数值模拟C/第十四届中国海洋（岸）工程学术讨论会论文集.呼和浩特：中国海洋工程协会，2 0 0 9：118-12 2.9 MEYERHOF G G.Ultimate bearing capacity of footings onsand layer overlying clayJJ.Canadian Geotechnical Journal,tem developmentR.National Renewable Energy Lab.(NREL),Golden,CO(United States),2009.11NB/T101052018，海上风电场工程风电机组基础设计规范S.12孔德森，刘一，邓美旭，等.非均质土中海上风电单桩基础动力响应特性J.工程科学学报，2 0 2 1，4 3（5：7 10-7 19.13 SEIDEL M,KELMA S.Stochastic modelling of wind andWave induced loads on jacket pilesJ.Stahlbau,2012,81(9):705-710.14 LOPEZ-VILLALOBOS C A,RODRIGUEZ-HERNANDEZ O,MARTINEZ-ALVARADO O,et al.Effects of wind poWerspectrum analysis over resource assessmentJ.Renew-able Energy,2021,167:761-773.15 LUKASSEN L J,STEVENS R J A M,MENEVEAU C,et al.Modeling space-time correlations of velocity fluctuations inwind farmsJ).Wind energy,2018,21(7):474-487.基金项目：中海油能源发展股份有限公司重大科技专项（HFZXKT-JN2021-01)。作者简介：陈实（19 8 2 一），男，福建宁德人，高级工程师，2 0 0 4 年毕业于华南理工大学船舶与海洋工程专业，现在主要从事船舶与海洋工程结构物设计与研发工作。Email：c h e n s h i 1 c n o o c.c o m.c n收稿日期：2 0 2 3-0 7-131974,11(2):223-229.10 MEYERHOF G G,HANNA A M.Ultimate bearing capacityof foundations on layered soils under inclinJ.CanadianGeotechnical Journal,1978,15（4):56 5-57 2.11 DESAI C S,REESE L C.Analysis of circular footings on lay-ered soilsJ.Journal of Soil Mechanics&Foundations Div,1970(2):125-129.12王琦.自升式平台插桩物理模型实验及分析J.中国水运月刊,2 0 14，14(3)：36 4-36 6.13】南京水利科学研究院.中华人民共和国行业标准：土工试验规程SL237-1999M.北京：中国水利水电出版社，19 9 9.基金项目：中海石油（中国）有限公司科研项目“南海西部油田上产2 0 0 0 万方钻完井关键技术研究”（CNOOC-KJ135ZDXM38ZJ05ZJ）。作者简介：刘贤玉（19 8 7 一），男，广东湛江人，高级工程师，2 0 12 年毕业于中国石油大学（华东）工程力学专业，硕士，主要从事海上钻完井工艺相关研究工作。Email：l i u x y 19 c n o o c.c o m.c n收稿日期：2 0 2 3-0 7-31第4 9 卷第5期15