CN112668611B - 一种基于Kmeans和CEEMD-PE-LSTM的短期光伏发电功率预测方法 - Google Patents
一种基于Kmeans和CEEMD-PE-LSTM的短期光伏发电功率预测方法 Download PDFInfo
- Publication number
- CN112668611B CN112668611B CN202011424687.3A CN202011424687A CN112668611B CN 112668611 B CN112668611 B CN 112668611B CN 202011424687 A CN202011424687 A CN 202011424687A CN 112668611 B CN112668611 B CN 112668611B
- Authority
- CN
- China
- Prior art keywords
- days
- ceemd
- power generation
- kmeans
- photovoltaic 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.)
- Active
Links
- 238000010248 power generation Methods 0.000 title claims abstract description 37
- 238000000034 method Methods 0.000 title claims abstract description 19
- 238000013528 artificial neural network Methods 0.000 claims abstract description 6
- 238000007781 pre-processing Methods 0.000 claims abstract description 4
- 238000004422 calculation algorithm Methods 0.000 claims description 27
- 238000004088 simulation Methods 0.000 claims description 6
- 230000005855 radiation Effects 0.000 claims description 4
- 238000000354 decomposition reaction Methods 0.000 description 7
- 101001095088 Homo sapiens Melanoma antigen preferentially expressed in tumors Proteins 0.000 description 4
- 102100037020 Melanoma antigen preferentially expressed in tumors Human genes 0.000 description 4
- 238000012549 training Methods 0.000 description 4
- 238000007621 cluster analysis Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种基于Kmeans和CEEMD‑PE‑LSTM的短期光伏发电功率预测方法,包括以下步骤:S1、建立基于Kmeans和CEEMD‑PE‑LSTM的光伏发电功率预测模型;S2、对历史功率数据和影响因素数据进行数据预处理;对处理后的数据进行分析,选出待预测日的相似日;S3、对相似日和预测日的影响因素以及发电功率进行分解,分别得到n个固有模态分量以及一个残差量Res。S4、对S3得到的固有模态分量进行计算,将固有模态分量重构成高频分量和低频分量;S5、分别对高频、低频和残差量序列建立长短期记忆神经网络预测模型,预测各子序列功率;S6、将三种子序列预测结果叠加起来,得到光伏功率预测结果。本发明降低了时间序列复杂度,使得子序列具有一定规律和趋势,提高了功率预测模型的准确性。
Description
技术领域
本发明涉及一种光伏发电领域,更具体地,涉及一种基于Kmeans和CEEMD-PE-LSTM的短期光伏发电功率预测方法。
背景技术
通过精确预测光伏发电功率,可以降低光伏并网对电网有可能造成的损害,对电网的稳定运行有重要的意义。目前,光伏发电功率的预测方法大多都是结合天气影响因素和光伏输出功率进行分析预测。有文献中提出EMD-LSTM的光伏发电预测方法,运用EMD将原始序列分别为多个模态,结合LSTM预测的光伏发电功率精度有所提高,但是EMD存在模态混叠的问题。还有文献提出基于Kmeans-SVM的短期光伏发电功率预测模型,利用Kmeans算法对训练样本进行聚类分析找出预测日的相似日,结合SVM算法预测光伏发电功率。另外相关文献还提出一种将相似日和小波神经网络结合起来的光伏功率预测方法,先用灰色关联度算法选取相似日,结合小波神经网络预测光伏发电功率。
以上研究有考虑采用Kmeans算法选取相似日或者使用经验模态分解算法降低时间序列的复杂度,但没有将之结合起来建立预测模型,在预测精度和预测时间上有一点的缺陷,没有能达到理想中的预测效果。
发明内容
本发明针对现有技术中现有技术中光伏发电功率预测精度有待提高,预测时间过长的问题。提供一种基于Kmeans和CEEMD-PE-LSTM的短期光伏发电功率预测方法。
为实现上述目的,本发明采用以下技术方案:
一种基于Kmeans和CEEMD-PE-LSTM的短期光伏发电功率预测方法,包括以下步骤:
S1、建立光伏发电功率预测模型;
S2、对历史功率数据和影响因素数据进行数据预处理;对处理后的数据进行分析,选出待预测日的相似日;
S3、对相似日和预测日的影响因素以及发电功率进行分解,分别得到n个固有模态分量以及一个残差量Res;
S4、对S3得到的固有模态分量进行计算,将固有模态分量重构成高频分量和低频分量;
S5、分别对高频、低频和残差量序列建立长短期记忆神经网络预测模型,预测各子序列功率;
S6、将三种子序列预测结果叠加起来,得到光伏功率预测结果。
进一步地,在步骤S1中建立的光伏发电预测模型包括Kmeans算法、CEEMD算法、PE算法、LSTM算法。
进一步地,在步骤S2采用Kmeans算法对数据进行预处理。
进一步的,在步骤S3中运用CEEMD算法对相似日和预测日的影响因素以及发电功率进行分解。
进一步的,在步骤S4中使用PE算法将固有模态分量重构成高频分量和低频分量。
进一步的,在步骤S5和S6中采用LSTM算法预测各子序列功率并且将子序列预测结果叠加起来,得到光伏功率预测结果。
进一步的,相似日和预测日的影响因素包括辐射强度、温度、湿度、风速。
进一步的,用仿真实验来验证本方法的有效性和精确度;仿真实验中采集的数据间隔为15-30分钟。
进一步的,在仿真实验中选取晴天和雨天的历史数据作为初始样本。
本发明的有益效果为:相比于单一算法的预测结果,本发明建立多种算法结合的光伏发电预测模型,降低了时间序列的复杂度,又使得子序列具有一定规律和趋势,提高了功率预测模型的准确性和精度,减少了预测所需的时间,同时在不同的天气类型下,具有良好的稳定性和适应性。
附图说明
图1为光伏发电功率预测流程图;
图2为光伏发电功率CEEMD分解图;
图3为重构后的时间序列图;
图4为晴天预测结果图;
图5为雨天预测结果图。
具体实施方式
下面结合具体实施方式对本发明作进一步的说明。
一种基于Kmeans和CEEMD-PE-LSTM的短期光伏发电功率预测方法,包括以下步骤:
S1、建立光伏发电功率预测模型;
S2、对历史功率数据和影响因素数据进行数据预处理;对处理后的数据进行分析,选出待预测日的相似日;
S3、对相似日和预测日的影响因素以及发电功率进行分解,分别得到n个固有模态分量以及一个残差量Res;
S4、对S3得到的固有模态分量进行计算,将固有模态分量重构成高频分量和低频分量;
S5、分别对高频、低频和残差量序列建立长短期记忆神经网络预测模型,预测各子序列功率;
S6、将三种子序列预测结果叠加起来,得到光伏功率预测结果。
进一步地,在步骤S1中建立的光伏发电预测模型包括Kmeans算法、CEEMD算法、PE算法、LSTM算法。
进一步地,在步骤S2采用Kmeans算法对数据进行预处理。
进一步的,在步骤S3中运用CEEMD算法对相似日和预测日的影响因素以及发电功率进行分解。
进一步的,在步骤S4中使用PE算法将固有模态分量重构成高频分量和低频分量。
进一步的,在步骤S5和S6中采用LSTM算法预测各子序列功率并且将子序列预测结果叠加起来,得到光伏功率预测结果。
算例分析:
选取某小型光伏发电站的真实数据为例,采用2016年7月到2016年9月的数据作为样本,由于夜晚光伏发电基本为零,所以只考虑每天8:00到18:00的数据,采集间隔为15min的时间序列。将辐射强度、温度、湿度和风速等影响因素作为自变量,光伏发电功率作为因变量。选择晴天9月30日以及雨天9月28日作为预测日,在其余样本中选择晴天和雨天的历史数据分别作为两种天气下的初始训练样本。先采用Kmeans聚类找相似日,再进行CEEMD分解以及PE重构,最后送入LSTM预测模型中进行预测。
首先分别对两种类型天气的样本进行Kmeans聚类分析,将辐射强度、温度、湿度和风速作为输入来进行聚类分析。Kmeans聚类分析后,将与两种类型天气预测日同一类的数据分别作为各自的相似日样本,也就是训练样本,分类结果如下表所示。
得到训练集和预测集的数据后,对晴天和雨天的相似日和预测日数据中的自变量影响因素和因变量发电功率的时间序列进行CEEMD分解。以发电功率为例进行CEEMD分解。发电功率经过CEEMD分解得到7个IMF和1个Ref分量,发电功率分解结果如图3所示。按照频率从高到低排列,从频率最高的IMF1排到频率最低的IMF7,最后是残余项Ref。其他影响因素的分解结果也比较相似,
对分解后的时间序列进行PE计算,以计算出的排列熵值对发电功率子序列进行重构。分别计算出各子序列的排列熵值,并且进行归一化,如下表所示:
将子序列进行重构,PE值大于0.6,表明时间序列比较随机,PE值小于0.6,时间序列比较规则。所以,将IMF1到IMF2重构为高频序列,将IMF3到IMF7重构为低频信号。重构后的时间序列如图4所示。
将所有影响因素和发电功率的时间序列进行分解重构后,对高频、低频和残余量序列分别建立LSTM预测模型进行功率预测并且进行叠加,得到光伏功率预测结果。
为了验证本文模型对光伏发电功率预测的有效性,还另外建立了单独的LSTM模型和EMD-PE-LSTM模型,将预测结果来进行比较。采用三种模型分别预测9月30日(晴天)以及9月28日(雨天),预测结果分别如图4、图5所示。
可以看出,晴天天气下,发电功率从始而终都很平滑,三种模型的预测结果都与实际值很接近,预测精度都很高。雨天天气下,发电功率波动较大,三种模型的预测精度明显低于晴天的预测精度,预测值与实际值有较大误差,但CMEED-PE-LSTM模型在三种模型的预测结果中效果最好。
分别计算晴天和雨天MAPE和RMSE,得到误差指标如下表所示。
晴天预测精度整体比雨天预测精度高,这是因为,晴天的光伏发电功率比较平滑,没有大的波动,预测效果很好;雨天的光伏发电功率较大波动,预测难度较高,随机性大,在有波动的地方,预测不准,从而影响整体预测精度。在两种天气条件下,本文模型的MAPE和RMSE值均最小。晴天条件下,CEEMD-PE-LSTM模型的MAPE比LSTM模型的预测精度高2%左右,比EMD-PE-LSTM模型的预测精度也所提高,RMSE为2.23kW,在三种模型中也是最小的;雨天CEEMD-PE-LSTM模型的MAPE为5.10%,比LSTM、EMD-PE-LSTM分别低1.35%和4.98%,RMSE为5.56kW,比其他两种模型的都低。综合看本文模型可以很好预测光伏发电功率,不管是光伏发电功率平稳时还是有较大波动时,预测精度都有了明显提高。
以上所述实施例仅是为充分说明本发明而所举的较佳的实施例,本发明的保护范围不限于此。本技术领域的技术人员在本发明基础上所作的等同替代或变换,均在本发明的保护范围之内。
Claims (4)
1.一种基于Kmeans和CEEMD-PE-LSTM的短期光伏发电功率预测方法,其特征在于,包括以下步骤:
S1、建立光伏发电功率预测模型;
S2、采用Kmeans算法对历史功率数据和影响因素数据进行数据预处理;对处理后的数据进行分析,选出待预测日的相似日;
S3、运用CEEMD算法对相似日和预测日的影响因素以及发电功率进行分解,分别得到n个固有模态分量以及一个残差量Res;
S4、对S3得到的固有模态分量进行计算,使用PE算法将固有模态分量重构成高频分量和低频分量;
S5、分别对高频、低频和残差量序列建立长短期记忆神经网络预测模型,预测各子序列功率;
S6、将三种子序列预测结果叠加起来,得到光伏功率预测结果。
2.根据权利要求1所述的一种基于Kmeans和CEEMD-PE-LSTM的短期光伏发电功率预测方法,其特征在于,相似日和预测日的影响因素包括辐射强度、温度、湿度、风速。
3.根据权利要求1所述的一种基于Kmeans和CEEMD-PE-LSTM的短期光伏发电功率预测方法,其特征在于,用仿真实验来验证本方法的有效性和精确度;仿真实验中采集的数据间隔为15-30分钟。
4.根据权利要求3所述的一种基于Kmeans和CEEMD-PE-LSTM的短期光伏发电功率预测方法,其特征在于,在仿真实验中选取晴天和雨天的历史数据作为初始样本。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011424687.3A CN112668611B (zh) | 2020-12-08 | 2020-12-08 | 一种基于Kmeans和CEEMD-PE-LSTM的短期光伏发电功率预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011424687.3A CN112668611B (zh) | 2020-12-08 | 2020-12-08 | 一种基于Kmeans和CEEMD-PE-LSTM的短期光伏发电功率预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112668611A CN112668611A (zh) | 2021-04-16 |
CN112668611B true CN112668611B (zh) | 2024-02-02 |
Family
ID=75401527
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011424687.3A Active CN112668611B (zh) | 2020-12-08 | 2020-12-08 | 一种基于Kmeans和CEEMD-PE-LSTM的短期光伏发电功率预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112668611B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112906995A (zh) * | 2021-04-23 | 2021-06-04 | 河北工业大学 | 基于emd-fe-lstm及迭代误差修正的光伏功率短期预测方法 |
CN113205226B (zh) * | 2021-05-28 | 2022-10-14 | 河北工业大学 | 一种结合注意力机制以及误差修正的光伏功率预测方法 |
CN116167465A (zh) * | 2023-04-23 | 2023-05-26 | 杭州经纬信息技术股份有限公司 | 基于多变量时间序列集成学习的太阳辐照度预测方法 |
CN116468469A (zh) * | 2023-04-24 | 2023-07-21 | 北京化工大学 | 一种基于改进emd-lstm的网约车需求预测方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107273426A (zh) * | 2017-05-18 | 2017-10-20 | 四川新网银行股份有限公司 | 一种基于深度语义路径搜索的短文本聚类方法 |
CN109002915A (zh) * | 2018-07-13 | 2018-12-14 | 福州大学 | 基于Kmeans-GRA-Elman模型的光伏电站短期功率预测方法 |
CN110380424A (zh) * | 2019-06-19 | 2019-10-25 | 湖南工业大学 | 电网电压不平衡下级联statcom改进正负序双环控制策略 |
CN110738363A (zh) * | 2019-09-30 | 2020-01-31 | 华中科技大学 | 一种光伏发电功率预测模型及其构建方法和应用 |
CN111241466A (zh) * | 2020-01-15 | 2020-06-05 | 上海海事大学 | 一种基于深度学习的船舶流量预测方法 |
CN111504635A (zh) * | 2020-04-21 | 2020-08-07 | 哈尔滨理工大学 | 基于差分进化概率神经网络的行星齿轮故障诊断方法 |
CN111695290A (zh) * | 2020-05-14 | 2020-09-22 | 天津大学 | 一种适用于变化环境下的短期径流智能预报混合模型方法 |
CN111967660A (zh) * | 2020-08-05 | 2020-11-20 | 河海大学 | 一种基于svr的超短期光伏预测残差修正方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20170286838A1 (en) * | 2016-03-29 | 2017-10-05 | International Business Machines Corporation | Predicting solar power generation using semi-supervised learning |
-
2020
- 2020-12-08 CN CN202011424687.3A patent/CN112668611B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107273426A (zh) * | 2017-05-18 | 2017-10-20 | 四川新网银行股份有限公司 | 一种基于深度语义路径搜索的短文本聚类方法 |
CN109002915A (zh) * | 2018-07-13 | 2018-12-14 | 福州大学 | 基于Kmeans-GRA-Elman模型的光伏电站短期功率预测方法 |
CN110380424A (zh) * | 2019-06-19 | 2019-10-25 | 湖南工业大学 | 电网电压不平衡下级联statcom改进正负序双环控制策略 |
CN110738363A (zh) * | 2019-09-30 | 2020-01-31 | 华中科技大学 | 一种光伏发电功率预测模型及其构建方法和应用 |
CN111241466A (zh) * | 2020-01-15 | 2020-06-05 | 上海海事大学 | 一种基于深度学习的船舶流量预测方法 |
CN111504635A (zh) * | 2020-04-21 | 2020-08-07 | 哈尔滨理工大学 | 基于差分进化概率神经网络的行星齿轮故障诊断方法 |
CN111695290A (zh) * | 2020-05-14 | 2020-09-22 | 天津大学 | 一种适用于变化环境下的短期径流智能预报混合模型方法 |
CN111967660A (zh) * | 2020-08-05 | 2020-11-20 | 河海大学 | 一种基于svr的超短期光伏预测残差修正方法 |
Non-Patent Citations (1)
Title |
---|
基于改进EMD-ARIMA的光伏发电系统短期功率预测;仇琦;杨兰;丁旭;董志强;苏然;郑凌蔚;;电力科学与工程(第08期);第46-54页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112668611A (zh) | 2021-04-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112668611B (zh) | 一种基于Kmeans和CEEMD-PE-LSTM的短期光伏发电功率预测方法 | |
CN110610280B (zh) | 一种电力负荷短期预测方法、模型、装置及系统 | |
Raza et al. | An ensemble framework for day-ahead forecast of PV output power in smart grids | |
CN109272156B (zh) | 一种超短期风电功率概率预测方法 | |
CN107292446B (zh) | 一种基于考虑分量关联性小波分解的混合风速预测方法 | |
Dou et al. | Hybrid model for renewable energy and loads prediction based on data mining and variational mode decomposition | |
Rodríguez et al. | Forecasting intra-hour solar photovoltaic energy by assembling wavelet based time-frequency analysis with deep learning neural networks | |
Yin et al. | A novel asexual-reproduction evolutionary neural network for wind power prediction based on generative adversarial networks | |
Wang et al. | Improved prediction method of PV output power based on optimised chaotic phase space reconstruction | |
CN115034485A (zh) | 一种基于数据空间的风电功率区间预测方法及装置 | |
Wang et al. | A novel wind power prediction model improved with feature enhancement and autoregressive error compensation | |
Sun et al. | Short-term power load prediction based on VMD-SG-LSTM | |
CN111950752A (zh) | 光伏电站发电量预测方法、装置、系统及其存储介质 | |
Zhang et al. | Prediction method of line loss rate in low‐voltage distribution network based on multi‐dimensional information matrix and dimensional attention mechanism‐long‐and short‐term time‐series network | |
CN112348236A (zh) | 用于智能用电监测终端的异常日负荷需求预测系统与方法 | |
CN115713044B (zh) | 一种多工况切换下的机电设备剩余寿命分析方法和装置 | |
CN116995670A (zh) | 基于多模式分解和多分支输入的光伏功率超短期预测方法 | |
CN111539558A (zh) | 一种采用优化极限学习机的电力负荷预测方法 | |
CN116596129A (zh) | 一种电动汽车充电场站短期负荷预测模型构建方法 | |
CN113761777B (zh) | 一种基于hp-ovmd的超短期光伏功率预测方法 | |
CN111126645A (zh) | 一种基于数据挖掘技术和改进支持向量机的风电预测算法 | |
CN113919600A (zh) | 一种居民负荷超短期预测方法 | |
CN114298408A (zh) | 一种基于ceemd-lstm-mlr的短期电力负荷预测方法 | |
CN114386666A (zh) | 一种基于时空相关性的风电场短期风速预测方法 | |
Huang et al. | Scenario Reduction of Power Systems with Renewable Generations Using Improved Time-GAN |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB03 | Change of inventor or designer information |
Inventor after: Yu Huijun Inventor after: Li Bingchen Inventor after: Liu Jingyu Inventor before: Yu Huijun Inventor before: Liu Jingyu |
|
CB03 | Change of inventor or designer information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |