CN107330183A - 一种基于运行数据的风电利用率计算方法 - Google Patents
一种基于运行数据的风电利用率计算方法 Download PDFInfo
- Publication number
- CN107330183A CN107330183A CN201710514679.XA CN201710514679A CN107330183A CN 107330183 A CN107330183 A CN 107330183A CN 201710514679 A CN201710514679 A CN 201710514679A CN 107330183 A CN107330183 A CN 107330183A
- Authority
- CN
- China
- Prior art keywords
- mrow
- power
- msub
- wind
- data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Wind Motors (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
本发明公开了一种基于运行数据的风电利用率计算方法,所述方法包括以下步骤:1)基于风机运行原理对异常数据进行初步筛选剔除;2)采用参数优化的四分位法进一步剔除异常数据;3)采用聚类算法深入剔除异常数据;4)采用分段Hermite插值进行缺失数据重构;5)建立相关型号风电机组的理论功率计算模型及风电场理论功率恢复;6)风电利用率相关指标计算。本发明方法基于风电场SCADA(Supervisory Control And Data Acquisition,数据采集与监视控制***)***记录的风机运行数据,结合有效的数据质量控制手段提出了风电场理论功率恢复算法,这为风电利用率的计算提供了一种合理的途径,也为风电场运行状况评估及风电、电网的规划发展提供了科学的依据。
Description
技术领域
本发明属于风电发配电技术领域,涉及一种基于运行数据的风电利用率计算方法。
背景技术
随着风电的快速发展,其间歇性和随机性的特点导致了风电弃风限电现象日益严峻,严重影响到风电场真实出力情况的评价以及对风电、电网发展的科学规划。弃风限电是基于电网的运行原理,在风电大发期间,为保证电网的安全稳定运行,在一定程度上限制风电出力的情况。如何通过技术手段恢复弃风限电期间的风电理论出力,改进并提高风电场理论功率恢复精度,进而计算出准确的弃风情况及风电利用情况,对于风电利用率的可靠计算,进而合理评价风电场/群运行的安全性和经济性、科学衡量风电消纳情况等实际问题具有重大的理论研究价值和应用指导意义。
风电理论功率恢复是风电利用率计算中的关键步骤。目前,国内风电行业对计算风电理论功率还没有形成标准。现有的计算方法一般是直接采用风机的出厂功率曲线作为理论功率计算曲线,但由于出厂功率曲线是在严格的试验条件下得到的,易受天气状况、地形和风机自身出力特性变化等影响,因而计算误差较大。
发明内容
本项发明利用无弃风条件下风电机组的实际运行数据,通过异常数据剔除、缺失数据重构等数据质量控制手段获取有效的风速、功率、风轮转速等实际运行数据。在此基础上,分别对不同型号的风电机组建立相应理论功率计算模型。然后,根据风速测量数据,分别恢复特定型号每台风机的理论功率。最后,通过累加得到风电场的理论功率,进而实现风电利用率的准确计算及风电利用情况的合理评价。基于运行数据的风电场理论功率恢复方法将为风电利用率的计算提供一种有效途径,也将为风电场运行情况评估和风电、电网的规划发展提供合理手段和科学依据。
本发明的具体技术方案是:
一种基于运行数据的风电利用率计算方法,其特征在于,所述方法包括以下步骤:
1)基于风机运行原理对异常数据进行初步筛选剔除;
2)采用参数优化的四分位法进一步剔除异常数据;
3)采用聚类算法深入剔除异常数据;
4)采用分段Hermite插值进行缺失数据重构;
5)建立相关型号风电机组的理论功率计算模型及风电场理论功率恢复;
6)风电利用率相关指标计算。
进一步,所述步骤1)包括:
S101:根据风机运行原理将其划分为5个运行区间:
区间1:0≤V<Vcut_in,P=0;
区间2:
区间3:
区间4:Vrated≤V<Vcut_out,P=Prated<Pmax;
区间5:V≥Vcut_out,P=0;
其中,V为测量风速,Vcut_in为切入风速,为额定转速对应风速,Vrated为额定风速,Vcut_out为切出风速;ωr为风轮转速,为最低风轮转速,为额定风轮转速,为最大风轮转速;P为有功功率,为额定转速对应的有功功率,Prated为额定有功功率,Pmax为最大有功功率;
S102:同一时间段下,基于相同采样周期获取V、ωr和P的时间序列,定义数据点(V,ωr,P),基于风机运行原理得到如下数据初筛剔除原则:
a.当0≤V<Vcut_in时,机组空转未并网,P=0;当V≥Vcut_out时,风轮叶片顺桨且风轮低速旋转或刹车,机组离网切机,P=0;当V≠0,机组发生故障或人为干预时刹车,ωr=0,机组离网切机。根据以上三种情形剔除所有超限的异常数据点;
b.当时,风机并网发电,根据风机运行特性,采用可变风速间隔εi将风速区间划分为k个间隔,则风速间隔可表示为Vi min≤Vi<Vi max(Vi min=Vcut_in+∑εi-1,Vi max=Vcut_in+∑εi, i=1,2,… k,ε0=0);
已知最优叶尖速比λopt时,风机的理论最优转速不同 Vi min≤Vi<Vi max内,有最小理论最优转速则最小理论最优出力也即(为最优风轮转矩系数,为最大风能利用系数);风机的最大理论最优转速则最大理论最优出力也即因此,Vi min≤Vi<Vi max时, 实际操作过程中,时,其中ξ′、ξ″为计及波动性的风轮转速上下限缩放系数;相应地,风机有功功率输出为其中η′、η″为计及波动性的风机有功功率上下限缩放系数。对风速V进行合理划分并在不同Vi内,初筛剔除超限异常数据点;
c.当时,有实际操作时,其中,α′、α″和κ′、κ″分别为该区间内风轮转速和有功输出的计及波动性的上下限缩放系数。在该风速区间内,可根据上述范围要求对超限数据点进行初筛剔除。
d.当Vrated≤V<Vcut_out时,P=Prated<Pmax。实际操作时,Prated-γ′≤P<Prated+γ″<Pmax,其中,β′、β″和γ′、γ″分别为该区间内风轮转速和有功输出的计及波动性的上下限缩放系数。在该风速区间内,可根据此范围对超限数据点进行初筛剔除。
进一步,步骤2)中所述四分位法为:根据采集到的V、ωr和P数据,建立三维散点图,对上述三维数据运用四分位法得到某维数据的有效取值范围[Fl,Fu]=[Q1-μIqr,Q3+μIqr],有效取值范围外的数据视为异常值进行剔除。其中,F1、Fu为数据有效取值范围的上下界限;沿单维坐标方向对该维数据值排序,然后依据数据点个数将其均匀划分为四部分,由此得到三个四分位数Q1、Q2、 Q3;Iqr=Q3-Q1是四分位距,μ为四分位系数,通常取1.5;
所述步骤2)中的参数优化为:为了保证四分位系数μ取值的合理性,采用粒子群算法对μ的取值进行寻优,寻优过程如下:
S201:在V、ωr和P组成的三维数据空间中,在V、ωr和P方向上分别设定μV、和μP,其取值范围均为[a,b],组成三维空间中的点在三维空间中随机产生m个粒子,这些粒子代表μ的取值;
S202:初始化一组μ值的位置和速度,在第t次寻优迭代的过程中,μ的位置表示为速度表示为
S203:得到一组μ值后,就可以确定各维坐标方向有效数据的取值范围;然后,对未更新前的某型号风机理论功率计算模型进行检测,将三维有效数据中的风速值带入模型,得到该模型下的理论功率,计算该模型下的理论功率与三维有效数据中实测功率值的相关性r;
定义相关性系数
其中,Pmk为第k个样本的实际功率;Ppk为第k个样本的理论恢复功率;和分别为实际功率Pmk和理论功率Ppk的期望值;N为样本个数。
S204:在每一次迭代中,μ通过跟踪两个极值来更新本身的速度和位置,一个极值是μ本身迄今为止搜索到的最优解,称为个体极值,表示为:另一个极值是到目前为止找到的最优解,称为全局极值,表示为:
在第t+1次迭代计算时,μi_labe根据以下规则来更新自己的速度和位置:
vi_labe(t+1)=αvi_labe(t)+c1r1(μi_labe(t)-xi_labe(t))+c2r2(μg_labe(t)-xi_labe(t))
xi_labe(t+1)=xi_labe(t)+vi_labe(t+1)
其中,i=1,2,…,m;下标labe可分别代表V、ωr和P;α为惯性权重,其主要是作用是为了权衡全局搜索和局部搜索的能力,α(α<0.8)较小时,其局部搜索能力较强,而α(α>1.2)较大时,其全局搜索能力较强,且总是搜寻新的区域;c1,c2是两个学习因子,通常二者均取为2;r1,r2是两个均匀分布在[0,1] 之间的随机数。
S205:μ每更新一次,即依据S203得到一个r值,当r值达到最大值或满足要求时,迭代达到最优,迭代结束。
进一步,步骤3)中对于由V、ωr和P组成的三维数据,根据机组的运行情况初始化产生kclus个聚类中心,通过计算所有数据点到聚类中心的距离对数据进行重新分类,将数据划分到离它最近的那个聚类中心所在的类中,并根据划分好的类产生新的聚类中心;如此反复进行,直到满足距离评价指标将同一个聚类中到聚类中心的距离超过该聚类中所有点到聚类中心距离平均值的n倍的数据视为异常数据,并进行删除;通常情况下, n取2.5,3.5,4.5等,根据不同的情况选取合适的n值。
进一步,所述步骤4)中:对于非连续的缺失点,直接运用缺失点两端的节点,建立三次Hermite插值多项式,进行缺失值的填补;对于连续的缺失点,采用分段三次Hermite插值法,先由缺失数据段两端的节点填补最中心的缺失值,将缺失数据分成两段再依次进行最中心缺失数据的填补。
进一步,所述步骤5)中:根据风机型号建立相应理论功率计算模型,通过累计得到风电场理论功率。风机理论功率计算模型的建立按如下步骤进行:
S501:采用Bin方法将处理过的数据按风速间隔τm/s得到NBin个区间划分;
S502:求出各个小区间内风速和功率的平均值,得到对应的点(Vi,Pi),i=1,2,…,NBin;
S503:采用最小二乘法对所有点(Vi,Pi)进行曲线拟合,得到某种型号风机的单机理论功率计算模型为
其中,Prated为常数,是风机的理论功率最大值;P为功率;V为风速;aj (j=0,1,…,6)为多项式系数。
为了保证模型的有效性,采用滑动窗口法更新数据;记窗口长度即模型所用的数据长度为L,上一模型建模点为xt-T,模型更新周期为T,则当前建模点为 xt,建模数据为从建模点xt开始至向后长度为L的数据。
为了评价理论功率恢复精度,采用相关性系数、理论功率恢复准确率和理论功率恢复均方根误差作为衡量理论功率恢复精度的指标;相关性系数的定义参考 S203中相关性系数的定义,它反映了两个变量之间的相关性,理论功率恢复准确率和理论功率恢复均方根误差能够反映理论功率恢复情况。
定义功率恢复准确率为
定义功率恢复均方根误差为
其中,Pmk为第k个样本的实际功率,Ppk为第k个样本的理论功率,N为样本个数,Pcap为风电场额定装机容量。
进一步,所述步骤6)中,对风电场的实测功率进行积分可以得到风电场实际发电量;根据各风机实测风速值和不同型号风机理论功率计算模型得到的各风机的理论功率并累加得到风电场的理论功率,然后,通过积分得到相应的风电场理论发电量;上述内容用公式表达如下
其中,Gprac为风电场实际发电量,t0为初始时刻,tinte为积分时间,Pprac(V(t)) 为风电场实际功率,Gtheo为风电场理论发电量,Ptheo(V(t))为风电场理论功率。
定义弃风率ρ为
其中,Gtheo为理论发电量,Gprac为实际发电量;
定义最优风电利用率η为
其中,G′theo为每月理论发电量,G′full为每月满发电量;
定义统计时间段内的最大极限风电利用率为
ηmax=max(η)
定义统计时间段内的最小极限风电利用率为
ηmin=min(η)
根据定义,利用上述计算得到的风电场理论发电量、实际发电量和额定装机容量等,可以分析并评价风电场的风电利用情况。
与现有技术相比,本发明具有以下有益效果:
(1)本发明的方法以风电场SCADA***记录的运行数据为基础,依据各型号风电机组的历史运行数据,对弃风限电等引发的异常数据进行剔除,对缺失数据进行填补,避免了机组检修、故障、限电等情况下对风电利用率计算带来的偏差等不利影响。如图2所示为三维原始采样数据,图3为三维原始数据经过异常数据剔除后的效果图,图4为风速、功率在相应单维坐标方向经过缺失值填补后组合成的二维散点图。通过这3幅图可以看出采用一系列数据质量控制手段进行异常数据剔除和缺失值填补后地处理效果。
(2)本发明方法基于不同型号风机的历史运行数据,对不同型号风机单独建立相应理论功率计算模型。在此基础上,恢复每台风机的理论功率并进行累加,从而得到风电场理论功率,准确度更高。图5为根据风机型号建立的单机理论功率计算模型,其数学表达是如下
其中,Prated为常数,是风机的理论功率最大值;P为功率;V为风速;aj (j=0,1,…,6)为多项式系数。
(3)本发明在风电场理论功率计算的基础上,给出了衡量风电利用率的计算指标:相关性系数、理论功率恢复准确率和理论功率恢复均方根误差,合理解决了风电利用率的评价问题。
(4)本发明的方法基于风电场SCADA***的记录数据,提出的风电场理论功率计算方法为风电利用率的计算提供了一种有效的途径,也为风电场运行情况的评估、风电与电网的规划发展等提供了科学的依据。
附图说明
图1所示为风电利用率计算方法的实施步骤示意图。
图2为采集到的单机原始风速、功率、风轮转速的三维数据。
图3为采集到的单机原始风速、功率的二维数据。
图4为原始数据经过异常数据剔除后的风速、功率、风轮转速的三维数据。
图5为原始数据经过异常数据剔除后的风速、功率二维数据效果图。
图6为风速、功率在单一方向经过缺失值填补后组合成的二维散点图。
图7为利用Bin方法得到的风电单机理论功率计算模型。
图8为采用滑动窗口法的建模数据更新机制。
具体实施方式
下面结合附图,对本发明进一步详细说明。
某一风电场拥有30台1.5MW的风电机组,采样周期为5min/点,研究数据为2015年1月1日至2015年12月31日采集到的数据。
如图1所示为风电利用率计算方法的实施步骤示意图,图2为采集到的单机原始风速、功率、风轮转速数据,图3为采集到的单机原始风速、功率数据,一种基于运行数据的风电利用率计算方法具体包括以下步骤:
步骤一:基于风机运行原理的异常数据初筛剔除,异常数据占数据总量的比例越小,采用数据预处理方法进行异常数据剔除的准确度越高。根据风机的运行机理,可对异常运行数据进行初步筛选和剔除。根据风机的运行机理,将风机的运行区间及运行特征分为如下5个区间:
区间1:0≤V<Vcut_in,P=0;
区间2:
区间3:
区间4:Vrated≤V<Vcut_out,P=Prated<Pmax;
区间5:V≥Vcut_out,P=0。
然后,分别采集相同时间段下相同采样周期的V、ωr和P的时间序列,通过风机运行原理分析得到如下数据剔除原则:
a.当0≤V<Vcut_in时,机组空转未并网,P=0;当V≥Vcut_out时,风轮叶片顺桨且风轮低速旋转或刹车,机组离网切机,P=0;当V≠0,机组发生故障或人为干预时刹车,ωr=0,机组离网切机。根据以上三种情形剔除所有超限的异常数据点;
b.当时,风机并网发电,根据风机运行特性,采用可变风速间隔εi将风速区间划分为k个间隔,则风速间隔可表示为Vi min≤Vi<Vi max(Vi min=Vcut_in+∑εi-1,Vi max=Vcut_in+∑εi, i=1,2,… k,ε0=0);
已知最优叶尖速比λopt时,风机的理论最优转速不同 Vi min≤Vi<Vi max内,有最小理论最优转速则最小理论最优出力也即(为最优风轮转矩系数,为最大风能利用系数);风机的最大理论最优转速则最大理论最优出力也即因此,Vi min≤Vi<Vi max时, 实际操作过程中,时,其中ξ′、ξ″为计及波动性的风轮转速上下限缩放系数;相应地,风机有功功率输出为其中η′、η″为计及波动性的风机有功功率上下限缩放系数。对风速V进行合理划分并在不同Vi内,初筛剔除超限异常数据点;
c.当时,有实际操作时,其中,α′、α″和κ′、κ″分别为该区间内风轮转速和有功输出的计及波动性的上下限缩放系数。在该风速区间内,可根据上述范围要求对超限数据点进行初筛剔除。
d.当Vrated≤V<Vcut_out时,P=Prated<Pmax。实际操作时,Prated-γ′≤P<Prated+γ″<Pmax,其中,β′、β″和γ′、γ″分别为该区间内风轮转速和有功输出的计及波动性的上下限缩放系数。在该风速区间内,可根据此范围对超限数据点进行初筛剔除。
综上,根据上述剔除原则对异常数据进行初步筛选和剔除,提高有效数据在总数据量中的占比。
步骤二:采用参数优化的四分位法进一步对异常数据进行剔除。
经过数据初筛,保留的数据点中有效数据比例大大增加,但是依然存在可观的异常数据点。进一步地,基于三维散点图,对三维数据进行四分位法异常数据剔除;同时,采用粒子群算法对四分位参数进行优化。
根据采集到的V、ωr和P数据,建立三维散点图,对三维数据运用四分位法得到数据的有效取值范围[Fl,Fu]=[Q1-μIqr,Q3+μIqr],并将有效取值范围外的数据视为异常值进行剔除;其中,F1、Fu为数据有效取值范围的上下界限;沿单维坐标方向对该维数据值排序,然后依据数据点个数将其均匀划分为四部分,由此得到三个四分位数Q1、Q2、Q3;Iqr=Q3-Q1是四分位距,μ为四分位系数,通常取1.5;
为了保证四分位系数μ取值的合理性,采用粒子群算法对μ的取值进行寻优,寻优过程如下:
在V、ωr和P组成的三维数据空间中,在V、ωr和P方向上分别设定μV、μωr和μP,其取值范围均为[a,b],组成三维空间中的点μ=(μV,μωr,μP),在三维空间中随机产生m个粒子,这些粒子代表μ的取值;
初始化一组μ值的位置和速度,在第t次寻优迭代的过程中,μ的位置表示为:速度表示
得到一组μ值后,就可以确定各维坐标方向有效数据的取值范围;然后,对未更新前的某型号风机理论功率计算模型进行检测,将三维有效数据中的风速值带入模型,得到该模型下的理论功率,计算该模型下的理论功率与三维有效数据中实测功率值的相关性r;
定义相关性系数
其中,Pmk为第k个样本的实际功率;Ppk为第k个样本的理论恢复功率;和分别为实际功率Pmk和理论功率Ppk的期望值;N为样本个数。
在每一次迭代中,μ通过跟踪两个极值来更新本身的速度和位置,一个极值是μ本身迄今为止搜索到的最优解,称为个体极值,表示为:另一个极值是到目前为止找到的最优解,称为全局极值,表示为:
在第t+1次迭代计算时,μi_labe根据以下规则来更新自己的速度和位置:
vi_labe(t+1)=αvi_labe(t)+c1r1(μi_labe(t)-xi_labe(t))+c2r2(μg_labe(t)-xi_labe(t))
xi_labe(t+1)=xi_labe(t)+vi_labe(t+1)
其中,i=1,2,…,m;下标labe可分别代表V、ωr和P;α为惯性权重,其主要是作用是为了权衡全局搜索和局部搜索的能力,α(α<0.8)较小时,其局部搜索能力较强,而α(α>1.2)较大时,其全局搜索能力较强,且总是搜寻新的区域;c1,c2是两个学习因子,通常二者均取为2;r1,r2是两个均匀分布在[0,1] 之间的随机数。
μ每更新一次,就可以根据前述相关性系数定义得到一个r值,当r值达到最大值或满足要求时,迭代达到最优,迭代结束。
步骤三:采用聚类算法深入剔除异常数据,对于由V、ωr和P组成的三维数据,根据机组的运行情况初始化产生4个聚类中心,通过计算所有数据点到聚类中心的距离对数据进行重新分类,将数据划分到离它最近的那个聚类中心所在的类中,并根据划分好的类产生新的聚类中心;如此反复进行,直到满足距离评价指标将同一个聚类中到聚类中心的距离超过该聚类中所有点到聚类中心距离平均值的n倍的数据视为异常数据,并进行删除;通常情况下,n可以取2.5,3.5,4.5等,可以根据不同的情况选取合适的n值。原始数据经过异常数据剔除后的效果图,如图4、5所示,图4为原始数据经过异常数据剔除后的风速、功率、风轮转速的三维数据,图5为原始数据经过异常数据剔除后的风速、功率二维数据效果图。
步骤四:采用分段Hermite插值进行缺失数据重构,对于非连续的缺失点,可以直接运用缺失点两端的节点,建立三次Hermite插值多项式,进行缺失值的填补;对于连续的缺失点,采用分段三次Hermite插值法,先由缺失数据段两端的节点填补最中心的缺失值,将缺失数据分成两段再依次进行最中心缺失数据的填补。如图6所示,为风速、功率在相应坐标维度经过缺失值填补后组合成的二维散点图。
其中,单一坐标维度填补数据所构造插值多项式的方法为:对于时间序列 X=[x1,x2,x3,…,xn]中互异的节点xk,xk+1,给定的函数值f(xi)=yi,f′(xi)=y′i, i=k,k+1。则其存在唯一的Hermite插值多项式H3(x)满足
采用基函数法,令
H3(x)=αk(x)yk+αk+1(x)yk+1+βk(x)y′k+βk+1(x)y′k+1
其中,αk(x)、αk+1(x)、βk(x)、βk+1(x)是关于节点xk及xk+1的三次Hermite 插值基函数,它们应分别满足条件
αk(xk)=1,αk(xk+1)=0,α′k(xk)=α′k(xk+1)=0;
αk+1(xk)=0,αk+1(xk+1)=1,α′k+1(xk)=α′k+1(xk+1)=0;
βk(xk)=βk(xk+1)=0,β′k(xk)=1,β′k(xk+1)=0;
βk+1(xk)=βk+1(xk+1)=0,β′k+1(xk)=0,β′k+1(xk+1)=1.
解得
步骤五:建立某型号风机理论功率计算模型及风电场理论功率恢复。如图7 所示,根据风机型号建立相应理论功率计算模型。在此基础上,恢复每台风机的理论功率,并通过累加得到风电场理论功率。
采用Bin方法将处理过的数据按风速间隔τm/s得到NBin个区间划分;
求出各个小区间内风速和功率的平均值,得到对应的点(Vi,Pi),i=1,2,…, NBin;
采用最小二乘法对通过点(Vi,Pi)的曲线进行拟合,得到某种型号风机的单机理论功率的计算模型为
其中,Prated为常数,是风机的理论功率最大值;P为功率;V为风速;aj (j=0,1,…,6)为多项式系数。
为了保证模型的有效性,需要定期对模型进行更新,这里将通过更新数据来完成对模型的更新,采用滑动窗口法更新数据;记窗口长度即模型所用的数据长度为L,上一模型建模点为xt-T,模型更新周期为T,则当前建模点为xt,建模数据为从建模点xt开始至向后长度为L的数据,更新机制如图8所示。
为了评价理论功率恢复精度,采用相关性系数、理论功率恢复准确率和理论功率恢复均方根误差作为衡量理论功率恢复精度的指标。其中,相关性系数可以反映两个变量之间的相关性,理论功率恢复准确率和理论功率恢复均方根误差能够反映理论功率恢复情况。
定义功率恢复准确率为
定义功率恢复均方根误差为
其中,Pmk为第k个样本的实际功率,Ppk为第k个样本的理论功率,N为样本个数,Pcap为风电场额定装机容量。
步骤六:计算风电利用率相关指标,对风电场的实测功率进行积分可以得到风电场实际发电量;根据各风机实测风速值和不同型号风机理论功率计算模型得到的各风机的理论功率并累加得到风电场的理论功率,然后,通过积分得到相应的风电场理论发电量;可用公式表达如下
其中,Gprac为风电场实际发电量,t0为初始时刻,tinte为积分时间,Pprac(V(t)) 为风电场实际功率,Gtheo为风电场理论发电量;Ptheo(V(t))为风电场理论功率。
定义弃风率ρ为
其中,Gtheo为理论发电量,Gprac为实际发电量。
定义最优风电利用率η为
其中,G′theo为每月理论发电量,G′full为每月满发电量;
定义统计时间段内的最大极限风电利用率为
ηmax=max(η)
定义统计时间段内的最小极限风电利用率为
ηmin=min(η)
根据定义,利用上述计算得到的风电场理论发电量、实际发电量和额定装机容量等,可以分析并评价风电场的风电利用情况。
以上对本发明的技术方案进行了详细说明。显然,本发明并不局限于所描述的实施例。基于本发明中的实施例,熟悉本技术领域的人员还可据此做出多种变化,但任何与本发明等同或相类似的变化都属于本发明保护的范围。
Claims (7)
1.一种基于运行数据的风电利用率计算方法,其特征在于,所述方法包括以下步骤:
1)基于风机运行原理对异常数据进行初步筛选剔除;
2)采用参数优化的四分位法进一步剔除异常数据;
3)采用聚类算法深入剔除异常数据;
4)采用分段Hermite插值进行缺失数据重构;
5)建立相关型号风电机组的理论功率计算模型及风电场理论功率恢复;
6)风电利用率相关指标计算。
2.根据权利要求1所述方法,其特征在于,所述步骤1)包括以下步骤:
S101:根据风机运行原理将其划分为5个运行区间:
区间1:0≤V<Vcut_in,P=0;
区间2:
区间3:
区间4:Vrated≤V<Vcut_out,P=Prated<Pmax;
区间5:V≥Vcut_out,P=0;
其中,V为测量风速,Vcut_in为切入风速,为额定风轮转速对应风速,Vrated为额定风速,Vcut_out为切出风速;ωr为风轮转速,为最低风轮转速,为额定风轮转速,为最大风轮转速;P为有功功率,为额定风轮转速对应的有功功率,Prated为额定有功功率,Pmax为最大有功功率;
所述区间1风机处于启动阶段;区间2风机处于最大风能跟踪阶段;区间3风机处于额定转速到额定功率的过渡阶段;区间4风机处于额定风速以上阶段;区间5风速过大风机处于切出阶段;
S102:同一时间段下,基于相同采样周期获取V、ωr和P的时间序列,定义数据点(V,ωr,P),基于风机运行原理得到如下数据初筛剔除原则:
a.当0≤V<Vcut_in时,机组空转未并网,P=0;当V≥Vcut_out时,风轮叶片顺桨且风轮低速旋转或刹车,机组离网切机,P=0;当V≠0,机组发生故障或人为干预时刹车,ωr=0,机组离网切机;根据以上三种情形剔除所有超限的异常数据点;
b.当时,风机并网发电,根据风机运行特性,采用可变风速间隔εi将风速区间划分为k个间隔,则风速间隔可表示为Vi min≤Vi<Vi max(Vi min=Vcut_in+∑εi-1,Vi max=Vcut_in+∑εi,i=1,2,…k,ε0=0);
已知最优叶尖速比λopt时,风机的理论最优转速不同Vi min≤Vi<Vi max内,有最小理论最优转速则最小理论最优出力也即(为最优风轮转矩系数,为最大风能利用系数);风机的最大理论最优转速则最大理论最优出力也即因此,Vi min≤Vi<Vi max时, 实际操作过程中,时,其中ξ′、ξ″为计及波动性的风轮转速上下限缩放系数;相应地,风机有功功率输出为其中η′、η″为计及波动性的风机有功功率上下限缩放系数;对风速V进行合理划分并在不同Vi内,初筛剔除超限异常数据点;
c.当时,有实际操作时,其中,α′、α″和κ′、κ″分别为该区间内风轮转速和有功输出的计及波动性的上下限缩放系数;在该风速区间内,根据上述范围要求对超限数据点进行初筛剔除;
d.当Vrated≤V<Vcut_out时,P=Prated<Pmax;实际操作时,Prated-γ′≤P<Prated+γ″<Pmax,其中,β′、β″和γ′、γ″分别为该区间内风轮转速和有功输出的计及波动性的上下限缩放系数;在该风速区间内,根据此范围对超限数据点进行初筛剔除。
3.根据权利要求1所述方法,其特征在于,步骤2)中所述四分位法为:根据采集到的V、ωr和P数据,建立三维散点图,对上述三维数据运用四分位法得到某维数据的有效取值范围[Fl,Fu]=[Q1-μIqr,Q3+μIqr],有效取值范围外的数据视为异常值进行剔除;其中,F1、Fu为数据有效取值范围的上下界限;沿单维坐标方向对该维数据值排序,然后依据数据点个数将其均匀划分为四部分,由此得到三个四分位数Q1、Q2、Q3;Iqr=Q3-Q1是四分位距,μ为四分位系数,一般可取1.5;
所述步骤2)中的参数优化为:为了保证四分位系数μ取值的合理性,采用粒子群算法对μ的取值进行寻优,寻优过程如下:
S201:在V、ωr和P组成的三维数据空间中,在V、ωr和P方向上分别设定μV、和μP,其取值范围均为[a,b],组成三维空间中的点在三维空间中随机产生m个粒子,这些粒子代表μ的取值;
S202:初始化一组μ值的位置和速度,在第t次寻优迭代的过程中,μ的位置表示为:速度表示
S203:得到一组μ值后,就可以确定各维坐标方向有效数据的取值范围;然后,对未更新前的某型号风机理论功率计算模型进行检测,将三维有效数据中的风速值带入模型,得到该模型下的理论功率,计算该模型下的理论功率与三维有效数据中实测功率值的相关性r;
定义相关性系数
<mrow>
<mi>r</mi>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>m</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>m</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>p</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>p</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<msqrt>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>m</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>m</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>p</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>p</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mfrac>
</mrow>
其中,Pmk为第k个样本的实际功率;Ppk为第k个样本的理论恢复功率;和分别为实际功率Pmk和理论功率Ppk的期望值;N为样本个数;
S204:在每一次迭代中,μ通过跟踪两个极值来更新本身的速度和位置,一个极值是μ本身迄今为止搜索到的最优解,称为个体极值,表示为:另一个极值是到目前为止找到的最优解,称为全局极值,表示为:
在第t+1次迭代计算时,μi_labe根据以下规则来更新自己的速度和位置:
vi_labe(t+1)=αvi_labe(t)+c1r1(μi_labe(t)-xi_labe(t))+c2r2(μg_labe(t)-xi_labe(t))
xi_labe(t+1)=xi_labe(t)+vi_labe(t+1)
其中,i=1,2,…,m;下标labe可分别代表V、ωr和P;α为惯性权重,其作用是为了权衡全局搜索和局部搜索的能力,α(α<0.8)较小时,其局部搜索能力较强,而α(α>1.2)较大时,其全局搜索能力较强,且总是搜寻新的区域;c1,c2是两个学习因子,二者均取为2;r1,r2是两个均匀分布在[0,1]之间的随机数;
S205:μ每更新一次,即依据S203得到一个r值,当r值达到最大值或满足要求时,迭代达到最优,迭代结束。
4.根据权利要求1所述方法,其特征在于,所述步骤3)中对于由V、ωr和P组成的三维数据,根据机组的运行情况初始化产生kclus个聚类中心,通过计算所有数据点到聚类中心的距离对数据进行重新分类,将数据划分到离它最近的那个聚类中心所在的类中,并根据划分好的类产生新的聚类中心;如此反复进行,直到满足距离评价指标将同一个聚类中到聚类中心的距离超过该聚类中所有点到聚类中心距离平均值的n倍的数据视为异常数据,并进行删除;n取2.5,3.5,4.5,根据不同的情况选取合适的n值。
5.根据权利要求1所述方法,其特征在于,所述步骤4)中:对于非连续的缺失点,直接运用缺失点两端的节点,建立三次Hermite插值多项式,进行缺失值的填补;对于连续的缺失点,采用分段三次Hermite插值法,先由缺失数据段两端的节点填补最中心的缺失值,将缺失数据分成两段再依次进行最中心缺失数据的填补。
6.根据权利要求1所述方法,其特征在于,所述步骤5)中:根据风机型号建立相应理论功率计算模型,通过累计得到风电场理论功率;风机理论功率计算模型的建立按如下步骤进行:
S501:采用Bin方法将处理过的数据按风速间隔τm/s得到NBin个区间划分;
S502:求出各个小区间内风速和功率的平均值,得到对应的点(Vi,Pi),i=1,2,…,NBin;
S503:采用最小二乘法对所有点(Vi,Pi)进行曲线拟合,得到某种型号风机的单机理论功率计算模型为
<mrow>
<mi>P</mi>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>,</mo>
<mn>0</mn>
<mo>&le;</mo>
<mi>V</mi>
<mo><</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>c</mi>
<mi>u</mi>
<mi>t</mi>
<mo>_</mo>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
<msup>
<mi>V</mi>
<mn>6</mn>
</msup>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<msup>
<mi>V</mi>
<mn>5</mn>
</msup>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<msup>
<mi>V</mi>
<mn>4</mn>
</msup>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>3</mn>
</msub>
<msup>
<mi>V</mi>
<mn>3</mn>
</msup>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>4</mn>
</msub>
<msup>
<mi>V</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>5</mn>
</msub>
<mi>V</mi>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>6</mn>
</msub>
<mo>,</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>c</mi>
<mi>u</mi>
<mi>t</mi>
<mo>_</mo>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>&le;</mo>
<mi>V</mi>
<mo><</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>r</mi>
<mi>a</mi>
<mi>t</mi>
<mi>e</mi>
<mi>d</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>r</mi>
<mi>a</mi>
<mi>t</mi>
<mi>e</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>r</mi>
<mi>a</mi>
<mi>t</mi>
<mi>e</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>&le;</mo>
<mi>V</mi>
<mo><</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>c</mi>
<mi>u</mi>
<mi>t</mi>
<mo>_</mo>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,Prated为常数,是风机的理论功率最大值;P为功率;V为风速;aj(j=0,1,…,6)为多项式系数;
采用滑动窗口法更新数据;记窗口长度即模型所用的数据长度为L,上一模型建模点为xt-T,模型更新周期为T,则当前建模点为xt,建模数据为从建模点xt开始至向后长度为L的数据;
采用相关性系数、理论功率恢复准确率和理论功率恢复均方根误差作为衡量理论功率恢复精度的指标;相关性系数的定义参考S203中相关性系数的定义,它反映了两个变量之间的相关性;理论功率恢复准确率和理论功率恢复均方根误差能够反映理论功率恢复情况;
定义功率恢复准确率为
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>a</mi>
<mi>c</mi>
<mi>c</mi>
<mi>u</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>1</mn>
<mo>-</mo>
<msqrt>
<mrow>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>m</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>p</mi>
<mi>k</mi>
</mrow>
</msub>
</mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>c</mi>
<mi>a</mi>
<mi>p</mi>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>&times;</mo>
<mn>100</mn>
<mi>%</mi>
</mrow>
定义功率恢复均方根误差为
<mrow>
<mi>N</mi>
<mi>R</mi>
<mi>M</mi>
<mi>S</mi>
<mi>E</mi>
<mo>=</mo>
<msqrt>
<mrow>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>m</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>p</mi>
<mi>k</mi>
</mrow>
</msub>
</mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>c</mi>
<mi>a</mi>
<mi>p</mi>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>&times;</mo>
<mn>100</mn>
<mi>%</mi>
</mrow>
其中,Pmk为第k个样本的实际功率,Ppk为第k个样本的理论功率,N为样本个数,Pcap为风电场额定装机容量。
7.根据权利要求1所述方法,其特征在于,所述步骤6)中:对风电场的实测功率进行积分可以得到风电场实际发电量;根据各风机实测风速值和不同型号风机理论功率计算模型得到的各风机的理论功率并累加得到风电场的理论功率,然后,通过积分得到相应的风电场理论发电量;上述内容用公式表达如下
<mrow>
<msub>
<mi>G</mi>
<mrow>
<mi>p</mi>
<mi>r</mi>
<mi>a</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>=</mo>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<mrow>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msub>
<mi>t</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
<mi>t</mi>
<mi>e</mi>
</mrow>
</msub>
</mrow>
</msubsup>
<msub>
<mi>P</mi>
<mrow>
<mi>p</mi>
<mi>r</mi>
<mi>a</mi>
<mi>c</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
<mrow>
<msub>
<mi>G</mi>
<mrow>
<mi>t</mi>
<mi>h</mi>
<mi>e</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>=</mo>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<mrow>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msub>
<mi>t</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
<mi>t</mi>
<mi>e</mi>
</mrow>
</msub>
</mrow>
</msubsup>
<msub>
<mi>P</mi>
<mrow>
<mi>t</mi>
<mi>h</mi>
<mi>e</mi>
<mi>o</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
其中,Gprac为风电场实际发电量,t0为初始时刻,tinte为积分时间,Pprac(V(t))为风电场实际功率,Gtheo为风电场理论发电量,Ptheo(V(t))为风电场理论功率。
定义弃风率ρ为
<mrow>
<mi>&rho;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>G</mi>
<mrow>
<mi>t</mi>
<mi>h</mi>
<mi>e</mi>
<mi>o</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>p</mi>
<mi>r</mi>
<mi>a</mi>
<mi>c</mi>
</mrow>
</msub>
</mrow>
<msub>
<mi>G</mi>
<mrow>
<mi>t</mi>
<mi>h</mi>
<mi>e</mi>
<mi>o</mi>
</mrow>
</msub>
</mfrac>
</mrow>
其中,Gtheo为理论发电量,Gprac为实际发电量;
定义最优风电利用率η为
<mrow>
<mi>&eta;</mi>
<mo>=</mo>
<mfrac>
<msubsup>
<mi>G</mi>
<mrow>
<mi>t</mi>
<mi>h</mi>
<mi>e</mi>
<mi>o</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
<msubsup>
<mi>G</mi>
<mrow>
<mi>f</mi>
<mi>u</mi>
<mi>l</mi>
<mi>l</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
</mfrac>
</mrow>
其中,G′theo为每月理论发电量,G′full为每月满发电量;
定义统计时间段内的最大极限风电利用率为
ηmax=max(η)
定义统计时间段内的最小极限风电利用率为
ηmin=min(η)
根据定义,利用上述计算得到的风电场理论发电量、实际发电量和额定装机容量,分析并评价风电场的风电利用情况。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710514679.XA CN107330183B (zh) | 2017-06-29 | 2017-06-29 | 一种基于运行数据的风电利用率计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710514679.XA CN107330183B (zh) | 2017-06-29 | 2017-06-29 | 一种基于运行数据的风电利用率计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107330183A true CN107330183A (zh) | 2017-11-07 |
CN107330183B CN107330183B (zh) | 2020-04-17 |
Family
ID=60197325
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710514679.XA Active CN107330183B (zh) | 2017-06-29 | 2017-06-29 | 一种基于运行数据的风电利用率计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107330183B (zh) |
Cited By (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108052963A (zh) * | 2017-12-01 | 2018-05-18 | 北京金风慧能技术有限公司 | 风电功率预测建模的数据筛选方法、装置及风力发电机组 |
CN108448625A (zh) * | 2018-04-09 | 2018-08-24 | 中能电力科技开发有限公司 | 一种基于数据驱动的风机切入风速计算方法 |
CN108536958A (zh) * | 2018-04-09 | 2018-09-14 | 中能电力科技开发有限公司 | 一种基于功率曲线健康状态分级的风机实时评价方法 |
CN108695845A (zh) * | 2018-05-28 | 2018-10-23 | 华润电力风能(汕头潮南)有限公司 | 一种风电机组可利用率的计算方法及*** |
CN108710689A (zh) * | 2018-05-22 | 2018-10-26 | 中国华能集团清洁能源技术研究院有限公司 | 一种风电场气象数据管理***及方法 |
CN108763584A (zh) * | 2018-06-11 | 2018-11-06 | 北京天泽智云科技有限公司 | 一种风功率曲线散点过滤的方法及其*** |
CN108847682A (zh) * | 2018-06-15 | 2018-11-20 | 中国电力科学研究院有限公司 | 一种基于运行数据的光伏发电站并网性能评估方法及*** |
CN108932580A (zh) * | 2018-06-05 | 2018-12-04 | 浙江运达风电股份有限公司 | 基于数据建模的风电机组变桨轴承磨损监测和预警方法 |
CN108985490A (zh) * | 2018-06-14 | 2018-12-11 | 内蒙古大学 | 基于滑动窗口特征数据拟合的风力发电机功率曲线修正方法 |
CN109871661A (zh) * | 2019-03-26 | 2019-06-11 | 西安交通大学 | 一种叶端定时信号丢失的识别补全方法 |
CN110206681A (zh) * | 2018-02-28 | 2019-09-06 | 北京金风科创风电设备有限公司 | 自适应地调整风力发电机组的转速的方法和装置 |
CN111027835A (zh) * | 2019-12-02 | 2020-04-17 | 国网宁夏电力有限公司经济技术研究院 | 电网规划决策数据噪音识别与降噪方法及*** |
CN111209914A (zh) * | 2019-12-18 | 2020-05-29 | 中电投电力工程有限公司 | 一种异常风功率数据的剔除方法 |
CN111275570A (zh) * | 2020-01-13 | 2020-06-12 | 浙江大学 | 一种基于迭代式统计与假设检验的风电机组功率异常值检测方法 |
CN111522808A (zh) * | 2020-04-29 | 2020-08-11 | 贵州电网有限责任公司 | 一种风电机组异常运行数据处理方法 |
CN111860956A (zh) * | 2020-06-18 | 2020-10-30 | 浙江运达风电股份有限公司 | 基于发电量可利用率的风电机组损失电量在线计算方法 |
CN111881617A (zh) * | 2020-07-02 | 2020-11-03 | 上海电气风电集团股份有限公司 | 数据处理方法、风力发电机组的性能评估方法及*** |
CN112085115A (zh) * | 2020-09-16 | 2020-12-15 | 华北电力大学 | 高维空间异常数据优化识别方法 |
CN112234655A (zh) * | 2020-07-30 | 2021-01-15 | 国网甘肃省电力公司电力科学研究院 | 一种送端电网安全稳定的源网协调调峰优化方法 |
CN112347655A (zh) * | 2020-11-17 | 2021-02-09 | 国网青海省电力公司 | 一种基于机组运行性能评估的风电场理论功率计算方法 |
WO2021027011A1 (zh) * | 2019-08-14 | 2021-02-18 | 北京天泽智云科技有限公司 | 提高风电***数据质量的方法及装置 |
CN112784215A (zh) * | 2021-01-21 | 2021-05-11 | 中国三峡新能源(集团)股份有限公司 | 基于实时功率遥测量的风电场理论发电量计算方法 |
CN113095602A (zh) * | 2021-05-27 | 2021-07-09 | 中国水利水电科学研究院 | 一种以出力平稳为目标的流域风电场群规划方法 |
CN113569399A (zh) * | 2021-07-20 | 2021-10-29 | 华能四平风力发电有限公司 | 风电机组运行数据处理方法、***、设备及存储介质 |
CN114033631A (zh) * | 2021-11-08 | 2022-02-11 | 浙江运达风电股份有限公司 | 一种风电机组风能利用系数在线辨识方法 |
CN115983663A (zh) * | 2022-11-04 | 2023-04-18 | 广西电网有限责任公司 | 基于新能源的发电分析方法及相关设备 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104252649A (zh) * | 2014-09-25 | 2014-12-31 | 东南大学 | 一种计及多风电场相关性的区域风电出力预测方法 |
CN105160060A (zh) * | 2015-07-17 | 2015-12-16 | 中国电力科学研究院 | 一种基于实际功率曲线拟合的风电场理论功率确定方法 |
US20160025071A1 (en) * | 2014-07-28 | 2016-01-28 | State Grid Corporation Of China | Method of computing theoretical power of wind farm based on sample wind turbine method |
CN105930933A (zh) * | 2016-04-26 | 2016-09-07 | 华北电力科学研究院有限责任公司 | 风电场理论功率曲线确定方法及装置 |
CN105986961A (zh) * | 2016-04-28 | 2016-10-05 | 华北电力大学 | 一种变速变桨风力机功率优化控制方法 |
CN106894950A (zh) * | 2017-03-09 | 2017-06-27 | 华电电力科学研究院 | 一种基于风电机组传动特性的功率特性验证方法 |
-
2017
- 2017-06-29 CN CN201710514679.XA patent/CN107330183B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160025071A1 (en) * | 2014-07-28 | 2016-01-28 | State Grid Corporation Of China | Method of computing theoretical power of wind farm based on sample wind turbine method |
CN104252649A (zh) * | 2014-09-25 | 2014-12-31 | 东南大学 | 一种计及多风电场相关性的区域风电出力预测方法 |
CN105160060A (zh) * | 2015-07-17 | 2015-12-16 | 中国电力科学研究院 | 一种基于实际功率曲线拟合的风电场理论功率确定方法 |
CN105930933A (zh) * | 2016-04-26 | 2016-09-07 | 华北电力科学研究院有限责任公司 | 风电场理论功率曲线确定方法及装置 |
CN105986961A (zh) * | 2016-04-28 | 2016-10-05 | 华北电力大学 | 一种变速变桨风力机功率优化控制方法 |
CN106894950A (zh) * | 2017-03-09 | 2017-06-27 | 华电电力科学研究院 | 一种基于风电机组传动特性的功率特性验证方法 |
Non-Patent Citations (2)
Title |
---|
朱倩雯等: "《风电场输出功率异常数据识别与重构方法研究》", 《电力***保护与控制》 * |
赵永宁等: "《风电场弃风异常数据簇的特征及处理方法》", 《电力***自动化》 * |
Cited By (36)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108052963A (zh) * | 2017-12-01 | 2018-05-18 | 北京金风慧能技术有限公司 | 风电功率预测建模的数据筛选方法、装置及风力发电机组 |
CN110206681B (zh) * | 2018-02-28 | 2020-05-05 | 北京金风科创风电设备有限公司 | 自适应地调整风力发电机组的转速的方法和装置 |
CN110206681A (zh) * | 2018-02-28 | 2019-09-06 | 北京金风科创风电设备有限公司 | 自适应地调整风力发电机组的转速的方法和装置 |
CN108448625A (zh) * | 2018-04-09 | 2018-08-24 | 中能电力科技开发有限公司 | 一种基于数据驱动的风机切入风速计算方法 |
CN108536958A (zh) * | 2018-04-09 | 2018-09-14 | 中能电力科技开发有限公司 | 一种基于功率曲线健康状态分级的风机实时评价方法 |
CN108448625B (zh) * | 2018-04-09 | 2020-09-11 | 中能电力科技开发有限公司 | 一种基于数据驱动的风机切入风速计算方法 |
CN108536958B (zh) * | 2018-04-09 | 2021-11-05 | 中能电力科技开发有限公司 | 一种基于功率曲线健康状态分级的风机实时评价方法 |
CN108710689A (zh) * | 2018-05-22 | 2018-10-26 | 中国华能集团清洁能源技术研究院有限公司 | 一种风电场气象数据管理***及方法 |
CN108695845A (zh) * | 2018-05-28 | 2018-10-23 | 华润电力风能(汕头潮南)有限公司 | 一种风电机组可利用率的计算方法及*** |
CN108932580A (zh) * | 2018-06-05 | 2018-12-04 | 浙江运达风电股份有限公司 | 基于数据建模的风电机组变桨轴承磨损监测和预警方法 |
CN108763584A (zh) * | 2018-06-11 | 2018-11-06 | 北京天泽智云科技有限公司 | 一种风功率曲线散点过滤的方法及其*** |
CN108763584B (zh) * | 2018-06-11 | 2021-11-02 | 北京天泽智云科技有限公司 | 一种风功率曲线散点过滤的方法及其*** |
CN108985490A (zh) * | 2018-06-14 | 2018-12-11 | 内蒙古大学 | 基于滑动窗口特征数据拟合的风力发电机功率曲线修正方法 |
CN108847682B (zh) * | 2018-06-15 | 2022-10-04 | 中国电力科学研究院有限公司 | 一种基于运行数据的光伏发电站并网性能评估方法及*** |
CN108847682A (zh) * | 2018-06-15 | 2018-11-20 | 中国电力科学研究院有限公司 | 一种基于运行数据的光伏发电站并网性能评估方法及*** |
CN109871661A (zh) * | 2019-03-26 | 2019-06-11 | 西安交通大学 | 一种叶端定时信号丢失的识别补全方法 |
CN109871661B (zh) * | 2019-03-26 | 2020-12-25 | 西安交通大学 | 一种叶端定时信号丢失的识别补全方法 |
WO2021027011A1 (zh) * | 2019-08-14 | 2021-02-18 | 北京天泽智云科技有限公司 | 提高风电***数据质量的方法及装置 |
CN111027835A (zh) * | 2019-12-02 | 2020-04-17 | 国网宁夏电力有限公司经济技术研究院 | 电网规划决策数据噪音识别与降噪方法及*** |
CN111209914A (zh) * | 2019-12-18 | 2020-05-29 | 中电投电力工程有限公司 | 一种异常风功率数据的剔除方法 |
CN111275570A (zh) * | 2020-01-13 | 2020-06-12 | 浙江大学 | 一种基于迭代式统计与假设检验的风电机组功率异常值检测方法 |
CN111522808A (zh) * | 2020-04-29 | 2020-08-11 | 贵州电网有限责任公司 | 一种风电机组异常运行数据处理方法 |
CN111522808B (zh) * | 2020-04-29 | 2023-07-28 | 贵州电网有限责任公司 | 一种风电机组异常运行数据处理方法 |
CN111860956A (zh) * | 2020-06-18 | 2020-10-30 | 浙江运达风电股份有限公司 | 基于发电量可利用率的风电机组损失电量在线计算方法 |
CN111881617A (zh) * | 2020-07-02 | 2020-11-03 | 上海电气风电集团股份有限公司 | 数据处理方法、风力发电机组的性能评估方法及*** |
CN111881617B (zh) * | 2020-07-02 | 2024-03-26 | 上海电气风电集团股份有限公司 | 数据处理方法、风力发电机组的性能评估方法及*** |
CN112234655A (zh) * | 2020-07-30 | 2021-01-15 | 国网甘肃省电力公司电力科学研究院 | 一种送端电网安全稳定的源网协调调峰优化方法 |
CN112085115A (zh) * | 2020-09-16 | 2020-12-15 | 华北电力大学 | 高维空间异常数据优化识别方法 |
CN112347655A (zh) * | 2020-11-17 | 2021-02-09 | 国网青海省电力公司 | 一种基于机组运行性能评估的风电场理论功率计算方法 |
CN112784215B (zh) * | 2021-01-21 | 2022-04-15 | 中国三峡新能源(集团)股份有限公司 | 基于实时功率遥测量的风电场理论发电量计算方法 |
CN112784215A (zh) * | 2021-01-21 | 2021-05-11 | 中国三峡新能源(集团)股份有限公司 | 基于实时功率遥测量的风电场理论发电量计算方法 |
CN113095602A (zh) * | 2021-05-27 | 2021-07-09 | 中国水利水电科学研究院 | 一种以出力平稳为目标的流域风电场群规划方法 |
CN113569399A (zh) * | 2021-07-20 | 2021-10-29 | 华能四平风力发电有限公司 | 风电机组运行数据处理方法、***、设备及存储介质 |
CN114033631A (zh) * | 2021-11-08 | 2022-02-11 | 浙江运达风电股份有限公司 | 一种风电机组风能利用系数在线辨识方法 |
CN115983663A (zh) * | 2022-11-04 | 2023-04-18 | 广西电网有限责任公司 | 基于新能源的发电分析方法及相关设备 |
CN115983663B (zh) * | 2022-11-04 | 2023-09-12 | 广西电网有限责任公司 | 基于新能源的发电分析方法及相关设备 |
Also Published As
Publication number | Publication date |
---|---|
CN107330183B (zh) | 2020-04-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107330183A (zh) | 一种基于运行数据的风电利用率计算方法 | |
Flores et al. | Application of a control algorithm for wind speed prediction and active power generation | |
CN103390202B (zh) | 一种基于相似数据选取的光伏电站输出功率预测方法 | |
Shi et al. | Hybrid forecasting model for very-short term wind power forecasting based on grey relational analysis and wind speed distribution features | |
CN102663251B (zh) | 基于计算流体力学模型的风电场功率物理预测方法 | |
CN110880789B (zh) | 一种风电、光伏联合发电***经济调度方法 | |
CN103683274B (zh) | 区域中长期风电发电量概率预测方法 | |
CN105160060B (zh) | 一种基于实际功率曲线拟合的风电场理论功率确定方法 | |
CN102184337B (zh) | 一种受气象信息影响的新能源发电能力动态组合分析方法 | |
CN105046374A (zh) | 一种基于核极限学习机模型的功率区间预测方法 | |
CN106875033A (zh) | 一种基于动态自适应的风电集群功率预测方法 | |
CN102419394B (zh) | 一种预测分辨率可变的风光功率预测方法 | |
CN106351793A (zh) | 用于改进风力发电的***和方法 | |
CN106505631B (zh) | 智能风电风功率预测*** | |
CN102663513A (zh) | 利用灰色关联度分析的风电场功率组合预测建模方法 | |
Verma et al. | Markov models based short term forecasting of wind speed for estimating day-ahead wind power | |
CN103138256A (zh) | 一种新能源电力消纳全景分析***及方法 | |
CN104978608A (zh) | 一种风电功率预测装置及预测方法 | |
CN107909211A (zh) | 基于模糊c均值聚类算法的风场等值建模及优化控制方法 | |
CN101794996A (zh) | 风电场出力实时预测方法 | |
CN103902837A (zh) | 一种基于经验Copula函数进行风速预测的方法 | |
CN103473621A (zh) | 风电场短期功率预测方法 | |
Shetty et al. | Comparison of modeling methods for wind power prediction: a critical study | |
CN108985490A (zh) | 基于滑动窗口特征数据拟合的风力发电机功率曲线修正方法 | |
CN104657584A (zh) | 一种基于Lorenz***的风速预测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |