CN102636706B - 一种电网中参数误差支路的辨识方法 - Google Patents
一种电网中参数误差支路的辨识方法 Download PDFInfo
- Publication number
- CN102636706B CN102636706B CN201210061638.7A CN201210061638A CN102636706B CN 102636706 B CN102636706 B CN 102636706B CN 201210061638 A CN201210061638 A CN 201210061638A CN 102636706 B CN102636706 B CN 102636706B
- Authority
- CN
- China
- Prior art keywords
- branch road
- residual
- error
- parameter
- branches
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title abstract description 13
- 230000002093 peripheral effect Effects 0.000 claims abstract 5
- 238000005259 measurement Methods 0.000 claims description 25
- 238000012850 discrimination method Methods 0.000 claims description 12
- 238000001514 detection method Methods 0.000 description 16
- 239000011159 matrix material Substances 0.000 description 12
- 238000011109 contamination Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 230000011664 signaling Effects 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 239000004615 ingredient Substances 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Abstract
Description
技术领域
本发明涉及一种电网中参数误差支路的辨识方法,属于电力系统运行和控制技术领域。
背景技术
现代化的调度系统要求能迅速、准确而全面地掌握电力系统的实际运行状态,因此能量管理系统(energy management system, EMS)得到了广泛的应用。在EMS中,状态估计是重要的组成部分,其结果直接影响电网调度的智能化分析与决策,状态估计领域所涉及的状态估计、不良数据的检测和辨识以及结线分析,都是建立在电力系统的开关遥信量均是正确、所给定的网络元件参数也均是正确的基础上。可见,一个正确的网络结构和元件参数值对调度中心是很重要的。它不仅是进行状态估计的依据,也是实现电网的自动化(如自动发电控制、经济调度、自动电压/无功控制和稳定安全分析等)所必不可少的。因此,如何检测辨识出网络中存在误差甚至是错误的参数,从而可以为调度人员和数据库提供一个完整、可靠的具有一定精度的参数值奠定基础,则显得非常有意义。解决这个问题的有效的方法是设计相应的软件,检测辨识存在误差甚至错误的网络参数以便进行修正更新。
数据库中所得到的电力系统网络的参数带有误差的原因一般有:
一、有时因缺少实测参数量而直接采用设计参数或参数测量条件与实际运行条件差别较大,这都会使给定的参数值与实际运行中的元件参数有差别;
二、不准确的制造数据(如在变压器模型中铁耗习惯上地被忽略)以及缺乏线路长度估计。对于超过200km的线路,如果采用集中参数模型则误差预计超过1%;
三、网架结构的改变并未适当地在数据库中得到更新(如一段架空线路已被电缆取代);
四、对于温度(尤其是串联电阻)或者环境条件(尤其是并联电抗)的依赖,实际运行中的元件参数因环境变化等原因而局部地、缓慢地发生着变化;
五、在未通知控制中心的情况下进行局部的分接头变更——遥信错误,调度中心对运行中的自动调压变压器的分接头位置和补偿电容器的组数掌握得不一定确切。
由于参数的误差和不确定性会给状态估计、不良数据辨识及其他高级安全分析程序带来不利的影响,错误的分接头信息和不正确的参数值可能产生以下后果:
(1)状态估计可信度的显著的降低,由此,可能误导其他高级应用程序(如安全评估);
(2)由于电网参数存在误差甚至出现错误参数,正确的量测量则可能会被辨识为坏数据;
(3)调度员对于状态估计结果失去信任。
基于以上原因,电网中参数误差支路的准确辨识就显得尤为重要。传统的参数误差支路辨识方法直接利用残差进行比较,最大值者为存在参数误差的支路,但由于当存在单个参数误差的情况下,可疑支路检测出来不止一个,即除了本身的参数误差支路的残差呈现出超过检测门槛值的特征外,还有一些准确参数支路的残差也超过门槛值,这种“残差污染”现象的结果使参数误差支路被模糊,导致辨识的困难。另外,在多个参数误差支路的情况下,情况更为严重,即多个参数误差相互作用下,导致在部分的或全部参数误差支路的残差接近于正常残差的现象,引起所谓“残差淹没”现象。就是说,残差淹没使参数误差支路不能呈现任何残差特征。残差淹没一般发生在强相关的同类型测点间,例如同为有功测点或无功测点之间。可以认为,不同类型的测点间,或者虽为同类测点,但在系统内的电气距离较大或弱相关时,事实上不会发生残差淹没。此外,实际上全部参数误差支路测点的残差都淹没的情况是比较少的,多数情况是部分参数误差支路点的残差被淹没。但是,不管属于哪种情况,都使参数误差支路的辨识不成功。这就有必要研究新的参数误差支路辨识方法。
发明内容
本发明所要解决的技术问题在于克服现有技术的不足,提供一种电网中参数误差支路的辨识方法,能够在离线情况下有效而准确地逐一辨识出网络中存在参数误差的支路。
本发明的电网中参数误差支路的辨识方法,包括以下步骤:
步骤A、从电网SCADA(supervisory control and data acquisition,数据采集与监控)系统采集到的量测数据中,选取某一断面的量测数据,此断面不包含不良数据和拓扑错误;
步骤B、利用所选取的量测数据,对各节点的电压幅值和相角进行离线状态估计;利用状态估计结果计算各状态量的量测残差,并进而求得电网中各支路的标准化残差;
,
步骤D、逐一对各可疑支路按照以下方法进行加权残差率绝对值求和,得到各可疑支路的加权残差率绝对值和:
步骤D1、对本支路的残差率绝对值赋予较大的权值,并同时寻找出其余可疑支路中与本支路首末节点直接相连的所有周边支路;
步骤E、选取加权残差率绝对值和最大的可疑支路作为参数误差支路。
进一步地,步骤E之后还包括:
步骤F、对步骤E中选取的参数误差支路进行参数估计,并修正参数误差;然后重新执行步骤B—步骤E,直至再无可疑支路。
优选地,所述较大的权值取值为1,较小的权值的取值范围为0.2—0.3,更小的权值的取值范围为0.1以下。
相比现有技术,本发明具有以下有益效果:
一、本发明方法能够有效避免“残差污染”和“残差淹没”现象所导致的误检,漏检,提高辨识准确率;
二、本发明基于离线取样计算,进行基于加权残差率绝对值的参数误差支路辨识,同时修正参数误差后,并在同一断面下重复进行。其易于与现有状态估计程序相结合,对程序的改变极小。
附图说明
图1为本发明方法的流程示意图;
图2为支路量测距离概念示意图;
图3为不同距离支路上参数误差对量测估计误差的影响示意图;
图4为标准正态分布图。
具体实施方式
下面结合附图对本发明的技术方案进行详细说明:
参数误差支路的辨识是现代电力系统状态估计的重要功能之一,其目的在于检测辨识存在误差甚至错误的参数,为提高参数估计的准确性以及状态估计可靠性奠定基础。常规状态估计程序采用加权最小二乘算法,迭代求得最优估计状态量,选取收敛标准可取基准电压幅值的10-4~10-6。
而状态估计误差方差阵是:
量测估计误差是:
量测估计误差方差阵是:
在状态真值附近线性化:
相应地,量测估计误差方差阵为:
则:
以及
或
则标准化形式的残差方程如下:
由标准化残差方程可知:
这意味着,在正常量测条件下的标准化残差,也是0均值的正态分布的随机变量。因此,标准化残差协方差矩阵的对角元素就是相应标准化残差的方差,有:
可知上式右端矩阵的对角元素均为1,故有:
的标准正态分布随机变量。
因此,在正常量测条件下,下述事件的概率为
于是,标准化残差门槛值可定为:
法在检测辨识可疑参数方面存在着某些问题,在单个参数误差的情况下,可疑数据检测出来不止一个,即除了参数误差支路据点的残差呈现出超过检测门槛值的特征外,还有一些准确参数的正常测点的残差也超过门槛值,这种现象称为残差污染。残差污染的结果使参数误差支路被模糊,导致辨识参数误差的困难。
在多个参数误差的情况下,情况更为严重,会引起所谓“残差淹没”情况。所谓残差淹没,是指多个参数误差相互作用的结果,导致在部分的或全部参数误差支路上的残差接近于正常残差的现象。就是说,残差淹没使参数误差支路不能呈现任何残差特性。因此,用法不能把它们作为可疑支路检测出来,造成漏检的结果。残差淹没一般发生在强相关的同类型测点间,此外,实际上全部参数误差支路的残差都淹没的情况是比较少的,多数情况是部分参数误差支路的残差淹没。但是,不管属于哪种情况,都使参数误差的检测不成功。这就有必要研究新的参数误差检测方法。
本发明的思路是首先将大于检测门槛值的支路作为可疑支路;然后根据参数误差对不同距离的支路量测残差的影响不同的原理,对各可疑支路进行加权残差率绝对值求和,最后选取加权残差率绝对值和最大的可疑支路作为参数误差支路。
参数误差相关分析如下:
一、支路量测距离的概念:如图2所示,若1支路为在电网中选取的某条待辨识支路,则定义为距离1支路,与距离1支路首末节点直接相连的所有周边支路定为距离2支路,与所有距离2支路首末节点直接相连的不包括距离2支路的所有外层支路定义为距离3支路。
二、参数误差对不同距离的支路量测残差的影响不同,即如果某支路存在参数误差,在不同距离的支路上,存在参数误差的量测估计量误差与无参数误差的量测估计量误差之比各不相同,图3显示了不同距离支路上参数误差对量测估计误差的影响,具体表现为:当参数误差从5%增值至25%时,通过对存在参数误差时的量测估计误差和正常的量测估计误差之比可以看出,距离1支路的误差比明显偏大,而距离2支路的误差比相对较小,距离3支路的误差比则更小,而距离空间上更远的支路4、5、6的误差比已经非常小。
三、在与所有距离3支路首末节点直接相连的不包括距离3支路的更外层支路上,残差的变化已经非常不明显,即可说明参数误差对支路功率量测残差的影响仅是局部影响。
因此,可以根据以上结论,逐一对可疑支路进行加权残差率绝对值求和,选取加权残差率绝对值和最大的可疑支路作为参数误差支路。
具体的,本发明方法如图1所示,按照以下步骤执行:
(1)初始化,包括:对节点次序优化、形成节点导纳矩阵、设置门槛值、状态量设置初值、分配内存。
(2)在SCADA系统定时采样的测量数据库中,选取某一特定的量测断面,因为是离线采样、计算,所以此断面要求应当不包含量测坏数据及拓扑错误。因此,在以下状态估计计算中,当残差超过阈值时必然表明存在参数误差。量测断面中包括:节点电压幅值、发电机有功功率和无功功率、负荷有功功率和无功功率、线路首末两端的有功功率和无功功率。
(3)进行离线状态估计,采用常规的最小二乘算法,状态量为各节点的电压幅值和相角(除平衡节点相角为零),迭代计算状态修正量,与门槛值进行比较(为基准电压幅值的10-4~10-5),状态修正量绝对值最大者小于给定的门槛值时结束计算,得到此量测断面下的状态估计量。
(5)可疑集检查:大于检测门槛值的支路作为可疑支路纳入可疑集S,并计算残差率绝对值,其中,的计算公式为,、分别为第i个量测量的残差和量测值,,为量测量总数。如果S为空集,则表明不存在参数误差支路;如果S不为空集,则表明电网某支路存在参数误差。
(6)将残差率绝对值按照由大到小的顺序排列在可疑集中,按照如下方法逐一对列入可疑集中的支路进行加权残差率绝对值求和:
步骤一:在可疑集中,将对应于本支路的距离1支路残差率绝对值赋予较大的权值(本具体实施方式中其值为1),并同时寻找出距离2支路残差率绝对值;
步骤三:将距离3支路残差率绝对值赋予更小的权值(本具体实施方式中其取值范围为0.1以下),最后将距离1支路、距离2支路、距离3支路的所有残差率绝对值进行加权求和,计算得到可疑集中某支路的加权残差率绝对值和。
(7)按照如上方法计算可疑集中的所有支路的加权残差率绝对值和,进行比较,取值最大者为存在参数误差的支路。
(8)若网络中存在多个参数误差支路,由于参数误差对残差影响是局部影响,因此,其相互影响的程度较小。当辨识出某条参数误差支路后,进行参数估计,并修正参数误差。返回(3)重新检测辨识网络中所有的参数误差支路,直至可疑集S为空集时结束。
Claims (5)
1.一种电网中参数误差支路的辨识方法,其特征在于,包括以下步骤:
步骤A、从电网SCADA系统采集到的量测数据中,选取某一断面的量测数据,此断面不包含不良数据和拓扑错误;
步骤B、利用所选取的量测数据,对各节点的电压幅值和相角进行离线状态估计;利用状态估计结果计算各状态量的量测残差,并进而求得电网中各支路的标准化残差;
步骤D、逐一对各可疑支路按照以下步骤D1、D2和D3进行加权残差率绝对值求和,得到各可疑支路的加权残差率绝对值和:
步骤E、选取加权残差率绝对值和最大的可疑支路作为参数误差支路。
2.如权利要求1所述电网中参数误差支路的辨识方法,其特征在于,步骤E之后还包括:
步骤F、对步骤E中选取的参数误差支路进行参数估计,并修正参数误差;然后重新执行步骤B—步骤E,直至再无可疑支路。
3.如权利要求1或2所述电网中参数误差支路的辨识方法,其特征在于,所述较大的权值取值为1,较小的权值的取值范围为0.2—0.3,更小的权值的取值范围为0.1以下。
4.如权利要求1或2所述电网中参数误差支路的辨识方法,其特征在于,所述预设标准化残差阈值为2.81。
5.如权利要求1或2所述电网中参数误差支路的辨识方法,其特征在于,步骤B中所述对各节点的电压幅值和相角进行离线状态估计,采用最小二乘算法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210061638.7A CN102636706B (zh) | 2012-03-12 | 2012-03-12 | 一种电网中参数误差支路的辨识方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210061638.7A CN102636706B (zh) | 2012-03-12 | 2012-03-12 | 一种电网中参数误差支路的辨识方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102636706A CN102636706A (zh) | 2012-08-15 |
CN102636706B true CN102636706B (zh) | 2014-02-19 |
Family
ID=46621175
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210061638.7A Expired - Fee Related CN102636706B (zh) | 2012-03-12 | 2012-03-12 | 一种电网中参数误差支路的辨识方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102636706B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103823140B (zh) * | 2014-02-27 | 2016-05-18 | 华北电力大学 | 基于道路-回路方程的电网拓扑错误辨识系统及其方法 |
CN103941130B (zh) * | 2014-05-13 | 2017-02-15 | 南方电网科学研究院有限责任公司 | 一种基于多量测断面的可疑支路识别方法 |
CN104809339B (zh) * | 2015-04-22 | 2018-01-05 | 华南理工大学 | 基于有功平衡的scada系统电网有功坏数据识别方法 |
CN108022046A (zh) * | 2017-12-05 | 2018-05-11 | 国网江西省电力有限公司景德镇供电分公司 | 一种电力系统数据质量评估方法、存储介质及设备 |
CN107958334A (zh) * | 2017-12-05 | 2018-04-24 | 国网江西省电力有限公司景德镇供电分公司 | 一种用于电力行业对电力数据进行分析检验的方法 |
CN111931318B (zh) * | 2020-07-07 | 2022-09-27 | 中国电力科学研究院有限公司 | 一种基于图计算的供电路径分析方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101593977A (zh) * | 2009-04-17 | 2009-12-02 | 哈尔滨工业大学 | 电力网络支路参数估计方法 |
CN101635457A (zh) * | 2009-05-14 | 2010-01-27 | 国家电网公司 | 一种基于状态估计残差参数灵敏度的电网参数估计方法 |
CN102136105A (zh) * | 2011-03-22 | 2011-07-27 | 哈尔滨工业大学 | 基于相量测量单元的电力网络新息图参数估计方法 |
CN102280877A (zh) * | 2011-07-25 | 2011-12-14 | 清华大学 | 一种多量测断面的电力系统不良支路参数辨识方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103154934A (zh) * | 2010-06-11 | 2013-06-12 | Abb研究有限公司 | 检测状态估计网络模型数据误差 |
-
2012
- 2012-03-12 CN CN201210061638.7A patent/CN102636706B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101593977A (zh) * | 2009-04-17 | 2009-12-02 | 哈尔滨工业大学 | 电力网络支路参数估计方法 |
CN101635457A (zh) * | 2009-05-14 | 2010-01-27 | 国家电网公司 | 一种基于状态估计残差参数灵敏度的电网参数估计方法 |
CN102136105A (zh) * | 2011-03-22 | 2011-07-27 | 哈尔滨工业大学 | 基于相量测量单元的电力网络新息图参数估计方法 |
CN102280877A (zh) * | 2011-07-25 | 2011-12-14 | 清华大学 | 一种多量测断面的电力系统不良支路参数辨识方法 |
Non-Patent Citations (2)
Title |
---|
杨滢等.集成于EMS中的参数估计软件的开发与应用.《电网技术》.2006,第30卷(第4期), |
集成于EMS中的参数估计软件的开发与应用;杨滢等;《电网技术》;20060228;第30卷(第4期);第43-49页 * |
Also Published As
Publication number | Publication date |
---|---|
CN102636706A (zh) | 2012-08-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103324847B (zh) | 电力系统动态不良数据检测与辨识方法 | |
CN102636706B (zh) | 一种电网中参数误差支路的辨识方法 | |
CN103944165B (zh) | 一种大电网参数辨识估计方法 | |
CN107016236B (zh) | 基于非线性量测方程的电网假数据注入攻击检测方法 | |
CN103326358B (zh) | 基于同步相角测量装置的电力系统动态状态估计方法 | |
CN104836223B (zh) | 电网参数错误与不良数据协同辨识与估计方法 | |
CN103454559B (zh) | 一种配电网单相接地故障区段定位方法及定位装置 | |
CN102175922B (zh) | 基于pmu量测数据的电力线路参数的辨识与估计方法 | |
CN110289613A (zh) | 基于灵敏度矩阵的配电网拓扑识别与线路参数辨识方法 | |
CN104483575A (zh) | 用于非侵入式电力监测的自适应负荷事件检测方法 | |
CN107370150B (zh) | 基于同步相量量测的电力系统状态估计不良数据处理方法 | |
CN113078630B (zh) | 一种基于实时量测数据的低压配电网拓扑辨识方法 | |
CN106786531B (zh) | 一种配电终端优化配置方法 | |
CN106127047B (zh) | 一种基于延森-香农距离的电力系统恶意数据检测方法 | |
CN103413044A (zh) | 一种基于变电站量测信息的电力系统局部拓扑估计方法 | |
CN102005757A (zh) | 一种基于量测生成树的电网拓扑错误的辨识方法 | |
CN107453484A (zh) | 一种基于wams信息的scada数据校准方法 | |
CN103066591B (zh) | 一种基于实时测量的电网参数偏差识别方法 | |
CN105406471A (zh) | 电网不良数据辨识与估计方法 | |
CN103324858A (zh) | 配电网三相潮流状态估计方法 | |
CN109193635A (zh) | 一种基于自适应稀疏回归方法的配电网拓扑结构重建方法 | |
CN105183938A (zh) | 电网不良数据辨识与估计方法 | |
CN104537581A (zh) | 采用模糊相似度匹配的电压暂降源在线定位方法 | |
CN105322539A (zh) | 一种配电网scada系统电压数据修正方法 | |
CN116736212B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140219 |
|
CF01 | Termination of patent right due to non-payment of annual fee |