甘 進 吳衛(wèi)國 李旭坤 王 旭

(武漢理工大學交通學院 武漢 430070)

?

海上風電打樁船運動響應分析*

甘 進 吳衛(wèi)國 李旭坤 王 旭

(武漢理工大學交通學院 武漢 430070)

針對雙體打樁船在海上打樁施工中的安全性問題，建立了打樁施工過程中風浪流作用下施工船及錨鏈的耦合系統(tǒng)數(shù)值分析模型，通過計算得到了頻域及時域下船體的運動響應.根據(jù)對施工環(huán)境條件及打樁垂直度要求，計算了打樁船多種施工工況下的運動響應，并對雙體打樁船施工安全性進行了評價.

海上風電；打樁船；垂直度；運動響應；數(shù)值計算；AQWA

0 引 言

風力發(fā)電是近年來世界各國普遍關注的可再生能源開發(fā)項目之一，發(fā)展速度非?？?1997～2004年，全球風電裝機容量年均增長率達26.1%.目前全球風電裝機容量已達5 000萬kW，相當于47座標準核電站.經(jīng)過連續(xù)多年的高速增長，我國風電裝機容量已居世界第1位[1].目前我國正在大力推動海上風電發(fā)展，將從以陸上風電開發(fā)為主向陸上和海上風電全面開發(fā)轉變，目標是成為海上風電大國.

滿足為海上風電施工的不同需求，國家及各類涉風電公司研發(fā)和建造了多型用于海洋風電施工安裝的先進裝備[2].其中風電安裝打樁船作為大型海上施工船舶，其安全性備受關注.傳統(tǒng)的打樁作業(yè)流程為：(1)鋼管樁浮運到指定施工地點；(2)鋼管樁立管作業(yè)以及打定位樁，構建臨時施工平臺，保證樁體的垂直度；(3)在施工平臺上進行打樁作業(yè).采用雙體船進行打樁相比傳統(tǒng)打樁方法而言，其省去了打定位樁以及搭建施工平臺的步驟，縮減了施工周期與成本.與普通單體打樁船相比，雙體打樁船具有較好的穩(wěn)性和操縱性，能較好的保證樁體的垂直度要求.海上打樁施工規(guī)范要求，打樁船施工期間樁體的垂直度偏差不得超過0.3%.樁體的垂直度主要受到船體運動的影響，如果打樁船在波浪中的運動響應過大，則不能進行正常打樁施工.在何種海況條件下可以進行施工作業(yè)，需要進一步計算.雙體船在波浪中的運動響應的計算原理與普通的單體船計算原理是一樣的，都是基于勢流理論求解規(guī)則波下船舶六自由度的幅值響應算子RAO，再通過線性系統(tǒng)的分析方法以及頻率特性得到脈沖響應函數(shù)并進行傅里葉變換，就可以得到船舶在規(guī)則或者不規(guī)則波下的時域運動響應.當前，國內外學者大多采用這種方法開展船舶波浪運動響應的研究.祁祺等[3]利用AQWA對FPSO系泊系統(tǒng)運動響應進行了數(shù)值模擬，給出各自由度方向上的一、二階的運動響應結果， 通過對比分析得到了較為合理的系泊方式，為工程應用提供了參考依據(jù).陳鵬等[4]計算了半潛式平臺水動力性能及系泊性能的計算分析，發(fā)現(xiàn)當環(huán)境力方向在60°附近時，平臺的系泊纜繩張力和水平偏移最大，需要在系泊方案設計中引起重視.Lee等[5]基于脈沖響應函數(shù)理論，創(chuàng)建局部非線性切片理論，計算了雙體船迎浪與隨浪狀態(tài)下的波浪動載荷，并與模型試驗做了對比.

文中以海上風電安裝打樁船“華爾辰”號為研究對象，建立了船體和錨泊耦合系統(tǒng)模型，開展了耦合系統(tǒng)在不同波浪條件下的運動響應分析，以滿足規(guī)范為前提，確定了可施工的海況，為海上施工的順利開展提供了決策依據(jù).

1 打樁船數(shù)值分析模型

雙體打樁船“華爾辰”號的船舶主尺度見表1，根據(jù)打樁船的型線圖、總布置圖等資料建立雙體打樁船幾何模型，離散為面元模型，見圖1.AQWA基于勢流理論和三維二階輻射/繞射切片理論計算水動力學系數(shù)[6].施工中所受環(huán)境力參數(shù)見表2.

表1 雙體打樁船“華爾辰”號主尺參數(shù)

圖1 打樁船網(wǎng)格劃分

參數(shù)數(shù)值波浪有義波高H1/3/m譜峰周期Tp/s角度/(°)0.843.490,90,180風速度/(m·s-1)風譜角度/(°)6Ochi&shin0,-90,-180流速度/(m·s-1)角度/(°)0.50,90,180

系泊纜繩的布置形式，見圖2.系泊纜繩由錨鏈和錨繩構成，參數(shù)見表3.

圖2 系泊錨鏈的水平布置示意圖

頻域下幅值響應算子的計算因受網(wǎng)格限制和實際工程的需要，取0°～180°間7個浪向角，即每隔30°取一個浪向角，頻率取0.1～0.294 Hz中的20個進行計算.

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

2 運動響應頻域分析

海上風電施工中樁體的垂直度是保證海上風機基礎施工質量的重點，故將影響打樁垂直度最敏感的船體縱搖與橫搖響應作為計算分析對象.打樁船運動響應頻域分析結果見圖3.

由圖3a)可知，由于浪向角和船都是對稱的緣故，計算得到的運動響應也基本呈對稱分布，在波浪橫向入射，即入射角為90°時，出現(xiàn)橫搖響應的峰值，且頻率為0.068 Hz(圓頻率0.43 rad/s)，根據(jù)微幅波理論

(1)

極淺水情況下有

(2)

(3)

此時波長約為137.3 m，為低頻情況下需避免的海況.入射波角度為90°為響應出現(xiàn)最大的位置也是合理的.由圖3b)可知對于縱搖響應，順浪與斜浪的時候響應最明顯，而橫浪的響應相比就為小值，也是合理的.縱搖響應在0.068 Hz(圓頻率0.43 rad/s)出現(xiàn)峰值，與橫搖響應相同，在頻率為0.128 Hz(圓頻率0.804 rad/s)每一浪向均出現(xiàn)尖點，且相較峰值來說不小，甚至有浪向為最大，這樣的海況也是要避免的.

圖3 打樁船運動響應頻域分析結果

圖3c)～f)提取的均是附加質量矩陣與阻尼系數(shù)矩陣的主對角線上的值，對于橫搖與縱搖分別是A44,A55,B44,B55.在低頻和高頻段，縱搖的附加質量A55小于橫搖附加質量A44，在中頻段橫搖附加質量增加，出現(xiàn)峰值于0.204 Hz(圓頻率1.28 rad/s).對阻尼系數(shù)而言，低頻段時，縱搖阻尼系數(shù)增長很快，之后又迅速下降，最后趨于一個較小的值，而橫搖阻尼系數(shù)則一直很小，高頻段有兩個頻率發(fā)生了突變，由于文中研究的是雙體船，長寬比較小，橫穩(wěn)性好，縱搖時的阻尼系數(shù)與附加質量較橫搖時大的并不多，也反映出雙體打樁船在保持垂直度時的優(yōu)勢.再對比橫搖與縱搖，橫搖時的附加質量與輻射阻尼變化規(guī)律類似，附加質量要比輻射阻尼大得多.縱搖時，附加質量大體上是隨著頻率增大而呈先下降再增長的趨勢，阻尼則正好相反.

3 運動響應時域分析

由于本船在施工中鋼管樁是直立狀態(tài)，重心很高，在波浪、流與風入射方向相反時，不論是對縱搖或是橫搖都是最不利的情況，所以根據(jù)頻域計算得出的結論，即0°與90°存在對垂直度影響最不利的部分，針對這2種工況進行時域的模擬.施工中允許的垂直度偏差為0.3%，即0.17°，在時域模擬中，橫搖、縱搖角度的幅值不大于0.17°為評價標準.本節(jié)對3種環(huán)境載荷工況下船舶的運動響應進行了分析.圖4～6分別為環(huán)境載荷工況I、II、III的運動響應計算結果.

環(huán)境載荷工況I時，見圖4a)，波流沿船艏向船艉方向入射，而風作用與波流相反，由圖4b)和圖4c)可知，縱搖的運動響應是遠大于橫搖的，幅值為0.051°，此時縱搖響應最大且縱搖幅值遠小于0.17°，滿足打樁垂直度的要求.

環(huán)境載荷工況II時，見圖5a)，風浪和流均沿船舶橫向入射，且風和波流方向相反，此時為對橫搖響應最為不利的工況，由圖5b)～c)可知，橫搖響應幅值為0.278°，是超過了要求的0.17°的施工要求，可見此種海況條件不利于施工要求和安全，應盡量避免.

圖4 環(huán)境載荷工況Ⅰ的運動響應計算結果

圖5 環(huán)境載荷工況Ⅱ的運動響應計算結果

圖6 環(huán)境載荷工況Ⅲ的運動響應計算結果

環(huán)境載荷工況III時，見圖6a)，該工況與環(huán)境載荷工況I類似，只是環(huán)境載荷作用的方向相反，此時仍然是縱搖不利的工況，最大縱搖角依然在安全的范圍之內，而橫搖響應則極小，可忽略不計.

4 結 論

1) 根據(jù)海上風電打樁船“華爾辰”號的實際情況，計算了在給定波浪條件下的該船的運動響應水平.根據(jù)計算結果可知，雙體打樁船在布置時應盡量使船艏迎浪，避免波浪橫向或斜向入射.并且在施工時，應盡量避開波長大，周期長的波浪條件進行施工作業(yè).

2) 文中所采用的基于水動力軟件分析AQWA的海上風電安裝船運動響應數(shù)值分析方法，能用于評估海上風電施工方案的合理性和安全性；計算設計施工工況下的打樁垂直度以及分析惡劣工況下的打樁船運動響應等，為海上風電工程施工提供技術保障.

[1]劉林,葛旭波,張義斌,等.我國海上風電發(fā)展現(xiàn)狀及分析[J].能源技術經(jīng)濟,2012,24(3):66-72.

[2]甘進,王旭,葉毅，等.一種海洋工程輔助施工裝置：ZL201420701483.3[P].2015-04-08.

[3]祁祺,張濤,文攀,等.基于AQWA的FPSO系泊系統(tǒng)響應數(shù)值模擬[J].艦船科學技術,2011,33(12):14-18.

[4]陳鵬,馬駿,黃進浩,等.基于AQWA的半潛式平臺水動力分析及系泊性能計算分析[J].船海工程,2013,42(3):44-47.

[5]LEE C M, JONES H D, CURPHEY R M. Predication of motions and hydrodynamic loads of catama-arns[J].Marine Technology,1973,10(4):392-405.

[6]CHAN H S. Prediction of motion and wave loads of twin-hull ships[J].Marine Structures,1993(6):75-102.

Motion Response Analysis of Offshore Pile Driving Barge for Wind Power Construction

GAN Jin WU Weiguo LI Xukun WANG Xu

(SchoolofTransportation,WuhanUniversityofTechnology,Wuhan430070,China)

In view of the safety problems of pile driving catamaran in pile driving at sea, an analysis model is built to simulate the behaviors of the coupled system between hull and mooring under the condition of wind, wave and current. The motion responses in both frequency and time domains are obtained. According to the requirements of environmental conditions and piling perpendicular degree in construction, the values of motion responses for the pile driving catamaran are calculated under many kinds of working conditions, and the construction safety of the pile driving catamaran is evaluated.

offshore wind power; drilling barge; perpendicularity; motion response; numerical calculation; AQWA

2016-08-12

*國家自然科學基金項目(51409201)、中央高?；究蒲袠I(yè)務費專項(2015-zy-025)資助

U674.32 doi:10.3963/j.issn.2095-3844.2016.05.014

甘進(1982- )：男，博士，講師，主要研究領域為海洋工程結構安全可靠性