CN104184144B - 一种用于多电压等级电网模型的抗差状态估计方法 - Google Patents
一种用于多电压等级电网模型的抗差状态估计方法 Download PDFInfo
- Publication number
- CN104184144B CN104184144B CN201410450017.7A CN201410450017A CN104184144B CN 104184144 B CN104184144 B CN 104184144B CN 201410450017 A CN201410450017 A CN 201410450017A CN 104184144 B CN104184144 B CN 104184144B
- Authority
- CN
- China
- Prior art keywords
- measurement
- matrix
- sigma
- weight
- state estimation
- 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 73
- 238000005259 measurement Methods 0.000 claims abstract description 168
- 239000011159 matrix material Substances 0.000 claims abstract description 118
- 238000012937 correction Methods 0.000 claims abstract description 10
- 230000009466 transformation Effects 0.000 claims abstract description 5
- 238000004364 calculation method Methods 0.000 claims description 44
- 230000008569 process Effects 0.000 claims description 37
- 230000008030 elimination Effects 0.000 claims description 23
- 238000003379 elimination reaction Methods 0.000 claims description 23
- 239000013585 weight reducing agent Substances 0.000 claims description 16
- 238000004458 analytical method Methods 0.000 claims description 11
- 238000004422 calculation algorithm Methods 0.000 claims description 9
- FFBHFFJDDLITSX-UHFFFAOYSA-N benzyl N-[2-hydroxy-4-(3-oxomorpholin-4-yl)phenyl]carbamate Chemical compound OC1=C(NC(=O)OCC2=CC=CC=C2)C=CC(=C1)N1CCOCC1=O FFBHFFJDDLITSX-UHFFFAOYSA-N 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 238000012163 sequencing technique Methods 0.000 claims description 3
- 238000000638 solvent extraction Methods 0.000 claims description 3
- 241000284466 Antarctothoa delta Species 0.000 claims description 2
- 238000011161 development Methods 0.000 description 4
- 230000018109 developmental process Effects 0.000 description 4
- 238000000354 decomposition reaction Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000005192 partition Methods 0.000 description 2
- 230000011664 signaling Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000010937 topological data analysis Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及一种用于多电压等级电网模型的抗差状态估计方法,所述方法包括读取电网模型和SCADA量测;进行网络拓扑分析,将物理电网模型转化成母线‑支路计算模型;建立有功加权雅可比矩阵Ra ‑1/2Ba和无功计算加权雅可比矩阵Rr ‑1/2Br,计算有功系数矩阵A和无功系数矩阵B;利用正交变换获得有功系数矩阵和无功系数矩阵的上三角矩阵L和L1,进行状态估计迭代求解;计算量测残差vi,选择量测类型基准值Si,计算遥测估计值误差计算量测权因子,判断是否有量测处于降权区或淘汰区,没有则状态估计计算结束;有则根据量测权因子修正量测权重矩阵Ra ‑1和Rr ‑1,重新生成加权雅可比矩阵Ra ‑1/2Ba和Rr ‑1/2Br,计算A和B,返回步骤(4)。
Description
技术领域
本发明涉及一种抗差状态估计方法,具体讲涉及一种用于多电压等级电网模型的抗差状态估计方法。
背景技术
随着智能电网的建设和发展,电网规模越来越大,电网建模更加精细化,对电力系统在线分析计算,尤其是状态估计计算带来了新的挑战,如何从计算精度和计算速度两个方面对状态估计算法进行提升,快速准确的计算出大电网全网实时运行状态,为其他各类分析应用提供准确可靠的数据源,已经成为大电网状态估计计算所面临的新问题。因此,具有更强的不良数据辨识能力、更高计算精度和计算速度的抗差状态估计模型与算法已经成为电力系统研究中的一个热点方向。
现有的M估计通过等价权原理化为加权最小二乘状态估计的形式,引入权函数对量测权重进行灵活调整,继承了最小二乘的求解形式,在迭代计算过程中,根据残差动态调节量测权重,具有一定的抗差效果。但是,随着我国大规模可再生能源的集约化开发,特高压交、直流输电技术的快速发展,以及智能电网一体化运行的深入开展,各级调度中心对在线分析和仿真计算的精度提出了更高的要求,近年来已有部分省级调度中心试验将省、地电网模型拼接,并基于拼接后的大电网模型进行在线分析和调度员培训模拟。省、地电网模型拼接后电网模型电压等级跨度非常大,不同电压等级量测量的数值及其满刻度值往往也相差数倍甚至数十倍。抗差最小二乘估计中常用的权函数有Huber权函数、Hampel权函数、Turkey权函数、丹麦权函数、IGG权函数等。现有的基于IGG权函数法的抗差状态估计采用量测残差绝对值作为判断是否进行降权处理的依据,对不同电压等级的具有不同特点的量测没有进行区分,其对低电压等级设备量测或小容量设备量测的坏数据辨识能力有限。
另外,在抗差状态估计数学模型进行求解迭代过程中,需要根据权重调整情况重新建立加权雅可比矩阵及正交分解,因此,基于等价权原理的抗差最小二乘状态估计的计算效率较传统加权最小二乘方法有所下降,限制了该算法在在线分析领域中的实用化进程。
综上所述,要将抗差最小二乘状态估计算法更好的应用到实际电力系统中去,还需要对抗差状态估计算法进行深入研究。
发明内容
针对现有技术的不足,本发明提供一种用于多电压等级电网模型的抗差状态估计方法,在IGG分段权函数中引入量测类型基准值,实现对不同电压等级量测的区分,对处于降权区的量测权因子计算公式进行修正,整体提高抗差状态估计对不同电压等级量测坏数据的抗差性能。
本发明的目的是采用下述技术方案实现的:
一种用于多电压等级电网模型的抗差状态估计方法,其改进之处在于,所述方法包括
(1)读取电网模型和SCADA量测;
(2)进行网络拓扑分析,将物理电网模型转化成母线-支路计算模型;
(3)建立有功加权雅可比矩阵Ra -1/2Ba和无功计算加权雅可比矩阵Rr -1/2Br,计算有功系数矩阵A和无功系数矩阵B;
(4)利用正交变换获得有功系数矩阵和无功系数矩阵的上三角矩阵L和L1,进行状态估计迭代求解;
(5)计算量测残差vi,选择量测类型基准值Si,计算遥测估计值误差
(6)计算量测权因子,判断是否有量测处于降权区或淘汰区,没有则状态估计计算结束;有则根据量测权因子修正量测权重矩阵Ra -1和Rr -1,重新生成加权雅可比矩阵Ra -1/2Ba和Rr -1/2Br,计算A和B,返回步骤(4)。
优选的,所述步骤(1)包括从调度控制系统模型库获取电网模型,从SCADA应用获取实时遥信遥测数据。
优选的,所述步骤(2)包括根据电网模型电气联接关系和开关/刀闸运行状态,进行厂站母线拓扑分析,形成闭合开关/刀闸相连的结点集合,分配计算母线编号;根据投运支路联接的计算母线集合,进行电气岛划分,判断电气岛带电状态。
优选的,所述步骤(3)包括
(3.1)根据节点相关支路数进行节点排序优化,采用支路电抗形成导纳矩阵,建立有功量测的雅可比矩阵Ba;形成导纳矩阵,建立无功量测的雅可比矩阵Br;
(3.2)获取系统量测初始权重设置,建立有功加权雅可比矩阵Ra -1/2Ba和无功计算加权雅可比矩阵Rr -1/2Br;其中,Ra -1和Rr -1分别为有功量测权重矩阵和无功量测权重矩阵;
(3.3)根据公式 和 分别计算有功系数矩阵A和计算无功系数矩阵B,其中,V0取系统参考节点电压。
优选的,所述步骤(4)包括
(4.1)利用正交变换分别获得A和B的因子表L和L1,置迭代计数器l=1;
(4.2)按当前状态量电压幅值V(l)和电压相角θ(l)计算有功量测ha,根据公式计算a(l);其中,a(l)为有功自由矢量,za为有功量测矢量,ha为用状态量V(l)和θ(l)表达的有功量测;
(4.3)利用A的因子表L由公式AΔθ(l)=a(l)解出电压相角修正矢量Δθ(l),判断|Δθ(l)|max是否达到收敛精度;若不收敛,则计算θ(l+1)=θ(l)+Δθ(l),若收敛,则有功收敛标志位1;
(4.4)按当前状态量电压幅值V(l)和电压相角θ(l)计算无功量测hr,根据公式计算b(l);其中,b(l)为无功自由矢量,zr为无功量测矢量,hr为用状态量V(l)和θ(l)表达的无功量测;
(4.5)利用B的因子表L1由公式BΔV(l)=b(l)解出电压幅值修正矢量ΔV(l),判断|ΔV(l)|max是否达到收敛精度;若不收敛,则计算V(l+1)=V(l)+ΔV(l),若收敛,则无功收敛标志位1;
(4.6)判断有功和无功收敛标志位是否都为1,为1则状态估计迭代收敛,否则迭代计数器l=l+1;判断是否达到最大迭代次数,达到最大迭代次数,则状态估计不收敛,计算结束;未达到最大迭代次数则返回步骤(4.2)。
进一步地,采用Givens正交三角化的列超前并行消去算法求解L和L1,采用4进程并行模式,所述步骤(4.1)包括
(4.1.1)记加权雅可比矩阵Ra -1/2Ba为BW,Givens变换将BW的某一元素bil零化,可以看成矩阵BW的两行左乘一个二阶平面旋转阵,即 其中βi,βj为BW矩阵的i行和j行,i,j∈{1,…,m},i≠j,c=sinθij,s=cosθij,θij表示旋转的角度;由βi'中的bil’=0,可确定c和s;βj'的元素和βi'的其余元素为:
b'jt=cbjt+sbit,(t=1,…n)
b'it=-sbjt+cbit,(t=1,…,n,t≠l);
(4.1.2)BW矩阵分块存储,BW矩阵储存在进程0中,将BW矩阵按行分解为4块,序号为0、1、2、3,将第1、2、3块分别发送给进程1、2和3,进程0保存并处理第0块BW矩阵数据;
(4.1.3)进程0从BW矩阵取一行数据,经消去操作后,将此行数据发送给进程1,进程1完成相应的列消去操作后,再将此行数组发送给进程2,以此类推;进程0将数据发送后可再次从BW矩阵取下一行数据,进行消去操作,并重复上述过程,以实现并行计算;
(4.1.4)进程3完成列消去操作后,将此行数组保存到L矩阵当中;
(4.1.5)全部消去完成后,进程3将最终生成的L矩阵发送给进程0。
优选的,所述步骤(5)包括
(5.1)根据公式νi=zi-hi(V,θ)计算各量测残差νi;
(5.2)根据量测电压等级和量测设备类型选择量测类型基准值Si,根据Si计算量测估计误差
优选的,所述步骤(6)包括
根据改进IGG权函数 计算所有量测权因子,其中,为保权区,为降权区,为淘汰区;判断是否有量测处于降权区或淘汰区,没有则状态估计计算结束;有则根据量测权因子修正量测权重矩阵Ra -1和Rr -1,重新生成有功加权雅可比矩阵Ra -1/2Ba和无功计算加权雅可比矩阵Rr -1/2Br,计算有功系数矩阵A和计算无功系数矩阵B,返回步骤(4),直至没有量测处于降权区或淘汰区,则状态估计结束。
进一步地,所述公式 的变换过程如下:
(9.1)IGG权函数 采用量测残差绝对值对权函数进行分区,引入量测类型基准值Si,计算遥测数据估计值误差vi/Si,对权函数进行分区,区分了不同电压等级的量测;其中,k1,k2∈R,0<k1<k2,σ为电力系统遥测量的标准误差,|ν|为量测残差绝对值;
(9.2)IGG权函数 通过等价权化作用于加权最小二乘法目标函数时,处于降权区的量测的目标函数分量为Ji(x)=rik1σ|v|,修正原有权函数中处于降权区的量测权因子为
(3)通过(9.1)和(9.2)得到改进的IGG权函数为 其中,为遥测数据估计值误差,遥测数据估计值误差门槛k1,k2取值根据实际系统量测情况设置,k1大于0.03,k2大于k1。
与现有技术比,本发明的有益效果为:
(1)使用遥测估计值相对误差对不同电压等级量测进行区分。将量测满刻度值(因各量测满刻度值难以统计,实际系统使用量测类型基准值代替量测满刻度值)引入到抗差状态估计的权函数中,不同电压等级的量测具有相同的IGG权函数形式。但由于各量测满刻度值不同,其分区策略也不同,实现了对不同电压等级量测的区分,在现有抗差状态估计抗差性能的基础上,提高了对低电压等级量测坏数据的辨识能力和抗差能力,整体提高抗差状态估计的计算精度。
(2)充分利用各级调度中心自动化系统的多核CPU配置高性能服务器的优势,使用Givens正交三角化的列超前并行消去算法和矩阵向量并行乘法对状态估计修正方程的求解进行加速,从而显著缩短每次调节权重后加权雅可比矩阵的正交变换时间,提高抗差状态估计的计算速度。
附图说明
图1为本发明提供的一种用于多电压等级电网模型的抗差状态估计方法流程图;
图2为本发明提供的并行Givens变化程序设计。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步的详细说明。
本发明一种用于多电压等级电网模型的抗差状态估计方法,具体为基于等价权原理的抗差最小二乘状态估计保留了加权最小二乘状态估计的求解形式。改进IGG权函数法抗差最小二乘状态估计先利用加权最小二乘状态估计进行求解,根据WLS的结果计算量测估计值误差,确定量测全因子,对量测权重进行修正,重新建立加权雅可比矩阵、进行正交变换和快速分解迭代求解,直至无量测权重调整为止,状态估计计算结束。
1、抗差最小二乘状态估计数学模型
电力系统状态估计的量测方程是:
z=h(x)+v (1)
式中:z为量测向量,h(x)为量测方程,v为量测误差向量,x为状态向量。
由于各量测量的精度不同,对不同量测取不同权重Ri,精度高的取权重大些,精度低的取权重小些,当状态量的估计值为最优时,目标函数J(x)最小。
加权最小二乘法(WLS)目标函数为:
式中:ri为量测i的权重(即方差的倒数),vi残差v的第i个分量,m为量测个数,σ为标准差。
式(2)的迭代修正公式如下:
(HTR-1H)Δx(k)=HTR-1v(k) (3)
抗差估计的形式如下:
式中:vi为残差v的第i个分量;Pi为先验权;Ψ函数是待选择的,ui=vi/(σiv(li));v(li)为已知函数,li T为R-1/2H矩阵的行矢量;H为m×n阶雅可比矩阵,m,n分别表示量测量及状态量的个数,R-1/2H=[l1,l2,...,lm]T则表示加权雅可比矩阵。将式(4)整理成矩阵形式:
HTR-1Qv=0 (5)
式(6)的最小二乘解如下:
(HTR-1Q(k)H)Δx(k)=HTR-1Q(k)v(k) (7)
适当地选择Ψ函数可使对x的M估计具有较高的效率及较强的抗差能力。Ψ函数确定后,等价权Q(k)也就确定了,式(6)便是修正量求解方程。本专利抗差最小二乘状态估计的等价权函数采用IGG权函数的形式。
2、IGG权函数
本方法在利用等加权原理调节量测量权重时选用IGG权函数形式作为抗差最小二乘状态估计的权函数形式。IGG法的权函数分为保权区、降权区和淘汰区。在状态估计过程中,充分利用主体可靠信息,使其保持原权值不变。对可疑量测值依其可信度分别降权处理,对显著异常量测值则取权值为0,相当于删除该量测。
电力系统遥测量的标准误差σ大约为量测正常范围的0.5%~2.0%,误差大于±3σ的量测就可以认为是不良数据,在正常量测条件下,误差大于±3σ的量测值出现的概率仅为0.27%,实际系统中,可根据量测冗余度及全网量测整体精度调整不良数据检测标准。
设定两个界限k1σ和k2σ(k1,k2∈R,0<k1<k2),当|v|≤k1σ时,量测处于保权区,权因子为1,其目标函数分量Ji(x)=rivi 2保留了最小二乘的形式;当k1σ<|v|≤k2σ时,量测处于降权区,需要降低量测的观测权,其目标函数分量为Ji(x)=rik1σ|v|,小于最小二乘形式Ji(x)=rivi 2,达到量测降权目的;对于|v|>k2σ的情况,这些量测不应作为观测信息,其权因子取0,相当于删除该量测重新进行状态估计。权因子取值形式如下:
当式(8)的权函数作用于加权最小二乘法(WLS)目标函数时,处于降权区的量测的目标函数分量为Ji(x)=rik1σ|v|,为残差绝对值的一次函数,一定程度上起到降权目的,但其目标函数分量仍然随残差绝对值的增大而增大,说明在此区间内,量测坏数据的污染依然随残差绝对值增大而增大。另外,式(8)所示的IGG权函数未对不同电压等级量测进行区分。对此,本方法使用改进的IGG权函数。
3、改进IGG权函数
(1)由于不同电压等级、不同电气设备的量测量的满刻度值不同,甚至相差几倍至数十倍,其测量精度及误差也就不同。通常,较低电压等级的量测值较小,其坏数据的绝对残差也较小。如公式(8)中所示的利用量测残差绝对值作为分区界限,不能很好的处理不同量测特点的坏数据,尤其对低电压等级设备量测坏数据辨识能力有限。为此,本专利改进IGG权函数法抗差最小二乘状态估计,在权函数中引入量测类型基准值Si,对不同电压等级量测进行区分;
量测类型基准值通常是指该量测量的满刻度值。为使计算方便,线路功率量测、零阻抗支路量测、串联补偿器功率量测及母线电压量测可根据其电压等级选取相应的基准值;发电机功率量测及变压器功率量测则可使用设备额定容量作为量测基准值。
利用量测类型基准值Si计算遥测数据估计值误差,根据估计值误差对IGG权函数重新进行分段,增加了不同电压等级不良数据辨识的区分度。
(2)将权函数中处于降权区的量测权因子修改为此方法具有更强的抗差性能。
改进后权函数如下所示。
其中,为遥测数据估计值误差。遥测数据估计值误差门槛k1,k2取值根据实际系统量测情况设置,通常情况,k1大于0.03,k2大于k1。
式(9)与式(8)相比,有两点改进:第一、引入量测类型基准值Si对权函数重新进行分段,虽然不同电压等级的量测具有相同的权函数形式,但其权函数分段策略却不同,实现对不同电压等级量测的区分。对低电压等级量测,其具有较小的量测类型基准值,其处于降权区的量测的残差绝对值较小,符合调度自动化系统对不同电压等级量测精度的要求。第二、处于降权区的量测权因子形式与式(8)不同,利用式(9)的权函数,代入加权最小二乘法(WLS)目标函数,则处于降权区的量测的目标函数分量为常数:Ji(x)=rik1 2σ2,即,当量测处于降权区时,其残差污染不再随残差增大而扩大,有效抑制了残差污染。
如图1所示,本发明一种用于多电压等级电网模型的抗差状态估计方法流程如下:
步骤1:读取电网模型和SCADA实时数据。从调度控制系统模型库获取电网模型,从SCADA应用获取实时遥信遥测数据。
步骤2:进行网络拓扑分析,通过拓扑分析将物理电网模型转化成母线-支路计算模型。根据电网模型电气联接关系和开关/刀闸运行状态,进行厂站母线拓扑分析,形成闭合开关/刀闸相连的结点集合(即计算母线),分配计算母线编号;根据投运支路联接的计算母线集合,进行电气岛划分,判断电气岛带电状态。
步骤3:根据量测量、量测权重及导纳矩阵分别建立有功加权雅可比矩阵Ra -1/2Ba和无功计算加权雅可比矩阵Rr -1/2Br,计算有功系数矩阵A和无功系数矩阵B;
步骤3-1:根据节点相关支路数进行节点排序优化,用支路电抗形成导纳矩阵,建立有功量测的雅可比矩阵Ba;形成正常导纳矩阵,建立无功量测的雅可比矩阵Br;
步骤3-2:获取系统量测初始权重设置,建立有功加权雅可比矩阵Ra -1/2Ba和无功计算加权雅可比矩阵Rr -1/2Br,其中Ra -1和Rr -1分别为有功量测权重矩阵和无功量测权重矩阵;
步骤3-3:根据公式 和 分别计算有功系数矩阵A和计算无功系数矩阵B,其中V0取系统参考节点电压。
步骤4:利用正交变换获得有功和无功迭代上三角矩阵L和L1,进行状态估计迭代求解;
步骤4-1:利用Givens正交三角化的列超前并行消去算法分别获得A和B的因子表L和L1,置迭代计数器l=1。以4进程并行模式为例说明,所述步骤4-1包括:
步骤4-1-1:记加权雅可比矩阵Ra -1/2Ba为BW,Givens变换将BW的某一元素bil零化,可以看成矩阵BW的两行左乘一个二阶平面旋转阵,即 其中βi,βj为BW矩阵的i行和j行,i,j∈{1,…,m},i≠j,c=sinθij,s=cosθij,θij表示旋转的角度。由βi'中的bil’=0,可以确定c和s。βj'的元素和βi'的其余元素为:
b'jt=cbjt+sbit,(t=1,…n) (11)
b'it=-sbjt+cbit,(t=1,…,n,t≠l) (12)
步骤4-1-2:BW矩阵分块存储。BW矩阵储存在进程0中,将BW矩阵按行分解为4块,序号为0、1、2、3,将第1、2、3块分别发送给进程1、2和3,进程0保存并处理第0块BW矩阵数据。
步骤4-1-3:进程0从BW矩阵取一行数据,经消去操作后,将此行数据发送给进程1,进程1完成相应的列消去操作后,再将此行数组发送给进程2,以此类推。进程0将数据发送后可再次从BW矩阵取下一行数据,进行消去操作,并重复上述过程,以实现并行计算。
步骤4-1-4:进程3完成列消去操作后,将此行数组保存到L矩阵当中。
步骤4-1-5:全部消去完成后,进程3将最终生成的L矩阵发送给进程0。
步骤4-2:按当前状态量电压幅值V(l)和电压相角θ(l)计算有功量测ha,根据公式计算a(l)。其中:a(l)为有功自由矢量,za为有功量测矢量,ha为用状态量V(l)和θ(l)表达的有功量测;
步骤4-3:利用A的因子表L由公式AΔθ(l)=a(l)解出电压相角修正矢量Δθ(l),判断|Δθ(l)|max是否达到收敛精度;若不收敛,则计算θ(l+1)=θ(l)+Δθ(l),若收敛,则有功收敛标志位1;
步骤4-4:按当前状态量电压幅值V(l)和电压相角θ(l)计算无功量测hr,根据公式计算b(l)。其中:b(l)为无功自由矢量,zr为无功量测矢量,hr为用状态量V(l)和θ(l)表达的无功量测;
步骤4-5:利用B的因子表L1由公式BΔV(l)=b(l)解出电压幅值修正矢量ΔV(l),判断|ΔV(l)|max是否达到收敛精度;若不收敛,则计算V(l+1)=V(l)+ΔV(l),若收敛,则无功收敛标志位1;
步骤4-6:判断有功和无功收敛标志位是否都为1,若都是1则状态估计迭代收敛,否则迭代计数器l=l+1;判断是否达到最大迭代次数,若达到最大迭代次数,则状态估计不收敛,计算结束;未达到最大迭代次数则返回步骤4-2。
步骤5:计算遥测估计值误差,根据式(9)计算各量测权因子,若有量测权重调整,则修正量测权重矩阵,重新建立加权雅可比矩阵,重新计算A和B,返回步骤(4);若无量测权重调整,则状态估计结束。
步骤5-1:根据公式νi=zi-hi(V,θ)计算各量测残差νi;
步骤5-2:根据量测电压等级和量测设备类型选择量测类型基准值Si,具体Si取值可参考2013年8月7日国家电网公司国家电力调度控制中心印发《智能电网调度控制系统实用化验收办法(试行)》,根据Si计算量测估计误差|νi|/Si;
步骤6:根据改进IGG权函数 计算所有量测权因子,其中,为保权区,为降权区,为淘汰区;判断是否有量测处于降权区或淘汰区,没有则状态估计计算结束;有则根据量测权因子修正量测权重矩阵Ra -1和Rr -1,重新生成有功加权雅可比矩阵Ra -1/2Ba和无功计算加权雅可比矩阵Rr -1/2Br,计算有功系数矩阵A和计算无功系数矩阵B,返回步骤(4),直至没有量测处于降权区或淘汰区,则状态估计结束。
最后应当说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,所属领域的普通技术人员参照上述实施例依然可以对本发明的具体实施方式进行修改或者等同替换,这些未脱离本发明精神和范围的任何修改或者等同替换,均在申请待批的本发明的权利要求保护范围之内。
Claims (9)
1.一种用于多电压等级电网模型的抗差状态估计方法,其特征在于,所述方法包括
(1)读取电网模型和SCADA量测;
(2)进行网络拓扑分析,将物理电网模型转化成母线‐支路计算模型;
(3)建立有功加权雅可比矩阵Ra -1/2Ba和无功计算加权雅可比矩阵Rr -1/2Br,计算有功系数矩阵A和无功系数矩阵B;
(4)利用正交变换获得有功系数矩阵和无功系数矩阵的上三角矩阵L和L1,进行状态估计迭代求解;
(5)计算量测残差vi,选择量测类型基准值Si,计算遥测估计值误差
(6)计算量测权因子,判断是否有量测处于降权区或淘汰区,没有则状态估计计算结束;有则根据量测权因子修正量测权重矩阵Ra -1和Rr -1,重新生成加权雅可比矩阵Ra -1/2Ba和Rr -1/2Br,计算A和B,返回步骤(4)。
2.如权利要求1所述的一种用于多电压等级电网模型的抗差状态估计方法,其特征在于,所述步骤(1)包括从调度控制系统模型库获取电网模型,从SCADA应用获取实时遥信遥测数据。
3.如权利要求1所述的一种用于多电压等级电网模型的抗差状态估计方法,其特征在于,所述步骤(2)包括根据电网模型电气联接关系和开关/刀闸运行状态,进行厂站母线拓扑分析,形成闭合开关/刀闸相连的结点集合,分配计算母线编号;根据投运支路联接的计算母线集合,进行电气岛划分,判断电气岛带电状态。
4.如权利要求1所述的一种用于多电压等级电网模型的抗差状态估计方法,其特征在于,所述步骤(3)包括
(3.1)根据节点相关支路数进行节点排序优化,采用支路电抗形成导纳矩阵,建立有功量测的雅可比矩阵Ba;形成导纳矩阵,建立无功量测的雅可比矩阵Br;
(3.2)获取系统量测初始权重设置,建立有功加权雅可比矩阵Ra -1/2Ba和无功计算加权雅可比矩阵Rr -1/2Br;其中,Ra -1和Rr -1分别为有功量测权重矩阵和无功量测权重矩阵;
(3.3)根据公式 和 分别计算有功系数矩阵A和计算无功系数矩阵B,其中,V0取系统参考节点电压。
5.如权利要求1所述的一种用于多电压等级电网模型的抗差状态估计方法,其特征在于,所述步骤(4)包括
(4.1)利用正交变换分别获得A和B的因子表L和L1,置迭代计数器l=1;
(4.2)按当前状态量电压幅值V(l)和电压相角θ(l)计算有功量测ha,根据公式计算a(l);其中,a(l)为有功自由矢量,za为有功量测矢量,ha为用状态量V(l)和θ(l)表达的有功量测;
(4.3)利用A的因子表L由公式AΔθ(l)=a(l)解出电压相角修正矢量Δθ(l),判断|Δθ(l)|max是否达到收敛精度;若不收敛,则计算θ(l+1)=θ(l)+Δθ(l),若收敛,则有功收敛标志位1;
(4.4)按当前状态量电压幅值V(l)和电压相角θ(l)计算无功量测hr,根据公式计算b(l);其中,b(l)为无功自由矢量,zr为无功量测矢量,hr为用状态量V(l)和θ(l)表达的无功量测;
(4.5)利用B的因子表L1由公式BΔV(l)=b(l)解出电压幅值修正矢量ΔV(l),判断|ΔV(l)|max是否达到收敛精度;若不收敛,则计算V(l+1)=V(l)+ΔV(l),若收敛,则无功收敛标志位1;
(4.6)判断有功和无功收敛标志位是否都为1,为1则状态估计迭代收敛,否则迭代计数器l=l+1;判断是否达到最大迭代次数,达到最大迭代次数,则状态估计不收敛,计算结束;未达到最大迭代次数则返回步骤(4.2)。
6.如权利要求5所述的一种用于多电压等级电网模型的抗差状态估计方法,其特征在于,采用Givens正交三角化的列超前并行消去算法求解L和L1,采用4进程并行模式,所述步骤(4.1)包括
(4.1.1)记加权雅可比矩阵Ra -1/2Ba为BW,Givens变换将BW的某一元素bil零化,可以看成矩阵BW的两行左乘一个二阶平面旋转阵,即 其中βi,βj为BW矩阵的i行和j行,i,j∈{1,…,m},i≠j,c=sinθij,s=cosθij,θij表示旋转的角度;由βi'中的bil’=0,可确定c和s;βj'的元素和βi'的其余元素为:
b'jt=cbjt+sbit,(t=1,…n)
b′′it=-sbjt+cbit,(t=1,…,n,t≠l);
(4.1.2)BW矩阵分块存储,BW矩阵储存在进程0中,将BW矩阵按行分解为4块,序号为0、1、2、3,将第1、2、3块分别发送给进程1、2和3,进程0保存并处理第0块BW矩阵数据;
(4.1.3)进程0从BW矩阵取一行数据,经消去操作后,将此行数据发送给进程1,进程1完成相应的列消去操作后,再将此行数组发送给进程2,以此类推;进程0将数据发送后可再次从BW矩阵取下一行数据,进行消去操作,并重复上述过程,以实现并行计算;
(4.1.4)进程3完成列消去操作后,将此行数组保存到L矩阵当中;
(4.1.5)全部消去完成后,进程3将最终生成的L矩阵发送给进程0。
7.如权利要求1所述的一种用于多电压等级电网模型的抗差状态估计方法,其特征在于,所述步骤(5)包括
(5.1)根据公式νi=zi-hi(V,θ)计算各量测残差νi;
(5.2)根据量测电压等级和量测设备类型选择量测类型基准值Si,根据Si计算量测估计误差
8.如权利要求1所述的一种用于多电压等级电网模型的抗差状态估计方法,其特征在于,所述步骤(6)包括
根据改进IGG权函数 计算所有量测权因子,其中, 为保权区, 为降权区, 为淘汰区,遥测数据估计值误差门槛k1,k2取值根据实际系统量测情况设置,k2大于k1;判断是否有量测处于降权区或淘汰区,没有则状态估计计算结束;有则根据量测权因子修正量测权重矩阵Ra -1和Rr -1,重新生成有功加权雅可比矩阵Ra -1/2Ba和无功计算加权雅可比矩阵Rr -1/2Br,计算有功系数矩阵A和计算无功系数矩阵B,返回步骤(4),直至没有量测处于降权区或淘汰区,则状态估计结束。
9.如权利要求8所述的一种用于多电压等级电网模型的抗差状态估计方法,其特征在于,所述公式 的变换过程如下:
(6.1)IGG权函数 采用量测残差绝对值对权函数进行分区,引入量测类型基准值Si,计算遥测数据估计值误差vi/Si,对权函数进行分区,区分了不同电压等级的量测;其中,k1,k2∈R,0<k1<k2,σ为电力系统遥测量的标准误差,|ν|为量测残差绝对值;
(6.2)IGG权函数 通过等价权化作用于加权最小二乘法目标函数时,处于降权区的量测的目标函数分量为Ji(x)=rik1σ|v|,修正原有权函数中处于降权区的量测权因子为
(6.3)通过(6.1)和(6.2)得到改进的IGG权函数为 其中,为遥测数据估计值误差,遥测数据估计值误差门槛k1,k2取值根据实际系统量测情况设置,k1大于0.03,k2大于k1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410450017.7A CN104184144B (zh) | 2014-09-05 | 2014-09-05 | 一种用于多电压等级电网模型的抗差状态估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410450017.7A CN104184144B (zh) | 2014-09-05 | 2014-09-05 | 一种用于多电压等级电网模型的抗差状态估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104184144A CN104184144A (zh) | 2014-12-03 |
CN104184144B true CN104184144B (zh) | 2016-08-24 |
Family
ID=51964984
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410450017.7A Active CN104184144B (zh) | 2014-09-05 | 2014-09-05 | 一种用于多电压等级电网模型的抗差状态估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104184144B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104992010B (zh) * | 2015-06-25 | 2018-02-13 | 国电南瑞科技股份有限公司 | 一种基于拓扑分区的多断面联合参数估计方法 |
CN105512502B (zh) * | 2016-01-13 | 2018-04-17 | 重庆大学 | 一种基于残差归一化的权函数最小二乘状态估计方法 |
CN107611962B (zh) * | 2017-08-24 | 2020-06-16 | 中国科学院深圳先进技术研究院 | 电网系统支路搜索方法、系统及电子设备 |
CN107657019B (zh) * | 2017-09-26 | 2020-05-12 | 广州供电局有限公司 | 电网系统的网络拓扑结构获取方法和系统 |
EP3579366A1 (en) * | 2018-06-08 | 2019-12-11 | Siemens Aktiengesellschaft | Method and apparatus for control of an electric grid |
CN110943445A (zh) * | 2018-09-21 | 2020-03-31 | 中国电力科学研究院有限公司 | 一种柔直及交流混合电网抗差状态估计方法和装置 |
CN110417009A (zh) * | 2019-07-29 | 2019-11-05 | 天津大学 | 基于不同采样周期数据的配电网混合抗差状态估计方法 |
CN111062610B (zh) * | 2019-12-16 | 2022-07-29 | 国电南瑞科技股份有限公司 | 一种基于信息矩阵稀疏求解的电力系统状态估计方法和系统 |
CN112383048B (zh) * | 2020-10-20 | 2024-07-09 | 重庆大学 | 考虑公专用配电变压器测量特点的中压馈线电网三相状态估计方法 |
CN112749369A (zh) * | 2021-01-19 | 2021-05-04 | 东方电子股份有限公司 | 一种基于Givens正交变换的电力系统状态估计方法 |
CN113139288A (zh) * | 2021-04-21 | 2021-07-20 | 广东电网有限责任公司广州供电局 | 基于gpu二次加速的配电网分布式抗差状态估计方法及装置 |
CN113991722B (zh) * | 2021-10-08 | 2024-05-31 | 国电南瑞科技股份有限公司 | 多馈入直流短路比计算方法、装置、存储介质及计算设备 |
CN116365555B (zh) * | 2023-05-10 | 2023-08-11 | 国网江苏省电力有限公司苏州供电分公司 | 一种计及局部三相不对称的电网状态估计方法和系统 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6253120B1 (en) * | 1997-04-28 | 2001-06-26 | Seiko Seiki Kabushiki Kaisha | Position and/or force controlling apparatus using sliding mode decoupling control |
CN101383511B (zh) * | 2008-10-10 | 2010-08-04 | 清华大学 | 基于数据采集系统量测数据的电力系统状态估计方法 |
CN101599643B (zh) * | 2009-04-23 | 2011-04-27 | 清华大学 | 一种基于指数型目标函数的电力系统抗差状态估计方法 |
-
2014
- 2014-09-05 CN CN201410450017.7A patent/CN104184144B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN104184144A (zh) | 2014-12-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104184144B (zh) | 一种用于多电压等级电网模型的抗差状态估计方法 | |
CN107453357B (zh) | 一种基于分层求解的配电网状态估计方法 | |
CN101599643B (zh) | 一种基于指数型目标函数的电力系统抗差状态估计方法 | |
CN110516912B (zh) | 一种配电台区户变关系的识别方法 | |
CN106505557B (zh) | 一种遥测错误辨识方法及装置 | |
WO2022021726A1 (zh) | 一种基于pmu的电力系统状态估计性能评价方法 | |
CN110543720B (zh) | 基于sdae-elm伪量测模型的状态估计方法 | |
CN104836223A (zh) | 电网参数错误与不良数据协同辨识与估计方法 | |
CN103236691B (zh) | 基于复仿射数学理论的三相不平衡潮流计算方法 | |
CN110070282A (zh) | 一种基于综合关联度的低压台区线损影响因素分析方法 | |
WO2016078477A1 (zh) | 一种变电站三相线性广义状态估计方法 | |
CN107038292A (zh) | 一种基于自适应多变量非参数核密度估计的多风电场出力相关性建模方法 | |
CN109088407B (zh) | 基于深度信念网络伪量测建模的配电网状态估计方法 | |
CN109711435A (zh) | 一种基于遗传算法的支持向量机在线电压稳定性监测方法 | |
CN111460374A (zh) | 一种考虑节点差异性的配电网d-pmu优化配置方法 | |
CN105406471A (zh) | 电网不良数据辨识与估计方法 | |
CN103034787A (zh) | 一种微电网状态估计的方法 | |
CN107634516A (zh) | 一种基于灰色‑马尔可夫链的配网状态估计方法 | |
CN111064180A (zh) | 基于ami潮流匹配的中压配电网拓扑检测与辨识方法 | |
CN108132441A (zh) | 储能电池模组荷电状态的运行范围确定方法及装置 | |
CN105048473B (zh) | 变电站自动电压控制中低压无功设备电压灵敏度统计方法 | |
CN109829246B (zh) | 一种基于参数综合可疑度的线路参数辨识方法 | |
CN111091223B (zh) | 一种基于物联网智能感知技术的配变短期负荷预测方法 | |
CN106372440B (zh) | 一种并行计算的配电网自适应抗差状态估计方法及装置 | |
CN104239716B (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 |