楊 鋒，高季章，張金接，曾 迪

（1.中國水利水電科學(xué)研究院 工程抗震研究中心，北京 100038；2.北京中水科水電科技開發(fā)有限公司，北京 100038）

1 研究背景

風(fēng)電是可再生能源開發(fā)的主要組成部分，近幾年得到了較快的發(fā)展，我國東部沿海的經(jīng)濟(jì)發(fā)展和電網(wǎng)特點(diǎn)與歐洲類似，適于大規(guī)模發(fā)展海上風(fēng)電，但該區(qū)域海床表層多是淤泥質(zhì)地層（深度＞10m），目前建設(shè)的近海風(fēng)機(jī)基本均采用樁基基礎(chǔ)結(jié)構(gòu)[1]。近年來，我國研究人員對海上風(fēng)機(jī)樁基礎(chǔ)-塔架結(jié)構(gòu)體系的動力特性與環(huán)境荷載下的響應(yīng)分析均開展了一定的研究。嚴(yán)根華等[2]考慮基礎(chǔ)結(jié)構(gòu)的不同沖刷深度，進(jìn)行了土基海流耦合條件下的海上風(fēng)機(jī)塔柱支撐結(jié)構(gòu)動力特性研究；劉志強(qiáng)[3]進(jìn)行了海上風(fēng)機(jī)三樁導(dǎo)管架基礎(chǔ)-塔架在風(fēng)荷載、冰荷載和地震荷載下的動力響應(yīng)分析；陳法波[4]計算分析了風(fēng)、浪和地震等荷載下的基礎(chǔ)結(jié)構(gòu)動力響應(yīng)；陳小波[5]進(jìn)行了基于SPH理論的近海單樁風(fēng)機(jī)結(jié)構(gòu)體系與波浪動力相互作用研究。以上研究均從理論上進(jìn)行了計算分析，但缺乏現(xiàn)場結(jié)構(gòu)體系的振動監(jiān)測驗證，造成目前設(shè)計階段仍重點(diǎn)進(jìn)行基礎(chǔ)結(jié)構(gòu)的動力校核，難以基于動力響應(yīng)分析進(jìn)行優(yōu)化設(shè)計，并據(jù)此提出相應(yīng)的振動控制措施。

為進(jìn)行環(huán)境荷載下風(fēng)機(jī)樁基礎(chǔ)-塔架結(jié)構(gòu)體系振動特性的深入研究，在江蘇響水近海2MW試驗風(fēng)機(jī)的基礎(chǔ)、塔架不同位置安裝多套振動監(jiān)測儀器，進(jìn)行了長期的樁基礎(chǔ)-塔架結(jié)構(gòu)體系振動監(jiān)測。本文依據(jù)實(shí)測相關(guān)數(shù)據(jù)，分析了結(jié)構(gòu)體系振動與風(fēng)、浪等環(huán)境荷載的相關(guān)性，并與動力有限元數(shù)值模擬計算結(jié)果進(jìn)行了對比分析，驗證了數(shù)值模擬分析的可行性。

2 監(jiān)測儀器布置

針對動力響應(yīng)分析工作，為獲取相關(guān)基礎(chǔ)資料，進(jìn)行了風(fēng)、浪和結(jié)構(gòu)振動的實(shí)時監(jiān)測。

2.1 振動監(jiān)測儀器風(fēng)機(jī)基礎(chǔ)-塔架結(jié)構(gòu)體系在波浪、風(fēng)荷載作用下振動相對較小、頻率分布范圍大，因此現(xiàn)場監(jiān)測儀器采用FBAES-T型高精度地震加速度計。

江蘇響水近海試驗風(fēng)機(jī)采用八樁混凝土高樁承臺基礎(chǔ)，根據(jù)基礎(chǔ)-塔架結(jié)構(gòu)體系的特點(diǎn)，現(xiàn)場在鋼管樁基礎(chǔ)、塔架的不同關(guān)鍵位置共埋設(shè)了4支振動監(jiān)測儀，其中在距離混凝土平臺底面0m、2m深的樁基內(nèi)分別安裝了1支，塔架內(nèi)安裝了2支（分別位于塔架底部和頂部），并采用無線傳輸實(shí)時傳輸監(jiān)測數(shù)據(jù)?，F(xiàn)場儀器的布設(shè)如圖1。

2.2 波浪觀測儀器采用“浪龍”儀器進(jìn)行實(shí)時波浪觀測，可觀測波浪的波高、周期及波向等多個波浪參數(shù)，現(xiàn)場儀器的布設(shè)如圖2。

圖1 現(xiàn)場振動監(jiān)測儀器安裝位置

圖2 現(xiàn)場“浪龍”儀器布設(shè)

2.3 測風(fēng)儀器響水近海試驗風(fēng)機(jī)東南1km位置，建有1座70m高的海上測風(fēng)塔，在70m、60m等不同高度安裝有多套測風(fēng)儀器，本文采用現(xiàn)場實(shí)測的實(shí)時風(fēng)資源數(shù)據(jù)。

3 實(shí)測數(shù)據(jù)分析

響水試驗風(fēng)機(jī)振動監(jiān)測時段為2011年3月—2012年8月，同期進(jìn)行了波浪和風(fēng)速現(xiàn)場觀測。觀測時間較長，數(shù)據(jù)量大，本文選擇有代表性的時段進(jìn)行詳細(xì)、重點(diǎn)分析。根據(jù)實(shí)測資料，選擇代表性時間段為2012年8月1—3日（“達(dá)維”臺風(fēng)）和2011年7月26日—8月7日（較大陣風(fēng)）。

3.1 振動監(jiān)測數(shù)據(jù)現(xiàn)場所有振動監(jiān)測數(shù)據(jù)的整理統(tǒng)計結(jié)果如表1，代表性時間段數(shù)據(jù)如圖3。觀測期內(nèi)ZD2、ZD3、ZD4位置的最大振動峰值加速度均出現(xiàn)在2011年7月30日和2012年8月2日，在選擇的代表性時間段；ZD4位置的各方向振動加速度最大，ZD1位置其次，ZD3和ZD2加速度相對最小，各位置Z方向的加速度相差不大。

表1 不同監(jiān)測部位觀測期內(nèi)的最大加速度

3.2 波浪觀測數(shù)據(jù)參考風(fēng)機(jī)附近區(qū)域的實(shí)測完整年海浪的平均波高和最大波高變化統(tǒng)計（如表2）、平均周期和最大周期變化統(tǒng)計（如表3），結(jié)合風(fēng)機(jī)位置實(shí)測波浪資料，2012年8月2—3日“達(dá)維”臺風(fēng)經(jīng)過風(fēng)機(jī)位置時的1h有效波高最大達(dá)3.22m，最大峰值波高達(dá)5.38m。

3.3 實(shí)測風(fēng)數(shù)據(jù)代表性時間段的風(fēng)速變化如圖4所示，風(fēng)速在一段時間能變化較大，具有一定的特殊性（風(fēng)速大、均方差大），分別是較大的陣風(fēng)和臺風(fēng)時刻，為海上風(fēng)機(jī)必須考慮的工況。

3.4 振動響應(yīng)與環(huán)境動荷載相關(guān)性2011年7月26日—8月7日代表性連續(xù)時段內(nèi)最大振動加速度出現(xiàn)在陣風(fēng)較大時刻，此時均方差為3.4，最大振動加速度（Gxy）為0.256g，最大波浪浪高為2.88m；風(fēng)速、均方差與塔架頂部（ZD4位置）振動加速度存在一定的相關(guān)性，均方差較大時，振動加速度相對較大。風(fēng)機(jī)振動、波浪數(shù)據(jù)和現(xiàn)場實(shí)測的風(fēng)速、均方差（標(biāo)準(zhǔn)差）等參數(shù)進(jìn)行統(tǒng)計分析結(jié)果圖5。

圖3 代表性時間段的實(shí)測振動數(shù)據(jù)

表2 現(xiàn)場實(shí)測海浪波高變化

表3 現(xiàn)場實(shí)測海浪周期變化

圖4 不同代表性時段的相關(guān)時刻70m高度實(shí)測風(fēng)速

2011年3月—2012年8月整體觀測時段內(nèi)不同方差、風(fēng)速對應(yīng)的振動加速度平均值如圖6。該圖顯示，實(shí)測風(fēng)速的均方差與塔架頂部振動加速度存在明顯的相關(guān)關(guān)系，隨著均方差的增大，風(fēng)機(jī)塔架頂部的振動加速度也相應(yīng)增加。

圖5 峰值加速度與風(fēng)速均方差、浪高相關(guān)性

圖6 觀測時段內(nèi)均方差、風(fēng)速平均值與振動加速度的關(guān)系

4 數(shù)值模擬

根據(jù)實(shí)測數(shù)據(jù)海上風(fēng)機(jī)基礎(chǔ)-塔架結(jié)構(gòu)體系的相關(guān)分析，獲得了一些基本規(guī)律，但都是定性的。結(jié)合江蘇響水試驗風(fēng)機(jī)的相關(guān)參數(shù)，本文建立了動力響應(yīng)分析計算模型，模擬了不同工況的實(shí)測風(fēng)、浪脈動時程曲線，進(jìn)行基礎(chǔ)-塔架結(jié)構(gòu)體系動力響應(yīng)計算分析，將計算結(jié)果與實(shí)測結(jié)果對比分析，討論建立動力模型計算分析的可靠性。

4.1 風(fēng)、浪荷載模擬采用Kaimal風(fēng)速功率譜模型，順風(fēng)向脈動風(fēng)荷載直接利用現(xiàn)場實(shí)測的10min風(fēng)速均方差（標(biāo)準(zhǔn)差）、10min平均風(fēng)速進(jìn)行模擬計算，并采用實(shí)測的10min瞬時最大風(fēng)速、最小風(fēng)速（采樣頻率：0.5Hz）進(jìn)行修正。Kaimal功率譜模型的計算公式為[6]：

式中：L、V、σu分別為湍流長度、測量高度的平均風(fēng)速和風(fēng)速標(biāo)準(zhǔn)差。

采用改進(jìn)BROTSCHNEIDER譜（B-M波浪譜），隨機(jī)波浪直接利用實(shí)測的有效波高、有效周期進(jìn)行模擬計算，并用實(shí)測的最大波高進(jìn)行校正。B-M波浪譜模型的計算公式為[7]：

式中：h、ω0分別為實(shí)測有效波高和有效周期。

表4 試驗風(fēng)機(jī)附近位置“達(dá)維”臺風(fēng)的實(shí)測風(fēng)速

表5 試驗風(fēng)機(jī)位置“達(dá)維”臺風(fēng)時刻的實(shí)測波浪

圖7 2012年8月2日“達(dá)維”臺風(fēng)經(jīng)過時瞬時風(fēng)荷載和隨機(jī)波浪力的模擬時程曲線

模擬計算采用Matlab編程實(shí)現(xiàn)，2012年8月2日的“達(dá)維”臺風(fēng)經(jīng)過試驗風(fēng)機(jī)位置時的實(shí)測風(fēng)浪數(shù)據(jù)如表4—5所示，脈動時程模擬計算結(jié)果如圖7。

4.2 動力有限元模型建立本文考慮土-樁-上部結(jié)構(gòu)相互作用時采用“m法”，建立模型時，具體假設(shè)樁周圍的土作為一系列非耦合的彈簧連接在樁和遠(yuǎn)場土之間；彈簧的材料特性可以是線彈性的，也可以是非線性的，彈簧串聯(lián)或并聯(lián)施加在樁上；樁基模擬成一系列的離散單元，來模擬樁-土共同作用。以響水2MW試驗風(fēng)機(jī)八樁混凝土高樁承臺基礎(chǔ)為例，建立動力響應(yīng)計算分析模型，風(fēng)機(jī)的參數(shù)如表6所示[8]。

表6 響水2MW試驗風(fēng)機(jī)的參數(shù)

模型的具體參數(shù)為：鋼管樁直徑為1.40m（壁厚22mm），樁長約為66.5m，底端高程為-64.00m，頂端高程1.80m；8根樁在承臺底部沿10m直徑的圓周均勻分布，斜度為7∶1，樁頂埋入承臺深度2.0m；承臺直徑14.00m，厚度3.0～4.5m，頂高程為4.30m。

4.3 風(fēng)浪聯(lián)合作用響應(yīng)分析選擇不同代表性時刻的實(shí)測風(fēng)和波浪數(shù)據(jù)，分別進(jìn)行大風(fēng)浪、中風(fēng)浪、小風(fēng)浪及臺風(fēng)風(fēng)浪等不同工況的結(jié)構(gòu)響應(yīng)分析。具體選擇的工況如表7所示。

表7 風(fēng)浪聯(lián)合作用工況的選取

圖8 波浪工況4作用下基礎(chǔ)-塔架結(jié)構(gòu)不同位置的加速度時程

風(fēng)浪工況4作用下，試驗風(fēng)機(jī)現(xiàn)場儀器不同布置位置的動力響應(yīng)如下圖8所示，各種不同工況的計算結(jié)果與現(xiàn)場實(shí)測如表8所示，ZD2、ZD3、ZD4三個位置的有限元計算值與實(shí)測值相關(guān)性分析如圖9所示。

表8 不同風(fēng)浪工況下的計算值與監(jiān)測值 （單位：m/s）

圖9 ZD2、ZD3、ZD4三個位置的有限元計算值與現(xiàn)場監(jiān)測值相關(guān)性分析

由以上計算結(jié)果與實(shí)測值對比分析可見得出如下結(jié)論。

（1）ZD1位置儀器實(shí)測結(jié)果較計算分析結(jié)果明顯偏大，實(shí)測值為計算值的5～10倍，與該儀器距離2m的ZD2實(shí)測值也相差較大，且實(shí)測結(jié)果均為這一規(guī)律，主要原因在于：ZD1儀器安裝位置為現(xiàn)場波面的變化區(qū)域，波浪荷載直接作用在該位置樁基上，對樁基有明顯的拍擊作用，而本文計算時按脈動荷載輸入隨機(jī)波浪荷載，未考慮波浪拍擊樁基造成的樁基礎(chǔ)振動；同時，距離ZD1位置2m的ZD2觀測結(jié)果明顯減小，且與計算值接近，也說明ZD1的觀測結(jié)果主要為波浪拍擊樁基儀器位置產(chǎn)生的局部振動。

（2）風(fēng)浪荷載聯(lián)合下，上部的ZD2、ZD3、ZD4三個監(jiān)測點(diǎn)現(xiàn)場實(shí)測值與計算動力響應(yīng)相差不大（同一量級），并具有較好的相關(guān)性，相關(guān)系數(shù)大于0.9，動力有限元模擬計算結(jié)果較好反映了響水試驗風(fēng)機(jī)結(jié)構(gòu)體系不同工況的動力響應(yīng)。

（3）計算結(jié)果與現(xiàn)場實(shí)測結(jié)果存在一定的差異，主要原因在于：計算分析時未考慮風(fēng)機(jī)運(yùn)行時風(fēng)機(jī)機(jī)組自身振動的影響，可能造成部分高頻振動缺失；計算時風(fēng)荷載直接按照水平推力系數(shù)加載，風(fēng)荷載變化與振動計算結(jié)果對應(yīng)性好，但實(shí)際風(fēng)機(jī)運(yùn)行過程中，偏航、啟動等過程相對風(fēng)速瞬時變化有延遲效應(yīng)，可能造成計算結(jié)果比實(shí)測結(jié)果偏大或偏小。

5 結(jié)論

（1）結(jié)構(gòu)體系的實(shí)測振動加速度與實(shí)測風(fēng)速的均方差存在明顯的相關(guān)關(guān)系，隨著均方差的增大，風(fēng)機(jī)塔架頂部的振動加速度也相應(yīng)增加；（2）通過現(xiàn)場實(shí)測風(fēng)、波浪荷載的時程曲線模擬，動力有限元計算結(jié)果較好反映了響水試驗風(fēng)機(jī)不同工況的動力響應(yīng)，可用于該區(qū)域海上風(fēng)機(jī)結(jié)構(gòu)體系的驗證分析與優(yōu)化設(shè)計，也可為類似工程提供參考；（3）波浪荷載對基礎(chǔ)-塔架結(jié)構(gòu)體系的振動也有一定的影響，特別是樁基位置，但本文采用波浪的加載方式與實(shí)際波浪對樁基的作用存在一定的差異，有待進(jìn)一步開展波浪荷載的加載方法研究。

[1]楊鋒，邢占清，等.近海風(fēng)機(jī)基礎(chǔ)結(jié)構(gòu)型式研究[J].水利水電技術(shù)，2009，40（9）：35-38.

[2]嚴(yán)根華，古華，等.土基流固耦合條件下海上風(fēng)電場塔架支撐結(jié)構(gòu)動力特性初探[J].中國工程科學(xué)，2010，12（11）：49-52.

[3]劉志強(qiáng).海上風(fēng)力發(fā)電支撐體系在環(huán)境荷載作用下動力響應(yīng)分析[D].大連：大連理工大學(xué)，2009.

[4]陳法波.海上風(fēng)機(jī)結(jié)構(gòu)動力反應(yīng)分析[D].大連：大連理工大學(xué)，2010.

[5]陳小波.近海風(fēng)機(jī)結(jié)構(gòu)體系環(huán)境荷載及動力響應(yīng)研究[D].大連：大連理工大學(xué)，2011.

[6]王耀南，孫春順，等.用實(shí)測風(fēng)速校正的短期風(fēng)速仿真研究[J].中國電機(jī)工程學(xué)報，2008，28（11）：94-100.

[7]竺艷蓉，謝峻，等.各種波浪譜在海洋工程中適用性的研究[J].海洋學(xué)報，1995，17（6）：126-131.

[8]周建華，張金接，邢占清.超靜孔隙水壓力影響的海上風(fēng)機(jī)樁基承載力高應(yīng)變檢測方法研究[J].中國水利水電科學(xué)研究院學(xué)報，2012，10（1）：17-22.