干货!数据分析报告之风资源分析(2)

本文是《风资源分析》系列的第二篇,深入探讨风能玫瑰图、bin统计图的绘制及数据分析。通过风能玫瑰组图分析风能的月份变化,使用add_subplot和windrose库结合实现。bin图用于风速段内变量统计,通过划分风速区间、计算均值并作图。同时,介绍bin数据导出和按月份对比湍流的方法,以及日分布散点图的时间戳处理技巧。
摘要由CSDN通过智能技术生成

这篇文章紧接着上一篇《干货!数据分析报告之风资源分析(2)》,由于分析的角度和涉及的内容太多了,因此分为1和2两篇。
在上一篇中,主要介绍了绘图前期的数据预处理,可视化展示的直方图,威布尔分布图,散点图和玫瑰图,接下来继续进行分析。

风能玫瑰组图

在平时的应用中,多是采用单个的玫瑰图来进行分析,在本次的分析中,我用了风能玫瑰组图,按照不同的月份来进行玫瑰图区分,看风能的月份变化情况。

import matplotlib.cm as cm
import numpy as np

label=['E','','60','','30','','N','',
       '330','','300','','W','','240','',
       '210','','S','','150','','120','']  #设置玫瑰图标签

#筛选出风速风向数据进行分析
dk=datr.loc[:,('Time','WS_main_Mean','WD_main_Mean')]  #提取出风速风向数据
dk=dk.set_index('Time')  #将时间设置成索引,稍后按月份进行画图
dk.columns=['ws','wd']

month=['2019-7','2019-8','2019-9','2019-10','2019-11','2019-12',
       '2020-3','2020-4','2020-5','2020-6','2020-7']
plt.rcParams['figure.dpi']=100
fig = plt.figure(figsize=(15,9),edgecolor='grey',frameon=False)
fig.subplots_adjust(left=0.03,bottom=0.04,right=0.97,top=0.95,wspace=0.17,hspace=0.3)
for i in range(0,len(month)):
    ax = fig.add_subplot(3,4,i+1, projection='windrose') #添加组图
    wd=dk[month[i]]['wd']
    ws=dk[month[i]]['ws']
    bins=np.arange(0,max(ws),2)  #设置风速间隔为2m/s
    ax.contourf(wd,ws,bins=bins,cmap=cm.hot_r,nsector=360)    
    ax.set_xticks(np.linspace(0,2*np.pi,24,endpoint=False))  #设置刻度个数,与label对应
    ax.set_xticklabels(label,fontsize=12)
    ax.set_yticklabels('')
    ax.grid(linestyle='-.',color='black',lw=0.15,zorder=0)
    ax.set_title(month[i]+&
自己开发的风资源分析工具包WindAnalysis-WindAnalysis风数据分析工具包教程-V1.4.pdf 本帖最后由 He_Challen 于 2017-9-6 14:40 编辑 由于工作的原因,今年项目开始转型风电项目,在慢慢上手的过程中发现,风电所涉及的软件清一色北欧的,好不好用只有用了才知道。因为仅是为前期风电开发做技术分析,老外的软件一个是不容易上手,二是操作复杂。随下决心自己开发一套专门用于项目前期的风资源分析工具包。就这样开始而一发不可收拾,从最开始的结构搭建、输出设计便沉迷此中两个月,推出的前三个版本都不太稳定,要么是兼容不好,要么是数据处理的时逻辑顺序有问题,总之在最初的三个版本在大量项目的测风数据的测试下暴漏出一堆又一堆的BUG。说实话,中途曾想过放弃,一个人孤军奋战实在是太孤独难耐了,多年工作环境造就的内心还是比较强大的,最终还是坚持了下来。在飞机上、动车上、出差的酒店里、办公桌前开始了一遍又一遍的调试修改,度过了一个又一个难免的夜晚。最终完成的兼容性和稳定性都可靠的V.1.4.1版本,经反复测试没有问题后,将这个版本作为目前能完成的最终的版本发出来供同行们使用,方便工作和分析。下面对工具包中的WindAnalysis1和WindAnalysis2的功能做个介绍,过一阵闲了编个教程发出来供大家使用。WindAnalysis1工具包能够对获取的整个测风数据构建dateset结构体,根据时间序列进行综合整理分析,通过运行可以获得如下分析结果:a.不同高度风速、风向、温度、压强的时间序列分布图; 风速、风向、温度、压强.jpg b.整个测风数据质量判断,及质量分析图; 测风数据质量评估.jpg c.不同高度湍流强度按照风速的分布、各风速对应的湍流强度与其平均湍流强度的分布图; 湍流分布.jpg d.不同高度月平均风速分布图; 月平均风速.jpg e.不同高度日平均风速分布图; 日平均风速.jpg f.不同高度风速频率分布直方图; 风频分布.jpg g.不同高度风速风向玫瑰图; 风向、风能玫瑰图.jpg h.风切变拟合和计算; 风切变拟合.jpg i.风切变系数随月分布图; 月风切变.jpg WindAnalysis2为针对特定高度H处的风资源进行详细分析,包括:a.测风时间序列上风速、湍流偏离测风周期内平均值的偏离程度; 风速、湍流时间序列分布.jpg b.风速的威布尔分布拟合和参数计算; 威布尔分布.jpg c.威布尔分布拟合的误差和相关系数R2的计算分析; 拟合误差分析.jpg d.风切变拟合和切边系数计算; 风切变拟合.jpg e.指定轮毂高度处的平均风速推算及威布尔分布拟合; 轮毂高度处威布尔分布.jpg f.根据选型风机的参数,绘制功率曲线和推力系数曲线; 功率特性曲线.jpg 不仅限于以上figure图文件的生成,还能够估算出指定轮毂高度hub(hub>H)测风塔处的发电量,在Command Window窗口中输出计算结果,作为风资源分析的参考。 计算结果.png WindAnalysis风数据分析工具包教程-V1.4.pdf WindAnalysis1-V1.4.1.zip WindAnalysis2-V1.4.1.zip -------------------------------------------------------------------
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值