CN112347655A - 一种基于机组运行性能评估的风电场理论功率计算方法 - Google Patents
一种基于机组运行性能评估的风电场理论功率计算方法 Download PDFInfo
- Publication number
- CN112347655A CN112347655A CN202011290260.9A CN202011290260A CN112347655A CN 112347655 A CN112347655 A CN 112347655A CN 202011290260 A CN202011290260 A CN 202011290260A CN 112347655 A CN112347655 A CN 112347655A
- Authority
- CN
- China
- Prior art keywords
- data
- power
- time
- wind turbine
- abnormal 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
Images
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/06—Wind turbines or wind farms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
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)
Abstract
一种基于机组运行性能评估的风电场理论功率计算方法,涉及新能源发电技术领域,其结构为:数据识别阶段,对在数据采集或通信等环节发生故障或受到干扰导致的异常数据进行识别;数据修复阶段,针对数据识别阶段识别出来异常数据的特点,分三个类别对新能源运行数据实现异常数据重构;非样板机分组及各组出力比例系数识别阶段,动态识别每组的出力比例系数;风电场的理论功率计算阶段,根据式统计时间内理论发电功率的计算式计算风电场的理论功率。本发明的有益效果在于:易操作,建模简单,并且能够准确识别各组非样板机和样板机之间的出力比例系数,从而有效提高风电场理论功率的计算精度。
Description
技术领域
本发明涉及新能源发电技术领域,具体涉及一种基于机组运行性能评估的风电场理论功率计算方法。
背景技术
随着风力发电技术的逐渐成熟,风电可再生能源的关注度大大提高,然而由于风电场大规模的快速无序投运以及电网架构建设的相对滞后,造成电网对风电接纳能力受限,使得弃风现象时有发生。因此开展风电场理论发电功率的研究十分必要,它能够对风电场的弃风电量进行科学、准确的评估,缓解网厂矛盾,有利于调度部门评估风电场的实际出力水平,为科学统计风电场限电电量提供参考。
目前,风电场的理论功率计算主要采用样板机法,即选择性能良好、运行稳定和故障相对较少的单机设备作为样板机,通过其实际功率推算出风电场的理论功率,但是这种方法由于默认样板机组和非样板机组的出力比例系数为1,而未考虑受地理位置、排布方式以及气象等因素所造成的样板机和非样板机出力比例系数的不确定性,使得原有样板机法计算的理论功率会出现较大误差。
发明内容
由于地理位置、排布方式以及气象等因素所造成的样板机和非样板机出力比例系数的不确定性,使得原有样板机法计算的理论功率会出现较大误差,为了克服现有技术存在的缺陷,本发明提出了一种基于机组运行性能评估的风电场理论功率计算方法。
本发明提供一种基于机组运行性能评估的风电场理论功率计算方法,具体包括:数据识别阶段、数据修复阶段、非样板机分组及各组出力比例系数识别阶段和风电场的理论功率计算阶段。
数据识别阶段。对在数据采集或通信等环节发生故障或受到干扰导致的异常进行识别,步骤如下:
步骤1:根据空缺值识别准则,对所有新能源风电机组时序功率数据中某个时间点或时间段出现的空值进行识别,并将该异常数据剔除。
步骤2:根据越限值识别准则,对所有新能源风电机组时序功率数据中超出合理范围的越限值进行识别,并将该异常数据剔除。
步骤3:根据不刷新值识别准则,对所有新能源风电机组时序功率数据中连续一段时间不刷新的数据值进行识别,并将该异常数据剔除。
步骤4:根据畸变数据识别准则,对所有新能源风电机组时序功率数据中产生的畸变数据进行识别,并将该异常数据剔除。
数据修复阶段。针对步骤1-4识别出来异常数据的特点,分三个类别,包括单点异常、多点异常和连续异常,对新能源运行数据实现异常数据重构,步骤如下:
步骤5:对于新能源风电机组时序功率数据中的单点异常数据,采用常量/均值的方法对其进行填充。
步骤6:对于新能源风电机组时序功率数据中不超4h的多点异常数据,采用基于ARMA模型预测的重构方法对异常数据进行重构,并建立平均可信度γp计算加权重构值时刻k正向预测值和反向预测值从而提高k时刻重构结果的可靠性;对于超过4h的多点异常则作为连续异常数据,数据重构过程转至步骤7。
步骤7:对于新能源风电机组时序功率数据中的连续异常数据,采用基于序列时延相关性的重构方法首先找出与新能源风电机组Wy时延相关性最大的风电机组Wx(x=1,2,...,N,N为新能源风电机组的个数),利用Wx的出力曲线对Wy的连续异常运行数据进行数据重构。
非样板机分组及各组出力比例系数识别阶段。该阶段基于机组运行性能对新能源电站中同型号的非样板机进行分组后,再动态识别每组的出力比例系数,从而更准确的掌握非样板机的相对出力情况,步骤如下:
步骤8:通过对风电场中第i台非样板机的基本出力性能特征向量Wi和样板机的基本出力性能平均特征向量的Pearson相关系数Ris以0.1为一个区间单位进行划分,并将每个区间单位的非样板机划分为一组,记为第g组,其中非样板机数量为共计20组。
步骤9:通过计算k型号风机第g组非样板机在统计时间内的实际出力平均值序列与样板机在统计时间内的实际出力平均值序列的功率差率得出第g组非样板机相对于样板机出力的比例系数为以此动态识别各组的出力比例系数。
风电场的理论功率计算阶段。
本发明的有益效果在于:考虑到新能源风电机组在数据采集或通信等环节发生故障或受到干扰而导致的数据异常,建立了不同类型的异常数据识别准则对新能源风电机组的时序出力数据进行异常数据识别与剔除;考虑到剔除异常数据后会破坏数据的完整性以及影响数据的可用性,采用不同的方法分别对单点、多点和连续异常数据进行重构,从而有效减小后续对新能源电站理论功率计算带来的误差;考虑到原有样板机法在计算理论功率时未能准确计算样板机和非样板机之间的出力比例系数,基于风电机组运行性能评估的非样板机分组方法,相较于考虑地形和尾流效应的分组方法,具有易操作,建模简单,并且能够准确识别各组非样板机和样板机之间的出力比例系数,从而有效提高风电场理论功率的计算精度。
附图说明
图1为本发明风电场的理论功率计算流程框图。
具体实施方式
实施例1,如图1所示,本发明提供一种基于机组运行性能评估的风电场理论功率计算方法,其包括四个阶段,数据识别阶段、数据修复阶段、非样板机分组及各组出力比例系数识别阶段和风电场的理论功率计算阶段;包括以下步骤如下:
步骤1:根据空缺值识别准则,对所有新能源风电机组时序功率数据中某个时间点或时间段出现的空缺值进行识别,并将该异常数据剔除;
步骤2:根据越限值识别准则,对所有新能源风电机组时序功率数据中超出合理范围的越限值进行识别,并将该异常数据剔除;
步骤3:根据不刷新值识别准则,对所有新能源风电机组时序功率数据中连续一段时间不刷新的数据值进行识别,并将该异常数据剔除;
步骤4:根据畸变数据识别准则,对所有新能源风电机组时序功率数据中产生的畸变数据进行识别,并将该异常数据剔除;
步骤5:对于新能源风电机组时序功率数据中的单点异常数据,采用常量/均值的方法对其进行填充;
步骤6:对于新能源风电机组时序功率数据中不超4h的多点异常数据,采用基于ARMA模型预测的重构方法对异常数据进行重构,并建立平均可信度γp计算加权重构值时刻k正向预测值和反向预测值从而提高k时刻重构结果的可靠性;对于超过4h的多点异常则作为连续异常数据,数据重构过程转至步骤7;
步骤7:对于新能源风电机组时序功率数据中的连续异常数据,采用基于序列时延相关性的重构方法首先找出与新能源风电机组Wy时延相关性最大的风电机组Wx(x=1,2,...,N,N为新能源风电机组的个数),利用Wx的出力曲线对Wy的连续异常运行数据进行数据重构;
步骤8:通过对风电场中第i台非样板机的基本出力性能特征向量Wi和样板机的基本出力性能平均特征向量的Pearson相关系数Ris以0.1为一个区间单位进行划分,并将每个区间单位的非样板机划分为一组,记为第g组,其中非样板机数量为共计20组;
步骤9:通过计算k型号风机第g组非样板机在统计时间内的实际出力平均值序列与样板机在统计时间内的实际出力平均值序列的功率差率得出第g组非样板机相对于样板机出力的比例系数为以此动态识别各组的出力比例系数;
对风电场的理论功率改进计算具体包括以下步骤:
对所有新能源风电机组的时序出力数据X={x1,x2,...,xn}进行异常数据识别及重构。
不同类型异常数据的识别准则如下:
(1)空缺值的识别
若存在xi满足:
则判断xi存在数据丢失。
(2)越限值的识别
若存在xi满足下式关系,则判断xi超出正常数值范围。
xi>xh|xi<xl (2)
式中,xh,xl分别为正常数值范围上限和下限。
(3)不刷新值的识别
若数据连续一段时间不刷新,则除去第一个数据点外均辨识为不刷新值。
(4)畸变数据的识别
相邻风电机组的功率通常具有相似趋势,故基于相邻风电机组出力空间相关性,利用N个与风电机组A具有空间相关性的相邻风电机组的功率数据XM={xM1,xM2,...,xMn}(M=1,2,…,N),对观测风电机组的功率数据XA={xA1,xA2,...,xAn}进行校验,识别异常数据。具体如下:基于式(3)对功率进行归一化,得到观测风电机组的功率归一化数据XA *={xA1 *,xA2 *,...,xAn *}和相邻风电机组的功率归一化数据XM *={xM1 *,xM2 *,...,xMn *},二者在每个采样点的差值为式(4)。
eAi-Mi=|xAi *-xMi *| (4)
风电机组A和M的功率差值和平均值分别为:
EA-M={eA1-M1,eA2-M2,…,eAn-Mn} (5)
因为不同风电机组、不同时间段获得的功率数据异常情况不同,为得到较高可信度的识别结果,对于已剔除数据缺失、不刷新以及数据越限的待识别数据集,需计算观测风电机组与相邻风电机组归一化后输出功率差的平均值和标准差。由于不能事先确定相邻风电机组出力数据的正确性,所以令根据莱以特准则判断异常数据,即当时,则认为xi为畸变数据。
不同异常数据的重构方法如下:
(1)单点异常
对于单点异常数据可以采用填充法对异常数据进行工程化重构,即选择一个经验常量或者按照某种规则对异常数据进行替换,该过程简单省时,且具有可观的精确性。
1)常量填充:基于人工经验,对单点异常数据采用同一个指定的常量值进行填充。
2)均值填充:可采用数据序列的均值或异常数据前后数据的均值对单点异常数据进行填充。
(2)多点异常
对于多点异常,由于异常数据两侧均存在完好的实际运行数据,因此可以根据两侧实际运行数据,采用基于ARMA模型预测的重构方法对异常数据进行重构。
Pt=β1Pt-1+β2Pt-2+…+βpPt-p+Zt (7)
Zt=εt+λ1εt-1+λ2εt-2+…+λ1εt-1 (8)
式中,βi为自回归参数;λi为滑动平均参数;εi为预测误差。
2)对于第p步预测结果,如果预测误差小于误差阈值ε0时,则认为该预测结果可信。通过对大量完好运行数据进行预测,确定数据序列第p步预测的平均可信度γp:
式中,Count(·)为次数统计;Ep为第p步预测误差;N为总次数。
基于自身出力规律的重构方法考虑到ARMA的预测效果,一般重构时间长度限制在4h以内。
(3)连续异常
随着异常数据时间尺度增加,对于基于ARMA模型预测的重构方法,其正确率将会显著降低。由于各类新能源运行数据间存在一定的相关性,因此对于连续异常数据,可以采用基于序列时延相关性的重构方法对异常数据进行重构。
1)时延相关性
假设有两个时间序列X={x1,x2,…,xn}和Y={y1,y2,…,yn},Y相对于X延迟l时刻的相关系数R(l)计算公式为:
取最大值时对应的延迟为(最大时延相关点),如果(为相关性阈值),则X和Y具有时延相关性。(可采用三点预测探查法,从最有可能出现最大时延相关点的地方设置第1个探查点(采用三角法确定),在该探查点左右两侧分别以“几何渐进”的方式设置其他探查点。
2)基于新能源风电机组出力时延相关性的重构
假设缺失数据时间段为[tm,tn],待重构风电机组为W0,其余风电机组为W1,W2,…。取各风电机组在时段[tm-t0,tm]内的出力数据。首先计算W0与W1在该时段内出力的最大时延相关点和相应的时延相关系数,若R(l1')=max(R(l1'),R(l2'),…),则W0与W1相对应的最大时延相关点为l1',,相应的时延相关系数为R(l1');之后再按此分别计算W0与W2,W3…在该时段内出力的最大时延相关点和相应的时延相关系数;若有W0与W1的时延相关性系数最大,则利用W1的出力曲线对W0的历史运行数据进行重构。即W0和W1的出力关系可用线性回归拟合为:
对风电场非样板机的分组及各组出力比例系数的动态识别和理论功率的计算如下:
(1)基于机组运行性能评估的非样板机分组
首先,建立风电机组的基本发电性能指标(非限电情况),如下:
a)月最大功率w1
月最大功率w1(kW)为观察日前30天风电机组实发功率Pr={p1,p2,…,pn}中的最大值,即:
w1=max{pi},i=1,2,…,n (13)
式中,n为样本采样点数。
b)月平均发电功率w2
式中,Pr,j,t为观察日前30天第j日第t个采样点的实发功率(kW)。
c)平均发电功率年变化率w3
机组平均发电功率年变化率反映机组发电水平随时间变化的趋势:
式中,Pr,av,y为机组观察日前第一个365天内的平均实际发电功率(kW);P'r,av,y为机组观察日前第二个365天内的平均实际发电功率(kW);P'r,j,t为观察日前第一个365天内第j日第t个采样点的实发功率(kW);P”r,j,t为观察日前第二个365天内第j日第t个采样点的实发功率(kW);M为每日采样点个数。
d)月满发时数w4
月满发时数w4(h)为观察日前30天风电机组实发电量与额定功率之比,可用于不同装机容量发电***的比较,其计算式如下:
式中,Pn为观察日前30天风电机组额定功率(kW),Er为观察日前30天的实发电量。
利用风电机组基本发电性能指标,建立统计时间内风电机组基本出力性能特征向量W={w1,w2,w3,w4}。设k型号风机第i台非样板机在统计时间内的基本出力性能特征向量为Wi k={wi1,wi2,wi3,wi4},k型号样板机在统计时间内的基本出力性能平均特征向量为其中,Mk为k型号样板机全风电场的总数量;为k型号风机第m台样板机在统计时间内的基本出力性能特征向量。
设以0.1为一个区间单位对Ris进行划分,通常Ris∈[0.1×(g-11),0.1×(g-10)),g为整数,且g∈[1,20]。将每个区间单位的非样板机划分为一组,设第g组非样板机数量为(为了提高计算精度,可以减少区间单位的长度)。
(2)各组出力比例系数的动态识别
设第g组非样板机在统计时间内的实际出力平均值序列为:
设k型号风机样板机在统计时间内的实际出力平均值序列为:
式中,n为在统计时间内的采样点个数。
因此,第g组非样板机相对于样板机出力的比例系数可取为理想情况下,当风电场内所有机组的理论出力都相等时,此时g=10,非样板机相对于样板机出力的比例系数均为1,电站的理论功率计算公式与目前理想的计算公式完全相符。
(3)风电场的理论功率计算如下:
利用样板机的实际出力,分别计算各组非样板机的理论发电功率,最终得到新能源电站在统计时间内理论发电功率的计算式:
Claims (5)
1.一种基于机组运行性能评估的风电场理论功率计算方法,其特征在于:包括四个阶段,数据识别阶段、数据修复阶段、非样板机分组及各组出力比例系数识别阶段和风电场的理论功率计算阶段;
数据识别阶段,对在数据采集或通信等环节发生故障或受到干扰导致的异常数据进行识别,步骤如下:
步骤1:根据空缺值识别准则,对所有新能源风电机组时序功率数据中某个时间点或时间段出现的空缺值进行识别,并将该异常数据剔除;
步骤2:根据越限值识别准则,对所有新能源风电机组时序功率数据中超出合理范围的越限值进行识别,并将该异常数据剔除;
步骤3:根据不刷新值识别准则,对所有新能源风电机组时序功率数据中连续一段时间不刷新的数据值进行识别,并将该异常数据剔除;
步骤4:根据畸变数据识别准则,对所有新能源风电机组时序功率数据中产生的畸变数据进行识别,并将该异常数据剔除;
数据修复阶段,针对步骤1-4识别出来异常数据的特点,分三个类别,包括单点异常、多点异常和连续异常,对新能源运行数据实现异常数据重构,步骤如下:
步骤5:对于新能源风电机组时序功率数据中的单点异常数据,采用常量/均值的方法对其进行填充;
步骤6:对于新能源风电机组时序功率数据中不超4h的多点异常数据,采用基于ARMA模型预测的重构方法对异常数据进行重构,并建立平均可信度γp计算加权重构值时刻k正向预测值和反向预测值从而提高k时刻重构结果的可靠性;对于超过4h的多点异常则作为连续异常数据,数据重构过程转至步骤7;
步骤7:对于新能源风电机组时序功率数据中的连续异常数据,采用基于序列时延相关性的重构方法首先找出与新能源风电机组Wy时延相关性最大的风电机组Wx,x=1,2,...,N,N为新能源风电机组的个数,利用Wx的出力曲线对Wy的连续异常运行数据进行数据重构;
非样板机分组及各组出力比例系数识别阶段,基于机组运行性能对新能源电站中同型号的非样板机进行分组后,再动态识别每组的出力比例系数,从而更准确的掌握非样板机的相对出力情况,步骤如下:
步骤8:通过对风电场中第i台非样板机的基本出力性能特征向量Wi和样板机的基本出力性能平均特征向量的Pearson相关系数Ris以0.1为一个区间单位进行划分,并将每个区间单位的非样板机划分为一组,记为第g组,其中非样板机数量为共计20组;
步骤9:通过计算k型号风机第g组非样板机在统计时间内的实际出力平均值序列与样板机在统计时间内的实际出力平均值序列的功率差率得出第g组非样板机相对于样板机出力的比例系数为以此动态识别各组的出力比例系数;
风电场的理论功率计算阶段,步骤如下:
2.根据权利要求1所述的一种基于机组运行性能评估的风电场理论功率计算方法,其特征在于:所述异常数据识别具体包括:
空缺值的识别:设X={x1,x2,...,xn}为风电机组在统计时间内的实际功率,若存在xi满足:
则判断xi存在数据丢失;
越限值的识别:若存在xi满足下式关系,则判断xi超出正常数值范围:
xi>xh|xi<x1 (2)
式中,xh,x1分别为正常数值范围上限和下限;
不刷新值的识别:若数据连续一段时间不刷新,则除去第一个数据点外均辨识为不刷新值;
畸变数据的识别:相邻风电机组的功率通常具有相似趋势,故基于相邻风电机组出力空间相关性,利用N个与风电机组A具有空间相关性的相邻风电机组的功率数据XM={xM1,xM2,...,xMn},M=1,2,…,N,对观测风电机组的功率数据XA={xA1,xA2,...,xAn}进行校验,识别异常数据,具体如下:
基于式(3)对功率进行归一化,得到观测风电机组的功率归一化数据XA *={xA1 *,xA2 *,...,xAn *}和相邻风电机组的功率归一化数据XM *={xM1 *,xM2 *,...,xMn *},二者在每个采样点的差值为式;
eAi-Mi=|xAi *-xMi *| (4)
风电机组A和M的功率差值和平均值分别为:
EA-M={eA1-M1,eA2-M2,…,eAn-Mn} (5)
3.根据权利要求1所述的一种基于机组运行性能评估的风电场理论功率计算方法,其特征在于:所述异常数据重构具体包括:
单点异常:对于单点异常数据可以采用填充法对异常数据进行工程化重构,即选择一个经验常量或者按照某种规则对异常数据进行替换填充;常量填充:基于人工经验,对单点异常数据采用同一个指定的常量值进行填充;均值填充:可采用数据序列的均值或异常数据前后数据的均值对单点异常数据进行填充;
多点异常:对于多点异常,由于异常数据两侧均存在完好的实际运行数据,因此可以根据两侧实际运行数据,采用基于ARMA模型预测的重构方法对异常数据进行重构;
Pt=β1Pt-1+β2Pt-2+…+βpPt-p+Zt (7)
Zt=εt+λ1εt-1+λ2εt-2+…+λ1εt-1 (8)
式中,βi为自回归参数;λi为滑动平均参数;εi为预测误差;
其次对于第p步预测结果,如果预测误差小于误差阈值ε0时,则认为该预测结果可信;通过对大量完好运行数据进行预测,确定数据序列预测的平均可信度γp:
式中,Count(·)为次数统计;Ep为数据序列预测误差;N为总次数;
连续异常:随着异常数据时间尺度增加,对于基于ARMA模型预测的重构方法,其正确率将会显著降低,由于风电机组运行数据间存在一定的相关性,因此对于连续异常数据,可以采用基于序列时延相关性的重构方法对异常数据进行重构;
时延相关性,假设有两个时间序列X={x1,x2,…,xn}和Y={y1,y2,…,yn},Y相对于X延迟l时刻的相关系数R(l)计算公式为:
基于新能源风电机组出力时延相关性的重构,假设缺失数据时间段为[tm,tn],待重构风电机组为W0,其余风电机组为W1,W2,…,取各风电机组在时段[tm-t0,tm]内的出力数据,首先计算W0与W1在该时段内出力的最大时延相关点和相应的时延相关系数,若R(l1')=max(R(l1'),R(l2'),…),则W0与W1相对应的最大时延相关点为l1',相应的时延相关系数为R(l1');之后再按此分别计算W0与W2,W3…在该时段内出力的最大时延相关点和相应的时延相关系数;若有W0与W1的时延相关性系数最大,则利用W1的出力曲线对W0的历史运行数据进行重构,即W0和W1的出力关系可用线性回归拟合为:
4.根据权利要求1所述的一种基于机组运行性能评估的风电场理论功率计算方法,其特征在于:所述非样板机分组,首先,在非限电情况下,建立风电机组的基本发电性能指标;
月最大功率w1(kW)为观察日前30天内风电机组实发功率Pr={p1,p2,…,pn}中的最大值,即:
w1=max{pi},i=1,2,…,n (13)
式中,n为样本采样点数;
月平均发电功率w2
式中,Pr,j,t为观察日前30天内第j日第t个采样点的实发功率(kW);
平均发电功率年变化率w3,机组平均发电功率年变化率反映机组发电水平随时间变化的趋势:
式中,Pr,av,y为机组观察日前第一个365天内的平均实际发电功率(kW);P′r,av,y为机组观察日前第二个365天内的平均实际发电功率(kW);P′r,j,t为观察日前第一个365天内第j日第t个采样点的实发功率(kW);P″r,j,t为观察日前第二个365天内第j日第t个采样点的实发功率(kW);M为每日采样点个数;
月满发时数w4(h)为观察日前30天内风电机组实发电量与额定功率之比,可用于不同装机容量发电***的比较,其计算式如下:
式中,Pn为观察日前30天内风电机组额定功率(kW),Er为观察日前30天内的实发电量;
利用风电机组基本发电性能指标,建立统计时间内风电机组基本出力性能特征向量W={w1,w2,w3,w4},设k型号风机第i台非样板机在统计时间内的基本出力性能特征向量为k型号样板机在统计时间内的基本出力性能平均特征向量为其中,Mk为k型号样板机全风电场的总数量;为k型号风机第m台样板机在统计时间内的基本出力性能特征向量;
各组出力比例系数的动态识别,设第g组非样板机在统计时间内的实际出力平均值序列为:
设k型号风机样板机在统计时间内的实际出力平均值序列为:
式中,n为在统计时间内的采样点个数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011290260.9A CN112347655B (zh) | 2020-11-17 | 2020-11-17 | 一种基于机组运行性能评估的风电场理论功率计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011290260.9A CN112347655B (zh) | 2020-11-17 | 2020-11-17 | 一种基于机组运行性能评估的风电场理论功率计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112347655A true CN112347655A (zh) | 2021-02-09 |
CN112347655B CN112347655B (zh) | 2023-03-17 |
Family
ID=74364096
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011290260.9A Active CN112347655B (zh) | 2020-11-17 | 2020-11-17 | 一种基于机组运行性能评估的风电场理论功率计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112347655B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115983663A (zh) * | 2022-11-04 | 2023-04-18 | 广西电网有限责任公司 | 基于新能源的发电分析方法及相关设备 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104112236A (zh) * | 2014-05-29 | 2014-10-22 | 国家电网公司 | 风电场发电功率的计算方法 |
CN104124685A (zh) * | 2014-07-28 | 2014-10-29 | 国家电网公司 | 基于样板风机法的风电场理论功率计算方法 |
CN105930933A (zh) * | 2016-04-26 | 2016-09-07 | 华北电力科学研究院有限责任公司 | 风电场理论功率曲线确定方法及装置 |
CN105989236A (zh) * | 2015-03-03 | 2016-10-05 | 南京理工大学 | 一种基于空间相关性的台风过境风电功率预测方法 |
CN107330183A (zh) * | 2017-06-29 | 2017-11-07 | 华北电力大学 | 一种基于运行数据的风电利用率计算方法 |
CN107944622A (zh) * | 2017-11-21 | 2018-04-20 | 华北电力大学 | 基于连续时段聚类的风电功率预测方法 |
CN109960778A (zh) * | 2017-12-26 | 2019-07-02 | 北京金风慧能技术有限公司 | 计算风电场理论功率的方法和装置 |
-
2020
- 2020-11-17 CN CN202011290260.9A patent/CN112347655B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104112236A (zh) * | 2014-05-29 | 2014-10-22 | 国家电网公司 | 风电场发电功率的计算方法 |
CN104124685A (zh) * | 2014-07-28 | 2014-10-29 | 国家电网公司 | 基于样板风机法的风电场理论功率计算方法 |
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 |
CN105989236A (zh) * | 2015-03-03 | 2016-10-05 | 南京理工大学 | 一种基于空间相关性的台风过境风电功率预测方法 |
CN105930933A (zh) * | 2016-04-26 | 2016-09-07 | 华北电力科学研究院有限责任公司 | 风电场理论功率曲线确定方法及装置 |
CN107330183A (zh) * | 2017-06-29 | 2017-11-07 | 华北电力大学 | 一种基于运行数据的风电利用率计算方法 |
CN107944622A (zh) * | 2017-11-21 | 2018-04-20 | 华北电力大学 | 基于连续时段聚类的风电功率预测方法 |
CN109960778A (zh) * | 2017-12-26 | 2019-07-02 | 北京金风慧能技术有限公司 | 计算风电场理论功率的方法和装置 |
Non-Patent Citations (4)
Title |
---|
刘永前等: "基于分段支持向量机的风电机组理论功率计算研究", 《太阳能学报》 * |
张东英等: "风电场有功功率异常运行数据重构方法", 《电力***自动化》 * |
王运等: "基于样板风机法的风电有功控制策略优化", 《宁夏电力》 * |
郭海思等: "风电机组理论发电功率计算及对比优化", 《电力与能源》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115983663A (zh) * | 2022-11-04 | 2023-04-18 | 广西电网有限责任公司 | 基于新能源的发电分析方法及相关设备 |
CN115983663B (zh) * | 2022-11-04 | 2023-09-12 | 广西电网有限责任公司 | 基于新能源的发电分析方法及相关设备 |
Also Published As
Publication number | Publication date |
---|---|
CN112347655B (zh) | 2023-03-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111680820B (zh) | 分布式光伏电站故障诊断方法和装置 | |
CN109308571B (zh) | 配电线路线变关系检测方法 | |
CN108400595B (zh) | 一种考虑新能源出力相关性的电压暂降随机预估方法 | |
CN109842373A (zh) | 基于时空分布特性的光伏阵列故障诊断方法及装置 | |
CN106570790B (zh) | 一种计及风速数据分段特性的风电场出力数据修复方法 | |
CN106682763B (zh) | 一种用于大量样本数据的电力负荷优化预测方法 | |
CN110570122A (zh) | 一种计及风速季节特性及集电***元件故障的海上风电场可靠性评估方法 | |
CN111654392A (zh) | 基于互信息的低压配电网拓扑识别方法及*** | |
CN112952828A (zh) | 一种用于低压配电网的用户台区分析方法 | |
CN111900713A (zh) | 网源协调下考虑负荷和风电随机性多场景输电网规划方法 | |
CN111313403A (zh) | 一种基于马尔可夫随机场的低压配电***网络拓扑识别方法 | |
CN113937764A (zh) | 一种低压配电网高频量测数据处理与拓扑识别方法 | |
CN112347655B (zh) | 一种基于机组运行性能评估的风电场理论功率计算方法 | |
CN113866552B (zh) | 一种基于机器学习的中压配网用户用电异常诊断方法 | |
CN109447512B (zh) | 基于均匀设计的大电网可靠性评估方法 | |
CN111221811A (zh) | 一种基于集抄***的低压配电网络线路参数估计方法 | |
CN108460177B (zh) | 一种大规模多状态串并联***的可靠性近似计算方法 | |
CN112001531B (zh) | 基于有效载荷能力的风电短期运行容量可信度评估方法 | |
CN107204616B (zh) | 基于自适应稀疏伪谱法的电力***随机状态估计方法 | |
CN113505909A (zh) | 一种面向短期风电功率趋势预测的误差补偿方法 | |
CN112182499A (zh) | 一种基于时序电量数据的低压配电网拓扑结构辨识方法 | |
CN116992328A (zh) | 基于无监督学习的低压台区5g高频数据拓扑识别方法 | |
CN113627655B (zh) | 一种配电网灾前故障场景模拟预测方法及装置 | |
CN114169118B (zh) | 考虑分布式电源出力相关性的配电网拓扑结构辨识方法 | |
CN117882261A (zh) | 预测电网频率 |
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 |