CN105552904B - 基于双线性化的多区域电网全分布式抗差状态估计方法 - Google Patents
基于双线性化的多区域电网全分布式抗差状态估计方法 Download PDFInfo
- Publication number
- CN105552904B CN105552904B CN201610067806.1A CN201610067806A CN105552904B CN 105552904 B CN105552904 B CN 105552904B CN 201610067806 A CN201610067806 A CN 201610067806A CN 105552904 B CN105552904 B CN 105552904B
- Authority
- CN
- China
- Prior art keywords
- region
- formula
- node
- variable
- stage
- 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
- 238000000034 method Methods 0.000 title claims abstract description 57
- 238000005259 measurement Methods 0.000 claims abstract description 39
- 238000006243 chemical reaction Methods 0.000 claims abstract description 5
- 238000004891 communication Methods 0.000 claims abstract description 5
- 239000011159 matrix material Substances 0.000 claims description 17
- 238000002347 injection Methods 0.000 claims description 12
- 239000007924 injection Substances 0.000 claims description 12
- 230000003416 augmentation Effects 0.000 claims description 6
- 239000000203 mixture Substances 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 abstract description 3
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 206010068052 Mosaicism Diseases 0.000 description 1
- 241001282153 Scopelogadus mizolepis Species 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 210000003765 sex chromosome Anatomy 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 238000004804 winding Methods 0.000 description 1
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
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05F—SYSTEMS FOR REGULATING ELECTRIC OR MAGNETIC VARIABLES
- G05F1/00—Automatic systems in which deviations of an electric quantity from one or more predetermined values are detected at the output of the system and fed back to a device within the system to restore the detected quantity to its predetermined value or values, i.e. retroactive systems
- G05F1/66—Regulating electric power
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B23/00—Testing or monitoring of control systems or parts thereof
- G05B23/02—Electric testing or monitoring
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B23/00—Testing or monitoring of control systems or parts thereof
- G05B23/02—Electric testing or monitoring
- G05B23/0205—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults
- G05B23/0208—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterized by the configuration of the monitoring system
- G05B23/0213—Modular or universal configuration of the monitoring system, e.g. monitoring system having modules that may be combined to build monitoring program; monitoring system that can be applied to legacy systems; adaptable monitoring system; using different communication protocols
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B23/00—Testing or monitoring of control systems or parts thereof
- G05B23/02—Electric testing or monitoring
- G05B23/0205—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults
- G05B23/0218—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults
- G05B23/0243—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults model based detection method, e.g. first-principles knowledge model
- G05B23/0254—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults model based detection method, e.g. first-principles knowledge model based on a quantitative model, e.g. mathematical relationships between inputs and outputs; functions: observer, Kalman filter, residual calculation, Neural Networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- 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)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Automation & Control Theory (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Power Engineering (AREA)
- Radar, Positioning & Navigation (AREA)
- Algebra (AREA)
- Operations Research (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Electromagnetism (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Complex Calculations (AREA)
Abstract
本发明涉及基于双线性化的多区域电网全分布式抗差状态估计方法,属于电力***运行和控制技术领域。该方法分为三阶段:第一阶段是对量测的预处理,得到中间状态变量的预估计;第二阶段把中间状态变量进行非线性转换,得到新的量测,即各支路两端电压相角差以及电压对数幅值之和的两倍。第三阶段基于这些新的量测,估计出最终的状态变量。第二阶段的非线性转换,各区可以并行求解,不需要与其他区域通信。第一阶段与第三阶段的状态估计,针对的是线性化的***,每个控制区利用本地与邻居的通信,采用交替方向乘子法,求得本地状态量。本发明适用于大规模多区域状态估计,能够很好地保护各区域数据隐私,具有很高的敏捷性、灵活性与鲁棒性。
Description
技术领域
本发明属于电力***运行和控制技术领域,特别涉及一种基于双线性化的多区域电网全分布式抗差状态估计方法。
背景技术
电力***的状态估计就是以数据采集与监视控制***采集的实时量测为数据输入源,剔除其中的不良数据之后,结合电网模型,按照特定的状态估计模型,对电网的状态量进行最优估算。状态估计在能量管理***中处于极为重要的地位,是后续一系列高级应用的重要基石;但是,随着***规模增大,现有的集中式的管理***,将面临如下挑战:
(1)海量信息问题:巨大的网络规模极有可能造成通信拥塞和信息处理瓶颈;
(2)可维护问题:互联区域数众多,异动频繁,控制中心难以实时地维护整个电网的全局模型;
(3)隐私性问题:不同控制区可能隶属于不同的运营主体,由于商业机密,控制中心难以采集各控制区的所有信息。
上述挑战都将促使集中式的状态估计变革为全分布式架构,决策机制由单一模式变革为自治模式。全分布式架构不需要协调层来对各区域进行集中管理协调,各区域完全自治,可以并行计算本区域的子问题,并与相邻区域交互边界信息,即可准确获得整个***的状态变量,同时抑制坏数据带来的影响。
发明内容
本发明的目的是为克服已有技术的不足之处,提出一种基于双线性化的多区域电网全分布式抗差状态估计方法,采用本发明提出的方法所进行的状态估计方法,能够很好地保护各子区域的数据隐私,各区域之间只需交互极少量信息,即可协同获得电力***状态,并且能抑制坏数据带来的影响,具有很高的敏捷性与鲁棒性。
本发明提出一种基于双线性化的多区域电网全分布式抗差状态估计方法,其特征在于,该方法包括以下步骤:
1)将电力***分成多个控制区,一个控制区为一条或若干条母线和挂接在上面的发电机;设给定的电力***被划分为R个区域,R为正整数;定义每个控制区中所有节点的集合为节点j表示该区域中任意一节点,mj:=|Mj|为包含节点j的区域数量;
2)对划分后的电力***的每个控制区域建立由目标函数和约束条件组成的分布式第一阶段线性模型,区域a和b为***中任意两个相邻的控制区域,区域a其节点集Na={1,2,4…},区域b节点集Nb={3,5…};
2-1)设任意控制区域a目标函数为最小化全网残差平方和以及保证坏数据稀疏性所引入的惩罚项如式(1):
其中,为区域a的目标函数,ya为区域a的中间状态变量,由{Ua,i,Ka,ij,La,ij}组成,La,ij=Va,iVa,jsinθa,ij,θa,ij=θa,i-θa,j,Va,i和Va,j分别为区域a节点i和节点j的电压幅值,θa,i和θa,j分别为区域a节点i和节点j的电压相角,为惩罚项;为区域a第一阶段坏数据向量,za为区域a的量测包括电压幅值量测、支路有功量测、支路无功量测、节点注入有功量测和注入无功量测,Ba为区域a的第一阶段量测矩阵,其元素由式(2)-(6)量测方程确定:
其中,和是区域a节点i的注入有功和无功量测,和是区域a支路ij的i端的有功和无功量测,是区域a节点i电压量测值的平方,gij为支路ij的电导,gsh,i为节点i的并联电导,bij为支路ij的电纳,bsh,i为节点i的并联电纳,εP,εQ,εU分别为有功量测、无功量测、电压幅值平方的误差项;λ为常数;
2-2)约束约束条件包括:区域a的零注入等式约束,如式(7)所示:
Ea代表区域a的零注入量测矩阵;
为保证相邻区域边界状态量一致所引入的约束,如式(8)所示:
其中,Ka,ij为区域a线路ij两端电压幅值以及相角差之余弦的乘积,La,ij为区域a线路ij两端电压幅值以及相角差之正弦的乘积,Гa,b为区域a和b之间的联络线集合,Δa为a的相邻区域的集合;
3)令交替方向乘子法的迭代下标t=0;给定交替方向乘子法的收敛标准ε∈R+;采用交替方向乘子法对步骤2)建立的分布式第一阶段线性模型进行求解;具体包括:
3-1)更新控制区域a的中间状态变量,如式(9)所示:
ρf∈R+为交替方向乘子法的罚因子,其中R+为正实数集,为增广信息矩阵,I为单位阵,为第t步迭代中区域a的辅助中间状态变量,由与{Ua,Ka,La}所对应的辅助中间状态变量组成,初始值选为为区域a的第一阶段辅助矩阵;
3-2)更新控制区域a的坏数据变量,如式(10)所示:
其中代表阈值算符,相当于对向量中的每个分量做如式(11)运算:
其中l代表分量序号;
3-3)控制区a与相邻控制区b交互边界辅助变量K和L,对相邻控制区a和b之间联络线上的辅助变量K和L求均值,如式(12):
3-4)更新区域a的辅助中间状态变量,如式(13):
其中,均为的分量;
3-5)判断第一阶段中交替方向乘子法是否收敛:
计算原始残差向量如式(14):
计算对偶残差向量如式(15):
计算总残差如式(16):
若δf,t≥εf,则t:=t+1,返回步骤3-1);反之,说明交替方向乘子法已收敛,进行步骤4);
4)对控制区域a的辅助中间状态变量进行第二阶段非线性变换,从而得到分布式第三阶段线性模型所需要的量测,如式(17)-(19):
αa,i=lnUa,i (47)
其中αa,i为区域a节点i的电压幅值自然对数的两倍,αa,ij为区域a支路ij两端αa,i和αa,j之和,θa,ij为区域a的支路ij两端的电压相角差;
5)对控制区a域建立由目标函数和约束条件组成的分布式第三阶段线性模型;
5-1)目标函数如下,为最小化全网残差平方和以及保证坏数据稀疏性所引入的惩罚项
其中为区域a的目标函数,xa为区域a的最终状态变量包括节点电压幅值和相角,为区域第二阶段坏数据向量,为惩罚项;为区域a的经过非线性转换后的量测,包括{αa,i,θa,ij,αa,ij},Ca为区域a的第三阶段量测矩阵,其元素由式(20)、(21)量测方程决定:
αa,ij=αa,i+αa,j+εα (50)
θa,ij=θa,i-θa,j+εθ (51)
其中,εα,εθ是第三阶段量测αa,ij,θa,ij的误差项;
5-2)约束条件是用以保证相邻区域在边界处一致的等式约束,如式(22):
其中,xa,i为区域a节点i的状态变量,为区域a的边界节点集;
6)令交替方向乘子法的迭代下标t=0;给定交替方向乘子法的收敛标准ε∈R+;采用交替方向乘子法对步骤5)的分布式第三阶段线性模型进行求解;
6-1)更新区域a的最终状态变量,如式(23)所示:
ρs∈R+为交替方向乘子法的罚因子,,为第二阶段增广信息矩阵,为第t步迭代中区域a的辅助最终状态变量,初值为0;
6-2)更新区域a的坏数据变量,如式(24)所示:
6-3)控制区a利用本地和与邻居的通信,区域a与相邻区域b交互边界节点的状态变量最新估计值,对区域a任意边界节点i上的状态变量求均值,如式(25)所示:
其中,Mi为包含节点i的区域集合,mi为包含节点i的区域数;
6-4)更新区域a的辅助状态变量,如式(26)所示:
4-5)判断第三阶段中交替方向乘子法是否收敛,计算原始残差向量如式(27)所示:
计算对偶残差向量,如式(28)、(29)所示:
计算总残差,如式(30)所示:
若δs,t≥ε,则t:=t+1,返回步骤4-1);反之,说明交替方向乘子法已收敛,则区域a输出xa为各区域a的电力***状态变量,包括各区域a的电压幅值和相角,并结束计算。
本发明提出的一种基于双线性化的多区域电网全分布式抗差状态估计方法,其优点是:
(1)可扩展性:经过分区后,该方法各区域的子问题规模十分小,因此可以应对大规模***;
(2)可维护性:整体而言,该方法只需各区域与相邻区域的协同迭代,不需要控制中心进行协调、处理,不需要维护庞大的集中模型;
(3)私密性:各区域只需在第一阶段和第三阶段中与相邻区域交互边界信息,因此本区域内的数据隐私得到很好的保护。
附图说明
图1是本方法第一阶段分解示意图。
图2是本方法第三阶段分解示意图。
具体实施方式
本发明提出一种基于双线性化的多区域电网全分布式抗差状态估计方法,下面结合具体实施例进一步详细说明如下:
本发明提出一种基于双线性化的多区域电网全分布式抗差状态估计方法,其特征在于,该方法包括以下步骤:
1)将电力***分成多个控制区,一个控制区为一条或若干条母线和挂接在上面的发电机;假设给定的电力***被划分为R个区域,R为正整数;区域a、区域b为任意两个相邻区域,定义每个控制区中所有节点的集合为节点j表示该区域中任意一节点,mj:=|Mj|为包含节点j的区域数量;
2)对划分后的电力***的每个控制区域建立由目标函数和约束条件组成的分布式第一阶段线性模型,如图1所示,区域a和b为***中任意两个相邻的控制区域,区域a包含节点1、2、4等,其节点集Na={1,2,4…},区域b包含节点3、5等,其节点集Nb={3,5…};(以下具体步骤描述以区域a表示任意一个区域,每个区域均为相同处理)。
2-1)设任意控制区域a目标函数为最小化全网残差平方和以及保证坏数据稀疏性所引入的惩罚项如式(1):
其中,为区域a的目标函数,ya为区域a的中间状态变量,由{Ua,i,Ka,ij,La,ij}组成,La,ij=Va,iVa,jsinθa,ij,θa,ij=θa,i-θa,j,Va,i和Va,j分别为区域a节点i和节点j的电压幅值,θa,i和θa,j分别为区域a节点i和节点j的电压相角,为惩罚项;为区域a第一阶段坏数据向量,za为区域a的量测包括电压幅值量测、支路有功量测、支路无功量测、节点注入有功量测和注入无功量测,Ba为区域a的第一阶段量测矩阵,其元素由式(2)-(6)量测方程确定:
其中,和是区域a节点i的注入有功和无功量测,和是区域a支路ij的i端的有功和无功量测,是区域a节点i电压量测值的平方,gij为支路ij的电导,gsh,i为节点i的并联电导,bij为支路ij的电纳,bsh,i为节点i的并联电纳,εP,εQ,εU分别为有功量测、无功量测、电压幅值平方的误差项;λ为常数,一般取值为1.34e-2;
2-2)约束约束条件包括:区域a的零注入等式约束,如式(7)所示:
Ea代表区域a的零注入量测矩阵;
为了保证相邻区域边界状态量一致所引入的约束,如式(8)所示:
其中,Ka,ij为区域a线路ij两端电压幅值以及相角差之余弦的乘积,La,ij为区域a线路ij两端电压幅值以及相角差之正弦的乘积,Гa,b为区域a和b之间的联络线集合,Δa为a的相邻区域的集合;如图1所示,在联络线1-3上,为了保证相邻区域边界状态量一致所引入的约束为Ka,13=Kb,13,La,13=Lb,13。
3)令交替方向乘子法的迭代下标t=0;给定交替方向乘子法的收敛标准ε∈R+;采用交替方向乘子法对步骤2)建立的分布式第一阶段线性模型进行求解;具体包括:
3-1)更新控制区域a的中间状态变量,如式(9)所示:
ρf∈R+为交替方向乘子法的罚因子,一般取值为1,其中R+为正实数集,为增广信息矩阵,I为单位阵,为第t步迭代中区域a的辅助中间状态变量,由与{Ua,Ka,La}所对应的辅助中间状态变量组成,初始值选为为区域a的第一阶段辅助矩阵;
3-2)更新控制区域a的坏数据变量,如式(10)所示:
其中代表阈值算符,相当于对向量中的每个分量做如式(11)运算:
其中l代表分量序号;
3-3)区域a与相邻控制区b交互边界信息(辅助变量K和L),对相邻控制区b之间联络线上的辅助变量K和L求均值,如式(12):
3-4)更新区域a的辅助中间状态变量,如式(13):
其中,均为的分量;
3-5)判断第一阶段中交替方向乘子法是否收敛:
计算原始残差向量如式(14):
计算对偶残差向量如式(15):
计算总残差如式(16):
若δf,t≥εf,则t:=t+1,返回步骤3-1);反之,说明交替方向乘子法已收敛,进行步骤4);
4)对控制区域a的辅助中间状态变量进行第二阶段非线性变换,从而得到分布式第三阶段线性模型所需要的量测值,如式(17)-(19):
αa,i=lnUa,i (77)
其中αa,i为区域a节点i的电压幅值自然对数的两倍,αa,ij为区域a支路ij两端αa,i和αa,j之和,θa,ij为区域a的支路ij两端的电压相角差;
5)对控制区域a建立由目标函数和约束条件组成的分布式第三阶段线性模型,以图2为例,该***包含5个节点两个区域,区域a包含节点1、2、3、4,区域b包含节点3、5,节点3为两个区域共享的节点。
5-1)目标函数如下,为最小化全网残差平方和以及保证坏数据稀疏性所引入的惩罚项
其中为区域a的目标函数,xa为区域a的最终状态变量包括节点电压幅值和相角,为区域第二阶段坏数据向量,为惩罚项;为区域a的经过非线性转换后的量测,包括{αa,i,θa,ij,αa,ij},Ca为区域a的第三阶段量测矩阵,其元素由式(20)、(21)量测方程决定:
αa,ij=αa,i+αa,j+εα (80)
θa,ij=θa,i-θa,j+εθ (81)
其中,εα,εθ是第三阶段量测αa,ij,θa,ij的误差项;
5-2)约束条件是为了保证相邻区域在边界处一致的等式约束,如式(22):
其中,xa,i为区域a节点i的状态变量,为区域a的边界节点集;以图2为例,为了保证相邻区域在边界处一致,在节点3上所引入的等式约束为xa,3=xb,3.
6),令交替方向乘子法的迭代下标t=0;给定交替方向乘子法的收敛标准ε∈R+;采用交替方向乘子法对步骤5)的分布式第三阶段线性模型进行求解;
6-1)更新区域a的最终状态变量,如式(23)所示:
ρs∈R+为交替方向乘子法的罚因子,一般取为为第二阶段增广信息矩阵,为第t步迭代中区域a的辅助最终状态变量,初值为0;
6-2)更新区域a的坏数据变量,如式(24)所示:
6-3)控制区a利用本地和与邻居的通信,区域a与相邻区域b交互边界节点的状态变量最新估计值,对区域a边界节点上的状态变量求均值,设区域a的任一边界节点i,如式(25)所示:
其中,Mi为包含节点i的区域集合,mi为包含节点i的区域数;
6-4)更新区域a的辅助状态变量,如式(26)所示:
4-5)判断第三阶段中交替方向乘子法是否收敛,如式(27)所示:
计算原始残差向量如下:
计算对偶残差向量,如式(23)所示:
计算总残差,如式(23)所示:
若δs,t≥ε,则t:=t+1,返回步骤4-1);反之,说明交替方向乘子法已收敛,则区域a输出xa,
为区域a的电力***状态变量,包括区域a的电压幅值和相角,并结束计算。
Claims (1)
1.一种基于双线性化的多区域电网全分布式抗差状态估计方法,其特征在于,该方法包括以下步骤:
1)将电力***分成多个控制区,一个控制区为一条或若干条母线和挂接在上面的发电机;设给定的电力***被划分为R个区域,R为正整数;定义每个控制区中所有节点的集合为节点j表示该区域中任意一节点,mj:=|Mj|为包含节点j的区域数量;
2)对划分后的电力***的每个控制区域建立由目标函数和约束条件组成的分布式第一阶段线性模型,区域a和b为***中任意两个相邻的控制区域,区域a其节点集Na={1,2,4…},区域b节点集Nb={3,5…};
2-1)设任意控制区域a目标函数为最小化全网残差平方和以及保证坏数据稀疏性所引入的惩罚项如式(1):
其中,为区域a的目标函数,ya为区域a的中间状态变量,由{Ua,i,Ka,ij,La,ij}组成,Ka,ij=Va,iVa,jcosθa,ij,La,ij=Va,iVa,jsinθa,ij,θa,ij=θa,i-θa,j,Va,i和Va,j分别为区域a节点i和节点j的电压幅值,θa,i和θa,j分别为区域a节点i和节点j的电压相角,为惩罚项;为区域a第一阶段坏数据向量,za为区域a的量测包括电压幅值量测、支路有功量测、支路无功量测、节点注入有功量测和注入无功量测,Ba为区域a的第一阶段量测矩阵,其元素由式(2)-(6)量测方程确定:
其中,和是区域a节点i的注入有功和无功量测,和是区域a支路ij的i端的有功和无功量测,是区域a节点i电压量测值的平方,gij为支路ij的电导,gsh,i为节点i的并联电导,bij为支路ij的电纳,bsh,i为节点i的并联电纳,εP,εQ,εU分别为有功量测、无功量测、电压幅值平方的误差项;λ为常数;
2-2)约束约束条件包括:区域a的零注入等式约束,如式(7)所示:
Ea代表区域a的零注入量测矩阵;
为保证相邻区域边界状态量一致所引入的约束,如式(8)所示:
其中,Ka,ij为区域a线路ij两端电压幅值以及相角差之余弦的乘积,La,ij为区域a线路ij两端电压幅值以及相角差之正弦的乘积,Гa,b为区域a和b之间的联络线集合,Δa为a的相邻区域的集合;
3)令交替方向乘子法的迭代下标t=0;给定交替方向乘子法的收敛标准ε∈R+;采用交替方向乘子法对步骤2)建立的分布式第一阶段线性模型进行求解;具体包括:
3-1)更新控制区域a的中间状态变量,如式(9)所示:
ρf∈R+为交替方向乘子法的罚因子,其中R+为正实数集,为增广信息矩阵,I为单位阵,为第t步迭代中区域a的辅助中间状态变量,由与{Ua,Ka,La}所对应的辅助中间状态变量组成,初始值选为 为区域a的第一阶段辅助矩阵;
3-2)更新控制区域a的坏数据变量,如式(10)所示:
其中代表阈值算符,相当于对向量中的每个分量做如式(11)运算:
其中l代表分量序号;
3-3)控制区a与相邻控制区b交互边界辅助变量K和L,对相邻控制区a和b之间联络线上的辅助变量K和L求均值,如式(12):
3-4)更新区域a的辅助中间状态变量,如式(13):
其中,均为的分量;
3-5)判断第一阶段中交替方向乘子法是否收敛:
计算原始残差向量如式(14):
计算对偶残差向量如式(15):
计算总残差如式(16):
若δf,t≥εf,则t:=t+1,返回步骤3-1);反之,说明交替方向乘子法已收敛,进行步骤4);
4)对控制区域a的辅助中间状态变量进行第二阶段非线性变换,从而得到分布式第三阶段线性模型所需要的量测,如式(17)-(19):
αa,i=lnUa,i (17)
其中αa,i为区域a节点i的电压幅值自然对数的两倍,αa,ij为区域a支路ij两端αa,i和αa,j之和,θa,ij为区域a的支路ij两端的电压相角差;
5)对控制区a域建立由目标函数和约束条件组成的分布式第三阶段线性模型;
5-1)目标函数如下,为最小化全网残差平方和以及保证坏数据稀疏性所引入的惩罚项
其中为区域a的目标函数,xa为区域a的最终状态变量包括节点电压幅值和相角,为区域第二阶段坏数据向量,为惩罚项;为区域a的经过非线性转换后的量测,包括{αa,i,θa,ij,αa,ij},Ca为区域a的第三阶段量测矩阵,其元素由式(20)、(21)量测方程决定:
αa,ij=αa,i+αa,j+εα (20)
θa,ij=θa,i-θa,j+εθ (21)
其中,εα,εθ是第三阶段量测αa,ij,θa,ij的误差项;
5-2)约束条件是用以保证相邻区域在边界处一致的等式约束,如式(22):
s.t.xa,i=xb,i,
其中,xa,i为区域a节点i的状态变量,为区域a的边界节点集;
6)令交替方向乘子法的迭代下标t=0;给定交替方向乘子法的收敛标准ε∈R+;采用交替方向乘子法对步骤5)的分布式第三阶段线性模型进行求解;
6-1)更新区域a的最终状态变量,如式(23)所示:
ρs∈R+为交替方向乘子法的罚因子,为第二阶段增广信息矩阵,为第t步迭代中区域a的辅助最终状态变量,初值为0;
6-2)更新区域a的坏数据变量,如式(24)所示:
6-3)控制区a利用本地和与邻居的通信,区域a与相邻区域b交互边界节点的状态变量最新估计值,对区域a任意边界节点i上的状态变量求均值,如式(25)所示:
其中,Mi为包含节点i的区域集合,mi为包含节点i的区域数;
6-4)更新区域a的辅助状态变量,如式(26)所示:
4-5)判断第三阶段中交替方向乘子法是否收敛,计算原始残差向量如式(27)所示:
计算对偶残差向量,如式(28)、(29)所示:
计算总残差,如式(30)所示:
若δs,t≥ε,则t:=t+1,返回步骤4-1);反之,说明交替方向乘子法已收敛,则区域a输出xa为各区域a的电力***状态变量,包括各区域a的电压幅值和相角,并结束计算。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610067806.1A CN105552904B (zh) | 2016-01-30 | 2016-01-30 | 基于双线性化的多区域电网全分布式抗差状态估计方法 |
US15/331,876 US10222815B2 (en) | 2016-01-30 | 2016-10-23 | Method and device for estimating state of power system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610067806.1A CN105552904B (zh) | 2016-01-30 | 2016-01-30 | 基于双线性化的多区域电网全分布式抗差状态估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105552904A CN105552904A (zh) | 2016-05-04 |
CN105552904B true CN105552904B (zh) | 2018-02-02 |
Family
ID=55831918
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610067806.1A Active CN105552904B (zh) | 2016-01-30 | 2016-01-30 | 基于双线性化的多区域电网全分布式抗差状态估计方法 |
Country Status (2)
Country | Link |
---|---|
US (1) | US10222815B2 (zh) |
CN (1) | CN105552904B (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106887844B (zh) * | 2017-04-20 | 2019-05-21 | 清华大学 | 一种基于稀疏恢复的快速分解配电网三相状态估计方法 |
CN107508284B (zh) * | 2017-08-15 | 2020-05-19 | 华北电力大学 | 计及电气互联的微电网分布式优化调度方法 |
CN108767844B (zh) * | 2018-04-25 | 2021-06-04 | 上海大学 | 数据注入攻击下网络化多区域电力***的自适应状态估计方法 |
CN109255171B (zh) * | 2018-08-29 | 2023-09-05 | 深圳十沣科技有限公司 | 一种自动判定数值模拟计算收敛的方法 |
CN109120011B (zh) * | 2018-09-29 | 2019-12-13 | 清华大学 | 一种考虑分布式电源的分布式配电网拥塞调度方法 |
CN109301818B (zh) * | 2018-10-17 | 2021-07-13 | 国网河南省电力公司电力科学研究院 | 一种考虑调峰、调压需求的广域分布式储能***调度方法 |
CN109494711B (zh) * | 2018-10-24 | 2022-09-20 | 华北电力大学 | 一种多区域并行的全分布式状态估计方法 |
CN109462234B (zh) * | 2018-11-27 | 2022-05-24 | 国家电网有限公司 | 联络线扩展区域状态估计方法及装置 |
CN110391935B (zh) * | 2019-07-11 | 2022-04-15 | 南瑞集团有限公司 | 计及信息物理耦合特性和信息扰动的量测度评估方法及*** |
CN110601233B (zh) * | 2019-09-30 | 2023-02-21 | 国家电网公司西北分部 | 一种电力***中储能电站的调峰调度方法 |
CN110796385A (zh) * | 2019-11-04 | 2020-02-14 | 北京交通大学 | 电力***状态估计方法、装置、设备及存储介质 |
CN110969292A (zh) * | 2019-11-22 | 2020-04-07 | 全球能源互联网研究院有限公司 | 基于拓扑图的电力***状态测算方法、装置以及电子设备 |
CN112255457B (zh) * | 2020-09-22 | 2022-06-07 | 天津电气科学研究院有限公司 | 适用于自动准同期装置的相角差测量方法 |
CN112993989B (zh) * | 2021-03-05 | 2022-12-16 | 广东电网有限责任公司广州供电局 | 一种主动配电***抗差状态估计数据处理方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6690739B1 (en) * | 2000-01-14 | 2004-02-10 | Shou Yee Mui | Method for intersymbol interference compensation |
DE102008042356A1 (de) * | 2008-09-25 | 2010-04-08 | Carl Zeiss Smt Ag | Projektionsbelichtungsanlage mit optimierter Justagemöglichkeit |
WO2012061674A2 (en) * | 2010-11-04 | 2012-05-10 | Siemens Corporation | Stochastic state estimation for smart grids |
US9983670B2 (en) * | 2012-09-14 | 2018-05-29 | Interaxon Inc. | Systems and methods for collecting, analyzing, and sharing bio-signal and non-bio-signal data |
CN103336904A (zh) * | 2013-07-08 | 2013-10-02 | 国家电网公司 | 一种基于分段线性权因子函数的抗差状态估计方法 |
CN103413053B (zh) * | 2013-08-21 | 2016-09-14 | 国家电网公司 | 一种基于内点法的电力***抗差状态估计方法 |
CN104701849B (zh) * | 2015-03-02 | 2017-10-31 | 清华大学 | 一种主动配电网全分布式自律电压控制方法 |
-
2016
- 2016-01-30 CN CN201610067806.1A patent/CN105552904B/zh active Active
- 2016-10-23 US US15/331,876 patent/US10222815B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
US20170220521A1 (en) | 2017-08-03 |
US10222815B2 (en) | 2019-03-05 |
CN105552904A (zh) | 2016-05-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105552904B (zh) | 基于双线性化的多区域电网全分布式抗差状态估计方法 | |
Yang et al. | Distributed optimal dispatch of virtual power plant via limited communication | |
CN110826818B (zh) | 多个巡检员对多个站点进行巡检任务规划和路径设计方法 | |
Zhao et al. | Optimal placement of PMUs and communication links for distributed state estimation in distribution networks | |
Kong et al. | Three-stage distributed state estimation for AC-DC hybrid distribution network under mixed measurement environment | |
CN107301482A (zh) | 一种基于城市综合管廊的配电网中压线路路径规划方法 | |
CN104617577B (zh) | 一种主动配电网全分布式自律经济调度方法 | |
Li et al. | Selection of China’s imported grain distribution centers in the context of the Belt and Road initiative | |
CN106921156A (zh) | 一种基于多采样周期混合量测的主动配电网状态估计方法 | |
CN108695902A (zh) | 一种梯级水库群生态-发电动态互馈调控方法 | |
Mhanna et al. | Component-based dual decomposition methods for the OPF problem | |
CN103632313A (zh) | 一种基于pmu数据的电网动态可观方法 | |
CN106026086B (zh) | 一种电网中运行状态的动态估计方法 | |
CN106019077B (zh) | 一种电流型行波测距装置优化布点方法 | |
CN103077480A (zh) | 一种电力***的安全校核方法 | |
CN104701849A (zh) | 一种主动配电网全分布式自律电压控制方法 | |
CN109888773A (zh) | 一种电力***多区域分布式状态评估方法 | |
CN106155908B (zh) | 一种pas应用性能评估方法及*** | |
CN108376983A (zh) | 一种基于边界融合的电网状态估计***及方法 | |
CN104463437A (zh) | 一种基于多目标决策的电力通信接入网设计评价方法 | |
Hu et al. | Credibility-based distributed frequency estimation for plug-in electric vehicles participating in load frequency control | |
CN106443276A (zh) | 一种交流高压多回输电线路无线电干扰计算方法及*** | |
CN107919678A (zh) | 交直流混合配电网的供电恢复方法及*** | |
CN109494711A (zh) | 一种多区域并行的全分布式状态估计方法 | |
CN102142682A (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 |