一种基于因子-主属性模型的中长期电力负荷预测方法
【专利摘要】本发明公开了一种基于因子-主属性模型的中长期电力负荷预测方法，所述方法包括：步骤1：建立(n-1)个影响因子指标X1-Xn-1和1个预测对象Xn的原始矩阵，采用Z标准化，对影响因子指标数据X1-Xn-1和Xn进行预处理，将指标数据无量纲化，得到指标矩阵Aoxn；步骤2：根据因子分析方法，确定选择的公因子，并计算相应的因子得分，建立因子预测模型步骤3：根据主属性算法，筛选得到主属性m1,…,mr，主属性模型为其中v是特征向量；步骤4：根据公式建立因子-主属性的中长期电力负荷预测模型，计算得到归一化后的预测变量值，经过Z标准化公式变换，重新计算得到相应的原变量实际值。
【专利说明】—种基于因子-主属性模型的中长期电力负荷预测方法
【技术领域】
[0001]本发明属于电力负荷预测方法研究领域，涉及一种受多因子影响的复杂问题预测方法，尤其涉及一种基于因子-主属性模型的中长期电力负荷预测方法。
【背景技术】
[0002]预测学是通过数学模型的研究预测并控制事物发展的不确定性，找到事物发展规律的方法和理论研究的一门学科。其应用广泛，在各个行业和领域占据着重要的地位，特别是在工业现代化发展时期，对于合理规划和决策发挥了很大作用。
[0003]预测学的理论和方法众多，可以从定性和定量的角度进行分类。定性预测方法主观性强，主要凭借的是人的经验和分析能力，适用于缺乏历史统计资料或者研究对象的趋势变化较大的情况。定量预测方法对统计资料的要求和依赖性高，计算量大，是从数量上分析把握客观事物发展变化趋势。这两类方法分别有很多相应的预测方法，而通常都是来自这两类方法中的单一方法被单独应用到各个预测问题上，对研究对象的研究要么是偏重性质的分析，要么是偏重数据结构的体现，缺乏相互结合。
[0004]随着经济、社会和科技的快速发展，一般待预测对象的背景越来越复杂，常常受诸多因素影响，并且这些影响来源甚广，影响的范围也不同，所以传统的通过单个变量进行拟合预测的方法已经不适合现在日益多元化的预测问题。
[0005]传统的电力负荷中长期预测一般采用的是单一的经济预测方法，研究数据局限于研究对象本身单一的数据，或者仅仅涉及单一的影响因素。随着深入研究，发现电力负荷预测值受多方面影响，而由中长期负荷的特点，负荷的中长期预测与经济、政策以及负荷自身发展规律密切相关，最近研究发现，还与气候因素有关。面对现代电力发展提出的准确性和可解释性的预测新要求，中长期负荷预测需要适应这种要求。
[0006]针对上述问题，本发明提出一种基于因子-主属性模型的中长期电力负荷预测方法。该方法分别从影响因素和数据特点出发，根据预测对象的指标矩阵和预测对象电力负荷值矩阵数值关系构建因子-主属性预测模型。应用该模型可以实现对有多因素影响的电力负荷做出合理准确预测，为电网合理规划、调度以及供电企业的基建和供电等决策提供参考。
【发明内容】
[0007]本发明的目的在于，提供一种基于因子-主属性模型的中长期电力负荷预测方法，对于受多因素影响的电力负荷，兼顾外在影响因素和数据变化特征，做出准确的预测，并且对预测结果进行很好的解释，而且能提高预测的准确度。
[0008]为了实现上述目的，本发明提出的技术方案是，一种基于因子-主属性模型的中长期电力负荷预测方法，其特征是所述方法包括下列步骤:
[0009]步骤1:建立(η-1)个影响因子指标X1-Xlri和I个预测对象Xn的原始矩阵，采用Z标准化，对影响因子指标数据X1-Xlri和Xn进行预处理，将指标数据无量纲化，得到指标矩阵Atjxn,矩阵中的ο代表数据的时间维度上的个数。其中，对指标数据进行Z标准化是为了将不同单位的指标数据无量纲化，使得这些指标数据具有相对可比性。Z标准化公式为
χ.? 1|.異其中μ i为变量的期望，σ为标准差。所述η-1个影响因子指标X1-Xlri根据
Cr
预测对象的特点进行选择，通常是根据专家法进行选择，所述预测对象Xn为全社会年用电量。
[0010]步骤2:根据因子分析方法，确定选择的公因子，并计算相应的因子得分，最后建立因子预测模型X = 4 X /? + In XF2+-*- + Iik X Fk =
[0011]步骤3:根据主属性算法，筛选得到主属性HI1，…，πν，得到主属性模型为X ---(/IilH-----1-Wf),其中V是特征向量。
V
I
[0012]步骤4:根据公式11 = 0(/.,6+…+ +…+i?r),建立因子-主属性的中
V
长期电力负荷预测模型，计算得到归一化后的预测变量值，经过Z标准化公式变换，重新计算得到相应的原变量实际值。其中α与β是因子预测模型和主属性模型算法的载荷因子。
[0013]所述步骤2包括如下子步骤:
[0014]子步骤201:建立原方程组X = AY,其中A即为步骤I中得到的指标矩阵Atjxn，计算X的协方差矩阵的特征值，将矩阵A各列按特征值大小排序；计算各个变量在总方差中的比例并进行排列，按照需要的方差贡献率大小选择相应的公因子个数，如大于90%的方差贡献率时，排列在90%方差贡献率之前的变量个数k，即为相应的公因子个数k。这里的公因子不是直接由指标矩阵得到的某一个指标，而是抽象出来的一个概念。
[0015]令f为公因子，爲=丨私,Og = ^jT1 ? Xi为矩阵A的特征值，a*,,为步骤I中相应A矩阵的元素；运用主成分分析法建立因子分析方程组
A = aufi + anfz + aJ'i + …+ eufi+SΛ\ = /| + U、, /-, 4- + -" + Un f, + ￡
_一—，其中k为步骤I得到的公因子个数，ei(i = l，…，..*'? = + --:./: + +*'■+./* + S
12)为误差，可以忽略，求解因子载荷，根据主成分分析法原理计算，得到相应k个主成分的
卜丨?12…aU
QQ“》 Q
因子载荷矩阵Φ，即.=::: ?
Lg ?1 …h
[0016]子步骤202:将因子载荷矩阵Φ进行旋转，得到因子载荷矩阵L。
[0017]其中，因子旋转方式分为正交旋转和斜交旋转两种。通常选用正交旋转方式，即坐标轴在旋转过程中始终保持垂直，新生成的因子可保持不相关；采用最大方差法进行正交旋转，得到旋转后在k个主成分上各个指标因子相应的正交旋转因子载荷，
[0018]经过正交旋转方式旋转的因子载荷记为IijQ = I, 2，…，P ； j = I, 2，…，k),因子载荷矩阵即为L ;
[0019]子步骤203:根据子步骤202中得到的正交旋转方式旋转后因子载荷矩阵L和经过Z标准化的原始变量，采用未加权的最小二乘法计算因子得分系数，如式
F = (VL)其中η为得分系数，X为Z标准化的原始指标变量，L’为L的转置矩阵。
[0020]所述步骤3包括如下子步骤:
[0021]子步骤301:将预测对象的指标矩阵Atjxn进行奇异值分解SVD (Singular ValueDecomposit1n)，得至丨J
[0022]Aoxn = Uoxo ΣοΧηντηΧη
[0023]其中矩阵Utjxtj中列向量相互正交，Stjxn是对角矩阵，VTnXn中列向量相互正交。根据SVD得到特征值为λ i，具体操作可以应用matlab中的奇异值分解函数实现，特征向量为
Vi, { Vi, I ^ i ^ η}构成矩阵丫以…|?i =^,其中的σ ^是奇异值，μ j构成矩阵Um。。对角矩阵1。&是全体奇异值从大到小的排列；
[0024]子步骤302:将排列在前面的r大的奇异值近似地描述矩阵为:
?UoxXxrVlI，经过变换可以得到A-V衝^UkoxXxr=Ahmr?所述Aht^是H的一个简洁替代表示，表示相乘后得到的一个矩阵；所述r表示的是奇异值按大小排序。
[0025]其中，一般的，最大的排列在前10 %的特征值占据全体特征值的99 %。所以，可以使用最大的r个特征值近似地描述矩阵Aoxn ；
[0026]子步骤303:根据奇异值分解后的AtjxnVn^矩阵计算主属性量化值λ i，根据λ i进行排序，值最大的r列即是选择的最为重要的r列，即可得到r个主属性，记作Hi1,…，；
[0027]其中主属性指的是对预测变量的值有重要影响的主要属性；主属性量化值指的是对主属性衡量的一种计量。
[0028]子步骤304:建立主属性模型I = !(% +…+W,)，其中V为预测变量在子步骤302
V
得到的矩阵\xr中的相应值。
[0029]本发明在因子分析的基础上，构建了因子-主属性预测模型；该模型不但可以描述预测对象的外在因素和内在结构因素的影响，同时，本发明对于受多因素影响的对象的预测具有较高的准确性和可解释性，能够对基于该预测对象的其他问题做出更好的决策。
【专利附图】
【附图说明】
[0030]图1是基于因子-主属性模型的预测方法流程图；
【具体实施方式】
[0031]下面结合附图，对本预测方法的实施例作详细说明。应该强调的是，下述说明仅仅是示例性的，而不是为了限制本发明的范围及其应用。
[0032]本申请使用基于因子-主属性模型的预测方法对中长期电力负荷的实例进行预测，所选用的实例对象中长期电力负荷用年负荷量(单位:亿千瓦时)表示，使用的指标数据可以在中华人民共和国国家统计局查到。另外，本申请从经济和消费角度选择对中长期电力负荷产生影响的6个方面相应指标，实际应用情况不限于这6个指标，往往更加复杂，这里仅仅是示例。
[0033]步骤1:选择国内生产总值(X1)、第一产业生产总值(X2)、第二产业生产总值(X3)、第三产业生产总值(X4)、能源消费增长率(X5)、电力消费增长率(X6)作为中长期电力负荷的6个影响因子指标，采用全社会年用电量(X7)表示中长期电力负荷，建立相应的指标体系；其中，这些指标中国内生产总值X1、第一产业生产总值X2、第二产业生产总值X3、第三产业生产总值X4的单位都是亿元，能源消费增长率X5是根据能源消费量(单位:万吨标准煤)计算增长得到的百分比，电力消费增长率X6是根据电力消费(单位:亿千瓦时)计算增长得到的百分比。本实例采集2001年-2011年的相应指标数据进行计算，建立6个指标和I
?09655,17 15781.3 49512,3 44361.6 3 3 9.3 }4682,5Il
120332.69 J 6537 53896.8 49898.9 6 11M 16385.82个预测对象的原始矩阵為=.-■■ - ■ - N
472881.56 47486.2 220412,8 205205 7.112.1 47022 J
一 、X — μ,、一
[0034]通过Z标准化公式Zij其中为变量的期望，σ为标
CF
准差，对指标数据进行标准化预处理，将指标数据无量纲化，得到指标矩阵:U3!96 -1.12120 -1.15368 -1.10911-1.17365-0.82570-1.334211
' _ -1.04425 >1.04925 -1,0771! -1.00673-0.54218-0.07877-1.17394 J > 7..?, * * * I O
*..■*N?麵
*.,**?I*I
1.85160 1.89703 !.83077 1.86476 -0,28491 0.01086 1.70864 J
[0035]步骤2:根据因子分析方法，确定选择的公因子，并计算相应的因子得分，建立因子预测模型X=Iil X忒+ /,2 X為+…+4 X 。
[0036]这一过程包括如下子步骤:
[0037]子步骤201:构造原方程组X = AY,其中A即为步骤I中得到的指标矩阵Atjxn，计算X的协方差矩阵的特征值，将矩阵A各列按特征值大小排序；计算各个变量在总方差中的比例并进行排列，按照需要的方差贡献率大小进行选择相应的公因子个数，如大于90%的方差贡献率时，排列在90%方差贡献率之前的变量个数k，即为相应的公因子个数k。这里的公因子不是直接由指标矩阵得到的某一个指标，而是抽象出来的一个概念。
[0038]令f为公因子，fi = fi / E，e# = ^ciii ? Ai为矩阵A的特征值，
为步骤I中相应A矩阵的元素；运用主成分分析法构造因子分析的方程组
X1 = U1 J\ + QvJ1 + ai:J- +-- + OlJk +￡
Xj =a,,+-- f\+a^f\+...+ an f, + ε
' 二一,其中Mi = 1，…，11)为误差，可以忽略；可以求
Λ = +11,,?/?.+"-+(InJl +G
解各个指标在k个主成分上的因子载荷矩阵Φ。本例通过主成分法计算，应用特征值等于I为条件，计算得到主成分个数k = 2。根据主成分分析法原理计算，得到相应这两个主成
-0.991 0.132""
0.989 0,136
0.990 0.137
分上的因子载荷矩阵Φ= 0.991 0.127 ;具体计算可以通过SPSS软件实现。可以看到各
-0.473 0.818
-0.387 0.867
—0.975 0.201 —
个指标在这两个主成分上的数值分布还不够明显，需要进行因子载荷的旋转。
[0039]子步骤202:将因子载荷矩阵Φ进行旋转。
[0040]其中，旋转方式分为正交旋转和斜交旋转两种。通常选用正交旋转方式，即坐标轴在旋转过程中始终保持垂直，新生成的因子可保持不相关；
[0041]本实例采用最大方差法进行正交旋转，得到旋转后在2个主成分上各个指标因子
'0.984 -0.176'
0.983 -0.171
0.985 -0.171
相应的正交旋转因子载荷，此时的正交旋转因子载荷矩阵i= 0.982 -0.181 ο
-0,202 0,923
-0.105 0.943
_ 0,990 4,105—
[0042]子步骤203:根据上一步骤中得到的正交旋转方式旋转后因子载荷矩阵L和经过Z标准化的原始变量，采用未加权的最小二乘法计算因子得分系数，如式F = = ,其中η为得分系数，X为Z标准化的原始指标变量，L'表示L的转置矩阵。本例通过SPSS软件进行相应的未加权的最小二乘法计算。
[0043]计算得到Ii1= (0.206, 0.207, 0.207, 0.205, 0.077, 0.102, 0.217)，n2 = (0.025，0.027，0.028，0.021，0.536，0.562，0.068)，则相应的
'f； =0.206X, + 0.207X, +0.207X, +0.205X4 +0,077J5 +0,102^+0.217^_
F2 =0.025^,+0.027,r, +0,028^,+0.02IX4 +0.536^, +0.562^ +0.068’ 1 1
?~1.39520-L29H79
-1.16594-0.52118
-0.7351H13/092
-0.473 fi51.42070
?0.4 hm0.44964
1,-,7)是标准化后的矩阵A12x7中的相应数值，则得到# =| -0.134470.55579 0
10.264560.44063
10.29304-1.54703
0,58048-1.08096
I 1,299610.11177
[1.885250,15851
[0044]为了简便，本例以2011年的因子得分作为2012年因子得分
Pi和A的值，根据因子预测模型X =ιΛ4'七in>a\ + -"+iik戍得至ij
Xf12 = Iu XF1+ Iia XF2= 0,99χ 1.88525 + (-0,105) x 0.15851 = 1.84975395 ο
[0045]步骤3:根据主属性算法，筛选得到主属性Hi1，…，πν，得到主属性模型为
I
M = —(mx+ ….￥mr) ο
V
[0046]这一过程包括如下子步骤:
[0047]子步骤301:将预测对象的指标矩阵Anx7进行奇异值分解SVD (Singular ValueDecomposit1n)，得到
[0048]Alix7 — Ulixil Σ 11X7V 7X7
[0049]其中矩阵Unxn中列向量相互正交，矩阵Σ11Χ7是对角矩阵，矩阵Vt7x7是V7x7的转置，即特征向量矩阵的转置矩阵，其中列向量相互正交。根据SVD得到特征值为Ai，特征向量为V J, { V J, I ^ i ^ 11}构成矩阵νηΧη，|σ；=+Α%1，其中的σ J是奇异值，μ i
构成矩阵Um。。对角矩阵2。&是全体奇异值从大到小的排列；
[0050]计算得到
"0,41970,14730,0046-0.1568-0.20560.17750,8389 —
0.41910.1455-0.0646-0,39930.7615-0.2329-0,0736
0.41890.15580.02960.11950.0127().800()-0.3S0N
V.7 = 0.41990.1390-0,0071-0.3937-0,6 Π 7-0.3598-03817/Λ' ，
-0.24840.6655-0.70240.0142-0,03500.02480,0000
-0.25510.65370.7063-0.08920.0265-0.00690.0002
^ 0.41200.20790.05210.79910.0405-0.3797-0.0033 _
【权利要求】
1.一种基于因子-主属性模型的中长期电力负荷预测方法，其特征是所述方法包括: 步骤1:建立(η-1)个影响因子指标X1-Xlri和I个预测对象Xn的原始矩阵，采用Z标准化，对影响因子指标数据X1-Xlri和Xn进行预处理，将指标数据无量纲化，得到指标矩阵A_，矩阵中的ο代表数据的时间维度上的个数；其中，对指标数据进行Z标准化以便将不同单位的指标数据无量纲化；Z标准化公式为Zi =?，其中UiS变量的期望，σ为标准差；所σ述η-1个影响因子指标X1-Xlri根据预测对象的特点按照专家法进行选择，所述预测对象Xn为全社会年用电量；步骤2:根据因子分析方法，确定选择的公因子，并计算相应的因子得分，建立因子预测模型I = 4 X /;+/,: X 為 + …+ Iik X Fk;步骤3:根据主属性算法，筛选得到主属性HI1，…，!V主属性模型为I = !(叫+...+ %)?
V其中V是特征向量；步骤4:根据公式P = <1(4#； +…+ //；) +於丄化+…+ mr),建立因子-主属性的中长期F电力负荷预测模型，计算得到归一化后的预测变量值，经过Z标准化公式变换，重新计算得到相应的原变量实际值；其中α与β是因子预测模型和主属性模型算法的载荷因子。
2.根据权利要求1所述的方法，其特征是所述步骤2包括如下子步骤:子步骤201:建立原方程组X = ΑΥ，其中A为步骤I中得到的指标矩阵Α_，运用主成
X1 =^rllZl +^2/2 + ^13/3 + …+ 沒
x7 =^1/,* +仏7/; +α71/; + “*+β，￡./;.+ε分分析法建立因子分析方程组」.a 政她,其中ε = 1，…，
Λ =αΛ+ aUtfi + a,J'i +***+itUtfk + e12)为误差可以忽略，f为公因子，k为步骤I得到的公因子个数，求解因子载荷，根据主成^ll aU …?I*'
龜Q■ , ? Qr分分析法原理计算，得到相应k个主成分的因子载荷矩阵Φ，即Φ=.21.22f ；
αΛ …aI 子步骤202:将因子载荷矩阵φ进行旋转，得到因子载荷矩阵L ； 采用最大方差法进行正交旋转，得到旋转后在k个主成分上各个指标因子相应的正交旋转因子载荷，正交旋转因子载荷矩阵即为L ；子步骤203:根据步骤202中得到的正交旋转方式旋转后因子载荷矩阵L和经过Z标准化的原始变量，采用未加权的最小二乘法计算因子得分系数，# = (IZI)H = 其中η为得分系数，X为Z标准化的原始指标变量，L'为L的转置矩阵。
3.根据权利要求1所述的方法，其特征是所述步骤3包括如下子步骤: 子步骤301:将预测对象的指标矩阵Atjxn进行奇异值分解，得到Atjxn = υοΧοΣοΧηντηΧη ； 其中矩阵Um。中列向量相互正交，Stjxn是对角矩阵，VTnXn是Vn5ta的转置，其中列向量相互正交；根据SVD得到特征值Ai，特征向量Vi，(Vi, I ^ i ^nj构成矩阵乂_，fI]1i =νλΓ*μ? =^Avi ，其中的σ i是奇异值，μ 1构成矩阵U#。;对角矩阵Stjxn是全体奇异值从大到小的排列；子步骤302:将前!■大的奇异值近似地描述矩阵为:A卿?经过变换可以得至IJH.,.-1d 二 Ah^f;所述r表示的是奇异值按大小排序，r = 7； 子步骤303:根据奇异值分解后的AtjxnVn^矩阵计算主属性量化值Xi，根据^进行排序，选择值最大的两个Xi，其对应的原始变量为国内生产总值Xl和电力消费增长率X6，值最大的r列即是选择的最为重要的r列，即可得到r个主属性，记作Hi1,…，mr ；子步骤304:建立主属性模型X = ^iml +…+ %)，其中v为预测变量在子步骤302得到
V的矩阵VnA中的相应值。
【文档编号】G06Q50/06GK104200283SQ201410428808
【公开日】2014年12月10日 申请日期:2014年8月27日 优先权日:2014年8月27日
【发明者】李国栋, 刘琳, 黄琳华, 李凯, 宋志新, 李小龙 申请人:华北电力大学, 国家电网公司, 国网新疆电力公司信息通信公司