CN105589998A - 一种光伏输出功率超短期预测方法 - Google Patents
一种光伏输出功率超短期预测方法 Download PDFInfo
- Publication number
- CN105589998A CN105589998A CN201510980687.4A CN201510980687A CN105589998A CN 105589998 A CN105589998 A CN 105589998A CN 201510980687 A CN201510980687 A CN 201510980687A CN 105589998 A CN105589998 A CN 105589998A
- Authority
- CN
- China
- Prior art keywords
- power output
- photovoltaic power
- screening
- historical data
- value
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
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)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种光伏输出功率超短期预测方法,包括下述步骤:A、使用经验模态分解法对光伏输出功率历史数据进行分解,得到本征模函数和剩余分量;B、对所有本征模函数和剩余分量分别建立最小二乘支持向量机预测模型,进行光伏输出功率超短期预测,得到预测值;C、将步骤B最后所得所有预测值相加,得到最终预测结果。利用本发明的光伏输出功率超短期预测方法,能够充分利用支持向量机良好的回归特性、小样本适用性、不易陷入局部最优特性、计算速度快特性以及经验模态分解法对非线性数据良好的平稳化处理能力。
Description
技术领域
本发明涉及电力系统分析技术领域,尤其涉及一种新能源电力系统分析技术。
背景技术
能源和环境问题促使以光伏为代表的清洁能源获得了广泛关注。由于光照和光伏发电的自身特性,使得光伏发电呈现波动性和间歇性的特点。光伏大规模并网后会对电网的安全、经济和可靠运行产生多方面影响,因此必须获得光伏出力的超短期预测数值,从而调节常规电源实时出力曲线,达到优化调度。
但是,现有技术中,对于光伏出力超短期预测的技术方案较少,技术效果也具有较大局限,例如使用统计学原理的模型精度较低。造成此结果的一个原因是光伏输出功率具有强烈的随机波动特性,规律难以把握。
发明内容
鉴于此,本发明的目的在于克服现有技术的缺点,提出一种光伏输出功率超短期预测方法,首先使用经验模态分解法对光伏出力历史序列进行平稳化预处理,使得序列规律更容易被学习,然后使用具有良好回归特性的最小二乘支持向量机进行预测,从而获得更高预测精度。
为了实现此目的,本发明采取的技术方案为如下。
一种光伏输出功率超短期预测方法,包括下述步骤:
A、使用经验模态分解法对光伏输出功率历史数据进行分解,得到本征模函数和剩余分量;
B、对所有本征模函数和剩余分量分别建立最小二乘支持向量机预测模型,进行光伏输出功率超短期预测,得到预测值;
C、将步骤B最后所得所有预测值相加,得到最终预测结果。
步骤A中使用经验模态分解法进行分解包括:
A1、识别光伏输出功率历史数据序列x(t)中所有极大值点并使用三次样条插值法拟合,得到上包络线eup(t),
A2、识别光伏输出功率历史数据序列x(t)中所有极小值点并使用三次样条插值法拟合,得到下包络线elow(t),
A3、根据计算均值包络线即上下包络线的平均值m1(t),
A4、根据h1(t)=x(t)-m1(t)计算h1(t),判断h1(t)是否满足筛分终止条件,若满足,则h1(t)为第一个本征模函数,记为c1(t)=h1(t),若h1(t)不满足筛分终止条件,则将h1(t)视为新的序列x(t),从步骤A1开始重新筛分,直到选出满足筛分终止条件的本征模函数,记为c1(t)=h1(t),
A5、根据r1(t)=x(t)-c1(t)计算r1(t),若r1(t)极值点数目大于等于2个,则将r1(t)作为新的序列x(t),从步骤A1开始重新筛分,若r1(t)极值点数目不足2个,则停止分解,r1(t)作为剩余分量予以保留,
最终,光伏输出功率历史数据序列被分解为本征模函数和剩余分量,表示如下:
其中ci(t)为第i个本征模函数,rN(t)为剩余分量。
步骤B所述的最小二乘支持向量机预测模型,它的回归方式为:
其中:为模型输出的出力预测值,ω为权向量,为映射函数,为训练样本和测试样本,b为偏置,最小二乘支持向量机中使用径向基核函数,共包含两个未知参数需要进行寻优:正则化参数γ和平方带宽σ2;
使用遗传方法,得到γ和σ2的最佳取值。
其中,所述筛分终止条件为:
设此时进行第k次筛分,则当σ小于给定数值σmax时,筛分终止:
其中h(j,k-1)(t)和h(j,k)(t)是两个连续筛分的结果序列,n为光伏输出功率历史数据序列x(t)样本点个数,
σmax取值为0.1。
通过采用本发明的光伏输出功率超短期预测方法,能够获得以下有益技术效果:首先采用经验模态分解法对光伏出力历史序列进行预处理,降低其波动性和随机性,使其分量变得更平稳更容易被学习。然后使用最小二乘支持向量机模型作为主体预测模型,充分利用了支持向量机良好的回归特性、小样本适用性、不易陷入局部最优特性、计算速度快特性等。以上特性使本发明所提模型适宜应用于光伏超短期预测领域,并且实际算例证明,该模型具有较高精度以及相对其他模型而言的较大优势。
附图说明
图1为本发明具体实施方式中光伏输出功率超短期预测方法的流程示意图。
图2为本发明具体实施方式中光伏出力样本数据经经验模态分解后的结果示意图。
图3本发明具体实施方式中EMD-LSSVM-GA法预测结果示意图。
图4本发明具体实施方式中LSSVM-GA法预测结果示意图。
具体实施方式
下面结合附图,对本发明作详细说明。
以下公开详细的示范实施例。然而,此处公开的具体结构和功能细节仅仅是出于描述示范实施例的目的。
然而,应该理解,本发明不局限于公开的具体示范实施例,而是覆盖落入本公开范围内的所有修改、等同物和替换物。在对全部附图的描述中,相同的附图标记表示相同的元件。
同时应该理解,如在此所用的术语“和/或”包括一个或多个相关的列出项的任意和所有组合。另外应该理解,当部件或单元被称为“连接”或“耦接”到另一部件或单元时,它可以直接连接或耦接到其他部件或单元,或者也可以存在中间部件或单元。此外,用来描述部件或单元之间关系的其他词语应该按照相同的方式理解(例如,“之间”对“直接之间”、“相邻”对“直接相邻”等)。
如图1所示,本发明具体实施方式中公开了一种光伏输出功率超短期预测方法,它包括下述步骤:
使用经验模态分解法对光伏输出功率历史数据进行分解,得到本征模函数和剩余分量;
B、对所有本征模函数和剩余分量分别建立最小二乘支持向量机预测模型,进行光伏输出功率超短期预测,得到预测值;
C、将步骤B最后所得所有预测值相加,得到最终预测结果。
因此,本发明的光伏输出功率超短期预测方法首先采用经验模态分解法对光伏出力历史序列进行预处理,降低其波动性和随机性,使其分量变得更平稳更容易被学习。然后使用最小二乘支持向量机模型作为主体预测模型,充分利用了支持向量机良好的回归特性、小样本适用性、不易陷入局部最优特性、计算速度快特性等。
特别地,在具体实施方式中,步骤A中使用经验模态分解法进行分解包括:
A1、识别光伏输出功率历史数据序列x(t)中所有极大值点并使用三次样条插值法拟合,得到上包络线eup(t),
A2、识别光伏输出功率历史数据序列x(t)中所有极小值点并使用三次样条插值法拟合,得到下包络线elow(t),
A3、根据计算均值包络线即上下包络线的平均值m1(t),
A4、根据h1(t)=x(t)-m1(t)计算h1(t),判断h1(t)是否满足筛分终止条件,若满足,则h1(t)为第一个本征模函数,记为c1(t)=h1(t),若h1(t)不满足筛分终止条件,则将h1(t)视为新的序列x(t),从步骤A1开始重新筛分,直到选出满足筛分终止条件的本征模函数,记为c1(t)=h1(t),
A5、根据r1(t)=x(t)-c1(t)计算r1(t),若r1(t)极值点数目大于等于2个,则将r1(t)作为新的序列x(t),从步骤A1开始重新筛分,若r1(t)极值点数目不足2个,则停止分解,r1(t)作为剩余分量予以保留,
最终,光伏输出功率历史数据序列被分解为本征模函数和剩余分量,表示如下:
其中ci(t)为第i个本征模函数,rN(t)为剩余分量。
另外,步骤B所述的最小二乘支持向量机预测模型,它的回归方式为:
其中:为模型输出的出力预测值,ω为权向量,为映射函数,为训练样本和测试样本,b为偏置,最小二乘支持向量机中使用径向基核函数,共包含两个未知参数需要进行寻优:正则化参数γ和平方带宽σ2;
使用遗传方法,得到γ和σ2的最佳取值。
特别地,所述筛分终止条件为:
设此时进行第k次筛分,则当σ小于给定数值σmax时,筛分终止:
其中h(j,k-1)(t)和h(j,k)(t)是两个连续筛分的结果序列,n为光伏输出功率历史数据序列x(t)样本点个数,
σmax取值为0.1。
下面结合以一个更加具体的示例对本发明作进一步的描述:
步骤1:获取样本数据。获取中国西南地区某60kW光伏电站2012年的出力历史数据,该数据采样间隔为15min。选取2012年10月1号~3号每日08:00至17:45的出力数据,共计120个。
步骤2:使用经验模态分解法对步骤1获取的样本数据进行分解,得到5个分量:4个本征模函数和1个剩余分量。分解结果见图2。
步骤3:构造训练样本和测试样本。
设序列{xi1,xi2,xi3,...xin}(n为数据个数)为分解出的lp个分量中的第i个分量。为满足建模要求,将序列{xi1,xi2,xi3,...xin}转化为符合LSSVM(最小二乘向量机)输入格式的矩阵形式(Xt,Yt),其中Xt={xt-m,xt-m+1,...,xt-1},Yt=xt,t=1,2,...,n。m为滑动窗口大小,取为3,表示根据前m个出力值推断第m+1个出力值大小。本发明实例中,lp为5,n为120。此时共得到117组(Xt,Yt)。将前85组(Xt,Yt)作为训练样本,后32组(Xt,Yt)作为测试样本。
步骤4:使用遗传方法进行参数寻优,为5个分量分别建立最小二乘支持向量机预测模型。遗传方法寻优的具体设置为:种群规模为200,最大进化代数为500,交叉概率为0.4,变异概率为0.01,参数γ的搜索范围为[1,20],σ2的搜索范围为[1,100]。
步骤5:结果合成。使用步骤4寻优得到的参数构建最小二乘支持向量机预测模型,使用步骤3构造的训练样本和测试样本进行预测,每个分量各得到一个预测结果。将5个结果相加,得到最终预测结果,如图2所示。
步骤6:误差评价。以上方法简称为EMD-LSSVM-GA法。将经验模态分解环节去掉后,单纯使用LSSVM-GA法对同样样本进行预测,预测结果见图3。比较两种方法的误差大小。选取均方根误差(Root-mean-squareError,RMSE)作为衡量预测结果的指标。RMSE的计算方法为:
其中,s为预测样本数,Pj‘为第j’个预测样本的实测值,P’j’为第j’个预测样本的预测值,Cap为光伏电站容量。
经计算,EMD-LSSVM-GA法的RMSE为2.38%,LSSVM-GA法的RMSE为3.34%。因此,本发明提出的方法具有较高预测精度。同时观察图形也可看出,EMD-LSSVM-GA法的预测结果能很好地跟踪出力变化趋势,在平缓和拐角处都有较好的拟合能力。LSSVM-GA模型在拐点处出现几次较大误差,容易导致决策人员的错误判断,且滞后现象较严重。因此EMD-LSSVM-GA法有较高应用价值。
需要说明的是,上述实施方式仅为本发明较佳的实施方案,不能将其理解为对本发明保护范围的限制,在未脱离本发明构思前提下,对本发明所做的任何微小变化与修饰均属于本发明的保护范围。
Claims (4)
1.一种光伏输出功率超短期预测方法,包括下述步骤:
A、使用经验模态分解法对光伏输出功率历史数据进行分解,得到本征模函数和剩余分量;
B、对所有本征模函数和剩余分量分别建立最小二乘支持向量机预测模型,进行光伏输出功率超短期预测,得到预测值;
C、将步骤B最后所得所有预测值相加,得到最终预测结果。
2.根据权利要求1中所述的光伏输出功率超短期预期方法,其特征在于,步骤A中使用经验模态分解法进行分解包括:
A1、识别光伏输出功率历史数据序列x(t)中所有极大值点并使用三次样条插值法拟合,得到上包络线eup(t),
A2、识别光伏输出功率历史数据序列x(t)中所有极小值点并使用三次样条插值法拟合,得到下包络线elow(t),
A3、根据计算均值包络线即上下包络线的平均值m1(t),
A4、根据h1(t)=x(t)-m1(t)计算h1(t),判断h1(t)是否满足筛分终止条件,若满足,则h1(t)为第一个本征模函数,记为c1(t)=h1(t),若h1(t)不满足筛分终止条件,则将h1(t)视为新的序列x(t),从步骤A1开始重新筛分,直到选出满足筛分终止条件的本征模函数,记为c1(t)=h1(t),
A5、根据r1(t)=x(t)-c1(t)计算r1(t),若r1(t)极值点数目大于等于2个,则将r1(t)作为新的序列x(t),从步骤A1开始重新筛分,若r1(t)极值点数目不足2个,则停止分解,r1(t)作为剩余分量予以保留,
最终,光伏输出功率历史数据序列被分解为本征模函数和剩余分量,表示如下:
其中ci(t)为第i个本征模函数,rN(t)为剩余分量。
3.根据权利要求1中所述的光伏输出功率超短期预测方法,其特征在于,步骤B所述的最小二乘支持向量机预测模型,它的回归方式为:
其中:为模型输出的出力预测值,ω为权向量,为映射函数,为训练样本和测试样本,b为偏置,最小二乘支持向量机中使用径向基核函数,共包含两个未知参数需要进行寻优:正则化参数γ和平方带宽σ2;
使用遗传方法,得到γ和σ2的最佳取值。
4.根据权利要求2中所述的光伏输出功率超短期预测方法,其特征在于,所述筛分终止条件为:
设此时进行第k次筛分,则当σ小于给定数值σmax时,筛分终止:
其中h(j,k-1)(t)和h(j,k)(t)是两个连续筛分的结果序列,n为光伏输出功率历史数据序列x(t)样本点个数,
σmax取值为0.1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510980687.4A CN105589998A (zh) | 2015-12-23 | 2015-12-23 | 一种光伏输出功率超短期预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510980687.4A CN105589998A (zh) | 2015-12-23 | 2015-12-23 | 一种光伏输出功率超短期预测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105589998A true CN105589998A (zh) | 2016-05-18 |
Family
ID=55929576
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510980687.4A Pending CN105589998A (zh) | 2015-12-23 | 2015-12-23 | 一种光伏输出功率超短期预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105589998A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106656032A (zh) * | 2016-08-29 | 2017-05-10 | 华北电力大学(保定) | 针对雾霾积灰的光伏发电功率输出减少率估计方法 |
CN109891518A (zh) * | 2016-10-28 | 2019-06-14 | 皇家飞利浦有限公司 | 时间敏感的风险模型计算 |
CN110516844A (zh) * | 2019-07-25 | 2019-11-29 | 太原理工大学 | 基于emd-pca-lstm的多变量输入光伏功率预测方法 |
CN112651537A (zh) * | 2019-10-10 | 2021-04-13 | 国网河北省电力有限公司 | 一种光伏发电超短期功率预测方法及系统 |
CN113361803A (zh) * | 2021-06-28 | 2021-09-07 | 嘉兴国电通新能源科技有限公司 | 基于生成对抗网络的超短期光伏功率预测方法 |
CN113469442A (zh) * | 2021-07-05 | 2021-10-01 | 国网山东省电力公司东营供电公司 | 基于经验模态分解的光伏功率超短期预测方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103049798A (zh) * | 2012-12-05 | 2013-04-17 | 浙江大学城市学院 | 一种应用于光伏发电系统的短期发电功率预测方法 |
CN103390902A (zh) * | 2013-06-04 | 2013-11-13 | 国家电网公司 | 一种基于最小二乘法的光伏电站超短期功率预测方法 |
-
2015
- 2015-12-23 CN CN201510980687.4A patent/CN105589998A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103049798A (zh) * | 2012-12-05 | 2013-04-17 | 浙江大学城市学院 | 一种应用于光伏发电系统的短期发电功率预测方法 |
CN103390902A (zh) * | 2013-06-04 | 2013-11-13 | 国家电网公司 | 一种基于最小二乘法的光伏电站超短期功率预测方法 |
Non-Patent Citations (3)
Title |
---|
CUI DU等: "The Real-time Forecast of Photovoltaic Output Based on Empirical Mode Decomposition", 《PROCEEDINGS OF THE 2015 2ND INTERNATIONAL WORKSHOP ON MATERIALS ENGINEERING AND COMPUTER SCIENCES:ADVANCES IN COMPUTER SCIENCE RESEARCH》 * |
J TERRADAS等: "Application of statistical techniques to the analysis of solar coronal oscillations", 《ASTROPHYSICAL JOURNAL》 * |
杜翠等: "气象数据弱相关的光伏出力短期预测", 《现代电力》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106656032A (zh) * | 2016-08-29 | 2017-05-10 | 华北电力大学(保定) | 针对雾霾积灰的光伏发电功率输出减少率估计方法 |
CN106656032B (zh) * | 2016-08-29 | 2018-08-07 | 华北电力大学(保定) | 针对雾霾积灰的光伏发电功率输出减少率估计方法 |
CN109891518A (zh) * | 2016-10-28 | 2019-06-14 | 皇家飞利浦有限公司 | 时间敏感的风险模型计算 |
CN110516844A (zh) * | 2019-07-25 | 2019-11-29 | 太原理工大学 | 基于emd-pca-lstm的多变量输入光伏功率预测方法 |
CN112651537A (zh) * | 2019-10-10 | 2021-04-13 | 国网河北省电力有限公司 | 一种光伏发电超短期功率预测方法及系统 |
CN113361803A (zh) * | 2021-06-28 | 2021-09-07 | 嘉兴国电通新能源科技有限公司 | 基于生成对抗网络的超短期光伏功率预测方法 |
CN113469442A (zh) * | 2021-07-05 | 2021-10-01 | 国网山东省电力公司东营供电公司 | 基于经验模态分解的光伏功率超短期预测方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105589998A (zh) | 一种光伏输出功率超短期预测方法 | |
Sun et al. | Short-term wind power forecasting based on VMD decomposition, ConvLSTM networks and error analysis | |
Miller et al. | Framework for analyzing transformative technologies in life cycle assessment | |
De Giorgi et al. | Comparison of strategies for multi-step ahead photovoltaic power forecasting models based on hybrid group method of data handling networks and least square support vector machine | |
CN103279807B (zh) | 一种恶劣天气下电网静态风险评估方法 | |
CN105425123A (zh) | 一种紫外成像与红外成像协同检测电力设备故障的方法及系统 | |
Xu et al. | A decision tree-based on-line preventive control strategy for power system transient instability prevention | |
CN104601109A (zh) | 一种电-图模型的光伏热斑效应检测方法 | |
Azzopardi et al. | Decision support system for ranking photovoltaic technologies | |
CN103077300A (zh) | 一种基于二型模糊逻辑的分布式光伏电源发电量预测方法 | |
CN106845562B (zh) | 光伏组件的故障监测系统及数据处理方法 | |
Zhang et al. | Intelligent systems for power system dynamic security assessment: Review and classification | |
CN103886223B (zh) | 功率预测方法和系统 | |
CN104332998A (zh) | 一种电力系统直流紧急功率调制改善频率安全的控制性能量化评价指标计算方法 | |
Tannahill et al. | Modeling of system of systems via data analytics—Case for “Big Data” in SoS | |
Weckesser et al. | An improved on-line contingency screening for power system transient stability assessment | |
Rasheed et al. | Asymmetric relationship between competitive industrial performance, renewable energy, industrialization, and carbon footprint: Does artificial intelligence matter for environmental sustainability? | |
CN110942180B (zh) | 基于xgboost算法的工业设计匹配服务方预测方法 | |
Yamaguchi et al. | Real-time prediction for future profile of car travel based on statistical data and greedy algorithm | |
CN102761888B (zh) | 一种基于特征选择的传感网络异常检测方法和装置 | |
Acquaye et al. | Green supply chain management: the development of supply chain carbon maps | |
Siddula et al. | Optimal Placement of Hybrid Wind‐Solar System Using Deep Learning Model | |
Miraftabzadeh et al. | Exploring the Synergy of Artificial Intelligence in Energy Storage Systems for Electric Vehicles | |
Xie et al. | Walmart Sale Forecasting Model Based On LSTM And LightGBM | |
CN110310033B (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 |