CN106374462A - 一种基于单时刻量测数据的戴维南等值参数解析的电力系统运行状态的监测方法 - Google Patents
一种基于单时刻量测数据的戴维南等值参数解析的电力系统运行状态的监测方法 Download PDFInfo
- Publication number
- CN106374462A CN106374462A CN201610877210.8A CN201610877210A CN106374462A CN 106374462 A CN106374462 A CN 106374462A CN 201610877210 A CN201610877210 A CN 201610877210A CN 106374462 A CN106374462 A CN 106374462A
- Authority
- CN
- China
- Prior art keywords
- delta
- node
- matrix
- formula
- thevenin
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 58
- 238000012544 monitoring process Methods 0.000 title claims abstract description 19
- 239000011159 matrix material Substances 0.000 claims abstract description 53
- 238000004458 analytical method Methods 0.000 claims description 12
- 238000005259 measurement Methods 0.000 claims description 11
- 230000005611 electricity Effects 0.000 claims description 5
- 241000208340 Araliaceae Species 0.000 claims description 4
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims description 4
- 235000003140 Panax quinquefolius Nutrition 0.000 claims description 4
- 235000008434 ginseng Nutrition 0.000 claims description 4
- 238000002347 injection Methods 0.000 claims description 4
- 239000007924 injection 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 2
- 238000004364 calculation method Methods 0.000 abstract description 16
- 238000000205 computational method Methods 0.000 description 5
- 230000008901 benefit Effects 0.000 description 3
- 238000005070 sampling Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002360 preparation method 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
-
- 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)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
一种基于单时刻量测数据的戴维南等值参数解析的电力系统运行状态的监测方法,本发明涉及基于单时刻量测数据的戴维南等值参数解析的电力系统运行状态的监测方法。本发明的目的是为了解决现有戴维南等值参数计算方法精度低和计算量大的问题。具体步骤为:步骤一、根据实时的网络拓扑结构情况,从数据库里在线调用相应的网络阻抗矩阵参数;步骤二、利用各个负荷节点和PV节点的电压和电流量测信息以及网络阻抗矩阵参数计算矩阵C;步骤三、对矩阵C求逆,并根据解析公式计算出相应节点的戴维南等值参数。本发明用于电力系统运行状态的监测。
Description
技术领域
本发明涉及基于单时刻量测数据的戴维南等值参数解析的电力系统运行状态的监测方法。
背景技术
近些年,由于电力负荷的快速增长,在世界范围内,电力网络逐渐运行于安全域的边缘,因此由电压崩溃而引起的大面积停电事故时有发生,造成了很大的经济损失。随着间歇性新能源的大量接入,电网运行环境将更加恶化,这一趋势引起了各国政府和电力公司的高度重视。为了保证电力系统的安全稳定运行,为用户提供优质可靠的电能,电力调度中心需要实时监控电力系统的运行情况。在各种电力系统的运行状态在线辨评估方法中,基于戴维南等值参数的阻抗模比值法因其简单快速而成为主要方法之一。因此,研究和发明一种有效的戴维南等值参数计算方法,使基于戴维南等值的阻抗模比值法能够更好地为电力系统的运行状态监视和控制提供服务,是实现电力系统安全稳定运行的一项重要-任务。
目前为止,戴维南等值参数计算方法已经引起电力行业的高度重视,其方法主要有以下几种:
1)基于多时刻量测数据的计算方法:该方法假定在多个采样时刻内(或在离线状态下的多个系统仿真状态)系统的等值参数不变,通过等值点的电压变化与电流变化的比值计算等值参数(使用开路电压和短路电流计算等值参数的方法是这一类方法中离线计算的一个特例)。这类方法一般需要保证在采样过程中,系统侧波动远远小于负荷侧波动,且两者均不会对系统运行状态产生明显影响。实测数据很难满足这些要求,而离线系统仿真的数据虽然可以满足要求,但需要多次潮流计算,在大规模电网的监控上有一定难度。
2)基于单时刻量测数据的计算方法:该方法使用单一时刻下的量测数据通过解析的方法计算戴维南等值参数。但目前为止,这一类方法在解析过程中都对公式有一定程度的近似处理,在精度上仍有提高空间。
发明内容
本发明的目的是为了解决现有戴维南等值参数计算方法精度低和计算量大的问题,而提出一种基于单时刻量测数据的戴维南等值参数解析的电力系统运行状态的监测方法。
一种基于单时刻量测数据的戴维南等值参数解析的电力系统运行状态的监测方法具体步骤为:
步骤一、根据实时的网络拓扑结构情况,从数据库里在线调用相应的网络阻抗矩阵参数;
步骤二、利用各个负荷节点和PV节点的电压和电流量测信息以及网络阻抗矩阵参数计算矩阵C;
步骤三、对矩阵C求逆,并根据解析公式计算出相应节点的戴维南等值参数。
本发明的有益效果为:
相较于多时刻量测的方法(在线量测),本方法不受网络侧扰动影响,可以准确的计算出戴维南等值参数;相较于多时刻量测的方法(离线潮流计算),本方法在对电网内多个负荷节点计算戴维南等值参数的情况下,可以节约大量的潮流计算量。以3120节点算例为例,对3120节点网络内的2314个负荷节点计算戴维南等值参数,本方法可以节约近97%的计算时间。
本发明所提出的方法具有明显的解析性,不需要进行多时刻下的数据量测(或者潮流计算)即可实现戴维南等值阻抗的快速计算,有效解决了传统方法的技术性难题;直接利用网络阻抗矩阵数据和单时间断面的节点量测信息,即可实现戴维南等值阻抗的快速计算,具有原理明晰、计算结果准确、计算速度快和操作简单等优点;而且影响戴维南等值参数的因素也可以解析地被识别。大量仿真结果表明:本发明所提出的方法计算速度快、线性度好、计算结果准确;与被公认的基于多时刻量测数据的方法相比,本发明所提出的方法不受网络扰动影响,能够快速地、解析地计算各个节点的戴维南等值参数。
表2给出了不同系统运行状态下本发明方法与多时刻量测数据计算方法的结果。通过对比以看出,本发明所提出的方法计算结果快速准确。与基于多时刻量测数据计算方法相比,本发明所提出的方法,避免了对多时刻量测数据的苛刻要求,可以解析地计算各节点的戴维南等值参数,具有概念明晰、原理简单、结果准确性高、计算速度快等优点。
附图说明
图1是本发明流程图;
图2是本发明具体实施方式例所用IEEE-4节点系统的拓扑结构图;
图3是IEEE-4节点网络导纳矩阵。
具体实施方式
具体实施方式一:结合图1说明本实施方式,本实施方式的一种基于单时刻量测数据的戴维南等值参数解析的电力系统运行状态的监测方法具体步骤为:
步骤一、根据实时的网络拓扑结构情况,从数据库里在线调用相应的网络阻抗矩阵参数;
步骤二、利用各个负荷节点和PV节点的电压和电流量测信息以及网络阻抗矩阵参数计算矩阵C;
步骤三、对矩阵C求逆,并根据解析公式计算出相应节点的戴维南等值参数。
具体实施方式二:本实施方式与具体实施方式一不同的是:所述步骤二中利用各个负荷节点和PV节点的电压和功率量测信息以及网络阻抗矩阵参数计算矩阵C;具体过程为:
利用单一时刻下各个负荷节点和PV节点的电压和电流量测信息以及网络阻抗矩阵参数计算矩阵C,矩阵C是由和分别对和的偏导数子矩阵所构成,矩阵C如公式(1)所示
式中,AL为负荷节点集合;AG为未饱和发电机节点集合;AGΤ饱和发电机节点与联络节点集合;为负荷节点集合的有功矩阵;为负荷节点集合的无功矩阵;为未饱和发电节点集合的有功矩阵;为饱和发电机节点与联络节点集合的有功矩阵;为饱和发电机节点与联络节点集合的无功矩阵;为负荷节点集合的电压幅值;为负荷节点集合的电压相角;为未饱和发电节点集合的电压相角;为饱和发电机节点与联络节点集合的电压幅值;为饱和发电机节点与联络节点集合的电压相角;
公式(1)的子矩阵中各个元素可由公式(2)-(5)计算得出
式中,Pi、Qi、Ui、δi分别为节点i处的有功功率、无功功率、电压幅值和电压相角;Gij和Bij分别为节点i与节点j之间线路的互电导和互电纳;Gii和Bii分别为节点i的自电导和自电纳;Uj为节点j处的电压幅值;δij为节点i与节点j之间线路的电压相角;δj为节点j处的电压相角;i为第i个节点;j为第j个节点;i∈{AL,AG,AGΤ};j∈{AL,AG,AGΤ}。
其它步骤及参数与具体实施方式一相同。
具体实施方式三:本实施方式与具体实施方式一或二不同的是:所述步骤三中对矩阵C求逆,具体过程为:
对矩阵C求逆,得到如公式(6)所示的关系
式中和为零矩阵,则和的表达式改写为如公式(7)和(8)所示的形式:
其中,和分别是C-1的子矩阵的对角矩阵。
其它步骤及参数与具体实施方式一或二相同。
具体实施方式四:本实施方式与具体实施方式一至三之一不同的是:所述步骤三中根据解析公式计算出相应节点的戴维南等值参数;具体过程为:
通过公式(7)-(8)得到如公式(9)所示的计算戴维南等值阻抗的解析公式
式中,m为虚数单位;A、B、kL为中间参数。
其它步骤及参数与具体实施方式一至三之一相同。
具体实施方式五:本实施方式与具体实施方式一至四之一不同的是:所述
则戴维南等值参数通过公式(10)计算得出
式中,为负荷节点集合的注入电流相量;为负荷节点集合的电压幅值相量;
其它步骤及参数与具体实施方式一至四之一相同。
采用以下实施例验证本发明的有益效果:
实施例一:
本实施例一种基于单时刻量测数据的戴维南等值参数解析的电力系统运行状态的监测方法具体是按照以下步骤制备的:
以IEEE 4节点系统为例进行仿真和分析,具体步骤和结果如下:
1)假设IEEE 4节点系统在各种拓扑结构下的网络阻抗矩阵已经离线得到,本发明只需根据当前运行的网络拓扑结构(网络结构完好或有个别线路停运)情况,直接从数据库中调用相应的网络阻抗矩阵数据,因此节余了其它方法需要计算网络矩阵或雅可比矩阵的计算时间和计算量。本实施例子所使用的IEEE 4节点系统拓扑结构与网络导纳矩阵数据分别如图2和图3所示。
2)利用潮流计算的数据来代替各负荷节点的实时量测数据,包括各重负荷节点的电压和电流量测信息,IEEE 4节点系统各个节点数据如表1所示,表1是利用潮流计算的IEEE-4节点系统各节点数据;。利用公式(1)-(5)计算矩阵C。
表1
节点号 | 有功注入 | 无功注入 | 电压幅值 | 电压相角 |
1 | -30 | -18 | 0.983826 | -2.40308 |
2 | -55 | -13 | 1.026451 | -5.29563 |
3 | 50 | 18.93217 | 1.1 | 0.326293 |
4 | 36.70733 | 9.938253 | 1.05 | 0 |
3)对矩阵C求逆,提取负荷节点子矩阵的对角项,通过公式(6)-(10)计算得到各个负荷节点的戴维南等值参数,其结果如表2所示,表2是本专利中方法与基于多时刻数据方法得出的戴维南等值参数;
表2
以上结果与基于多时刻量测数据计算方法的结果相吻合。为了进一步证实本发明所提方法的有效性,表2给出了不同系统运行状态下本发明方法与多时刻量测数据计算方法的结果。通过对比以看出,本发明所提出的方法计算结果快速准确。与基于多时刻量测数据计算方法相比,本发明所提出的方法,避免了对多时刻量测数据的苛刻要求,可以解析地计算各节点的戴维南等值参数,具有概念明晰、原理简单、结果准确性高、计算速度快等优点。
本发明还可有其它多种实施例,在不背离本发明精神及其实质的情况下,本领域技术人员当可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。
Claims (5)
1.一种基于单时刻量测数据的戴维南等值参数解析的电力系统运行状态的监测方法,其特征在于:一种基于单时刻量测数据的戴维南等值参数解析的电力系统运行状态的监测方法具体步骤为:
步骤一、根据实时的网络拓扑结构情况,从数据库里在线调用相应的网络阻抗矩阵参数;
步骤二、利用各个负荷节点和PV节点的电压和电流量测信息以及网络阻抗矩阵参数计算矩阵C;
步骤三、对矩阵C求逆,并根据解析公式计算出相应节点的戴维南等值参数。
2.根据权利要求1所述一种基于单时刻量测数据的戴维南等值参数解析的电力系统运行状态的监测方法,其特征在于:所述步骤二中利用各个负荷节点和PV节点的电压和功率量测信息以及网络阻抗矩阵参数计算矩阵C;具体过程为:
利用单一时刻下各个负荷节点和PV节点的电压和电流量测信息以及网络阻抗矩阵参数计算矩阵C,矩阵C是由和分别对和的偏导数子矩阵所构成,矩阵C如公式(1)所示
式中,AL为负荷节点集合;AG为未饱和发电机节点集合;AGΤ饱和发电机节点与联络节点集合;为负荷节点集合的有功矩阵;为负荷节点集合的无功矩阵;为未饱和发电节点集合的有功矩阵;为饱和发电机节点与联络节点集合的有功矩阵;为饱和发电机节点与联络节点集合的无功矩阵;为负荷节点集合的电压幅值;为负荷节点集合的电压相角;为未饱和发电节点集合的电压相角;为饱和发电机节点与联络节点集合的电压幅值;为饱和发电机节点与联络节点集合的电压相角;
公式(1)的子矩阵中各个元素可由公式(2)-(5)计算得出
式中,Pi、Qi、Ui、δi分别为节点i处的有功功率、无功功率、电压幅值和电压相角;Gij和Bij分别为节点i与节点j之间线路的互电导和互电纳;Gii和Bii分别为节点i的自电导和自电纳;Uj为节点j处的电压幅值;δij为节点i与节点j之间线路的电压相角;δj为节点j处的电压相角;i为第i个节点;j为第j个节点;i∈{AL,AG,AGΤ};j∈{AL,AG,AGΤ}。
3.根据权利要求2所述一种基于单时刻量测数据的戴维南等值参数解析的电力系统运行状态的监测方法,其特征在于:所述步骤三中对矩阵C求逆,具体过程为:
对矩阵C求逆,得到如公式(6)所示的关系
式中和为零矩阵,则和的表达式改写为如公式(7)和(8)所示的形式:
其中,和分别是C-1的子矩阵的对角矩阵。
4.根据权利要求3所述一种基于单时刻量测数据的戴维南等值参数解析的电力系统运行状态的监测方法,其特征在于:所述步骤三中根据解析公式计算出相应节点的戴维南等值参数;具体过程为:
通过公式(7)-(8)得到如公式(9)所示的计算戴维南等值阻抗的解析公式
式中,m为虚数单位;A、B、kL为中间参数。
5.根据权利要求4所述一种基于单时刻量测数据的戴维南等值参数解析的电力系统运行状态的监测方法,其特征在于:所述
则戴维南等值参数通过公式(10)计算得出
式中,为负荷节点集合的注入电流相量;为负荷节点集合的电压幅值相量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610877210.8A CN106374462B (zh) | 2016-10-08 | 2016-10-08 | 一种基于单时刻量测数据的戴维南等值参数解析的电力系统运行状态的监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610877210.8A CN106374462B (zh) | 2016-10-08 | 2016-10-08 | 一种基于单时刻量测数据的戴维南等值参数解析的电力系统运行状态的监测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106374462A true CN106374462A (zh) | 2017-02-01 |
CN106374462B CN106374462B (zh) | 2018-07-20 |
Family
ID=57895546
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610877210.8A Expired - Fee Related CN106374462B (zh) | 2016-10-08 | 2016-10-08 | 一种基于单时刻量测数据的戴维南等值参数解析的电力系统运行状态的监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106374462B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103647284A (zh) * | 2013-10-22 | 2014-03-19 | 上海交通大学 | 一种解决单时间断面问题的电压稳定预测方法 |
CN104778367A (zh) * | 2015-04-20 | 2015-07-15 | 山东大学 | 基于单一状态断面的广域戴维南等值参数在线计算方法 |
CN104917177A (zh) * | 2015-06-26 | 2015-09-16 | 四川大学 | 一种电力系统戴维南等值参数在线识别方法及其装置 |
-
2016
- 2016-10-08 CN CN201610877210.8A patent/CN106374462B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103647284A (zh) * | 2013-10-22 | 2014-03-19 | 上海交通大学 | 一种解决单时间断面问题的电压稳定预测方法 |
CN104778367A (zh) * | 2015-04-20 | 2015-07-15 | 山东大学 | 基于单一状态断面的广域戴维南等值参数在线计算方法 |
CN104917177A (zh) * | 2015-06-26 | 2015-09-16 | 四川大学 | 一种电力系统戴维南等值参数在线识别方法及其装置 |
Non-Patent Citations (2)
Title |
---|
付正鑫: "基于广域测量技术戴维南等值参数的计算及在电压稳定性分析中的应用", 《中国优秀硕士学位论文全文数据库(电子期刊)工程科技Ⅱ辑》 * |
罗华伟等: "电网戴维南等值参数的快速计算", 《中国电机工程学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN106374462B (zh) | 2018-07-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Usman et al. | Validation of a PMU‐based fault location identification method for smart distribution network with photovoltaics using real‐time data | |
CN101383511B (zh) | 基于数据采集系统量测数据的电力系统状态估计方法 | |
CN101599643B (zh) | 一种基于指数型目标函数的电力系统抗差状态估计方法 | |
CN103248043B (zh) | 一种基于同步相角测量装置的电力系统多区域分布式状态估计方法 | |
CN102157938B (zh) | 电力系统电压稳定薄弱节点在线识别方法 | |
CN110417011A (zh) | 一种基于互信息与迭代随机森林的在线动态安全评估方法 | |
CN104007321B (zh) | 一种无功补偿电容器组的参数在线辨识方法 | |
Tripathy et al. | Sparse S‐transform for location of faults on transmission lines operating with unified power flow controller | |
WO2017016019A1 (zh) | 一种电力系统广域同调辨识方法及其装置 | |
CN112018784B (zh) | 一种基于同步相量测量数据的次同步谐振溯源方法 | |
CN104749488A (zh) | 一种基于连续数据窗的直流线路时域故障测距方法 | |
CN106228459A (zh) | 基于蒙特卡洛的等值可靠性评估方法 | |
CN106203889A (zh) | 非健全信息下基于最小二乘拟合法的配电网损耗评估方法 | |
CN107167689B (zh) | 考虑联络变压器变比的概率短路灵敏度计算方法 | |
CN109901022A (zh) | 基于同步量测数据的配电网区域定位方法 | |
CN103972889B (zh) | 一种配电线路阻抗在线辨识方法 | |
Guillen et al. | Data‐driven short‐circuit detection and location in microgrids using micro‐synchrophasors | |
CN105203877B (zh) | 能够消除量测误差影响的单回输电线路零序参数辨识方法 | |
CN104240036B (zh) | 一种基于临界系统等效阻抗的暂态电压稳定量化评估方法 | |
CN106374462A (zh) | 一种基于单时刻量测数据的戴维南等值参数解析的电力系统运行状态的监测方法 | |
CN104934969B (zh) | 一种电力线路参数的计算方法 | |
Gunadin et al. | Steady state stability assessment using extreme learning machine based on modal analysis | |
CN110007195B (zh) | 基于电流不平衡量的scada不良数据剔除方法及系统 | |
CN108089058A (zh) | 一种抗差配电网线路正序参数的辨识方法 | |
CN102798751A (zh) | 一种新型电压稳定性检测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180720 |