寧楠,范俊秋,徐潘宇馳,程培軍,林盛振,謝敏
(1.貴州電網(wǎng)有限責(zé)任公司貴安供電局,貴州 貴安 550003;2.廣東電網(wǎng)有限責(zé)任公司天河供電局,廣東 廣州 510075 3.華南理工大學(xué),廣東 廣州 510640)
高比例新能源并入的電力系統(tǒng)具有良好的應(yīng)用前景[1]。與此同時(shí),新能源自身固有的波動(dòng)性、隨機(jī)性和間歇性[2]等特點(diǎn),也為電力系統(tǒng)的安全經(jīng)濟(jì)調(diào)度過(guò)程帶來(lái)了復(fù)雜的不確定性沖擊。為了準(zhǔn)確應(yīng)對(duì)高維新能源并網(wǎng)后引入的調(diào)度風(fēng)險(xiǎn)問(wèn)題,有必要對(duì)高維新能源發(fā)電的相關(guān)性建模以及相應(yīng)的系統(tǒng)動(dòng)態(tài)經(jīng)濟(jì)調(diào)度問(wèn)題進(jìn)行研究[3-6]。
隨著新能源的并網(wǎng)容量進(jìn)一步增長(zhǎng),多個(gè)新能源電場(chǎng)基于可再生能源類型特點(diǎn)、所處地理位置和氣象因素的相似性[7],存在高維變量出力相關(guān)性較強(qiáng)的特征。Copula函數(shù)在新能源出力建模中得到了廣泛的應(yīng)用[8-11]。但是,當(dāng)系統(tǒng)中接入的新能源維度增加時(shí),單一的Copula函數(shù)難以對(duì)高維變量展開建模分析。藤結(jié)構(gòu)的pair-Copula模型[12-13]使用二元Copula函數(shù)作為基礎(chǔ),根據(jù)多種逐層合并的方式,對(duì)應(yīng)形成多元Copula的聯(lián)合分布函數(shù),能夠有效提高建模的準(zhǔn)確性。
為順利求解含高維新能源的動(dòng)態(tài)經(jīng)濟(jì)調(diào)度問(wèn)題,已有研究通常引入隨機(jī)規(guī)劃理論,從而解決經(jīng)濟(jì)調(diào)度模型中新能源出力的間歇、波動(dòng)以及隨機(jī)問(wèn)題[14]。典型方案有機(jī)會(huì)約束規(guī)劃法[15-17]、場(chǎng)景分析法[18-20]等。以場(chǎng)景分析為基礎(chǔ)的隨機(jī)規(guī)劃模型,主要將抽樣生成的若干誤差場(chǎng)景作為紐帶,不斷使隨機(jī)優(yōu)化向確定優(yōu)化進(jìn)行轉(zhuǎn)換[21-22],但其應(yīng)用的有效性受限于場(chǎng)景數(shù)量的選擇。分布式優(yōu)化技術(shù)通過(guò)合理分解建立在場(chǎng)景分析基礎(chǔ)上的經(jīng)濟(jì)調(diào)度問(wèn)題,盡可能全面地考慮了各個(gè)場(chǎng)景的貢獻(xiàn)[23]。但是,已有研究中對(duì)于高維新能源間的出力相關(guān)性考慮不夠精細(xì)。
針對(duì)以上研究存在的不足,本文首先基于動(dòng)態(tài)C藤pair-Copula理論建立模型,能夠準(zhǔn)確描述高維新能源間的出力相關(guān)特性。然后,基于場(chǎng)景分析與多學(xué)科協(xié)同優(yōu)化的思想,設(shè)計(jì)及求解含高維新能源的動(dòng)態(tài)經(jīng)濟(jì)調(diào)度場(chǎng)景協(xié)同優(yōu)化模型,從而有效提升引入高維相依新能源的協(xié)同優(yōu)化經(jīng)濟(jì)調(diào)度的計(jì)算效率。
為充分考慮并網(wǎng)新能源發(fā)電機(jī)組的隨機(jī)特性給系統(tǒng)運(yùn)行產(chǎn)生的影響,本文基于動(dòng)態(tài)C藤pair-Copula函數(shù),形成高維新能源的聯(lián)合出力分布模型,并通過(guò)抽樣技術(shù)完成各新能源出力隨機(jī)誤差場(chǎng)景生成過(guò)程。
Copula思想與模型以其靈活、實(shí)用的特性得到金融分析建模以及隨機(jī)新能源相關(guān)性分析等領(lǐng)域的廣泛關(guān)注與使用。Copula技術(shù)通過(guò)整合多維變量的邊緣分布函數(shù),以及與之相適應(yīng)的單個(gè)Copula函數(shù),能夠方便地獲得其聯(lián)合分布數(shù)學(xué)模型進(jìn)行模擬。為確保高維變量之間的相依特性獲得更加精確的刻畫,pair-Copula函數(shù)在單純的Copula建模思想基礎(chǔ)上,進(jìn)一步對(duì)高維變量進(jìn)行精細(xì)化建模。pair-Copula以二元Copula分布為基礎(chǔ),采取逐層合并的方式來(lái)構(gòu)造多元高維Copula聯(lián)合分布函數(shù)。圖1給出了pair-Copula函數(shù)的兩種常用結(jié)構(gòu)圖,即C-vine(C藤)和D-vine(D藤)。該方法能夠引入多類型的Copula分布,有效提升模型的適用性和異構(gòu)性[24]。假若存在與其他所有新能源發(fā)電均具有較強(qiáng)相關(guān)性的電源,則適宜采用C藤結(jié)構(gòu),否則應(yīng)采用D藤結(jié)構(gòu)。
圖1 基于pair-Copula方法的兩種藤Copula結(jié)構(gòu)
側(cè)重于研究的實(shí)際性,本文通過(guò)對(duì)五類常用的靜態(tài)Copula函數(shù)以及三類動(dòng)態(tài)Copula函數(shù)加以有機(jī)結(jié)合,構(gòu)建基于C藤pair-Copula函數(shù)的高維新能源出力動(dòng)態(tài)模型。文中模型采用的三類動(dòng)態(tài)Copula函數(shù)分別是時(shí)變gaussian Copula函數(shù)、時(shí)變t-Copula函數(shù)與時(shí)變clayton Copula函數(shù)[25-26]。動(dòng)態(tài)Copula函數(shù)的時(shí)變性,指的是其表達(dá)式中具備相應(yīng)的時(shí)變參數(shù)。具體公式如下所示:
時(shí)變gaussian Copula函數(shù)的參數(shù):
(1)
時(shí)變t-Copula函數(shù)的參數(shù):
(2)
時(shí)變clayton Copula函數(shù)的參數(shù):
(3)
式(3)中,ρct為時(shí)變clayton Copula函數(shù)的時(shí)變參數(shù),ωρc、βρc、αρc為ρct的相關(guān)系數(shù)。
1.1 高維新能源的出力相關(guān)性建模
建立n維新能源出力相關(guān)性模型的具體流程如下:
Step 1:求取n維原始變量的邊緣分布函數(shù)。首先,讀取多個(gè)新能源電場(chǎng)的歷史出力數(shù)據(jù)X,并進(jìn)行相應(yīng)的標(biāo)幺化;其次,利用核密度估計(jì)法,求出邊緣概率分布F(x);最后,基于各變量間的線性相關(guān)系數(shù)由大及小的順序原則,重新排列這n個(gè)隨機(jī)變量,并記排序后的隨機(jī)變量為{x1,x2,…,xn},x1即為根節(jié)點(diǎn)。同時(shí),令ui=Fi(xi)。
Step 2:選取最優(yōu)Copula函數(shù)。利用與經(jīng)驗(yàn)Copula函數(shù)的歐氏距離最小,對(duì)該層根節(jié)點(diǎn)與其余節(jié)點(diǎn)間的Copula函數(shù)進(jìn)行選??;此外,相應(yīng)的Copula參數(shù)將基于極大似然函數(shù)法求得。
Step 3:求取下一層節(jié)點(diǎn)的pair-Copula序列。通過(guò)結(jié)合該層節(jié)點(diǎn)的序列,以及下一層節(jié)點(diǎn)的最優(yōu)Copula函數(shù),從而根據(jù)式(4)求取下層各節(jié)點(diǎn)pair-Copula序列。
(4)
式(4)中,Ci,j|k為二元Copula分布函數(shù)。此外,vj代表向量v中的第j個(gè)變量,v-j代表向量v去掉vj后的向量。
Step 4:重復(fù)Step 2和Step 3,直至最后一層節(jié)點(diǎn)的建模結(jié)束,即完成了描述n維新能源出力相關(guān)性的高維動(dòng)態(tài)模型構(gòu)建。
基于場(chǎng)景采樣分析技術(shù),對(duì)動(dòng)態(tài)C藤pair-Copula模型展開采樣,具體過(guò)程如下:
Step 1:生成n個(gè)不僅服從[0,1]均勻分布,而且各自獨(dú)立的隨機(jī)數(shù)wi(i=1,2,…n);
Step 2:設(shè)n維采樣序列為{x1,x2,…,xn},令w1=x1,w2=F(x2|x1),…,wn=F(xn|x1,…,xn-1),根據(jù)式(5)求取n維采樣序列。
(5)
基于上述過(guò)程對(duì)新能源出力的高維動(dòng)態(tài)模型進(jìn)行采樣,獲得采樣結(jié)果后,在場(chǎng)景削減的過(guò)程中引入基于概率距離的快速前代消除思想,有助于后續(xù)完整模型的構(gòu)建。
2.1 目標(biāo)函數(shù)
為了保證電網(wǎng)調(diào)度結(jié)果能夠在安全穩(wěn)定運(yùn)行的情況下兼顧節(jié)能環(huán)保問(wèn)題,以發(fā)電總?cè)己某杀九c廢氣排放成本之和最小建立目標(biāo)函數(shù)。
minF=F1+F2
(6)
式中,F(xiàn)為系統(tǒng)的總運(yùn)行成本,F(xiàn)1為系統(tǒng)的燃油消耗成本,F(xiàn)2為系統(tǒng)的廢氣排放成本;A2,i、A1,i、A0,i表示常規(guī)機(jī)組i的燃油消耗成本系數(shù);PGi(t)代表常規(guī)機(jī)組i在時(shí)段t內(nèi)的有功功率;Cp代表常規(guī)機(jī)組廢氣排放的環(huán)保處理單價(jià);B2,i、B1,i、B0,1代表常規(guī)機(jī)組i的廢氣排放系數(shù);T代表調(diào)度周期,取為24h;N表示常規(guī)機(jī)組數(shù)量。
2.2 約束條件
模型的約束條件如式(7)至式(12)所示。
(7)
式(7)為系統(tǒng)運(yùn)行中的功率平衡約束。其中PRj(t)表示新能源發(fā)電機(jī)組j在時(shí)段t內(nèi)的有功功率;PLoad(t)表示系統(tǒng)在第t時(shí)段的負(fù)荷總和;s為場(chǎng)景區(qū)分標(biāo)志,預(yù)測(cè)場(chǎng)景下有s=0,隨機(jī)誤差場(chǎng)景下則為s≠0;Nr為新能源發(fā)電機(jī)組的數(shù)量。
(8)
(9)
式(8)為系統(tǒng)運(yùn)行中的向上、向下旋轉(zhuǎn)備用約束,式(9)為常規(guī)機(jī)組出力的上限、下限約束。式中,PGi,max、PGi,min分別代表常規(guī)機(jī)組i出力的上限以及下限;α表示負(fù)荷預(yù)測(cè)偏差率,一般取5%~10%。
(10)
式(10)為系統(tǒng)運(yùn)行中的機(jī)組出力爬坡約束。式中,rui、rdi分別代表常規(guī)機(jī)組i向上、向下出力爬坡率,ΔT表示單位調(diào)度時(shí)段。
(11)
式(11)為系統(tǒng)運(yùn)行中的線路有功潮流約束。式中,Pmn(t)為以節(jié)點(diǎn)m、n為兩個(gè)端點(diǎn)的輸電線路mn在時(shí)段t內(nèi)的有功功率;Pmn,max為輸電線路mn的有功功率最大限制值;kl表示關(guān)鍵輸電斷面k所包含的線路條數(shù);Pk,j(t)為在時(shí)段t內(nèi),關(guān)鍵輸電斷面k內(nèi)的第j條線路所輸送的有功功率;Pcut,max(k)表示關(guān)鍵輸電斷面k的有功潮流最大限制值。
(12)
式(12)為系統(tǒng)運(yùn)行中的場(chǎng)景功率轉(zhuǎn)移約束。式中,ΔPGi(t)為該場(chǎng)景內(nèi)的常規(guī)機(jī)組i于時(shí)段t的可轉(zhuǎn)移功率;Ghyd為水電機(jī)組集合,Ei,hyd為第i臺(tái)水電機(jī)組的水電量日最大限制值;Ggas為燃?xì)鈾C(jī)組集合,Ei,gas為第i臺(tái)燃?xì)鈾C(jī)組的氣電量日最大限制值。
多學(xué)科協(xié)同優(yōu)化方法充分融合吸收了分布式優(yōu)化技術(shù)以及二級(jí)優(yōu)化技術(shù)的核心思想[27-28],能夠針對(duì)性地為具有高維特征的模型設(shè)計(jì)高效的多級(jí)優(yōu)化框架,從而將高維、復(fù)雜的原始模型解耦成相對(duì)低維、簡(jiǎn)化的優(yōu)化模型,高效地降低了原始模型的求解難度,也順利解決了“維數(shù)災(zāi)”等難題,是求解復(fù)雜系統(tǒng)優(yōu)化問(wèn)題的有力工具。其本質(zhì)在于,按照適應(yīng)于實(shí)際復(fù)雜系統(tǒng)優(yōu)化問(wèn)題的標(biāo)準(zhǔn)對(duì)原模型進(jìn)行合理分解,從而簡(jiǎn)化原模型,獲得單獨(dú)的主優(yōu)化問(wèn)題和部分相對(duì)簡(jiǎn)單的子優(yōu)化問(wèn)題,并在主優(yōu)化級(jí)對(duì)所有子優(yōu)化問(wèn)題的最優(yōu)解進(jìn)行場(chǎng)景分析與優(yōu)化調(diào)整,最終求得原始模型的最優(yōu)解。
3.1 基于多學(xué)科協(xié)同優(yōu)化的高維場(chǎng)景解耦動(dòng)態(tài)經(jīng)
濟(jì)調(diào)度模型
具體的建模展開如下:
(1)主優(yōu)化級(jí)模型
主優(yōu)化級(jí)模型的目標(biāo)函數(shù),即為式(6)中原系統(tǒng)級(jí)模型的目標(biāo)函數(shù)。主優(yōu)化級(jí)模型的約束條件不僅涵蓋式(7)至式(12)中所有只涉及預(yù)測(cè)場(chǎng)景的部分,還包含主優(yōu)化級(jí)的一致性約束,如式(13)所示:
(13)
(2)子優(yōu)化級(jí)模型
(14)
s.t.式(7),s≠0
如上所述,基于MCO技術(shù)與場(chǎng)景分析方法的高維場(chǎng)景解耦動(dòng)態(tài)經(jīng)濟(jì)調(diào)度模型建模計(jì)算框架如圖2所示。
圖2 基于MCO的高維場(chǎng)景解耦動(dòng)態(tài)經(jīng)濟(jì)調(diào)度模型計(jì)算結(jié)構(gòu)Fig.2 Computational structure of high-dimensional scene decoupling dynamic economic dispatch model based on MCO
基于動(dòng)態(tài)C藤pair-Copula函數(shù)對(duì)高維新能源進(jìn)行建模,并結(jié)合場(chǎng)景分析法與多學(xué)科協(xié)同優(yōu)化方法對(duì)電網(wǎng)動(dòng)態(tài)經(jīng)濟(jì)調(diào)度問(wèn)題展開求解。求解步驟流程如圖3所示,詳細(xì)描述如下:
圖3 基于場(chǎng)景分析的電網(wǎng)動(dòng)態(tài)經(jīng)濟(jì)調(diào)度多學(xué)科協(xié)同優(yōu)化算法流程圖Fig.3 Flow chart of multidisciplinary collaborative optimization algorithm for power grid dynamic economic dispatch based on scenario analysis
Step 1:主優(yōu)化級(jí)模型進(jìn)行數(shù)據(jù)初始化,并傳送到子優(yōu)化級(jí)模型,同時(shí)令迭代次數(shù)k=1。
|(F(k+1)-F(k))/F(k+1)|≤ε
(15)
式(15)中,F(xiàn)(k)、F(k+1)為第k次以及第k+1次迭代求解得到的主優(yōu)化級(jí)模型目標(biāo)函數(shù)值;ε為迭代收斂的精度,設(shè)置為10-3。
仿真系統(tǒng)采用標(biāo)準(zhǔn)IEEE 118節(jié)點(diǎn)系統(tǒng)進(jìn)行建模分析。選擇10個(gè)裝機(jī)容量均為75MW的實(shí)際風(fēng)力發(fā)電場(chǎng)Wi(i=1,2,…,10),并且所有風(fēng)電場(chǎng)全部落在美國(guó)加利福尼亞州的某一片區(qū)中。風(fēng)電場(chǎng)接入原仿真系統(tǒng)發(fā)電機(jī)節(jié)點(diǎn)位置的對(duì)應(yīng)關(guān)系如表1所示。
表1 風(fēng)電場(chǎng)接入系統(tǒng)位置對(duì)應(yīng)關(guān)系Tab.1 Position correspondence of wind farm access system
另外,電力系統(tǒng)在調(diào)度當(dāng)日的最大負(fù)荷為6000MW,電網(wǎng)負(fù)荷預(yù)測(cè)偏差率α為5%。如圖4所示,為部分風(fēng)力發(fā)電場(chǎng)的預(yù)測(cè)出力曲線??梢钥吹剑诟呔S新能源接入系統(tǒng)的情景下,其功率輸出的相關(guān)關(guān)系復(fù)雜性較高。
圖4 部分風(fēng)電場(chǎng)的預(yù)測(cè)出力曲線Fig.4 Predictive outputs curves of partial wind farms
4.1.1 高維新能源的出力相關(guān)性建模
基于1.2節(jié)中的建模流程,取采樣時(shí)長(zhǎng)為1年、采樣間隔為15分鐘,以獲得10個(gè)風(fēng)電場(chǎng)的功率輸出歷史數(shù)據(jù),從而建立風(fēng)電場(chǎng)之間的出力相關(guān)性模型。如圖5所示,即為風(fēng)電場(chǎng)歷史出力數(shù)據(jù)間的線性相關(guān)系數(shù)熱度圖。
圖5 風(fēng)電場(chǎng)出力線性相關(guān)系數(shù)Fig.5Linear correlation coefficients of wind power outputs
由圖5中相關(guān)性熱度圖的整體情況可見(jiàn),風(fēng)電場(chǎng)2的出力與其余風(fēng)電場(chǎng)出力間的線性相關(guān)系數(shù)較高,也即相關(guān)性更強(qiáng)。因此,選擇建模根節(jié)點(diǎn)為風(fēng)電場(chǎng)2的出力序列,基于動(dòng)態(tài)C藤結(jié)構(gòu)進(jìn)行pair-Copula建模。如表2所示,即為基于動(dòng)態(tài)C藤結(jié)構(gòu)建模所得到的風(fēng)電場(chǎng)出力模型最優(yōu)Copula函數(shù)結(jié)果。
表2 風(fēng)電場(chǎng)出力相關(guān)性模型的最優(yōu)Copula函數(shù)Tab.2 Optimal Copula function of wind farm output correlation model
如圖6所示,為時(shí)變t-Copula函數(shù)C36|297101458的參數(shù)變化過(guò)程。由圖6可見(jiàn),所建模型充分利用了動(dòng)態(tài)Copula函數(shù)中相關(guān)參數(shù)的時(shí)變性,因此能更完備地描述C藤結(jié)構(gòu)的動(dòng)態(tài)特征。
圖6 時(shí)變t-Copula的參數(shù)變化圖Fig.6 Parametervariation diagram of time varying t-Copula
為驗(yàn)證所建模型的有效性,計(jì)算10維風(fēng)電場(chǎng)出力的靜態(tài)pair-Copula以及動(dòng)態(tài)pair-Copula模型的歐氏距離和AIC[29-31]值,如表3所示。同時(shí),模型的K(z)檢驗(yàn)結(jié)果由圖7展示。
表3 不同pair-Copula模型的歐氏距離和AIC值Tab.3 Euclidean distance and AIC value of different pair-Copula models
圖7 不同pair-Copula模型的K(z)檢驗(yàn)Fig.6 K(z) test of different pair-Copula models
查看表3與圖7的結(jié)果可見(jiàn),本文所提模型對(duì)比采取靜態(tài)Copula函數(shù)的pair-Copula模型,其 AIC值以及與經(jīng)驗(yàn)Copula分布的歐氏距離更小。因此,基于動(dòng)態(tài)C藤pair-Copula函數(shù)的本文模型擬合效果更佳,同時(shí)對(duì)高維新能源實(shí)際出力概率分布的表達(dá)也更為精準(zhǔn)。另外,IEEE118節(jié)點(diǎn)系統(tǒng)不同區(qū)域的常規(guī)機(jī)組出力情況與全體風(fēng)電機(jī)組各時(shí)段的出力曲線,如圖8所示。
圖8 IEEE118節(jié)點(diǎn)系統(tǒng)的機(jī)組出力曲線Fig.8 Unit power output curves of IEEE 118 node system
4.1.2 高維新能源的動(dòng)態(tài)相關(guān)關(guān)系影響分析
基于是否考慮高維新能源的動(dòng)態(tài)相關(guān)關(guān)系,分析其對(duì)于調(diào)度結(jié)果帶來(lái)的影響。采用靜態(tài)C藤 pair-Copula模型和不考慮風(fēng)電場(chǎng)出力相關(guān)性的模型作為對(duì)比模型。然后,每個(gè)模型基于抽樣技術(shù)生成10個(gè)隨機(jī)誤差場(chǎng)景。如表4所示,即為利用文中第2小節(jié)提出的MCO方法求解三種模型所得的結(jié)果。
表4 三種模型的結(jié)果對(duì)比Tab.4 Result comparison between three methods
由表4可知,對(duì)比不考慮高維新能源相關(guān)性以及只考慮靜態(tài)C藤pair-Copula函數(shù)的情況,計(jì)及高維動(dòng)態(tài)C藤pair-Copula函數(shù)的參數(shù)所具備的時(shí)變特性進(jìn)行建模的系統(tǒng)總運(yùn)行成本最低,經(jīng)濟(jì)調(diào)度結(jié)果最優(yōu)。因此,考慮高維動(dòng)態(tài)C藤pair-Copula函數(shù)進(jìn)行建模,能夠使得所建經(jīng)濟(jì)調(diào)度模型符合實(shí)際環(huán)境,從而將為含有高維新能源的電力系統(tǒng)動(dòng)態(tài)經(jīng)濟(jì)調(diào)度模型提供更為精確的設(shè)計(jì)依據(jù)。
本文基于動(dòng)態(tài)Copula函數(shù)中參數(shù)所具備的時(shí)變特性,以及C藤pair-Copula函數(shù)的核心框架,構(gòu)建了高維新能源出力的概率模型,該模型計(jì)及高維新能源出力之間的相關(guān)特性,對(duì)高維相依新能源結(jié)構(gòu)的整體與局部特性進(jìn)行了準(zhǔn)確地刻畫。其次,本文立足于多學(xué)科協(xié)同優(yōu)化模型的思想,結(jié)合場(chǎng)景分析法求解含高維新能源接入電力系統(tǒng)的動(dòng)態(tài)經(jīng)濟(jì)調(diào)度問(wèn)題,不僅能夠滿足系統(tǒng)級(jí)優(yōu)化模型在每次迭代中可行域存在的要求,而且可以反映出系統(tǒng)級(jí)問(wèn)題最優(yōu)解與誤差場(chǎng)景級(jí)問(wèn)題最優(yōu)解間隱含的耦合關(guān)系。最后,基于IEEE118節(jié)點(diǎn)系統(tǒng)仿真驗(yàn)證結(jié)果表明:動(dòng)態(tài)C藤pair-Copula模型能夠完備地描述出高維相依新能源的實(shí)際結(jié)構(gòu)特征,同時(shí)文中應(yīng)用的基于場(chǎng)景分析的電網(wǎng)動(dòng)態(tài)經(jīng)濟(jì)調(diào)度多學(xué)科協(xié)同優(yōu)化算法具備降維求解能力強(qiáng)的優(yōu)勢(shì)。