CN105760961A - 光伏输出功率预测方法和系统 - Google Patents
光伏输出功率预测方法和系统 Download PDFInfo
- Publication number
- CN105760961A CN105760961A CN201610131283.2A CN201610131283A CN105760961A CN 105760961 A CN105760961 A CN 105760961A CN 201610131283 A CN201610131283 A CN 201610131283A CN 105760961 A CN105760961 A CN 105760961A
- Authority
- CN
- China
- Prior art keywords
- prediction
- sample
- training
- item
- photovoltaic
- 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
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000012549 training Methods 0.000 claims abstract description 118
- 230000006870 function Effects 0.000 claims description 27
- 238000004458 analytical method Methods 0.000 claims description 19
- 238000012706 support-vector machine Methods 0.000 claims description 15
- 238000000354 decomposition reaction Methods 0.000 claims description 11
- 238000010276 construction Methods 0.000 claims description 9
- 238000012360 testing method Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 abstract description 17
- 230000008859 change Effects 0.000 abstract description 11
- 238000013528 artificial neural network Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 230000008901 benefit Effects 0.000 description 5
- 230000007547 defect Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 238000013519 translation Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 210000002569 neuron Anatomy 0.000 description 3
- 238000010248 power generation Methods 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000005286 illumination Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 241000282414 Homo sapiens Species 0.000 description 1
- 101001095088 Homo sapiens Melanoma antigen preferentially expressed in tumors Proteins 0.000 description 1
- 102100037020 Melanoma antigen preferentially expressed in tumors Human genes 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 238000002790 cross-validation Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000002803 fossil fuel Substances 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- 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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- 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)
- Computer Hardware Design (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- Development Economics (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Game Theory and Decision Science (AREA)
- Evolutionary Computation (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种光伏输出功率预测方法和系统,其中,方法包括:对光伏功率数据按天气类型进行分类,并将预设范围内的光伏功率数据分解为若干个随机项分量和趋势项分量,对随机项分量和趋势项分量,分别构造预测模型并进行训练,根据训练后的预测模型分别获取待预测的光伏功率的随机项分量的第一预测值和趋势项分量的第二预测值,将第一预测值和第二预测值进行叠加求和,获得待预测的光伏输出功率的预测功率,考虑了随机项的非平稳性影响预测效果,提高了对突变天气下光伏输出功率预测的精度。
Description
技术领域
本发明涉及新能源发电和智能电网技术领域,特别是涉及一种光伏输出功率预测方法和系统。
背景技术
随着社会的发展,化石燃料面临枯竭,环境的问题也日益突出。为了解决这一难题,人类开始关注可再生能源,其中绿色无污染的太阳能就成为大家关注的焦点。作为一种有效利用太阳能的方式,光伏发电已经取得了显著的发展。但是考虑到光伏发电受到诸多气候因素的影响,功率变化的随机性较强,一旦功率发生较大变化,必须投入大量备用设备进行削峰填谷,同时输出功率的扰动严重时可能影响电网的安全稳定运行,因此对光伏输出功率的预测研究就显得尤为必要。
考虑实际光伏输出功率的非平稳性,直接对其进行预测的误差较大。现有的光伏输出功率预测方法有人工神经网络(ArtificialNeuralNetworks,ANN)、支持向量机(SupportVectorMachines,SVM)。但是在实际应用过程中,上述方法对突变天气下的光伏输出功率的预测效果不佳。一方面,上述方法数据分析的复杂度较高,另一方面,上述方法对突变天气下的光伏输出功率预测的精度较低。综上所述,现有技术的功率预测效果较差。
发明内容
基于此,有必要针对功率预测效果较差的问题,提供一种光伏输出功率预测方法和系统。
一种光伏输出功率预测方法,包括以下步骤:
根据获取光伏功率数据时的天气特征将光伏功率数据划分为不同天气类型对应的光伏功率数据,将预设范围内的光伏功率数据分解为若干个随机项分量和趋势项分量;
根据随机项分量构造第一训练样本和第一预测样本,根据所述第一训练样本和第一预测样本建立随机项分量的第一预测模型,并根据所述第一训练样本对所述第一预测模型进行训练;
根据趋势项分量构造第二训练样本和第二预测样本,根据所述第二训练样本和第二预测样本建立趋势项分量的第二预测模型,并根据所述第二训练样本对所述第二预测模型进行训练;
将预测样本的输入导入训练后的第一预测模型和训练后的第二预测模型,分别获取待预测的光伏功率的随机项分量的第一预测值和趋势项分量的第二预测值,将第一预测值和第二预测值进行叠加求和,获得待预测的光伏输出功率的预测功率。
一种光伏输出功率预测系统,包括:
分解模块,用于根据获取光伏功率数据时的天气特征将光伏功率数据划分为不同天气类型对应的光伏功率数据,将预设范围内的光伏功率数据分解为若干个随机项分量和趋势项分量;
第一训练模块,用于根据随机项分量构造第一训练样本和第一预测样本,根据所述第一训练样本和第一预测样本建立随机项分量的第一预测模型,并根据所述第一训练样本对所述第一预测模型进行训练;
第二训练模块,用于根据趋势项分量构造第二训练样本和第二预测样本,根据所述第二训练样本和第二预测样本建立趋势项分量的第二预测模型,并根据所述第二训练样本对所述第二预测模型进行训练;
预测模块,用于将预测样本的输入导入训练后的第一预测模型和训练后的第二预测模型,分别获取待预测的光伏功率的随机项分量的第一预测值和趋势项分量的第二预测值,将第一预测值和第二预测值进行叠加求和,获得待预测的光伏输出功率的预测功率。
上述光伏输出功率预测方法和系统,通过对光伏功率数据按天气类型进行分类,并将预设范围内的光伏功率数据分解为若干个随机项分量和趋势项分量,对随机项分量和趋势项分量,分别构造预测模型并进行训练,根据训练后的预测模型分别获取待预测的光伏功率的随机项分量的第一预测值和趋势项分量的第二预测值,将第一预测值和第二预测值进行叠加求和,获得待预测的光伏输出功率的预测功率,考虑了随机项的非平稳性影响预测效果,提高了对突变天气下光伏输出功率预测的精度。
附图说明
图1为本发明的光伏输出功率预测方法流程图;
图2为三种天气类型下的光伏输出功率预测模型框图;
图3为为小波分析效果示意图;
图4为BP神经网络结构示意图;
图5为晴天光伏输出功率短期预测结果;
图6为阴雨天光伏输出功率短期预测结果;
图7为多云天光伏输出功率短期预测结果;
图8为本发明的光伏输出功率预测系统的结构示意图。
具体实施方式
下面结合附图对本发明的光伏输出功率预测方法和系统的实施例进行描述。
图1为本发明的光伏输出功率预测方法流程图。如图1所示,所述光伏输出功率预测方法可包括如下步骤:
S1,根据获取光伏功率数据时的天气特征将光伏功率数据划分为不同天气类型对应的光伏功率数据,将预设范围内的光伏功率数据分解为若干个随机项分量和趋势项分量;
将预设范围内的光伏功率数据分解为若干个随机项分量和趋势项分量的步骤可包括:
对波预设范围内的光伏功率数据进行小波分析,获取离散细节信号和离散逼近信号;
将所述离散细节信号设为随机项分量,将所述离散逼近信号设为趋势项分量。
S2,根据随机项分量构造第一训练样本和第一预测样本,根据所述第一训练样本和第一预测样本建立随机项分量的第一预测模型,并根据所述第一训练样本对所述第一预测模型进行训练;
具体地,根据随机项分量构造第一训练样本和第一预测样本的步骤可包括:
根据预设时间内获取的光伏功率数据构造第一训练样本和第一预测样本的第一输入项;
根据光伏功率预测的影响因子的历史数据构造第一训练样本的第二输入项,根据光伏功率预测的影响因子的预测数据构造第一预测样本的第二输入项;
根据预设时间内获取的光伏功率数据的随机项分量构造第一训练样本和第一预测样本的输出项;
根据所述第一输入项、第一训练样本的第二输入项和输出项构造第一训练样本,根据所述第一输入项、第一预测样本的第二输入项和输出项构造第一预测样本。
根据所述第一训练样本和第一预测样本建立随机项分量的第一预测模型的步骤可包括:
设置若干个弱学习器,初始化所述第一训练样本的样本权重,根据初始化的第一训练样本对各个弱学习器进行训练,获取各个弱学习器的预测函数;
采用预设的测试集对所述预测函数进行测试,得到各个弱学习器的预测误差,根据所述预测误差更新所述样本权重;
重复根据初始化的第一训练样本对各个弱学习器进行训练的步骤,直到达到预设的训练次数,或所述预测误差小于预设的误差期望值;
将各个弱学习器的预测函数进行集成,得到随机项分量的预测函数。
S3,根据趋势项分量构造第二训练样本和第二预测样本,根据所述第二训练样本和第二预测样本建立趋势项分量的第二预测模型,并根据所述第二训练样本对所述第二预测模型进行训练;
根据所述第二训练样本对所述第二预测模型进行训练的步骤可包括:
根据所述第二训练样本设定所述第二预测模型的模型参数;其中,所述模型参数包括SVM的核函数;
建立支持向量机模型,根据所述模型参数确定支持向量机的最优核宽和惩罚因子;
将所述最优核宽和惩罚因子代入SVM的核函数建立所述第二预测模型。
S4,将预测样本的输入导入训练后的第一预测模型和训练后的第二预测模型,分别获取待预测的光伏功率的随机项分量的第一预测值和趋势项分量的第二预测值,将第一预测值和第二预测值进行叠加求和,获得待预测的光伏输出功率的预测功率。
下面根据一个实施例对本发明的技术方案进行进一步阐述。
在一个实施例中,按照天气状况可将光伏功率数据分成晴天、阴雨天和多云天三种天气类型对应的光伏功率数据。首先,考虑到各季节的日出、日落时间不同,为保证数据均有值,可取每天8:00-17:00共10小时的数据进行训练和预测。其次,不同天气类型下的光照强度存在较大差异,光伏输出功率也存在明显差别。因此为了更准确的对光伏功率进行预测,可将光伏功率数据按照天气类型分为晴天、阴雨天、多云天这三种情况,并分别建模。三种天气类型下的光伏输出功率预测模型框图如图2所示。
再者,为避免突变数据或坏数据对预测结果的影响,可对光伏输出功率数据进行预处理,可使用数据纵向对比法,其计算公式如下:
若满足 则:
其中,L(i,t)为第i天第t时刻的光伏功率数据;L(i-1,t)为第i-1天第t时刻的光伏功率数据,即L(i,t)前一天同时刻的光伏功率数据;L(i+1,t)为第i+1天第t时刻的光伏功率数据,即L(i,t)后一天同时刻的光伏功率数据;α(t)和β(t)为预先设定的两个阈值。
对已完成坏数据剔除的原始光伏输出功率数据采用小波分析分解,小波分析分解的具体步骤可包括:
①设x(t)是平方可积信号,即x(t)∈L2(R),则x(t)的连续小波变换定义为:
式中,称为变换后的基本小波,a为尺度参数或伸缩因子,即拉伸压缩基本小波,b为位移参数或平移因子,即平移基本小波,t为时间,ψ(t)为基本小波。
②将a和b同时离散化,即令得到x(t)的离散小波变换(discretewavelettransform,DWT):
为第j次伸缩变换的初始伸缩因子,b0为伸缩后进行平移变换的初始平移因子。该变换既节约计算,又克服了连续小波变换后的信息冗余。
③将信号x(t)正交投影到空间Vj和Wj,能对应得到分辨力j下x(t)的离散逼近信号cj(t)和离散细节信号dj(t)。令j由零逐级增大,便可最终可将信号分解为:
④由此可将波动性较大的光伏功率数据分解为一系列的随机项分量和趋势项分量,具体的分解效果如图3所示。
可选择与预测日同天气类型前5天每时刻光伏功率构造样本输入,对各随机项分量和趋势项分量采用相同方法进行处理,分析并选取光伏功率预测的影响因素(例如,温度、光照强度等)的历史数据和预测数据作为补充输入变量,以所对应预测日的随机项分量和趋势项分量数据作为输出,构造出训练样本和预测样本并归一化;样本输入影响因素的构成如表1所示。
表1样本输入影响因素构成
由此输入向量可表示为:
X(i,t)=[L(i-1,t),L(i-2,t),L(i-3,t),L(i-4,t),L(i-5,t),T(i,t),T(i-1,t),T(i-2,t),T(i-3,t),T(i-4,t),T(i-5,t),S(i,t),S(i-1,t),S(i-2,t),S(i-3,t),S(i-4,t),S(i-5,t)],i>5,j=1,…,M,输出向量为y(i,t)=L(i,t);可对上述输入向量进行归一化,归一化公式可为:
为归一化的输入向量,xmin和xmax分别为x(i)的最小值和最大值,x(i)为输入向量的第i个元素。可通过算法优势互补建立组合预测模型。其中,趋势项可采用SVM算法,随机项可采用BP算法进行预测处理,并可将集成学习引入随机项的预测模型。具体步骤如下:
设定模型参数,其中小波分析采用db4小波基,分解尺度为3;SVM采用RBF核函数;集成学习中弱学习器个数为3,每个学习器的最大训练次数为10,对晴天、阴雨天以及多云天的训练终止条件ξ分别设定为10%、15%、25%,BP神经网络的隐含层为10,其他参数采用默认值;
经过5折交叉验证确定的网格搜索确定支持向量机的最佳参数σ和C。
RBF核函数计算公式具体为:
K(x,xi)=exp(-||x-xi||2/σ2),其中K(x,xi)表示SVM模型的核函数;x表示归一化后的输入向量;xi表示支持向量;σ表示高斯核的核宽,C表示惩罚因子。
对于趋势项分量采用SVM模型进行预测的步骤如下:
假设含有l个训练样本的训练集样本对{(xi,yi),i=1,2,…,l},其中,xi(xi∈Rd)是第i个训练样本的输入列向量,yi∈R为对应的输出值。
设在高维特征空间中建立的线性回归函数为
式中,为原空间映射到高维特征空间的非线性映射函数,ω为权重,b为阈值。
ε-线性不敏感损失函数可以表示为:
式中,f(x)为回归函数返回的预测值,y为对应的真实值。即表示f(x)与y之间的差别小于等于预设值ε,则损失等于0。
引入松弛变量ξ、ξ*,并将上述寻找ω和b的问题用数学语言描述,即
式中,C为惩罚因子,ε规定了回归函数的误差要求,l是训练样本的个数。
根据Mercer定理,定义如下核函数:
是将原空间映射到高维特征空间的非线性映射函数。将拉格朗日乘子αi、引入,并化为对偶形式,得到最优解,最终得到回归函数为:
式中,αi、不为0,对应的样本点即为支持向量。
BP神经网络即采用误差反向传播算法的多层前向神经网络。其重要特点是信号前向传递,误差反向传播。前向传递中,输入信号从输入层经隐含层逐层处理,直至输出层。相邻两层神经元完全互连,不相邻层无连接。每一层的神经元状态只影响下一层神经元状态,如果输出层得不到期望输出,则转入反向传播,根据预测误差调整网络的权值和阀值,使BP神经网络预测输出不断逼近预测输出。其网络结构如图4所示。
考虑到随机项的非平稳性影响预测效果以及BP算法的固有缺点,本发明对随机项的预测模型引入集成学习。具体步骤如下:
给定训练样本S={(x1,y1),...,(xm,ym)},其中,xi∈x,yi∈y,i=1,...,m,弱学习器的个数为T,每个弱学习器的最大训练次数为K。x对应训练样本的输入,y对应相应的输出。
初始化样本权重;
D1(i)=1/m,i=1,...,m(11)
利用样本权重Dt构造弱学习器;得到第t个弱学习器的预测函数表示为ht,其中ht:X→Y;
对弱学习器进行训练,评估弱学习器的训练误差,训练误差可表示为:
εt=∑(|yi-ht(i)|Dt(i))(12)
按照预训练误差εt重新分配各弱学习器的权重:
更新样本权重;
式中,ξ为误差期望,Zt为归一化因子,以使
若尚未达到最大迭代次数K或未满足弱学习器训练误差小于误差期望这一终止条件,则按照公式(13)调整训练样本的样本权重,返回继续执行,否则训练结束,退出循环迭代。
将每个弱学习器的预测函数进行集成,输出最终的预测函数为:
为评价预测效果的好坏,采用相对误差(relativeerror,RE)和平均绝对值百分误差(meanabsolutepercenterror,MAPE)这2个误差指标对预测结果进行误差分析。其值越小,效果越好。指标的表达式为:
式中,ytrue表示真实值,ypre表示预测值,N表示数据总数。
为了说明本发明的预测效果,分别对三种天气类型下的光伏输出功率进行短期预测,其预测效果如图5-图7所示,预测结果及指标分析如表2-表4所示。
表2晴天预测误差统计
表3阴雨天预测误差统计
表4多云天预测误差统计
由图5-图7所示的预测效果和表2-表4所示的指标分析可以明显发现,本发明方法预测结果的趋势更贴近实际光伏功率,预测误差明显较小。综上所述,本发明可较高的实现光伏输出功率短期预测,具有较高的工程应用价值。
针对突变天气下光伏输出功率的预测误差较大的问题,本发明提出了一种小波分析和集成学习的光伏输出功率短期预测方法。首先按照天气状况将光伏输出功率数据分成晴天、阴雨天和多云天三种类型,保证了样本数据的一致性;并采用小波分析将以上三种光伏输出功率数据分解为一系列的随机项分量和趋势项分量,降低了数据的波动性;然后,通过算法优势互补建立组合预测模型。其中,趋势项采用SVM算法,随机项采用BP算法进行预测处理;并考虑到随机项的非平稳性和BP算法的固有缺点,将集成学习引入随机项的预测模型,进一步提高了模型的预测效果;最后将各分量的预测结果进行叠加得到最终预测结果。本发明方法从数据预处理和模型优化两个方面进行处理,有效的提高了模型的预测精度,可较好的运用于工程实际问题。
与已有技术相比,本发明的有益效果体现在:
1、本发明引入小波分解降低数据的波动性,原理简单,结果准确,有效提高了预测精度;
2、本发明一方面针对单一算法预测突变天气的光伏输出功率时误差较大的问题,通过算法优势互补建立组合预测模型;另一方面考虑到随机项的非平稳性影响预测效果以及BP算法的固有缺点,对随机项的预测模型引入集成学习,通过对具有一定预测精度的弱学习器组合形成具有很高预测精度的强学习器,进一步提升了预测模型的性能,提高了对突变天气下光伏输出功率预测的精度。
图8为本发明的光伏输出功率预测系统的结构示意图。如图8所示,所述光伏输出功率预测系统可包括:
分解模块10,用于根据获取光伏功率数据时的天气特征将光伏功率数据划分为不同天气类型对应的光伏功率数据,将预设范围内的光伏功率数据分解为若干个随机项分量和趋势项分量;
第一训练模块20,用于根据随机项分量构造第一训练样本和第一预测样本,根据所述第一训练样本和第一预测样本建立随机项分量的第一预测模型,并根据所述第一训练样本对所述第一预测模型进行训练;
第二训练模块30,用于根据趋势项分量构造第二训练样本和第二预测样本,根据所述第二训练样本和第二预测样本建立趋势项分量的第二预测模型,并根据所述第二训练样本对所述第二预测模型进行训练;
预测模块40,用于将预测样本的输入导入训练后的第一预测模型和训练后的第二预测模型,分别获取待预测的光伏功率的随机项分量的第一预测值和趋势项分量的第二预测值,将第一预测值和第二预测值进行叠加求和,获得待预测的光伏输出功率的预测功率。
其中,所述分解模块可包括:
小波分析单元,用于对波预设范围内的光伏功率数据进行小波分析,获取离散细节信号和离散逼近信号;
设置单元,用于将所述离散细节信号设为随机项分量,将所述离散逼近信号设为趋势项分量。
所述第一训练模块可包括:
第一构造单元,用于根据预设时间内获取的光伏功率数据构造第一训练样本和第一预测样本的第一输入项;
第二构造单元,用于根据光伏功率预测的影响因子的历史数据构造第一训练样本的第二输入项,根据光伏功率预测的影响因子的预测数据构造第一预测样本的第二输入项;
第三构造单元,用于根据预设时间内获取的光伏功率数据的随机项分量构造第一训练样本和第一预测样本的输出项;
第四构造单元,用于根据所述第一输入项、第一训练样本的第二输入项和输出项构造第一训练样本,根据所述第一输入项、第一预测样本的第二输入项和输出项构造第一预测样本。
针对突变天气下光伏输出功率的预测误差较大的问题,本发明提出了一种小波分析和集成学习的光伏输出功率短期预测方法。首先按照天气状况将光伏输出功率数据分成晴天、阴雨天和多云天三种类型,保证了样本数据的一致性;并采用小波分析将以上三种光伏输出功率数据分解为一系列的随机项分量和趋势项分量,降低了数据的波动性;然后,通过算法优势互补建立组合预测模型。其中,趋势项采用SVM算法,随机项采用BP算法进行预测处理;并考虑到随机项的非平稳性和BP算法的固有缺点,将集成学习引入随机项的预测模型,进一步提高了模型的预测效果;最后将各分量的预测结果进行叠加得到最终预测结果。本发明方法从数据预处理和模型优化两个方面进行处理,有效的提高了模型的预测精度,可较好的运用于工程实际问题。
与已有技术相比,本发明的有益效果体现在:
1、本发明引入小波分解降低数据的波动性,原理简单,结果准确,有效提高了预测精度;
2、本发明一方面针对单一算法预测突变天气的光伏输出功率时误差较大的问题,通过算法优势互补建立组合预测模型;另一方面考虑到随机项的非平稳性影响预测效果以及BP算法的固有缺点,对随机项的预测模型引入集成学习,通过对具有一定预测精度的弱学习器组合形成具有很高预测精度的强学习器,进一步提升了预测模型的性能,提高了对突变天气下光伏输出功率预测的精度。
本发明的光伏输出功率预测系统与本发明的光伏输出功率预测方法一一对应,在上述光伏输出功率预测方法的实施例阐述的技术特征及其有益效果均适用于光伏输出功率预测系统的实施例中,特此声明。
以上所述实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。
Claims (10)
1.一种光伏输出功率预测方法,其特征在于,包括以下步骤:
根据获取光伏功率数据时的天气特征将光伏功率数据划分为不同天气类型对应的光伏功率数据,将预设范围内的光伏功率数据分解为若干个随机项分量和趋势项分量;
根据随机项分量构造第一训练样本和第一预测样本,根据所述第一训练样本和第一预测样本建立随机项分量的第一预测模型,并根据所述第一训练样本对所述第一预测模型进行训练;
根据趋势项分量构造第二训练样本和第二预测样本,根据所述第二训练样本和第二预测样本建立趋势项分量的第二预测模型,并根据所述第二训练样本对所述第二预测模型进行训练;
将预测样本的输入导入训练后的第一预测模型和训练后的第二预测模型,分别获取待预测的光伏功率的随机项分量的第一预测值和趋势项分量的第二预测值,将第一预测值和第二预测值进行叠加求和,获得待预测的光伏输出功率的预测功率。
2.根据权利要求1所述的光伏输出功率预测方法,其特征在于,将预设范围内的光伏功率数据分解为若干个随机项分量和趋势项分量的步骤包括:
对波预设范围内的光伏功率数据进行小波分析,获取离散细节信号和离散逼近信号;
将所述离散细节信号设为随机项分量,将所述离散逼近信号设为趋势项分量。
3.根据权利要求1所述的光伏输出功率预测方法,其特征在于,根据随机项分量构造第一训练样本和第一预测样本的步骤包括:
根据预设时间内获取的光伏功率数据构造第一训练样本和第一预测样本的第一输入项;
根据光伏功率预测的影响因子的历史数据构造第一训练样本的第二输入项,根据光伏功率预测的影响因子的预测数据构造第一预测样本的第二输入项;
根据预设时间内获取的光伏功率数据的随机项分量构造第一训练样本和第一预测样本的输出项;
根据所述第一输入项、第一训练样本的第二输入项和输出项构造第一训练样本,根据所述第一输入项、第一预测样本的第二输入项和输出项构造第一预测样本。
4.根据权利要求1所述的光伏输出功率预测方法,其特征在于,根据所述第一训练样本和第一预测样本建立随机项分量的第一预测模型的步骤包括:
设置若干个弱学习器,初始化所述第一训练样本的样本权重,根据初始化的第一训练样本对各个弱学习器进行训练,获取各个弱学习器的预测函数;
采用预设的测试集对所述预测函数进行测试,得到各个弱学习器的预测误差,根据所述预测误差更新所述样本权重;
重复根据初始化的第一训练样本对各个弱学习器进行训练的步骤,直到达到预设的训练次数,或所述预测误差小于预设的误差期望值;
将各个弱学习器的预测函数进行集成,得到随机项分量的预测函数。
5.根据权利要求4所述的光伏输出功率预测方法,其特征在于,根据所述预测误差更新所述样本权重的步骤包括:
根据如下公式更新所述样本权重:
其中,
式中,Dt(i)为更新前的样本权重,Dt+1(i)为更新后的样本权重,ξ为误差期望指,Zt为归一化因子,ht(i)为第t个弱学习器的预测函数,t为弱学习器的个数,εt为预测误差,yi为第一训练样本的输出项,i为当前训练次数。
6.根据权利要求1所述的光伏输出功率预测方法,其特征在于,根据所述第二训练样本对所述第二预测模型进行训练的步骤包括:
根据所述第二训练样本设定所述第二预测模型的模型参数;其中,所述模型参数包括SVM的核函数;
建立支持向量机模型,根据所述模型参数确定支持向量机的最优核宽和惩罚因子;
将所述最优核宽和惩罚因子代入SVM的核函数建立所述第二预测模型。
7.根据权利要求1所述的光伏输出功率预测方法,其特征在于,还包括以下步骤:
计算所述第一预测值的第一相对误差和第一平均绝对值百分误差,计算所述第二预测值的第二相对误差和第二平均绝对值百分误差;
根据所述第一相对误差和第一平均绝对值百分误差对所述第一预测值进行误差分析,根据所述第二相对误差和第二平均绝对值百分误差对所述第二预测值进行误差分析。
8.一种光伏输出功率预测系统,其特征在于,包括:
分解模块,用于根据获取光伏功率数据时的天气特征将光伏功率数据划分为不同天气类型对应的光伏功率数据,将预设范围内的光伏功率数据分解为若干个随机项分量和趋势项分量;
第一训练模块,用于根据随机项分量构造第一训练样本和第一预测样本,根据所述第一训练样本和第一预测样本建立随机项分量的第一预测模型,并根据所述第一训练样本对所述第一预测模型进行训练;
第二训练模块,用于根据趋势项分量构造第二训练样本和第二预测样本,根据所述第二训练样本和第二预测样本建立趋势项分量的第二预测模型,并根据所述第二训练样本对所述第二预测模型进行训练;
预测模块,用于将预测样本的输入导入训练后的第一预测模型和训练后的第二预测模型,分别获取待预测的光伏功率的随机项分量的第一预测值和趋势项分量的第二预测值,将第一预测值和第二预测值进行叠加求和,获得待预测的光伏输出功率的预测功率。
9.根据权利要求8所述的光伏输出功率预测系统,其特征在于,所述分解模块包括:
小波分析单元,用于对波预设范围内的光伏功率数据进行小波分析,获取离散细节信号和离散逼近信号;
设置单元,用于将所述离散细节信号设为随机项分量,将所述离散逼近信号设为趋势项分量。
10.根据权利要求8所述的光伏输出功率预测系统,其特征在于,所述第一训练模块包括:
第一构造单元,用于根据预设时间内获取的光伏功率数据构造第一训练样本和第一预测样本的第一输入项;
第二构造单元,用于根据光伏功率预测的影响因子的历史数据构造第一训练样本的第二输入项,根据光伏功率预测的影响因子的预测数据构造第一预测样本的第二输入项;
第三构造单元,用于根据预设时间内获取的光伏功率数据的随机项分量构造第一训练样本和第一预测样本的输出项;
第四构造单元,用于根据所述第一输入项、第一训练样本的第二输入项和输出项构造第一训练样本,根据所述第一输入项、第一预测样本的第二输入项和输出项构造第一预测样本。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610131283.2A CN105760961A (zh) | 2016-03-08 | 2016-03-08 | 光伏输出功率预测方法和系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610131283.2A CN105760961A (zh) | 2016-03-08 | 2016-03-08 | 光伏输出功率预测方法和系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105760961A true CN105760961A (zh) | 2016-07-13 |
Family
ID=56332791
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610131283.2A Pending CN105760961A (zh) | 2016-03-08 | 2016-03-08 | 光伏输出功率预测方法和系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105760961A (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106803135A (zh) * | 2017-01-17 | 2017-06-06 | 广东工业大学 | 一种光伏发电系统输出功率的预测方法及装置 |
CN106815659A (zh) * | 2017-01-20 | 2017-06-09 | 国网浙江省电力公司电力科学研究院 | 一种基于混合模型的超短期太阳辐射预测方法及其装置 |
CN107133692A (zh) * | 2017-04-16 | 2017-09-05 | 信阳师范学院 | 基于加号模型的光伏发电功率预测方法 |
CN107358323A (zh) * | 2017-07-13 | 2017-11-17 | 上海交通大学 | 一种短期光伏发电的功率预测方法 |
CN107545327A (zh) * | 2017-08-23 | 2018-01-05 | 国网辽宁省电力有限公司经济技术研究院 | 基于支持向量机的光伏发电短期输出功率综合预测方法 |
CN108256693A (zh) * | 2018-02-11 | 2018-07-06 | 阳光电源股份有限公司 | 一种光伏发电功率预测方法、装置及系统 |
CN110598896A (zh) * | 2019-07-26 | 2019-12-20 | 陕西省水利电力勘测设计研究院 | 基于预测误差校正的光伏功率预测方法 |
CN110738363A (zh) * | 2019-09-30 | 2020-01-31 | 华中科技大学 | 一种光伏发电功率预测模型及其构建方法和应用 |
CN111168569A (zh) * | 2020-01-08 | 2020-05-19 | 华中科技大学 | 一种磨削材料去除量预测方法、装置、设备及存储介质 |
CN113159102A (zh) * | 2021-02-24 | 2021-07-23 | 山东大学 | 一种多时间尺度的光伏功率预测方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105023070A (zh) * | 2015-08-12 | 2015-11-04 | 河海大学常州校区 | 一种光伏系统输出功率预测方法 |
CN105117975A (zh) * | 2015-08-06 | 2015-12-02 | 国家电网公司 | 一种光伏输出功率值的分频预测方法 |
-
2016
- 2016-03-08 CN CN201610131283.2A patent/CN105760961A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105117975A (zh) * | 2015-08-06 | 2015-12-02 | 国家电网公司 | 一种光伏输出功率值的分频预测方法 |
CN105023070A (zh) * | 2015-08-12 | 2015-11-04 | 河海大学常州校区 | 一种光伏系统输出功率预测方法 |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106803135A (zh) * | 2017-01-17 | 2017-06-06 | 广东工业大学 | 一种光伏发电系统输出功率的预测方法及装置 |
CN106815659B (zh) * | 2017-01-20 | 2020-09-18 | 国网浙江省电力公司电力科学研究院 | 一种基于混合模型的超短期太阳辐射预测方法及其装置 |
CN106815659A (zh) * | 2017-01-20 | 2017-06-09 | 国网浙江省电力公司电力科学研究院 | 一种基于混合模型的超短期太阳辐射预测方法及其装置 |
CN107133692A (zh) * | 2017-04-16 | 2017-09-05 | 信阳师范学院 | 基于加号模型的光伏发电功率预测方法 |
CN107358323A (zh) * | 2017-07-13 | 2017-11-17 | 上海交通大学 | 一种短期光伏发电的功率预测方法 |
CN107545327A (zh) * | 2017-08-23 | 2018-01-05 | 国网辽宁省电力有限公司经济技术研究院 | 基于支持向量机的光伏发电短期输出功率综合预测方法 |
CN108256693A (zh) * | 2018-02-11 | 2018-07-06 | 阳光电源股份有限公司 | 一种光伏发电功率预测方法、装置及系统 |
CN108256693B (zh) * | 2018-02-11 | 2024-02-13 | 阳光电源股份有限公司 | 一种光伏发电功率预测方法、装置及系统 |
CN110598896A (zh) * | 2019-07-26 | 2019-12-20 | 陕西省水利电力勘测设计研究院 | 基于预测误差校正的光伏功率预测方法 |
CN110738363A (zh) * | 2019-09-30 | 2020-01-31 | 华中科技大学 | 一种光伏发电功率预测模型及其构建方法和应用 |
CN110738363B (zh) * | 2019-09-30 | 2022-04-01 | 华中科技大学 | 一种光伏发电功率预测方法 |
CN111168569A (zh) * | 2020-01-08 | 2020-05-19 | 华中科技大学 | 一种磨削材料去除量预测方法、装置、设备及存储介质 |
CN113159102A (zh) * | 2021-02-24 | 2021-07-23 | 山东大学 | 一种多时间尺度的光伏功率预测方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105760961A (zh) | 光伏输出功率预测方法和系统 | |
Li et al. | Smart wind speed forecasting approach using various boosting algorithms, big multi-step forecasting strategy | |
CN111860982A (zh) | 一种基于vmd-fcm-gru的风电场短期风电功率预测方法 | |
CN109146162B (zh) | 一种基于集成循环神经网络的概率风速预测方法 | |
CN112100911B (zh) | 一种基于深度bilstm的太阳辐射预测方法 | |
CN109242212A (zh) | 一种基于变模式分解与长短记忆网络的风电预测方法 | |
CN104463349A (zh) | 一种基于多时段综合相似日的光伏发电功率预测方法 | |
CN111027775A (zh) | 基于长短期记忆网络的梯级水电站发电量预测方法 | |
CN110443417A (zh) | 基于小波变换的多模型集成负荷预测方法 | |
CN112434891A (zh) | 基于wcnn-alstm的太阳辐照度时间序列的预测方法 | |
CN111553510A (zh) | 一种短期风速预测方法 | |
CN111144644A (zh) | 基于变分异方差高斯过程回归的短期风速预测方法 | |
CN113111592A (zh) | 一种基于emd-lstm的短期风电功率预测方法 | |
CN113935513A (zh) | 一种基于ceemdan的短期电力负荷预测方法 | |
CN116345555A (zh) | 一种基于cnn-isca-lstm模型的短期光伏发电功率预测方法 | |
Yang et al. | A novel method of wind speed prediction by peephole LSTM | |
CN113837434A (zh) | 太阳能光伏发电预测方法、装置、电子设备和存储介质 | |
CN116632820A (zh) | 基于序列延拓vmd与风速波动特性的风电功率预测方法 | |
CN118095891A (zh) | 考虑源荷气象特性解耦的有源配电网净负荷预测方法及系统 | |
CN118040678A (zh) | 一种短期海上风电功率组合预测方法 | |
CN118035728A (zh) | 一种基于双重信号分解和双网络通道的短期光伏功率集成预测方法 | |
CN118134284A (zh) | 一种基于多阶段注意力机制的深度学习风电功率预测方法 | |
CN117113054A (zh) | 一种基于图神经网络和Transformer的多元时间序列预测方法 | |
CN116613732A (zh) | 一种基于shap值选择策略的多元负荷预测方法及系统 | |
CN113779861A (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: 20160713 |