高小童，秦志龍

（1.山東省膠東調(diào)水工程棘洪灘水庫管理處，山東　青島　266001；2.國網(wǎng)山東省電力公司青島供電公司，山東　青島　266002）

基于混合Copula函數(shù)的含風(fēng)電場電力系統(tǒng)可靠性評估

高小童1，秦志龍2

（1.山東省膠東調(diào)水工程棘洪灘水庫管理處，山東青島266001；2.國網(wǎng)山東省電力公司青島供電公司，山東青島266002）

隨著風(fēng)力發(fā)電的發(fā)展，風(fēng)電場風(fēng)速之間的相關(guān)性對電力系統(tǒng)可靠性的影響越來越大?；诨旌螩opula函數(shù)建立一種多元相關(guān)性風(fēng)速的模型，并在含多個風(fēng)電場的IEEE-RTS系統(tǒng)可靠性評估中應(yīng)用該模型，通過實例分析，在描述風(fēng)電場風(fēng)速相關(guān)結(jié)構(gòu)方面，提出的混合Copula函數(shù)法比單一Copula函數(shù)法更合適，更有優(yōu)勢。

風(fēng)速相關(guān)性；混合Copula函數(shù)；電力系統(tǒng)可靠性；蒙特卡洛仿真

0　引言

風(fēng)能是一種清潔的可再生能源，大力開發(fā)風(fēng)電能源是合理調(diào)整電源結(jié)構(gòu)的需要，也是國內(nèi)能源發(fā)展戰(zhàn)略的重要組成部分［1］。文獻［1-4］在研究含風(fēng)電場風(fēng)速模型過程中僅僅考慮了風(fēng)速隨機性，沒有計及風(fēng)電場風(fēng)速之間相關(guān)性。隨著風(fēng)電場的密集建設(shè)，距離較近風(fēng)電場之間的風(fēng)速相關(guān)性研究是確保含多個風(fēng)電場電力系統(tǒng)可靠性的基礎(chǔ)［5-11］。文獻［5-8］提出了一種產(chǎn)生線性相關(guān)性風(fēng)速的方法，但是該方法不能模擬風(fēng)速之間的非線性相關(guān)性；文獻 ［9-11］均僅將單一Copula函數(shù)理論運用到含風(fēng)電場的電力系統(tǒng)可靠性評估中。與計及風(fēng)速相關(guān)性比較，不計風(fēng)速相關(guān)性會造成對含風(fēng)電場發(fā)電系統(tǒng)可靠性樂觀估計［8，11］；與混合Copula函數(shù)法相比，單一Copula函數(shù)將高估風(fēng)速相關(guān)性對電力系統(tǒng)可靠性的影響。

提出混合Copula函數(shù)理論來模擬風(fēng)電場風(fēng)速相關(guān)性結(jié)構(gòu)，并應(yīng)用于含風(fēng)電場發(fā)電系統(tǒng)可靠性評估。

1　基于混合Copula函數(shù)的風(fēng)速相關(guān)性模型

1.1單一Copula函數(shù)

阿基米德Copula函數(shù)具有構(gòu)造簡單、對稱性和可結(jié)合性良好等特性［11-13］。

常用的阿基米德Copula函數(shù)有Clayton Copula，Gumbel Copula和Frank Copula函數(shù)［11-13］，表達式分別為

式中：α為相關(guān)參數(shù)；（u1，…，uN）為N元隨機變量。

Clayton Copula、Gumbel Copula密度函數(shù)具有非對稱性的相關(guān)結(jié)構(gòu)，無法表達變量間對稱的相關(guān)關(guān)系；Frank Copula函數(shù)具有對稱性的相關(guān)結(jié)構(gòu)，無法表達變量間非對稱的相關(guān)關(guān)系。Clayton Copula函數(shù)適用于描述分布下尾部的變化，但Gumbel Copula并不適用；Clayton Copula函數(shù)不適用于描述分布上尾部的變化，而Gumbel Copula適用；Frank Copula對描述分布上、下尾相關(guān)性的變化均不適用［13］。

Clayton Copula、Gumbel Copula和Frank Copula函數(shù)，在描述風(fēng)電場風(fēng)速間相關(guān)結(jié)構(gòu)方面各有不同特點［13］，某個單一Copula函數(shù)很難表達風(fēng)電場風(fēng)速多樣性、差異性的特點。

1.2混合Copula函數(shù)

在實際應(yīng)用中應(yīng)采用混合Copula函數(shù)來表達風(fēng)電場的相關(guān)風(fēng)速數(shù)據(jù)，混合后的Copula函數(shù)為［14］

式中：Cn（u1，u2；αn）為單一 Copula函數(shù)；λn為權(quán)重系數(shù)，滿足 0≤λn≤1；αn為相關(guān)參數(shù)［14］。不同的Copula函數(shù)組合，在描述風(fēng)電場風(fēng)速結(jié)構(gòu)時的效果也有所不同。歐式距離檢驗法、P-P圖檢驗法、K-S檢驗法是檢驗Copula函數(shù)擬合效果的幾種常用方法［15］。

選用 Gumbel、Frank、Clayton 3種 Copula函數(shù)的線性組合構(gòu)造反映風(fēng)電場風(fēng)速相關(guān)性的混合Copula函數(shù)［13-15］，能夠反映風(fēng)電場風(fēng)速變化的各種情況。

2　相關(guān)風(fēng)速產(chǎn)生方法

2.1單一Copula函數(shù)模型中的分步參數(shù)估計

利用分步參數(shù)估計法估計邊緣分布函數(shù)和Copula函數(shù)中的未知參數(shù)［8，11］。

首先要確定風(fēng)速的邊緣分布，采用非參數(shù)核密度分布模型來表達風(fēng)速邊緣分布［16］。非參數(shù)核密度分布密度函數(shù)為

非參數(shù)核密度分布累積函數(shù)為

式中：n為樣本數(shù)；Xi為樣本值；h為帶寬；K （·）為核函數(shù)；Φ（·）為單變量標(biāo)準(zhǔn)正態(tài)累積分布函數(shù)。

得出邊緣分布的參數(shù)后，用擬牛頓法得到相關(guān)性參數(shù)α的估計值［8，11］。

式中：c為單一Copula密度函數(shù)；n為樣本數(shù)；Xi為樣本值；h為帶寬；N為Copula函數(shù)的元數(shù)。

2.2混合Copula函數(shù)權(quán)重系數(shù)

混合Copula函數(shù)權(quán)重系數(shù)通過理論Copula函數(shù)和經(jīng)驗Copula函數(shù)之間的最短歐式距離來選取。經(jīng)驗Copula函數(shù)為［8，11-12］

理論Copula函數(shù)與經(jīng)驗Copula函數(shù)之間的最短歐式距離為［8，11-12］

用歐式距離最小的理論求解選取混合Copula函數(shù)權(quán)重系數(shù)λn。

2.3多維風(fēng)電場相關(guān)風(fēng)速

描述風(fēng)速相關(guān)結(jié)構(gòu)的最優(yōu)混合Copula函數(shù)確定后，采用一種合成法產(chǎn)生具有相關(guān)性的多維風(fēng)速樣本，步驟為：

1）產(chǎn)生U（0，1）中的一個隨機數(shù)t。

2）確定i∈｛1，2，…N｝使得

3）產(chǎn)生Ci（u1，u2）的一個隨機數(shù)對（u1，u2）。

4）重復(fù)上述3步T次，得出服從C（u1，u2）的T個隨機數(shù)對。

3　算例研究

選用的風(fēng)速數(shù)據(jù)為山東青島地區(qū)兩風(fēng)電場風(fēng)速歷史數(shù)據(jù)。

3.1選取最優(yōu)混合Copula函數(shù)［8，11］

采用分步參數(shù)估計法得到單一Copula函數(shù)的參數(shù)，如表1所示。

表1　單一Copula函數(shù)參數(shù)估計值

結(jié)合表1中Copula函數(shù)的參數(shù)α，由公式（10）可以求解出混合Copula函數(shù)的權(quán)重系數(shù)λ1=0.05，λ2=0.17，λ3=0.78。

選取K-S檢驗法、歐式距離檢驗法對多個風(fēng)電場相關(guān)風(fēng)速的擬合效果進行檢驗［15］，如表2所示。

表2　風(fēng)電場風(fēng)速的擬合效果檢驗方法

由表2可知，與單一Copula分布函數(shù)相比，混合Copula分布函數(shù)具有最短歐式距離。根據(jù)K-S檢驗（H值），能通過K-S檢驗僅僅只有混合Copula函數(shù)（H=0），其他單一Copula函數(shù)（H=1）都沒有通過K-S檢驗。因此相比單一Copula分布函數(shù)，混合Copula分布函數(shù)更合適描述2個風(fēng)電場風(fēng)速的相關(guān)結(jié)構(gòu)。

3.2含多個風(fēng)電場的發(fā)電系統(tǒng)可靠性評估

將同一地區(qū)兩個容量為160 MW風(fēng)電場接入到IEEE-RTS可靠性測試系統(tǒng)［17］中。由風(fēng)速與輸出功率之間的函數(shù)計算出風(fēng)電機組的輸出功率［6］。風(fēng)電機組的切入、切出、額定速度分別為4.0 m/s、20.0 m/s、11.1 m/s，其強迫停運率為0.05［8，11］。

IEEE-RTS系統(tǒng)在2個風(fēng)電場接入前后的系統(tǒng)可靠性指標(biāo)如表3所示。缺電時間期望（LOLE）和電量不足期望（LOEE）［8，11］是評估發(fā)電系統(tǒng)可靠性的2個重要指標(biāo)。其中，情況1為不考慮風(fēng)電場風(fēng)速之間的相關(guān)性；情況2為考慮風(fēng)電場風(fēng)速之間的相關(guān)性，并基于單一Copula函數(shù)法建立風(fēng)電場風(fēng)速相關(guān)性模型；情況3為考慮風(fēng)電場風(fēng)速之間的相關(guān)性，并基于混合Copula函數(shù)法建立風(fēng)電場風(fēng)速相關(guān)性模型。

表3　IEEE-RTS系統(tǒng)在2個風(fēng)電場接入前后的系統(tǒng)可靠性指標(biāo)

由表3可知，在2個風(fēng)電場接入IEEE-RTS系統(tǒng)后，其系統(tǒng)可靠性明顯有所提高；發(fā)電系統(tǒng)中的風(fēng)電場之間如不計風(fēng)速相關(guān)性會造成含多個風(fēng)電場發(fā)電系統(tǒng)可靠性的評估過于樂觀。與此同時，混合Copula函數(shù)法和單一Copula函數(shù)法分別應(yīng)用于含多個風(fēng)電場發(fā)電系統(tǒng)可靠性評估中，與混合Copula函數(shù)法相比，單一Copula函數(shù)法也會高估風(fēng)速之間相關(guān)性對系統(tǒng)可靠性產(chǎn)生的影響。

研究不同裝機容量的2個風(fēng)電場接入IEEERTS系統(tǒng)后其可靠性水平。針對不同容量的風(fēng)電場，混合Copula函數(shù)法和單一Copula函數(shù)法分別應(yīng)用于含多個風(fēng)電場發(fā)電系統(tǒng)可靠性評估，其系統(tǒng)可靠性水平如表4、表5所示，兩者之間的相對誤差如表6所示。

表4　可靠性指標(biāo)（混合Copula函數(shù)法）

表5　可靠性指標(biāo)（單一Copula函數(shù)法）

表6　可靠性指標(biāo)間的相對誤差

由表4～6可知，隨著風(fēng)電場裝機容量的增大，應(yīng)用單一Copula函數(shù)法與混合Copula函數(shù)法得出系統(tǒng)可靠性指標(biāo)之間的相對誤差增大。綜上所述隨著越來越多風(fēng)電場加入電力系統(tǒng)中，建立精確的風(fēng)電場風(fēng)速相關(guān)模型越來越重要。

4　結(jié)語

基于混合Copula函數(shù)建立了一種多元相關(guān)性風(fēng)速模型，并在含多個風(fēng)電場的發(fā)電系統(tǒng)可靠性評估中應(yīng)用該模型。通過對多個風(fēng)電場加入IEEE-RTS系統(tǒng)進行算例分析，與混合Copula函數(shù)法相比，單一Copula函數(shù)法高估了風(fēng)速的相關(guān)性對系統(tǒng)可靠性的影響。由此可見，風(fēng)電場風(fēng)速相關(guān)性模型的精確程度是其系統(tǒng)可靠性評估計算中的關(guān)鍵因素，且所提出的混合Copula函數(shù)法在描述風(fēng)電場風(fēng)速相關(guān)性方面比單一Copula函數(shù)法更合適，更有優(yōu)勢。所建立的精確風(fēng)速相關(guān)性模型為以后青島地區(qū)乃至山東地區(qū)風(fēng)電場的規(guī)劃打下了理論基礎(chǔ)。

［1］吳林偉，張建華，劉若溪.考慮風(fēng)電機組故障的風(fēng)電場可靠性模型及其應(yīng)用［J］.電力系統(tǒng)自動化，2012，36（16）：31-35.

［2］劉威，趙淵，周家啟，等.計及風(fēng)電場的發(fā)輸配電系統(tǒng)可靠性評估［J］.電網(wǎng)技術(shù)，2008，32（13）：69-74.

［3］吳昊，張焰，劉波.考慮風(fēng)電場影響的發(fā)輸電系統(tǒng)可靠性評估［J］.電力系統(tǒng)保護與控制，201l，39（4）：36-41.

［4］張碩，李庚銀，周明.含風(fēng)電場的發(fā)輸電系統(tǒng)可靠性評估［J］.中國電機工程學(xué)報，2010，30（7）：8-14.

［5］XIE Kaigui，BILLITON R.Considering wind speed correlation of WECS in reliability evaluation using the time-shifting technique ［J］.Electric Power System Research，2009，79（4）：687-693.

［6］WANGDEEW，BILLITONR.Consideringload-carrying capability and wind speed correlation of WECS in generation adequacy as sessment［J］.IEEE Trans on Power Systems，2006，21（3）：734-741.

［7］FEIJOO A E，CIDRAS J，DORNELAS J L G.Wind speed simulation in wind farms forestay-state security assessment of electri cal power system［J］.IEEE Trans on Energy Conversion，1999，14（4）：1 582-1 588.

［8］秦志龍.計及相關(guān)性的含風(fēng)電場和光伏電站電力系統(tǒng)可靠性評估［D］.重慶：重慶大學(xué)，2013.

［9］GEORGE Papaefthymiou，DOROTA Kurowicka.Using Copula for modeling stochastic dependence in power system uncertainly analysis［J］.IEEE Trans on Power System，2009，24（1）：40-49.

［10］徐玉琴，張林浩.考慮風(fēng)速相關(guān)性的風(fēng)電接入能力分析［J］.可再生能源，2014，32（2）：201-206.

［11］秦志龍，李文沅，熊小伏.考慮風(fēng)速相關(guān)性的發(fā)輸電系統(tǒng)可靠性評估［J］.電力系統(tǒng)自動化，2013，37（16）：47-52.

［12］韋艷華，張世英.Copula理論及其在金融分析上的應(yīng)用［M］.北京：清華大學(xué)出版社，2008.

［13］莊丹琴.Copula函數(shù)的選擇及其在金融分析中的若干應(yīng)用［D］.馬鞍山：安徽工業(yè)大學(xué)，2010.

［14］季峰，蔡興國，王俊.基于混合Copula函數(shù)的風(fēng)電功率相關(guān)性分析［J］.電力系統(tǒng)自動化，2014，38（2）：1-6.

［15］王小紅，周步祥，傅利，等.基于綜合Copula函數(shù)的風(fēng)電功率相關(guān)性分析及其在無功優(yōu)化中的應(yīng)用［J］.可再生能源，2014，32 （8）：1 111-1 117.

［16］QIN Zhilong，LI Wenyuan，XIONG Xiaofu.Estimating wind speed probability distribution using kernel density method［J］. Electric Power Systems Research，2011，81（12）：2 139-2 146.

［17］BILLITON R，LI Wenyuan.Reliability assessment of electrical power systems using Monte Carlo methods［M］.New York：Plenum Press，1994.

Reliability Evaluation of Power System with Wind Farms Based on Mix-Copula Function

GAO Xiaotong1，QIN Zhilong2
（1.Jihongtan Reservoir Management Office of Jiaodong Water Diversion Project in Shandong Province，Qingdao 266001，China；2.State Grid Qingdao Power Supply Company，Qingdao 266002，China）

With the development of wind power，the relationship between the wind speed of wind farms has more and more influence on power system reliability.Based on mix-Copula functions，a model of multi-dimensional correlation wind speed is presented，and is applied to the IEEE-RTS reliability test system containing multiple wind farms.Through the example analysis，the mixed Copula function is more appropriate and has more advantages than single copula function for the description of the wind speed dependence structure.

wind speed correlation；mix-copula function；power system reliability；Monte Carlo simulation

TM732

A

1007－9904（2016）08－0025－04

2016-04-21

高小童（1985），女，工程師，從事水利、風(fēng)電等清潔能源研究工作；秦志龍（1982），男，工程師，從事電網(wǎng)調(diào)度與控制的相關(guān)工作。