點(diǎn)吸收式波浪能俘獲裝置結(jié)構(gòu)優(yōu)化方法研究*

王 鑫 李大鳴 王兵振 韓林生

點(diǎn)吸收式波浪能俘獲裝置結(jié)構(gòu)優(yōu)化方法研究*

王 鑫1, 2李大鳴1①王兵振2韓林生2

(1. 天津大學(xué)水利工程仿真與安全國家重點(diǎn)實(shí)驗(yàn)室 天津 300072; 2. 國家海洋技術(shù)中心 天津 300112)

根據(jù)布設(shè)海域波浪能資源特征, 對(duì)點(diǎn)吸收式波浪能裝置結(jié)構(gòu)進(jìn)行優(yōu)化, 可有效提高裝置的能量俘獲效率。本文以威海褚島北部海域?yàn)檠b置布設(shè)目標(biāo)海域, 在對(duì)該海域波浪資源特征進(jìn)行統(tǒng)計(jì)分析的基礎(chǔ)上, 計(jì)算得到裝置的直徑, 同時(shí)利用數(shù)值軟件計(jì)算出多組工況下點(diǎn)吸收裝置吃水深度和裝置固有周期的對(duì)應(yīng)關(guān)系, 并利用統(tǒng)計(jì)學(xué)方法得出裝置固有周期隨吃水深度的變化規(guī)律, 進(jìn)而分析得出裝置在該海域的最佳設(shè)計(jì)吃水深度, 為點(diǎn)吸收波浪能裝置結(jié)構(gòu)優(yōu)化設(shè)計(jì)提供了新的思路。該方法對(duì)于其他技術(shù)類型波浪能裝置的結(jié)構(gòu)優(yōu)化設(shè)計(jì)同樣具有借鑒意義。

波浪能; 點(diǎn)吸收; 結(jié)構(gòu)優(yōu)化; 吃水深度

點(diǎn)吸收式波浪能裝置利用自由浮動(dòng)的浮子隨波浪上下振蕩俘獲波浪能量, 通過浮子連接的液壓裝置或機(jī)械裝置將俘獲的波浪能轉(zhuǎn)換機(jī)械能, 再通過電機(jī)將機(jī)械能轉(zhuǎn)換為電能。點(diǎn)吸收式波浪能裝置采用相對(duì)簡(jiǎn)便的錨泊式固定, 易于進(jìn)行從裝機(jī)功率幾千瓦到幾百千瓦的系列化設(shè)計(jì), 在解決邊遠(yuǎn)海島、各類海上功能設(shè)施供電問題方面具有廣泛的應(yīng)用前景。

由美國海洋電力技術(shù)公司研發(fā)的波浪發(fā)電浮標(biāo)(Power Buoy)為點(diǎn)吸收式波浪能發(fā)電裝置的典型代表。目前有單機(jī)功率40kW和150kW兩種產(chǎn)品, 已基本實(shí)現(xiàn)了產(chǎn)品化。國內(nèi), 中國科學(xué)院廣州能源研究所研制與航標(biāo)形成一體化結(jié)構(gòu)的10W、60W、100W系列化點(diǎn)吸收波浪能發(fā)電裝置已初步具備產(chǎn)品化條件(王立國等, 2013; 王坤林等, 2017; 盛松偉等, 2017, 2019); 山東大學(xué)劉延俊等(2017)在點(diǎn)吸收裝置的液壓系統(tǒng)、海洋電力技術(shù)系統(tǒng)優(yōu)化等方面取得了豐碩的研究成果。此外, 集美大學(xué)、中國海洋大學(xué)等科研單位也在點(diǎn)吸收式波浪能裝置的研發(fā)方面開展了大量有益的工作(羅續(xù)業(yè)等, 2014; 麻常雷等, 2017)。

有關(guān)點(diǎn)吸收式波浪能發(fā)電裝置能量俘獲裝置的外形, 業(yè)內(nèi)學(xué)者已開展了大量的研究工作, 分析對(duì)比了長方形、球形、水平圓柱形、垂直圓柱形等多種外形結(jié)構(gòu)。研究結(jié)果表明, 垂直圓柱形對(duì)于波浪能的吸收和轉(zhuǎn)換最為有利(劉秋林, 2016)。本文即針對(duì)垂直圓柱形點(diǎn)吸收式波浪能裝置的結(jié)構(gòu)優(yōu)化方法開展研究。

1 目標(biāo)海域波浪資源特征分析

目標(biāo)海域的波浪資源特征是開展裝置外形優(yōu)化設(shè)計(jì)的依據(jù)。本文選擇威海褚島北部海域?yàn)檠b置布設(shè)目標(biāo)海域, 并利用布設(shè)于該海域的浮標(biāo)觀測(cè)數(shù)據(jù)對(duì)該海域的波浪資源特征進(jìn)行分析, 連續(xù)觀測(cè)時(shí)間為2015年1月24日到12月24日, 經(jīng)統(tǒng)計(jì)分析, 出現(xiàn)概率較高的波高和波周期如表1所示。

表1 2015年有效波高-平均周期頻率分布

Tab.1 Effective wave height-average periodic frequency distribution in 2015

可見, 該海域全年以1m以下的小波浪為主, 波周期主要集中在3—5s之間, 其中以3—4s周期波浪最為集中, 低功率的點(diǎn)吸收式波浪能開發(fā)利用技術(shù), 正是適合于該目標(biāo)海域的波浪資源特點(diǎn)的技術(shù)形式。

2 點(diǎn)吸收波浪能裝置結(jié)構(gòu)關(guān)鍵技術(shù)參數(shù)分析

點(diǎn)吸收式波浪能裝置符合弗汝德-克雷洛夫(Froude-Krylov, 簡(jiǎn)稱F-K)條件, 即入射波動(dòng)場(chǎng)的波壓強(qiáng)分布不因潛體的存在而改變【何光宇, 2015 #49;毛來鋒, 2013 #54】, 則作用在潛體上的波浪力一般表達(dá)式如公式(1):

式中,為作用在潛體上的波浪力,F為未受擾動(dòng)的入射波壓強(qiáng)對(duì)浮體的作用力,為反映附加質(zhì)量效應(yīng)的繞射系數(shù), 通常繞射系數(shù)需要通過模型試驗(yàn)給予確定。

F可分解為入射波作用在整個(gè)潛體上的水平作用力和垂直作用力。對(duì)于垂直圓柱俘體, 在某一時(shí)刻浮體承受的波浪力的垂直作用力為F(何光宇等, 2015; 楊紹輝等, 2016a, b; Tagliabue, 2019)【何光宇, 2015 #49;Borgarino, 2012 #39】:

式中,1()為零一階第一類貝塞爾函數(shù), 浮體的半徑為, 浮體高為l, 吃水深度為d, 其中為波浪作用到浮體的時(shí)間, 水深為,C為垂直繞射系數(shù)1.02, ρ為海水密度1.025×103kg/m3, g為重力加速度9.8m/s2,為波高,為波數(shù)。當(dāng)= 0時(shí),F等于波浪力的振幅:

由波浪力垂直分量的計(jì)算公式(2)可見, 在特定海域水深條件下, 裝置浮子的直徑、入水深度兩個(gè)關(guān)鍵設(shè)計(jì)參數(shù)是裝置獲能效率的主要決定因素。

3 確定裝置浮子直徑

波浪能發(fā)電裝置一般包含三個(gè)能量轉(zhuǎn)換部分: 波浪能俘獲系統(tǒng)、機(jī)械傳動(dòng)系統(tǒng)、發(fā)電機(jī)。分析當(dāng)前波浪能開發(fā)利用技術(shù)的研究水平, 波浪能俘獲系統(tǒng)的能量轉(zhuǎn)換效率一般約為40%, 機(jī)械傳動(dòng)系統(tǒng)的工作效率可達(dá)到70%, 發(fā)電機(jī)的合理工作效率為90%, 則當(dāng)前技術(shù)水平下波浪能發(fā)電裝置的整體能量理想的轉(zhuǎn)換效率約為:

以低裝機(jī)功率的裝置為研究對(duì)象, 設(shè)研究目標(biāo)的裝機(jī)功率為3kW。發(fā)電裝置輸入端的波浪能功率in(單位: kW)需要達(dá)到:

利用常用的波浪能量簡(jiǎn)化公式計(jì)算波浪能裝置所處波浪場(chǎng)的波浪功率密度Pd (單位: kW/m)

根據(jù)波浪資源的分析統(tǒng)計(jì)結(jié)果, 褚島海域平均波周期主要分布在3—5s的區(qū)間內(nèi), 以1m以下小波浪為主, 大部分情況下波高不超過2m。因此選擇出現(xiàn)概率較高的波高和波周期為裝置的設(shè)計(jì)輸入, 即波浪能裝置的設(shè)計(jì)波高為1m, 波浪平均周期為4s, 這種波況條件下的波浪能功率密度為:

對(duì)于水平截面為圓型的垂直圓柱點(diǎn)吸收波浪能俘獲裝置來說, 裝機(jī)3kW的波浪俘獲裝置的直徑(單位: m)約為:

在此, 將圓整為6.0m。

4 裝置最優(yōu)吃水深度分析

點(diǎn)吸收波浪能裝置是通過裝置浮子的垂蕩運(yùn)動(dòng)來獲取能量, 即對(duì)于垂蕩運(yùn)動(dòng), 每個(gè)裝置都有一個(gè)固有周期, 當(dāng)波浪的周期與該周期相同時(shí), 裝置的垂蕩運(yùn)動(dòng)幅值達(dá)到最大, 在該條件下, 波浪能裝置的俘獲效率達(dá)到最大。首先要研究得出裝置吃水深度與裝置固有周期的對(duì)應(yīng)關(guān)系; 進(jìn)而在資源分析的基礎(chǔ)上, 確定目標(biāo)海域出現(xiàn)頻率相對(duì)較高的波周期, 以此作為裝置固有周期的設(shè)計(jì)目標(biāo), 并通過固有周期與裝置吃水深度的對(duì)應(yīng)關(guān)系計(jì)算出裝置的最佳吃水深度。本文對(duì)點(diǎn)吸收裝置固有周期與裝置浮子吃水深度的對(duì)應(yīng)關(guān)系分析借助AQWA軟件, 計(jì)算前文分析得出的6m直徑浮子不同吃水深度對(duì)應(yīng)的裝置固有周期, 進(jìn)而分析得出二者的對(duì)應(yīng)變化關(guān)系。

4.1 計(jì)算方法驗(yàn)證

利用AQWA計(jì)算張亮等(2015a)的一個(gè)例子, 將本文采用計(jì)算方法得出的結(jié)果與算例結(jié)果進(jìn)行比較, 以驗(yàn)證本文算法的合理性。

計(jì)算為一種圓柱形點(diǎn)吸式波浪能裝置, 如圖1所示。俘獲裝置的外徑為6m, 內(nèi)孔直徑為2m, 高度為1.6m, 靜止?fàn)顟B(tài)下吃水深度為0.68m。裝置的質(zhì)量為17.517t。

圖1 計(jì)算的點(diǎn)吸收裝置外形尺寸

在利用軟件建立模型的過程中, 將全局坐標(biāo)系的-平面與靜水面重合,軸與波浪傳播的方向平行, 方向與波浪傳播方向相反;軸與鉛錘方向平行, 豎直向上; 坐標(biāo)原點(diǎn)位于水線面的形心處。模型水下深度根據(jù)考察對(duì)象的具體數(shù)值進(jìn)行設(shè)定, 水線以上的高度設(shè)定為1.0m。模型的質(zhì)量按模型的排水量估算, 為保證運(yùn)動(dòng)的穩(wěn)定性, 將模型的重心設(shè)定在浮心下面100mm處。

在此基礎(chǔ)上, 利用ANSYS mesh模塊生成俘獲裝置表面網(wǎng)格單元, 裝置網(wǎng)格劃分情況如圖2所示, 最大單元尺寸設(shè)置為0.2m, 共計(jì)9920個(gè)單元, 網(wǎng)格類型采用程序默認(rèn)設(shè)(De O Falc?o, 2009; Lopes, 2009; Penalba, 2017)。最大允許波頻率為1.14Hz, 計(jì)算過程中采用規(guī)則波, 入射波浪的周期范圍為2.4—5s。初始入水深度0.5m, 俘獲裝置初始靜平衡狀態(tài)處于豎直位置。水線面將浮子模型分為上下兩部分, 只有水下部分參與水動(dòng)力繞射計(jì)算。

圖2 計(jì)算網(wǎng)格

在實(shí)際工作過程中, 點(diǎn)吸收裝置以垂直方向的垂蕩運(yùn)動(dòng)為主, 在導(dǎo)柱的限制下其他方向的往復(fù)振蕩和擺動(dòng)的幅度很小。因此, 在計(jì)算模型中為限制繞三個(gè)坐標(biāo)軸的搖擺運(yùn)動(dòng), 特別是繞軸的搖擺運(yùn)動(dòng), 設(shè)置了繞3個(gè)坐標(biāo)軸搖擺運(yùn)動(dòng)的較大數(shù)值的阻尼。

波浪能俘獲裝置的垂蕩運(yùn)動(dòng)響應(yīng)幅值比的計(jì)算結(jié)果如圖3所示。圖中, 橫軸為入射波浪的周期, 縱軸為裝置垂蕩運(yùn)動(dòng)響應(yīng)幅值比。由圖3可知, 裝置在波周期為3.1s時(shí)的垂蕩響應(yīng)幅值最大, 即裝置的固有周期為3.1s。

圖3 算例裝置振蕩幅值與周期關(guān)系

文獻(xiàn)的計(jì)算結(jié)果如圖4所示, 該計(jì)算結(jié)果同樣表明裝置在波周期為3.1s時(shí)垂蕩運(yùn)動(dòng)響應(yīng)達(dá)到最大, 且對(duì)照垂蕩運(yùn)動(dòng)響應(yīng)隨波浪周期的變化趨勢(shì)與本文計(jì)算的結(jié)果, 兩者近似相同(張亮等, 2015a, b; 鄭雄波, 2016)。這說明本文建立的計(jì)算模型合理可行, 應(yīng)用該方法進(jìn)行圓柱型點(diǎn)吸收式裝置的波浪能轉(zhuǎn)換特性計(jì)算可有效支撐浮子吃水深度與裝置固有周期對(duì)應(yīng)關(guān)系的分析。

圖4 文獻(xiàn)中算例的計(jì)算結(jié)果

注：引自張亮等, 2015a

4.2 分析確定裝置的最優(yōu)吃水深度

考察了初始吃水深度對(duì)圓柱形點(diǎn)吸收式波浪能俘獲裝置固有周期的影響, 計(jì)算過程中裝置的外徑均為6m, 參照國內(nèi)點(diǎn)吸收波浪能裝置的設(shè)計(jì)經(jīng)驗(yàn), 水面以上部分的高度均取為1m (林江波, 2006; 史宏達(dá)等, 2011, 2017)。開展以下9組計(jì)算, 得到不同浮子吃水深度下的裝置固有周期, 計(jì)算結(jié)果分布要能夠覆蓋目標(biāo)海域波周期主要的分布區(qū)間, 即3—5s。計(jì)算工況相關(guān)參數(shù)如表2, 各工況中裝置浮心位置均位于水下部分的中心, 重心位置位于浮心位置正下方, 兩者相距0.1m。計(jì)算結(jié)果如圖5、表3所示。

表2 裝置工況相關(guān)參數(shù)

Tab.2 The parameters of operation condition for the device

圖5 工況1至9周期與響應(yīng)振幅系數(shù)關(guān)系

表3 裝置吃水深度和裝置固有周期對(duì)應(yīng)表

Tab.3 Corresponding values between the draft and natural period of the device

根據(jù)表3, 利用插值法繪制裝置浮子吃水深度和固有周期對(duì)應(yīng)關(guān)系曲線(圖6), 橫軸為吃水深度, 縱軸為固有周期。

圖6 裝置浮子吃水深度和固有周期對(duì)應(yīng)關(guān)系曲線

根據(jù)對(duì)目標(biāo)海域波浪能資源特征的分析, 該海域全年3—5s周期的波浪出現(xiàn)的概率最大, 其中以3—4s周期波浪最為集中, 裝置固有周期設(shè)計(jì)為4s可最大限度的獲得波浪能量, 根據(jù)分析得出的浮子吃水深度和固有周期的對(duì)應(yīng)關(guān)系, 可得出在目標(biāo)海域波浪能資源條件下, 裝置的最優(yōu)吃水深度為2.44m。

至此, 根據(jù)目標(biāo)海域波浪能資源特性, 對(duì)圓柱形浮子點(diǎn)吸收波浪能轉(zhuǎn)換裝置的外形關(guān)鍵技術(shù)參數(shù)的優(yōu)化設(shè)計(jì)完成, 既直徑6m, 固有周期4s, 最優(yōu)吃水深度2.44m。

5 結(jié)論

本文在充分掌握目標(biāo)海域波浪能資源特性的情況下, 研究了對(duì)點(diǎn)吸收波浪能俘獲裝置結(jié)構(gòu)進(jìn)行優(yōu)化的方法, 并重點(diǎn)關(guān)注點(diǎn)吸收裝置浮子直徑和吃水深度的分析和確定方法。通過研究可得出以下結(jié)論: 根據(jù)弗汝德-克雷洛夫定律分析, 在所處水深確定的情況下, 點(diǎn)吸收波浪能裝置浮子的直徑和吃水深度是其關(guān)鍵參數(shù), 是裝置能量俘獲能力的重要決定因素。對(duì)點(diǎn)吸收波浪能裝置結(jié)構(gòu)的優(yōu)化方法研究, 即是對(duì)兩個(gè)參數(shù)在目標(biāo)海域環(huán)境條件下最優(yōu)值的研究確定。通過計(jì)算可以得到浮子吃水深度和浮子固有周期之間的對(duì)應(yīng)關(guān)系; 分析裝置目標(biāo)布置海域的波浪資源特征, 可得到出現(xiàn)概率最高, 也即“蘊(yùn)含波浪能量”最多的波浪周期值, 該值即為點(diǎn)吸收裝置的設(shè)計(jì)目標(biāo)固有周期, 與之對(duì)應(yīng)浮子吃水深度即為裝置的最優(yōu)吃水深度, 最終得出裝置外形結(jié)構(gòu)的優(yōu)化設(shè)計(jì)方案。該方法可應(yīng)用于大多數(shù)圓柱形點(diǎn)吸收裝置的外形結(jié)構(gòu)優(yōu)化設(shè)計(jì)工作。

王立國, 游亞戈, 張亞群等, 2013. 波浪能發(fā)電裝置動(dòng)力攝取系統(tǒng)研究進(jìn)展. 機(jī)床與液壓, 41(1): 165—168, 162

王坤林, 盛松偉, 葉 寅等, 2017. 波浪能裝置中液壓發(fā)電系統(tǒng)Boost變換機(jī)理及控制策略. 電力系統(tǒng)自動(dòng)化, 41(12): 173—178

史宏達(dá), 曲 娜, 曹飛飛等, 2017. 振蕩浮子波能發(fā)電裝置浮子運(yùn)動(dòng)性能的試驗(yàn)研究. 中國海洋大學(xué)學(xué)報(bào), 47(6): 124—130

史宏達(dá), 劉 棟, 劉 臻, 2011. 一種新型波能發(fā)電裝置的運(yùn)動(dòng)響應(yīng)計(jì)算. 中國海洋大學(xué)學(xué)報(bào), 41(10): 111—116

劉延俊, 賀彤彤, 2017. 波浪能利用發(fā)展歷史與關(guān)鍵技術(shù). 海洋技術(shù)學(xué)報(bào), 36(4): 76—81

劉秋林, 2016. 點(diǎn)吸收浮子陣列的波能轉(zhuǎn)換特性研究. 北京: 清華大學(xué)博士學(xué)位論文, 23—67

楊紹輝, 何宏舟, 陳 滬等, 2016a. 陣列筏式波浪能發(fā)電系統(tǒng)設(shè)計(jì)與試驗(yàn)研究. 機(jī)械工程學(xué)報(bào), 52(11): 57—62

楊紹輝, 何宏舟, 李 暉等, 2016b. 點(diǎn)吸收式波浪能發(fā)電技術(shù)的研究現(xiàn)狀與展望. 海洋技術(shù)學(xué)報(bào), 35(3): 8—16

何光宇, 楊紹輝, 何宏舟等, 2015. 陣列式波浪能發(fā)電裝置的水動(dòng)力分析. 水力發(fā)電學(xué)報(bào), 34(2): 118—124

張 亮, 國 威, 王樹齊, 2015a. 一種點(diǎn)吸式波浪能裝置水動(dòng)力性能優(yōu)化. 哈爾濱工業(yè)大學(xué)學(xué)報(bào), 47(7): 117—121

張 亮, 國 威, 鄭雄波, 2015b. 圓柱浮子垂蕩運(yùn)動(dòng)性能的CFD模擬與流場(chǎng)分析. 華中科技大學(xué)學(xué)報(bào)(自然科學(xué)版), 43(12): 65—70

林江波, 2006. 浮子式海浪發(fā)電船的動(dòng)態(tài)分析與仿真. 秦皇島: 燕山大學(xué)碩士學(xué)位論文, 20—35

羅續(xù)業(yè), 夏登文, 2014. 海洋可再生能源開發(fā)利用戰(zhàn)略研究報(bào)告. 北京: 海洋出版社, 72—81

鄭雄波, 2016. 兩類點(diǎn)吸式波能裝置水動(dòng)力特性研究. 哈爾濱: 哈爾濱工程大學(xué)博士學(xué)位論文, 42—56

盛松偉, 王坤林, 吝紅軍等, 2019. 100 kW鷹式波浪能發(fā)電裝置“萬山號(hào)”實(shí)海況試驗(yàn). 太陽能學(xué)報(bào), 40(3): 709—714

盛松偉, 張亞群, 王坤林等, 2017. 波浪能發(fā)電裝置能量轉(zhuǎn)換系統(tǒng)多級(jí)負(fù)載試驗(yàn). 可再生能源, 35(2): 311—316

麻常雷, 夏登文, 王海峰, 2017. 國內(nèi)外海洋能進(jìn)展及前景展望研究. 北京: 海洋出版社, 56—70, 121—132

De O Falc?o A F, 2009. Wave energy utilization: a review of the technologies. Renewable and Sustainable Energy Reviews, 14(3): 899—918

Lopes M F P, Hals J, Gomes R P F, 2009. Experimental and numerical investigation of non-predictive phase-control strategies for a point-absorbing wave energy converter. Ocean Engineering, 36(5): 386—402

Penalba M, Touzón I, Lopez-Mendia J, 2017. A numerical study on the hydrodynamic impact of device slenderness and array size in wave energy farms in realistic wave climates. Ocean Engineering, 142: 224—232

Tagliabue A, Izzo L, Mella M, 2019. Absorbed weak polyelectrolytes: impact of confinement, topology, and chemically specific interactions on ionization, conformation free energy, counterion condensation, and absorption equilibrium. Journal of Polymer Science Part B: Polymer Physics, 57(9): 491—510

THE STRUCTURAL OPTIMIZATION OF POINT ABSORPTION IN WAVE ENERGY CAPTURE DEVICE

WANG Xin1, 2, LI Da-Ming1, WANG Bing-Zhen2, HAN Lin-Sheng2

(1. State Key Laboratory of Hydraulic Engineering Simulation and Safety, Tianjin University, Tianjin 300072, China; 2. National Ocean Technology Center, Tianjin 300112, China)

According to the characteristics of wave energy resources in a sea area, the energy capture efficiency can be effectively improved by optimizing the structure of the point absorption device. Based on the statistical analysis of the characteristics of wave resources in the north area of Chudao Island, Weihai, Shandong, East China, the diameter of the device, the corresponding relationship between the draft of point absorption device, and the natural period of the device under multiple working conditions were calculated using numerical software, and the variation of the natural period of the device with the draft was understood by using a statistical method. The optimal design draft of the device in the sea area was determined, which provides a new idea for the structural optimization design of the point absorbing wave energy device. This method can also be used for reference in the structural optimization design of other types of wave energy devices.

wave energy; point absorption; structural optimization; draft

* 國家重點(diǎn)研發(fā)計(jì)劃海洋環(huán)境安全保障專項(xiàng), 2016YFC1401204號(hào)。王 鑫, 博士研究生, E-mail: wxtjnk022@126.com

李大鳴, 教授, E-mail: lidaming@tju.edu.cn

2019-09-22,

2020-01-19

TM614

10.11693/hyhz20190900179