CN107169881A - 引入pmu支路电流相量的电力系统静态状态估计方法 - Google Patents
引入pmu支路电流相量的电力系统静态状态估计方法 Download PDFInfo
- Publication number
- CN107169881A CN107169881A CN201710416590.XA CN201710416590A CN107169881A CN 107169881 A CN107169881 A CN 107169881A CN 201710416590 A CN201710416590 A CN 201710416590A CN 107169881 A CN107169881 A CN 107169881A
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- pmu
- mfrac
- node
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000005259 measurement Methods 0.000 claims abstract description 65
- 239000011159 matrix material Substances 0.000 claims description 27
- 230000005540 biological transmission Effects 0.000 claims description 8
- 230000005611 electricity Effects 0.000 claims description 8
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 239000013589 supplement Substances 0.000 abstract description 2
- 238000002347 injection Methods 0.000 description 9
- 239000007924 injection Substances 0.000 description 9
- 238000005516 engineering process Methods 0.000 description 4
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 240000008042 Zea mays Species 0.000 description 1
- 235000005824 Zea mays ssp. parviglumis Nutrition 0.000 description 1
- 235000002017 Zea mays subsp mays Nutrition 0.000 description 1
- 235000005822 corn Nutrition 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000012948 formulation analysis Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种引入PMU支路电流相量的电力系统静态状态估计方法。本发明首先根据系统中已知的数据将量测方程列写出来。接着将PMU装置根据系统经过计算后达到可观的要求配置到系统中,然后将PMU得到的数据应用到系统中,本发明以PMU数据为主,辅以经过量测变换后的SCADA数据和PMU伪量测数据及类PMU数据对PMU可观测不足点的数据进行补充,使系统达到可观,从而使系统的估计精度得到提高。本发明可以在配置较少PMU数据的情况下得到较精确的数据,充分利用PMU量测数据,可以更好的进行状态估计。
Description
技术领域
本发明涉及电力系统中的状态估计方法,具体涉及的是引入PMU支路电流相量的电力系统静态状态估计方法。
背景技术
随着电网规模的不断扩大,电力系统中的网络结构日趋复杂,电力系统调度中心的自动化水平也需要逐步向高级发展,能量管理系统(Energy Management System,EMS)得到了广泛应用。其中状态估计是能量管理系统的核心功能模块,基于电网的结构、参数以及实时量测,为其他高级应用软件提供一个可靠而完整的电力系统实时数据库,是电力系统运行、控制和安全评估等工作的基础。
在基于全球定位系统(Global Positioning System,GPS)的相量量测装置(Phasor Measurement Unit,PMU)量测技术应用于电力系统之前,状态估计的测量数据一直来自于数据采集与监视控制系统(Supervisory Control And Data Acquisition,SCADA)系统。而PMU量测可以直接测量节点电压和相角,一旦电网所有节点均配置有PMU,状态估计的问题将大大简化,基于SCADA的状态估计将失去意义。
当前利用PMU量测的状态估计模型一般或者只利用PMU量测进行线性估计,或者在SCADA量测的基础上增加PMU电压幅值和相角量测进行非线性估计,前者在当前PMU不可观测配置的情况下难以进行,后者没有利用PMU支路电流相量量测。支路电流幅值量测在输电网状态估计中应用较少,原因是:①电流幅值量测不易进行有功和无功解耦;②电流幅值量测不能提供方向信息,会产生多解性问题;③线路电纳较小时,利用电流幅值量测形成雅克比矩阵可能会由于相应的雅克比举证元素较小而产生数值问题。总之,在输电网状态估计中直接利用电流幅值量测存在不便,而支路电流的相角量测还无法直接利用。
发明内容
本发明的目的是提供一种计及PMU支路电流相量的状态估计方法。利用PMU电压相量和电流相量量测计算出相关支路的等值功率量测值,然后在状态估计中加以利用,有利于提高状态估计的精度,从而保证电网能够更加稳定的运行。
本发明解决其技术问题所采用的技术方案是:这种计及PMU支路电流的状态估计方法如下:
一、读取EMS系统中的电网数据,包括:输电线支路号、节点编号、变压器变比、开关位置、网络参数;
二、读取系统中的实时量测值,包括:节点有功量测、节点无功量测,支路有功量测、支路无功量测,列写系统量测方程如下:
z=h(x)+v;
目标函数方程:J(x)=[z-h(x)]TR-1[z-h(x)];
三、初始化,包括:对状态量设置初值,节点PMU进行配置,获得节点导纳矩阵,设置迭代门限;
四、输入试量测量z,包括SCADA系统中节点有功量测、节点无功量测,支路有功量测、支路无功量测;PMU装置测得的节点电压幅值和相角,节点支路电流幅值和相角;
五、将电流相量变换为支路潮流,利用PMU电压相量和电流相量量测计算出相关的等值功率量测值,然后再状态估计中加以利用:
其中,Pij和Qij分别为节点i与节点j之间的支路ij的有功功率和无功功率,Ii为节点i的电流,和分别为节点i的电压相角和电流相角;Pij和Qij的权值按照如下误差传递公式计算;
其中, 为等值有功量测和等值无功量测的误差方差, 分别为电压幅值、电流幅值、电压相角、电流相角量测误差的标准差;和
将电流相量变换为PMU相关节点电压相量,利用PMU节点电压相量和支路电流相量,推算出相关节点的电压相量,为:
其中,为节点j的相电压,为支路ij的电流向量,Yi0为节点i的对地电容,Yij为节点导纳矩阵的元素;的权值按照如下误差传递公式计算;
其中,为等值电压向量量测的误差方差, 和分别为电压幅值Ui、电压相角θ、支路电流Iij和支路电流相角对应的标准差;的权值表示为
转化的支路功率和电压相量代入到量测方程中,进行迭代;
检验目标函数J(x)=[z-h(x)]TR-1[z-h(x)]值是否小于收敛标准,若满足则状态估计结束,得到可观测下的状态估计误差协方差和估计误差总方差:
估计误差协方差:
估计误差总方差:
若不满足返回步骤五继续进行迭代;
六、将此估算结果送入系统数据库中,供电力系统调度中心的工作人员随时估计出电力系统当前的运行状态,保证系统正常运行。
与现代技术相比,本发明的有益效果在于:
1、本发明提供的电力系统状态估计方法,原理可靠,和传统以SCADA数据为主的基于同步相量状态估计相比具有对PMU数据利用率高,符合电力系统中PMU装置数量逐渐增多的发展规律;同时将PMU支路电流相量变换为支路潮流以及相关电压相量,更有利于估计精度的提高。本发明提供的方法的估计精度较高,估计稳定性好,可以满足电力系统状态估计要求,确保电力系统的安全、稳定、经济运行。
2、本发明可以缓解扩展卡尔曼滤波在电力系统运行正常时,状态估计结果较为准确,但在电力系统的参数或运行状态发生了变化的时候,导致滤波初期估计协方差下降太快,从而滤波不稳定甚至发散,造成估计结果的不准确等现象,可以在电力系统发生波动时获得较好的状态估计值,为电力系统的稳定运行提供实时数据的保障。
附图说明
图1是引入PMU支路电流相量的电力系统静态状态估计方法的流程图;
图2是本发明采用的π型等值电路图。
具体实施方式
下面结合附图对本发明做进一步的详细说明:
结合图1、图2所示,这种引入PMU支路电流相量的电力系统静态状态估计方法如下:
将网络中已知的数据输入其中,包括结线信息状态。
对系统进行初始化,为进行新的一次状态估计做准备。
输入SCADA系统中测得的量测量,同时将PMU装置测得的数据同时输入到系统中,使得SCADA系统的数据和PMU的数据进行混合。
根据PMU配置的数量和位置判断系统是否达到可观,如果是则进行模型的假设和列写目标函数。
在上述系统中,给出了系统支路的参数、网络的结线和量测系统,非线性的量测方程改写为:
z=h(x)+v (1)
已知量测矢量z以后,状态估计矢量为了使目标函数:
J(x)=[z-h(x)]TR-1[z-h(x)] (2)
达到最小x的值。式中z为量测量矢量,h(x)为量测量函数矢量,x=e+jf为状态矢量,v为量测误差,R-1为系统量测权重矩阵,对x进行线性化后得:
h(x)≈h(x0)+H(x0)△x (3)
H(x)是m×n阶量测矢量的雅克比矩阵。
最终得到目标函数为:
J(x)=[△z-H(x0)△x]TR-1[△z-H(x0)△x] (4)
对其进行配方分析得:
(H-1R-1H)△x=HTR-1(z-h(x)) (5)
为量测量的雅克比矩阵,HTR-1H为状态估计的信息矩阵。
线性状态估计方程为:
PMU设备装配系统将是一个渐进的过程,因此,在相当长的一段时间内,PMU量测对系统还将是不完全可观测的,在这种情况下需要利用SCADA量测补充PMU不可观测的区域,进行线性状态估计具体过程如下:
使用上标O表示PMU量测可观测区域,上标U表示PMU量测不可观测区域,则节点电压方程可表示为:
其中:和分别为可观测和不可观测区域的节点电压矩阵;和分别为可观测和不可观测区域的节点注入电流矩阵;E为单位矩阵;Yi OO为节点注入可观自导纳矩阵,Yi OU为节点注入可观与不可观互导纳矩阵,Yi UO为节点注入不可观和可观互导纳矩阵,Yi UU为节点注入不可观自导纳矩阵;为支路可观电流矩阵,为支路可观自导纳矩阵,为支路可观与不可观互导纳矩阵。
对不可观测区域只采用了节点电压和节点注入电流的理由:
1)节点的注入量是决定系统状态的因变量,而支路电流(或功率)只是描述系统状态的从变量,采用节点电压和节点注入电流已完全可以描述不可观测区域的状态了。
2)不可观测区域的节点电压和节点注入电流可以方便的由状态估计结果得到。
3)对PMU不可观测区域补偿节点电压和节点注入电流量测后,已达到对系统的完全可观测性。
4)不可观测区域的节点注入电流矩阵表示为:
其中,Pi和Q i为注入的有功和无功功率,为的共轭阵;这里要用到节点电压的实部和虚部,而电压的实部和虚部要由压的幅值和相角求出。
如果对不可观测区域使用SCADA量测则得不到电压相角,也就无法求出节点注入电流,本发明对节点电流公式进行变型以便利用SCADA量测。
设节点等效注入可观导纳矩阵Yi U表示为:
将式(9)带入下式,有
整理得:
变型后的节点电流方程为:
把该节点电流方程作为线性计算的线性方程,带入式(6)进行线性状态估。
进行非线性状态估计中,分别引入PMU电压量测信息和相角量测信息、PMU功率量测信息以及PMU电流量测信息进行非线性状态估计。
引入PMU电压量测信息和相角量测信息进行非线性状态估计具体过程为:
1)直接将PMU节点电压相量量测方程加入到非线性状态估计中,PMU电压和相角量测方程为:
其中,和分别为节点i的PMU电压量测信息和相角量测息,Ui和θi为节点i的电压和相角,对应的雅克比矩阵的元素为1,其他的元素均为0;这种模型与传统状态估计模型中母线电压幅值测量值的用法完全相同,当在某条母线配置PMU时,量测方程中增加了上述2个方程以后,量测雅可比矩阵增2行,且每行只有一个取值为1的非零元素。
2)引入功率相角差;
对于含有PMU相角量测的大电网非线性状态估计中,引入新量测量θij,其中节i与节点j之间的θij=θi-θj,θj为节点j的相角;则有
这种对雅可比矩阵的修改非常简单,可以直接用于状态估计中。采用相角差的方法时,将不存在参考点的问题,直接量测(或计算)相角差都不依靠参考点。由于θij为节点之间的相对相角差,可以消除由于PMU和SCADA选取平衡节点选取不同,重新归算产生的误差,还可以消除PMU量测设备的由于共模因素引起的电压相角误差。θij主要约束节点间的相角差的作用。
以上两种方法中,对于归算误差或共模误差较小的PMU相角量测,可以直接用方法1),对于误差较大的相角量测,使用方法2)。
引入PMU功率量测信息进行非线性状态估计具体过程为:
其中,Pi m和分别为节点i的有功功率和无功功率量测信息,Pi和Qi分别为节点i的有功功率和无功功率;PMU功率量测按节点注入功率和支路功率的类型,其对应的雅克比矩阵系数与SCADA量测功率量对应的系数一致。
直接引入支路存在困难,原因在于各支路电流是电压相角θi和θj的函数,电流向量随着各个节点电压的相角变化而变化,一般情况下,各个节点相角θi和θj不能简单地考虑为0,无法用PQ解耦计算,给计算带来了困难。
引入PMU电流量测信息进行非线性状态估计具体过程分为以下两种情况:
将电流量转化为支路潮流,有
其中,Pij和Qij分别为节点i与节点j之间的支路ij的有功功率和无功功率,Ii为节点i的电流,和分别为节点i的电压相角和电流相角;Pij和Qij的权值按照如下误差传递公式计算;
其中, 为等值有功量测和等值无功量测的误差方差, 分别为电压幅值、电流幅值、电压相角、电流相角量测误差的标准差;和
将电流相量变换为PMU相关节点电压相量,利用PMU节点电压相量和支路电流相量,推算出相关节点的电压相量,为:
其中,为节点j的相电压,为支路ij的电流向量,Yi0为节点i的对地电容,Yij为节点导纳矩阵的元素;的权值按照如下误差传递公式计算;
其中,为等值电压向量量测的误差方差, 和分别为电压幅值Ui、电压相角θ、支路电流Iij和支路电流相角对应的标准差;的权值表示为
将此估算结果送入系统数据库中,供电力系统调度中心的工作人员随时估计出电力系统当前的运行状态,保证系统正常运行。
基于同步相量的状态估计,是在各个节点加入了电压、电流相直接可测的PMU量测量,由于技术和经济原因,不能任何节点都安装PMU装置,所以本发明是一种基于SCADA和PMU混合量测的静态状态估计方法,进而获得更高的估计精度和估计速度,通过将加入PMU装置的PMU数据和SCADA系统的数据进行融合,并且利用支路电流来得到支路功率以及相邻节点电压相量,使得电网系统达到可观。这样可以在配置较少PMU数据的情况下得到较精确的数据,可以更好的进行状态估计。
综上,本发明是在电力系统中加入PMU支路电流相量后进行静态状态估计的一种方法,是一种在电网中加入PMU装置后进行电力系统状态估计的方法。本发明首先根据系统中已知的数据将量测方程列写出来。接着将PMU装置根据系统经过计算后达到可观的要求配置到系统中,然后将PMU得到的数据应用到系统中,本发明以PMU数据为主,辅以经过量测变换后的SCADA数据和PMU伪量测数据及类PMU数据对PMU可观测不足点的数据进行补充,使系统达到可观,从而使系统的估计精度得到提高。本发明易于在原有系统中配置少量PMU的基础上增加PMU数量后的状态估计。
Claims (3)
1.引入PMU支路电流相量的电力系统静态状态估计方法,其特征在于:
一、读取EMS系统中的电网数据;
二、读取系统中的实时量测值,列写系统量测方程如下:
z=h(x)+v;
目标函数方程:J(x)=[z-h(x)]TR-1[z-h(x)];
式中z为量测矢量,h(x)为量测量函数矢量,v为量测误差,R-1为系统量测权重矩阵;
三、初始化,包括:对状态量设置初值,节点PMU进行配置,获得节点导纳矩阵,设置迭代门限;
四、输入量测矢量z,包括SCADA系统中节点有功量测、节点无功量测,支路有功量测和支路无功量测;PMU装置测得的节点电压幅值和相角,节点支路电流幅值和相角;
五、将电流相量变换为支路潮流,利用PMU电压相量和电流相量量测计算出相关的等值功率量测值,然后再状态估计中加以利用:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<msub>
<mi>I</mi>
<mi>i</mi>
</msub>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>&theta;</mi>
<msub>
<mi>u</mi>
<mi>i</mi>
</msub>
</msub>
<mo>-</mo>
<msub>
<mi>&theta;</mi>
<msub>
<mi>i</mi>
<mi>i</mi>
</msub>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<msub>
<mi>I</mi>
<mi>i</mi>
</msub>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>&theta;</mi>
<msub>
<mi>u</mi>
<mi>i</mi>
</msub>
</msub>
<mo>-</mo>
<msub>
<mi>&theta;</mi>
<msub>
<mi>i</mi>
<mi>i</mi>
</msub>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,Pij和Qij分别为节点i与节点j之间的支路ij的有功功率和无功功率,Ui为节点i的电压,Ii为节点i的电流,和分别为节点i的电压相角和电流相角;Pij和Qij的权值按照如下误差传递公式计算:
<mrow>
<msub>
<mi>R</mi>
<msub>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<msub>
<mi>&sigma;</mi>
<msub>
<mi>u</mi>
<mi>i</mi>
</msub>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>I</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<msub>
<mi>&sigma;</mi>
<msub>
<mi>i</mi>
<mi>i</mi>
</msub>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&theta;</mi>
<msub>
<mi>u</mi>
<mi>i</mi>
</msub>
</msub>
</mrow>
</mfrac>
<msub>
<mi>&sigma;</mi>
<msub>
<mi>&theta;</mi>
<msub>
<mi>u</mi>
<mi>i</mi>
</msub>
</msub>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&theta;</mi>
<msub>
<mi>i</mi>
<mi>i</mi>
</msub>
</msub>
</mrow>
</mfrac>
<msub>
<mi>&sigma;</mi>
<msub>
<mi>&theta;</mi>
<msub>
<mi>i</mi>
<mi>i</mi>
</msub>
</msub>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>R</mi>
<msub>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<msub>
<mi>&sigma;</mi>
<msub>
<mi>u</mi>
<mi>i</mi>
</msub>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>I</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<msub>
<mi>&sigma;</mi>
<msub>
<mi>i</mi>
<mi>i</mi>
</msub>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&theta;</mi>
<msub>
<mi>u</mi>
<mi>i</mi>
</msub>
</msub>
</mrow>
</mfrac>
<msub>
<mi>&sigma;</mi>
<msub>
<mi>&theta;</mi>
<msub>
<mi>u</mi>
<mi>i</mi>
</msub>
</msub>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&theta;</mi>
<msub>
<mi>i</mi>
<mi>i</mi>
</msub>
</msub>
</mrow>
</mfrac>
<msub>
<mi>&sigma;</mi>
<msub>
<mi>&theta;</mi>
<msub>
<mi>i</mi>
<mi>i</mi>
</msub>
</msub>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,为等值有功量测和等值无功量测的误差方差,分别为电压幅值、电流幅值、电压相角、电流相角量测误差的标准差;和
将电流相量变换为PMU相关节点电压相量,利用PMU节点电压相量和支路电流相量,推算出相关节点的电压相量,为:
<mrow>
<msub>
<mover>
<mi>U</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>j</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mover>
<mi>I</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>Y</mi>
<mrow>
<mi>i</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>Y</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>U</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<mo>-</mo>
<msub>
<mi>Y</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,为节点j的相电压,为节点i的相电压,为支路ij的电流向量,Yi0为节点i的对地电容,Yij为节点导纳矩阵的元素;的权值按照如下误差传递公式计算;
<mrow>
<msub>
<mi>R</mi>
<msub>
<mover>
<mi>u</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>j</mi>
</msub>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mover>
<mi>U</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>j</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<msub>
<mi>&sigma;</mi>
<msub>
<mi>u</mi>
<mi>i</mi>
</msub>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mover>
<mi>U</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>j</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>I</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<msub>
<mi>&sigma;</mi>
<msub>
<mi>i</mi>
<mi>i</mi>
</msub>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mover>
<mi>U</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>j</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&theta;</mi>
<msub>
<mi>u</mi>
<mi>i</mi>
</msub>
</msub>
</mrow>
</mfrac>
<msub>
<mi>&sigma;</mi>
<msub>
<mi>&theta;</mi>
<msub>
<mi>u</mi>
<mi>i</mi>
</msub>
</msub>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mover>
<mi>U</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>j</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&theta;</mi>
<msub>
<mi>i</mi>
<mi>i</mi>
</msub>
</msub>
</mrow>
</mfrac>
<msub>
<mi>&sigma;</mi>
<msub>
<mi>&theta;</mi>
<msub>
<mi>i</mi>
<mi>i</mi>
</msub>
</msub>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,为等值电压向量量测的误差方差,和分别为电压幅值Ui、电压相角θ、支路电流Iij和支路电流相角对应的标准差;的权值表示为
转化的支路功率和电压相量代入到量测方程中,进行迭代;
检验目标函数J(x)=[z-h(x)]TR-1[z-h(x)]值是否小于收敛标准,若满足则状态估计结束,得到可观测下的状态估计误差协方差和估计误差总方差:
估计误差协方差:
估计误差总方差:
若不满足返回步骤五继续进行迭代;
六、将此估算结果送入系统数据库中,供电力系统调度中心的工作人员随时估计出电力系统当前的运行状态,保证系统正常运行。
2.根据权利要求1所述的引入PMU支路电流相量的电力系统静态状态估计方法,其特征在于:所述的电网数据包括输电线支路号、节点编号、变压器变比、开关位置和网络参数。
3.根据权利要求1所述的引入PMU支路电流相量的电力系统静态状态估计方法,其特征在于:所述的实时量测值包括节点有功量测、节点无功量测,支路有功量测、支路无功量测。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710416590.XA CN107169881A (zh) | 2017-06-06 | 2017-06-06 | 引入pmu支路电流相量的电力系统静态状态估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710416590.XA CN107169881A (zh) | 2017-06-06 | 2017-06-06 | 引入pmu支路电流相量的电力系统静态状态估计方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107169881A true CN107169881A (zh) | 2017-09-15 |
Family
ID=59824430
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710416590.XA Pending CN107169881A (zh) | 2017-06-06 | 2017-06-06 | 引入pmu支路电流相量的电力系统静态状态估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107169881A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107665112A (zh) * | 2017-09-29 | 2018-02-06 | 北京紫优能源科技有限公司 | 一种电力控制流程的组态方法及装置 |
CN110543720A (zh) * | 2019-09-03 | 2019-12-06 | 北京交通大学 | 基于sdae-elm伪量测模型的状态估计方法 |
CN113659562A (zh) * | 2021-07-01 | 2021-11-16 | 国网江苏省电力有限公司电力科学研究院 | 可控移相器安装位置确定方法、系统、介质及计算设备 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130262001A1 (en) * | 2012-03-27 | 2013-10-03 | Hongbo Sun | State Estimation for Power System Using Hybrid Measurements |
CN103606113A (zh) * | 2013-12-02 | 2014-02-26 | 东北石油大学 | 基于pmu装置的电力系统静态状态估计方法 |
CN103840452A (zh) * | 2014-03-04 | 2014-06-04 | 国家电网公司 | 一种引入pmu量测信息的大电网状态估计方法 |
-
2017
- 2017-06-06 CN CN201710416590.XA patent/CN107169881A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130262001A1 (en) * | 2012-03-27 | 2013-10-03 | Hongbo Sun | State Estimation for Power System Using Hybrid Measurements |
CN103606113A (zh) * | 2013-12-02 | 2014-02-26 | 东北石油大学 | 基于pmu装置的电力系统静态状态估计方法 |
CN103840452A (zh) * | 2014-03-04 | 2014-06-04 | 国家电网公司 | 一种引入pmu量测信息的大电网状态估计方法 |
Non-Patent Citations (2)
Title |
---|
薛辉等: ""基于PMU量测数据和SCADA数据融合的电力系统状态估计方法"", 《电网技术》 * |
赵红噶等: ""计及PMU支路电流相量的状态估计模型"", 《电力系统自动化》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107665112A (zh) * | 2017-09-29 | 2018-02-06 | 北京紫优能源科技有限公司 | 一种电力控制流程的组态方法及装置 |
CN110543720A (zh) * | 2019-09-03 | 2019-12-06 | 北京交通大学 | 基于sdae-elm伪量测模型的状态估计方法 |
CN113659562A (zh) * | 2021-07-01 | 2021-11-16 | 国网江苏省电力有限公司电力科学研究院 | 可控移相器安装位置确定方法、系统、介质及计算设备 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100486073C (zh) | 考虑负荷电压特性的非线性动态状态估计算法 | |
CN103326358B (zh) | 基于同步相角测量装置的电力系统动态状态估计方法 | |
CN103972884A (zh) | 一种电力系统状态估计方法 | |
CN103840452B (zh) | 一种引入pmu量测信息的大电网状态估计方法 | |
CN106921156A (zh) | 一种基于多采样周期混合量测的主动配电网状态估计方法 | |
CN103399209B (zh) | 一种特高压双极直流输电线路工频参数测量方法 | |
CN101651343B (zh) | 基于混合动态仿真的电力系统模型校核方法及系统 | |
CN107843810A (zh) | 一种基于状态估计的有源配电网故障区段在线定位方法 | |
CN101505061B (zh) | 一种基于时域仿真的戴维南等值参数跟踪的计算方法 | |
CN106443246A (zh) | 基于pmu量测数据的小干扰稳定参数的在线辨识方法 | |
CN107577870A (zh) | 基于同步相量量测的配电网电压功率灵敏度鲁棒估计方法 | |
CN107169881A (zh) | 引入pmu支路电流相量的电力系统静态状态估计方法 | |
CN103886193B (zh) | 一种电力系统模糊自适应抗差估计方法 | |
CN103606113A (zh) | 基于pmu装置的电力系统静态状态估计方法 | |
CN108400592A (zh) | 一种计及潮流约束的配电网区间状态估计算法 | |
CN103869171A (zh) | 超高压同塔四回交流/双回双极直流线路零序参数测量法 | |
CN106501612B (zh) | 基于非均匀传输线的双回部分同塔线路零序参数测量方法 | |
CN105588984B (zh) | 一种混压双极直流输电线路零序参数精确测量方法 | |
CN103760400A (zh) | 基于广域量测系统的电压稳定在线监测方法 | |
Xue et al. | A dynamic state estimation method with PMU and SCADA measurement for power systems | |
CN106786536B (zh) | 考虑外网扩展Ward等值的状态估计方法 | |
CN105223436A (zh) | 一种同塔双回交流输电线路参数测量和计算方法 | |
CN104124684A (zh) | 基于广域量测信息的电力系统电压稳定性在线评估方法 | |
CN106526424B (zh) | 一种输电线路单相接地故障参数识别方法 | |
CN113702706A (zh) | 一种基于电力电子变换器的电网阻抗测量方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170915 |