CN103280783A - 基于临界机组对的故障临界切除时间计算方法 - Google Patents
基于临界机组对的故障临界切除时间计算方法 Download PDFInfo
- Publication number
- CN103280783A CN103280783A CN2013101692734A CN201310169273A CN103280783A CN 103280783 A CN103280783 A CN 103280783A CN 2013101692734 A CN2013101692734 A CN 2013101692734A CN 201310169273 A CN201310169273 A CN 201310169273A CN 103280783 A CN103280783 A CN 103280783A
- Authority
- CN
- China
- Prior art keywords
- critical
- unit
- delta
- fault
- equ
- 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
技术领域
本发明涉及电力系统稳定与控制技术领域,特别涉及一种基于临界机组对的故障下临界切除时间计算方法。
背景技术
临界切除时间(Critical Clearing Time,CCT)是指电力系统在发生故障的情况下,能够保持系统稳定的最晚切除时间。故障临界切除时间是表征系统暂态稳定性的常用指标,CCT越大,则表明系统抵御外界冲击的能力越强,它为运行人员提供了重要的操作依据。因此暂态稳定研究中的一个重要的内容就是如何快速而准确地计算CCT。由于大电网下运行方式的增多以及预想事故集的增大,对于CCT的计算要耗费大量的时间,因此提高单个CCT的计算速度,从而缩短整个预想事故集的CCT计算时间尤为重要。在本领域中,通常采用时域仿真方法、扩展等面积法,以及暂态能量函数方法等方法来计算CCT的值。但是,时域仿真方法由于计算时间过长而无法适应在线分析的要求;扩展等面积法虽能提高计算速度,但依赖于机群分群的正确性,而对于临界机群的鉴别,目前仍没有很好的被解决;暂态能量函数方法虽较时域仿真方法的计算速度快很多,但该方法本身对系统信息量的需求是极大的,在迭代计算过程中,需要获取所有发电机的暂态过程数据。
本发明借助临界机组对信息来计算故障后系统的CCT,避免了对于临界机群的鉴别问题,从而在在获得精确计算结果的同时,缩短了计算时间,并减轻了电网运行人员的工作量。
发明内容
本发明提供了一种基于临界机组对的故障临界切除时间计算方法,其包含以下步骤:一种基于临界机组对的故障临界切除时间计算方法,包括:
1)从广域相量测量系统采集电力系统运行中的各个参量数据,形成数据文件;
2)针对调度部门提供的预想事故集中的每一故障,基于所述数据文件以最大切除时间进行电力系统暂态数值仿真,生成仿真数据;
3)根据仿真数据选取临界机组对;
4)判断临界机组对的稳定性,并最终判断所述电力系统的稳定性;
5)根据稳定性的判断结果将故障分为:失稳故障,不失稳故障;
6)如果所述故障为失稳故障,则计算临界机组对的临界切除时间;
7)计算所述电力系统的临界切除时间。
本发明的有益效果如下:在大型电力系统中,针对故障集中每一故障计算系统临界切除时间时,采用临界机组对进行计算,避开了识别发电机临界机群这一难题;且无论包含多少台发电机的电力系统,只需要很少量的几台发电机信息,减少了繁重的数据处理过程;不需要机组凝聚等值,因此计算公式简单,计算速度快;尤其对于超大型的电力系统而言,这种效果将更加明显。
附图说明
图1为临界机组对暂态过程功角曲线示意图
图2为10机39节点系统网络图
具体实施方式
下面结合实例来说明本发明的具体实施步骤。本实施在新英格兰10机39节点系统中进行,该系统的接线图如图2所示。
步骤1:从电力调度部门获取当前电网运行状态下的预想事故集,并从预想事故集中提取一个故障。假设所获取的故障为节点3处发生三相接地短路,tclmax=0.27秒切除故障。
步骤:2:选取临界机组对
选取临界机组对的方法为:
若记Ωcr为临界机群集合,Ωnon-cr为剩余机群集合,则对于任意的i∈Ωcr,j∈Ωnon-cr,称机组对(i,j)为“临界机组对”。临界机组对选取步骤如下:
(2)选择转子功角δ最大的n1个机组及转子角速度Δω最大的n2个机组作为临界机的候选集合Ωcr,转子功角δ最小的m1个机组及转子角速度Δω最小的m2个机组作为非临界机的候选集合Ωnon-cr,由此形成临界机组对集合Ω的初始判别:
Ω={(i,j)|i∈Ωcr,j∈Ωnon-cr}
其中,n1、n2、m1、m2的取值遵循如下规律:当系统中机组数目为个位数时,n1、n2可都取为2,m1、m2取为1;当系统中机组数目达十位数时,n1、n2可都取为2-3,m1、m2取为1-2;当系统中机组数目达百位数时,n1、n2可都取为3-4,m1、m2取为2-3。
(3)检测Ωcr和Ωnon-cr是否有交集?若有,则设系统中有N台发电机,将N台发电机在此时刻的值由大到小排序,依次编号为1、2、……、N。若对应的编号i满足i<N/2(N/2若为小数,则取整),将Gi归属集合Ωcr。反之,若i>N/2,则将Gi归属集合Ωnon-cr。
(4)检测Ωcr和Ωnon-cr是否有空集?若有,则对剩余机群,重复步骤(2)和(3),直至Ωcr和Ωnon-cr都不为空集。
(5)将集合Ωcr和Ωnon-cr中的元素按照排列组合组成临界机组对,若Ωcr中含有n个元素,Ωnon-cr中含有m个元素,则临界机组对为n×m对,即Ω中有n×m个元素。
步骤1中的故障,3个工频周期后即0.32秒时,发电机组的各个变量信息如表1所示
表10.32s时刻发电机功角、角速度、角加速度值
因此初始临界机集合为Ωcr={34,38},非临界机集合为Ωnon-cr={39}。两集合没有交集,因此最终确定的临界机组对集合为Ω={34-39,38-39}。
步骤3:判断临界机组对的稳定性,系统稳定性
判断临界机组对的稳定性的方法为:
图1所示为任意机组对(i,j)故障前、故障中及故障后机组对功角曲线示意图,机组对(i,j)故障切除后加速面积为:
式中,
Mi、Mj——机组i、j的转子惯量时间常数(s);
Pmi、Pmj——机组i、j的机械功率(p.u.);
δij0——故障前机组对稳态工作时对应的相位角(rad);
δijc——故障切除时刻机组对相位角(rad)。
机组对(i,j)故障切除后减速面积为为:
式中,
——机组i、j的故障后电磁功率(p.u.);
δiju——故障后不稳态平衡点处机组对的相位角。
机组对稳定判据为:若A<C,机组对(i,j)是稳定的;若A>C,机组对(i,j)是失稳的;若A=C,机组对(i,j)临界稳定。
系统稳定性的判别方法为:
实际上,系统中任意“临界机组对”的稳定性是能够代替整个系统稳定性的问题的,即(1)若临界机组对(i,j)失稳,根据稳定性的定义可知,系统一定是失稳的;(2)若临界机组对(i,j)稳定,系统一定是稳定的。
因此任一临界机组对的稳定性即可代表全系统的稳定性。只要正确的选择了一对临界机组对,即可判断全系统的稳定性。为了防止临界机组对选择错误而使结果误判,在前面的选择临界机组对时,选择了几组,因此全系统的稳定性判断方法应为:若几组临界机组对都稳定,则系统是稳定的;若有一组临界机组对失稳,则系统是失稳的。
因此针对步骤1中的故障,利用公式(a)(b)计算第一个临界机组对34-39的加速面积A,减速面积C,比较A与C的大小,判断结果为临界机组对34-39失稳,则系统为失稳。已不需要再计算临界机组对38-39的稳定性。
步骤4:系统临界切除时间计算
按照公式
计算实际临界切除时间tcr,并用下式校验
|tcr-tcl|<ε(d)
若满足公式(d),则计算结束。若不满足,则以这个tcr为新的切除时间tcl,再次仿真获得机组对实际故障轨迹,若失稳,仍采用公式(c)计算新的临界切除时间;若稳定,采用公式下面(e)计算新的临界切除时间。
重复以上求解过程,直至相邻两次tcr小于允许误差值ε为止。
结合图1及公式(c)(d)(e),计算临界机组对34-39临界切除时间过程如表2所示
表2临界切除时间迭代过程
则系统在该故障下的临界切除时间为0.22秒。
尽管上面已经通过对本发明的具体实施例的描述对本发明进行了披露,但是,应该理解,上述的所有实施例和示例均是示例性的,而非限制性的。本领域的技术人员可在所附权利要求的精神和范围内设计对本发明的各种修改、改进或者等同物。这些修改、改进或者等同物也应当被认为包括在本发明的保护范围内。
Claims (6)
1.一种基于临界机组对的故障临界切除时间计算方法,其特征在于包括以下步骤:
1)从广域相量测量系统采集电力系统运行中的各个参量数据,形成数据文件;
2)针对调度部门提供的预想事故集中的每一故障,基于所述数据文件以最大切除时间进行电力系统暂态数值仿真,生成仿真数据;
3)根据仿真数据选取临界机组对;
4)判断临界机组对的稳定性,并最终判断所述电力系统的稳定性;
5)根据稳定性的判断结果将故障分为:失稳故障,不失稳故障;
6)如果所述故障为失稳故障,则计算临界机组对的临界切除时间;
7)计算所述电力系统的临界切除时间。
2.如权利要求1所述的基于临界机组对的故障临界切除时间计算方法,其特征在于:
步骤3)中所述的选取临界机组对进一步包括:
(2)选择故障切除后某时刻t1转子功角δ最大的n1个机组及转子角速度Δω最大的n2个机组作为临界机的候选集合Ωcr,转子功角δ最小的m1个机组及转子角速度Δω最小的m2个机组作为非临界机的候选集合Ωnon-cr,由此形成临界机组对集合Ω的初始判别:
Ω={(i,j)|i∈Ωcr,j∈Ωnon-cr}
其中,Ωcr为临界机群集合,Ωnon-cr为剩余机群集合,则对于任意的i∈Ωcr,j∈Ωnon-cr,机组对(i,j)为临界机组对;
(3)检测Ωcr和Ωnon-cr是否有交集,若有,则设系统中有N台发电机,将N台发电机在t1时刻的值由大到小排序,依次编号为1、2、……、N;若对应的编号i满足i<N/2(N/2若为小数,则取整),将Gi归属集合Ωcr;反之,若i>N/2,则将Gi归属集合Ωnon-cr,其中,Gi代表第i台发电机;
(4)检测Ωcr和Ωnon-cr是否有空集,若有,则对剩余机群,重复步骤(2)和(3),直至Ωcr和Ωnon-cr都不为空集;
(5)将集合Ωcr和Ωnon-cr中的元素按照排列组合组成临界机组对,若Ωcr中含有n个元素,Ωnon-cr中含有m个元素,则临界机组对为n×m对,即Ω中有n×m个元素。
3.如权利要求1所述的基于临界机组对的故障临界切除时间计算方法,其特征在于:
步骤4)中所述的判断临界机组对的稳定性的方法进一步包括:
(1)临界机组对(i,j)故障切除后加速面积通过以下公式计算:
其中,
A—加速面积;
Mi、Mj——机组i、j的转子惯量时间常数;
Pmi、Pmj——机组i、j的机械功率;
δij0——故障前机组对稳态工作时对应的相位角;
δijc——故障切除时刻机组对相位角;
(2)临界机组对(i,j)故障切除后减速面积为:
其中,
C—减速面积;
δiju——故障后不稳态平衡点处机组对的相位角;
(3)其中,若A<C,机组对(i,j)是稳定的;若A>C,机组对(i,j)是失稳的;若A=C,机组对(i,j)临界稳定。
4.如权利要求1所述的基于临界机组对的故障临界切除时间计算方法,其特征在于:
步骤4)中所述的判断系统稳定性的步骤进一步包括:
若临界机组对(i,j)失稳,则系统为失稳的;若临界机组对(i,j)稳定,则系为稳定。
5.如权利要求1所述的基于临界机组对的故障临界切除时间计算方法,其特征在于:
步骤6)中所述的计算临界机组对的临界切除时间进一步包括:
(1)通过以下公式计算临界切除时间tcr
式中,δijc1——前一次仿真故障切除时间对应的功角;
公式(a)中求得的是实际临界功角δijcr,通过暂态数值仿真数据即可确定与δijcr对应的实际临界切除时间tcr值;
(2)使用以下公式校验:
|tcr-tcl|<ε(b)
式中,tcr—本次求得的临界切除时间;
tcl—前次仿真求得的临界切除时间;
ε—误差精度;
(3)若上述校验公式(b)成立,则计算结束;若不成立,则以tcr代替切除时间tcl,再次仿真获得机组对实际故障轨迹,若失稳,仍采用公式(a)计算新的临界切除时间;若稳定,则采用以下公式计算新的临界切除时间。
直至相邻两次tcr小于允许误差值ε为止。
6.如权利要求1所述的基于临界机组对的故障临界切除时间计算方法,其特征在于:
步骤7)中所述的计算全系统的临界切除时间的步骤进一步包括:
使用临界机组对的临界切除时间表示全系统的临界切除时间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310169273.4A CN103280783B (zh) | 2013-05-09 | 2013-05-09 | 基于临界机组对的故障临界切除时间计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310169273.4A CN103280783B (zh) | 2013-05-09 | 2013-05-09 | 基于临界机组对的故障临界切除时间计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103280783A true CN103280783A (zh) | 2013-09-04 |
CN103280783B CN103280783B (zh) | 2015-07-22 |
Family
ID=49063259
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310169273.4A Expired - Fee Related CN103280783B (zh) | 2013-05-09 | 2013-05-09 | 基于临界机组对的故障临界切除时间计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103280783B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107069712A (zh) * | 2017-04-01 | 2017-08-18 | 华北电力大学 | 基于李雅普诺夫的电力系统临界切除时间计算系统及方法 |
CN108199418A (zh) * | 2017-12-29 | 2018-06-22 | 湖南大学 | 一种互联电网暂态稳定紧急控制切机方案制定方法 |
CN110333444A (zh) * | 2019-07-11 | 2019-10-15 | 中国石油大学(华东) | 一种基于电动机铭牌参数的暂降临界切除时间计算方法 |
CN111969924A (zh) * | 2020-08-18 | 2020-11-20 | 四川大学 | 一种交流接触器自适应动作控制方法 |
CN113381401A (zh) * | 2021-05-31 | 2021-09-10 | 国网宁夏电力有限公司电力科学研究院 | 基于归一化能量并行计算临界切除时间的方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101751518A (zh) * | 2009-12-21 | 2010-06-23 | 江西省电力公司 | 基于集群的故障临界切除时间并行计算方法 |
CN102570422A (zh) * | 2012-01-05 | 2012-07-11 | 东北电力大学 | 根据实测轨迹计算电力系统极限切除时间的方法 |
CN102855382A (zh) * | 2012-07-09 | 2013-01-02 | 中国电力科学研究院 | 一种电力系统三相短路故障临界切除时间的在线求取方法 |
-
2013
- 2013-05-09 CN CN201310169273.4A patent/CN103280783B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101751518A (zh) * | 2009-12-21 | 2010-06-23 | 江西省电力公司 | 基于集群的故障临界切除时间并行计算方法 |
CN102570422A (zh) * | 2012-01-05 | 2012-07-11 | 东北电力大学 | 根据实测轨迹计算电力系统极限切除时间的方法 |
CN102855382A (zh) * | 2012-07-09 | 2013-01-02 | 中国电力科学研究院 | 一种电力系统三相短路故障临界切除时间的在线求取方法 |
Non-Patent Citations (2)
Title |
---|
LU FANG,YU JI-LAI: "Transient stability analysis with equal area criterion directly used to a non-equivalent generator pair", 《INTERNATIONAL CONFERENCE ON POWER ENGINEERING, ENERGY AND ELECTRICAL DRIVES, 2009. POWERENG 2009》, 20 March 2009 (2009-03-20), pages 386 - 389, XP031453878 * |
王守相等: "事故临界切除时区计算及其在事故扫描中的应用", 《电网技术》, vol. 27, no. 10, 31 October 2003 (2003-10-31), pages 72 - 77 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107069712A (zh) * | 2017-04-01 | 2017-08-18 | 华北电力大学 | 基于李雅普诺夫的电力系统临界切除时间计算系统及方法 |
CN108199418A (zh) * | 2017-12-29 | 2018-06-22 | 湖南大学 | 一种互联电网暂态稳定紧急控制切机方案制定方法 |
CN108199418B (zh) * | 2017-12-29 | 2021-03-19 | 湖南大学 | 一种互联电网暂态稳定紧急控制切机方案制定方法 |
CN110333444A (zh) * | 2019-07-11 | 2019-10-15 | 中国石油大学(华东) | 一种基于电动机铭牌参数的暂降临界切除时间计算方法 |
CN110333444B (zh) * | 2019-07-11 | 2023-10-20 | 中国石油大学(华东) | 一种基于电动机铭牌参数的暂降临界切除时间计算方法 |
CN111969924A (zh) * | 2020-08-18 | 2020-11-20 | 四川大学 | 一种交流接触器自适应动作控制方法 |
CN111969924B (zh) * | 2020-08-18 | 2022-03-29 | 四川大学 | 一种交流接触器自适应动作控制方法 |
CN113381401A (zh) * | 2021-05-31 | 2021-09-10 | 国网宁夏电力有限公司电力科学研究院 | 基于归一化能量并行计算临界切除时间的方法及系统 |
CN113381401B (zh) * | 2021-05-31 | 2022-06-07 | 国网宁夏电力有限公司电力科学研究院 | 基于归一化能量并行计算临界切除时间的方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN103280783B (zh) | 2015-07-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103280783B (zh) | 基于临界机组对的故障临界切除时间计算方法 | |
CN102761125B (zh) | 基于响应的电力系统暂态稳定性实时闭环系统及控制方法 | |
CN103746370B (zh) | 一种风电场可靠性建模方法 | |
CN108092304B (zh) | 基于相轨迹与相对动能变化率的电力系统暂态稳定性判别方法 | |
CN102510071B (zh) | 一种电网系统紧急控制方法及装置 | |
CN103500254B (zh) | 一种考虑地理分区边界优化的电网关键断面自动搜索方法 | |
CN103346558A (zh) | 一种基于广域测量系统的电力系统功角暂态失稳判别方法 | |
CN101794998A (zh) | 基于多机电力系统单一发电机电磁功率简洁表达方式的暂态稳定在线分析方法 | |
CN102445660B (zh) | 基于灰色Verhulst模型的发电机功角预测方法 | |
CN102185325A (zh) | 基于网络量测信息的电力系统暂态稳定性定量评价方法 | |
CN103473478A (zh) | 基于能量函数的电网暂态稳定性评估方法 | |
CN103632043A (zh) | 一种基于实测响应信息的电力系统主导失稳模式识别方法 | |
CN103995948B (zh) | 一种基于多项式模型的振荡中心电压预测方法 | |
CN102510072B (zh) | 一种电网系统暂态失稳判别方法 | |
CN105071380A (zh) | 一种电力系统广域同调辨识方法及其装置 | |
CN102570448B (zh) | 基于wams的电力系统分群及稳定裕度评估系统及方法 | |
CN102723710B (zh) | 一种用于电力系统暂态稳定分析的临界机群识别方法 | |
CN109995023A (zh) | 基于谱聚类区域划分的电力系统稳定性的判断方法和系统 | |
CN105021952A (zh) | 基于广域相角量测信息的大电网故障定位方法 | |
CN109617061A (zh) | 一种不对称故障情况下稳定裕度的快速评估方法 | |
CN108054768A (zh) | 基于主成分分析的电力系统暂态稳定评估方法 | |
CN105676013A (zh) | 广域测量系统的机组运行状态监视功能的测试系统及方法 | |
CN110969214B (zh) | 一种基于支持向量机综合模型的暂态安全域在线构建方法 | |
CN104505827A (zh) | 基于响应信息的复杂电力系统闭环控制方法 | |
CN104809512B (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: 20150722 Termination date: 20210509 |
|
CF01 | Termination of patent right due to non-payment of annual fee |