-
全潜式浮式风机是丁红岩等综合3种传统的结构型式单立柱(Spar-buoy)式、张力腿式(Tension-leg platform)和半潜式(Semi-submersible platform)基础的优点提出一种新型浮式风机[12],适合超过60 m的深海区域。全潜式浮式风机的拖航系统由基础、风机、拖船和缆绳组成。全潜式基础由立柱、水平浮箱、浮筒、斜撑和水平浮箱组成,基础结构如图1所示,尺寸如表1和表2所示。
表 1 全潜式浮式风机参数
Table 1. Parameters of submersible floating OWT
参数 值 参数 值 浮筒直径/m 11 水平浮箱高度/m 4 浮筒高度/m 15 水平浮箱宽度/m 6.5 浮筒中轴线间距/m 50 水平撑杆高度/m 3 立柱高度/m 22 排水体积/m3 12 811 立柱直径/m 8.3 基础质量/t 4 812 撑杆直径/m 2 表 2 DTU 10 MW风机参数
Table 2. Parameters of DTU 10 MW OWT
参数 值 叶轮直径/m 178.3 叶片长度/m 86.366 轮毂高度/m 119 叶轮质量/t 227.962 机舱质量/t 446.036 塔筒质量/t 628.442 -
本文采用Moses(Mulyi-operational Structural Engineering Simulator)软件[5]进行数值模拟计算,是一种将建模与结构应力分析整合到一起的编程软件。其主要功能包括海上安装模拟、水动力分析、稳性校核、结构分析校核等多个方面。计算时Moses可以实现模拟环境荷载并进行应力分析。
在Moses中建立的模型拖船——全潜式浮式风机拖航系统如图2所示,由拖船-龙须缆-风机组成,其中风机的主要结构可划分为全潜式浮式基础、塔筒和机头。风机和拖船通过Y形的龙须缆相连,龙须缆一端系在拖船尾部,另外两端分别与全潜式浮式风机前端立柱相连。风机整体吃水为5.15 m,基础底座提供的浮力可以承担起整个风机结构的重量。以塔筒底面的几何形心为原点,x轴为垂直正对风机叶片方向,z轴为垂直海面向上,y轴按右手定则确定,则上部风机的重心为(−0.3,0.0,85.5),考虑全潜式基础后,以基础底面的几何形心为原点,重心位置为(−0.3,0,32)。
本文主要研究海洋中风荷载对全潜式浮式风机拖航运动稳定性的影响。采用NPD风谱生成时程风速序列进行模拟,张延涛等[13]对渤海风速进行采集与分析,得到某天渤海10 min内平均风速分布在0~18.5 m/s。为研究风速对全潜式浮式风机拖航运动响应的影响,得到不同风速下较为明显的运动响应差别,同时考虑到特殊天气的影响,这里选取渤海某天10 min内的,水面以上0~15.3 m高度范围内的平均风速,即10.3 m/s、20.6 m/s和30.9 m/s 3种不同高度下的风速进行模拟。选取5组不同工况开展全潜式浮式风机的拖航仿真,如表3所示。全潜式浮式风机的基础为漂浮式基础,而漂浮式风机多适用于水深超过50 m的深海区域[14],本文为研究风速和风向对全潜式浮式风机在深海拖航时的影响,在数值模拟过程中采用控制变量的方法,研究不同风速下全潜式浮式风机的运动响应,以实际情况为基础设置环境条件,拖航水深为100 m。由于实际工程中全部拖航过程可持续24 h以上,数值模拟全过程成本较大,同时一定环境条件下全潜式浮式风机的拖航运动响应时程曲线具有周期性变化规律,因此选择拖航从静止到开始运动的初始阶段中,稳定的1 000 s观察风速和风向对拖航影响程度。以全潜式浮式风机基础底面中点为原点,拖缆点高度同吃水高度为5.2 m。其中工况H1、H2、H3为一组研究不同风速对全潜式浮式风机拖航运动稳性的影响,工况H2、H4、H5为一组研究不同风向对风机拖航运动稳性的影响。
表 3 拖航工况设置
Table 3. Towing condition settings
工况 风速/(m·s−2) 有义波高/m 波浪周期/s 风向/(°) H1 10.3 1 8 180 H2 20.6 1 8 180 H3 30.9 1 8 180 H4 20.6 1 8 90 H5 20.6 1 8 0 -
风是一种不确定的自然环境条件,风压[15]和风力可用规定时段内的平均风速来计算。风荷载[16]的基本表达式为:
风压:
$$ q = \dfrac{1}{2}{\rho _{\mathrm{a}}}U_{T,Z}^2 $$ (1) 式中:
q ——基本风压(Pa);
ρa ——空气密度(kg/m3);
UT,Z ——T平均时间内的高度Z处对应的平均风速(m/s)。
和一般建筑结构荷载计算不同,此处基本风压指的是T平均时间内高度Z处的风压大小。
风力:
$$ {F_{\mathrm{W}}} = q\displaystyle \sum\limits_1^n {{C_{{Z}}}} {C_{\mathrm{S}}}{A_n} $$ (2) 式中:
CZ——受风结构高度系数;
CS——构件形状系数;
An ——受风部件的迎风面积(m2)。
由于全潜式浮式风机是高耸结构,风速随着高度的变化而变化,为准确计算风力,不同的高度需要采用不同的高度系数,规范API RP 2SK[17]给出的基于1 min平均风速的风力高度系数转换关系如表4所示,Moses中可以根据需要设置不同规范中的风荷载定义方式进行计算,本文选择国际上广泛认可的API规范作为依据进行计算。
表 4 风力高度系数CZ
Table 4. Height coefficient CZ
水面以上面积中心高度h/m CZ 0<h≤15.3 1.00 15.3<h≤30.5 1.18 30.5<h≤46.0 1.31 46.0<h≤61.0 1.40 61.0<h≤76.0 1.47 由阵风引起的脉动风荷载也同样需要引起注意,当风场的主要能量周期和拖航系统的震荡周期相近时,脉动风可能对拖航系统产生较大的震荡[17],甚至有可能会引起风和结构物的共振。由此引入了风谱的概念。根据API RP 2A[18]推荐,本文使用NPD谱进行计算:
海平面以上z米处的1 h平均风速Uz为:
$$ {{{U}}_{{z}}}= {{{U}}_{10}}\left[ {1 + {C{\mathrm{ln}}} \left( {\dfrac{{{z}}}{{10}}} \right)} \right] $$ (3) $$ {{C}} = 0.057\;3\sqrt {1 + 0.15{{{U}}_{10}}} $$ (4) 式中:
Uz ——海平面以上z米的1 h平均风速(m/s);
U10 ——海平面以上10 m处的1 h平均风速(m/s)。
NPD谱表述了风中某一质点纵向风速能量密度波动:
$$ {S_{{\rm{NP}}{{\rm{D}}_f}}}= \dfrac{{320{{\left( {\dfrac{{{{{U}}_{0}}}}{{10}}} \right)}^2}{{\left( {\dfrac{{{z}}}{{10}}} \right)}^{0.45}}}}{{{{\left( {1 + {{{f}}^{0.468}}} \right)}^{3.561}}}} $$ (5) 式中:
${S_{{\rm{NP}}{{\rm{D}}_f}}} $ ——频率f的能量谱密度(m2/s);
f ——频率(Hz)。
-
浮体在海洋中受力情况复杂,受到惯性力、附加力、回复力和环境荷载等的共同影响,国内外一般基于势流理论解决此类问题[19]。时域势流理论可以通过建立时域控制方程的方法计算浮体在任意波浪击力下的运动响应。目前常采用间接时域法,引入脉冲频谱来建立浮体时域运动方程[20]:
$$ (M + m)\ddot {X_t} + \int\limits_{ - \infty }^t {K_{t - \tau }\dot {X_\tau }{d_\tau }} + {K_{\mathrm{C}}}{X_t} = F $$ (6) 式中:
M ——浮体质量矩阵;
m ——附加质量矩阵;
Xt ——浮体运动位移;
KC——回复力矩阵;
Kt ——时延函数;
F ——波浪力。
Moses中根据纽马克β假设对式进行计算,式(6)可简化为:
$$ M\ddot X + B\dot X + KX = F $$ (7) 根据拉格朗日中值定理及纽马克假设,t时刻的速度可以表示为:
$$ {\dot X_{t + 1}} = {\dot X_i} + (1 - \delta ){\ddot X_t}\Delta t + \delta {\ddot X_{i + 1}}\Delta t $$ (8) 将一段时间内的时间段按照一定步长进行离散,根据拉格朗日中值定理并进行泰勒级数展开,得:
$$ {X_{t + 1}} = {X_t} + \dot X\Delta t + \left[ {\left( {\dfrac{1}{2} - \beta } \right){{\ddot X}_i} + \beta {{\ddot X}_{t + 1}}} \right]\Delta {t^2} $$ (9) 式中,$ \delta $和β是和精度、稳定性有关的参数,Moses中默认$ \delta $=0.5,保证积分过程中振幅保持不变,取β=0.25为平均加速度法。经计算,可以得到全潜式浮式风机在模拟工况中的运动响应时程曲线。
Effect of Wind Loads on Towing Response of Submersible Floating OWT
-
摘要:
目的 为了解决浮式风机拖航运输的安全问题,对全潜式浮式基础风机的拖航运动特性展开研究。 方法 针对全潜式浮式风机拖航过程,通过Moses软件建立了“全潜式浮式风机-拖缆绳-拖船”一体化拖航系统数值模型,设置了五组不同的海况对拖航过程进行计算。文章研究风速和风浪方向对全潜式浮式风机拖航运动响应的影响,并根据分析结果,提出有助于全潜式浮式风机拖航安全的建议。 结果 结果表明:相比于横摇运动,全潜式浮式风机的纵摇和垂荡运动受风速影响较为明显,风速越大,运动幅度越大;当风浪方向为90°时,风机的横摇运动受风浪方向影响较大,风浪方向为0°时,风机的纵摇和垂荡运动受风浪方向影响较大。 结论 在保证拖航过程在拖航窗口期的前提下,实际工程中建议对全潜式浮式风机采用逆风拖航的方式,以降低风机振动和共振的风险。同时,应尽量避免风向与拖航方向垂直的情况,两者垂直容易引发较大程度的横摇运动,使拖航过程较为危险。逆风拖航时,在超过平均海况的较大风速30.9 m/s下,全潜式浮式风机的运动响应及拖缆力仍能满足安全拖航的要求。 Abstract:Introduction In order to solve the safety problem of towing transportation of floating OWT( offshore wind turbine ), the towing characteristics of fully submersible floating OWT are studied. Method The model of the fully submersible floating OWT towing system was established through Moses software. For the towing process, five groups of sea conditions were set up for simulation, the numerical model calculations were carried out, the influence of wind speed and wind wave direction on towing response of fully submersible floating OWT was studied. According to the analysis results, some suggestions were put forward to help the towing safety of the fully submersible floating OWT. Result The results show that, compared with rolling, the pitch and heave of the fully submersible floating OWT are more obviously affected by the wind speed, and the greater the wind speed, the greater the motion amplitude; when the wind wave direction is 90°, the rolling of OWT is greatly influenced by the wind wave direction, and when the wind wave direction is 0°, the pitch and heave of the OWT is greatly influenced by the wind wave direction. Conclusion On the premise of ensuring that the towing process is within the towing window period, it is suggested that the fully submersible floating OWT should be towed against the wind in practical engineering to reduce the risk of vibration and resonance of OWT. At the same time, we should try our best to avoid the situation that the wind direction is perpendicular to the towing direction, which will easily lead to a large degree of rolling and make the towing process more dangerous. When towing against the wind, the motion response and towing force of the fully submersible floating OWT can still meet the requirements of safe towing at a large wind speed of 30.9 m/s exceeding the average sea condition. -
表 1 全潜式浮式风机参数
Tab. 1. Parameters of submersible floating OWT
参数 值 参数 值 浮筒直径/m 11 水平浮箱高度/m 4 浮筒高度/m 15 水平浮箱宽度/m 6.5 浮筒中轴线间距/m 50 水平撑杆高度/m 3 立柱高度/m 22 排水体积/m3 12 811 立柱直径/m 8.3 基础质量/t 4 812 撑杆直径/m 2 表 2 DTU 10 MW风机参数
Tab. 2. Parameters of DTU 10 MW OWT
参数 值 叶轮直径/m 178.3 叶片长度/m 86.366 轮毂高度/m 119 叶轮质量/t 227.962 机舱质量/t 446.036 塔筒质量/t 628.442 表 3 拖航工况设置
Tab. 3. Towing condition settings
工况 风速/(m·s−2) 有义波高/m 波浪周期/s 风向/(°) H1 10.3 1 8 180 H2 20.6 1 8 180 H3 30.9 1 8 180 H4 20.6 1 8 90 H5 20.6 1 8 0 表 4 风力高度系数CZ
Tab. 4. Height coefficient CZ
水面以上面积中心高度h/m CZ 0<h≤15.3 1.00 15.3<h≤30.5 1.18 30.5<h≤46.0 1.31 46.0<h≤61.0 1.40 61.0<h≤76.0 1.47 -
[1] 黄俊. 海上风电基础特点及中国海域的适用性分析 [J]. 风能, 2020(2): 36-40. DOI: 10.3969/j.issn.1674-9219.2020.02.013. HUANG J. Analysis of basic characteristics of offshore wind power and applicability of China waters [J]. Wind energy, 2020(2): 36-40. DOI: 10.3969/j.issn.1674-9219.2020.02.013. [2] 陈易人, 姚靳羽, 李明轩, 等. 带月池驳船式浮式风机水动力性能 [J]. 上海交通大学学报: 1-35. [2024-01-02]. DOI: 10.16183/j.cnki.jsjtu.2022.521. CHEN Y R, YAO J Y, LI M X, et al. Hydrodynamic performance of a barge-type floating offshore wind turbine with moonpool [J]. Journal of Shanghai Jiaotong University, 1-35. [2024-01-02]. DOI: 10.16183/j.cnki.jsjtu.2022.521. [3] DYMARSKI P, DYMARSKI C, CIBA E. Stability analysis of the floating offshore wind turbine support structure of cell spar type during its installation [J]. Polish maritime research, 2019, 26(4): 109-116. DOI: 10.2478/pomr-2019-0072. [4] HAN Y Q, LE C H, DING H Y, et al. Stability and dynamic response analysis of a submerged tension leg platform for offshore wind turbines [J]. Ocean engineering, 2017, 129(6): 68-82. DOI: 10.1016/j.oceaneng.2016.10.048. [5] 张浦阳, 石建超, 丁红岩, 等. 海上风电复合筒型基础结构浮运分析 [J]. 太阳能学报, 2014, 35(11): 2313-2319. DOI: 10.3969/j.issn.0254-0096.2014.11.037. ZHANG P Y, SHI J C, DING H Y, et al. Floating analysis of composite bucket foundation of offshore wind turbine [J]. Acta energiae solaris sinica, 2014, 35(11): 2313-2319. DOI: 10.3969/j.issn.0254-0096.2014.11.037. [6] 韩彦青, 丁红岩, 张浦阳, 等. 基于多体动力学的浮式风机拖航中的环境影响分析 [J]. 天津大学学报(自然科学与工程技术版), 2017, 50(10): 1055-1061. DOI: 10.11784/tdxbz20160 5059. HAN Y Q, DING H Y, ZHANG P Y, et al. Environmental effects on towing of floating wind turbine based on multibody dynamics [J]. Journal of Tianjin University (Science and Technology Edition), 2017, 50(10): 1055-1061. DOI: 10.11784/tdxbz20160 5059. [7] DING H Y, HAN Y Q, LE C H, et al. Dynamic analysis of a floating wind turbine in wet tows based on multi-body dynamics [J]. Journal of renewable and sustainable energy, 2017, 9(3): 033301. DOI: 10.1063/1.4982742. [8] 丁红岩, 韩彦青, 张浦阳, 等. 全潜式浮式风机基础在不同风况下的动力特性研究 [J]. 振动与冲击, 2017, 36(6): 201-206, 228. DOI: 10.13465/j.cnki.jvs.2017.06.031. DING H Y, HAN Y Q, ZHANG P Y, et al. Dynamic analysis of the submersible foundation for floating wind turbine in different wind conditions [J]. Journal of vibration and shock, 2017, 36(6): 201-206, 228. DOI: 10.13465/j.cnki.jvs.2017.06.031. [9] 倪道俊, 肖瑶瑶. 海上风力发电复合筒型基础拖航稳性研究 [J]. 南方能源建设, 2021, 8(4): 26-31. DOI: 10.16516/j.gedi.issn2095-8676.2021.04.004. NI D J, XIAO Y Y. Research on towing stability of composite bucket foundation for offshore wind power generation [J]. Southern energy construction, 2021, 8(4): 26-31. DOI: 10.16516/j.gedi.issn2095-8676.2021.04.004. [10] 任灏, 马兆荣, 李聪, 等. 海上风电多筒导管架基础湿拖过程稳性控制研究 [J]. 南方能源建设, 2021, 8(增刊1): 65-69. DOI: 10.16516/j.gedi.issn2095-8676.2021.S1.010. REN H, MA Z R, LI C, et al. Research on stability control of offshore wind power multi jacket foundation during wet towing [J]. Southern energy construction, 2021, 8(Suppl.1): 65-69. DOI: 10.16516/j.gedi.issn2095-8676.2021.S1.010. [11] 蔡彦枫, 王海龙, 周川, 等. 基于长期实测资料的风速测量相关推测方法对比 [J]. 南方能源建设, 2017, 4(3): 97-102, 106. DOI: 10.16516/j.gedi.issn2095-8676.2017.03.018. CAI Y F, WANG H L, ZHOU C, et al. Comparison of measure-correlation-predict algorithms in offshore wind power assessment with multi-year observation of automatic weather stations [J]. Southern energy construction, 2017, 4(3): 97-102, 106. DOI: 10.16516/j.gedi.issn2095-8676.2017.03.018. [12] 乐丛欢, 滕丽霞, 李彦娥, 等. 新型全潜式浮式风机在不同风况下的动力响应分析 [J]. 中国海洋大学学报, 2019, 49(增刊1): 120-127. DOI: 10.16441/j.cnki.hdxb.20180406. LE C H, TENG L X, LI Y E, et al. Dynamic response analysis of a new type submersible floating wind turbines under different wind conditions [J]. Journal of Ocean University of China, 2019, 49(Suppl.1): 120-127. DOI: 10.16441/j.cnki.hdxb.20180406. [13] 张延涛, 吕柏呈, 武文华, 等. 基于现场实测的渤海风速特性研究 [J]. 海洋工程, 2020, 38(1): 147-153. DOI: 10.16483/j.issn.1005-9865.2020.01.015. ZHANG Y T, LÜ B C, WU W H, et al. Wind characteristics analysis of Bohai Bay based on field measurement [J]. The ocean engineering, 2020, 38(1): 147-153. DOI: 10.16483/j.issn.1005-9865.2020.01.015. [14] 陈嘉豪, 裴爱国, 马兆荣, 等. 海上漂浮式风机关键技术研究进展 [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. [15] 李聪, 马兆荣, 刘晋超, 等. 海上测风塔整体耦合风荷载计算 [J]. 南方能源建设, 2015, 2(3): 86-90. DOI: 10.16516/j.gedi.issn2095-8676.2015.03.017. LI C, MA Z R, LIU J C, et al. Full-coupled wind load calculation of offshore wind measurement mast [J]. Southern energy construction, 2015, 2(3): 86-90. DOI: 10.16516/j.gedi.issn2095-8676.2015.03.017. [16] 李敏生, 王振华. 中国输电线路规范的风荷载计算比较 [J]. 南方能源建设, 2018, 5(3): 89-93. DOI: 10.16516/j.gedi.issn2095-8676.2018.03.014. LI M S, WANG Z H. Comparison of wind load calculation for China transmission codes [J]. Southern energy construction, 2018, 5(3): 89-93. DOI: 10.16516/j.gedi.issn2095-8676.2018.03.014. [17] 杜齐鲁, 黄海龙, 周益人, 等. 海洋工程试验中API谱特性的风模拟 [J]. 水利水运工程学报, 2016(1): 17-22. DOI: 10.16198/j.cnki.1009-640X.2016.01.003. DU Q L, HUANG H L, ZHOU Y R, et al. API simulation spectrum in ocean engineering based on physical model tests [J]. Hydro-science and engineering, 2016(1): 17-22. DOI: 10.16198/j.cnki.1009-640X.2016.01.003. [18] 美国石油学会. API RP 2A-WSD-2014海上固定平台规划、设计和建造的推荐作法(工作应力设计法): API RP 2A-WSD-2014 [S]. 华盛顿: 美国石油学会, 2014. American Petroleum Institute. Planning, designing and constructing fixed offshore platforms-working stress design: API RP 2A-WSD-2014 [S]. Washington: API, 2014. [19] 扈喆, 张晓莹, 李妍, 等. 改进的浮体运动响应间接时域计算方法 [J]. 哈尔滨工程大学学报, 2021, 42(10): 1482-1489. DOI: 10.11990/jheu.202007010. HU Z, ZHANG X Y, LI Y, et al. Improved indirect time domain method for calculating the motion response of floating bodies [J]. Journal of Harbin Engineering University, 2021, 42(10): 1482-1489. DOI: 10.11990/jheu.202007010. [20] CUMMINS W E. The impulse response function and ship motions [J]. Schiffstechnik, 1962, 9(1): 101-109.