一种电力数据补全方法
【技术领域】
本发明涉及电力技术领域,具体涉及一种电力数据补全方法。
【背景技术】
电力负荷预测作为电力系统重要工作内容之一,预测的精度高低会对经 济调度、实时控制、运行计划及发展规划等方面产生很大影响。随着电力市 场进一步开放、市场机制更加成熟、市场开放更加全面,电力负荷预测将发 挥更大作用。传统的负荷预测方式例如回归预测法、时间序列法等方法在实 时性和准确性方面已经不能满足目前的需求,尤其是在多重因素的影响下面 临巨大挑战。
电力系统作为城市生命线工程的重要环节,对于城市的发展建设和居民 生活保障都发挥着重要作用,电力系统的安全和稳定是一切社会活动的基础。 一方面随着国家科技进步、社会现代化进程的加快,工程建设、企业生产和 国民生活对电能的依赖性更强,电力消耗巨大带来的需求量进一步增加,同 时对于供电质量以及可靠性安全性等方面也会有更高的要求。因此,通过智 能电网建设可以显著提升电力系统安全稳定运行水平以及供电的可靠性,这 样不仅有利于智能化水平和用户互动化水平的提升,在缓解能源压力以及资 源有效优化的配置等方面同样受益良多。负荷预测作为现阶段智能电网建设 中的一项重要工作,精确的负荷预测对于保证电力系统的安全、稳定、经济 运行起到举足轻重的作用,其工程价值以及理论意义也有深远的影响。与此 同时,随着电力市场化交易改革进一步深入,交易规模也在进一步扩大,其 中在现货市场、输配电价以及增量配电等等领域稳中有进,从政策层面上对 市场的配套设施、辅助服务也提供一定的帮助,电力市场改革所带来的巨大 经济效益以及社会效益会集中体现,这一系列的机构设立和政策出台标志着有竞争力的市场结构以及市场体系逐步形成,市场的力量也将为电力行业带 来新的活力与生机。
电力负荷预测逐步成为研究人员关注的重点方向,其水平的高低甚至可 以作为衡量相关电力企业的管理水平是否达到现代化的标准,不同时间周期 负荷预测的精度对电网稳定安全经济运行起到重要的保障。尤其在电力事业 发展速度空前的今天,用电管理的方式逐步与市场对接、用电管理走向市场, 电力负荷预测问题的解决已经成为重要的任务,由此可见,有必要针对电力 负荷预测工作展开深入探讨和研究。
【发明内容】
为解决前述问题,本发明提供了一种电力数据补全方法,以准确地将缺 失的电力数据补充完整。
为了达到上述目的,本发明采用如下技术方案:
建立GABP神经网络;
利用缺失的电力负荷数据之前的电力负荷数据对GABP神经网络进行训 练,得到第一GABP神经网络,利用第一GABP神经网络对缺失的电力负荷 数据进行计算,得到第一补全数据;
利用缺失的电力负荷数据之后的电力负荷数据对GABP神经网络进行训 练,得到第二GABP神经网络,利用第二GABP神经网络对缺失的电力负荷 数据进行计算,得到第二补全数据;
建立适应度函数,利用适应度函数对第一补全数据和第二补全数据进行 计算,得到第三补全数据;
建立PSO-SVR模型,并训练所述PSO-SVR模型,利用训练好的PSO-SVR 模型对所述第三补全数据进行计算,得到完整的电力负荷数据。
进一步地,所述建立GABP神经网络具体包括如下步骤:
建立BP神经网络模型,所述BP神经网络模型包括输入层、隐含层和输 出层,隐含层的第i个神经元的输入为:
其中,x
p为输入样本数据,t
p为输出样本数据,w
ij表示隐含层第i个神 经元与输入层第j个神经元的连接权值;θ
i表示隐含层中第i个神经元阈值;
表示输入层第j个神经元的输入值,
表示输入层第j个神经元的输出值,M 代表输入层神经元个数;
隐含层第i个神经元输出为:
其中,f(x)为传递函数;
隐含层所含神经元输出的结果与连接权值相乘后,将所有数值相加输入 到输出层,所述输出层中,第k个神经元的输入为:
其中, w
ki为输出层第k个神经元与隐含层第i个神经元的连接权值,
表示输入层 第i个神经元的输出值,θ
k表示隐含层中第k个神经元阈值;
所述输出层中将输入数值带入所述传递函数中,得到所述输出层的输出 结果为
作为优选,建立适应度函数,利用适应度函数对第一补全数据和第二补 全数据进行计算,得到第三补全数据具体包括如下步骤:
依据适应度函数分别求解不同种群的适应度,所述适应度函数为
其中,K为输出层中神经元个数,M’为输入样本数 据总数,y
pi为BP神经网络的预期输出,y
oi表示BP神经网络的实际输出;
再依据以下两条公式计算个体被选中的概率pi以及种群累计概率ppi:
其中,M代表种群中全部个体的数量;
对种群累计概率高的群体中的个体进行随机配对,在配对个体中随机选 择交叉点的位置,将两个个体的部分结构进行替换,重组而生成新的个体;
将新的个体在变异点上的基因值按照变异概率取反,所述取反为1变0, 0变1,所述变异率的取值范围在0.0001~0.1之间。
作为优选,所述最PSO-SVR模型为
其中,y为PSO-SVR模型的输出数据,x为PSO-SVR模型的输入数据,l表示维度的总 数,K(x
i,x)表示输入数据到高维特征空间的非线性映射,η
i表示第i个维度 上的参数η,h表示待定常数,训练所述PSO-SVR模型为训练参数η
i和参数h。
本发明所提供的方法具有如下有益效果:
(1)本申请所提供的方法,采用双向的GABP神经网络算法,相较于单 向的补全方法,不仅利用到缺失点之前的负荷数据,还利用到缺失点之后的 数据,因此使得补全的数据更加精准。
(2)本申请所提供的方法,引入PSO-SVR模型,将两个方向的计算结 果融合,取代了传统的加权求和方式,这样将缺失之后的数据变化趋势同样 考虑进去,进一步保证了数据补全的准确性。
本发明的这些特点和优点将会在下面的具体实施方式以及附图中进行详 细的揭露。本发明最佳的实施方式或手段将结合附图来详尽表现,但并非是 对本发明技术方案的限制。另外,在每个下文和附图中出现的这些特征、要 素和组件是具有多个,并且为了表示方便而标记了不同的符号或数字,但均 表示相同或相似构造或功能的部件。
【附图说明】
下面结合附图对本发明作进一步说明:
图1为本发明实施例的流程图;
图2本发明实施例中补全数据的流程图。
【具体实施方式】
下面结合本发明实施例的附图对本发明实施例的技术方案进行解释和说 明,但下述实施例仅为本发明的优选实施例,并非全部。基于实施方式中的 实施例,本领域技术人员在没有做出创造性劳动的前提下所获得其他实施例, 都属于本发明的保护范围。
在本说明书中引用的“一个实施例”或“实例”或“例子”意指结合实 施例本身描述的特定特征、结构或特性可被包括在本专利公开的至少一个实 施例中。短语“在一个实施例中”在说明书中的各位置的出现不必都是指同 一个实施例。
实施例:
一种电力负荷预测方法,包括如下步骤:
以0到30min内任意时长为时间间隔,采集电力负荷数据;
收集若干天内的电力负荷数据,判断若干天内的电力负荷数据是否缺失, 若电力负荷数据缺失,则进行补全,具体步骤如下:
建立GABP神经网络,具体包括:建立BP神经网络模型,BP神经网络 模型包括输入层、隐含层和输出层,隐含层的第i个神经元的输入为:
其中,x
p为输入样本数据,t
p为输出样本数据,w
ij表示隐含层第i个神 经元与输入层第j个神经元的连接权值;θ
i表示隐含层中第i个神经元阈值;
表示输入层第j个神经元的输入值,
表示输入层第j个神经元的输出值,M 代表输入层神经元个数;
隐含层第i个神经元输出为:
其中,f(x)为传递函数,常用 的传递函数包括purelin函数、tansig函数等;
隐含层所含神经元输出的结果与连接权值相乘后,将所有数值相加输入 到输出层,输出层中,第k个神经元的输入为:
其中,w
ki为输出层第k个神经元与隐含层第i个神经元的连接权值,
表示输入层第i 个神经元的输出值,θ
k表示隐含层中第k个神经元阈值;
输出层中将输入数值带入传递函数中,得到输出层的输出结果为
利用缺失的电力负荷数据之前的电力负荷数据对GABP神经网络进行训 练,得到第一GABP神经网络,利用第一GABP神经网络对缺失的电力负荷 数据进行计算,得到第一补全数据,这一步骤的主要目的是挖掘负荷数据缺 失点的前向趋势信息,并获得前向GABP神经网络的基础;
利用缺失的电力负荷数据之后的电力负荷数据对GABP神经网络进行训 练,得到第二GABP神经网络,利用第二GABP神经网络对缺失的电力负荷 数据进行计算,得到第二补全数据,这一步骤使用反向排序的加载数据来训 练模型。该步骤可以挖掘负载数据缺失点后的数据变化趋势,从而获得基于 反向GABP神经网络的负载数据互补结果。
建立适应度函数,利用适应度函数对第一补全数据和第二补全数据进行 计算,得到第三补全数据,具体步骤包括:依据适应度函数分别求解不同种 群的适应度,不同的种群分别对应利用第一补全数据、第二补全数据取值的 集合,而取值的集合为运算的中间步骤,是根据输入来确定,并非固定的值 或固定的范围,适应度函数为
其中,K为输出层中 神经元个数,M’为输入样本数据总数,y
pi为BP神经网络的预期输出,即BP 神经网络的预测值,y
oi表示BP神经网络的实际输出,即训练BP神经网络的 数据里用来检验预测结果的值;
再依据以下两条公式计算个体被选中的概率pi以及种群累计概率ppi:
对种群累计概率高的群体中的个体进行随机配对,随机变化中,适应度 函数值较高的个体会有更大概率遗传到下一代,这个概率与适应度函数值呈 正比例关系,在配对个体中随机选择交叉点的位置,将两个个体的部分结构 进行替换,重组而生成新的个体,替换的部分结构为个体编码是0、1组成的 长度为l的字符串,将变异点的字符串相互交换位置,变异点即随机选择的点 的字符串;
将新的个体在变异点上的基因值按照变异概率取反,取反为1变0,0变 1,变异率的取值范围在0.0001~0.1之间。本领域技术人员均知晓,在遗传算 法中,配对、交叉、变异的操作是完全随机的,类似于物种进化过程,通过 反复的运算规则取得最终效果。
建立并训练PSO-SVR模型,PSO-SVR模型为
其 中,y为PSO-SVR模型的输出数据,x为PSO-SVR模型的输入数据,l表示 维度的总数,K(x
i,x)表示输入数据到高维特征空间的非线性映射,本领域技 术人员均知晓,非线性映射过程通过一个核函数运算实现,核函数的形式可 以是多样的,最简单的形式为求x
i与x两个特征向量的内积,η
i表示第i个维 度上的参数η,h表示待定常数,训练PSO-SVR模型为利用第一补全数据和第 二补全数据训练参数η
i和参数h,利用训练好的PSO-SVR模型对第三补全数 据进行计算,得到完整的电力负荷数据。
PSO-SVR模型的作用是将第一GABP神经网络和第二GABP神经网络结 合起来,对正向数据和反向数据产生影响。为了实现数据结构的双向组合, 需要考虑因素。一般组合模型在权重分配时通常采用简单的线性加权或平均 法。这种简单的处理方法使得权重分布和计算方法不合理,导致预测结果不 准确。在PSO-SVR模型中,输入数据为第一GABP神经网络负载数据完成结 果和第二GABP神经网络负载数据完成结果,即第一补全数据和第二补全数 据,输出数据为实际负载数据值,通过PSO-SVR模型计算最优权重,进一步 保证了数据补全的准确性。
获得完整的电力负荷数据以后,建立最小二乘支持向量机模型,训练最 小二乘支持向量机模型,最小二乘支持向量机模型为
其中,y为最小二乘支持向量机模型的输出数据,x为最小二乘支持向量机模 型的输入数据,具体为
格式的若干天内的电力负荷数据,其中,m 为若干天内的电力负荷数据的天数,t为采集数据时间间隔的时长,l表示维 度的总数,K'(x
i,x)表示输入数据到高维特征空间的非线性映射,本领域技术 人员均知晓,非线性映射过程通过一个核函数运算实现,核函数的形式可以 是多样的,最简单的形式为求x
i与x两个特征向量的内积,λ
i表示第i个维度上的参数λ,q表示待定常数,训练最小二乘支持向量机模型为利用若干天内 完整的电力负荷数据训练参数λ
i和参数。
建立容积卡尔曼滤波预测模型,利用若干天内完整的电力负荷数据训练 容积卡尔曼滤波预测模型,容积卡尔曼滤波预测模型为
其中,下标k表示时刻,x
k为k时刻 的电力负荷数据,y
k为k时刻的电力负荷预测数据,均采用格式为1*n的若干 天内的电力负荷数据,n为若干天内的电力负荷数据。y
k=H(x
k)+v
k,H(x
k)为 非线性函数,是基于测量数据拟合得到的,输入不同的测量数据拟合出来的 具体函数关系不同;v
k为量测噪声,有不同的模型进行表示,例如正态分布,
为待训练的增益矩阵,a、b、c分别为状态方程的参数,状态 方程为x
k=ax
k-1+bx
k-2+cx
k-3+d。
建立灰色神经网络模型,训练灰色神经网络模型,具体包括如下步骤:
建立灰色神经网络结构、参数a’、b’、c’以及初始权值,灰色神经网 络包括LA层、LB层、LC层以及LD层,初始权值为:
ω11=a',ω21=-y1(0),ω22=u1.....ω2n=un-1
ω31=ω32=...=ω3n=1+ea't
其中,ω11,ω21,ω22,.....ω2n,ω31,ω32,......ω3n为灰色神经网络的权值,
参数a’为LA层的参数,计算公式为:a'=ω11t;
参数c’为LC层的参数,计算公式为:c1'=b'ω21,c2'=y2(t)b'ω22,....cn'=yn(t)b'w2n, 其中,t为输入,y2(t)至yn(t)为第n个节点对t时刻的预测值;
LD层具有参数d’,计算公式为:d'=ω31c'1+ω32c'2+...+ω3nc'n-θy1
采用最小二成支持向量机模型的计算结果和容积卡尔曼滤波预测模型的 计算结果对灰色神经网络进行训练,计算预测误差,根据误差调整权值,误 差的比较基准根据业务需要而定,比如浙江省电力市场运行规则所引导的盈 利目标,不断变化;
利用训练好的最小二乘支持向量机模型对若干天内的电力负荷数据进行 计算,得到第一预测结果;
利用训练好的容积卡尔曼滤波预测模型对若干天内的电力负荷数据进行 计算,得到第二预测结果;
利用训练好的灰色神经网络模型对第一预测结果和第二预测结果进行计 算,得到最终预测结果,此步骤中,采用最终预测方程进行计算,得到最终 预测结果,最终预测方程为:
本实施例所提供的技术方案,发挥了最小二乘支持向量机预测算法泛化 能力强的优点,在数据量较少的情况下依然能进行准确预测,同时发挥了容 积卡尔曼滤波预测算法实时性强、具有一定矫正能力的特点,并且,采用灰 色神经网络对两种算法的预测结果进行结合,有效实现算法之间的优势互补, 进一步提高负荷预测的精度。
以上,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此, 熟悉该本领域的技术人员应该明白本发明包括但不限于附图和上面具体实施 方式中描述的内容。任何不偏离本发明的功能和结构原理的修改都将包括在 权利要求书的范围中。