CN105303471A - 电力系统状态的双曲余弦型最大指数平方抗差估计方法 - Google Patents
电力系统状态的双曲余弦型最大指数平方抗差估计方法 Download PDFInfo
- Publication number
- CN105303471A CN105303471A CN201510850198.7A CN201510850198A CN105303471A CN 105303471 A CN105303471 A CN 105303471A CN 201510850198 A CN201510850198 A CN 201510850198A CN 105303471 A CN105303471 A CN 105303471A
- Authority
- CN
- China
- Prior art keywords
- sigma
- exp
- cosh
- node
- rsqb
- 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.)
- Granted
Links
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了属于电力系统调度自动化领域的电力系统状态的双曲余弦型最大指数平方抗差估计方法。该方法是基于抗差性好、计算效率高的双曲余弦型最大指数平方抗差状态估计;包括如下步骤:步骤A.提取电力系统节点注入有功和无功、支路有功和无功功率,以及节点电压幅值参数;以此建立双曲余弦型最大指数平方抗差状态估计模型;步骤B.利用原-对偶内点算法,对所述双曲余弦型最大指数平方抗差状态估计模型求解。本发明有效抑制包括一致性不良数据在内的多个不良数据,显示了良好的抗差性,并有很高的计算效率,得到广泛应用。
Description
技术领域
本发明属于电力系统调度自动化领域。特别涉及一种电力系统状态的双曲余弦型最大指数平方抗差估计方法。
背景技术
电力系统状态估计是能量管理系统的基础和核心。现在几乎每一个大型调度中心都安装了状态估计器,状态估计已成为电网安全运行的基石。自1970国外学者首次提出状态估计以来,人们对状态估计的研究和应用已经有40多年的历史了,这期间涌现出了各种各样的状态估计方法。
目前,在国内外应用最为广泛的状态估计是加权最小二乘法(Weightedleastsquares,WLS)。WLS模型简洁,求解容易,但是其抗差性很差。为了增强抗差性,一般有两种方法。第一种是在WLS估计之后加入不良数据辨识环节,例如最大正则化残差检验法(LNR)或估计辨识方法等;另一种是采用抗差状态估计方法。目前,国内外学者已经提出的抗差状态估计方法(Robuststateestimation)包括加权最小绝对值估计(Weightedleastabsolutevalue,WLAV)、非二次准则法(QL、QC等)以合格率最大为目标的状态估计(Maximumnormalmeasurementrate,MNMR)等。但是这些抗差状态估计方法的估计性能仍有待提高。
发明内容
本发明的目的是提出一种电力系统状态的双曲余弦型最大指数平方抗差估计方法,其特征在于,该方法是基于抗差性好、计算效率高的双曲余弦型最大指数平方抗差状态估计;包括如下步骤:
步骤A.提取电力系统节点注入有功和无功、支路有功和无功功率,以及节点电压幅值参数;以此建立双曲余弦型最大指数平方抗差状态估计模型;
步骤B.利用原-对偶内点算法,对所述双曲余弦型最大指数平方抗差状态估计模型求解。
所述步骤A建立双曲余弦型最大指数平方抗差状态估计模型为: s.t.g(x)=0,r=z-h(x),其中:z∈Rm为量测矢量,包括节点注入有功和无功、支路有功和无功以及节点电压幅值量测;x∈Rn为状态矢量,包括节点电压幅值和平衡节点除外的其他各个节点相角;h:Rn→Rm为由状态矢量到量测矢量的非线性映射;ri为残差矢量r的第i个元素;g(x):Rn→Rc为零注入功率等式约束;wi为第i个量测量的权重,σ0和σ1为窗宽参数。
所述步骤B包括:
步骤B1:令x为平启动状态变量;选择r(0)=λ(0)=π(0)=0;设置收敛判据ε=10-6,置迭代计数器k=0;
步骤B2:求解修正方程,得到[dxTdrTdλTdπT];
步骤B3,修正变量
步骤B4:判断是否收敛,若max(dx)<ε,则转步骤B6,否则进入步骤B5;
步骤B5:令迭代计数器k=k+1,进入步骤B2;
步骤B6:输出最优解,结束。
所述步骤B2包括:
步骤B21:形成量测方程以及零注入功率约束对应的雅克比矩阵及形成量测方程以及零注入功率约束对应的海森矩阵▽2h(x)及▽2g(x);
其中h(x)为状态矢量到量测矢量的映射,即为量测估计值;z=h(x)+ε为量测方程,z和x如步骤A所示,ε∈Rm为量测误差向量;h(x)的具体表达式介绍如下:对于节点i的电压幅值量测,vi=vi;对于注入有功量测,j∈Ni意思为与节点i相联的所有节点;Gij,Bij分别为节点导纳矩阵第i行第j列的实部和虚步,vi和vj分别为节点i和节点j的电压幅值;θij为节点i和节点j的相角差;对于注入无功量测,对于支路ij的有功量测,其中gsi为支路ij的首端对地电导,其中gij和bij分别为支路ij的串联电导和串联电纳;对于支路ij的无功量测, 其中bsi为支路ij的首端对地电纳。
g(x)=0为零注入功率约束;
步骤B22:引入拉格朗日函数
式中:λ∈Rc及π∈Rm为拉格朗日乘子矢量;计算
式中,W为Rm×m的对角阵,其对角元素为
步骤B23:求解方程 得到[dxTdrTdλTdπT]。
本发明的有益效果对电力系统状态采用双曲余弦型最大指数平方抗差状态估计,在估计过程中可有效抑制包括一致性不良数据在内的多个不良数据,显示了良好的抗差性,并具有很高的计算效率,非常适宜于实际工程应用。
具体实施方式
本发明提出一种电力系统状态的双曲余弦型最大指数平方抗差状态估计方法,下面结合实施例详细描述本发明。
本发明是基于抗差性好、计算效率高的双曲余弦型最大指数平方抗差状态估计(Hyperboliccosinemaximumexponentialsquarestateestimation,COSH-MES)方法,包括下列步骤:
步骤A:提取电力系统节点注入有功和无功、支路有功和无功功率,以及节点电压幅值参数;以此建立双曲余弦型最大指数平方抗差状态估计模型;具体地,本发明提出的COSH-MES的模型如下所示
s.t.g(x)=0(2)
r=z-h(x)(3)
式中:z∈Rm为量测矢量,常包括节点注入有功和无功、支路有功和无功以及节点电压幅值量测等;x∈Rn为包括节点电压幅值和相角的状态矢量(平衡节点相角除外);h:Rn→Rm为由状态矢量到量测矢量的非线性映射;ri是残差矢量r的第i个元素;g(x):Rn→Rc为零注入功率等式约束;wi为第i个量测量的权重,σ0和σ1为窗宽参数。
步骤B:利用原-对偶内点算法,对所述双曲余弦型最大指数平方抗差状态估计模型求解。
(1)COSH-MES模型的求解方法
注意到COSH-MES模型(1)~(3)是一个含有等式约束和不等式约束的最优化问题,适宜用原-对偶内点算法进行求解。为使本领域技术人员更好地理解本发明,首先给出详细的推导过程如下:
引入拉格朗日函数
式中:λ∈Rc及π∈Rm为拉格朗日乘子矢量。
为取得最优值,根据KKT条件,可得
式中:
以上方程由牛顿法求解可得
[▽2g(x)λ-▽2h(x)π]dx+GTdλ-HTdπ=-Lx(9)
Gdx=-Lλ(10)
-Hdx-dr=-Lπ(11)
式(12)可表示为矩阵形式,为
Wdr+dπ=-Lr(13)
式中:W为Rm×m的对角阵,其对角元素为
根据式(9)、(10)、(11)及(13),可得修正方程为
求解式(14)可得[dxTdrTdλTdπT],则迭代即可持续进行。
(2)COSH-MES模型的求解步骤
在介绍COSH-MES模型的求解推导过程之后,发明人将求解步骤归纳如下:
步骤B1:进行初始化,令x为平启动状态变量;选择r(0)=λ(0)=π(0)=0;设置收敛判据ε=10-6,置迭代计数器k=0。
具体地,令x(0)∈Rn代表由所有节点电压幅值和相角组成的的平启动状态变量(参考节点相角除外);选择r(0)=λ(0)=π(0)=0,其中λ∈Rc及π∈Rm为拉格朗日乘子矢量,m为量测量的个数,而c为零注入功率约束的个数;置迭代计数器k=0。
步骤B2:求解式(14)修正方程,以完成对变量的修正,得到[dxTdrTdλTdπT]。
步骤B3:修正变量为:
步骤B4:判断是否收敛,若max(dx)<ε,则转步骤B6,否则进入步骤B5;
步骤B5:令迭代计数器k=k+1,进入步骤B2;以及
步骤B6:输出最优解,结束。
实施例
设定利用IEEE标准系统检验基于原-对偶内点算法的COSH-MES的性能。试验采用全量测,量测值通过在潮流计算的结果上叠加白噪声(均值为0,标准差为τ)来获得。对于电压量测,取τV=0.005p.u.;对于功率量测,取τPQ=1MW/MVar。测试环境为PC机,CPU为Intel(R)Core(TM)i3M370、主频为2.40GHz、内存2.00GB。
1.抗差性能的比较
发明人将本发明的COSH-MES与其他状态估计器进行比较,来测试COSH-MES的抗差性。
在IEEE-14系统上设置4个一致性不良数据(P1-2、Q1-2、P1、Q1)。所设置的不良量测值以及量测量的正确值如表1所示。
表1COSH-MES对IEEE14系统一致性不良数据的辨识
作为对比,首先用广为应用的WLS进行估计,并用LNR进行不良数据的辨识(简记为WLS+LNR)。首次辨识的结果为:10个量测量的标准化残差大于门槛值(3.0),这10个量测量被认为是可疑数据;其中标准化残差最大的量测量为P2-1,删去该量测后重新运行WLS;此时发现P2的标准化残差最大。以上过程循环4次,4个良好的量测量被LNR误认为是可疑数据而被删去,但真正的不良数据仍然存在。可见,WLS+LNR不能辨识一致性不良数据。
应用COSH-MES方法的估计结果如表1所示。可以发现,即使量测量中存在一致性不良数据,COSH-MES的估计值与真值也可很好地吻合。在IEEE其他系统的多次试验也表明COSH-MES在估计的过程中可以自动抑制不良数据,具有良好的抗差性。
2.计算效率的比较
发明人为了进行效率比较,在正常量测条件下分别对四种状态估计器WLS、WLAV、MNMR以及COSH-MES进行了测试,其中后三种属于抗差状态估计器。在试验中,WLS采用牛顿法求解,其他三种状态估计采用内点法求解;且MNMR采用两阶段法,即第一阶段进行WLS估计,第二阶段将WLS的估计值作为MNMR估计的初值进行计算。
共进行50次仿真试验,状态估计收敛时的迭代次数以及平均计算耗时如表2所示。由表2可见,在这四种状态估计器中,WLS的计算效率最高;而在后三种抗差状态估计器中,COSH-MES的计算效率最高;而且随着系统规模的增大,COSH-MES的迭代次数以及计算耗时增长的很缓慢,因而COSH-MES适用于实际的大规模系统的估计。
表2四种状态估计器的迭代次数以及计算耗时
综上所述,本发明提出的COSH-MES在估计过程中可有效抑制包括一致性不良数据在内的多个不良数据,显示了良好的抗差性,并具有很高的计算效率,非常适宜于实际工程应用。
Claims (4)
1.一种电力系统状态的双曲余弦型最大指数平方抗差估计方法,其特征在于,该方法是基于抗差性好、计算效率高的双曲余弦型最大指数平方抗差状态估计;包括如下步骤:
步骤A.提取电力系统节点注入有功和无功、支路有功和无功功率,以及节点电压幅值参数;以此建立双曲余弦型最大指数平方抗差状态估计模型;
步骤B.利用原-对偶内点算法,对所述双曲余弦型最大指数平方抗差状态估计模型求解。
2.根据权利要求1所述电力系统状态的双曲余弦型最大指数平方抗差估计方法,其特征在于,所述步骤A建立双曲余弦型最大指数平方抗差状态估计模型为: s.t.g(x)=0,r=z-h(x),其中:z∈Rm为量测矢量,包括节点注入有功和无功、支路有功和无功以及节点电压幅值量测;x∈Rn为状态矢量,包括节点电压幅值和平衡节点除外的其他各个节点相角;h:Rn→Rm为由状态矢量到量测矢量的非线性映射;ri为残差矢量r的第i个元素;g(x):Rn→Rc为零注入功率等式约束;wi为第i个量测量的权重,σ0和σ1为窗宽参数。
3.根据权利要求1所述电力系统状态的双曲余弦型最大指数平方抗差估计方法,其特征在于,所述步骤B包括:
步骤B1:令x为平启动状态变量;选择r(0)=λ(0)=π(0)=0;设置收敛判据ε=10-6,置迭代计数器k=0;
步骤B2:求解修正方程,得到[dxTdrTdλTdπT];
步骤B3,修正变量
步骤B4:判断是否收敛,若max(dx)<ε,则转步骤B6,否则进入步骤B5;
步骤B5:令迭代计数器k=k+1,进入步骤B2;
步骤B6:输出最优解,结束。
4.根据权利要求1所述电力系统状态的双曲余弦型最大指数平方抗差估计方法,其特征在于,所述步骤B2包括:
步骤B21:形成量测方程以及零注入功率约束对应的雅克比矩阵及形成量测方程以及零注入功率约束对应的海森矩阵及
其中h(x)为状态矢量到量测矢量的映射,即为量测估计值;z=h(x)+ε为量测方程,z和x如步骤A所示,ε∈Rm为量测误差向量;h(x)的具体表达式介绍如下:对于节点i的电压幅值量测,vi=vi;对于注入有功量测,j∈Ni意思为与节点i相联的所有节点;Gij,Bij分别为节点导纳矩阵第i行第j列的实部和虚步,vi和vj分别为节点i和节点j的电压幅值;θij为节点i和节点j的相角差;对于注入无功量测,对于支路ij的有功量测,其中gsi为支路ij的首端对地电导,其中gij和bij分别为支路ij的串联电导和串联电纳;对于支路ij的无功量测,其中bsi为支路ij的首端对地电纳,
g(x)=0为零注入功率约束;
步骤B22:引入拉格朗日函数
式中:λ∈Rc及π∈Rm为拉格朗日乘子矢量;计算
式中,W为Rm×m的对角阵,其对角元素为
步骤B23:求解方程 得到[dxTdrTdλTdπT]。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510850198.7A CN105303471B (zh) | 2015-11-27 | 2015-11-27 | 电力系统状态的双曲余弦型最大指数平方抗差估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510850198.7A CN105303471B (zh) | 2015-11-27 | 2015-11-27 | 电力系统状态的双曲余弦型最大指数平方抗差估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105303471A true CN105303471A (zh) | 2016-02-03 |
CN105303471B CN105303471B (zh) | 2019-05-03 |
Family
ID=55200701
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510850198.7A Active CN105303471B (zh) | 2015-11-27 | 2015-11-27 | 电力系统状态的双曲余弦型最大指数平方抗差估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105303471B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107516164A (zh) * | 2017-08-11 | 2017-12-26 | 华北电力大学 | 一种反余切导数型抗差状态估计方法 |
CN111507591A (zh) * | 2020-04-07 | 2020-08-07 | 山东科技大学 | 电力系统状态确定方法、装置、计算机介质及存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003143756A (ja) * | 2001-10-30 | 2003-05-16 | Toshiba Corp | 電力系統の縮約方法及び処理プログラム記録媒体 |
CN102801162A (zh) * | 2012-08-23 | 2012-11-28 | 清华大学 | 一种两阶段线性加权最小二乘电力系统状态估计方法 |
CN102868157A (zh) * | 2012-09-11 | 2013-01-09 | 清华大学 | 一种基于最大指数绝对值目标函数的抗差状态估计方法 |
-
2015
- 2015-11-27 CN CN201510850198.7A patent/CN105303471B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003143756A (ja) * | 2001-10-30 | 2003-05-16 | Toshiba Corp | 電力系統の縮約方法及び処理プログラム記録媒体 |
CN102801162A (zh) * | 2012-08-23 | 2012-11-28 | 清华大学 | 一种两阶段线性加权最小二乘电力系统状态估计方法 |
CN102868157A (zh) * | 2012-09-11 | 2013-01-09 | 清华大学 | 一种基于最大指数绝对值目标函数的抗差状态估计方法 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107516164A (zh) * | 2017-08-11 | 2017-12-26 | 华北电力大学 | 一种反余切导数型抗差状态估计方法 |
CN111507591A (zh) * | 2020-04-07 | 2020-08-07 | 山东科技大学 | 电力系统状态确定方法、装置、计算机介质及存储介质 |
CN111507591B (zh) * | 2020-04-07 | 2021-03-19 | 山东科技大学 | 电力系统状态确定方法、装置、计算机介质及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN105303471B (zh) | 2019-05-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102868157B (zh) | 一种基于最大指数绝对值目标函数的抗差状态估计方法 | |
CN102801162B (zh) | 一种两阶段线性加权最小二乘电力系统状态估计方法 | |
Schenato et al. | Bayesian linear state estimation using smart meters and pmus measurements in distribution grids | |
CN101599643B (zh) | 一种基于指数型目标函数的电力系统抗差状态估计方法 | |
CN107016489A (zh) | 一种电力系统抗差状态估计方法和装置 | |
CN107038292A (zh) | 一种基于自适应多变量非参数核密度估计的多风电场出力相关性建模方法 | |
CN105514977A (zh) | 一种电力系统状态的双曲余弦型抗差状态估计方法 | |
CN105512502B (zh) | 一种基于残差归一化的权函数最小二乘状态估计方法 | |
CN103413053B (zh) | 一种基于内点法的电力系统抗差状态估计方法 | |
CN103886193B (zh) | 一种电力系统模糊自适应抗差估计方法 | |
CN104836223A (zh) | 电网参数错误与不良数据协同辨识与估计方法 | |
CN105514978B (zh) | 一种混合整数二次规划形式的抗差状态估计方法 | |
CN105490269A (zh) | 基于wams量测的多区域电力系统状态估计方法及系统 | |
CN105305440B (zh) | 电力系统状态的双曲余弦型最大指数绝对值抗差估计方法 | |
CN105303471A (zh) | 电力系统状态的双曲余弦型最大指数平方抗差估计方法 | |
CN102427227B (zh) | 考虑零注入约束的修正快速解耦电力系统状态估计方法 | |
CN104239716A (zh) | 一种基于参数偏差灵敏度的电网设备参数识别与估计方法 | |
CN104393588A (zh) | 一种基于最小均方误差的最优pmu布置方法 | |
CN106372440A (zh) | 一种并行计算的配电网自适应抗差状态估计方法及装置 | |
CN109255541B (zh) | 一种基于最小二乘和一乘之和的配电网抗差状态估计方法 | |
CN104252571B (zh) | 基于多预测‑校正内点法的wlav抗差状态估计方法 | |
CN103838962B (zh) | 一种计及pmu量测的分步线性状态估计方法 | |
CN106159941B (zh) | 一种考虑实际量测误差传递特性的电力系统状态估计方法 | |
CN106022968A (zh) | 一种双曲正切型抗差状态估计方法 | |
CN103593565B (zh) | 结构风险最小化的加权最小二乘电力系统状态估计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |