一種用于多電壓等級(jí)電網(wǎng)模型的抗差狀態(tài)估計(jì)方法
【專利摘要】本發(fā)明涉及一種用于多電壓等級(jí)電網(wǎng)模型的抗差狀態(tài)估計(jì)方法，所述方法包括讀取電網(wǎng)模型和SCADA量測(cè)；進(jìn)行網(wǎng)絡(luò)拓?fù)浞治觯瑢⑽锢黼娋W(wǎng)模型轉(zhuǎn)化成母線-支路計(jì)算模型；建立有功加權(quán)雅可比矩陣Ra-1/2Ba和無(wú)功計(jì)算加權(quán)雅可比矩陣Rr-1/2Br，計(jì)算有功系數(shù)矩陣A和無(wú)功系數(shù)矩陣B；利用正交變換獲得有功系數(shù)矩陣和無(wú)功系數(shù)矩陣的上三角矩陣L和L1，進(jìn)行狀態(tài)估計(jì)迭代求解；計(jì)算量測(cè)殘差vi，選擇量測(cè)類型基準(zhǔn)值Si，計(jì)算遙測(cè)估計(jì)值誤差計(jì)算量測(cè)權(quán)因子，判斷是否有量測(cè)處于降權(quán)區(qū)或淘汰區(qū)，沒(méi)有則狀態(tài)估計(jì)計(jì)算結(jié)束；有則根據(jù)量測(cè)權(quán)因子修正量測(cè)權(quán)重矩陣Ra-1和Rr-1，重新生成加權(quán)雅可比矩陣Ra-1/2Ba和Rr-1/2Br，計(jì)算A和B，返回步驟(4)。
【專利說(shuō)明】—種用于多電壓等級(jí)電網(wǎng)模型的抗差狀態(tài)估計(jì)方法

【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及一種抗差狀態(tài)估計(jì)方法，具體講涉及一種用于多電壓等級(jí)電網(wǎng)模型的抗差狀態(tài)估計(jì)方法。

【背景技術(shù)】
[0002]隨著智能電網(wǎng)的建設(shè)和發(fā)展，電網(wǎng)規(guī)模越來(lái)越大，電網(wǎng)建模更加精細(xì)化，對(duì)電力系統(tǒng)在線分析計(jì)算，尤其是狀態(tài)估計(jì)計(jì)算帶來(lái)了新的挑戰(zhàn)，如何從計(jì)算精度和計(jì)算速度兩個(gè)方面對(duì)狀態(tài)估計(jì)算法進(jìn)行提升，快速準(zhǔn)確的計(jì)算出大電網(wǎng)全網(wǎng)實(shí)時(shí)運(yùn)行狀態(tài)，為其他各類分析應(yīng)用提供準(zhǔn)確可靠的數(shù)據(jù)源，已經(jīng)成為大電網(wǎng)狀態(tài)估計(jì)計(jì)算所面臨的新問(wèn)題。因此，具有更強(qiáng)的不良數(shù)據(jù)辨識(shí)能力、更高計(jì)算精度和計(jì)算速度的抗差狀態(tài)估計(jì)模型與算法已經(jīng)成為電力系統(tǒng)研究中的一個(gè)熱點(diǎn)方向。
[0003]現(xiàn)有的M估計(jì)通過(guò)等價(jià)權(quán)原理化為加權(quán)最小二乘狀態(tài)估計(jì)的形式，引入權(quán)函數(shù)對(duì)量測(cè)權(quán)重進(jìn)行靈活調(diào)整，繼承了最小二乘的求解形式，在迭代計(jì)算過(guò)程中，根據(jù)殘差動(dòng)態(tài)調(diào)節(jié)量測(cè)權(quán)重，具有一定的抗差效果。但是，隨著我國(guó)大規(guī)?？稍偕茉吹募s化開發(fā)，特高壓交、直流輸電技術(shù)的快速發(fā)展，以及智能電網(wǎng)一體化運(yùn)行的深入開展，各級(jí)調(diào)度中心對(duì)在線分析和仿真計(jì)算的精度提出了更高的要求，近年來(lái)已有部分省級(jí)調(diào)度中心試驗(yàn)將省、地電網(wǎng)模型拼接，并基于拼接后的大電網(wǎng)模型進(jìn)行在線分析和調(diào)度員培訓(xùn)模擬。省、地電網(wǎng)模型拼接后電網(wǎng)模型電壓等級(jí)跨度非常大，不同電壓等級(jí)量測(cè)量的數(shù)值及其滿刻度值往往也相差數(shù)倍甚至數(shù)十倍?？共钭钚《斯烙?jì)中常用的權(quán)函數(shù)有Huber權(quán)函數(shù)、Hampel權(quán)函數(shù)、Turkey權(quán)函數(shù)、丹麥權(quán)函數(shù)、IGG權(quán)函數(shù)等?，F(xiàn)有的基于IGG權(quán)函數(shù)法的抗差狀態(tài)估計(jì)采用量測(cè)殘差絕對(duì)值作為判斷是否進(jìn)行降權(quán)處理的依據(jù)，對(duì)不同電壓等級(jí)的具有不同特點(diǎn)的量測(cè)沒(méi)有進(jìn)行區(qū)分，其對(duì)低電壓等級(jí)設(shè)備量測(cè)或小容量設(shè)備量測(cè)的壞數(shù)據(jù)辨識(shí)能力有限。
[0004]另外，在抗差狀態(tài)估計(jì)數(shù)學(xué)模型進(jìn)行求解迭代過(guò)程中，需要根據(jù)權(quán)重調(diào)整情況重新建立加權(quán)雅可比矩陣及正交分解，因此，基于等價(jià)權(quán)原理的抗差最小二乘狀態(tài)估計(jì)的計(jì)算效率較傳統(tǒng)加權(quán)最小二乘方法有所下降，限制了該算法在在線分析領(lǐng)域中的實(shí)用化進(jìn)程。
[0005]綜上所述，要將抗差最小二乘狀態(tài)估計(jì)算法更好的應(yīng)用到實(shí)際電力系統(tǒng)中去，還需要對(duì)抗差狀態(tài)估計(jì)算法進(jìn)行深入研究。


【發(fā)明內(nèi)容】

[0006]針對(duì)現(xiàn)有技術(shù)的不足，本發(fā)明提供一種用于多電壓等級(jí)電網(wǎng)模型的抗差狀態(tài)估計(jì)方法，在IGG分段權(quán)函數(shù)中引入量測(cè)類型基準(zhǔn)值，實(shí)現(xiàn)對(duì)不同電壓等級(jí)量測(cè)的區(qū)分，對(duì)處于降權(quán)區(qū)的量測(cè)權(quán)因子計(jì)算公式進(jìn)行修正，整體提高抗差狀態(tài)估計(jì)對(duì)不同電壓等級(jí)量測(cè)壞數(shù)據(jù)的抗差性能。
[0007]本發(fā)明的目的是采用下述技術(shù)方案實(shí)現(xiàn)的:
[0008]一種用于多電壓等級(jí)電網(wǎng)模型的抗差狀態(tài)估計(jì)方法，其改進(jìn)之處在于，所述方法包括
[0009](I)讀取電網(wǎng)模型和SCADA量測(cè)；
[0010](2)進(jìn)行網(wǎng)絡(luò)拓?fù)浞治?，將物理電網(wǎng)模型轉(zhuǎn)化成母線-支路計(jì)算模型；
[0011 ] (3)建立有功加權(quán)雅可比矩陣Ra_1/2Ba和無(wú)功計(jì)算加權(quán)雅可比矩陣Rr_1/2Br，計(jì)算有功系數(shù)矩陣A和無(wú)功系數(shù)矩陣B ；
[0012](4)利用正交變換獲得有功系數(shù)矩陣和無(wú)功系數(shù)矩陣的上三角矩陣L和L1，進(jìn)行狀態(tài)估計(jì)迭代求解；
[0013](5)計(jì)算量測(cè)殘差Vi,選擇量測(cè)類型基準(zhǔn)值Si,計(jì)算遙測(cè)估計(jì)值誤差
[0014](6)計(jì)算量測(cè)權(quán)因子，判斷是否有量測(cè)處于降權(quán)區(qū)或淘汰區(qū)，沒(méi)有則狀態(tài)估計(jì)計(jì)算結(jié)束；有則根據(jù)量測(cè)權(quán)因子修正量測(cè)權(quán)重矩陣RJ1和R,，重新生成加權(quán)雅可比矩陣Ra_1/2Ba和Rr-l72Br,計(jì)算A和B,返回步驟(4)。
[0015]優(yōu)選的，所述步驟(I)包括從調(diào)度控制系統(tǒng)模型庫(kù)獲取電網(wǎng)模型，從SCADA應(yīng)用獲取實(shí)時(shí)遙信遙測(cè)數(shù)據(jù)。
[0016]優(yōu)選的，所述步驟(2)包括根據(jù)電網(wǎng)模型電氣聯(lián)接關(guān)系和開關(guān)/刀閘運(yùn)行狀態(tài)，進(jìn)行廠站母線拓?fù)浞治?，形成閉合開關(guān)/刀閘相連的結(jié)點(diǎn)集合，分配計(jì)算母線編號(hào)；根據(jù)投運(yùn)支路聯(lián)接的計(jì)算母線集合，進(jìn)行電氣島劃分，判斷電氣島帶電狀態(tài)。
[0017]優(yōu)選的，所述步驟⑶包括
[0018](3.1)根據(jù)節(jié)點(diǎn)相關(guān)支路數(shù)進(jìn)行節(jié)點(diǎn)排序優(yōu)化，采用支路電抗形成導(dǎo)納矩陣，建立有功量測(cè)的雅可比矩陣Ba ;形成導(dǎo)納矩陣，建立無(wú)功量測(cè)的雅可比矩陣；
[0019](3.2)獲取系統(tǒng)量測(cè)初始權(quán)重設(shè)置，建立有功加權(quán)雅可比矩陣Ra_1/2Ba和無(wú)功計(jì)算加權(quán)雅可比矩陣R,_1/2B,;其中，R；1和Rf1分別為有功量測(cè)權(quán)重矩陣和無(wú)功量測(cè)權(quán)重矩陣；
[0020](3.3)根據(jù)公式 J = F04KUWA)]和 5 = 2KU' ΚΛ-Α)]分別計(jì)算有功系數(shù)矩陣A和計(jì)算無(wú)功系數(shù)矩陣B，其中，V0取系統(tǒng)參考節(jié)點(diǎn)電壓。
[0021]優(yōu)選的，所述步驟(4)包括
[0022](4.1)利用正交變換分別獲得A和B的因子表L和L1，置迭代計(jì)數(shù)器1 = 1;
[0023](4.2)按當(dāng)前狀態(tài)量電壓幅值V(1)和電壓相角Θ (1)計(jì)算有功量測(cè)ha，根據(jù)公式
αα> =K2(-Ba)f K1Iza計(jì)算a⑴;其中，a(1)為有功自由矢量，Za為有功量測(cè)矢量，ha為用狀態(tài)量V(1)和Θ (1)表達(dá)的有功量測(cè)；
[0024](4.3)利用A的因子表L由公式ΑΛ Θ (1) = 解出電壓相角修正矢量Δ Θ⑴，判斷I Λ Θ (1) Imax是否達(dá)到收斂精度；若不收斂，則計(jì)算Θ _ = θ (1)+Δ Θ (1)，若收斂，則有功收斂標(biāo)志位I ;
[0025](4.4)按當(dāng)前狀態(tài)量電壓幅值Va)和電壓相角θ ω計(jì)算無(wú)功量測(cè)匕，根據(jù)公式b11' =K 計(jì)算b(1);其中，b(1)為無(wú)功自由矢量，Zr為無(wú)功量測(cè)矢量，hr為用狀態(tài)量Va)和θ⑴表達(dá)的無(wú)功量測(cè)；
[0026](4.5)利用B的因子表L1由公式BAV(1) = b(1)解出電壓幅值修正矢量Λ V⑴，判斷I λ Va) Imax是否達(dá)到收斂精度；若不收斂，則計(jì)算V(1+1) =V(1)+AVa)，若收斂，則無(wú)功收斂標(biāo)志位I ;
[0027](4.6)判斷有功和無(wú)功收斂標(biāo)志位是否都為1，為I則狀態(tài)估計(jì)迭代收斂，否則迭代計(jì)數(shù)器I = 1+1 ;判斷是否達(dá)到最大迭代次數(shù)，達(dá)到最大迭代次數(shù)，則狀態(tài)估計(jì)不收斂，計(jì)算結(jié)束；未達(dá)到最大迭代次數(shù)則返回步驟(4.2)。
[0028]進(jìn)一步地,采用Givens正交三角化的列超前并行消去算法求解L和L1,采用4進(jìn)程并行模式，所述步驟(4.1)包括
[0029](4.1.1)記加權(quán)雅可比矩陣Ra_1/2Ba為Bw，Givens變換將Bw的某一元素bn零化，

可以看成矩陣Bw的兩行左乘一個(gè)二階平面旋轉(zhuǎn)陣，即

【權(quán)利要求】
1.一種用于多電壓等級(jí)電網(wǎng)模型的抗差狀態(tài)估計(jì)方法，其特征在于，所述方法包括 (1)讀取電網(wǎng)模型和SCADA量測(cè)； (2)進(jìn)行網(wǎng)絡(luò)拓?fù)浞治?，將物理電網(wǎng)模型轉(zhuǎn)化成母線-支路計(jì)算模型； (3)建立有功加權(quán)雅可比矩陣Ra_1/2Ba和無(wú)功計(jì)算加權(quán)雅可比矩陣R,/2B?計(jì)算有功系數(shù)矩陣A和無(wú)功系數(shù)矩陣B; (4)利用正交變換獲得有功系數(shù)矩陣和無(wú)功系數(shù)矩陣的上三角矩陣L和L1,進(jìn)行狀態(tài)估計(jì)迭代求解； (5)計(jì)算量測(cè)殘差Vi，選擇量測(cè)類型基準(zhǔn)值Si,計(jì)算遙測(cè)估計(jì)值誤差P; (6)計(jì)算量測(cè)權(quán)因子，判斷是否有量測(cè)處于降權(quán)區(qū)或淘汰區(qū)，沒(méi)有則狀態(tài)估計(jì)計(jì)算結(jié)束；有則根據(jù)量測(cè)權(quán)因子修正量測(cè)權(quán)重矩陣R,和R,，重新生成加權(quán)雅可比矩陣1^_1/283和Rfv2Bp計(jì)算A和B，返回步驟(4)。
2.如權(quán)利要求1所述的一種用于多電壓等級(jí)電網(wǎng)模型的抗差狀態(tài)估計(jì)方法，其特征在于，所述步驟(I)包括從調(diào)度控制系統(tǒng)模型庫(kù)獲取電網(wǎng)模型，從SCADA應(yīng)用獲取實(shí)時(shí)遙信遙測(cè)數(shù)據(jù)。
3.如權(quán)利要求1所述的一種用于多電壓等級(jí)電網(wǎng)模型的抗差狀態(tài)估計(jì)方法，其特征在于，所述步驟(2)包括根據(jù)電網(wǎng)模型電氣聯(lián)接關(guān)系和開關(guān)/刀閘運(yùn)行狀態(tài)，進(jìn)行廠站母線拓?fù)浞治?，形成閉合開關(guān)/刀閘相連的結(jié)點(diǎn)集合，分配計(jì)算母線編號(hào)；根據(jù)投運(yùn)支路聯(lián)接的計(jì)算母線集合，進(jìn)行電氣島劃分，判斷電氣島帶電狀態(tài)。
4.如權(quán)利要求1所述的一種用于多電壓等級(jí)電網(wǎng)模型的抗差狀態(tài)估計(jì)方法，其特征在于，所述步驟(3)包括 (3.1)根據(jù)節(jié)點(diǎn)相關(guān)支路數(shù)進(jìn)行節(jié)點(diǎn)排序優(yōu)化，采用支路電抗形成導(dǎo)納矩陣，建立有功量測(cè)的雅可比矩陣Ba ;形成導(dǎo)納矩陣，建立無(wú)功量測(cè)的雅可比矩陣； (3.2)獲取系統(tǒng)量測(cè)初始權(quán)重設(shè)置，建立有功加權(quán)雅可比矩陣Ra_1/2Ba和無(wú)功計(jì)算加權(quán)雅可比矩陣R,_1/2B,;其中，R；1和Rf1分別為有功量測(cè)權(quán)重矩陣和無(wú)功量測(cè)權(quán)重矩陣；
(3.3)根據(jù)公式/i = K)[(-5 Y K1 (-Bu)]和S = JV [(-5 V K ' (-Br)]分別計(jì)算有功系數(shù)矩陣A和計(jì)算無(wú)功系數(shù)矩陣B，其中，Vtl取系統(tǒng)參考節(jié)點(diǎn)電壓。
5.如權(quán)利要求1所述的一種用于多電壓等級(jí)電網(wǎng)模型的抗差狀態(tài)估計(jì)方法，其特征在于，所述步驟(4)包括 (4.1)利用正交變換分別獲得A和B的因子表L和L1，置迭代計(jì)數(shù)器1 = 1; (4.2)按當(dāng)前狀態(tài)量電壓幅值Vq)和電壓相角θ⑴計(jì)算有功量測(cè)ha，根據(jù)公式計(jì)算a(1);其中，a(1)為有功自由矢量，Za為有功量測(cè)矢量，ha為用狀態(tài)量V(1)和Θ (1)表達(dá)的有功量測(cè)； (4.3)利用A的因子表L由公式ΑΛθ(1) =a(1)解出電壓相角修正矢量Λ θ(1)，判斷Δ Θ (1) |max是否達(dá)到收斂精度；若不收斂，則計(jì)算Θ (1+1) = θ (1)+Δ Θ (1)，若收斂，則有功收斂標(biāo)志位I ; (4.4)按當(dāng)前狀態(tài)量電壓幅值Vq)和電壓相角θ⑴計(jì)算無(wú)功量測(cè)tv根據(jù)公式,([ζ,.-?"1/#)]計(jì)算b(1);其中，b⑴為無(wú)功自由矢量，％為無(wú)功量測(cè)矢量，hr為用狀態(tài)量νω和θ (1)表達(dá)的無(wú)功量測(cè)； (4.5)利用B的因子表L1由公式BAVa) = b(1)解出電壓幅值修正矢量Λν(1)，判斷AV(1)L.是否達(dá)到收斂精度；若不收斂，則計(jì)算V(1+1) = V(1)+AVa)，若收斂，則無(wú)功收斂標(biāo)志位I ; (4.6)判斷有功和無(wú)功收斂標(biāo)志位是否都為1，為I則狀態(tài)估計(jì)迭代收斂，否則迭代計(jì)數(shù)器I = 1+1 ;判斷是否達(dá)到最大迭代次數(shù)，達(dá)到最大迭代次數(shù)，則狀態(tài)估計(jì)不收斂，計(jì)算結(jié)束；未達(dá)到最大迭代次數(shù)則返回步驟(4.2)。
6.如權(quán)利要求5所述的一種用于多電壓等級(jí)電網(wǎng)模型的抗差狀態(tài)估計(jì)方法，其特征在于，采用Givens正交三角化的列超前并行消去算法求解L和L1，采用4進(jìn)程并行模式，所述步驟(4.1)包括 (4.1.1)記加權(quán)雅可比矩陣Ra_1/2Ba為Bw，Givens變換將Bw的某一元素bu零化，可以


(c s\( β \ ( β'Λ看成矩陣Bw的兩行左乘一個(gè)二階平面旋轉(zhuǎn)陣，S卩— ==〗，其中βρ?為~矩


\s 0■八 A y \Pi)陣的 i 行和 j 行，i, j ^ {I, ---,ml, i ^ j,c = sin Θ u，s = cos θ Θ。表示旋轉(zhuǎn)的角度；由β/中的bu ’ =0，可確定(3和8;0/的元素和β/的其余元素為:
b，Jt = cbJt+sbit, (t = I, *..η)
b，it = -sbJt+cbit, (t = I,...，n, t 關(guān) I)；(4.1.2) Bw矩陣分塊存儲(chǔ)，Bw矩陣儲(chǔ)存在進(jìn)程O中，將Bw矩陣按行分解為4塊，序號(hào)為O、1、2、3，將第1、2、3塊分別發(fā)送給進(jìn)程1、2和3，進(jìn)程O保存并處理第O塊Bw矩陣數(shù)據(jù)；(4.1.3)進(jìn)程OWBw矩陣取一行數(shù)據(jù)，經(jīng)消去操作后，將此行數(shù)據(jù)發(fā)送給進(jìn)程1，進(jìn)程I完成相應(yīng)的列消去操作后，再將此行數(shù)組發(fā)送給進(jìn)程2，以此類推；進(jìn)程O將數(shù)據(jù)發(fā)送后可再次從Bw矩陣取下一行數(shù)據(jù)，進(jìn)行消去操作，并重復(fù)上述過(guò)程，以實(shí)現(xiàn)并行計(jì)算； (4.1.4)進(jìn)程3完成列消去操作后，將此行數(shù)組保存到L矩陣當(dāng)中； (4.1.5)全部消去完成后，進(jìn)程3將最終生成的L矩陣發(fā)送給進(jìn)程O。
7.如權(quán)利要求1所述的一種用于多電壓等級(jí)電網(wǎng)模型的抗差狀態(tài)估計(jì)方法，其特征在于，所述步驟(5)包括 (5.1)根據(jù)公式Vi = Z1-hi(V，Θ )計(jì)算各量測(cè)殘差Vi ; (5.2)根據(jù)量測(cè)電壓等級(jí)和量測(cè)設(shè)備類型選擇量測(cè)類型基準(zhǔn)值Si，根據(jù)Si計(jì)算量測(cè)估、山、口^ I計(jì)祆差]r。
8.如權(quán)利要求1所述的一種用于多電壓等級(jí)電網(wǎng)模型的抗差狀態(tài)估計(jì)方法，其特征在于，所述步驟(6)包括根據(jù)改進(jìn)IGG權(quán)函數(shù)
計(jì)算所有量測(cè)權(quán)因子，其中， 為保權(quán)區(qū)，矣為降權(quán)區(qū)，為淘汰區(qū)；判斷是否有量測(cè)處于降權(quán)區(qū)或淘汰區(qū)，沒(méi)有則狀態(tài)估計(jì)計(jì)算結(jié)束；有則根據(jù)量測(cè)權(quán)因子修正量測(cè)權(quán)重矩陣R,和R,，重新生成有功加權(quán)雅可比矩陣Ra_1/2Ba和無(wú)功計(jì)算加權(quán)雅可比矩陣R,/2Bp計(jì)算有功系數(shù)矩陣A和計(jì)算無(wú)功系數(shù)矩陣B，返回步驟(4)，直至沒(méi)有量測(cè)處于降權(quán)區(qū)或淘汰區(qū)，則狀態(tài)估計(jì)結(jié)束。
9.如權(quán)利要求8所述的一種用于多電壓等級(jí)電網(wǎng)模型的抗差狀態(tài)估計(jì)方法，其特征在
于，所述公式
的變換過(guò)程如下: (9.1)166權(quán)函數(shù)
采用量測(cè)殘差絕對(duì)值對(duì)權(quán)函數(shù)進(jìn)行分區(qū)，
引入量測(cè)類型基準(zhǔn)值Si，計(jì)算遙測(cè)數(shù)據(jù)估計(jì)值誤差Vi/Si，對(duì)權(quán)函數(shù)進(jìn)行分區(qū)，區(qū)分了不同電壓等級(jí)的量測(cè)；其中，kl，k2 e R，0〈kl〈k2，σ為電力系統(tǒng)遙測(cè)量的標(biāo)準(zhǔn)誤差，| ν |為量測(cè)殘差絕對(duì)值；
(9.2) IGG權(quán)函數(shù)
通過(guò)等價(jià)權(quán)化作用于加權(quán)最小二乘法目標(biāo)
函數(shù)時(shí)，處于降權(quán)區(qū)的量測(cè)的目標(biāo)函數(shù)分量為Ji (X) = Tik1 O |v|,修正原有權(quán)函數(shù)中處于降權(quán)區(qū)的量測(cè)權(quán)因子為.^~ *
K(3)通過(guò)(9.1)和(9.2)得到改進(jìn)的IGG權(quán)函數(shù)為
其中，￥為遙測(cè)數(shù)據(jù)估計(jì)值誤差，遙測(cè)數(shù)據(jù)估計(jì)值誤差門檻k1; k2取值根據(jù)實(shí)際系統(tǒng)量測(cè)情況設(shè)置，K大于0.03，k2大于k1D
【文檔編號(hào)】G06F19/00GK104184144SQ201410450017
【公開日】2014年12月3日 申請(qǐng)日期:2014年9月5日 優(yōu)先權(quán)日:2014年9月5日
【發(fā)明者】徐杰, 郭凌旭, 張志君, 徐晟 , 范廣民, 羅雅迪, 趙昆, 鄒煜, 郎燕生, 賈育培, 王磊, 王少芳, 竇成龍 申請(qǐng)人:國(guó)家電網(wǎng)公司, 中國(guó)電力科學(xué)研究院, 國(guó)網(wǎng)天津市電力公司