CN102427227B - 考虑零注入约束的修正快速解耦电力系统状态估计方法 - Google Patents
考虑零注入约束的修正快速解耦电力系统状态估计方法 Download PDFInfo
- Publication number
- CN102427227B CN102427227B CN 201110316552 CN201110316552A CN102427227B CN 102427227 B CN102427227 B CN 102427227B CN 201110316552 CN201110316552 CN 201110316552 CN 201110316552 A CN201110316552 A CN 201110316552A CN 102427227 B CN102427227 B CN 102427227B
- Authority
- CN
- China
- Prior art keywords
- node
- zero
- vector
- phase angle
- iteration
- 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
- 238000002347 injection Methods 0.000 title claims abstract description 46
- 239000007924 injection Substances 0.000 title claims abstract description 46
- 238000012937 correction Methods 0.000 title claims abstract description 42
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000004364 calculation method Methods 0.000 claims abstract description 12
- 239000000243 solution Substances 0.000 claims abstract description 9
- 239000011159 matrix material Substances 0.000 claims description 30
- 238000005259 measurement Methods 0.000 claims description 21
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000000205 computational method Methods 0.000 claims description 3
- 230000003760 hair shine Effects 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 3
- 238000004088 simulation Methods 0.000 abstract description 2
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及一种考虑零注入约束的修正快速解耦电力系统状态估计方法,属于电力系统调度自动化与电力系统仿真技术领域。本发明方法包括:建立最小二乘状态估计模型,并按照通常的快速解耦估计对该状态估计模型进行迭代求解,在每一次迭代中,非零注入节点的节点电压幅值和相角按照快速解耦法计算出的电压幅值和相角的修正量进行修正,而零注入节点的节点电压幅值和相角不取通常快速解耦估计的计算结果,而是根据由零注入等式约束确立的零注入节点状态变量与非零注入节点状态变量之间的关系得到。本发明的整个计算过程与传统的状态估计计算流程非常相似,实现方便,同时能够保证零注入节点的注入功率严格为0,状态估计结果严格满足潮流方程。
Description
技术领域
本发明涉及一种考虑零注入约束的修正快速解耦电力系统状态估计方法,属于电力系统调度自动化与电力系统仿真技术领域。
背景技术
电力系统状态估计是电力系统能量管理系统的关键基础模块。在实际的电力系统中,存在许多既不挂接发电机也不挂接负荷的零注入节点。在电力系统状态估计的计算结果中,这些零注入节点的节点注入功率应当严格为0,否则,电力系统状态估计的计算结果将不能严格满足潮流方程,这将导致电力系统状态估计的计算结果与调度员潮流的计算结果有偏差,给电力系统的其他高级应用带来很大的不便。
目前处理零注入节点的一般做法是设置权重很大的零注入节点功率伪量测,以保证状态估计结果中零注入节点的注入功率较小。这是一种近似的方法,无法使得零注入节点的注入功率严格为0。实际上,目前国内外零注入节点注入功率较大的问题相当严重,研究简单高效的能够保证零注入节点的注入功率严格为0的电力系统状态估计求解方法非常重要。
发明内容
本发明的目的是提出一种考虑零注入约束的修正快速解耦电力系统状态估计方法,可以使用本发明提出的方法求解含零注入等式约束的电力系统状态估计模型,以保证状态估计结果完全满足潮流方程。
本发明提出的考虑零注入约束的修正快速解耦电力系统状态估计方法,包括以下步骤:
(1)建立一个含等式约束的电力状态估计模型:
s.t c(x)=0
等式约束为:使零注入节点的节点注入功率为0,以c(x)=0表示,其中x是电力系统的状态变量,包括零注入节点的电压幅值向量Vz、零注入节点的电压相角向量θz、非零注入节点的电压幅值向量Vn和非零注入节点的电压相角向量θn,zi是第i号量测的实时测量值,hi(x)是第i号量测的实时测量方程,m为量测个数;
(2)根据电力系统当前的拓扑结构和网络参数,形成电力系统当前的节点导纳矩阵,并计算以下系数矩阵F
其中矩阵Gzz和Bzz分别是节点导纳矩阵中零注入节点对应的对角子阵的实部和虚部,矩阵Gzn和Bzn分别是节点导纳矩阵中零注入节点与非零注入节点的交叉非对角子阵对应的实部和虚部;
(3)根据电力系统当前的拓扑结构和网络参数,形成快速解耦估计有功迭代的系数矩阵A和无功迭代的系数矩阵R,分别为:
其中Ha,Hr分别是电力系统状态估计的有功迭代和无功迭代中使用的定常量测雅可比矩阵,Wa,Wr分别为电力系统状态估计的有功量测和无功量测的权重矩阵,上标T表示转置;
(4)设置电力系统状态估计的节点电压幅值的初值,记为V(0),节点电压相角的初值为θ(0),并设置迭代次数,记为k=0;
(5)第k次迭代得到节点电压幅值向量V(k)和电压相角向量θ(k),将θ(k)中零注入节点的电压相角向量和非零注入节点的电压相角向量分别记为和将V(k)中零注入节点的电压幅值向量和非零注入节点的电压幅值向量分别记为和
其中矩阵F为步骤(2)中的计算结果,Φ为复数理论中用直角坐标表示的复数到用极坐标表示的复数的转换映射;Φ-1表示Φ的逆映射。是非零注入节点的状态变量,由非零注入节点的电压幅值向量和相角向量组成,计算得到的是零注入节点的状态变量向量,由零注入节点的电压幅值向量和相角向量组成,由向量就可以得到和
Φ的表达式是:
Φ-1的表达式是:
e=Vcosθ
f=Vsinθ
其中e,f是用直角坐标表示的节点电压的实部和虚部,V,θ是用极坐标表示的节点电压的幅值和相角;。
(8)求解如下的快速解耦估计有功迭代修正方程,得到第k次迭代节点相角向量的修正量Δθ(k)
AΔθ(k)=α(k)
(9)根据设定的收敛精度ε,对上述修正量Δθ(k)进行判断,若max|Δθ(k)|<ε,且无功迭代已收敛,则状态估计计算完成;若max|Δθ(k)|<ε,且无功迭代未收敛,则标记当前有功迭代已收敛,并进行步骤(12);若max|Δθ(k)|≥ε,则进行步骤(10);
(11)根据上述第k+1次迭代的电压相角变量利用步骤(6)的方法,重新计算零注入节点对应的电压幅值向量相角子向量,记为和
(13)求解以下的快速解耦估计无功迭代修正方程,得到第k次迭代节点电压幅值向量的修正量ΔV(k)
RΔV(k)=β(k);
(14)根据设定的收敛精度ε,对上述修正量ΔV(k)进行判断,若max|ΔV(k)|<ε,且有功迭代已收敛,则状态估计计算完成;若max|ΔV(k)|<ε,且有功迭代未收敛,则标记当前无功迭代已收敛,并进行步骤(15),若max|ΔV(k)|≥ε,则进行步骤(5);
本发明提出的一种考虑零注入约束的修正快速解耦电力系统状态估计方法,其优点是:
1、本发明方法可以保证状态估计结果完全满足潮流方程,计算结果没有功率不平衡量,状态估计结果与调度员潮流结果完全一致。
2、本发明方法计算速度与现有的大权重法状态估计程序相当,但大权重法状态估计的计算结果不能完全满足潮流方程。本发明方法的计算速度远远快于现有的其他能够保证状态估计结果严格满足潮流方程的方法。
3、本方法数值稳定性好于现有的任何状态估计解法,收敛非常可靠。
4、本发明方法与现在广泛使用的传统状态估计算法兼容性非常好,只需要很小的程序改动即可实现,实现起来很方便。
具体实施方式
本发明提出的考虑零注入约束的修正快速解耦电力系统状态估计方法,包括以下步骤:
(1)建立一个含等式约束的电力状态估计模型:
s.t c(x)=0
等式约束为:使零注入节点的节点注入功率为0,以c(x)=0表示,其中x是电力系统的状态变量,包括零注入节点的电压幅值向量Vz、零注入节点的电压相角向量θz、非零注入节点的电压幅值向量Vn和非零注入节点的电压相角向量θn,zi是第i号量测的实时测量值,hi(x)是第i号量测的实时测量方程,m为量测个数;
(2)根据电力系统当前的拓扑结构和网络参数,形成电力系统当前的节点导纳矩阵,并计算以下系数矩阵F
其中矩阵Gzz和Bzz分别是节点导纳矩阵中零注入节点对应的对角子阵的实部和虚部,矩阵Gzn和Bzn分别是节点导纳矩阵中零注入节点与非零注入节点的交叉非对角子阵对应的实部和虚部;
(3)根据电力系统当前的拓扑结构和网络参数,形成快速解耦估计有功迭代的系数矩阵A和无功迭代的系数矩阵R,分别为:
其中Ha,Hr分别是电力系统状态估计的有功迭代和无功迭代中使用的定常量测雅可比矩阵,Wa,Wr分别为电力系统状态估计的有功量测和无功量测的权重矩阵,上标T表示转置;
(4)设置电力系统状态估计的节点电压幅值的初值,记为V(0),节点电压相角的初值为θ(0),并设置迭代次数,记为k=0;
(5)第k次迭代得到节点电压幅值向量V(k)和电压相角向量θ(k),将θ(k)中零注入节点的电压相角向量和非零注入节点的电压相角向量分别记为和将V(k)中零注入节点的电压幅值向量和非零注入节点的电压幅值向量分别记为和
通过下式计算零注入节点的状态变量向量
其中矩阵F为步骤(2)中的计算结果,Φ为复数理论中用直角坐标表示的复数到用极坐标表示的复数的转换映射;Φ-1表示Φ的逆映射。是非零注入节点的状态变量,由非零注入节点的电压幅值向量和相角向量组成,计算得到的是零注入节点的状态变量向量,由零注入节点的电压幅值向量和相角向量组成,由向量就可以得到和
Φ的表达式是:
Φ-1的表达式是:
e=Vcosθ
f=Vsinθ
其中e,f是用直角坐标表示的节点电压的实部和虚部,V,θ是用极坐标表示的节点电压的幅值和相角;。
(8)求解如下的快速解耦估计有功迭代修正方程,得到第k次迭代节点相角向量的修正量Δθ(k)
AΔθ(k)=α(k)
(9)根据设定的收敛精度ε,对上述修正量Δθ(k)进行判断,若max|Δθ(k)|<ε,且无功迭代已收敛,则状态估计计算完成;若max|Δθ(k)|<ε,且无功迭代未收敛,则标记当前有功迭代已收敛,并进行步骤(12);若max|Δθ(k)|≥ε,则进行步骤(10);通常取ε=0.0001。
(13)求解以下的快速解耦估计无功迭代修正方程,得到第k次迭代节点电压幅值向量的修正量ΔV(k)
RΔV(k)=β(k);
(14)根据设定的收敛精度ε,对上述修正量ΔV(k)进行判断,若max|ΔV(k)|<ε,且有功迭代已收敛,则状态估计计算完成;若max|ΔV(k)|<ε,且有功迭代未收敛,则标记当前无功迭代已收敛,并进行步骤(15),若max|ΔV(k)|≥ε,则进行步骤(5);
Claims (1)
1.一种考虑零注入约束的修正快速解耦电力系统状态估计方法,其特征在于该方法包括以下步骤:
(1)建立一个含等式约束的电力状态估计模型:
等式约束为:使零注入节点的节点注入功率为0,以c(x)=0表示,其中x是电力系统的状态变量,包括零注入节点的电压幅值向量Vz、零注入节点的电压相角向量θz、非零注入节点的电压幅值向量Vn和非零注入节点的电压相角向量θn,zi是第i号量测的实时测量值,hi(x)是第i号量测的实时测量方程,m为量测个数;
(2)根据电力系统当前的拓扑结构和网络参数,形成电力系统当前的节点导纳矩阵,并计算以下系数矩阵F
其中矩阵Gzz和Bzz分别是节点导纳矩阵中零注入节点对应的对角子阵的实部和虚部,矩阵Gzn和Bzn分别是节点导纳矩阵中零注入节点与非零注入节点的交叉非对角子阵对应的实部和虚部;
(3)根据电力系统当前的拓扑结构和网络参数,形成快速解耦估计有功迭代的系数矩阵A和无功迭代的系数矩阵R,分别为:
其中Ha,Hr分别是电力系统状态估计的有功迭代和无功迭代中使用的定常量测雅可比矩阵,Wa,Wr分别为电力系统状态估计的有功量测和无功量测的权重矩阵,上标T表示转置;
(4)设置电力系统状态估计的节点电压幅值的初值,记为V(0),节点电压相角的初值为θ(0),并设置迭代次数,记为k=0;
(5)第k次迭代得到节点电压幅值向量V(k)和电压相角向量θ(k),将θ(k)中零注入节点的电压相角向量和非零注入节点的电压相角向量分别记为和将V(k)中零注入节点的电压幅值向量和非零注入节点的电压幅值向量分别记为和
其中矩阵F为步骤(2)中的计算结果,Φ为复数理论中用直角坐标表示的复数到用极坐标表示的复数的转换映射;Φ-1表示Φ的逆映射,是非零注入节点的状态变量,由非零注入节点的电压幅值向量和相角向量组成,计算得到的是零注入节点的状态变量向量,由零注入节点的电压幅值向量和相角向量组成,由向量就可以得到 和
Φ的表达式是:
Φ-1的表达式是:
e=V cosθ
f=V sinθ
其中e,f是用直角坐标表示的节点电压的实部和虚部,V,θ是用极坐标表示的节点电压的幅值和相角;
(8)求解如下的快速解耦估计有功迭代修正方程,得到第k次迭代节点相角向量的修正量Δθ(k)
AΔθ(k)=α(k);
(9)根据设定的收敛精度ε,对上述修正量Δθ(k)进行判断,若max|Δθ(k)|<ε,且无功迭代已收敛,则状态估计计算完成;若max|Δθ(k)|<ε,且无功迭代未收敛,则标记当前有功迭代已收敛,并进行步骤(12);若max|Δθ(k)|≥ε,则进行步骤(10);通常取ε=0.0001;
(13)求解以下的快速解耦估计无功迭代修正方程,得到第k次迭代节点电压幅值向量的修正量ΔV(k)
RΔV(k)=β(k);
(14)根据设定的收敛精度ε,对上述修正量ΔV(k)进行判断,若max|ΔV(k)|<ε,且有功迭代已收敛,则状态估计计算完成;若max|ΔV(k)|<ε,且有功迭代未收敛,则标记当前无功迭代已收敛,并进行步骤(15),若max|ΔV(k)|≥ε,则进行步骤(5);
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110316552 CN102427227B (zh) | 2011-10-18 | 2011-10-18 | 考虑零注入约束的修正快速解耦电力系统状态估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110316552 CN102427227B (zh) | 2011-10-18 | 2011-10-18 | 考虑零注入约束的修正快速解耦电力系统状态估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102427227A CN102427227A (zh) | 2012-04-25 |
CN102427227B true CN102427227B (zh) | 2013-08-21 |
Family
ID=45961177
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110316552 Active CN102427227B (zh) | 2011-10-18 | 2011-10-18 | 考虑零注入约束的修正快速解耦电力系统状态估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102427227B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102842908B (zh) * | 2012-09-11 | 2014-11-12 | 河海大学 | 含多变压器支路的配电网三相解耦潮流计算方法 |
CN105811394B (zh) * | 2014-12-30 | 2019-03-08 | 国家电网公司 | 一种用于电力系统潮流近似计算的电压估算方法 |
CN107579516B (zh) * | 2017-09-13 | 2020-08-25 | 国电南瑞科技股份有限公司 | 一种提高电力系统状态估计计算速度的方法 |
CN107994567B (zh) * | 2017-11-22 | 2020-09-01 | 华北电力大学 | 一种广义快速分解状态估计方法 |
CN107957974B (zh) * | 2017-11-29 | 2020-12-18 | 广西大学 | 一种提取极坐标潮流方程雅可比矩阵的方法 |
CN110021928B (zh) * | 2018-11-30 | 2023-05-23 | 南京南瑞继保电气有限公司 | 一种考虑温度影响的状态估计解耦方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6545482B1 (en) * | 2000-09-22 | 2003-04-08 | Nxtphase Technology Srl | Monitoring wide area dynamic swing recordings on an electrical power system |
CN101325336A (zh) * | 2008-07-25 | 2008-12-17 | 江西电力调度中心 | 一种基于带约束状态估计的多区域电网潮流匹配方法 |
CN102185308A (zh) * | 2010-03-19 | 2011-09-14 | 清华大学 | 一种考虑零注入量测等式约束的电力系统状态估计方法 |
-
2011
- 2011-10-18 CN CN 201110316552 patent/CN102427227B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6545482B1 (en) * | 2000-09-22 | 2003-04-08 | Nxtphase Technology Srl | Monitoring wide area dynamic swing recordings on an electrical power system |
CN101325336A (zh) * | 2008-07-25 | 2008-12-17 | 江西电力调度中心 | 一种基于带约束状态估计的多区域电网潮流匹配方法 |
CN102185308A (zh) * | 2010-03-19 | 2011-09-14 | 清华大学 | 一种考虑零注入量测等式约束的电力系统状态估计方法 |
Non-Patent Citations (1)
Title |
---|
兰华,李积捷.电力系统状态估计算法的研究现状和展望.《继电器》.2007,第35卷(第10期), * |
Also Published As
Publication number | Publication date |
---|---|
CN102427227A (zh) | 2012-04-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102427229B (zh) | 基于修正牛顿法的带零注入约束的电力系统状态估计方法 | |
CN102427227B (zh) | 考虑零注入约束的修正快速解耦电力系统状态估计方法 | |
CN103840452B (zh) | 一种引入pmu量测信息的大电网状态估计方法 | |
CN102831315B (zh) | 一种电力系统状态估计量测方程的精确线性化方法 | |
CN101505061B (zh) | 一种基于时域仿真的戴维南等值参数跟踪的计算方法 | |
CN105322546B (zh) | 交直流解耦混合潮流算法 | |
CN103532137B (zh) | 一种三相四线低压配电网的状态估计方法 | |
CN106099920A (zh) | 一种基于参数估计的现代输电网虚假数据攻击方法 | |
CN102185308B (zh) | 一种考虑零注入量测等式约束的电力系统状态估计方法 | |
CN106786493A (zh) | 一种多馈入直流相互作用因子的实用计算方法 | |
CN103413053A (zh) | 一种基于内点法的电力系统抗差状态估计方法 | |
CN105512502A (zh) | 一种基于残差归一化的权函数最小二乘状态估计方法 | |
CN103400303A (zh) | 基于内网常规Ward等值电路的外网静态等值方法 | |
CN104182644A (zh) | 一种融合分布式电源特性的配电网状态估计方法 | |
CN105071387A (zh) | 基于复数域标幺化理论的配电网快速分解状态估计方法 | |
CN106786536B (zh) | 考虑外网扩展Ward等值的状态估计方法 | |
CN108460547A (zh) | 热电联供系统热电联合状态估计方法及系统 | |
CN101788608B (zh) | 一种独立三绕组变压器电抗参数的估计方法 | |
CN105514978B (zh) | 一种混合整数二次规划形式的抗差状态估计方法 | |
CN103825270B (zh) | 一种配电网三相状态估计雅可比矩阵常数化的处理方法 | |
CN103001220B (zh) | 一种基于最优匹配回路功率的配电网拓扑结构辨识方法 | |
CN103838962A (zh) | 一种计及pmu量测的分步线性状态估计方法 | |
CN102545216B (zh) | 用于电力系统暂态稳定仿真的发电机节点电压预测方法 | |
CN103986158A (zh) | 一种分布式电源配电网潮流计算方法 | |
CN104022501B (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |