(国网新疆电力有限公司经济技术研究院，新疆 乌鲁木齐 830049)

0 引 言

随着生态文明建设的持续深入，电力作为国家核心战略能源，公众对其的关注度日益增强。目前，中国的经济实力在全球经济带中的作用突出。为了保证经济平稳、有序、协调发展，中国的电力电网规划就需要与建设同步协调发展，而中长期的电力电量预测作为系统规划及运行中重要的组成部分，对电力电网制定战略发展规划、营销策略、合理资源配置和促进经济社会发展等方面具有不可替代的重要意义[1-2]。为适应当前新形势、新环境下的电力电量变化，就要求中长期电量预测在分析历史变化的基础上，能够有效自适应规律的变化，分析模拟出预测对象与相关因素的关系。一个良好的预测模型应当能够做出这种变化规律描述[3]。

中长期电量预测的对象往往都是以地区级为单位，其时间序列长度也是以年度逐年变化。一般传统的预测方法主要以趋势外推法为主，常用的方法包括时间序列法和灰色理论等[4-5]。时间序列法侧重于目标自身的特征变化规律而忽视了相关因素的影响作用。传统的灰色模型GM(1，1)主要用于波动变化不大，变化趋势体现为指数增长的样本数据，但随着中国经济社会发展受全球化影响，已进入新形势、新环境和新常态，这种预测模型的方式已不再能够满足精度需求。文献[6]中提出了基于粒子群优化传统GM(1，1)模型的预测方法。文献[7]将偏最小二乘回归分析理论用于对电力负荷的中长期预测中。文献[8]提出了一种多变量残差修正进行多变量预测的方法，但是该方法的预测精度难以把控。文献[9]利用不同维度的灰色预测模型训练样本，结合神经网络进行预测。文献[10]组合优化了5种不同预测方法的实际应用性，提出了一种优化的组合预测模型。目前，国内主要的预测方向主要集中在短期电力负荷电量的预测，在中长期的电力电量预测中尚有空缺[11]。

每个地区的电力总量是由不同产业电量及居民用电组成，不同的用电主体总是呈现出不同用电特质，相应的用电规律也有所不同，但都具有各自的规律。通过挖掘分析不同产业以及GDP对用电总量的影响，考虑多方面的关联因素，改进传统MGM(1，n)预测模型中背景值的计算方法，使预测模型能够自适应当前新形势下的变化规律，取得了较好的预测结果。

1 改进的多变量灰色预测模型

1.1 灰色预测模型

目前常用的回归分析、神经网络等预测模型主要是应用于多数据样本的模型预测中。对于小样本数据因其规律变化难以突显，因而无法应用。为解决小样本数据的研究，邓聚龙教授提出了灰色理论[12]，灰色预测模型就是以基于小样本数据的规律特性为主要研究对象。为弱化小样本数据的随机特性，灰色预测模型通常采用相邻累加使数据产生出有规律的趋势，从而生成新的白色数据序列。

灰色预测模型一般分为两类方程：一类是离散型的灰色微分方程；另一类是连续型的灰色微分方程。预测模型主要是通过数据信息的分析、比对，得到主变量的变化规律完成预测。在两类方程选用中，主要取决于数据类型及特征。对于不适宜建模的数据样本，需要对样本进行必要的预处理，使其符合模型的基本要求从而进行预测[13-14]。灰色模型的构建原理如图1所示。

1.2 改进的多变量灰色预测模型

多变量灰色预测模型MGM(1，m)是灰色理论中的一种，其中m表示变量的个数。该模型能够较好地体现出预测变量体系中各因素量之间相互制约和相互发展的关联性，可以从整体的角度统一描述[15]。将该模型用于中长期电力电量预测中，相对于仅考虑单因素的GM(1，1)而言，有更高预测精度；同时，MGM预测模型可优化的方向很多。下面针对MGM背景值进行了优化改进。改进的MGM背景值计算，不再采取传统的相邻数据累加求均值的计算方法，而是基于省级电力电量特性及敏感分析引入权重系数和规律重构计算公式，强化了背景值对相邻2个节点数据规律特征的突显，使得改进后的预测方法计算出的背景值能够适应当前电量变化的规律趋势。

图1 灰色模型构建原理

基于MGM(1，m)模型的矩阵形式为

(1)

建立预测模型为

(2)

式中：

(3)

解式(1)得对应的时间响应函数为

X(1)(t)=eA(t-1)[X(1)(1)+A-1B]-A-1B

(4)

将式(4)进行离散化处理：

(j=1,2,…,m;k=2,3,…,n)

(5)

通过最小二乘法进行模型估计，计算各参数估计值：

=(PTP)-1PT(Q1,Q2,…,Qm)

(6)

式中：

(7)

(8)

预测模型的时间相量表示为

(9)

还原得到初始序列目标模型为

=X(1)(k)-X(1)(k-1)(k=2,3,…,n)

(10)

改进的多变量灰色模型的预测流程主要分为一次累加、背景值计算、最小二乘法相关参数值计算、时间相应表示、还原初始向量从而得到预测结果，具体流程如图2所示。

图2 改进的MGM预测流程

1.3 模型检验

预测模型结果的检验方法主要有3种[16]，分别是残差检验、关联度检验以及后验差检验。

1)残差检验

残差检验主要是基于平均相对误差的一种检验方法。首先，需要先进行预测值的绝对误差计算；然后，根据上一步结果计算出平均相对误差；最后，以计算结果检验预测模型的预测结果。所研究项目采用此法进行检验。

2)关联度检验

关联度检验首先需要进行关联度系数η(t)的计算，计算公式为

(11)

(12)

其中，当r计算结果的绝对值大于0.6，则表示预测模型检验合格，否则检验不合格。

3)后验差检验

后验差检验是以基于原始数列S1和绝对误差序列的标准S2为对象，通过S1和S2计算出相应的方差比C和小误差概率P，并将计算结果与预测精度等级表进行对应检测预测模型的预测结果，其计算公式为：

(13)

(14)

(15)

P=P{| (Δ(0)(t)-Δ(0)|<0.674 5S1}

(16)

预测精度等级如表1所示。

表1 预测精度等级

2 仿真验证与结果分析

2.1 数据样本

仿真验证在Matlab环境下进行，模型样本数据选用某省2010—2018年全社会用电总量、第一产业增长值、第二产业增长值、第三产业增长值以及居民生活用电量5项数据。以2010—2016年的样本数据作为原始数据进行模型的训练和参数的拟合计算，用2017—2108年的数据作为检验样本，验证改进MGM预测模型的预测精度和效果。表2为原始数据样本。

表2 原始数据样本

从以上历史数据中可以发现，随着经济社会的快速发展，该省全社会用电总量呈现逐年上升的趋势，同时对应的居民用电量也成相关性的逐年上升趋势。相应地可以发现，不同产业增长值与全社会用电总量的总体趋势虽然有一定的相关性，但是不同产业总会出现与全社会用电总量不相匹配的变化趋势，比如，2013年第一产业的增长值不是上升的趋势，而是出现了明显的下降，此后却又是极速的上升；第二产业在2015—2016年期间放缓了增长速率的同时，出现了一定程度的下降；第三产业在2014年的增长值不是上升的趋势，而是出现了明显的下降，此后却又迅速恢复了增长趋势：因此，充分说明随着经济社会的不断发展以及全球性的突发事件、全国性的突发事件影响，会使得各个相关因素与全社会用电总量之间的相关关系出现一定难以预测的变化。此时就要求相应的预测模型要能够在一定程度自适应这种变化，做出准确的预测，达到规划要求的预测精度，对规划工作提供更为准确的参考。

2.2 结果分析与比较

首先，运用改进的MGM预测模型进行拟合预测，预测结果如表3、表4和图3、图4所示。

表3 改进MGM预测模型化的拟合结果

表4 改进MGM预测模型的预测结果

图3 改进MGM的预测结果

图4 改进MGM的相对误差

通过改进MGM预测模型的预测结果可以看出，预测过程中出现程度不同的拟合误差，但平均相当误差在1%以内，与比较样本数据的预测误差在0.6%以内，预测结果随着模型训练不断提高，最终的预测结果精度能够满足当前社会经济环境下规划精度的要求。

利用相同的历史样本数据，采用传统MGM预测模型进行预测比较，同样在Matlab中实现，预测结果如表5、表6和图5、图6所示。

表5 传统MGM预测模型的拟合结果

表6 传统MGM预测模型的预测结果

图5 MGM的预测结果

图6 MGM的相对误差

由预测结果可见，传统MGM预测模型的平均相对误差只能控制在30%以内，相对较大。虽然随着年份的增长预测模型训练的成熟，预测结果的误差不断下降，在2018年的预测结果表现出的相对误差仅0.85%，因此该方式虽然可以进行一定程度上的预测，但预测精度难以满足现行环境下规划精度的要求。

将改进MGM预测结果与传统MGM预测结果进行比较，如图7所示。

图7 改进MGM和传统MGM的预测结果比较

由图7可得，改进MGM预测结果与真实样本数据仅存在0.6%以内的误差，表现出较为精确的预测精度，可以应对当前复杂环境下的规划预测需求。而传统的MGM预测结果，虽然随着模型的预测加深在后期表现出较为不错的预测精度，但整体而言精度偏差相对较大，不能自适应当前复杂多变的环境变化。

3 结 语

基于某省级全社会用电总量的电量数据，分析了其相关影响因素的变化发展趋势。通过分析因素间的相关关系以及当前经济社会发展需求，改进了传统MGM(1，m)模型中背景值的修正公式，使得改进的模型能够对电量数据构造出自适应的样本序列，最终实验预测某省2017—2018年的全社会用电总量的预测结果，表现出平均误差1%以内的精确结果。最后，得出以下结论：1)改进的MGM预测模型能够表现出较好的自适应能力，可以得到较高的精度结果。2)改进的MGM预测模型虽然依旧遵循灰色模型的传统规律，在中期预测结果中表现良好，但随着时间推移会呈现误差上升的趋势。然而改进后的MGM预测模型可以将误差控制在1%以内，依然能够满足规划精度的要求。