CN102944720A - 考虑动态负荷下的电力系统短路电流衰减计算方法 - Google Patents
考虑动态负荷下的电力系统短路电流衰减计算方法 Download PDFInfo
- Publication number
- CN102944720A CN102944720A CN2012104744492A CN201210474449A CN102944720A CN 102944720 A CN102944720 A CN 102944720A CN 2012104744492 A CN2012104744492 A CN 2012104744492A CN 201210474449 A CN201210474449 A CN 201210474449A CN 102944720 A CN102944720 A CN 102944720A
- Authority
- CN
- China
- Prior art keywords
- short
- generator
- motor
- dynamic load
- 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.)
- Granted
Links
Images
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及一种考虑动态负荷的电力系统短路电流衰减计算方法,根据动态负荷无功需求对动态和静态负荷进行恰当分配,根据发电机和电动机磁链(电势)不突变原理进行发电机和电动机电势的衰减计算,以发电机和电动机短路前瞬间的磁链(电势)作为磁链衰减计算的初始值,通过列写电路方程(线性代数方程组),根据磁链(电势)的衰减求解得到考虑负荷动态特性的短路电流衰减情况。此计算方法能够用于考虑动态负荷的电力系统短路电流衰减特性研究,并且考虑动态负荷后较未考虑动态负荷时所得计算结果更符合电网实际短路电流情况。
Description
技术领域
本发明涉及一种短路电流计算方法,特别涉及一种考虑动态负荷下的电力系统短路电流衰减计算方法。
背景技术
短路是电力系统的严重故障。电力系统短路电流中包含短路电流周期分量和非周期分量。一般意义上的短路电流计算(短路电流实用计算)主要是指短路电流周期分量计算,本说明的短路电流计算也指短路电流周期分量的计算。
根据目前的短路电流计算标准,对于输电网的短路电流计算,可以不考虑负荷的动态特性,即将负荷等效成恒阻抗、恒功率或恒电流的静态模型,国内电力行业采用的常规短路电流计算在基于潮流计算时即将负荷等效为上述静态模型。
然而,对于负荷的不同处理方式,可能会导致不同的短路电流计算结论和短路电流超标判断。特别地,对于大负荷中心的输电网进行短路电流计算时,当故障切除时间较短时,动态负荷对三相短路电流周期分量的影响较大,值得考虑。有研究表明,如果在仿真计算时忽略负荷的动态特性,一旦发生短路故障,现场的故障录波值比仿真计算值明显偏大,这样就给电力系统的安全运行带来隐患。
同时,国内电力行业的常规短路电流计算只能计算短路瞬间短路电流的周期分量,不能计算短路电流周期分量的衰减情况,事实上,短路电流在断路器跳闸前是存在一定程度衰减的。因此,断路器开断时短路电流是否已衰减到常规短路电流计算所得结果,需要在考虑负荷动态特性的同时计及短路电流衰减进行计算。
发明内容
本发明是针对电力系统常规短路电流计算不考虑动态负荷、不考虑短路电流的衰减,给电力系统带来很大的安全隐患的问题,提出了一种考虑动态负荷下的电力系统短路电流衰减计算方法,该方法基于发电机和电动机磁链(电势)不突变原理,根据动态负荷无功需求对动态和静态负荷进行恰当分配,通过列写电路方程(线性代数方程组),根据电势的衰减求解得到考虑负荷动态特性的短路电流衰减情况,使计算结果更加接近实际情况。
本发明的技术方案为:一种考虑动态负荷下的电力系统短路电流衰减计算方法,具体包括如下步骤:
1)在已知潮流分布的情况下对动态和静态负荷进行分配,假定母线的负荷功率是 ,有功部分P L ,按照静态有功负荷与动态有功负荷的比例分为和P,即;无功部分,先根据母线电压V及动态有功负荷P确定动态负荷需要的无功Q,即Q=f(P,V),然后总的无功负荷减去动态负荷需要的无功Q后剩余部分即为静态负荷的无功,即;
2)动态负荷无功需求计算,设定电力系统中的动态负荷电动机作为动态负荷模型,根据已知负荷节点的端电压V,分配给动态负荷的有功功率P,设定电动机等值电路计算动态负荷无功功率Q,动态负荷为P+jQ;
4)对预设故障点,根据磁链(电势)不突变原理进行发电机和电动机电势的衰减计算,即当t=0时,电动机短路后暂态电势的衰减初始值等于电动机短路前瞬间暂态电势,发电机短路后d、q轴次暂态电势等于发电机短路前瞬间d、q轴次暂态电势;5)编程解网络方程(即节点导纳矩阵),根据发电机和电动机电势的衰减得到短路后各时刻短路电流周期分量的衰减情况。
所述步骤4)中对于单鼠笼感应电动机,短路后暂态电势的衰减计算公式为:
发电机短路后d轴次暂态电势的衰减计算公式:
式中,为发电机短路前瞬间q轴次暂态电势;为发电机短路前瞬间定子绕组d轴电流;分别为发电机的稳态电抗、暂态电抗和次暂态电抗;为基于节点导纳矩阵计算的短路后瞬间发电机机端电压与机端电流的商;分别为暂态分量和次暂态分量衰减的时间常数。
所述步骤5)当短路点与发电机或电动机直接相连时,
电动机反馈的短路电流周期分量的衰减公式如下:
发电机提供的短路电流周期分量的衰减公式如下:
其中,,为发电机短路前瞬间q轴次暂态电势;为发电机短路前瞬间定子绕组d轴电流,分别为发电机的稳态电抗、暂态电抗和次暂态电抗,为发电机机端到短路点的电抗,与基于节点导纳矩阵计算的短路后瞬间发电机机端电压与机端电流的商一致,分别为暂态分量和次暂态分量衰减的时间常数。
本发明的有益效果在于:本发明为考虑动态负荷特性的电力系统短路电流衰减计算方法,该方法考虑电力系统中实际存在的动态负荷,将动态负荷等效为感应电动机,根据电势的衰减求解得到发电机和电动机提供的短路电流衰减情况,此计算方法能够用于考虑动态负荷的电网短路电流衰减特性研究,并且考虑动态负荷后较未考虑动态负荷时所得计算结果更符合电网实际短路电流情况。
附图说明
图1为负荷分配示意图;
图2为感应电动机T形等值电路图;
图3为本发明考虑动态负荷下的电力系统短路电流衰减计算方法实施例A厂500kV母线三相短路时短路电流周期分量衰减对比示意图;
图4为本发明考虑动态负荷下的电力系统短路电流衰减计算方法实施例B站500kV母线三相短路时短路电流周期分量衰减对比示意图;
图5为本发明考虑动态负荷下的电力系统短路电流衰减计算方法实施例C站500kV母线三相短路时短路电流周期分量衰减对比示意图;
图6为本发明考虑动态负荷下的电力系统短路电流衰减计算方法实施例C站220kV母线三相短路时短路电流周期分量衰减对比示意图;
图7为本发明考虑动态负荷下的电力系统短路电流衰减计算方法实施例D站220kV母线1-3-4段三相短路时短路电流周期分量衰减对比示意图;
图8为本发明考虑动态负荷下的电力系统短路电流衰减计算方法实施例B站220kV母线三相短路时短路电流周期分量衰减对比示意图。
具体实施方式
考虑动态负荷下的电力系统短路电流衰减计算方法根据磁链(电势)不突变原理,以发电机和电动机短路前瞬间的磁链(电势)作为衰减计算的初值,通过列写电路方程(线性代数方程组),用计算机求解得到考虑负荷动态特性的短路电流衰减情况,即考虑动态负荷的电力系统短路电流衰减实用计算方法。
考虑动态负荷下的电力系统短路电流衰减计算方法的具体实施步骤如下:
(1)在已知潮流分布的情况下对动态和静态负荷进行分配。假定某母线的负荷功率是,欲将该负荷分出一部分作为动态负荷,如图1所示负荷分配示意图。对于有功部分,直接按照静态与动态负荷的比例分为(静态负荷)和P(动态负荷),;而对于无功部分,先根据母线电压V及动态有功负荷P确定动态负荷需要的无功Q,即Q=f(P,V),然后总的无功负荷减去动态负荷需要的无功Q后剩余部分即为静态负荷的无功,即,这样考虑动态负荷后的系统潮流与原潮流一致。
(2)动态负荷无功需求计算过程。电力系统中的动态负荷包含同步电动机和异步电动机,其中以异步电动机居多,单鼠笼感应电动机是最常用的一种模型。下面以单鼠笼感应电动机作为动态负荷模型,计算其无功需求。
已知负荷节点的端电压V,分配给动态负荷的有功功率P,以及感应电动机T形等值电路(如图2所示)的相关参数,其中分别为考虑配电网经参数修正后的电动机定子绕组电阻、漏抗,为定子、转子绕组间的互感抗,分别为转子绕组电阻、漏抗。
感应电动机需要的无功功率Q的计算过程如式(1)至式(5)所示。
式中:
(2)
(4)对预设故障点,根据电势不突变原理进行发电机和电动机电势的衰减初始值计算。对于发电机,以d轴为例,短路后d轴次暂态电势的衰减计算如式(7)所示。
式中,为发电机短路前瞬间q轴次暂态电势;为发电机短路前瞬间定子绕组d轴电流;分别为发电机的稳态电抗、暂态电抗和次暂态电抗;为基于节点导纳矩阵计算的短路后瞬间发电机机端电压与机端电流的商;分别为暂态分量和次暂态分量衰减的时间常数。
对于单鼠笼感应电动机,短路后暂态电势的衰减计算如式(8)所示。
(5)编程解网络方程(即节点导纳矩阵),根据发电机和电动机电势的衰减得到短路后各时刻短路电流周期分量的衰减情况。
作为特例,当短路点与发电机或电动机直接相连时,发电机提供的短路电流周期分量的衰减如式(9)所示,电动机反馈的短路电流周期分量的衰减如(10)所示,这将与求解网络方程所得结果一致。
计算短路电流的衰减可用上述解节点电压方程(为线性代数方程)的方法。该方法所得结果可与一般认为更加精确(同时也更复杂、计算量更大)的机电暂态仿真计算结果相比较判断其正确性。机电暂态仿真计算方法,通过发电机的微分方程与网络代数方程联立求解得到暂态过程中的电压、电流、角度、频率等变量。但由于机电暂态计算过程复杂,运算量大,一般只用来进行电力系统稳定性计算,而不在电力企业日常进行的短路电流计算中使用。
例1这里选择我国某大都市电网实际数据作为本发明的应用举例,应用本发明提出的实用计算新方法对该电网典型500kV和220kV母线三相短路时进行考虑负荷动态特性的短路电流周期分量的衰减计算(初始状态基于潮流),并与机电暂态的计算结果进行对比,如表1所示A厂500kV母线三相短路时短路电流周期分量衰减对比表和图3至图8所示短路电流周期分量衰减对比图。
表1
短路时间/s | 实用计算新方法/kA | 机电暂态/kA | 相差率/% |
0.00 | 63.52 | 63.53 | -0.02 |
0.01 | 60.06 | 61.17 | -1.81 |
0.02 | 57.23 | 57.95 | -1.24 |
0.03 | 54.90 | 56.15 | -2.23 |
0.04 | 52.95 | 54.44 | -2.74 |
0.05 | 51.32 | 52.96 | -3.10 |
0.06 | 49.92 | 51.63 | -3.31 |
0.07 | 48.73 | 50.44 | -3.39 |
0.08 | 47.70 | 49.36 | -3.36 |
0.09 | 46.80 | 48.36 | -3.23 |
0.10 | 46.01 | 47.44 | -3.01 |
从以上厂站的短路电流周期分量衰减计算来看,实用计算新方法得到的结果和机电暂态的计算结果比较接近。在研究的时间范围0-0.1s内,相差率在4%以内,其中短路电流的初值(0s时刻)两种计算方法一致,最大负偏差为-3.71%(出现在C站500kV母线短路后0.07s),最大正偏差为1.65%(出现在D站220kV母线1-3-4段短路后0.1s)。
例2为了进一步检验本发明提出的考虑负荷动态特性的短路电流衰减实用计算方法,现将举例1表1中A厂500kV母线三相短路时短路电流周期分量的衰减计算结果与该厂该母线实际发生过的一次短路事故的现场录波数据进行对比分析。该事故为三相金属性接地短路故障,保护在40ms后动作切除故障。
由于仿真计算和录波数据基于该电网在一年内的不同的运行方式,所以两种情况下短路电流数据的绝对值没有可比性,但其衰减率(即衰减变化趋势)具有一定的可比性。现将表1中用本发明方法计算所得的短路电流衰减情况与现场录波衰减数据做一对比,同时列出计算时未考虑动态负荷(即负荷全部等效为静态负荷)的情况,如表2所示A厂500kV母线三相短路时短路电流周期分量衰减率对比表。
表2
短路时间/s | 动态负荷/% | 静态负荷/% | 录波数据/% |
0.01 | 5.45 | 4.49 | 4.85 |
0.02 | 9.90 | 7.90 | 13.56 |
0.03 | 13.57 | 10.53 | 14.57 |
0.04 | 16.64 | 12.58 | 19.42 |
从表2可以看出,从短路后0.02s开始,本发明方法所得的考虑动态负荷后的短路电流周期分量衰减率比未考虑动态负荷的情况更接近实际故障录波分解得到的短路电流周期分量衰减情况,即本发明所述方法所得的短路电流周期分量的衰减速度和实际故障录波数据所揭示的情况更接近。
通过应用举例1和2的计算分析,可以说明本发明提出的考虑负荷动态特性的短路电流衰减实用计算方法能够用于考虑动态负荷的电网短路电流衰减特性研究,并且考虑动态负荷后较未考虑动态负荷时所得计算结果更符合电网实际短路电流情况。
Claims (3)
1.一种考虑动态负荷下的电力系统短路电流衰减计算方法,其特征在于,具体包括如下步骤:
1)在已知潮流分布的情况下对动态和静态负荷进行分配,假定母线的负荷功率是 ,有功部分P L ,按照静态有功负荷与动态有功负荷的比例分为和P,即;无功部分,先根据母线电压V及动态有功负荷P确定动态负荷需要的无功Q,即Q=f(P,V),然后总的无功负荷减去动态负荷需要的无功Q后剩余部分即为静态负荷的无功,即;
2)动态负荷无功需求计算,设定电力系统中的动态负荷电动机作为动态负荷模型,根据已知负荷节点的端电压V,分配给动态负荷的有功功率P,设定电动机等值电路计算动态负荷无功功率Q,动态负荷为P+jQ;
4)对预设故障点,根据磁链(电势)不突变原理进行发电机和电动机磁链(电势)的衰减初始值计算,即当t=0时,电动机短路后暂态电势等于电动机短路前瞬间暂态电势,发电机短路后d、q轴次暂态电势分别等于发电机短路前瞬间d、q轴次暂态电势;
5)编程解网络方程(即节点导纳矩阵),根据发电机和电动机电势的衰减得到短路后各时刻短路电流周期分量的衰减情况。
3.根据权利要求1所述考虑动态负荷下的电力系统短路电流衰减计算方法,其特征在于,所述步骤5)当短路点与发电机或电动机直接相连时,
电动机反馈的短路电流周期分量的衰减公式如下:
发电机提供的短路电流周期分量的衰减公式如下:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210474449.2A CN102944720B (zh) | 2012-11-21 | 2012-11-21 | 考虑动态负荷下的电力系统短路电流衰减计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210474449.2A CN102944720B (zh) | 2012-11-21 | 2012-11-21 | 考虑动态负荷下的电力系统短路电流衰减计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102944720A true CN102944720A (zh) | 2013-02-27 |
CN102944720B CN102944720B (zh) | 2015-02-25 |
Family
ID=47727680
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210474449.2A Expired - Fee Related CN102944720B (zh) | 2012-11-21 | 2012-11-21 | 考虑动态负荷下的电力系统短路电流衰减计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102944720B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103698586A (zh) * | 2014-01-15 | 2014-04-02 | 重庆大学 | 一种确定含双馈感应发电机三相短路电流的磁链解析方法 |
WO2015035892A1 (zh) * | 2013-09-10 | 2015-03-19 | 国家电网公司 | 一种双馈异步发电机组暂态短路电流的计算方法 |
CN106410783A (zh) * | 2016-04-14 | 2017-02-15 | 上海电力学院 | 一种基于坐标变换的电力系统短路电流衰减计算方法 |
CN106990326A (zh) * | 2017-05-10 | 2017-07-28 | 上海电力学院 | 电力系统短路电流直流分量计算方法 |
CN113595070A (zh) * | 2021-07-28 | 2021-11-02 | 国网湖南省电力有限公司 | 基于电网正常运行时潮流计算的短路转移阻抗计算方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101957398A (zh) * | 2010-09-16 | 2011-01-26 | 河北省电力研究院 | 一种基于机电与电磁暂态混合仿真技术检测计算电网一次时间常数的方法 |
CN102646981A (zh) * | 2012-04-19 | 2012-08-22 | 上海电力学院 | 电力系统短路电流衰减实用计算方法 |
-
2012
- 2012-11-21 CN CN201210474449.2A patent/CN102944720B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101957398A (zh) * | 2010-09-16 | 2011-01-26 | 河北省电力研究院 | 一种基于机电与电磁暂态混合仿真技术检测计算电网一次时间常数的方法 |
CN102646981A (zh) * | 2012-04-19 | 2012-08-22 | 上海电力学院 | 电力系统短路电流衰减实用计算方法 |
Non-Patent Citations (2)
Title |
---|
尹建华: "考虑感应电动机负荷电磁暂态特性的系统短路电流计算方法", 《南方电网技术》 * |
曹炜等: "基于PSS/E考虑动态负荷的短路电流计算方法", 《上海电力学院学报》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015035892A1 (zh) * | 2013-09-10 | 2015-03-19 | 国家电网公司 | 一种双馈异步发电机组暂态短路电流的计算方法 |
CN103698586A (zh) * | 2014-01-15 | 2014-04-02 | 重庆大学 | 一种确定含双馈感应发电机三相短路电流的磁链解析方法 |
CN103698586B (zh) * | 2014-01-15 | 2016-09-07 | 重庆大学 | 一种确定含双馈感应发电机三相短路电流的磁链解析方法 |
CN106410783A (zh) * | 2016-04-14 | 2017-02-15 | 上海电力学院 | 一种基于坐标变换的电力系统短路电流衰减计算方法 |
CN106410783B (zh) * | 2016-04-14 | 2019-02-22 | 上海电力学院 | 一种基于坐标变换的电力系统短路电流衰减计算方法 |
CN106990326A (zh) * | 2017-05-10 | 2017-07-28 | 上海电力学院 | 电力系统短路电流直流分量计算方法 |
CN106990326B (zh) * | 2017-05-10 | 2019-05-24 | 上海电力学院 | 电力系统短路电流直流分量计算方法 |
CN113595070A (zh) * | 2021-07-28 | 2021-11-02 | 国网湖南省电力有限公司 | 基于电网正常运行时潮流计算的短路转移阻抗计算方法 |
CN113595070B (zh) * | 2021-07-28 | 2023-05-16 | 国网湖南省电力有限公司 | 基于电网正常运行时潮流计算的短路转移阻抗计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102944720B (zh) | 2015-02-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Lindenmeyer et al. | An induction motor parameter estimation method | |
Zhang et al. | An improved synchronverter model and its dynamic behaviour comparison with synchronous generator | |
CN102944720B (zh) | 考虑动态负荷下的电力系统短路电流衰减计算方法 | |
Kamel | Effect of wind generation system types on Micro-Grid (MG) fault performance during both standalone and grid connected modes | |
WO2019007354A1 (zh) | 一种考虑不同转差的双馈风力发电系统短路电流计算方法 | |
CN104699996B (zh) | 基于双馈风电机组等值的电力系统非对称短路工频电气量计算方法 | |
CN106990326A (zh) | 电力系统短路电流直流分量计算方法 | |
CN111460750B (zh) | 一种发电机、变压器保护整定值自动生成系统和方法 | |
Jabr et al. | Effects of main and leakage flux saturation on the transient performances of doubly-fed wind driven induction generator | |
Abdel‐Akher et al. | Unbalanced distribution power‐flow model and analysis of wind turbine generating systems | |
Zhang et al. | Improved schemes for traditional current‐based phase selectors in wind power systems | |
Kamel et al. | Effect of load type on standalone micro grid fault performance | |
Suvorov et al. | A validation approach for short‐circuit currents calculation in large‐scale power systems | |
CN115270597A (zh) | 配电网继电保护适应性仿真测试方法、仿真模型和装置 | |
Shao et al. | Voltage stability analysis of multi-infeed HVDC systems using small-signal stability assessment | |
Ouyang et al. | Short-circuit current of doubly fed induction generator under partial and asymmetrical voltage drop | |
Qi et al. | Generic EMT modeling method of Type-4 wind turbine generators based on detailed FRT studies | |
Aboul-Seoud et al. | Dynamic modeling of induction motor loads for transient voltage stability studies | |
CN102608498B (zh) | 新能源输电线路故障选相方法 | |
Gong et al. | Research on influence and resolution of the relay protections with electric vehicle charging station integrating into distribution network | |
CN113437732B (zh) | 一种光伏发电并网联络线纵联保护方法及系统 | |
Liu et al. | A FLEXIBLE REAL TIME NETWORK MODEL FOR EVALUATING HVDC SYSTEMS'IMPACT ON AC PROTECTION PERFORMANCE | |
Mora et al. | The effect of induction generators on the transient stability of a laboratory electric power system | |
Zhang et al. | Equivalent modeling method of induction motor contribution to short-circuit current | |
CN104899677B (zh) | 扰动恢复中500kV辐射状供电网络负荷最大无功需求快估方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150225 Termination date: 20201121 |