王瑞慶，王宏福

（安陽師范學(xué)院計算機與信息工程學(xué)院，安陽 455000）

電價序列的高階矩波動特征

王瑞慶，王宏福

（安陽師范學(xué)院計算機與信息工程學(xué)院，安陽 455000）

電價的分布特性是電力市場風(fēng)險管理和電力金融產(chǎn)品定價的重要依據(jù)。在對電力市場現(xiàn)貨電價的變動規(guī)律綜合分析的基礎(chǔ)上，使用條件方差、條件偏度、條件峰度來刻畫電價序列的二階矩、三階矩、四階矩時變特征，使用正弦函數(shù)和負(fù)荷平方來刻畫電價序列的多重周期及其與負(fù)荷之間的相關(guān)性，建立了一個基于正態(tài)分布概率密度函數(shù)的Gram-Charlier展開的多周期GARCH-M模型。對PJM電力市場歷史數(shù)據(jù)的分析表明：系統(tǒng)負(fù)荷對電價均值具有顯著的影響，電價序列具有周、半月、月、雙月、季、半年等多重周期和二階矩波動集聚性，二階矩、三階矩和四階矩具有明顯的時變特征，三階矩和四階矩的尖峰跳躍具有同步性。

電力價格；電力負(fù)荷；高階矩；波動集聚；廣義條件異方差均值模型

電價是電力市場供求平衡時形成的出清價，不僅受氣候、系統(tǒng)負(fù)荷、發(fā)電成本、可用發(fā)電容量、輸電網(wǎng)絡(luò)阻塞等客觀因素的影響，還受到市場交易規(guī)則、參與者的競價策略及其對價格的心理反應(yīng)等主觀因素的影響，這些因素使得準(zhǔn)確的電價預(yù)測較為困難。當(dāng)前預(yù)測方法主要包括通過模擬電力市場競爭規(guī)則來預(yù)測市場出清電價的長期預(yù)測方法和依據(jù)大量歷史數(shù)據(jù)建立反映電價變化規(guī)律數(shù)學(xué)模型的短期預(yù)測方法[1]。

神經(jīng)網(wǎng)絡(luò)對非確定性、非精確性規(guī)律具有自適應(yīng)能力，能有效處理多變量和非線性問題[2-10]。但由于神經(jīng)網(wǎng)絡(luò)方法的參數(shù)調(diào)整不夠靈活，學(xué)習(xí)速度較慢，在實際應(yīng)用中遇到了困難。

時間序列方法需要的歷史數(shù)據(jù)相對較少，能準(zhǔn)確反映歷史電價變化的連續(xù)性，常用的有自回歸滑動平均（auto-regressive and moving average）和帶外生變量的ARMA（ARMA with exogenous variables，ARMAX）模型。文獻(xiàn)[11]建立了一個以負(fù)荷作為外生解釋變量的ARMAX模型；文獻(xiàn)[12]注意到大多數(shù)電價序列均是非平穩(wěn)隨機過程的特點，提出一種基于自回歸積分滑動平均ARIMA（autoregressive integrated moving average）模型的電價預(yù)測方法；文獻(xiàn)[13]建立了一個基于ARMA模型的分時段電價預(yù)測模型，使得對價格飛升（price spikes）的預(yù)測準(zhǔn)確度得到較大提高。文獻(xiàn)[14-15]分別建立了將誤差校正、小波變換與ARIMA相結(jié)合的電價預(yù)測模型。文獻(xiàn)[16]將電價分解成工作日和周末兩個時間序列，通過移動平均法和離散傅里葉變換將這兩個時間序列分成趨勢分量、周期分量和隨機分量三個組成部分，然后采用移動平均法、外推法和最小二乘支持向量機對這三個組成部分進(jìn)行預(yù)測，進(jìn)一步提高了電價預(yù)測的準(zhǔn)確性。但這些模型均基于電價服從常方差正態(tài)分布的假設(shè)，不能有效處理電價序列的高階矩波動。

使用條件方差、條件偏度和條件峰度來刻畫電價序列的二階矩（方差）、三階矩（偏度）和四階矩（峰度）的時變特征，使用正弦函數(shù)和負(fù)荷平方來刻畫電價序列的多重周期性及其與負(fù)荷之間的相關(guān)性，建立一個基于正態(tài)分布概率密度函數(shù)的Gram-Charlier展開的多周期廣義自回歸條件異方差均值GARCH-M（general autoregressive conditional heteroskedasticity in mean）模型。該模型待估參數(shù)少、計算速度快，可同時處理電價序列的趨勢變化、高階矩波動、多重周期及其與負(fù)荷的相關(guān)性。對PJM電力市場歷史數(shù)據(jù)的分析表明：系統(tǒng)負(fù)荷對電價均值具有顯著影響，電價序列具有周、半月、月、雙月、季、半年等多重周期和二階矩波動集聚性，二階矩、三階矩和四階矩具有明顯時變特征，三階矩和四階矩的尖峰跳躍具有同步性。

1 模型及求解方法

1.1 多周期GARCH-M模型

電價預(yù)測模型可看作一個多輸入單輸出系統(tǒng)，輸出為當(dāng)期電價，輸入為系統(tǒng)負(fù)荷、參與者的報價策略、燃料價格、季節(jié)、氣候等影響因素。考慮到系統(tǒng)負(fù)荷和電價在各個電力市場中均是公開信息，因此本文選擇當(dāng)期負(fù)荷、歷史負(fù)荷與歷史電價作為輸入，并使用時間序列分析的GARCH-M模型對這一系統(tǒng)進(jìn)行描述。設(shè)pt表示t期現(xiàn)貨電價，dt2表示t期系統(tǒng)負(fù)荷平方，εt表示t期殘差，則描述電價變化的多周期GARCH-M模型可表述為

式中：B表示滯后算子；It-1表示t-1期的可用信息集合；ht表示εt的條件二階矩（方差）；zt表示歸一化殘差；u、v、p、q分別為均值方程中、pt及εt的滯后階數(shù)；rh、sh分別為二階矩方程中ht和εt2的滯后階數(shù)；m為電價序列一年內(nèi)的變化周期數(shù)；（ft）描述電價序列的趨勢和季節(jié)性變化；dwkd是一個描述工作日和休息日用電差異的虛擬變量（工作日取值為1、休息日取值為0）。α=（α0，α1，α2，α11，…，α1m，α21，…，α2m）、γ=（γ0，…，γu）、φ=（φ1，…，φp）、κ=（κ1，…，κq）、θ=（θ0，…，θ）v、β=（β0，β11，…，β1rh，β21，…，β2sh）為待估參數(shù)。使用多個正弦函數(shù)可允許電價序列存在多個周期，每個周期的幅度和峰值位置分別由α1i和α2i描述。為保證二階矩為正，要求β0〉0，β1i，β2j≥0，?i∈[1，rh]，j∈[1，sh]。

1.2 參數(shù)估計

在對模型進(jìn)行參數(shù)估計時，需要對殘差的分布做出假定。設(shè)zt的條件概率密度函數(shù)符合正態(tài)分布密度函數(shù)的Gram-Charlier展開（在四階矩截斷），則εt的條件概率密度函數(shù)可表示為[17]2

式中：ηt和λt分別表示εt的條件三階矩和條件四階矩；δ=（δ0，δ11，…，δ1rη，δ21，…，δ2sη）表示三階矩方程的待估參數(shù)；τ=（τ0，τ11，…，τ1rλ，τ21，…，τ2sλ）表示四階矩方程的待估參數(shù)。為保證四階矩為正，要求τ0〉0，τ1i，τ2j≥0，?i∈[1，rλ]，j∈[1，sλ]。

若記ξ=（α，φ，θ，γ，κ，β，δ，τ），則條件對數(shù)似然函數(shù)可表示為

式中：lt（ξ）=ln g（εt|It-1）表示t期觀察值的條件對數(shù)似然函數(shù)。通過最大化L（ξ），可以獲得模型參數(shù)的估計值ξ^。由于L（ξ）是關(guān)于ξ的非線性函數(shù)，初始值的選取關(guān)系到L（ξ）能否收斂到全局最大值。為了增加估計結(jié)果的精度，本文采用逐次逼近的求解方法，即先求解簡化模型，然后以其結(jié)果作為求解復(fù)雜模型的初始值。

1.3 模型檢驗

Nyblom統(tǒng)計量可用于檢驗?zāi)Ｐ偷姆€(wěn)定性[18]，其漸近分布只依賴于待估參數(shù)個數(shù)。Nyblom統(tǒng)計量WN可表示為

Nyblom統(tǒng)計量也可用于檢驗單個系數(shù)的穩(wěn)定性，相應(yīng)于第k個待估參數(shù)的Nyblom統(tǒng)計量WN，k為

Cramer-Von Mises統(tǒng)計量可用于檢驗殘差是否符合正態(tài)分布的Gram-Charlier展開。若記Gram-Charlier展開的累積分布函數(shù)為FN（z），殘差的累積分布函數(shù)為F（z），則Cramer-Von Mises統(tǒng)計量WCVM為

2 PJM電力市場實證研究

本文的研究樣本來源于美國PJM電力市場2007年6月1日至2010年9月9日的日平均現(xiàn)貨電價和負(fù)荷[19]，樣本總數(shù)為1197。表1給出了日平均現(xiàn)貨電價和負(fù)荷序列的描述性統(tǒng)計結(jié)果。從表1可看出，兩個樣本的三階矩和四階矩都顯著異于正態(tài)，現(xiàn)貨電價和負(fù)荷序列均呈現(xiàn)明顯的右偏形態(tài)，同時具有較為明顯的尖峰胖尾特征。J-B統(tǒng)計量非常顯著，說明樣本期間內(nèi)現(xiàn)貨電價和負(fù)荷的分布具有非正態(tài)性。

表1 樣本數(shù)據(jù)的描述性統(tǒng)計結(jié)果Tab.1Descriptive statistical results of the sample datas

對樣本數(shù)據(jù)的趨勢圖、相關(guān)系數(shù)及偏相關(guān)系數(shù)進(jìn)行分析，可將多周期GARCH-M模型具體化為（m，p，q，u，v，rh，sh，rη，sη，rλ，sλ），取值分別為（52，1，3，1，0，1，1，1，1，1，1）。表2給出了該模型的極大似然估計結(jié)果。對表2的數(shù)據(jù)進(jìn)行分析，可以得到如下結(jié)論。

（1）Cramer-Von Mises統(tǒng)計量和各參數(shù)的Nyblom統(tǒng)計量均小于99%置信限的臨界值，表明截止到四階矩的正態(tài)概率密度函數(shù)的Gram-Charlier展開能很好描述殘差的實際分布，各參數(shù)估計值穩(wěn)定有效，可用該模型進(jìn)行預(yù)測及分析。

（2）參數(shù)α1i，α2i，i∈（2，4，6，12，24，52）的t統(tǒng)計量均在95%置信水平上顯著，說明電價序列存在周、半月、月、雙月、季、半年等多周期效應(yīng)，且周周期的峰值幅度較大。

（3）參數(shù)α2的t統(tǒng)計量在99%置信水平上顯著，表明工作日和休息日的系統(tǒng)負(fù)荷對電價均值的影響存在較大差異。

（4）參數(shù)γ0和γ1的t統(tǒng)計量均在99%置信水平上顯著，表明對電價均值具有顯著影響。但在均值方程中包含項后，dwkd的系數(shù)α2的符號由正變負(fù)，表明dwkd和間存在一定替代效應(yīng)。

表2 GARCH-M模型的估計結(jié)果Tab.2Estimation results of GARCH-M model

（5）二階矩方程中的參數(shù)β11的t統(tǒng)計量在99%置信水平上顯著，表明二階矩波動具有較強的記憶性，前期波動對本期波動的影響呈衰減趨勢，電價序列存在明顯波動集聚現(xiàn)象。

（6）二階矩方程中的參數(shù)β21的t統(tǒng)計量在99%置信水平上顯著，表明外部沖擊對二階矩的波動具有加強作用。β11和β21之和接近于1，表明電價序列可能存在共積GARCH效應(yīng)。

（7）三階矩方程中的參數(shù)δ21的t統(tǒng)計量在90%的置信水平上顯著，表明外部沖擊對三階矩的波動具有加強作用。δ11的t統(tǒng)計量在90%的置信水平上不顯著，說明電價序列的三階矩不存在明顯的波動持續(xù)性，這一現(xiàn)象不同于傳統(tǒng)金融市場的收益序列。圖1給出了電價序列的條件三階矩，從圖中可見，在樣本期間該序列表現(xiàn)出明顯的右偏特征（三階矩絕大部分集中在0.2左右），其波動不具有明顯的集聚現(xiàn)象。

圖1 電價序列的三階矩Fig.1Third moment of electricity price series

（8）四階矩方程中的參數(shù)τ21的t統(tǒng)計量在90%置信水平上顯著，表明外部沖擊對四階矩波動具有加強作用，但強度要弱于三階矩和二階矩（|τ21|〈|δ21|〈|β21|）。τ11的t統(tǒng)計量在90%置信水平上不顯著，說明電價序列的四階矩不存在明顯波動集聚現(xiàn)象，如圖2。觀察圖1與圖2可發(fā)現(xiàn)，三階矩和四階矩尖峰跳躍具有明顯的同步時變特征。

圖2 電價序列的四階矩Fig.2Fourth moment of electricity price series

3 結(jié)語

在對電力市場現(xiàn)貨電價的影響因素和變化規(guī)律分析的基礎(chǔ)上，使用條件方差、條件偏度和條件峰度來刻畫電價序列的二階矩、三階矩和四階矩的時變特征，使用正弦函數(shù)和負(fù)荷平方來刻畫電價序列的多重周期性及其與負(fù)荷之間的相關(guān)性，建立了一個基于正態(tài)分布密度函數(shù)的Gram-Charlier展開的多周期GARCH-M模型。該模型待估參數(shù)少、計算速度快，可同時處理電價序列的趨勢變化、多重周期、高階矩波動及其與負(fù)荷之間的相關(guān)性。對PJM電力市場2007年6月1日至2010年9月9日的歷史數(shù)據(jù)的分析表明：系統(tǒng)負(fù)荷對電價均值具有顯著影響，電價序列具有周、半月、月、雙月、季、半年等多重周期和二階矩波動集聚性，二階矩、三階矩和四階矩具有明顯的時變特征，三階矩和四階矩的尖峰跳躍具有同步性。但虛擬變量dwkd和負(fù)荷平方dt2之間的替代性導(dǎo)致dwkd系數(shù)為負(fù)，因此如何更好地處理兩者之間的關(guān)系是下一步要解決的主要問題。

[1]陳思杰，周浩（Chen Sijie，Zhou Hao）．電力市場電價預(yù)測方法綜述（Electricity price forecast methods of electricity market）[J]．繼電器（Relay），2006，34（11）：54-60．

[2]Szkuta B R，Sanabria L A，Dillon T S.Electricity price short-term forecasting using artificial neural networks[J]. IEEE Trans on Power Systems，1999，14（3）：851-857.

[3]Catalao J P S，Mariano S J P S，Mendes V M F，et al. Short-term electricity prices forecasting in a competitive market：a neural network approach[J].Electric Power System Research，2007，77（10）：1297-1304.

[4]Guo Jau-Jia，Luh Peter B.Selecting input factors for clusters of Gaussian radial basis function networks to improve market clearing price prediction[J].IEEE Trans on Power Systems，2003，18（2）：665-672.

[5]張顯，王錫凡，陳芳華，等（Zhang Xian，Wang Xifan，Chen Fanghua，et al）．分時段短期電價預(yù)測（Short-term electricity price forecasting based on period-decoupled price sequence）[J]．中國電機工程學(xué)報（Proceedings of the CSEE），2005，25（15）：1-6．

[6]Hong Y-Y，Hsiao C-Y.Locational marginal price forecasting in deregulated electricity markets using artificial intelligence[J].IEE Proceedings-Generation，Transmission and Distribution，2002，149（5）：621-626.

[7]涂啟玉，張茂林（Tu Qiyu，Zhang Maolin）．小波神經(jīng)網(wǎng)絡(luò)預(yù)測電價的新改進(jìn)（Electricity price forecasting using wavelet neural networks optimized by GA）[J].電力系統(tǒng)及其自動化學(xué)報（Proceedings of the CSU-EPSA），2011，23（2）：157-160．

[8]Zhang Li，Luh Peter B.Neural network-based market clearing price prediction and confidence interval estimation with an improved extended Kalman filter method[J]. IEEE Trans on Power Systems，2005，20（1）：59-66.

[9]Amjady N.Day-ahead price forecasting of electricity markets by a new fuzzy neural network[J].IEEE Trans on Power Systems，2006，21（2）：887-896.

[10]Fan S，Mao C，Chen L.Next-day electricity price forecasting using a hybrid network[J].IET Generation，Transmission and Distribution，2007，1（1）：176-182.

[11]王瑞慶，王晛，李渝曾（Wang Ruiqing，Wang Xian，Li Yuzeng）．基于時間序列ARMAX模型的短期電價預(yù)測方法（Short-term electricity price forecasting Methods based on ARMAX models）[J]．華東電力（East China Electric Power），2009，37（5）：727-730．

[12]Contreras J，Espinola R，Nogales F J，et al.ARIMA models to predict next-day electricity prices[J].IEEE Trans on Power Systems，2003，18（3）：1014-1020.

[13]Crespo Cuaresma Jesus，Hlouskova Jaroslava，Kossmeier Stephan，et al.Forecasting electricity spot-prices using linear univariate time-series models[J].Applied Energy，2004，77（1）：87-106.

[14]周明，嚴(yán)正，倪以信，等（Zhou Ming，Yan Zheng，Ni Yixin，et al）．含誤差預(yù)測校正的ARIMA電價預(yù)測新方法（A novel ARIMA approach on electricity price forecasting with the improvement of predicted error）[J]．中國電機工程學(xué)報（Proceedings of the CSEE），2004，24（12）：63-68．

[15]Conejo A J，Plazas M A，Espinola R，et al.Day-ahead electricity price forecasting using the wavelet transform and ARIMA models[J].IEEE Trans on Power Systems，2005，20（2）：1035-1042.

[16]熊高峰，韓鵬，聶坤凱（Xiong Gaofeng，Han Peng，Nie Kunkai）．時間序列分解在短期電價分析與預(yù)測中的應(yīng)用（Application of time series decomposition to short-term electricity price analysis and forecasting）[J]．電力系統(tǒng)及其自動化學(xué)報（Proceedings of the CSU-EPSA），2011，23 （3）：95-100．

[17]Angel Leon，Gonzalo Rubio，Gregorio Serna.Autoregressive conditional volatility，skewness and kurtosis[J].The Quarterly Review of Economics and Finance，2005，45（2-3）：599-618.

[18]Nyblom J.Testing for the constancy of parameters over time[J].Journal of the American Statistical Association，1989，84（405）：223-230.

[19]PJM Energy Market.Day-ahead energy markets[EB/OL]. http：//www.pjm.com/markets-and-operations/energy/dayahead.aspx，2012.

Higher Moments Volatility Characteristics of Electricity Price Series

WANG Rui-qing，WANG Hong-fu
（School of Computer&Information Engineering，Anyang Normal University，Anyang 455000，China）

The distribution properties of electricity prices are the important information for the risk management of electricity markets and the pricing of electricity financial derivatives.With comprehensive consideration of the changing rules of the electricity spot price，a multiple periodic GARCH-M model based on Gram-Charlier series expansion of the normal probability density function is proposed，in which the second，third and fourth moments of electricity price series are described by conditional variance，skewness and kurtosis，and the multiple periods and the relationship among load and spot price are also taken into account by sine function and system load squares.The numerical example based on the historical data of the PJM market shows that the system loads have a significant effect on the mean electricity prices，there exists second moment volatility clustering and weekly，semi-monthly，monthly，bimonthly，quarterly and semi-annual periods，and the second，third and fourth moments manifest the clear time-varying characteristics and there also exists synchronization between the peak jumping of third and fourth moments.

electricity price；electric load；higher moments；volatility clustering；GARCH-M model

TM73

A

1003-8930（2013）04-0058-05

王瑞慶（1965—），男，博士，教授，研究方向為電力運籌與控制、電力市場風(fēng)險管理、數(shù)據(jù)挖掘與智能計算。Email：ayrqwang@163.com

2012-01-16；

2012-04-02

河南省教育廳自然科學(xué)研究計劃項目（2010B120002）

王宏福（1979—），男，碩士，講師，研究方向為知識管理、數(shù)據(jù)挖掘與智能計算。Email：anytc123@163.com