高等电力系统稳态分析第三章电力系统状态估计
电力系统状态估计概述
电力系统状态估计研究综述摘要：电力系统状态估计是当代电力系统能量管理系统(EMS)的重要组成部分。
本文介绍了电力系统状态估计的概念、数学模型，阐述了状态估计的必要性及其作用，系统介绍了状态估计的研究现状，最后对状态估计的研究方向进行了展望。
关键词：电力系统；状态估计；能量管理系统0引言状态估计是当代电力系统能量管理系统(EMS)的重要组成部分，尤其在电力市场环境中发挥更重要的作用。
它是将可用的冗余信息(直接量测值及其他信息) 转变为电力系统当前状态估计值的实时计算机程序和算法。
准确的状态估计结果是进行后续工作(如安全分析、调度员潮流和最优潮流等)必不可少的基础。
随着电力市场的发展，状态估计的作用更显重要⑴o状态估计的理论研究促进了工程应用，而状态估计软件的工程应用也推动了状态估计理论的研究和发展。
迄今为止，这两方面都取得了大量成果。
然而，状态估计领域仍有不少问题未得到妥善解决，随着电力系统规模的不断扩大，电力工业管理体制向市场化迈进，对状态估计有了新要求，各种新技术和新理论不断涌现，为解决状态估计的某些问题提供了可能。
本文就电力系统状态估计的研究现状和进一步的研究方向进行了综合阐述。
1电力系统状态估计的概念1.1电力系统状态估计的基本定义状态估计也被称为滤波，它是利用实时量测系统的冗余度来提高数据精度，自动排除随机干扰所引起的错误信息，估计或预报系统的运行状态(或轨迹)o 状态估计作为近代计算机实时数据处理的手段，首先应用于宇宙飞船、卫星、导弹、潜艇和飞机的追踪、导航和控制中。
它主要使用了六十年代初期由卡尔曼、布西等人提出的一种递推式数字滤波方法，该方法既节约内存，又大大降低了每次估计的计算量[2,4]o电力系统状态估计的研究也是由卡尔曼滤波开始。
但根据电力系统的特点，即状态估计主要处理对象是某一时间断面上的高维空间(网络)问题，而且对量测误差的统计知识又不够清楚，因此便于采用基于统计学的估计方法如最小方差估计、极大验后估计、极大似然估计等方法，目前很多电力系统实际采用的状态估计算法是最小二乘法。
高等电力系统稳态分析 第三章 电力系统状态估计
二、电力系统状态估计－必要性
电力系统需要随时监视系统的运行状态 需要提供调度员所关心的所有数据 测量所有关心的量是不经济的，也是不 可能的，需要利用一些测量量来推算其 它电气量 由于误差的存在，直接测量的量不甚可 靠，甚至有坏数据
三、状态估计的作用
降低量测系统投资，少装测点 计算出未测量的电气量 利用量测系统的冗余信息，提高量测数 据的精度
对角元随测量量的增多而减小，亦即测量越多 时，估计越准确。 测量量的测量值与估计值的差，称为残差r， 表达式为：
ˆ Hx v Hx ˆ r zz
[I H(HT R 1H)1 HT R 1 ]v Wv
式中W称为残差灵敏度矩阵，表示残差与测量 误差之间的关系
一、最小二乘原理
令
J ( x) 0 2500x 3 3400x 5740 x x 1.36x 2.296 0 x 0.9852 x2,3 0.4926 j1.445
3
二、例题
状态的估计值x=0.9852 量测的估计值： 电流I=x=0. 9852 p.u.=0.9852A 电压U=Rx=0.9852p.u.=9. 852V 有功P=Rx2=0.9706p.u.=9.706W 量测的残差值： 电流残差νI=1.05-
由于通常测量误差的均值为零，所以估 计误差的均值为
ˆ ) (HTR 1H)1 HTR 1 E( v) 0 E (x x
在工程中往往以估计误差的协方差阵来 衡量状态量的估计值与真值间的差异， 估计误差的协方差阵为
T 1 1 T 1 T 1 T 1 T 1 1 T
T ˆ ˆ c E[(x x)(x x) ]
电力系统稳态分析-第三章
P jQ U
( R jX )
U
PR QX U
(电压降落纵分量）
U
PX QR （电压降落横分量） U
j0
0
（取 U
Ue
U
）
U 1
已知首端电压和首端功率时：
P1 R Q 1 X U1
(电压降落纵分量）
U 1
P1 X Q 1 R （电压降落横分量） U1
电压损耗：
U1 P1 R Q 1 X U1
I1
R jX
I2
U1
j
B 2
j
B 2
U2
中等长度线路的集中参 数等值电路
(2) 电力线路上的功率损耗
①阻抗支路上的功率损耗
P
P2 Q 2
2
2
U2
2
2
R
Q
ቤተ መጻሕፍቲ ባይዱ
P2 Q 2 U2
2
2
2
X
②导纳支路上的功率损耗
Py G U Q y 1 2
2
2
jBU
2
G
jB
2
电力系统稳态分析
第三章 简单电力网络的分析计算
第二节 电力线路和变压器运行状况的计算和分析
主讲 周任军
1、电力线路的分析计算
(1) 电力线路上的电压降落和电压损耗
U 1 P1 R Q 1 X U1 (电压降落纵分量）
P1 X Q 1 R U 1 （电压降落横分量） U1
1端流向 2 端； 2 端流向 1端；
即有功功率总是从电压相位超前的点流向电压相位滞后的点。
2、变压器运行情况分析
电力系统状态估计概述
电力系统状态估计研究综述摘要：电力系统状态估计是当代电力系统能量管理系统(EMS)的重要组成部分。
本文介绍了电力系统状态估计的概念、数学模型，阐述了状态估计的必要性及其作用，系统介绍了状态估计的研究现状，最后对状态估计的研究方向进行了展望。
关键词：电力系统；状态估计；能量管理系统0 引言状态估计是当代电力系统能量管理系统(EMS)的重要组成部分, 尤其在电力市场环境中发挥更重要的作用。
它是将可用的冗余信息(直接量测值及其他信息)转变为电力系统当前状态估计值的实时计算机程序和算法。
准确的状态估计结果是进行后续工作(如安全分析、调度员潮流和最优潮流等)必不可少的基础。
随着电力市场的发展，状态估计的作用更显重要[1]。
状态估计的理论研究促进了工程应用，而状态估计软件的工程应用也推动了状态估计理论的研究和发展。
迄今为止，这两方面都取得了大量成果。
然而，状态估计领域仍有不少问题未得到妥善解决，随着电力系统规模的不断扩大，电力工业管理体制向市场化迈进，对状态估计有了新要求，各种新技术和新理论不断涌现，为解决状态估计的某些问题提供了可能。
本文就电力系统状态估计的研究现状和进一步的研究方向进行了综合阐述。
1 电力系统状态估计的概念1.1电力系统状态估计的基本定义状态估计也被称为滤波，它是利用实时量测系统的冗余度来提高数据精度，自动排除随机干扰所引起的错误信息，估计或预报系统的运行状态(或轨迹)。
状态估计作为近代计算机实时数据处理的手段，首先应用于宇宙飞船、卫星、导弹、潜艇和飞机的追踪、导航和控制中。
它主要使用了六十年代初期由卡尔曼、布西等人提出的一种递推式数字滤波方法，该方法既节约内存，又大大降低了每次估计的计算量[2,4]。
电力系统状态估计的研究也是由卡尔曼滤波开始。
但根据电力系统的特点，即状态估计主要处理对象是某一时间断面上的高维空间(网络)问题，而且对量测误差的统计知识又不够清楚，因此便于采用基于统计学的估计方法如最小方差估计、极大验后估计、极大似然估计等方法，目前很多电力系统实际采用的状态估计算法是最小二乘法。
电力系统状态估计
电 力 工 程 系
North China Electric Power University
Department of Electrical Engineering
电力工程系
Department of Electrical Engineering
Thanks
Http\\ \ee
∆x = x − x0
H(x0 ) = ∂h(x) ∂x x=x0 ，H是m×n阶量测向量的雅可比矩阵
电 力 工 程 系
North China Electric Power University
Department of Electrical Engineering
二、状态估计的数学模型及算法
J (x) = [∆z − H(x0 )∆x]T R−1[∆z − H(x0 )∆x]
P =Vi2g −ViVj g cosθij −ViVjbsin θij ij
v是量测误差；z和v都是随机变量，是均值为零，方差为σ2 的正态分布随机向量 。 电 力 工 程 系
North China Electric Power University
Department of Electrical Engineering
∆z = z − h(x0 )
要使J(x)极小，应有 ： ∆ x = [HT (x0 )R−1H(x0 )]−1 HT (x0 )R−1∆z
^
迭代公式 ： x = [HT (x ∆
l+ ^ (l+1) ^ (l )
^
^ (l )
)R−1H(x )]−1 HT (x )R−1[z − h(x )]
^ (l )
North China Electric Power University
3 电力系统状态估计算法
量测雅克比矩阵 信息矩阵 H R H 状态估计误差方差阵
T 1
h( x ) H x
ˆ(k ) x x
T T 1 ˆ x x ˆ E x x H R H
1
量测估计误差
ˆ) H(x ˆ )x ˆ H(x ˆ )( x x ˆ) ˆ z h( x zz
收敛条件
ˆi ( l ) x
max i
x
ˆ (l ) ) J ( x ˆ ( l 1) ) J J(x ˆ (l ) x
三个收敛条件 任选其一即可
概述
WLS
FDSE
变换量测
比较
1
示例
ˆ (k ) ) z ( k ) z h( x
R 1 H H T R 1z ( k ) x ˆ ( k 1) x ˆ ( k ) x ˆ (k ) x
WLS与FDSE 求解方法的 区别？
l＝l+1
计算b；求解△ v( l ) v 否 是
Q-V迭代
vi( l )
max
v ?
KQ=1 v(l+1)=v(l)+△v(l) 否
KQ=0 KP=0? 是
概述
WLS
FDSE
变换量测
比较
示例
WLS与FDSE的区别
算法 求解方式 方程维数 系数矩阵 WLS 同时求解 θ 和 v n=na+nr 变化的 FDSE 分别求解 θ 和 v na和nr 常数
采用PQ分解法求解潮流 的思想，将有功和无功 解耦以及雅克比矩阵常 数化的方法用在加权最 小二乘法中，形成了快 速分解状态估计算法。
电力系统状态估计的原理
电力系统状态估计的原理
电力系统状态估计是指对电力系统的各个分量进行在线监测，并通过对监测数据的处理和分析，对电力系统的状态进行估计的技术。
电力系统状态估计的原理主要包括以下几个方面：
1.电力系统模型：电力系统状态估计需要建立电力系统的数学模型，包括线路参数、节点电压、母线注入功率等参数。
通常使用潮流方程来描述电力系统的运行情况。
2.测量数据：通过电力系统中的传感器和测量设备，获取电压、电流、功率、功角等各个分量的实时测量数据。
这些数据是电力系统状态估计的基础。
3.潮流方程求解：根据电力系统的模型和测量数据，可以建立潮流方程组，并利用数值方法求解潮流方程组，得到所有节点的电压、相角和功率等信息。
4.数据处理：将测量数据与潮流方程求解结果进行比对和匹配，通过误差最小化的方法，对电力系统状态进行修正和估计。
常用的方法有最小二乘法、卡尔曼滤波和最大似然估计等。
5.状态量调整：根据估计结果，对电力系统中的状态量进行调整。
比如，根据估计的电压值，调整变压器的调压装置，使得电压保持在合适的范围内。
6.结果评估：对估计结果进行评估，分析估计的准确性和可靠性。
如果发现估计结果与测量数据的差异较大，可能需要重新调整模型或校准测量设备。
综上所述，电力系统状态估计的原理主要是建立电力系统模型，获取实时测量数据，通过潮流方程求解和数据处理，对电力系统状态进行估计和调整，以实现对电力系统运行状态的实时监测和评估。
3节电力系统状态估计(WLS算法)
3节电系统状态估计报告【任务说明】:闭合的开关:打开的开关:打开的刀闸:线路:负荷G:发电机:母线:连接线(没有阻抗) Unit2Unit13节点系统主接线图任务：1、采用最小二乘状态估计算法,所有量测的权重都取1.0,编写状态估计程序(C/Matlab)。
2、按量测类型，列出量测方程（每一类写出一个方程）3、画出程序流程4、提交源程序，程序中每个函数的作用5、提交计算的输出结果(屏幕拷贝)系统参数：功率基值：100MW电压基值：230 kV线路阻抗参数（标么值）：线路量测(流出母线为正)：母线电压量测：负荷量测(流出母线为正)：发电量测(流入母线为正)：注：量测存在误差【数据预处理】首先根据基值将已知的量测值均转换为标幺值，并将功率值转换为流入量，得到如下数据：线路导纳参数（标么值）：线路注入功率量测(标幺值)：负荷点注入功率量测(标幺值)：发电机节点注入量测(流入母线为正)：发电机量测真值unit2 0.88-j0.0424 0.8892-j0.0424unit3 0.23+j0.24 0.2304+j0.2378母线电压量测（标幺值）：母线电压量测真值(幅值/角度)1 1.0087 1.0130/02 1.0198 1.0242/3.233 1.0281 1.0281/1.82【量测方程】选择节点1的电压相角为参考，为0度，以vi表示误差值。
1）节点1电压量测方程：Vi=Vi+v1即1.0087=V1+v12）1-3支路1号节点处注入有功功率功率：P ij=V i2g ij-V i V j(g ij cos+b ij sin)+v20.613=V12g13-V1V3(g13cos+b13sin)+v2即0.613=-1.6171V12-V1V3(-1.6171cos +13.698sin)+v2 3）1号节点注入功率：P i=V i2G ii +G ij cos+B ij sin+v3P1=V12G11+G1j cos+B1j sin+v3即-1.11=3.5613V12+V1V2(-1.9442cos -10.5107sin)+V1V3(-1.6171 cos -13.698 sin)+v3【流程图】【计算结果】其中iterations 为迭代次数，可见本例的迭代次数为4，收敛较快，状态估计得到的节点1、2、3电压分别为：234.0144444444444444444444444444444444444444444444【程序说明】遥测数据给定V 0，，k=0计算H(V (k),)和h(V (k),)A=H T R -1H, b=H T R -1(Z-h)求解A X=b,得Xk=k+1X (k+1)=X (k)+XNmax|X|<Y结束1、计算h矩阵的函数cal_hfunction h=cal_h(V,th0,B,G) %其中，V为节点电压估计值，th0为节点电压相角估计%值，B为节点电导矩阵，G为节点电纳矩阵b=-B; %线路电导矩阵g=-G; %线路电纳矩阵P=zeros(3,1); %初始化，节点注入功率Q=zeros(3,1);PP=zeros(3,3); %线路注入功率QQ=PP;th=[0;th0]; %节点1的电压相角为0for i=1:3P_P=0;Q_Q=0;for j=1:3if(j~=i)P_P=P_P+V(i)*V(j)*(G(i,j)*cos(th(i)-th(j))+B(i,j)*sin(th(i)-th(j)));Q_Q=Q_Q+V(i)*V(j)*(G(i,j)*sin(th(i)-th(j))-B(i,j)*cos(th(i)-th(j)));PP(i,j)=(V(i)^2)*g(i,j)-V(i)*V(j)*(g(i,j)*cos(th(i)-th(j))+b(i,j)*sin(th(i)-th(j)));QQ(i,j)=-(V(i)^2)*b(i,j)-V(i)*V(j)*(g(i,j)*sin(th(i)-th(j))-b(i,j)*cos(th(i)-th(j)));endendP(i)=(V(i)^2)*G(i,i)+P_P;Q(i)=-(V(i)^2)*B(i,i)+Q_Q;endVV=[V(1);V(2);V(3)];h=[P;Q;PP(1,2);PP(2,3);PP(3,1);QQ(1,2);QQ(2,3);QQ(3,1);PP(1,3);PP(2,1);P P(3,2);QQ(1,3);QQ(2,1);QQ(3,2);VV];2、计算H矩阵的函数cal_HHfunction H=cal_HH(V,th0,G,B,P,Q) %其中，P，Q为根据电压估计值计算得到的节点%注入电压b=-B;g=-G;PV=zeros(3,3); %节点注入功率对电压幅值的偏导数QV=zeros(3,3);Pth=zeros(3,3); %节点注入功率对电压相角的偏导数Qth=zeros(3,3);PPV=zeros(3,3); %P ij对V j的偏导数QQV=zeros(3,3); %Q ij对V j的偏导数PPth=zeros(3,3); %P ij对th j的偏导数QQth=zeros(3,3); %Q ij对th j的偏导数PPV1=zeros(3,3); %P ij对V i的偏导数QQV1=zeros(3,3); %Q ij对V i的偏导数PPth1=zeros(3,3); %P ij对th i的偏导数QQth1=zeros(3,3); %Q ij对th i的偏导数VV=eye(3);Vth=zeros(3,2);th=[0;th0];for i=1:3for j=1:3if (i~=j)PV(i,j)=V(i)*(G(i,j)*cos(th(i)-th(j))+B(i,j)*sin(th(i)-th(j)));QV(i,j)=V(i)*(G(i,j)*sin(th(i)-th(j))-B(i,j)*cos(th(i)-th(j)));PPV(i,j)=-V(i)*(g(i,j)*cos(th(i)-th(j))+b(i,j)*sin(th(i)-th(j)));QQV(i,j)=-V(i)*(g(i,j)*sin(th(i)-th(j))-b(i,j)*cos(th(i)-th(j)));PPV1(i,j)=2*V(i)*g(i,j)-V(j)*(g(i,j)*cos(th(i)-th(j))+b(i,j)*sin(th(i)-th(j)));QQV1(i,j)=-2*V(i)*b(i,j)-V(j)*(g(i,j)*sin(th(i)-th(j))-b(i,j)*cos(th(i)-th(j)));if (j~=1)Pth(i,j)=V(i)*V(j)*(G(i,j)*sin(th(i)-th(j))-B(i,j)*cos(th(i)-th(j)));Qth(i,j)=-V(i)*V(j)*(G(i,j)*cos(th(i)-th(j))+B(i,j)*sin(th(i)-th(j)));PPth(i,j)=-V(i)*V(j)*(g(i,j)*sin(th(i)-th(j))-b(i,j)*cos(th(i)-th(j)));QQth(i,j)=-V(i)*V(j)*(-g(i,j)*cos(th(i)-th(j))-b(i,j)*sin(th(i)-th(j)));endif(i~=1)PPth1(i,j)=V(i)*V(j)*(g(i,j)*sin(th(i)-th(j))-b(i,j)*cos(th(i)-th(j))); QQth1(i,j)=-V(i)*V(j)*(g(i,j)*cos(th(i)-th(j))+b(i,j)*sin(th(i)-th(j)));endelsePV(i,j)=(G(i,i)*(V(i)^2)+P(i))/V(i);QV(i,j)=(Q(i)-(V(i)^2)*B(i,i))/V(i);if (j~=1)Pth(i,j)=-B(i,i)*(V(i)^2)-Q(i);Qth(i,j)=P(i)-(V(i)^2)*G(i,i);endendendendH=[[PV,Pth(:,2:3)];[QV,Qth(:,2:3)];...[PPV1(1,2),PPV(1,2),0,PPth(1,2),0;...0,PPV1(2,3),PPV(2,3),PPth1(2,3),PPth(2,3);...PPV(3,1),0,PPV1(3,1),0,PPth1(3,1)];...[QQV1(1,2),QQV(1,2),0,QQth(1,2),0;...0,QQV1(2,3),QQV(2,3),QQth1(2,3),QQth(2,3);...QQV(3,1),0,QQV1(3,1),0,QQth1(3,1)];...[PPV1(1,3),0,PPV(1,3),0,PPth(1,3);...PPV(2,1),PPV1(2,1),0,PPth1(2,1),0;...0,PPV(3,2),PPV1(3,2),PPth(3,2),PPth1(3,2)];...[QQV1(1,3),0,QQV(1,3),0,QQth(1,3);...QQV(2,1),QQV1(2,1),0,QQth1(2,1),0;...0,QQV(3,2),QQV1(3,2),QQth(3,2),QQth1(3,2)];...[VV,Vth]];3、主程序calculate_all.m文件format longG=[3.5613,-1.9442,-1.6171;...-1.9442,3.0993,-1.1551;...-1.6171,-1.1551,2.7722]; %B为节点电导矩阵B=[-24.2087,10.5107,13.698;...10.5107,-20.295,9.7843;...13.698,9.7843,-23.4832]; %G为节点电纳矩阵P=[-1.11;0.88;0.23]; %节点注入功率量测值Q=[-0.135;-0.0424;0.24];PP=[0.613;-0.24;-0.459]; %线路1-2,2-3,3-1注入功率在首端的量测值QQ=[-0.012;0.066;-0.165];PP1=[0.467;-0.6;0.24]; %线路1-3,2-1,3-2注入功率在首端的量测值QQ1=[0.148;-0.024;-0.072];V=[1.0087;1.0198;1.0281]; %节点电压幅值量测值R=diag(ones(21,1)); %权重都取为1Z=[P;Q;PP;QQ;PP1;QQ1;V]; %量测值矩阵V0=[1;1;1]; %初值th0=[0;0];delta=100;iterations=0; %迭代次数while delta>0.000001iterations=iterations+1;h=cal_h(V0,th0,B,G); %计算h矩阵H=cal_HH(V0,th0,G,B,h(1:3,1),h(4:6,1)); %计算H矩阵A=H'*inv(R)*H;b=H'*inv(R)*(Z-h);d=A\b; %求解修正值delta=max(abs(d));V0=V0+d(1:3,1); %修正估计值th0=th0+d(4:5,1);enditerationsV0=V0*230; %转换为有名值th0=th0*180/pi; %转换为度for i=1:3j=num2str(i);v=num2str(V0(i));show1=strcat('The voltage magnitude of node ',j,' is', v,' kV');disp(show1);endfor i=1:2j=num2str(i+1);th=num2str(th0(i));show1=strcat('The phase angle of node ',j,' is ',th,' degrees');disp(show1);end。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
一、最小二乘原理
• 假设测量函数线性
n
hi(x) hijxj j1
(i1,2, ,m)
• 则状态量的值x与测量值 z间的关系为
zHxv
• 式中：H为m*n矩阵。 • 按最小二乘法建立目标函数
J(x)(zH)Tx (zH)x • 极值条件 J (x) 0
x
一、最小二乘原理
• 加权（提高精度）
JW (x)(zH)TW x(zH)x
– 测量误差的方差为
11
K
2 i
rii
c1zi
c2(F)2
一、测量方程
– 测量误差的方差阵
2 1
R
2 2
2 m
二、电力系统状态的可观察性
• 必要但非充分条件：雅可比矩阵的秩等 于n。 h(x) H(x) x xx0
• 有冗余度的目的是提高测量系统的可靠 性和提高状态估计的精确度。
第三章 电力系统状态估 计
State Estimation
参考书籍
• 《电力系统状态估计》于尔铿
第一节 概述
一、什么是状态估计
• 环境噪声使理想的运动方程无法精确求 解。
• 测量系统的随机误差，使测量向量不能 直接通过理想的测量方程求出状态真值。
• 通过统计学的方法加以处理以求出对状 态向量的估计值。这种方法，称为状态 估计。
• 动态估计与静态估计
二、电力系统状态估计－必要性
• 电力系统需要随时监视系统的运行状态 • 需要提供调度员所关心的所有数据 • 测量所有关心的量是不经济的，也是不
可能的，需要利用一些测量量来推算其 它电气量 • 由于误差的存在，直接测量的量不甚可 靠，甚至有坏数据
三、状态估计的作用
• 降低量测系统投资，少装测点
evi2
/
2
2 i
p(vi )
2ቤተ መጻሕፍቲ ባይዱ
2 i
• 方差越大表示误差大的概率增大
一、测量方程
– 用协方差表示不同时刻测量数据误差之间均 值的相关度
Ri vi(tk)vi(tkmt) m
• 通常 m0时，Ri 0 ；当 m0，Ri i2表
示不同时间的测量之间是不相关的，一般情况下， 不同测量的误差之间也是不相关的。
• 计算出未测量的电气量 • 利用量测系统的冗余信息，提高量测数
据的精度
– 独立测量量的数目与状态量数目之比，成为 冗余度。
四、状态估计的流程
结构信息 测量信息
修正输入
预过滤 假设模型 估计计算
检测 BD
识别
有 无 结束
四、状态估计与潮流计算的关系
• 潮流计算是状态估计的一个特例 • 状态估计用于处理实时数据，或者有冗
j i
j i
Q i f i ( G ie j j B if j j) e i ( G if j j B ie j j) ( 2 5 )
j i
j i
• 支路潮流
一、测量方程
S ij P i j jiQ j U i 2 y * i oU i(U * i U *j)y * ij (2 6 )
• 定义测量量向量为 z，待求的系统状态
量为 x，通过网络方程可以从估计的状
态量 xˆ ，求出估计的计算值 zˆ，如果测 量有误差，则计算值 zˆ与实际值 z之间 有误差 zzˆ，称为残差向量。
• 求出的状态量不可能使残差向量为零， 但可以得到一个使残差平方和为最小的 状态估计值。
第二节 电力系运行状态的表征与可观 察性
余的矛盾方程的场合 • 潮流计算用于无冗余矛盾方程的场合 • 两者的求解算法不同 • 在线应用中，潮流计算在状态估计的基
础上进行，也就是说，由状态估计提供 经过加工处理过的熟数据，作为潮流计 算的原始数据。
四、状态估计与潮流计算的关系
n节 点注 入量
潮流 计算
n节 点电 压
网络 参数
潮流计算
m维 测量 量
2
除1外加上所有节点的电压 同1
模值
3N-2
3
每条支路两侧的有功、无功 式(2-6)、式(2-7) 4M
潮流
4
除3外，再加所有节点的电 式(2-6)、式(2- 4M+N
压模值
7)、式(2-9)
5
完全的测量系统
式(2-4)~式(2-7)、4M+3N 式(2-9)
一、测量方程
• 节点注入功率方程式
P i e i ( G ie j j B if j j) f i ( G if j j B ie j j) ( 2 4 )
（在状态估计中，平衡节点的电压模值 也是测量值，需要当作状态量，只有平 衡节点电压相角可以确定）
一、测量方程
• 五种基本测量方式(N为节点数、M为支 路数)
测量方式 Z的分量
方程式h(x)
z的维数
1
平衡节点电压模值除平衡节 式(2-4)、式(2- 2N-1
点外所有节点的注入功率 5)、式(2-9)
• W为一适当选择的加权正定阵
• 假设W=R-1，R为测量误差方差阵
• 于是目标函数可以写成
J(x)(zH)Tx R 1(zH)x
S ji P ji jj Q i U 2 jy * j oU j(U *j U * i)y * ij (2 7 )
• 电压实部、虚部和模值、相角的关系
i
arctfai n ei
(28)
uiei2fi2 (29)
一、测量方程
• 数学模型不完整
• 测量系统的系统误差
• 随机误差
– 随机误差的概率密度函数
一、测量方程
• 测量矢量：z=[z1,z2,…,zm]T, m维 • 测量误差矢量： ν=[ν1, ν2,…, νm]T, m维 • 测量函数：h(x)=[h1(x),h2(x),…,hm(x)]T • 状态量：x＝ [x 1, x 2,…, x n]T, n维 • 对于N节点的系统，状态量数目为n=2N-1
测量 噪声
估计 算法
n节 点电 压
网络 参数
状态估计
四、状态估计与潮流计算的关系
Vi,Pi,Qi,Ii Pij,Qij,Iij
状态估计
模拟操作： 开关操作 出力调整 负荷调整 分接头调整
Vi,Pi,Qi 潮流计算
四、状态估计基本思路
• 电力系统的测量量一般包括支路功率、 节点注入功率、节点电压模值等；状态 变量是各节点的电压模值和相角。
• 保证可观性是测量点布置的最低要求。
三、坏数据的可检测和可辨识性
• 可检测：可以判断系统中是否有坏数据 • 可辨识：若有坏数据，可以找出谁是坏
数据 • 量测冗余度越大，坏数据的可检测和可
辨识性越好 • 例：
一杆秤称重，不可检测、不可辨识 两杆秤称重，可检测、不可辨识 三杆秤称重，可检测、可辨识
第三节 最小二乘估计