洪 雯，沈 青，魯愛斌，鄭 雪

(武漢電力職業(yè)技術學院，湖北 武漢 430072)

第二次工業(yè)革命后，電能使用愈發(fā)廣泛，電力極大促進生產(chǎn)力的快速發(fā)展，推動了社會的進步。電力系統(tǒng)在確保安全穩(wěn)定運行的前提下，最主要任務就是盡可能經(jīng)濟、有效地滿足用戶的需求，保證持續(xù)不間斷的供電。而電力系統(tǒng)重要環(huán)節(jié)之一的輸電線路經(jīng)?？缭缴交鹨装l(fā)地段，若山火發(fā)生在輸電線路附近，山火產(chǎn)生的污穢容易引起空氣放電，會導致多條輸電線路跳閘，此外跳閘線路的潮流將轉(zhuǎn)移到其它線路，可能造成其它線路潮流接連越限，進而發(fā)生連鎖故障，導致整個系統(tǒng)崩潰，造成長時間的大面積停電事故。

安全校正控制的功能是采取一系列合理措施，在故障初期在短時間內(nèi)進盡可能解決電網(wǎng)中的各種潮流越限情況，迅速將全網(wǎng)潮流控制在安全范圍內(nèi)，以防止連鎖性的事故出現(xiàn)和惡化，將潮流越限的波及范圍和經(jīng)濟損失控制在最低程度。因此，山火導致多條線路跳閘后，電網(wǎng)能否安全穩(wěn)定運行，主要取決于安全校正策略是否正確、合理。

目前，已有的安全校正策略主要是靈敏度法[1-3]和優(yōu)化規(guī)劃法[4-5]。靈敏度法通過計算出控制變量對支路潮流變化的關系，然后根據(jù)支路潮流的越限量反推出控制變量的調(diào)整量來達到消除支路潮流越限的目的，便于操作實施，但是無法處理系統(tǒng)運行中的一些安全約束問題，容易造成原本處于平衡狀態(tài)的機組出現(xiàn)潮流越限的狀況，也可能在調(diào)整一些支路的潮流滿足約束條件后有其他的支路出現(xiàn)了潮流越限的情況。優(yōu)化規(guī)劃法通過求解數(shù)學模型得出控制方案，因為可以設置復雜的約束條件和合理的目標函數(shù)，故而有較好的安全性和經(jīng)濟性。

針對山火故障導致的多條輸電線路跳閘進而出現(xiàn)的潮流越限問題，本文采用優(yōu)化規(guī)劃類算法，根據(jù)當前系統(tǒng)狀態(tài)進行在線安全校正，以期在短時間內(nèi)恢復系統(tǒng)正常運行，避免連鎖故障的發(fā)生。

1 安全校正策略的優(yōu)化變量

本文選擇優(yōu)化規(guī)劃類方法制定安全校正策略。安全校正策略的優(yōu)化模型對電力系統(tǒng)進行了合理的簡化，將發(fā)電機和負荷的參數(shù)調(diào)整問題轉(zhuǎn)化為建立、求解非線性的連續(xù)和二元混合規(guī)劃問題。首先基于啟發(fā)式搜索原理確定優(yōu)化變量，然后以負荷總有功切除量最小為目標函數(shù)，將系統(tǒng)運行中的常規(guī)和頻率約束作為約束條件，建立優(yōu)化模型，并基于改進二元粒子群算法對該模型進行優(yōu)化計算，求解得出滿足約束條件且目標函數(shù)最小的安全校正方案。

當某電網(wǎng)出現(xiàn)支路潮流越限時，則需要進行發(fā)電機有功和無功功率調(diào)整和切負荷等操作來消除越限。對于實際電網(wǎng)，節(jié)點數(shù)眾多，若想逐一調(diào)整所有的發(fā)電機和負荷顯然是一件既費時又費力的事情。針對這種情況，本文采用啟發(fā)式搜索來確定優(yōu)化變量。

所謂啟發(fā)式搜索指的是在狀態(tài)空間中對系統(tǒng)中所有可能的優(yōu)化變量進行評估，得到優(yōu)化變量中最好的一部分，這部分最好的變量即為下一輪的搜索范圍，直到找到目標，通過這樣的方法，能避開大多數(shù)沒有意義的搜索過程，節(jié)約搜索時間，提高搜索速度。

本文基于負荷轉(zhuǎn)移因子，篩選出對潮流越限的支路、效果最好的發(fā)電機和負荷節(jié)點，將這些節(jié)點作為優(yōu)化變量，大大縮減了優(yōu)化變量的數(shù)量，縮短了優(yōu)化時間。

1.1 負荷轉(zhuǎn)移因子

負荷轉(zhuǎn)移系數(shù)是指發(fā)電機和負荷節(jié)點的網(wǎng)絡注入功率的變化引起支路潮流變化的比例系數(shù)。該系數(shù)越大，在該節(jié)點進行單位出力的調(diào)節(jié)，過載支路上潮流變化趨勢越大，這樣就能保證在調(diào)節(jié)量盡可能小的情況下達到所需的效果。

定義負荷轉(zhuǎn)移因子為母線注入電流變化時，電網(wǎng)中每一條線路上電流的變化，即：

(1)

式中：ΔIk表示母線注入電流的變化量；ΔIij表示母線k電流的變化在線路ij上所引起的電流變化量。

根據(jù)基爾霍夫定理1的推導，可得網(wǎng)絡方程可表示為U=ZI，其中，Z為阻抗矩陣，進一步得到母線電壓變化量與母線注入電流的關系：

(2)

進一步得到負荷轉(zhuǎn)移因子：

(3)

負荷轉(zhuǎn)移因子表示的是節(jié)點k負荷出現(xiàn)功率變化時, 在支路ij上所導致的功率變化。其中，Zik、Zjk為阻抗矩陣的元素；zij為線路ij的阻抗。

計算出每一個節(jié)點分別對每一條越限支路的負荷轉(zhuǎn)移因子后，從中篩選出對每一條越限支路影響較大的節(jié)點的集合，這些發(fā)電機節(jié)點即為對該條越限支路調(diào)節(jié)效果最好的一類發(fā)電機。

1.2 連續(xù)型變量

采用啟發(fā)式搜索的方法，基于負荷轉(zhuǎn)移系數(shù)得到對潮流越限支路潮流調(diào)節(jié)效果最大的部分發(fā)電機。由于發(fā)電機的有功功率是通過調(diào)整原動機出力、無功功率是通過調(diào)整勵磁參數(shù)實現(xiàn)的，故可以實現(xiàn)連續(xù)調(diào)節(jié)。

由于切除負荷只能有功和無功同時切除，如果切負荷后，不相應調(diào)整發(fā)電機的有功和無功功率，系統(tǒng)有功和無功就無法同時滿足平衡，可能出現(xiàn)頻率或者電壓失去穩(wěn)定的問題，因此將發(fā)電機的有功和無功均作為連續(xù)型優(yōu)化變量。

確定連續(xù)型變量的步驟如下：

1)首先基于1.1節(jié)的方法計算出所有發(fā)電機對待調(diào)節(jié)支路A的負荷轉(zhuǎn)移因子(Load Shift Factor, LSF)并從大到小排序；

2)然后篩選出對該待調(diào)節(jié)支路潮流調(diào)節(jié)效果最好的發(fā)電機節(jié)點xA；

3)重復步驟1)和步驟2)，依次得到對第i個待調(diào)節(jié)支路潮流調(diào)節(jié)效果最好的發(fā)電機節(jié)點xi；

4)對這些發(fā)電機節(jié)點取并集即得到對所有待調(diào)節(jié)支路調(diào)節(jié)效果最優(yōu)的發(fā)電機節(jié)點的集合，將這些發(fā)電機節(jié)點的有功、無功功率作為連續(xù)型變量x。

1.3 0-1二元變量

由于負荷無法實現(xiàn)連續(xù)調(diào)節(jié)，只能切或者不切，所以在切負荷控制中，對于負荷的控制是一個0-1優(yōu)化的問題。并且不一定每個負荷節(jié)點都只有一個負荷，對于有多個負荷的負荷節(jié)點，每一個負荷都是獨立的，應當分別考慮，而不是全部切除或保留。

確定0-1二元型變量的步驟如下：

1)首先基于1.1節(jié)的方法計算出所有負荷對待調(diào)節(jié)支路A的負荷轉(zhuǎn)移因子(Load Shift Factor, LSF)并從大到小排序；

2)然后篩選出對該待調(diào)節(jié)支路潮流調(diào)節(jié)效果最好的負荷節(jié)點yA；

3)重復步驟1)和步驟2)，依次得到對第i個待調(diào)節(jié)支路潮流調(diào)節(jié)效果最好的負荷節(jié)點yi；

4)對這些負荷節(jié)點取并集即得到對所有待調(diào)節(jié)支路調(diào)節(jié)效果最優(yōu)的負荷節(jié)點的集合yk；

5)由于文中負荷節(jié)點在實際中指變電站，而負荷指變電站的每條出線，因為變電站出線上的功率無法實現(xiàn)連續(xù)調(diào)節(jié)，故將變電站集合yk中每條出線的切除與否作為獨立的0-1二元型優(yōu)化變量y，1表示保留該條出線，0表示斷開該條出線。

1.4 優(yōu)化變量的初值

優(yōu)化變量初值的選取對整個優(yōu)化結果有較大的影響，一個好的初始值能讓優(yōu)化過程更為快捷，取得更好的優(yōu)化效果。

在優(yōu)化中，通常是隨機生成的一組數(shù)據(jù)，將該數(shù)據(jù)賦給變量作為初始值。但是在電力系統(tǒng)中，采用這種方法隨機性太強，如果某一次優(yōu)化的初始值不好，甚至無法得出可行解。故而在本文中，以系統(tǒng)故障前優(yōu)化變量的值作為初始值。經(jīng)過多次試驗證明，這樣不僅更快捷，并且得到的最終結果優(yōu)化變量的改變量也更小。

2 安全校正策略的目標函數(shù)

為了盡可能滿足對用戶的供電，保持供電的連續(xù)性，故而希望安全校正得出的方案中負荷有功總切除量最小。設置目標函數(shù)為：

式中：n為二元型變量個數(shù)；PLi為二元型變量i的有功功率；yi為二元型變量i的值，0表示該負荷被切除。

3 安全校正策略的約束條件

為了保證安全校正的安全性，故設置以下約束條件：

1)節(jié)點電壓約束：Vmin≤V1,2…n≤Vmax

2)發(fā)電機有功上下限約束：Pgmin≤Pg1,2…m≤Pgmax

3)發(fā)電機無功上下限約束：Qgmin≤Qg1,2…m≤Qgmax

4)負荷有功上下限約束：PLDmin≤PLD1,2…n≤PLDmax

5)支路潮流約束：PLine≤PLinemax

6)發(fā)電機視在功率約束：Sgmin≤Sg1,2…m≤Sgmax

7)功率因數(shù)約束即發(fā)電機的功率因數(shù)應該保持使發(fā)電機正常工作的一定范圍之內(nèi)，不得低于設定值：η1,2…m≥ηmin。

在本文中，以上約束的實現(xiàn)方法是通過在目標函數(shù)中加入懲罰函數(shù)部分。其原理是，通過在懲罰函數(shù)前加一個非常大的系數(shù)，該系數(shù)通常為目標函數(shù)的數(shù)量級倍數(shù)，可使懲罰函數(shù)的數(shù)量級遠遠大于目標函數(shù)本身的值。一旦發(fā)生越限，目標函數(shù)就會變的相當大，而安全校正方案的目標函數(shù)是負荷有功總切除量最小，這一組解自然就會被排除。

4 優(yōu)化模型的求解方法

在實際電力系統(tǒng)中，負荷是不可能實現(xiàn)連續(xù)調(diào)節(jié)的，所以模型采用改進二元粒子群優(yōu)化算法(Improved Binary Particle Swarm Optimizer, IBPSO)求解二元變量最優(yōu)規(guī)劃問題。

IBPSO基于兩個假設：

1)對于0-1二元變量y，由于優(yōu)化過程中y的極值點yopt未知，假定yopt取0或1的先驗概率均為0.5；

2)在尋找的過程中，個體極值Pbest和全局極值Gbest對最佳值的判定是獨立的。

記Pbest作出正確判斷的概率為pp，Gbest作出正確判斷的概率為pg，由于Pbest和Gbest是粒子群迭代過程的歷史最優(yōu)值，它們分別發(fā)現(xiàn)全局最優(yōu)解yopt的概率超過平均值，因此pp>0.5，pg>0.5。并且，為促使粒子群分散搜索，避免陷入局部最優(yōu)，要求pp>pg。取pg=0.8，pp=0.7。應用概率論中的條件概率Bayes公式，引入?yún)?shù)α和β，記：

(5)

通過式(5)計算得α=0.9，β=0.7，迭代次數(shù)為q+1時，二元粒子位置向量yq+1的0-1決策的依據(jù)如下。取r=rand，有：

模型基于二元粒子群優(yōu)化算法，通過調(diào)用電力系統(tǒng)分析綜合程序(Power System Analysis Software Package，簡稱PSASP)進行潮流計算來進行安全校正優(yōu)化模型的求解。安全校正優(yōu)化模型的整體求解流程圖如圖1所示。

圖1 安全校正優(yōu)化模型的整體求解流程圖

5 算理驗證

5.1 故障設置

采用華中某省電網(wǎng)的潮流數(shù)據(jù)進行算例驗證，對該電網(wǎng)應用上節(jié)的結構脆弱性計算方法，綜合2018年的該電網(wǎng)交流線路、母線、兩繞組變壓器和三繞組變壓器參數(shù)，和該電網(wǎng)網(wǎng)架數(shù)據(jù)，得到該電網(wǎng)網(wǎng)線路脆弱性排序結果。在得到的排序結果中篩選出500 kV等級的交流線路，根據(jù)數(shù)值大小進行排序，并節(jié)選了排序靠前的部分結果列入表1中。數(shù)值越接近1，表明該條線路退出對該電網(wǎng)網(wǎng)的結構穩(wěn)定影響越大；數(shù)值越接近0，則該條線路退出運行對該電網(wǎng)網(wǎng)結構穩(wěn)定性影響越小。

表1 某電網(wǎng)500 kV交流線路脆弱性排序表(節(jié)選)

山火雖有突發(fā)性，但均發(fā)生在植被茂密地區(qū)和干燥、高溫天氣。在上表中從上往下觀察，由于艷山紅500母線所處的地理位置在山區(qū)附近，所以本研究假設序號為2和3的兩條線路因山火跳閘。本研究假想的山火導致輸電線路跳閘的故障形式如表2所示。

5.2 優(yōu)化變量

由于不同的故障類型，需要優(yōu)化的優(yōu)化變量不同。所以針對本研究的假想故障，本節(jié)基于上節(jié)所述方法，選取該故障形式下對應的優(yōu)化變量。

5.2.1 待調(diào)節(jié)支路的選取

斷開上述兩條線路后，重新計算全網(wǎng)潮流。對于所有的約束條件，此時只有支路潮流越限，節(jié)點電壓未發(fā)生越限。其中發(fā)現(xiàn)，故障線路跳閘后，有4條線路的功率顯著增加，故選取這4條線路為待調(diào)節(jié)支路，即應該選取對這4條線路潮流影響較大的發(fā)電機和負荷作為優(yōu)化變量。這4條線路如表3所示。

表3 待調(diào)節(jié)支路表達式表 p.u.

5.2.2 優(yōu)化變量的選取

根據(jù)第1節(jié)所述方法得到優(yōu)化變量如表4所示，其中負荷轉(zhuǎn)移因子越接近1，表示影響程度越大。

表4 優(yōu)化變量表(節(jié)選)

5.2.3 優(yōu)化變量的初值

連續(xù)型優(yōu)化變量的初值(以及上下限)如表5所示。0-1二元型優(yōu)化變量的初值如表6所示。值得注意的是，文中負荷節(jié)點在實際中指變電站，而負荷指變電站的每條出線。假設每個變電站有4條出線，則相當于負荷節(jié)點的功率只能以1/4為為最小切除單位進行切除。

將以上發(fā)電機和負荷節(jié)點的的有功和無功功率作為優(yōu)化對象進行優(yōu)化。

表5 連續(xù)型優(yōu)化變量的初值表 p.u.

表6 0-1二元型優(yōu)化變量的初值表 p.u.

5.3 安全校正方案

基于優(yōu)化模型和求解方法，得出安全校正方案。表7列出了發(fā)電機的安全校正方案，表8列出了負荷的安全校正方案。

5.4 安全校正方案校核

根據(jù)優(yōu)化出的安全校正方案修改相應的發(fā)電機和負荷功率，重新計算全網(wǎng)潮流，得到安全校正方案對故障后越限值較大的4條線路的調(diào)整效果，也就是待調(diào)節(jié)支路優(yōu)化前后的有功值變化情況，如表9所示。

從表9中可以看出，按照優(yōu)化得到的安全校正方案調(diào)整后，該電網(wǎng)的潮流均降至安全穩(wěn)定范圍內(nèi)，不再有發(fā)生連鎖故障的隱患，安全校正方案的有效性得到驗證。

表7 發(fā)電機安全校正方案表 p.u.

表8 負荷安全校正方案表 p.u.

表9 待調(diào)節(jié)支路的優(yōu)化效果表

5.5 結果分析

通過對比表5和表7中優(yōu)化變量的優(yōu)化結果和原始數(shù)據(jù)，發(fā)現(xiàn)優(yōu)化后發(fā)電機的功率都有所減少，三個負荷節(jié)點中，湘黔東01的負荷全部切除，湘黔東02的負荷全部保留了。之所以優(yōu)化出的結果是切機和切負荷，是因為這些節(jié)點對待調(diào)節(jié)支路的調(diào)節(jié)效果均為正相關，即減小這些節(jié)點的功率，可以使得待調(diào)節(jié)支路的潮流相應減少。反之，若某些節(jié)點對待調(diào)節(jié)支路的調(diào)節(jié)效果為負相關，則可能需要通過增發(fā)功率才可以相應減小待調(diào)節(jié)支路的潮流。

通過對比安全校正方案調(diào)整前后的潮流結果，發(fā)現(xiàn)優(yōu)化后，系統(tǒng)的潮流相比故障后都顯著減少。對比優(yōu)化后該電網(wǎng)潮流和原始潮流，發(fā)現(xiàn)優(yōu)化后的潮流較為貼合原始潮流數(shù)據(jù)，只是略有小幅增加，但也沒有超過潮流極限值，由此驗證了優(yōu)化模型求得的安全校正方案的合理性和正確性。

6 結 語

本文詳細介紹了一種安全校正策略的優(yōu)化模型，該模型將復雜的發(fā)電機和負荷的參數(shù)調(diào)整問題轉(zhuǎn)化為建立、求解非線性的連續(xù)和二元混合規(guī)劃問題，通過限定粒子群算法中粒子的搜索域來實現(xiàn)對優(yōu)化變量的參數(shù)調(diào)節(jié)范圍的約束，并通過求解數(shù)學模型得出控制方案。從根本上避免了越限和一些不符合實際的結果，也去除了大量無效的搜索范圍，使得在更短的時間可以得到合理有效的安全校正方案。