石陸豐，程建生，段金輝，隋元松

（解放軍理工大學(xué)野戰(zhàn)工程學(xué)院，南京 210007）

0 引 言

近年來隨著能源枯竭問題的日益嚴(yán)峻，海上風(fēng)能的開發(fā)受到重視，并已逐漸由淺海向深海發(fā)展。浮基風(fēng)電平臺是深海風(fēng)電機(jī)組的基礎(chǔ)結(jié)構(gòu)。由于浮基風(fēng)電系統(tǒng)重心高，受到風(fēng)、浪、流及冰荷載等的影響，受力情況復(fù)雜，使得浮基風(fēng)電平臺的設(shè)計面臨很多技術(shù)上挑戰(zhàn)和理論應(yīng)用的難點。目前研究較多的浮基風(fēng)電平臺主要有 Spar結(jié)構(gòu)基礎(chǔ)[1]、張力腿結(jié)構(gòu)基礎(chǔ)( tension leg platform，TLP)[2,3]、半潛式結(jié)構(gòu)基礎(chǔ)和浮筒( 或駁船) 基礎(chǔ)。

半潛式浮基風(fēng)電平臺由于具有水深適用范圍廣、安裝工作量少、總成本低等優(yōu)點，受到了各國研究人員的廣泛關(guān)注。Zambrano T[4]等人開發(fā)了一種支撐3個風(fēng)力機(jī)的半潛式平臺MiniFloat，采用基于墨西哥的暴風(fēng)模型傅里葉譜進(jìn)行了深海風(fēng)電浮式基礎(chǔ)結(jié)構(gòu)的運動響應(yīng)分析，通過采用龍格-庫塔法得到了平臺結(jié)構(gòu)在六個自由度上的響應(yīng)，并求解出作用在結(jié)構(gòu)上的風(fēng)浪荷載和錨鏈承受力。而后美國又開發(fā)了單風(fēng)力機(jī)半潛式平臺WindFloat[5～7]，并針對平臺的流體和結(jié)構(gòu)動力特性分別展開了研究。Thomas Zambrano等人[8]基于Wamit軟件平臺計算作用在多風(fēng)機(jī)半潛平臺上的波浪力，并分析了波浪作用下半潛平臺的六個自由度運動響應(yīng)。唐友剛[9]針對所設(shè)計的水深為60m的600kW水平軸半潛式海上風(fēng)力發(fā)電機(jī)，使用有限元法建立了包括風(fēng)機(jī)、塔架、浮式基礎(chǔ)和系泊系統(tǒng)的數(shù)學(xué)模型，對波浪載荷作用下的浮式基礎(chǔ)動力響應(yīng)進(jìn)行了頻域分析。研究發(fā)現(xiàn)，平臺的運動會導(dǎo)致風(fēng)力機(jī)功率的波動，引起電力輸出的不穩(wěn)定[10]。因此，在海上浮式風(fēng)力機(jī)的設(shè)計中，平臺的運動特性應(yīng)作為主導(dǎo)因素考慮。

為了提升半潛式浮基平臺的穩(wěn)性，采用懸掛壓載物降低平臺重心的方法設(shè)計了一種半潛式浮基風(fēng)電平臺，并應(yīng)用ADINA軟件，基于CFD方法，采用邊界造波法、ALE動網(wǎng)格技術(shù)以及網(wǎng)格變疏增大數(shù)值黏性的消波法構(gòu)建了數(shù)值波浪水池，對平臺在波浪中的運動進(jìn)行數(shù)值模擬。

1 浮基風(fēng)電平臺概念化設(shè)計

受風(fēng)電機(jī)組的影響，深海浮式基礎(chǔ)系統(tǒng)較普通石油鉆井平臺而言傾覆力矩較大，為了保證風(fēng)機(jī)發(fā)電的穩(wěn)定和整個系統(tǒng)部件免受疲勞損傷，設(shè)計在平臺的下部懸掛大重量的壓載物來降低整個系統(tǒng)的重心以抑制整個平臺在各種荷載作用下的晃動。此壓載物設(shè)計總質(zhì)量 2400t，體積 1000m3，尺寸為 10m×10m×10m的立方體，材料為防腐鋼筋混凝土。浮基風(fēng)電平臺概念化設(shè)計見圖1、2，設(shè)計的主要參數(shù)見表1。

圖1 深海浮基風(fēng)電平臺概念化設(shè)計

圖2 深海浮基風(fēng)電平臺模型

表1 浮基風(fēng)電平臺概念化設(shè)計的主要參數(shù)

2 波浪水池的數(shù)學(xué)模型

2.1 流體模型

流體控制方程用非定常不可壓縮雷諾平均 Navier-Stokes（RANS）方程求解。為了跟蹤固體邊界和流體邊界的運動狀況，很好地反映耦合界面的網(wǎng)格變化，采用任意拉格朗日-歐拉（Arbitrary Lagrangian-Eulerian description）坐標(biāo)系。ALE描述下的連續(xù)性方程和動量方程為：

2.2 造波與消波

采用線性波數(shù)值造波方法實現(xiàn)造波。線性波理論是最常用的波浪理論之一，其表達(dá)形式簡單，使用方便，線性波主要用于小振幅波。計算中采用自由面位移的理論公式為：

波浪水質(zhì)點在各方向的速度為：

色散關(guān)系為：

式中：H——波高，h——靜水深，λ——波長，k=2π/λ——波數(shù)，g——重力加速度，ω=2π/T——圓頻率，T——波浪周期，c——波速。依據(jù)上述公式可知，給定模擬波浪的周期、波高和靜水深，便可以得到入射邊界上邊界點的速度值，將其賦予邊界進(jìn)行造波。

利用ADINA軟件建立計算所用的三維幾何模型，入射波波高H=6m，波長λ=100m ，水深h=200m，根據(jù)色散關(guān)系式（7）求得波浪周期T=6s，同時，λ/2＜h＜∞，屬于深水波。通過對數(shù)值水池某時刻的波形和某位置的波幅歷時時程進(jìn)行監(jiān)控。圖3所示為采用本文數(shù)值方法所造波在t=48s時刻波面變形的三維形態(tài)。在距離入口100m、200m處設(shè)置了兩個觀測點進(jìn)行波面的歷時時程監(jiān)控，兩監(jiān)測點波面歷時時程數(shù)值解與解析解對比（見圖4）。從圖4可知，所構(gòu)建的三維波浪數(shù)值模型波面的數(shù)值解和理論解吻合很好。消波區(qū)采用網(wǎng)格變疏增大黏性的數(shù)值黏性消波方法，即對消波區(qū)段人為地定義較大的黏性來增大消波區(qū)域段的阻尼，達(dá)到消波的目的。

圖3 t=48s時刻線性波三維波面

圖4 100m、200m監(jiān)測點處波面歷時曲線對比

2.3 浮基風(fēng)電平臺運動與流體的耦合

浮基風(fēng)電平臺的六個自由度控制方程表達(dá)式如下：

式中：下標(biāo)c——浮基風(fēng)電平臺的轉(zhuǎn)動中心，m——平臺的質(zhì)量，——速度矢量，Ic——慣性矩張量，——角速度矢量，和——作用在平臺上的力和力矩矢量。

當(dāng)作用在平臺上的力和力矩已知時，其運動的速度、角速度及位置、姿態(tài)等可以通過式（8）、（9）求解積分得到。

3 數(shù)值模型驗證

采用 ADINA軟件進(jìn)行數(shù)值波浪水池的構(gòu)建及驗證和深海風(fēng)電浮式基礎(chǔ)在波浪作用下的運動響應(yīng)研究。依據(jù)第20屆ITTC半潛平臺原型1:70的尺寸建立了數(shù)值計算模型，進(jìn)行波浪作用下三維模型的數(shù)值模擬，并將計算結(jié)果通過分析整理，與模型試驗結(jié)果比較，以驗證構(gòu)建的數(shù)值波浪水池的可靠性。模型結(jié)構(gòu)見圖5。

數(shù)值波浪水池計算區(qū)域范圍為1000m×400m×200m，浮體位于區(qū)域中間，吃水深度為24m，計算區(qū)域分3部分，波浪生成區(qū)長400m，工作區(qū)長200m，消波區(qū)長400m，見圖6。

采用四面體網(wǎng)格，整個流固耦合區(qū)域網(wǎng)格總數(shù)（包括流體和固體）為509586個，計算區(qū)域網(wǎng)格劃分情況見圖7。

圖5 ITTC半潛平臺數(shù)值模型

圖6 ITTC平臺模型分析區(qū)域

圖7 計算區(qū)域網(wǎng)格劃分

水池左側(cè)為入射邊界，波浪水池側(cè)壁和水池底部結(jié)構(gòu)邊界條件都設(shè)置為固壁邊界，與半潛平臺接觸的流體表面定義為FSI邊界(即流固耦合邊界)，而自由表面邊界設(shè)置在整個流體模型的上表面處。

初始條件：風(fēng)電浮式基礎(chǔ)平臺靜止，不受外力的作用；壓力設(shè)為靜水壓；時間步長取為波浪周期的1/200。針對下述波浪環(huán)境條件，將ADINA軟件計算得到的數(shù)值模擬計算結(jié)果與在哈爾濱工程大學(xué)深水池所做試驗的結(jié)果進(jìn)行對比（見圖8，表2）：

表2 規(guī)則波作用下的數(shù)值模擬工況

圖8 ITTC半潛平臺試驗?zāi)Ｐ?/p>

從圖9可見，ITTC試驗平臺的首搖響應(yīng)試驗值與計算值有較大差異，這是因為理論分析時，浪向是不會改變的，但模型試驗時，由于ITTC試驗平臺的艏搖，實際浪向和波浪力也發(fā)生了改變，試驗結(jié)果必然會與理論分析結(jié)果有所差異。但總體來說，平臺各自由度 RAO通過數(shù)值計算所得的結(jié)果與試驗結(jié)果變化規(guī)律基本相似，吻合程度達(dá)到一定的要求，從而證明了構(gòu)建的三維數(shù)值波浪水池是可靠的，數(shù)值計算所采用的計算方法能夠很好地模擬規(guī)則波作用下平臺的運動狀態(tài)。

圖9 浪向角180°時計算結(jié)果與試驗結(jié)果對比

4 水動力響應(yīng)分析

利用上述方法建立的數(shù)值波浪水池對概念設(shè)計的半潛式浮基風(fēng)電平臺在波浪作用下的動力響應(yīng)進(jìn)行數(shù)值模擬研究，從波高、波浪周期兩方面對浮基平臺進(jìn)行動力響應(yīng)分析。

4.1 模型建立

首先建立計算所需的幾何模型，對結(jié)構(gòu)進(jìn)行網(wǎng)格劃分，均采用四面體網(wǎng)格，除6根連接壓載物的桿件網(wǎng)格尺寸為0.5m外，其余結(jié)構(gòu)網(wǎng)格長度均取為1.0m，輸出的整體結(jié)構(gòu)有限元模型見圖10。其中節(jié)點數(shù)為36516個，單元數(shù)為155681個。

4.2 計算區(qū)域及網(wǎng)格劃分

計算區(qū)域大小取為600m×200m×200m，即是三維波浪水池的3部分尺寸為波浪生成區(qū)、工作區(qū)和消波區(qū)，其長度分別為150m、150m和300m，浮體距入口200m處，水深為200m。

在造波區(qū)和工作區(qū)，x方向網(wǎng)格長度取為2m，y方向網(wǎng)格長度取為4m，Z方向為5m，底邊為20m，消波區(qū)表面網(wǎng)格X向為10m，Z向為10m，其中工作區(qū)中流固耦合界面網(wǎng)格尺寸為0.5m，連桿部分為0.25m。整個流固耦合計算模型網(wǎng)格點數(shù)為109427點，單元數(shù)為624021個，計算區(qū)域網(wǎng)格劃分情況見圖11。

4.3 不同波浪周期情況下浮式基礎(chǔ)橫搖、縱搖運動分析

分別對波高5m，波浪周期分別為4s、6s和10s，對應(yīng)波長分別為25m、56m和156m時的浮基風(fēng)電平臺進(jìn)行動力響應(yīng)分析。橫搖響應(yīng)時程曲線見圖12，縱搖響應(yīng)曲線見圖13。

圖10 風(fēng)電浮基整體結(jié)構(gòu)有限元模型

圖11 流固耦合分析中流體模塊

圖12 不同波浪周期時橫搖響應(yīng)時程曲線

圖13 不同波浪周期時的縱搖響應(yīng)時程曲線

由圖12、13可知，波浪周期為4s時，浮式基礎(chǔ)的縱搖角度很小，在所選時間段內(nèi)，最大縱搖角峰值達(dá)到 0.5°，而橫搖角更小，最大橫搖角達(dá)到 0.1°；波浪周期為 6s時，浮式基礎(chǔ)的縱搖角度比波浪周期為4s時的縱搖角度稍有增加，在所選時間段內(nèi)，最大縱搖角達(dá)到1.6°，而橫搖角表現(xiàn)出相同的增大規(guī)律，最大橫搖角度達(dá)到0.21°；波浪周期為10s時，浮式基礎(chǔ)的縱搖增加幅度相比于4s情況較大些，在所選時間段內(nèi)，最大縱搖角度達(dá)到5.4°，橫搖角變化增加幅度不大，最大橫搖角度達(dá)到0.65°。

4.4 不同波高情況下浮式基礎(chǔ)橫搖、縱搖運動分析

取波浪周期為10s，浪高分別為3m、5m和7m的風(fēng)電浮基進(jìn)行動力響應(yīng)分析，通過數(shù)值模擬得到3種浪高情況下風(fēng)電浮基的橫搖、縱搖響應(yīng)時程曲線，見圖14、15。

圖14 不同波高時的橫搖響應(yīng)時程曲線

圖15 不同波高時的縱搖響應(yīng)時程曲線

從圖 14、15中浮基風(fēng)電平臺橫搖、縱搖響應(yīng)時程曲線中可知，浮基的橫搖、縱搖角度整體上隨波高的增大而增大，并且波高3m時風(fēng)電浮基最大縱搖角度達(dá)到3.8°，波高5m時浮基縱搖角度如上述為5.4°，波高為7m時浮基的縱搖角度達(dá)到7.8°，3m、5m和7m波高情況下浮基橫搖角度峰值分別為0.52°、0.65°和 0.79°。

5 結(jié) 語

利用懸掛壓載物降低平臺重心的方法設(shè)計了一種半潛式浮基風(fēng)電平臺，并基于CFD方法，應(yīng)用ADINA軟件建立了數(shù)值波浪水池，對第20屆ITTC半潛平臺在規(guī)則波作用下的動力響應(yīng)進(jìn)行數(shù)值模擬，模擬結(jié)果與試驗值吻合較好，驗證了所構(gòu)建的三維數(shù)值波浪水池的可靠性。

利用構(gòu)建的數(shù)值波浪水池對新設(shè)計的半潛式浮基風(fēng)電平臺在規(guī)則波作用下的運動進(jìn)行數(shù)值模擬，給出了浮基風(fēng)電平臺在不同波高、波浪周期情況下橫搖、縱搖時程曲線。通過數(shù)值計算結(jié)果分析，得出結(jié)論：在所選海況下，所設(shè)計的半潛式浮基風(fēng)電平臺的縱搖角最大值為7.8°，略超出單個風(fēng)機(jī)系統(tǒng)在正常作業(yè)情況下所要求的系統(tǒng)俯仰角運動范圍（國外著名學(xué)者Thomas Zambrano等人[11]提出在正常作業(yè)情況下，深海風(fēng)電場單個風(fēng)機(jī)系統(tǒng)的平均俯仰角應(yīng)在±6°范圍內(nèi)，而動態(tài)俯仰角在±15°范圍內(nèi)）。本數(shù)值計算結(jié)果將有助于了解浮基風(fēng)電平臺塔架頂端的運動、位移和受力，為海上風(fēng)機(jī)系統(tǒng)的優(yōu)化設(shè)計提供合理的荷載參數(shù)。

[1] 楊雄文，樊洪海. Spar平臺結(jié)構(gòu)型式及總體性能分析[J]. 石油礦場機(jī)械，2008，37(5): 32-35.

[2] 楊雄文，樊洪海. TLP平臺結(jié)構(gòu)型式及總體性能分析[J]. 石油機(jī)械，2008，36(5): 70-73.

[3] 徐萬海，曾曉輝，吳應(yīng)湘，等. 深水張力腿平臺與系泊系統(tǒng)的耦合動力響應(yīng)[J]. 振動與沖擊，2009, 28(2): 145-150.

[4] Zambrano T, Mac Cready T, Kiceniuk T, Jr, et al. Dynamic Modeling of Deepwater Offshore Wind Turbine Structures in Gulf of Mexico Storm Conditions[A]. 25th International Conference on Offshore Mechanics and Arctic Engineering[C], Hamburg,2006. 1-6.

[5] Bye A, Erbrich C, Rognlien B, et al Geotechnical Design of Bucket Foundations[A]. Offshore Technology Conference[C],Houston, 1995.

[6] Rasmussen J L, Feld T, S rensen P H. Bucket Foundation for Offshore Wind Farms Comparison of Simplified Model and FE Calculations[A]. OWEMES 2000[C], Sicilia, 2000.

[7] Feld T. Design Procedures for Bucket Foundations—A New Innovative Foundation Concept Applied to Offshore Wind Turbines[A]. Offshore Wind Energy Conference[C], Brussels, 2001.

[8] Thomas Zambrano, Tyler MacCready and Taras Kiceniuk Jr. et al. Dynamic Modeling of Deepwater Offshore Wind Turbine Structures in Gulf of Mexico storm condition[A]. Proceedings of OMAE 2006 25th International Conference on Offshore Mechanics and Arctic Engineering[C]. Hamburg, ASME, 2006: OMAE2006-92029.

[9] Tang Yougang, Hu Jun, Liu Liqin. Study on the Dynamic Response for Floating Foundation of Offshore Wind Turbine[A]. The ASME 2011 30th International Conference on Ocean, Offshore and Arctic Engineering[C], Rotterdam, 2011, 19-24.

[10] 葉小嶸，張 亮，吳海濤，等. 平臺運動對海上浮式風(fēng)機(jī)的氣動性能影響研究[J]. 華中科技大學(xué)學(xué)報：自然科學(xué)版，2012,40(3): 123-126.

[11] J.M. Jonkman, M.L. Buhl Jr. Development and Verification of a Fully Coupled Simulator for Offshore Wind Turbines[A]. The 45th AIAA Aerospace Science Meeting and Exhibit, Wind Energy Symposium[C]. Nevada, 2007: NREL/CP-500-40979.