CN105590027A - 光伏功率异常数据的识别方法 - Google Patents
光伏功率异常数据的识别方法 Download PDFInfo
- Publication number
- CN105590027A CN105590027A CN201510948258.9A CN201510948258A CN105590027A CN 105590027 A CN105590027 A CN 105590027A CN 201510948258 A CN201510948258 A CN 201510948258A CN 105590027 A CN105590027 A CN 105590027A
- Authority
- CN
- China
- Prior art keywords
- photovoltaic power
- data
- abnormal data
- probability
- power
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Photovoltaic Devices (AREA)
Abstract
本发明提供一种光伏功率异常数据的识别方法,包括:利用光伏电站辐照度和光伏功率实测数据,拟合辐照度-光伏功率Copula函数参数;根据Copula函数描述辐照度和光伏功率这两个随机变量的相关关系建立概率功率曲线;通过辐照度-光伏功率散点图,总结异常数据点特征建立异常数据判别准则;基于Copula函数和异常数据判别准则,识别筛选异常数据,建立新数据集;以及若有异常数据被识别,在剔除异常数据后重复以上步骤,对新数据集继续进行识别;若无,则直接利用判别准则和概率功率曲线,识别原始数据集中异常数据。本方法适用于各类光伏电站光伏功率异常数据识别,具有普遍性,且可解决原始数据异常数据比例高的问题。
Description
技术领域
本发明属于电力系统新能源发电领域。尤其涉及一种在光伏功率数据含有高比例异常数据的条件下,光伏功率异常数据的识别方法。
背景技术
准确可信的光伏功率时间序列数据是光伏发电性能分析、功率预测等工作的基础。然而很多光伏电站现场采集的功率数据质量较差,极大阻碍了这些数据的信息挖掘和深化应用。产生光伏功率异常数据的原因众多,如通信故障、设备异常、人为限电等,其中人为限电导致的光伏功率数据异常问题在我国尤为严重。高比例的光伏功率异常数据会对光伏功率和光资源波动性研究、光伏功率与辐照度、温度等因素之间真实规律的提取造成不利影响,直接利用现场数据会降低光伏功率预测模型的精度和有效性,也会对光伏电站的运行管理、电网的运营调度产生不利影响。因此需要进行光伏功率高比例异常数据识别。
进行光伏功率异常数据识别的难点在于准确把握光伏功率数据的强随机性、分散性等特性及我国光伏数据含高比例异常数据的实际情况。现有方法多依赖对组件状态的检测,如根据组件温度和工作电压关系筛选异常数据,需要对工作电压进行检测,但是现阶段我国缺少光伏组件层级的数据,且在光伏发电性能分析和光伏功率预测中更关注场站层面的异常数据。
然而,现有方法不能完全适应我国光伏实际情况;对光伏功率分布的某个概率密度函数假定,如利用3-sigma原则进行异常数据识别,按这类假定,各个辐照度区间内光伏功率的分布规律是相互独立的,但是实际上若辐照度和光伏功率是两个具有关联性的随机变量,各个辐照度区间内的功率分布规律也不是独立的,若进行独立处理不符合实际规律,使得在实际应用中对异常数据的识别能力有限;简单异常识别规则的制定,采用规则识别的方法,一方面过分依赖于经验规律,另一方面未考虑光伏功率与主要影响因素的相关关系,使得在实际应用中识别效果欠佳;较少针对高比例异常数据问题,限电、设备故障等因素导致光伏功率数据含异常数据比例高,在实际应用中,高比例异常数据会导致统计分析结果偏离真实情况,使得异常数据误识别率高。
发明内容
综上所述,确有必要提供一种能够适应我国数据以光伏场站级数据类型为主、异常数据比例高特点的异常数据识别方法。
一种光伏功率异常数据的识别方法,包括以下步骤:利用光伏电站辐照度和光伏功率实测数据,拟合辐照度-光伏功率Copula函数参数;根据Copula函数描述辐照度和光伏功率这两个随机变量的相关关系建立概率功率曲线;通过辐照度-光伏功率散点图,总结异常数据点特征建立异常数据判别准则;基于Copula函数和异常数据判别准则,识别筛选异常数据,建立新数据集;以及若有异常数据被识别,在剔除异常数据后重复以上步骤,对新数据集继续进行识别;若无,则直接利用判别准则和概率功率曲线,识别原始数据集中异常数据。
相对于现有技术,本发明将辐照度和光伏功率作为两个随机变量,利用Copula函数描述两随机变量的相关关系,建立概率功率曲线,根据实测数据的观察建立异常数据判别准则,并针对异常数据比例高的特点,通过反复进行减少原数据集中异常数据对辐照度-功率之间真实规律发现的干扰。该方法易于操作、具有普适性,不需要复杂的优化计算,切合我国光伏数据以场站级数据为主、异常数据比例高等实际特点。
附图说明
图1为本发明提供的光伏功率异常数据的识别方法的流程图。
图2为本发明提供的光伏功率异常数据的识别方法的流程框图。
图3为甘肃某光伏电站辐照度-功率散点及异常数据类型示意图。
具体实施方式
下面根据说明书附图并结合具体实施例对本发明的技术方案进一步详细表述。
请参阅图1及图2,本发明提供的光伏功率异常数据的识别方法包括以下步骤:
步骤S10,利用光伏电站辐照度和光伏功率实测数据,拟合辐照度-光伏功率Copula函数参数;
步骤S20,根据Copula函数描述辐照度和光伏功率这两个随机变量的相关关系建立概率功率曲线;
步骤S30,通过辐照度-光伏功率散点图,总结异常数据点特征建立异常数据判别准则;
步骤S40,基于Copula函数和异常数据判别准则,识别筛选异常数据,建立新数据集;以及
步骤S50,若有异常数据被识别,在剔除异常数据后跳转至步骤S10,对新数据集继续进行识别;若无,则直接利用判别准则和概率功率曲线,识别原始数据集中异常数据。
在步骤S10中,所述辐照度-光伏功率Copula函数参数可通过以下步骤得到:
步骤S11,获取光伏电站辐照度和功率实测数据,进行数据归一化操作,筛选已标记的错误数据;
步骤S12,利用统计分析得到光伏功率P的累积概率分布函数和辐照度R的累积概率分布函数;
步骤S13,结合辐照度-光伏功率散点图的观察,确定选取的Copula函数类型;
步骤S14,利用和,得到唯一的Copula函数C连接R和R,,拟合Copula函数参数。
在步骤S11中,已标记的错误数据指的是在数据采集过程中已被系统识别并标记为错误数据的数据。
在步骤S12中,常用的Copula函数有ClaytonCopula,GaussianCopula,FrankCopula函数等。
可以理解,上述所述的辐照度-光伏功率Copula函数参数的获取方法仅仅为具体的实施例及表达方式,还可以根据实际需要进行选择。
在步骤S20中,所述概率功率曲线可通过以下方法计算:
步骤S21,给定辐照度累计概率分布取值,确定光伏功率累计概率分布的条件概率分布函数。
步骤S22,设功率曲线的置信概率为,即有的数据落于概率区间内,的数据落于区间外。
步骤S23,设置信区间的不对称系数为,计算置信区间上下边界的分位概率,,表示数据点高于上边界的概率为,低于下边界的概率为。
步骤S24,利用光伏功率累计概率分布的条件概率分布函数计算,对应的分位数,。
步骤S25,通过光伏功率累积概率分布函数求逆计算不同辐照度r取值下光伏功率的上下边界取值和,形成概率功率曲线中的上下两条曲线。
在步骤S21中,光伏功率累计概率分布的条件概率分布函数可通过以下公式计算:
(1);
其中,
在步骤S22中,可通过以下公式计算:
(2)
在步骤S23中,由于光伏功率异常数据的分布不一定均匀,因此设置信区间的不对称系数为,,可通过以下公式计算:
(3)
(4)
当时,置信概率区间对称,当时,置信概率区间向上偏移。
在步骤S24中,分位数,可通过以下公式计算:
(5)
(6)。
可以理解,上述概率功率曲线的获得仅为具体的实施例或表达方式,还可以根据Copula函数描述辐照度和光伏功率这两个随机变量的相关关系及实际需要进行选择,只要能够得出概率功率的分布即可。
在步骤S30中,所述异常数据判别准则可通过观察辐照度-光伏功率散点图,总结异常数据点特征获得,以我国甘肃省某光伏电站实测数据为例,请参阅图3,图3是我国甘肃某光伏电站辐照度-实测功率散点图及异常数据类型示意图,对应的异常数据判别准则如下:
表1光伏功率异常数据判别准则
其中,和是时刻t的总辐照度和光伏功率,和是对应的概率功率曲线上下边界,单位持续时间为T,、、是整数,是允许误差,和分别是需要设定的参数。
在步骤S40中,只有在概率功率曲线以外,且满足异常数据判别准则的数据点才被识别为异常数据。
在步骤S50中,若步骤S40中有异常数据被识别,则在剔除识别的异常数据点后跳转至步骤S10,对新数据集进行步骤S10-S50操作;若步骤S40中无异常数据被识别,则直接利用判别准则和最终形成的概率功率曲线对原始数据集进行异常数据识别。
本发明提供的光伏功率异常数据的识别方法,通过将辐照度和光伏功率作为两个随机变量,利用Copula函数描述两随机变量的相关关系,建立概率功率曲线,根据实测数据的观察建立异常数据判别准则,并针对异常数据比例高的特点,通过反复进行Copula函数参数拟合、概率功率曲线建立、异常数据识别、新数据集建立的操作减少原数据集中异常数据对辐照度-功率之间真实规律发现的干扰。该方法易于操作、具有普适性,不需要复杂的优化计算,切合我国光伏数据以场站级数据为主、异常数据比例高等实际特点。
另外,本领域技术人员还可在本发明精神内作其它变化,当然这些依据本发明精神所作的变化,都应包含在本发明所要求保护的范围内。
Claims (9)
1.一种光伏功率异常数据的识别方法,包括以下步骤:
利用光伏电站辐照度和光伏功率实测数据,拟合辐照度-光伏功率Copula函数参数;
根据Copula函数描述辐照度和光伏功率这两个随机变量的相关关系建立概率功率曲线;
通过辐照度-光伏功率散点图,总结异常数据点特征建立异常数据判别准则;
基于Copula函数和异常数据判别准则,识别筛选异常数据,建立新数据集;以及
若有异常数据被识别,在剔除异常数据后重复以上步骤,对新数据集继续进行识别;若无,则直接利用判别准则和概率功率曲线,识别原始数据集中异常数据。
2.如权利要求1所述的光伏功率异常数据的识别方法,其特征在于,所述Copula函数参数通过以下步骤获得:
获取光伏电站辐照度和功率实测数据,进行数据归一化操作,筛选已标记的错误数据;
利用统计分析得到光伏功率P的累积概率分布函数和辐照度R的累积概率分布函数;
结合辐照度-光伏功率散点图的观察,确定选取的Copula函数类型;
利用和,得到唯一的Copula函数C连接光伏功率P和辐照度R,,拟合Copula函数参数。
3.如权利要求2所述的光伏功率异常数据的识别方法,其特征在于,所述Copula函数为ClaytonCopula,GaussianCopula,FrankCopula函数中的一种。
4.如权利要求2所述的光伏功率异常数据的识别方法,其特征在于,概率功率曲线获取方法,包括以下步骤:
给定辐照度累计概率分布取值,确定光伏功率累计概率分布的条件概率分布函数;
设功率曲线的置信概率为,即有的数据落于概率区间内,的数据落于区间外;
设置信区间的不对称系数为,计算置信区间上下边界的分位概率,,其中表示数据点高于上边界的概率,表示数据点低于下边界的概率;
利用光伏功率累计概率分布的条件概率分布函数计算,对应的分位数,;
通过光伏功率累积概率分布函数求逆计算不同辐照度r取值下光伏功率的上下边界取值和,形成概率功率曲线中的上下两条曲线。
5.如权利要求4所述的光伏功率异常数据的识别方法,其特征在于,光伏功率累计概率分布的条件概率分布函数通过下式获取:
。
6.如权利要求4所述的光伏功率异常数据的识别方法,其特征在于,所述的通过下式获取:
。
7.如权利要求4所述的光伏功率异常数据的识别方法,其特征在于,,通过下式获取:
;
。
8.如权利要求4所述的光伏功率异常数据的识别方法,其特征在于,分位数,通过下式获取:
;
。
9.如权利要求1所述的光伏功率异常数据的识别方法,其特征在于,在概率功率曲线以外,且满足异常数据判别准则的数据点被识别为异常数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510948258.9A CN105590027A (zh) | 2015-12-17 | 2015-12-17 | 光伏功率异常数据的识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510948258.9A CN105590027A (zh) | 2015-12-17 | 2015-12-17 | 光伏功率异常数据的识别方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105590027A true CN105590027A (zh) | 2016-05-18 |
Family
ID=55929603
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510948258.9A Pending CN105590027A (zh) | 2015-12-17 | 2015-12-17 | 光伏功率异常数据的识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105590027A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106169910A (zh) * | 2016-07-20 | 2016-11-30 | 国网青海省电力公司 | 基于群搜索算法的光伏电池参数辨识方法 |
CN107229824A (zh) * | 2017-05-22 | 2017-10-03 | 华北电力科学研究院有限责任公司 | 光伏电站发电单元功率曲线建模方法及装置 |
CN107274021A (zh) * | 2017-06-16 | 2017-10-20 | 南京国电南自电网自动化有限公司 | 一种光伏功率预测系统干扰数据处理系统及方法 |
CN110389949A (zh) * | 2019-07-23 | 2019-10-29 | 华北电力大学 | 一种光伏阵列数据清洗方法 |
CN110555220A (zh) * | 2018-05-31 | 2019-12-10 | 中国电力科学研究院有限公司 | 一种光电转换模型的率定方法及系统 |
CN112085258A (zh) * | 2020-08-13 | 2020-12-15 | 国网上海市电力公司 | 一种基于大数据技术的区域光伏发电量异常实时监测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103559561A (zh) * | 2013-11-13 | 2014-02-05 | 上海电气集团股份有限公司 | 一种光伏电站辐照度的超短期预测方法 |
CN105134484A (zh) * | 2015-07-28 | 2015-12-09 | 国家电网公司 | 风电功率异常数据点的识别方法 |
-
2015
- 2015-12-17 CN CN201510948258.9A patent/CN105590027A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103559561A (zh) * | 2013-11-13 | 2014-02-05 | 上海电气集团股份有限公司 | 一种光伏电站辐照度的超短期预测方法 |
CN105134484A (zh) * | 2015-07-28 | 2015-12-09 | 国家电网公司 | 风电功率异常数据点的识别方法 |
Non-Patent Citations (3)
Title |
---|
叶家宇等: "基于概率分布的异常数据发现与识别算法", 《计算机应用与软件》 * |
朱倩雯等: "风电场输出功率异常数据识别与重构方法研究", 《电力系统保护与控制》 * |
李强等: "基于Copula的我国台湾和韩国股票市场相关性研究", 《管理工程学报》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106169910A (zh) * | 2016-07-20 | 2016-11-30 | 国网青海省电力公司 | 基于群搜索算法的光伏电池参数辨识方法 |
CN106169910B (zh) * | 2016-07-20 | 2018-12-25 | 国网青海省电力公司 | 基于群搜索算法的光伏电池参数辨识方法 |
CN107229824A (zh) * | 2017-05-22 | 2017-10-03 | 华北电力科学研究院有限责任公司 | 光伏电站发电单元功率曲线建模方法及装置 |
CN107229824B (zh) * | 2017-05-22 | 2020-03-13 | 华北电力科学研究院有限责任公司 | 光伏电站发电单元功率曲线建模方法及装置 |
CN107274021A (zh) * | 2017-06-16 | 2017-10-20 | 南京国电南自电网自动化有限公司 | 一种光伏功率预测系统干扰数据处理系统及方法 |
CN110555220A (zh) * | 2018-05-31 | 2019-12-10 | 中国电力科学研究院有限公司 | 一种光电转换模型的率定方法及系统 |
CN110555220B (zh) * | 2018-05-31 | 2022-10-25 | 中国电力科学研究院有限公司 | 一种光电转换模型的率定方法及系统 |
CN110389949A (zh) * | 2019-07-23 | 2019-10-29 | 华北电力大学 | 一种光伏阵列数据清洗方法 |
CN110389949B (zh) * | 2019-07-23 | 2022-12-16 | 华北电力大学 | 一种光伏阵列数据清洗方法 |
CN112085258A (zh) * | 2020-08-13 | 2020-12-15 | 国网上海市电力公司 | 一种基于大数据技术的区域光伏发电量异常实时监测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105590027A (zh) | 光伏功率异常数据的识别方法 | |
CN106643765A (zh) | 一种采集异常维修时长计算方法 | |
CN103869192B (zh) | 智能电网线损检测方法和系统 | |
CN105184394B (zh) | 基于配电网cps在线数据挖掘的优化控制方法 | |
CN110555785B (zh) | 一种月度计划安全稳定校核方法及系统 | |
CN104794544A (zh) | 基于智能算法的分布式光伏发电防窃电监测方法 | |
CN103927695A (zh) | 基于自学习复合数据源的风电功率超短期预测方法 | |
CN101465544B (zh) | 一种用于发电侧控制性能标准考核的实时采集与监控系统 | |
CN103631681A (zh) | 一种在线修复风电场异常数据的方法 | |
CN113011481B (zh) | 一种基于决策树算法的电能表功能异常评估方法及系统 | |
CN105303469A (zh) | 线损异常原因数据挖掘分析的方法与系统 | |
CN104052058A (zh) | 一种基于马尔可夫链蒙特卡罗法的系统谐波概率评估方法 | |
CN104281779A (zh) | 一种异常数据判定与处理方法及装置 | |
CN105488592A (zh) | 一种光伏电站发电量预测的方法 | |
CN102521080A (zh) | 电力用户用电信息采集系统的计算机数据修复方法 | |
CN105141253A (zh) | 基于光伏出力曲线斜率的光伏窃电识别方法 | |
CN108448721A (zh) | 电力系统主站监控信号与标准信号自动匹配的方法 | |
CN103996079A (zh) | 一种基于条件概率的风电功率加权预测方法 | |
CN103887813A (zh) | 基于风功率预测不确定度的风电系统运行的控制方法 | |
CN104123588A (zh) | 智能配电网综合评价系统及方法 | |
CN106780125A (zh) | 一种基于月均用电量的采集异常紧急程度计算方法 | |
CN104574221B (zh) | 一种基于损失电量特征参数的光伏电站运行状态辨识方法 | |
CN109146265A (zh) | 基于动态规划时间序列相似性的三相不平衡源溯源方法 | |
CN101534010B (zh) | 基于微扰的割集电压稳定域局部边界求解方法 | |
CN104376174A (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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160518 |