首页 / 专利库 / 风能 / 风能利用系数 / 一种计及风资源利用率预测评估的方法

一种计及资源利用率预测评估的方法

阅读:989发布:2020-05-18

专利汇可以提供一种计及资源利用率预测评估的方法专利检索,专利查询,专利分析的服务。并且本 发明 一种计及 风 资源利用率预测评估的方法,涉及 风 力 发动机 的监控领域,该方法步骤是:构建风 电机 组功率五参数模型;灰狼优化 算法 进行功率模型参数识别;等间隔划分风速段,采用非参数核 密度 估计法建立各风速段的风资源利用率预测误差概率分布;进行计及风资源利用率区间估计,求得风电机组在给定 置信度 下的风资源利用率的置信区间;预测评估风资源利用率的 波动 程度,本发明克服了 现有技术 所存在的没有进行风资源利用率波动程度的预测评估的 缺陷 。,下面是一种计及资源利用率预测评估的方法专利的具体信息内容。

1.一种计及资源利用率预测评估的方法,其特征在于:将数据样本等间隔划分风速段,采用非参数核密度估计法,分别确定各风速段的预测误差ξ的概率分布F(ξ),再根据预测误差概率分布F(ξ)求得给定置信度下的风资源利用率Cpa的置信区间,进而预测评估风资源利用率的波动程度,具体步骤如下:
第一步,辨识风电机组功率五参数模型中的最优解五个参数:
根据实际风速v,引入灰狼优化算法辨识风电机组功率五参数模型的参数,具体操作过程如下:
(1.1)步,风电机组功率五参数模型:
设计风电机组功率五参数模型如下公式(1)所示,
公式(1)中, 为通过五参数模型计算出的功率预测值,v为由观测得到的实际风速,X=[a,b,c,d,g]为风电机组功率五参数模型的特定参数向量,a为预期的最大响应,b为斜率因子,c为过渡位置参数,d为最小响应,g为不对称参数;
(1.2)步,在给定范围内产生灰狼优化算法的随机灰狼种群:
初始化狼群Y,狼群中的灰狼数量为M,设置灰狼优化算法中的最大迭代次数为N,设定上述(1.1)步中的风电机组功率五参数的上下界,在给定上下界范围内,即在种群中灰狼数量为M的训练集以及参数的上下界范围内产生随机灰狼种群,该随机灰狼种群的特定参数向量的矩阵Y表示如公式(2)所示,
Y=[X1,X2,…,Xi,…,XM]T   (2),
公式(2)中,Xi为第i个个体特定参数位置向量,Xi=[a,b,c,d,g]都是由上述(1.1)步中的a,b,c,d,g五个参数组成,其中,1≤i≤M,M为狼群中的灰狼数量;
(1.3)步,计算每个个体特定参数位置向量的适应度值:
使用如下公式(3)所示均方根误差作为适应度函数,由此计算每个个体的特定参数位置向量适应度值F(Xi),
公式(3)中,vk为给定风速,其中1≤k≤n,n为给定风速样本个数,Pm(Xi,vk)为上述(1.2)中的第i个个体特定参数位置向量的确定的五参数模型功率,Pa(vk)为风速vk对应的风电机组实际功率;
(1.4)步,选择适应度值最好的3个个体特定参数位置向量:
适应度值F(Xi)的函数值越小,个体越符合要求,根据此原则选择确定上述(1.3)步计算适应度值F(Xi)最好的3个特定参数个体位置向量,分别记为Xα、Xβ、Xδ;
(1.5)步,更新适应度值最好的个体特定参数位置向量:
上述(1.4)步中的最好的3个个体特定参数位置向量中的Xα由如下公式(4-1)~(4-3)得到:
D=|C×XP(t)-X(t)|   (4-1),
X(t+1)=XP(t)-A×D   (4-2),
公式(4-1)中,D为个体特定参数位置向量与最优解的距离,t为当前迭代次数,C是与随机数r1有关的参数,C=2×r1,r1是[0,1]的随机数,Xp为最优解的个体特定参数位置向量,X(t)表示当前个体特定参数位置向量,
公式(4-2)中,A是与随机数r2有关的参数,A=2h×r2-h,收敛因子h随迭代次数从2线性递减到0,r2是[0,1]的随机数,
公式(4-3)中,N为最大迭代次数,
上述(1.4)步中的Xα带领Xβ和Xδ指导整个种群逼近最优解的数学描述如下,
Dα=|C1×Xα(t)-X(t)|   (5-1),
Dβ=|C2×Xβ(t)-X(t)|   (5-2),
Dδ=|C3×Xδ(t)-X(t)|   (5-3),
公式(5-1)~(5-3)中,C1、C2和C3为随机向量,Dα、Dβ、Dδ分别定义了Xα、Xβ、Xδ的前进步长和方向;
更新候选个体特定参数位置向量Xω的位置由以下公式确定,
X1=Xα-A1×Dα   (6-1),
X2=Xβ-A2×Dβ   (6-2),
X3=Xδ-A3×Dδ   (6-3),
X(t+1)=(X1+X2+X3)/3   (6-4),
公式(6-1)~(6-4)中,A1、A2、A3为随机向量,X(t+1)为更新候选个体特定参数位置向量Xω的位置;
(1.6)步,更新参数h、A、C:
随迭代进行不断更新参数a、A、C;
(1.7)步,如果当前迭代次数t小于最大迭代次数N,转到上述(1.3)步:
(1.8)步,若当前迭代次数t达到最大迭代次数N,输出当前适应度值F(Xi)最好的特定参数个体位置向量Xα的最终位置向量,其中的五个参数即为最优解;
至此,辨识了风电机组功率五参数模型中的最优解五个参数,即a为预期的最大响应,b为斜率因子,c为过渡位置参数,d为最小响应,g为不对称参数五个参数;
第二步,确定风资源利用率的预测误差ξ:
(2.1)步,采用风能利用系数Cp表征风资源利用率,Cp定义方法如下公式(7)所示,
3
公式(7)中,Cp为风资源利用率,ρ为空气密度,单位kg/m ,P为风电机组实际有功功率,单位W;S为叶轮扫风面积,单位m2;v为由观测得到的实际风速,单位m/s,
(2.2)步,采用预测误差统计分析方法,确定风资源利用率的预测误差ξ:
风资源利用率的预测误差ξ的确定方法如下公式(8)所示,
公式(8)中,Cpa为风资源利用率实际值,将观测得到的实际风速v及其对应的风电机组实际有功功率P代入公式(7),再根据此时的空气密度与叶轮扫风面积即可确定风资源利用率实际值,Cpm为风资源利用率预测值,根据实际风速由上述公式(1)求得风电机组功率预测值,将该风电机组功率预测值代入上述公式(7)来确定;Cpmax为理论上的风资源利用率最大值,根据贝兹理论,取0.593,
至此确定风资源利用率的预测误差ξ;
第三步,确定各风速段的风资源利用率Cp的预测误差ξi及其概率分布F(ξi):
所述各风速段为Di,其中i=1,2,···,λ,λ为风速段的所分段数,
采用非参数核密度估计法建立各风速段的风资源利用率Cp的预测误差ξi的概率分布F(ξi),
(3.1)等间隔划分风速段:
将数据样本等间隔划分成多个风速段,分别进行概率分布统计,设定每一风速段的段长为Δv,整体风速波动范围为[vl,vh],则所分某风速段Di由如下公式(9)所示:
Di=[vl+(i-1)Δv,vl+iΔv]   (9),
公式(9)中,i=1,2,···,λ,λ为风速段的所分段数,
λ=[(vh-vl)/Δv]+1
(3.2)求取各风速段的风资源利用率Cp的预测误差ξi的概率密度函数f(ξi),并绘出风资源利用率预测误差概率密度曲线:
采用非参数核密度估计法求取各风速段的风资源利用率Cp的预测误差ξi的概率密度函数f(ξi),并绘出风资源利用率预测误差概率密度曲线,操作如下:
取标准高斯核函数K(x)作为非参数核密度估计函数,表达式为公式(10),
对于某一风速段Di的风资源利用率Cp的预测误差ξi,其概率密度函数f(ξi)如下公式(11)所示,
公式(11)中,Ti为样本总数;L为窗宽;K(x)为标准高斯核函数;ξm为误差样本,并依次绘出风资源利用率预测误差概率密度曲线,横坐标为风资源利用率Cp的预测误差ξi,纵坐标为概率密度函数f(ξi);
(3.3)确定各风速段的风资源利用率Cp的预测误差ξi的概率分布F(ξi):
采用交叉验证法,选择各风速段的最佳窗宽L值,求得各风速段的风资源利用率Cp的预测误差ξi的概率分布F(ξi);
第四步,确定风电机组在给定置信度下的风资源利用率波动区间:
根据上述第三步所得各风速段的风资源利用率Cp的预测误差ξi的概率分布F(ξi),确定风电机组在给定置信度下的风资源利用率波动区间的过程如下,
利用上述第三步所得各风速段的风资源利用率Cp的预测误差ξi的概率分布F(ξi),做出给定置信度1-α下计及风资源利用率实际值Cpa的置信区间为公式(12)所示,
Cpa∈[Cpl,Cpu]1-α   (12),
公式(12)中,Cpl为风资源利用率置信下限;Cpu为风资源利用率置信上限;1-α为给定置信度,有1-α=P{CplCpl=Cpm+G(α/2)×Cpmax   (13),
Cpu=Cpm+G(1-α/2)×Cpmax   (14),
公式(13)和公式(14)中,G(η)为各风速段的风资源利用率Cp的预测误差ξi的概率分布F(ξi)的反函数,其中η如下公式(15)所示,
P{ξ≤G(η)}=η   (15),
由此,确定风电机组在给定置信度下的风资源利用率波动区间为[Cpl,Cpu];
第五步,预测评估风资源利用率的波动程度:
(5.1)步,对于某一风资源利用率的预测值Cpm,找到其对应的实际风速,判断其所在风速段;
(5.2)步,求得预测误差ξ的α/2和1-α/2对应点:
根据上述第三步中已经求得的各风速段的风资源利用率预测误差概率密度f(ξi),查找该风资源利用率预测值所处风速段的风资源利用率Cp预测误差ξi的概率密度曲线,求得预测误差ξ的α/2和1-α/2对应点;
(5.3)步,根据上述第四步中给出的风资源利用率置信区间的置信下限Cpl和风资源利用率置信上限Cpu的计算式公式(13)和公式(14),求出上述(5.2)步中所述风资源利用率Cp预测值所处风速段的区间估计的边界值,进而给出给定置信度1-α下该风资源利用率Cp预测值的置信区间;
(5.4)步,同理求得各风速段的所有风资源利用率Cp预测值的置信区间,分别连接各个风资源利用率Cp预测值置信区间上下限,形成全风速段的风资源利用率Cp估计曲线;
至此,根据全风速段的风资源利用率Cp估计曲线预测评估风资源利用率的波动程度。

说明书全文

一种计及资源利用率预测评估的方法

技术领域

[0001] 本发明的技术方案涉及发动机的监控领域,具体地说是一种计及风资源利用率预测评估的方法。

背景技术

[0002] 风资源利用率反映风电机组性能,风资源利用率预测评估是风电机组状态监测和风电功率预测的有效手段之一。风能的间歇性、随机性和波动性使得风电机组发电性能偏离设计指标、风资源利用率稳定性差,这给电网的安全运行和经济调度带来极大挑战。
[0003] 目前,研究风资源利用率预测评估的文献较少,一般的风资源评估方法很少对风资源利用率进行评估,并且仅局限于确定性的评估。CN105005581A公开了一种风电场风资源数据处理方法及系统,该技术仅进行风速风向等数据的监测,未对风资源的利用率及风电机组健康状态进行监测评估;CN104331621B公开了一种风资源计算方法,该方法并未充分考虑测风地以及风电机组的运行状态,只是根据单一指标进行简单的风资源等级划分;CN103696913B公开了一种实时监测风电机组运行时的风能利用偏差的方法,利用偏差的方法实时监测风能利用偏差,但未对风能利用率的波动程度进行评估。总之,现有的技术方案没有进行风资源利用率波动程度的预测评估。

发明内容

[0004] 本发明所要解决的技术问题是:提供一种计及风资源利用率预测评估的方法,该方法将数据样本等间隔划分风速段,采用非参数核密度估计法,分别求得各风速段的预测误差ξ的概率分布F(ξ),再根据预测误差概率分布F(ξ)求得给定置信度下的风资源利用率Cpa的置信区间,进而预测评估风资源利用率的波动程度,克服了现有技术所存在的没有进行风资源利用率波动程度的预测评估的缺陷
[0005] 本发明解决该技术问题所采用的技术方案是:一种计及风资源利用率预测评估的方法,将数据样本等间隔划分风速段,采用非参数核密度估计法,分别确定各风速段的预测误差ξ的概率分布F(ξ),再根据预测误差概率分布F(ξ)求得给定置信度下的风资源利用率Cpa的置信区间,进而预测评估风资源利用率的波动程度,具体步骤如下:
[0006] 第一步,辨识风电机组功率五参数模型中的最优解五个参数:
[0007] 根据实际风速v,引入灰狼优化算法辨识风电机组功率五参数模型的参数,具体操作过程如下:
[0008] (1.1)步,风电机组功率五参数模型:
[0009] 设计风电机组功率五参数模型如下公式(1)所示,
[0010]
[0011] 公式(1)中, 为通过五参数模型计算出的功率预测值,v为由观测得到的实际风速,X=[a,b,c,d,g]为风电机组功率五参数模型的特定参数向量,a为预期的最大响应,b为斜率因子,c为过渡位置参数,d为最小响应,g为不对称参数;
[0012] (1.2)步,在给定范围内产生灰狼优化算法的随机灰狼种群:
[0013] 初始化狼群Y,狼群中的灰狼数量为M,设置灰狼优化算法中的最大迭代次数为N,设定上述(1.1)步中的风电机组功率五参数的上下界,在给定上下界范围内,即在种群中灰狼数量为M的训练集以及参数的上下界范围内产生随机灰狼种群,该随机灰狼种群的特定参数向量的矩阵Y表示如公式(2)所示,
[0014] Y=[X1,X2,…,Xi,…,XM]T  (2),
[0015] 公式(2)中,Xi为第i个个体特定参数位置向量,Xi=[a,b,c,d,g]都是由上述(1.1)步中的a,b,c,d,g五个参数组成,其中,1≤i≤M,M为狼群中的灰狼数量;
[0016] (1.3)步,计算每个个体特定参数位置向量的适应度值:
[0017] 使用如下公式(3)所示均方根误差作为适应度函数,由此计算每个个体的特定参数位置向量适应度值F(Xi),
[0018]
[0019] 公式(3)中,vk为给定风速,其中1≤k≤n,n为给定风速样本个数,Pm(Xi,vk)为上述(1.2)中的第i个个体特定参数位置向量的确定的五参数模型功率,Pa(vk)为风速vk对应的风电机组实际功率;
[0020] (1.4)步,选择适应度值最好的3个个体特定参数位置向量:
[0021] 适应度值F(Xi)的函数值越小,个体越符合要求,根据此原则选择确定上述(1.3)步计算适应度值F(Xi)最好的3个特定参数个体位置向量,分别记为Xα、Xβ、Xδ;
[0022] (1.5)步,更新适应度值最好的个体特定参数位置向量:
[0023] 上述(1.4)步中的最好的3个个体特定参数位置向量中的Xα由如下公式(4-1)~(4-3)得到:
[0024] D=|C×XP(t)-X(t)|  (4-1),
[0025] X(t+1)=XP(t)-A×D  (4-2),
[0026]
[0027] 公式(4-1)中,D为个体特定参数位置向量与最优解的距离,t为当前迭代次数,C是与随机数r1有关的参数,C=2×r1,r1是[0,1]的随机数,Xp为最优解的个体特定参数位置向量,X(t)表示当前个体特定参数位置向量,
[0028] 公式(4-2)中,A是与随机数r2有关的参数,A=2h×r2-h,收敛因子h随迭代次数从2线性递减到0,r2是[0,1]的随机数,
[0029] 公式(4-3)中,N为最大迭代次数,
[0030] 上述(1.4)步中的Xα带领Xβ和Xδ指导整个种群逼近最优解的数学描述如下,[0031] Dα=|C1×Xα(t)-X(t)|  (5-1),
[0032] Dβ=|C2×Xβ(t)-X(t)|  (5-2),
[0033] Dδ=|C3×Xδ(t)-X(t)|  (5-3),
[0034] 公式(5-1)~(5-3)中,C1、C2和C3为随机向量,Dα、Dβ、Dδ分别定义了Xα、Xβ、Xδ的前进步长和方向;
[0035] 更新候选个体特定参数位置向量Xω的位置由以下公式确定,
[0036] X1=Xα-A1×Dα  (6-1),
[0037] X2=Xβ-A2×Dβ  (6-2),
[0038] X3=Xδ-A3×Dδ  (6-3),
[0039] X(t+1)=(X1+X2+X3)/3  (6-4),
[0040] 公式(6-1)~(6-4)中,A1、A2、A3为随机向量,X(t+1)为更新候选个体特定参数位置向量Xω的位置;
[0041] (1.6)步,更新参数h、A、C:
[0042] 随迭代进行不断更新参数a、A、C;
[0043] (1.7)步,如果当前迭代次数t小于最大迭代次数N,转到上述(1.3)步:
[0044] (1.8)步,若当前迭代次数t达到最大迭代次数N,输出当前适应度值F(Xi)最好的特定参数个体位置向量Xα的最终位置向量,其中的五个参数即为最优解;
[0045] 至此,辨识了风电机组功率五参数模型中的最优解五个参数,即a为预期的最大响应,b为斜率因子,c为过渡位置参数,d为最小响应,g为不对称参数五个参数;
[0046] 第二步,确定风资源利用率的预测误差ξ:
[0047] (2.1)步,采用风能利用系数Cp表征风资源利用率,Cp定义方法如下公式(7)所示,[0048]
[0049] 公式(7)中,Cp为风资源利用率,ρ为空气密度,单位kg/m3,P为风电机组实际有功功率,单位W;S为叶轮扫风面积,单位m2;v为由观测得到的实际风速,单位m/s,
[0050] (2.2)步,采用预测误差统计分析方法,确定风资源利用率的预测误差ξ:
[0051] 风资源利用率的预测误差ξ的确定方法如下公式(8)所示,
[0052]
[0053] 公式(8)中,Cpa为风资源利用率实际值,将观测得到的实际风速v及其对应的风电机组实际有功功率P代入公式(7),再根据此时的空气密度与叶轮扫风面积即可确定风资源利用率实际值,Cpm为风资源利用率预测值,根据实际风速由上述公式(1)求得风电机组功率预测值,将该风电机组功率预测值代入上述公式(7)来确定;Cpmax为理论上的风资源利用率最大值,根据贝兹理论,取0.593,
[0054] 至此确定风资源利用率的预测误差ξ;
[0055] 第三步,确定各风速段的风资源利用率Cp的预测误差ξi及其概率分布F(ξi):
[0056] 所述各风速段为Di,其中i=1,2,···,λ,λ为风速段的所分段数,
[0057] 采用非参数核密度估计法建立各风速段的风资源利用率Cp的预测误差ξi的概率分布F(ξi),
[0058] (3.1)等间隔划分风速段:
[0059] 将数据样本等间隔划分成多个风速段,分别进行概率分布统计,设定每一风速段的段长为Δv,整体风速波动范围为[vl,vh],则所分某风速段Di由如下公式(9)所示:
[0060] Di=[vl+(i-1)Δv,vl+iΔv]  (9),
[0061] 公式(9)中,i=1,2,···,λ,λ为风速段的所分段数,
[0062] λ=[(vh-vl)/Δv]+1
[0063] (3.2)求取各风速段的风资源利用率Cp的预测误差ξi的概率密度函数f(ξi),并绘出风资源利用率预测误差概率密度曲线:
[0064] 采用非参数核密度估计法求取各风速段的风资源利用率Cp的预测误差ξi的概率密度函数f(ξi),并绘出风资源利用率预测误差概率密度曲线,操作如下:
[0065] 取标准高斯核函数K(x)作为非参数核密度估计函数,表达式为公式(10),[0066]
[0067] 对于某一风速段Di的风资源利用率Cp的预测误差ξi,其概率密度函数f(ξi)如下公式(11)所示,
[0068]
[0069] 公式(11)中,Ti为样本总数;L为窗宽;K(x)为标准高斯核函数;ξm为误差样本,[0070] 并依次绘出风资源利用率预测误差概率密度曲线,横坐标为风资源利用率Cp的预测误差ξi,纵坐标为概率密度函数f(ξi);
[0071] (3.3)确定各风速段的风资源利用率Cp的预测误差ξi的概率分布F(ξi):
[0072] 采用交叉验证法,选择各风速段的最佳窗宽L值,求得各风速段的风资源利用率Cp的预测误差ξi的概率分布F(ξi);
[0073] 第四步,确定风电机组在给定置信度下的风资源利用率波动区间:
[0074] 根据上述第三步所得各风速段的风资源利用率Cp的预测误差ξi的概率分布F(ξi),确定风电机组在给定置信度下的风资源利用率波动区间的过程如下,
[0075] 利用上述第三步所得各风速段的风资源利用率Cp的预测误差ξi的概率分布F(ξi),做出给定置信度1-α下计及风资源利用率实际值Cpa的置信区间为公式(12)所示,[0076] Cpa∈[Cpl,Cpu]1-α  (12),
[0077] 公式(12)中,Cpl为风资源利用率置信下限;Cpu为风资源利用率置信上限;1-α为给定置信度,有1-α=P{Cpl
[0078] Cpl=Cpm+G(α/2)×Cpmax  (13),
[0079] Cpu=Cpm+G(1-α/2)×Cpmax  (14),
[0080] 公式(13)和公式(14)中,G(η)为各风速段的风资源利用率Cp的预测误差ξi的概率分布F(ξi)的反函数,其中η如下公式(15)所示,
[0081] P{ξ≤G(η)}=η  (15),
[0082] 由此,确定风电机组在给定置信度下的风资源利用率波动区间为[Cpl,Cpu];
[0083] 第五步,预测评估风资源利用率的波动程度:
[0084] (5.1)步,对于某一风资源利用率的预测值Cpm,找到其对应的实际风速,判断其所在风速段;
[0085] (5.2)步,求得预测误差ξ的α/2和1-α/2对应点:
[0086] 根据上述第三步中已经求得的各风速段的风资源利用率预测误差概率密度f(ξi),查找该风资源利用率预测值所处风速段的风资源利用率Cp预测误差ξi的概率密度曲线,求得预测误差ξ的α/2和1-α/2对应点;
[0087] (5.3)步,根据上述第四步中给出的风资源利用率置信区间的置信下限Cpl和风资源利用率置信上限Cpu的计算式公式(13)和公式(14),求出上述(5.2)步中所述风资源利用率Cp预测值所处风速段的区间估计的边界值,进而给出给定置信度1-α下该风资源利用率Cp预测值的置信区间;
[0088] (5.4)步,同理求得各风速段的所有风资源利用率Cp预测值的置信区间,分别连接各个风资源利用率Cp预测值置信区间上下限,形成全风速段的风资源利用率Cp估计曲线;
[0089] 至此,根据全风速段的风资源利用率Cp估计曲线预测评估风资源利用率的波动程度。
[0090] 上述一种计及风资源利用率预测评估的方法,所述非参数核密度估计法的英文为nonparametric kernel density estimation,缩写为NKDE;交叉验证法的英文为cross validation,缩写为CV;灰狼优化算法的英文缩写为GWO,均是本技术领域公知的方法。
[0091] 本发明的有益效果是:与现有技术相比,本发明具有以下突出的实质性特点和显著进步:
[0092] (1)本发明提出将数据样本等间隔划分风速段,采用非参数核密度估计法,分别求得各风速段的预测误差ξ的概率分布F(ξ),再根据预测误差概率分布F(ξ)求得给定置信度下的风资源利用率Cpa的置信区间,进而预测评估风资源利用率的波动程度,克服了现有的技术方案没有进行风资源利用率波动程度的预测评估的缺陷。
[0093] (2)CN105005581A公开了一种风电场风资源数据处理方法及系统,仅进行风速风向等数据的监测,未对风资源的利用率及风电机组健康状态进行监测评估;CN104331621B公开了一种风资源计算方法,并未充分考虑测风地以及风电机组的运行状态,只是根据单一指标进行简单的风资源等级划分,未对风资源利用率波动程度进行预测评估;CN103696913B公开了一种实时监测风电机组运行时的风能利用偏差的方法,利用偏差的方法实时监测风能利用偏差,但未对未来风能利用率的波动程度进行评估;总之,上述现有技术没有对风资源利用率波动程度进行预测评估,仅局限于风资源利用率确定性的评估。
[0094] (3)风资源利用率反映风电机组性能,风资源利用率预测评估是风电机组状态监测和风电功率预测的有效手段之一。风能的间歇性、随机性和波动性使得风电机组发电性能偏离设计指标、风资源利用率稳定性差,这给电网的安全运行和经济调度带来极大挑战。
[0095] 风力发电的不确定性导致风资源利用率的不确定性,本发明一种计及风资源利用率预测评估的方法,合理评估风电机组健康状态及其风资源利用率,并在此基础上分析风资源利用率的波动程度,实现科学预测风电功率,这些是大规模风电并网必须予以解决的关键问题。附图说明
[0096] 下面结合附图和实施例对本发明进一步说明。
[0097] 图1是实施例1中处理后的风速-功率散点图。
[0098] 图2是实施例1中灰狼优化算法模型参数辨识结果的风速-功率曲线图。
[0099] 图3(a)是实施例1中风电机组风资源利用率Cp的实际值测试图。
[0100] 图3(b)是实施例1中风电机组风资源利用率Cp的预测值测试图。
[0101] 图4是实施例1中风电机组风速分布特征柱状图。
[0102] 图5是实施例1中不同风速下的风资源利用率Cp预测误差散点图。
[0103] 图6是实施例1中不同风速段风资源利用率Cp预测误差ξ概率密度f(ξ)曲线图。
[0104] 图7(a)是实施例1中风速段4的风资源利用率Cp概率分布及80%置信区间曲线图。
[0105] 图7(b)是实施例1中风速段7的风资源利用率Cp概率分布及60%置信区间曲线图。
[0106] 图8(a)是实施例1中风资源利用率Cp实际值测试图。
[0107] 图8(b)是实施例1中80%置信区间(1-α)=0.8的风资源利用率Cp实际值测试图。
[0108] 图8(c)是实施例160%置信区间(1-α)=0.6的风资源利用率Cp实际值测试图。
[0109] 图9是本发明方法操作流程示意图。

具体实施方式

[0110] 图9所示实施例表明,本发明方法操作流程是:构建风电机组功率五参数模型→灰狼优化算法进行功率模型参数识别→等间隔划分风速段,采用非参数核密度估计法建立各风速段的风资源利用率预测误差概率分布→进行计及风资源利用率区间估计,求得风电机组在给定置信度下的风资源利用率的置信区间→预测评估风资源利用率的波动程度。
[0111] 实施例
[0112] 本实施例一种计及风资源利用率预测评估的方法,选用的2MW风电机组的技术参数如下表1所示:
[0113] 表1.2MW风电机组的技术参数
[0114]
[0115] 统计某典型风电场SCADA系统中的实际风速v、风电机组实际有功功率P全年运行数据,进行数据预处理,绘制成P-v散点图,如图1本实施例中处理后的风速-功率散点图所示。筛选后的数据被随机分为2个子集,都分别包括训练样本和测试样本。第1个子集用于风电机组功率模型参数辨识及模型评价,第2个子集用于风资源利用率预测、误差概率分布统计及置信区间估计。
[0116] 本实施例一种计及风资源利用率预测评估的方法,具体步骤如下:
[0117] 第一步,辨识风电机组功率五参数模型中的最优解五个参数:
[0118] 根据实际风速v,引入灰狼优化算法辨识风电机组功率五参数模型的参数,具体操作过程如下:
[0119] (1.1)步,风电机组功率五参数模型:
[0120] 设计风电机组功率五参数模型如下公式(1)所示,
[0121]
[0122] 公式(1)中, 为通过五参数模型计算出的功率预测值,v为由观测得到的实际风速,X=[a,b,c,d,g]为风电机组功率五参数模型的特定参数向量,a为预期的最大响应,b为斜率因子,c为过渡位置参数,d为最小响应,g为不对称参数;
[0123] 设定,上述风电机组功率五参数模型中a,b,c,d,g的上下限分别为[1950,2050],[-20,20],[0,20],[-20,20],[0,2],精度取小数点后四位;
[0124] (1.2)步,在给定范围内产生灰狼优化算法的随机灰狼种群:
[0125] 初始化狼群Y,狼群中的灰狼数量为M,设置灰狼优化算法中的最大迭代次数为N,设定上述(1.1)步中的风电机组功率五参数的上下界,在给定上下界范围内,即在种群中灰狼数量为M的训练集以及参数的上下界范围内产生随机灰狼种群,在本实施例中,初始种群数量M=50,最大迭代次数N=200,该随机灰狼种群的特定参数向量的矩阵Y表示如公式(2)所示,
[0126] Y=[X1,X2,…,Xi,…,XM]T  (2),
[0127] 公式(2)中,Xi为第i个个体特定参数位置向量,Xi=[a,b,c,d,g]都是由上述(1.1)步中的a,b,c,d,g五个参数组成,其中,1≤i≤M,M为狼群中的灰狼数量;
[0128] (1.3)步,计算每个个体特定参数位置向量的适应度值:
[0129] 使用如下公式(3)所示均方根误差作为适应度函数,由此计算每个个体的特定参数位置向量适应度值F(Xi),
[0130]
[0131] 公式(3)中,vk为给定风速,其中1≤k≤n,n为给定风速样本个数,Pm(Xi,vk)为上述(1.2)中的第i个个体特定参数位置向量的确定的五参数模型功率,Pa(vk)为风速vk对应的风电机组实际功率;
[0132] (1.4)步,选择适应度值最好的3个个体特定参数位置向量:
[0133] 适应度值F(Xi)的函数值越小,个体越符合要求,根据此原则选择确定上述(1.3)步计算适应度值F(Xi)最好的3个特定参数个体位置向量,分别记为Xα、Xβ、Xδ;
[0134] (1.5)步,更新适应度值最好的个体特定参数位置向量:
[0135] 上述(1.4)步中的最好的3个个体特定参数位置向量中的Xα由如下公式(4-1)~(4-3)得到:
[0136] D=|C×XP(t)-X(t)|  (4-1),
[0137] X(t+1)=XP(t)-A×D  (4-2),
[0138]
[0139] 公式(4-1)中,D为个体特定参数位置向量与最优解的距离,t为当前迭代次数,C是与随机数r1有关的参数,C=2×r1,r1是[0,1]的随机数,Xp为最优解的个体特定参数位置向量,X(t)表示当前个体特定参数位置向量,
[0140] 公式(4-2)中,A是与随机数r2有关的参数,A=2h×r2-h,收敛因子h随迭代次数从2线性递减到0,r2是[0,1]的随机数,
[0141] 公式(4-3)中,N为最大迭代次数,
[0142] 上述(1.4)步中的Xα带领Xβ和Xδ指导整个种群逼近最优解的数学描述如下,[0143] Dα=|C1×Xα(t)-X(t)|  (5-1),
[0144] Dβ=|C2×Xβ(t)-X(t)|  (5-2),
[0145] Dδ=|C3×Xδ(t)-X(t)|  (5-3),
[0146] 公式(5-1)~(5-3)中,C1、C2和C3为随机向量,Dα、Dβ、Dδ分别定义了Xα、Xβ、Xδ的前进步长和方向;
[0147] 更新候选个体特定参数位置向量Xω的位置由以下公式确定,
[0148] X1=Xα-A1×Dα  (6-1),
[0149] X2=Xβ-A2×Dβ  (6-2),
[0150] X3=Xδ-A3×Dδ  (6-3),
[0151] X(t+1)=(X1+X2+X3)/3  (6-4),
[0152] 公式(6-1)~(6-4)中,A1、A2、A3为随机向量,X(t+1)为更新候选个体特定参数位置向量Xω的位置;
[0153] (1.6)步,更新参数h、A、C:
[0154] 随迭代进行不断更新参数a、A、C;
[0155] (1.7)步,如果当前迭代次数t小于最大迭代次数N,转到上述(1.3)步:
[0156] (1.8)步,若当前迭代次数t达到最大迭代次数N,输出当前适应度值F(Xi)最好的特定参数个体位置向量Xα的最终位置向量,其中的五个参数即为最优解;
[0157] 至此,辨识了风电机组功率五参数模型中的最优解五个参数,即a为预期的最大响应,b为斜率因子,c为过渡位置参数,d为最小响应,g为不对称参数五个参数,本实施例中灰狼优化算法模型参数辨识结果的风速-功率曲线图如图2所示,经上述步骤(1.3)步~(1.8)步操作,某次运算得到的最优解[a,b,c,d,g]=[1995.4289,-19.8367,11.9052,0.5638,0.1601];
[0158] 第二步,确定风资源利用率的预测误差ξ:
[0159] (2.1)步,采用风能利用系数Cp表征风资源利用率,Cp定义方法如下公式(7)所示,[0160]
[0161] 公式(7)中,Cp为风资源利用率,ρ为空气密度,单位kg/m3,P为风电机组实际有功功率,单位W;S为叶轮扫风面积,单位m2;v为由观测得到的实际风速,单位m/s,
[0162] (2.2)步,采用预测误差统计分析方法,确定风资源利用率的预测误差ξ:
[0163] 风资源利用率的预测误差ξ的确定方法如下公式(8)所示,
[0164]
[0165] 公式(8)中,Cpa为风资源利用率实际值,将观测得到的实际风速v及其对应的风电机组实际有功功率P代入公式(7),再根据此时的空气密度与叶轮扫风面积即可确定风资源利用率实际值,Cpm为风资源利用率预测值,根据实际风速由上述公式(1)求得风电机组功率预测值,将该风电机组功率预测值代入上述公式(7)来确定;Cpmax为理论上的风资源利用率最大值,根据贝兹理论,取0.593,
[0166] 图3(a)为本实施例中风电机组风资源利用率Cp的实际值测试图,图3(b)为本实施例计算后的中风电机组风资源利用率Cp的预测值测试图,如图3(a)与图3(b)所示,由于预测方法、模型精度及自然环境的差异,风电机组Cp的实际值与预测值并不完全相符。实际运行的风电场中,风速服从威布尔概率分布,
[0167] 图4是本实施例中该风电机组风速分布特征柱状图。由该图可知,该风电机组所捕获的风速主要集中于3m/s和13m/s之间,因此主要对此风速区间的风资源利用率预测误差进行研究。此区间风速主要集中于vC到vN段,对应于功率特性曲线的线性放大区,此时较小的风速变化就会引起较大的功率波动,从而引起较大的Cp预测误差。
[0168] 图5是本实施例中不同风速下的风资源利用率Cp预测误差散点图。由图5可知,预测误差总体上在零附近波动,随着风速增大,波动逐渐减小,表2统计了各风速段内Cp预测误差的均值μ和标准差σ。
[0169] 表2.各风速段内风资源利用率的预测误差ξ的均值和标准差
[0170]
[0171] 由表2可见,各风速段的Cp预测误差平有明显差别,因此此处以Δv=1m/s把风速3~13m/s的风速区间等分为10个风速段;
[0172] 至此确定风资源利用率的预测误差ξ;
[0173] 第三步,确定各风速段的风资源利用率Cp的预测误差ξi及其概率分布F(ξi):
[0174] 所述各风速段为Di,其中i=1,2,···,λ,λ为风速段的所分段数,
[0175] 采用非参数核密度估计法建立各风速段的风资源利用率Cp的预测误差ξi的概率分布F(ξi),
[0176] (3.1)等间隔划分风速段:
[0177] 将数据样本等间隔划分成多个风速段,分别进行概率分布统计,设定每一风速段的段长为Δv,整体风速波动范围为[vl,vh],则所分某风速段Di由如下公式(9)所示:
[0178] Di=[vl+(i-1)Δv,vl+iΔv]  (9),
[0179] 公式(9)中,i=1,2,···,λ,λ为风速段的所分段数,
[0180] λ=[(vh-vl)/Δv]+1
[0181] (3.2)求取各风速段的风资源利用率Cp的预测误差ξi的概率密度函数f(ξi),并绘出风资源利用率预测误差概率密度曲线:
[0182] 采用非参数核密度估计法求取各风速段的风资源利用率Cp的预测误差ξi的概率密度函数f(ξi),并绘出风资源利用率预测误差概率密度曲线,操作如下:
[0183] 取标准高斯核函数K(x)作为非参数核密度估计函数,表达式为公式(10),[0184]
[0185] 对于某一风速段Di的风资源利用率Cp的预测误差ξi,其概率密度函数f(ξi)如下公式(11)所示,
[0186]
[0187] 公式(11)中,Ti为样本总数;L为窗宽;K(x)为标准高斯核函数;ξm为误差样本,[0188] 并依次绘出风资源利用率预测误差概率密度曲线,横坐标为风资源利用率Cp的预测误差ξi,纵坐标为概率密度函数f(ξi),根据本实施例上述所述第(3.2)步分别计算各风速段内每一点的概率密度,采用三次样条插值求取各段Cp预测误差的概率密度曲线,如图6本实施例中不同风速段风资源利用率Cp预测误差ξ概率密度f(ξ)曲线图所示,风速段9和10的预测误差较小,其概率分布主要集中在零值附近,说明达到额定风速以后Cp值变化较小,预测值与实际值更为接近,而其它风速段的预测误差较为分散,说明Cp值波动范围较大,准确预测更为困难;
[0189] (3.3)确定各风速段的风资源利用率Cp的预测误差ξi的概率分布F(ξi):
[0190] 采用交叉验证法,选择各风速段的最佳窗宽L值,求得各风速段的风资源利用率Cp的预测误差ξi的概率分布F(ξi);
[0191] 第四步,确定风电机组在给定置信度下的风资源利用率波动区间:
[0192] 根据上述第三步所得各风速段的风资源利用率Cp的预测误差ξi的概率分布F(ξi),确定风电机组在给定置信度下的风资源利用率波动区间的过程如下,
[0193] 利用上述第三步所得各风速段的风资源利用率Cp的预测误差ξi的概率分布F(ξi),做出给定置信度1-α下计及风资源利用率实际值Cpa的置信区间为公式(12)所示,[0194] Cpa∈[Cpl,Cpu]1-α  (12),
[0195] 公式(12)中,Cpl为风资源利用率置信下限;Cpu为风资源利用率置信上限;1-α为给定置信度,有1-α=P{Cpl
[0196] Cpl=Cpm+G(α/2)×Cpmax  (13),
[0197] Cpu=Cpm+G(1-α/2)×Cpmax  (14),
[0198] 公式(13)和公式(14)中,G(η)为各风速段的风资源利用率Cp的预测误差ξi的概率分布F(ξi)的反函数,其中η如下公式(15)所示,
[0199] P{ξ≤G(η)}=η  (15),
[0200] 由此,确定风电机组在给定置信度下的风资源利用率波动区间为[Cpl,Cpu];
[0201] 第五步,预测评估风资源利用率的波动程度:
[0202] (5.1)步,对于某一风资源利用率的预测值Cpm,找到其对应的实际风速,判断其所在风速段;
[0203] (5.2)步,求得预测误差ξ的α/2和1-α/2对应点:
[0204] 根据上述第三步中已经求得的各风速段的风资源利用率预测误差概率密度f(ξi),查找该风资源利用率预测值所处风速段的风资源利用率Cp预测误差ξi的概率密度曲线,求得预测误差ξ的α/2和1-α/2对应点;
[0205] (5.3)步,根据上述第四步中给出的风资源利用率置信区间的置信下限Cpl和风资源利用率置信上限Cpu的计算式公式(13)和公式(14),求出上述(5.2)步中所述风资源利用率Cp预测值所处风速段的区间估计的边界值,进而给出给定置信度1-α下该风资源利用率Cp预测值的置信区间;
[0206] (5.4)步,同理求得各风速段的所有风资源利用率Cp预测值的置信区间,分别连接各个风资源利用率Cp预测值置信区间上下限,形成全风速段的风资源利用率Cp估计曲线;
[0207] 至此,根据全风速段的风资源利用率Cp估计曲线预测评估风资源利用率的波动程度。
[0208] 由预测误差概率密度曲线经过折算可得Cp的概率密度曲线,选取连续545个测试样本分别进行置信度为0.8和0.6的风资源利用率评估;
[0209] 图7(a)是本实施例中风速段4的风资源利用率Cp概率分布及80%置信区间曲线图,为本实施例第98个测点所在的风速段4的Cp概率密度曲线及80%置信区间,其预测值为0.2937,置信区间为[0.1707,0.4659],实际值为0.2432;图7(b)是本实施例中风速段7的风资源利用率Cp概率分布及60%置信区间曲线图,为本实施例第484个测点所在风速段7的Cp概率密度曲线及60%置信区间,其预测值为0.3155,置信区间为[0.1798,0.4102],实际值为0.2879。两测点都恰好被包含在置信区间内,表明区间估计可靠,
[0210] 由图7(a)和图7(b)可知,基于非参数估计预测误差概率分布求得的置信区间并不关于预测值对称;不同风速段的Cp概率分布有明显差别,这也说明本发明采用风速分段的非参数估计方法是合理的;
[0211] 求出测试样本中所有Cp预测值在给定置信度下的置信区间,分别连接上下限形成全风速段的置信区间,结果如图8(b)本实施例中80%置信区间(1-α)=0.8的风资源利用率Cp实际值测试图和图8(c)本实施例中60%置信区间(1-α)=0.6的风资源利用率Cp实际值测试图所示,
[0212] 由图8(a)本实施例中风资源利用率Cp实际值测试图、图8(b)和图8(c)可知:给定置信度,本发明所用方法得到的置信区间能够包含大部分实际值;置信度越低,置信区间宽度越小,因为置信区间越小,其包含实际值的概率就越小。此外,由于原始数据中的弃风、停机和欠功率运行等因素,所以Cp置信区间的估计范围相对保守,主要受数据预处理技术水平的限制。
[0213] 上述一种计及风资源利用率预测评估的方法,所述非参数核密度估计法的英文为nonparametric kernel density estimation,缩写为NKDE;交叉验证法的英文为cross validation,缩写为CV;灰狼优化算法的英文缩写为GWO,均是本技术领域公知的方法。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈