一种基于分解滤波的网络等值算法制造技术

技术编号:13835200 阅读:99 留言:0更新日期:2016-10-15 15:09
本发明专利技术公开了一种基于分解滤波的网络等值算法,用两组第一次测量的电压向量观测数据、电流向量观测数据求取递推参数估计的初值,在递推最小二乘算法中加入遗忘因子,对矩阵P进行UD分解,将附加遗忘因子的递推最小二乘算法公式中的P(k)进行UD分解,并在递推过程中实时修正U(k)、D(k),而不是修正P(k),从而能保证P(k)的非负定性。应用附加遗忘因子的UD分解滤波算法进行戴维南等值参数估计,从而实现系统的网络等值。本发明专利技术的有益效果是将递推算法中的协方差矩阵P进行UD分解,避免了递推过程中对P的修正,有效地保证了其非负定性,提高了参数估计的数值稳定性。

【技术实现步骤摘要】

本专利技术属于电力
,涉及一种基于分解滤波的网络等值算法
技术介绍
随着电力系统规模的不断扩大,现在电力系统常常运行在稳定的边界。近年来国内外发生了多起电力系统电压崩溃事故,电力系统电压问题引起了国内外学者的普遍关注。与此同时,同步相量测量装置(PMU)在电力系统中的应用日益广泛;它能够高精度地直接测量节点电压相量和相关支路电流相量,使得根据系统实时运行状态监视其电压稳定情况成为可能。基于PMU的测量数据估计系统等值模型参数是实现电力系统电压稳定实时监测的关键技术。假设任一时刻从系统中的某一条负荷母线看进去,可以将整个系统等值为一个由该母线和其余外部系统的戴维南等值电源母线组成的2节点系统,如图1所示;其中,外部系统的戴维南等值电势和等值阻抗Zeq阻抗需要根据负荷侧母线电压、电流相量的测量值进行实时估计。目前,应用较为广泛的参数估计方法是递推最小二乘算法。上述线性方程组中未知参数相量θ的最小二乘估计的递推算法公式为: θ ^ ( k ) = θ ^ ( k - 1 ) + G ( k ) [ y ( k ) - φ T ( k ) θ ^ ( k - 1 ) ] G ( k ) = P ( k ) φ ( k ) = P ( k - 1 ) φ ( k ) [ I + φ T ( k ) P ( k - 1 ) φ ( k ) ] - 1 P ( k ) = [ I - G ( k ) φ T ( k ) ] P ( k - 1 ) ]]>式中,k是迭代次数,G(k)是第k次迭代的估计误差校正系数矩阵,P(k)是误差的协方差矩阵。可以证明本文档来自技高网...

【技术保护点】
一种基于分解滤波的网络等值算法,其特征在于按照以下步骤进行:步骤1:用两组第一次测量的电压向量观测数据、电流向量观测数据求取递推参数估计的初值U(0)和D(0);θ^(0)=[Φ0W0Φ0T]-1Φ0W0Z0]]>U(0)D(0)UT(0)=P(0)=[Φ0W0Φ0T]-1]]>式中,对矩阵P(0)进行UD分解即可得到U(0)和D(0);其中,和是输入的前2组电流向量观测数据;和是输入的前2组电压向量观测数据;j为虚数单位,T上标表示矩阵转置,ρ称为遗忘因子;步骤2:用当前的电压相量电流相量得到φT(k)和y(k);φT(k)=10-Ir(k)Ii(k)01-Ii(k)-Ir(k),]]>y(k)=[Vr(k) Vi(k)]T;其中,k是递推次数,和为第k次递推时输入的电压和电流向量,Vr(k)和Vi(k)分别为输入的电压相量的实部和虚部,Ir(k)和Ii(k)分别为输入的电流相量的实部和虚部;步骤3:求取f(k)和g(k):f(k)=UT(k‑1)φ(k),g(k)=D(k‑1)f(k);步骤4:求取β(k):β(k)=ρI+fT(k)g(k)其中I是单位矩阵;步骤5:求取G(k):G(k)=U(k‑1)g(k)β‑1(k)步骤6:进行戴维南等值参数的估计:θ^(k)=θ^(k-1)+G(k)[y(k)-φT(k)θ^(k-1)]]]>式中,y(k)=[Vr(k) Vi(k)]T;步骤7:设将其进行UD分解,得到H(k)和步骤8:求取U(k)和D(k),计算公式如下:U(k)=U(k-1)H(k),D(k)=ρ-1D‾(k)]]>步骤9:返回步骤2,直到没有新的测量数据为止。...

【技术特征摘要】
1.一种基于分解滤波的网络等值算法,其特征在于按照以下步骤进行:步骤1:用两组第一次测量的电压向量观测数据、电流向量观测数据求取递推参数估计的初值U(0)和D(0); θ ^ ( 0 ) = [ Φ 0 W 0 Φ 0 T ] - 1 Φ 0 W 0 Z 0 ]]> U ( 0 ) D ( 0 ) U T ( 0 ) = P ( 0 ) = [ Φ 0 W 0 Φ 0 T ] - 1 ]]>式中,对矩阵P(0)进行UD分解即可得到U(0)和D(0);其中,和是输入的前2组电流向量观测数据;和是输入的前2组电压向量观测数据;j为虚数单位,T上标表示矩阵转置,ρ称为遗忘因子;步骤2:用当前的电压相量电流相量得到φT(k)和y(k); φ T ( k ) = 1 0 - I r ( k ) I i ( k ) 0 1 - I i ( k ) - I r ( k ) , ]]>y(k)=[Vr(k) Vi(k)]T;其中,k是递推次数,和为第k次递推时输入的电压和电流向量,Vr(k)...

【专利技术属性】
技术研发人员:王颖陈亮邵华张章张倩茅习朋荆志朋齐晓光张丽洁康伟任志刚李树水贺春光凌云鹏安佳坤高立坡朱士加赵建华申永鹏翟广心张欣悦赵阳胡诗尧唐帅孟繁华
申请(专利权)人:国家电网公司国网河北省电力公司经济技术研究院
类型:发明
国别省市:北京;11

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1