一种风电机组功率曲线异常散点识别系统及方法与流程

1.本发明属于风力发电技术领域，特别涉及一种风电机组功率曲线异常散点识别系统及方法。


背景技术：

2.在利用机组运行数据，依据《gbt 18451.2
‑
2012风力发电机组功率特性测试》标准绘制功率曲线过程中，在机组的功率散点图中会出现很多异常的散点，如图2所示。其原因可分为三种：第一，由于电网调度指令，风电场实行限功率措施；第二，在所计算的时段内，机组有启/停机的动作，在计算平均功率的过程中，导致不能反映真实的情况，拉低了功率值；第三，风速风向变化大或者湍流等的影响，风速仪测得的数据与机组实际面对的风况不符。
3.目前，在剔除这类数据的过程中，有根据桨距角、发电机转速和转矩等之间的对应关系判断机组是否出力异常来判断异常散点的，这种方法需要借助额外的参量，增加了数据的维度，并且需要掌握机组内在的控制逻辑，通用性不强。还有通过于均值距离判别(avdc)、局部异常因子(lof)以及考虑噪声的基于密度的聚类(dbscan)等方法识别异常散点，达到剔除异常点，保证功率曲线准确性的目的，但该方法所涉及的算法较复杂，计算量大，效率低。


技术实现要素：

4.本发明的目的在于提供一种风电机组功率曲线异常散点识别系统及方法，解决了现有算法复杂、计算量大及效率低的问题。
5.本发明是通过以下技术方案来实现：
6.一种风电机组功率曲线异常散点识别方法，包括以下步骤：
7.s1、根据风电机组的设计功率曲线，确定每个风速中心点的功率波动范围；
8.s2、根据采集的某时间段内的风速和功率数据，建立风速
‑
功率散点数据集，对风速
‑
功率散点数据集进行分风速段处理，形成多个分风速段的风速
‑
功率散点集合；
9.s3、对风速
‑
功率散点集合中额定风速前的功率数据进行异常散点识别及剔除；对风速
‑
功率散点集合中额定风速后的功率数据异常散点识别及剔除；
10.对剔除完成后的功率数据取中位数，作为该风速区间的参考功率，进而依据s1中对应风速段的功率波动范围及经验值得出该风速段的上限功率值和下限功率值，进而得到该风速中心点的功率散点变化范围；
11.s4、依据得出的风速中心点的功率散点变化范围，对风速
‑
功率散点数据集进行自适应加窗处理，识别并剔除异常散点数据；
12.s5、得到剔除掉异常散点的风速
‑
功率数据，进而得到风电机组实测功率曲线。
13.进一步，s1具体包括以下步骤：
14.s1.1、对设计功率曲线上的功率提升段内的风速段进行分风速段处理，得到若干
个风速中心点；
15.s1.2、确定出每个风速中心点的上限风速和下限风速，据风电机组的设计功率曲线得到每个风速中心点对应的风速中心点功率、上限风速功率和下限风速功率；
16.s1.3、根据风速中心点功率、上限风速功率和下限风速功率计算出该风速中心点在上下限风速内的功率变化率。
17.进一步，s1.3具体为：某风速中心点的风速为x，风速中心点功率为px，上限风速对应的上限风速功率为px

，下限风速对应的下限风速功率为px
‑
，则上限风速相较风速中心点的功率变化率为
18.下限风速相较风速中心点的功率变化率为
19.进一步，s2具体为：以一定的时间间隔t，按时间序列计算平均风速和平均功率，进而得到风速
‑
功率散点数据集[wsn，pn]，n为总的点的个数；wsn代表风速数据，pn代表功率数据；
[0020]
以间隔wst划分风速段，wst的值与风速间隔对应，形成分风速段的风速
‑
功率散点集合[twsx，tpx]，x为对应的风速中心点的值，twsx代表风速中心点为xm/s的风速的集合，tpx代表风速中心点为xm/s的功率的集合。
[0021]
进一步，s3中异常散点识别及剔除具体为：将风速
‑
功率散点集合中的功率值由小到大排序，删除一部分高数据和一部分低数据。
[0022]
进一步，删除一部分高数据和一部分低数据，具体为：将最低的后10％的数据和最高的前2％的数据删除。
[0023]
进一步，s3中，对剔除完成后的功率数据取中位数，取中位数作为该风速区间的参考功率spower
x
；考虑到散点的波动，在确定的vx
‑
和vx

的基础上增加10％作为功率的上下限值，即：
[0024]
spower
x
‑
＝spower
x
*[1
‑
(vx
‑
)
‑
10％]
[0025]
spower
x
＝spower
x
*[1 (vx

) 10％]
[0026]
式中：
[0027]
spower
x
‑
—风速中心为x m/s风速段的功率下限值；
[0028]
spower
x
—风速中心为x m/s风速段的功率上限值。
[0029]
进一步，s4具体为：依据得出的x m/s风速中心点的功率散点变化范围(spower
x
‑
，spower
x
]，以该限值对风速
‑
功率散点集合[tws
x
，tp
x
]中的数据进行识别和剔除，超出(spower
x
‑
，spower
x
]中的认为是异常点，将其剔除。
[0030]
一种风电机组功率曲线异常散点识别系统，包括数据采集模块、分风速段处理模块、异常散点识别模块、功率值计算模块、自适应加窗处理模块和风电机组实测功率曲线生成模块；
[0031]
数据采集模块，用于采集某时间段内的风速和功率数据，建立风速
‑
功率散点数据集；
[0032]
分风速段处理模块，用于对风速
‑
功率散点数据集进行分风速段处理，形成多个分风速段的风速
‑
功率散点集合；
[0033]
异常散点识别模块，用于对风速
‑
功率散点集合中额定风速前的功率数据进行异常散点识别及剔除，对风速
‑
功率散点集合中额定风速后的功率数据进行异常散点识别及剔除；
[0034]
功率值计算模块，用于对剔除完成后的功率数据取中位数，作为该风速区间的参考功率，进而依据对应风速段的功率波动范围及经验值得出该风速段的上限功率值和下限功率值，进而得到该风速中心点的功率散点变化范围；
[0035]
自适应加窗处理模块，用于根据风速中心点的功率散点变化范围，对风速
‑
功率散点数据集进行自适应加窗处理，识别并剔除异常散点数据；
[0036]
风电机组实测功率曲线生成模块，用于根据得到的剔除掉异常散点的风速
‑
功率数据生成风电机组实测功率曲线。
[0037]
与现有技术相比，本发明具有以下有益的技术效果：
[0038]
本发明根据风电机组的设计功率曲线，通过计算风速中心点对应的功率上下限范围，确定该风速区间的功率波动范围，相当于自适应加窗的过程，进而将功率散点数据中的异常数据剔除，减少了异常数据的干扰，保证了功率曲线的准确性。该方法简单实用，大大减小了固定比例加窗筛选数据的误差，准确性高且计算量小，具有很大实际意义。
附图说明
[0039]
图1是设计功率曲线图；
[0040]
图2是采样间隔为1s的风速和功率的原始数据，图2(a)为风速随时间变化的原始数据，图2(b)为功率随时间变化的原始数据；
[0041]
图3是以时间间隔t处理后的风速
‑
功率数据，图3(a)为风速随时间变化的平均数据，图3(b)为功率随时间变化的平均数据；
[0042]
图4为异常功率散点及加窗过程示意图；
[0043]
图5是固定加窗后的示意图；
[0044]
图6是本发明的一种风电机组功率曲线异常散点识别方法处理流程图；
[0045]
图7为异常功率散点及加窗过程示意图。
具体实施方式
[0046]
下面结合具体的实施例对本发明做进一步的详细说明，所述是对本发明的解释而不是限定。
[0047]
如图6所示，本发明公开了一种风电机组功率曲线异常散点识别方法，包括以下步骤：
[0048]
s1.根据图1所示的机组的设计功率曲线，确定不同风速中心点功率的波动范围，该步骤的目的是确定在不同风速中心点功率的变化范围，为后续筛选异常的离散功率点提供依据。
[0049]
s2.参考《gbt 18451.2
‑
2012风力发电机组功率特性测试》标准，以一定的时间间隔t，按时间序列计算平均风速和平均功率，进而得到风速
‑
功率散点数据集，如图4中的所有功率散点。
[0050]
s3.得到风速
‑
功率散点数据集[wsn，pn](n为总的点的个数)后，以间隔wst划分风
速段，wst的值与表1中的风速间隔对应(例如0.5m/s)，形成分风速段的风速
‑
功率散点集合[twsx，tpx](x为对应的风速中心点的值，例如3、3.5、4等)。例如，[tws5，tp5]中包含是数据情况是，风速集合tws5在(4.75，5.25]范围内及其对应的功率值集合tp5。
[0051]
s4.对风速
‑
功率散点集合[twsx，tpx]中额定风速前的每一组进行异常散点识别及剔除处理。
[0052]
1)得出风速中心点的功率散点变化范围(spower
i
‑
，spower
i
]，以该限值对风速
‑
功率散点集合[tws
i
，tp
i
]中的数据进行识别和剔除，超出(spower
i
‑
，spower
i
]中的认为是异常点，将其剔除。
[0053]
对风速
‑
功率散点集合[twsx，tpx]中额定风速前的每一组进行上述操作，得出剔除掉异常散点的风速
‑
功率散点集合。
[0054]
上述步骤1)—5)中提到的10％和2％等数据，可以根据实际处理效果进行优化。
[0055]
由于功率曲线中功率提升段的功率变化是非线性的，每一段的变化率不一样，因此，在原设计功率曲线的基础上得到功率变化率进而确定功率上下限值的过程是一种自适应加窗过程，可以更加准确的识别异常散点，同时也不会像图5中那样固定限值导致散点图出现锯齿状，结果更加异常且失准。
[0056]
s5.风速
‑
功率散点集合[twsx，tpx]中额定风速后的异常散点识别及剔除。由于额定风速后，机组的功率值是确定的，因此可以根据机组正常的功率波动设置一个固定的变化率，一般取5％(经验值)即可。
[0057]
s6.根据s4和s5的方法，得到最终正常的风速
‑
功率散点的数据[nwsx，npx]绘制功率曲线。
[0058]
以下给一个具体的实施例来说明：
[0059]
s1.以某2mw级别机型的机组为例说明，设计的数据如下表所示：
[0060]
根据该表格建立如图1所示的机组的设计功率曲线，此处的风速间隔是0.5m/s，可通过软件计算或者插值的方法得到风速间隔更小的数据，例如0.25m/s。
[0061]
如图1，由于在功率提升段(3m/s
‑
10m/s范围)功率曲线的斜率是变化的，因此，按照风速段(表1中以0.5m/s为间隔举例)划分风速中心点，计算在风速中心点上下限风速内功率的波动范围。
[0062]
表1设计功率曲线波动范围计算
[0063]
[0064][0065]
根据设计的功率曲线，计算功率曲线波动范围，此处以5m/s风速中心的数据举例分析。
[0066][0067]
此时，得到的v5
‑
＝|p5
‑
(p5
‑
)|/p5＝0.1445、v5

＝|p5
‑
(p5

)|/p5＝0.1678即为5m/s的风速中心处，功率的波动范围。
[0068]
当到达额定功率以后，即使风速再增加，功率也会稳定在一个范围内，不会再继续上升，因此，这一部分的异常数据剔除就相对简单，在额定功率的基础上设定一个固定的波动范围就可以，这里选择5％，这个值可以根据实际的情况调整。
[0069]
s2.以某机组运行约1个月的数据举例
[0070]
图3是在图2的基础上，每t时间分别对风速和功率做平均(t可以是2～10分钟)。例如，原始数据是1s采集一个数，得到图2的风速
‑
时间图(见图2(a))和功率
‑
时间图(见图2(b))。但是，鉴于风速的波动性太大，且数据量太多，因此，按时间顺序，每2分钟对风速和功率做一个平均，这样，平抑了风速的波动，同时减少了数据量，得到了图3所示的风速
‑
时间图(见图3(a))和功率
‑
时间图(见图3(b))。这是《gbt 18451.2
‑
2012风力发电机组功率特性测试》标准中的通用方法。
[0071]
把图2转换成横坐标是风速，纵坐标是功率的形式就可以得到功率散点图，如图4所示，对应的数据集就是[wsn，pn]。
[0072]
其中，wsn—t时间间隔平均处理后的风速，ws(wind speed)；
[0073]
pn—t时间间隔平均处理后的功率，p(power)。
[0074]
s3.以间隔0.5m/s划分风速段，与表1中的风速间隔对应，形成分风速段的风速
‑
功率散点集合[twsx，tpx](x为对应的风速中心点的值，例如3、3.5、4等)。例如，[tws5，tp5]中包含是数据情况，以图4黑色方框内的5m/s为例，风速中心点为5m/s的风速
‑
功率散点集合[tws5，tp5]，其中，tws5—风速中心点为5m/s的风速的集合；tp5—风速中心点为5m/s的功率的集合。
[0075]
s4.风速中心点为5m/s的风速
‑
功率散点集合[tws5，tp5]为例进行处理说明。
[0076]
(1)对风速
‑
功率散点集合[tws5，tp5]，按功率的值进行由小到大的排序；
[0077]
(2)排完序的功率值，将最低的后10％的数据和最高的前2％的数据删除；
[0078]
(3)删除完成后的数据，取其中位数作为该风速区间的参考功率spower5；
[0079]
(4)考虑到散点的波动，在步骤s1的表1中确定的v5
‑
和v5

的基础上增加10％
[0080]
(该数值可根据具体情况进行优化)作为功率的上下限值，即：spower5
‑
＝spower5*[1
‑
(v5
‑
)
‑
10％]
[0081]
spower5

＝spower5*[1 (v5

) 10％]
[0082]
式中：
[0083]
spower5
‑
—风速中心为5m/s风速段的功率下限值；
[0084]
spower5

—风速中心为5m/s风速段的功率上限值。
[0085]
(5)得出5m/s风速中心点的功率散点变化范围(spower5
‑
，spower5

]，以该限值对风速
‑
功率散点集合[tws5，tp5]中的数据进行识别和剔除，超出(spower5
‑
，spower5

]中的认为是异常点，将其剔除。
[0086]
s5.风速
‑
功率散点集合[twsx，tpx]中额定风速后的异常散点识别及剔除。由于额定风速后，机组的功率值是确定的，因此可以根据机组正常的功率波动设置一个固定的变化率，一般取5％(经验值)即可，如图7中横向方框内所示，在这个范围，机组到达额定功率，即使风速上升，功率也不会线性增加，而是维持在一个比较稳定的状态。
[0087]
s6.根据s4和s5的方法，得到最终正常的风速
‑
功率散点的数据[nwsx，npm]，进而绘制功率曲线，其中，nwsx—剔除完异常数据的风速的集合；npx—剔除完异常数据的功率的集合。
[0088]
本发明还公开了一种风电机组功率曲线异常散点识别系统，包括数据采集模块、分风速段处理模块、异常散点识别模块、功率值计算模块、自适应加窗处理模块和风电机组实测功率曲线生成模块；
[0089]
数据采集模块，用于采集某时间段内的风速和功率数据，建立风速
‑
功率散点数据集；
[0090]
分风速段处理模块，用于对风速
‑
功率散点数据集进行分风速段处理，形成多个分风速段的风速
‑
功率散点集合；
[0091]
异常散点识别模块，用于对风速
‑
功率散点集合中额定风速前的功率数据进行异常散点识别及剔除，对风速
‑
功率散点集合中额定风速后的功率数据进行异常散点识别及剔除；
[0092]
功率值计算模块，用于对剔除完成后的功率数据取中位数，作为该风速区间的参考功率，进而依据对应风速段的功率波动范围及经验值得出该风速段的上限功率值和下限
功率值，进而得到该风速中心点的功率散点变化范围；
[0093]
自适应加窗处理模块，用于根据风速中心点的功率散点变化范围，对风速
‑
功率散点数据集进行自适应加窗处理，识别并剔除异常散点数据；
[0094]
风电机组实测功率曲线生成模块，用于根据得到的剔除掉异常散点的风速
‑
功率数据生成风电机组实测功率曲线。