張亮，趙玉娜，馬勇，張學(xué)偉，荊豐梅

(哈爾濱工程大學(xué)船舶工程學(xué)院，黑龍江哈爾濱150001)

隨著化石能源的日益枯竭，海上風(fēng)能作為可再生清潔能源已取得了令人矚目的成就［1-4］，它具有風(fēng)速高且穩(wěn)定和能量收益高等優(yōu)勢(shì)［5］。目前，世界上已建海上風(fēng)場(chǎng)近40多處，單個(gè)風(fēng)機(jī)的裝機(jī)容量已高達(dá)5 MW［6］。與近海風(fēng)資源相比，深海風(fēng)資源具有能量密度相對(duì)較高以及開發(fā)區(qū)域廣闊的優(yōu)點(diǎn)，因而深海浮式風(fēng)電場(chǎng)的建設(shè)成為必然發(fā)展趨勢(shì)［7］。海上漂浮式風(fēng)力機(jī)系統(tǒng)主要由上部風(fēng)機(jī)、支撐塔架、浮式基礎(chǔ)和錨泊系統(tǒng)組成。浮式基礎(chǔ)作為塔架和風(fēng)力機(jī)的支撐平臺(tái)，它的性能對(duì)于浮式風(fēng)機(jī)的性能有著重要的影響［8］。Spar式平臺(tái)水線面面積小、吃水深，其底部較大的壓載可降低重心，使平臺(tái)具有良好的穩(wěn)性［9］。國(guó)際上，浮式風(fēng)機(jī)支撐平臺(tái)的研究(平臺(tái)設(shè)計(jì)以及風(fēng)機(jī)系統(tǒng)耦合分析方面)已經(jīng)取得了巨大的進(jìn)步，此外挪威國(guó)家石油公司于2009年在挪威西南海域建成了Spar式海上浮式示范風(fēng)力機(jī)Hywind［10］。國(guó)內(nèi)對(duì)Spar平臺(tái)在海上油氣開采方面的研究以及應(yīng)用已經(jīng)比較成熟。但是與海上油氣平臺(tái)相比，風(fēng)力機(jī)葉輪運(yùn)轉(zhuǎn)引入的巨大氣動(dòng)載荷如推力及扭矩對(duì)Spar風(fēng)力機(jī)平臺(tái)的水動(dòng)力性能有很大影響，上部葉輪氣動(dòng)載荷的引入無(wú)疑增大了Spar式風(fēng)力機(jī)平臺(tái)設(shè)計(jì)以及研究難度，因此國(guó)內(nèi)對(duì)于Spar式浮式風(fēng)機(jī)的研究仍處于起步階段。本文設(shè)計(jì)Spar式浮式風(fēng)機(jī)支撐平臺(tái)，并運(yùn)用頻域水動(dòng)力方法對(duì)其水動(dòng)力性能進(jìn)行數(shù)值模擬，研究有無(wú)垂蕩板、工作水深以及重心高度對(duì)Spar型浮式風(fēng)機(jī)運(yùn)動(dòng)性能影響。

1 風(fēng)力機(jī)及Spar平臺(tái)參數(shù)

1.1 風(fēng)力機(jī)選型

在設(shè)計(jì)海上浮式風(fēng)機(jī)系統(tǒng)時(shí)，首要選擇風(fēng)力機(jī)并確定其參數(shù)。由于NREL5MW水平軸風(fēng)力機(jī)及塔架參數(shù)齊全，本文選用它，其具體參數(shù)見(jiàn)表1。

表15MW水平軸風(fēng)力機(jī)性能參數(shù)Table 1 Parameters for NREL 5MW wind turbine

1.2 Spar式風(fēng)力機(jī)平臺(tái)及垂蕩板設(shè)計(jì)

本文根據(jù)5 MW風(fēng)機(jī)的氣動(dòng)載荷(推力，扭矩等)，考慮風(fēng)力機(jī)系統(tǒng)的穩(wěn)性等方面，經(jīng)過(guò)多次的優(yōu)化設(shè)計(jì)，最終確定Spar式風(fēng)力機(jī)平臺(tái)的主尺度如表2，其結(jié)構(gòu)示意圖見(jiàn)1(a)。

表2 無(wú)垂蕩板Spar風(fēng)力機(jī)平臺(tái)主尺度Table 2 Main parameters for Spar platform without heave-damping plate

圖1 Spar式風(fēng)力機(jī)及垂蕩板結(jié)構(gòu)示意圖Fig.1 Structure diagrams of Spar-type wind turbine

根據(jù)Spar式風(fēng)力機(jī)平臺(tái)水動(dòng)力特性，設(shè)計(jì)了直徑大于Spar柱體直徑兩倍的垂蕩板。垂蕩板模型見(jiàn)圖1(b)。帶垂蕩板的Spar風(fēng)力機(jī)平臺(tái)的具體參數(shù)見(jiàn)表3。

表3 帶垂蕩板Spar風(fēng)力機(jī)平臺(tái)主尺度Table 3 Main parameters for Spar platform with heave-damping plate

2 垂蕩板阻尼分析及其對(duì)Spar式風(fēng)力機(jī)平臺(tái)水動(dòng)力性能的影響

本文中的Spar型風(fēng)機(jī)基礎(chǔ)可視為圓筒，垂蕩周期的簡(jiǎn)化計(jì)算公式如下［11］:

式中:h為吃水，CM為慣性系數(shù)。

從式(1)中可以看出，平臺(tái)基礎(chǔ)吃水越大平臺(tái)垂蕩的固有周期也就越大。通常在海上浮式風(fēng)機(jī)基礎(chǔ)設(shè)計(jì)中，為了避免浮式基礎(chǔ)與波浪產(chǎn)生共振，盡量使基礎(chǔ)垂蕩的固有周期大于波浪的固有周期，以滿足風(fēng)機(jī)安全工作的要求?，F(xiàn)有的Classic Spar一般安裝在波浪周期較小的海域，采用大吃水來(lái)獲得較大的固有周期，而這種做法對(duì)于波浪周期較大的海域則不再適用。為了改善較大波浪周期情況下Spar平臺(tái)的垂蕩運(yùn)動(dòng)性能，一些學(xué)者提出了圓形垂蕩板的應(yīng)用。垂蕩板可以明顯的提高Spar平臺(tái)的垂蕩性能，因此可以不再依賴改變Spar尺寸來(lái)滿足垂蕩要求。

2.1 垂蕩板及Spar式風(fēng)機(jī)平臺(tái)阻尼計(jì)算

滕斌等［12］應(yīng)用勢(shì)流理論計(jì)算了Spar平臺(tái)垂蕩板的水動(dòng)力系數(shù)并與試驗(yàn)結(jié)果進(jìn)行了對(duì)比，分析發(fā)現(xiàn)真實(shí)流體中輻射阻尼在總阻尼中所占的百分比很小，而粘性阻尼是由垂蕩板底部的渦引起的，隨著板的運(yùn)動(dòng)狀態(tài)和波浪的運(yùn)動(dòng)特性的加劇而增加，是垂蕩板的主要阻尼。

本文在計(jì)算垂蕩板對(duì)平臺(tái)運(yùn)動(dòng)響應(yīng)的影響時(shí)，將其簡(jiǎn)化為粘性阻尼加載在整個(gè)平臺(tái)結(jié)構(gòu)上。垂蕩板的縱蕩和垂蕩的粘性阻尼系數(shù)按照下式計(jì)算:

而垂蕩板的縱搖運(yùn)動(dòng)阻尼系數(shù)按照下式計(jì)算:

式中:D為構(gòu)建基準(zhǔn)長(zhǎng)度;Cd為Morison系數(shù);L為結(jié)構(gòu)總長(zhǎng)度;d1，d2分別為結(jié)構(gòu)上下邊緣吃水的深度。由于該垂蕩板結(jié)構(gòu)對(duì)稱，因此縱蕩與橫蕩阻尼系數(shù)相同，縱搖與橫搖阻尼系數(shù)相同。垂蕩板艏搖運(yùn)動(dòng)是一個(gè)很小的量，因此給艏搖阻尼系數(shù)一個(gè)足夠的值即可。Cd的取值是垂蕩板水動(dòng)力性能計(jì)算的關(guān)鍵因素。在哈爾濱工程大學(xué)船模拖曳水池對(duì)該Spar式浮式風(fēng)機(jī)進(jìn)行了模型試驗(yàn)，通過(guò)自由衰減實(shí)驗(yàn)得到無(wú)因次衰減系數(shù)，通過(guò)轉(zhuǎn)化可得到Cd為5，加入到式(2)和(3)中進(jìn)行阻尼系數(shù)的計(jì)算。

運(yùn)用頻域水動(dòng)力分析對(duì)Spar式風(fēng)力機(jī)平臺(tái)進(jìn)行建模及數(shù)值模擬，可得Spar式風(fēng)力機(jī)平臺(tái)在六自由度上的勢(shì)流阻尼曲線如圖2所示。其中，橫(縱)蕩、垂蕩及橫(縱)搖勢(shì)流阻尼給出的是無(wú)因次化后的曲線。由于艏搖勢(shì)流阻尼較小，沒(méi)有進(jìn)行無(wú)因次化，給出了該項(xiàng)系數(shù)的數(shù)值曲線。由圖中可知，六自由度勢(shì)流阻尼均呈先增大后減小的趨勢(shì)。其中，橫(縱)蕩、橫(縱)搖以及艏搖在波頻范圍內(nèi)(1.5 rad/s附近)出現(xiàn)峰值，而垂蕩勢(shì)流阻尼在低頻范圍內(nèi)(0.25 rad/s)出現(xiàn)峰值，之后急劇減小，在波頻范圍內(nèi)響應(yīng)值幾乎為零。

圖2 Spar基礎(chǔ)六自由度勢(shì)流阻尼Fig.2 Potential flow damping of each motion

2.2 垂蕩板對(duì)Spar式風(fēng)力機(jī)平臺(tái)運(yùn)動(dòng)性能的影響

為分析垂蕩板對(duì)Spar式風(fēng)力機(jī)平臺(tái)運(yùn)動(dòng)性能的影響，本文在 0°浪向角下，波浪周期范圍取為 0～125 s，涵蓋了所有常見(jiàn)的波浪周期，對(duì)帶垂蕩板和不帶垂蕩板的Spar風(fēng)力機(jī)平臺(tái)模型進(jìn)行水動(dòng)力性能數(shù)值計(jì)算。由于浮式風(fēng)力機(jī)基礎(chǔ)的艏搖運(yùn)動(dòng)響應(yīng)值很小，在本文以后的內(nèi)容中將予以忽略。其余自由度運(yùn)動(dòng)響應(yīng)計(jì)算結(jié)果如圖3所示。

圖3 2種Spar模型運(yùn)動(dòng)響應(yīng)比較Fig.3 Response comparison between two types of Spar turbine platform

從圖3中可以看出，垂蕩板對(duì)風(fēng)力機(jī)平臺(tái)在平面內(nèi)的運(yùn)動(dòng)，即縱(橫)蕩及縱(橫)搖的運(yùn)動(dòng)幅值以及固有周期幾乎沒(méi)有影響，而對(duì)垂蕩方向運(yùn)動(dòng)幅值和固有周期有很大影響。具體來(lái)說(shuō)，引入垂蕩板后，垂蕩運(yùn)動(dòng)運(yùn)動(dòng)曲線較為平坦，無(wú)陡峭的峰值出現(xiàn)，且引入垂蕩板后垂蕩響應(yīng)幅值的峰值是未引入垂蕩板時(shí)垂蕩響應(yīng)幅值峰值的36%，即垂蕩響應(yīng)峰值減小至原來(lái)的1/3。由此可以看出，垂蕩板大幅度地減小了平臺(tái)在垂蕩方向的運(yùn)動(dòng)幅值。而且從圖3(c)可以看出，垂蕩板稍微增大了垂蕩運(yùn)動(dòng)的固有周期，這是由于垂蕩板的周圍區(qū)域內(nèi)會(huì)產(chǎn)生一定的漩渦進(jìn)而增加了垂蕩阻尼，導(dǎo)致周期增大。因此，在設(shè)計(jì)浮式風(fēng)機(jī)基礎(chǔ)時(shí)，柱體底部增加垂蕩板是很有必要的，能夠明顯改善平臺(tái)的垂蕩運(yùn)動(dòng)性能，使其具有更好的垂蕩水動(dòng)力性能。

3 工作水深對(duì)Spar式風(fēng)力機(jī)平臺(tái)運(yùn)動(dòng)性能的影響

由于Spar式平臺(tái)與海底必須存在一定的距離，以便于系泊系統(tǒng)的設(shè)計(jì)，因此它可以用于不同水深的海域。工作水深對(duì)Spar式風(fēng)力機(jī)平臺(tái)的運(yùn)動(dòng)響應(yīng)有一定的影響，計(jì)算不同水深下平臺(tái)的水動(dòng)力性能對(duì)浮式風(fēng)機(jī)支撐平臺(tái)的設(shè)計(jì)有著重要的意義。本文選取200、280、500、1 000 m4 個(gè)水深對(duì) Spar式風(fēng)力機(jī)平臺(tái)的水動(dòng)力響應(yīng)進(jìn)行計(jì)算，計(jì)算結(jié)果如圖4所示。表4對(duì)比了不同水深的三自由度最大幅值。

表4 2種Spar模型運(yùn)動(dòng)響應(yīng)比較Table 4 Response comparison between two types of Spar turbine platform

圖4 不同水深Spar運(yùn)動(dòng)響應(yīng)比較Fig.4 Response comparison of Spar turbine platform in different operating ocean depths

結(jié)合圖4和表4，可以看出各自由度運(yùn)動(dòng)響應(yīng)幅值隨水深的增加而遞減，縱蕩和縱搖受水深影響不大。垂蕩響應(yīng)幅值在水深小于500 m時(shí)受水深影響較大，而當(dāng)水深大于500 m(是平臺(tái)設(shè)計(jì)吃水4倍)時(shí)垂蕩響應(yīng)幅值變化不大。這是由于Spar式風(fēng)力機(jī)平臺(tái)吃水較深，適合在深水海域(水深是平臺(tái)設(shè)計(jì)吃水4倍的海域)工作。

4 重心高度對(duì)Spar式風(fēng)力機(jī)平臺(tái)運(yùn)動(dòng)性能的影響

一般情況下，Spar平臺(tái)的重心都在浮心之下以保證其穩(wěn)性，通過(guò)調(diào)節(jié)壓載可以改變Spar平臺(tái)的重心位置。為分析裝載壓載水之后Spar式風(fēng)力機(jī)平臺(tái)的重心位置對(duì)平臺(tái)的水動(dòng)力性能的影響。本文選取了一系列不同重心高度分別計(jì)算平臺(tái)運(yùn)動(dòng)響應(yīng)，研究平臺(tái)重心高度對(duì)運(yùn)動(dòng)響應(yīng)的影響。選取的重心高度分別為-70、-77.4、-77.6、-77.8、-79、-80、-82、-83、-84、-85、-86、-87、-88、-90 m(水線面處高度為 0 m)。將其按1～14N排列，設(shè)為橫坐標(biāo)，將每個(gè)重心高度對(duì)應(yīng)的共振幅值設(shè)為縱坐標(biāo)，做成對(duì)比如圖5所示。

圖5 不同重心下各自由度運(yùn)動(dòng)響應(yīng)Fig.5 Motion response v.s.height of center of gravity

從圖5中可以看出重心位置對(duì)垂蕩運(yùn)動(dòng)響應(yīng)沒(méi)有影響?？v(橫)蕩及縱(橫)搖的運(yùn)動(dòng)響應(yīng)隨重心位置變化的趨勢(shì)基本相同，當(dāng)重心位置位于-77.4～-82 m之間時(shí)，這4個(gè)自由度上的運(yùn)動(dòng)響應(yīng)基本變化很小。而當(dāng)重心位置大于-82 m之后，運(yùn)動(dòng)響應(yīng)幅值急劇減小，這4個(gè)自由度上的運(yùn)動(dòng)均在重心在水線面下-82 m處出現(xiàn)最大值，最大值分別達(dá)到6.19 m和28.3°，這說(shuō)明Spar平臺(tái)重心在水線面以下82 m處共振最為劇烈。重心位置達(dá)到水線面以下90 m時(shí)，運(yùn)動(dòng)響應(yīng)幅值又急劇增大。將重心位置對(duì)平臺(tái)的設(shè)計(jì)吃水進(jìn)行無(wú)量綱化處理之后，可以看出，Spar式風(fēng)力機(jī)平臺(tái)重心高度設(shè)置在0～0.8倍的設(shè)計(jì)重浮心距離范圍內(nèi)時(shí)，Spar平臺(tái)的水動(dòng)力性能良好。

5 結(jié)論

本文針對(duì)NREL5MW風(fēng)力機(jī)設(shè)計(jì)了Spar式風(fēng)力機(jī)平臺(tái)，并采用頻域水動(dòng)力分析方法對(duì)Spar式風(fēng)力機(jī)平臺(tái)的水動(dòng)力性能進(jìn)行數(shù)值模擬，分析平臺(tái)有無(wú)縱蕩板、平臺(tái)的重心高度以及平臺(tái)的工作水深對(duì)平臺(tái)水動(dòng)力性能的影響，研究結(jié)果表明:

1)垂蕩板可顯著提高平臺(tái)的垂蕩性能，可以降低平臺(tái)垂蕩的響應(yīng)幅值峰值到無(wú)垂蕩板時(shí)峰值的1/3，但對(duì)于橫(縱)蕩及橫(縱)搖無(wú)明顯影響。

2)隨著工作水深的增加，橫(縱)蕩、橫(縱)搖影響較小;當(dāng)工作水深小于4倍平臺(tái)設(shè)計(jì)吃水時(shí)，水深對(duì)垂蕩響應(yīng)幅值影響較大，而當(dāng)工作水深大于4倍平臺(tái)設(shè)計(jì)吃水時(shí)，此影響變得不明顯，這主要是因?yàn)镾par式風(fēng)力機(jī)平臺(tái)適合在深水區(qū)域(水深大于4倍平臺(tái)設(shè)計(jì)吃水)工作。

3)重心高度的變化幾乎對(duì)垂蕩的運(yùn)動(dòng)無(wú)影響;隨著平臺(tái)重心的降低，平臺(tái)的橫(縱)蕩、橫(縱)搖的響應(yīng)幅值變化較大，當(dāng)加壓載之后的平臺(tái)重心低于設(shè)計(jì)重心0～0.8倍的設(shè)計(jì)重浮心距離范圍內(nèi)時(shí)，Spar平臺(tái)的水動(dòng)力性能較好。因此，實(shí)際操作中應(yīng)注意控制壓載水量，如果加壓載之后的平臺(tái)重心高度與設(shè)計(jì)重心高度偏差超出上述范圍，平臺(tái)的橫(縱)蕩、橫(縱)搖運(yùn)動(dòng)響應(yīng)峰值將會(huì)增大5倍以上，嚴(yán)重降低平臺(tái)的水動(dòng)力性能。

［1］郭越，王占坤.中歐海上風(fēng)電產(chǎn)業(yè)發(fā)展比較［J］.中外能源，2011，16(3):26-30.GUO Yue，WANG Zhankun.A comparison between European and Chinese offshore wind power industries［J］.Sino-Global Energy，2011，16(3):26-30.

［2］賀德馨.我國(guó)風(fēng)工程研究現(xiàn)狀與展望［J］.力學(xué)與實(shí)踐，2002，24(4):10-19.HE Dexin.The status and vista of wind engineering studies in China［J］.Mechanics in Engineering，2002，24(4):10-19.

［3］李曉燕，余志.海上風(fēng)力發(fā)電進(jìn)展［J］.太陽(yáng)能學(xué)報(bào)，2004，25(1):78-84.LI Xiaoyan，YU Zhi.Development of offshore wind power［J］.Acta Energiae Solaris Sinica，2004，25(1):78-84.

［4］張亮，吳海濤，荊豐梅，等.海上漂浮式風(fēng)力機(jī)研究進(jìn)展及發(fā)展趨勢(shì)［J］.海洋技術(shù)，2010，29(4):122-125.ZHANG Liang，WU Haitao，JING Fengmei，et al.Study on offshore wind turbine and its development［J］.Ocean Technology，2010，29(4):122-125.

［5］WANG Yi，DUAN Menglan，SHANG Jinghong.Application of an abandoned jacket for an offshore structure base of wind turbine in Bohai heavy ice conditions［C］//Proceedings of the 19th International Offshore and Polar Engineering Conference.Osaka，Japan，2009:384-389.

［6］TWIDEL J，GAUDIOSI G.Offshore wind power［M］.Brentwood:Multi-science Publishing Co.，Ltd，2009:2-14.

［7］張德.風(fēng)能資源數(shù)值模擬及其在中國(guó)能源評(píng)估中的應(yīng)用研究［D］.蘭州:蘭州大學(xué)，2009:14-36.ZHNAG De.Study on wind energy numerical simulation and its application to wind energy resources assessment in China［D］.Lanzhou:Lanzhou Unviersity，2009:14-36.

［8］葉小嶸，張亮，吳海濤，等.平臺(tái)運(yùn)動(dòng)對(duì)海上浮式風(fēng)機(jī)的氣動(dòng)性能影響研究［J］.華中科技大學(xué)學(xué)報(bào)，2012，40(3):123-126.YE Xiaorong，ZHANG Liang，WU Haitao，et al.Influence of platform motion response on aerodynamic performance of floating offshore wind turbine［J］.Journal of Huazhong University of Science and Technology，2012，40(3):123-126.

［9］WITHEE J E.Fully coupled dynamic analysis of a floating wind turbine system［M］.Cambridge:Massachusetts Institute of Technology，2004.

［10］NIELSEN F G，HANSON T D，SKAARE B.Integrated dynamic analysis of floating offshore wind turbines［C］//Proceedings of 25th International Conference on Offshore Mechanics and Arctic Engineering，Hamburg，Germany，2006.

［11］THIAGARAJAN K P.Influence of heave plate geometry on the heave response of classic spars［C］//Proceedings of 21st International Conference on Offshore Mechanics and Artic Engineering.Oslo，Norway，2002:621-627.

［12］滕斌，鄭苗子，姜?jiǎng)俪?，?Spar平臺(tái)垂蕩板水動(dòng)力系數(shù)計(jì)算與分析［J］.海洋工程，2010，8:1-4.TENG Bin，ZHENG Miaozi，JIANG Shengchao，et al.Calculation and analysis of the hydrodynamic coefficients of heave-plates of Spar platform［J］.The Ocean Engineering，2010，8:1-4.