八、电力系统状态估计.pptVIP

?版权所有 于是快速解耦法的雅可比矩阵为 ?版权所有 再用Ri－1计算出信息矩阵 式中：A和B可由线路参数和加权值求出，其表示式为 有功测量向量与无功测量向量写成标幺值后分别为 ?版权所有 于是迭代方程式为 ?版权所有 从初值电压开始的迭代结果如下： 迭代序号l θ2（rad） θ3（rad） U2 U3 1 0.001490 -0.001600 1.048792 1.052866 2 0.001388 -0.001431 1.048850 1.052718 3 0.001354 -0.001449 1.048847 1.052726 最后，按状态量的估计值，计算出三个节点注入功率估计值的标幺值如下 P1＋jQ1＝0.11051－j0.24572 P2＋jQ2＝0.20439－j0.24346 P3＋jQ3＝－0.3149－j0.4925 ?版权所有 结 束 待续… ?版权所有 在电力系统中，一般取权重为各量测量方差的倒数，即 ，这样 最后达到 其中 代表状态量x的估计值 ?版权所有 对上面的加权最小二乘法，写成矩阵形式，得状态估计的目标函数： 即在给定量测向量z之后，状态估计量 是使目标函数 达到最小的x值。 式中R是以 为对角元素的m?m阶量测误差方差阵。表示量测权重，式的含意即是使量测量加权残差平方和为最小。 ?版权所有 （三）基本加权最小二乘法状态估计 加权最小二乘法状态估计的目标函数： 由于h(x)为x的非线性函数，无法直接计算 ，需要用迭代的方法求解。 先假定状态量初值为x(0)，使h(x)在x(0)处线性化，并用泰勒级数在x(0)附近展开h(x)，并略去二阶以上项： h(x)=h(x(0))+H(x(0))Δx 式中： Δx= x-x(0)，H(x(0))是函数向量h(x)的雅可比矩阵，其元素为 ?版权所有 取Δz=z- h(x(0))，展开J(x) ，得 上式中第一项与Δx无关，因此，要使目标函数最小，第二项应为0，从而有： 展开 其中： ?版权所有 只有当x(0)充分接近 时泰勒级数略去高数项后才能是足够近似的。应用上式作逐次迭代，可以得到 。若以(l)表示迭代序号，上面两式可以写成： 由此得到： 按上两式进行迭代修正，直到目标函数 接近于最小为止 （3） ?版权所有 收敛判据可以是下三项中任意一项： ?版权所有 经过l次迭代满足收敛标准时，求得 ，即为最优状态估计值 。此时量测量的估计值是 几个概念 状态估计的误差为 ，可得 测量误差：v = z-h(x) 残差：量测量与量测估计值之差 。 状态估计误差方差阵： ?版权所有 状态估计误差方差阵： 其中， 由于真值x是未知的，近似用 代替估计误差方差阵中的x，有 称HTR-1H 为信息矩阵（gain matrix） ?版权所有 加权最小二乘法估计步骤 从状态量的初值计算测量函数向量h(x(0))和雅可比矩阵H(x(0))。 由测量z和h(x(0))计算残差z-h(x(l))和目标函数J (x(l)),并用雅可比矩阵H(x(l))计算信息矩阵［HTR-1H］和向量HTR-1[z-h(x(l))]。 解方程求取状态修正量Δx(l)，并取其中绝对值最大值max|Δxi(l)| 检查是否达到收敛标准 若未达到收敛标准，修改状态量x(l+1)=x(l)+Δx(l)，继续迭代计算，直到收敛为止。 将计算结果送入不良数据检测于辨识入口 ?版权所有 ?版权所有 （四）关于H矩阵 HTR-1H一般为稀疏矩阵，所以可用稀疏矩阵技巧进行求解。由前述可得 或写成 A阵是n?n的对称稀疏矩阵，它的结构与导纳矩阵不一样，是取决于网络结构与测点的布置。 对线路，不论在线路哪一侧，也不论是有功或无功，只要有一个测量就能出现aij元素 对节点i的有功或无功注入的测量值，不仅与节点i的状态量有关，而且还与同节点i有直接连接的相邻节点的状态量有关。 节点i的电压测量值仅在H阵i列有非零元素，在A阵中也只影响相应的i行对角元 ?版权所有 对于图2-5所示的例子，在H阵中，相应于节点i注入测量的行（设为m行）的i列以及与i相关的各节点（如i、j、k）的列均为非零元素，即hme、hmi、hmj、hmk为非零元素，即相应的H阵为 可以看出，相应这一测量值，在A阵（下三角）中将使aie、aje、aj