CN109599870A - 一种抑制短路电流的联络线电抗快速计算方法 - Google Patents
一种抑制短路电流的联络线电抗快速计算方法 Download PDFInfo
- Publication number
- CN109599870A CN109599870A CN201811356521.5A CN201811356521A CN109599870A CN 109599870 A CN109599870 A CN 109599870A CN 201811356521 A CN201811356521 A CN 201811356521A CN 109599870 A CN109599870 A CN 109599870A
- Authority
- CN
- China
- Prior art keywords
- node
- branch
- circuit current
- network
- short circuit
- 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
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
- H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了属于电力系统调度自动化技术领域的一种抑制短路电流的联络线电抗快速计算方法。该方法包括以下步骤:计算电力系统中节点的短路电流,选取短路水平过高需要进行短路电流限制的节点;以该节点作为支路一端的节点,增加一条新的支路,在原有系统的节点电抗矩阵的基础上利用支路追加法形成新网络的节点电抗矩阵;在故障分量网络中计算新网络下该节点的短路电流,在追加支路中获得限制该节点短路电流的最优支路。
Description
技术领域
本发明属于电力系统调度自动化技术领域,具体涉及一种抑制短路电流的联络线电抗快速计算方法。
背景技术
城市电网在500kV网架不断发展和完善的推动下,逐步形成了主网为500kV环网,区域电网为220kV电网的网络构造。在此基础上,推进实行以500kV变电站为核心的电网分区运行。城市电网采用分区运行方式后,能够大幅提高系统输电能力,有效控制并降低网损,便于系统调整运行方式以适应电网发展,降低事故预案编制的难度。可见,实施分区方案是城市电网发展的必经阶段,可以显著促进电网的安全稳定运行。但是,对于具有非独立供电可靠性的分区,当分区内发生故障时如主变N-2等,分区内的主变可能会发生过载情况,需要考虑投入分区联络线以获取其他分区的功率支援,此时,由于联络线的投入会导致分区互联,分区内的厂站短路电流会上升,有可能超出该电压等级下的遮断电流。为了满足系统安全稳定运行的要求,需要采取一定的短路电流限制措施。
发明内容
本发明的目的是提出一种抑制短路电流的联络线电抗快速计算方法,其特征在于,包括以下步骤:
步骤A.基于系统网络参数形成节点导纳矩阵,对节点导纳矩阵求逆得到描述网络的节点阻抗矩阵,求的系统各节点的短路电流,获取短路水平最高的节点;
步骤B.基于以上计算结果,由该节点引入一条支路,支路电抗为Xk,利用支路追加法推导新网络的节点电抗矩阵,计算此时该节点的短路电流,在能够选支路中获得最优支路。
所述步骤A包括:
步骤A1:任一网络用节点阻抗矩阵表示的节点电压方程为U=ZI,I和U分别为:用矩阵表示如下:
式(1)中,电压相量U为网络各节点对地电压;电流相量I为网络外部向各节点的注入电流;系数矩阵即节点阻抗矩阵Z;
上述矩阵的矩阵对角元素Zii称为自阻抗,表达为:
即当除节点i外所有其他节点注入电流均为零时,节点i电压与注入电流之比,也就是从节点i看进网络的等效阻抗;
非对角元素Zij称为互阻抗,表示为:
即当除节点i外所有其他节点注入电流均为零时,节点j电压与节点i电流之比,或者说节点i单位电流在节点j引起的电压值;显然Zij是不会为零的,即节点阻抗矩阵总是满阵;
步骤A2:基于步骤A1形成的节点阻抗矩阵,假设某一点f发生短路故障时,短路点f的接地阻抗为Zf=0;该网络中只有短路点f有注入电流-if,即if由f点流向地;故节点电压方程为
短路点电压故障分量为:
由此求得短路电流:
当Zf=0时,
考虑到在输电网规划中,关心的是电力系统中有功潮流的分布,而不需要计算各节点的电压幅值,且对于超高压电力网,线路电阻比线路电抗小得多,此时可忽略所有的电阻,则式(7)简化为:
所述步骤B包括如下步骤:
步骤B1:原输电系统网络的节点阻抗矩阵为X=B0 -1,支路k两端的节点为l,m;这里的支路是指两节点间各线路的并联,线路是支路的一个元件;当支路k闭合时,相当于在l,m间增加一条电抗为xk的线路,形成新的网络;新网络的节点电抗矩阵为X′,并且设网络变化前后节点注入功率不变;其节点注入电流相量I和电压相量U,由于注入功率不变,故满足:
U=X′I (9)
原网络的注入电流表示为:
式中:ek为支路k的节点-支路关联矢量,仅在l,m位置对应元素为+1和-1,其余元素为0。
由原网络知:
节点l,m间的电压为:
将式(11)代入式(12),得:
由此得到:
将式(14)代入到式(11)中,得:
将式(15)与(9)比较,得:
从而节点电抗矩阵的变化量为:
将xk换为xlm,ek换为Alm,则式(16)为:
实际追加的支路其xk为正值时追加的是电抗支路,xk为负值时追加的是电容支路。由公式(8)、(18)可知,若要减小节点短路电流,则X′要增大,应当追加电容支路;
步骤B2:由新的节点阻抗矩阵可得短路电流最高节点此时的自阻抗为X′ff,代入公式(8)可得:
选择不同的支路与xlm,即可得到取其最小值I′fmin即为所求最优解。
本发明的有益效果是实际电网经常会由于网架结构不合理、负荷增长而出现短路电流超过最大值的情况,从而影响系统和设备的安全稳定运行。此时电网公司需要在合适的位置新增联络线以降低短路电流。本发明给出了一种可以快速计算得到联络线电阻的方法,可有效降低短路电流。
附图说明
图1为制短路电流的联络线电抗快速计算方法流程图。
图2为单支路等值网络示意图。
图3为14节点系统结构示意图。
具体实施方式
本发明的目的是提出一种抑制短路电流的联络线电抗快速计算方法,下面结合附图和实施例予以说明。
图1所示为制短路电流的联络线电抗快速计算方法流程图。图中所示步骤为:
步骤A.基于系统网络参数形成节点导纳矩阵,对节点导纳矩阵求逆得到描述网络的节点阻抗矩阵,求的系统各节点的短路电流,获取短路水平最高的节点;
步骤B.基于以上计算结果,由该节点引入一条支路,支路电抗为Xk,利用支路追加法推导新网络的节点电抗矩阵,计算此时该节点的短路电流,在能够选支路中获得最优支路(如图2所示)。
以下面实施例详细描述本发明,
图3所示为14节点系统结构示意图。本次算例基于Matlab的14节点系统进行的;为了简化计算过程,在计算过程中有以下计算条件:
1)假设系统有无限大的容量,短路后,系统母线电压维持不变;
2)在计算短路电流时,只考虑发电机、线路等的电抗而忽略其电阻,且不考虑变压器支路的影响即令变压器变比为1,计算过程中的参数均为标幺值;
3)短路电流计算过程中采取的短路形式为三相短路,因为单项短路与两相短路的短路电流均小于三相短路电流。
步骤A1:基于系统的网络参数与拓扑结构形成节点导纳矩阵Y=G+jB,Y是14×14阶的方阵,其中G为电导,B为电纳,令G=0得Y=jB,对Y求导得节点电抗矩阵X=Y-1,取X的对角元即可得各节点的自阻抗Xii=diag(X)。
在对节点导纳矩阵求逆的中,由于导纳阵为奇异阵无法求逆,故去掉系统的平衡节点1,得到13×13阶的节点导纳阵Y,求得的节点自阻抗为:Xii=[0.0496 0.1495 0.09540.0869 0.2485 0.2242 0.4003 0.2336 0.2793 0.3140 0.3788 0.3116 0.3656]T;
步骤A2:代入公式(8)可得节点短路电流:
If=[20.17 6.69 10.49 11.51 4.02 4.46 2.50 4.28 3.58 3.18 2.67 3.212.73]T,
Ifmax=If2=20.17,因此节点2为短路水平最高的节点,针对该节点进行支路追加。
步骤B1:原14节点系统共有20条支路,现以节点2作为支路一端节点新增一条支路,另一端可选节点为节点6、节点7、节点8、节点9、节点10、节点11、节点12、节点13、节点14九种情况。
以节点6为支路末端追加支路,支路电抗Xlm取为0.2。代入公式(18)得即可得新的节点电抗矩阵,取其对角元得以下节点自阻抗:
Xii=[0.0492 0.1488 0.0912 0.0818 0.1431 0.2059 0.3821 0.2039 0.23960.2459 0.2761 0.2197 0.3129]T,
步骤B2:代入公式(8)可得节点短路电流:
If=[20.31 6.72 10.96 12.23 6.98 4.86 2.62 4.90 4.17 4.07 3.62 4.553.19]T;
I′f2=20.31,短路电流增大,应投入电容支路,取Xlm=-0.2。
分别以节点6、节点7、节点8、节点9、节点10、节点11、节点12、节点13、节点14为末端节点,重复步骤B1,B2得到节点2短路电流如下表所示:
表1节点2路电流对比
末端节点 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | 13 | 14 |
短路电流I<sub>f2</sub> | 23.18 | 12.81 | 20.44 | 33.94 | 21.28 | 20.83 | 20.58 | 20.86 | 20.56 |
由上表可知,If2的最小值If2min=12.81,在末端节点为节点7时取得。由以上结果可知,通过支路追加法可以降低节点2的短路电流且最优支路为节点2至节点7支路。
Claims (3)
1.一种抑制短路电流的联络线电抗快速计算方法,其特征在于,包括以下步骤:
步骤A.基于系统网络参数形成节点导纳矩阵,对节点导纳矩阵求逆得到描述网络的节点阻抗矩阵,求得系统各节点的短路电流,获取短路水平最高的节点;
步骤B.基于以上计算结果,由该节点引入一条支路,支路电抗为Xk,利用支路追加法推导新网络的节点阻抗矩阵,计算此时该节点的短路电流,在能够选取的支路中获得最优支路。
2.根据权利要求1所述一种抑制短路电流的联络线电抗快速计算方法,其特征在于,所述步骤A包括:
步骤A1:任一网络用节点阻抗矩阵表示的节点电压方程为U=ZI,I和U分别为:
用矩阵表示如下:
式(1)中,电压相量U为网络各节点对地电压;电流相量I为网络外部向各节点的注入电流;系数矩阵即节点阻抗矩阵Z;
上述矩阵的矩阵对角元素Zii称为自阻抗,表达为:
即当除节点i之外所有其他节点注入电流均为零时,节点i电压与注入电流之比,也就是从节点i看进网络的等效阻抗;
非对角元素Zij称为互阻抗,表示为:
即当除节点i外所有其他节点注入电流均为零时,节点j电压与节点i电流之比,或者说节点i单位电流在节点j引起的电压值;显然Zij是不会为零的,即节点阻抗矩阵总是满阵;
步骤A2:基于步骤A1形成的节点阻抗矩阵,某一点f发生短路故障时,短路点f的接地阻抗为Zf=0;该网络中只有短路点f有注入电流-if,即if由f点流向地;故节点电压方程为:
短路点电压故障分量为:
由此求得短路电流:
当Zf=0时,
考虑到在输电网规划中,关心的是电力系统中有功潮流的分布,而不需要计算各节点的电压幅值,且对于超高压电力网,线路电阻比线路电抗小得多,此时忽略所有的电阻,则式(7)简化为,
3.根据权利要求1所述一种抑制短路电流的联络线电抗快速计算方法,其特征在于,所述步骤B包括如下步骤:
步骤B1:原输电系统网络的节点阻抗矩阵为X=B0 -1,支路k两端的节点为l,m;这里的支路是指两节点间各线路的并联,线路是支路的一个元件;当支路k闭合时,相当于在l,m间增加一条电抗为xk的线路,形成新的网络;新网络的节点电抗矩阵为X′,并且设网络变化前后节点注入功率不变;其节点注入电流相量I和电压相量U,由于注入功率不变,故满足:
U=X′I (9)
原网络的注入电流表示为:
式中:ek为支路k的节点-支路关联矢量,仅在l,m位置对应元素为+1和-1,其余元素为0,
由原网络知:
节点l,m间的电压为:
将式(11)代入式(12),得:
由此得到:
将式(14)代入到式(11)中,得:
将式(15)与(9)比较,得:
从而节点电抗矩阵的变化量为:
将xk换为xlm,ek换为Alm,则式(16)为:
实际支路中xk为正值时是电抗支路,xk为负值时是电容支路,由公式(8)、(18)可知,若要减小节点短路电流,若要增大X′,应当增加电容支路;
步骤B2:由新的节点阻抗矩阵求得短路电流最大的节点,此时的自阻抗为X′ff,代入公式(8)得:
选择不同的支路与xlm,即得到取其最小值I′fmin即为所求最优解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811356521.5A CN109599870A (zh) | 2018-11-14 | 2018-11-14 | 一种抑制短路电流的联络线电抗快速计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811356521.5A CN109599870A (zh) | 2018-11-14 | 2018-11-14 | 一种抑制短路电流的联络线电抗快速计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109599870A true CN109599870A (zh) | 2019-04-09 |
Family
ID=65958160
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811356521.5A Pending CN109599870A (zh) | 2018-11-14 | 2018-11-14 | 一种抑制短路电流的联络线电抗快速计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109599870A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112751338A (zh) * | 2020-12-29 | 2021-05-04 | 广东电网有限责任公司电力调度控制中心 | 一种限制电网短路电流的电抗布置方法及相关装置 |
CN113447803A (zh) * | 2021-06-08 | 2021-09-28 | 国网上海市电力公司 | 断路器开断能力校核的短路电流计算电压系数取值方法 |
CN113904313A (zh) * | 2021-09-30 | 2022-01-07 | 广东电网有限责任公司 | 一种含可控移相器的电力系统的短路电流抑制方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8264238B1 (en) * | 2008-06-13 | 2012-09-11 | Mks Instruments, Inc. | Method for calibrating a broadband voltage/current probe |
CN103219719A (zh) * | 2013-03-21 | 2013-07-24 | 国家电网公司 | 一种利用线路出串运行限制短路电流的分析方法 |
CN104331619A (zh) * | 2014-10-30 | 2015-02-04 | 国网上海市电力公司 | 一种含逆变器接口电源电力系统的短路电流计算方法 |
CN105974271A (zh) * | 2016-07-25 | 2016-09-28 | 山东理工大学 | 一种配电线路故障端口的节点阻抗矩阵计算方法 |
CN108521141A (zh) * | 2018-03-30 | 2018-09-11 | 河海大学 | 一种计及风电场电压分布特性的短路电流计算方法 |
-
2018
- 2018-11-14 CN CN201811356521.5A patent/CN109599870A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8264238B1 (en) * | 2008-06-13 | 2012-09-11 | Mks Instruments, Inc. | Method for calibrating a broadband voltage/current probe |
CN103219719A (zh) * | 2013-03-21 | 2013-07-24 | 国家电网公司 | 一种利用线路出串运行限制短路电流的分析方法 |
CN104331619A (zh) * | 2014-10-30 | 2015-02-04 | 国网上海市电力公司 | 一种含逆变器接口电源电力系统的短路电流计算方法 |
CN105974271A (zh) * | 2016-07-25 | 2016-09-28 | 山东理工大学 | 一种配电线路故障端口的节点阻抗矩阵计算方法 |
CN108521141A (zh) * | 2018-03-30 | 2018-09-11 | 河海大学 | 一种计及风电场电压分布特性的短路电流计算方法 |
Non-Patent Citations (2)
Title |
---|
冯树海 等: ""适用于拓扑调整的短路电流快速计算方法"", 《电力系统保护与控制》 * |
郭杰: ""广州电网短路电流特性及运行对策研究"", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112751338A (zh) * | 2020-12-29 | 2021-05-04 | 广东电网有限责任公司电力调度控制中心 | 一种限制电网短路电流的电抗布置方法及相关装置 |
CN113447803A (zh) * | 2021-06-08 | 2021-09-28 | 国网上海市电力公司 | 断路器开断能力校核的短路电流计算电压系数取值方法 |
CN113447803B (zh) * | 2021-06-08 | 2024-05-28 | 国网上海市电力公司 | 断路器开断能力校核的短路电流计算电压系数取值方法 |
CN113904313A (zh) * | 2021-09-30 | 2022-01-07 | 广东电网有限责任公司 | 一种含可控移相器的电力系统的短路电流抑制方法及装置 |
CN113904313B (zh) * | 2021-09-30 | 2023-06-23 | 广东电网有限责任公司 | 一种含可控移相器的电力系统的短路电流抑制方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Liu et al. | An efficient power flow algorithm for distribution systems with polynomial load | |
CN108923396B (zh) | 一种多端柔性直流电网的短路电流快速计算方法 | |
CN103887810B (zh) | 基于短路比动态跟踪的直流系统持续换相失败判断方法 | |
CN109599870A (zh) | 一种抑制短路电流的联络线电抗快速计算方法 | |
CN107492908A (zh) | 一种基于分布式电源接入的馈线优化规划方法 | |
CN107272645B (zh) | 中性点经电阻接地的光伏电站并网故障模型与分析方法 | |
CN107508322B (zh) | 一种考虑光伏电站并网的电力系统短路故障计算方法 | |
CN108400609A (zh) | 适应于调相机的抑制多馈入直流换相失败选址方法 | |
CN105514970B (zh) | 一种基于图论的upfc监控子站布点算法 | |
CN107292481A (zh) | 一种基于节点重要度的电网关键节点评估方法 | |
CN109066612B (zh) | 基于机电暂态和电磁暂态结合的操作过电压计算方法 | |
CN104899396A (zh) | 一种修正系数矩阵的快速分解法潮流计算方法 | |
CN106786607B (zh) | 一种确定upfc工程安装位置和安装容量的方法 | |
CN108667027A (zh) | 含柔性直流的交直流系统潮流转移搜索与量化分析方法 | |
CN105160148A (zh) | 一种交直流电网连锁故障关键线路辨识方法 | |
CN106849120B (zh) | 多馈入系统中无功补偿装置的布点方法及系统 | |
CN103235198B (zh) | 一种同塔双回输电线路的导线相序排列的确定方法 | |
CN109802418A (zh) | 适用于特高压直流接入的局部电网电磁暂态模型等值方法 | |
CN107800137B (zh) | 一种基于mmc的upfc机电暂态仿真方法 | |
CN107037321B (zh) | 一种小电流接地电力系统的单相接地的稳态计算方法 | |
Halamay et al. | Feasibility study of a high-voltage DC & AC multi-circuit hybrid transmission line | |
CN113270868B (zh) | 用于交流电气化铁路列车供电系统的动态潮流计算方法 | |
CN115714404A (zh) | 治理配电网络电流不平衡的换相软开关系统及评价方法 | |
CN112531695B (zh) | 一种电能质量分支自适应补偿方法及系统 | |
CN109193660B (zh) | 一种多端直流配电网实时仿真装置及其仿真方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
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: 20190409 |