CN105930933B - 风电场理论功率曲线确定方法及装置 - Google Patents
风电场理论功率曲线确定方法及装置 Download PDFInfo
- Publication number
- CN105930933B CN105930933B CN201610266284.8A CN201610266284A CN105930933B CN 105930933 B CN105930933 B CN 105930933B CN 201610266284 A CN201610266284 A CN 201610266284A CN 105930933 B CN105930933 B CN 105930933B
- Authority
- CN
- China
- Prior art keywords
- data
- power
- wind speed
- acquisition system
- single machine
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 55
- 230000002159 abnormal effect Effects 0.000 claims abstract description 96
- 230000005856 abnormality Effects 0.000 claims abstract description 12
- 238000012360 testing method Methods 0.000 claims description 23
- 241001269238 Data Species 0.000 claims description 14
- 238000009826 distribution Methods 0.000 claims description 13
- 238000010586 diagram Methods 0.000 claims description 12
- 238000010276 construction Methods 0.000 claims description 8
- 238000003064 k means clustering Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000007689 inspection Methods 0.000 claims description 4
- 239000004744 fabric Substances 0.000 claims description 3
- 230000008030 elimination Effects 0.000 description 14
- 238000003379 elimination reaction Methods 0.000 description 14
- 230000005611 electricity Effects 0.000 description 7
- 230000006978 adaptation Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 125000006850 spacer group Chemical group 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Wind Motors (AREA)
Abstract
本发明公开了一种风电场理论功率曲线确定方法及装置,该方法包括:获取风电场各风机在预设时间内的历史数据,历史数据包括风速和功率;基于历史数据构建单机风速功率散点图;剔除单机风速功率散点图中的异常数据;针对剔除异常数据后的单机风速功率散点数据,分别计算以预设风速间隔划分的各个区间对应的数据集合的聚类中心,并按聚类中心构建单机等值功率曲线;基于单机等值功率曲线,对所剔除的异常数据点中的风速进行线性插值计算对应的理论功率数据,并按该理论功率数据重构单机理论功率数据;根据各风机理论功率数据得到风电场理论功率曲线。本发明仅使用风速和功率两项数据能最大限度排除干扰数据,还原真实风电场理论功率曲线。
Description
技术领域
本发明涉及风电场技术领域,尤其涉及一种风电场理论功率曲线确定方法及装置。
背景技术
近年来,随着风电高速发展,我国弃风限电问题也越趋严重。风电场的理论功率数据是准确研究风电功率预测、客观评价弃风限电情况和合理规划风电容量配置的重要基础,因此有必要对风电历史运行数据中的限电以及故障等异常数据进行还原重构,以最大限度地反映风电场真实理论出力水平。
对于风电场理论功率的研究主要基于历史数据回归重构的角度进行开展,包括标杆风机法、机舱风速法等。然而,风电场实际运行过程中影响数据质量的因素有很多,由于条件约束、外部干扰和设备故障等原因,实际运行数据中充斥着大量异常数据,严重破坏了各因素之间的对应关系,造成回归重构结果存在偏差,精度不足。
实际运行中,风电场除了因限电带来大量数据外,众多因素会造成各数据之间的对应关系发生畸变,如较小的风速(切入风速以下)对应着较大的风电功率输出(甚至达到额定功率),因此不能仅依据风机是否处于限电状态来判断风机运行数据是否正常,而且状态信息也同时存在一些异常数据。
针对上述问题,目前尚未提出有效的解决方案。
发明内容
本发明提供了一种风电场理论功率曲线确定方法及装置,以至少解决现有技术中大量异常数据的存在导致通过历史数据回归重构得到的风电场理论功率存在偏差的问题。
根据本发明的一个方面,提供了一种风电场理论功率曲线确定方法,包括:获取风电场中各风机在预设时间内的历史数据,其中所述历史数据包括:风速和功率;基于所述历史数据构建单机风速功率散点图;剔除所述单机风速功率散点图中的异常数据;针对剔除异常数据后的单机风速功率散点数据,按照以预设风速间隔划分的区间,分别计算各个区间对应的数据集合的聚类中心,并按照所述聚类中心构建单机等值功率曲线;基于所述单机等值功率曲线,对所剔除的异常数据点中的风速进行线性插值计算得到对应的理论功率数据,并根据所述理论功率数据重构单机在所述预设时间内的理论功率数据;根据各风机在所述预设时间内的理论功率数据得到风电场在所述预设时间内的理论功率曲线。
在一个实施例中,剔除所述单机风速功率散点图中的异常数据,包括:利用罗曼诺夫斯基准则剔除所述单机风速功率散点图中的异常数据;以及利用聚类分析法剔除所述单机风速功率散点图中剩余的异常数据。
在一个实施例中,利用罗曼诺夫斯基准则剔除所述单机风速功率散点图中沿第一方向分布的异常数据,包括:在所述单机风速功率散点图中,按照预设功率间隔将风机功率范围划分为p个功率区间,并按功率区间将风速功率散点数据划分为p个数据集合X1,X2,…Xp;将数据集合Xi内的受怀疑数据xh去掉,计算该数据集合Xi内不包括受怀疑数据xh的风速数据的算数平均值其中,1≤i≤p;计算该数据集合Xi内不包括受怀疑数据xh的风速数据的标准偏差s1;根据第一预设显著水平α1和该数据集合Xi内包含的风速数据量m,在t分布表中查出检验系数K(α1,m);若则确定受怀疑数据xh为异常数据并进行剔除,在该数据集合Xi内除受怀疑数据xh之外的风速数据中选定新的受怀疑数据,重复执行计算算术平均值、标准偏差和检验系数的步骤,直到剔除该数据集合Xi内所有的异常数据;若则确定该数据集合Xi内除受怀疑数据xh之外的风速数据均不是异常数据。
在一个实施例中,采用以下公式计算该数据集合Xi内不包括受怀疑数据xh的风速数据的算数平均值 其中,xi表示数据集合Xi内的第i个风速数据;采用以下公式计算该数据集合Xi内不包括受怀疑数据xh的风速数据的标准偏差s1:
在一个实施例中,利用罗曼诺夫斯基准则剔除所述单机风速功率散点图中沿第二方向分布的异常数据,包括:在所述单机风速功率散点图中,按照预设风速间隔将风机风速范围划分为q个风速区间,并按风速区间将风速功率散点数据划分为q个数据集合Y1,Y2,…Yq;将数据集合Yi内的受怀疑数据yh去掉,计算该数据集合Yi内不包括受怀疑数据yh的功率数据的算数平均值其中,1≤i≤q;计算该数据集合Yi内不包括受怀疑数据yh的功率数据的标准偏差s2;根据第二预设显著水平α2和该数据集合Yi内包含的功率数据量n,在t分布表中查出检验系数K(α2,n);若则确定受怀疑数据yh为异常数据并进行剔除,在该数据集合Yi内除受怀疑数据yh之外的功率数据中选定新的受怀疑数据,重复执行计算算术平均值、标准偏差和检验系数的步骤,直到剔除该数据集合Yi内所有的异常数据;若则确定该数据集合Yi内除受怀疑数据yh之外的功率数据均不是异常数据。
在一个实施例中,采用以下公式计算该数据集合Yi内不包括受怀疑数据yh的功率数据的算数平均值 其中,yi表示数据集合Yi内的第i个功率数据;采用以下公式计算该数据集合Yi内不包括受怀疑数据yh的功率数据的标准偏差s2:
在一个实施例中,利用聚类分析法剔除所述单机风速功率散点图中剩余的异常数据,包括:在所述单机风速功率散点图中,按照预设风速间隔将风机风速范围划分为q个风速区间,并按风速区间将风速功率散点数据划分为q个数据集合Y1,Y2,…Yq;基于k均值聚类算法将数据集合Yi内的风速功率散点数据划分为k类数据C1,C2,…Ck,并计算各个聚类的聚类中心w1,w2,…wk;确定聚类中心w1,w2,…wk中的最大聚类中心wmax以及对应的聚类Cmax;确定与最大聚类中心wmax的距离大于预设值的所有聚类,并将确定的聚类作为异常数据进行剔除。
在一个实施例中,采用以下公式对所剔除的异常数据点中的风速进行线性插值计算得到对应的理论功率数据P:其中,(xi,yi),(xj,yj)表示所述单机等值功率曲线中相邻的两点,v表示所剔除的异常数据点中的风速,且v∈[xi,xj]。
根据本发明的另一个方面,提供了一种风电场理论功率曲线确定装置,包括:获取单元,用于获取风电场中各风机在预设时间内的历史数据,其中所述历史数据包括:风速和功率;第一构建单元,用于基于所述历史数据构建单机风速功率散点图;剔除单元,用于剔除所述单机风速功率散点图中的异常数据;第二构建单元,用于针对剔除异常数据后的单机风速功率散点数据,按照以预设风速间隔划分的区间,分别计算各个区间对应的数据集合的聚类中心,并按照所述聚类中心构建单机等值功率曲线;重构单元,用于基于所述单机等值功率曲线,对所剔除的异常数据点中的风速进行线性插值计算得到对应的理论功率数据,并根据所述理论功率数据重构单机在所述预设时间内的理论功率数据;计算单元,用于根据各风机在所述预设时间内的理论功率数据得到风电场在所述预设时间内的理论功率曲线。
在一个实施例中,所述剔除单元包括:第一剔除子单元,用于利用罗曼诺夫斯基准则剔除所述单机风速功率散点图中的异常数据;以及第二剔除子单元,用于利用聚类分析法剔除所述单机风速功率散点图中剩余的异常数据。
在一个实施例中,所述剔除单元包括:第一划分模块,用于在所述单机风速功率散点图中,按照预设功率间隔将风机功率范围划分为p个功率区间,并按功率区间将风速功率散点数据划分为p个数据集合X1,X2,…Xp;第一计算模块,用于将数据集合Xi内的受怀疑数据xh去掉,计算该数据集合Xi内不包括受怀疑数据xh的风速数据的算数平均值其中,1≤i≤p;第二计算模块,用于计算该数据集合Xi内不包括受怀疑数据xh的风速数据的标准偏差s1;查找模块,用于根据第一预设显著水平α1和该数据集合Xi内包含的风速数据量m,在t分布表中查出检验系数K(α1,m);第一确定模块,用于在的情况下,则确定受怀疑数据xh为异常数据并进行剔除,在该数据集合Xi内除受怀疑数据xh之外的风速数据中选定新的受怀疑数据,重复执行计算算术平均值、标准偏差和检验系数的步骤,直到剔除该数据集合Xi内所有的异常数据;第二确定模块,用于在的情况下,则确定该数据集合Xi内除受怀疑数据xh之外的风速数据均不是异常数据。
在一个实施例中,所述第一计算模块具体用于采用以下公式计算该数据集合Xi内不包括受怀疑数据xh的风速数据的算数平均值 其中,xi表示数据集合Xi内的第i个风速数据;所述第二计算模块具体用于采用以下公式计算该数据集合Xi内不包括受怀疑数据xh的风速数据的标准偏差s1:
在一个实施例中,所述第二剔除子单元包括:第二划分模块,用于在所述单机风速功率散点图中,按照预设风速间隔将风机风速范围划分为q个风速区间,并按风速区间将风速功率散点数据划分为q个数据集合Y1,Y2,…Yq;第三计算模块,用于基于k均值聚类算法将数据集合Yi内的风速功率散点数据划分为k类数据C1,C2,…Ck,并计算各个聚类的聚类中心w1,w2,…wk;第三确定模块,用于确定聚类中心w1,w2,…wk中的最大聚类中心wmax以及对应的聚类Cmax;第四确定模块,用于确定与最大聚类中心wmax的距离大于预设值的所有聚类,并将确定的聚类作为异常数据进行剔除。
通过本发明的风电场理论功率曲线确定方法及装置,利用异常数据剔除和重构,最大限度地排除干扰数据,尽可能还原真实的风电场理论功率曲线;同时,该方法仅使用了风速和功率两项数据,实施简单,适应于各类复杂环境。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的限定。在附图中:
图1是本发明实施例的风电场理论功率曲线确定方法的流程图;
图2是本发明实施例的风电场理论功率曲线确定装置的结构框图;
图3是本发明实施例的基于异常数据剔除和重构的风电场理论功率曲线确定方法的示意图;
图4是本发明实施例的风电场实际功率示意图;
图5是本发明实施例的001#风机实际功率示意图;
图6是本发明实施例的001#风机实际风速示意图;
图7是本发明实施例的001#风机的风速功率散点图;
图8是本发明实施例的001#风机的等值功率曲线示意图;
图9是本发明实施例的风电场实际功率和理论功率曲线示意图。
具体实施方式
下面结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明的保护范围。
本发明实施例提供了一种风电场理论功率曲线确定方法,图1是本发明实施例的风电场理论功率曲线确定方法的流程图,如图1所示,该方法包括:
步骤S101,获取风电场中各风机在预设时间内的历史数据,其中历史数据包括:风速和功率。该预设时间可以是几个月、一年或者几年。历史数据中还可以包括切入风速和切出风速,其中,切入风速为风机并网发电的启动风速,切出风速为风机并网发电的最大风速。
步骤S102,基于上述历史数据构建单机风速功率散点图。
本发明对于单机风速功率散点图的横纵坐标不做限制,即以风速为横坐标,功率为纵坐标,或者,以功率为横坐标,风速为纵坐标。较优的,单机风速功率散点图以风速为横坐标,以功率为纵坐标。
步骤S103,剔除单机风速功率散点图中的异常数据。
步骤S104,针对剔除异常数据后的单机风速功率散点数据,按照以预设风速间隔划分的区间,分别计算各个区间对应的数据集合的聚类中心,并按照聚类中心构建单机等值功率曲线。对于散点图,数据集合中的数据可以用坐标表示,聚类中心实际是数据集合中的某一个数据,也可用坐标表示,将各个聚类中心连接成曲线,则得到等值功率曲线。
步骤S105,基于单机等值功率曲线,对所剔除的异常数据点中的风速进行线性插值计算得到对应的理论功率数据,并根据该理论功率数据重构单机在预设时间内的理论功率数据。具体的,对异常数据点中的风速进行插值计算得到对应的准确的功率数据,从而可以得到较为准确的单机理论功率数据。
步骤S106,根据各风机在预设时间内的理论功率数据得到风电场在预设时间内的理论功率曲线。其中,可以先对各风机的理论功率数据进行求和,得到风电场的理论功率数据,然后以时间为横坐标,功率为纵坐标建立风电场理论功率曲线。
通过上述方法,利用异常数据剔除和重构,最大限度地排除干扰数据,尽可能还原真实的风电场理论功率曲线;同时,该方法仅使用了风速和功率两项数据,实施简单,适应于各类复杂环境。
在一个实施例中,步骤S103中剔除单机风速功率散点图中的异常数据可以包括:利用罗曼诺夫斯基准则剔除单机风速功率散点图中的异常数据;以及利用聚类分析法剔除单机风速功率散点图中剩余的异常数据。
利用罗曼诺夫斯基准则进行异常数据剔除,可以包括:剔除横向分布差异的数据和/或剔除纵向分布差异的数据,二者方法类似,横向是指在单机风速功率散点图中以纵坐标为基准划分区间或数据集合,纵向是指在单机风速功率散点图中以横坐标为基准划分区间或数据集合。但是利用罗曼诺夫斯基准则并不能剔除全部异常数据,例如不能剔除限电数据,因此,本发明还利用聚类分析法进行剩余异常数据的剔除,从而尽可能剔除全部异常数据。针对单机风速功率散点图中数据点的分布特点,利用聚类分析法主要是纵向上的异常数据剔除。
当然,也可以先利用聚类分析法剔除单机风速功率散点图中的异常数据,再利用罗曼诺夫斯基准则剔除剩余的异常数据。
在一个实施例中,利用罗曼诺夫斯基准则剔除单机风速功率散点图中沿第一方向分布的异常数据,可以包括以下步骤:
(1)在单机风速功率散点图中,按照预设功率间隔将风机功率范围划分为p个功率区间,并按功率区间将风速功率散点数据划分为p个数据集合X1,X2,…Xp。优选的,预设功率间隔可以是某一固定值,例如5kW。预设功率间隔也可以是不等的间隔,例如在合理的情况下按2kW、3kW、4kW逐个递增,或者以2kW、3kW和5kW这三个间隔进行循环。
(2)将数据集合Xi内的受怀疑数据xh去掉,计算该数据集合Xi内不包括受怀疑数据xh的风速数据的算数平均值其中,1≤i≤p。具体的,可以采用以下公式计算该算数平均值 其中,m表示数据集合Xi内包含的风速数据量,xi表示数据集合Xi内的第i个风速数据。
(3)计算该数据集合Xi内不包括受怀疑数据xh的风速数据的标准偏差s1。具体的,可以采用以下公式计算该标准偏差s1:其中,m表示数据集合Xi内包含的风速数据量,xi表示数据集合Xi内的第i个风速数据。
(4)根据第一预设显著水平α1和该数据集合Xi内包含的风速数据量m,在t分布表中查出检验系数K(α1,m)。其中,显著水平是估计总体参数落在某一区间内,可能犯错误的概率,显著水平可以预先设定,一般取值为0.05;t分布表中α1和m数列交叉的数即为相应的检验系数。
(5)若则确定受怀疑数据xh为异常数据并进行剔除,在该数据集合Xi内除受怀疑数据xh之外的风速数据中选定新的受怀疑数据,重复执行计算算术平均值、标准偏差和检验系数的步骤,直到剔除该数据集合Xi内所有的异常数据;若则确定该数据集合Xi内除受怀疑数据xh之外的风速数据均不是异常数据。
本实施例中,按照预设功率间隔划分区间和数据集合,根据之前对横向和纵向的含义说明,可知,如果单机风速功率散点图以功率作为纵坐标,则本实施例属于横向异常数据剔除;如果单机风速功率散点图以功率作为横坐标,则本实施例属于纵向异常数据剔除。本实施例中使用罗曼诺夫斯基准则剔除异常数据,能够较为精确地剔除大部分异常数据,保证最终得到的风电场理论功率曲线的准确性较高。
在另一个实施例中,利用罗曼诺夫斯基准则剔除单机风速功率散点图中沿第二方向分布的异常数据,可以包括以下步骤:
(1)在单机风速功率散点图中,按照预设风速间隔将风机风速范围划分为q个风速区间,并按风速区间将风速功率散点数据划分为q个数据集合Y1,Y2,…Yq。优选的,预设风速间隔可以是某一固定值,例如1m/s。预设风速间隔也可以是不等的间隔,例如在合理的情况下按1m/s、2m/s、3m/s逐个递增,或者以1m/s、2m/s和3m/s这三个间隔进行循环。
(2)将数据集合Yi内的受怀疑数据yh去掉,计算该数据集合Yi内不包括受怀疑数据yh的功率数据的算数平均值其中,1≤i≤q。具体的,可以采用以下公式计算该算数平均值 其中,n表示数据集合Yi内包含的功率数据量,yi表示数据集合Yi内的第i个功率数据。
(3)计算该数据集合Yi内不包括受怀疑数据yh的功率数据的标准偏差s2。具体的,可以采用以下公式计算该标准偏差s2:其中,n表示数据集合Yi内包含的功率数据量,yi表示数据集合Yi内的第i个功率数据。
(4)根据第二预设显著水平α2和该数据集合Yi内包含的功率数据量n,在t分布表中查出检验系数K(α2,n)。其中,显著水平是估计总体参数落在某一区间内,可能犯错误的概率,显著水平可以预先设定,一般取值为0.05;t分布表中α2和n数列交叉的数即为相应的检验系数。
(5)若则确定受怀疑数据yh为异常数据并进行剔除,在该数据集合Yi内除受怀疑数据yh之外的功率数据中选定新的受怀疑数据,重复执行计算算术平均值、标准偏差和检验系数的步骤,直到剔除该数据集合Yi内所有的异常数据;若则确定该数据集合Yi内除受怀疑数据yh之外的功率数据均不是异常数据。
本实施例中,按照预设风速间隔划分区间和数据集合,根据之前对横向和纵向的含义说明,可知,如果单机风速功率散点图以风速作为纵坐标,则本实施例属于横向异常数据剔除;如果单机风速功率散点图以风速作为横坐标,则本实施例属于纵向异常数据剔除。本实施例中使用罗曼诺夫斯基准则剔除异常数据,能够较为精确地剔除大部分异常数据,保证最终得到的风电场理论功率曲线的准确性较高。
需要说明的是,为了尽可能的剔除所有异常数据,提高最终得到的理论功率曲线的正确性,可以同时使用上述两个实施例中描述的利用罗曼诺夫斯基准则剔除异常数据的方法,即同时使用横向异常数据剔除和纵向异常数据剔除,执行顺序不区分先后。
利用聚类分析法剔除单机风速功率散点图中剩余的异常数据,可以包括以下步骤:
(1)在单机风速功率散点图中,按照预设风速间隔将风机风速范围划分为q个风速区间,并按风速区间将风速功率散点数据划分为q个数据集合Y1,Y2,…Yq。优选的,预设风速间隔可以是某一固定值,例如1m/s。预设风速间隔也可以是不等的间隔,例如在合理的情况下按1m/s、2m/s、3m/s逐个递增,或者以1m/s、2m/s和3m/s这三个间隔进行循环。
(2)基于k均值聚类算法将数据集合Yi内的风速功率散点数据划分为k类数据C1,C2,…Ck,并计算各个聚类的聚类中心w1,w2,…wk。
具体的,k均值聚类算法的步骤如下:
①为每个聚类确定一个初始聚类中心,这样就有k个初始聚类中心;
②将各个聚类中的数据按照最小距离原则分配到与其距离最小的初始聚类中心所在的聚类;
③将每个聚类中的数据的均值作为新的聚类中心;
④重复步骤②、③,直到聚类中心不再变化;
⑤流程结束,得到k个聚类及其对应的k个聚类中心。
(3)确定聚类中心w1,w2,…wk中的最大聚类中心wmax以及对应的聚类Cmax。
(4)确定与最大聚类中心wmax的距离大于预设值的所有聚类,并将确定的聚类作为异常数据进行剔除。其中,该预设值的取值可以为0.2Pr,Pr为风机的额定功率。
本实施例中使用聚类分析法剔除异常数据,能够较为精确地剔除大部分异常数据,保证最终得到的风电场理论功率曲线的准确性较高。
步骤S105中,可以采用以下公式对所剔除的异常数据点中的风速进行线性插值计算得到对应的理论功率数据P:
其中,(xi,yi),(xj,yj)表示单机等值功率曲线中相邻的两点,v表示所剔除的异常数据点中的风速,且v∈[xi,xj]。
基于同一发明构思,本发明实施例还提供了一种风电场理论功率曲线确定装置,可以用于实现上述实施例所描述的方法。由于该装置解决问题的原理与上述方法相似,因此该装置的实施可以参见上述方法的实施,重复之处不再赘述。以下所使用的,术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的***较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。
图2是本发明实施例的风电场理论功率曲线确定装置的结构框图,如图2所示,该装置包括:获取单元21、第一构建单元22、剔除单元23、第二构建单元24、重构单元25和计算单元26,下面对该结构进行具体说明。
获取单元21,用于获取风电场中各风机在预设时间内的历史数据,其中历史数据包括:风速和功率;
第一构建单元22,用于基于上述历史数据构建单机风速功率散点图;
剔除单元23,用于剔除单机风速功率散点图中的异常数据;
第二构建单元24,用于针对剔除异常数据后的单机风速功率散点数据,按照以预设风速间隔划分的区间,分别计算各个区间对应的数据集合的聚类中心,并按照聚类中心构建单机等值功率曲线;
重构单元25,用于基于单机等值功率曲线,对所剔除的异常数据点中的风速进行线性插值计算得到对应的理论功率数据,并根据该理论功率数据重构单机在预设时间内的理论功率数据;
计算单元26,用于根据各风机在预设时间内的理论功率数据得到风电场在预设时间内的理论功率曲线。
通过上述装置,利用异常数据剔除和重构,最大限度地排除干扰数据,尽可能还原真实的风电场理论功率曲线;同时,该装置仅使用了风速和功率两项数据,实施简单,适应于各类复杂环境。
在一个实施例中,剔除单元23可以包括:第一剔除子单元,用于利用罗曼诺夫斯基准则剔除单机风速功率散点图中的异常数据;以及第二剔除子单元,用于利用聚类分析法剔除单机风速功率散点图中剩余的异常数据。
在一个实施例中,第一剔除子单元可以包括:
第一划分模块,用于在单机风速功率散点图中,按照预设功率间隔将风机功率范围划分为p个功率区间,并按功率区间将风速功率散点数据划分为p个数据集合X1,X2,…Xp;
第一计算模块,用于将数据集合Xi内的受怀疑数据xh去掉,计算该数据集合Xi内不包括受怀疑数据xh的风速数据的算数平均值其中,1≤i≤p;
第二计算模块,用于计算该数据集合Xi内不包括受怀疑数据xh的风速数据的标准偏差s1;
第一查找模块,用于根据第一预设显著水平α1和该数据集合Xi内包含的风速数据量m,在t分布表中查出检验系数K(α1,m);
第一确定模块,用于在的情况下,则确定受怀疑数据xh为异常数据并进行剔除,在该数据集合Xi内除受怀疑数据xh之外的风速数据中选定新的受怀疑数据,重复执行计算算术平均值、标准偏差和检验系数的步骤,直到剔除该数据集合Xi内所有的异常数据;
第二确定模块,用于在的情况下,则确定该数据集合Xi内除受怀疑数据xh之外的风速数据均不是异常数据。
具体的,第一计算模块可以采用以下公式计算该数据集合Xi内不包括受怀疑数据xh的风速数据的算数平均值
第二计算模块可以采用以下公式计算该数据集合Xi内不包括受怀疑数据xh的风速数据的标准偏差s1:
在一个实施例中,第一剔除子单元还可以包括:
第二划分模块,用于在单机风速功率散点图中,按照预设风速间隔将风机风速范围划分为q个风速区间,并按风速区间将风速功率散点数据划分为q个数据集合Y1,Y2,…Yq;
第三计算模块,用于将数据集合Yi内的受怀疑数据yh去掉,计算该数据集合Yi内不包括受怀疑数据yh的功率数据的算数平均值其中,1≤i≤q;
第四计算模块,用于计算该数据集合Yi内不包括受怀疑数据yh的功率数据的标准偏差s2;
第二查找模块,用于根据第二预设显著水平α2和该数据集合Yi内包含的功率数据量n,在t分布表中查出检验系数K(α2,n);
第三确定模块,用于在的情况下,则确定受怀疑数据yh为异常数据并进行剔除,在该数据集合Yi内除受怀疑数据yh之外的功率数据中选定新的受怀疑数据,重复执行计算算术平均值、标准偏差和检验系数的步骤,直到剔除该数据集合Yi内所有的异常数据;
第四确定模块,用于在的情况下,则确定该数据集合Yi内除受怀疑数据yh之外的功率数据均不是异常数据。
具体的,第三计算模块可以采用以下公式计算该数据集合Yi内不包括受怀疑数据yh的功率数据的算数平均值
第四计算模块可以采用以下公式计算该数据集合Yi内不包括受怀疑数据yh的功率数据的标准偏差s2:
第二剔除子单元可以包括:
第三划分模块,用于在单机风速功率散点图中,按照预设风速间隔将风机风速范围划分为q个风速区间,并按风速区间将风速功率散点数据划分为q个数据集合Y1,Y2,…Yq;
第五计算模块,用于基于k均值聚类算法将数据集合Yi内的风速功率散点数据划分为k类数据C1,C2,…Ck,并计算各个聚类的聚类中心w1,w2,…wk;
第五确定模块,用于确定聚类中心w1,w2,…wk中的最大聚类中心wmax以及对应的聚类Cmax;
第六确定模块,用于确定与最大聚类中心wmax的距离大于预设值的所有聚类,并将确定的聚类作为异常数据进行剔除。
重构单元25可以采用以下公式对所剔除的异常数据点中的风速进行线性插值计算得到对应的理论功率数据P:
其中,(xi,yi),(xj,yj)表示单机等值功率曲线中相邻的两点,v表示所剔除的异常数据点中的风速,且v∈[xi,xj]。
当然,上述单元或模块划分只是一种示意划分,本发明并不局限于此。只要能实现本发明的目的的模块划分,均应属于本发明的保护范围。
为了对上述风电场理论功率曲线确定方法及装置进行更为清楚的解释,下面结合具体的实施例来进行说明,然而值得注意的是该实施例仅是为了更好地说明本发明,并不构成对本发明不当的限定。
图3是本发明实施例的基于异常数据剔除和重构的风电场理论功率曲线确定方法的示意图。在本实施例中,散点图以风速为横坐标,功率为纵坐标。
(1)获得风电场中各风机的风速全年历史数据v1,v2,…vn和功率全年历史数据P1,P2,…Pn(其中n为风机数量),其中包括切入风速vr和切出风速vc。
(2)基于各风机历史数据构建以风速vi为横坐标和以功率Pi为纵坐标的单机风速功率散点图。
(3)利用横向罗曼诺夫斯基准则剔除单机风速功率散点图中横向分布差异的数据。具体如下:
①在单机风速功率散点图中,以PΔ为间隔单位将风机功率范围[0,Pe]划分为p个区间,并按功率区间将风速功率数据划分为p个数据集合X1,X2,…Xp,之后针对各个数据集合,利用罗曼诺夫斯基准则进行异常数据剔除。
②将数据集合Xi内的受怀疑数据xh(xh可选取与均值偏差最大的数据)去掉,计算数据集合Xi内不包括xh的风速数据的算数平均值1≤i≤p。算数平均值的具体计算公式参见上述描述,此处不再赘述。
③计算数据集合Xi内不包括xh的风速数据的标准偏差s。标准偏差s的具体计算公式参见上述描述,此处不再赘述。
④根据预设显著水平α(一般取0.05)和数据集合Xi内包含的风速数据量m,在t分布表中查出检验系数K(α,m)。
⑤若则可确定xh为异常数据并进行剔除,重复执行步骤②、③、④,直到剔除该数据集合Xi内所有的异常数据;若则可确定剩余风速数据(该数据集合Xi内除受怀疑数据xh之外的风速数据)均不是异常数据。
(4)利用纵向罗曼诺夫斯基准则剔除单机风速功率散点图中纵向分布差异的数据,其具体实现方法与步骤(3)类似,此处不再赘述。
(5)利用纵向聚类分析法剔除异常数据。
①在单机风速功率散点图中,以vΔ为间隔单位将风机风速范围[0,vc]划分为q个区间,并按功率区间将风速功率数据划分为q个数据集合Y1,Y2,…Yq,之后针对各个数据集合,利用聚类分析法进行异常数据剔除。
②基于k均值聚类算法将数据集合Yi内的风速功率散点数据划分为k类数据C1,C2,…Ck,并计算各个聚类的聚类中心w1,w2,…wk,聚类中心即与数据集合内各点的距离和最小的那个数据。
③找出w1,w2,…wk中的最大聚类中心wmax以及对应的聚类Cmax。
④找出所有与最大聚类中心wmax的距离大于预设值T(一般取0.2Pr,Pr为风机的额定功率)的聚类,并将这些聚类视为异常数据进行剔除。
(6)针对剔除异常数据后的单机风速功率散点数据,按照划分的q个区间,分别计算各个区间对应的数据集合的聚类中心(x1,y1),(x2,y2),…,(xq,yq),并按照聚类中心构建单机等值功率曲线。
(7)基于构建的单机等值功率曲线,对所剔除的异常数据点中的风速v进行线性插值计算得到对应的理论功率数据P,并根据计算的理论功率数据P重构单机全年理论功率数据P′i。功率数据P的具体计算公式参见上述描述,此处不再赘述。
(8)将各风机全年理论功率数据求和得到风电场全年理论功率曲线。
以某地区电网内一座风电场为例进行计算分析,对上述风电场理论功率曲线确定方法进行解释说明。
该风电场并网装机容量99.5MW,共39台风机,其中2MW风机11台,2.5MW风机13台,3MW风机15台。
(1)获取该风电场各风机2015年1-8月的实际功率和实际风速数据,风电场实际功率如图4所示,1号(001#)风机实际功率如图5所示,001#风机实际风速如图6所示。
(2)以风速vi为横坐标,以功率Pi为纵坐标构建各单机的风速功率散点图,001#风机的风速功率散点图如图7所示。
(3)选取PΔ为20kW,vΔ为0.5m/s,分别利用罗曼诺夫斯基准则和k均值聚类剔除异常数据,并构建单机等值功率曲线。001#风机的等值功率曲线如图8所示,与图6相比,去除了大量的异常数据,其中灰色的曲线为001#风机的等值功率曲线。
(4)重构风机2015年1-8月的理论功率曲线,进而得到风电场2015年1-8月的理论功率曲线,如图9所示,为风电场1月1日至1月11日的实际功率和理论功率曲线示意图,其中,实线表示理论功率,虚线表示实际功率。
(5)弃风电量评估
根据该风电场2015年1-8月的实际功率和理论功率曲线,可知该风电场1-8月实际发电量为1.26亿千瓦时,理论发电量为1.88亿千瓦时,弃风电量为0.62亿千瓦时。
综上所述,本发明提出了一种基于异常数据剔除和重构的风电场理论功率曲线确定方法及装置,该方案将异常数据剔除和重构应用到风电场理论功率估算中,通过剔除异常数据来还原风速理论功率曲线,并以此回归重构异常数据点对应的理论功率数值来构建风电场理论功率曲线,最大限度地排除干扰数据,尽可能地还原真实的风电场理论功率曲线。
流程图中或在此以其他方式描述的任何过程或方法描述可以被理解为,表示包括一个或更多个用于实现特定逻辑功能或过程的步骤的可执行指令的代码的模块、片段或部分,并且本发明的优选实施方式的范围包括另外的实现,其中可以不按所示出或讨论的顺序,包括根据所涉及的功能按基本同时的方式或按相反的顺序,来执行功能,这应被本发明的实施例所属技术领域的技术人员所理解。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (9)
1.一种风电场理论功率曲线确定方法,其特征在于,包括:
获取风电场中各风机在预设时间内的历史数据,其中所述历史数据包括:风速和功率;
基于所述历史数据构建单机风速功率散点图;
剔除所述单机风速功率散点图中的异常数据;
剔除所述单机风速功率散点图中的异常数据,包括:
利用罗曼诺夫斯基准则剔除所述单机风速功率散点图中的异常数据;以及
利用聚类分析法剔除所述单机风速功率散点图中剩余的异常数据;
利用罗曼诺夫斯基准则剔除所述单机风速功率散点图中沿第一方向分布的异常数据,包括:
在所述单机风速功率散点图中,按照预设功率间隔将风机功率范围划分为p个功率区间,并按功率区间将风速功率散点数据划分为p个数据集合X1,X2,…Xp;
将数据集合Xi内的受怀疑数据xh去掉,计算该数据集合Xi内不包括受怀疑数据xh的风速数据的算数平均值其中,1≤i≤p;
计算该数据集合Xi内不包括受怀疑数据xh的风速数据的标准偏差s1;
根据第一预设显著水平α1和该数据集合Xi内包含的风速数据量m,在t分布表中查出检验系数K(α1,m);
若则确定受怀疑数据xh为异常数据并进行剔除,在该数据集合Xi内除受怀疑数据xh之外的风速数据中选定新的受怀疑数据,重复执行计算算术平均值、标准偏差和检验系数的步骤,直到剔除该数据集合Xi内所有的异常数据;
若则确定该数据集合Xi内除受怀疑数据xh之外的风速数据均不是异常数据;
针对剔除异常数据后的单机风速功率散点数据,按照以预设风速间隔划分的区间,分别计算各个区间对应的数据集合的聚类中心,并按照所述聚类中心构建单机等值功率曲线;
基于所述单机等值功率曲线,对所剔除的异常数据点中的风速进行线性插值计算得到对应的理论功率数据,并根据所述理论功率数据重构单机在所述预设时间内的理论功率数据;
根据各风机在所述预设时间内的理论功率数据得到风电场在所述预设时间内的理论功率曲线。
2.根据权利要求1所述的方法,其特征在于,
采用以下公式计算该数据集合Xi内不包括受怀疑数据xh的风速数据的算数平均值 其中,xi表示数据集合Xi内的第i个风速数据;
采用以下公式计算该数据集合Xi内不包括受怀疑数据xh的风速数据的标准偏差s1:
3.根据权利要求1所述的方法,其特征在于,利用罗曼诺夫斯基准则剔除所述单机风速功率散点图中沿第二方向分布的异常数据,包括:
在所述单机风速功率散点图中,按照预设风速间隔将风机风速范围划分为q个风速区间,并按风速区间将风速功率散点数据划分为q个数据集合Y1,Y2,…Yq;
将数据集合Yi内的受怀疑数据yh去掉,计算该数据集合Yi内不包括受怀疑数据yh的功率数据的算数平均值其中,1≤i≤q;
计算该数据集合Yi内不包括受怀疑数据yh的功率数据的标准偏差s2;
根据第二预设显著水平α2和该数据集合Yi内包含的功率数据量n,在t分布表中查出检验系数K(α2,n);
若则确定受怀疑数据yh为异常数据并进行剔除,在该数据集合Yi内除受怀疑数据yh之外的功率数据中选定新的受怀疑数据,重复执行计算算术平均值、标准偏差和检验系数的步骤,直到剔除该数据集合Yi内所有的异常数据;
若则确定该数据集合Yi内除受怀疑数据yh之外的功率数据均不是异常数据。
4.根据权利要求3所述的方法,其特征在于,
采用以下公式计算该数据集合Yi内不包括受怀疑数据yh的功率数据的算数平均值 其中,yi表示数据集合Yi内的第i个功率数据;
采用以下公式计算该数据集合Yi内不包括受怀疑数据yh的功率数据的标准偏差s2:
5.根据权利要求1所述的方法,其特征在于,利用聚类分析法剔除所述单机风速功率散点图中剩余的异常数据,包括:
在所述单机风速功率散点图中,按照预设风速间隔将风机风速范围划分为q个风速区间,并按风速区间将风速功率散点数据划分为q个数据集合Y1,Y2,…Yq;
基于k均值聚类算法将数据集合Yi内的风速功率散点数据划分为k类数据C1,C2,…Ck,并计算各个聚类的聚类中心w1,w2,…wk;
确定聚类中心w1,w2,…wk中的最大聚类中心wmax以及对应的聚类Cmax;
确定与最大聚类中心wmax的距离大于预设值的所有聚类,并将确定的聚类作为异常数据进行剔除。
6.根据权利要求1所述的方法,其特征在于,采用以下公式对所剔除的异常数据点中的风速v进行线性插值计算得到对应的理论功率数据P:
其中,(xi,yi),(xj,yj)表示所述单机等值功率曲线中相邻的两点,v表示所剔除的异常数据点中的风速,且v∈[xi,xj]。
7.一种风电场理论功率曲线确定装置,其特征在于,包括:
获取单元,用于获取风电场中各风机在预设时间内的历史数据,其中所述历史数据包括:风速和功率;
第一构建单元,用于基于所述历史数据构建单机风速功率散点图;
剔除单元,用于剔除所述单机风速功率散点图中的异常数据;
所述剔除单元包括:
第一剔除子单元,用于利用罗曼诺夫斯基准则剔除所述单机风速功率散点图中的异常数据;以及
第二剔除子单元,用于利用聚类分析法剔除所述单机风速功率散点图中剩余的异常数据;
所述第一剔除子单元包括:
第一划分模块,用于在所述单机风速功率散点图中,按照预设功率间隔将风机功率范围划分为p个功率区间,并按功率区间将风速功率散点数据划分为p个数据集合X1,X2,…Xp;
第一计算模块,用于将数据集合Xi内的受怀疑数据xh去掉,计算该数据集合Xi内不包括受怀疑数据xh的风速数据的算数平均值其中,1≤i≤p;
第二计算模块,用于计算该数据集合Xi内不包括受怀疑数据xh的风速数据的标准偏差s1;
查找模块,用于根据第一预设显著水平α1和该数据集合Xi内包含的风速数据量m,在t分布表中查出检验系数K(α1,m);
第一确定模块,用于在的情况下,则确定受怀疑数据xh为异常数据并进行剔除,在该数据集合Xi内除受怀疑数据xh之外的风速数据中选定新的受怀疑数据,重复执行计算算术平均值、标准偏差和检验系数的步骤,直到剔除该数据集合Xi内所有的异常数据;
第二确定模块,用于在的情况下,则确定该数据集合Xi内除受怀疑数据xh之外的风速数据均不是异常数据;
第二构建单元,用于针对剔除异常数据后的单机风速功率散点数据,按照以预设风速间隔划分的区间,分别计算各个区间对应的数据集合的聚类中心,并按照所述聚类中心构建单机等值功率曲线;
重构单元,用于基于所述单机等值功率曲线,对所剔除的异常数据点中的风速进行线性插值计算得到对应的理论功率数据,并根据所述理论功率数据重构单机在所述预设时间内的理论功率数据;
计算单元,用于根据各风机在所述预设时间内的理论功率数据得到风电场在所述预设时间内的理论功率曲线。
8.根据权利要求7所述的装置,其特征在于,
所述第一计算模块具体用于采用以下公式计算该数据集合Xi内不包括受怀疑数据xh的风速数据的算数平均值 其中,xi表示数据集合Xi内的第i个风速数据;
所述第二计算模块具体用于采用以下公式计算该数据集合Xi内不包括受怀疑数据xh的风速数据的标准偏差s1:
9.根据权利要求7所述的装置,其特征在于,所述第二剔除子单元包括:
第二划分模块,用于在所述单机风速功率散点图中,按照预设风速间隔将风机风速范围划分为q个风速区间,并按风速区间将风速功率散点数据划分为q个数据集合Y1,Y2,…Yq;
第三计算模块,用于基于k均值聚类算法将数据集合Yi内的风速功率散点数据划分为k类数据C1,C2,…Ck,并计算各个聚类的聚类中心w1,w2,…wk;
第三确定模块,用于确定聚类中心w1,w2,…wk中的最大聚类中心wmax以及对应的聚类Cmax;
第四确定模块,用于确定与最大聚类中心wmax的距离大于预设值的所有聚类,并将确定的聚类作为异常数据进行剔除。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610266284.8A CN105930933B (zh) | 2016-04-26 | 2016-04-26 | 风电场理论功率曲线确定方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610266284.8A CN105930933B (zh) | 2016-04-26 | 2016-04-26 | 风电场理论功率曲线确定方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105930933A CN105930933A (zh) | 2016-09-07 |
CN105930933B true CN105930933B (zh) | 2019-05-31 |
Family
ID=56837132
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610266284.8A Active CN105930933B (zh) | 2016-04-26 | 2016-04-26 | 风电场理论功率曲线确定方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105930933B (zh) |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107247997A (zh) * | 2017-05-18 | 2017-10-13 | 北京唐浩电力工程技术研究有限公司 | 一种风电场风机电量分析方法 |
CN107330183B (zh) * | 2017-06-29 | 2020-04-17 | 华北电力大学 | 一种基于运行数据的风电利用率计算方法 |
CN107527057B (zh) * | 2017-09-07 | 2020-03-31 | 国能日新科技股份有限公司 | 一种风速功率异常数据剔除方法及装置 |
CN107843745B (zh) * | 2017-09-20 | 2020-02-07 | 新疆金风科技股份有限公司 | 风力发电机组的风速仪的故障诊断方法和装置、存储介质 |
CN109751195B (zh) * | 2017-11-01 | 2020-05-15 | 中国船舶重工集团海装风电股份有限公司 | 一种风力发电机功率曲线的获取方法及装置 |
CN109960778B (zh) * | 2017-12-26 | 2023-06-27 | 北京金风慧能技术有限公司 | 计算风电场理论功率的方法和装置 |
CN108287320A (zh) * | 2018-02-01 | 2018-07-17 | 安徽创世科技股份有限公司 | 一种电池电量检测优化方法 |
CN108763584B (zh) * | 2018-06-11 | 2021-11-02 | 北京天泽智云科技有限公司 | 一种风功率曲线散点过滤的方法及其*** |
CN110334138B (zh) * | 2019-03-15 | 2023-02-24 | 贵州航天林泉电机有限公司 | 一种基于matlab的数据一致性的分析方法 |
CN110334310B (zh) * | 2019-07-01 | 2023-04-07 | 北京源清慧虹信息科技有限公司 | 一种风机功率曲线的确定方法及确定装置 |
CN110826904B (zh) * | 2019-11-01 | 2022-08-02 | 三一重能股份有限公司 | 风机的数据处理方法、装置、处理设备及可读存储介质 |
CN111563543B (zh) * | 2020-04-26 | 2023-07-04 | 国网冀北电力有限公司电力科学研究院 | 一种风电机组的风速-发电功率数据的清洗方法及装置 |
CN112347655B (zh) * | 2020-11-17 | 2023-03-17 | 国网青海省电力公司 | 一种基于机组运行性能评估的风电场理论功率计算方法 |
CN112800103B (zh) * | 2020-12-09 | 2023-03-03 | 华能陕西发电有限公司 | 一种scada数据挖掘的风电机组功率曲线的获取***及方法 |
CN112950403A (zh) * | 2021-01-29 | 2021-06-11 | 上海电气风电集团股份有限公司 | 风机多工况数据异常处理的方法、装置及计算机可读存储介质 |
CN113536198B (zh) * | 2021-07-13 | 2022-12-13 | 中国华能集团清洁能源技术研究院有限公司 | 一种风电机组功率曲线异常散点识别***及方法 |
CN115170820B (zh) * | 2022-05-13 | 2023-08-01 | 中铁西北科学研究院有限公司 | 一种应用于数据曲线过渡阶段的特征提取及界限识别方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103291544A (zh) * | 2013-06-21 | 2013-09-11 | 华北电力大学 | 数字化风电机组功率曲线绘制方法 |
CN105160060A (zh) * | 2015-07-17 | 2015-12-16 | 中国电力科学研究院 | 一种基于实际功率曲线拟合的风电场理论功率确定方法 |
CN105389634A (zh) * | 2015-12-01 | 2016-03-09 | 广东智造能源科技研究有限公司 | 一种组合式短期风电功率预测***及方法 |
-
2016
- 2016-04-26 CN CN201610266284.8A patent/CN105930933B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103291544A (zh) * | 2013-06-21 | 2013-09-11 | 华北电力大学 | 数字化风电机组功率曲线绘制方法 |
CN105160060A (zh) * | 2015-07-17 | 2015-12-16 | 中国电力科学研究院 | 一种基于实际功率曲线拟合的风电场理论功率确定方法 |
CN105389634A (zh) * | 2015-12-01 | 2016-03-09 | 广东智造能源科技研究有限公司 | 一种组合式短期风电功率预测***及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105930933A (zh) | 2016-09-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105930933B (zh) | 风电场理论功率曲线确定方法及装置 | |
CN101592538B (zh) | 一种基于实测数据风电场稳态输出功率的计算方法 | |
CN105628977B (zh) | 风速仪修正系数的标定方法、装置和*** | |
CN105160060B (zh) | 一种基于实际功率曲线拟合的风电场理论功率确定方法 | |
CN102682222B (zh) | 一种基于风电波动规律的连续潮流计算方法 | |
US20230012218A1 (en) | Method and apparatus for detecting yaw-to-wind abnormality, and device and storage medium thereof | |
CN104156892A (zh) | 一种有源配电网电压跌落仿真与评估方法 | |
CN109784563B (zh) | 一种基于虚拟测风塔技术的超短期功率预测方法 | |
CN106548256A (zh) | 一种风电场时空动态相关性建模的方法和*** | |
CN109167387A (zh) | 风场风电功率预测方法 | |
DK2981711T3 (en) | PROCEDURE FOR EFFICIENCY MONITORING OF A WINDOW ENERGY PARK | |
CN108629520B (zh) | 一种微气象环境下的高压输电线路运行状态评估方法 | |
CN103902837A (zh) | 一种基于经验Copula函数进行风速预测的方法 | |
CN103475021A (zh) | 一种基于统计模型的风电场弃风电量确定方法 | |
CN110648249A (zh) | 年度电力平衡测算方法、装置及设备 | |
CN101414317A (zh) | 一种处理风电场静态功率等值分散性问题的等效风速法 | |
CN110390343B (zh) | 一种空间气象数据的订正方法及*** | |
CN103277250B (zh) | 基于标杆机组的大型风电基地弃风电量的评估方法 | |
CN107229824A (zh) | 光伏电站发电单元功率曲线建模方法及装置 | |
Di Piazza et al. | Environmental data processing by clustering methods for energy forecast and planning | |
CN106684905A (zh) | 一种考虑风电预测不确定性的风电场动态等值方法 | |
CN103440405A (zh) | 基于实测数据的风电场稳态建模方法及*** | |
CN105825002A (zh) | 一种基于动态灰关联分析法的风电场动态等值建模方法 | |
CN105279384A (zh) | 一种基于风力机机舱风速的来流风速计算方法及装置 | |
CN105356462A (zh) | 风电场谐波预测评估方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |