馬 靜，張俁妤，馬 偉，王增平

（華北電力大學(xué) 新能源電力系統(tǒng)國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京 102206）

0 引言

電力系統(tǒng)拓?fù)浞治鰯?shù)據(jù)主要由開關(guān)量組成，然而在開關(guān)量采集過(guò)程中，由于設(shè)備或網(wǎng)絡(luò)等原因，難免造成信息損壞或偏差，致使拓?fù)湫畔⒊霈F(xiàn)錯(cuò)誤［1-4］。

目前，國(guó)內(nèi)外學(xué)者針對(duì)電力系統(tǒng)拓?fù)溴e(cuò)誤辨識(shí)提出了諸多方法，如殘差法［4-5］、規(guī)則法［6-7］、最小信息損失法［8］、信息圖法［9-11］、轉(zhuǎn)移潮流法［12-14］等。 這些方法或?qū)顟B(tài)估計(jì)與拓?fù)浞治鱿嘟Y(jié)合，在狀態(tài)估計(jì)的過(guò)程中排除拓?fù)溴e(cuò)誤信息；或利用不良數(shù)據(jù)的特點(diǎn)，單獨(dú)進(jìn)行拓?fù)溴e(cuò)誤辨識(shí)，可有效辨識(shí)拓?fù)溴e(cuò)誤并修改不良數(shù)據(jù)。但傳統(tǒng)方法數(shù)據(jù)來(lái)源較單一，大多僅通過(guò)數(shù)據(jù)采集與監(jiān)視控制SCADA（Supervisory Control And Data Acquisition）系統(tǒng)，采用潮流來(lái)進(jìn)行分析，無(wú)法保證信息來(lái)自同一時(shí)間斷面［15］，且計(jì)算過(guò)程較復(fù)雜。

隨著智能電網(wǎng)的發(fā)展，同步相量測(cè)量技術(shù)及其相應(yīng)的廣域測(cè)量系統(tǒng)WAMS（Wide-Area Measurement System）［16-18］為電力系統(tǒng)拓?fù)浞治鎏峁┝诵碌募夹g(shù)手段。為此，本文借助相量測(cè)量單元PMU（Phasor Measurement Unit）提供的同步電氣量信息，提出了一種基于道路-回路方程的拓?fù)溴e(cuò)誤辨識(shí)方法。該方法首先利用網(wǎng)絡(luò)的拓?fù)浼s束與基爾霍夫定律建立道路-回路方程，計(jì)算支路電流；在此基礎(chǔ)上，考慮各種拓?fù)溴e(cuò)誤及不良數(shù)據(jù)并存的情況，根據(jù)計(jì)算值與量測(cè)值間的差值對(duì)拓?fù)湫畔⒑椭冯娏髁繙y(cè)值進(jìn)行校驗(yàn)；最后，通過(guò)IEEE 9與IEEE 39節(jié)點(diǎn)系統(tǒng)對(duì)該方法進(jìn)行了仿真驗(yàn)證，結(jié)果表明即使存在諸多拓?fù)溴e(cuò)誤及不良數(shù)據(jù)，也能實(shí)現(xiàn)準(zhǔn)確、可靠的網(wǎng)絡(luò)拓?fù)浔孀R(shí)。

1 道路-回路方程

1.1 道路矩陣［19］

電網(wǎng)中，道路是有方向的樹支，表示電流流向，故道路以樹的頂點(diǎn)為起點(diǎn)，經(jīng)支路流向樹的根。道路矩陣T（l×n階，l為支路數(shù)，n為節(jié)點(diǎn)數(shù)）中元素為：

以圖1所示的電網(wǎng)為例，選取支路L9作為連支，其他支路為樹支，并以母線B1、B3、B6作為樹的頂點(diǎn)，母線B2作為樹的根，道路用箭頭標(biāo)示。

圖1 IEEE 9節(jié)點(diǎn)系統(tǒng)道路-回路標(biāo)記圖Fig.1 Road-loop mark map of IEEE 9-bus system

由式（1）可得道路矩陣如式（2）所示。其中，與連支對(duì)應(yīng)的部分TL為零矩陣。

根據(jù)基本回路各樹支電流與支路電流方向，列寫回-支關(guān)聯(lián)矩陣B（L×l階矩陣，L為基本回路數(shù)），規(guī)定連支正方向?yàn)榛芈氛较?，則B元素為：

由式（3）可得，圖1中的回-支關(guān)聯(lián)矩陣B為：

1.2 道路-回路方程

支路電流由兩部分組成：一是節(jié)點(diǎn)注入電流對(duì)本支路的貢獻(xiàn)，二是回路電流對(duì)回路中支路的貢獻(xiàn)。若節(jié)點(diǎn)注入電流為IN，連支電流為IL，則支路電流為I′b（“′”代表其為計(jì)算值），其可由道路-回路方程表示為：

對(duì)于節(jié)點(diǎn)數(shù)目較多，且接線復(fù)雜大電網(wǎng)的局部網(wǎng)絡(luò)，道路-回路方程依然適用。假設(shè)節(jié)點(diǎn)i為局部網(wǎng)絡(luò)內(nèi)道路的起點(diǎn)，而支路k與支路t分別以節(jié)點(diǎn)i為終點(diǎn)和起點(diǎn)，則支路t的電流為：

其中，iNi為節(jié)點(diǎn)i的注入電流；Tt表示道路矩陣T的第t行；Bt和Bk分別表示回-支關(guān)聯(lián)矩陣B的第t行和第k列。

由式（6）可得，對(duì)于局部網(wǎng)絡(luò)道路的起點(diǎn)i，其等值節(jié)點(diǎn)注入電流i′Ni等于以該節(jié)點(diǎn)為終點(diǎn)的支路電流ilk與該節(jié)點(diǎn)外加電流iNi之和，減去支路k回路電流BkIL的結(jié)果。根據(jù)式（6）依次算出局部網(wǎng)絡(luò)起點(diǎn)的等值注入電流，與其他內(nèi)部節(jié)點(diǎn)的注入電流即可構(gòu)成局部網(wǎng)絡(luò)的IN，代入式（5）可得局部網(wǎng)絡(luò)中支路電流 I′b。

2 拓?fù)溴e(cuò)誤與不良數(shù)據(jù)對(duì)支路電流的影響

根據(jù)道路-回路方程求出支路電流后，再利用下述原則檢測(cè)系統(tǒng)拓?fù)溴e(cuò)誤和不良數(shù)據(jù)。

原則1：不存在拓?fù)溴e(cuò)誤與不良數(shù)據(jù)時(shí)，任一條支路電流的計(jì)算誤差其中，ε1取最大支路電流絕對(duì)值的 0.3%［20］，即

原則2：節(jié)點(diǎn)注入電流之和滿足KCL，即其中，kr為由實(shí)際測(cè)量裝置的允許誤差確定的可靠系數(shù)［20］（本文取1）；ns為求和節(jié)點(diǎn)數(shù)；0.2%為PMU量測(cè)量的相對(duì)誤差［21］；為求和節(jié)點(diǎn)注入電流最大絕對(duì)值。

原則3：若支路電流值則判斷支路為閉合，否則為斷開。 其中，ε3取40%的支路額定電流［22］，即 ε3=0.04ILN。

根據(jù)上述原則，利用道路-回路方程進(jìn)行拓?fù)溴e(cuò)誤辨識(shí)和不良數(shù)據(jù)的分析。在求解過(guò)程中，將電源電流和負(fù)荷電流視為節(jié)點(diǎn)注入電流。同時(shí)，計(jì)及支路對(duì)地導(dǎo)納，并將等值π型支路對(duì)地導(dǎo)納的電流均歸算為支路末端母線的注入電流。

2.1 拓?fù)溴e(cuò)誤的影響分析

電力網(wǎng)絡(luò)拓?fù)溴e(cuò)誤，一般分為開關(guān)狀態(tài)閉合而實(shí)際斷開、開關(guān)狀態(tài)斷開而實(shí)際閉合，以及雙母接線判斷錯(cuò)誤3種情況。

設(shè)節(jié)點(diǎn)i為拓?fù)溴e(cuò)誤支路k的末端節(jié)點(diǎn)，iΔk為該支路對(duì)地導(dǎo)納電流。將T中節(jié)點(diǎn)、IN中的節(jié)點(diǎn)電流以及支路按照距離道路根部的遠(yuǎn)近程度排序。若道路中節(jié)點(diǎn)i排在m個(gè)節(jié)點(diǎn)之后，則支路k排在m-1條支路之后。

（1）開關(guān)狀態(tài)閉合而實(shí)際斷開。

當(dāng)某條支路開關(guān)狀態(tài)閉合而實(shí)際斷開時(shí)，求取該支路末端節(jié)點(diǎn)注入電流將增大iΔk，此時(shí)道路-回路方程可表示為：

由式（7）可知，支路開關(guān)狀態(tài)為閉合而實(shí)際為斷開的情況下，拓?fù)溴e(cuò)誤支路的電流小于ε3，且道路中該支路之后的支路電流均不滿足原則1，且電流差值與拓?fù)溴e(cuò)誤支路的對(duì)地導(dǎo)納電流有關(guān)。

（2）開關(guān)狀態(tài)斷開而實(shí)際閉合。

當(dāng)某支路開關(guān)狀態(tài)為斷開而實(shí)際為閉合時(shí)，該支路末端節(jié)點(diǎn)注入電流將減小 iΔk，由式（8）得，拓?fù)渲烦霈F(xiàn)錯(cuò)誤后支路電流均不滿足原則1，且差值與拓?fù)溴e(cuò)誤支路電流及其對(duì)地電流有關(guān)。

若閉合連支的拓?fù)錉顟B(tài)為斷開，矩陣B受影響，導(dǎo)致回路中支路電流不滿足原則1，該差值與連支電流有關(guān)；此外，連支末端節(jié)點(diǎn)注入電流將受到影響，且道路中位于回路后的支路電流出現(xiàn)較小差值，仍不滿足原則1，該差值與連支對(duì)地導(dǎo)電流有關(guān)。

（3）雙母接線判斷錯(cuò)誤。

對(duì)于雙母連接的廠站出線，若母線連接判斷錯(cuò)誤：設(shè)節(jié)點(diǎn)i與節(jié)點(diǎn)u之間可以進(jìn)行倒閘操作，其出線為k，且兩母線間有s條支路。開關(guān)狀態(tài)為支路k與節(jié)點(diǎn)u連接，而實(shí)際支路k與母線i連接，根據(jù)式（9）可知，兩母線間的支路電流計(jì)算值與量測(cè)值之間均出現(xiàn)差值，且差值與支路k對(duì)地導(dǎo)納電流及其所在回路電流有關(guān)。

2.2 不良數(shù)據(jù)的影響分析

為了進(jìn)一步提高拓?fù)溴e(cuò)誤辨識(shí)能力，還需有效識(shí)別不良數(shù)據(jù)，并對(duì)其進(jìn)行修正。由于道路-回路方程中包含連支電流和節(jié)點(diǎn)電流，因此需先驗(yàn)證節(jié)點(diǎn)注入電流是否正確。首先，根據(jù)原則2，若節(jié)點(diǎn)注入電流出現(xiàn)不良數(shù)據(jù)，則節(jié)點(diǎn)注入電流之和大于等于ε2。根據(jù)式（7），出現(xiàn)不良數(shù)據(jù)節(jié)點(diǎn)的支路電流計(jì)算值與量測(cè)值之間均存在差值，差值恰好為不良數(shù)據(jù)誤差且可根據(jù)原則2修正節(jié)點(diǎn)注入電流。

若連支電流出現(xiàn)不良數(shù)據(jù)，由式（10）可得，回路中的支路電流計(jì)算值與量測(cè)值之間均出現(xiàn)差值，差值恰好也為不良數(shù)據(jù)誤差（其中，Δi為不良數(shù)據(jù)的誤差，Bk為B的第k列）。

若樹支電流出現(xiàn)不良數(shù)據(jù)，不影響計(jì)算值的求取，可直接利用計(jì)算值對(duì)量測(cè)值進(jìn)行修改。

3 基于道路-回路方程的拓?fù)溴e(cuò)誤辨識(shí)方法

根據(jù)第2節(jié)中的分析，可進(jìn)行開關(guān)量與電氣量一致性校驗(yàn)：首先對(duì)拓?fù)湫畔⑦M(jìn)行錯(cuò)誤辨識(shí)，再對(duì)不良數(shù)據(jù)進(jìn)行修正。具體步驟如下。

步驟 1：根據(jù)電力網(wǎng)絡(luò)形成 T、B、IN、IL，利用道路-回路方程求取I′b，判斷所有支路電流是否均滿足原則1與原則3。

步驟2：若出現(xiàn)差值，判斷是否有拓?fù)溴e(cuò)誤。

a.若某連支斷開，但其電流量測(cè)值大于ε3，且該回路中的支路電流不滿足原則1，則該連支實(shí)際為閉合；若某樹支斷開，但其電流計(jì)算值大于ε3，且道路上該樹支末端節(jié)點(diǎn)后的全部支路電流均不滿足原則1，則該樹支實(shí)為閉合。

b.若支路閉合，但其計(jì)算值小于ε3，且系統(tǒng)道路中該支路后的全部支路電流不滿足原則1，則該支路實(shí)為斷開。

c.若支路閉合，其計(jì)算值小于ε3，但將其支路對(duì)地導(dǎo)納電流減去后，節(jié)點(diǎn)電流之和增大，該支路拓?fù)湫畔⒄_，支路輕載。

d.對(duì)于可進(jìn)行倒閘操作的支路，若其計(jì)算值與量測(cè)值大于ε3且閉合，但兩母線間的支路電流不滿足原則1，則該支路母線連接判斷錯(cuò)誤。

根據(jù)上述方法檢測(cè)支路拓?fù)?，出現(xiàn)錯(cuò)誤進(jìn)行修改，更新 T、B、IN、IL，計(jì)算 I′b，再次判斷是否存在拓?fù)溴e(cuò)誤。

步驟3：無(wú)拓?fù)溴e(cuò)誤或拓?fù)溴e(cuò)誤已修正后，判斷是否存在不良數(shù)據(jù)。

a.求取系統(tǒng)節(jié)點(diǎn)注入電流之和，若不滿足原則2，且系統(tǒng)道路中某一電源或負(fù)荷節(jié)點(diǎn)后的支路電流均不滿足原則1，則該節(jié)點(diǎn)注入電流出現(xiàn)不良數(shù)據(jù)，根據(jù)原則2或根據(jù)差值進(jìn)行修正。

b.若某一回路中的支路電流不滿足原則1，則該回路連支電流出現(xiàn)不良數(shù)據(jù)，根據(jù)差值對(duì)連支電流進(jìn)行修正。

c.個(gè)別樹支電流不滿足原則1，即可用計(jì)算出的樹支電流對(duì)量測(cè)值進(jìn)行修正。

對(duì)不良數(shù)據(jù)修正后，繼續(xù)計(jì)算I′b，直至所有支路電流滿足原則1。

步驟4：拓?fù)溴e(cuò)誤辨識(shí)與不良數(shù)據(jù)修正完成。

4 算例分析

4.1 IEEE 9節(jié)點(diǎn)系統(tǒng)

基于PSCAD/EMTDC仿真軟件搭建了圖1所示的IEEE 9節(jié)點(diǎn)系統(tǒng)模型，并假設(shè)母線B2與B7可以倒閘操作。其中，系統(tǒng)電源100 MV·A、220 kV，支路對(duì)地導(dǎo)納均為1.623×10-4S，則ILN=0.223+j0.138 kA，ε3=0.00893+j0.00556 kA。節(jié)點(diǎn)注入電流如式（11）所示，回路電流為iL9=-0.43628+j0.216385 kA，依式（5）可得各支路電流計(jì)算結(jié)果如表1所示。

表1 IEEE 9節(jié)點(diǎn)系統(tǒng)中支路電流Table 1 Branch currents of IEEE 9-bus systemkA

由表1可以看出，系統(tǒng)各支路電流均滿足原則1，即不存在拓?fù)湫畔㈠e(cuò)誤和不良數(shù)據(jù)。

（1）開關(guān)閉合而實(shí)際斷開的情況。

以圖1所示系統(tǒng)中支路L6為例，假設(shè)支路兩端的斷路器均閉合，但L6實(shí)際處于斷線狀態(tài)。當(dāng)系統(tǒng)中支路均閉合時(shí)，其 T、B 矩陣分別如式（2）、（4）所示，節(jié)點(diǎn)注入電流如式（12）所示。

依據(jù)式（5），各支路電流計(jì)算結(jié)果如表2所示。ε1和 ε2分別取值為 ε1=（2.793+j1.324）×10-3kA，ε2=（1.676+j0.9734）×10-2kA。 由表 2可知，支路 L6的量測(cè)值與計(jì)算值均小于ε3，且其末端節(jié)點(diǎn)后支路電流不滿足原則1，因此判斷支路L6實(shí)為斷開。

修改支路L6的拓?fù)湫畔⒑?，系統(tǒng)道路如圖2所示，根據(jù)式（1）、（3）形成矩陣 T、B，節(jié)點(diǎn)注入電流為式（13），各支路電流計(jì)算結(jié)果如表3所示。分析表3可知，拓?fù)湫畔⑴c量測(cè)量均正確。

表2 支路均閉合情況下各支路電流Table 2 Branch currents when all branches are closed kA

圖2 L6斷開后系統(tǒng)道路-回路標(biāo)記圖Fig.2 Road-loop mark map with L6opened

表3 L6斷開后各支路電流Table3 Branch currents when L6is opened kA

（2）雙母線接線判斷錯(cuò)誤情況。

若系統(tǒng)拓?fù)浣Y(jié)構(gòu)如圖3所示，根據(jù)定義形成矩陣T、B，節(jié)點(diǎn)注入電流為式（14），各支路電流計(jì)算結(jié)果如表4所示。

圖3 L6與B2相連后道路-回路標(biāo)記圖Fig.3 Road-loop mark map with L6andB2connected

表4 L6與B2相連后計(jì)算結(jié)果Table 4 Calculative results when L6and B2are connected kA

分析表4可知，支路L6符合原則3，但L6之后可進(jìn)行倒閘操作的兩母線B2和B7間的支路L8電流不符合原則1，其差值與支路L6所在回路以及該線路對(duì)地導(dǎo)納電流差接近，判斷L6母線連接拓?fù)湫畔㈠e(cuò)誤。修改拓?fù)湫畔?，更新后的矩陣T、B以及節(jié)點(diǎn)注入電流如式（2）、（4）、（11）所示，各支路電流計(jì)算結(jié)果如表1所示，可知拓?fù)湫畔⑴c量測(cè)值均正確。

4.2 IEEE 39節(jié)點(diǎn)系統(tǒng)

基于PSCAD/EMTDC仿真軟件搭建IEEE 39節(jié)點(diǎn)系統(tǒng)進(jìn)一步驗(yàn)證本方法的正確性和有效性。系統(tǒng)結(jié)構(gòu)與道路標(biāo)記如圖4所示，其中“▲”箭頭標(biāo)記為連支，其余支路為樹支。取支路的最小額定值為0.2693+j0.1617 kA，則 ε3=0.01077+j0.00647 kA。

對(duì)于節(jié)點(diǎn)數(shù)目較多、接線復(fù)雜的網(wǎng)絡(luò)，可進(jìn)行局部拓?fù)溴e(cuò)誤辨識(shí)，本文僅以圖4所標(biāo)記的2個(gè)局部網(wǎng)絡(luò)（區(qū)域1和區(qū)域2）為例。

圖4 IEEE 39節(jié)點(diǎn)系統(tǒng)道路-回路標(biāo)記圖Fig.4 Road-loop mark map of IEEE 39-bus system

（1）多處不良數(shù)據(jù)存在的情況。

以節(jié)點(diǎn)為終點(diǎn)的局部網(wǎng)絡(luò)（區(qū)域1）如圖5所示。 根據(jù)式（1）、（3）形成矩陣 T、B，并根據(jù)式（6）分析，節(jié)點(diǎn)注入電流可用節(jié)點(diǎn)其余支路電流之和與該節(jié)點(diǎn)注入電流等值（如在回路中，除去其所在回路電流影響），則該局部網(wǎng)絡(luò)等值節(jié)點(diǎn)注入電流如式（15）所示，且連支電流 i26=0.6675+j0.9609kA，i27=-（0.2943+j0.2441）kA，各支路電流計(jì)算結(jié)果如表5所示。

圖5 以節(jié)點(diǎn)為終點(diǎn)的局部網(wǎng)絡(luò)道路-回路標(biāo)記圖Fig.5 Partial road-loop mark map with Bus as end

根據(jù)式（15）與表 5，ε1和 ε2分別取值為 ε1=（6.452+j9.652）×10-3kA、ε2=（4.3+j6.4）×10-2kA，節(jié)點(diǎn)總的輸出電流，全部支路電流均滿足原則3，無(wú)拓?fù)溴e(cuò)誤。

同時(shí)，根據(jù)局部網(wǎng)絡(luò)KCL，與節(jié)點(diǎn)注入電流之和的絕對(duì)差值大于ε2，支路19、20、25電流均不滿足原則1。其中，節(jié)點(diǎn)既不是負(fù)荷節(jié)點(diǎn)，也不是電源節(jié)點(diǎn)，則判斷節(jié)點(diǎn)注入電流出現(xiàn)差值?？紤]連支電流也可能有錯(cuò)誤，不能直接用KCL修正節(jié)點(diǎn)注入電流。比較造成各線路電流差值的因素，可得支路19差值減去支路25差值，即為節(jié)點(diǎn)注入電流的誤差值，運(yùn)用該差值對(duì)注入電流進(jìn)行修正得i12=-（0.6614+j0.5166）kA，各支路電流計(jì)算結(jié)果如表6所示。

表5 圖5網(wǎng)絡(luò)計(jì)算結(jié)果Table 5 Calculative results of fig.5 kA

表6 節(jié)點(diǎn)12注入電流不良數(shù)據(jù)修正后計(jì)算結(jié)果Table 6 Calculative results after bad current data of Bus 12 are corrected kA

由表6可以看出，與連支26相關(guān)的樹支電流均不滿足原則1，且局部網(wǎng)絡(luò)節(jié)點(diǎn)注入電流仍不滿足原則2，但兩者差值相等。初步判斷支路26電流錯(cuò)誤，運(yùn)用支路19的電流差值對(duì)其修正得i26=1.33495+j1.92183 kA。然后，運(yùn)用修正后的量測(cè)值進(jìn)行計(jì)算，結(jié)果如表7所示，電氣量均滿足原則1—3，無(wú)拓?fù)溴e(cuò)誤與不良數(shù)據(jù)。

表7 支路26電流不良數(shù)據(jù)修正后計(jì)算結(jié)果Table 7 Calculative results after bad current data of Bus 26 are corrected kA

對(duì)于疑似不良數(shù)據(jù)，可通過(guò)擴(kuò)大區(qū)域范圍，加入受不良數(shù)據(jù)影響的支路，進(jìn)一步提高拓?fù)溴e(cuò)誤辨識(shí)與不良數(shù)據(jù)修正的準(zhǔn)確性。

（2）開關(guān)狀態(tài)斷開而實(shí)際閉合的情況。

圖6為IEEE 39節(jié)點(diǎn)系統(tǒng)中局部網(wǎng)絡(luò)（區(qū)域2）。圖中，支路29斷開。根據(jù)定義形成矩陣T、B且該局部網(wǎng)絡(luò)節(jié)點(diǎn)注入電流如式（16）所示，各支路電流計(jì)算結(jié)果如表8所示。

圖6 IEEE 39節(jié)點(diǎn)系統(tǒng)中局部網(wǎng)絡(luò)道路-回路標(biāo)記圖Fig.6 Partial road-loop mark map of IEEE 39-bus system

表8 圖6網(wǎng)絡(luò)計(jì)算結(jié)果Table 8 Calculative results of fig.6 kA

根據(jù)圖6可知節(jié)點(diǎn)為道路起點(diǎn)，與本區(qū)域其他節(jié)點(diǎn)無(wú)關(guān)，需要單獨(dú)進(jìn)行KCL校驗(yàn)。該區(qū)域ε1=（7.988+j6.904）×10-3kA；對(duì)于節(jié)點(diǎn)，ε2=（4.01+j3.45）×10-3kA，其余節(jié)點(diǎn)求和閾值 ε2=（0.64+j1.59）×10-2kA。分析表8，可知連支29的量測(cè)值大于ε3，與拓?fù)湫畔⒉环?，且該回路中的支路電流均不滿足原則1，因此，判斷支路29拓?fù)溴e(cuò)誤。將矩陣B修改后，各支路電流計(jì)算結(jié)果如表9所示。

表9 支路29電流不良數(shù)據(jù)修正后計(jì)算結(jié)果Table 9 Calculative results after bad current data of Bus 29 are corrected kA

從表9可以看出，其支路電流符合原則1、3，且該局部網(wǎng)絡(luò)節(jié)點(diǎn)注入電流符合原則2。此時(shí)，該局部網(wǎng)絡(luò)中不存在拓?fù)溴e(cuò)誤和不良數(shù)據(jù)。

5 結(jié)論

本文提出一種基于道路-回路方程的拓?fù)溴e(cuò)誤辨識(shí)方法。該方法利用網(wǎng)絡(luò)拓?fù)浼s束和基爾霍夫定律，借助PMU提供的同步電氣量信息，建立道路-回路方程，計(jì)算支路電流。在考慮多種拓?fù)溴e(cuò)誤及不良數(shù)據(jù)的情況下，根據(jù)計(jì)算值與量測(cè)值間的差值對(duì)拓?fù)湫畔⒓爸冯娏髁繙y(cè)值進(jìn)行校驗(yàn)。基于IEEE 9和IEEE 39節(jié)點(diǎn)系統(tǒng)的仿真結(jié)果均表明，該方法不受拓?fù)溴e(cuò)誤與不良數(shù)據(jù)雙重因素的影響，辨識(shí)結(jié)果的準(zhǔn)確性及可靠性高。

［1］CLEMENT S K A，COSTA A S.Topology error identification using lagrange multipliers［J］.IEEE Trans on Power Systems，1998，13（2）：347-353.

［2］SINGH H，ALVARADO F L.Networktopologydetermination using least absolute value state estimation［J］.IEEE Trans on Power Systems，1995，10（3）：1159-1165.

［3］ABUR A，KIM H，CELIK M.Identifying the unknown circuit breaker status’s in power networks［J］.IEEE Trans on Power Systems，1995，10（4）：2029-2036.

［4］趙俊博，張葛祥，黃彥全.含新能源電力系統(tǒng)狀態(tài)估計(jì)研究現(xiàn)狀和展望［J］.電力自動(dòng)化設(shè)備，2014，34（5）：7-20，34.ZHAO Junbo，ZHANG Gexiang，HUANG Yanquan.Statusand prospect of state estimation for power system containing renewable energy［J］.Electric Power Automation Equipment，2014，34（5）：7-20，34.

［5］張伯明，王世纓，相年德.電力系統(tǒng)動(dòng)態(tài)狀態(tài)估計(jì)中不正常事件的處理［J］.中國(guó)電機(jī)工程學(xué)報(bào)，1993，13（3）：52-58.ZHANG Boming，WANG Shiying，XIANG Niande.Anomaly processing in power system dynamic state estimation［J］.Proceedings of the CSEE，1993，13（3）：52-58.

［6］SINGH N，GLAVITSCH H.Detection and identification of topological errors in online power system analysis［J］.IEEE Trans on Power Systems，1991，6（1）：324-331.

［7］黃滔，盧建剛，張輝.廣東省調(diào)提高狀態(tài)估計(jì)計(jì)算精度的措施［J］.電網(wǎng)技術(shù)，2004，28（16）：78-81.HUANG Tao，LU Jiangang，ZHANG Hui.Measures to improve state estimation calculation accuracy in Guangdong provincial power dispatching centre［J］.Power System Technology，2004，28（16）：78-81.

［8］孫宏斌，高峰，張伯明.最小信息損失狀態(tài)估計(jì)在拓?fù)溴e(cuò)誤辨識(shí)中的應(yīng)用［J］.中國(guó)電機(jī)工程學(xué)報(bào)，2005，25（18）：1-5 SUN Hongbin，GAO Feng，ZHANG Boming.Application of minimum information loss based state estimation to topology error identification［J］.Proceedings of the CSEE，2005，25（18）：1-5.

［9］周蘇荃，柳焯.新息圖法拓?fù)溴e(cuò)誤辨識(shí)［J］.電力系統(tǒng)自動(dòng)化，2000，24（4）：23-27.ZHOU Suquan，LIU Zhuo.An innovation graph approach totopology error identification［J］.Automation of Electric Power Systems，2000，24（4）：23-27.

［10］周蘇荃，柳焯.新息圖法識(shí)別多重網(wǎng)絡(luò)結(jié)構(gòu)動(dòng)態(tài)變化［J］.中國(guó)電機(jī)工程學(xué)報(bào)，2001，21（10）：67-72.ZHOU Suquan，LIU Zhuo.Identification of multiple dynamic network structure changes by the innovation graph approach［J］.Proceedings of the CSEE，2001，21（10）：67-72.

［11］周蘇荃，李昌，柳焯.利用回路新息相角差代數(shù)和識(shí)別拓?fù)浣Y(jié)構(gòu)變化［J］.電力自動(dòng)化設(shè)備，2002，22（6）：26-28.ZHOU Suquan，LI Chang，LIU Zhuo.Indentification of topology change by algebraic sum of innovation angle differences in a loop［J］.Electric Power Automation Equipment，2002，22（6）：26-28.

［12］陳艷波，周京陽(yáng)，于爾鏗，等.轉(zhuǎn)移潮流法拓?fù)溴e(cuò)誤辨識(shí)［J］.電力系統(tǒng)自動(dòng)化，2010，34（1）：20-24.CHEN Yanbo，ZHOU Jingyang，YU Erkeng，et al.Transfer power flow approach to topology error identification［J］.Automation of Electric Power Systems，2010，34（1）：20-24.

［13］陳艷波，何光宇，周京陽(yáng)，等.基于改進(jìn)轉(zhuǎn)移潮流法的拓?fù)溴e(cuò)誤辨識(shí)方法［J］.電網(wǎng)技術(shù)，2012，36（3）：95-100.CHEN Yanbo，HE Guangyu，ZHOU Jingyang，et al.An improved power flow transfer approach with enhanced ability to identify topology error and bad data［J］.Power System Technology，2012，36（3）：95-100.

［14］蘇慧玲，李揚(yáng).基于準(zhǔn)穩(wěn)態(tài)功率轉(zhuǎn)移分布因子的電力系統(tǒng)復(fù)雜網(wǎng)絡(luò)特性分析［J］.電力自動(dòng)化設(shè)備，2013，33（9）：47-53.SU Huiling，LI Yang.Analysis of complex network characteristics based on quasi-steady PTDF for power system［J］.Electric Power Automation Equipment，2013，33（9）：47-53.

［15］薛輝，賈清泉，王寧，等.基于PMU量測(cè)數(shù)據(jù)和SCADA數(shù)據(jù)融合的電力系統(tǒng)狀態(tài)估計(jì)方法［J］.電網(wǎng)技術(shù)，2008，32（14）：44-49.XUE Hui，JIA Qingquan，WANG Ning，et al.A novel power system state estimation method based on merging PMU-measured data into SCADA data［J］.Power System Technology，2008，32（14）：44-49.

［16］肖晉宇，謝小榮，李建，等.電網(wǎng)廣域動(dòng)態(tài)安全監(jiān)測(cè)系統(tǒng)及其動(dòng)態(tài)模擬試驗(yàn)［J］.電網(wǎng)技術(shù)，2004，28（6）：5-9.XIAO Jinyu，XIE Xiaorong，LI Jian，et al.Wide-area dynamic security monitoring system for power network and its dynamic simulation［J］.Power System Technology，2004，28（6）：5-9.

［17］蔡運(yùn)清，汪磊，MORISO K，等.廣域保護(hù)（穩(wěn)控）技術(shù)的現(xiàn)狀及展望［J］.電網(wǎng)技術(shù)，2004，28（8）：20-25.CAI J Y，WANG Lei，MORISO K，et al.Current status and prospect of wide-area protection（dynamic stability control） technologies［J］.Power System Technology，2004，28（8）：20-25.

［18］PHADKE A G，THORP J S，KARIMI K J，et al.State estimation with phasor measurements［C］∥Power Industry Computer Application Conference PICA 85.San Fransisco，USA：［s.n.］，1985：233-241.

［19］張伯明，陳壽孫，嚴(yán)正，等.高等電力系統(tǒng)分析［M］.2版.北京：清華大學(xué)出版社，2007.

［20］于爾鏗.電力系統(tǒng)狀態(tài)估計(jì)［M］.北京：水利電力出版社，1985：85-89.

［21］熊小伏，劉曉放.基于WAMS的繼電保護(hù)靜態(tài)特性監(jiān)視及其隱藏故障診斷［J］.電力系統(tǒng)自動(dòng)化，2009，33（9）：6-10.XIONG Xiaofu，LIU Xiaofang.Monitoring of static characteristics and diagnosis of hidden failures in protection system based on WAMS［J］.Automation of Electric Power Systems，2009，33（9）：6-10.

［22］陳建國(guó)，徐長(zhǎng)寶，高吉普，等.數(shù)字化變壓器保護(hù)電子式互感器的應(yīng)用與研究［J］.南方電網(wǎng)技術(shù)，2012，6（5）：105-109.

（continued on page 96）（continued from page 88）CHEN Jianguo，XU Changbao，GAO Jipu，et al.Application and research of the electronic transformer for the digital transformer protection［J］.Southern Power System Technology，2012，6（5）：105-109.