馬騫，鄧卓明，吳云亮，李鵬，姚文峰，郭知非

(1.中國南方電網(wǎng)電力調(diào)度控制中心，廣州510663；2. 南方電網(wǎng)科學(xué)研究院，廣州510663)

0 引言

受端電網(wǎng)的電能往往由其外部遠距離的發(fā)電廠通過高壓輸電線傳輸獲取，并且對應(yīng)著密集的工商業(yè)負(fù)荷，換流站、輸電線路以及異步電動機、電弧爐、降溫壓縮電機等負(fù)荷的無功功率需求量巨大[1]。在負(fù)荷側(cè)安裝一定數(shù)量和容量的靜止同步補償器(static synchronous compensator, STATCOM)做無功補償，不僅能在電網(wǎng)穩(wěn)定運行時維持電壓穩(wěn)定、提高供電質(zhì)量、降低電能損耗，同時在發(fā)生較為嚴(yán)重的短路故障時，能提供堅強的無功支撐，快速恢復(fù)電壓至安全水平，避免甩負(fù)荷。目前在經(jīng)濟條件相對較好的廣東電網(wǎng)，已經(jīng)逐漸在負(fù)荷端如東莞、水鄉(xiāng)配置了STATCOM[2]。STATCOM的規(guī)劃和控制也獲得了越來越多的關(guān)注。文獻[3]提出了在多直流饋入系統(tǒng)各母線對多回直流逆變站電壓支撐強度的指標(biāo)，并用于南方電網(wǎng)STATCOM配置方案的研究。文獻[4]描述并復(fù)現(xiàn)了HVDC系統(tǒng)在故障恢復(fù)過程中發(fā)生“二次換相失敗”的現(xiàn)象，提出了根據(jù)故障的嚴(yán)重程度差異化配置STATCOM無功補償模式。STATCOM作為一種新興的電力電子無功補償裝置，盡管有著極其優(yōu)秀的無功補償性能(反應(yīng)時間短，補償容量靈活多變)，但是其造價仍然高昂，規(guī)劃配置STATCOM必須精打細(xì)算，避免電網(wǎng)投資的浪費[5]。

目前電網(wǎng)規(guī)劃人員對STATCOM的規(guī)劃配置大多采用枚舉試錯的方式，通過大量的仿真試驗獲得滿足安全約束的規(guī)劃方案。然而，大量的仿真試錯不僅工作量巨大、耗費大量的時間，而且獲得的規(guī)劃方案也僅僅是一個相對較好、能滿足安全約束的方案，并不是最經(jīng)濟的方案。并且現(xiàn)代電力系統(tǒng)模型越來越復(fù)雜，規(guī)模也越來越大，組合爆炸的問題將導(dǎo)致枚舉無法進行。

因此，國內(nèi)外一些文獻提出了數(shù)學(xué)化的規(guī)劃方案。在規(guī)劃STATCOM的配置時，需要考慮STATCOM的安裝地點以及安裝容量，即對應(yīng)定點和定容的規(guī)劃。對此國內(nèi)外已有相應(yīng)的研究成果。文獻[6 - 7]采用基于動態(tài)無功—電壓靈敏度的電壓控制敏感因子和多饋入相互作用因子，計算獲得電壓控制敏感點，以便確定動態(tài)無功源的布點。在此基礎(chǔ)獲得候選節(jié)點后，文獻[8]進一步采用基于分解的多目標(biāo)進化算法確定動態(tài)無功源的容量。文獻[9]指定動態(tài)無功—電壓靈敏度最大的節(jié)點為候選節(jié)點，然后指定STATCOM的容量為對應(yīng)節(jié)點補償?shù)纳舷拗?。文獻[10]利用軌跡靈敏度計算和奇異值分解確定動態(tài)無功源的布點，然后通過時域仿真和非線性優(yōu)化之間的交替迭代獲得規(guī)劃容量。盡管上述的方法能夠獲得相對較好的、工程上能用的規(guī)劃方案，然而，因為它們把動態(tài)無功源布點及容量完全割裂開來，獲得的規(guī)劃方案只是布點和容量二層的分別最優(yōu)，并不是它們整體的最優(yōu)方案。

為此，本文分別采用整數(shù)變量和連續(xù)變量表征STATCOM的安裝地點和安裝容量，把整體的規(guī)劃問題解析為數(shù)學(xué)上的混合整數(shù)規(guī)劃模型，以便獲得布點和容量規(guī)劃的整體最優(yōu)方案。為了提高求解效率，本文采取了凸松弛技術(shù)[11 - 13]，借助STATCOM布點的某些特征，關(guān)于整數(shù)變量整體凸化原始的混合整數(shù)規(guī)劃模型，然后對凸模型實施一種緊松弛。經(jīng)算例驗算，本文的混合整數(shù)規(guī)劃模型和凸松弛技術(shù)能夠獲得STATCOM定點和容量的最優(yōu)規(guī)劃方案，有效緩解電網(wǎng)短期電壓失穩(wěn)的問題。同時，相比割裂布點和容量的傳統(tǒng)規(guī)劃方法，本文的方法獲得的規(guī)劃方案更優(yōu)，即規(guī)劃的費用更小。本文的創(chuàng)新在于提出一個通用的混合整數(shù)規(guī)劃模型和凸松弛求解方案，避免大規(guī)模枚舉和試錯仿真工作，并且獲得最經(jīng)濟的規(guī)劃方案。同時，得益于電網(wǎng)元器件的模型的通用性，本文提出的優(yōu)化模型可以實施封裝，后續(xù)可以根據(jù)用戶需要對模型進行修改，從而推廣到其他的電網(wǎng)。該模型的通用性將可以節(jié)省大量仿真試錯的工作量和時間。

1 STATCOM規(guī)劃的混合整數(shù)規(guī)劃模型

1.1 目標(biāo)函數(shù)

本文目標(biāo)函數(shù)J(·)定義為STATCOM規(guī)劃方案的投資費用，包括與基建相關(guān)的固定安裝費用(例如征地補償)以及與容量相關(guān)的容量費用(例如設(shè)備購置)2個部分。STATCOM總體的規(guī)劃投資費用由固定安裝費用和容量費用疊加[14]。

J(·)=vT[c1e+c2·(QCmax-QLmax)]

(1)

式中：v為STATCOM安裝地點的二進制列向量(即所有元素皆為二進制數(shù)值)，其中第i個元素的數(shù)值為0，代表候選節(jié)點i不安裝STATCOM，無需記及固定安裝費用和容量費用，第i個元素的數(shù)值為1，代表候選節(jié)點i安裝STATCOM，需要記及固定安裝費用和容量費用；c1為固定安裝費用，萬元；c2單位容量費用，萬元/Mvar；e為單位列向量；QCmax=[qc1,qc2,qc3,…,qci]為STATCOM的容性輸出限值列向量；QLmax=[ql1,ql2,ql3, …,qli]為STATCOM的感性輸出限值列向量；其中第i個元素qci和qli分別對應(yīng)候選節(jié)點i的容性容量(正值)和感性容量限值(負(fù)值)。需要注意的是，QCmax和QLmax均為可參與優(yōu)化的變量，而不是固定不變的參數(shù)。

在極小化目標(biāo)函數(shù)(即規(guī)劃投資費用)的過程中，各個站址STATCOM的配置容量(即QCmax-QLmax)會逐漸下降直至收斂，極小化目標(biāo)函數(shù)能夠限值各個站址的配置容量。

1.2 電壓和功角安全約束

根據(jù)南方電網(wǎng)安全穩(wěn)定導(dǎo)則，電壓安全約束為故障清除1 s時間內(nèi)，節(jié)點電壓恢復(fù)至0.75 p.u.以上。

Vl.k(t)≥0.75,t∈[tc+Td,tf]

(2)

式中：l和k分別為節(jié)點編號和故障編號；t為仿真時刻；tc和tf分別為故障清除時刻和仿真終止時刻；Td為持續(xù)時間1 s；Vl.k(t)為節(jié)點電壓。

功角安全約束則必須保證所有發(fā)電機功角與慣量中心轉(zhuǎn)子角的功角差小于給定的數(shù)值120 °。

|δg.k(t)-δCOI.k(t)|≤δmax

(3)

式中：g為發(fā)電機編號；δg.k(t)為發(fā)電機功角；δmax為給定的閾值120 °；δCOI.k(t)為慣量中心轉(zhuǎn)子角。

(4)

式中：Mg為發(fā)電機慣量；NG為發(fā)電機數(shù)量。

此外，約束條件還應(yīng)包括動態(tài)設(shè)備的微分代數(shù)方程、系統(tǒng)的潮流方程、變量上下限約束以及邊界約束等等。因為這些方程都比較常規(guī)且繁瑣，此處不一一列舉，具體的方程可參考文獻[9,15]。其中，發(fā)電機采用經(jīng)典四階模型，勵磁系統(tǒng)采用經(jīng)典四階模型，異步電機負(fù)荷采用經(jīng)典三階模型，STATCOM采用增益-時間-常數(shù)模型。

1.3 STATCOM容量約束

在實際工程中，各個站址由于場地、檢修等限制，可能存在數(shù)值不相等的受限容量。該限制可以通過控制變量QCmax和QLmax的上下限約束(不等式約束)表達。

(5)

1.4 混合整數(shù)規(guī)劃模型

STATCOM規(guī)劃問題的混合整數(shù)規(guī)劃模型可采用簡約形式表示，具體如式(6a)—(6e)所示。

(6a)

(6b)

0=Gk(xk(t),yk(t),u,v)

(6c)

0≤Hk(xk(t),yk(t),u)

(6d)

v∈{0, 1},t∈[t0,tf]

(6e)

式中：xk(t)、yk(t)、u、v分別為狀態(tài)變量、代數(shù)變量、連續(xù)決策變量(即規(guī)劃安裝容量)以及二進制決策變量(即規(guī)劃安裝布點)的列向量；t0為仿真起始時刻；在方程方面，式(6b) 為微分方程，包含發(fā)電機及其勵磁系統(tǒng)、異步電機負(fù)荷以及STATCOM的動態(tài)；式(6c)為受端電網(wǎng)的潮流方程，其中STATCOM以注入功率的形式掛在對應(yīng)節(jié)點上；式(6d)包含所有變量(狀態(tài)變量、代數(shù)變量和決策變量)的上下限約束以及1.2節(jié)提到的電壓和功角安全約束。關(guān)于式(6b)提到的設(shè)備模型和微分方程可參考文獻[9,15]。

2 混合整數(shù)規(guī)劃模型的凸松弛

凸松弛方法是一種針對和處理混合整數(shù)規(guī)劃問題的處理方案。經(jīng)過凸松弛處理，混合整數(shù)規(guī)劃的求解效率可以極大提高。凸松弛包含2個步驟：凸化以及松弛。

凸化過程需要列舉二進制決策變量的所有可能組合，然后混合整數(shù)規(guī)劃模型整體關(guān)于二進制決策變量做線性疊加。

線性疊加可以使得規(guī)劃模型關(guān)于二進制變量整體凸化。因為STATCOM以注入功率掛在潮流方程對應(yīng)節(jié)點，該功率存在與否是由決策變量v決定的。因此潮流方程(6c)是直接與v相關(guān)的。此外，在1.1節(jié)介紹過目標(biāo)函數(shù)也是與v直接相關(guān)的。除此以外，微分方程、上下限約束和安全約束都沒有與v直接相關(guān)。因此只需對式(6a)和(6c)做關(guān)于二進制變量的線性疊加即可。

本文借用一個虛擬4節(jié)點網(wǎng)絡(luò)為例子詳細(xì)介紹凸化。假定在該虛擬的4節(jié)點網(wǎng)絡(luò)中，節(jié)點1—3為負(fù)荷節(jié)點，并被選為STATCOM的候選節(jié)點，節(jié)點4是平衡節(jié)點。因為在此案例共3個候選節(jié)點，每個節(jié)點安裝或不安裝(1/0)共2種情況，共有8種布點方案，具體如表1所示。

表1 虛擬案例的8種布點方案Tab.1 Eight plans of locations in the virtual case

為表1中每種組合引入替代變量αw。 依據(jù)線性疊加的計算，節(jié)點1—3的潮流方程可分別凸化為：

(7a)

(7b)

(7c)

式中Gl.k(xl.k,yl.k,ul,vl.w)和vl.w分別為節(jié)點l的潮流方程函數(shù)以及vw的具體值。

由表1和式(7a)—(7c)可知，因為整數(shù)決策變量共8種組合，凸化需要枚舉8次，引入8個替代變量，才能完成線性疊加。然而，α1-α8實際并非獨立變量，可以做一定的簡化。在對應(yīng)的潮流方程上分別引入變量θ1.1=α1+α2+α3+α5，θ1.2=α4+α6+α7+α8，θ2.1=α1+α2+α4+α6，θ2.2=α3+α5+α7+α8，θ3.1=α1+α3+α4+α7， 以及θ3.2=α2+α5+α6+α8， 其中θl.i代表節(jié)點l的第i個θ變量。式(7a)—(7c)可簡化為：

0=G1.k(x1.k,y1.k,u1, 0)θ1.1+G1.k(x1.k,y1.k,u1, 1)θ1.2

(8a)

0=G2.k(x2.k,y2.k,u2, 0)θ2.1+G2.k(x2.k,y2.k,u2, 1)θ2.2

(8b)

0=G3.k(x3.k,y3.k,u3, 0)θ3.1+G3.k(x3.k,y3.k,u3, 1)θ3.2

(8c)

并且同樣滿足θ1.1+θ1.2=1，θ2.1+θ2.2=1以及θ3.1+θ3.2=1的約束條件，且θl.i為二進制變量。這是因為STATCOM以注入功率掛在負(fù)荷節(jié)點上，節(jié)點l安裝STATCOM與否只影響節(jié)點l列寫潮流方程，并不影響其他節(jié)點列寫相應(yīng)的潮流方程，因此節(jié)點間列寫潮流方程是相互獨立的。

同理，對于目標(biāo)函數(shù)J(·)， 假定Jl(ul,vl.w)是節(jié)點l規(guī)劃安裝STATCOM的費用，則目標(biāo)函數(shù)可凸化如式(9)所示。

(9)

松弛過程則把上述的替代變量θl.i從二進制的取值范圍松弛至連續(xù)區(qū)間[0,1]，混合整數(shù)規(guī)劃模型整體轉(zhuǎn)換為連續(xù)規(guī)劃問題，從而避免了直接求解混合整數(shù)規(guī)劃問題。在優(yōu)化領(lǐng)域，基于梯度的算法或求解器都能比較高效地求解連續(xù)規(guī)劃問題，因此提高了STATCOM規(guī)劃的效率。同時，經(jīng)文獻[11 - 13]的驗證，上述關(guān)于整數(shù)變量的凸化及松弛是一個緊松弛過程，可以獲得滿足邊界條件約束的規(guī)劃方案。

3 整體規(guī)劃流程

依據(jù)上述的混合整數(shù)規(guī)劃模型以及凸松弛方法，本文歸納了STATCOM的整體規(guī)劃流程如下。

1)采用三相接地短路故障N-2掃描受端電網(wǎng)薄弱區(qū)域，獲得電網(wǎng)在規(guī)劃STATCOM前的短期電壓失穩(wěn)的故障集；

2)以安裝費用最小為目標(biāo)函數(shù)，建立全網(wǎng)電壓安全約束的STATCOM布點及容量規(guī)劃的規(guī)劃模型式(6a)—(6e)；

3)枚舉二進制布點變量v的所有組合，采用固定的列向量vw(w= 1, 2,…,NW)代表對應(yīng)的每1種規(guī)劃方案；

5)關(guān)于二進制變量，整體做線性疊加以凸化混合整數(shù)規(guī)劃模型；

6)松弛替代變量θl.i至連續(xù)區(qū)間[0,1]；

7)在通用代數(shù)建模系統(tǒng)(general algebraic modeling system，GAMS)上調(diào)用CONOPT求解器[16]求解經(jīng)過凸松弛處理后的連續(xù)優(yōu)化模型，獲得STATCOM布點及容量規(guī)劃方案；

8)采用PSD-BPA軟件仿真校驗規(guī)劃方案是否滿足暫態(tài)電壓、功角穩(wěn)定約束。

4 算例分析

本文在2個系統(tǒng)上實施了上述的規(guī)劃方法。系統(tǒng)A是IEEE標(biāo)準(zhǔn)39節(jié)點測試網(wǎng)絡(luò)，包含有10發(fā)電機、39節(jié)點以及17異步電機。系統(tǒng)B是廣東省級電網(wǎng)，包含有181發(fā)電機、1 009節(jié)點以及361異步電機。本章的算例計算是在一臺高性能服務(wù)器上完成的，它具有16組Dell M620刀片機，每個刀片又配有2個Xeon E5- 2650處理器(其中每個處理器都是8核8線程)和128 GB的運行內(nèi)存。計算平臺分別采用GAMS 24.4.6做優(yōu)化計算[16]，MATLAB 2015a 做編程，以及PSD-BPA做時域仿真[17]。在上述2個測試系統(tǒng)上進行N-2掃描，分別獲得6個短期電壓失穩(wěn)的故障組成故障集。以下STATCOM的布點和容量規(guī)劃將解決故障集導(dǎo)致的短期電壓失穩(wěn)問題。

首先以系統(tǒng)A為例子，把本文的規(guī)劃方法對比傳統(tǒng)的布點及容量完全割裂的規(guī)劃方法，即先通過靈敏度計算確定安裝地點再確定安裝容量。經(jīng)過計算，表2列舉了本文的布點和容量統(tǒng)一規(guī)劃方案，共需6臺STATCOM，配置總?cè)萘繛?86.52 Mvar。表3列舉了傳統(tǒng)的布點及容量完全割裂的規(guī)劃方案，同樣需要6臺STATCOM共338.68 Mvar。由2個表格可看出，雖然兩個方案需要的STATCOM數(shù)量即布點數(shù)相同，但本文的布點和容量統(tǒng)一規(guī)劃方案需要的容量和費用遠小于傳統(tǒng)的規(guī)劃方法，因此更加經(jīng)濟，是安裝地點和安裝容量整體的最優(yōu)結(jié)果；而傳統(tǒng)的做法僅僅只是安裝地點和安裝容量的分別最優(yōu)。

表2 本文的布點和容量統(tǒng)一規(guī)劃方案Tab.2 Unified plan of location and capacity in this paper Mvar

STATCOM的規(guī)劃難點在于混合整數(shù)規(guī)劃。有學(xué)者指出，如果全體負(fù)荷母線均為STATCOM備選安裝點，只規(guī)劃STATCOM容量，即可避免混合整數(shù)規(guī)劃問題。然而，該方法只能最小化容量的規(guī)劃，卻無法記及STATCOM的固定安裝費用c1。STATCOM總體投資費用變得很大，因為幾乎所有母線都規(guī)劃安裝STATCOM。

同樣以系統(tǒng)A為例，如果所有負(fù)荷點均為STATCOM的規(guī)劃節(jié)點，以此極小化STATCOM的規(guī)劃容量。該規(guī)劃方案列舉在圖1中，在17個負(fù)荷

表3 傳統(tǒng)的布點及容量完全割裂的規(guī)劃方案Tab.3 Conventional separated plan of location and capacity Mvar

節(jié)點的規(guī)劃容量共219.28 Mvar，略小于表2即本文方案的286.52 Mvar。但是，圖1的規(guī)劃方案中，規(guī)劃容量被分散到17臺STATCOM，數(shù)量較大，因此基建費用較高。對應(yīng)的STATCOM規(guī)劃費用是32 716.52 萬元，遠高于表2的方案20 691.087 萬元。因此在實際運行中，南方電網(wǎng)較多情況是集中安裝少量的大容量STATCOM以減少基建費用，提高投資的效率。

圖1 STATCOM在所有負(fù)荷節(jié)點的規(guī)劃容量Fig.1 Capacities of STATCOM at all of the loads

為了驗證本文凸松弛方法的正確性，系統(tǒng)A作為公共可用的IEEE測試系統(tǒng)，其算例用于案例的比較和復(fù)現(xiàn)。傳統(tǒng)的布點及容量完全割裂的規(guī)劃方法或者標(biāo)準(zhǔn)求解器，在涉及大型實際電網(wǎng)案例時會面臨巨大的求解壓力；相反，它們可以快捷地求解系統(tǒng)A類的測試系統(tǒng)，并提供相應(yīng)的計算結(jié)果，作為凸松弛方法的對比。在另一方面，系統(tǒng)B作為一個實際省級電網(wǎng)，其算例可驗證凸松弛方法求解STATCOM布點及容量規(guī)劃問題的高效性。在此算例上，凸松弛方法在可接受的時間框架內(nèi)可獲取STATCOM規(guī)劃問題的方案，然而像SBB這種求解器卻不能完成，足以驗證凸松弛方法的高效性。

如圖2所示，系統(tǒng)B在故障集下節(jié)點電壓恢復(fù)緩慢，大量節(jié)點在故障清除1 s后仍然不能恢復(fù)0.75 p.u.以上，短期電壓失穩(wěn)。采用本文凸松弛算法計算獲得的STATCOM規(guī)劃方案如表4所示。在該規(guī)劃方案下，如圖3所示，系統(tǒng)B能夠恢復(fù)短期電壓安全。

圖2 系統(tǒng)B在故障集下的電壓軌跡Fig.2 Voltages of system B under contingency set

表4 系統(tǒng)B的STATCOM規(guī)劃方案Tab.4 Planning of STATCOM for system B

圖3 含STATCOM規(guī)劃的系統(tǒng)B在故障集下的電壓Fig.3 Voltages of system B under contingency set with planning of STATCOM

針對混合整數(shù)規(guī)劃，目前存在成熟的商業(yè)求解器SBB。SBB求解器是GAMS平臺上其中一種較為流行和通用的混合整數(shù)非線性規(guī)劃求解器。SBB結(jié)合了求解混合整數(shù)線性規(guī)劃的分支定界(branch and bound，B&B)方法以及某些標(biāo)準(zhǔn)的非線性規(guī)劃求解器(可供用戶選擇)如MINOS、SNOPT。在求解混合整數(shù)非線性規(guī)劃問題時，SBB先對整數(shù)規(guī)劃部分實施分支定界，求取一個相對較好的整數(shù)規(guī)劃結(jié)果，然后再基于用戶選擇的非線性規(guī)劃求解器對非線性規(guī)劃部分實施求解，從而獲得混合整數(shù)非線性規(guī)劃的結(jié)果。筆者在從事混合整數(shù)非線性規(guī)劃問題研究的時候，先后比較過BARON、DICOPT、SBB等混合整數(shù)非線性規(guī)劃求解器，最后發(fā)現(xiàn)SBB求解器對于本文的STATCOM布點及容量規(guī)劃問題具有較好的求解效果。因此，文章最后將本文的凸松弛方法和SBB求解器的結(jié)果進行對比，以便驗證本文方法的有效性和高效性。如表5展示了本文規(guī)劃方法與SBB求解器的結(jié)果對比。

表5 本文方法和SBB求解器計算結(jié)果的比較Tab.5 Comparison of the results between the method in this paper and the SBB solver

在規(guī)模比較小的系統(tǒng)A中，凸松弛方法與SBB求解器的STATCOM整體規(guī)劃費用基本一致，相差無幾，因此本文的規(guī)劃方案是可靠的。并且，本文凸松弛方法的計算時間遠短于SBB求解器，說明凸松弛方法的求解效率更高。面對大規(guī)模系統(tǒng)B的STATCOM規(guī)劃，SBB求解器需要求解大型混合整數(shù)規(guī)劃問題，在花費較長時間后仍然不能獲得STATCOM的規(guī)劃方案。整數(shù)規(guī)劃問題可以理解為組合優(yōu)化問題，而組合優(yōu)化問題的規(guī)模隨著控制變量數(shù)量的增加呈現(xiàn)指數(shù)增長，極其容易引發(fā)數(shù)學(xué)上的指數(shù)爆炸問題。對于實際1 009節(jié)點系統(tǒng)B的規(guī)劃問題，由于節(jié)點數(shù)量龐大，其控制變量(包括STATCOM的規(guī)劃容量和安裝地點)的數(shù)量也相應(yīng)非常大。一般而言，當(dāng)控制變量的數(shù)量大于20，組合的數(shù)量已經(jīng)超過百萬。如此龐大的組合數(shù)量，在上述分支定界的過程要存儲很多分支和葉子節(jié)點，花費很多內(nèi)存空間。分支定界法的求解效率基本上由值界方法決定，若界估計不好，在極端情況下將與窮舉搜索沒多大區(qū)別，在有限的時間內(nèi)不能計算所有的組合。因此，面對1 009節(jié)點系統(tǒng)B，SBB求解器不能獲得最優(yōu)解，其本質(zhì)原因是變量的數(shù)量太大，同時組合數(shù)量指數(shù)式增長，出現(xiàn)了整數(shù)規(guī)劃的指數(shù)爆炸問題。但凸松弛方法不需要直接求解混合整數(shù)優(yōu)化，只需要求解連續(xù)的優(yōu)化問題。許多經(jīng)典的基于梯度的優(yōu)化方法，如內(nèi)點法，可以輕松獲得最優(yōu)解，因此凸松弛方法在系統(tǒng)B仍然得到了最優(yōu)的規(guī)劃方案。因此，本文的規(guī)劃方法更適用于南方電網(wǎng)的STATCOM規(guī)劃。

5 結(jié)語

本文針對受端電網(wǎng)STATCOM的布點和容量規(guī)劃問題，采用了布點和容量統(tǒng)一規(guī)劃的混合整數(shù)規(guī)劃模型，并提出了一種針對混合整數(shù)規(guī)劃的凸松弛方法。該方法首先關(guān)于整數(shù)變量整體凸化混合整數(shù)規(guī)劃模型，然后對凸模型實施一種緊松弛，避免直接求解混合整數(shù)規(guī)劃。

在IEEE標(biāo)準(zhǔn)39節(jié)點以及廣東電網(wǎng)的算例計算表明，該方法比傳統(tǒng)的割裂布點和容量的規(guī)劃方法的STATCOM投資費用更小。布點和容量同時規(guī)劃的必要性得到了驗證。相比商業(yè)化的求解器SBB，本文凸松弛方法更適用于在大型受端電網(wǎng)如廣東電網(wǎng)的STATCOM規(guī)劃問題，因求解效率更高，計算時間更短，容易獲得大規(guī)模規(guī)劃問題的最優(yōu)解。