CN113722375B - 一种基于多算法优化的双层混合光伏预测方法及计算机设备 - Google Patents
一种基于多算法优化的双层混合光伏预测方法及计算机设备 Download PDFInfo
- Publication number
- CN113722375B CN113722375B CN202110913414.3A CN202110913414A CN113722375B CN 113722375 B CN113722375 B CN 113722375B CN 202110913414 A CN202110913414 A CN 202110913414A CN 113722375 B CN113722375 B CN 113722375B
- Authority
- CN
- China
- Prior art keywords
- photovoltaic
- random
- time sequence
- data
- layer
- 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
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 45
- 238000005457 optimization Methods 0.000 title claims abstract description 25
- 238000012549 training Methods 0.000 claims description 38
- 230000005855 radiation Effects 0.000 claims description 35
- 238000004364 calculation method Methods 0.000 claims description 20
- 238000012937 correction Methods 0.000 claims description 19
- 238000010801 machine learning Methods 0.000 claims description 15
- 238000013528 artificial neural network Methods 0.000 claims description 8
- 238000007621 cluster analysis Methods 0.000 claims description 8
- 238000012216 screening Methods 0.000 claims description 8
- 101001095088 Homo sapiens Melanoma antigen preferentially expressed in tumors Proteins 0.000 claims description 6
- 102100037020 Melanoma antigen preferentially expressed in tumors Human genes 0.000 claims description 6
- 238000001556 precipitation Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 238000010845 search algorithm Methods 0.000 claims description 6
- 230000006870 function Effects 0.000 claims description 4
- 230000007787 long-term memory Effects 0.000 claims description 4
- 230000006403 short-term memory Effects 0.000 claims description 4
- 230000005484 gravity Effects 0.000 claims description 3
- 238000007477 logistic regression Methods 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 230000015654 memory Effects 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 5
- 238000005265 energy consumption Methods 0.000 description 5
- 238000004590 computer program Methods 0.000 description 3
- 238000000354 decomposition reaction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000010220 Pearson correlation analysis Methods 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000306 recurrent effect Effects 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/24—Querying
- G06F16/245—Query processing
- G06F16/2458—Special types of queries, e.g. statistical queries, fuzzy queries or distributed queries
- G06F16/2465—Query processing support for facilitating data mining operations in structured databases
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/24—Querying
- G06F16/245—Query processing
- G06F16/2458—Special types of queries, e.g. statistical queries, fuzzy queries or distributed queries
- G06F16/2474—Sequence data queries, e.g. querying versioned data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/044—Recurrent networks, e.g. Hopfield networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- 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
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate 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
- 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
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Strategic Management (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Probability & Statistics with Applications (AREA)
- General Health & Medical Sciences (AREA)
- Fuzzy Systems (AREA)
- Tourism & Hospitality (AREA)
- Water Supply & Treatment (AREA)
- Entrepreneurship & Innovation (AREA)
- Computing Systems (AREA)
- Primary Health Care (AREA)
- Biomedical Technology (AREA)
- Public Health (AREA)
- Artificial Intelligence (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Biophysics (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Evolutionary Computation (AREA)
- Molecular Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Photovoltaic Devices (AREA)
Abstract
本发明公开了一种基于多算法优化的双层混合光伏预测方法及计算机设备，属于光伏预测技术领域。现有的光伏出力预测方案，准确度不高，无法实现短期光伏功率准确预测。本发明的一种基于多算法优化的双层混合光伏预测方法，对获取的数据信息进行分类和筛选，并将其划分为：时序特征量、随机特征量；并将光伏历史输出功率数据分成时序分量和随机分量，并根据时序特征量形成时序层、随机层；然后叠加时序层和随机层，得到光伏出力预测模型，对光伏出力进行准确预测。本发明能够综合考虑随机因素和时序因素对模型的影响，有效提高模型预测的准确度，实现光伏出力的短期准确预测，以增强电力系统整体稳定性。
Description
技术领域
本发明涉及一种基于多算法优化的双层混合光伏预测方法及计算机设备，属于光伏预测技术领域。
背景技术
当前，利用太阳能的主要问题是其时空规律性不明显，太阳辐射量的大小可以分为两部分，第一部分表示时序部分，第二部分表示随机部分。时序分量是由地球的自转和公转运动决定的，随机分量是由天气等随机因素决定的，这就导致光伏出力也存在时序性和确定性两部分。太阳能的间歇性对光伏电站输出功率的影响很大，增加了电力系统能量调度的难度。
自上个世纪八十年代人工神经网络诞生以来，越来越多的学者将其应用到光伏功率预测的研究上。如今，随着时间序列神经网络，如循环神经网络(RNN)的出现，利用时间序列预测已成为光伏功率预测的一大热门方法。然而，现有的时间序列预测方法基本上都是将光伏出力直接看作时间序列进行预测，而未专门考虑天气等随机因素对光伏出力的影响。
为了提高光伏短期的预测精度，很多学者选择对气象数据进行分类预测，即天气分型。这种分型预测方法有效地减小了相似日下的光伏短期预测误差，从整体上提高了预测精度。然而，常规的天气分型方法通常采取降水，风速，辐射量等随机因素进行分析，缺乏对时序分量的考虑，效果欠佳。
因此，现有的光伏出力预测方案，没有综合考虑随机因素和时序因素对模型的影响，导致模型预测的准确度不高，无法实现短期光伏功率准确预测，进而不利于电网管理人员合理调度和安排电能消纳，影响电力系统整体稳定性。
发明内容
针对现有技术的缺陷，本发明的目的在于提供一种能兼顾考虑随机因素和时序因素的模型，对获取的数据信息进行分类和筛选，并将其划分为：时序特征量、随机特征量；并将光伏历史输出功率数据分成时序分量和随机分量，并根据时序特征量形成时序层、随机层；然后叠加时序层和随机层，得到光伏出力预测模型，对光伏出力进行预测；有效提高模型预测的准确度，有助于电网管理人员合理调度和安排电能消纳，以增强电力系统整体稳定性的基于多算法优化的双层混合光伏预测方法及计算机设备。
为实现上述目的，本发明的技术方案为：
一种基于多算法优化的双层混合光伏预测方法，包括以下步骤：
步骤1：获取某地天文数据、地理数据、气象数据、气象数据以及光伏历史输出功率数据；
步骤2：对步骤1中获取的数据信息进行分类和筛选，并将其划分为：时序特征量、随机特征量；
步骤3：对步骤1的气象数据，利用修正大气清晰度指数k’T和水平面直射比Bd，运用聚类分析算法，将其划分为天气类型1，天气类型2，天气类型3；
并将天气类型2细分为天气类型2-1和天气类型2-2；
步骤4：利用集合经验模态分解EEMD方法，并根据皮尔逊相关性模型将步骤1中的光伏历史输出功率数据分成时序分量和随机分量；
所述皮尔逊相关性模型，用于构建某地光伏输出功率与某地天文、地理、气象指标的关联性，并能根据光伏历史输出功率各信号分量的关联性大小，把各信号分量划分为时序分量或随机分量；
所述关联性为各信号分量与步骤2中时序特征量的相关系数；
步骤5：结合步骤3中天气类型，构建时序训练模型，对步骤4中的时序分量进行训练，形成滑动窗口，得到时序层；
所述时序训练模型为经过引力搜索算法GSA优化的长短期记忆人工神经网络LSTM模型；
步骤6：结合步骤3中天气类型，构建随机训练模型，对步骤4中的随机分量进行训练,得到随机层；
所述随机训练模型为经过贝叶斯优化的机器学习模型；
步骤7：叠加步骤5中时序层和步骤6中随机层，得到光伏出力预测模型；
所述光伏出力预测模型能够根据某地某时的天文数据、地理数据、气象数据，得到某时间段的光伏输出功率。
由于现有的光伏短期预测方法在天气分型上大多忽略了时序因素对模型的影响；现有的时间序列光伏短期预测方法在原理设计上大多忽略了随机部分对模型的影响。
因此本发明综合考虑随机因素和时序因素对模型的影响，对获取的数据信息进行分类和筛选，并将其划分为：时序特征量、随机特征量；并将光伏历史输出功率数据分成时序分量和随机分量，并根据时序特征量形成时序层、随机层；然后叠加时序层和随机层，得到光伏出力预测模型，对光伏出力进行准确预测；有效提高模型预测的准确度，实现光伏出力的短期准确预测，有助于电网管理人员合理调度和安排电能消纳，以增强电力系统整体稳定性。
进一步，本发明兼顾了光伏出力的时间连续性和随机波动性，有效提高预测准确度，可以适用于大多数光伏出力预测情景，预测速度快，适用性强。
作为优选技术措施：
所述步骤1中，获取的数据特征量包括年，月，日，时，本站气压，环境温度，相对湿度，降水，平均风速，水平能见度，总云量，赤纬角，时角，修正大气清晰度指数，地外总辐射，水平面总辐射，散射辐射，法向辐射，反射辐射和日照百分率。
作为优选技术措施：
所述步骤2中对获取的数据进行处理，包括对获取的某地天文、地理和气象数据进行分类，筛选；
其中时序特征量包括年、月、日、时、赤纬角、时角、地外总辐射I0；
所述随机特征量包括本站气压，环境温度，相对湿度，降水，平均风速，水平能见度，总云量，修正大气清晰度指数，水平面总辐射，散射辐射，法向辐射，反射辐射和日照百分率；
同时筛除异常值，其包括筛除太阳高度角小于10°的数据。
作为优选技术措施：
所述步骤3中的天气类型1为恶劣天气，天气类型2为转折天气，天气类型3为晴空天气，天气类型2-1为晴转阴，天气类型2-2为阴转晴。
作为优选技术措施：
所述步骤3中聚类分析算法为二阶聚类和逻辑回归；
所述修正大气清晰度指数k’T，其计算公式为：
式中，kT为修正前的大气清晰度指数，m为大气质量；
所述修正前的大气清晰度指数，其计算公式为：
式中，I表示获取的气象数据中的水平面总辐射，γ为日地距离订正系数，ESC为太阳常数，取值为1367±7W/m2，δ为赤纬角，φ和ω分别为纬度和时角；
所述赤纬角，其计算公式为：
式中，n为所求日期在一年中的日序数；
所述日地距离订正系数，其计算公式为：
式中，n为所求日期在一年中的日序数。
作为优选技术措施：
所述步骤4中，信号分量的个数为12个；
选取与步骤2中时序特征量相关系数最大的信号分量作为光伏历史输出功率数据的时序分量，其余11个信号分量作为光伏历史输出功率数据的随机分量。
作为优选技术措施：
所述步骤5中的引力搜索算法GSA算法选取的超参数变量为训练迭代次数和训练长度；其适应度函数为预测值与实际值的平均绝对百分比误差MAPE。
作为优选技术措施：
所述训练迭代次数为14次，训练长度为6天；
所述平均绝对百分比误差MAPE的计算公式为：
其中N为测试集样本数量，Pf为预测值，Po为实际值。
作为优选技术措施：
所述步骤6中的经过贝叶斯优化的机器学习算法迭代次数为50次，所述随机层为机器学习算法的预测结果，机器学习算法包括装袋Bagging和高斯过程回归GaussianProcess Regression,GPR。
作为优选应用技术措施：
一种计算机设备，包括一个或多个处理器；
存储装置，用于存储一个或多个程序；
当所述一个或多个程序被所述一个或多个处理器执行时，使得所述一个或多个处理器实现上述的一种基于多算法优化的双层混合光伏预测方法。
与现有技术相比，本发明具有以下有益效果：
本发明经过不断探索以及试验，对获取的数据信息进行分类和筛选，并将其划分为：时序特征量、随机特征量；并将光伏历史输出功率数据分成时序分量和随机分量，并根据时序特征量形成时序层、随机层；然后叠加时序层和随机层，得到光伏出力预测模型，对光伏出力进行准确预测。
本发明能够综合考虑随机因素和时序因素对模型的影响，有效提高模型预测的准确度，实现光伏出力的短期准确预测，有助于电网管理人员合理调度和安排电能消纳，以增强电力系统整体稳定性。
附图说明
图1为本发明的方法流程图。
图2为本发明时序特征量与不同信号分量imf之间的相关性示意图。
图3为本发明滑动窗口LSTM示意图；
图4为本发明不同天气类型与直射比、修正清晰度指数的关系示意图。
图5为本发明预测结果和实际值对比示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
相反，本发明涵盖任何由权利要求定义的在本发明的精髓和范围上做的替代、修改、等效方法以及方案。进一步，为了使公众对本发明有更好的了解，在下文对本发明的细节描述中，详尽描述了一些特定的细节部分。对本领域技术人员来说没有这些细节部分的描述也可以完全理解本发明。
一种基于多算法优化的双层混合光伏预测方法，包括以下步骤：
步骤1：获取某地天文数据、地理数据、气象数据、气象数据以及光伏历史输出功率数据；
步骤2：对步骤1中获取的数据信息进行分类和筛选，并将其划分为：时序特征量、随机特征量；
步骤3：对步骤1的气象数据，利用修正大气清晰度指数k’T和水平面直射比Bd，运用聚类分析算法，将其划分为天气类型1，天气类型2，天气类型3；
并将天气类型2细分为天气类型2-1和天气类型2-2；
步骤4：利用集合经验模态分解EEMD方法，并根据皮尔逊相关性模型将步骤1中的光伏历史输出功率数据分成时序分量和随机分量；
所述皮尔逊相关性模型，用于构建某地光伏输出功率与某地天文、地理、气象指标的关联性，并能根据光伏历史输出功率各信号分量的关联性大小，把各信号分量划分为时序分量或随机分量；
所述关联性为各信号分量与步骤2中时序特征量的相关系数；
步骤5：结合步骤3中天气类型，构建时序训练模型，对步骤4中的时序分量进行训练，形成滑动窗口，得到时序层；
所述时序训练模型为经过引力搜索算法GSA优化的长短期记忆人工神经网络LSTM模型；
步骤6：结合步骤3中天气类型，构建随机训练模型，对步骤4中的随机分量进行训练,得到随机层；
所述随机训练模型为经过贝叶斯优化的机器学习模型；
步骤7：叠加步骤5中时序层和步骤6中随机层，得到光伏出力预测模型；
所述光伏出力预测模型能够根据某地某时的天文数据、地理数据、气象数据，得到某时间段的光伏输出功率。
本发明经过不断探索以及试验，对获取的数据信息进行分类和筛选，并将其划分为：时序特征量、随机特征量；并将光伏历史输出功率数据分成时序分量和随机分量，并根据时序特征量形成时序层、随机层；然后叠加时序层和随机层，得到光伏出力预测模型，对光伏出力进行准确预测。
本发明能够综合考虑随机因素和时序因素对模型的影响，有效提高模型预测的准确度，实现光伏出力的短期准确预测，有助于电网管理人员合理调度和安排电能消纳，以增强电力系统整体稳定性。本发明一种最佳实施例：
一种基于多算法优化的双层混合光伏预测方法，包括以下步骤：
步骤1：获取当地天文、地理和气象数据；
步骤2：对获取的数据进行分类和筛选；
步骤3：利用修正大气清晰度指数k’T和水平面直射比Bd，运用聚类分析算法，将气象数据中的历史天气样本划分为天气类型1，天气类型2，天气类型3，并将天气类型2细分为天气类型2-1和天气类型2-2；
步骤4：利用集合经验模态分解(EEMD)并依据各信号分量(imf分量)与各天文、地理、气象指标的皮尔逊相关性将历史功率数据分成时序分量和随机分量；
步骤5：利用经过引力搜索算法(GSA)优化的长短期记忆人工神经网络(LSTM)对各个天气类型的时序分量进行训练，形成滑动窗口，得到时序层；
步骤6：利用经过贝叶斯优化的机器学习算法对各个天气类型的随机分量进行训练,得到随机层；
步骤7：叠加时序层和随机层，得到光伏出力预测结果。
本发明获取数据特征量的一种具体实施例：
获取的特征量共有20个，分别为年，月，日，时，本站气压，环境温度，相对湿度，降水，平均风速，水平能见度，总云量，赤纬角，时角，修正大气清晰度指数，地外总辐射，水平面总辐射，散射辐射，法向辐射，反射辐射和日照百分率。
本发明对获取的数据进行处理的一种具体实施例：
对获取的数据进行处理，包括对的获取当地天文、地理和气象数据进行分类，筛选。其中将年、月、日、时、赤纬角、时角、地外总辐射I0定义为时序特征量，其他特征量定义为随机特征量，筛除异常值并筛除太阳高度角小于10°的数据。
本发明天气类型划分的一种具体实施例：
天气类型1为恶劣天气，天气类型2为转折天气，天气类型3为晴空天气，天气类型2-1为晴转阴，天气类型2-2为阴转晴。所述步骤3中的聚类分析算法指二阶聚类和逻辑回归。
本发明功率数据划分的一种具体实施例：
利用EEMD将功率数据共分为12个imf分量，经皮尔逊相关性分析后，选取与每个时序特征量相关系数最大的imf分量作为历史功率数据的时序分量，其余特征量为随机分量。
本发明时序分量进行训练的一种具体实施例：
GSA算法选取的超参数变量为训练迭代次数和训练长度。适应度函数选取为预测值与实际值的平均绝对百分比误差(MAPE),计算公式为：
其中N为测试集样本数量，Pf为预测值，Po为实际值。
经过GSA优化的LSTM算法每轮迭代次数为14次，训练集长度选择为6天，所述时序层为滑动窗口LSTM的结果。
本发明随机分量进行训练的一种具体实施例：
经过贝叶斯优化的机器学习算法迭代次数为50次，所述随机层为机器学习算法的预测结果，机器学习算法包括装袋树(Bagging)和高斯过程回归(Gaussian ProcessRegression,GPR)。
本发明修正大气清晰度指数的计算实施例：
修正大气清晰度指数k’T，其计算公式为：
式中，kT为修正前的大气清晰度指数，m为大气质量。
清晰度指数表示大气的透明程度，与天气状况及太阳辐射密切相关，其公式为：
大气层外水平面上太阳辐射量I0：
其中，ISC、γ、δ分别为太阳常数、日地距离订正系数和赤纬角，φ、ω分别为纬度和时角，其计算公式分别如下：
ISC＝1367±7W/m2 (17)
式中，n为所求日期在一年中的日序数。
然而清晰度指数不仅与气象条件相关，而且与天空中太阳位置有关。为了降低太阳高度角对清晰度指数的影响，对其进行修正如下：
其中，k’T是修正后的清晰度指数，m是大气质量。
如图1所示，应用本发明的一种实施例：
1)获取当地天文、地理和气象数据；
2)对获取的数据进行分类和筛选；
3)利用修正大气清晰度指数k’T和水平面直射比Bd，运用聚类分析算法，将历史天气样本划分为天气类型1，天气类型2，天气类型3，并将天气类型2细分为天气类型2-1和天气类型2-2，
如图2和表1所示，在不同特征量下，不同天气类型表现出不同的特点。总体而言，天气类型1为恶劣天气，天气类型2为转折天气，天气类型3为晴空天气，天气类型2-1为晴转阴，天气类型2-2为阴转晴。
表1:不同天气的类型的特征量平均值
4)利用集合经验模态分解(EEMD)并依据各信号分量(imf分量)与各天文、地理、气象指标的皮尔逊相关性将历史功率数据分成时序分量和随机分量，如图2所示，选取的时序分量为imf2、imf7、imf10、imf11、imf12，选取的随机分量为imf1、imf3、imf4、imf5、imf6、imf8、imf9。
5)利用经过GSA优化的长短期记忆人工神经网络(LSTM)对各个天气类型的时序分量进行训练，形成滑动窗口，得到时序层，可参见图3。
滑动窗口是建立时序模型中的一种数据滚动方式，其工作原理为：每次都选择下个时刻之前的66组数据，即时间跨度为6天的数据进行训练。这种训练方式比起传统的整体训练方式，训练时间短，迭代效果好。在针对时序模型的训练中，优势更为明显。
6)利用经过贝叶斯优化的机器学习算法对各个天气类型的随机分量进行训练,得到随机层，具体使用的算法如表2所示。
表2：不同天气类型使用的机器学习算法
7)叠加时序层和随机层，得到光伏出力预测结果。
所述叠加指最终所得时序层结果和随机层结果直接相加。
8)将经过优化的双层混合预测方法所得预测结果与其他方法所得结果相比，如图5以及表3所示，易见其准确度有显著优势。
表3：不同预测方法的误差指标
本发明方法中的天气分型，特征量选择相关参数计算过程及方式如下：
kT为修正前的大气清晰度指数，m为大气质量。
所述修正前的大气清晰度指数，参见图4，其计算公式为：
式中，I表示获取的气象数据中的水平面总辐射，γ为日地距离订正系数，ESC为太阳常数，取值为1367±7W/m2，δ为赤纬角，φ和ω分别为纬度和时角。
所述赤纬角，其计算公式为：
式中，n为所求日期该天在一年中的日序数。
所述日地距离订正系数，其计算公式为：
式中，n为所求日期该天在一年中的日序数。
修正大气清晰度指数k’T，其计算公式为：
评价预测准确度的指标有归一化均方根误差(NRMSE)和平均绝对百分比误差(MAPE),计算公式为：
本发明的计算机设备实施例：
一种计算机设备，其包括：
一个或多个处理器；
存储装置，用于存储一个或多个程序；
当所述一个或多个程序被所述一个或多个处理器执行时，使得所述一个或多个处理器实现上述的一种基于多算法优化的双层混合光伏预测方法。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
最后应当说明的是：以上实施例仅用以说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解：依然可以对本发明的具体实施方式进行修改或者等同替换，而未脱离本发明精神和范围的任何修改或者等同替换，其均应涵盖在本发明的权利要求保护范围之内。
Claims (10)
1.一种基于多算法优化的双层混合光伏预测方法，其特征在于，包括以下步骤：
步骤1：获取某地天文数据、地理数据、气象数据以及光伏历史输出功率数据；
步骤2：对步骤1中获取的数据信息进行分类和筛选，并将其划分为：时序特征量、随机特征量；
步骤3：对步骤1的气象数据，利用修正大气清晰度指数和水平面直射比/>，运用聚类分析算法，将其划分为天气类型1，天气类型2，天气类型3；
并将天气类型2细分为天气类型2-1和天气类型2-2；
步骤4：利用集合经验模态分解EEMD方法，并根据皮尔逊相关性模型将步骤1中的光伏历史输出功率数据分成时序分量和随机分量；
所述皮尔逊相关性模型，用于构建某地光伏输出功率与某地天文、地理、气象指标的关联性，并能根据光伏历史输出功率各信号分量的关联性大小，把各信号分量划分为时序分量或随机分量；
所述关联性为各信号分量与步骤2中时序特征量的相关系数；
步骤5：结合步骤3中天气类型，构建时序训练模型，对步骤4中的时序分量进行训练，形成滑动窗口，得到时序层；
所述时序训练模型为经过引力搜索算法GSA优化的长短期记忆人工神经网络LSTM模型；
步骤6：结合步骤3中天气类型，构建随机训练模型，对步骤4中的随机分量进行训练,得到随机层；
所述随机训练模型为经过贝叶斯优化的机器学习模型；
步骤7：叠加步骤5中时序层和步骤6中随机层，得到光伏出力预测模型；
所述光伏出力预测模型能够根据某地某时的天文数据、地理数据、气象数据，得到某时间段的光伏输出功率。
2.根据权利要求1所述的一种基于多算法优化的双层混合光伏预测方法，其特征在于，
所述步骤1中，获取的数据特征量包括年，月，日，时，本站气压，环境温度，相对湿度，降水，平均风速，水平能见度，总云量，赤纬角，时角，修正大气清晰度指数，地外总辐射，水平面总辐射，散射辐射，法向辐射，反射辐射和日照百分率。
3.根据权利要求2所述的一种基于多算法优化的双层混合光伏预测方法，其特征在于，
所述步骤2中对获取的数据进行处理，包括对获取的某地天文、地理和气象数据进行分类，筛选；
其中时序特征量包括年、月、日、时、赤纬角、时角、地外总辐射I 0；
所述随机特征量包括本站气压，环境温度，相对湿度，降水，平均风速，水平能见度，总云量，修正大气清晰度指数，水平面总辐射，散射辐射，法向辐射，反射辐射和日照百分率；
同时筛除异常值，其包括筛除太阳高度角小于10°的数据。
4.根据权利要求1所述的一种基于多算法优化的双层混合光伏预测方法，其特征在于，
所述步骤3中的天气类型1为恶劣天气，天气类型2为转折天气，天气类型3为晴空天气，天气类型2-1为晴转阴，天气类型2-2为阴转晴。
5.根据权利要求1所述的一种基于多算法优化的双层混合光伏预测方法，其特征在于，
所述步骤3中聚类分析算法为二阶聚类和逻辑回归；
所述修正大气清晰度指数，其计算公式为：
(1)
式中，k T为修正前的大气清晰度指数，m为大气质量；
所述修正前的大气清晰度指数，其计算公式为：
(2)
(3)
式中，I表示获取的气象数据中的水平面总辐射，γ为日地距离订正系数，ESC为太阳常数，取值为1367±7W/m2，δ为赤纬角，φ和ω分别为纬度和时角；
所述赤纬角，其计算公式为：
(4)
式中，n为所求日期在一年中的日序数；
所述日地距离订正系数，其计算公式为：
(5)
式中，n为所求日期在一年中的日序数。
6.根据权利要求1所述的一种基于多算法优化的双层混合光伏预测方法，其特征在于，
所述步骤4中，信号分量的个数为12个；
选取与步骤2中时序特征量相关系数最大的信号分量作为光伏历史输出功率数据的时序分量，其余11个信号分量作为光伏历史输出功率数据的随机分量。
7.根据权利要求1所述的一种基于多算法优化的双层混合光伏预测方法，其特征在于，
所述步骤5中的引力搜索算法GSA算法选取的超参数变量为训练迭代次数和训练长度；其适应度函数为预测值与实际值的平均绝对百分比误差MAPE。
8.根据权利要求7所述的一种基于多算法优化的双层混合光伏预测方法，其特征在于，
所述训练迭代次数为14次，训练长度为6天；
所述平均绝对百分比误差MAPE的计算公式为：
(6)
其中为测试集样本数量，P f为预测值，P o为实际值。
9.根据权利要求1所述的一种基于多算法优化的双层混合光伏预测方法，其特征在于，
所述步骤6中的经过贝叶斯优化的机器学习算法迭代次数为50次，所述随机层为机器学习算法的预测结果，机器学习算法包括装袋Bagging和高斯过程回归GPR。
10.一种计算机设备，其特征在于，包括一个或多个处理器；
存储器，用于存储一个或多个程序；
当所述一个或多个程序被所述一个或多个处理器执行时，使得所述一个或多个处理器实现如权利要求1-9任一所述的一种基于多算法优化的双层混合光伏预测方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110913414.3A CN113722375B (zh) | 2021-08-10 | 2021-08-10 | 一种基于多算法优化的双层混合光伏预测方法及计算机设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110913414.3A CN113722375B (zh) | 2021-08-10 | 2021-08-10 | 一种基于多算法优化的双层混合光伏预测方法及计算机设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113722375A CN113722375A (zh) | 2021-11-30 |
| CN113722375B true CN113722375B (zh) | 2023-10-17 |
Family
ID=78675370
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110913414.3A Active CN113722375B (zh) | 2021-08-10 | 2021-08-10 | 一种基于多算法优化的双层混合光伏预测方法及计算机设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113722375B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115062856A (zh) * | 2022-06-28 | 2022-09-16 | 广东电网有限责任公司 | 基于电网私有云的新能源厂站中期出力预测方法和系统 |
| CN117175569B (zh) * | 2023-09-06 | 2024-03-22 | 国网上海市电力公司 | 一种基于精细化天气分型的光伏预测方法和系统 |
| CN116881743B (zh) * | 2023-09-06 | 2024-01-19 | 国网上海市电力公司 | 一种用于光伏功率预测的气象转折天气识别方法和系统 |
| CN117148875B (zh) * | 2023-10-30 | 2024-01-16 | 余姚市宏宇输变电工程有限公司 | 基于储能配合的光伏板转角控制方法、系统及存储介质 |
| CN117272851B (zh) * | 2023-11-23 | 2024-02-02 | 太原理工大学 | 一种盐碱地光伏板下植物接收光量的建模预测方法 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108923738A (zh) * | 2018-08-09 | 2018-11-30 | 上海电力学院 | 一种基于天气类型判别的双跟踪光伏发电系统控制方法 |
| CN109948281A (zh) * | 2019-03-29 | 2019-06-28 | 上海电力学院 | 基于天气类型有效识别和组合预测的直散分离建模方法 |
| CN110601204A (zh) * | 2019-10-14 | 2019-12-20 | 国网辽宁省电力有限公司盘锦供电公司 | 基于随机变量状态时序模拟的光伏并网系统概率潮流分析方法 |
| CN111242355A (zh) * | 2020-01-06 | 2020-06-05 | 中国电力科学研究院有限公司 | 一种基于贝叶斯神经网络的光伏概率预测方法及系统 |
| CN111585310A (zh) * | 2020-04-07 | 2020-08-25 | 国网冀北电力有限公司 | 一种分布式电源出力预测的方法和装置 |
| CN112348271A (zh) * | 2020-11-12 | 2021-02-09 | 华北电力大学 | 基于vmd-ipso-gru的短期光伏功率预测方法 |
| CN112507793A (zh) * | 2020-11-05 | 2021-03-16 | 上海电力大学 | 一种超短期光伏功率预测方法 |
| CN112668806A (zh) * | 2021-01-17 | 2021-04-16 | 中国南方电网有限责任公司 | 一种基于改进随机森林的光伏功率超短期预测方法 |
-
2021
- 2021-08-10 CN CN202110913414.3A patent/CN113722375B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108923738A (zh) * | 2018-08-09 | 2018-11-30 | 上海电力学院 | 一种基于天气类型判别的双跟踪光伏发电系统控制方法 |
| CN109948281A (zh) * | 2019-03-29 | 2019-06-28 | 上海电力学院 | 基于天气类型有效识别和组合预测的直散分离建模方法 |
| CN110601204A (zh) * | 2019-10-14 | 2019-12-20 | 国网辽宁省电力有限公司盘锦供电公司 | 基于随机变量状态时序模拟的光伏并网系统概率潮流分析方法 |
| CN111242355A (zh) * | 2020-01-06 | 2020-06-05 | 中国电力科学研究院有限公司 | 一种基于贝叶斯神经网络的光伏概率预测方法及系统 |
| CN111585310A (zh) * | 2020-04-07 | 2020-08-25 | 国网冀北电力有限公司 | 一种分布式电源出力预测的方法和装置 |
| CN112507793A (zh) * | 2020-11-05 | 2021-03-16 | 上海电力大学 | 一种超短期光伏功率预测方法 |
| CN112348271A (zh) * | 2020-11-12 | 2021-02-09 | 华北电力大学 | 基于vmd-ipso-gru的短期光伏功率预测方法 |
| CN112668806A (zh) * | 2021-01-17 | 2021-04-16 | 中国南方电网有限责任公司 | 一种基于改进随机森林的光伏功率超短期预测方法 |
Non-Patent Citations (5)
| Title |
|---|
| An LSTM Short-Term Solar Irradiance Forecasting Under Complicated Weather Conditions;Yunjun Yu 等;《IEEE Access》;第7卷;145651 - 145666 * |
| 一种新型天气分型方法及其在光伏功率预测中的应用;李芬 等;《上海交通大学学报》;第55卷(第10期);1510-1519 * |
| 基于CEEMDAN-GSA-LSTM和SVR的光伏功率短期区间预测;李芬 等;《上海交通大学学报》;1-22 * |
| 基于渗透率的区域配电网分布式光伏并网消纳能力分析基于渗透率的区域配电网分布式光伏并网消纳能力分析;赵波 等;《电力系统自动化》;105-111 * |
| 提高光伏消纳能力的配电网联络开关优化配置;王晓港;;电力系统及其自动化学报(第09期);134-140 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113722375A (zh) | 2021-11-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113722375B (zh) | 一种基于多算法优化的双层混合光伏预测方法及计算机设备 | |
| Delle Monache et al. | Kalman filter and analog schemes to postprocess numerical weather predictions | |
| Carneiro et al. | Review on photovoltaic power and solar resource forecasting: current status and trends | |
| CN111079980B (zh) | 基于自适应分类策略和混合优化svr的光功率预测方法 | |
| CN110110912B (zh) | 一种光伏功率多模型区间预测方法 | |
| Fan et al. | A new approach to forecasting typhoon frequency over the western North Pacific | |
| CN110909919A (zh) | 融合注意力机制的深度神经网络模型的光伏功率预测方法 | |
| CN104463349A (zh) | 一种基于多时段综合相似日的光伏发电功率预测方法 | |
| CN111815027A (zh) | 一种光伏场站发电功率预测方法及系统 | |
| CN115936177A (zh) | 一种基于神经网络的光伏输出功率预测方法及系统 | |
| CN117613883A (zh) | 发电功率预测方法、装置、计算机设备及存储介质 | |
| Ogbulezie et al. | A review of regression models employed for predicting diffuse solar radiation in North-Western Africa | |
| CN114399081A (zh) | 一种基于天气分类的光伏发电功率预测方法 | |
| CN107256437B (zh) | 一种光伏电站超短期辐照度预测方法及系统 | |
| CN117374956A (zh) | 一种针对综合能源站光伏发电的短期预测方法 | |
| CN116029419A (zh) | 基于深度学习的长期新能源日均发电功率预测方法及系统 | |
| CN116128168A (zh) | 基于因果扩充卷积和Autoformer的气象预测法 | |
| Sulaiman et al. | Long-term solar irradiance forecasting using multilinear predictors | |
| Hsieh et al. | The influence of large-scale radiation anomalies on tropical cyclone frequency | |
| CN111967652B (zh) | 一种双层协同实时校正光伏预测方法 | |
| Dong et al. | An improved prediction approach on solar irradiance of photovoltaic Power Station | |
| CN115900936A (zh) | 基于天气分型的水平面总辐射和散射辐射预测方法 | |
| CN116484998A (zh) | 基于气象相似日的分布式光伏电站功率预测方法及系统 | |
| CN115238967A (zh) | 一种结合云图和相邻电站集群的光伏功率预测方法及装置 | |
| Sahani et al. | Precise single step and multistep short-term photovoltaic parameters forecasting based on reduced deep convolutional stack autoencoder and minimum variance multikernel random vector functional network |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |