一种考虑组件温度的光伏出力短期预测方法
【技术领域】
[0001] 本发明设及一种光伏出力短期预测方法，尤其是设及一种考虑组件温度的光伏出 力短期预测方法。
【背景技术】
[0002] 光伏发电相对于传统能源发电而言，具有一定的间歇性、波动性、周期性等特点， 该些特点会给电网的运行带来较大的影响。从实际来说，对电力调度部口而言，能够较准确 的对光伏发电功率进行预测，可W提早得到光伏电站输出功率的数据及变化规律。可W提 高电网的稳定性、可靠性和经济性。
[0003] 研究表明，晶体娃太阳能电池的最大输出功率与组件温度之间存在负相关性，即 组件温度越低，最大输出功率越高；组件温度越高，最大输出功率越低。也可W说，当组件 温度不断上升时，短路电流不断变大，开路电压逐渐减小，输出功率值也不断减小。
[0004] 而目前对光伏发电系统输出功率预测，用环境温度代替光伏组件实际温度，或者 不考虑组件温度的影响，必然会造成明显的误差。因此本发明考虑组件温度对光伏出力的 影响，可W进一步提高光伏出力的精确度。
【发明内容】
[0005] 本发明的上述技术问题主要是通过下述技术方案得W解决的：
[0006] 一种考虑组件温度的光伏出力短期预测方法，其特征在于，包括W下步骤：
[0007] 步骤1，选取相邻日影响组件温度的天气因素，对预测日组件温度进行预测，具体 包括运用粒子群进行算法参数寻优W及最小二乘支持向量机算法对组件温度进行预测两 个子步骤；
[000引其中，所述运用粒子群进行算法参数寻优的具体方法是：
[0009] 步骤1. 11，粒子初始化，设置相关参数；
[0010] 步骤1. 12,评价和计算每个粒子的适应值；
[0011] 步骤1. 13,计算每个粒子的个体最优位置和群体最优位置，并计算每次迭代时粒 子个体历史最优位置的高斯扰动值；
[0012] 步骤1. 14,按照下式对粒子速度和位置进行更新，并计算新的粒子适应度值；
[001引Vi, j (t+1) = wVi, j (t) +(3江1虹，j (t)甘2gausSi, j (t) -Xi, j (t)]
[0014] +C2r3[Pg，j(t)-Xi,j(t)]
[0015] Xi j(t+1)=Xi j(t)+Vi j(t+l)
[0016] gausSij(t) =r4gaussian(y, 5^)
[0017] 其中，w为惯性权重，Cl、C2为加速因子，r1、f2、r3、表示（；0, 1)区间服从均匀分 布的随机数，Vu (t)为粒子i在第t次迭代时的速度，Pu (t)为粒子i在第t次迭代时的 历史最优位置，gaussu(t)为粒子i在第t迭代时产生的高斯扰动，y为均值，S2为方差， xu(t)为粒子i在第t次迭代时的位置，是第t次迭代时种群的最优位置；
[0018] 步骤1. 15,若算法满足终止条件，则停止迭代，输出全局最优粒子位置及其适应度 值；否则转到步骤1. 12继续循环；
[0019] 所述最小二乘支持向量机算法对组件温度进行预测的具体方法是：
[0020] 步骤1. 21，对收集的气象数据和组件温度数据进行归一化处理，形成训练样本矩 阵；
[002U步骤1. 22,确定LSSVM模型所需的运行参数，设定GDPSO算法的参数；
[0022] 步骤1. 23,计算各个粒子当前的适应值和高斯扰动项的值，更新粒子的速度和位 置；
[0023] 步骤1. 24,计算更新后粒子的新适应值，并完成对粒子个体和群体的最优解替 换；
[0024] 步骤1.25,若达到最大迭代次数，则结束优化，输出优化参数；否则，返回步骤 1. 23继续循环；
[0025] 步骤1. 26,将优化的参数赋给LSSVM，得到优化的预测模型，完成预测；
[0026] 步骤2,运用滴权法选取预测日的相似日；
[0027] 步骤3,分析影响光伏出力在不同时刻在出力中的变异程度，结合统计方法 佑DPSO-LSSVM)对光伏出力进行短期预测。
[002引在上述的一种考虑组件温度的光伏出力短期预测方法，所述的步骤2中，具体的 操作步骤如下：
[0029] 步骤2. 1，W每天7:00-20:00中每60min为一个基准点，提取预测日气象特征数 值，构成预测日指标向量
[0030] X〇=[X0(7)，x0(；8)，…，x0(20)，y0(7)，y0(；8)，…，y0(20)，z0(7)，z0(；8)，…，z0(20)]; 将与预测日接近的14天作为历史日，提取历史日的相对应指标构成历史日指标向量：
[0031] Xj=[xj(7),Xj(8), ???,Xj(20),yj(7),yj(8), ???,y0(20),Zj(7),Zj(8), ???,Zj(20)], 其中，Xj表示组件温度，yj表示湿度，zj表示福照强度，j= 1，2，…，14 ;
[0032] 步骤2. 2,求解Xj.与X。差值绝对值构成差值矩阵，并对差值矩阵进行归一化处理， 得到预测日与历史日的指标数据差值归一化矩阵；
[0033] 步骤2. 3,由差值矩阵确定各个基值点的权重参数，将权重与各历史天的指标差值 归一化矩阵相应点相乘，并求取加权和作为相似误差；由相似误差从小到大排序，选出相似 误差最小的5天作为预测日的相似日。
[0034] 在上述的一种考虑组件温度的光伏出力短期预测方法，所述的步骤3中，具体的 操作步骤如下：
[0035] 步骤3. 1，在相似日中选取与预测时刻对应的所需时刻的组件温度、湿度W及福照 强度等指标数据；
[0036] 步骤3. 2,构成矩阵，利用滴权法程序分析此时刻每个指标所提供的信息权重；
[0037] 步骤3. 3,将得到的权重与其对应的相似日中各指标原始数据进行相乘，得到新的 加权输入矩阵；
[003引步骤3. 4,将加权输入矩阵和输出功率作为输入数据，对此时刻的光伏输出功率进 行预测。
[0039] 因此，本发明具有如下优点；1.本发明提出的考虑组件温度影响因素前提下，利 用滴权法与GDPSD-LSSVM算法进行光伏功率短期预测，具有误差小，预测效果良好的优点。【附图说明】
[0040] 图1是本发明中运用运用粒子群寻优的最小二乘支持向量机算法对组件温度进 行预测的流程图。
【具体实施方式】
[0041] 下面通过实施例，并结合附图，对本发明的技术方案作进一步具体的说明。
[0042] 实施例；
[0043] 下面是采用本发明的具体实施例的步骤，如下：
[0044] 一、组件温度预测方法
[0045] 运用粒子群寻优的最小二乘支持向量机算法对组件温度进行预测W备后续的输 出功率预测提供数据支持。
[0046] 1)通过机器学习，针对光伏发电系统，根据所提供的训练样本，寻找系统输入和 输出之间的未知依赖关系。本发明采用最小二乘支持向量机（LSSVM)算法进行机器学习。 LSSVM优化指标采用的平方项，运用等式约束代替标准支持向量机的不等式约束，将标准 支持向量机中的求解二次规划问题转化为求解线性问题，可简化计算复杂度，提高计算速 度。LSSVM采用结构风险最小化准则，使用核函数解决非线性回归问题，假设（x^yi)，i= 1，2，"'N是给定的数据集合，且满足条件XGr，yGR，对于非线性的样本输入向量Xi，可 W通过预先定义的映射函数抑.r)将其映射到高维特征空间，在新的高维特征空间中，可将 其当成线性函数求解，因为原来的非线性函数经过映射函数的变换己经成为高维特征空 间中的线性函数，此时决策函数可W表示为：
[0047] /(.V)二U??刮.V) + A
[0048] 其中，W是权重向量，b是偏移量。采用结构风险最小化原则寻求权重向量和偏移 量的值，LSSVM的优化问题可W表示为
[0化1] 其中，丫为正则化参数，曰1表示误差。从上式可W看出，最小二乘支持向量机直接 将损失函数定义为误差的平方和，将传统的标准支持向量机中的不等式约束转换成了等式 约束。
[0052]引入拉格朗日乘子a1，可从隐下式转化为；
[0053]
[0054] 根据KKT条件可得：
[0 化 5]
[0056] 消去W和e河W得到线性方程：
[0057]
[0058] 其中，a*= [a。02,…ajT, 1*= [1，1，…U了为狀1 维列向Y= [y^y^wyj T，QGR胃，并且a二例-V,)?如r,) =A'(和y,)。非线性回归的预测模型的表达式可W写为；
[0059]
[0060] 其中，为拉格朗日乘子，K(Xi，Xj.)为核函数，满足Mercer条件。因为RBF核函数 在非线性函数预测方面表现出很好的特性，本发明选用RBF核函数作为核函数进行预测：
[00W] 2)由LSSVM原理可知，在使用最小二乘支持向量机建模过程中，惩罚参数和核参 数对模型性能和泛化能力有很大的影响，本发明采用效果较好的粒子群算法来优化对参数 的选择。通过对标准粒子群优化算法（PSO)算法分析可知，标准PSO算法很容易陷入局部 最优，导致算法进化后期停滞不前，为了解决该个问题，本发明采用高斯扰动的粒子群优化 算法（GDPSO)进行参数寻优。此算法在粒子运动过程中增加高斯扰动项，丰富种群的多样 性，提高粒子寻找到最优可行解的概率。高斯扰动的范围会随着迭代的进行逐渐减小，并根 据不同粒子上的不同维度信息进行不同程度的扰动，具有一定的动态随机性。将高斯扰动 项加入速度迭代公式中，得到
[0062] Vi,j(t+1) =wVi,j(t) +(3江1 虹,j(t) +r2gausSi,j(t) -Xi,j(t)]
[0063] +C2r3 [pg_j(t)-Xij(t)]
[0064] Xi j(t+1) =Xi j(t)+Vij(t+l)
[0065]gauss；j(t) =r4gaussian(y,5
[0066] 其中，w为惯性权重，Cl、C2为加速因子，r1、r2、rs、r4表示化1)区间服从均匀分 布的随机数，Vu(t)为粒子i在第t次迭代时的速度，Pu(t)为粒子i在第t次迭代时的 历史最优位置，gaussu(t)为粒子i在第t迭代时产生的高斯扰动，y为均值，S2为方差， xu(t)为粒子i在第t次迭代时的位置，