田縱橫 劉海波 喻 飛 劉 爽 甘 進(jìn)*

(武漢理工大學(xué)交通學(xué)院1) 武漢 430063) (長江勘測規(guī)劃設(shè)計研究院2) 武漢 430010)

0 引 言

水面浮式光伏電站，將太陽能發(fā)電技術(shù)與水庫、湖泊等水面結(jié)合起來，建成水面上的光伏電站，與陸地光伏電站相比，具有節(jié)約用地、發(fā)電效率高和充分利用水面資源實(shí)現(xiàn)“漁光互補(bǔ)”等多方面優(yōu)勢，但是由于其工作環(huán)境的特殊性，對浮體結(jié)構(gòu)的承載能力提出了更高的要求.

風(fēng)浪是湖泊中重要的水動力因素，同時也是水面光伏浮體結(jié)構(gòu)安全性評估需要考慮的主要環(huán)境載荷，由于內(nèi)陸湖泊及小型水庫的風(fēng)浪環(huán)境遠(yuǎn)不如海洋環(huán)境惡劣，利用現(xiàn)有的海浪譜分析方法進(jìn)行水面光伏浮體波浪載荷短、長期預(yù)報，顯然存在不足，尋找到適用于水面光伏浮體的湖泊及小型水庫風(fēng)浪載荷至關(guān)重要.目前，國內(nèi)外對湖泊的風(fēng)浪形成和發(fā)展規(guī)律早已展開研究,計算風(fēng)浪要素時大多數(shù)采用半理論、半經(jīng)驗(yàn)法[1]，主要有美國的SMB法、太湖(淺水湖泊)風(fēng)浪計算公式、小風(fēng)區(qū)風(fēng)浪公式、福建莆田實(shí)驗(yàn)站法等.

水面光伏電站是近年來才迅猛發(fā)展的新能源項目，日本的水上漂浮式光伏發(fā)電技術(shù)保持著世界領(lǐng)先地位，韓國和印度也在積極研發(fā)和推廣該技術(shù)[2].國內(nèi)相關(guān)研究還處于起步階段，潘霄等[3]針對單個光伏浮式基礎(chǔ)結(jié)構(gòu)進(jìn)行了靜態(tài)結(jié)構(gòu)強(qiáng)度分析，何嘯[4]提出了一種海洋漂浮式光伏發(fā)電裝置，研究了海浪影響下的浮體運(yùn)動特性及光伏發(fā)電性能.由于水面光伏浮體陣列結(jié)構(gòu)龐大，作業(yè)環(huán)境及載荷復(fù)雜，尚未有水面光伏浮體陣列結(jié)構(gòu)安全性評估方法的研究，也無相關(guān)行業(yè)規(guī)范，因此，需要借鑒文獻(xiàn)[5-7]相關(guān)海上浮式平臺研究方法和成果展開研究.

光伏浮體結(jié)構(gòu)安全性評估主要包括三個方面：①基于三維勢流理論的頻域內(nèi)運(yùn)動響應(yīng)及波浪載荷的分析；②錨固狀態(tài)下鋼絲繩系泊力計算和可靠性分析；③風(fēng)浪聯(lián)合作用工況下的浮體結(jié)構(gòu)強(qiáng)度計算和優(yōu)化.基于此，段艷麗等[8]針對系泊狀態(tài)下的半潛式平臺波浪載荷展開研究，通過數(shù)值分析軟件ANSYS得到了平臺的波浪力和波浪彎矩.薄景富[9]介紹了設(shè)計波法在半潛式平臺整體響應(yīng)和結(jié)構(gòu)強(qiáng)度計算中的應(yīng)用.肖越等[10]利用數(shù)值方法研究了錨泊浮體在風(fēng)、浪、流聯(lián)合作用下的運(yùn)動響應(yīng)和錨鏈線張力變化規(guī)律.曹劍鋒等[11]采用ANSYS/AQWA軟件對極大型浮體水動力響應(yīng)、波浪載荷及總體強(qiáng)度進(jìn)行了計算分析.因此，本文擬結(jié)合相關(guān)風(fēng)浪載荷規(guī)范和船舶與海洋工程領(lǐng)域結(jié)構(gòu)強(qiáng)度直接計算方法，對水面光伏浮體的運(yùn)動響應(yīng)及結(jié)構(gòu)安全展開分析.

1 分析模型及載荷選取

1.1 數(shù)值模型

按照某光伏電站項目的水面光伏浮體建立數(shù)值分析模型，浮體結(jié)構(gòu)主要分為主浮體、連接浮體和光伏組件三部分，見圖1.主浮體和連接浮體均采用高聚乙烯材料，浮體間的連接方式采用螺栓鉸接.

圖1 水面光伏浮體結(jié)構(gòu)圖

1.2 風(fēng)載荷

光伏結(jié)構(gòu)基本風(fēng)壓的計算應(yīng)參照文獻(xiàn)[5]并結(jié)合光伏浮體設(shè)計壽命進(jìn)行選取.以安徽省六安市附近某小型水庫水面光伏電站設(shè)計為例，水面光伏浮體設(shè)計壽命為25年，根據(jù)規(guī)范要求，需采用10年一遇和25年一遇基本風(fēng)壓進(jìn)行風(fēng)載強(qiáng)度校核，全國各城市重現(xiàn)期基本風(fēng)壓計算見式(1)，經(jīng)基本風(fēng)壓公式(2)折算，需校核風(fēng)速分別為18和21 m/s.

(1)

(2)

式中：x100為100年一遇基本風(fēng)壓，取0.4 kN/m2;x10為10年一遇基本風(fēng)壓, 取0.2 kN/m2;R為重現(xiàn)期；w0為基本風(fēng)壓；ρ為空氣密度；v0為需計算校核風(fēng)速.

1.3 波浪載荷

水面光伏浮體工作環(huán)境為內(nèi)陸湖泊或小型水庫的水灣或近岸水域，風(fēng)是產(chǎn)生波浪運(yùn)動的最主要原因.由于廣泛應(yīng)用于海洋結(jié)構(gòu)領(lǐng)域的設(shè)計波法，未考慮到風(fēng)載荷對湖面波浪的直接影響，且計算所需要的湖泊或水庫的實(shí)測風(fēng)浪數(shù)據(jù)難以獲得，故單一使用設(shè)計波法進(jìn)行波浪要素的確定并不合理，因此，本文采用內(nèi)陸湖泊或水庫風(fēng)浪計算公式與設(shè)計波法相結(jié)合的方式，對浮體結(jié)構(gòu)在不同風(fēng)浪條件下的危險波浪工況進(jìn)行計算分析.

本文計算的水面光伏電站位于丘陵環(huán)繞的有限風(fēng)區(qū)水域，有效風(fēng)區(qū)長度為1 000 m.考慮到波浪的波長與風(fēng)區(qū)長度和水域大小密切相關(guān)，有限的風(fēng)區(qū)長度和水域不利于大波長的產(chǎn)生，故參考文獻(xiàn)[5]中莆田-鶴地試驗(yàn)站公式，對校核風(fēng)速下內(nèi)陸湖泊或水庫的平均波高和波長進(jìn)行相關(guān)計算，其平均波高和波長計算公式分別為

(3)

Lm=0.038 6vD0.5g-0.5

(4)

式中:Hm為平均波高；d為平均水深；v為校核風(fēng)速；D為風(fēng)區(qū)長度；Lm為平均波長.綜上可得，10年一遇和25年一遇基本風(fēng)速下的平均波高分別為0.27和0.32 m，平均波長分別為7和8.2 m.

2 水動力計算分析

2.1 頻域運(yùn)動響應(yīng)及危險波工況

根據(jù)內(nèi)陸湖泊和小型水庫波浪特點(diǎn)，水面光伏浮體的破壞形式多為小范圍內(nèi)風(fēng)浪載荷作用下的局部區(qū)域破壞，當(dāng)浮體陣列組合尺寸接近或等于波長時，將可能發(fā)生最嚴(yán)重中拱或中垂現(xiàn)象.結(jié)合10年一遇和25年一遇風(fēng)速下的波長，選取4×4光伏浮體組合進(jìn)行研究，4×4模型長8 m、寬7.5 m，水動力模型見圖2.

圖2 4×4組合水動力模型圖

采用ANSYS/AQWA數(shù)值分析軟件對4×4光伏浮體組合模型進(jìn)行浮體運(yùn)動和波浪彎矩響應(yīng)分析時，應(yīng)首先開展各種響應(yīng)傳遞函數(shù)的預(yù)報，以掌握浮體結(jié)構(gòu)在波浪中響應(yīng)的基本特征.

在光伏浮體運(yùn)動響應(yīng)數(shù)值計算分析時，按照DNV規(guī)范中設(shè)計波法要求，設(shè)置七個波浪入射方向：0°，15°，30°，45°，60°，75°，90°.每個浪向和單位波高下，波浪頻率取 23個，包括了內(nèi)陸湖泊和水庫的主要能量成分和浮體陣列的主要響應(yīng)區(qū)域：0.1～8 rad/s.4×4光伏浮體組合模型在15 m水深條件下，不同波浪入射角的六自由度運(yùn)動響應(yīng)傳遞函數(shù)，見圖3.

圖3 六自由度運(yùn)動響應(yīng)傳遞函數(shù)

考慮到垂向波浪彎矩、水平波浪彎矩及風(fēng)載荷，是影響光伏浮體陣列組合結(jié)構(gòu)安全性的主要因素，風(fēng)浪聯(lián)合作用時可能會產(chǎn)生最危險中垂工況，因此，針對4×4光伏浮體組合模型的垂向波浪彎矩和水平波浪彎矩進(jìn)行詳細(xì)分析，以確定最危險工況時的波浪頻率，見圖4.

圖4 波浪彎矩RAO

由圖4可知，波浪入射角度為0°和90°時，沿X方向和Y方向的垂向波浪彎矩最大，波浪頻率均為0.63 Hz，波浪入射角度為60°時，水平波浪彎矩達(dá)到最大，波浪頻率為1.12 Hz.

結(jié)合文獻(xiàn)[5]要求的風(fēng)載荷、內(nèi)陸湖泊或水庫風(fēng)浪計算公式確定的平均波高和波長，以及頻域內(nèi)計算分析所得到的波浪頻率和相位，可以得到用以錨泊分析及結(jié)構(gòu)強(qiáng)度計算的六種危險波工況，見表1.

表1 4×4光伏浮體組合模型危險波工況

2.2 時域錨泊強(qiáng)度分析

針對100 kW光伏浮體陣列實(shí)際尺寸建立錨泊數(shù)值分析模型，模型長28 m、寬24.12 m，由14×14個光伏浮體單元組成.應(yīng)用AQWA-DRIFT開展風(fēng)浪聯(lián)合作用下的時域計算，計算工況為表1所確定的六種危險波工況，錨固系統(tǒng)為18根鋼絲繩對稱布置的形式，具體布置方案及風(fēng)浪入射角見圖5.

圖5 100 kW光伏陣列錨泊布置圖

在水面光伏浮體錨泊方案設(shè)計的初級階段，運(yùn)用時域耦合的數(shù)值分析方法對錨泊系統(tǒng)進(jìn)行動力分析，研究在4×4光伏浮體陣列發(fā)生局部破壞的六種工況時，設(shè)計的100kW光伏浮體陣列錨泊系統(tǒng)是否穩(wěn)定，鋼絲繩的應(yīng)力是否在許用應(yīng)力范圍之內(nèi)以及能否滿足特定的系泊要求.下面給出時域耦合分析得到的各工況下鋼絲繩最大平均張力的統(tǒng)計數(shù)值，見表2.

表2 不同工況下鋼絲繩張力計算結(jié)果

由表2可知，在工況4時鋼絲繩張力達(dá)到最大值16.8 kN.圖6為該工況下鋼絲繩張力的時域響應(yīng)曲線圖.通過光伏浮體對角拉伸試驗(yàn)，測得單個耳板對接強(qiáng)度大致為10～15 kN，實(shí)際錨泊時采用單根錨繩連接多個耳板的連接方式，每個耳板實(shí)際承載力遠(yuǎn)小于16.8 kN.故而，該錨固方案滿足水面光伏浮體結(jié)構(gòu)強(qiáng)度要求.

圖6 工況4鋼絲繩張力時域響應(yīng)曲線圖

3 浮體陣列結(jié)構(gòu)強(qiáng)度校核

對于大尺度水面光伏浮體陣列來說，由于內(nèi)陸湖泊或小型水庫波浪載荷的實(shí)際作用區(qū)域有限，風(fēng)浪載荷對其破壞位置一般為小范圍內(nèi)的局部區(qū)域，見圖7.為提高設(shè)計階段光伏浮體陣列結(jié)構(gòu)強(qiáng)度的計算效率，本文針對100 kW光伏浮體陣列最危險的三個典型區(qū)域進(jìn)行強(qiáng)度校核, 波浪載荷是將水動力求解的一階波浪速度勢代入伯努利方程求得壓力分布后,通過AQWA-WAVE模塊以壓力載荷的形式加載到有限元模型濕表面上，風(fēng)載荷按照文獻(xiàn)[5]的要求施加在光伏板壁面上.

圖7 典型危險區(qū)域圖

結(jié)構(gòu)分析時使用的典型區(qū)域模型為4×4光伏浮體結(jié)構(gòu)模型，見圖8.根據(jù)典型區(qū)域的位置特點(diǎn)，分別對三種區(qū)域采用不同邊界條件：①號區(qū)域限制A點(diǎn)耳板和B、C端的x，y，z方向位移；②號區(qū)域限制D、E、F端的x，y，z位移；③號區(qū)域四周剛性固定.采用ANSYS軟件進(jìn)行結(jié)構(gòu)強(qiáng)度計算時，利用beam connection近似替代浮體與浮體之間的螺栓連接，有利于節(jié)省大量模型處理和計算分析的時間.

圖8 4×4光伏浮體結(jié)構(gòu)計算模型圖

針對4×4光伏浮體組合結(jié)構(gòu)，分別對三個典型區(qū)域展開六個危險波工況下的結(jié)構(gòu)強(qiáng)度計算分析，風(fēng)速方向與波浪入射方向保持一致.不同工況下的三個典型區(qū)域應(yīng)力與位移計算結(jié)果見表3～5.

表3 ①號區(qū)域不同工況下浮體強(qiáng)度計算結(jié)果

表4 ②號區(qū)域不同工況下浮體強(qiáng)度計算結(jié)果

表5 ③號區(qū)域不同工況下浮體強(qiáng)度計算結(jié)果

對比各工況下不同區(qū)域的應(yīng)力大小，可以確定①號區(qū)域?yàn)槿N區(qū)域中最危險區(qū)域，這也是其處于100 kW光伏浮體陣列的邊角位置決定的.由于0°方向風(fēng)載荷在光伏浮體上的受風(fēng)投影面積較小，故波浪載荷與風(fēng)浪聯(lián)合載荷作用下的應(yīng)力、位移分布相差不大，以波浪載荷作用為主.工況3和工況4的風(fēng)浪入射角為60°，主要考慮水平波浪彎矩與風(fēng)載荷聯(lián)合作用對浮體陣列結(jié)構(gòu)的影響，從表中計算結(jié)果分析，①號區(qū)域和②號區(qū)域風(fēng)載荷影響較大，③號區(qū)域由于處于陣列內(nèi)部，結(jié)構(gòu)較強(qiáng)，主要以局部動壓力產(chǎn)生的應(yīng)力集中為主.工況5和工況6 是垂向波浪彎矩和風(fēng)載荷聯(lián)合作用產(chǎn)生的最危險中垂工況，由于90°方向風(fēng)載荷與各排光伏板受風(fēng)投影面積達(dá)到最大，考慮風(fēng)載荷作用后，三個區(qū)域位置處的光伏浮體陣列結(jié)構(gòu)應(yīng)力和變形顯著增加，其中，①號區(qū)域工況6的應(yīng)力和位移達(dá)到所有工況的最大值，分別為16.03 MPa和0.085 m，見圖9.

圖9 ①號區(qū)域工況6應(yīng)力和位移圖

在結(jié)構(gòu)強(qiáng)度計算分析的基礎(chǔ)上，本文探索了水面光伏浮體陣列在不同工況下的應(yīng)力集中位置和結(jié)構(gòu)破壞形式規(guī)律，發(fā)現(xiàn)處于迎浪邊角位置的光伏浮體陣列區(qū)域結(jié)構(gòu)最弱，陣列內(nèi)發(fā)生應(yīng)力集中的位置主要為耳板與主浮體連接處、連接浮體與主浮體搭接處和主浮體底板角隅處.圖10為①號和②號區(qū)域工況6載荷下連接浮體與主浮體搭接處局部應(yīng)力圖.

圖10 ①號和②號區(qū)域工況6應(yīng)力和位移圖

考慮到①號區(qū)域在工況6的風(fēng)浪載荷下，光伏浮體陣列應(yīng)力和位移變形較大，不利于光伏浮體正常作業(yè)，故針對邊角區(qū)域光伏浮體實(shí)施輔助加強(qiáng)措施，使用長角鋼和螺栓將浮體陣列左側(cè)的耳板串聯(lián)，以及將下側(cè)的主浮體支撐平臺串聯(lián)，在減小浮體陣列的垂向位移變形的同時，將系泊力通過鋼架結(jié)構(gòu)均勻的傳遞到多個浮體結(jié)構(gòu)上，有利于顯著提高浮體陣列結(jié)構(gòu)的安全性.表6為①號區(qū)域改良前后應(yīng)力和位移的對比情況，圖11為改良后浮體陣列的應(yīng)力和位移圖.

表6 ①號區(qū)域改良設(shè)計效果對比

圖11 ①號區(qū)域改良方案應(yīng)力和位移圖

4 浮體陣列結(jié)構(gòu)安全性評估流程

本文針對大尺度水面光伏浮體陣列結(jié)構(gòu)的特殊性，結(jié)合內(nèi)陸湖泊和水庫風(fēng)浪經(jīng)驗(yàn)計算公式和設(shè)計波法，提出一種在缺乏內(nèi)陸湖泊或水庫風(fēng)浪資料的情況下，適用于設(shè)計階段的水面光伏浮體結(jié)構(gòu)錨固強(qiáng)度及結(jié)構(gòu)安全性評估流程，其路線圖見圖12.

圖12 水面光伏浮體結(jié)構(gòu)安全性評估方法流程圖

5 結(jié) 論

1) 100 kW水面光伏浮體陣列錨固方案及結(jié)構(gòu)強(qiáng)度滿足設(shè)計要求，建議進(jìn)一步加強(qiáng)耳板結(jié)構(gòu)，提高耳板承載能力.

2) 由于內(nèi)陸湖泊和水庫波浪載荷的特殊性，大尺度水面光伏浮體陣列破壞位置多為迎浪處的局部邊角區(qū)域，內(nèi)部浮體區(qū)域較為安全.

3) 水面光伏浮體陣列內(nèi)應(yīng)力集中位置多發(fā)生在耳板與主浮體連接處、連接浮體與主浮體搭接處和主浮體底板角隅處.