CN107947160A - 一种基于矩阵变换的双线性wlav状态估计方法 - Google Patents
一种基于矩阵变换的双线性wlav状态估计方法 Download PDFInfo
- Publication number
- CN107947160A CN107947160A CN201711223404.7A CN201711223404A CN107947160A CN 107947160 A CN107947160 A CN 107947160A CN 201711223404 A CN201711223404 A CN 201711223404A CN 107947160 A CN107947160 A CN 107947160A
- Authority
- CN
- China
- Prior art keywords
- mrow
- matrix
- node
- mtd
- formula
- 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
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种基于矩阵变换的双线性WLAV状态估计方法,其特征包括:1、获得电力系统的网络参数;2、得到节点导纳矩阵;3、从电力系统中获取遥测数据;4、计算各量测量的权重;5、计算雅克比矩阵;6、通过坐标系线性变换的方法实现对连接长短支路的杠杆点注入量测的抗差估计。本发明能降低杠杆量测点坏数据对状态估计精度的影响,从而能得到更精确的状态估计结果。
Description
技术领域
本发明涉及一种基于矩阵变换的双线性WLAV状态估计方法,属于电力系统状态估计领域。
背景技术
电力系统状态估计是能量管理系统中的核心基础模块,准确的状态估计结果能为现代电网的实时分析、控制提供保障。抗差估计能够有效利用正确数据并抑制异常数据,当量测系统含有不良数据或者量测误差偏离原分布模式时,仍然能够得出精确的状态量估计值。WLAV(Weighted LeastAbsolute Value)是一种被广泛研究的抗差状态估计方法,该方法能够有效抑制一般坏数据,提高状态估计精度,但非线性内点法的计算效率不能满足工程实际需求。变量代换方法,将状态估计问题分解为两步线性化过程,提出了双线性WLAV抗差估计,该方法与传统的WLAV状态估计相比,精度略有降低,但运算效率有显著提升。
对杠杆点坏数据缺乏抗差能力是WLAV抗差估计的主要缺陷之一,对WLAV状态估计进行改进得到的双线性WLAV状态估计依然存在此问题。对此的一种处理方法为:在状态估计前进行杠杆点检测,并剔除相关量测,但此方法会导致局部量测冗余度降低,影响电网可观测性。
发明内容
本发明为克服现有技术中双线性WLAV状态估计对杠杆量测坏数据缺乏抗差能力这一问题,提供一种基于矩阵变换的双线性WLAV状态估计方法,以期能降低杠杆量测点坏数据对状态估计精度的影响,从而能得到更精确的状态估计结果。
本发明为实现上述目的,采用如下技术方案:
本发明一种基于矩阵变换的双线性WLAV状态估计方法的特点是按如下步骤进行:
步骤1、获得电力系统的网络参数,包括:输电线路的支路号、首端节点和末端节点编号、串联电阻、串联电抗、并联电导、并联电纳、变压器变比和阻抗;
步骤2、初始化,包括:对状态量x设置初值,所述状态量x包括节点电压的幅值V和相角θ;根据所述电力系统的网络参数得到节点导纳矩阵;
步骤3、从所述电力系统中获取遥测数据z,包括:电压幅值、节点注入有功功率、节点注入无功功率、线路首端有功功率、线路首端无功功率、线路末端有功功率以及线路末端无功功率;所述遥测数据z即为电力系统状态估计的量测量;
步骤4、根据量测仪表的精度,计算各量测量的权重;
步骤5、取中间变量y,令:
式(1)中,n为系统的节点个数;Vi,Vj分别为节点i、节点j的电压幅值,θij为连接节点i和节点j的支路两端相角差,θij=θi-θj,θi、θj分别为节点i、节点j的电压相角;
步骤6、根据所述节点导纳矩阵计算雅克比矩阵C;
步骤7、根据所述雅克比矩阵C,计算对角阵Sr和矩阵Q,所述对角阵Sr用于对雅克比矩阵C的行向量进行尺度缩放,所述矩阵Q用于对雅克比矩阵进行旋转变换;
步骤8、利用内点法求解式(2)得到中间变量y:
式(2)中,是对y进行旋转后得到的矩阵,w为权重矩阵,是对角元素为各量测量权重的对角矩阵,l,u为松弛变量,且l,u≥0;
步骤9、对中间变量y进行非线性变换,获得如式(3)所示变换后的中间变量y′:
步骤10、根据所述变换后的中间变量y′,利用内点法求得中间变量x′=[lnV θ]T;步骤11、利用式(4)得到所述状态量x,从而实现电力系统状态估计:
与现有技术相比,本发明的有益效果体现在:
1、本发明通过对雅克比矩阵进行尺度缩放及坐标系线性变换的方法实现对杠杆量测的抗差估计,与传统的双线性WLAV状态估计相比,对杠杆量测坏数据有更好的抗差能力。
2、考虑到对杠杆点坏数据缺乏抗差能力是WLAV抗差估计的主要缺陷之一,本发明通过雅克比矩阵的尺度缩放来实现双线性WLAV状态估计方法对长短支路导致的支路杠杆点量测抗差能力的改进,进一步通过坐标系线性变换的方法实现对连接长短支路的杠杆点注入量测的抗差估计,该方法在不剔除杠杆量测坏数据的情况下改进对其的抗差能力,不会影响电网可观测性。
附图说明
图1为本发明的流程图;
图2为本发明采用的线路π形等值电路图。
具体实施方式
下面结合附图对发明的技术方案进行详细说明:
本实施例中,如图1所示,一种基于矩阵变换的双线性WLAV状态估计方法是按如下步骤进行:
步骤1、获得电力系统的网络参数,包括:输电线路的支路号、首端节点和末端节点编号、串联电阻、串联电抗、并联电导、并联电纳、变压器变比和阻抗;
步骤2、初始化,包括:对状态量x设置初值,状态量x包括节点电压的幅值V和相角θ;根据电力系统的网络参数得到节点导纳矩阵;
步骤3、从电力系统中获取遥测数据z,包括:电压幅值、节点注入有功功率、节点注入无功功率、线路首端有功功率、线路首端无功功率、线路末端有功功率以及线路末端无功功率;遥测数据z即为电力系统状态估计的量测量;
电力系统状态估计属于非线性优化问题,其量测方程可表示为:
z=h(x)+ε (1)
式(1)中,h(x)为非线性函数;ε为量测误差向量。
步骤4、根据量测仪表的精度,计算各量测量的权重;
步骤5、取中间变量y,令:
式(2)中,n为系统的节点个数;Vi,Vj分别为节点i、节点j的电压幅值,θij为连接节点i和节点j的支路两端相角差,θij=θi-θj,θi、θj分别为节点i、节点j的电压相角;
步骤6、根据节点导纳矩阵计算雅克比矩阵C;
结合电力系统的特性,通过引入中间变量可将非线性方程式(1)分解为:
z=Cy+ε (3)
y′=f(y) (4)
y′=Dx′+εy′ (5)
x=t(x′) (6)
式(3)-式(6)中,y,y′,x′为中间变量;f(y)和t(x′)为非线性函数;C和D为常系数矩阵;εy′为中间变量y′的误差向量。
状态估计量测量包括节点电压量测、节点功率量测、支路有功和无功量测,在极坐标系下,其计算方程为:
式(7)中,Pi,Qi分别为节点i的注入有功及无功功率;Pij,Qij分别为支路ij上的有功和无功功率;Gij,Bij分别为节点i和节点j之间的互电导和互电纳;Gii,Bii分别为节点i的自电导和自电纳;gij,bij,yc分别为图2所示的支路ij的π形等效电路的串联电导、串联电纳和对地电纳。
由上式可得雅克比矩阵C的各元素值如表1所示。
表1雅克比矩阵C
步骤7、根据雅克比矩阵C,计算对角阵Sr和矩阵Q,对角阵Sr用于对雅克比矩阵C的行向量进行尺度缩放,矩阵Q用于对雅克比矩阵进行旋转变换;
步骤8、利用内点法求解式(8)得到中间变量y:
式(8)中,是对y进行旋转后得到的矩阵,w为权重矩阵,是对角元素为各量测量权重的对角矩阵,l,u为松弛变量,且l,u≥0;
构造出式(8)的拉格朗日函数:
式(9)中,α,β,η为拉格朗日乘子。
由KKT(Karush-Kuhn-Tucker)条件可得如下KKT方程:
式(10)中,A,B,L,U分别为以α,β,l,u为对角元素的对角阵,μ为扰动因子,e为单位列向量。应用牛顿法解如上KKT方程,可得修正方程:
-Srdη-Srdα=-Ll (13)
Srdη-Srdβ=-Lu (14)
取初值η=0,α=β=w则Ll=Lu=0,联立式(13)-(16)可得出如下结果:
将式(17)带入式(11)-(12)可得到修正方程的矩阵形式:
式(18)中
步骤9、对中间变量y进行非线性变换,获得如式(20)所示变换后的中间变量y′:
步骤10、根据变换后的中间变量y′,利用内点法求得中间变量x′=[lnV θ]T;
式(21)中,I为单位矩阵;L为节点关联矩阵;Lr为去掉平衡节点的节点关联矩阵。
中间变量y′的权重矩阵为:
Wy′=F-TCTWCF-1 (22)
式(22)中,F为非线性函数f(y)对y求偏导所得的雅克比矩阵。
由式(22)可知,中间变量y′的权重矩阵Wy′为实对称矩阵,由矩阵理论分析可知,所有实对称矩阵均正交相似于由其特征值组成的对角矩阵。
用主轴化方法将其转化为标准型矩阵如下:
式(23)中,Λ为以矩阵Wy′特征值为对角元素的对角阵,H为矩阵线性无关的特征向量组成的正交阵。
式(5)的WLS目标函数为:
min(y′-Dx)TWy′(y′-Dx) (24)
将式(23)带入式(24):
min(HTy′-HTDx)TΛ(HTy′-HTDx) (25)
理论证明矩阵变换前后两者二次型函数值是相等的,即由式(24)和式(25)所解得的状态量x是相同的。将变换后对角阵代替原非对角权重矩阵,建立基于式(5)的WLAV状态估计模型为:
利用内点法求解式(26)时,修正方程为:
迭代求解式(27)即可获得分步模型中变量x′的WLAV解。
步骤11、利用式(28)得到状态量x,从而实现电力系统状态估计:
Claims (1)
1.一种基于矩阵变换的双线性WLAV状态估计方法,其特征是按如下步骤进行:
步骤1、获得电力系统的网络参数,包括:输电线路的支路号、首端节点和末端节点编号、串联电阻、串联电抗、并联电导、并联电纳、变压器变比和阻抗;
步骤2、初始化,包括:对状态量x设置初值,所述状态量x包括节点电压的幅值V和相角θ;根据所述电力系统的网络参数得到节点导纳矩阵;
步骤3、从所述电力系统中获取遥测数据z,包括:电压幅值、节点注入有功功率、节点注入无功功率、线路首端有功功率、线路首端无功功率、线路末端有功功率以及线路末端无功功率;所述遥测数据z即为电力系统状态估计的量测量;
步骤4、根据量测仪表的精度,计算各量测量的权重;
步骤5、取中间变量y,令:
式(1)中,n为系统的节点个数;Vi,Vj分别为节点i、节点j的电压幅值,θij为连接节点i和节点j的支路两端相角差,θij=θi-θj,θi、θj分别为节点i、节点j的电压相角;
步骤6、根据所述节点导纳矩阵计算雅克比矩阵C;
步骤7、根据所述雅克比矩阵C,计算对角阵Sr和矩阵Q,所述对角阵Sr用于对雅克比矩阵C的行向量进行尺度缩放,所述矩阵Q用于对雅克比矩阵进行旋转变换;
步骤8、利用内点法求解式(2)得到中间变量y:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mi>min</mi>
</mtd>
<mtd>
<mrow>
<msup>
<mi>w</mi>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>+</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>S</mi>
<mi>r</mi>
</msub>
<mi>z</mi>
<mo>-</mo>
<msub>
<mi>S</mi>
<mi>r</mi>
</msub>
<mi>C</mi>
<mi>Q</mi>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mo>+</mo>
<msub>
<mi>S</mi>
<mi>r</mi>
</msub>
<mi>l</mi>
<mo>-</mo>
<msub>
<mi>S</mi>
<mi>r</mi>
</msub>
<mi>u</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式(2)中,是对y进行旋转后得到的矩阵,w为权重矩阵,是对角元素为各量测量权重的对角矩阵,l,u为松弛变量,且l,u≥0;
步骤9、对中间变量y进行非线性变换,获得如式(3)所示变换后的中间变量y′:
步骤10、根据所述变换后的中间变量y′,利用内点法求得中间变量x′=[lnV θ]T;
步骤11、利用式(4)得到所述状态量x,从而实现电力系统状态估计:
<mrow>
<mi>x</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mi>ln</mi>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>)</mo>
</mrow>
</mrow>
</msup>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>&theta;</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711223404.7A CN107947160A (zh) | 2017-11-28 | 2017-11-28 | 一种基于矩阵变换的双线性wlav状态估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711223404.7A CN107947160A (zh) | 2017-11-28 | 2017-11-28 | 一种基于矩阵变换的双线性wlav状态估计方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107947160A true CN107947160A (zh) | 2018-04-20 |
Family
ID=61949543
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711223404.7A Pending CN107947160A (zh) | 2017-11-28 | 2017-11-28 | 一种基于矩阵变换的双线性wlav状态估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107947160A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110232199A (zh) * | 2018-11-30 | 2019-09-13 | 南京南瑞继保电气有限公司 | 一种考虑矩阵尺度伸缩的双线性wlav抗差状态估计方法 |
CN112039198A (zh) * | 2020-07-22 | 2020-12-04 | 清华大学 | 基于dane的电网拓扑结构关键特征提取方法和装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102801162A (zh) * | 2012-08-23 | 2012-11-28 | 清华大学 | 一种两阶段线性加权最小二乘电力系统状态估计方法 |
CN103886193A (zh) * | 2014-03-13 | 2014-06-25 | 河海大学 | 一种电力系统模糊自适应抗差估计方法 |
CN105303269A (zh) * | 2015-11-27 | 2016-02-03 | 华北电力大学 | 一种消除杠杆点的最优变换方法 |
CN106709195A (zh) * | 2016-12-30 | 2017-05-24 | 河海大学 | 计及等式约束的双线性wlav状态估计方法 |
CN105119276B (zh) * | 2015-08-18 | 2018-09-21 | 河海大学 | 一种基于admm的电力系统wlav抗差估计方法 |
-
2017
- 2017-11-28 CN CN201711223404.7A patent/CN107947160A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102801162A (zh) * | 2012-08-23 | 2012-11-28 | 清华大学 | 一种两阶段线性加权最小二乘电力系统状态估计方法 |
CN103886193A (zh) * | 2014-03-13 | 2014-06-25 | 河海大学 | 一种电力系统模糊自适应抗差估计方法 |
CN105119276B (zh) * | 2015-08-18 | 2018-09-21 | 河海大学 | 一种基于admm的电力系统wlav抗差估计方法 |
CN105303269A (zh) * | 2015-11-27 | 2016-02-03 | 华北电力大学 | 一种消除杠杆点的最优变换方法 |
CN106709195A (zh) * | 2016-12-30 | 2017-05-24 | 河海大学 | 计及等式约束的双线性wlav状态估计方法 |
Non-Patent Citations (1)
Title |
---|
厉超等: "基于变量代换内点法的加权最小绝对值抗差状态估计", 《电力系统自动化》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110232199A (zh) * | 2018-11-30 | 2019-09-13 | 南京南瑞继保电气有限公司 | 一种考虑矩阵尺度伸缩的双线性wlav抗差状态估计方法 |
CN110232199B (zh) * | 2018-11-30 | 2023-12-19 | 南京南瑞继保电气有限公司 | 一种考虑矩阵尺度伸缩的双线性wlav抗差状态估计方法 |
CN112039198A (zh) * | 2020-07-22 | 2020-12-04 | 清华大学 | 基于dane的电网拓扑结构关键特征提取方法和装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102801162B (zh) | 一种两阶段线性加权最小二乘电力系统状态估计方法 | |
CN104092212B (zh) | 一种基于pmu量测的电力系统多区域分布式状态估计方法 | |
CN103944165B (zh) | 一种大电网参数辨识估计方法 | |
CN105958470B (zh) | 一种基于双线性凸优化理论电力系统双线性抗差估计方法 | |
CN107577870A (zh) | 基于同步相量量测的配电网电压功率灵敏度鲁棒估计方法 | |
CN104778367B (zh) | 基于单一状态断面的广域戴维南等值参数在线计算方法 | |
CN103886193B (zh) | 一种电力系统模糊自适应抗差估计方法 | |
TWI479159B (zh) | 電壓穩定度即時估測方法 | |
CN103279676B (zh) | 一种基于变量代换的电力系统wlav抗差估计方法 | |
CN106026086B (zh) | 一种电网中运行状态的动态估计方法 | |
CN103020726A (zh) | 面向全pmu量测的抗差状态估计方法 | |
CN107947160A (zh) | 一种基于矩阵变换的双线性wlav状态估计方法 | |
CN105512502A (zh) | 一种基于残差归一化的权函数最小二乘状态估计方法 | |
CN101788608B (zh) | 一种独立三绕组变压器电抗参数的估计方法 | |
CN104182644A (zh) | 一种融合分布式电源特性的配电网状态估计方法 | |
CN107104430A (zh) | 一种配电网双解耦状态估计方法 | |
CN106709195A (zh) | 计及等式约束的双线性wlav状态估计方法 | |
CN107167707A (zh) | 一种基于参数未知的同杆双回线故障测距方法及装置 | |
CN103838962B (zh) | 一种计及pmu量测的分步线性状态估计方法 | |
CN106159941B (zh) | 一种考虑实际量测误差传递特性的电力系统状态估计方法 | |
CN109217304A (zh) | 基于wams系统量测量的线路参数辨识方法 | |
CN104252571B (zh) | 基于多预测‑校正内点法的wlav抗差状态估计方法 | |
CN102280877A (zh) | 一种多量测断面的电力系统不良支路参数辨识方法 | |
CN109829246A (zh) | 一种基于参数综合可疑度的线路参数辨识方法 | |
CN105811394B (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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180420 |