基于電網(wǎng)分區(qū)和類(lèi)擴(kuò)展變量松弛的電力系統(tǒng)潮流優(yōu)化方法
【專(zhuān)利摘要】本發(fā)明公開(kāi)了一種基于電網(wǎng)分區(qū)和類(lèi)擴(kuò)展變量松弛的電力系統(tǒng)潮流優(yōu)化方法,方法包括分解電力系統(tǒng),得到多個(gè)相互獨(dú)立的子系統(tǒng);協(xié)調(diào)各個(gè)子系統(tǒng),得到基于電網(wǎng)分區(qū)的電力系統(tǒng)分解協(xié)調(diào)模型;采用輔助問(wèn)題原理對(duì)電力系統(tǒng)分解協(xié)調(diào)模型進(jìn)行分布式并行最優(yōu)潮流計(jì)算;采用類(lèi)擴(kuò)展變量?jī)?nèi)點(diǎn)法分別求解各子系統(tǒng),得到電力系統(tǒng)的潮流最優(yōu)解。本發(fā)明符合大型互聯(lián)電力系統(tǒng)數(shù)據(jù)分散、結(jié)構(gòu)分布的特點(diǎn);有效地降低了子系統(tǒng)的維數(shù);各個(gè)分區(qū)的優(yōu)化計(jì)算只需要本區(qū)域的數(shù)據(jù)和邊界節(jié)點(diǎn)數(shù)據(jù),不需要區(qū)域外詳細(xì)的數(shù)學(xué)模型,避免了數(shù)據(jù)傳輸?shù)钠款i,加快了優(yōu)化速度;子系統(tǒng)優(yōu)化問(wèn)題采用的類(lèi)擴(kuò)展變量?jī)?nèi)點(diǎn)法擴(kuò)大了優(yōu)化問(wèn)題的解空間,提高了收斂性。
【專(zhuān)利說(shuō)明】
基于電網(wǎng)分區(qū)和類(lèi)擴(kuò)展變量松弛的電力系統(tǒng)潮流優(yōu)化方法
技術(shù)領(lǐng)域
[0001] 本發(fā)明設(shè)及電力系統(tǒng)優(yōu)化技術(shù)領(lǐng)域,尤其是一種基于電網(wǎng)分區(qū)和類(lèi)擴(kuò)展變量松弛 的電力系統(tǒng)潮流優(yōu)化方法。
【背景技術(shù)】
[0002] 基于電網(wǎng)分區(qū)的分解協(xié)調(diào)并行優(yōu)化算法的基本思路是設(shè)法將大系統(tǒng)優(yōu)化問(wèn)題有 效分解為一系列子問(wèn)題,并進(jìn)行分布式求解和協(xié)調(diào)優(yōu)化,最早應(yīng)用于日發(fā)電計(jì)劃優(yōu)化和分 布式最優(yōu)潮流計(jì)算中。求解主要基于巧中考慮:①保留每個(gè)區(qū)域的獨(dú)立性來(lái)進(jìn)行總的優(yōu)化問(wèn) 題的計(jì)算;②W分布方式或并行方式求解最優(yōu)潮流問(wèn)題,提高計(jì)算速度;③使區(qū)域調(diào)度員能 夠協(xié)調(diào)地調(diào)度電網(wǎng),而無(wú)需知道其他區(qū)域的運(yùn)行數(shù)據(jù)。
[0003] 近些年,基于電網(wǎng)分區(qū)的分解協(xié)調(diào)并行優(yōu)化算法來(lái)有了一定的研究,除采用并行 處理方式的智能算法外,其實(shí)現(xiàn)方法目前大致有=類(lèi):第一類(lèi)方法先基于最優(yōu)化問(wèn)題的 Karush K化n Tucker化KT)條件得到全網(wǎng)的系統(tǒng)修正方程,再運(yùn)用對(duì)角加邊模型、近似牛頓 方向等矩陣?yán)碚搶?duì)該系統(tǒng)修正方程進(jìn)行解禪。但該類(lèi)方法計(jì)算效率受制于子網(wǎng)間的禪合 度、區(qū)域分解的數(shù)目,且存在協(xié)調(diào)層負(fù)擔(dān)重等問(wèn)題,更適用于同地分布式并行計(jì)算的場(chǎng)景。 第二類(lèi)方法通過(guò)對(duì)外部網(wǎng)絡(luò)進(jìn)行等值,進(jìn)行各分區(qū)優(yōu)化的獨(dú)立求解,并引入外層協(xié)調(diào)環(huán)節(jié) 實(shí)現(xiàn)全網(wǎng)的優(yōu)化。該類(lèi)方法在進(jìn)行本地計(jì)算時(shí)需要相鄰子網(wǎng)的運(yùn)行參數(shù)W進(jìn)行外網(wǎng)等值, 影響了子網(wǎng)優(yōu)化計(jì)算的獨(dú)立性,且該等值參數(shù)的計(jì)算和更新相對(duì)繁瑣。第=類(lèi)方法基于拉 格朗日松弛理論,將子網(wǎng)間的禪合約束松弛至目標(biāo)函數(shù)中并采用一定的分解算法實(shí)現(xiàn)其可 分性。W基于輔助問(wèn)題原理(Auxi 1 iary Prob 1 em Princ ip 1 e,APP)的分解算法為例,該經(jīng)典 方法通過(guò)構(gòu)造輔助函數(shù),將原鞍點(diǎn)問(wèn)題分化為可獨(dú)立求解的極值問(wèn)題,有效地解決了集中 式并行優(yōu)化算法帶來(lái)的問(wèn)題。
[0004] 非線(xiàn)性?xún)?nèi)點(diǎn)法具有較好的收斂性和較快的計(jì)算速度,已經(jīng)成為求解電力系統(tǒng)優(yōu)化 問(wèn)題的有力工具,但其解空間并不令人滿(mǎn)意。類(lèi)擴(kuò)展變量?jī)?nèi)點(diǎn)法通過(guò)在同種類(lèi)別的不等式 約束方程中增加相同的類(lèi)擴(kuò)展變量,并在目標(biāo)函數(shù)中用罰因子迫使類(lèi)擴(kuò)展變量趨向0實(shí)現(xiàn) 約束條件的弱化,擴(kuò)大了優(yōu)化問(wèn)題的解空間,提高了收斂性。
【發(fā)明內(nèi)容】
[0005] 本發(fā)明針對(duì)現(xiàn)有技術(shù)的不足,提出一種基于電網(wǎng)分區(qū)和類(lèi)擴(kuò)展變量松弛的電力系 統(tǒng)潮流優(yōu)化方法,符合大型互聯(lián)電力系統(tǒng)數(shù)據(jù)分散、結(jié)構(gòu)分布的特點(diǎn);有效地降低了子系統(tǒng) 的維數(shù),使得各子系統(tǒng)的優(yōu)化問(wèn)題規(guī)模較小,提高了全網(wǎng)最優(yōu)潮流的計(jì)算速度;各個(gè)分區(qū)的 優(yōu)化計(jì)算只需要本區(qū)域的數(shù)據(jù)和邊界節(jié)點(diǎn)數(shù)據(jù),不需要區(qū)域外詳細(xì)的數(shù)學(xué)模型,避免了數(shù) 據(jù)傳輸?shù)钠款i,加快了優(yōu)化速度;子系統(tǒng)優(yōu)化問(wèn)題采用的類(lèi)擴(kuò)展變量?jī)?nèi)點(diǎn)法擴(kuò)大了優(yōu)化問(wèn) 題的解空間,提高了收斂性。
[0006] 為了實(shí)現(xiàn)上述發(fā)明目的,本發(fā)明提供W下技術(shù)方案:
[0007] -種基于電網(wǎng)分區(qū)和類(lèi)擴(kuò)展變量松弛的電力系統(tǒng)潮流優(yōu)化方法,所述方法包括如 下步驟:
[0008] 步驟1.分解電力系統(tǒng),得到多個(gè)相互獨(dú)立的子系統(tǒng);
[0009] 步驟2.協(xié)調(diào)各個(gè)所述子系統(tǒng),得到基于電網(wǎng)分區(qū)的電力系統(tǒng)分解協(xié)調(diào)模型;
[0010] 步驟3.采用輔助問(wèn)題原理對(duì)所述電力系統(tǒng)分解協(xié)調(diào)模型進(jìn)行分布式并行最優(yōu)潮 流計(jì)算;
[0011] 步驟4.采用類(lèi)擴(kuò)展變量?jī)?nèi)點(diǎn)法分別求解各所述子系統(tǒng),得到所述電力系統(tǒng)的潮流 最優(yōu)解。
[0012] 進(jìn)一步地,所述步驟1之前,包括:
[0013] 確定所述電力系統(tǒng)的最優(yōu)潮流數(shù)學(xué)模型:
[0014]
(1)
[0015] 式(1)中:f (X)為目標(biāo)函數(shù),X為自變量構(gòu)成的列向量;h(x)和g(x)分別為等式約束 和不等式約束條件;g和容分別為不等式約束條件g(x)的下限值與上限值。
[0016] 進(jìn)一步地,所述步驟1包括:
[0017] 1-1.根據(jù)所述電力系統(tǒng)的地理分布,在聯(lián)絡(luò)線(xiàn)處確定邊界節(jié)點(diǎn);
[0018] 1-2.判斷所述邊界節(jié)點(diǎn)是否為有源節(jié)點(diǎn);
[0019] 若是,則進(jìn)入步驟1-3;
[0020] 若否,則進(jìn)入步驟1-4;
[0021] 1-3.在所述聯(lián)絡(luò)線(xiàn)上增加虛擬的無(wú)源節(jié)點(diǎn),將身上無(wú)源節(jié)點(diǎn)作為新的邊界節(jié)點(diǎn), 進(jìn)入步驟1-4;
[0022] 1-4.將所述電力系統(tǒng)在每個(gè)所述邊界節(jié)點(diǎn)處分解,得到多個(gè)相互獨(dú)立的所述子系 統(tǒng)。
[0023] 進(jìn)一步地,所述步驟2包括:
[0024] 2-1.將所述邊界節(jié)點(diǎn)參數(shù)作為協(xié)調(diào)變量;
[0025] 2-2.各個(gè)所述子系統(tǒng)之間通過(guò)交換所述協(xié)調(diào)變量進(jìn)行協(xié)調(diào),得到基于電網(wǎng)分區(qū)的 電力系統(tǒng)分解協(xié)調(diào)模型。
[00%] 進(jìn)一步地,所述步驟2-2包括:
[0027] a.各個(gè)所述子系統(tǒng)之間通過(guò)交換所述協(xié)調(diào)變量進(jìn)行協(xié)調(diào),得到基于電網(wǎng)分區(qū)的所 述
[0028] 電力系統(tǒng)分解協(xié)調(diào)模型的目標(biāo)函數(shù)minf (X):
[0029] minf (X) =fi(xi)+f2(X2)... .+fn(Xn) (2)
[0030] 式(2)中,Xl、X2...Xn分別為所述各個(gè)子系統(tǒng);fl(Xl)、f2(X2)...fn(Xn)分別為各所 述子系統(tǒng)的目標(biāo)函數(shù);
[0031] b.確定所述電力系統(tǒng)分解協(xié)調(diào)模型的邊界約束。
[0032] 進(jìn)一步地,所述步驟3包括:
[0033] 3-1.根據(jù)所述輔助問(wèn)題原理中的核函數(shù),構(gòu)造所述電力系統(tǒng)分解協(xié)調(diào)模型的輔助 函數(shù);
[0034] 3-2.將所述電力系統(tǒng)分解協(xié)調(diào)模型求解的主問(wèn)題轉(zhuǎn)化為求解所述輔助問(wèn)題的鞍 點(diǎn)問(wèn)題;
[0035] 3-3.根據(jù)設(shè)置的闊值對(duì)所述輔助問(wèn)題進(jìn)行迭代,直到循環(huán)結(jié)束,得到所述電力系 統(tǒng)分解協(xié)調(diào)模型的分布式并行最優(yōu)潮流計(jì)算結(jié)果。
[0036] 進(jìn)一步地,所述步驟4包括:
[0037] 4-1.在各所述子系統(tǒng)中的同類(lèi)別的不等式約束方程中增加相同的類(lèi)擴(kuò)展變量,且 所述類(lèi)擴(kuò)展變量按上限和下限分別設(shè)置;
[0038] 4-2.在各所述子系統(tǒng)的目標(biāo)函數(shù)中增加帶罰因子的擴(kuò)展變量二次項(xiàng),擴(kuò)大不等式 約束的解空間;
[0039] 4-3.將對(duì)數(shù)障礙函數(shù)化不等式約束作為等式約束,用拉格朗日乘子法處理等式約 束,得到引入類(lèi)擴(kuò)展變量后的拉格朗日函數(shù);
[0040] 4-4.推導(dǎo)KKT條件,得到所述引入類(lèi)擴(kuò)展變量后的拉格朗日函數(shù)的修正方程組;
[0041] 4-5.求解所述修正方程,得到所述電力系統(tǒng)的潮流最優(yōu)解。
[0042] 進(jìn)一步地,所述步驟4-5包括:
[0043] (1)初始化所述方程組:置迭代次數(shù)為0,給定自變量、松弛因子、拉格朗日乘子、類(lèi) 擴(kuò)展變量的初值;
[0044] (2)計(jì)算目標(biāo)函數(shù)和約束的雅可比矩陣及海森矩陣,形成線(xiàn)性系統(tǒng);
[0045] (3)求解所述線(xiàn)性系統(tǒng),得到每次迭代的修正量;
[0046] (4)計(jì)算原變量仿射步長(zhǎng)和對(duì)偶變量仿射步長(zhǎng);
[0047] (5)更新各變量及類(lèi)擴(kuò)展變量;
[0048] (6)根據(jù)設(shè)定的收斂標(biāo)準(zhǔn),判斷最優(yōu)潮流是否收斂;
[0049] 若是,則停止計(jì)算,得到所述電力系統(tǒng)的潮流最優(yōu)解;
[0050] 若否,則返回步驟(2)。
[0051] 本發(fā)明提出一種將輔助問(wèn)題原理(APP)和類(lèi)擴(kuò)展變量?jī)?nèi)點(diǎn)法相結(jié)合的電力系統(tǒng)并 行最優(yōu)潮流算法,利用APP建立起基于電網(wǎng)分區(qū)的分解協(xié)調(diào)模型,利用類(lèi)擴(kuò)展變量?jī)?nèi)點(diǎn)法來(lái) 處理各子系統(tǒng)的優(yōu)化問(wèn)題,從而在降低計(jì)算維數(shù)、縮短計(jì)算時(shí)間的同時(shí)提高算法收斂性。 IEEE 118節(jié)點(diǎn)網(wǎng)絡(luò)的仿真分析驗(yàn)證了該算法的合理性和有效性。其有益效果為:
[0052] (1)本發(fā)明提出的方法符合大型互聯(lián)電力系統(tǒng)數(shù)據(jù)分散、結(jié)構(gòu)分布的特點(diǎn);
[0053] (2)該方法有效地降低了子系統(tǒng)的維數(shù),使得各子系統(tǒng)的優(yōu)化問(wèn)題規(guī)模較小,提高 了全網(wǎng)最優(yōu)潮流的計(jì)算速度;
[0054] (3)各個(gè)分區(qū)的優(yōu)化計(jì)算只需要本區(qū)域的數(shù)據(jù)和邊界節(jié)點(diǎn)數(shù)據(jù),不需要區(qū)域外詳 細(xì)的數(shù)學(xué)模型,避免了數(shù)據(jù)傳輸?shù)钠款i,加快了優(yōu)化速度;
[0055] (4)子系統(tǒng)優(yōu)化問(wèn)題采用的類(lèi)擴(kuò)展變量?jī)?nèi)點(diǎn)法擴(kuò)大了優(yōu)化問(wèn)題的解空間,提高了 收斂性。
【附圖說(shuō)明】
[0056] 圖1為本發(fā)明的基于電網(wǎng)分區(qū)的電力系統(tǒng)的潮流并行優(yōu)化方法的流程示意圖;
[0057] 圖2為本發(fā)明的W二分區(qū)為例中的電力系統(tǒng)分解前的結(jié)構(gòu)示意圖;
[0058] 圖3為本發(fā)明的W二分區(qū)為例中的電力系統(tǒng)分解后的結(jié)構(gòu)示意圖;
[0059] 圖4為本發(fā)明的具體應(yīng)用例的計(jì)算流程圖。
【具體實(shí)施方式】
[0060] 下面結(jié)合附圖對(duì)本發(fā)明進(jìn)行詳細(xì)描述,本部分的描述僅是示范性和解釋性,不應(yīng) 對(duì)本發(fā)明的保護(hù)范圍有任何的限制作用。
[0061] 如圖1所示,本發(fā)明提供一種基于電網(wǎng)分區(qū)和類(lèi)擴(kuò)展變量松弛的電力系統(tǒng)潮流優(yōu) 化方法,包括如下步驟:
[0062] 步驟1.分解電力系統(tǒng),得到多個(gè)相互獨(dú)立的子系統(tǒng);
[0063] 步驟2.協(xié)調(diào)各個(gè)子系統(tǒng),得到基于電網(wǎng)分區(qū)的電力系統(tǒng)分解協(xié)調(diào)模型;
[0064] 步驟3.采用輔助問(wèn)題原理對(duì)電力系統(tǒng)分解協(xié)調(diào)模型進(jìn)行分布式并行最優(yōu)潮流計(jì) 算;
[0065] 步驟4.采用類(lèi)擴(kuò)展變量?jī)?nèi)點(diǎn)法分別求解各子系統(tǒng),得到電力系統(tǒng)的潮流最優(yōu)解。
[0066] 其中,步驟1之前,包括:
[0067] 確定電力系統(tǒng)的最優(yōu)潮流數(shù)學(xué)模型:
[0068]
(1)
[0069] 式(1)中:f (X)為目標(biāo)函數(shù),X為自變量構(gòu)成的列向量;h(x)和g(x)分別為等式約束 和不等式約束條件;g和蓉分別為不等式約束條件g(x)的下限值與上限值。
[0070] 其中,步驟1包括:
[0071 ] 1-1.根據(jù)電力系統(tǒng)的地理分布,在聯(lián)絡(luò)線(xiàn)處確定邊界節(jié)點(diǎn);
[0072] 1-2.判斷邊界節(jié)點(diǎn)是否為有源節(jié)點(diǎn);
[0073] 若是,則進(jìn)入步驟1-3;
[0074] 若否,則進(jìn)入步驟1-4;
[0075] 1-3.在聯(lián)絡(luò)線(xiàn)上增加虛擬的無(wú)源節(jié)點(diǎn),將身上無(wú)源節(jié)點(diǎn)作為新的邊界節(jié)點(diǎn),進(jìn)入 步驟1-4;
[0076] 1-4.將電力系統(tǒng)在每個(gè)邊界節(jié)點(diǎn)處分解,得到多個(gè)相互獨(dú)立的子系統(tǒng)。
[0077] 其中,步驟2包括:
[0078] 2-1.將邊界節(jié)點(diǎn)參數(shù)作為協(xié)調(diào)變量;
[0079] 2-2.各個(gè)子系統(tǒng)之間通過(guò)交換協(xié)調(diào)變量進(jìn)行協(xié)調(diào),得到基于電網(wǎng)分區(qū)的電力系統(tǒng) 分解協(xié)調(diào)模型。
[0080] 其中,步驟2-2包括:
[0081] a.各個(gè)子系統(tǒng)之間通過(guò)交換協(xié)調(diào)變量進(jìn)行協(xié)調(diào),得到基于電網(wǎng)分區(qū)的電力系統(tǒng)分 解協(xié)調(diào)模型的目標(biāo)函數(shù)minf (X):
[0082] minf (X) =fi(xi)+f2(X2)... .+fn(Xn) (2)
[0083] 式(2)中,xl、x2...xn分別為各個(gè)子系統(tǒng);fl(xl)、f 2(x2)...f n(xn)分別為各子 系統(tǒng)的目標(biāo)函數(shù);
[0084] b.確定電力系統(tǒng)分解協(xié)調(diào)模型的邊界約束。
[00化]其中,步驟3包括:
[0086] 3-1.根據(jù)輔助問(wèn)題原理中的核函數(shù),構(gòu)造電力系統(tǒng)分解協(xié)調(diào)模型的輔助函數(shù);
[0087] 3-2.將電力系統(tǒng)分解協(xié)調(diào)模型求解的主問(wèn)題轉(zhuǎn)化為求解輔助問(wèn)題的鞍點(diǎn)問(wèn)題;
[0088] 3-3.根據(jù)設(shè)置的闊值對(duì)輔助問(wèn)題進(jìn)行迭代,直到循環(huán)結(jié)束,得到電力系統(tǒng)分解協(xié) 調(diào)模型的分布式并行最優(yōu)潮流計(jì)算結(jié)果。
[0089] 其中,步驟4包括:
[0090] 4-1.在各子系統(tǒng)中的同類(lèi)別的不等式約束方程中增加相同的類(lèi)擴(kuò)展變量,且類(lèi)擴(kuò) 展變量按上限和下限分別設(shè)置;
[0091] 4-2.在各子系統(tǒng)的目標(biāo)函數(shù)中增加帶罰因子的擴(kuò)展變量二次項(xiàng),擴(kuò)大不等式約束 的解空間;
[0092] 4-3.將對(duì)數(shù)障礙函數(shù)化不等式約束作為等式約束,用拉格朗日乘子法處理等式約 束,得到引入類(lèi)擴(kuò)展變量后的拉格朗日函數(shù);
[009引4-4.推導(dǎo)KKT條件,得到引入類(lèi)擴(kuò)展變量后的拉格朗日函數(shù)的修正方程組;
[0094] 4-5.求解修正方程,得到電力系統(tǒng)的潮流最優(yōu)解。
[00巧]其中,步驟4-5包括:
[0096] (1)初始化方程組:置迭代次數(shù)為0,給定自變量、松弛因子、拉格朗日乘子、類(lèi)擴(kuò)展 變量的初值;
[0097] (2)計(jì)算目標(biāo)函數(shù)和約束的雅可比矩陣及海森矩陣,形成線(xiàn)性系統(tǒng);
[0098] (3)求解線(xiàn)性系統(tǒng),得到每次迭代的修正量;
[0099] (4)計(jì)算原變量仿射步長(zhǎng)和對(duì)偶變量仿射步長(zhǎng);
[0100] (5)更新各變量及類(lèi)擴(kuò)展變量;
[0101] (6)根據(jù)設(shè)定的收斂標(biāo)準(zhǔn),判斷最優(yōu)潮流是否收斂;
[0102] 若是,則停止計(jì)算,得到電力系統(tǒng)的潮流最優(yōu)解;
[0103] 若否,則返回步驟(2)。
[0104] 本發(fā)明W二分區(qū)為例介紹大電網(wǎng)的分區(qū)及其并行優(yōu)化模型,如下:
[0105] 1、最優(yōu)潮流數(shù)學(xué)模型
[0106] 最優(yōu)潮流計(jì)算在數(shù)學(xué)上可歸結(jié)為有約束的非線(xiàn)性規(guī)劃問(wèn)題,其數(shù)學(xué)模型一般為:
[0107]
(1)
[0108] 式中:f (X)為目標(biāo)函數(shù),最優(yōu)潮流的目標(biāo)函數(shù)通常為系統(tǒng)總的發(fā)電費(fèi)用(發(fā)出有功 無(wú)功的成本)最小,或系統(tǒng)總的線(xiàn)路有功功率損耗最小;X為自變量構(gòu)成的列向量,表示發(fā)電 機(jī)有功、無(wú)功出力W及節(jié)點(diǎn)電壓幅值角度等控制變量或狀態(tài)變量;h(x)和g(x)分別表示為 等式約束和不等式約束條件。
[0109] 2、電網(wǎng)分區(qū)
[0110] 最優(yōu)潮流計(jì)算的復(fù)雜程度與電網(wǎng)規(guī)模密切相關(guān),為降低問(wèn)題求解的難度,可W考 慮根據(jù)電網(wǎng)的實(shí)際地理分布,在某些聯(lián)絡(luò)線(xiàn)處將整個(gè)電網(wǎng)分解為多個(gè)相對(duì)獨(dú)立的區(qū)域,分 別進(jìn)行計(jì)算。本發(fā)明采用分解協(xié)調(diào)方法,即通過(guò)"復(fù)制"邊界節(jié)點(diǎn)將一個(gè)大系統(tǒng)直接分解成 多個(gè)子系統(tǒng),各個(gè)子系統(tǒng)相對(duì)獨(dú)立,W邊界節(jié)點(diǎn)的參數(shù)為協(xié)調(diào)變量,每個(gè)分區(qū)的解只和自己 的內(nèi)部變量和邊界變量有關(guān),分區(qū)之間僅依靠交換邊界節(jié)點(diǎn)數(shù)據(jù)進(jìn)行協(xié)調(diào)。下面W二分區(qū) 為例介紹大電網(wǎng)的分區(qū)及其并行優(yōu)化過(guò)程:
[0111] 如圖2所示,電網(wǎng)S的兩個(gè)子系統(tǒng)乂1,)(2通過(guò)邊界Xb相連,相連的支路稱(chēng)為聯(lián)絡(luò)線(xiàn)支 路。分解后如圖3所示,"復(fù)制"邊界節(jié)點(diǎn)Xb,進(jìn)行網(wǎng)絡(luò)分離。用XbiQ = I,2)表示各子系統(tǒng)的 "復(fù)制"邊界節(jié)點(diǎn),C康示分離后的每一個(gè)分區(qū),用E表示所有分區(qū)的集合,即E=促} Q = 1, 2)。
[0112] 可W把每一個(gè)邊界節(jié)點(diǎn)看作一個(gè)虛擬的發(fā)電機(jī)節(jié)點(diǎn),Xbi為第i分區(qū)虛擬發(fā)電機(jī)節(jié) 點(diǎn)的參數(shù)變量。虛擬發(fā)電機(jī)彌補(bǔ)了因?yàn)閷?duì)電網(wǎng)的直接分解而造成的潮流不平衡,同時(shí)也使 得Cl形成了一個(gè)相對(duì)獨(dú)立的子系統(tǒng)。
[0113] 要使分解后的系統(tǒng)同原系統(tǒng)等值,應(yīng)規(guī)定Xbl和Xb2在電氣上為同一點(diǎn),具有相同的 電氣量,即
舊等,即Xbl = Xb2。其中,Ch和 氣C =1,2 )表示邊界節(jié)點(diǎn)的虛擬發(fā)電機(jī)注入功率;F%,:郝0'h (f二I,2)表示電壓幅值和相 角在進(jìn)行分區(qū)時(shí),若實(shí)際電網(wǎng)的邊界節(jié)點(diǎn)為有源節(jié)點(diǎn),應(yīng)先在聯(lián)絡(luò)線(xiàn)支路上增加虛擬的無(wú) 源節(jié)點(diǎn),取代原來(lái)的邊界節(jié)點(diǎn),然后再進(jìn)行分區(qū)。同時(shí),當(dāng)分區(qū)內(nèi)的約束與區(qū)外變量有關(guān)時(shí), 可W將運(yùn)個(gè)變量所在的節(jié)點(diǎn)看作邊界節(jié)點(diǎn),使分區(qū)約束只與本分區(qū)內(nèi)的W及邊界節(jié)點(diǎn)的變 量有關(guān),與區(qū)外的變量無(wú)關(guān)。
[0114] 運(yùn)樣,原目標(biāo)函數(shù)f(x)可W寫(xiě)成各分區(qū)的目標(biāo)函數(shù)之和,記為fl(Xl)+f2(X2),Xl和 X2分別表示2個(gè)區(qū)域的內(nèi)部參數(shù)變量。全網(wǎng)的等式和不等式約束分解為各區(qū)域的約束。對(duì)于 邊界節(jié)點(diǎn),由于子系統(tǒng)間在邊界節(jié)點(diǎn)上應(yīng)具有相同的電氣量,應(yīng)該增加一個(gè)等式約束條件。
[0115] 采用上述分解協(xié)調(diào)法后,得到2分區(qū)等值最優(yōu)潮流分解協(xié)調(diào)模型如下:
[0116] 目標(biāo)函數(shù):
[0117] minf (X) =fi(xi)+f2(X2) (2)
[011引邊界約束:
[0119] hi(xi) =0,h2(X2) =0
[0120]
(3)
[0121] 目(X) =Xb 廣 Xb2 = 0
[01。]式中,義=[財(cái),別1,町別2^,。佔(zhàn))為分區(qū)。的目標(biāo)函數(shù),目^)為全局約束。
[0123] 對(duì)于多分區(qū)的電力系統(tǒng),其分解協(xié)調(diào)模型建立的主要思想可用W上2分區(qū)系統(tǒng)來(lái) 說(shuō)明。
[0124] 3、基于輔助問(wèn)題原理的并行最優(yōu)潮流計(jì)算方法
[0125] 不考慮子系統(tǒng)內(nèi)部的約束條件,僅考慮由于分區(qū)產(chǎn)生的邊界點(diǎn)約束,此時(shí)根據(jù)式 (2)、式(3),優(yōu)化問(wèn)題的拉格朗日函數(shù)可W寫(xiě)為:
[0126] L(x,人)=fi(xi)+f2(X2)+<^,目(X)〉 (4)
[0127] 為了提高算法的收斂性,采用增廣拉格朗日方法,即在原拉格朗日函數(shù)上增加一 個(gè)二次項(xiàng):
[012引
巧)
[0129] 式中,A為拉格朗日乘子向量,C為懲罰因子,是一個(gè)常數(shù)。
[0130] 可W看出,引入二次項(xiàng)并不會(huì)影響最終的計(jì)算結(jié)果,當(dāng)最后迭代收斂時(shí),子系統(tǒng)間 邊界節(jié)點(diǎn)電氣量趨向一致,此時(shí)二次項(xiàng)的值趨向0。由于0(x)是和2個(gè)區(qū)域變量都相關(guān)的函 數(shù),且加入的二次項(xiàng)本身是不可分解的,從而破壞了兩個(gè)子系統(tǒng)之間的可分性,使得優(yōu)化計(jì) 算無(wú)法直接在2個(gè)分區(qū)內(nèi)各自獨(dú)立地完成。在運(yùn)種情況下,引入輔助問(wèn)題原理(APP)來(lái)解決 運(yùn)個(gè)問(wèn)題。
[0131] 3.1輔助問(wèn)題原理(APP)
[0132] APP的思想是,通過(guò)構(gòu)造核函數(shù),將求解LU, A)的鞍點(diǎn)問(wèn)題轉(zhuǎn)化成求解一系列迭代 形式的輔助鞍點(diǎn)問(wèn)題,而每一次迭代的輔助問(wèn)題均為可分解的問(wèn)題。
[0133] W求極小值為例,說(shuō)明輔助算法原理。假設(shè)Ei(X)可微,且LipscMtz連續(xù),函數(shù)E2 (X)不一定可微,對(duì)于主問(wèn)題mi址i(x)+E2(x)構(gòu)造一輔助問(wèn)題:minG(x)+eE2(x),可W證明若 存在X氣吏得
[0134] G'(X*)=巧'iU*) (6)
[0135] 成立,則X*也成為主問(wèn)題的解,因此,構(gòu)造輔助函數(shù)G(x,v),其中,W(X)為核函數(shù):
[0136] G(x,v) = W(x)+<eE'i(v)-W'(v),x〉 (7)
[0137] 則主問(wèn)題可W轉(zhuǎn)化為求解W下輔助問(wèn)題的鞍點(diǎn)問(wèn)題,即:
[013引 xk" = arg minW (x)+<eE'I(^k)-W'(yk) ,x〉+eE2(x) (8)
[0139] 式中,argmin表示使輔助問(wèn)題W (x)+<eE' i(xk)-W' (yk) ,x〉+巧2(x)取得最小值時(shí) X的值。當(dāng)||xk-xk"| I小于指定闊值時(shí)停止,否則繼續(xù)迭代。
[0140] 3.2基于A(yíng)PP的并行最優(yōu)潮流計(jì)算方法
[0141] 對(duì)于主問(wèn)題LU,人),令
[0142] L(x,A)=f(x)+J(x,A) (9)
[0143] f (x) =fi(xi)+f2(x2) (10)
[0144]
(11)
[0145] 其中,J(x,A)可微,f (X)不一定可微。構(gòu)造輔助函數(shù)為:
[0146]
(12)
[0147] 井由U責(zé)雲(yún)取搶!造的核函數(shù),構(gòu)造為;
[0148] (13)
[0149] e在求極小值和極大值問(wèn)題過(guò)程中分別取E = Ei = I和e = E2;b為正數(shù)。
[0150] 根據(jù)輔助問(wèn)題原理的兩層算法模型,在第k次迭代中,求解LU,A)的鞍點(diǎn)問(wèn)題相當(dāng) 于交替求解W下極小值和極大值問(wèn)題:
[0151] 解極小值問(wèn)題:
[0152]
(14)
[0153] 解極大值問(wèn)題:
[0154] 人 k+i 二人 k+E2 目(xk+i) (15)
[0156]
[0155] 很明顯,式(14)中不會(huì)再出現(xiàn)不同變量間的交叉乘積,可被進(jìn)一步分解。對(duì)于2分 區(qū)系統(tǒng),將f(x),0(x)和X代入式(14)、式(15),令e=l+b,整理得:
[0157]
[0158] 不難看出,式(16)可W由各分區(qū)并行求解式(18)夾連現(xiàn):
[0159]
(1汾
[0160] 式中,i = i ,2; i = l時(shí),q = l; i = 2時(shí),q = -l。
[0161] 4、子系統(tǒng)優(yōu)化問(wèn)題的類(lèi)擴(kuò)展變量?jī)?nèi)點(diǎn)法求解
[0162] 前面將一個(gè)大網(wǎng)絡(luò)的優(yōu)化問(wèn)題分解為2個(gè)子系統(tǒng)的優(yōu)化問(wèn)題,并針對(duì)各個(gè)子系統(tǒng) 形成了各自的目標(biāo)函數(shù)及約束條件,W子系統(tǒng)Cl為例,其子優(yōu)化問(wèn)題及其約束為:
[0163]
[0164] 為了提高算法的收斂性,運(yùn)里采用類(lèi)擴(kuò)展變量?jī)?nèi)點(diǎn)法求解子系統(tǒng)的優(yōu)化問(wèn)題。類(lèi) 擴(kuò)展變量?jī)?nèi)點(diǎn)法通過(guò)在同種類(lèi)別的不等式約束方程中增加相同的類(lèi)擴(kuò)展變量,并在目標(biāo)函 數(shù)中用罰因子迫使類(lèi)擴(kuò)展變量的平方和趨向0,從而擴(kuò)大了不等式約束的解空間。當(dāng)原問(wèn)題 可行時(shí),該算法可W收斂到原問(wèn)題的最優(yōu)解;當(dāng)原問(wèn)題無(wú)解時(shí),可W自動(dòng)到更大的可行域內(nèi) 尋優(yōu),快速得到近似解,從計(jì)算結(jié)果中可W方便地得到調(diào)整的措施。
[01化]4.1類(lèi)擴(kuò)展變量的引入
[0166] 將式(19)中不等式約束按類(lèi)別增加擴(kuò)展變量,考慮到同一類(lèi)約束從上限和下限松 弛的需求通常相互獨(dú)立,類(lèi)擴(kuò)展變量需按上限和下限分別設(shè)置。為使優(yōu)化問(wèn)題的解盡量向 原設(shè)定的陥估吿:嫌-#曰拉尿,掀由墻巾n惡巧巧子的擴(kuò)展變量二次項(xiàng)。則式(19)變?yōu)椋?br>[0167] (20)
[0168] 式(20)中,n為設(shè)置擴(kuò)展變量的不等式約束的種類(lèi),a、e為類(lèi)擴(kuò)展變量列向量,分別 代表下擴(kuò)展變量和上擴(kuò)展變量,維度均為n; ak、&分別表示Cl、e的第k個(gè)元素;Ck為第k個(gè)類(lèi)擴(kuò) 展變量二次項(xiàng)的罰因子,T為不等式約束和類(lèi)擴(kuò)展變量的聯(lián)系矩陣,為r Xn維矩陣,形式如 下所示:
[0169]
[0170] 若式(20)中第i個(gè)不等式約束中增加了第j個(gè)類(lèi)擴(kuò)展變量,則Tl,J為1,T第i行的其 它元素均為0。
[0171] 4.2子系統(tǒng)優(yōu)化問(wèn)題的求解方案
[0172] 根據(jù)內(nèi)點(diǎn)法的核屯、思想,用對(duì)數(shù)障礙函數(shù)化不等式約束為等式約束,用拉格朗日 乘子法處理等式約束,與式(20)對(duì)應(yīng)引入類(lèi)擴(kuò)展變量后的拉格朗日函數(shù)為:
[0180] Fl=-▽巧 i(xi,孔 i)ZL_1T (26)[0181] F2 = -▽巧 i(xi,孔 i)WU_1t (27)[0182] F3 = 2C 巧 Tzl-It (28)[0183] F4 = 2C-TW^T (29)[0184] Lx,=Lx+Vjcgi(xi,孔i)(Lz,-Lw,) (30)[0185] LZ=La^TLz' (31)[0186] Le ' =L葉 Tlw ' (32)[0187] Lz'=L_i(ZLz+Li) (33)[018 引 Lw ' =U-I(WLw-Lu) (34)[0189] 方程組求解步驟如下:[0190] (1)初始化:置迭代次數(shù)為0,給定自變量、松弛因子、拉格朗日乘子、類(lèi)擴(kuò)展變量的
[0173]
[0174]
[0175]
[0176]
[0177]
[017 引
[0179] 初值;
[0191] (2)計(jì)算目標(biāo)函數(shù)和約束的雅可比矩陣及海森矩陣,形成式(22)~(24)所示的線(xiàn) 性系統(tǒng);
[0192] (3)求解式(22)~(24),得到每次迭代的修正量;
[01 93] (4 H+當(dāng)原巧音化射擊長(zhǎng)Agfn巧對(duì)儲(chǔ)巧章仿射步長(zhǎng)Aafd :
[0194] (游)
[0195] :
[0196] 巧技)
[0197] 式(36)中,對(duì)偶間隙P和仿射對(duì)偶間隙Paf為:
[019 引 P = ITz-UV (37)
[0199] Paf= (1+入a印 A laf)T(Z+入afd A Zaf)-(U+入a印 A Uaf)T(w+入afd A Waf) (38)
[0200] (5)按下式更新各變量及類(lèi)擴(kuò)展變量:
[0201]
(39)
[0202] (6)根據(jù)設(shè)定的收斂標(biāo)準(zhǔn),判斷最優(yōu)潮流是否收斂。若收斂則停止計(jì)算,否則轉(zhuǎn)步 驟2。
[0203] 如圖4所示,具體應(yīng)用例5WIEEE 118節(jié)點(diǎn)測(cè)試系統(tǒng)為例,對(duì)本發(fā)明討論的計(jì)算方 法進(jìn)行驗(yàn)證和分析。其中,對(duì)系統(tǒng)進(jìn)行二分區(qū),節(jié)點(diǎn)23,68,69作為分界點(diǎn)。測(cè)試系統(tǒng)的基本 數(shù)據(jù)信息如表1所示。
[0204] 夫IT腳,R 11 8韋T占娜陸系統(tǒng)其太掛據(jù)
[0205]
[0206]
[0207] W系統(tǒng)網(wǎng)損最小為目標(biāo)函數(shù),為了比較算法的計(jì)算效果,采用W下=種方式對(duì)測(cè) 試系統(tǒng)進(jìn)行最優(yōu)潮流計(jì)算。
[020引(1)電網(wǎng)未分區(qū),采用非線(xiàn)性原-對(duì)偶內(nèi)點(diǎn)法計(jì)算。
[0209] (2)電網(wǎng)二分區(qū),采用非線(xiàn)性原-對(duì)偶內(nèi)點(diǎn)法計(jì)算。
[0210] (3)電網(wǎng)二分區(qū),采用類(lèi)擴(kuò)展變量?jī)?nèi)點(diǎn)法計(jì)算。
[OW]其中,從邊界差值0(x)<O.O3和迭代次數(shù)>50作為全系統(tǒng)優(yōu)化收斂的判據(jù),Wp< 1D-5和M h(x) M w<2D-4作為子系統(tǒng)優(yōu)化的收斂判據(jù)。優(yōu)化結(jié)果如表2所示。
[0212] 表2優(yōu)化結(jié)果對(duì)比
[0213]
[0214] 由表2可W看出,在計(jì)算時(shí)間上,采用了分區(qū)算法后有明顯下降。在方法1中,電網(wǎng) 未分區(qū)時(shí),盡管迭代次數(shù)較少,但每次迭代耗時(shí)過(guò)多,故總時(shí)間仍遠(yuǎn)高于分區(qū)算法。在方法3 中,不僅采用分區(qū)算法節(jié)省了總的計(jì)算時(shí)間,而且由于類(lèi)擴(kuò)展變量的引入,擴(kuò)大了優(yōu)化問(wèn)題 的解空間,使得優(yōu)化后的發(fā)電費(fèi)用相比方法2有所下降,而總迭代次數(shù)卻無(wú)顯著變化。
[0215] 本發(fā)明一種基于電網(wǎng)分區(qū)和類(lèi)擴(kuò)展變量松弛的電力系統(tǒng)潮流優(yōu)化方法,其有益效 果為:
[0216] (1)本發(fā)明提出的方法符合大型互聯(lián)電力系統(tǒng)數(shù)據(jù)分散、結(jié)構(gòu)分布的特點(diǎn);
[0217] (2)該方法有效地降低了子系統(tǒng)的維數(shù),使得各子系統(tǒng)的優(yōu)化問(wèn)題規(guī)模較小,提高 了全網(wǎng)最優(yōu)潮流的計(jì)算速度;
[0218] (3)各個(gè)分區(qū)的優(yōu)化計(jì)算只需要本區(qū)域的數(shù)據(jù)和邊界節(jié)點(diǎn)數(shù)據(jù),不需要區(qū)域外詳 細(xì)的數(shù)學(xué)模型,避免了數(shù)據(jù)傳輸?shù)钠款i,加快了優(yōu)化速度;
[0219] (4)子系統(tǒng)優(yōu)化問(wèn)題采用的類(lèi)擴(kuò)展變量?jī)?nèi)點(diǎn)法擴(kuò)大了優(yōu)化問(wèn)題的解空間,提高了 收斂性。
[0220] 最后應(yīng)當(dāng)說(shuō)明的是:W上實(shí)施例僅用W說(shuō)明本發(fā)明的技術(shù)方案而非對(duì)其限制,盡 管參照上述實(shí)施例對(duì)本發(fā)明進(jìn)行了詳細(xì)的說(shuō)明,所屬領(lǐng)域的普通技術(shù)人員應(yīng)當(dāng)理解:依然 可W對(duì)本發(fā)明的【具體實(shí)施方式】進(jìn)行修改或者等同替換,而未脫離本發(fā)明精神和范圍的任何 修改或者等同替換,其均應(yīng)涵蓋在本發(fā)明的權(quán)利要求范圍當(dāng)中。
【主權(quán)項(xiàng)】
1. 一種基于電網(wǎng)分區(qū)和類(lèi)擴(kuò)展變量松弛的電力系統(tǒng)潮流優(yōu)化方法,其特征在于,所述 方法包括如下步驟: 步驟1.分解電力系統(tǒng),得到多個(gè)相互獨(dú)立的子系統(tǒng); 步驟2.協(xié)調(diào)各個(gè)所述子系統(tǒng),得到基于電網(wǎng)分區(qū)的電力系統(tǒng)分解協(xié)調(diào)模型; 步驟3.采用輔助問(wèn)題原理對(duì)所述電力系統(tǒng)分解協(xié)調(diào)模型進(jìn)行分布式并行最優(yōu)潮流計(jì) 算; 步驟4.采用類(lèi)擴(kuò)展變量?jī)?nèi)點(diǎn)法分別求解各所述子系統(tǒng),得到所述電力系統(tǒng)的潮流最優(yōu) 解。2. 如權(quán)利要求1所述的方法,其特征在于,所述步驟1之前,確定所述電力系統(tǒng)的最優(yōu)潮 流數(shù)學(xué)模型:(1) 式(1)中:f(x)為目標(biāo)函數(shù),X為自變量構(gòu)成的列向量;h(x)和g(x)分別為等式約束和不 等式約束條件;S和穿分別為不等式約束條件g(x)的下限值與上限值。3. 如權(quán)利要求2所述的方法,其特征在于,所述步驟1包括: 1-1.根據(jù)所述電力系統(tǒng)的地理分布,在聯(lián)絡(luò)線(xiàn)處確定邊界節(jié)點(diǎn); 1-2.判斷所述邊界節(jié)點(diǎn)是否為有源節(jié)點(diǎn); 若是,則進(jìn)入步驟1-3; 若否,則進(jìn)入步驟1-4; 1-3.在所述聯(lián)絡(luò)線(xiàn)上增加虛擬的無(wú)源節(jié)點(diǎn),將身上無(wú)源節(jié)點(diǎn)作為新的邊界節(jié)點(diǎn),進(jìn)入 步驟1-4; 1- 4.將所述電力系統(tǒng)在每個(gè)所述邊界節(jié)點(diǎn)處分解,得到多個(gè)相互獨(dú)立的所述子系統(tǒng)。4. 如權(quán)利要求3所述的方法,其特征在于,所述步驟2包括: 2- 1.將所述邊界節(jié)點(diǎn)參數(shù)作為協(xié)調(diào)變量; 2- 2.各個(gè)所述子系統(tǒng)之間通過(guò)交換所述協(xié)調(diào)變量進(jìn)行協(xié)調(diào),得到基于電網(wǎng)分區(qū)的電力 系統(tǒng)分解協(xié)調(diào)模型。5. 如權(quán)利要求4所述的方法,其特征在于,所述步驟2-2包括: a. 各個(gè)所述子系統(tǒng)之間通過(guò)交換所述協(xié)調(diào)變量進(jìn)行協(xié)調(diào),得到基于電網(wǎng)分區(qū)的所述電 力系統(tǒng)分解協(xié)調(diào)模型的目標(biāo)函數(shù)minf (X): minf (x)=fl(xi)+f2(x2)··· .+fn(xn) (2) 式(2)中,xl、x2...xn分別為所述各個(gè)子系統(tǒng);fl(Xl)、f2(X2)...fn(Xn)分別為各所述子 系統(tǒng)的目標(biāo)函數(shù); b. 確定所述電力系統(tǒng)分解協(xié)調(diào)模型的邊界約束。6. 如權(quán)利要求5所述的方法,其特征在于,所述步驟3包括: 3- 1.根據(jù)所述輔助問(wèn)題原理中的核函數(shù),構(gòu)造所述電力系統(tǒng)分解協(xié)調(diào)模型的輔助函 數(shù); 3-2.將所述電力系統(tǒng)分解協(xié)調(diào)模型求解的主問(wèn)題轉(zhuǎn)化為求解所述輔助問(wèn)題的鞍點(diǎn)問(wèn) 題; 3- 3.根據(jù)設(shè)置的闊值對(duì)所述輔助問(wèn)題進(jìn)行迭代,直到循環(huán)結(jié)束,得到所述電力系統(tǒng)分 解協(xié)調(diào)模型的分布式并行最優(yōu)潮流計(jì)算結(jié)果。7. 如權(quán)利要求6所述的方法,其特征在于,所述步驟4包括: 4- 1.在各所述子系統(tǒng)中的同類(lèi)別的不等式約束方程中增加相同的類(lèi)擴(kuò)展變量,且所述 類(lèi)擴(kuò)展變量按上限和下限分別設(shè)置; 4-2.在各所述子系統(tǒng)的目標(biāo)函數(shù)中增加帶罰因子的擴(kuò)展變量二次項(xiàng),擴(kuò)大不等式約束 的解空間; 4-3 .將對(duì)數(shù)障礙函數(shù)化不等式約束作為等式約束,用拉格朗日乘子法處理等式約束, 得到引入類(lèi)擴(kuò)展變量后的拉格朗日函數(shù); 4-4.推導(dǎo)KKT條件,得到所述引入類(lèi)擴(kuò)展變量后的拉格朗日函數(shù)的修正方程組; 4-5.求解所述修正方程,得到所述電力系統(tǒng)的潮流最優(yōu)解。8. 如權(quán)利要求7所述的方法,其特征在于,所述步驟4-5包括: (1) 初始化所述方程組:置迭代次數(shù)為0,給定自變量、松弛因子、拉格朗日乘子、類(lèi)擴(kuò)展 變量的初值; (2) 計(jì)算目標(biāo)函數(shù)和約束的雅可比矩陣及海森矩陣,形成線(xiàn)性系統(tǒng); (3) 求解所述線(xiàn)性系統(tǒng),得到每次迭代的修正量; (4) 計(jì)算原變量仿射步長(zhǎng)和對(duì)偶變量仿射步長(zhǎng); (5) 更新各變量及類(lèi)擴(kuò)展變量; (6) 根據(jù)設(shè)定的收斂標(biāo)準(zhǔn),判斷最優(yōu)潮流是否收斂; 若是,則停止計(jì)算,得到所述電力系統(tǒng)的潮流最優(yōu)解; 若否,則返回步驟(2)。
【文檔編號(hào)】H02J3/06GK105846429SQ201610363732
【公開(kāi)日】2016年8月10日
【申請(qǐng)日】2016年5月26日
【發(fā)明人】文晶, 丁平, 安寧, 湯涌, 王毅, 李芳 , 陳興雷, 黃彥浩, 楊小煜, 蔡靖, 孫璐, 何蕾, 李木, 李木一, 李文臣, 趙敏, 徐希望
【申請(qǐng)人】中國(guó)電力科學(xué)研究院, 國(guó)網(wǎng)遼寧省電力有限公司電力科學(xué)研究院, 國(guó)家電網(wǎng)公司