基于相似日聚类的微电网光伏出力空间预测方法
【专利摘要】本发明针对现有相似日聚类预测光伏出力的不足，提出一种基于相似日聚类的微电网光伏出力空间预测方法，属于新能源出力预测领域，包含步骤：利用训练期内微电网中输入输出光伏电站的发电数据，建立空间关系来预测预测期内输出电站的发电数据；对输入电站训练期和预测期内的发电数据进行粒子化处理；对粒子化后的输入电站训练期内的数据进行SOM神经网络聚类训练；利用SOM神经网络对粒子化后的输入电站预测期进行仿真分类；找出相似日；利用相似日数据，结合最小二乘法多项式拟合，得出发电数据预测值。与现有方法相比，本发明简化了研究对象，促进了聚类过程，显著提高了预测精度。
【专利说明】
基于相似日聚类的微电网光伏出力空间预测方法
技术领域
[0001] 本发明属于微电网光伏出力预测技术领域，具体涉及一种基于相似日聚类的微电 网光伏出力空间预测方法。
【背景技术】
[0002] 作为智能电网重要组成部分的微电网，是由分布式电源、储能、负荷、能量转换装 置和监控、保护装置组合而成的小型发配电系统，其既可以与外部电网并网运行，也可以孤 立运行，具有自我控制、保护和管理功能。
[0003] 太阳能光伏发电是微电网中常见的一种微电源发电形式。作为一种无污染的清洁 能源，一种可持续的能源替代方式，太阳能光伏发电已得到人们越来越多的重视，其已作为 世界上许多国家开辟新能源的一条重要途径。
[0004] 为实现微电网的自我控制、保护和管理功能，对作为微电源的光伏发电系统出力 进行准确有效的预测必不可少。然而光伏发电系统受环境因素影响明显，发电功率的变化 存在着间歇性、波动性和随机性，准确有效的预测变得困难和复杂。
[0005] 目前，有关使用S0M神经网络进行相似日聚类来进行光伏出力预测的国内外研究 中，相似日的聚类大多直接依据天气环境因素，而从空间角度直接依据与目标电站相关联 的其它电站的发电量非常罕见，如文献《基于天气类型聚类识别的光伏系统短期无辐照度 发电预测模型研究》（出处：中国电机工程学报，2011年第31卷34期28-35页）就是直接依据 天气环境因素来进行聚类。直接依据天气环境因素，存在影响因素复杂、获取数据困难和检 测精度不高等问题。另外，针对发电量数据变动范围大，影响S0M神经网络分类性能的情况， 提出了一种粒子化数据处理方法，对发电量数据进行粒子化处理，处理后的数据量纲清晰、 层次分明，便于S0M神经网络分类。
【发明内容】
[0006] 本发明的目的在于，针对上述现有技术的不足，避免考虑复杂的天气环境因素，从 空间角度直接依据微电网内相关联光伏电站的发电量数据，并采用粒子化数据处理方法， 对发电量数据进行粒子化处理，采用S0M神经网络进行相似日聚类，结合最小二乘法多项式 拟合，对预测期内目标输出电站的发电数据进行预测。空间角度简化了研究对象，粒子化法 促进了聚类过程，相似日聚类显著提高了预测精度，为基于相似日的微电网光伏出力预测 提供了较高的技术参考价值。
[0007] 为实现上述目的，本发明采用的技术方案是：
[0008] 提供一种基于相似日聚类的微电网光伏出力空间预测方法，其特征在于，包含以 下步骤：
[0009] 1)以微电网内的η个光伏电站为研究对象(其中p个作为输入电站，q个作为输出电 站，n = p+q)，n个光伏电站的发电数据（日发电量)分为训练期和预测期，训练期内输入电站 和输出电站的发电数据均已知，预测期内输入电站的发电数据已知而输出电站的发电数据 未知，利用训练期的发电数据，从空间角度建立输入电站和输出电站之间的关系，利用这种 关系则可预测预测期内输出电站的发电数据；
[0010] 2)对P个输入电站训练期和预测期内的发电数据进行粒子化处理，具体过程如下：
[0011] 201)将p个输入电站训练期和预测期内的发电数据(w维），按各个输入电站时间顺 序分别归一化到[0，a]，a为大于0的整数，具体使用下式(1):
[0013] 其中，m= 1，2,…，p，i = 1，2,…，w，xmi为第m个光伏电站的第i个发电数据，Xmi*为 Xmi归一化后对应的数据，Xm max和Xm min分别为Xmi中最大值和最小值；
[0014] 202)将Xmi1f确到整数位或小数点后1位，具体使用下式(2):
[0015] f(xmi*)=R0UND(xmi*,N) (2)
[0016] 其中，R0UNDO为四舍五入运算，其中N=0，1，当N为0时表示精确到整数位，N为1时 表示精确到小数点后1位；
[0017] 3)对粒子化后的p个输入电站训练期内的数据进行S0M神经网络聚类训练，具体过 程如下：
[0018] 301)将粒子化后的p个输入电站训练期内的数据作为输入向量X;
[0019] 302)初始化S0M神经网络的连接权值向量W为随机值，W表示S0M输入层神经元和输 出层神经元之间的全连接；
[0020] 303)将输入向量X和连接权值向量W分别进行归一化处理，具体使用下式(3):
[0022]其中，| | · | |表示2-范数运算；
[0023] 304)计算归一化后输入向量Y和归一化后连接权值向量F之间的欧氏距离d;
[0024] 305)选出d值最小对应的输出层神经元为获胜神经元，对获胜神经元及其邻域内 的神经元权值进行调整；
[0025] 306)判断是否达到预先训练结束要求，否则循环执行，得出最终训练分类；
[0026] 4)利用步骤3)中达到预先训练结束要求的S0M神经网络对粒子化后的输入电站预 测期进行仿真分类；
[0027] 5)找出步骤4)中仿真分类和步骤3)中训练分类相同的类别，这些类别所对应的日 期则可作为相似日；
[0028] 6)利用相似日粒子化前的发电数据，结合最小二乘法多项式拟合，得出预测期内 输出电站的发电数据预测值，具体过程如下：
[0029] 601)利用与预测期为相似日的训练期内输入输出电站粒子化前的发电数据，使用 最小二乘法多项式拟合，求出多项式中的系数，所使用的最小二乘法多项式拟合原理是:对 于给定的一组数据{(幻，^)，（」=1，2，~，8)}，其中8为相似日的个数，若拟合的曲线为 7 = ? (X)，则第j个点的误差距离为F(Xj)-^，所有点的误差距离平方和
的最小值对应的系数，则可得到拟合曲线y=F(X);
[0030] 602)利用拟合曲线y = F(x)，输入预测期内输入电站粒子化前的发电数据，则可求 出预测期内输出电站的发电数据预测值。
[0031] 本发明采用的以上技术方案与现有技术相比，可达到以下技术效果：
[0032] (1)预测微电网内的光伏出力，从空间角度直接依据与目标光伏电站相关联的其 它光伏电站的出力，简化了研究对象，避免了复杂的天气环境因素；
[0033] (2)S0M神经网络是一种无监督、自学习网络，无需外界给出评价函数，能根据数据 自身特征进行分类，具有适应性强、自稳定性好、处理数据量大、拓扑有序等特点；
[0034] (3)针对发电量数据变动范围大，影响S0M神经网络分类性能的情况，提出了一种 粒子化数据处理方法，对发电量数据进行粒子化处理，处理后的数据量纲清晰、层次分明， 便于S0M神经网络分类；
[0035] (4)利用粒子化后的数据，进行S0M神经网络相似日聚类，结合最小二乘法多项式 拟合，显著提高了预测精度，为基于相似日的微电网光伏出力预测提供了较高的技术参考 价值。
【附图说明】
[0036] 图1是本发明所述的一种基于相似日聚类的微电网光伏出力空间预测方法总体框 架示意图；
[0037] 图2是其中一个输入电站A训练期和预测期内的发电数据进行粒子化处理前的数 据图；
[0038] 图3是其中一个输入电站A训练期和预测期内的发电数据进行粒子化处理后的数 据图；
[0039] 图4是S0M神经网络对输入电站训练期内的数据进行分类的结果图；
[0040] 图5是分出的9个类中每个类的个数图。
【具体实施方式】
[0041] 下面结合附图和实施例对本发明作详细说明，需要说明的是，该具体实施例只用 来解释本发明，而非对本发明范围的限定。
[0042] 图1是本发明所述的一种基于相似日聚类的微电网光伏出力空间预测方法总体框 架示意图，如图1所示，一种基于相似日聚类的微电网光伏出力空间预测方法具体包含以下 步骤：
[0043] 1)以微电网内的η个光伏电站为研究对象(其中p个作为输入电站，q个作为输出电 站，n = p+q)，n个光伏电站的发电数据（日发电量)分为训练期和预测期，训练期内输入电站 和输出电站的发电数据均已知，预测期内输入电站的发电数据已知而输出电站的发电数据 未知，利用训练期的发电数据，从空间角度建立输入电站和输出电站之间的关系，利用这种 关系则可预测预测期内输出电站的发电数据。
[0044] 本实施例中，η为7，p为6，q为1，训练期为2013年1月1日至2014年11月30日，预测期 为2014年12月1日。
[0045] 2)对p个输入电站训练期和预测期内的发电数据进行粒子化处理，具体过程如下：
[0046] 201)将p个输入电站训练期和预测期内的发电数据(w维），按各个输入电站时间顺 序分别归一化到[0，a]，a为大于0的整数，具体使用下式(1):
[0048] 其中，m= 1，2,…，p，i = 1，2,…，w，xmi为第m个光伏电站的第i个发电数据，XmZ为 Xmi归一化后对应的数据，Xm max和Xm min分别为Xmi中最大值和最小值。
[0049] 本实施例中，w为700，a为10。
[0050] 202)将Xmi1f确到整数位或小数点后1位，具体使用下式(2):
[0051] f(xmi*)=ROUND(xmi*,N) (2)
[0052] 其中，ROUND 〇为四舍五入运算，其中N=0，1，当N为0时表示精确到整数位，N为1时 表示精确到小数点后1位。
[0053] 本实施例中，因粒子化归一区间为[0，10]，区间间隔范围大，可将Xm，精确到整数 位，即N为0。
[0054]图2所示是其中一个输入电站A训练期和预测期内的发电数据进行粒子化处理前 的数据图，图3所示是其对应的粒子化处理后的数据图。可见处理后的数据量纲清晰、层次 分明，便于S0M神经网络分类。
[0055] 3)对粒子化后的p个输入电站训练期内的数据进行S0M神经网络（自组织映射网 络)聚类训练，具体过程如下：
[0056] 301)将粒子化后的p个输入电站训练期内的数据作为输入向量X。
[0057] 302)初始化SOM神经网络的连接权值向量W为随机值，W表示SOM输入层神经元和输 出层神经元之间的全连接。
[0058] 303)将输入向量X和连接权值向量W分别进行归一化处理，具体使用下式(3):
[0060]其中，| | · | |表示2-范数运算。
[0061 ] 304)计算归一化后输入向量Y和归一化后连接权值向量F之间的欧氏距离d。
[0062] 305)选出d值最小对应的输出层神经元为获胜神经元，对获胜神经元及其邻域内 的神经元权值进行调整。
[0063] 306)判断是否达到预先训练结束要求，否则循环执行，得出最终训练分类。
[0064] 本实施例中，SOM神经网络输出层神经元个数为9,即SOM神经网络可将输入电站训 练期内的发电数据分成9个类。训练次数为100,当SOM神经网络训练达到100次时训练结束。 图4所示为SOM神经网络对输入电站训练期内的数据进行分类的结果图。图5所示为分出的9 个类中每个类的个数图。
[0065] 4)利用步骤3)中达到预先训练结束要求的SOM神经网络对粒子化后的输入电站预 测期进行仿真分类。
[0066] 5)找出步骤4)中仿真分类和步骤3)中训练分类相同的类别，这些类别所对应的日 期则可作为相似日。
[0067] 本实施例中，对输入电站预测期的仿真分类为1，类别为1的相似日天数为74,如图 5所示。为保证训练分类相同的类别所对应的相似日天数不至于太少(小于10)，要有充足的 样本数据或者减小SOM神经网络分类的数目。
[0068] 6)利用相似日粒子化前的发电数据，结合最小二乘法多项式拟合，得出预测期内 输出电站的发电数据预测值，具体过程如下：
[0069] 601)利用与预测期为相似日的训练期内输入输出电站粒子化前的发电数据，使用 最小二乘法多项式拟合，求出多项式中的系数，所使用的最小二乘法多项式拟合原理是:对 于给定的一组数据{(幻，^)，（」=1，2，~， 8)}，其中8为相似日的个数，若拟合的曲线为7 = ?
(X)，则第j个点的误差距离为F(Xj)-^，所有点的误差距离平方和
最小值对应的系数，则可得到拟合曲线y = F(x)。
[0070] 602)利用拟合曲线y = F(x)，输入预测期内输入电站粒子化前的发电数据，则可求 出预测期内输出电站的发电数据预测值。
[0071] 本实施例中，对于6个输入电站来说，求取目标输出电站发电量数据预测值的多项 式表达式为y = am+a2X2+a3X3+a4X4+a5X5+a6X6，其中X1-X6为预测期内已知的6个输入电站的 发电数据，y为预测期内未知的输出电站的发电数据， ai-a6为系数。求出的最小二乘法多项 式的系数如表1所示。
[0072]预测性能参数采用绝对百分误差APE，其具体公式为APE=(|y-t|/t)X100%，其 中y为预测值，t为实际值。预测值y和实际值t以及APE如表2所示(精确到小数点后四位）。 [0073]本实施例还增加了对比试验，发电数据不经过粒子化处理直接进行S0M神经网络 分类后，使用最小二乘法多项式拟合预测，效果是：S0M分类时间加长且预测APE值偏大，预 测性能明显低于本发明所提方法。
[0074]表 1
[0076]以上所述仅为本发明的优选实施例，并不用于限制本发明。凡在本发明的思想和 原理内做的任何轻易想到的修改、改进或替换等，均应包含在本发明的保护范围内。
【主权项】
1. 一种基于相似日聚类的微电网光伏出力空间预测方法，其特征在于，包含以下步骤： 1) 以微电网内的η个光伏电站为研究对象(其中P个作为输入电站，q个作为输出电站，η = P+q)，n个光伏电站的发电数据（日发电量)分为训练期和预测期，训练期内输入电站和输 出电站的发电数据均已知，预测期内输入电站的发电数据已知而输出电站的发电数据未 知，利用训练期的发电数据，从空间角度建立输入电站和输出电站之间的关系，利用这种关 系则可预测预测期内输出电站的发电数据； 2) 对p个输入电站训练期和预测期内的发电数据进行粒子化处理，具体过程如下： 201) 将p个输入电站训练期和预测期内的发电数据(w维），按各个输入电站时间顺序分 别归一化到[〇，a]，a为大于0的整数，具体使用下式（1):其中，m=l，2,…，p，i = l，2,…，w，xmi为第m个光伏电站的第i个发电数据，XmZ为xmi归一 化后对应的数据，Xm max和Xm min分别为Xmi中最大值和最小值； 202) 将Xmi1f确到整数位或小数点后1位，具体使用下式(2): f(xmi*) = ROUND (xmi*,N) (2) 其中，ROUND()为四舍五入运算，其中N = 0，1，当N为0时表示精确到整数位，N为1时表 示精确到小数点后1位； 3) 对粒子化后的p个输入电站训练期内的数据进行S0M神经网络聚类训练，具体过程如 下： 301) 将粒子化后的p个输入电站训练期内的数据作为输入向量X; 302) 初始化S0M神经网络的连接权值向量W为随机值，W表示S0M输入层神经元和输出层 神经元之间的全连接； 303) 将输入向量X和连接权值向量W分别讲行丨]3-化々卜理，具体使用下式(3):其中，11 · 11表示2-范数运算； 304) 计算归一化后输入向量Y和归一化后连接权值向量w之间的欧氏距离d; 305) 选出d值最小对应的输出层神经元为获胜神经元，对获胜神经元及其邻域内的神 经元权值进行调整； 306) 判断是否达到预先训练结束要求，否则循环执行，得出最终训练分类； 4) 利用步骤3)中达到预先训练结束要求的S0M神经网络对粒子化后的输入电站预测期 进行仿真分类； 5) 找出步骤4)中仿真分类和步骤3)中训练分类相同的类别，这些类别所对应的日期则 可作为相似日； 6) 利用相似日粒子化前的发电数据，结合最小二乘法多项式拟合，得出预测期内输出 电站的发电数据预测值，具体过程如下： 601)利用与预测期为相似日的训练期内输入输出电站粒子化前的发电数据，使用最小 二乘法多项式拟合，求出多项式中的系数，所使用的最小二乘法多项式拟合原理是:对于给 定的一组数据{(Xj，yj)，（ j = l，2，···，s)}，其中s为相似日的个数，若拟合的曲线为y = F(x)， 则第j个点的误差距离为F(Xj)-yp所有点的误差距离平方和为：t|>(七)-A]2,求出 -A]2的最小值对应的系数，则可得到拟合曲线y=F( X); l 602)利用拟合曲线y = F(x)，输入预测期内输入电站粒子化前的发电数据，则可求出预 测期内输出电站的发电数据预测值。
【文档编号】G06Q50/06GK105868916SQ201610222346
【公开日】2016年8月17日
【申请日】2016年4月5日
【发明人】付青, 单英浩, 耿炫
【申请人】中山大学