CN110232199B - 一种考虑矩阵尺度伸缩的双线性wlav抗差状态估计方法 - Google Patents
一种考虑矩阵尺度伸缩的双线性wlav抗差状态估计方法 Download PDFInfo
- Publication number
- CN110232199B CN110232199B CN201811454693.6A CN201811454693A CN110232199B CN 110232199 B CN110232199 B CN 110232199B CN 201811454693 A CN201811454693 A CN 201811454693A CN 110232199 B CN110232199 B CN 110232199B
- Authority
- CN
- China
- Prior art keywords
- equation
- measurement
- matrix
- bilinear
- power
- 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.)
- Active
Links
- 239000011159 matrix material Substances 0.000 title claims abstract description 67
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000005259 measurement Methods 0.000 claims abstract description 48
- 239000013598 vector Substances 0.000 claims description 23
- 238000003491 array Methods 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000012937 correction Methods 0.000 claims description 9
- 230000009466 transformation Effects 0.000 claims description 8
- 230000005540 biological transmission Effects 0.000 claims description 4
- 230000009977 dual effect Effects 0.000 claims description 4
- 230000014509 gene expression Effects 0.000 claims description 4
- 238000001514 detection method Methods 0.000 abstract description 2
- 230000000694 effects Effects 0.000 description 7
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- 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)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Human Resources & Organizations (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Public Health (AREA)
- Power Engineering (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Complex Calculations (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公布了一种考虑矩阵尺度伸缩的双线性WLAV抗差状态估计方法。双线性WLAV抗差估计通过引入中间变量,将电力系统非线性的量测方程转变为两步线性化过程,简化了模型,降低了估计时间。但当电网中出现长短支路情况,会导致杠杆量测和信息矩阵病态,影响双线性WLAV抗差估计的数值稳定性和坏数据检测辨识。本发明通过对雅克比矩阵进行尺度伸缩,可以有效消除长短支路导致的杠杆点量测,并降低信息矩阵的条件数,提高数值稳定性,具有工程应用价值。
Description
技术领域
发明涉及一种考虑矩阵尺度伸缩的双线性WLAV抗差状态估计方法,属于电力系统运行和控制技术领域。
背景技术
加权最小绝对值(Weighted LeastAbsolute Values,WLAV)估计是误差数据服从拉普拉斯分布下的最优估计,传统基于原对偶内点法(Primal-Dual IPM,PDIPM)的WLAV抗差估计存在2个不足:一是模型复杂,迭代系数矩阵维数大,需要计算海森矩阵,计算时间远大于加权最小二乘法(Weighted Least Squares,WLS);二是对于杠杆点坏数据和多相关一致坏数据的抗差效果差。为提高WLAV抗差估计的计算速度,一般通过计及高阶项,并利用高阶内点法修正障碍参数,提高单次迭代下修正量的幅值和方向精度,从而有效减少抗差估计的迭代次数并提高估计精度。但由于高阶项的计算,该类方法单次迭代的计算时间增大,无法有效提高整体的计算速度,且仍然无法检测出杠杆点坏数据。
双线性WLAV抗差估计通过引入中间变量,实现状态估计的非线性量测方程分步线性化。将内点法引入分步线性化后的量测方程进行WLAV抗差估计,无需形成海森矩阵,同时减小了库恩-塔克(Karush-Kuhn-Tucker,KKT)方程组中系数矩阵的维数。双线性WLAV方法中迭代方程系数矩阵中不包含海森矩阵,且稀疏度高,矩阵阶数小,有效提高了WLAV抗差估计的计算效率,计算时间远小于传统的基于非线性量测方程的PDIPM,并且具有良好的抗差性能。
可是当电网出现长短支路时,系统中便会存在杠杆量测并造成信息矩阵病态,而双线性WLAV抗差估计对杠杆量测并没有抗差效果,影响了该方法的数值稳定性以及对坏数据的检测辨识,限制了双线性WLAV抗差估计的工程应用。本发明通过对雅克比矩阵进行尺度伸缩,将杠杆量测转化为非杠杆量测,实现了对杠杆点的抗差效果,解决了双线性WLAV抗差估计存在的缺陷,具有工程应用价值。
发明内容
发明目的:双线性WLAV抗差估计在电网存在长短支路时,无法检测辨别杠杆点量测,且信息矩阵的条件数相当大,造成状态估计结果严重偏离真实值,而本发明提出的一种考虑矩阵尺度伸缩的双线性WLAV抗差状态估计方法能够解决现有技术存在的缺陷。
技术方案:本发明提出一种考虑矩阵尺度伸缩的双线性WLAV抗差状态估计方法,包括以下步骤:
(1)获得电力系统的网络参数和量测量。网络参数包括:母线编号、名称、补偿电容,输电线路的支路号、首端节点和末端节点编号、串联电阻、串联电抗、并联电导、并联电纳、变压器变比和阻抗;量测量z包括:节点电压幅值、支路首端有功功率和无功功率、支路末端有功功率和无功功率;
(2)程序初始化。包括:设置迭代精度、最大迭代次数、原对偶内点法的松弛变量、拉格朗日乘子和罚因子初值,并利用(1)中的网络参数形成节点导纳矩阵;
(3)设置中间变量:
y=[Mij Lij Oi]T。
式中:Mij=UiUjcosθij,Lij=UiUjsinθij,Oi=Ui 2,Ui和Uj表示节点i和j对应的电压幅值,θij表示节点i和j之间的相角差;
(4)对雅克比矩阵进行尺度伸缩,并利用中间变量建立如下的分步状态估计模型:
y'=f(y)
x=t(x')
式中:x为待求的状态量;y、y′和x′为中间变量,f(y)和t(x′)为非线性函数;εz和εy′分别为量测量z和中间变量y′的误差向量;C和D分别表示两步线性方程各自的雅克比矩阵;Sr和Sc分别表示常系数雅克比矩阵C行向量和列向量进行尺度伸缩的对角阵,其中Sr和Sc的对角元素分别取C中行和列的最大值;Sr'和Sc'分别表示雅克比矩阵D行向量和列向量进行尺度伸缩的对角阵,其中Sr'和Sc'的对角元素分别取D中行和列的最大值;
(5)针对第一步线性方程,建立WLAV抗差估计模型:
式中:w为量测量z的m维权重向量,l和u为松弛变量;
(6)根据内点法,建立上式的拉格朗日函数:
式中:η、α和β为拉格朗日乘子;
(7)根据拉格朗日函数L,列出库恩-塔克方程,并进行泰勒展开,获得修正方程;
(8)将(7)中求解得到的y代入第二步线性方程,进行和第一步线性方程同样的操作,解出状态量x,输出结果。
有益效果:本发明提供一种考虑矩阵尺度伸缩的双线性WLAV抗差状态估计方法,通过对双线性WLAV抗差模型中两步线性过程的雅克比矩阵进行伸缩变换,将杠杆量测转化为非杠杆量测,不仅实现了对杠杆点的抗差效果,提高了估计结果的精度,还降低了信息矩阵的条件数,提高了数值稳定性,具有工程应用价值。
附图说明
图1为本发明方法流程图。
图2为本发明方法所应用的算例系统,采用的是IEEE14节点系统。
具体实施方式
下面结合附图对发明的技术流程进行详细说明:
如图1所示,本发明具体包括如下步骤:
1.获得电力系统的网络参数和量测量。网络参数包括:母线编号、名称、补偿电容,输电线路的支路号、首端节点和末端节点编号、串联电阻、串联电抗、并联电导、并联电纳、变压器变比和阻抗;量测量z包括:节点电压幅值、支路首端有功功率和无功功率、支路末端有功功率和无功功率。
2.程序初始化。包括:设置迭代精度、最大迭代次数、原对偶内点法的松弛变量、拉格朗日乘子和罚因子初值,并利用(1)中的网络参数形成节点导纳矩阵。
3.设置中间变量:
y=[Mij Lij Oi]T。
式中:Mij=UiUjcosθij,Lij=UiUjsinθij,Oi=Ui 2,Ui和Uj表示节点i和j对应的电压幅值,θij表示节点i和j之间的相角差。
4.对雅克比矩阵进行尺度伸缩,并利用中间变量建立如下的分步状态估计模型,具体过程如下:
电力系统状态估计的量测方程为:
z=h(x)+ε
式中:x为状态量(维数n=2N-1,N为节点数);z为量测量(维数m,m>n);h为m维非线性量测函数;ε为m维量测误差。
电力系统支路功率方程表达式为:
式中:Pij和Qij分别为线路i-j上始端有功功率和无功功率,其方向规定:由i流向j为正,由j流向i为负;g为线路电导,b为线路电纳,yc为线路对地导纳值。
采用第3步的中间变量后,支路功率方程可以化为如下格式:
Pij=Oig-Mijg-Lijb
Qij=-Oi(b+yc)+Mijb-Lijg
易看出上式中支路功率与变量y为线性关系,若将电压量测量设为平方的格式,则所有量测量与y均为线性关系。
设电网共有k条支路,可知y的维数为N+2k,令m>n,则中间变量的维数小于量测量维数,建立如下的量测量-中间变量函数方程:
z=Cy+εz
式中:εz为量测量z的误差向量,C表示该线性方程的雅克比矩阵。
通过一步计算即可得到y的加权最小二乘解:
y=(CTWC)-1CTWz
式中:W表示量测量z的权重矩阵。
对y做适当的变换如下:
同时对状态量x中电压幅值取其对数,将其化为:
可得y′与x′的函数表达式为:
式中:I为单位阵,F为支路节点关联矩阵,FR为删除平衡节点的节点支路关联矩阵。y′与x′为线性关系,D表示雅克比矩阵,则x′的最小二乘解为:
x′=(DTWy′D)DTWy′y′
式中:Wy′为y′对应的权重矩阵。求出x′后,对其中的对数电压部分取对数即可获得电压幅值的估计值。
上面即为双线性WLAV抗差估计模型,而本发明对上述两步线性方程的雅克比矩阵C和D进行尺度伸缩,以消除长短支路导致的杠杆点量测,最后建立考虑矩阵尺度伸缩的双线性WLAV抗差状态估计模型:
y'=f(y)
x=t(x')
式中:x为待求的状态量;y、y′和x′为中间变量,f(y)和t(x′)为非线性函数;εz和εy′分别为量测量z和中间变量y′的误差向量;C和D分别表示两步线性方程各自的雅克比矩阵;Sr和Sc分别表示常系数雅克比矩阵C行向量和列向量进行尺度伸缩的对角阵,其中Sr和Sc的对角元素分别取C中行和列的最大值;Sr'和Sc'分别表示雅克比矩阵D行向量和列向量进行尺度伸缩的对角阵,其中Sr'和Sc'的对角元素分别取D中行和列的最大值。
5.以第一步线性方程为例,建立本发明的WLAV抗差估计模型,具体过程如下:
电力系统中WLAV状态估计的数学模型为:
min wT|ε|
s.t.ε=z-h(x)
式中:w=[w1,w2,…,wm]T为权重向量,σi为第i个量测量误差的标准差。
上式的目标函数中有绝对值量,不满足连续可微条件,因此无法利用拉格朗日法进行求解。添加不等式约束,代替上述模型中等式约束:
添加松弛变量2l和2u,将上式中的不等式约束转化为等式约束:
再将上式中两个等式约束相加减可得新的等式约束:
代入最原始的方程中,获得目标函数连续可微的WLAV抗差估计模型化简为:
而本发明的双线性WLAV抗差状态估计模型对雅克比矩阵进行了伸缩变换,依据第4步中考虑矩阵尺度伸缩的双线性WLAV抗差状态估计模型,以第一步线性方程为例,建立本发明方法的WLAV抗差估计模型,如下:
6.根据内点法,建立上式的拉格朗日函数:
式中:η、α和β为拉格朗日乘子;
7.根据拉格朗日函数L,列出库恩-塔克方程,并进行泰勒展开,获得修正方程,过程如下:
拉格朗日函数L的库恩-塔克等式为:
式中:(A,L,B,U)分别是以(α,l,β,u)为对角元素的对角阵;μ为扰动因子。
将库恩-塔克等式泰勒展开,保留一阶项,获得修正方程如下:
-dη-dα=-Ll
dη-dβ=-Lu
取初值η=0,α=β=w,则Ll=Lu=0,代入库恩-塔克等式中可以得到:
进而对应的迭代方程可以简化为:
式中:中间变量/>将求解得到的dη回代,即可获得其他变量的修正量。
8.将第7步中求解得到的y代入第二步线性方程,建立本发明方法的WLAV抗差估计模型,仿照第一步线性方程的处理过程,解出状态量x,输出结果。
上述方法具体步骤如下:
(1)获得电力系统的网络参数和量测量。网络参数包括:母线编号、名称、补偿电容,输电线路的支路号、首端节点和末端节点编号、串联电阻、串联电抗、并联电导、并联电纳、变压器变比和阻抗;量测量z包括:节点电压幅值、支路首端有功功率和无功功率、支路末端有功功率和无功功率;
(2)程序初始化。包括:设置迭代精度、最大迭代次数、原对偶内点法的松弛变量、拉格朗日乘子和罚因子初值,并利用(1)中的网络参数形成节点导纳矩阵;
(3)设置中间变量:
y=[Mij Lij Oi]T。
式中:Mij=UiUjcosθij,Lij=UiUjsinθij,Oi=Ui 2,Ui和Uj表示节点i和j对应的电压幅值,θij表示节点i和j之间的相角差;
(4)对雅克比矩阵进行尺度伸缩,并利用中间变量建立如下的分步状态估计模型:
y'=f(y)
x=t(x')
式中:x为待求的状态量;y、y′和x′为中间变量,f(y)和t(x′)为非线性函数;εz和εy′分别为量测量z和中间变量y′的误差向量;C和D分别表示两步线性方程各自的雅克比矩阵;Sr和Sc分别表示常系数雅克比矩阵C行向量和列向量进行尺度伸缩的对角阵,其中Sr和Sc的对角元素分别取C中行和列的最大值;Sr'和Sc'分别表示雅克比矩阵D行向量和列向量进行尺度伸缩的对角阵,其中Sr'和Sc'的对角元素分别取D中行和列的最大值;
(5)针对第一步线性方程,建立WLAV抗差估计模型:
式中:w为量测量z的m维权重向量,l和u为松弛变量;
(6)根据内点法,建立上式的拉格朗日函数:
式中:η、α和β为拉格朗日乘子;
(7)根据拉格朗日函数L,列出库恩-塔克方程,并进行泰勒展开,获得修正方程;
(8)将(7)中求解得到的y代入第二步线性方程,进行和第一步线性方程同样的操作,解出状态量x,输出结果。
本发明提出的一种考虑矩阵尺度伸缩的双线性WLAV抗差状态估计方法,通过对双线性WLAV抗差估计中两步线性方程的雅克比矩阵进行尺度伸缩,将杠杆量测转化为非杠杆量测,实现杠杆点坏数据抗差效果,提高了估计精度,并且降低了信息矩阵的条件数,提高了数值稳定性。
下面介绍本发明的实施例:
本发明采用图2所示的IEEE14节点的标准算例,分析本发明提出的方法对杠杆点坏数据的抗差效果。当系统中含有杠杆量测坏数据时,表1展示了双线性WLAV抗差估计和本发明提出方法的状态估计结果,表2为两种方法的信息矩阵条件数及迭代次数。
表1 IEEE14节点系统杠杆量测坏数据状态估计结果
表2信息矩阵条件数及迭代次数
由表1可以看出,通过雅克比矩阵的尺度伸缩,可以实现杠杆点坏数据的抗差效果,相比双线性WLAV具有更高的精度。由表2发现,本发明提出的方法信息矩阵的条件数大幅度下降,数值稳定性明显提高,而迭代次数几乎不变,相比之下有更好的工程应用价值。
最后应该说明的是,结合上述实施例仅说明本发明的技术方案而非对其限制。所属领域的普通技术人员应当理解到,本领域技术人员可以对本发明的具体实施方式进行修改或者等同替换,但这些修改或变更均在申请待批的权利要求保护范围之中。
Claims (2)
1.一种考虑矩阵尺度伸缩的双线性WLAV抗差状态估计方法,其特征在于,所述方法具体包括如下步骤:
步骤(1)、获得电力系统的网络参数和量测量,网络参数包括:母线编号、名称、补偿电容,输电线路的支路号、首端节点和末端节点编号、串联电阻、串联电抗、并联电导、并联电纳、变压器变比和阻抗;量测量z包括:节点电压幅值、支路首端有功功率和无功功率、支路末端有功功率和无功功率;
步骤(2)、程序初始化,包括:设置迭代精度、最大迭代次数、原对偶内点法的松弛变量、拉格朗日乘子和罚因子初值,并利用步骤(1)中的网络参数形成节点导纳矩阵;
步骤(3)、设置中间变量:
y=[Mij Lij Oi]T
式中:Mij=UiUjcosθij,Lij=UiUjsinθij,Oi=Ui 2,Ui和Uj表示节点i和j对应的电压幅值,θij表示节点i和j之间的相角差;
步骤(4)、对雅克比矩阵进行尺度伸缩,并利用中间变量建立如下的分步状态估计模型,具体过程如下:
电力系统状态估计的量测方程为:
z=h(x)+ε
式中:x为状态量;z为量测量,维数为m,m>n;h为m维非线性量测函数;ε为m维量测误差;
电力系统支路功率方程表达式为:
式中:Pij和Qij分别为线路i-j上始端有功功率和无功功率,其方向规定:由i流向j为正,由j流向i为负;g为线路电导,b为线路电纳,yc为线路对地导纳值;
采用步骤(3)的中间变量后,支路功率方程化为如下格式:
Pij=Oig-Mijg-Lijb
Qij=-Oi(b+yc)+Mijb-Lijg
上式中支路功率与变量y为线性关系,若将电压量测量设为平方的格式,则所有量测量与y均为线性关系;
设电网共有k条支路,已知y的维数为N+2k,令m>n,则中间变量的维数小于量测量维数,建立如下的量测量-中间变量函数方程:
z=Cy+εz
式中:εz为量测量z的误差向量,C表示该线性方程的雅克比矩阵;
通过一步计算即可得到y的加权最小二乘解:
y=(CTWC)-1CTWz
式中:W表示量测量z的权重矩阵;
对y做适当的变换如下:
同时对状态量x中电压幅值取其对数,将其化为:
可得y′与x′的函数表达式为:
式中:I为单位阵,F为支路节点关联矩阵,FR为删除平衡节点的节点支路关联矩阵,y′与x′为线性关系,D表示雅克比矩阵,则x′的最小二乘解为:
x′=(DTWy′D)DTWy′y′
式中:Wy′为y′对应的权重矩阵,求出x′后,对其中的对数电压部分取对数即可获得电压幅值的估计值;
对线性方程的雅克比矩阵C和D进行尺度伸缩,以消除长短支路导致的杠杆点量测,最后建立考虑矩阵尺度伸缩的双线性WLAV抗差状态估计模型:
y'=f(y)
x=t(x')
式中:x为待求的状态量;y、y′和x′为中间变量,f(y)和t(x′)为非线性函数;εz和εy′分别为量测量z和中间变量y′的误差向量;C和D分别表示两步线性方程各自的雅克比矩阵;Sr和Sc分别表示常系数雅克比矩阵C行向量和列向量进行尺度伸缩的对角阵,其中Sr和Sc的对角元素分别取C中行和列的最大值;Sr'和Sc'分别表示雅克比矩阵D行向量和列向量进行尺度伸缩的对角阵,其中Sr'和Sc'的对角元素分别取D中行和列的最大值;
步骤(5)、针对步骤(4)中的第一步线性方程,建立WLAV抗差估计模型:
式中:w为量测量z的m维权重向量,l和u为松弛变量;
步骤(6)、根据内点法,建立上式的拉格朗日函数:
式中:η、α和β为拉格朗日乘子;
步骤(7)、根据拉格朗日函数L,列出库恩-塔克方程,并进行泰勒展开,获得修正方程;
步骤(8)、将步骤(7)中求解得到的y代入步骤(4)中的第二步线性方程,进行和第一步线性方程同样的操作,解出状态量x,输出结果。
2.根据权利要求1所述的一种考虑矩阵尺度伸缩的双线性WLAV抗差状态估计方法,其特征在于,所述步骤(7)具体包括如下内容:
拉格朗日函数L的库恩-塔克等式为:
式中:A,L,B,U分别是以α,l,β,u为对角元素的对角阵;μ为扰动因子;
将库恩-塔克等式泰勒展开,保留一阶项,获得修正方程如下:
-dη-dα=-Ll
dη-dβ=-Lu
取初值η=0,α=β=w,则Ll=Lu=0,代入库恩-塔克等式中能够得到:
进而对应的迭代方程能够简化为:
式中:中间变量/>将求解得到的dη回代,即可获得其他变量的修正量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811454693.6A CN110232199B (zh) | 2018-11-30 | 2018-11-30 | 一种考虑矩阵尺度伸缩的双线性wlav抗差状态估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811454693.6A CN110232199B (zh) | 2018-11-30 | 2018-11-30 | 一种考虑矩阵尺度伸缩的双线性wlav抗差状态估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110232199A CN110232199A (zh) | 2019-09-13 |
CN110232199B true CN110232199B (zh) | 2023-12-19 |
Family
ID=67861919
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811454693.6A Active CN110232199B (zh) | 2018-11-30 | 2018-11-30 | 一种考虑矩阵尺度伸缩的双线性wlav抗差状态估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110232199B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111191426B (zh) * | 2020-01-08 | 2023-03-14 | 山东大学 | 一种基于Gaussian程序输出文件的数据提取及处理方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102801162A (zh) * | 2012-08-23 | 2012-11-28 | 清华大学 | 一种两阶段线性加权最小二乘电力系统状态估计方法 |
CN103279676A (zh) * | 2013-06-07 | 2013-09-04 | 河海大学 | 一种基于变量代换的电力系统wlav抗差估计方法 |
CN105119276A (zh) * | 2015-08-18 | 2015-12-02 | 河海大学 | 一种基于admm的电力系统wlav抗差估计方法 |
CN105303269A (zh) * | 2015-11-27 | 2016-02-03 | 华北电力大学 | 一种消除杠杆点的最优变换方法 |
WO2016058248A1 (zh) * | 2014-10-15 | 2016-04-21 | 国家电网公司 | 一种基于双线性凸优化理论电力系统双线性抗差估计方法 |
CN105958470A (zh) * | 2014-10-20 | 2016-09-21 | 国家电网公司 | 一种基于双线性凸优化理论电力系统双线性抗差估计方法 |
CN106709195A (zh) * | 2016-12-30 | 2017-05-24 | 河海大学 | 计及等式约束的双线性wlav状态估计方法 |
CN107016489A (zh) * | 2017-03-09 | 2017-08-04 | 中国电力科学研究院 | 一种电力系统抗差状态估计方法和装置 |
CN107947160A (zh) * | 2017-11-28 | 2018-04-20 | 国家电网公司 | 一种基于矩阵变换的双线性wlav状态估计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102982379B (zh) * | 2012-11-16 | 2016-08-03 | 清华大学 | 基于集合论估计模型的电力系统状态估计方法及系统 |
-
2018
- 2018-11-30 CN CN201811454693.6A patent/CN110232199B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102801162A (zh) * | 2012-08-23 | 2012-11-28 | 清华大学 | 一种两阶段线性加权最小二乘电力系统状态估计方法 |
CN103279676A (zh) * | 2013-06-07 | 2013-09-04 | 河海大学 | 一种基于变量代换的电力系统wlav抗差估计方法 |
WO2016058248A1 (zh) * | 2014-10-15 | 2016-04-21 | 国家电网公司 | 一种基于双线性凸优化理论电力系统双线性抗差估计方法 |
CN105958470A (zh) * | 2014-10-20 | 2016-09-21 | 国家电网公司 | 一种基于双线性凸优化理论电力系统双线性抗差估计方法 |
CN105119276A (zh) * | 2015-08-18 | 2015-12-02 | 河海大学 | 一种基于admm的电力系统wlav抗差估计方法 |
CN105303269A (zh) * | 2015-11-27 | 2016-02-03 | 华北电力大学 | 一种消除杠杆点的最优变换方法 |
CN106709195A (zh) * | 2016-12-30 | 2017-05-24 | 河海大学 | 计及等式约束的双线性wlav状态估计方法 |
CN107016489A (zh) * | 2017-03-09 | 2017-08-04 | 中国电力科学研究院 | 一种电力系统抗差状态估计方法和装置 |
CN107947160A (zh) * | 2017-11-28 | 2018-04-20 | 国家电网公司 | 一种基于矩阵变换的双线性wlav状态估计方法 |
Non-Patent Citations (4)
Title |
---|
一种处理杠杆量测的配电网抗差估计算法;李慧等;《电力系统自动化》;20050315(第03期);第35-39页 * |
基于多中心-校正内点法的加权最小绝对值状态估计;刘小波;《贵州电力技术》;20140929(第09期);第32-35页 * |
基于多预测校正内点法的WLAV抗差状态估计;颜全椿 等;《电网技术》;20130805;第2194-2200页 * |
考虑零注入约束的双线性WLAV状态估计;沈甜甜等;《河海大学学报(自然科学版)》;20171125(第06期);第90-95页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110232199A (zh) | 2019-09-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2019056753A1 (zh) | 一种分布式光伏电站集群的动态等值建模方法 | |
CN109888773B (zh) | 一种电力系统多区域分布式状态评估方法 | |
CN103279676B (zh) | 一种基于变量代换的电力系统wlav抗差估计方法 | |
CN108804386B (zh) | 一种电力系统负荷裕度的并行化计算方法 | |
CN107749627B (zh) | 基于改进匹配追踪的智能配电网潮流雅可比矩阵估计方法 | |
CN109617080B (zh) | 基于改进的雅可比矩阵的直角坐标牛顿法潮流计算方法 | |
CN110232199B (zh) | 一种考虑矩阵尺度伸缩的双线性wlav抗差状态估计方法 | |
CN116706921B (zh) | 基于hhl算法的量子牛顿-拉夫逊法潮流计算方法和系统 | |
CN115453193B (zh) | 基于pqm、ttu和sm量测数据协同的配电网谐波状态估计方法 | |
CN111355236A (zh) | 一种计及中性点电压变量的配电网三相潮流计算方法 | |
CN108918930A (zh) | 一种负荷分解中功率信号自适应重构方法 | |
WO2016058248A1 (zh) | 一种基于双线性凸优化理论电力系统双线性抗差估计方法 | |
CN106410811B (zh) | 首次迭代小阻抗支路端点改变雅可比矩阵的潮流计算方法 | |
CN110148946B (zh) | 一种基于辅助因子两步求解的孤岛微电网潮流计算方法 | |
CN117112996A (zh) | 一种基于自适应尺度分解的海洋气象数据同化方法 | |
CN109255541B (zh) | 一种基于最小二乘和一乘之和的配电网抗差状态估计方法 | |
CN109061544B (zh) | 一种电能计量误差估计方法 | |
CN111049144A (zh) | 基于定变化率矩阵的有源配电网潮流计算pv节点处理方法 | |
CN113890042B (zh) | 基于有源配电网的潮流计算方法 | |
CN115169246A (zh) | 基于Direct优化和双样本点更新的可靠性分析方法 | |
CN112636325B (zh) | 一种配电网抗差状态计算方法、装置、终端及存储介质 | |
CN111505379B (zh) | 一种基于奇异值分解的配电网同步相量量测数据压缩方法 | |
CN114372036A (zh) | 电力系统的状态估计方法、装置、设备和计算机存储介质 | |
CN110336288B (zh) | 基于矩阵运算提取雅可比元素的配电网三相潮流计算方法 | |
CN110021928B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |