1　引言 　　電力系統(tǒng)短路故障是一個(gè)復(fù)雜的電磁暫態(tài)過程，應(yīng)用對稱分量法分析計(jì)算三相對稱電力系統(tǒng)的不對稱短路或非對稱斷相故障時(shí)，總是將故障端口處三相不對稱的支路阻抗參數(shù)等值轉(zhuǎn)換成三個(gè)互相獨(dú)立的正、負(fù)、零序三相對稱的電流源，從而將難于求解的三相不對稱故障的計(jì)算，轉(zhuǎn)換成易于求解的三個(gè)互相獨(dú)立而又對稱的三相電路網(wǎng)絡(luò)的計(jì)算。 　　工程實(shí)用計(jì)算中，視各序網(wǎng)絡(luò)為線性網(wǎng)絡(luò)，如用下標(biāo)（S）表示序列，n（S）為各序網(wǎng)絡(luò)的獨(dú)立節(jié)點(diǎn)數(shù)，應(yīng)用對稱分量法計(jì)算故障時(shí)，需解算獨(dú)立節(jié)點(diǎn)數(shù)為n（S）的各序大型電路網(wǎng)絡(luò)，計(jì)算工作量相當(dāng)大。在很多情況下我們感興趣的可能是大型網(wǎng)絡(luò)中的某個(gè)局部網(wǎng)絡(luò)，其關(guān)聯(lián)節(jié)點(diǎn)數(shù)r（S）遠(yuǎn)小于n（S），或者為了加快故障的計(jì)算速度，或者受所用計(jì)算機(jī)容量的限制，常常需要將大型網(wǎng)絡(luò)節(jié)點(diǎn)方程縮減成預(yù)定節(jié)點(diǎn)數(shù)目的等值節(jié)點(diǎn)方程。本文提出了將REI等值法［1］用于多節(jié)點(diǎn)配電系統(tǒng)短路電流計(jì)算，對擴(kuò)展的REI等值網(wǎng)絡(luò)應(yīng)用局部節(jié)點(diǎn)等值方程進(jìn)行了研究，使多節(jié)點(diǎn)配電系統(tǒng)短路電流計(jì)算的數(shù)學(xué)模型得到簡化，短路電流計(jì)算程序更簡捷、實(shí)用。 [b]2　局部節(jié)點(diǎn)等值方程 [/b]　　局部節(jié)點(diǎn)等值方程是將大型網(wǎng)絡(luò)節(jié)點(diǎn)分為兩大類：保留節(jié)點(diǎn)r（S）及消去節(jié)點(diǎn)e（S），e（S）＝n（S）－r（S）。計(jì)算時(shí)，在各序網(wǎng)中先對保留節(jié)點(diǎn)編號，后對消去節(jié)點(diǎn)編號，則其節(jié)點(diǎn)導(dǎo)納方程Y（S）U（S）＝I（S）按保留節(jié)點(diǎn)、消去節(jié)點(diǎn)進(jìn)行分塊，可得 　 ?。?） 　　式中　Y（S）為節(jié)點(diǎn)導(dǎo)納矩陣；U（S）為節(jié)點(diǎn)電壓列矢量；I（S）為節(jié)點(diǎn)注入電流列矢量；Yrr（S）包括保留節(jié)點(diǎn)的自導(dǎo)納和互導(dǎo)納；Yre（S）、Yer（S）為保留節(jié)點(diǎn)與消去節(jié)點(diǎn)之間的互導(dǎo)納；Yee（S）包括消去節(jié)點(diǎn)的自導(dǎo)納和互導(dǎo)納；Ur（S）是保留節(jié)點(diǎn)電壓列矢量；Ue（S）為消去節(jié)點(diǎn)電壓列矢量；Ir（S）為保留節(jié)點(diǎn)注入電流列矢量；Ie（S）為消去節(jié)點(diǎn)注入電流列矢量。 　　由于在等值縮減的前后，Ur（S）應(yīng)保持不變，因此應(yīng)用高斯消去法消去e（S）個(gè)節(jié)點(diǎn)后得保留節(jié)點(diǎn)的等值導(dǎo)納方程為 　?。╕rr（S）－Yre（S）Y－1ee（S）Yer（S））Ur（S）＝Ir（S）－Yre（S）Y－1ee（S）Ie（S） 　?。?） 　　簡記為　YNrr（S）Ur（S）＝INr（S）　　　　　　　　　　　　　　　　　　　　　　　　　　?。?） 　　式中　YNrr（S）為保留節(jié)點(diǎn)的等值導(dǎo)納矩陣；INr（S）為等值電流源［2］。 　　YNrr（S）＝Y(jié)rr（S）－Yre（S）Y－1ee（S）Yer（S）　　　　　　　　　　　　?。?） 　　INr（S）＝Ir（S）－Yre（S）Y－1ee（S）Ie（S）　　　　　　　　　　　　?。?） 　　由式（4）和式（5）可知，當(dāng)被消節(jié)點(diǎn)是有源節(jié)點(diǎn)，即Ie（S）不等于零時(shí)，則原網(wǎng)絡(luò)Y（S）對應(yīng)的保留節(jié)點(diǎn)的子塊Yrr（S）及保留節(jié)點(diǎn)的電源Ir（S）均需修正；如被消節(jié)點(diǎn)是無源節(jié)點(diǎn)即Ie（S）等于零時(shí)，則勿需修正保留節(jié)點(diǎn)的電流源。實(shí)用工程計(jì)算中，故障前認(rèn)為系統(tǒng)處于三相對稱運(yùn)行狀態(tài)，所以負(fù)序網(wǎng)的Ie（2）、零序網(wǎng)的Ie（o）均為零。 　　如果應(yīng)用式（4）與式（5）將Yrr（S）修正成YNrr（S），將Ir（S）修正成INr（S），對多節(jié)點(diǎn)配電網(wǎng)絡(luò)短路電流計(jì)算將會出現(xiàn)兩種情況：（1）按局部節(jié)點(diǎn)等值方程編計(jì)算程序計(jì)算YNrr（S）與INr（S）時(shí)需輸入全網(wǎng)的原始數(shù)據(jù)，工作量很大；（2）在只有微機(jī)的條件下計(jì)算具有幾百個(gè)節(jié)點(diǎn)的配電網(wǎng)絡(luò)時(shí)由于計(jì)算機(jī)容量的限制，分塊前全網(wǎng)節(jié)點(diǎn)導(dǎo)納矩陣Y（S）就無法形成，因而無法分塊，程序無法進(jìn)行。 　　如果將REI等值法與局部節(jié)點(diǎn)等值法結(jié)合運(yùn)用，將使多節(jié)點(diǎn)配電網(wǎng)絡(luò)的短路電流計(jì)算問題得到進(jìn)一步簡化。 [b]3　將REI等值法用于短路電流計(jì)算 [/b]　　REI等值法為P.Dimo等人提出并應(yīng)用于發(fā)達(dá)國家電力系統(tǒng)潮流計(jì)算與安全分析［3］。它將電網(wǎng)的節(jié)點(diǎn)分為兩組：應(yīng)保留的節(jié)點(diǎn)與要消去的節(jié)點(diǎn)。首先將要消去的節(jié)點(diǎn)中的有源節(jié)點(diǎn)按其相關(guān)性質(zhì)歸并為若干組，每組用一個(gè)虛擬的等價(jià)有源節(jié)點(diǎn)來代替，并通過一個(gè)無損耗的虛構(gòu)網(wǎng)絡(luò)（REI網(wǎng)絡(luò)）與其相聯(lián)。虛擬有源節(jié)點(diǎn)上的有功、無功注入功率是該組有源節(jié)點(diǎn)有功與無功功率的代數(shù)和。在接入REI網(wǎng)絡(luò)與虛擬等價(jià)節(jié)點(diǎn)后，原來的有源節(jié)點(diǎn)就變成無源節(jié)點(diǎn)了，然后將所有要消去的無源節(jié)點(diǎn)用常規(guī)的方法消去。等值變化過程如圖1所示。 　　 圖1　REI等值網(wǎng)絡(luò)簡化過程 　 　Fig.1　Simplified processing of REI equivalent networ 　REI網(wǎng)絡(luò)中各個(gè)導(dǎo)納的數(shù)值用如下方法確定： 　　對每個(gè)要消去的有源節(jié)點(diǎn)，其注入電流關(guān)系式為 　[img=199,25]http://zszl.cepee.com/cepee_kjlw_pic/files/wx/zgdjgcxb/zgdj2000/0004/image4/65-2.gif[/img]　（6） 　　式中　[img=11,22]http://zszl.cepee.com/cepee_kjlw_pic/files/wx/zgdjgcxb/zgdj2000/0004/image4/65-3.gif[/img][sub]k[/sub]為要消去的有源節(jié)點(diǎn)注入電流；[img=18,24]http://zszl.cepee.com/cepee_kjlw_pic/files/wx/zgdjgcxb/zgdj2000/0004/image4/65-4.gif[/img]為原始網(wǎng)絡(luò)有源節(jié)點(diǎn)注入功率的共軛值；[img=23,24]http://zszl.cepee.com/cepee_kjlw_pic/files/wx/zgdjgcxb/zgdj2000/0004/image4/65-5.gif[/img]為基本潮流解的節(jié)點(diǎn)復(fù)電壓共扼值。 　[img=149,50]http://zszl.cepee.com/cepee_kjlw_pic/files/wx/zgdjgcxb/zgdj2000/0004/image4/65-6.gif[/img]?。?） 　　式中　[img=11,22]http://zszl.cepee.com/cepee_kjlw_pic/files/wx/zgdjgcxb/zgdj2000/0004/image4/65-3.gif[/img][sub]R[/sub]為虛擬有源節(jié)點(diǎn)的注入電流。 　　在構(gòu)造REI網(wǎng)絡(luò)時(shí)參數(shù)應(yīng)保持原始網(wǎng)絡(luò)各有源節(jié)點(diǎn)的注入功率不變，因此 　[img=264,27]http://zszl.cepee.com/cepee_kjlw_pic/files/wx/zgdjgcxb/zgdj2000/0004/image4/65-7.gif[/img]?。?） 　　式中　Yk為REI等值網(wǎng)絡(luò)中，節(jié)點(diǎn)G與要消去節(jié)點(diǎn)之間的導(dǎo)納。 　　為滿足無損網(wǎng)的條件，則 　[img=89,38]http://zszl.cepee.com/cepee_kjlw_pic/files/wx/zgdjgcxb/zgdj2000/0004/image4/66-1.gif[/img]?。?） 　[img=220,53]http://zszl.cepee.com/cepee_kjlw_pic/files/wx/zgdjgcxb/zgdj2000/0004/image4/66-2.gif[/img]　（10） 　[img=145,27]http://zszl.cepee.com/cepee_kjlw_pic/files/wx/zgdjgcxb/zgdj2000/0004/image4/66-3.gif[/img]?。?1） 　　式中[img=51,24]http://zszl.cepee.com/cepee_kjlw_pic/files/wx/zgdjgcxb/zgdj2000/0004/image4/66-4.gif[/img]　分別為虛擬有源節(jié)點(diǎn)R的注入功率和復(fù)電壓；YR為虛擬有源節(jié)點(diǎn)R與節(jié)點(diǎn)G之間的導(dǎo)納。 　　式（8）中的UG可以是任意的，取UG＝0，于是REI網(wǎng)絡(luò)的構(gòu)造就變成了唯一的，由此得 　[img=269,49]http://zszl.cepee.com/cepee_kjlw_pic/files/wx/zgdjgcxb/zgdj2000/0004/image4/66-5.gif[/img]　?。?）′ 　[img=265,26]http://zszl.cepee.com/cepee_kjlw_pic/files/wx/zgdjgcxb/zgdj2000/0004/image4/66-6.gif[/img]?。?1）′ 　　構(gòu)造出REI網(wǎng)絡(luò)后要消去那些不感興趣的節(jié)點(diǎn)，假定擴(kuò)展了REI網(wǎng)絡(luò)后其網(wǎng)絡(luò)導(dǎo)納矩陣是Y，以下標(biāo)E表示要消去的節(jié)點(diǎn)集，以下標(biāo)I表示要保留的節(jié)點(diǎn)集，于是Y可寫成： 　[img=112,43]http://zszl.cepee.com/cepee_kjlw_pic/files/wx/zgdjgcxb/zgdj2000/0004/image4/66-7.gif[/img]?。?2） 　　消去節(jié)點(diǎn)集E中所有的節(jié)點(diǎn)，得到由節(jié)點(diǎn)集Ⅰ中節(jié)點(diǎn)所組成的簡化網(wǎng)絡(luò)，如圖1（c）所示。經(jīng)過外部等值后的節(jié)點(diǎn)導(dǎo)納矩陣為 　　YIIN＝Y(jié)II－YIEYEE－1YEI　　　　　　　　　　　　　　　?。?3） 　　等值電流源為 　　INI＝II－YIEYEE－1IE　　　　　　　　　　　　　　　?。?4） 　　保留節(jié)點(diǎn)的等值導(dǎo)納方程簡記為 　　YNIIUI＝INI　　　　　　　　　　　　　　　　（15） 　　將這一理論應(yīng)用到短路電流計(jì)算中，我們看到，式（13）～（15）與式（3）～（5）分別對應(yīng)一致，不同的是接入REI網(wǎng)絡(luò)與虛擬等價(jià)節(jié)點(diǎn)后，原來的有源節(jié)點(diǎn)就變成了無源節(jié)點(diǎn)，將虛擬等價(jià)節(jié)點(diǎn)R作為保留節(jié)點(diǎn)，則IE均為零。式（14）成為 　　INI＝II　　　　　　　　　　　　　　　　（16） 　　由式（8）與式（11）可見，REI網(wǎng)絡(luò)的參數(shù)與網(wǎng)絡(luò)的運(yùn)行參數(shù)UK有關(guān)，在潮流計(jì)算與安全分析中，只在基本運(yùn)行方式下才滿足和原網(wǎng)絡(luò)相互等值的關(guān)系。當(dāng)系統(tǒng)運(yùn)行情況偏離于基本運(yùn)行方式時(shí)，如果仍保持REI網(wǎng)絡(luò)參數(shù)不變，就會出現(xiàn)誤差，但在短路電流計(jì)算中，由于采用實(shí)用算法，取UK為K點(diǎn)短路前的額定電壓，其標(biāo)么值為1，則式（6）～（11）′可改變?yōu)槭剑?7）～（24）: 　[img=200,25]http://zszl.cepee.com/cepee_kjlw_pic/files/wx/zgdjgcxb/zgdj2000/0004/image4/66-8.gif[/img]　（17） 　[img=86,36]http://zszl.cepee.com/cepee_kjlw_pic/files/wx/zgdjgcxb/zgdj2000/0004/image4/66-9.gif[/img]?。?8） 　[img=67,24]http://zszl.cepee.com/cepee_kjlw_pic/files/wx/zgdjgcxb/zgdj2000/0004/image4/66-10.gif[/img]　（19） 　[img=87,37]http://zszl.cepee.com/cepee_kjlw_pic/files/wx/zgdjgcxb/zgdj2000/0004/image4/66-11.gif[/img]?。?0） 　[img=53,24]http://zszl.cepee.com/cepee_kjlw_pic/files/wx/zgdjgcxb/zgdj2000/0004/image4/66-12.gif[/img]　（21） 　　[img=61,26]http://zszl.cepee.com/cepee_kjlw_pic/files/wx/zgdjgcxb/zgdj2000/0004/image4/66-13.gif[/img]（22） 　[img=74,25]http://zszl.cepee.com/cepee_kjlw_pic/files/wx/zgdjgcxb/zgdj2000/0004/image4/66-14.gif[/img]?。?3） 　　[img=65,24]http://zszl.cepee.com/cepee_kjlw_pic/files/wx/zgdjgcxb/zgdj2000/0004/image4/66-15.gif[/img]（24） 　　式（18）～（24）構(gòu)造短路電流計(jì)算用的REI等值網(wǎng)絡(luò)很容易實(shí)現(xiàn)。在擴(kuò)展了REI網(wǎng)絡(luò)后采用局部節(jié)點(diǎn)等值法，則不需修正保留節(jié)點(diǎn)的電流源II，對保留節(jié)點(diǎn)等值導(dǎo)納矩陣YII的修正，可根據(jù)其特點(diǎn)用簡便的方法進(jìn)行。 [b]4　保留節(jié)點(diǎn)等值導(dǎo)納矩陣的分析與修正 [/b]　　按式（13）修正保留節(jié)點(diǎn)等值導(dǎo)納矩陣，對多節(jié)點(diǎn)配電網(wǎng)絡(luò)將出現(xiàn)第2節(jié)中所說的兩種情況。分析修正前后的保留節(jié)點(diǎn)導(dǎo)納矩陣YII與YNII之間的異同可以發(fā)現(xiàn)：①兩矩陣階數(shù)相同；②用Yij表示YII中的互導(dǎo)納，YNij表示YNII中的互導(dǎo)納，則YNij＝Y(jié)ij（i≠j），即保留節(jié)點(diǎn)相互之間的阻抗參數(shù)不變；③自導(dǎo)納YNii＝Y(jié)ii（i≠m），表示與節(jié)點(diǎn)i連接的其它節(jié)點(diǎn)沒有一個(gè)被消去；④自導(dǎo)納YNmm≠Ymm（m≠i），表示與節(jié)點(diǎn)m連接的節(jié)點(diǎn)中至少有一個(gè)被消去。 　　分析表明YNII對于YII，其中只有與被消去的節(jié)點(diǎn)有關(guān)聯(lián)的節(jié)點(diǎn)m對應(yīng)的YNmm與Ymm不相等。而YII中屬于節(jié)點(diǎn)m性質(zhì)的節(jié)點(diǎn)是極其有限的，因?yàn)閷?shí)際的配電網(wǎng)絡(luò)?？稍谀骋荒妇€處將網(wǎng)絡(luò)分為被保留部分與被消去部分。故將YII修正為YNII，只需根據(jù)網(wǎng)絡(luò)結(jié)構(gòu)逐個(gè)按定義求出YNm1m1～YNmnmn即可。 [b]5　算例 [/b]　　山西省天脊煤化工集團(tuán)供電系統(tǒng)有兩路110 kV進(jìn)線；兩臺分別為30 MW與16 MW的自備發(fā)電機(jī)；十幾臺2 MW的異步電機(jī)；十個(gè)較大的變電站；及180多臺需計(jì)算短路電流的用電設(shè)備。經(jīng)過最大可能化簡，尚有170多個(gè)節(jié)點(diǎn)，200多條支路需計(jì)算短路電流，才能滿足繼電保護(hù)整定計(jì)算與合理選擇電氣設(shè)備等安全可靠供電管理的要求。而該單位用于供電管理的計(jì)算機(jī)受容量限制在Windows 95或98環(huán)境下運(yùn)行不能計(jì)算多于約100個(gè)節(jié)點(diǎn)的網(wǎng)絡(luò)。應(yīng)用本文提出的方法，我們?yōu)樵撓到y(tǒng)編制的短路電流計(jì)算軟件，在現(xiàn)有的計(jì)算機(jī)上運(yùn)行良好，滿足了該單位供電管理的要求。 　　5.1　算法步驟 　?。?）將圖2所示網(wǎng)絡(luò)在適當(dāng)?shù)哪妇€處分解為兩（或多個(gè)）部分，以一個(gè)約100節(jié)點(diǎn)的網(wǎng)作為內(nèi)部系統(tǒng)，其余作為外部系統(tǒng)。內(nèi)部系統(tǒng)的節(jié)點(diǎn)作為保留節(jié)點(diǎn)。編節(jié)點(diǎn)號時(shí)先編保留節(jié)點(diǎn)中無源節(jié)點(diǎn)，然后編保留節(jié)點(diǎn)中的電源節(jié)點(diǎn)，再編邊界節(jié)點(diǎn)，最后編準(zhǔn)備消去的節(jié)點(diǎn)，如圖2所示。圖中電抗參數(shù)為實(shí)際元件算出的標(biāo)么值，為說明問題，略去了低壓網(wǎng)應(yīng)計(jì)算的電阻參數(shù)。節(jié)點(diǎn)13，14為邊界節(jié)點(diǎn)。 　?。?）將準(zhǔn)備消去的外部系統(tǒng)中具有相關(guān)性質(zhì)（如電源、大型異步電機(jī)等）的有源節(jié)點(diǎn)歸為一組，如圖2中節(jié)點(diǎn)15～17、18～20。用以上方法求出每一組的YR與YK，構(gòu)造出REI網(wǎng)，如圖3所示，得 　 　圖2　算例配電網(wǎng)接線圖 　　 Fig.2　Distribution network sample 　 　圖3　外部系統(tǒng)有源節(jié)點(diǎn)構(gòu)造REI網(wǎng) 　　 Fig.3　Make up supply mains nodels of outside 　　 part system into REI network 　　SR1＝S15＋S16＋S17，經(jīng)阻抗X13～R1接到邊界節(jié)點(diǎn)13上，得 　　[img=294,35]http://zszl.cepee.com/cepee_kjlw_pic/files/wx/zgdjgcxb/zgdj2000/0004/image4/67-1.gif[/img][img=194,37]http://zszl.cepee.com/cepee_kjlw_pic/files/wx/zgdjgcxb/zgdj2000/0004/image4/67-2.gif[/img]　節(jié)點(diǎn)15～17成為無源節(jié)點(diǎn)，可歸入消去節(jié)點(diǎn)集E中，用同樣的方法可對有源節(jié)點(diǎn)18～20進(jìn)行修正。帶有電源的節(jié)點(diǎn)13，14稱為邊界節(jié)點(diǎn)，將其歸入保留節(jié)點(diǎn)集I，其等值網(wǎng)絡(luò)如圖4所示。 　　 圖4　保留節(jié)點(diǎn)等值網(wǎng)絡(luò) 　 　Fig.4　Reservation nodes equivalent network 　?。?）按圖4等值網(wǎng)絡(luò)輸入原始數(shù)據(jù)，計(jì)算出短路電流已足夠精確。此條件下形成的導(dǎo)納矩陣只修正了外部系統(tǒng)電源，而未修正外部系統(tǒng)導(dǎo)納，用YN1II表示。YN1II與YNII的區(qū)別與第4節(jié)中分析相同，只有邊界點(diǎn)的自導(dǎo)納不相等。如要求更精確的短路電流值，需將YN1II修正為YNII，即將邊界節(jié)點(diǎn)的自導(dǎo)納YN1mm修正為YNmm。在算例中邊界點(diǎn)13的導(dǎo)納為YN11313。修正方法是先按圖2求出Y1313，然后按節(jié)點(diǎn)13周圍的結(jié)構(gòu)求出[img=168,32]http://zszl.cepee.com/cepee_kjlw_pic/files/wx/zgdjgcxb/zgdj2000/0004/image4/67-3.gif[/img][img=115,34]http://zszl.cepee.com/cepee_kjlw_pic/files/wx/zgdjgcxb/zgdj2000/0004/image4/67-4.gif[/img]用此方法逐一求出Y[sup]N[/sup][sub]mm[/sub]，用Y[sup]N[/sup][sub]II[/sub]作為計(jì)算過程應(yīng)用的保留節(jié)點(diǎn)導(dǎo)納矩陣。 　　配電網(wǎng)絡(luò)每一支路在任意運(yùn)行方式下，任意處、任何類型的短路電流均可按上述方法計(jì)算。 　　5.2　精度分析 　　以算例中13點(diǎn)發(fā)生三相短路為例說明計(jì)算精度，該點(diǎn)故障時(shí)受異步電機(jī)供出短路電流影響最大。 　?。?）按常規(guī)短路電流計(jì)算方法計(jì)算 [img=297,58]http://zszl.cepee.com/cepee_kjlw_pic/files/wx/zgdjgcxb/zgdj2000/0004/image4/67-5.gif[/img] 　?。?）按本文提出的短路電流計(jì)算方法計(jì)算 [img=316,49]http://zszl.cepee.com/cepee_kjlw_pic/files/wx/zgdjgcxb/zgdj2000/0004/image4/67-6.gif[/img] 　?。?）將節(jié)點(diǎn)導(dǎo)納矩陣YN1II修正為YNII時(shí)Id＝2.63 　　5.3　程序框圖 [b]6　結(jié)論 [/b]　?。?） 多節(jié)點(diǎn)配電系統(tǒng)可以采用本文提出的方法在適當(dāng)?shù)哪妇€處分解為內(nèi)外部系統(tǒng)，依次將每一部分作為內(nèi)部系統(tǒng)完成短路電流計(jì)算。 　?。?） REI等值網(wǎng)的作用是將外部系統(tǒng)中的電源等值聯(lián)接到邊界節(jié)點(diǎn)上，邊界節(jié)點(diǎn)歸入保留節(jié)點(diǎn)。原有源節(jié)點(diǎn)變成無源節(jié)點(diǎn)，使短路電流計(jì)算精度明顯提高。 [b]參考文獻(xiàn)： [/b]　　［1］　諸駿偉.電力系統(tǒng)分析（上冊）［M］.北京：水利電力出版社，1995，11：119～124. 　?。?］　米麟書，江世芳.電網(wǎng)繼電保護(hù)的計(jì)算機(jī)輔助分析［M］.北京：水利電力出版社，1995，11：25～44. 　?。?］　Wu F F, Narasimhamurthi N. Necessary Condition for REI Reduction to be Exact［C］. IEEE PES Winter Meeting, 1979：A79 065～4. 　?。?］　陳亞民.電力系統(tǒng)計(jì)算程序及其實(shí)現(xiàn)［M］.北京：水利電力出版社，1995，11：11～146.