CN116167508A - 一种基于气象因素分解的短期光伏出力快速预测方法及系统 - Google Patents
一种基于气象因素分解的短期光伏出力快速预测方法及系统 Download PDFInfo
- Publication number
- CN116167508A CN116167508A CN202310057695.6A CN202310057695A CN116167508A CN 116167508 A CN116167508 A CN 116167508A CN 202310057695 A CN202310057695 A CN 202310057695A CN 116167508 A CN116167508 A CN 116167508A
- Authority
- CN
- China
- Prior art keywords
- fish
- value
- photovoltaic output
- model
- meteorological
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 56
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 37
- 241000251468 Actinopterygii Species 0.000 claims abstract description 89
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 38
- 238000005457 optimization Methods 0.000 claims abstract description 36
- 238000012706 support-vector machine Methods 0.000 claims abstract description 35
- 238000004364 calculation method Methods 0.000 claims abstract description 32
- 230000002431 foraging effect Effects 0.000 claims abstract description 28
- 238000012549 training Methods 0.000 claims abstract description 24
- 230000006870 function Effects 0.000 claims description 41
- 235000013305 food Nutrition 0.000 claims description 39
- 230000008569 process Effects 0.000 claims description 17
- 238000010606 normalization Methods 0.000 claims description 13
- 238000005295 random walk Methods 0.000 claims description 12
- 238000004590 computer program Methods 0.000 claims description 7
- 239000011159 matrix material Substances 0.000 claims description 6
- 241001465754 Metazoa Species 0.000 claims description 3
- 230000006399 behavior Effects 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 230000009467 reduction Effects 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 238000012804 iterative process Methods 0.000 claims 1
- 239000000126 substance Substances 0.000 claims 1
- 238000010248 power generation Methods 0.000 description 9
- 238000003860 storage Methods 0.000 description 6
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
Images
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
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- 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
- 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)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- General Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Biophysics (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- Operations Research (AREA)
- Primary Health Care (AREA)
- Life Sciences & Earth Sciences (AREA)
- Quality & Reliability (AREA)
- Biomedical Technology (AREA)
- Entrepreneurship & Innovation (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于气象因素分解的短期光伏出力快速预测方法及系统,预测方法包括:(1)提出一种气象因素相关性分解系数计算方法,实现对光伏出力气象因素分量的快速分解计算;(2)建立改进支持向量机的预测模型对气象因素数值进行预测,采用一种一致性测度函数来选择模型的训练样本对模型进行训练,并基于改进鱼群觅食优化算法优化模型参数;(3)最后,基于优化完成的模型,快速预测下一时刻的气象因素值,计算其对应的光伏出力分量,并利用气象因素相关性分解系数计算总光伏出力。本发明实现对光伏出力气象因素分量的快速分解计算;同样也能基于预测的个光伏出力气象因素分量反求总光伏出力,实现光伏短期出力的快速预测,提高光伏电站的实际运行效率。
Description
技术领域
本发明属于光伏出力预测技术领域,尤其涉及一种基于气象因素分解的短期光伏出力快速预测方法及系统。
背景技术
随着世界经济以及全球化的快速发展,能源也越来越不能完全满足现代人们的日常生产以及生活使用需求。因此如何发展可再生清洁能源,对于整个世界各国的经济发展意义重大。太阳能发电作为可再生清洁能源之一,其性能清洁环保无污染,生产成本低,而且光伏发电系统的正常运行以及维护的主要工作量较小。然而,光伏发电装机容量的迅速增大也给电网的安全稳定运行带来了一些问题。光伏发电是一种光生伏特效应,光的质量将影响光伏出力,而光的质量会受温度、相对湿度、风速、云层厚度等气象因素的影响特别大。可能会给其他电网设备造成冲击、带来电压波动,不利于电力系统的稳定正常运行。当天气变化时,光伏出力会出现间歇、波动较大的情况。当光伏装机容量达到一定规模并入电网后,会对整个电网的稳定运行造成较大的冲击,加大了电网调度部门协调安排发电的难度。
因此,当大规模光伏并网时,同时保证电网对光伏发电功率的消纳,需要对光伏发电功率进行预测,尤其是针对不确定性的天气的光伏出力预测。对光伏出力预测有利于提高电力系统的安全稳定运行和电网调度部门的调度决策,包括安排发电计划、调整发电功率的峰谷差、电力系统的潮流优化和安排设备检修等。对大型光伏电站在不同天气下发电功率的精准预测,能为电网调度部门提供数据依据,有利于电网调度部门统筹安排常规发电厂与光伏电站发电的协调配合,对电力系统安全稳定运行有积极的影响。
发明内容
本发明针对现有技术中的不足,提供一种基于气象因素分解的短期光伏出力快速预测方法及系统。本发明基于改进支持向量机模型对气象因素进行快速预测,计算其对应的光伏出力分量,并利用气象因素相关性分解系数计算总光伏出力总量,实现光伏短期出力的快速预测,提高光伏电站的实际运行效率。
为实现上述目的,本发明采用以下技术方案:
一种基于气象因素分解的短期光伏出力快速预测方法,其特征在于,包括:
提出一种将气象因素相关性分解系数计算方法,实现对光伏出力气象因素分量的快速分解计算;
建立改进支持向量机的预测模型对气象因素数值进行预测,采用一种一致性测度函数来选择模型的训练样本对模型进行训练,并基于改进鱼群觅食优化算法优化模型参数,
最后,基于优化完成的模型,快速预测下一时刻的气象因素值,计算其对应的光伏出力分量,并利用气象因素相关性分解系数计算总光伏出力总量。
提出一种分解光伏出力的气象因素相关性系数计算方法,可以根据此系数将光伏出力进行分解,计算方法见式(1)
式中,Ppv表示光伏出力总量,Ppv Z、Ppv W和Ppv R分别表示光伏出力基于太阳辐照度、环境温度和相对湿度这三种气象因素的气象分量,λZ、λW和λR分别表示气象因素太阳辐照度、环境温度和相对湿度对光伏出力的相关性权重因子,其计算方法见式(2),
式中,kZ、kW和kR分别表示气象因素太阳辐照度、环境温度和相对湿度对光伏出力的相关性系数,其计算方法见式(3),
式中,Q表示计算所用的样本数量,Zq *表示第q个样本中气象因素Z归一化处理后的数值,表示第q个样本中光伏的输出功率经过式(4)归一化处理后的数值,YZ *和YP *分别表示对应样气象因素Z和本的光伏的输出Ppv的样本平均值。
式中,Ppv *表示经过式归一化处理后的某时刻光伏输出功率数值,Ppv为某时刻光伏输出功率实际值,Ppv min和Ppv max分别为光伏电站出力最小值和最大值。
气象因素太阳辐照度、环境温度和相对湿度其归一化计算方式见式(5)~(7),
式中,Z*、W*和R*分别为太阳辐照度、环境温度和相对湿度这三种气象因素经过归一化处理后的数值,Z、W和R分别为太阳辐照度、环境温度和相对湿度这三种气象因素的实际数值,Zmin和Zmax分别为太阳辐照度的最小值和最大值。
同理可求出λW和λR。
采用一种改进支持向量机相的预测模型,该模型运用一种一致性测度函数,从大量的历史气象-光伏出力数据中提取与预测样本一致性较高的数据创建训练样本,训练支持向量机预测模型,最后预测下一时刻的气象因素数值,
具体预测过程如下
步骤1,基于一致性理论,采用一种一致性测度函数,选取一致性较高的训练样本;
步骤2,利用一种改进鱼群觅食优化算法,求得气象因素对应的支持向量机相模型中的最优惩罚系数和核函数参数
步骤3,构造样本对(Xt,Y t),Xt={xt-m,xt-m+1,…,xt-1}为模型输入,Yt=xt即为模型的输出,m为决定预测值的历史数据个数,即模型输入向量的维数;
步骤4,若是单步预测则将最新实际的xt加入历史序列末端,若是进行多步预测则将预测出来的xt加入历史序列末端,确定输入为Xt+1,Xt+1={xt-m+1,xt-m+2,…,xt}对xt+1预测,利用上面的方法,最终得到xt+1的值。
步骤5,每预测完一步,用最新的实际数据或者预测值更新历史数据,再对下一时刻的值进行预测。
再根据式(1),可以快速计算出t+1时刻光伏出力总值Ppv t+1。
其中所述改进支持向量机模型具体建立方式见下,
对于非线性回归问题,给定训练样本集{xn,yn},n=1、2、…、N,其中输入量xn∈Rn,对应的输出量yn∈R,N为样本数。通过非线性函数将xn,映射到高维特征空间,再进行回归运算,支持向量回归表达式见式(9):
式中,αn,α′n,为大于零的拉格朗日乘子,对应的向量为支持向量,其值根据式(10)来确定。xn为训练样本中的输入分量,x为测试样本中的输入分量,xn、x∈Rn。b为偏置值,是一个常数。
式中,c为惩罚系数;β1和β2是松弛变量,其值大于零;ε为允许误差系数,yn为样本输入xn对应实际值,f(xn)为模型输出值。
ψ(xn,x)为核函数,其作用是给定低维空间的xn和x,直接计算出在高维空间的内积值,其具体函数表达式见式(11)。
式中,σ为核函数参数。
在运用模型对气象因素进行预测时,需要确定模型的2个参数值,即惩罚系数C和核函数参数σ。为提高模型的预测精度采用一种改进鱼群觅食优化算法,对这2个参数进行优化选择。
所述一致性测度函数计算方法见式(12)
式中,A={a1,a2,…,aD}为待测样本序列,B={b1,b2,…,bD}为待选样本序列,
上式函数的值域为[0 1],τ(A,B)越接近1,则待测样本序列A与待与选样本序列B的一致性越高。
若满足式(13),则选取序列B为改进支持向量机模型的训练样本。
τ(A,B)≥γ (13)
式中,γ为一致性系数基准值,其取值范围为[0 1]。
采用一种改进鱼群觅食优化算法,对改进的支持向量机模型中的2个参数值,即惩罚系数c和核函数参数σ进行寻优计算。
改进鱼群觅食优化是一种新的群体智能优化算法,其模拟动物界中鱼群觅食过程和行为,主要包括5个步骤,具体如下。
步骤1:参数设置,设置鱼群优化算法的种群数量L,最大迭代次数S,变量个数H;变量的上边界OU=[oU 1,oU 2,…],变量的下边界OD=[oD 1,oD 2,…];
步骤2:种群位置的初始化,种群包括鱼群和食物2种生物,食物位置U和鱼儿个体位置V由式(14)和式(15)
式中:U和V分别为食物和鱼儿个体的位置矩阵;uln和vln分别为第l个食物和鱼儿个体对应的第h个参数值;l=1,2,…,L,h=1,2,…,H。
鱼个体和食物初始位置U0n和v0n根据式(16)计算
式中:r1和r2为随机数;oU h和oD h分别为第h个变量的上、下边界。
步骤3:精英鱼的选择,精英鱼是指每次迭代过程所有鱼群中最优的鱼个体,可由适应度值E进行判断,可以根据式(17)计算。
式中:f(x)和y分别为关于输入x的改进支持向量机模型预测值关于输入x的和实际值。
每次迭代过程中鱼儿的适应度矩阵E为:
因此,初代精英鱼即为算法首次迭代过程中具有最优适应度值的鱼个体。
步骤4:优化迭代,在算法中,每个食物的位置移动受精英个体和鱼个体的影响。第s代食物的位置ul s为:
式中:R1 s为依据精英鱼个体位置的第s代食物个体的随机游走;R2 s为依据鱼群位置的第s代食物的随机游走。
为更有利于鱼群觅食,算法中会根据式(20)逐渐缩小随机游走的范围,
式中:J随机游走范围缩小比例,可由式(21)计算得出。
式中:s#为当前迭代次数位置系数,计算方法见式(22)
步骤5:最优鱼儿(精英)的选择,鱼儿在捕猎食物之后,为了更有利于新的觅食,鱼儿将调整位置。在优化算法中,适应度值大于鱼儿的这类食物将被捕获。在捕获到食物后,鱼儿将根据捕获到的食物的位置更新其位置,即:
vl s=ul s,f(ul s)>f(vl s) (23)
在每次迭代过程中,鱼儿都会根据式(23)更新其位置和对应的适应度,进而可确定在该次迭代中最优的鱼儿个体。当达到算法终止的条件时,即可获得最优的鱼儿个体(即精英),此时其值就是改进支持向量机模型的最优惩罚系数c和核函数参数σ。
一种基于气象因素分解的短期光伏出力快速预测系统,其特征在于,包括:
光伏出力气象分量分解模块,用于确定分解光伏出力的气象因素相关性系数;
基于改进支持向量机模型的预测模块,用于快速预测下一时刻的气象因素值,并根据气象因素预测值计算光伏总出力预测值;
模型训练样本选择模块,用于选择与待测数据组一致性较高的历史数据作为训练样本数据;
模型优化模块,用于基于一种改进鱼群觅食优化算法,对改进的支持向量机模型中的惩罚系数c和核函数参数σ进行寻优计算。
一种电子设备,其特征在于,包括:存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行计算机程序时,实现如权利要求1-5任一项所述的基于气象因素分解的短期光伏出力快速预测方法。
本发明的有益效果是:
1、针对光伏出力受气象因素影响较大的问题,本发明一种气象因素相关性分解系数计算方法,实现对光伏出力气象因素分量的快速分解计算;同样也能基于预测的个光伏出力气象因素分量反求总光伏出力;
2、建立一种改进支持向量机的预测模型对气象因素数值进行预测,采用一种一致性测度函数来选择模型的训练样本对模型进行训练,并基于改进鱼群觅食优化算法优化模型参数,实现光伏短期出力的快速预测,提高光伏电站的实际运行效率。
附图说明
图1是本发明提出的一种基于气象因素分解的短期光伏出力快速预测方法的总体流程图。
图2本发明所建立改进支持向量机模型结构图。
图3本发明所采用改进鱼群觅食优化算法的流程图。
具体实施方式
下面结合附图和实施例对本发明作进一步详细的说明。
如图1所示一种基于气象因素分解的短期光伏出力快速预测方法及系统,首先提出一种气象因素相关性分解系数计算方法,实现对光伏出力气象因素分量的快速分解计算;其次建立一种如图2所示改进支持向量机的预测模型对气象因素数值进行预测,采用一种一致性测度函数来选择模型的训练样本对模型进行训练,并采用一种如图3所示基于改进鱼群觅食优化算法优化模型参数;最后,基于优化完成的模型,快速预测下一时刻的气象因素值,计算其对应的光伏出力分量,并利用气象因素相关性分解系数计算总光伏出力,实现光伏短期出力的快速预测,提高光伏电站的实际运行效率。
本发明所提出的一种面向交直流混合电网的电力系统无功优化方法,具体实
施过程如下:
一、提出一种分解光伏出力的气象因素相关性系数计算方法,可以根据此系数将光伏出力进行分解,计算方法见式(1)
Ppv=λZ·Ppv Z+λW·Ppv W+λR·Ppv R (1)
式中,Ppv表示光伏出力总量,Ppv Z、Ppv W和Ppv R分别表示光伏出力基于太阳辐照度、环境温度和相对湿度这三种气象因素的气象分量,λZ、λW和λR分别表示气象因素太阳辐照度、环境温度和相对湿度对光伏出力的相关性权重因子,其计算方法见式(2),
式中,kZ、kW和kR分别表示气象因素太阳辐照度、环境温度和相对湿度对光伏出力的相关性系数,其计算方法见式(3),
式中,Q表示计算所用的样本数量,Zq *表示第q个样本中气象因素Z归一化处理后的数值,Ppv * q表示第q个样本中光伏的输出功率经过式(4)归一化处理后的数值,YZ *和YP *分别表示对应样气象因素Z和本的光伏的输出Ppv的样本平均值。
相关系数k的取值范围为(-1,1),其中k>0为正相关,表明一组变量随另一组变量的增大(减小)一致的发生增大(减小),k<0为负相关,表明一组变量随另一组变量的增大(减小)相反的发生减小(增大),其中的绝对值越接近1则相关程度越大、越接近0则相关程度越小。
式中,Ppv *表示经过式归一化处理后的某时刻光伏输出功率数值,Ppv为某时刻光伏输出功率实际值,Ppv min和Ppv max分别为光伏电站出力最小值和最大值。
气象因素太阳辐照度、环境温度和相对湿度其归一化计算方式见式(5)~(7),
式中,Z*、W*和R*分别为太阳辐照度、环境温度和相对湿度这三种气象因素经过归一化处理后的数值,Z、W和R分别为太阳辐照度、环境温度和相对湿度这三种气象因素的实际数值,Zmin和Zmax分别为太阳辐照度的最小值和最大值。
湿度是用来形容大气干燥程度的物理量,一般自身就是一个相对百分数。
同理可求出λW和λR。
二、采用一种改进支持向量机相的预测模型,该模型运用一种一致性测度函数,从大量的历史气象-光伏出力数据中提取与预测样本一致性较高的数据创建训练样本,训练支持向量机预测模型,最后预测下一时刻的气象因素数值,
具体预测过程如下
步骤1,基于一致性理论,采用一种一致性测度函数,选取一致性较高的训练样本;
步骤2,利用一种改进鱼群觅食优化算法,求得气象因素对应的支持向量机相模型中的最优惩罚系数和核函数参数
步骤3,构造样本对(Xt,Y t),Xt={xt-m,xt-m+1,…,xt-1}为模型输入,Yt=xt即为模型的输出,m为决定预测值的历史数据个数,即模型输入向量的维数;
步骤4,若是单步预测则将最新实际的xt加入历史序列末端,若是进行多步预测则将预测出来的xt加入历史序列末端,确定输入为Xt+1,Xt+1={xt-m+1,xt-m+2,…,xt}对xt+1预测,利用上面的方法,最终得到xt+1的值。
步骤5,每预测完一步,用最新的实际数据或者预测值更新历史数据,再对下一时刻的值进行预测。
步骤6,根据预测得到t+1时刻的气象因素太阳辐照度、环境温度和相对湿度归一化后的值Z* t+1、W* t+1和R* t+1数值,根据式(8)计算t+1时刻光伏出力的各气象因素分量Ppv Z t+1、Ppv W t+1和Ppv R t+1,
式中,Z* t、W* t和R* t分别为气象因素太阳辐照度、环境温度和相对湿度归一化后的值,Ppv Z t、Ppv W t和Ppv R t分别为t时刻光伏出力的太阳辐照度、环境温度和相对湿度分量值。
再根据式(1),可以快速计算出t+1时刻光伏出力总值Ppv t+1。
其中所述改进支持向量机的结构如图2所示,具体建立方式见下,
对于非线性回归问题,给定训练样本集{xn,yn},n=1、2、…、N,其中输入量xn∈Rn,对应的输出量yn∈R,N为样本数。通过非线性函数将xn,映射到高维特征空间,再进行回归运算,支持向量回归表达式见式(9):
式中,αn,α′n,为大于零的拉格朗日乘子,对应的向量为支持向量,其值根据式(10)来确定。xn为训练样本中的输入分量,x为测试样本中的输入分量,xn、x∈Rn。b为偏置值,是一个常数。
式中,c为惩罚系数;β1和β2是松弛变量,其值大于零;ε为允许误差系数,yn为样本输入xn对应实际值,f(xn)为模型输出值。
ψ(xn,x)为核函数,其作用是给定低维空间的xn和x,直接计算出在高维空间的内积值,其具体函数表达式见式(11)。
式中,σ为核函数参数。
在运用模型对气象因素进行预测时,需要确定模型的2个参数值,即惩罚系数C和核函数参数σ。为提高模型的预测精度采用一种改进鱼群觅食优化算法,对这2个参数进行优化选择。
三、所述一致性测度函数计算方法见式(12)
式中,A={a1,a2,…,aD}为待测样本序列,B={b1,b2,…,bD}为待选样本序列,
上式函数的值域为[0 1],τ(A,B)越接近1,则待测样本序列A与待与选样本序列B的一致性越高。
若满足式(13),则选取序列B为改进支持向量机模型的训练样本。
τ(A,B)≥γ (13)
式中,γ为一致性系数基准值,其取值范围为[0 1]。
四、采用一种改进鱼群觅食优化算法,对改进的支持向量机模型中的2个参数值,即惩罚系数c和核函数参数σ进行寻优计算。
改进鱼群觅食优化是一种新的群体智能优化算法,其模拟动物界中鱼群觅食过程和行为,主要包括5个步骤,具体如下。
步骤1:参数设置,设置鱼群优化算法的种群数量L,最大迭代次数S,变量个数H;变量的上边界OU=[oU 1,oU 2,…],变量的下边界OD=[oD 1,oD 2,…];
步骤2:种群位置的初始化,种群包括鱼群和食物2种生物,食物位置U和鱼儿个体位置V由式(14)和式(15)
式中:U和V分别为食物和鱼儿个体的位置矩阵;uln和vln分别为第l个食物和鱼儿个体对应的第h个参数值;l=1,2,…,L,h=1,2,…,H。
鱼个体和食物初始位置U0n和v0n根据式(16)计算
式中:r1和r2为随机数;oU h和oD h分别为第h个变量的上、下边界。
步骤3:精英鱼的选择,精英鱼是指每次迭代过程所有鱼群中最优的鱼个体,可由适应度值E进行判断,可以根据式(17)计算。
式中:f(x)和y分别为关于输入x的改进支持向量机模型预测值关于输入x的和实际值。
每次迭代过程中鱼儿的适应度矩阵E为:
因此,初代精英鱼即为算法首次迭代过程中具有最优适应度值的鱼个体。
步骤4:优化迭代,在算法中,每个食物的位置移动受精英个体和鱼个体的影响。第s代食物的位置ul s为:
式中:R1 s为依据精英鱼个体位置的第s代食物个体的随机游走;R2 s为依据鱼群位置的第s代食物的随机游走。
为更有利于鱼群觅食,算法中会根据式(20)逐渐缩小随机游走的范围,
式中:J随机游走范围缩小比例,可由式(21)计算得出。
式中:s#为当前迭代次数位置系数,计算方法见式(22)
步骤5:最优鱼儿(精英)的选择,鱼儿在捕猎食物之后,为了更有利于新的觅食,鱼儿将调整位置。在优化算法中,适应度值大于鱼儿的这类食物将被捕获。在捕获到食物后,鱼儿将根据捕获到的食物的位置更新其位置,即:
vl s=ul s,f(ul s)>f(vl s) (23)
在每次迭代过程中,鱼儿都会根据式(23)更新其位置和对应的适应度,进而可确定在该次迭代中最优的鱼儿个体。当达到算法终止的条件时,即可获得最优的鱼儿个体(即精英),此时其值就是改进支持向量机模型的最优惩罚系数c和核函数参数σ。
在另一实施例中,本发明提出了一种计算机可读存储介质,存储有计算机程序,计算机程序使计算机执行如上所述的基于气象因素分解的短期光伏出力快速预测方法。
在另一实施例中,本发明提出了一种电子设备,包括:存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,处理器执行计算机程序时,实现如上所述的基于气象因素分解的短期光伏出力快速预测方法。
在本申请所公开的实施例中,计算机存储介质可以是有形的介质,其可以包含或存储以供指令执行系统、装置或设备使用或与指令执行系统、装置或设备结合使用的程序。计算机存储介质可以包括但不限于电子的、磁性的、光学的、电磁的、红外的、或半导体系统、装置或设备,或者上述内容的任何合适组合。计算机存储介质的更具体示例会包括基于一个或多个线的电气连接、便携式计算机盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦除可编程只读存储器(EPROM或快闪存储器)、光纤、便捷式紧凑盘只读存储器(CD-ROM)、光学储存设备、磁储存设备、或上述内容的任何合适组合。
本领域普通技术人员可以意识到,结合本申请所公开的实施例描述的各示例的单元及算法步骤,能够以电子硬件或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本申请的范围。
以上仅是本发明的优选实施方式,本发明的保护范围并不仅局限于上述实施例,凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理前提下的若干改进和润饰,应视为本发明的保护范围。
Claims (7)
1.一种基于气象因素分解的短期光伏出力快速预测方法,其特征是:包括下列步骤:
(1)提出一种气象因素相关性分解系数计算方法,实现对光伏出力气象因素分量的快速分解计算;
(2)建立改进支持向量机的预测模型对气象因素数值进行预测,采用一种一致性测度函数来选择模型的训练样本对模型进行训练,并基于改进鱼群觅食优化算法优化模型参数;
(3)最后,基于优化完成的模型,快速预测下一时刻的气象因素值,计算其对应的光伏出力分量,并利用气象因素相关性分解系数计算总光伏出力。
2.根据权利要求1所述的基于气象因素分解的短期光伏出力快速预测方法,其特征是:
步骤(1)所述的气象因素相关性分解系数计算方法,见式(1)
Ppv=λZ·Ppv Z+λW·Ppv W+λR·Ppv R (1)
式中,Ppv表示光伏出力总量,Ppv Z、Ppv W和Ppv R分别表示光伏出力基于太阳辐照度、环境温度和相对湿度这三种气象因素的气象分量;λZ、λW和λR分别表示气象因素太阳辐照度、环境温度和相对湿度对光伏出力的相关性权重因子,其计算方法见式(2),
式中,kZ、kW和kR分别表示气象因素太阳辐照度、环境温度和相对湿度对光伏出力的相关性系数,其计算方法见式(3),
式中,Q表示计算所用的样本数量,Zq *表示第q个样本中气象因素Z归一化处理后的数值,Ppv * q表示第q个样本中光伏的输出功率经过式(4)归一化处理后的数值,YZ *和YP *分别表示对应样气象因素Z和本的光伏的输出Ppv的样本平均值;
式中,Ppv *表示经过式归一化处理后的某时刻光伏输出功率数值,Ppv为某时刻光伏输出功率实际值,Ppv min和Ppv max分别为光伏电站出力最小值和最大值;
气象因素太阳辐照度、环境温度和相对湿度其归一化计算方式见式(5)~(7),
式中,Z*、W*和R*分别为太阳辐照度、环境温度和相对湿度这三种气象因素经过归一化处理后的数值,Z、W和R分别为太阳辐照度、环境温度和相对湿度这三种气象因素的实际数值,Zmin和Zmax分别为太阳辐照度的最小值和最大值;
同理可求出λW和λR。
3.根据权利要求1所述的基于气象因素分解的短期光伏出力快速预测方法,其特征是:步骤(2)包括:
采用一种改进支持向量机相的预测模型,该模型运用一种一致性测度函数,从大量的历史气象-光伏出力数据中提取与预测样本一致性较高的数据创建训练样本,训练支持向量机预测模型,最后预测下一时刻的气象因素数值;
具体预测过程如下
步骤1,基于一致性理论,采用一种一致性测度函数,选取一致性较高的训练样本;
步骤2,利用一种改进鱼群觅食优化算法,求得气象因素对应的支持向量机相模型中的最优惩罚系数和核函数参数
步骤3,构造样本对(Xt,Yt),Xt={xt-m,xt-m+1,…,xt-1}为模型输入,Yt=xt即为模型的输出,m为决定预测值的历史数据个数,即模型输入向量的维数;
步骤4,若是单步预测则将最新实际的xt加入历史序列末端,若是进行多步预测则将预测出来的xt加入历史序列末端,确定输入为Xt+1,Xt+1={xt-m+1,xt-m+2,…,xt}对xt+1预测,利用上面的方法,最终得到xt+1的值;
步骤5,每预测完一步,用最新的实际数据或者预测值更新历史数据,再对下一时刻的值进行预测;
步骤6,根据预测得到t+1时刻的气象因素太阳辐照度、环境温度和相对湿度归一化后的值Z* t+1、W* t+1和R* t+1数值,根据式(8)计算t+1时刻光伏出力的各气象因素分量Ppv Z t+1、Ppv W t+1和Ppv R t+1,
式中,Z* t、W* t和R* t分别为气象因素太阳辐照度、环境温度和相对湿度归一化后的值,Ppv Z t、Ppv W t和Ppv R t分别为t+1时刻光伏出力的太阳辐照度、环境温度和相对湿度分量值;
再根据式(1),可以快速计算出t+1时刻光伏出力总值Ppv t+1;
其中所述改进支持向量机模型具体建立方式见下:
对于非线性回归问题,给定训练样本集{xn,yn},n=1、2、…、N,其中输入量xn∈Rn,对应的输出量yn∈R,N为样本数;通过非线性函数将xn,映射到高维特征空间,再进行回归运算,支持向量回归表达式见式(9):
式中,αn,α′n,为大于零的拉格朗日乘子,对应的向量为支持向量,其值根据式(10)来确定;xn为训练样本中的输入分量,x为测试样本中的输入分量,xn、x∈Rn;b为偏置值,是一个常数;
式中,c为惩罚系数;β1和β2是松弛变量,其值大于零;ε为允许误差系数,yn为样本输入xn对应实际值,f(xn)为模型输出值;
ψ(xn,x)为核函数,其作用是给定低维空间的xn和x,直接计算出在高维空间的内积值,其具体函数表达式见式(11):
式中,σ为核函数参数。
5.根据权利要求3所述的基于气象因素分解的短期光伏出力快速预测方法,其特征是:
采用一种改进鱼群觅食优化算法,对改进的支持向量机模型中的2个参数值,即惩罚系数c和核函数参数σ进行寻优计算;
改进鱼群觅食优化是一种新的群体智能优化算法,其模拟动物界中鱼群觅食过程和行为,主要包括5个步骤,具体如下:
步骤1:参数设置,设置鱼群优化算法的种群数量L,最大迭代次数S,变量个数H;变量的上边界OU=[oU 1,oU 2,…],变量的下边界OD=[oD 1,oD 2,…];
步骤2:种群位置的初始化,种群包括鱼群和食物2种生物,食物位置U和鱼儿个体位置V由式(14)和式(15)
式中:U和V分别为食物和鱼儿个体的位置矩阵;uln和vln分别为第l个食物和鱼儿个体对应的第h个参数值;l=1,2,…,L,h=1,2,…,H;
鱼个体和食物初始位置U0n和v0n根据式(16)计算
式中:r1和r2为随机数;oU h和oD h分别为第h个变量的上、下边界;
步骤3:精英鱼的选择,精英鱼是指每次迭代过程所有鱼群中最优的鱼个体,可由适应度值E进行判断,可以根据式(17)计算;
式中:f(x)和y分别为关于输入x的改进支持向量机模型预测值关于输入x的和实际值;
每次迭代过程中鱼儿的适应度矩阵E为:
因此,初代精英鱼即为算法首次迭代过程中具有最优适应度值的鱼个体;
步骤4:优化迭代,在算法中,每个食物的位置移动受精英个体和鱼个体的影响;第s代食物的位置ul s为:
式中:R1 s为依据精英鱼个体位置的第s代食物个体的随机游走;R2 s为依据鱼群位置的第s代食物的随机游走;
为更有利于鱼群觅食,算法中会根据式(20)逐渐缩小随机游走的范围,
式中:J随机游走范围缩小比例,可由式(21)计算得出:
式中:s#为当前迭代次数位置系数,计算方法见式(22)
步骤5:最优鱼儿(精英)的选择,鱼儿在捕猎食物之后,为了更有利于新的觅食,鱼儿将调整位置;在优化算法中,适应度值大于鱼儿的这类食物将被捕获;在捕获到食物后,鱼儿将根据捕获到的食物的位置更新其位置,即:
vl s=ul s,f(ul s)>f(vl s) (23)
在每次迭代过程中,鱼儿都会根据式(23)更新其位置和对应的适应度,进而可确定在该次迭代中最优的鱼儿个体;当达到算法终止的条件时,即可获得最优的鱼儿个体(即精英),此时其值就是改进支持向量机模型的最优惩罚系数c和核函数参数σ。
6.一种基于气象因素分解的短期光伏出力快速预测系统,其特征是:包括:
光伏出力气象分量分解模块,用于确定分解光伏出力的气象因素相关性系数;
基于改进支持向量机模型的预测模块,用于快速预测下一时刻的气象因素值,并根据气象因素预测值计算光伏总出力预测值;
模型训练样本选择模块,用于选择与待测数据组一致性较高的历史数据作为训练样本数据;
模型优化模块,用于基于一种改进鱼群觅食优化算法,对改进的支持向量机模型中的惩罚系数c和核函数参数σ进行寻优计算。
7.一种电子设备,其特征在于,包括:存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行计算机程序时,实现如权利要求1-6任一项所述的基于气象因素分解的短期光伏出力快速预测方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310057695.6A CN116167508B (zh) | 2023-01-16 | 2023-01-16 | 一种基于气象因素分解的短期光伏出力快速预测方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310057695.6A CN116167508B (zh) | 2023-01-16 | 2023-01-16 | 一种基于气象因素分解的短期光伏出力快速预测方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116167508A true CN116167508A (zh) | 2023-05-26 |
CN116167508B CN116167508B (zh) | 2024-03-15 |
Family
ID=86412792
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310057695.6A Active CN116167508B (zh) | 2023-01-16 | 2023-01-16 | 一种基于气象因素分解的短期光伏出力快速预测方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116167508B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117154851A (zh) * | 2023-10-24 | 2023-12-01 | 国网浙江省电力有限公司宁波供电公司 | 一种基于空间域的光伏输出控制方法、电子设备及介质 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018045689A1 (zh) * | 2016-09-07 | 2018-03-15 | 江南大学 | 光伏二极管箝位型三电平逆变器的决策树svm故障诊断方法 |
CN110059891A (zh) * | 2019-04-28 | 2019-07-26 | 湖北工业大学 | 一种基于vmd-svm-wsa-gm组合模型的光伏电站输出功率预测方法 |
CN110909958A (zh) * | 2019-12-05 | 2020-03-24 | 国网江苏省电力有限公司南通供电分公司 | 一种计及光伏并网功率的短期负荷预测方法 |
CN111079980A (zh) * | 2019-11-22 | 2020-04-28 | 天合云能源互联网技术(杭州)有限公司 | 基于自适应分类策略和混合优化svr的光功率预测方法 |
CN112348271A (zh) * | 2020-11-12 | 2021-02-09 | 华北电力大学 | 基于vmd-ipso-gru的短期光伏功率预测方法 |
CN112561178A (zh) * | 2020-12-21 | 2021-03-26 | 国网甘肃省电力公司电力科学研究院 | 一种考虑微气象因素的分布式光伏出力预测方法 |
CN114004424A (zh) * | 2021-12-21 | 2022-02-01 | 南方电网科学研究院有限责任公司 | 一种光伏功率预测方法、系统、设备及存储介质 |
CN114819365A (zh) * | 2022-05-05 | 2022-07-29 | 合肥工业大学 | 一种基于lsgan-dcn网络的短期光伏功率预测方法 |
CN115347571A (zh) * | 2022-10-17 | 2022-11-15 | 国网江西省电力有限公司电力科学研究院 | 基于迁移学习的光伏发电功率短期预测方法及装置 |
CN115423143A (zh) * | 2022-07-15 | 2022-12-02 | 国网安徽省电力有限公司经济技术研究院 | 一种基于回归分析和lssvm的源荷一体化预测方法 |
-
2023
- 2023-01-16 CN CN202310057695.6A patent/CN116167508B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018045689A1 (zh) * | 2016-09-07 | 2018-03-15 | 江南大学 | 光伏二极管箝位型三电平逆变器的决策树svm故障诊断方法 |
CN110059891A (zh) * | 2019-04-28 | 2019-07-26 | 湖北工业大学 | 一种基于vmd-svm-wsa-gm组合模型的光伏电站输出功率预测方法 |
CN111079980A (zh) * | 2019-11-22 | 2020-04-28 | 天合云能源互联网技术(杭州)有限公司 | 基于自适应分类策略和混合优化svr的光功率预测方法 |
CN110909958A (zh) * | 2019-12-05 | 2020-03-24 | 国网江苏省电力有限公司南通供电分公司 | 一种计及光伏并网功率的短期负荷预测方法 |
CN112348271A (zh) * | 2020-11-12 | 2021-02-09 | 华北电力大学 | 基于vmd-ipso-gru的短期光伏功率预测方法 |
CN112561178A (zh) * | 2020-12-21 | 2021-03-26 | 国网甘肃省电力公司电力科学研究院 | 一种考虑微气象因素的分布式光伏出力预测方法 |
CN114004424A (zh) * | 2021-12-21 | 2022-02-01 | 南方电网科学研究院有限责任公司 | 一种光伏功率预测方法、系统、设备及存储介质 |
CN114819365A (zh) * | 2022-05-05 | 2022-07-29 | 合肥工业大学 | 一种基于lsgan-dcn网络的短期光伏功率预测方法 |
CN115423143A (zh) * | 2022-07-15 | 2022-12-02 | 国网安徽省电力有限公司经济技术研究院 | 一种基于回归分析和lssvm的源荷一体化预测方法 |
CN115347571A (zh) * | 2022-10-17 | 2022-11-15 | 国网江西省电力有限公司电力科学研究院 | 基于迁移学习的光伏发电功率短期预测方法及装置 |
Non-Patent Citations (2)
Title |
---|
朱梅梅;苏建徽;陈智慧;: "基于EEMD和IPSO的SVM短期光伏出力预测", 电气工程学报, no. 04 * |
葛乐;陆文伟;袁晓冬;周前;: "基于改进相似日和ABC-SVM的光伏电站功率预测", 太阳能学报, no. 03 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117154851A (zh) * | 2023-10-24 | 2023-12-01 | 国网浙江省电力有限公司宁波供电公司 | 一种基于空间域的光伏输出控制方法、电子设备及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN116167508B (zh) | 2024-03-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107766990B (zh) | 一种光伏电站发电功率的预测方法 | |
CN109802430B (zh) | 一种基于LSTM-Attention网络的风电电网控制方法 | |
Mahmoud et al. | An advanced approach for optimal wind power generation prediction intervals by using self-adaptive evolutionary extreme learning machine | |
Liu et al. | An improved photovoltaic power forecasting model with the assistance of aerosol index data | |
Yang et al. | A weather-based hybrid method for 1-day ahead hourly forecasting of PV power output | |
CN112733462B (zh) | 一种结合气象因素的超短期风电场功率预测方法 | |
CN111008728A (zh) | 一种用于分布式光伏发电系统短期出力的预测方法 | |
CN107766986A (zh) | 泄漏积分型回声状态网络在线学习光伏功率预测方法 | |
CN104978613A (zh) | 一种考虑组件温度的光伏出力短期预测方法 | |
CN113554466B (zh) | 一种短期用电量预测模型构建方法、预测方法和装置 | |
CN114970952B (zh) | 一种考虑环境因素的光伏出力短期预测方法及系统 | |
CN110866633B (zh) | 一种基于svr支持向量回归的微电网超短期负荷预测方法 | |
CN114021420A (zh) | 一种分布式光伏超短期发电功率预测方法及系统 | |
CN114462718A (zh) | 基于时间滑动窗口的cnn-gru风电功率预测方法 | |
CN116167508B (zh) | 一种基于气象因素分解的短期光伏出力快速预测方法及系统 | |
CN112149905A (zh) | 一种基于小波变换和小波神经网络的光伏电站短期功率预测方法 | |
CN116070769A (zh) | 一种超短期风电场功率多步区间预测模块化方法及其设备 | |
CN116341728A (zh) | 一种基于数据驱动的超短期光伏输出功率预测方法 | |
Meng et al. | A new PV generation power prediction model based on GA-BP neural network with artificial classification of history day | |
Yang et al. | Kalman filter photovoltaic power prediction model based on forecasting experience | |
Alharbi et al. | Short-term wind speed and temperature forecasting model based on gated recurrent unit neural networks | |
CN112215478B (zh) | 光储电站的功率协调控制方法和装置、存储介质 | |
Zhang et al. | Real-time prediction of solar radiation based on online sequential extreme learning machine | |
Jun-Ma et al. | Photovoltaic power generation prediction based on MEA-BP neural network | |
CN112036672B (zh) | 一种基于迭代订正的新能源发电超短期功率预测方法及系统 |
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 |