CN107016489A - 一种电力***抗差状态估计方法和装置 - Google Patents
一种电力***抗差状态估计方法和装置 Download PDFInfo
- Publication number
- CN107016489A CN107016489A CN201710137890.4A CN201710137890A CN107016489A CN 107016489 A CN107016489 A CN 107016489A CN 201710137890 A CN201710137890 A CN 201710137890A CN 107016489 A CN107016489 A CN 107016489A
- Authority
- CN
- China
- Prior art keywords
- node
- measurement
- power
- state variable
- power system
- 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
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000005259 measurement Methods 0.000 claims abstract description 226
- 239000011159 matrix material Substances 0.000 claims abstract description 97
- 238000002347 injection Methods 0.000 claims abstract description 44
- 239000007924 injection Substances 0.000 claims abstract description 44
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 21
- 238000004364 calculation method Methods 0.000 claims abstract description 17
- 238000001514 detection method Methods 0.000 claims abstract description 12
- 230000008859 change Effects 0.000 claims abstract description 5
- 238000012937 correction Methods 0.000 claims description 22
- 230000017105 transposition Effects 0.000 claims description 6
- 230000011218 segmentation Effects 0.000 abstract description 2
- 230000006870 function Effects 0.000 description 30
- 238000004590 computer program Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 238000012545 processing Methods 0.000 description 5
- 230000009466 transformation Effects 0.000 description 4
- 230000005611 electricity Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000005304 joining Methods 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- JEIPFZHSYJVQDO-UHFFFAOYSA-N ferric oxide Chemical compound O=[Fe]O[Fe]=O JEIPFZHSYJVQDO-UHFFFAOYSA-N 0.000 description 1
- 230000009931 harmful effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
Classifications
-
- 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
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
-
- 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
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Entrepreneurship & Innovation (AREA)
- Theoretical Computer Science (AREA)
- Development Economics (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Educational Administration (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Power Engineering (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明实施例提供一种基于分段非线性权函数的电力***抗差状态估计方法和装置,方法包括构建电力***的节点导纳矩阵、量测量量测方程以及零注入等式约束;计算电力***的雅可比矩阵和量测量残差,得到各量测量的权重;更新状态变量,并判断迭代算法是否达到收敛阈值,若是则得到最终的电力***的节点电压幅值和相位。本发明将不良数据检测与辨识与状态估计同时进行,在状态估计计算的同时,能够改变量测量权重,在迭代中不断降低不良数据权重,抑制不良数据带来的影响,有效的减小了计算量,提高工作效率;且直接将量测量的残差带入分段非线性权函数,避免标准化残差的计算,能够大大提升计算速度,节省计算时间,十分适合工程应用。
Description
技术领域
本发明涉及电力***状态估计技术,具体涉及一种基于分段非线性权函数的电力***抗差状态估计方法和装置。
背景技术
传统的电力***状态估计应用最为广泛的是加权最小二乘估计算法,其模型简单,计算方法方便,在量测误差服从理想的正态分布的假设下能够得到方差最小且无偏的估计结果。
然而在电力***量测中,除了随机误差和***误差外,还不可避免的存在一些与真值相差很大的不良数据,这些不良数据会对估计结果带来严重的不良影响。混在量测数据中的不良数据会使加权最小二乘法的估计结果失去可靠性和精确性。为提高在线状态估计结果的准确性和实用性,在实际应用中,通常是在状态估计程序后加入不良数据的检测和辨识环节,利用残差检验等方法去检验不良数据。然而,这种方式需要状态估计和不良数据检测辨识两个环节多次的循环检验,其工作量巨大,计算时间长,并且由于残差检验方法自身固有的缺陷,辨识可靠性不够高。
现有的抗差状态估计方法包括加权最小绝对值(Weighted least absoluteabsolute,WLAV)和基于变量代换内点法的加权最小绝对值抗差状态估计等,其求解过程复杂,计算量大,无法应用到大型电力***的实时估计中。
发明内容
为了克服上述现有技术的不足,本发明提供一种基于分段非线性权函数的电力***抗差状态估计方法和装置,首先采集网络参数和SCADA***参数,并根据采集的网络参数和SCADA***参数确定节点导纳矩阵、量测量量测方程以及零注入等式约束;然后以节点电压幅值和相位为状态变量计算雅可比矩阵和量测量残差,并通过分段非线性权函数得到各量测量的权重;最后通过牛顿迭代法更新状态变量,并进行收敛判断,进而得到状态估计结果。本发明能满足计算效率要求,又具备良好抗差性能。
为了实现上述发明目的,本发明采取如下技术方案:
本发明提供一种基于分段非线性权函数的电力***抗差状态估计方法,所述方法包括:
根据数据采集与监视控制***SCADA采集到的电力***的网络参数和的节点与支路运行参数,构建电力***的节点导纳矩阵、量测量量测方程以及零注入等式约束;
以电力***的节点电压幅值和相位作为状态变量,并根据构建的电力***的节点导纳矩阵、量测量量测方程以及零注入等式约束计算电力***的雅可比矩阵和量测量残差,并通过分段非线性权函数得到各量测量的权重;
通过迭代算法更新状态变量,并根据状态变量的修正量判断迭代算法是否达到收敛阈值,若是则得到最终的电力***的节点电压幅值和相位。
所述根据采集的网络参数和SCADA***参数确定节点导纳矩阵、量测方程以及零注入等式约束包括:
所述节点导纳矩阵表示为:
其中,Yii表示节点导纳矩阵的对角线元素,Yij表示节点导纳矩阵的非对角线元素,yij表示节点i和节点j之间的支路阻抗zij的倒数。
所述根据采集的网络参数和SCADA***参数确定节点导纳矩阵、量测量量测方程以及零注入等式约束包括:
所述量测量量测方程表示为:
z=h(x)+r
其中,z表示量测量量测值,r表示量测量残差,x表示状态变量,h(x)为量测量计算方程,其包括节点功率量测方程、线路支路功率量测方程和变压器支路功率量测方程;
其中,节点功率量测方程表示为:
其中,Pi为节点i注入的有功功率,Qi为节点i注入的无功功率,vi为节点i的电压幅值,vj为节点j的电压幅值,θi为为节点i的电压相角,θj为节点j的电压相角,θij为表示节点i和节点j的电压相角差,Gij为节点导纳矩阵中对应节点i和j之间的互电导,Bij为节点导纳矩阵中对应节点i和j之间的互电纳,N为节点总数;
所述线路支路功率量测方程表示为:
其中,Pij为线路支路节点侧的有功功率,Qij为线路支路节点侧的无功功率,yc为线路支路对地电纳;
所述变压器支路功率量测方程表示为:
其中,Pij′为变压器支路节点i侧的有功功率,Qij′为变压器支路节点i侧的无功功率,Pji′为变压器支路节点j侧的有功功率,Qji′为变压器支路节点j侧的无功功率,K为变压器非标准变比,bT为变压器标准侧的电纳,节点i侧为标准侧,变比为1;节点j侧为非标准侧,变比为K。
所述根据采集的网络参数和SCADA***参数确定节点导纳矩阵、量测方程以及零注入等式约束包括:
所述零注入等式约束表示为:
c(x)=0
其中,c(x)为零注入节点的量测方程,表示节点注入有功、无功功率为零。
所述以数据采集与监视控制***的节点电压的幅值和相位作为状态变量,并根据构建的电力***的节点导纳矩阵、量测量量测方程以及零注入等式约束计算电力***的雅可比矩阵和量测量残差包括:
雅可比矩阵包含m个量测量计算方程对n个状态变量的偏导,其表示为:
其中,H(x)表示雅可比矩阵,其为m×n的矩阵;hm(x)表示第m个量测量计算方程,xn表示第n个状态变量;
对于节点电压幅值,其雅可比元素表示为:
对于节点注入有功功率和无功功率,其雅可比元素表示为:
其中,Gii为节点i的自电导,Bii为节点i的自电纳;
对于线路支路有功功率和无功功率以及变压器支路有功功率和无功功率,其雅可比元素表示为:
用l表示量测量索引,l=1,2,…,m,于是第l个量测量的残差rl表示为:
rl=zl-hl(xk)
其中,zl表示第l个量测量量测值;xk表示第k个状态变量,k=1,2,…,n;hl(xk)表示第l个量测量计算方程。
所述分段非线性权函数表示为:
其中,w(rl)表示第l个量测量的分段非线性权函数,a为不良数据的检测门槛值;
通过w(rl)得到第l个量测量的权重Rl,其表示为:
其中,为第l个量测量的初始权重,且σl为第l个量测量的标准差。
所述通过迭代算法更新状态变量,并根据状态变量的修正量判断迭代算法是否达到收敛阈值,若是则得到最终的电力***的节点电压幅值和相位包括:
采用牛顿迭代法得到状态变量的修正量,有:
其中,t表示迭代次数,Δx(t)为第t次迭代时状态变量的修正量,T表示转置,x(t)为第t次迭代时的状态变量,x(t+1)为第t+1次迭代时的状态变量,h(x(t))为第t次迭代时的状态变量对应的量测量计算方程,H(x(t))为第t次迭代时状态变量对应的雅可比矩阵,R为量测量的权重矩阵,且R=[R1,R2,…,Rl,…,Rm]T;
根据得到的Δx(t)更新状态变量x;
进行收敛判断,若满足max(|Δx(t)|)<εx,则表明收敛,结束迭代计算,输出状态评估结果;否则使迭代计数加一,重新进行迭代计算,其中εx表示收敛阈值。
本发明还提供一种基于分段非线性权函数的电力***抗差状态估计装置,所述装置包括:
构建模块,用于根据数据采集与监视控制***SCADA采集到的电力***的网络参数和的节点与支路运行参数,构建电力***的节点导纳矩阵、量测量量测方程以及零注入等式约束;
计算模块,用于以电力***的节点电压幅值和相位作为状态变量,并根据构建的电力***的节点导纳矩阵、量测量量测方程以及零注入等式约束计算电力***的雅可比矩阵和量测量残差,并通过分段非线性权函数得到各量测量的权重;
判断模块,用于通过迭代算法更新状态变量,并根据状态变量的修正量判断迭代算法是否达到收敛阈值,若是则得到最终的电力***的节点电压幅值和相位。
所述构建模块具体用于:
所述节点导纳矩阵表示为:
其中,Yii表示节点导纳矩阵的对角线元素,Yij表示节点导纳矩阵的非对角线元素,yij表示节点i和节点j之间的支路阻抗zij的倒数。
所述构建模块具体用于:
所述量测量量测方程表示为:
z=h(x)+r
其中,z表示量测量量测值,r表示量测量残差,x表示状态变量,h(x)为量测量计算方程,其包括节点功率量测方程、线路支路功率量测方程和变压器支路功率量测方程;
其中,节点功率量测方程表示为:
其中,Pi为节点i注入的有功功率,Qi为节点i注入的无功功率,vi为节点i的电压幅值,vj为节点j的电压幅值,θi为为节点i的电压相角,θj为节点j的电压相角,θij为表示节点i和节点j的电压相角差,Gij为节点导纳矩阵中对应节点i和j之间的互电导,Bij为节点导纳矩阵中对应节点i和j之间的互电纳,N为节点总数;
所述线路支路功率量测方程表示为:
其中,Pij为线路支路节点侧的有功功率,Qij为线路支路节点侧的无功功率,yc为线路支路对地电纳;
所述变压器支路功率量测方程表示为:
其中,Pij′为变压器支路节点i侧的有功功率,Qij′为变压器支路节点i侧的无功功率,Pji′为变压器支路节点j侧的有功功率,Qji′为变压器支路节点j侧的无功功率,K为变压器非标准变比,bT为变压器标准侧的电纳,节点i侧为标准侧,变比为1;节点j侧为非标准侧,变比为K。
所述构建模块具体用于:
所述零注入等式约束表示为:
c(x)=0
其中,c(x)为零注入节点的量测方程,表示节点注入有功、无功功率为零。
所述计算模块具体用于:
雅可比矩阵包含m个量测量计算方程对n个状态变量的偏导,其表示为:
其中,H(x)表示雅可比矩阵,其为m×n的矩阵;hm(x)表示第m个量测量计算方程,xn表示第n个状态变量;
对于节点电压幅值,其雅可比元素表示为:
对于节点注入有功功率和无功功率,其雅可比元素表示为:
其中,Gii为节点i的自电导,Bii为节点i的自电纳;
对于线路支路有功功率和无功功率以及变压器支路有功功率和无功功率,其雅可比元素表示为:
用l表示量测量索引,l=1,2,…,m,于是第l个量测量的残差rl表示为:
rl=zl-hl(xk)
其中,zl表示第l个量测量量测值;xk表示第k个状态变量,k=1,2,…,n;hl(xk)表示第l个量测量计算方程;
所述分段非线性权函数表示为:
其中,w(rl)表示第l个量测量的分段非线性权函数,a为不良数据的检测门槛值;
通过w(rl)得到第l个量测量的权重Rl,其表示为:
其中,为第l个量测量的初始权重,且σl为第l个量测量的标准差。
所述判断模块具体用于:
采用牛顿迭代法得到状态变量的修正量,有:
其中,t表示迭代次数,Δx(t)为第t次迭代时状态变量的修正量,T表示转置,x(t)为第t次迭代时的状态变量,x(t+1)为第t+1次迭代时的状态变量,h(x(t))为第t次迭代时的状态变量对应的量测量计算方程,H(x(t))为第t次迭代时状态变量对应的雅可比矩阵,R为量测量的权重矩阵,且R=[R1,R2,…,Rl,…,Rm]T;
根据得到的Δx(t)更新状态变量x;
进行收敛判断,若满足max(|Δx(t)|)<εx,则表明收敛,结束迭代计算,输出状态评估结果;否则使迭代计数加一,重新进行迭代计算,其中εx表示收敛阈值。
与最接近的现有技术相比,本发明提供的技术方案具有以下有益效果:
本发明将不良数据检测与辨识与状态估计同时进行,在状态估计计算的同时,能够改变量测量权重,在迭代中不断降低不良数据权重,抑制不良数据带来的影响,能够有效的减小计算量,提高工作效率;
与现有技术的权函数赋值相比,本发明直接将量测量的残差带入分段非线性权函数,避免标准化残差的计算,能够大大提升计算速度,节省计算时间,十分适合工程应用;
本发明采用分段非线性权函数,对于残差在检测门槛值以内的量测量,维持其权重,而对于残差在检测门槛值以外的量测量,通过指数型权函数降低其权重,使其接近于零而不等于零,能够有效的排除不良数据过多时,矩阵中对应零元素过多而出现非奇异矩阵的情况,使得状态估计能够有效收敛,极大的改善不良数据与量测误差共存时的数值稳定性。
附图说明
图1是本发明实施例中基于分段非线性权函数的电力***抗差状态估计方法流程图。
具体实施方式
下面结合附图对本发明作进一步详细说明。
本发明实施例提供一种基于分段非线性权函数的电力***抗差状态估计方法,如图1所示,该方法具体过程如下:
S101:根据数据采集与监视控制***SCADA采集到的电力***的网络参数和的节点与支路运行参数(节点与支路运行参数包括节点电压幅值、节点注入功率和支路潮流信息),构建电力***的节点导纳矩阵、量测量量测方程以及零注入等式约束;
S102:以S101中电力***的节点电压幅值和相位作为状态变量,并根据构建的电力***的节点导纳矩阵、量测量量测方程以及零注入等式约束计算电力***的雅可比矩阵和量测量残差,并通过分段非线性权函数得到各量测量的权重;
S103:通过迭代算法更新S102中状态变量,并根据状态变量的修正量判断迭代算法是否达到收敛阈值,若是则得到最终的电力***的节点电压幅值和相位。
上述S101中,根据数据采集与监视控制***SCADA采集到的电力***的网络参数和的节点与支路运行参数,构建电力***的节点导纳矩阵、量测量量测方程以及零注入等式约束具体如下:
1)节点导纳矩阵表示为:
其中,Yii表示节点导纳矩阵的对角线元素,Yij表示节点导纳矩阵的非对角线元素,yij表示节点i和节点j之间的支路阻抗zij的倒数。
2)根据采集的网络参数和SCADA***参数确定节点导纳矩阵、量测量量测方程以及零注入等式约束包括:
量测量量测方程表示为:
z=h(x)+r
其中,z表示量测量量测值,r表示量测量残差,x表示状态变量,h(x)为量测量计算方程,其包括节点功率量测方程、线路支路功率量测方程和变压器支路功率量测方程;
其中,节点功率量测方程表示为:
其中,Pi为节点i注入的有功功率,Qi为节点i注入的无功功率,vi为节点i的电压幅值,vj为节点j的电压幅值,θi为为节点i的电压相角,θj为节点j的电压相角,θij为表示节点i和节点j的电压相角差,Gij为节点导纳矩阵中对应节点i和j之间的互电导,Bij为节点导纳矩阵中对应节点i和j之间的互电纳,N为节点总数;
线路支路功率量测方程表示为:
其中,Pij为线路支路节点侧的有功功率,Qij为线路支路节点侧的无功功率,yc为线路支路对地电纳;
变压器支路功率量测方程表示为:
其中,Pij′为变压器支路节点i侧的有功功率,Qij′为变压器支路节点i侧的无功功率,Pji′为变压器支路节点j侧的有功功率,Qji′为变压器支路节点j侧的无功功率,K为变压器非标准变比,bT为变压器标准侧的电纳,节点i侧为标准侧,变比为1;节点j侧为非标准侧,变比为K。
3)零注入等式约束表示为:
c(x)=0
其中,c(x)为零注入节点的量测方程,表示节点注入有功、无功功率为零。
状态变量包括节点电压幅值和节点电压相位,量测量包括节点电压幅值、节点注入有功功率、节点注入无功功率、线路支路有功功率、线路支路无功功率、变压器支路有功功率和变压器支路无功功率。
上述S102中,以S101中电力***的节点电压幅值和相位作为状态变量,并根据构建的电力***的节点导纳矩阵、量测量量测方程以及零注入等式约束计算电力***的雅可比矩阵和量测量残差,并通过分段非线性权函数得到各量测量的权重具体过程如下:
(1)雅可比矩阵包含m个量测量计算方程对n个状态变量的偏导,其表示为:
其中,H(x)表示雅可比矩阵,其为m×n的矩阵;hm(x)表示第m个量测量计算方程,xn表示第n个状态变量;
1)对于节点电压幅值,其雅可比元素表示为:
2)对于节点注入有功功率和无功功率,其雅可比元素表示为:
其中,Gii为节点i的自电导,Bii为节点i的自电纳;
3)对于线路支路有功功率和无功功率以及变压器支路有功功率和无功功率,其雅可比元素表示为:
(2)用l表示量测量索引,l=1,2,…,m,于是第l个量测量的残差rl表示为:
rl=zl-hl(xk)
其中,zl表示第l个量测量量测值;xk表示第k个状态变量,k=1,2,…,n;hl(xk)表示第l个量测量计算方程;
(3)分段非线性权函数表示为:
其中,w(rl)表示第l个量测量的分段非线性权函数,a为不良数据的检测门槛值;
(4)通过w(rl)得到第l个量测量的权重Rl,其表示为:
其中,为第l个量测量的初始权重,且σl为第l个量测量的标准差。
上述S103中,通过迭代算法更新S102中的状态变量,并根据状态变量的修正量判断迭代算法是否达到收敛阈值,若是则得到最终的电力***的节点电压幅值和相位具体包括过程如下:
(1)采用牛顿迭代法得到状态变量的修正量,有:
其中,t表示迭代次数,Δx(t)为第t次迭代时状态变量的修正量,T表示转置,x(t)为第t次迭代时的状态变量,x(t+1)为第t+1次迭代时的状态变量,h(x(t))为第t次迭代时的状态变量对应的量测量计算方程,H(x(t))为第t次迭代时状态变量对应的雅可比矩阵,R为量测量的权重矩阵,且R=[R1,R2,…,Rl,…,Rm]T;
(2)根据得到的Δx(t)更新状态变量x;
(3)进行收敛判断,若满足max(|Δx(t)|)<εx,则表明收敛,结束迭代计算,输出状态评估结果;否则使迭代计数加一,重新进行迭代计算,其中εx表示收敛阈值。
本发明实施例同时还提供一种基于分段非线性权函数的电力***抗差状态估计装置,该装置包括构建模块、计算模块和判断模块,具体功能如下:
构建模块,用于根据数据采集与监视控制***SCADA采集到的电力***的网络参数和的节点与支路运行参数,构建电力***的节点导纳矩阵、量测量量测方程以及零注入等式约束;
计算模块,用于以电力***的节点电压幅值和相位作为状态变量,并根据构建的电力***的节点导纳矩阵、量测量量测方程以及零注入等式约束计算电力***的雅可比矩阵和量测量残差,并通过分段非线性权函数得到各量测量的权重;
判断模块,用于通过迭代算法更新状态变量,并根据状态变量的修正量判断迭代算法是否达到收敛阈值,若是则得到最终的电力***的节点电压幅值和相位。
上述的构建模块确定的节点导纳矩阵表示为:
其中,Yii表示节点导纳矩阵的对角线元素,Yij表示节点导纳矩阵的非对角线元素,yij表示节点i和节点j之间的支路阻抗zij的倒数。
上述的构建模块确定的量测量量测方程表示为:
z=h(x)+r
其中,z表示量测量量测值,r表示量测量残差,x表示状态变量,h(x)为量测量计算方程,其包括节点功率量测方程、线路支路功率量测方程和变压器支路功率量测方程;
其中,节点功率量测方程表示为:
其中,Pi为节点i注入的有功功率,Qi为节点i注入的无功功率,vi为节点i的电压幅值,vj为节点j的电压幅值,θi为为节点i的电压相角,θj为节点j的电压相角,θij为表示节点i和节点j的电压相角差,Gij为节点导纳矩阵中对应节点i和j之间的互电导,Bij为节点导纳矩阵中对应节点i和j之间的互电纳,N为节点总数;
其中,线路支路功率量测方程表示为:
其中,Pij为线路支路节点侧的有功功率,Qij为线路支路节点侧的无功功率,yc为线路支路对地电纳;
其中,变压器支路功率量测方程表示为:
其中,Pij′为变压器支路节点i侧的有功功率,Qij′为变压器支路节点i侧的无功功率,Pji′为变压器支路节点j侧的有功功率,Qji′为变压器支路节点j侧的无功功率,K为变压器非标准变比,bT为变压器标准侧的电纳,节点i侧为标准侧,变比为1;节点j侧为非标准侧,变比为K。
上述的构建模块确定的零注入等式约束表示为:
c(x)=0
其中,c(x)为零注入节点的量测方程,表示节点注入有功、无功功率为零。
状态变量包括节点电压幅值和节点电压相位,所述量测量包括节点电压幅值、节点注入有功功率、节点注入无功功率、线路支路有功功率、线路支路无功功率、变压器支路有功功率和变压器支路无功功率。
上述的计算模块以电力***的节点电压幅值和相位作为状态变量,并根据构建的电力***的节点导纳矩阵、量测量量测方程以及零注入等式约束计算电力***的雅可比矩阵和量测量残差,并通过分段非线性权函数得到各量测量的权重具体过程如下:
雅可比矩阵包含m个量测量计算方程对n个状态变量的偏导,其表示为:
其中,H(x)表示雅可比矩阵,其为m×n的矩阵;hm(x)表示第m个量测量计算方程,xn表示第n个状态变量;
对于节点电压幅值,其雅可比元素表示为:
对于节点注入有功功率和无功功率,其雅可比元素表示为:
其中,Gii为节点i的自电导,Bii为节点i的自电纳;
对于线路支路有功功率和无功功率以及变压器支路有功功率和无功功率,其雅可比元素表示为:
用l表示量测量索引,l=1,2,…,m,于是第l个量测量的残差rl表示为:
rl=zl-hl(xk)
其中,zl表示第l个量测量量测值;xk表示第k个状态变量,k=1,2,…,n;hl(xk)表示第l个量测量计算方程;
分段非线性权函数表示为:
其中,w(rl)表示第l个量测量的分段非线性权函数,a为不良数据的检测门槛值;
通过w(rl)得到第l个量测量的权重Rl,其表示为:
其中,为第l个量测量的初始权重,且σl为第l个量测量的标准差。
上述的判断模块通过迭代算法更新状态变量,并根据状态变量的修正量判断迭代算法是否达到收敛阈值,若是则得到最终的电力***的节点电压幅值和相位具体过程如下:
(1)采用牛顿迭代法得到状态变量的修正量,有:
其中,t表示迭代次数,Δx(t)为第t次迭代时状态变量的修正量,T表示转置,x(t)为第t次迭代时的状态变量,x(t+1)为第t+1次迭代时的状态变量,h(x(t))为第t次迭代时的状态变量对应的量测量计算方程,H(x(t))为第t次迭代时状态变量对应的雅可比矩阵,R为量测量的权重矩阵,且R=[R1,R2,…,Rl,…,Rm]T;
(2)根据得到的Δx(t)更新状态变量x;
(3)进行收敛判断,若满足max(|Δx(t)|)<εx,则表明收敛,结束迭代计算,输出状态评估结果;否则使迭代计数加一,重新进行迭代计算,其中εx表示收敛阈值。
本发明实施例首先在迭代计算开始进行时,各量测量权值使用设定的初始值,利用加权最小二乘法求得初值;其次根据量测残差用分段非线性权函数计算量更新权重,进行量测权值随残差变化的迭代,每次迭代前修改量测的权值,通过计算迭代方程更新状态量,并与收敛性判据作比较,直至收敛,得到状态估计结果。
本领域内的技术人员应明白,本申请的实施例可提供为方法、***、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
Claims (14)
1.一种电力***抗差状态估计方法,其特征在于,所述方法包括:
根据数据采集与监视控制***SCADA采集到的电力***的网络参数和的节点与支路运行参数,构建电力***的节点导纳矩阵、量测量量测方程以及零注入等式约束;
以电力***的节点电压幅值和相位作为状态变量,并根据构建的电力***的节点导纳矩阵、量测量量测方程以及零注入等式约束计算电力***的雅可比矩阵和量测量残差,并通过分段非线性权函数得到各量测量的权重;
通过迭代算法更新状态变量,并根据状态变量的修正量判断迭代算法是否达到收敛阈值,若是则得到最终的电力***的节点电压幅值和相位。
2.根据权利要求1所述的电力***抗差状态估计方法,其特征在于,所述节点导纳矩阵表示为:
其中,Yii表示节点导纳矩阵的对角线元素,Yij表示节点导纳矩阵的非对角线元素,yij表示节点i和节点j之间的支路阻抗zij的倒数。
3.根据权利要求1所述的电力***抗差状态估计方法,其特征在于,所述量测量量测方程表示为:
z=h(x)+r
其中,z表示量测量量测值,r表示量测量残差,x表示状态变量,h(x)为量测量计算方程,其包括节点功率量测方程、线路支路功率量测方程和变压器支路功率量测方程;
其中,节点功率量测方程表示为:
其中,Pi为节点i注入的有功功率,Qi为节点i注入的无功功率,vi为节点i的电压幅值,vj为节点j的电压幅值,θi为为节点i的电压相角,θj为节点j的电压相角,θij为表示节点i和节点j的电压相角差,Gij为节点导纳矩阵中对应节点i和j之间的互电导,Bij为节点导纳矩阵中对应节点i和j之间的互电纳,N为节点总数;
所述线路支路功率量测方程表示为:
其中,Pij为线路支路节点侧的有功功率,Qij为线路支路节点侧的无功功率,yc为线路支路对地电纳;
所述变压器支路功率量测方程表示为:
其中,Pij′为变压器支路节点i侧的有功功率,Qij′为变压器支路节点i侧的无功功率,Pji′为变压器支路节点j侧的有功功率,Qji′为变压器支路节点j侧的无功功率,K为变压器非标准变比,bT为变压器标准侧的电纳,节点i侧为标准侧,变比为1;节点j侧为非标准侧,变比为K。
4.根据权利要求1所述的电力***抗差状态估计方法,其特征在于,所述零注入等式约束表示为:
c(x)=0
其中,c(x)为零注入节点的量测方程,表示节点注入有功、无功功率为零。
5.根据权利要求3所述的电力***抗差状态估计方法,其特征在于,所述以数据采集与监视控制***的节点电压的幅值和相位作为状态变量,并根据构建的电力***的节点导纳矩阵、量测量量测方程以及零注入等式约束计算电力***的雅可比矩阵和量测量残差包括:
雅可比矩阵包含m个量测量计算方程对n个状态变量的偏导,其表示为:
其中,H(x)表示雅可比矩阵,其为m×n的矩阵;hm(x)表示第m个量测量计算方程,xn表示第n个状态变量;
对于节点电压幅值,其雅可比元素表示为:
对于节点注入有功功率和无功功率,其雅可比元素表示为:
其中,Gii为节点i的自电导,Bii为节点i的自电纳;
对于线路支路有功功率和无功功率以及变压器支路有功功率和无功功率,其雅可比元素表示为:
用l表示量测量索引,l=1,2,…,m,于是第l个量测量的残差rl表示为:
rl=zl-hl(xk)
其中,zl表示第l个量测量量测值;xk表示第k个状态变量,k=1,2,…,n;hl(xk)表示第l个量测量计算方程。
6.根据权利要求5所述的电力***抗差状态估计方法,其特征在于,所述通过分段非线性权函数得到各量测量的权重包括:
所述分段非线性权函数表示为:
其中,w(rl)表示第l个量测量的分段非线性权函数,a为不良数据的检测门槛值;
通过w(rl)得到第l个量测量的权重Rl,其表示为:
其中,为第l个量测量的初始权重,且σl为第l个量测量的标准差。
7.根据权利要求6所述的电力***抗差状态估计方法,其特征在于,所述通过迭代算法更新状态变量,并根据状态变量的修正量判断迭代算法是否达到收敛阈值,若是则得到最终的电力***的节点电压幅值和相位包括:
采用牛顿迭代法得到状态变量的修正量,有:
其中,t表示迭代次数,Δx(t)为第t次迭代时状态变量的修正量,T表示转置,x(t)为第t次迭代时的状态变量,x(t+1)为第t+1次迭代时的状态变量,h(x(t))为第t次迭代时的状态变量对应的量测量计算方程,H(x(t))为第t次迭代时状态变量对应的雅可比矩阵,R为量测量的权重矩阵,且R=[R1,R2,…,Rl,…,Rm]T;
根据得到的Δx(t)更新状态变量x;
进行收敛判断,若满足max(|Δx(t)|)<εx,则表明收敛,结束迭代计算,输出状态评估结果;否则使迭代计数加一,重新进行迭代计算,其中εx表示收敛阈值。
8.一种电力***抗差状态估计装置,其特征在于,所述装置包括:
构建模块,用于根据数据采集与监视控制***SCADA采集到的电力***的网络参数和的节点与支路运行参数,构建电力***的节点导纳矩阵、量测量量测方程以及零注入等式约束;
计算模块,用于以电力***的节点电压幅值和相位作为状态变量,并根据构建的电力***的节点导纳矩阵、量测量量测方程以及零注入等式约束计算电力***的雅可比矩阵和量测量残差,并通过分段非线性权函数得到各量测量的权重;
判断模块,用于通过迭代算法更新状态变量,并根据状态变量的修正量判断迭代算法是否达到收敛阈值,若是则得到最终的电力***的节点电压幅值和相位。
9.根据权利要求8所述的电力***抗差状态估计装置,其特征在于,所述构建模块具体用于:
所述节点导纳矩阵表示为:
其中,Yii表示节点导纳矩阵的对角线元素,Yij表示节点导纳矩阵的非对角线元素,yij表示节点i和节点j之间的支路阻抗zij的倒数。
10.根据权利要求8所述的电力***抗差状态估计装置,其特征在于,所述构建模块具体用于:
所述量测量量测方程表示为:
z=h(x)+r
其中,z表示量测量量测值,r表示量测量残差,x表示状态变量,h(x)为量测量计算方程,其包括节点功率量测方程、线路支路功率量测方程和变压器支路功率量测方程;
其中,节点功率量测方程表示为:
其中,Pi为节点i注入的有功功率,Qi为节点i注入的无功功率,vi为节点i的电压幅值,vj为节点j的电压幅值,θi为为节点i的电压相角,θj为节点j的电压相角,θij为表示节点i和节点j的电压相角差,Gij为节点导纳矩阵中对应节点i和j之间的互电导,Bij为节点导纳矩阵中对应节点i和j之间的互电纳,N为节点总数;
所述线路支路功率量测方程表示为:
其中,Pij为线路支路节点侧的有功功率,Qij为线路支路节点侧的无功功率,yc为线路支路对地电纳;
所述变压器支路功率量测方程表示为:
其中,Pij′为变压器支路节点i侧的有功功率,Qij′为变压器支路节点i侧的无功功率,Pji′为变压器支路节点j侧的有功功率,Qji′为变压器支路节点j侧的无功功率,K为变压器非标准变比,bT为变压器标准侧的电纳,节点i侧为标准侧,变比为1;节点j侧为非标准侧,变比为K。
11.根据权利要求8所述的电力***抗差状态估计装置,其特征在于,所述构建模块具体用于:
所述零注入等式约束表示为:
c(x)=0
其中,c(x)为零注入节点的量测方程,表示节点注入有功、无功功率为零。
12.根据权利要求10所述的电力***抗差状态估计装置,其特征在于,所述计算模块具体用于:
雅可比矩阵包含m个量测量计算方程对n个状态变量的偏导,其表示为:
其中,H(x)表示雅可比矩阵,其为m×n的矩阵;hm(x)表示第m个量测量计算方程,xn表示第n个状态变量;
对于节点电压幅值,其雅可比元素表示为:
对于节点注入有功功率和无功功率,其雅可比元素表示为:
其中,Gii为节点i的自电导,Bii为节点i的自电纳;
对于线路支路有功功率和无功功率以及变压器支路有功功率和无功功率,其雅可比元素表示为:
用l表示量测量索引,l=1,2,…,m,于是第l个量测量的残差rl表示为:
rl=zl-hl(xk)
其中,zl表示第l个量测量量测值;xk表示第k个状态变量,k=1,2,…,n;hl(xk)表示第l个量测量计算方程。
13.根据权利要求12所述的电力***抗差状态估计装置,其特征在于,所述分段非线性权函数表示为:
其中,w(rl)表示第l个量测量的分段非线性权函数,a为不良数据的检测门槛值;
通过w(rl)得到第l个量测量的权重Rl,其表示为:
其中,为第l个量测量的初始权重,且σl为第l个量测量的标准差。
14.根据权利要求13所述的电力***抗差状态估计装置,其特征在于,所述判断模块具体用于:
采用牛顿迭代法得到状态变量的修正量,有:
其中,t表示迭代次数,Δx(t)为第t次迭代时状态变量的修正量,T表示转置,x(t)为第t次迭代时的状态变量,x(t+1)为第t+1次迭代时的状态变量,h(x(t))为第t次迭代时的状态变量对应的量测量计算方程,H(x(t))为第t次迭代时状态变量对应的雅可比矩阵,R为量测量的权重矩阵,且R=[R1,R2,…,Rl,…,Rm]T;
根据得到的Δx(t)更新状态变量x;
进行收敛判断,若满足max(|Δx(t)|)<εx,则表明收敛,结束迭代计算,输出状态评估结果;否则使迭代计数加一,重新进行迭代计算,其中εx表示收敛阈值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710137890.4A CN107016489A (zh) | 2017-03-09 | 2017-03-09 | 一种电力***抗差状态估计方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710137890.4A CN107016489A (zh) | 2017-03-09 | 2017-03-09 | 一种电力***抗差状态估计方法和装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107016489A true CN107016489A (zh) | 2017-08-04 |
Family
ID=59439714
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710137890.4A Pending CN107016489A (zh) | 2017-03-09 | 2017-03-09 | 一种电力***抗差状态估计方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107016489A (zh) |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107994567A (zh) * | 2017-11-22 | 2018-05-04 | 华北电力大学 | 一种广义快速分解状态估计方法 |
CN108281967A (zh) * | 2018-01-04 | 2018-07-13 | 国网浙江省电力有限公司电力科学研究院 | 一种改善电力***潮流计算收敛性的算法启动方法 |
CN108879695A (zh) * | 2018-07-13 | 2018-11-23 | 海南电网有限责任公司电力科学研究院 | 一种直流潮流模型的建立与计算方法 |
CN109888773A (zh) * | 2019-02-25 | 2019-06-14 | 武汉大学 | 一种电力***多区域分布式状态评估方法 |
CN110232199A (zh) * | 2018-11-30 | 2019-09-13 | 南京南瑞继保电气有限公司 | 一种考虑矩阵尺度伸缩的双线性wlav抗差状态估计方法 |
CN110417009A (zh) * | 2019-07-29 | 2019-11-05 | 天津大学 | 基于不同采样周期数据的配电网混合抗差状态估计方法 |
CN110619105A (zh) * | 2019-10-14 | 2019-12-27 | 国网山东省电力公司潍坊供电公司 | 一种基于量测量和热平衡方程的输电线路温度估计方法 |
CN110781450A (zh) * | 2019-09-17 | 2020-02-11 | 广西电网有限责任公司电力科学研究院 | 一种配电馈线电压测量数据的不良数据检测***及方法 |
CN111711186A (zh) * | 2020-05-28 | 2020-09-25 | 中国南方电网有限责任公司 | 基于gpu并行计算的电力***pq分解状态估计方法 |
CN112383048A (zh) * | 2020-10-20 | 2021-02-19 | 重庆大学 | 考虑公专用配电变压器测量特点的中压馈线电网三相状态估计方法 |
CN112398117A (zh) * | 2020-09-24 | 2021-02-23 | 北京航空航天大学 | 一种引起线路负荷过载的虚假数据注入攻击构造及防御方法 |
CN112636325A (zh) * | 2020-11-11 | 2021-04-09 | 广东电网有限责任公司广州供电局 | 一种配电网抗差状态计算方法、装置、终端及存储介质 |
CN112865096A (zh) * | 2021-03-18 | 2021-05-28 | 国网上海市电力公司 | 一种考虑pmu量测相角偏差的配电网状态估计方法及*** |
CN112993989A (zh) * | 2021-03-05 | 2021-06-18 | 广东电网有限责任公司广州供电局 | 一种主动配电***抗差状态估计数据处理方法 |
CN113032718A (zh) * | 2021-03-29 | 2021-06-25 | 北京华大九天科技股份有限公司 | 解决电路仿真中牛顿迭代算法死循环的方法 |
CN113139288A (zh) * | 2021-04-21 | 2021-07-20 | 广东电网有限责任公司广州供电局 | 基于gpu二次加速的配电网分布式抗差状态估计方法及装置 |
CN113255959A (zh) * | 2021-04-09 | 2021-08-13 | 广东电网有限责任公司电力调度控制中心 | 一种电力***动态状态估计方法及*** |
CN113313339A (zh) * | 2021-03-26 | 2021-08-27 | 贵州大学 | 一种配电网动态状态估计算法 |
CN115470736A (zh) * | 2022-09-29 | 2022-12-13 | 长沙学院 | 适应储能电站变工况运行的电力***动态行为建模方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101219545B1 (ko) * | 2011-09-14 | 2013-01-09 | 주식회사 파워이십일 | 전력계통에서의 최적화 기법을 적용한 파라미터 추정 방법 |
CN104765962A (zh) * | 2015-04-08 | 2015-07-08 | 河海大学 | 一种计及温度变化的电力***状态估计方法 |
CN105512502A (zh) * | 2016-01-13 | 2016-04-20 | 重庆大学 | 一种基于残差归一化的权函数最小二乘状态估计方法 |
CN105514977A (zh) * | 2015-11-27 | 2016-04-20 | 华北电力大学 | 一种电力***状态的双曲余弦型抗差状态估计方法 |
WO2016078477A1 (zh) * | 2014-11-18 | 2016-05-26 | 国电南瑞科技股份有限公司 | 一种变电站三相线性广义状态估计方法 |
-
2017
- 2017-03-09 CN CN201710137890.4A patent/CN107016489A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101219545B1 (ko) * | 2011-09-14 | 2013-01-09 | 주식회사 파워이십일 | 전력계통에서의 최적화 기법을 적용한 파라미터 추정 방법 |
WO2016078477A1 (zh) * | 2014-11-18 | 2016-05-26 | 国电南瑞科技股份有限公司 | 一种变电站三相线性广义状态估计方法 |
CN104765962A (zh) * | 2015-04-08 | 2015-07-08 | 河海大学 | 一种计及温度变化的电力***状态估计方法 |
CN105514977A (zh) * | 2015-11-27 | 2016-04-20 | 华北电力大学 | 一种电力***状态的双曲余弦型抗差状态估计方法 |
CN105512502A (zh) * | 2016-01-13 | 2016-04-20 | 重庆大学 | 一种基于残差归一化的权函数最小二乘状态估计方法 |
Non-Patent Citations (1)
Title |
---|
闫丽梅等: "状态估计与不良数据检测方法及其在大庆油田电网中的应用", 《化工自动化及仪表》 * |
Cited By (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107994567A (zh) * | 2017-11-22 | 2018-05-04 | 华北电力大学 | 一种广义快速分解状态估计方法 |
CN107994567B (zh) * | 2017-11-22 | 2020-09-01 | 华北电力大学 | 一种广义快速分解状态估计方法 |
CN108281967A (zh) * | 2018-01-04 | 2018-07-13 | 国网浙江省电力有限公司电力科学研究院 | 一种改善电力***潮流计算收敛性的算法启动方法 |
CN108879695A (zh) * | 2018-07-13 | 2018-11-23 | 海南电网有限责任公司电力科学研究院 | 一种直流潮流模型的建立与计算方法 |
CN110232199A (zh) * | 2018-11-30 | 2019-09-13 | 南京南瑞继保电气有限公司 | 一种考虑矩阵尺度伸缩的双线性wlav抗差状态估计方法 |
CN110232199B (zh) * | 2018-11-30 | 2023-12-19 | 南京南瑞继保电气有限公司 | 一种考虑矩阵尺度伸缩的双线性wlav抗差状态估计方法 |
CN109888773B (zh) * | 2019-02-25 | 2022-08-05 | 武汉大学 | 一种电力***多区域分布式状态评估方法 |
CN109888773A (zh) * | 2019-02-25 | 2019-06-14 | 武汉大学 | 一种电力***多区域分布式状态评估方法 |
CN110417009A (zh) * | 2019-07-29 | 2019-11-05 | 天津大学 | 基于不同采样周期数据的配电网混合抗差状态估计方法 |
CN110781450A (zh) * | 2019-09-17 | 2020-02-11 | 广西电网有限责任公司电力科学研究院 | 一种配电馈线电压测量数据的不良数据检测***及方法 |
CN110781450B (zh) * | 2019-09-17 | 2022-12-16 | 广西电网有限责任公司电力科学研究院 | 一种配电馈线电压测量数据的不良数据检测方法及*** |
CN110619105A (zh) * | 2019-10-14 | 2019-12-27 | 国网山东省电力公司潍坊供电公司 | 一种基于量测量和热平衡方程的输电线路温度估计方法 |
CN111711186B (zh) * | 2020-05-28 | 2023-05-02 | 中国南方电网有限责任公司 | 基于gpu并行计算的电力***pq分解状态估计方法 |
CN111711186A (zh) * | 2020-05-28 | 2020-09-25 | 中国南方电网有限责任公司 | 基于gpu并行计算的电力***pq分解状态估计方法 |
CN112398117A (zh) * | 2020-09-24 | 2021-02-23 | 北京航空航天大学 | 一种引起线路负荷过载的虚假数据注入攻击构造及防御方法 |
CN112383048A (zh) * | 2020-10-20 | 2021-02-19 | 重庆大学 | 考虑公专用配电变压器测量特点的中压馈线电网三相状态估计方法 |
CN112636325A (zh) * | 2020-11-11 | 2021-04-09 | 广东电网有限责任公司广州供电局 | 一种配电网抗差状态计算方法、装置、终端及存储介质 |
CN112993989B (zh) * | 2021-03-05 | 2022-12-16 | 广东电网有限责任公司广州供电局 | 一种主动配电***抗差状态估计数据处理方法 |
CN112993989A (zh) * | 2021-03-05 | 2021-06-18 | 广东电网有限责任公司广州供电局 | 一种主动配电***抗差状态估计数据处理方法 |
CN112865096A (zh) * | 2021-03-18 | 2021-05-28 | 国网上海市电力公司 | 一种考虑pmu量测相角偏差的配电网状态估计方法及*** |
CN113313339A (zh) * | 2021-03-26 | 2021-08-27 | 贵州大学 | 一种配电网动态状态估计算法 |
CN113032718B (zh) * | 2021-03-29 | 2022-05-24 | 北京华大九天科技股份有限公司 | 解决电路仿真中牛顿迭代算法死循环的方法及装置 |
CN113032718A (zh) * | 2021-03-29 | 2021-06-25 | 北京华大九天科技股份有限公司 | 解决电路仿真中牛顿迭代算法死循环的方法 |
CN113255959A (zh) * | 2021-04-09 | 2021-08-13 | 广东电网有限责任公司电力调度控制中心 | 一种电力***动态状态估计方法及*** |
CN113139288A (zh) * | 2021-04-21 | 2021-07-20 | 广东电网有限责任公司广州供电局 | 基于gpu二次加速的配电网分布式抗差状态估计方法及装置 |
CN115470736A (zh) * | 2022-09-29 | 2022-12-13 | 长沙学院 | 适应储能电站变工况运行的电力***动态行为建模方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107016489A (zh) | 一种电力***抗差状态估计方法和装置 | |
CN104134999B (zh) | 基于多数据源的配电网量测有效性分析实用化计算方法 | |
CN104836223B (zh) | 电网参数错误与不良数据协同辨识与估计方法 | |
CN103944165B (zh) | 一种大电网参数辨识估计方法 | |
CN103532137B (zh) | 一种三相四线低压配电网的状态估计方法 | |
CN104184144B (zh) | 一种用于多电压等级电网模型的抗差状态估计方法 | |
CN107577870A (zh) | 基于同步相量量测的配电网电压功率灵敏度鲁棒估计方法 | |
CN103972884A (zh) | 一种电力***状态估计方法 | |
CN103886193B (zh) | 一种电力***模糊自适应抗差估计方法 | |
CN110417011A (zh) | 一种基于互信息与迭代随机森林的在线动态安全评估方法 | |
CN102403720A (zh) | 一种基于暂态电压安全裕度的超实时重合时序整定方法 | |
CN106228459A (zh) | 基于蒙特卡洛的等值可靠性评估方法 | |
CN105512502A (zh) | 一种基于残差归一化的权函数最小二乘状态估计方法 | |
CN103473602A (zh) | 一种电网理论线损数据预测***及预测方法 | |
CN108565852A (zh) | 一种三阶段递进的大电网电压稳定评估的故障筛选与排序方法 | |
CN103324858A (zh) | 配电网三相潮流状态估计方法 | |
CN106356840A (zh) | 基于同步相量量测的地区电力***状态估计方法及*** | |
CN104182644A (zh) | 一种融合分布式电源特性的配电网状态估计方法 | |
CN109829246B (zh) | 一种基于参数综合可疑度的线路参数辨识方法 | |
CN104866714A (zh) | 一种电力***自适应核密度抗差状态估计方法 | |
CN104617575A (zh) | 一种用于静态安全分析预想故障下潮流不可解评估方法 | |
CN114204560A (zh) | 一种中压配电网线路参数识别方法 | |
CN101788608A (zh) | 一种独立三绕组变压器电抗参数的估计方法 | |
CN111756031B (zh) | 一种电网运行趋势估算方法和*** | |
CN105977963B (zh) | 一种基于两阶段法的主配网一体化分布式状态估计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
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: 20170804 |