蒙宣伊，陽(yáng)航

（湘電風(fēng)能有限公司，湖南 湘潭 411100）

Spar式風(fēng)機(jī)基礎(chǔ)系統(tǒng)水動(dòng)力性能研究

蒙宣伊，陽(yáng)航

（湘電風(fēng)能有限公司，湖南 湘潭 411100）

本文基于三維水動(dòng)力學(xué)軟件Aqwa進(jìn)行了Spar式風(fēng)機(jī)基礎(chǔ)系統(tǒng)的水動(dòng)力性能研究。通過(guò)時(shí)域方法，研究系統(tǒng)在額定風(fēng)速工況下的運(yùn)動(dòng)響應(yīng)。計(jì)算時(shí)考慮風(fēng)浪聯(lián)合作用的影響。最后通過(guò)傅里葉變換，得到升沉、縱搖、縱蕩和錨鏈拉力響應(yīng)譜。

Spar式基礎(chǔ)；水動(dòng)力性能；響應(yīng)譜

海上風(fēng)能被公認(rèn)為是一種可以用來(lái)滿足能量增長(zhǎng)需求的可再生能源。相比海洋中其它可再生能源，比如潮汐能和波浪能，風(fēng)能的開(kāi)發(fā)及相關(guān)技術(shù)被認(rèn)為是成熟的，而且建設(shè)相當(dāng)好。其中大部分已建成并運(yùn)行的風(fēng)場(chǎng)主要是以固定式基礎(chǔ)形式，而且水深比較淺。對(duì)于每個(gè)可能建成的風(fēng)場(chǎng)來(lái)說(shuō)，其取決于波浪和風(fēng)特征、海床特性以及社會(huì)條件。在某一水深，選擇使用何種基礎(chǔ)時(shí)，主要考慮成本相關(guān)的問(wèn)題。

相比傳統(tǒng)固定式基礎(chǔ)，漂浮式基礎(chǔ)整體系統(tǒng)的性能研究是十分必要的，主要原因如下。

（1）它們的固有頻率非常低，通常會(huì)影響氣動(dòng)阻尼和穩(wěn)定性。

（2）對(duì)于半潛式和Spar來(lái)說(shuō)，它們的位移和旋轉(zhuǎn)運(yùn)動(dòng)會(huì)與機(jī)艙、葉輪的運(yùn)動(dòng)相互耦合。

（3）它們錨固在海床上的錨鏈系統(tǒng)必須包含在整體分析中。

Nielsen等對(duì)Spar基礎(chǔ)整體動(dòng)力分析進(jìn)行了研究。他們對(duì)Hywind的基礎(chǔ)進(jìn)行仿真，并將結(jié)果與縮比模型的試驗(yàn)結(jié)果進(jìn)行對(duì)比。Matsukuma和Utsunimiya采用多體動(dòng)力學(xué)理論對(duì)一種漂浮式基礎(chǔ)在恒定風(fēng)速下考慮葉輪旋轉(zhuǎn)時(shí)的運(yùn)動(dòng)響應(yīng)。Jokman等在OC3項(xiàng)目中對(duì)固定式和漂浮式基礎(chǔ)的結(jié)構(gòu)動(dòng)態(tài)響應(yīng)進(jìn)行了驗(yàn)證。Karmirad和Moan采用混合 aero-hydro-elastic時(shí)域方法進(jìn)行了一種Spar式基礎(chǔ)在極限情況下的結(jié)構(gòu)響應(yīng)研究。Karmake和Carlos等對(duì)采用長(zhǎng)期預(yù)報(bào)方法進(jìn)行了Spar和半潛式基礎(chǔ)的運(yùn)動(dòng)響應(yīng)對(duì)比分析。

本文給出了5MW風(fēng)機(jī)的Spar基礎(chǔ)的參數(shù)和外形，并進(jìn)行了錨泊系統(tǒng)的建模。采用時(shí)域方法，將風(fēng)機(jī)推力假定為一恒定力，計(jì)算系統(tǒng)在額定風(fēng)速情況下的運(yùn)動(dòng)響應(yīng)情況。通過(guò)傅里葉變換，得到各運(yùn)動(dòng)和錨鏈拉力的響應(yīng)譜。

1 Spar基礎(chǔ)模型和風(fēng)機(jī)系統(tǒng)

Spar式風(fēng)機(jī)系統(tǒng)包含Spar式基礎(chǔ)、塔筒、機(jī)艙和葉輪?；A(chǔ)吃水為120m，其中底段浮筒直徑為9.4m，長(zhǎng)度為108m，過(guò)渡段上部直徑為6.5m，長(zhǎng)度為8m?；A(chǔ)平臺(tái)的質(zhì)量（包含壓載）為7466330kg，重心高度為-45m。塔筒重量為347460kg，重心高度為海平面以上45m。葉輪重量和分布采用NREL5MW風(fēng)機(jī)模型，總重量為350000kg，重心高度為90m（MSL）。系統(tǒng)重心為相對(duì)平均海平面（0,0，-62.82）?；A(chǔ)和風(fēng)機(jī)模型參數(shù)如圖1和表1所示。

圖1 Spar基礎(chǔ)及錨泊系統(tǒng)

表1 風(fēng)機(jī)和基礎(chǔ)模型參數(shù)

進(jìn)行時(shí)域計(jì)算時(shí)，通過(guò)定義風(fēng)推力系數(shù)，將風(fēng)載簡(jiǎn)化為以某一推力作用在輪轂中心高度，即海平面以上90m。Spar基礎(chǔ)和錨泊系統(tǒng)如圖1所示。錨泊系統(tǒng)采用3根錨鏈互相成120度散布方式進(jìn)行設(shè)計(jì)，導(dǎo)纜孔的位置在水下70m，計(jì)算水深為320m。錨鏈直徑為0.09m，錨鏈總長(zhǎng)為902.2m。錨泊系統(tǒng)具體見(jiàn)圖1，具體參數(shù)如表2所示。

表2 錨泊系統(tǒng)參數(shù)

2 Spar基礎(chǔ)數(shù)值分析

2.1 工況設(shè)置

基礎(chǔ)結(jié)構(gòu)的時(shí)域分析是運(yùn)用Aqwa時(shí)域方法進(jìn)行的，在計(jì)算時(shí)考慮恒定風(fēng)推力對(duì)運(yùn)動(dòng)響應(yīng)的影響。在額定風(fēng)速工況下，風(fēng)速為11.2m/s，有義波高為3m，譜峰周期為10s，風(fēng)推力為720kN。計(jì)算所采用的波浪譜為Jonswap譜，其公式如式（1）所示：

式中：γ為升高因子，取值范圍為1～3.3；

σ為峰形參數(shù)，若ω≤ωp時(shí)，σ=0.07，若ω＞ωp時(shí)，σ=0.09；

計(jì)算步長(zhǎng)設(shè)置為1s，總時(shí)長(zhǎng)為3600s。

2.2 額定風(fēng)速工況下運(yùn)動(dòng)響應(yīng)分析

圖2～4分別為基礎(chǔ)在額定風(fēng)速情況下的升沉、縱搖、垂蕩和各錨鏈截面拉力的時(shí)歷曲線。對(duì)圖2～4來(lái)說(shuō)，藍(lán)線表示的是波浪作用下的響應(yīng)，而紅色為實(shí)際響應(yīng)曲線。圖2可以看出，基礎(chǔ)的實(shí)際響應(yīng)在風(fēng)浪聯(lián)合初始作用下產(chǎn)生一個(gè)約0.55m向下偏移。從圖3看出，在風(fēng)作用下，基礎(chǔ)實(shí)際有一個(gè)2.5度的初始傾角，而且實(shí)際運(yùn)動(dòng)響應(yīng)幅值是要大于僅考慮波浪作用時(shí)的響應(yīng)。圖4可以看出，基礎(chǔ)在風(fēng)作用下縱向會(huì)有8.78m左右的平均偏移。

圖2 額定風(fēng)速工況下升沉?xí)r歷曲線

圖3 額定風(fēng)速工況下縱搖時(shí)歷曲線

圖4 額定風(fēng)速工況下縱蕩時(shí)歷曲線

圖5給出了3個(gè)方向的錨鏈截面拉力的時(shí)歷曲線。其中1號(hào)錨鏈的方向?yàn)?度，因此該錨鏈所受的拉力要大于2、3號(hào)錨鏈。2、3號(hào)錨鏈由于基礎(chǔ)系統(tǒng)的對(duì)稱性，其所受的拉力時(shí)歷曲線基本重合，因此圖中用同一條曲線來(lái)表示。

2.3 響應(yīng)頻譜

圖5 額定風(fēng)速工況下各錨鏈截面拉力時(shí)歷曲線

對(duì)于仿真得到的運(yùn)動(dòng)和拉力時(shí)歷曲線，通過(guò)傅里葉變換，可以得到其響應(yīng)頻譜。圖6～8分別為升沉、縱搖和縱蕩的運(yùn)動(dòng)響應(yīng)頻譜。

圖6 升沉響應(yīng)頻譜

從圖6可以看出，基礎(chǔ)系統(tǒng)的升沉響應(yīng)頻譜存在兩峰值點(diǎn)，左側(cè)峰值處于系統(tǒng)升沉固有周期附近，右側(cè)峰值對(duì)應(yīng)頻率處于譜峰周期所對(duì)應(yīng)頻率附近。波浪所誘導(dǎo)的升沉響應(yīng)譜要略大于實(shí)際響應(yīng)頻譜。

圖7 縱搖響應(yīng)頻譜

從圖7可以看出，無(wú)論是否考慮風(fēng)推力，在譜峰周期附近有一個(gè)峰值點(diǎn)。在風(fēng)推力的作用下，基礎(chǔ)的實(shí)際縱搖響應(yīng)要大于僅僅考慮波浪作用情況下的響應(yīng)。

圖8 縱蕩響應(yīng)頻譜

從圖8可以看出，實(shí)際的縱蕩響應(yīng)要略小于僅波浪作用下的情況。

圖9為錨鏈張力響應(yīng)頻譜。從圖上可以看出錨鏈1的峰值是要大于2、3號(hào)錨鏈。對(duì)于2、3號(hào)錨鏈來(lái)說(shuō)，存在兩個(gè)峰值，其中一個(gè)處于譜峰周期附近，另一個(gè)出于系統(tǒng)固有周期附近。

圖9 錨鏈拉力響應(yīng)頻譜

3 結(jié)語(yǔ)

本文基于三維水動(dòng)力學(xué)軟件對(duì)Spar式風(fēng)機(jī)基礎(chǔ)系統(tǒng)的水動(dòng)力性能進(jìn)行了研究。從運(yùn)動(dòng)時(shí)歷來(lái)看，在額定風(fēng)速工況下運(yùn)動(dòng)響應(yīng)比較小，性能比較優(yōu)越。

在風(fēng)、浪聯(lián)合作用情況下，基礎(chǔ)的運(yùn)動(dòng)響應(yīng)譜存在雙峰值，其中一個(gè)位于系統(tǒng)固有頻率附近，另一個(gè)位于譜峰周期附近，這與實(shí)際比較吻合。錨鏈1的拉力響應(yīng)大于錨鏈2、3的響應(yīng)。

風(fēng)推力對(duì)基礎(chǔ)系統(tǒng)的響應(yīng)頻譜有影響。在風(fēng)推力作用下，升沉和縱蕩的實(shí)際響應(yīng)要略小于僅波浪誘導(dǎo)的響應(yīng)，而實(shí)際縱搖響應(yīng)要大于僅波浪誘導(dǎo)的縱搖響應(yīng)。

[1]Karimirad M., C. Michailides. V-shaped semisubmersible offshore wind turbine: An alternative concept for offshore wind technology[J]. Renewable Energy 85(2015)126-143.

[2]Breton SF, Moe G. Status, plans and technologies for offshore wind turbines in Europe and North America[J]. Renew Energy 2009;34(3):646-54.

[3]Shi W, Park H, Han J, Na S, Kim C. A study on the effect of different modeling parameters on the dynamic response of a jacket-type offshore wind turbine in the Korean Southwest Sea[J].Renew Energy 2013; 58:50-9.

[4]Shi W, Han J, Kim C, Lee D, Shin H, Park H. Feasibility study of offshore wind turbine substructures for southwest offshore wind farm project in Korea[J]. Renew Energy 2015; 74:406-13.

[5]Musial W, Butterfield S, Ram B. Energy from offshore wind[C]. Texas,USA: Offshore Technology Conference Houston; 2006. p. 1888-98.

[6]Myhr A, Bjerkseter C, ?gotnes A, Nygaard TA. Levelised cost of energy for offshore floating wind turbines in a life cycle perspective[J]. Renew Energy 2014; 66:714-28.

[7]Sun X, Huang D, Wu G. The current state of offshore wind energy technology development[J]. Energy 2012; 41:298-312.

[8]Nielsen, F. G., Hanson, T. D., and Skaare, B., 2006, Integrated Dynamic Analysis of Floating Offshore Wind Turbines[C].ASME’ Paper No. OMAE2006-92291.

[9]Matsukuma, H., and Utsunomiya, T., 2008, Motion Analysis of a Floating Offshore Wind Turbine Considering Rotor Rotation[J].IES J., Part A, 1(4), pp.268–279.

[10]Jonkman, J., and Musial, W., 2010, quot;Offshore Code Comparison Collaboration (OC3) for IEA Wind Task 23, Offshore Wind Technology and Deployment, National Renewable Energy Laboratory, Technical Report No. NREL/TP-5000-48191.

[11]Karimirad, M., and Moan, T., 2010, Extreme Structural Dynamic Response of a Spar Type Wind Turbine[C].’ASME’ Paper No.OMAE2010-20044.

[12]K.Bagbanci, D. Karmakar., C. Guedes Soares. Comparison of Spar and Semisubmersible Floater Concepts of Offshore Wind Turbines Using Long Term Analysis[J]. Journal of Offshore Mechanics and Arctic Engineering,2015,12, Vol.137.

[13]AQWA User’s Manual 2013. ANSYS.

[14]Design of Offshore Wind Turbine Structures. DNV-OS-J101. DNV.

P752 TK89

A

1671-0711（2017）12（上）-0162-04