專利名稱：：基于擴(kuò)展不確定度的電力系統(tǒng)狀態(tài)估計(jì)結(jié)果評價(jià)方法
技術(shù)領(lǐng)域：
：本技術(shù)發(fā)明屬于電力系統(tǒng)分析計(jì)算領(lǐng)域。
背景技術(shù)：
：在測點(diǎn)真值未知情況下，如何對電力系統(tǒng)狀態(tài)估計(jì)結(jié)果給出合理的評價(jià)方法，是進(jìn)行狀態(tài)估計(jì)研究首先并且必須要解決的問題。以往研究中僅僅給出了測點(diǎn)真值已知情況下的評價(jià)指標(biāo)，主要應(yīng)用于仿真研究，通過人為設(shè)定真值和量測值，進(jìn)行狀態(tài)估計(jì)研究，用范數(shù)和無窮范數(shù)對估計(jì)值和真值之間距離進(jìn)行衡量，即又=2]1、-^1，5"2=maxIx,.—I其中x,，r,分別為狀態(tài)估計(jì)結(jié)果:c的第/個(gè)分量的估計(jì)值和真值。&,S2分別用來衡量狀態(tài)估計(jì)結(jié)果整體偏差和最大偏差。由此可以看出，以往的研究中，被測量真值已知情況下，人們主要通過測量誤差來評價(jià)狀態(tài)估計(jì)結(jié)果的優(yōu)劣。而在實(shí)際應(yīng)用中，被測量真值是無法準(zhǔn)確獲知的未知量，因而誤差僅是理論上的概念，實(shí)際上是無法確切獲知的，使用起來具有局限性。如何在被測量真值未知情況下對狀態(tài)估計(jì)結(jié)果進(jìn)行合理評價(jià)，是目前狀態(tài)估計(jì)研究中的盲區(qū)，但卻是進(jìn)行狀態(tài)估計(jì)首先且必須要解決的重要問題。由于認(rèn)識到測量殘差概念的不足，國際計(jì)量學(xué)界在上世紀(jì)九十年代引入測量不確定度概念，并建議用采用測量不確定度對測量結(jié)果進(jìn)行評價(jià)。測量不確定度是定量說明測量結(jié)果質(zhì)量的一個(gè)參數(shù)，反映了可能的誤差分布范圍。不確定度只與測量條件有關(guān)而與測量結(jié)果無關(guān)，可以近似地理解為一定置信概率下的誤差限值。利用測量不確定度代替誤差表示測量結(jié)果，易于了解、便于評定，具有合理性和實(shí)用性。測量不確定度分為標(biāo)準(zhǔn)不確定度和擴(kuò)展不確定度兩類。標(biāo)準(zhǔn)不確定度是用標(biāo)準(zhǔn)偏差表示測量結(jié)果的不確定度，分為用A類計(jì)算方法得到的標(biāo)準(zhǔn)不確定度、用B類計(jì)算方法得到的標(biāo)準(zhǔn)不確定度和根據(jù)其它一些量值求出的合成標(biāo)準(zhǔn)不確定度。擴(kuò)展不確定度由合成標(biāo)準(zhǔn)不確定度^乘上覆蓋因子A:而得到。A:值一般為2，有時(shí)為3，取決于被測量的重要性、效益和風(fēng)險(xiǎn)。被測量真值分布的大部分可望包含在擴(kuò)展不確定度表示的區(qū)間內(nèi)，因而本發(fā)明提出的狀態(tài)估計(jì)結(jié)果評價(jià)方法是基于擴(kuò)展不確定度進(jìn)行的。人們可根據(jù)測量數(shù)據(jù)、有關(guān)材料和儀器的技術(shù)指標(biāo)、檢定證書和手冊中提供的數(shù)據(jù)來定量評定擴(kuò)展不確定度。我國于1999年，經(jīng)國家質(zhì)量技術(shù)監(jiān)督局批準(zhǔn)，頒布實(shí)施由全國法制計(jì)量技術(shù)委員會提出的《測量不確定度評定與表示》(JJF1059—1999)，依據(jù)此規(guī)范，不確定度的一般評定步驟為第l步建模；第2步分析影響測量的所有不確定度來源；第3步對符合A類評定的不確定度分量用A類評定方法，計(jì)算出A類標(biāo)準(zhǔn)不確定度；第4步將剩余的不確定度分量，分別用B類評定方法，計(jì)算出B類標(biāo)準(zhǔn)不確定度；第5歩求出合成標(biāo)準(zhǔn)不確定度；第6歩將合成標(biāo)準(zhǔn)不確定度^乘以置信因子k，求出擴(kuò)展不確定度f7-yb/c。本發(fā)明提出了一種電力系統(tǒng)狀態(tài)估計(jì)結(jié)果評價(jià)方法?；跀U(kuò)展不確定度，給出了正常測點(diǎn)和異常測點(diǎn)的定義，進(jìn)而構(gòu)建了測點(diǎn)正常率指標(biāo)和;72/3、/71/3指標(biāo)，利用這些指標(biāo)進(jìn)行狀態(tài)估計(jì)結(jié)果評價(jià)。同時(shí)，給出了衡量狀態(tài)估計(jì)結(jié)果合理性的具體方法。
發(fā)明內(nèi)容本發(fā)明的目的在于提供科學(xué)合理的電力系統(tǒng)狀態(tài)估計(jì)結(jié)果評價(jià)方法。本發(fā)明的特征在于所述方法是在調(diào)度中心的一個(gè)電力系統(tǒng)狀態(tài)估計(jì)結(jié)果評價(jià)計(jì)算機(jī)上依次按以下步驟實(shí)現(xiàn)的歩驟(l)，所述電力系統(tǒng)狀態(tài)估計(jì)結(jié)果評價(jià)計(jì)算機(jī)從能量管理系統(tǒng)中獲得電力系統(tǒng)的網(wǎng)絡(luò)參數(shù)，包括輸電支路的各支路號、首節(jié)點(diǎn)編號、末節(jié)點(diǎn)編號、串聯(lián)電阻、串聯(lián)電抗、并聯(lián)電導(dǎo)、并聯(lián)電納、以及變壓器的變比和阻抗；步驟(2)，從數(shù)據(jù)采集和監(jiān)視系統(tǒng)SCADA的網(wǎng)絡(luò)數(shù)據(jù)服務(wù)器獲得遙測量測值和遙信量測值，其中遙測量測值包括電壓幅值、發(fā)電機(jī)有功功率、發(fā)電機(jī)無功功率、負(fù)荷有功功率、負(fù)荷無功功率、支路首端有功功率、支路首端無功功率、支路末端有功功率、支路末端無功功率、支路電流幅值；遙信量測值包括開關(guān)、刀閘開合狀態(tài)以及變壓器分接頭位置，開關(guān)或刀閘閉合的量測值為i，斷開的量測值為o;步驟p)，構(gòu)造計(jì)算模型歩驟(3.1)，遍歷所述各開關(guān)、支路、變壓器、發(fā)電機(jī)、電容器以及負(fù)荷端口各節(jié)點(diǎn)，根據(jù)得到的各開關(guān)量把電網(wǎng)中通過閉合的開關(guān)或刀閘連接在一起的節(jié)點(diǎn)稱為計(jì)算用節(jié)點(diǎn)；步驟(3.2)，把歩驟(2)得到的各計(jì)算用節(jié)點(diǎn)的電壓幅值、發(fā)電機(jī)功率以及負(fù)荷功率的量測值按以下方法分配到對應(yīng)的所述計(jì)算用節(jié)點(diǎn)上若所述計(jì)算用節(jié)點(diǎn)上有多個(gè)電壓幅值量測值，則取平均值作為一個(gè)電壓幅只有所述計(jì)算用節(jié)點(diǎn)上所有的發(fā)電機(jī)和負(fù)荷都有功率量測時(shí)，才能把發(fā)電機(jī)功率減去負(fù)荷功率作為該計(jì)算用節(jié)點(diǎn)上的注入功率量測值；把各計(jì)算用節(jié)點(diǎn)上的電壓幅值、各計(jì)算用節(jié)點(diǎn)上的注入功率以及各支路上的潮流這些量測量定義為測點(diǎn)；歩驟(4)，從數(shù)據(jù)庫中獲取測量不確定度信息，其中包括測點(diǎn)序號、量測量、標(biāo)準(zhǔn)差以及擴(kuò)展不確定度c/;歩驟(5)，基于步驟(1)和步驟(2)得到的數(shù)據(jù)，運(yùn)行調(diào)度中心的n套狀態(tài)估計(jì)軟件，得到n組實(shí)時(shí)狀態(tài)估計(jì)結(jié)果V,x2,…,x"，所述實(shí)時(shí)狀態(tài)估計(jì)的狀態(tài)變量為計(jì)算用節(jié)點(diǎn)電壓的幅值v和相角e，所述n的值是任意預(yù)設(shè)的；歩驟(6)，根據(jù)步驟(5)得到的狀態(tài)估計(jì)結(jié)果x^，…，x"，分別計(jì)算出不同狀態(tài)結(jié)果下各測點(diǎn)的估計(jì)值，在狀態(tài)估計(jì)結(jié)果/下，各測點(diǎn)的估計(jì)值為/v,()、~()、&y()、"1,2…,"，計(jì)算用節(jié)點(diǎn)/的電壓幅值為，/v,(Z)=V;計(jì)算用節(jié)點(diǎn)/的有功注入功率為，盧計(jì)算用節(jié)點(diǎn)z'的無功注入功率為，=I/(Gysin(《-《)-^cos(《-;其中，/v,(/)、；^,()和/^(x"統(tǒng)一用^(x"表示，&()中&代表計(jì)算用節(jié)點(diǎn)/上的量測量g，g是計(jì)算用節(jié)點(diǎn)z的電壓幅值K或注入有功功率s、注入無功功率G;從計(jì)算用節(jié)點(diǎn)/到計(jì)算用節(jié)點(diǎn)y的支路/-7'上/側(cè)有功潮流為，=(v,"、-v,v《c一,-o-"v盡"<);從計(jì)算用節(jié)點(diǎn)/到計(jì)算用節(jié)點(diǎn)y的支路,'-y上/側(cè)無功潮流為，=v,VGysin(《-O",VAco);其中，(/)和^(/)統(tǒng)一用、(/)表示，1(/)中^代表從計(jì)算用節(jié)點(diǎn)/到計(jì)算用節(jié)點(diǎn)7'的支路/-7'上/側(cè)的量測量^，^是/側(cè)有功潮流《或無功潮流0/;從計(jì)算用節(jié)點(diǎn)/到計(jì)算用節(jié)點(diǎn)y的支路/-/上y側(cè)有功潮流為，從計(jì)算用節(jié)點(diǎn)/到計(jì)算用節(jié)點(diǎn)y的支路卜7上y側(cè)無功潮流為，=一(v/)2(A+sin(^—0+^cos(W;其中，、,(/)和/^(x"統(tǒng)一用/^(/)表示，K)中&代表從計(jì)算用節(jié)點(diǎn)/到計(jì)算用節(jié)點(diǎn)y的支路/-y'上y側(cè)的量測量g，g是y側(cè)有功潮流s,或無功潮流上述量測方程中，《，《分別為狀態(tài)估計(jì)結(jié)果^下計(jì)算用節(jié)點(diǎn)z'和計(jì)算用節(jié)點(diǎn)y的電壓幅值的估計(jì)值，6f，《分別為狀態(tài)估計(jì)結(jié)果/下計(jì)算用節(jié)點(diǎn)/和計(jì)算用節(jié)點(diǎn)y電壓相角的估計(jì)值，q為計(jì)算用節(jié)點(diǎn)/和計(jì)算用節(jié)點(diǎn)"司的電導(dǎo)，&為計(jì)算用節(jié)點(diǎn)/和計(jì)算用節(jié)點(diǎn)y間的電納，^為支路w的"形等值電路中計(jì)算用節(jié)點(diǎn)/和計(jì)算用節(jié)點(diǎn)乂的接地電納^;步驟(7)，對步驟(5)得到的所述狀態(tài)估計(jì)結(jié)果x',x2，…,？的合理性進(jìn)行評價(jià)歩驟(7.1)，對于測點(diǎn)g，g為上述測點(diǎn)g,，或&,的任意一種，在狀態(tài)估計(jì)結(jié)果/下，若有l(wèi)^(v)-zg|st/g，其中，4為該測點(diǎn)g的量測值，~(/)為該測點(diǎn)g的估計(jì)值，t/g為測點(diǎn)g的擴(kuò)展不確定度，則在所述狀態(tài)估計(jì)結(jié)果/下，測點(diǎn)g為正常測點(diǎn)，否則，所述測點(diǎn)g為異常測點(diǎn)；定義狀態(tài)估計(jì)結(jié)果/對應(yīng)的測點(diǎn)正常率為—/;*100%，其中，r為系統(tǒng)中總的測點(diǎn)數(shù)目，^為狀態(tài)估計(jì)結(jié)果？對應(yīng)的正常測點(diǎn)數(shù)目步驟(7.2)，在狀態(tài)估計(jì)結(jié)果/下，若有/個(gè)測點(diǎn)滿足|\(/)-&/s^，則定義狀態(tài)估計(jì)結(jié)果？的;/2/3指標(biāo)為義4=/〃*，/0,步驟(7.3)，在狀態(tài)估計(jì)結(jié)果/下，若有,個(gè)測點(diǎn)滿足|~(/)-4則定義狀態(tài)估計(jì)結(jié)果/的//1/3指標(biāo)為^，^步驟(7.4)，根據(jù)歩驟(7.1)，計(jì)算并比較所述各狀態(tài)估計(jì)結(jié)果？,x2,…,:c"對應(yīng)的測點(diǎn)正常率;71,72,…,^，測點(diǎn)正常率更大的狀態(tài)估計(jì)結(jié)果更合理；歩驟(7.5)，對于步驟(7.4)所述的n個(gè)測點(diǎn)正常率;/,;72,…,;7",若有/個(gè)狀態(tài)估計(jì)結(jié)果對應(yīng)的測點(diǎn)正常率相同，^w，則對于這/個(gè)狀態(tài)估計(jì)結(jié)果，繼續(xù)比較它們的&。指標(biāo)，/72/3指標(biāo)更大的狀態(tài)估計(jì)結(jié)果更合理；步驟(7.6)，對于步驟(7.5)所述的/個(gè)狀態(tài)估計(jì)結(jié)果中，若仍有m個(gè)狀態(tài)估計(jì)結(jié)果的72/3指標(biāo)相同，則繼續(xù)比較這附個(gè)狀態(tài)估計(jì)結(jié)果的;71/3指標(biāo)，仏/3指標(biāo)更大的狀態(tài)估計(jì)結(jié)果更合理；步驟(7.7)，對于步驟(7.6)所述的m個(gè)狀態(tài)估計(jì)結(jié)果中，若仍有w個(gè)狀態(tài)估計(jì)結(jié)果的/71/3指標(biāo)相同，w《m，則認(rèn)為該w個(gè)狀態(tài)估計(jì)結(jié)果的合理性相同。本發(fā)明通過構(gòu)造測點(diǎn)正常率/7指標(biāo)、"2/3指標(biāo)、％3指標(biāo)，給出了科學(xué)的電力系統(tǒng)狀態(tài)估計(jì)結(jié)果評價(jià)方法，可以對各實(shí)時(shí)狀態(tài)估計(jì)結(jié)果進(jìn)行合理性評價(jià)，并輸出合理性最高的狀態(tài)估計(jì)結(jié)果。該評價(jià)方法計(jì)算簡單，具有實(shí)用性和可推廣性的優(yōu)點(diǎn)。圖1是基于擴(kuò)展不確定度的電力系統(tǒng)狀態(tài)估計(jì)結(jié)果評價(jià)流程圖。圖2是支路的;r形等值電路圖。圖3是變壓器的;r形等值電路圖。圖4是算例中的四節(jié)點(diǎn)系統(tǒng)量測圖0，有功無功量測對，A，電壓量測。圖5是算例中的四節(jié)點(diǎn)系統(tǒng)計(jì)算模型圖。具體實(shí)施例方式本發(fā)明按照以下步驟進(jìn)行歩驟(l)，從能量管理系統(tǒng)中獲得電力系統(tǒng)的網(wǎng)絡(luò)參數(shù)，包括輸電支路的各支路號、首節(jié)點(diǎn)編號、末節(jié)點(diǎn)編號、串聯(lián)電阻、串聯(lián)電抗、并聯(lián)電導(dǎo)、并聯(lián)電納、以及變壓器的變比和阻抗。步驟(2)，從數(shù)據(jù)采集和監(jiān)視系統(tǒng)SCADA的網(wǎng)絡(luò)數(shù)據(jù)服務(wù)器獲得遙測量測值和遙信量測值，其中遙測量測值包括電壓幅值、發(fā)電機(jī)有功功率、發(fā)電機(jī)無功功率、負(fù)荷有功功率、負(fù)荷無功功率、支路首端有功功率、支路首端無功功率、支路末端有功功率、支路末端無功功率、支路電流幅值；遙信量測值包括開關(guān)、刀閘開合狀態(tài)以及變壓器分接頭位置，開關(guān)或刀閘閉合的量測值為1，斷開的量測值為0;所述SCADA系統(tǒng)采用輪詢方式采集數(shù)據(jù)，輪詢周期為510s。離線狀態(tài)下可以通過文件給出上述量測量，在線應(yīng)用時(shí)可以通過SCADA系統(tǒng)提供的API或使用TCP/IP通信協(xié)議從SCADA系統(tǒng)中獲得。步驟(3)，構(gòu)造計(jì)算模型歩驟(3.1)，把開關(guān)、支路、變壓器、發(fā)電機(jī)、電容器以及負(fù)荷端口定義為節(jié)點(diǎn)，遍歷所述這些節(jié)點(diǎn)，根據(jù)得到的各開關(guān)量把電網(wǎng)中通過閉合的開關(guān)或刀閘連接在一起的節(jié)點(diǎn)稱為計(jì)算用節(jié)點(diǎn)，下面所述節(jié)點(diǎn)均指計(jì)算用節(jié)點(diǎn)；步驟(3.2)，把步驟(2)得到的各計(jì)算用節(jié)點(diǎn)的電壓幅值、發(fā)電機(jī)功率以及負(fù)荷功率的量測值按以下方法分配到對應(yīng)的所述計(jì)算用節(jié)點(diǎn)上若所述計(jì)算用節(jié)點(diǎn)上有多個(gè)電壓幅值量測值，則取平均值作為一個(gè)電壓幅值量測值；只有所述計(jì)算用節(jié)點(diǎn)上所有的發(fā)電機(jī)和負(fù)荷都有功率量測時(shí)，才能把發(fā)電機(jī)功率減去負(fù)荷功率作為該計(jì)算用節(jié)點(diǎn)上的注入功率量測值；把各計(jì)算用節(jié)點(diǎn)上的電壓幅值、各計(jì)算用節(jié)點(diǎn)上的注入功率以及各支路上的潮流這些量測量定義為測點(diǎn)。歩驟(4)，從數(shù)據(jù)庫中獲取測量不確定度信息，其中包括測點(diǎn)序號、量測量、標(biāo)準(zhǔn)差以及擴(kuò)展不確定度t/。歩驟(5),基于歩驟(1)和步驟(2)得到的數(shù)據(jù)，運(yùn)行調(diào)度中心的n套狀態(tài)估計(jì)軟件，得到n組實(shí)時(shí)狀態(tài)估計(jì)結(jié)果;c'，;c2，…,x"，所述n的值是任意預(yù)設(shè)的。對于輸電網(wǎng)狀態(tài)估計(jì)而言，實(shí)時(shí)狀態(tài)估計(jì)的狀態(tài)變量為節(jié)點(diǎn)電壓的幅值v和相角步驟(6)，根據(jù)步驟(5)得到的狀態(tài)估計(jì)結(jié)果x',;c2，…,x"，分別計(jì)算出不同狀態(tài)結(jié)果下各測點(diǎn)的估計(jì)值，在狀態(tài)估計(jì)結(jié)果Z下，各測點(diǎn)的估計(jì)值為/v,()、Z^,(x"、1(/)、"1，2…，"，節(jié)點(diǎn)/的電壓幅值為，/v,(X;節(jié)點(diǎn)/的有功注入功率為，/^,(/)=2>>/(GyCOS(6^-《)+^sin(^-O);盧節(jié)點(diǎn)/的無功注入功率為，/i!e,(/)=^/(GySin(《-^cos(《-e/));其中，/v,(Z)、~,(^)和^(^)統(tǒng)一用^(^)表示，/^(/)中g(shù),代表節(jié)點(diǎn)z'上的量測量g，g是節(jié)點(diǎn)/的電壓幅值K或有功注入功率《、無功注入功率g;支路"_/上/側(cè)有功潮流為，支路"y上/側(cè)無功潮流為，1(/)=-—Gysin(《—e/)+v,，yC0-0;其中，(V)和/^(x"統(tǒng)一用/^(x"表示，/^(/)中&代表支路卜y上f側(cè)的量測量g，g是z'側(cè)有功潮流A或無功潮流仏；支路"y上7側(cè)有功潮流為，=(v/)2《-v>Xco《—o+A/A-O;支路"/上_/側(cè)無功潮流為，、(？)=-(v/)2(5yU+v,、Xsi《—《)+vyCO《-e/);其中，、,(？)和/^,(/)統(tǒng)一用\,(/)表示，^,(/)中g(shù),代表支路/-y'上y側(cè)的量測量g，g是y側(cè)有功潮流&或無功潮流仏；上述量測方程中，《，《分別為狀態(tài)估計(jì)結(jié)果/下節(jié)點(diǎn)/和節(jié)點(diǎn)/的電壓幅值的估計(jì)值，《，《分別為狀態(tài)估計(jì)結(jié)果z下節(jié)點(diǎn)/和節(jié)點(diǎn)y電壓相角的估計(jì)值，q為節(jié)點(diǎn)/和節(jié)點(diǎn)y間的電導(dǎo)，A為節(jié)點(diǎn)/和節(jié)點(diǎn)/間的電納，x為支路的;r形等值電路中節(jié)點(diǎn)/和節(jié)點(diǎn)y的接地電納^。支路zV的;r形等值電路如圖2所示。歩驟(7),對步驟(5)得到的狀態(tài)估計(jì)結(jié)果x、x2，…,x"的合理性進(jìn)行評價(jià)對于測點(diǎn)g，g為上述測點(diǎn)g,，^或^的任意一種，真值落在(Z「C/g,^+C/g)中概率為/，良口式中，f為系統(tǒng)真實(shí)狀態(tài)，4為該測點(diǎn)g的量測值，K為與置信概率p對應(yīng)的測點(diǎn)g的擴(kuò)展不確定度，可通過上面介紹的擴(kuò)展不確定度的評定方法獲取。對于測點(diǎn)g，在狀態(tài)估計(jì)結(jié)果/下，若有則在所述狀態(tài)估計(jì)結(jié)果/下，測點(diǎn)g為正常測點(diǎn)，否則，所述測點(diǎn)g為異常測占八"o設(shè)系統(tǒng)中總的測點(diǎn)數(shù)目為r,狀態(tài)估計(jì)結(jié)果z對應(yīng)的正常測點(diǎn)數(shù)目為y，則狀態(tài)估計(jì)結(jié)果/對應(yīng)的測點(diǎn)正常率為依次類推，可計(jì)算出狀態(tài)估計(jì)結(jié)果1^2,...,對應(yīng)的測點(diǎn)正常率分別為;/,/72,...，/7"。根據(jù)測點(diǎn)正常率對狀態(tài)估計(jì)結(jié)果做出評價(jià)，測點(diǎn)正常率越高，狀態(tài)估計(jì)結(jié)果越合理，即越接近系統(tǒng)的真實(shí)狀態(tài)。由于正常測點(diǎn)的定義與擴(kuò)展不確定^相關(guān)，而C/g又與一定的置信概率p相關(guān)，因而測點(diǎn)正常率是相應(yīng)于某一固定的置信概率/提出的。若狀態(tài)估計(jì)結(jié)果/，x2，…,x"中，有/個(gè)狀態(tài)估計(jì)結(jié)果的正常率相同，可認(rèn)為在置信概率p下衡量時(shí)，這/個(gè)狀態(tài)估計(jì)結(jié)果合理性相同。而在另一置信概率下衡量時(shí)，這/個(gè)狀態(tài)估計(jì)結(jié)果未必具有同等合理性。此時(shí)，需要引入/72/3指標(biāo)，仏/3指標(biāo)進(jìn)一步判別。對于上述的狀態(tài)估計(jì)結(jié)果/，若在總的r個(gè)測點(diǎn)中，有/個(gè)測點(diǎn)(/《^)滿3g則定義72,3指標(biāo)為《/3=//"100%同理，若在總的r個(gè)測點(diǎn)中，有f個(gè)測點(diǎn)(f、/)滿足W?！?丄t/3g則定義/71/3指標(biāo)為<=,/"100%對于測點(diǎn)正常率相同的/個(gè)狀態(tài)估計(jì)結(jié)果，通過比較72/3指標(biāo)進(jìn)一步判別狀態(tài)估計(jì)結(jié)果的合理性。72/3指標(biāo)更高的狀態(tài)估計(jì)結(jié)果更合理。如果這/個(gè)狀態(tài)估計(jì)結(jié)果中，仍有附個(gè)狀態(tài)估計(jì)結(jié)果的;/2/3指標(biāo)相同，附W，則繼續(xù)比較這m個(gè)狀態(tài)估計(jì)結(jié)果的/71/3指標(biāo)。71/3指標(biāo)更高的狀態(tài)估計(jì)結(jié)果更合理。如果這m個(gè)狀態(tài)估計(jì)結(jié)果中，仍有w個(gè)狀態(tài)估計(jì)結(jié)果的"。指標(biāo)相同，WSm，則認(rèn)為這w個(gè)狀態(tài)估計(jì)結(jié)果合理性相同。所述方法中當(dāng)支路中含變壓器時(shí)，支路兩端節(jié)點(diǎn)側(cè)的潮流的量測方程如下變壓器支路W上節(jié)點(diǎn)/側(cè)有功潮流為，變壓器支路&上節(jié)點(diǎn)z'側(cè)無功潮流為，=-如""、+vco《-o;變壓器支路zV上節(jié)點(diǎn)y側(cè)有功潮流為，變壓器支路上節(jié)點(diǎn)y'側(cè)無功潮流為，=-)、十會v,V-《)；上述量測方程中，K為變壓器的非標(biāo)準(zhǔn)變比，^為變壓器標(biāo)準(zhǔn)側(cè)的電納。變壓器的v形等值電路如圖3所示。以圖4所示的4節(jié)點(diǎn)系統(tǒng)為例，具體步驟說明如下步驟1:初始化從現(xiàn)有的能量管理系統(tǒng)中獲得電力系統(tǒng)的網(wǎng)絡(luò)參數(shù)。包括支路的串聯(lián)電阻、串聯(lián)電抗、并聯(lián)電導(dǎo)和并聯(lián)電納；變壓器的變比和阻抗。圖4所示的4節(jié)點(diǎn)系統(tǒng)支路參數(shù)如表1所示表l四節(jié)點(diǎn)系統(tǒng)支路參數(shù)<table>tableseeoriginaldocumentpage21</column></row><table>步驟2:獲得SCADA量測值某一個(gè)時(shí)刻4節(jié)點(diǎn)系統(tǒng)的量測數(shù)據(jù)如表2所示，其中遙信量的量測值為1表示開關(guān)閉合，0表示開關(guān)斷開，遙測量的量測值使用有名值有功功率的單位是MW,無功功率的單位是MVAR，電壓幅值的單位是kV。表2四節(jié)點(diǎn)系統(tǒng)量測數(shù)據(jù)<table>tableseeoriginaldocumentpage22</column></row><table>步驟4:從數(shù)據(jù)庫中獲取測點(diǎn)不確定度信息，如表4所示,表4四節(jié)點(diǎn)系統(tǒng)測點(diǎn)標(biāo)準(zhǔn)差以及擴(kuò)展不確定度測點(diǎn)序號量測量標(biāo)準(zhǔn)差擴(kuò)展不確定度U1Pi1.33.92Qi1.44.23P121.434,294Q121.424.265P131.293.876Q131.183.547P311.293.878Q311.243.729P321.293.8710Q321.444.3211P31.54.512Q31.434.2913P342.497.4714Q341.634.8915Vi1.113.3316v31.113.33步驟5:運(yùn)行狀態(tài)估計(jì)軟件，獲取狀態(tài)估計(jì)結(jié)果。運(yùn)行狀態(tài)估計(jì)軟件，獲取狀態(tài)估計(jì)結(jié)果x1、x2、x3、x4、jc5，如表中電壓的單位為kV，相角的單位為弧度。所示，表表5四節(jié)點(diǎn)系統(tǒng)狀態(tài)估計(jì)結(jié)果狀態(tài)估計(jì)結(jié)果？狀態(tài)估計(jì)結(jié)果x2狀態(tài)估計(jì)狀態(tài)估計(jì)結(jié)果x4狀態(tài)估計(jì)結(jié)果x5111.4636110.74976111.4967111.4749112.5207v2110.3402109.60363110.3731110.3541111.3933v3111.7018110.99001111.7291111.7117112.7166v4112.3304111.83313112.3254112.3286113.1516《00000《-0.40344-0.427315-0.40323-0.40178-0.39083《0.8134990.82324330.8003780.8049150.791826&7.9572718.17766387.8569937.933027.787569步驟6:計(jì)算所有量測量的估計(jì)值各量測量的量測方程如下節(jié)點(diǎn)/注入功率Pi，Qi:J。ye,支路/:/上節(jié)點(diǎn)/側(cè)潮流《=v,2g-v,"gcos《-v,v,sin《支路/:/上節(jié)點(diǎn)y側(cè)潮流=—、(6+Jc)+v,Vjgsin6>y+v,6cosPy變壓器支路&上節(jié)點(diǎn)/側(cè)潮流變壓器支路上節(jié)點(diǎn)/側(cè)潮流節(jié)點(diǎn)/電壓幅值上述量測方程中，v,，v,分別為節(jié)點(diǎn)i和節(jié)點(diǎn)j的電壓幅值，《為節(jié)點(diǎn)i和節(jié)點(diǎn)j電壓相角差，Gy，^為導(dǎo)納矩陣元素，g，b，尺為支路;r形模型中的參數(shù)，K為變壓器的非標(biāo)準(zhǔn)變比，6,為變壓器標(biāo)準(zhǔn)側(cè)的電納。將步驟5計(jì)算得到的狀態(tài)變量代入到量測方程中，計(jì)算出各量測量的估計(jì)值"/^)，求得的結(jié)果如表6所示。表6四節(jié)點(diǎn)系統(tǒng)量測量估計(jì)結(jié)果量測量量測值狀態(tài)估計(jì)結(jié)果x1狀態(tài)估計(jì)結(jié)果x2狀態(tài)估計(jì)結(jié)果x3狀態(tài)估計(jì)結(jié)果:c4狀態(tài)估計(jì)結(jié)果x518.7120.6322.1920.9520.6820.6734.0534.7535.2134.7234.6235.66p1241.7940.2441.7940.2540.1039.96Q1237.9637.5837.9637.6037,5138.22P13-19.1-19,62-19,60-19.30-19.42-19.29Q13-2.24-2.83-2.75-2.88-2.89-2.55P3117.5519.7119.7019.4019.5119.38Q3-10.74-10.84-10.74-10.81-10.79-11.39132.57130.73132.57129.45129,73129.30Q3242.6141.6142.2241.5541.5840.37P3-49.49-47.70-49.49-46.93-48.49-48.57Q3-41.93-41.22-41.93-41.14-41.11-42.01P34-193.22-198.15-201.76-195.78-197.73-197.26Q34-71.08-72.00-73.41-71.87-71.90-70.99V,112.15111.46110.75111.50111.47112.52v3110.99111.70110.99111.73111.71112.72步驟7:比較狀態(tài)估計(jì)結(jié)果的合理性首先，根據(jù)步驟3得到的測點(diǎn)量測值、步驟4得到的測點(diǎn)不確定度、步驟6得到的測點(diǎn)估計(jì)值，分別計(jì)算狀態(tài)估計(jì)結(jié)果？，x2,x3,x4,對應(yīng)的測點(diǎn)正常率7指標(biāo)，計(jì)算結(jié)果如表7所示。比較/7指標(biāo)，狀態(tài)估計(jì)結(jié)果^的合理性最差。x'乂,x4,x5的合理性比較轉(zhuǎn)到下一步。其次，改變擴(kuò)展不確定度取值為lc/時(shí)，計(jì)算出此時(shí)狀態(tài)估計(jì)結(jié)果5的;72/3指標(biāo)，計(jì)算結(jié)果如表7所示。比較;72/3指標(biāo)，x的合理性最好。、x4和？的合理性比較轉(zhuǎn)到下一步。最后，改變擴(kuò)展不確定度取值為^7時(shí)，計(jì)算出此時(shí)狀態(tài)估計(jì)結(jié)果xUW4乂的/^指標(biāo)，計(jì)算結(jié)果如表7所示。比較//1/3指標(biāo)，可以發(fā)現(xiàn)x3、x4和？的合理性順序?yàn)閄4>>x5。故最終所得合理性排序?yàn)?；c1>/>>x5>x2。表7四節(jié)點(diǎn)系統(tǒng)狀態(tài)估計(jì)結(jié)果的評價(jià)指標(biāo)^一狀態(tài)估計(jì)狀態(tài)估計(jì)狀態(tài)估計(jì)狀態(tài)估計(jì)狀態(tài)估計(jì)指標(biāo)結(jié)果x1結(jié)果x2結(jié)果x3結(jié)果x4結(jié)果x5/7100.0%93.8%100.0%100.0%100.0%72/3100.0%87.5%93.8%93.8%93.8%仏/362.5%68.8%62.5%68.8°/050.0%權(quán)利要求1.基于擴(kuò)展不確定度的電力系統(tǒng)狀態(tài)估計(jì)結(jié)果評價(jià)方法，其特征在于，所述方法是在調(diào)度中心的一個(gè)電力系統(tǒng)狀態(tài)估計(jì)結(jié)果評價(jià)計(jì)算機(jī)上依次按以下步驟實(shí)現(xiàn)的步驟(1)，所述電力系統(tǒng)狀態(tài)估計(jì)結(jié)果評價(jià)計(jì)算機(jī)從能量管理系統(tǒng)中獲得電力系統(tǒng)的網(wǎng)絡(luò)參數(shù)，包括輸電支路的各支路號、首節(jié)點(diǎn)編號、末節(jié)點(diǎn)編號、串聯(lián)電阻、串聯(lián)電抗、并聯(lián)電導(dǎo)、并聯(lián)電納、以及變壓器的變比和阻抗；步驟(2)，從數(shù)據(jù)采集和監(jiān)視系統(tǒng)SCADA的網(wǎng)絡(luò)數(shù)據(jù)服務(wù)器獲得遙測量測值和遙信量測值，其中遙測量測值包括電壓幅值、發(fā)電機(jī)有功功率、發(fā)電機(jī)無功功率、負(fù)荷有功功率、負(fù)荷無功功率、支路首端有功功率、支路首端無功功率、支路末端有功功率、支路末端無功功率、支路電流幅值；遙信量測值包括開關(guān)、刀閘開合狀態(tài)以及變壓器分接頭位置，開關(guān)或刀閘閉合的量測值為1，斷開的量測值為0；步驟(3)，構(gòu)造計(jì)算模型步驟(3.1)，遍歷所述各開關(guān)、支路、變壓器、發(fā)電機(jī)、電容器以及負(fù)荷端口各節(jié)點(diǎn)，根據(jù)得到的各開關(guān)量把電網(wǎng)中通過閉合的開關(guān)或刀閘連接在一起的節(jié)點(diǎn)稱為計(jì)算用節(jié)點(diǎn)；步驟(3.2)，把步驟(2)得到的各計(jì)算用節(jié)點(diǎn)的電壓幅值、發(fā)電機(jī)功率以及負(fù)荷功率的量測值按以下方法分配到對應(yīng)的所述計(jì)算用節(jié)點(diǎn)上若所述計(jì)算用節(jié)點(diǎn)上有多個(gè)電壓幅值量測值，則取平均值作為一個(gè)電壓幅值量測值；只有所述計(jì)算用節(jié)點(diǎn)上所有的發(fā)電機(jī)和負(fù)荷都有功率量測時(shí)，才能把發(fā)電機(jī)功率減去負(fù)荷功率作為該計(jì)算用節(jié)點(diǎn)上的注入功率量測值；把各計(jì)算用節(jié)點(diǎn)上的電壓幅值、各計(jì)算用節(jié)點(diǎn)上的注入功率以及各支路上的潮流這些量測量定義為測點(diǎn)；步驟(4)，從數(shù)據(jù)庫中獲取測量不確定度信息，其中包括測點(diǎn)序號、量測量、標(biāo)準(zhǔn)差以及擴(kuò)展不確定度U；步驟(5)，基于步驟(1)和步驟(2)得到的數(shù)據(jù)，運(yùn)行調(diào)度中心的n套狀態(tài)估計(jì)軟件，得到n組實(shí)時(shí)狀態(tài)估計(jì)結(jié)果x1，x2，…，xn，所述實(shí)時(shí)狀態(tài)估計(jì)的狀態(tài)變量為計(jì)算用節(jié)點(diǎn)電壓的幅值v和相角θ，所述n的值是任意預(yù)設(shè)的；步驟(6)，根據(jù)步驟(5)得到的狀態(tài)估計(jì)結(jié)果x1，x2，…，xn，分別計(jì)算出不同狀態(tài)結(jié)果下各測點(diǎn)的估計(jì)值，在狀態(tài)估計(jì)結(jié)果xk下，各測點(diǎn)的估計(jì)值為hVi(xk)、hPi(xk)、hQi(xk)、hPij(xk)、hQij(xk)、hPji(xk)、hQji(xk)，k＝1，2…，n，計(jì)算用節(jié)點(diǎn)i的電壓幅值為，hVi(xk)＝vik；計(jì)算用節(jié)點(diǎn)i的有功注入功率為，<mathsid="math0001"num="0001"><math><![CDATA[<mrow><msub><mi>h</mi><mi>Pi</mi></msub><mrow><mo>(</mo><msup><mi>x</mi><mi>k</mi></msup><mo>)</mo></mrow><mo>=</mo><munder><mi>&Sigma;</mi><mrow><mi>j</mi><mo>&Element;</mo><mi>i</mi></mrow></munder><msup><msub><mi>v</mi><mi>i</mi></msub><mi>k</mi></msup><msup><msub><mi>v</mi><mi>j</mi></msub><mi>k</mi></msup><mrow><mo>(</mo><msub><mi>G</mi><mi>ij</mi></msub><mi>cos</mi><mrow><mo>(</mo><msup><msub><mi>&theta;</mi><mi>i</mi></msub><mi>k</mi></msup><mi></mi><mo>-</mo><msup><msub><mi>&theta;</mi><mi>j</mi></msub><mi>k</mi></msup><mi></mi><mo>)</mo></mrow><mo>+</mo><msub><mi>B</mi><mi>ij</mi></msub><mi>sin</mi><mrow><mo>(</mo><msup><msub><mi>&theta;</mi><mi>i</mi></msub><mi>k</mi></msup><mi></mi><mo>-</mo><msup><msub><mi>&theta;</mi><mi>j</mi></msub><mi>k</mi></msup><mi></mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>;</mo></mrow>]]></math></maths>計(jì)算用節(jié)點(diǎn)i的無功注入功率為，<mathsid="math0002"num="0002"><math><![CDATA[<mrow><msub><mi>h</mi><mi>Qi</mi></msub><mrow><mo>(</mo><msup><mi>x</mi><mi>k</mi></msup><mo>)</mo></mrow><mo>=</mo><munder><mi>&Sigma;</mi><mrow><mi>j</mi><mo>&Element;</mo><mi>i</mi></mrow></munder><msup><msub><mi>v</mi><mi>i</mi></msub><mi>k</mi></msup><msup><msub><mi>v</mi><mi>j</mi></msub><mi>k</mi></msup><mrow><mo>(</mo><msub><mi>G</mi><mi>ij</mi></msub><mi>sin</mi><mrow><mo>(</mo><msup><msub><mi>&theta;</mi><mi>i</mi></msub><mi>k</mi></msup><mi></mi><mo>-</mo><msup><msub><mi>&theta;</mi><mi>j</mi></msub><mi>k</mi></msup><mi></mi><mo>)</mo></mrow><mo>-</mo><msub><mi>B</mi><mi>ij</mi></msub><mi>cos</mi><mrow><mo>(</mo><msup><msub><mi>&theta;</mi><mi>i</mi></msub><mi>k</mi></msup><mi></mi><mo>-</mo><msup><msub><mi>&theta;</mi><mi>j</mi></msub><mi>k</mi></msup><mi></mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>;</mo></mrow>]]></math></maths>其中，hVi(xk)、hPi(xk)和hQi(xk)統(tǒng)一用hgi(xk)表示，hgi(xk)中g(shù)i代表計(jì)算用節(jié)點(diǎn)i上的量測量g，g是計(jì)算用節(jié)點(diǎn)i的電壓幅值Vi或有功注入功率Pi、無功注入功率Qi；從計(jì)算用節(jié)點(diǎn)i到計(jì)算用節(jié)點(diǎn)j的支路i-j上i側(cè)有功潮流為，hPij(xk)＝(vik)2Gij-vikvjkGijcos(θik-θjk)-vikvjkBijsin(θik-θjk)；從計(jì)算用節(jié)點(diǎn)i到計(jì)算用節(jié)點(diǎn)j的支路i-j上i側(cè)無功潮流為，hQij(xk)＝-(vik)2(Bij+yc)-vikvjkGijsin(θik-θjk)+vikvjkBijcos(θik-θjk)；其中，hPij(xk)和hQij(xk)統(tǒng)一用hgij(xk)表示，hgij(xk)中g(shù)ij代表從計(jì)算用節(jié)點(diǎn)i到計(jì)算用節(jié)點(diǎn)j的支路i-j上i側(cè)的量測量g，g是i側(cè)有功潮流Pij或無功潮流Qij；從計(jì)算用節(jié)點(diǎn)i到計(jì)算用節(jié)點(diǎn)j的支路i-j上j側(cè)有功潮流為，hPji(xk)＝(vjk)2Gij-vikvjkGijcos(θik-θjk)+vikvjkBijsin(θik-θjk)；從計(jì)算用節(jié)點(diǎn)i到計(jì)算用節(jié)點(diǎn)j的支路i-j上j側(cè)無功潮流為，hQji(xk)＝-(vjk)2(Bij+yc)+vikvjkGijsin(θik-θjk)+vikvjkBijcos(θik-θjk)；其中，hPji(xk)和hQji(xk)統(tǒng)一用hgji(xk)表示，hgji(xk)中g(shù)ji代表從計(jì)算用節(jié)點(diǎn)i到計(jì)算用節(jié)點(diǎn)j的支路i-j上j側(cè)的量測量g，g是j側(cè)有功潮流Pji或無功潮流Qji；上述量測方程中，vik，vjk分別為狀態(tài)估計(jì)結(jié)果xk下計(jì)算用節(jié)點(diǎn)i和計(jì)算用節(jié)點(diǎn)j的電壓幅值的估計(jì)值，θik，θjk分別為狀態(tài)估計(jì)結(jié)果xk下計(jì)算用節(jié)點(diǎn)i和計(jì)算用節(jié)點(diǎn)j電壓相角的估計(jì)值，Gij為計(jì)算用節(jié)點(diǎn)i和計(jì)算用節(jié)點(diǎn)j間的電導(dǎo)，Bij為計(jì)算用節(jié)點(diǎn)i和計(jì)算用節(jié)點(diǎn)j間的電納，yc為支路i-j的π形等值電路中計(jì)算用節(jié)點(diǎn)i和計(jì)算用節(jié)點(diǎn)j的接地電納yc；步驟(7)，對步驟(5)得到的所述狀態(tài)估計(jì)結(jié)果x1，x2，…，xn的合理性進(jìn)行評價(jià)步驟(7.1)，對于測點(diǎn)g，g為上述測點(diǎn)gi，gij或gji的任意一種，在狀態(tài)估計(jì)結(jié)果xk下，若有|hg(xk)-Zg|≤Ug，其中，Zg為該測點(diǎn)g的量測值，hg(xk)為該測點(diǎn)g的估計(jì)值，Ug為測點(diǎn)g的擴(kuò)展不確定度，則在所述狀態(tài)估計(jì)結(jié)果xk下，測點(diǎn)g為正常測點(diǎn)，否則，所述測點(diǎn)g為異常測點(diǎn)；定義狀態(tài)估計(jì)結(jié)果xk對應(yīng)的測點(diǎn)正常率為ηk＝qk/r*100％，其中，r為系統(tǒng)中總的測點(diǎn)數(shù)目，qk為狀態(tài)估計(jì)結(jié)果xk對應(yīng)的正常測點(diǎn)數(shù)目；步驟(7.2)，在狀態(tài)估計(jì)結(jié)果xk下，若有sk個(gè)測點(diǎn)滿足<mathsid="math0003"num="0003"><math><![CDATA[<mrow><mo>|</mo><msub><mi>h</mi><mi>g</mi></msub><mrow><mo>(</mo><msup><mi>x</mi><mi>k</mi></msup><mo>)</mo></mrow><mo>-</mo><msub><mi>Z</mi><mi>g</mi></msub><mo>|</mo><mo>&le;</mo><mfrac><mn>2</mn><mn>3</mn></mfrac><msub><mi>U</mi><mi>g</mi></msub><mo>,</mo></mrow>]]></math>id="icf0003"file="A2009100894790005C1.tif"wi="33"he="8"top="177"left="144"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>sk≤qk，則定義狀態(tài)估計(jì)結(jié)果xk的η2/3指標(biāo)為<mathsid="math0004"num="0004"><math><![CDATA[<mrow><msubsup><mi>&eta;</mi><mrow><mn>2</mn><mo>/</mo><mn>3</mn></mrow><mi>k</mi></msubsup><mo>=</mo><msup><mi>s</mi><mi>k</mi></msup><mo>/</mo><mi>r</mi><mo>*</mo><mn>100</mn><mo>%</mo><mo>,</mo></mrow>]]></math></maths>步驟(7.3)，在狀態(tài)估計(jì)結(jié)果xk下，若有tk個(gè)測點(diǎn)滿足<mathsid="math0005"num="0005"><math><![CDATA[<mrow><mo>|</mo><msub><mi>h</mi><mi>g</mi></msub><mrow><mo>(</mo><msup><mi>x</mi><mi>k</mi></msup><mo>)</mo></mrow><mo>-</mo><msub><mi>Z</mi><mi>g</mi></msub><mo>|</mo><mo>&le;</mo><mfrac><mn>1</mn><mn>3</mn></mfrac><msub><mi>U</mi><mi>g</mi></msub><mo>,</mo></mrow>]]></math>id="icf0005"file="A2009100894790005C3.tif"wi="33"he="8"top="219"left="133"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/></maths>tk≤sk，則定義狀態(tài)估計(jì)結(jié)果xk的η1/3指標(biāo)為<mathsid="math0006"num="0006"><math><![CDATA[<mrow><msubsup><mi>&eta;</mi><mrow><mn>1</mn><mo>/</mo><mn>3</mn></mrow><mi>k</mi></msubsup><mo>=</mo><msup><mi>t</mi><mi>k</mi></msup><mo>/</mo><mi>r</mi><mo>*</mo><mn>100</mn><mo>%</mo><mo>,</mo></mrow>]]></math></maths>步驟(7.4)，根據(jù)步驟(7.1)，計(jì)算并比較所述各狀態(tài)估計(jì)結(jié)果x1，x2，…，xn對應(yīng)的測點(diǎn)正常率η1，η2，…，ηn，測點(diǎn)正常率更大的狀態(tài)估計(jì)結(jié)果更合理；步驟(7.5)，對于步驟(7.4)所述的n個(gè)測點(diǎn)正常率η1，η2，…，ηn，若有l(wèi)個(gè)狀態(tài)估計(jì)結(jié)果對應(yīng)的測點(diǎn)正常率相同，l≤n，則對于這l個(gè)狀態(tài)估計(jì)結(jié)果，繼續(xù)比較它們的η2/3指標(biāo)，η2/3指標(biāo)更大的狀態(tài)估計(jì)結(jié)果更合理；步驟(7.6)，對于步驟(7.5)所述的l個(gè)狀態(tài)估計(jì)結(jié)果中，若仍有m個(gè)狀態(tài)估計(jì)結(jié)果的η2/3指標(biāo)相同，m≤l，則繼續(xù)比較這m個(gè)狀態(tài)估計(jì)結(jié)果的η1/3指標(biāo)，η1/3指標(biāo)更大的狀態(tài)估計(jì)結(jié)果更合理；步驟(7.7)，對于步驟(7.6)所述的m個(gè)狀態(tài)估計(jì)結(jié)果中，若仍有w個(gè)狀態(tài)估計(jì)結(jié)果的η1/3指標(biāo)相同，w≤m，則認(rèn)為該w個(gè)狀態(tài)估計(jì)結(jié)果的合理性相同。2.根據(jù)權(quán)利要求1所述的基于擴(kuò)展不確定度的電力系統(tǒng)狀態(tài)估計(jì)結(jié)果評價(jià)方法，其特征在于，所述支路中含變壓器時(shí)，支路兩端的計(jì)算用節(jié)點(diǎn)側(cè)的潮流的量測方程如下從計(jì)算用節(jié)點(diǎn)/到計(jì)算用節(jié)點(diǎn)7的變壓器支路W上節(jié)點(diǎn)/側(cè)有功潮流為，從計(jì)算用節(jié)點(diǎn)z到計(jì)算用節(jié)點(diǎn)乂的變壓器支路上節(jié)點(diǎn)z側(cè)無功潮流為，1"4)=-如。、-《)；從計(jì)算用節(jié)點(diǎn)/到計(jì)算用節(jié)點(diǎn)7的變壓器支路i-y上節(jié)點(diǎn)/側(cè)有功潮流為，從計(jì)算用節(jié)點(diǎn)/到計(jì)算用節(jié)點(diǎn)j的變壓器支路w上節(jié)點(diǎn)y側(cè)無功潮流為，上述量測方程中，K為變壓器的非標(biāo)準(zhǔn)變比，Zv為變壓器標(biāo)準(zhǔn)側(cè)的電納。全文摘要基于擴(kuò)展不確定度的電力系統(tǒng)狀態(tài)估計(jì)結(jié)果評價(jià)方法屬于電力系統(tǒng)狀態(tài)分析計(jì)算領(lǐng)域，其特征在于，根據(jù)國際《測量不確定度評定與表示》，確定各測點(diǎn)的擴(kuò)展不確定度，基于該擴(kuò)展不確定度，構(gòu)造了測點(diǎn)正常率η指標(biāo)、η_2/3指標(biāo)、η_1/3指標(biāo)，再利用這些指標(biāo)對根據(jù)量測值進(jìn)行狀態(tài)估計(jì)計(jì)算得到的各個(gè)結(jié)果進(jìn)行評價(jià)，選擇這些指標(biāo)更大的狀態(tài)估計(jì)結(jié)果輸出。本發(fā)明具有合理性、實(shí)用性和可推廣性。文檔編號G06F17/50GK101615213SQ200910089479公開日2009年12月30日申請日期2009年7月21日優(yōu)先權(quán)日2009年7月21日發(fā)明者何光宇,鋒劉,梅生偉,王雅婷,董樹鋒申請人:清華大學(xué)