王 寶，馬燕如，楊 敏，胡光景，伍 飛，王 重

(1.國網(wǎng)安徽省電力有限公司經(jīng)濟技術(shù)研究院，安徽 合肥 230022；2.安徽明生電力投資集團有限公司，安徽 合肥 230022；3.國網(wǎng)安徽省電力有限公司，安徽 合肥 230022)

最大負荷預(yù)測是電力市場分析預(yù)測工作的核心內(nèi)容之一，是電網(wǎng)規(guī)劃與運行管理工作的重要基礎(chǔ)，對電網(wǎng)企業(yè)規(guī)劃、調(diào)控和交易等多個部門均起到重要作用[1,2]。然而最大負荷受多種因素影響，年度波動性大，大大增加了準確預(yù)測地區(qū)年度最大負荷的難度。

針對最大負荷預(yù)測的先前相關(guān)研究主要集中在日層面的短期負荷預(yù)測[3-5]，以及建筑物[6]、商業(yè)[7]等負荷預(yù)測。部分學者對最大負荷進行分解后開展預(yù)測，文獻[8]將最大負荷分解為基礎(chǔ)負荷和空調(diào)負荷，并從經(jīng)濟社會和氣象角度探索了地區(qū)空調(diào)負荷預(yù)測；文獻[9]將最大負荷分解為基礎(chǔ)負荷、非居民空調(diào)負荷和居民空調(diào)負荷，并基于經(jīng)濟社會與氣象指標開展了地區(qū)非居民和居民空調(diào)負荷預(yù)測，以上研究一方面缺乏有效的基礎(chǔ)負荷預(yù)測方法，另一方面未考慮電網(wǎng)供電能力維度因素，電網(wǎng)供電能力直接決定了最大負荷能否有效釋放。

本文結(jié)合最大負荷可分解特性，將年最大負荷分解為基礎(chǔ)負荷和空調(diào)負荷，并從經(jīng)濟增長、城鎮(zhèn)化進程、居民和服務(wù)業(yè)需求、氣象和電網(wǎng)供電能力5個維度系統(tǒng)總結(jié)基礎(chǔ)負荷和空調(diào)負荷的11個影響因素指標，為避免不同維度指標存在多重共線性而刪除部分指標的情況，采用主成分分析方法，分別構(gòu)建了基礎(chǔ)負荷和空調(diào)負荷預(yù)測模型，加總合成年度最大負荷預(yù)測模型。并以安徽省某地市為例，對預(yù)測過程和預(yù)測效果進行了分析說明。

一、最大負荷分解與多維度影響因素分析

1.最大負荷分解

最大負荷分解方法有最大負荷比較法[10,11]、基準負荷比較法[10,11]和基于溫度梯度法[2,12]等，本文采用常用的基準負荷比較法進行最大負荷分解。

4月、5月和10月氣溫總體適宜，可認為基本無空調(diào)負荷，基于4月、5月日負荷曲線計算春季基礎(chǔ)負荷曲線，基于10月日負荷曲線計算秋季基礎(chǔ)負荷曲線。假定某地區(qū)4月、5月共有工作日W1天，10月共有工作日W2天，對4月、5月各工作日相應(yīng)時點的負荷求平均得到春季基礎(chǔ)負荷曲線，對10月各工作日相應(yīng)時點的負荷求平均得到秋季基礎(chǔ)負荷曲線，如式(1)?？紤]到負荷的逐月自然增長，取春季基礎(chǔ)負荷曲線和秋季基礎(chǔ)負荷曲線平均作為夏季基礎(chǔ)負荷曲線，如式(2)。

(1)

式(1)中，PSBj和PFBj分別表示春季和秋季基礎(chǔ)負荷曲線中j時點的負荷，PSi,j和PFi,j分別表示4月、5月第i個工作日和10月第i個工作日在j時點的負荷。

(2)

式(2)中，PBj表示夏季基礎(chǔ)負荷曲線中j時點的負荷。

夏季最大負荷發(fā)生日負荷曲線與夏季基礎(chǔ)負荷曲線對應(yīng)時點之差即為夏季最大負荷發(fā)生日空調(diào)負荷曲線，空調(diào)負荷曲線上與最大負荷發(fā)生時點對應(yīng)的為最大負荷對應(yīng)的空調(diào)負荷，夏季基礎(chǔ)負荷曲線上與最大負荷發(fā)生時點對應(yīng)的為最大負荷對應(yīng)的基礎(chǔ)負荷。

2.多維度影響因素分析

基于上述最大負荷分解得出的基礎(chǔ)負荷和空調(diào)負荷，二者的影響因素差別較大，基礎(chǔ)負荷主要受經(jīng)濟社會發(fā)展方面的因素影響，而空調(diào)負荷主要受居民和服務(wù)業(yè)因素、氣象因素以及電網(wǎng)供電能力因素影響，因而有必要分別研究基礎(chǔ)負荷和空調(diào)負荷的多維度影響因素。

針對基礎(chǔ)負荷，其主要影響因素來自經(jīng)濟增長和城鎮(zhèn)化進程兩個維度。經(jīng)濟增長維度可選取GDP、工業(yè)增加值或服務(wù)業(yè)增加值三項指標中的一項，具體需要根據(jù)地區(qū)經(jīng)濟結(jié)構(gòu)的實際情況，對工業(yè)占絕對比重的地區(qū)，可選擇工業(yè)增加值指標，相反選擇服務(wù)業(yè)增加值指標；對工業(yè)和服務(wù)業(yè)比重相當?shù)牡貐^(qū)，可選擇GDP指標。城鎮(zhèn)化進程維度選擇城鎮(zhèn)化率指標。

針對空調(diào)負荷，取決于需求和供給兩方面共三個維度。在需求方面存在著兩類維度。一類維度是居民和服務(wù)業(yè)需求維度，主要包括人口、居民收入、空調(diào)擁有量、商業(yè)使用面積等指標。其中居民收入可采用城鄉(xiāng)居民人均可支配收入指標，城鄉(xiāng)居民人均可支配收入=城鎮(zhèn)居民人均可支配收入×城鎮(zhèn)化率+農(nóng)村居民人均可支配收入×(1-城鎮(zhèn)化率)，商業(yè)使用面積指標為歷史各年(商品房銷售面積-住宅銷售面積)的累計值，考慮銷售與投入使用存在一定的時間差，商業(yè)使用面積可按上年及之前歷史年份(商品房銷售面積-住宅銷售面積)的累計值計算。這一類維度決定了空調(diào)負荷正常的增長水平。另一類維度是氣象維度，主要包括最大負荷發(fā)生日最高氣溫和高溫累積效應(yīng)(可考慮最大負荷發(fā)生日前35 ℃以上高溫持續(xù)天數(shù))。該維度決定了正常的負荷增長能否釋放以及釋放程度，即達不到或超過正常的負荷增長水平。在供給方面主要是電網(wǎng)供電能力維度，由居民和服務(wù)業(yè)需求維度和氣象維度決定的空調(diào)負荷需求能否最終釋放取決于電網(wǎng)供電能力，該維度指標包括地區(qū)、城、農(nóng)網(wǎng)戶均配變?nèi)萘恐笜恕５貐^(qū)最大負荷多維度影響因素與指標如表1所示。

表1 地區(qū)最大負荷多維度影響因素與指標表

二、基于多維度與主成分分析的最大負荷分解預(yù)測模型

考慮到基礎(chǔ)負荷與空調(diào)負荷主要影響維度因素不同，因而有必要分別開展基礎(chǔ)負荷和空調(diào)負荷建模預(yù)測，最終加總合成最大負荷。由于影響空調(diào)負荷的維度因素較多，為避免直接建立回歸方程而存在多重共線性，本文采用主成分分析方法[13,14]。

1.基礎(chǔ)負荷預(yù)測模型

基礎(chǔ)負荷主要受經(jīng)濟增長和城鎮(zhèn)化進程兩個維度影響，可建立以基礎(chǔ)負荷為被解釋變量，經(jīng)濟增長和城鎮(zhèn)化率兩項指標作為解釋變量的多元線性回歸方程，經(jīng)濟增長指標需要進行可比價處理。

(3)

2.空調(diào)負荷預(yù)測模型

空調(diào)負荷受三個維度因素影響，可能的影響指標合計約9個，首先需要對地區(qū)空調(diào)負荷與各維度中的指標進行相關(guān)性分析，計算相關(guān)系數(shù)，如式(4)。在各維度中選擇相關(guān)性最強的關(guān)鍵指標，作為解釋變量，并形成解釋變量集X=[JFQXDW]，JF、QX和DW分別表示選取的居民和服務(wù)業(yè)需求維度指標、氣象維度指標以及電網(wǎng)供電能力維度指標。

(4)

居民和服務(wù)業(yè)需求、氣象以及電網(wǎng)供電能力三個維度是從不同角度解釋地區(qū)空調(diào)負荷增長的原因，但這三個維度中選取的指標之間存在相關(guān)性高的情況容易出現(xiàn)，從而導致直接基于這些指標建立多元線性回歸方程容易出現(xiàn)多重共線性，不得不剔除部分指標，從而降低了模型的實際意義。

本文采用主成分分析方法，先對選取的三個維度中的指標進行主成分提取，提取1～2個能夠高度反映原始指標的主成分，再以主成分作綜合解釋變量、以地區(qū)空調(diào)負荷為被解釋變量建立多元或一元線性回歸方程，既可消除直接建立回歸方程產(chǎn)生的多重共線性，還可以確定各原始指標對地區(qū)空調(diào)負荷影響程度大小。

首先對解釋變量集X中各指標進行如式(5)的標準化處理，然后計算xi序列形成的變量集x的相關(guān)系數(shù)矩陣R。

(5)

利用Matlab求得R的特征值Ri及其對應(yīng)的特征向量，如式(6)計算各特征值的貢獻度λi，確定貢獻度最大的1～2個特征值對應(yīng)的特征向量與變量集x各指標對應(yīng)相乘加總得到1～2個主成分。如：貢獻度最大的特征值對應(yīng)的特征向量為β=[β1,β2,…,βN]，則該特征值對應(yīng)的主成分yt=β1x1t+β2x2t+…+βNxNt。

(6)

式(6)中，N表示特征值數(shù)量，與變量集X和x指標個數(shù)相同。

以選取一個主成分yt為例，以該主成分為綜合解釋變量、以地區(qū)空調(diào)負荷為被解釋變量，建立回歸方程，并通過回推，可確定原始解釋變量集X中各指標對地區(qū)空調(diào)負荷影響大小，如式(7)。

(7)

3.最大負荷預(yù)測模型

將基礎(chǔ)負荷和空調(diào)負荷預(yù)測方程加總，合成形成地區(qū)最大負荷Ft預(yù)測方程，如式(8)。

Ft=FJt+FKt

(8)

三、實例分析

以安徽某地市為例，具體數(shù)據(jù)如表2所示，對上述分解預(yù)測模型進行闡述。

表2 某地區(qū)基礎(chǔ)負荷與空調(diào)負荷表 單位:萬千瓦

基于上述最大負荷分解方法，得到該地區(qū)歷史基礎(chǔ)負荷和最大負荷，隨著經(jīng)濟步入發(fā)展新常態(tài)，該地區(qū)基礎(chǔ)負荷增長呈現(xiàn)趨緩態(tài)勢，但空調(diào)負荷快速增長。

針對基礎(chǔ)負荷預(yù)測，該地區(qū)二產(chǎn)和三產(chǎn)經(jīng)濟比重較為接近，經(jīng)濟增長指標選擇GDP指標，城鎮(zhèn)化進程選擇常住人口城鎮(zhèn)化率指標。該地區(qū)基礎(chǔ)負荷與GDP和常住人口城鎮(zhèn)化率間相關(guān)系數(shù)分別高達0.998和0.999，圖1顯示了該地區(qū)空調(diào)負荷與常住人口城鎮(zhèn)化率間散點圖。

圖1 該地區(qū)基礎(chǔ)負荷與城鎮(zhèn)化率間散點圖

以GDP和常住人口城鎮(zhèn)化率為解釋變量、基礎(chǔ)負荷為被解釋變量構(gòu)建回歸方程，如式(9)。

FJt=-94.069+0.017JZt+508.929CHt

(9)

針對空調(diào)負荷預(yù)測，考慮統(tǒng)計局公布年度空調(diào)擁有量數(shù)據(jù)過于滯后，且空調(diào)擁有量和居民收入及人口相關(guān)性大，以及房地產(chǎn)銷售歷史數(shù)據(jù)無法全面獲取，且該地區(qū)居民用電比重明顯偏高，服務(wù)業(yè)對空調(diào)負荷影響明顯弱于居民，因而上述提出的空調(diào)負荷三個影響維度共9個指標簡化為7個。

表3給出了該地區(qū)空調(diào)負荷與三個維度中7個指標間的相關(guān)系數(shù)，相關(guān)性均較高，該地區(qū)空調(diào)負荷與農(nóng)網(wǎng)戶均配變?nèi)萘块g相關(guān)系數(shù)高于全市和城網(wǎng)，這與該地區(qū)農(nóng)村人口多、城鎮(zhèn)化率偏低，這幾年農(nóng)網(wǎng)升級改造力度大等實際情況是吻合的，因而電網(wǎng)供電能力維度選取農(nóng)網(wǎng)戶均配變?nèi)萘恐笜?。再加上居民和服?wù)業(yè)需求維度2個指標和氣象維度2個指標，共形成5個指標作為解釋變量。

表3 該地區(qū)空調(diào)負荷與三個維度指標間相關(guān)系數(shù)表

直接建立多元線性回歸方程，如式(10)發(fā)現(xiàn)，由于多重共線性，多項解釋變量前面的系數(shù)為負，經(jīng)濟意義不合理。

FKt=-82.352-0.001 29JF1t-1.344JF2t+29.856QX1t-5.342QX2t+122.302DW1t

(10)

式(10)中，JF1和JF2分別表示城鄉(xiāng)居民人均可支配收入和常住人口，QX1和QX2分別表示最大負荷發(fā)生日最高溫度和最大負荷發(fā)生前35 ℃以上高溫持續(xù)天數(shù)，DW1表示農(nóng)網(wǎng)戶均配變?nèi)萘俊?/p>

以下采用主成分分析方法，從5個原始指標中提取主成分。首先對各解釋變量進行標準化處理，經(jīng)Eviews軟件計算標準化后解釋變量的相關(guān)系數(shù)矩陣，再經(jīng)Matlab計算相關(guān)系數(shù)矩陣特征值及其特征向量，如表4所示。

表4 5個特征值及其對應(yīng)的特征向量表

第一特征值的貢獻率高達83.4%，能夠很好地反映各解釋變量信息，因而以該特征值及其對應(yīng)的特征向量形成的主成分yt=0.478x1t+0.469x2t+0.372x3t+0.440x4t+0.468x5t作為綜合解釋變量，圖2顯示了該地區(qū)空調(diào)負荷與綜合解釋變量散點圖，二者相關(guān)系數(shù)在0.99以上。

圖2 該地區(qū)空調(diào)負荷與綜合解釋變量(主成分)散點圖

以綜合解釋變量為解釋變量，以空調(diào)負荷作為被解釋變量，建立回歸方程，并進行回推還原得到地區(qū)空調(diào)負荷與三個維度5項原始指標間關(guān)系方程，如式(11)。

FKt=90.166+20.845yt

=90.166+20845×(0.478x1t+0.469x2t+0.372x3t+0.440x4t+0.468x5t)

=-838.45+0.0045JF1t+0.555JF2t+10.867Qx1t+3.213QX2t+18.872DW1t

(11)

從上述方程可看出，各指標對空調(diào)負荷均產(chǎn)生正向影響，與相關(guān)系數(shù)方向是一致的。在其他指標不變情況下，城鄉(xiāng)居民人均可支配收入每增加1 000元，空調(diào)負荷平均增加4.5萬千瓦；常住人口每增加10萬人，空調(diào)負荷平均增加5.5萬千瓦；最大負荷發(fā)生日最高溫度每增加1 ℃，空調(diào)負荷平均增加10.9萬千瓦；最大負荷發(fā)生前35 ℃以上高溫持續(xù)天數(shù)每增加1天，空調(diào)負荷平均增加3.2萬千瓦；農(nóng)網(wǎng)戶均配變?nèi)萘棵吭黾?.1千伏安/戶，空調(diào)負荷平均增加1.9萬千瓦。

將基礎(chǔ)負荷預(yù)測方程和空調(diào)負荷預(yù)測方程加總合成地區(qū)最大負荷預(yù)測方程，如式(12)。

Ft=FJt+FKt

=-93.252+0.017JZt+508.929CHt+0.0045JF1t+0.555JF2t+10.867QX1t+3.213QX2t+18.872DW1t

(12)

以影響基礎(chǔ)負荷和空調(diào)負荷的多維度指標2014～2018年實際值帶入模型，得到2014～2018年該地區(qū)基礎(chǔ)負荷、空調(diào)負荷和最大負荷回測值，如表5所示。經(jīng)與實際值比較，最大負荷回測的平均絕對百分誤差率在1%以內(nèi)。

表5 現(xiàn)模型回測效果表

表6給出了本文提出的最大負荷分解預(yù)測模型與基于主成分分析的最大負荷預(yù)測模型(不分解為基礎(chǔ)負荷和空調(diào)負荷)的預(yù)測效果對比，由于本文有針對性分別研究了基礎(chǔ)負荷和空調(diào)負荷多維度影響因素與建模，本文提出的分解預(yù)測方法預(yù)測效果明顯好于不分解方法。

表6 模型回測效果對比表

四、結(jié)論

最大負荷預(yù)測是電網(wǎng)規(guī)劃運行的重要基礎(chǔ)性工作，最大負荷預(yù)測準確與否直接關(guān)系到電網(wǎng)規(guī)劃工作能否有效滿足電網(wǎng)實際運行需要。針對最大負荷波動大、影響因素多樣，且本身具有可分解特性，本文從將最大負荷分解為基礎(chǔ)負荷和空調(diào)負荷的角度，系統(tǒng)總結(jié)了經(jīng)濟增長、城鎮(zhèn)化進程、居民和服務(wù)業(yè)需求、氣象和電網(wǎng)供電能力5個維度共11個最大負荷主要影響因素指標，考慮到各維度影響因素指標間存在多重共線性，構(gòu)建了基于多維度與主成分分析的年最大負荷分解預(yù)測模型，將分別構(gòu)建的基礎(chǔ)負荷和空調(diào)負荷預(yù)測方程加總合成形成最大負荷預(yù)測方程。實例表明本文提出的最大負荷分解預(yù)測模型能夠取得良好的預(yù)測效果，為電力市場分析預(yù)測和電網(wǎng)規(guī)劃人員提供了一種有效的年度最大負荷預(yù)測方法。