-
近十年来,随着全球能源向低碳能源转型,海上风电产业正蓬勃发展并形成规模化势头,中国在近年来也加速发展海上风电产业。截至2022年底,全球海上风电累计装机64.3 GW,中国累计海上风电装机占比49%,占据全球第一大的装机市场,并且2022年新增海上风电装机容量排名第一[1]。中国东南沿海年平均风速分布如图1所示。中国南海海域辽阔,深远海海上风能资源丰富,年平均风速可达9 m/s以上,具有开发深远海风电的巨大潜力。其中,广东省35~50 m水深的近海深水区海上风电规划容量在当前规划总容量的占比超过85%,而超过50 m水深海域的海上风电可开发容量更为巨大,深远海风电是广东省未来能源发展的重要方向[2]。
图 1 中国东南沿海年平均风速分布示意图
Figure 1. Annual average wind speed distribution map in the southeast coast of China
海上风机基础形式如图2所示。受技术和工程经济性限制,海上风机支撑结构形式也必然伴随着水深而发生变化,从固定式的支撑结构演变为漂浮式的支撑结构[3]。海上漂浮式风机极大地拓展了海上风电的应用范围,并且具有诸多的优势[4],例如:机位部署更加灵活、可在岸上完成整体组装、海上施工安装更加方便、对海床生态环境友好、可完全拆解与迁移、更经济地搭载更大功率的风电机组等。
相比于固定式风机基础,漂浮式风机基础受力状态更加复杂,浮式基础作为顺应式结构,对波浪荷载敏感,运动幅度较大,气动荷载-结构-水动力-系泊-控制系统的耦合作用更加明显[5],尤其对于大容量的风电机组,风机荷载导致的气动-运动耦合作用更加明显,解耦分离迭代的设计方法会给浮式风电设计造成较大的计算误差,尤其在目前缺乏规模化的浮式风电场建设经验情况下,荷载抗力系数和安全系数取值难以准确把握,容易造成设计过于保守或偏激。因此,需要采用一体化计算软件分析海上漂浮式风机系统耦合动力响应特性。
当前的海上漂浮式风机一体化时域计算软件可以大致划分为如下几类:一类从传统的海洋工程水动力软件发展而来,如Sima[6]、OrcaFlex[7]、Deeplines[8]、3DFLOAT[9]等;一类从传统的固定式海上风机计算分析软件发展而来,如FAST[10-11]、GH Bladed[12]、HAWC2[13]等;一类则来自于航空航天、车辆工程等多体机械软件仿真软件,如SIMPACK[14]、MSC- ADAMS[15]等。还有一些基于多个软件进行联合数值模拟,如AQWA-FAST[16]。近些年,部分学者针对浮式风机动力特点进行针对性的一体化软件开发,如DARwind[17-18]等。还有一些学者通过CFD技术对海上浮式风机进行数值计算[19-20],部分统计如表1所示。
表 1 海上漂浮式风机时域耦合计算程序统计表
Table 1. Statistical table of time-domain coupling calculation program for floating offshore wind turbines
名称 空气动力学 水动力学 结构动力学 系泊模型 FAST[10-11] (BEM/GDW)+DS+FVM PF+QD+(QTF)+MR+NA T: Mod/MB
P: RigidQS+FE+LM Bladed[12] (BEM/GDW)+DS MR+(IP+IWL) T: Mod/MB
P: RigidQS HAWC2[13] (BEM/GDW)+DS MR T: MB/FE
P: MB/FEFE 3Dfloat[9] BEM+FDT MR+(IWL) T: FE
P: FEFE Sima[6] (BEM/GDW)+DS PF+MR T: FE
P: FEFE DARwind[17-18] BEM PF+MR+OTF+IP+IWL+QD T: Mod/MB
P: RigidOS+LM 注:BEM-叶素动量理论;GDW-广义动态尾迹;DS-动态失速;FVM-自由涡;FDT-过滤后的动态推力;PF-势流理论;QD-二次拖曳力;OTF-二次传递函数;MR-Morison方程;NA-Newman近似;IP-瞬时位置;IWL-瞬时水位;T-风机;P-平台;Mod-模态;MB-多体动力学;FE-有限元;Q-准静态;LM-集中质量。 文章结合广东近海深水区某海上风电场,针对10 MW级的混凝土半潜式浮式风机开展一体化仿真计算,并且研究其在发电工况和极端工况下的耦合动力响应特性,以期揭示近海深水区大容量半潜型浮式风机的动力响应特性和关键设计因素。
-
OO-Star[21]风力机浮式基础由挪威Olav Olsen公司研发,属于半潜型平台,其基础底部是一个“星”型浮筒结构,上面布置一个中央立柱,外围设置3个立柱结构,由预应力钢筋混凝土材料构成,具体结构形式如图3所示。文章采用DTU 10 MW的风力发电机组,RNA总重为676 t,额定功率为10 MW,风机塔柱与OO-Star基础中央立柱连接,关于该风电机组和塔筒的详细数据可参考文献[22]。上述的漂浮式风机整体吃水为22 m,排水量2.350 9E+04 m3,典型参数见表2。
名称 数值 叶轮质量/t 230.72 机舱质量/t 446 塔柱质量/t 1 257 平台法兰顶距水面/m 11.0 基础质量(不含塔柱,系泊系统)/t 2.17E+04 重心高(相对水面)/m −15.22 平台横摇惯性矩(含压载)/(kg·m2) 9.43E+09 平台纵摇惯性矩(含压载)/(kg·m2) 9.43E+09 平台首摇惯性矩(含压载)/(kg·m2) 1.63E+010 设计吃水(含系泊)/m 22.0 总排水量/m3 2.350 9E+04 针对广东40~50 m的近海深水区,基于“刚度相似”准则,根据原始130 m水深母型系泊[21]设计出了适用于该近海深水区的浮式风力机系泊系统。该浅水系泊系统设有3根锚链,每根锚链之间的夹角为120°(如图4所示)。锚链直径为137 mm,长度为652 m。沿锚链长度方向距离导缆孔60 m处设置一个50 t配重块,然后在距离导缆孔134 m处,即锚链与海底的切点位置再增加一个10 t的配重,锚链破断张力为1.6992E+07 N,如图5所示。原型130 m水深的锚链系泊刚度和针对40~50 m的近海深水区改造后的锚链系泊刚度在纵荡方向的对比如图6所示,表明水深变浅后,系泊刚度的非线性特性大大增强。在平台大范围偏移下,将产生较大的系泊张力,容易引起系泊破断风险[23]。
-
为了更加准确地考虑漂浮式风机不同荷载和基础运动之间的耦合关系,文章采用一体化的海上漂浮式风机数值模型。一体化计算程序需要将风轮气动荷载、结构动态响应、基础平台受到的波浪和海流的水动力荷载、锚泊系统约束张力和控制系统作用力统一到一个数值模型中,并且尽可能地在同一个动力学方程中同时求解各个自由度的信息,从而尽可能地获得准确的计算结果。
本计算采用FAST(Fatigue, Aerodynamics, Structures and Turbulence)[10]软件进行上述的漂浮式风机一体化数值计算。该软件由美国可再生能源实验室(NREL)开发,是一款包含空气动力学、水动力学、柔性多体动力学以及机电控制系统为一体的,用于2~3桨叶水平轴风力发电机的载荷计算工具。软件综合了气动力模块、结构动力分析模块、水动力分析模块和控制器模块,可以执行整个浮式风机的气动-水动力-伺服-弹性耦合的时域耦合数值分析,软件各个功能模块的计算逻辑如图7所示。该软件气动荷载的计算主要采用BEM方法或GDW方法,水动力荷载采用势流方法求解大型结构体的水动力,并辅助以Morison方程修正细长体水动力荷载,结构动力学采用kane动力学方程,系泊模型包含悬链线、集中质量法和有限元模型等。
-
结合广东省近海水深区某海上风电场水文气象数据,该海域平均海平面为44 m,本研究数值模拟选取2种计算工况,分别为额定作业工况(轮毂处10 min平均风速为11.4 m/s)以及生存工况(50年一遇,轮毂处10 min平均风速为56.2 m/s)。每种工况对应的环境参数如表3所示。针对每种工况分别选取5个随机种子,模拟随机环境载荷,每个随机种子下模拟1 h,针对每种工况分别选取了3种环境载荷入射方向分别为0°、90°和45°。波浪采用随机波由三参数Jonswap谱定义,海流取为均匀流。
表 3 计算工况及环境参数
Table 3. Load cases and environmental parameters
名称 编号 方向/(°) 风 浪 流 轮毂处风速/(m·s−1) 风切变指数 有义波高/m 谱峰周期/s 谱峰因子 流速/(m·s−1) 发电工况 DLC 1a 0 11.4 0.067 1.75 7.5 3.3 0.46 DLC 1b 90 11.4 0.067 1.75 7.5 3.3 0.46 DLC 1c 45 11.4 0.067 1.75 7.5 3.3 0.46 生存工况 DLC 3a 0 56.2 0.096 11.53 14.93 3.3 1.19 DLC 3b 90 56.2 0.096 11.53 14.93 3.3 1.19 DLC 3c 45 56.2 0.096 11.53 14.93 3.3 1.19 文章根据上述10 MW半潜型浮式风力机的相关参数建立浮式风力机数值模型。利用SESAM软件[24]建立浮式基础的有限元模型,并求解浮式基础水动力系数,计算结果作为FAST软件水动力模块的输入文件,其气动计算、系泊计算、结构动力计算和风机控制算法均在FAST软件实现,按照设定的环境工况,进行数值模拟,最后对计算结果进行统计分析,计算流程如图8所示。
Integrated Calculation and Coupled Dynamic Analysis of a Concrete Semi-Submersible Floating Wind Turbine in Offshore Deep Water Area
-
摘要:
目的 文章旨在探索10 MW级半潜型浮式风机的一体化计算方法并分析其在40~50 m近海深水区的耦合动力响应特性。 方法 以10 MW混凝土半潜型浮式风机为例,构建一体化时域数值计算模型,统计分析其额定发电和极端工况下的平台动力响应和系泊张力特征。 结果 平台水平运动主要受波浪荷载、风荷载和系泊刚度特性的影响,最大水平运动和系泊张力发生在生存工况,垂向运动主要受波浪荷载影响,摇摆运动的均值主要受风荷载影响,上述浮式风机动力响应均满足设计指标。 结论 一体化数值计算方法较好地考虑浮式风机耦合动力特性,由于水深限制,近海深水区的海上漂浮式风机对水平运动约束和系泊非线性问题的优化更为重要,响应极值主要发生在极端工况,上述结论为此类漂浮式海上风机基础结构的研究与设计工作提供了一定的参考。 Abstract:Introduction The research aims to explore the integrated calculation method of a 10 MW semi-submersible floating wind turbine, and analyze its coupled dynamic response characteristics in 40~50 m offshore deep water areas. Method A 10 MW concrete semi-submersible floating wind turbine was taken as an example, and then numerical calculation was carried out by the integrated calculation method, and its coupled dynamic response under rated and survival conditions was statistically analyzed. Result The horizontal motion of the platform is mainly affected by the wave force, wind loading and mooring stiffness. The maximum value of motion and mooring tension occur in the survival condition, and the heave motion is mainly affected by the wave, but the mean value of the pitch/roll motion is mainly affected by the wind loading, all of which meet the design specification. Conclusion The integrated calculation method better considers the coupled dynamic behavior of floating wind turbines. Due to the limitation of water depth, the optimization of horizontal motion and mooring nonlinearity of offshore floating wind turbines is more important in offshore deep water areas, and the extreme response mainly occurs in survival conditions. The above conclusions provide an important reference for the research and design of the floating offshore wind turbines in offshore deep water area. -
表 1 海上漂浮式风机时域耦合计算程序统计表
Tab. 1. Statistical table of time-domain coupling calculation program for floating offshore wind turbines
名称 空气动力学 水动力学 结构动力学 系泊模型 FAST[10-11] (BEM/GDW)+DS+FVM PF+QD+(QTF)+MR+NA T: Mod/MB
P: RigidQS+FE+LM Bladed[12] (BEM/GDW)+DS MR+(IP+IWL) T: Mod/MB
P: RigidQS HAWC2[13] (BEM/GDW)+DS MR T: MB/FE
P: MB/FEFE 3Dfloat[9] BEM+FDT MR+(IWL) T: FE
P: FEFE Sima[6] (BEM/GDW)+DS PF+MR T: FE
P: FEFE DARwind[17-18] BEM PF+MR+OTF+IP+IWL+QD T: Mod/MB
P: RigidOS+LM 注:BEM-叶素动量理论;GDW-广义动态尾迹;DS-动态失速;FVM-自由涡;FDT-过滤后的动态推力;PF-势流理论;QD-二次拖曳力;OTF-二次传递函数;MR-Morison方程;NA-Newman近似;IP-瞬时位置;IWL-瞬时水位;T-风机;P-平台;Mod-模态;MB-多体动力学;FE-有限元;Q-准静态;LM-集中质量。 名称 数值 叶轮质量/t 230.72 机舱质量/t 446 塔柱质量/t 1 257 平台法兰顶距水面/m 11.0 基础质量(不含塔柱,系泊系统)/t 2.17E+04 重心高(相对水面)/m −15.22 平台横摇惯性矩(含压载)/(kg·m2) 9.43E+09 平台纵摇惯性矩(含压载)/(kg·m2) 9.43E+09 平台首摇惯性矩(含压载)/(kg·m2) 1.63E+010 设计吃水(含系泊)/m 22.0 总排水量/m3 2.350 9E+04 表 3 计算工况及环境参数
Tab. 3. Load cases and environmental parameters
名称 编号 方向/(°) 风 浪 流 轮毂处风速/(m·s−1) 风切变指数 有义波高/m 谱峰周期/s 谱峰因子 流速/(m·s−1) 发电工况 DLC 1a 0 11.4 0.067 1.75 7.5 3.3 0.46 DLC 1b 90 11.4 0.067 1.75 7.5 3.3 0.46 DLC 1c 45 11.4 0.067 1.75 7.5 3.3 0.46 生存工况 DLC 3a 0 56.2 0.096 11.53 14.93 3.3 1.19 DLC 3b 90 56.2 0.096 11.53 14.93 3.3 1.19 DLC 3c 45 56.2 0.096 11.53 14.93 3.3 1.19 -
[1] Global Wind Energy Council. Global wind report 2023 [R]. Brussels, Belgium: Global Wind Energy Council, 2023. [2] 王诗超, 刘嘉畅, 刘展志, 等. 海上风电产业现状及未来发展分析 [J]. 南方能源建设, 2023, 10(4): 103-112. DOI: 10.16516/j.gedi.issn2095-8676.2023.04.010. WANG S C, LIU J C, LIU Z Z, et al. Analysis of current situation and future development of offshore wind power industry [J]. Southern energy construction, 2023, 10(4): 103-112. DOI: 10.16516/j.gedi.issn2095-8676.2023.04.010. [3] CRUZ J, ATCHESON M. Floating offshore wind energy: the next generation of wind energy [M]. Cham: Springer, 2016. DOI: 10.1007/978-3-319-29398-1. [4] 陈嘉豪, 裴爱国, 马兆荣, 等. 海上漂浮式风机关键技术研究进展 [J]. 南方能源建设, 2020, 7(1): 8-20. DOI: 10.16516/j.gedi.issn2095-8676.2020.01.002. CHEN J H, PEI A G, MA Z R, et al. A review of the key technologies for floating offshore wind turbines [J]. Southern energy construction, 2020, 7(1): 8-20. DOI: 10.16516/j.gedi.issn2095-8676.2020.01.002. [5] CORDLE A, JONKMAN J. State of the art in floating wind turbine design tools [R]. Golden: National Renewable Energy Laboratory, 2011. [6] NIELSEN F G, HANSON T D, SKAARE B. Integrated dynamic analysis of floating offshore wind turbines [C]//25th International Conference on Offshore Mechanics and Arctic Engineering, Hamburg, Germany, June 4-9, 2006. Hamburg, Germany: ASME, 2006: 671-679. DOI: 10.1115/OMAE2006-92291. [7] THOMSEN J B, BERGUA R, JONKMAN J, et al. Modeling the TetraSpar floating offshore wind turbine foundation as a flexible structure in OrcaFlex and OpenFAST [J]. Energies, 2021, 14(23): 7866. DOI: 10.3390/en14237866. [8] BLONDEL F, BOISARD R, MILEKOVIC M, et al. Validation and comparison of aerodynamic modelling approaches for wind turbines [J]. Journal of physics:conference series, 2016, 753(2): 022029. DOI: 10.1088/1742-6596/753/2/022029. [9] MYHR A, HANSEN F E, MOELLER U, et al. A comparison of existing and conceptual designs for floating wind turbines [C]//Renewable Energy Research Conference 2010, July 6-August 6, 2010. [10] JONKMAN J M. Dynamics of offshore floating wind turbines—model development and verification [J]. Wind energy, 2009, 12(5): 459-492. DOI: 10.1002/we.347. [11] COULLING A J, GOUPEE A J, ROBERTSON A N, et al. Validation of a FAST semi-submersible floating wind turbine numerical model with DeepCwind test data [J]. Journal of renewable and sustainable energy, 2013, 5(2): 023116. DOI: 10.1063/1.4796197. [12] DNV. Wind turbine design software - Bladed [EB/OL]. [2023-06-05]. https://www.dnv.com/services/wind-turbine-design-software-bladed-3775. [13] LARSEN T J, HANSEN A M. How 2 HAWC2, the user's manual [R]. Roskilde, Denmark: Riso National Laboratory, 2007. [14] MATHA D, SCHLIPF M, CORDLE A, et al. Challenges in simulation of aerodynamics, hydrodynamics, and mooring-line dynamics of floating offshore wind turbines [R]. Golden: National Renewable Energy Laboratory, 2011. [15] WITHEE J W. Fully coupled dynamic analysis of a floating wind turbine system [D]. Cambridge: Massachusetts Institute of Technology, 2004. [16] YAO Y S, NING D Z, DENG S J, et al. Hydrodynamic investigation on floating offshore wind turbine platform integrated with porous shell [J]. Energies, 2023, 16(11): 4376. DOI: 10.3390/en16114376. [17] CHEN J H, HU Z Q, LIU G L, et al. Coupled aero-hydro-servo-elastic methods for floating wind turbines [J]. Renewable energy, 2019, 130: 139-153. DOI: 10.1016/j.renene.2018.06.060. [18] 陈嘉豪, 刘格梁, 胡志强. 海上浮式风机时域耦合程序原理及其验证 [J]. 上海交通大学学报, 2019, 53(12): 1440-1449. DOI: 10.16183/j.cnki.jsjtu.2019.12.006. CHEN J H, LIU G L, HU Z Q. Development and validation of a time-domain coupling simulation code for floating offshore wind turbines [J]. Journal of Shanghai jiaotong university, 2019, 53(12): 1440-1449. DOI: 10.16183/j.cnki.jsjtu.2019.12.006. [19] LI Y W, PAIK K J, XING T, et al. Dynamic overset CFD simulations of wind turbine aerodynamics [J]. Renewable energy, 2012, 37(1): 285-298. DOI: 10.1016/j.renene.2011.06.029. [20] TRAN T, KIM D, SONG J. Computational fluid dynamic analysis of a floating offshore wind turbine experiencing platform pitching motion [J]. Energies, 2014, 7(8): 5011-5026. DOI: 10.3390/en7085011. [21] PEGALAJAR-JURADO A, BREDMOSE H, BORG M, et al. State-of-the-art model for the LIFES50+ OO-Star wind floater semi 10MW floating wind turbine [J]. Journal of physics: conference series, 2018, 1104: 012024. DOI: 10.1088/1742-6596/1104/1/012024. [22] BAK C, ZAHLE F, BITSCHE R, et al. The DTU 10-MW reference wind turbine [R]. Fredericia, Denmark: DTU, 2013. [23] 李旭, 肖龙飞, 魏汉迪, 等. 系泊缆破断对半潜式平台水动力性能影响的试验研究 [J]. 海洋工程装备与技术, 2023, 10(1): 48-56. LI X, XIAO L F, WEI H D, et al. Experimental study on the influence of mooring line breaking on the hydrodynamic performance of semi-submersible platform [J]. Ocean engineering equipment and technology, 2023, 10(1): 48-56. [24] 张会良, 肖龙飞, 徐秀龙. 深水圆筒型钻井平台张紧式系泊系统设计 [J]. 海洋工程, 2021, 39(6): 10-18. DOI: 10.16483/j.issn.1005-9865.2021.06.002. ZHANG H L, XIAO L F, XU X L. Design of taut mooring system for a deep water cylindrical drilling platform [J]. The ocean engineering, 2021, 39(6): 10-18. DOI: 10.16483/j.issn.1005-9865.2021.06.002. [25] DNV. Standard for classification of wind turbine installation units: DNV—OS—J301 [S]. Oslo, Norway: DNV, 2011.