CN104917176B - 一种基于耦合支路介数的电网关键支路识别方法 - Google Patents
一种基于耦合支路介数的电网关键支路识别方法 Download PDFInfo
- Publication number
- CN104917176B CN104917176B CN201510358198.5A CN201510358198A CN104917176B CN 104917176 B CN104917176 B CN 104917176B CN 201510358198 A CN201510358198 A CN 201510358198A CN 104917176 B CN104917176 B CN 104917176B
- Authority
- CN
- China
- Prior art keywords
- mtd
- mrow
- msub
- branch road
- mtr
- 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
- 238000010168 coupling process Methods 0.000 title claims abstract description 36
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 36
- 230000008878 coupling Effects 0.000 title claims abstract description 27
- 238000000034 method Methods 0.000 title claims abstract description 20
- 239000011159 matrix material Substances 0.000 claims abstract description 59
- 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 3
- 230000005540 biological transmission Effects 0.000 abstract description 6
- 238000007796 conventional method Methods 0.000 abstract description 2
- 230000008569 process Effects 0.000 abstract description 2
- 230000005611 electricity Effects 0.000 description 6
- 230000007812 deficiency Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 241000272525 Anas platyrhynchos Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000037361 pathway Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 150000003839 salts Chemical group 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
- Navigation (AREA)
Abstract
本发明提供一种基于耦合支路介数的电网关键支路识别方法,包括以下步骤:建立边权连接矩阵、支路潮流权重矩阵、节点‑支路关联向量及节点电抗矩阵;断开支路Lij,若电网解列,则将支路Lij列入第一关键支路子集;若电网未解列,则将支路Lij列入备选集;计算耦合支路介数,并对支路进行排序,对于耦合支路介数大于0.5的支路,列入第二关键支路子集;得到电网关键支路集。本发明采用Floyd算法求取最短路径,而且求取了传输方向,优于常规方法;耦合支路介数是综合评估指标,既考虑了有功功率的绝对大小,还考虑了功率裕度的相对值,并且结合了前后两级故障,所得结果更加符合电网实际情况;耦合支路介数能够充分反映连锁故障发展演化过程中支路的关键性影响。
Description
技术领域
本发明属于电力***安全稳定分析技术领域,具体涉及一种基于耦合支路介数的电网关键支路识别方法。
背景技术
近年来,国内外电网中连锁故障时有发生,均造成重大的经济损失和严重的社会影响。连锁故障一般是由单一或多重故障扰动诱发而形成,具体而言,电网正常运行时支路传输一定的有功功率,当其中某条或某几条支路被切除而停运后,停运支路上的潮流就会转移到其他支路上,如果其他正常运行支路不能承载这部分转移潮流,就会因过负荷而停运,导致连锁故障发生。以往针对连锁故障的研究表明:电网中某些关键支路对于连锁故障的产生和发展具有推动作用。如何识别电网中的关键支路,对于正确认识连锁故障发展演化机理具有重要意义。
以往针对电网关键支路的研究,所提出的反映支路关键性影响的指标,主要存在两点不足:
(1)仅假设电能沿节点对间支路电抗值之和最小的路径输送,并不符合电力***实际运行情况。
(2)仅考虑支路自身的关键性影响因素,未从连锁故障角度出发,考虑其他相关支路的影响。
如何克服以上两点不足,提出能够充分反映支路关键性影响的指标,是需要深入研究的问题。
发明内容
为了克服上述现有技术的不足,本发明提供一种基于耦合支路介数的电网关键支路识别方法,耦合支路介数不仅考虑当前运行方式下的各支路潮流权重,同时还考虑了待评价支路的最大相关性支路的影响,用于电网关键支路的识别。
为了实现上述发明目的,本发明采取如下技术方案:
本发明提供一种基于耦合支路介数的电网关键支路识别方法,所述方法包括以下步骤:
步骤1:建立边权连接矩阵BQ、支路潮流权重矩阵W、节点-支路关联向量M及节点电抗矩阵X;
步骤2:以边权连接矩阵BQ为对象,断开支路Lij,若电网解列,则将支路Lij列入第一关键支路子集;若电网未解列,则将支路Lij列入备选集;
步骤3:计算备选集中支路Lij的耦合支路介数,并按照耦合支路介数从大到小的顺序将支路进行排序,对于耦合支路介数大于0.5的支路,列入第二关键支路子集;
步骤4:将第一关键支路子集和第二关键支路子集合并,得到电网关键支路集。
所述步骤1包括以下步骤:
步骤1-1:建立以支路电抗值为权重的边权连接矩阵BQ,有:
式(1)中,n为节点总数,且有:
式(2)中,BQij表示边权连接矩阵BQ的第i行、第j列元素;BQji表示边权连接矩阵BQ的第j行、第i列元素;xij为支路Lij的电抗值,i,j∈[1,n];
步骤1-2:建立包含支路有功功率的支路潮流权重矩阵W,有:
式(3)中,有:
式(4)中,Wij表示支路潮流权重矩阵W的第i行、第j列元素;Wji表示支路潮流权重矩阵W的第j行、第i列元素;Pij表示支路Lij的有功功率,Pmax表示支路Lij的额定有功功率;
步骤1-3:建立节点-支路关联向量M,有:
1)若电网中的节点与支路Lij无关,则节点-支路关联向量M对应元素取零;
2)若支路Lij的有功功率从节点i流向节点j,则节点i在节点-支路关联向量M中对应元素取1,节点j在则节点-支路关联向量M中对应元素-1;
步骤1-4:建立节点电抗矩阵X,有:
式(6)中,若i≠j,Xij表示节点i、j之间的互阻抗;若i=j,Xii表示节点i的自阻抗。
所述步骤2具体包括以下步骤:
步骤2-1:以边权连接矩阵BQ为对象,对断开支路Lij,令BQij=BQji=∞,之后从节点1开始,查找与节点1之间连接的节点,进而查找与节点1连接的所有节点,其他节点以此类推;
步骤2-2:若节点i不能遍历与节点i连接的所有节点,则表明支路Lij断开后电网解列,于是将支路Lij列入第一关键支路子集;
步骤2-3:若节点i可以遍历与节点i连接的所有节点,则表明支路Lij断开后电网未解列,则将支路Lij列入备选集。
所述步骤3中,以边权连接矩阵BQ为对象,计算备选集中支路Lij的耦合支路介数Sij,具体包括:
1)建立节点连接关系矩阵R,有:
式(7)中,n为节点总数,且有:
式(8)中,Rij表示节点连接关系矩阵R的第i行、第j列元素,Rji表示节点连接关系矩阵R的第j行、第i列元素,i,j∈[1,n];
从节点i开始遍历,找到所有可能到达节点j的路径,得到最短路径Z,最短路径Z由多条支路串联构成;
2)在最短路径Z上寻找与由节点i到节点j有功功率流向相同、且潮流权重最大的支路Lrs,潮流权重由支路潮流权重矩阵W得到;计算支路Lrs、Lij之间的潮流转移因子Drs-ij,有:
式(9)中,xrs表示支路Lrs的电抗值,xij表示支路Lij的电抗值;Mrs表示支路Lrs的节点支路关联向量,Mij表示支路Lij的节点支路关联向量,X表示节点电抗矩阵;
3)备选集中支路Lij的耦合支路介数Sij表示为:
Sij=Wij·Wrs·Drs-ij (10)
式(10)中,Wij表示支路Lij的潮流权重,Wrs表示支路Lrs的潮流权重,Wij和Wrs由支路潮流权重矩阵W得到。
与现有技术相比,本发明的有益效果在于:
1)本发明提供的基于耦合支路介数的电网关键支路识别方法,耦合支路介数不仅考虑当前运行方式下的各支路潮流权重,还考虑待评价支路Lij的最大相关性支路的影响,潮流转移因子越大,支路Lrs传输有功功率越限的风险越大;
2)同时考虑到支路Lrs的运行状况,支路Lrs的潮流权重Wrs越大,受支路Lij影响发生故障的风险越大,对下一级故障的影响也越大;
3)采用Floyd算法求取最短路径,而且求取了传输方向,优于常规方法;
4)耦合支路介数是综合评估指标,既考虑了有功功率的绝对大小,还考虑了功率裕度的相对值,并且结合了前后两级故障,所得结果更加符合电网实际情况;
5)耦合支路介数能够充分反映连锁故障发展演化过程中支路的关键性影响。
附图说明
图1是本发明实施例中基于耦合支路介数的电网关键支路识别方法流程图;
图2是本发明实施例中宁夏电网结构示意图。
具体实施方式
下面对本发明作进一步详细说明。
如图1,本发明提供一种基于耦合支路介数的电网关键支路识别方法,所述方法包括以下步骤:
步骤1:建立边权连接矩阵BQ、支路潮流权重矩阵W、节点-支路关联向量M及节点电抗矩阵X;
步骤2:以边权连接矩阵BQ为对象,断开支路Lij,若电网解列,则将支路Lij列入第一关键支路子集;若电网未解列,则将支路Lij列入备选集;
步骤3:计算备选集中支路Lij的耦合支路介数,并按照耦合支路介数从大到小的顺序将支路进行排序,对于耦合支路介数大于0.5的支路,列入第二关键支路子集;
步骤4:将第一关键支路子集和第二关键支路子集合并,得到电网关键支路集。
所述步骤1包括以下步骤:
步骤1-1:建立以支路电抗值为权重的边权连接矩阵BQ,有:
式(1)中,n为节点总数,且有:
式(2)中,BQij表示边权连接矩阵BQ的第i行、第j列元素;BQji表示边权连接矩阵BQ的第j行、第i列元素;xij为支路Lij的电抗值,i,j∈[1,n];
步骤1-2:建立包含支路有功功率的支路潮流权重矩阵W,有:
式(3)中,有:
式(4)中,Wij表示支路潮流权重矩阵W的第i行、第j列元素;Wji表示支路潮流权重矩阵W的第j行、第i列元素;Pij表示支路Lij的有功功率,Pmax表示支路Lij的额定有功功率;
步骤1-3:建立节点-支路关联向量M,有:
1)若电网中的节点与支路Lij无关,则节点-支路关联向量M对应元素取零;
2)若支路Lij的有功功率从节点i流向节点j,则节点i在节点-支路关联向量M中对应元素取1,节点j在则节点-支路关联向量M中对应元素-1;
步骤1-4:建立节点电抗矩阵X,有:
式(6)中,若i≠j,Xij表示节点i、j之间的互阻抗;若i=j,Xii表示节点i的自阻抗。
所述步骤2具体包括以下步骤:
步骤2-1:以边权连接矩阵BQ为对象,对断开支路Lij,令BQij=BQji=∞,之后从节点1开始,查找与节点1之间连接的节点,进而查找与节点1连接的所有节点,其他节点以此类推;
步骤2-2:若节点i不能遍历与节点i连接的所有节点,则表明支路Lij断开后电网解列,于是将支路Lij列入第一关键支路子集;
步骤2-3:若节点i可以遍历与节点i连接的所有节点,则表明支路Lij断开后电网未解列,则将支路Lij列入备选集。
所述步骤3中,以边权连接矩阵BQ为对象,计算备选集中支路Lij的耦合支路介数Sij,具体包括:
1)建立节点连接关系矩阵R,有:
式(7)中,n为节点总数,且有:
式(8)中,Rij表示节点连接关系矩阵R的第i行、第j列元素,Rji表示节点连接关系矩阵R的第j行、第i列元素,i,j∈[1,n];
从节点i开始遍历,找到所有可能到达节点j的路径,得到最短路径Z,最短路径Z由多条支路串联构成;
2)在最短路径Z上寻找与由节点i到节点j有功功率流向相同、且潮流权重最大的支路Lrs,潮流权重由支路潮流权重矩阵W得到;计算支路Lrs、Lij之间的潮流转移因子Drs-ij,有:
式(9)中,xrs表示支路Lrs的电抗值,xij表示支路Lij的电抗值;Mrs表示支路Lrs的节点支路关联向量,Mij表示支路Lij的节点支路关联向量,X表示节点电抗矩阵;
3)备选集中支路Lij的耦合支路介数Sij表示
Sij=Wij·Wrs·Drs-ij (10)
式(10)中,Wij表示支路Lij的潮流权重,Wrs表示支路Lrs的潮流权重,Wij和Wrs由支路潮流权重矩阵W得到。
实施例
以宁夏电网结构(如图2)对本发明进一步详细说明:
一、先将220kV和330kV电压等级支路归算到750kV电压等级,并将各节点统一编号如表1所示:
表1
针对表1宁夏电网节点编号,分别建立边权连接矩阵、支路潮流权重矩阵;
二、对图2所示宁夏电网,采用BFS算法对所有支路进行扫描,检测***间是否存在唯一联络线;结果显示,宁夏电网不存在唯一联络线。
三、采用Floyd算法确定支路Lij故障后由节点i到节点j的最短路径。以银川东变压器输电支路为例,银川东750kV节点编号为4,银川东330kV节点编号为5,采用Floyd算法,得到银川东变压器故障后,最短路径为银川东750kV(4)-太阳山750kV(8)-太阳山330kV(9)-古峰330kV(39)-盐州330kV(40)-蒋家南330kV(41)-银川东330kV(5)。
四、通过比较分析,最短路径4-8-9-39-40-41-5中潮流权重最大的为支路40-41,即盐州330kV-蒋家南330kV输电支路,计算相应的潮流转移因子D40-41-4-5=1。
五、计算支路41-5的耦合支路介数,并归一化为:
S4,5=W4,5·W40,41·D40,41-4.5=1
依次求取宁夏电网各条支路的耦合支路介数,结果如表2所示,仅列出耦合支路介数较大的10条支路:
表2
支路 | 电压等/kV | 支路潮/MW | 耦合支路介数 | 最大相关性支路潮/MW |
L4-5银川东变压器 | \ | 2487.7 | S(4,5)=1.00 | L40-41=2842 |
L40-41盐州-蒋家南 | 330 | 1196.8 | S(40,41)=1.00 | L4-5=3660.8 |
L5-41银川东-蒋家南 | 330 | 1416.3 | S(5,41)=0.38 | L4-5=3902.6 |
L39-40古峰-盐州 | 330 | 483.6 | S(39,40)=0.13 | L4-5=2954.3 |
L18-31六盘山-清水河 | 330 | 1120.6 | S(18,31)=0.11 | L15-16=1619.6 |
L8-9太阳山变压器 | \ | 1393.8 | S(8,9)=0.07 | L15-16=1625.0 |
L8-10太阳山-枣泉 | 750 | 1106.3 | S(8,10)=0.07 | L8-11=2488.1 |
L8-11太阳山-鸳鸯二 | 750 | 1381.7 | S(8,11)=0.07 | L8-10=2486.7 |
L2-26沙湖-贺兰山 | 750 | 1259.2 | S(2,26)=0.06 | L8-15=2478.2 |
L9-29太阳山-罗山 | 330 | 671.0 | S(9,29)=0.05 | L15-16=1449.4 |
由表2结果可知,排名前三的分别为银川东750kV站变压器支路、盐州-蒋家南330kV输电支路、蒋家南-银川东的330kV输电支路,均位于银川东站附近,可以称之为“关键区域”,在给定运行方式下应关注该地区各电气元件的运行状态。
最后应当说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,所属领域的普通技术人员参照上述实施例依然可以对本发明的具体实施方式进行修改或者等同替换,这些未脱离本发明精神和范围的任何修改或者等同替换,均在申请待批的本发明的权利要求保护范围之内。
Claims (3)
1.一种基于耦合支路介数的电网关键支路识别方法,其特征在于:所述方法包括以下步骤:
步骤1:建立边权连接矩阵BQ、支路潮流权重矩阵W、节点-支路关联向量M及节点电抗矩阵X;
步骤2:以边权连接矩阵BQ为对象,断开支路Lij,若电网解列,则将支路Lij列入第一关键支路子集;若电网未解列,则将支路Lij列入备选集;
步骤3:计算备选集中支路Lij的耦合支路介数,并按照耦合支路介数从大到小的顺序将支路进行排序,对于耦合支路介数大于0.5的支路,列入第二关键支路子集;
步骤4:将第一关键支路子集和第二关键支路子集合并,得到电网关键支路集;
所述步骤1包括以下步骤:
步骤1-1:建立以支路电抗值为权重的边权连接矩阵BQ,有:
<mrow>
<mi>B</mi>
<mi>Q</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>BQ</mi>
<mn>11</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>BQ</mi>
<mn>12</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>BQ</mi>
<mn>13</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>BQ</mi>
<mrow>
<mn>1</mn>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>BQ</mi>
<mn>21</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>BQ</mi>
<mn>22</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>BQ</mi>
<mrow>
<mn>2</mn>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>BQ</mi>
<mrow>
<mi>n</mi>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>BQ</mi>
<mrow>
<mi>n</mi>
<mn>2</mn>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>BQ</mi>
<mrow>
<mi>n</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式(1)中,n为节点总数,且有:
式(2)中,BQij表示边权连接矩阵BQ的第i行、第j列元素;BQji表示边权连接矩阵BQ的第j行、第i列元素;xij为支路Lij的电抗值,i,j∈[1,n];
步骤1-2:建立包含支路有功功率的支路潮流权重矩阵W,有:
<mrow>
<mi>W</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>W</mi>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>W</mi>
<mn>12</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>W</mi>
<mn>13</mn>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>W</mi>
<mrow>
<mn>1</mn>
<mi>n</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>W</mi>
<mn>21</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>W</mi>
<mn>22</mn>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>W</mi>
<mrow>
<mn>2</mn>
<mi>n</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>W</mi>
<mrow>
<mi>n</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>W</mi>
<mrow>
<mi>n</mi>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>W</mi>
<mrow>
<mi>n</mi>
<mi>n</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
式(3)中,有:
式(4)中,Wij表示支路潮流权重矩阵W的第i行、第j列元素;Wji表示支路潮流权重矩阵W的第j行、第i列元素;Pij表示支路Lij的有功功率,Pmax表示支路Lij的额定有功功率;
步骤1-3:建立节点-支路关联向量M,有:
1)若电网中的节点与支路Lij无关,则节点-支路关联向量M对应元素取零;
2)若支路Lij的有功功率从节点i流向节点j,则节点i在节点-支路关联向量M中对应元素取1,节点j在则节点-支路关联向量M中对应元素-1;
步骤1-4:建立节点电抗矩阵X,有:
<mrow>
<mi>X</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>X</mi>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>X</mi>
<mn>12</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>X</mi>
<mn>13</mn>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>X</mi>
<mrow>
<mn>1</mn>
<mi>n</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>X</mi>
<mn>21</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>X</mi>
<mn>22</mn>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>X</mi>
<mrow>
<mn>2</mn>
<mi>n</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>X</mi>
<mrow>
<mi>n</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>X</mi>
<mrow>
<mi>n</mi>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>X</mi>
<mrow>
<mi>n</mi>
<mi>n</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
式(6)中,若i≠j,Xij表示节点i、j之间的互阻抗;若i=j,Xii表示节点i的自阻抗。
2.根据权利要求1所述的基于耦合支路介数的电网关键支路识别方法,其特征在于:所述步骤2具体包括以下步骤:
步骤2-1:以边权连接矩阵BQ为对象,对断开支路Lij,令BQij=BQji=∞,之后从节点1开始,查找与节点1之间连接的节点,进而查找与节点1连接的所有节点,其他节点以此类推;
步骤2-2:若节点i不能遍历与节点i连接的所有节点,则表明支路Lij断开后电网解列,于是将支路Lij列入第一关键支路子集;
步骤2-3:若节点i可以遍历与节点i连接的所有节点,则表明支路Lij断开后电网未解列,则将支路Lij列入备选集。
3.根据权利要求1所述的基于耦合支路介数的电网关键支路识别方法,其特征在于:所述步骤3中,以边权连接矩阵BQ为对象,计算备选集中支路Lij的耦合支路介数Sij,具体包括:
1)建立节点连接关系矩阵R,有:
<mrow>
<mi>R</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>R</mi>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>R</mi>
<mn>12</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>R</mi>
<mn>13</mn>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>R</mi>
<mrow>
<mn>1</mn>
<mi>n</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>R</mi>
<mn>21</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>R</mi>
<mn>22</mn>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>R</mi>
<mrow>
<mn>2</mn>
<mi>n</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>R</mi>
<mrow>
<mi>n</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>R</mi>
<mrow>
<mi>n</mi>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>R</mi>
<mrow>
<mi>n</mi>
<mi>n</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
式(7)中,n为节点总数,且有:
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>j</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mi>j</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>j</mi>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mo>&NotEqual;</mo>
<mi>j</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
式(8)中,Rij表示节点连接关系矩阵R的第i行、第j列元素,Rji表示节点连接关系矩阵R的第j行、第i列元素,i,j∈[1,n];
从节点i开始遍历,找到所有可能到达节点j的路径,得到最短路径Z,最短路径Z由多条支路串联构成;
2)在最短路径Z上寻找与由节点i到节点j有功功率流向相同、且潮流权重最大的支路Lrs,潮流权重由支路潮流权重矩阵W得到;计算支路Lrs、Lij之间的潮流转移因子Drs-ij,有:
<mrow>
<msub>
<mi>D</mi>
<mrow>
<mi>r</mi>
<mi>s</mi>
<mo>-</mo>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>X</mi>
<mrow>
<mi>r</mi>
<mi>s</mi>
<mo>-</mo>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>/</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>r</mi>
<mi>s</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
<mo>-</mo>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>/</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
式(9)中,xrs表示支路Lrs的电抗值,xij表示支路Lij的电抗值;Mrs表示支路Lrs的节点支路关联向量,Mij表示支路Lij的节点支路关联向量,X表示节点电抗矩阵;
3)备选集中支路Lij的耦合支路介数Sij表示为:
Sij=Wij·Wrs·Drs-ij (10)
式(10)中,Wij表示支路Lij的潮流权重,Wrs表示支路Lrs的潮流权重,Wij和Wrs由支路潮流权重矩阵W得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510358198.5A CN104917176B (zh) | 2015-06-25 | 2015-06-25 | 一种基于耦合支路介数的电网关键支路识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510358198.5A CN104917176B (zh) | 2015-06-25 | 2015-06-25 | 一种基于耦合支路介数的电网关键支路识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104917176A CN104917176A (zh) | 2015-09-16 |
CN104917176B true CN104917176B (zh) | 2018-01-19 |
Family
ID=54085997
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510358198.5A Active CN104917176B (zh) | 2015-06-25 | 2015-06-25 | 一种基于耦合支路介数的电网关键支路识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104917176B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105335609B (zh) * | 2015-10-20 | 2018-06-29 | 国网天津市电力公司 | 一种评估送电断面功率转移比的方法 |
CN105303460A (zh) * | 2015-10-30 | 2016-02-03 | 国家电网公司 | 一种电网关键节点和关键支路的辨识方法 |
CN108667027B (zh) * | 2018-05-15 | 2020-04-24 | 四川大学 | 含柔性直流的交直流***潮流转移搜索与量化分析方法 |
CN109257207B (zh) * | 2018-08-28 | 2021-10-15 | 浙江工业大学 | 一种基于耦合链接介数差异指标的提高相互依存网络鲁棒性的方法 |
CN109449931B (zh) * | 2018-11-23 | 2021-11-16 | 中国电力科学研究院有限公司 | 结合加权潮流熵和全过程仿真的连锁故障集构建方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1976160A (zh) * | 2006-12-08 | 2007-06-06 | 浙江大学 | 一种大型电力***脆弱线路辨识的方法 |
CN103714398A (zh) * | 2013-12-30 | 2014-04-09 | 国家电网公司 | 一种基于抗灾变性指标扩展配电网的装置及方法 |
CN104156769A (zh) * | 2013-05-31 | 2014-11-19 | 贵州电网公司电力调度控制中心 | 电力***脆弱性评估方法 |
CN104466959A (zh) * | 2014-12-31 | 2015-03-25 | 广东电网有限责任公司电力调度控制中心 | 电力***关键线路辨识方法和*** |
-
2015
- 2015-06-25 CN CN201510358198.5A patent/CN104917176B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1976160A (zh) * | 2006-12-08 | 2007-06-06 | 浙江大学 | 一种大型电力***脆弱线路辨识的方法 |
CN104156769A (zh) * | 2013-05-31 | 2014-11-19 | 贵州电网公司电力调度控制中心 | 电力***脆弱性评估方法 |
CN103714398A (zh) * | 2013-12-30 | 2014-04-09 | 国家电网公司 | 一种基于抗灾变性指标扩展配电网的装置及方法 |
CN104466959A (zh) * | 2014-12-31 | 2015-03-25 | 广东电网有限责任公司电力调度控制中心 | 电力***关键线路辨识方法和*** |
Also Published As
Publication number | Publication date |
---|---|
CN104917176A (zh) | 2015-09-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104917176B (zh) | 一种基于耦合支路介数的电网关键支路识别方法 | |
CN105912835B (zh) | 一种计及风电场接入的电力线路脆弱性辨识模型建立方法 | |
CN106529791B (zh) | 一种电力***中支路重要度的评估方法 | |
CN110034581A (zh) | 风电并网下电力***的区间电气介数脆弱性评估方法 | |
CN107947151A (zh) | 一种电力***关键节点的识别方法 | |
CN104319785B (zh) | 一种基于源流路径电气剖分的风电***关键节点识别方法 | |
CN105048521B (zh) | 一种考虑多约束的受端电网直流接入能力计算方法 | |
CN103954885B (zh) | 基于分布参数的双回线故障单端定位***及定位方法 | |
CN104901308B (zh) | 电力***关键线路辨识方法 | |
CN104849614B (zh) | 小电流单相接地故障选线的判断方法 | |
CN107272645B (zh) | 中性点经电阻接地的光伏电站并网故障模型与分析方法 | |
CN111431173A (zh) | 考虑负荷互补特性的中压配电网环网结构确定方法及*** | |
CN105356494B (zh) | 多端vsc-hvdc并网***的可靠性计算方法 | |
CN104821578A (zh) | 计及可用输电能力的含大规模风电的输电***规划方法 | |
CN106779269A (zh) | 基于变异系数法的智能变电站二次***可靠性评估方法 | |
CN110350522A (zh) | 一种基于加权h指数的电力***脆弱线路辨识方法 | |
CN103904665B (zh) | 一种适用于双馈风机的电力***稳定器设计方法 | |
CN109375050A (zh) | 一种基于初始故障线路割集的电力***关键断面辨识方法 | |
CN102930140B (zh) | 一种同塔多回复杂输电线路传输方程的求解方法 | |
CN109768543B (zh) | 一种基于混合整数线性规划的弹性保底网架搜索建模方法 | |
CN107196307B (zh) | 一种变压器故障后电网有功潮流快速估算的方法 | |
CN104732302A (zh) | 一种基于免疫算法的多级电磁环网优化解环方法 | |
CN112052543B (zh) | 一种基于混合整数二阶锥规划的保底网架搜索建模方法 | |
CN105552880A (zh) | 基于状态枚举法的电力***典型故障集确定方法 | |
CN106777597A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |