专利名称：：基于扩展不确定度的电力系统状态估计结果评价方法
技术领域：
：本技术发明属于电力系统分析计算领域。
背景技术：
：在测点真值未知情况下，如何对电力系统状态估计结果给出合理的评价方法，是进行状态估计研究首先并且必须要解决的问题。以往研究中仅仅给出了测点真值已知情况下的评价指标，主要应用于仿真研究，通过人为设定真值和量测值，进行状态估计研究，用范数和无穷范数对估计值和真值之间距离进行衡量，即又=2]1、-^1，5"2=maxIx,.—I其中x,，r,分别为状态估计结果:c的第/个分量的估计值和真值。&,S2分别用来衡量状态估计结果整体偏差和最大偏差。由此可以看出，以往的研究中，被测量真值已知情况下，人们主要通过测量误差来评价状态估计结果的优劣。而在实际应用中，被测量真值是无法准确获知的未知量，因而误差仅是理论上的概念，实际上是无法确切获知的，使用起来具有局限性。如何在被测量真值未知情况下对状态估计结果进行合理评价，是目前状态估计研究中的盲区，但却是进行状态估计首先且必须要解决的重要问题。由于认识到测量残差概念的不足，国际计量学界在上世纪九十年代引入测量不确定度概念，并建议用采用测量不确定度对测量结果进行评价。测量不确定度是定量说明测量结果质量的一个参数，反映了可能的误差分布范围。不确定度只与测量条件有关而与测量结果无关，可以近似地理解为一定置信概率下的误差限值。利用测量不确定度代替误差表示测量结果，易于了解、便于评定，具有合理性和实用性。测量不确定度分为标准不确定度和扩展不确定度两类。标准不确定度是用标准偏差表示测量结果的不确定度，分为用A类计算方法得到的标准不确定度、用B类计算方法得到的标准不确定度和根据其它一些量值求出的合成标准不确定度。扩展不确定度由合成标准不确定度^乘上覆盖因子A:而得到。A:值一般为2，有时为3，取决于被测量的重要性、效益和风险。被测量真值分布的大部分可望包含在扩展不确定度表示的区间内，因而本发明提出的状态估计结果评价方法是基于扩展不确定度进行的。人们可根据测量数据、有关材料和仪器的技术指标、检定证书和手册中提供的数据来定量评定扩展不确定度。我国于1999年，经国家质量技术监督局批准，颁布实施由全国法制计量技术委员会提出的《测量不确定度评定与表示》(JJF1059—1999)，依据此规范，不确定度的一般评定步骤为第l步建模；第2步分析影响测量的所有不确定度来源；第3步对符合A类评定的不确定度分量用A类评定方法，计算出A类标准不确定度；第4步将剩余的不确定度分量，分别用B类评定方法，计算出B类标准不确定度；第5歩求出合成标准不确定度；第6歩将合成标准不确定度^乘以置信因子k，求出扩展不确定度f7-yb/c。本发明提出了一种电力系统状态估计结果评价方法。基于扩展不确定度，给出了正常测点和异常测点的定义，进而构建了测点正常率指标和;72/3、/71/3指标，利用这些指标进行状态估计结果评价。同时，给出了衡量状态估计结果合理性的具体方法。
发明内容本发明的目的在于提供科学合理的电力系统状态估计结果评价方法。本发明的特征在于所述方法是在调度中心的一个电力系统状态估计结果评价计算机上依次按以下步骤实现的歩骤(l)，所述电力系统状态估计结果评价计算机从能量管理系统中获得电力系统的网络参数，包括输电支路的各支路号、首节点编号、末节点编号、串联电阻、串联电抗、并联电导、并联电纳、以及变压器的变比和阻抗；步骤(2)，从数据采集和监视系统SCADA的网络数据服务器获得遥测量测值和遥信量测值，其中遥测量测值包括电压幅值、发电机有功功率、发电机无功功率、负荷有功功率、负荷无功功率、支路首端有功功率、支路首端无功功率、支路末端有功功率、支路末端无功功率、支路电流幅值；遥信量测值包括开关、刀闸开合状态以及变压器分接头位置，开关或刀闸闭合的量测值为i，断开的量测值为o;步骤p)，构造计算模型歩骤(3.1)，遍历所述各开关、支路、变压器、发电机、电容器以及负荷端口各节点，根据得到的各开关量把电网中通过闭合的开关或刀闸连接在一起的节点称为计算用节点；步骤(3.2)，把歩骤(2)得到的各计算用节点的电压幅值、发电机功率以及负荷功率的量测值按以下方法分配到对应的所述计算用节点上若所述计算用节点上有多个电压幅值量测值，则取平均值作为一个电压幅只有所述计算用节点上所有的发电机和负荷都有功率量测时，才能把发电机功率减去负荷功率作为该计算用节点上的注入功率量测值；把各计算用节点上的电压幅值、各计算用节点上的注入功率以及各支路上的潮流这些量测量定义为测点；歩骤(4)，从数据库中获取测量不确定度信息，其中包括测点序号、量测量、标准差以及扩展不确定度c/;歩骤(5)，基于步骤(1)和步骤(2)得到的数据，运行调度中心的n套状态估计软件，得到n组实时状态估计结果V,x2,…,x"，所述实时状态估计的状态变量为计算用节点电压的幅值v和相角e，所述n的值是任意预设的；歩骤(6)，根据步骤(5)得到的状态估计结果x^，…，x"，分别计算出不同状态结果下各测点的估计值，在状态估计结果/下，各测点的估计值为/v,()、~()、&y()、"1,2…,"，计算用节点/的电压幅值为，/v,(Z)=V;计算用节点/的有功注入功率为，卢计算用节点z'的无功注入功率为，=I/(Gysin(《-《)-^cos(《-;其中，/v,(/)、；^,()和/^(x"统一用^(x"表示，&()中&代表计算用节点/上的量测量g，g是计算用节点z的电压幅值K或注入有功功率s、注入无功功率G;从计算用节点/到计算用节点y的支路/-7'上/侧有功潮流为，=(v,"、-v,v《c一,-o-"v尽"<);从计算用节点/到计算用节点y的支路,'-y上/侧无功潮流为，=v,VGysin(《-O",VAco);其中，(/)和^(/)统一用、(/)表示，1(/)中^代表从计算用节点/到计算用节点7'的支路/-7'上/侧的量测量^，^是/侧有功潮流《或无功潮流0/;从计算用节点/到计算用节点y的支路/-/上y侧有功潮流为，从计算用节点/到计算用节点y的支路卜7上y侧无功潮流为，=一(v/)2(A+sin(^—0+^cos(W;其中，、,(/)和/^(x"统一用/^(/)表示，K)中&代表从计算用节点/到计算用节点y的支路/-y'上y侧的量测量g，g是y侧有功潮流s,或无功潮流上述量测方程中，《，《分别为状态估计结果^下计算用节点z'和计算用节点y的电压幅值的估计值，6f，《分别为状态估计结果/下计算用节点/和计算用节点y电压相角的估计值，q为计算用节点/和计算用节点"司的电导，&为计算用节点/和计算用节点y间的电纳，^为支路w的"形等值电路中计算用节点/和计算用节点乂的接地电纳^;步骤(7)，对步骤(5)得到的所述状态估计结果x',x2，…,？的合理性进行评价歩骤(7.1)，对于测点g，g为上述测点g,，或&,的任意一种，在状态估计结果/下，若有l^(v)-zg|st/g，其中，4为该测点g的量测值，~(/)为该测点g的估计值，t/g为测点g的扩展不确定度，则在所述状态估计结果/下，测点g为正常测点，否则，所述测点g为异常测点；定义状态估计结果/对应的测点正常率为—/;*100%，其中，r为系统中总的测点数目，^为状态估计结果？对应的正常测点数目步骤(7.2)，在状态估计结果/下，若有/个测点满足|\(/)-&/s^，则定义状态估计结果？的;/2/3指标为义4=/〃*，/0,步骤(7.3)，在状态估计结果/下，若有,个测点满足|~(/)-4则定义状态估计结果/的//1/3指标为^，^步骤(7.4)，根据歩骤(7.1)，计算并比较所述各状态估计结果？,x2,…,:c"对应的测点正常率;71,72,…,^，测点正常率更大的状态估计结果更合理；歩骤(7.5)，对于步骤(7.4)所述的n个测点正常率;/,;72,…,;7",若有/个状态估计结果对应的测点正常率相同，^w，则对于这/个状态估计结果，继续比较它们的&。指标，/72/3指标更大的状态估计结果更合理；步骤(7.6)，对于步骤(7.5)所述的/个状态估计结果中，若仍有m个状态估计结果的72/3指标相同，则继续比较这附个状态估计结果的;71/3指标，仏/3指标更大的状态估计结果更合理；步骤(7.7)，对于步骤(7.6)所述的m个状态估计结果中，若仍有w个状态估计结果的/71/3指标相同，w《m，则认为该w个状态估计结果的合理性相同。本发明通过构造测点正常率/7指标、"2/3指标、％3指标，给出了科学的电力系统状态估计结果评价方法，可以对各实时状态估计结果进行合理性评价，并输出合理性最高的状态估计结果。该评价方法计算简单，具有实用性和可推广性的优点。图1是基于扩展不确定度的电力系统状态估计结果评价流程图。图2是支路的;r形等值电路图。图3是变压器的;r形等值电路图。图4是算例中的四节点系统量测图0，有功无功量测对，A，电压量测。图5是算例中的四节点系统计算模型图。具体实施例方式本发明按照以下步骤进行歩骤(l)，从能量管理系统中获得电力系统的网络参数，包括输电支路的各支路号、首节点编号、末节点编号、串联电阻、串联电抗、并联电导、并联电纳、以及变压器的变比和阻抗。步骤(2)，从数据采集和监视系统SCADA的网络数据服务器获得遥测量测值和遥信量测值，其中遥测量测值包括电压幅值、发电机有功功率、发电机无功功率、负荷有功功率、负荷无功功率、支路首端有功功率、支路首端无功功率、支路末端有功功率、支路末端无功功率、支路电流幅值；遥信量测值包括开关、刀闸开合状态以及变压器分接头位置，开关或刀闸闭合的量测值为1，断开的量测值为0;所述SCADA系统采用轮询方式采集数据，轮询周期为510s。离线状态下可以通过文件给出上述量测量，在线应用时可以通过SCADA系统提供的API或使用TCP/IP通信协议从SCADA系统中获得。步骤(3)，构造计算模型歩骤(3.1)，把开关、支路、变压器、发电机、电容器以及负荷端口定义为节点，遍历所述这些节点，根据得到的各开关量把电网中通过闭合的开关或刀闸连接在一起的节点称为计算用节点，下面所述节点均指计算用节点；步骤(3.2)，把步骤(2)得到的各计算用节点的电压幅值、发电机功率以及负荷功率的量测值按以下方法分配到对应的所述计算用节点上若所述计算用节点上有多个电压幅值量测值，则取平均值作为一个电压幅值量测值；只有所述计算用节点上所有的发电机和负荷都有功率量测时，才能把发电机功率减去负荷功率作为该计算用节点上的注入功率量测值；把各计算用节点上的电压幅值、各计算用节点上的注入功率以及各支路上的潮流这些量测量定义为测点。歩骤(4)，从数据库中获取测量不确定度信息，其中包括测点序号、量测量、标准差以及扩展不确定度t/。歩骤(5),基于歩骤(1)和步骤(2)得到的数据，运行调度中心的n套状态估计软件，得到n组实时状态估计结果;c'，;c2，…,x"，所述n的值是任意预设的。对于输电网状态估计而言，实时状态估计的状态变量为节点电压的幅值v和相角步骤(6)，根据步骤(5)得到的状态估计结果x',;c2，…,x"，分别计算出不同状态结果下各测点的估计值，在状态估计结果Z下，各测点的估计值为/v,()、Z^,(x"、1(/)、"1，2…，"，节点/的电压幅值为，/v,(X;节点/的有功注入功率为，/^,(/)=2>>/(GyCOS(6^-《)+^sin(^-O);卢节点/的无功注入功率为，/i!e,(/)=^/(GySin(《-^cos(《-e/));其中，/v,(Z)、~,(^)和^(^)统一用^(^)表示，/^(/)中g,代表节点z'上的量测量g，g是节点/的电压幅值K或有功注入功率《、无功注入功率g;支路"_/上/侧有功潮流为，支路"y上/侧无功潮流为，1(/)=-—Gysin(《—e/)+v,，yC0-0;其中，(V)和/^(x"统一用/^(x"表示，/^(/)中&代表支路卜y上f侧的量测量g，g是z'侧有功潮流A或无功潮流仏；支路"y上7侧有功潮流为，=(v/)2《-v>Xco《—o+A/A-O;支路"/上_/侧无功潮流为，、(？)=-(v/)2(5yU+v,、Xsi《—《)+vyCO《-e/);其中，、,(？)和/^,(/)统一用\,(/)表示，^,(/)中g,代表支路/-y'上y侧的量测量g，g是y侧有功潮流&或无功潮流仏；上述量测方程中，《，《分别为状态估计结果/下节点/和节点/的电压幅值的估计值，《，《分别为状态估计结果z下节点/和节点y电压相角的估计值，q为节点/和节点y间的电导，A为节点/和节点/间的电纳，x为支路的;r形等值电路中节点/和节点y的接地电纳^。支路zV的;r形等值电路如图2所示。歩骤(7),对步骤(5)得到的状态估计结果x、x2，…,x"的合理性进行评价对于测点g，g为上述测点g,，^或^的任意一种，真值落在(Z「C/g,^+C/g)中概率为/，良口式中，f为系统真实状态，4为该测点g的量测值，K为与置信概率p对应的测点g的扩展不确定度，可通过上面介绍的扩展不确定度的评定方法获取。对于测点g，在状态估计结果/下，若有则在所述状态估计结果/下，测点g为正常测点，否则，所述测点g为异常测占八"o设系统中总的测点数目为r,状态估计结果z对应的正常测点数目为y，则状态估计结果/对应的测点正常率为依次类推，可计算出状态估计结果1^2,...,对应的测点正常率分别为;/,/72,...，/7"。根据测点正常率对状态估计结果做出评价，测点正常率越高，状态估计结果越合理，即越接近系统的真实状态。由于正常测点的定义与扩展不确定^相关，而C/g又与一定的置信概率p相关，因而测点正常率是相应于某一固定的置信概率/提出的。若状态估计结果/，x2，…,x"中，有/个状态估计结果的正常率相同，可认为在置信概率p下衡量时，这/个状态估计结果合理性相同。而在另一置信概率下衡量时，这/个状态估计结果未必具有同等合理性。此时，需要引入/72/3指标，仏/3指标进一步判别。对于上述的状态估计结果/，若在总的r个测点中，有/个测点(/《^)满3g则定义72,3指标为《/3=//"100%同理，若在总的r个测点中，有f个测点(f、/)满足W。—4丄t/3g则定义/71/3指标为<=,/"100%对于测点正常率相同的/个状态估计结果，通过比较72/3指标进一步判别状态估计结果的合理性。72/3指标更高的状态估计结果更合理。如果这/个状态估计结果中，仍有附个状态估计结果的;/2/3指标相同，附W，则继续比较这m个状态估计结果的/71/3指标。71/3指标更高的状态估计结果更合理。如果这m个状态估计结果中，仍有w个状态估计结果的"。指标相同，WSm，则认为这w个状态估计结果合理性相同。所述方法中当支路中含变压器时，支路两端节点侧的潮流的量测方程如下变压器支路W上节点/侧有功潮流为，变压器支路&上节点z'侧无功潮流为，=-如""、+vco《-o;变压器支路zV上节点y侧有功潮流为，变压器支路上节点y'侧无功潮流为，=-)、十会v,V-《)；上述量测方程中，K为变压器的非标准变比，^为变压器标准侧的电纳。变压器的v形等值电路如图3所示。以图4所示的4节点系统为例，具体步骤说明如下步骤1:初始化从现有的能量管理系统中获得电力系统的网络参数。包括支路的串联电阻、串联电抗、并联电导和并联电纳；变压器的变比和阻抗。图4所示的4节点系统支路参数如表1所示表l四节点系统支路参数<table>tableseeoriginaldocumentpage21</column></row><table>步骤2:获得SCADA量测值某一个时刻4节点系统的量测数据如表2所示，其中遥信量的量测值为1表示开关闭合，0表示开关断开，遥测量的量测值使用有名值有功功率的单位是MW,无功功率的单位是MVAR，电压幅值的单位是kV。表2四节点系统量测数据<table>tableseeoriginaldocumentpage22</column></row><table>步骤4:从数据库中获取测点不确定度信息，如表4所示,表4四节点系统测点标准差以及扩展不确定度测点序号量测量标准差扩展不确定度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:运行状态估计软件，获取状态估计结果。运行状态估计软件，获取状态估计结果x1、x2、x3、x4、jc5，如表中电压的单位为kV，相角的单位为弧度。所示，表表5四节点系统状态估计结果状态估计结果？状态估计结果x2状态估计状态估计结果x4状态估计结果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:计算所有量测量的估计值各量测量的量测方程如下节点/注入功率Pi，Qi:J。ye,支路/:/上节点/侧潮流《=v,2g-v,"gcos《-v,v,sin《支路/:/上节点y侧潮流=—、(6+Jc)+v,Vjgsin6>y+v,6cosPy变压器支路&上节点/侧潮流变压器支路上节点/侧潮流节点/电压幅值上述量测方程中，v,，v,分别为节点i和节点j的电压幅值，《为节点i和节点j电压相角差，Gy，^为导纳矩阵元素，g，b，尺为支路;r形模型中的参数，K为变压器的非标准变比，6,为变压器标准侧的电纳。将步骤5计算得到的状态变量代入到量测方程中，计算出各量测量的估计值"/^)，求得的结果如表6所示。表6四节点系统量测量估计结果量测量量测值状态估计结果x1状态估计结果x2状态估计结果x3状态估计结果:c4状态估计结果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:比较状态估计结果的合理性首先，根据步骤3得到的测点量测值、步骤4得到的测点不确定度、步骤6得到的测点估计值，分别计算状态估计结果？，x2,x3,x4,对应的测点正常率7指标，计算结果如表7所示。比较/7指标，状态估计结果^的合理性最差。x'乂,x4,x5的合理性比较转到下一步。其次，改变扩展不确定度取值为lc/时，计算出此时状态估计结果5的;72/3指标，计算结果如表7所示。比较;72/3指标，x的合理性最好。、x4和？的合理性比较转到下一步。最后，改变扩展不确定度取值为^7时，计算出此时状态估计结果xUW4乂的/^指标，计算结果如表7所示。比较//1/3指标，可以发现x3、x4和？的合理性顺序为X4>>x5。故最终所得合理性排序为；c1>/>>x5>x2。表7四节点系统状态估计结果的评价指标^一状态估计状态估计状态估计状态估计状态估计指标结果x1结果x2结果x3结果x4结果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%权利要求1.基于扩展不确定度的电力系统状态估计结果评价方法，其特征在于，所述方法是在调度中心的一个电力系统状态估计结果评价计算机上依次按以下步骤实现的步骤(1)，所述电力系统状态估计结果评价计算机从能量管理系统中获得电力系统的网络参数，包括输电支路的各支路号、首节点编号、末节点编号、串联电阻、串联电抗、并联电导、并联电纳、以及变压器的变比和阻抗；步骤(2)，从数据采集和监视系统SCADA的网络数据服务器获得遥测量测值和遥信量测值，其中遥测量测值包括电压幅值、发电机有功功率、发电机无功功率、负荷有功功率、负荷无功功率、支路首端有功功率、支路首端无功功率、支路末端有功功率、支路末端无功功率、支路电流幅值；遥信量测值包括开关、刀闸开合状态以及变压器分接头位置，开关或刀闸闭合的量测值为1，断开的量测值为0；步骤(3)，构造计算模型步骤(3.1)，遍历所述各开关、支路、变压器、发电机、电容器以及负荷端口各节点，根据得到的各开关量把电网中通过闭合的开关或刀闸连接在一起的节点称为计算用节点；步骤(3.2)，把步骤(2)得到的各计算用节点的电压幅值、发电机功率以及负荷功率的量测值按以下方法分配到对应的所述计算用节点上若所述计算用节点上有多个电压幅值量测值，则取平均值作为一个电压幅值量测值；只有所述计算用节点上所有的发电机和负荷都有功率量测时，才能把发电机功率减去负荷功率作为该计算用节点上的注入功率量测值；把各计算用节点上的电压幅值、各计算用节点上的注入功率以及各支路上的潮流这些量测量定义为测点；步骤(4)，从数据库中获取测量不确定度信息，其中包括测点序号、量测量、标准差以及扩展不确定度U；步骤(5)，基于步骤(1)和步骤(2)得到的数据，运行调度中心的n套状态估计软件，得到n组实时状态估计结果x1，x2，…，xn，所述实时状态估计的状态变量为计算用节点电压的幅值v和相角θ，所述n的值是任意预设的；步骤(6)，根据步骤(5)得到的状态估计结果x1，x2，…，xn，分别计算出不同状态结果下各测点的估计值，在状态估计结果xk下，各测点的估计值为hVi(xk)、hPi(xk)、hQi(xk)、hPij(xk)、hQij(xk)、hPji(xk)、hQji(xk)，k＝1，2…，n，计算用节点i的电压幅值为，hVi(xk)＝vik；计算用节点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>Σ</mi><mrow><mi>j</mi><mo>∈</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>θ</mi><mi>i</mi></msub><mi>k</mi></msup><mi></mi><mo>-</mo><msup><msub><mi>θ</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>θ</mi><mi>i</mi></msub><mi>k</mi></msup><mi></mi><mo>-</mo><msup><msub><mi>θ</mi><mi>j</mi></msub><mi>k</mi></msup><mi></mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>;</mo></mrow>]]></math></maths>计算用节点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>Σ</mi><mrow><mi>j</mi><mo>∈</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>θ</mi><mi>i</mi></msub><mi>k</mi></msup><mi></mi><mo>-</mo><msup><msub><mi>θ</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>θ</mi><mi>i</mi></msub><mi>k</mi></msup><mi></mi><mo>-</mo><msup><msub><mi>θ</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)统一用hgi(xk)表示，hgi(xk)中gi代表计算用节点i上的量测量g，g是计算用节点i的电压幅值Vi或有功注入功率Pi、无功注入功率Qi；从计算用节点i到计算用节点j的支路i-j上i侧有功潮流为，hPij(xk)＝(vik)2Gij-vikvjkGijcos(θik-θjk)-vikvjkBijsin(θik-θjk)；从计算用节点i到计算用节点j的支路i-j上i侧无功潮流为，hQij(xk)＝-(vik)2(Bij+yc)-vikvjkGijsin(θik-θjk)+vikvjkBijcos(θik-θjk)；其中，hPij(xk)和hQij(xk)统一用hgij(xk)表示，hgij(xk)中gij代表从计算用节点i到计算用节点j的支路i-j上i侧的量测量g，g是i侧有功潮流Pij或无功潮流Qij；从计算用节点i到计算用节点j的支路i-j上j侧有功潮流为，hPji(xk)＝(vjk)2Gij-vikvjkGijcos(θik-θjk)+vikvjkBijsin(θik-θjk)；从计算用节点i到计算用节点j的支路i-j上j侧无功潮流为，hQji(xk)＝-(vjk)2(Bij+yc)+vikvjkGijsin(θik-θjk)+vikvjkBijcos(θik-θjk)；其中，hPji(xk)和hQji(xk)统一用hgji(xk)表示，hgji(xk)中gji代表从计算用节点i到计算用节点j的支路i-j上j侧的量测量g，g是j侧有功潮流Pji或无功潮流Qji；上述量测方程中，vik，vjk分别为状态估计结果xk下计算用节点i和计算用节点j的电压幅值的估计值，θik，θjk分别为状态估计结果xk下计算用节点i和计算用节点j电压相角的估计值，Gij为计算用节点i和计算用节点j间的电导，Bij为计算用节点i和计算用节点j间的电纳，yc为支路i-j的π形等值电路中计算用节点i和计算用节点j的接地电纳yc；步骤(7)，对步骤(5)得到的所述状态估计结果x1，x2，…，xn的合理性进行评价步骤(7.1)，对于测点g，g为上述测点gi，gij或gji的任意一种，在状态估计结果xk下，若有|hg(xk)-Zg|≤Ug，其中，Zg为该测点g的量测值，hg(xk)为该测点g的估计值，Ug为测点g的扩展不确定度，则在所述状态估计结果xk下，测点g为正常测点，否则，所述测点g为异常测点；定义状态估计结果xk对应的测点正常率为ηk＝qk/r*100％，其中，r为系统中总的测点数目，qk为状态估计结果xk对应的正常测点数目；步骤(7.2)，在状态估计结果xk下，若有sk个测点满足<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>≤</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，则定义状态估计结果xk的η2/3指标为<mathsid="math0004"num="0004"><math><![CDATA[<mrow><msubsup><mi>η</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)，在状态估计结果xk下，若有tk个测点满足<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>≤</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，则定义状态估计结果xk的η1/3指标为<mathsid="math0006"num="0006"><math><![CDATA[<mrow><msubsup><mi>η</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)，根据步骤(7.1)，计算并比较所述各状态估计结果x1，x2，…，xn对应的测点正常率η1，η2，…，ηn，测点正常率更大的状态估计结果更合理；步骤(7.5)，对于步骤(7.4)所述的n个测点正常率η1，η2，…，ηn，若有l个状态估计结果对应的测点正常率相同，l≤n，则对于这l个状态估计结果，继续比较它们的η2/3指标，η2/3指标更大的状态估计结果更合理；步骤(7.6)，对于步骤(7.5)所述的l个状态估计结果中，若仍有m个状态估计结果的η2/3指标相同，m≤l，则继续比较这m个状态估计结果的η1/3指标，η1/3指标更大的状态估计结果更合理；步骤(7.7)，对于步骤(7.6)所述的m个状态估计结果中，若仍有w个状态估计结果的η1/3指标相同，w≤m，则认为该w个状态估计结果的合理性相同。2.根据权利要求1所述的基于扩展不确定度的电力系统状态估计结果评价方法，其特征在于，所述支路中含变压器时，支路两端的计算用节点侧的潮流的量测方程如下从计算用节点/到计算用节点7的变压器支路W上节点/侧有功潮流为，从计算用节点z到计算用节点乂的变压器支路上节点z侧无功潮流为，1"4)=-如。、-《)；从计算用节点/到计算用节点7的变压器支路i-y上节点/侧有功潮流为，从计算用节点/到计算用节点j的变压器支路w上节点y侧无功潮流为，上述量测方程中，K为变压器的非标准变比，Zv为变压器标准侧的电纳。全文摘要基于扩展不确定度的电力系统状态估计结果评价方法属于电力系统状态分析计算领域，其特征在于，根据国际《测量不确定度评定与表示》，确定各测点的扩展不确定度，基于该扩展不确定度，构造了测点正常率η指标、η<sub>2/3</sub>指标、η<sub>1/3</sub>指标，再利用这些指标对根据量测值进行状态估计计算得到的各个结果进行评价，选择这些指标更大的状态估计结果输出。本发明具有合理性、实用性和可推广性。文档编号G06F17/50GK101615213SQ200910089479公开日2009年12月30日申请日期2009年7月21日优先权日2009年7月21日发明者何光宇,锋刘,梅生伟,王雅婷,董树锋申请人:清华大学