芦潮港海上风资源变化特性分析

doi:10.12096/j.2096-4528.pgt.22010

芦潮港海上风资源变化特性分析

张智伟¹, 张建平², 刘明², 纪海鹏³, 诸浩君¹, 周圣荻⁴

1.上海绿色环保能源有限公司, 上海市崇明区 200433

2.上海理工大学机械工程学院, 上海市杨浦区 200093

3.上海电力大学能源与机械工程学院, 上海市浦东新区 201306

4.上海海湾新能风力发电有限公司, 上海市崇明区 202156

Analysis on Variation Characteristics of Offshore Wind Resources in Luchao Port

ZHANG Zhiwei¹, ZHANG Jianping², LIU Ming², JI Haipeng³, ZHU Haojun¹, ZHOU Shengdi⁴

1.Shanghai Green Environmental Protection Energy Co. , Ltd. , Chongming District, Shanghai 200433, China

2.School of Mechanical Engineering, University of Shanghai for Science and Technology, Yangpu District, Shanghai 200093, China

3.College of Energy and Mechanical Engineering, Shanghai University of Electric Power, Pudong New Area, Shanghai 201306, China

4.Shanghai Gulf New Energy Wind Power Generation Co. , Ltd. , Chongming District, Shanghai 202156, China

收稿日期: 2022-01-19

基金资助:

国家自然科学基金项目. 12172228. 11572187
上海市自然科学基金项目. 22ZR1444400

Received: 2022-01-19

作者简介 About authors

张智伟(1978)，男，工程师，主要从事海上发电技术研究，dianjing@21cn.com；

张建平(1972)，男，博士，教授，主要从事能源安全低碳利用及智能化的研究，本文通信作者，jpzhanglzu@163.com；

刘明(1997)，男，硕士研究生，主要从事海上发电技术研究，mliuusst@163.com；

纪海鹏(1994)，男，硕士，工程师，主要从事海上发电技术研究，jihaipeng20171120@163.com；

诸浩君(1984)，男，工程师，主要从事海上发电技术研究，juleon1011@163.com；

周圣荻(1989)，女，工程师，主要从事海上发电技术研究，414639946@qq.com。

摘要

海上风资源特性对海上风电场的建设具有重大影响，是海上风机选型及排布的关键因素，因此对其进行研究十分必要。利用测风塔测量得到上海芦潮港地区风速数据，并分析了该地区风速、风向、风能密度分布及其变化特性。结果表明：在芦潮港地区海上风速分布描述上，Logistic函数具有最好的精度，Weibull函数很难描述单月风速分布；芦潮港地区风向主要来自正北与东北方向；随着芦潮港地区测风塔测风高度的增加，年平均海上风速逐步趋于8 m/s；芦潮港地区年平均风能密度不低于290 W/m²，具有较为充足的风能资源，可为风电场收益提供保障。

关键词： 海上风电场 ; 风资源特性 ; 海上风机 ; 风速 ; 风向 ; 风能密度

Abstract

The offshore wind resource character has a significant impact on the construction of offshore wind farms, it is a key factor in the selection and arrangement of offshore wind turbines. Therefore, it is necessary to study it. The wind speed data of Shanghai Luchao Port were measured and obtained by the anemometer tower, and the distribution and the variation characteristics of wind speed, wind direction and wind energy density in this area were given and analyzed. The results show that the Logistic function has the best accuracy in the description of offshore wind speed distribution in Luchao Port area, and the Weibull function is difficult to describe the monthly wind speed distribution. The wind direction in Luchao Port area mainly comes from due north and northeast. With the increase of anemometer tower height in Luchao Port area, the annual average offshore wind speed gradually tends to 8 m/s. The annual average wind energy density of Luchao Port area is no less than 290 W/m², which has sufficient wind energy resources to ensure the income of wind farm.

Keywords： offshore wind farm ; wind resource character ; offshore wind turbine ; wind speed ; wind direction ; wind energy density

PDF (4989KB) 元数据多维度评价相关文章导出 EndNote| Ris| Bibtex 收藏本文

本文引用格式

张智伟, 张建平, 刘明, 纪海鹏, 诸浩君, 周圣荻. 芦潮港海上风资源变化特性分析. 发电技术[J], 2022, 43(2): 260-267 DOI:10.12096/j.2096-4528.pgt.22010

ZHANG Zhiwei, ZHANG Jianping, LIU Ming, JI Haipeng, ZHU Haojun, ZHOU Shengdi. Analysis on Variation Characteristics of Offshore Wind Resources in Luchao Port. Power Generation Technology[J], 2022, 43(2): 260-267 DOI:10.12096/j.2096-4528.pgt.22010

0 引言

海上风资源特性对风电场的建设影响巨大，尤其是在海上风机选型及排布上具有决定性的意义^[1-3]。同时，由于风电场开发正逐步走向海上，因此对海上风资源特性开展研究十分必要^[4-6]。

近年来，研究者采用多种分布函数来描述风分布特征。杨正瓴等^[7]分析了季风特性对风电功率预测效果的影响，结果发现，改进空间相关性预测方法既适应季风区，又可用于我国北部非季风区冬季的风电功率预测。周齐等^[8]提出基于概率加权矩法(probability weighted moment method，PWMM)的三参数Weibull分布垂直外推方法，该方法可有效地体现平坦地形低速区的风速分布特征。周泽人等^[9]将混合Weibull分布与Weibull、Normal、Lognormal分布进行了对比，结果表明，混合Weibull分布更适用于较复杂的风速分布状况。Wang等^[10]提出了2种高斯混合分布模型，通过研究发现该模型可用于复杂地区的风速分布估算。姜海燕^[11]对比了Weibull、Rayleigh、Gamma、Lognormal分布在描述低风速场中的性能，结果发现，Weibull分布更接近实际风速频率分布。马契^[12]通过研究发现，两参数Weibull分布在概率密度较低的低风速区精度不高。Yu等^[13]将Nakagami与Rician分布用于评估渤海湾风资源，研究结果表明，渤海湾风速在12 m/s以下，主要在4~8 m/s，4、10月份的主要风速范围高于8、12月份，且夜间风速普遍高于白天。但陈楠等^[14]研究表明，在风资源评估过程中，陆上风资源评估方法在海上风资源评估应用时需要进行加检验。因此，在海上风资源评估研究较少的情况下，亟须对比分析不同分布函数描述海上风速分布的特点。

目前，学者们针对风速变化特性开展了一定的研究。Lin等^[15]通过模拟海风变化的特征发现，在地面以上10 m高处的风速被高估近1 m/s，而夜间近地层垂直风切变被低估。吴强等^[16]针对沿海丘陵地区的风速变化特性进行研究，结果发现，该地区的风速变化具有日周期性且受海陆风影响较大。封宇等^[17]应用机器学习算法，结合Monin-Obukhov相似理论对风速变化特性进行了研究，结果发现，海上风速沿垂直高度方向变化较小，切变系数较小。Sing等^[18]通过对某岛屿的风数据进行统计分析发现，风切变系数与温度之间相关性较好。Jeong等^[19]研究发现，HeMOSU-1海上风场风速的日变化、季节变化以及湍流强度均较小。靳双龙等^[20]通过研究不同区域风切变的特征及差异发现，下垫面平缓地区的风切变系数较小，而地形起伏大的地区风切变系数较大，且海滨地区的垂直风切变值明显大于内陆。

本文基于测风塔测量得到的上海芦潮港地区风速数据，分析了该地区不同高度的风速、风向变化以及风能密度分布特性，揭示了芦潮港海上风资源的变化特性，以期为芦潮港地区后期的海上风电场建设提供参考。

1 风速特性的描述

1.1 描述风速分布的函数

1）Weibull分布函数

两参数Weibull分布概率密度函数^[21]可表示为

f_{1} (v) = \frac{k}{c} {(\frac{v}{c})}^{k - 1} e x p [- {(\frac{v}{c})}^{k}]

(1)

式中： $v > 0$ ，为风速； $k$ 和 $c$ 分别为威布尔分布函数的形状参数和尺度参数。

2）Rayleigh分布函数

Rayleigh分布的概率密度函数^[11]可表示为

f_{2} (v) = \frac{v}{c_{1}^{2}} e x p (- \frac{v^{2}}{2 c_{1}^{2}})

(2)

式中 $c_{1}$ 为Rayleigh分布函数的尺度参数。

3）Gamma分布函数

Gamma分布的概率密度函数^[11]可表示为

f_{3} (v) = \frac{v^{k_{2} - 1}}{c_{2}^{k_{2}} Γ (k_{2})} e x p (- \frac{v}{c_{2}})

(3)

式中 $c_{2}$ 和 $k_{2}$ 分别为Gamma分布函数的尺度参数与形状参数。

4）Normal分布函数

Normal分布的概率密度函数^[22]可表示为

f_{4} (v) = \frac{1}{\tilde{σ} \sqrt[]{2 π}} e x p [- \frac{{(v - \tilde{μ})}^{2}}{2 {\tilde{σ}}^{2}}]

(4)

式中 $\tilde{μ}$ 和 $\tilde{σ}$ 分别为均值与标准方差。

5）Logistic分布函数

Logistic分布的概率密度函数^[22]可表示为

f_{5} (v) = \frac{e x p [- (v - μ) / σ]}{σ {1 + e x p [- {(v - μ) / σ]}}^{2}}

(5)

式中 $μ$ 和 $σ$ 分别为Logistic分布函数的位置参数与形状参数。

6）Log-Logistic分布函数

根据式(5)可得到Log-Logistic分布的概率密度函数^[22]：

f_{6} (v) = \frac{e x p (- μ / σ) v^{(1 / σ - 1)}}{σ {1 + [e x p {(- μ) v]}^{1 / σ}}^{2}}

(6)

1.2 切变指数

风切变指数 $γ$ 可表示为

γ = \frac{l g (v_{2} / v_{1})}{l g (h_{2} / h_{1})}

(7)

式中 $v_{1}$ 、 $v_{2}$ 分别为高度在 $h_{1}$ 、 $h_{2}$ 的风速。

1.3 风能密度

风能密度 $D_{w i n d}$ 计算公式如下：

D_{w i n d} = \frac{1}{2} \sum_{j = 1}^{q} ρ v_{j}^{3} t_{j}

(8)

式中： $q$ 为风速区间内风速数目； $ρ$ 为空气密度； $v_{j}$ 为第 $j$ 个风速区间的风速值； $t_{j}$ 为第 $j$ 个风速区间内风速发生时间。

2 风分布描述及分析

2.1 不同分布函数的风速分布

通过在上海芦潮港地区建立的海上测风塔测量得到了该地区在2004年10月到2005年9月的测风数据，测量高度分别在10、50、60、70 m，测量参数有风速、风向、气压以及气温，其中70 m高风速序列如图1所示。

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1 芦潮港海上风速序列

Fig. 1 Offshore wind speed series in Luchao port

采用6种分布函数对芦潮港2004年10月70 m高风速分布进行描述，结果如图2所示。可以看出：在单一月份情况下，与Gamma和Rayleigh分布函数相比，Normal、Log-Logistic以及Logistic分布能较好地描述风速概率密度分布情况；与其他分布函数相比，两参数Weibull分布虽然可以对单月风速概率密度分布进行描述，但是精确度并不高；Normal分布要差于Logistic分布。

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2 不同分布函数对2004年10月风速分布的描述

Fig. 2 Description of wind speed distribution in October 2004 by different distribution functions

2005年70 m高度上全年风速概率密度分布如图3所示，可以看出：在3~5 m/s区间内风速概率密度明显提高，说明该风速区间内的风速在全年占比相对较高；Logistic分布能够较好地描述芦潮港2005年风速分布情况，而Normal分布则略差。而Weibull分布由于描述效果更差，无法在该图中给出。

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3 2005年风速分布的描述

Fig. 3 Description of wind speed distribution in 2005

综合以上分析结果可知，虽然Weibull函数可用来描述单月风速概率密度，但不能很好地描述芦潮港2005年全年的风速概率密度分布，之所以产生这种情况，是因为各地风速分布情况均不同，尤其是陆地与海上风速分布的差异更为明显。由于目前获取的深远海域风速分布数据较少，因此，是否能用Weibull分布来描述深远海域风速分布有待研究。

2.2 不同风向风速分布

为了显示不同风向上风速分布情况，以2005年的风速分布为例，将不同风向上风速分布绘于图4，由于0°与60°这2个方向上风速概率密度分布具有代表性，因此本文未列出其他方向上风速概率密度分布。由图4可以看出：2005年各月风速分布主要呈现单、双峰的分布特性；低风速在全年各风向分布的占比较高。

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4 2005年不同风向风速分布

Fig. 4 Wind speed distribution of different wind directions in 2005

2.3 风向玫瑰图

为了清晰地观察到风向分布情况，将芦潮港地区2005年全年风向玫瑰图与2004年10月风向玫瑰图分别绘于图5。可以看出：芦潮港地区2005年全年风向主要来自正北方向，占所有风向比例为61.5%；芦潮港地区2004年10月份风向主要来自东北方向，45°~75°范围内的风速在该月占比最大，比例为22.5%。

图5

新窗口打开| 下载原图ZIP| 生成PPT

图5 风向分布玫瑰图

Fig. 5 Rose chart of wind direction distribution

3 风参数变化特性分析

3.1 风变化趋势

图6为指数函数拟合风速与高度变化之间的关系曲线。可以看出，随着高度的不断增加，风速先快速增大，在0~20 m高度内风速变化最为剧烈，当高度增至80 m处时，风速增速逐步变缓，最终风速接近8 m/s。

图6

新窗口打开| 下载原图ZIP| 生成PPT

图6 高度上的风速变化拟合曲线

Fig. 6 Fitting curve of wind speed change at height

根据图6分析可得出：在80 m高度以内，风速变化较为剧烈，不适宜安装大型风机；在80 m高度以上，风速随高度变化则较为平缓，风轮旋转域内风速变化不大。因此，芦潮港地区高度在80 m以上空域更有利于风机平稳发电。

为了探究不同高度之间风速变化相关性，将10、50、60、70 m高度之间风速相关性曲线分别绘于图7(a)—(f)。图7(a)—(f)中拟合直线的决定系数分别为0.807 1、0.766 4、0.745 5、0.991 2、0.977 4与0.992 6，可以看出：10 m高度风速主要分布在0~10 m/s，而其他高度风速则主要分布在0~15 m/s；从决定系数方面来看，高度相差越大，风速相关性越低。

图7

新窗口打开| 下载原图ZIP| 生成PPT

图7 不同高度之间风速相关性

Fig. 7 Wind speed correlation between different heights

此外，由图7(c)可知，10 m与70 m这2个高度上大部分风速呈线性相关。为了得到风向变化相关性，将这2个高度风向相关性曲线绘于图8。可以看出，很多情况下这2个高度风向相关性较差，这表明，在10 m与70 m高度内风向变化没有明显的规律。

图8

新窗口打开| 下载原图ZIP| 生成PPT

图8 10 m与70 m高度之间风向相关性

Fig. 8 Wind direction correlation between 10 m and 70 m height

3.2 不同高度间切变指数的变化

图9(a)—(f)给出了芦潮港地区2005年不同高度间切变指数变化情况。由图9(a)—(c)可以看出：在1—4月与11—12月，切变指数变化幅度与频率较高，且没有明显的规律，最高与最低切变指数分别可以达到2与-2，这是由于此时间段正处于春季与冬季的缘故；而其他月份切变指数随时间变化则较为平缓，说明夏季与秋季风速较稳定。

图9

新窗口打开| 下载原图ZIP| 生成PPT

图9 2005年不同高度之间切变指数变化

Fig. 9 Wind shear index changes between different heights in 2005

由图9(d)—(f)可知：在高度较为相近的情况下，风速切变指数随时间变化波动幅度相比于高度相差较多的情况下平缓许多，但可以观察到在2 000~2 500 h与4 000~5 000 h，也即3—4月与6—7月，切变指数变化仍有部分较大波动，特别是在4 700 h左右时，图9(e)、(f)中切变指数分别接近于-20与-40，从整个2005年来看，这种情况出现频率不高，一般是由台风引起的。

4 风能密度分布特性分析

为了直观显示各风速区间对应的风能密度情况，给出了芦潮港地区70 m高2005年全年与2004年10月份的风速与风能概率密度分布情况，如图10所示。

图10

新窗口打开| 下载原图ZIP| 生成PPT

图10 风速和风能密度分布

Fig. 10 Distributions of wind speed and wind energy density

由图10可见：2005年，在芦潮港地区5 m/s风速出现频率最高，3~14 m/s有效风速出现的频率为92.74%，对应的风能密度频率为99.98%，年平均风速为6.75 m/s，年平均风能密度为298.93 W/m²；2004年10月份芦潮港地区的平均风速为6.58 m/s，对于风速12 m/s与13 m/s，风能密度在该月中占比较大，平均风能密度为290.90 W/m²，风能密度较好，但风向分布不稳定。因此，风电场收益是否能够达到预期更取决于风机的排布。

5 结论

利用测风塔测量得到上海芦潮港地区风速数据，并分析了该地区风速、风向、风能密度分布及其变化特性，得到如下主要结论：

1）无论是全年还是单月情况下，Logistic分布在描述芦潮港地区风速分布上均具有较好的性能，同时在芦潮港地区不同风向上，风速分布具有单、双峰的特性，风向主要来自正北与东北方向；

2）芦潮港地区全年平均风速变化与高度之间存在指数关系，全年平均风速最高接近8 m/s，不同高度之间3、4、6、7月份切变指数变化较为剧烈；

3）芦潮港地区2005年与2004年10月份的平均风速均高于6.5 m/s，风能密度较好，全年平均风能密度不低于290 W/m²，有效风能在90%以上，风能较为丰富。

参考文献

原文顺序

文献年度倒序

文中引用次数倒序

被引期刊影响因子

[1]

， BAHAREH

， PENG

Z L

．

Assessment of long-term offshore wind energy potential in the south and southeast coasts of China based on a 55-year dataset

[J]．Energy，2021，224：120225． doi:10.1016/j.energy.2021.120225