可再生能源利用技术的发展可以克服化石能源带来的风险, 比如全球变暖、环境污染等, 作为化石燃料的替代资源, 波浪能在可再生能源中能量密度最高(Leijon et al, 2003)。我国当下发展中, 大力开发和利用海洋能资源对于我国实行可持续发展战略, 建设生态文明, 缓解常规能源紧张带来的能源问题具有重大战略意义(管轶, 2011)。

20世纪70年代, Tornkvist(1975)和Hulls(1977)就已经利用有限的大洋船报资料和浮标资料对全球海洋沿岸波浪能资源分布情况展开分析, 结果显示北大西洋、太平洋东北部、澳大利亚南部沿岸以及南非的西南部沿岸等海域为全球波浪能富集区域, 中国海域属于波浪能贫瘠区域。船报数据虽然比较准确, 但是其有空间和时间的限制, 一些极端海况下波浪数据会有遗漏(Gulev et al, 2003)。随着海洋探测技术的发展和大规模海洋调查的开展, 特别是海浪数值模式和卫星遥感的精确, 越来越多的卫星反演资料(Barstow et al, 1998)以及数值同化的再分析资料(Folley et al, 2009; Iglesias et al, 2010)等普遍被应用于波浪能的评估之中。

中国从20世纪80年代开展波浪能资源评估工作, 早期仅对部分海域展开分析(任建莉等, 2009), 全国性波浪能评估工作始于国家海洋局第二海洋研究所(王传崑等, 2009)。这些成果推进了我国波浪能评估工作的发展, 其后基于实测、数值模拟(郑崇伟等, 2011, 2012, 2013;宗芳伊等, 2014)以及再分析数据(万勇等, 2014;石洪源等, 2017)的波浪能评估工作相继开展。已有研究指出我国波浪能资源丰富地区为东海和南海, 考虑到波浪能资源的开发条件, 虽然东海和南海波能资源丰富, 但是该海域受台风影响较多(You et al, 2018), 波浪能装置在极端大浪和大风条件下生存力差, 遭受破坏风险高。因此, 我国波浪能的实际开发利用可以放置在渤海或者黄海某些条件较好的区域, 但是针对这两个海区的波能评估相对较少, 同时波浪能评估多采用深水评估公式, 但是黄渤海大部分水深相对较浅, 利用深水公式评估过渡区、浅水区波浪能资源会造成评估结果的不准。鉴于上述原因, 本文利用适用于各种水深条件的波浪能评估公式以及欧洲中期天气预报中心(European Centre for Medium-Range Weather Forecasting, ECMWF)提供的最新波浪数据集ECMWF Reanalysis v5(ERA5), 展开黄渤海波浪能资源评估工作。

1 数据资料 1.1 数据来源

ECMWF是全球几家最主要的再分析数据中心之一。ERA5是ECMWF第五代再分析数据集, 相比ERA-Interim(Dee et al, 2011), ERA5同化系统使用了若干专为再分析而开发的综合预报系统(IFS cycle 41r2) (Hersbach et al, 2016), 数据准确性得到了显著改善(Meng et al, 2018; Czernecki et al, 2019)。ERA5数据具有更高分辨率:从地面到0.01hPa(约80km)分为137个垂向σ层, 水平分辨率为31km(约0.25°)。ERA5提供的变量从ERA-Interim的100个增加到240个, 包括耦合波模型提供的波高和波向, 使用户能够更准确地分析过去的大气和海洋状态。许多研究人员将其用于大气和海洋研究, 如Bechtle等(2019)采用ERA5的风和波浪数据集, 进行风能和波能资源评估。Mahmoodi等(2019)基于ERA5再分析数据集对波斯湾波浪能量的时空特征展开分析。

1.2 插值方法

ERA5提供的数据时间分辨率为1h, 空间分辨率为0.25°, 在空间分辨率不高的情形下, 为了获得较高分辨率本文采用Kriging插值法处理数据集。Kriging插值法可将数学函数与指定数量的点或指定半径内的所有点进行拟合以确定每个位置的输出值(Burrough, 1986)。

2 评价方法 2.1 波浪能评估公式

波浪能资源的时空分布状况可以通过波浪能的特征量和评价指标进行评估。基于波浪谱的能量通量可以表述为式(1):

式中, P_w, ρ, g, f, h, θ, C_g(f, h)和S(f, θ)分别表示波能流密度、海水密度、重力加速度、频率、水深、波向、群速以及波谱密度。

群速可以写成:

这里, T是波浪周期, k是波数(k = 2π/L), L为波长。

对于自由波而言, 其波能流密度可以写成有效波高H_s和能量周期T_e的函数:

通过式(1)和(2), 可以得到

其中, 波长L的计算公式为: 。

因此, 式(4)可以写成:

该公式适用于各种水深条件, 对于深水区, 和tanh(kh)趋近于1, 因此波能流密度可以写成:

式(6)为波能评估中常用的公式, 本文做出如下证明, 利用式(5)除以公式(6), 得到一个新的函数:

其中, λ表示各种水深条件下的波能流密度对于深水区的波能流密度的比值。

以kh为未知数, 参考You(2008)给出的深浅水定义(kh＜0.1π为浅水, kh＞π为深水, 介于二者之间为过渡区), 得到的λ变化见图 1, 可以看出:浅水区域, 因为波浪的破碎, 利用深水公式计算会明显高估其波浪资源; 过渡区内, 使用深水公式部分海域是高估波浪能资源, 部分区域是低估波浪能资源, 最大低估量约为20%;深水区内, 也存在极小一部分低估区域。之前的研究中, 仅仅指出利用深水公式会高估浅水区域的波浪能资源(Mahmoodi et al, 2019), 却没有指出在过渡区内对波浪能资源的低估, 因此在波浪能评估中利用深水公式是不准确的。

考虑到水深对波浪能评估的影响, 本文利用Beji(2013)提出的新的色散方程μ_B和将计算结果精度进行提高的公式μ_*(Simarro et al, 2013), 把μ_B和μ_*代入公式(5)中得出新的波浪能评估公式(8):

其中, 而 在式中: 。

2.2 评价指标

本文的评价指标借鉴现有研究成果(郑崇伟等, 2012;万勇等, 2014;石洪源等, 2017), 选用可用度和富集度频率以及变异系数进行定量分析黄渤海域的波浪能资源。

在评定波浪能资源的过程中, 将大于2kW/m出现的频率视为可用度频率, 大于20kW/m出现的频率视为富集度频率, 计算黄渤海1980年1月-2018年12月逐1h的波能流密度, 对该海域波能流密度大于2kW/m和大于20kW/m出现的概率进行统计。

3 波浪能评估结果 3.1 黄渤海波浪要素的年特征

对中国近海的有效波高(H_s)、平均周期()、波能流密度(P_w)进行多年平均, 得出研究海域的波浪要素的分布规律(图)。从图 2中可以看出, 研究海域年均波高最大值为1.3m, 位于黄海东部地区。就海区而言, 渤海区域最小, 年均值不超过0.8m。年平均周期以及波能流密度的分布和有效波高H_s的分布类似, 黄海东部最大。从波能分布来看, 和前期研究结果类似, 黄渤海区域波能资源并不丰富。

3.2 波浪能季节分布特征

对波浪而言, 由于其与风速关系密切相关, 我国不同的季节风速差异巨大, 因此在研究中应当给出不同季节的波能分布情况。本研究分别统计了冬季(12-2月)、春季(3-5月)、夏季(6-8月)和秋季(9-11月)四个季节的波能流密度(P_w)见图 3。

由图 3可知, 我国周边海域的波能流密度分布具有季节性, 其中秋冬季节, 年均波能流密度较大; 春夏期间年均波能流密度较小。总的来看, 渤海四季变化明显, 整体仍低于黄海, 春夏两季波能流密度不足2kW/m, 秋冬最大值可超过4kW/m。黄海区域特征四季变化不大, 数值变化特征与渤海类似。在春季, 波能流密度的大值主要集中在:黄海南部附近海域(4-5kW/m); 夏季波能流密度大值区域分布与春季类似, 进入秋季, 波能流密度的大值区与春夏两季分布相同, 约6-8kW/m; 长江口附近波能流密度相对较大, 在10kW/m。对黄渤海而言, 远海的波能流密度大于近岸, 黄海东南部为波能流密度较大海区, 具有较大的开发潜力; 渤海中部及渤海海峡附近海域秋冬季节波能流密度相较其他季节而言比较丰富, 具有开发价值。在冬季, 黄海地区波能流密度大部分区域较大, 数值变化呈现出东-西走向。

3.3 新公式优越性分析

渤海是我国水深最浅的海域, 平均水深只有18m, 在本文中, 将渤海作为公式验证的区域。在渤海中, 分别利用公式(6)和公式(8)计算39年波能流密度, 将二者进行比值分析, 结果见图 4。在水深较浅的区域, 深水公式出现了明显的高估; 尤其在黄河入海口处, 最大比值接近140%, 在深水与浅水的过渡区, 比值在90%-100%之间; 在深水中, 深水公式与浅水公式比值基本为100%, 说明在深水区域, 二者计算得到的波能流密度基本相等。

3.4 能级频率计算结果

39年数据的能级频率统计结果见图 5。黄海大部分海域可用度频率在30%以上, 有些区域可用度频率达60%, 富集度频率大约在2%-5%之间, 主要集中在黄海东部海域; 渤海地区可用度频率最小, 最大不超过23%, 富集度出现的频率不超过3%。

黄海大部分海域波浪四季期间可用度频率都较高(图 6, 7)。春季, 黄海东南部可用度频率在50%以上, 其余海域较低, 渤海最小, 不足20%;富集度频率黄渤海域都较低, 最大值不足5%。夏季, 渤海整体不足10%, 黄海大部分海域可用度频率超过30%, 最大可达50%以上; 富集度频率有所增大, 但幅度不大, 最大富集度频率为5%;秋季, 可用度频率增大, 最大值仍出现在黄海东南部海域, 接近70%;富集度频率有了较大提高, 最大值可达8%, 出现在黄海东南部。冬季, 黄渤海域内可用度频率整体跃升, 渤海大部分海域在30%以上, 黄海大部分海域集中在50%以上, 渤海的富集度频率相较秋季而言有所下降, 最大值依旧不及5%, 黄海富集度频率上升明显。

3.5 变异系数计算结果

将中国海域1980-2018年共39年的波能流密度按季节和年度分别进行了变异系数的计算, 计算结果见图 8和图 9。按照季节来看, 春季, 黄海中部变异系数较小; 渤海中南部、黄海东南部海域变异系数较大; 夏季, 波能流密度变异性有增大的趋势, 波能最不稳定的季节; 秋季, 波能不稳定性延续夏季的趋势, 但是不稳定性减弱; 冬季, 波能趋于稳定, 其中黄海中部能稳定性最好(图 8)。

全年统计结果显示, 我国波能变异系数呈条带状分布, 总体上从南至北呈现增大的趋势, 但数值变化不大, 高值区最大不超过2.3, 低值区不低于0.5, 总体趋于稳定。沿海地区39年波能变异系数的低值区出现在黄海南部近长江口海域, 高值区位于渤海南部、黄海北部和渤海沿海地区(图 9)。

3.6 可利用波高结果

对于海洋能开发而言, 并不是所有波高都可以利用, 太小波浪带不动发电装置运转, 太大的波浪会造成发电装置的损坏, 例如台风浪, 因此本文选择有效波高介于1-4m之间的波浪开展研究(图 10)。春季, 渤海可用波浪出现频率最大值集中在渤海中部, 不超过20%。黄海分布呈现东南-西北走向, 大部分区域可用波浪出现频率超过20%, 最大不超过45%。夏季, 渤海可用波浪频率较低, 整体海域不超过10%;黄海大部分海域集中在20%以上, 最大值在40%左右。秋季, 可用波浪分布频率有所提高, 渤海还最大值分布在中部, 超过30%;黄海大部分海域超过30%, 最大值可超过50%。冬季, 可用波浪分布概率最大, 渤海大部分海域超过30%;黄海大部分海域超过50%, 最大处可达60%以上。黄海东南部海域四季都有开发价值, 渤海冬季具有开发价值。

3.7 重点开发利用区

根据刘首华等(2015)的定义, 波能流密度介于0-1kW/m属于波浪能贫乏区, 1-1.5kW/m属于可利用区, 1.5kW/m以上属于波浪能丰富区。结合上述分析, 本研究将成山头以东、渤海海峡外侧以及长江口外海海域作为重点开发利用区(图 11黑框处), 与前人分析相吻合(万勇等, 2014;刘首华等, 2015)。该区域的波能流密度年均可达3kW/m以上, 可用度频率超过30%, 富集度频率大于2%, 长江口外海甚至可达5%, 变异系数小于1.8, 说明稳定性好, 便于采集与利用; 可利用波高频率大多数时间在30%以上, 冬季甚至可达50%以上, 具有较大的开发价值, 且渤海海域部分受台风影响小, 发电装置不易受到破坏。在长江口外海海域的波浪能集中在5-9kW/m, 与万勇等(2015)计算结果基本一致, 刘首华等(2015)计算出成山头附近海域波能流密度10年平均近岸区域在1.5kW/m以上, 与本文的1.3kW/m相差不大, 冬季计算数值在3.5kW/m与本文的3kW/m略有差异, 说明旧公式对数值计算存在高估, 新公式在提高准确度方面有所进步。

4 结论

本文利用ECMWF最新的ERA5数据, 并利用新的波浪能评估公式, 对黄渤海海域1980-2018年间波浪能资源展开评估, 主要计算指标包括波浪能可用度频率、富集度频率、变异系数以及可利用波高占比等, 结果显示:

(1) 黄渤海区波浪能资源具有明显的季节性, 秋冬季节较高, 春夏季节较低, 冬季是波浪能资源开发的最佳季节。

(2) 波浪能富集区域主要集中在渤海海峡外侧、成山头东部以及长江口外海区域。

(3) 沿海地区39a波能变异系数的低值区出现在长江入海口北部, 高值区位于渤海南部、黄海北部和渤海沿海地区。

(4) 对黄渤海而言, 长江口外海区域、渤海海峡以及成山头附近海域波浪能资源最具有开发潜力。

综上所述, 本文利用新公式对原有波浪能计算方法进行改进, 对黄渤海海域的波浪能进行评估分析, 在近岸地区提高了精度, 使计算结果更加精确, 为今后在该海域的资源开发利用提供参考。