程定芳，任永建，陳正洪

(1.黃岡市氣象服務(wù)中心，湖北 黃岡 438000； 2.湖北省氣象服務(wù)中心，武漢 430074)

電力負荷受許多因素的共同影響，如經(jīng)濟發(fā)展水平、工農(nóng)業(yè)結(jié)構(gòu)、居民生活水平、電價及天氣氣候等，在極端天氣條件下天氣氣候的影響尤其不容忽視，這些因素對電力負荷的變化起著不同的影響.根據(jù)電力負荷、氣象、經(jīng)濟、社會等歷史數(shù)據(jù)，探索負荷歷史數(shù)據(jù)變化規(guī)律對未來負荷的影響，從而對未來的電力負荷進行科學(xué)的預(yù)測[1].

氣象因素是影響電力系統(tǒng)負荷預(yù)測的重要因素，研究氣象與負荷的關(guān)系對電力系統(tǒng)的負荷預(yù)測工作具有重要意義[2].因此，電力負荷與氣象條件的關(guān)系引起了眾多學(xué)者的關(guān)注，基于氣象條件的電力負荷預(yù)測做了許多卓有成效的工作.目前，國內(nèi)許多學(xué)者對城市電力負荷、用電量及其與氣象因子的關(guān)系及預(yù)測模型做了大量的研究.

張自銀等[3]利用電力負荷和氣象資料，定量分析了北京市夏季電力負荷逐日變率與主要氣象因子的關(guān)系，溫度、風(fēng)速和濕度與最大電力負荷關(guān)系顯著.廣東電力負荷具有非常明顯的線性增長趨勢，季節(jié)變化明顯，且在不同的季節(jié)與氣象因子的關(guān)系存在一定的差異[4].傅新姝和談建國[5]利用2010-2013年上海日最大電力負荷數(shù)據(jù)和同期氣象資料，分析了日最大負荷的時間變化特征及其與氣象因子的相關(guān)性.穆海振[6]基于盛夏期間日用電量與氣溫、降水、相對濕度等數(shù)據(jù)，探討了分位數(shù)回歸方法在日用電量預(yù)測中的應(yīng)用.葉殿秀等[7]利用電力負荷中分離出的氣象負荷與氣象因子采用逐步回歸方法建立預(yù)報模型，結(jié)果表明去除趨勢的氣象負荷預(yù)報模型的擬合效果更好.付曉輝等[8]利用相關(guān)分析法分析了華中電網(wǎng)氣象負荷與同期氣溫的相關(guān)關(guān)系，研究了氣象負荷隨氣溫變化的規(guī)律，并探討了相關(guān)關(guān)系的異同.江宇飛等[9]提出利用氣象指標體系進行電力負荷特性分析和預(yù)測的方法，可依據(jù)氣象條件得到不同的電力負荷預(yù)測結(jié)果.

隨著城市經(jīng)濟的迅猛發(fā)展，全球變暖背景下新能源的快速發(fā)展，使得電力負荷與眾多影響因子的關(guān)系變得更加復(fù)雜.由于電力負荷影響因素的復(fù)雜性以及負荷預(yù)測對電網(wǎng)系統(tǒng)的重要性，對電力負荷做出準確預(yù)測依然是電力和氣象工作者面臨的一項亟需解決的難題.盡管眾多學(xué)者在電力負荷與氣象因子的關(guān)系開展了大量的研究，但氣象因子對負荷預(yù)測影響的定量分析依然不足.隨著氣象部門數(shù)值天氣預(yù)報和智能網(wǎng)格預(yù)報技術(shù)的發(fā)展，天氣預(yù)報時間和空間的精細化程度有了很大的提高.目前逐15 min氣象因子預(yù)報已實現(xiàn)業(yè)務(wù)化，為氣象因子在電力負荷預(yù)測模型中廣泛應(yīng)用創(chuàng)造了很好的條件.在這一背景下，以華中電網(wǎng)某地區(qū)為研究對象，分析電力負荷變化特征及與氣象要素的關(guān)系，基于逐步回歸(SR)和神經(jīng)網(wǎng)絡(luò)(BP)方法建立預(yù)測模型，定量解析氣象因子對電力負荷預(yù)測的主要貢獻.

1主要計算方法介紹

1.1預(yù)測準確率

根據(jù)電網(wǎng)調(diào)度系統(tǒng)短期負荷預(yù)測考核細則，日預(yù)測準確率的計算依據(jù)為每日96點預(yù)測及實際數(shù)據(jù)，其中96點預(yù)測數(shù)據(jù)是指該日00:15至24：00、相鄰15 min的預(yù)測值(整點時刻值)，96點實際數(shù)據(jù)為該日00:15至24：00、相鄰15 min的實際值(整點時刻值)，具體計算公式如下：

(1)

(2)

式中,DR表示偏差率，LO觀測負荷值，LP預(yù)測負荷值;FA表示預(yù)測準確率，N表示日考核點數(shù)，一般取值96.

1.2氣象因子貢獻率

為進一步表征氣象因子在電網(wǎng)短期負荷預(yù)測中的影響，本文定義氣象因子貢獻率(meteorological contribution)，計算FF方法如下：

式中,FASR直接依據(jù)負荷數(shù)據(jù)利用逐步回歸方法得出的預(yù)測準確率，F(xiàn)ASR-M依據(jù)負荷數(shù)據(jù)和氣象數(shù)據(jù)利用逐步回歸方法得出的預(yù)測準確率，MCSR步回歸方法中氣象因子貢獻率.BP神經(jīng)網(wǎng)絡(luò)方法中氣象因子貢獻率計算方法與MCSR，標記為MCBP2相關(guān)性分析.

負荷數(shù)據(jù)選用華中電網(wǎng)某地區(qū)2018年6月每日96點電力負荷數(shù)據(jù).由于溫度對負荷的變化最為顯著[10-11]，氣象數(shù)據(jù)選用同期的溫度數(shù)據(jù)，也為每日96個值.圖1給出了6月平均負荷、溫度的日變化特征.可以看出，負荷具有明顯的雙波峰，從00∶00到5∶00左右負荷逐漸下降，然后逐漸增加，至11∶00負荷達到第一個峰值，之后負荷緩慢下降，到21∶00左右負荷達到第二個峰值，且明顯大于第一個峰值.而溫度的日變化則表現(xiàn)為單波峰特征，00∶00到6∶00溫度逐漸減小，而后逐漸增加，至14∶00左右溫度達到最大值.可見，00∶00到11∶00之間，負荷與溫度的變化存在較好的對應(yīng)關(guān)系，其他時段除受溫度的影響之外，還受到居民生活習(xí)慣、電價等因素的影響.

圖1 電力負荷和溫度的日變化曲線Fig.1 Daily variation of electric load and temperature

前人研究發(fā)現(xiàn)電力負荷除與當(dāng)日溫度相關(guān)關(guān)系顯著外，與前一日負荷的關(guān)系也比較密切[6,12].利用前一日負荷和當(dāng)日溫度建立預(yù)報模型可取的較好的預(yù)測效果[7].為此，本文計算了當(dāng)日負荷與前一日負荷、當(dāng)日溫度及前一日溫度之間的相關(guān)關(guān)系，結(jié)果如圖2所示.可以看出，當(dāng)日負荷與前一日負荷存在較好的相關(guān)關(guān)系，除6月30日外，其他當(dāng)日負荷與前一日負荷的相關(guān)均通過信度為0.001的顯著性水平，且大部分日期的相關(guān)關(guān)系在0.95以上，可見，前一日負荷對次日負荷的變化具有重要的影響.當(dāng)日溫度與負荷也存在很好的正相關(guān)，均通過了0.001的顯著性檢驗，但兩者的相關(guān)關(guān)系明顯小于前一日負荷與當(dāng)日負荷的相關(guān).由于氣象因子對負荷的影響具有累積效應(yīng)[13-14]，也計算了前一日負荷與當(dāng)日負荷之間的關(guān)系.除6月18日、6月29日和6月30日外，其他日期均為顯著的正相關(guān).

圖2 當(dāng)日負荷與前一日負荷(a)、當(dāng)日溫度(b)及前一日溫度(c)之間的相關(guān)關(guān)系，圖中實線為信度0.001的顯著性水平Fig.2 The correlation between the current load and the previous day load (a)，the current day temperature (b) and the previous day temperature (c)，The solid line in the figure is the significance level of 0.001

因此，以前一日負荷、當(dāng)日溫度、前一日溫度作為模型的輸入因子，當(dāng)日負荷作為模型的輸出因子.以逐步回歸、BP神經(jīng)網(wǎng)絡(luò)為基礎(chǔ)，采用滾動預(yù)測的方法，如利用6月1日-6月14日共14天96點負荷數(shù)據(jù)得到6月15日96點負荷數(shù)據(jù)，6月2日-6月15日共14天96點負荷數(shù)據(jù)得到6月16日96點負荷數(shù)據(jù)，依次類推，利用6月16日-6月29日的負荷數(shù)據(jù)得到6月30日的負荷數(shù)據(jù).

3預(yù)測結(jié)果分析

圖3給出利用BP和SR方法計算得到6月15日-6月30日的負荷預(yù)測準確率(FA).在BP方法中，6月16日FABP最高為93.99%，F(xiàn)ABP-M則提高至95.95%，引進氣象因子后FA提高了1.96%；6月30日FABP最低為64.68%，F(xiàn)ABP-M則提高至76.24%，引進氣象因子后FA提高了11.56%.綜合圖3(a)可見，引進氣象因子后，F(xiàn)A提高幅度在0.26%～13.41%.在SR方法中，6月15日FASR最高為96.09%，F(xiàn)ASR-M則提高至97.34%，引進氣象因子后FA提高了1.25%；6月30日FASR最低為65.89%，F(xiàn)ASR-M則提高至84.19%，引進氣象因子后FA提高了18.3%.綜合圖3(b)可見，引進氣象因子后，F(xiàn)A提高幅度在0.89%～18.3%.可以看出，不論是BP還是SR方法，氣象因子對負荷預(yù)測準確率的提高均有正的貢獻；且僅依靠歷史負荷單因子得出的FA越低的情況下，引進氣象因子后對FA改進的效果越明顯.

圖3 BP和SR方法得出的負荷預(yù)測準確率(單位：%)Fig.3 Accuracy of load forecasting obtained by BP and SR methods (unit:%)

為分析氣象因子對負荷預(yù)測的重要影響，選取BP方法中FA最高的6月16日和FA最差的6月30日，以及SR方法中FA最高的6月15日和FA最差的6月30日，分析預(yù)測負荷序列和實況負荷序列的差異.可以看出，6月15日和6月16日預(yù)測曲線與實況的對應(yīng)關(guān)系較好，對負荷上升和下降的轉(zhuǎn)折時間節(jié)點把握較好，疊加氣象因子后的預(yù)測序列更加接近實況.從影響因子的相關(guān)關(guān)系可以較好地印證(圖2)，6月15日的負荷與6月14日負荷和溫度、6月15日溫度存在較高的相關(guān)性，分別為0.98、0.81、0.82，6月16日的負荷與6月15日負荷和溫度、6月16日溫度相關(guān)性也較高，分別為0.94、0.61、0.74.6月30日而言，不論BP還是SR方法，都無法較好地預(yù)測出6∶00以后負荷序列的變化，單獨利用歷史負荷得出的預(yù)測序列與實況偏離較大，疊加氣象因子后預(yù)報偏差有所減小.從圖2中可以發(fā)現(xiàn)，6月30日負荷僅與當(dāng)日氣溫存在較好的正相關(guān)，而與6月29日負荷與溫度均呈負相關(guān)，與其他日期中負荷與影響因子的相關(guān)關(guān)系存在明顯的差異.

圖4 實況負荷與預(yù)測負荷序列的日變化曲線(單位：MW)Fig.4 Daily change of live load and predicted load series (unit: MW)

為進一步定量分析氣象因子對負荷預(yù)測的貢獻，圖5給出BP和SR方法中氣象因子對負荷預(yù)測的貢獻率.可以發(fā)現(xiàn)，氣象因子在兩種預(yù)測方法中對負荷預(yù)測準確率的提升均有正的貢獻，貢獻率分別為0.28%～17.87%和0.97%～17.78%.通過比較發(fā)現(xiàn)，利用歷史負荷序列單因子進行預(yù)測時，F(xiàn)A較高的日期，MC的貢獻率較低，6月15日，MCBP和MCSR分別為0.28%和1.29%；而利用歷史負荷序列單因子進行預(yù)測時，F(xiàn)A較低的日期，MC的貢獻率則較高，6月28和6月30日，MCBP和MCSR分別為17.54%、10.18%和17.87%、27.78%.通過上述分析可知，精細化氣象因子對短期負荷預(yù)測存在極其重大的影響，尤其是轉(zhuǎn)折天氣條件下.

圖5 BP和SR方法中氣象因子對負荷預(yù)測的貢獻率(單位：%)Fig.5 The contribution rate of meteorological factors to load forecasting in BP and SR methods (unit:%)

4小結(jié)與討論

本文定量分析了氣象因子對電力負荷預(yù)測影響，BP和SR預(yù)測模型中，氣象因子對FA的提高均有正的貢獻；且僅依靠歷史負荷單因子得出的FA越低的情況下，引進氣象因子后對FA改進的效果越明顯.氣象因子對負荷預(yù)測準確率的提升均有正的貢獻，貢獻率分別為0.28%～17.87%和0.97%～17.78%，尤其是轉(zhuǎn)折天氣條件下，精細化氣象因子對短期負荷預(yù)測的影響尤為重要.

本文在分析中僅選用了與負荷影響較大的溫度，而降水、風(fēng)速、相對濕度等要素對負荷也有一定的影響.同時需要進一步考慮工作日、節(jié)假日的影響，對夏季峰值和長假期間的變化有較好的擬合與預(yù)測[4].因此，在今后的與預(yù)測模型中需綜合考慮其他氣象因子、工作日、長假等因素的影響.