一种基于非参数核密度估计的风电功率数据修正方法
【专利摘要】本发明提出一种基于非参数核密度估计的风电功率数据修正方法，通过分析风速与功率的散点图，对实测风速小区间内的功率值进行置信度处理，将给定置信区间外的数据修正到置信区间之内。该方法无需考虑数据原本的分布情况，对数据的分布规律具有非常好的包容性。
【专利说明】
-种基于非参数核密度估计的风电功率数据修正方法
技术领域
[0001] 本发明属于电网风电预测领域，特别是设及到一种基于非参数核密度估计的风电 功率数据修正方法。
【背景技术】
[0002] 由于风电场原始采集的数据量庞大、数据种类较多，原始数据含有缺失、错误、重 复和噪声等各类数据问题，运对风电场各类数据挖掘、风资源和电场后评估和功率预测等 带来了严重的干扰。如果简单的抛弃缺失值、错误值，对数据挖掘的效果可能造成严重影 响，也改变了原始数据的真实有效性;比如错估风电场的年平均风速，低估或高估风电场的 出力等。数据修正是正确评估风电场情况的先决条件。
[0003] 目前主要依赖于《风电功率预测功能规范》和《风电场风能资源监测数据管理办 法》运两个文件中的方法对数据进行处理，但是上述提到的两个文件只是满足要求的一个 最低标准，对于更高精度的要求没有写出处理方式，同时不能处理风功率数据不一致的情 况。手动、半自动方式依赖于人的经验，不同的人修正的数据效果不一样，很难达到统一，同 时也需要大量的时间。
【发明内容】
[0004] 有鉴于此，本发明提出一种基于非参数核密度估计的风电功率数据修正方法，该 方法无需考虑数据原本的分布情况，对数据的分布规律具有非常好的包容性。
[0005] 为达到上述目的，本发明的技术方案是运样实现的：一种基于非参数核密度估计 的风电功率数据修正方法，包括：
[0006] 1)读取实测风速功率数据并划分为多个子区间；
[0007] 2)建立子区间的功率密度函数；
[000引3)计算各子区间实测风电功率满足给定置信水平的最短置信区间；
[0009] 4)对置信区间内的数据进行拟合，修正置信区间外的数据。
[0010] 进一步的，步骤1)所述划分子区间的方法为：
[0011] 101)将实测风速按从小到大顺序排列，W排序后的实测风速为横坐标，对应的实 测风电功率为纵坐标，形成历史实测风速-实测风电功率统计数据组；
[0012] 102) W-定的分辨率对实测风速进行等间隔划分，将实测风电功率数据划分为多 个子区间；
[OOK] 103)将每组数据分另赃行功率区间划分。
[0014] 更进一步的，步骤102)中，子区间数公式为：
[0015]
[0016] 区间公式为：
[0017] Di=[Pmin+(i-l) AP，Pmin+i AP]
[001 引 i = i,2,? ? ?,n
[0019]其中，Pmax为风电功率最大值，Pmin为最小值，A P为功率间隔;n为子区间数，Di为区 间。
[0020] 进一步的，步骤2)的具体方法为：
[0021] 统计各子区内的实测风电功率分布，分别建立各子区间实测风电功率的概率密度 函数；
[0022] 对于某一风电功率区间的风电功率数据，其概率密度函数为：
[0023]
[0024] 式中:KO为核函数;Xi为风电功率样本;h为带宽系数;N为样本数。
[0025] 更进一步的，所述核函数为高斯核函数。
[00%]进一步的，步骤3)的具体方法为：
[0027] 应用拉格朗日乘子法求取各子区间实测风电功率满足给定置信水平的最短置信 区间。
[0028] 进一步的，步骤4)的具体方法为：
[0029] 对置信区间内的数据进行拟合，拟合出一条风功率曲线，置信区间之外的数据根 据拟合出的曲线修正到置信区间之内。
[0030] 相对于现有技术，本发明所述的一种基于非参数核密度估计的风电功率数据修正 方法，主要有W下优势：
[0031] 本发明无需考虑数据原本的分布情况，对数据的分布规律具有非常好的包容性。 本发明保留可信程度较高的数据并将可信程度较低的风电功率数据修正到可信度的范围 内，提高了建模数据的有效性。对风电场各类数据挖掘、风资源和电场后评估和功率预测的 准确性的提高有重要意义。
【附图说明】
[0032] 图1是本发明的流程图。
[0033] 图2是本发明实施例中的核密度图。
[0034] 图3是本发明实施例中的最短置信区间图。
[0035] 图4是本发明实施例中的拟合图。
[0036] 图5是本发明实施例中的原始数据和修正结果对比图。
【具体实施方式】
[0037] 需要说明的是，在不冲突的情况下，本发明中的实施例及实施例中的特征可W相 互组合。
[0038] 下面将结合实施例W及附图来详细说明本发明。
[0039] 如图1所示为本发明的流程图。
[0040] 按照图1所示流程，首先，读取实测风速功率数据，将实测风速按从小到大顺序排 列，W排序后的实测风速为横坐标，对应的实测风电功率为纵坐标，形成历史实测风速-实 测风电功率统计数据组，W-定的分辨率对实测风速进行等间隔划分，将实测风电功率数 据划分为多个子区间；
[0041 ]设风电功率最大值为Pmax，最小值为Pmin，功率间隔为A P，则子区间数n为：
[0042]
[00创则所分区间化为：
[0044] Di= [Pmin+(i-l) A P，Pmin+i A 門
[0045] i = l，2, ? ? ?，n
[0046] 按照上述方法，将每组数据分别进行功率区间划分。
[0047] 然后，统计各子区内的实测风电功率分布，分别建立各子区间实测风电功率的概 率密度函数；
[0048] 对于某一风电功率区间的风电功率数据，其概率密度函数为：
[0049]
[0050] 式中：KO为核函数，本发明取高斯核函数;Xi为风电功率样本;h为带宽系数;N为 样本数。图2为核密度图。
[0051] 第=步，应用拉格朗日乘子法求取各子区间实测风电功率满足给定置信水平的最 短置信区间。
[0化2 ] 对于功率值X，给定置信水平a ( 0 < a < 1 )，使Xmin、Xmax满足P ( Xmin < X < Xmax ) = 1 -a， 并^ 且?巧5/S.Xmax~Xmin.^si^J'' :
[0化3] 求解过程如下：
[0化4]目标：min( Xmax-Xmin)
[0化5]约束条件：P(Xmin<X<Xmax) = l-a
[0056] 构造拉格朗日函数
[0057] L( Xmin , Xmax) - ( Xmax-Xmin) +人（P( Xmin<X<Xmax) -1+口）
[005引对该函数进行求解，得到Xmin和Xmax,即最短置信区间，图3中深黑色为原始数据，浅 灰色为标注的最短置信区间。
[0059] 最后，对置信区间内的数据进行拟合，拟合出一条风功率曲线，置信区间之外的数 据根据拟合出的曲线修正到置信区间之内。
[0060] 图4为对置信区间的数据进行拟合，拟合结果为浅灰色线。
[0061] 图5为原始数据和修正结果对比图，其中图5a为原始数据图，图5b为修正结果图。 可W看到，适用本发明的方法无需考虑数据原本的分布情况，对数据的分布规律具有非常 好的包容性。本发明保留可信程度较高的数据并将可信程度较低的风电功率数据修正到可 信度的范围内，提高了建模数据的有效性。对风电场各类数据挖掘、风资源和电场后评估和 功率预测的准确性的提高有重要意义。
[0062] W上所述仅为本发明的较佳实施例而已，并不用W限制本发明，凡在本发明的精 神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
【主权项】
1. 一种基于非参数核密度估计的风电功率数据修正方法，其特征在于，包括： 1) 读取实测风速功率数据并划分为多个子区间； 2) 建立子区间的功率密度函数； 3) 计算各子区间实测风电功率满足给定置信水平的最短置信区间； 4) 对置信区间内的数据进行拟合，修正置信区间外的数据。2. 根据权利要求1所述的一种基于非参数核密度估计的风电功率数据修正方法，其特 征在于，步骤1)所述划分子区间的方法为： 101) 将实测风速按从小到大顺序排列，以排序后的实测风速为横坐标，对应的实测风 电功率为纵坐标，形成历史实测风速-实测风电功率统计数据组； 102) 以一定的分辨率对实测风速进行等间隔划分，将实测风电功率数据划分为多个子 区间； 103) 将每组数据分别进行功率区间划分。3. 根据权利要求2所述的一种基于非参数核密度估计的风电功率数据修正方法，其特 征在于，步骤102)中，子区间数公式为：区间公式为： Di=[Pmin+(i-l) AP,Pmin+iAP] i = l，2, · · ·，n 其中，Pmax为风电功率最大值，Pmin为最小值，A P为功率间隔;n为子区间数，Di为区间。4. 根据权利要求1所述的一种基于非参数核密度估计的风电功率数据修正方法，其特 征在于，步骤2)的具体方法为： 统计各子区内的实测风电功率分布，分别建立各子区间实测风电功率的概率密度函 数； 对于某一风电功率冈间的风电功率数据，其概率密度函数为：式中:KO为核函数;X1为风电功率样本;h为带宽系数;N为样本数。5. 根据权利要求4所述的一种基于非参数核密度估计的风电功率数据修正方法，其特 征在于，所述核函数为高斯核函数。6. 根据权利要求1所述的一种基于非参数核密度估计的风电功率数据修正方法，其特 征在于，步骤3)的具体方法为： 应用拉格朗日乘子法求取各子区间实测风电功率满足给定置信水平的最短置信区间。7. 根据权利要求1所述的一种基于非参数核密度估计的风电功率数据修正方法，其特 征在于，步骤4)的具体方法为： 对置信区间内的数据进行拟合，拟合出一条风功率曲线，置信区间之外的数据根据拟 合出的曲线修正到置信区间之内。
【文档编号】G06F19/00GK106021972SQ201610644055
【公开日】2016年10月12日
【申请日】2016年8月8日
【发明人】向婕, 雍正, 吕建驰
【申请人】北京国能日新系统控制技术有限公司