CN104766142B - 基于eeac和轨迹灵敏度的暂态稳定约束最优潮流计算方法 - Google Patents
基于eeac和轨迹灵敏度的暂态稳定约束最优潮流计算方法 Download PDFInfo
- Publication number
- CN104766142B CN104766142B CN201510190459.7A CN201510190459A CN104766142B CN 104766142 B CN104766142 B CN 104766142B CN 201510190459 A CN201510190459 A CN 201510190459A CN 104766142 B CN104766142 B CN 104766142B
- Authority
- CN
- China
- Prior art keywords
- failure
- transient
- load flow
- transient stability
- stability
- 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
Abstract
本发明公开了基于EEAC和轨迹灵敏度的暂态稳定约束最优潮流计算方法，属于电力系统暂态分析的技术领域。所述方法：以常规最优潮流解为初始运行点计算暂态初值并设置预想故障集；对预想故障集中的各故障进行暂态时域仿真、EEAC等值操作，对故障进行分类，将满足故障稳定裕度要求的故障从预想故障集中剔除；计算不满足稳定裕度要求的故障的轨迹灵敏度，给出不同类型下各故障的裕度灵敏度，并根据各故障的裕度灵敏度构建暂态稳定性约束；将所述暂态稳定性约束嵌入常规最优潮流模型后得到含有简单暂态稳定约束的最优潮流模型，避免了只采用功角作为暂态稳定判据时无法量化暂态稳定性，减少了计算复杂性并解决了收敛困难的问题。
Description
技术领域
本发明公开了基于EEAC(Extended Equal Area Criterion，扩展等面积准则)和轨迹灵敏度的暂态稳定约束最优潮流计算方法，属于电力系统暂态分析的技术领域。
背景技术
随着国民经济的快速发展和人民生活水平的不断提高，电力工业迅猛发展，电力新技术、新措施逐步投入应用，给社会各方面带来了显著的效益。但是，由于电力体制改革的不断深入，电网运营商为了追求更高的经济效益，使电力系统的运行点不断逼近安全稳定边界。各区域电网间互联加强以及新技术的应用，导致电网规模和复杂性日益增大，使电力系统的安全稳定问题越来越突出。电力系统中任意元件发生故障，都有可能引起事故扩大。如果电网结构不合理，或者安全自动装置不够健全，都有可能使系统陷入稳定危机，造成稳定破坏或大面积停电，给国民经济带来重大损失。现代社会一旦发生大停电，其经济损失可能远远超过市场竞争可能带来的经济收益甚至危及国家和社会安全。
含暂态稳定约束的最优潮流(Transient Stability Constrained OptimalPower Flow,TSCOPF)可以兼顾系统的暂态稳定性和经济性，实现了电力系统暂态稳定离线预防控制，成为解决这一问题的有效手段。
TSCOPF抽象为数学模型后，是一个复杂的大规模非线性微分代数混合约束优化问题。目前计算的基本解法分为两类。
第一类方法为联立求解法，一般是将暂态微分方程差分化变为代数方程后加入最优潮流中联立求解。如文献一《暂态稳定约束下的最优潮流》(中国电机工程学报，2005年第25卷第12期第12页)所述。这类方法原理清晰，但是问题规模随着故障数目和系统规模的增大，呈几何级数增长，易出现维数灾的问题。由于差分化加入最优潮流中后，最优化问题庞大，容易出现收敛困难的问题。另外，这类方法的暂态稳定判据一般为功角阈值，文献二《暂态稳定预防控制和优化新进展》(电力系统自动化，2004年第28卷第10期)所述，该判据存在两个理论上的问题，一是角度稳定是否就意味着系统稳定，二是如何选取最大相对摇摆角。
第二类方法是分解迭代求解法，通常将TSCOPF问题分解为暂态稳定评估和最优潮流两个子问题，将暂态稳定分析得到的约束式嵌入最优潮流问题中迭代求解。这类方法可有效处理多故障问题、减少计算维数。如文献三《基于单机无穷大母线等值和轨迹灵敏度的暂态稳定约束最优潮流》(中国电机工程学报，2011年第31卷第13期)所述，这类解法可以定义系统暂态稳定裕度，可得到系统暂态稳定量化信息，便于优化。这类方法通常采用功角门槛值作为暂态稳定约束，没有量化稳定裕度。文献三定义了一般不稳定和稳定情况下的暂态稳定裕度，构造稳定性约束，极度失稳情况下仍采用功角门槛值构造暂态稳定性约束，采用二次规划算法进行迭代计算。最后的算例结果仅使多个故障达到稳定，并未满足一定的稳定性要求。由于各种随机因素的存在，必须留有一定的暂态稳定裕度。
发明内容
本发明所要解决的技术问题是针对上述背景技术的不足，提供了基于EEAC和轨迹灵敏度的暂态稳定约束最优潮流计算方法，将暂态稳定约束最优潮流分解为暂态稳定评估、灵敏度分析、最优潮流三个子问题进行迭代求解，解决了现有暂态稳定约束最优潮流计算复杂、难以收敛，以功角阈值作为稳定判据不准确以及选取最大相对摇摆角难的技术问题。
本发明为实现上述发明目的采用如下技术方案：
基于EEAC和轨迹灵敏度的暂态稳定约束最优潮流计算方法，包括如下步骤：
A、以常规最优潮流解为初始运行点计算暂态初值并设置预想故障集；
B、对预想故障集中的各故障进行暂态时域仿真、EEAC等值操作，对故障进行分类，将满足故障稳定裕度要求的故障从预想故障集中剔除，仅当预想故障集中所有故障都被剔除时直接输出最优解；
C、对于步骤B中不满足稳定裕度要求的故障：计算各故障的轨迹灵敏度，给出不同类型下各故障的裕度灵敏度，并根据各故障的裕度灵敏度构建暂态稳定性约束；
D、将所述暂态稳定性约束嵌入常规最优潮流模型后得到含有简单暂态稳定约束的最优潮流模型，重新计算常规潮流解，返回步骤A。
进一步的，作为所述暂态稳定约束最优潮流计算方法的优化方案，步骤C中所述暂态稳定性约束为：η为各故障稳定裕度构成的向量，为各故障裕度灵敏度构成的向量，u为系统控制变量构成的向量，ΔuT为系统控制变量在各故障下的变化量构成的向量，ηmin为各故障稳定裕度最小值构成的向量，
各故障稳定裕度构成的向量η包括：一般失稳故障的稳定裕度ηu、稳定故障的稳定裕度ηs、极度失稳故障的稳定裕度ηeu：
ηeu＝-Pa(tmin)＝-min{PmE(t)-PeE(t),t＞tcl}，
其中：tu为等值后系统达到DSP点的时刻，ME为发电机的等值惯性时间常数，ωE(tu)、δE(tu)分别为tu时刻发电机的转速、功角，tr为等值后系统首次达到FEP点的时刻，Pa(tr)、δE(tr)为分别为tr时刻的不平衡功率、发电机功角，tmin为不平衡功率达到最小值的时刻，Pa(tmin)为tmin时刻的不平衡功率值，PmE(t)、PeE(t)分别为等值后发电机在t时刻的机械功率、电磁功率，tcl为故障清除时刻。
再进一步的，作为所述暂态稳定约束最优潮流计算方法的优化方案，步骤C中所述一般失稳故障的稳定裕度ηu，在其暂态稳定性约束中增加校正因子λ：
其中：λ＝0.8～0.9，为系统控制变量在一般失稳故障下变化量构成的向量，ηumin为一般失稳故障的稳定裕度最小值。
进一步的，作为所述暂态稳定约束最优潮流计算方法的优化方案，步骤B中所述EEAC等值操作具体为：将对各故障进行暂态时域仿真得到的功角曲线进行EEAC等值，发电机等值后的模型为：
其中：δE、ωE、PmE、PeE分别为发电机的等值功角、等值转速、等值机械功率、等值电磁功率，Pa为等值后的不平衡功率。
再进一步的，作为所述暂态稳定约束最优潮流计算方法的优化方案，步骤B中所述对故障进行分类，采用如下判据：
一般失稳故障的判据：
稳定故障的判据：
极度失稳故障的判据：Pa(t)＞0,
其中：Pa(tu)、PmE(tu)、PeE(tu)分别为tu时刻的不平衡功率、发电机机械功率、发电机电磁功率，Pa(t)、ωE(t)分别为t时刻的不平衡功率、发电机转速，t0为故障开始时刻，PmE(tr)、PeE(tr)、ωE(tr)分别为tr时刻发电机的机械功率、电磁功率、转速，tend为暂态时域仿真结束时刻。
更进一步的，作为所述暂态稳定约束最优潮流计算方法的优化方案，步骤D中采用非线性原对偶内点法求解含有简单暂态稳定约束的最优潮流模型。
本发明采用上述技术方案，具有以下有益效果：
(1)采用EEAC法进行暂态稳定评估，考虑不平衡功率以及发电机相角构建了极度失稳、一般失稳和稳定情况下的暂态稳定判据，进而量化表示各故障的暂态稳定裕度，避免了只采用功角作为暂态稳定判据时无法量化暂态稳定性以及构建的暂态稳定约束不可靠的问题；
(2)针对各故障计算轨迹灵敏度，进而确定各故障的裕度灵敏度，再由轨迹灵敏度以及裕度灵敏度确定二次规划问题中暂态稳定约束的表达式，指导最优控制方向；
(3)本发明构建的将所述暂态稳定性约束嵌入常规最优潮流模型后得到含有简单暂态稳定约束的最优潮流模型，将一般失稳故障调整为稳定故障，将极度失稳故障调整为一般失稳故障，在以后的迭代中再将一般失稳故障调整为稳定故障，以满足给定的稳定性要求，实现系统由轻微失稳向稳定的调整，调整的过程中控制量变化较小，调整精确性较高且过度调整不明显；
(4)将暂态稳定约束最优潮流分解为暂态稳定评估、灵敏度分析、最优潮流三个子问题进行迭代求解，构建的各故障暂态稳定性约束嵌入最常规潮流计算后，其计算量与常规最优潮流计算几乎相同，减少了计算复杂性并解决了收敛困难的问题。
附图说明
图1是EEAC等值后一般失稳故障的等值功角曲线示意图；
图2是EEAC等值后稳定故障的等值功角曲线示意图；
图3是EEAC等值后极度失稳故障的等值功角曲线示意图；
图4是计算流程图。
具体实施方式
下面结合附图对发明的技术方案进行详细说明。本发明的思路是：将暂态稳定约束最优潮流分为三个子问题，包括暂态稳定评估问题、灵敏度分析问题和最优潮流问题，并将三个子问题迭代求解。如图4所示，具体按照以下步骤：
步骤1、基于EEAC对各预想故障进行暂态稳定评估：
步骤101、计算系统常规最优潮流，即计算不含暂态稳定约束的最优潮流，以此作为初始运行点。
建立发电机模型，负荷模型，计算暂态初值，对各故障进行暂态时域仿真。暂态稳定计算的数学模型为：
G(x(t),y(t),u)＝0 t∈T (4)
式(1)至式(4)中，x为系统状态量向量，包括发电机功角和转速等；y为系统代数量向量，包括系统节点电压的幅值和相角等；u为系统控制量向量，包括发电机有功出力和无功出力等；为故障发生前一时刻， 为故障发生时刻，tcl为故障清除时刻，tend为暂态仿真结束时刻；式(1)为潮流方程；式(2)为发电机暂态初值方程；式(3)和式(4)描述系统的暂态运行过程，为一组微分代数方程。
然后将仿真后功角曲线进行EEAC等值。如文献四《运动稳定性量化理论—非自治非线性多刚体系统的稳定性分析》(南京：江苏科学技术出版社，1999)所述，系统发电机等值后模型为：
式(5)中：ME为等值惯性时间常数；δE为等值功角；ωE为等值转速；PmE为等值机械功率；PeE为等值电磁功率；Pa为等值后的不平衡功率。
步骤102、EEAC等值后，将故障分为一般失稳、稳定、极度失稳三种情况。
在仿真时间内，若系统无DSP点，只有FEP点，则系统稳定；若系统达到DSP点，则系统一般失稳；若系统既无DSP点也无FEP点，则系统极度失稳。下面具体给出各稳定情况的判据以及裕度表达式：
(1)一般失稳故障的判据为：
式(6)中：tu为EEAC等值后系统达到DSP点的时刻，一般失稳故障的等值功角曲线示意图如图1所示，其稳定裕度定义为：
(2)系统稳定的判据为：
式(8)中：tr为系统EEAC等值后第一次达到FEP点的时刻；δE(tr)为系统EEAC等值后第一次达到FEP点时的等值功角，稳定故障的等值功角曲线示意图如图2所示，其稳定裕度可近似表达为：
式(9)中：δE(tu)是系统临界失稳时在tu时刻的等值功角。
(3)系统极度失稳的判据为：
由于极度失稳故障不存在DSP点和FEP点，极度失稳故障的等值功角曲线示意图如图3所示，所以将最小不平衡功率定义为其稳定裕度：
ηeu＝-Pa(tmin)＝-min{PmE(t)-PeE(t),t＞tcl} (11)
式(6)至式(11)中：Pa(tu)、PmE(tu)、PeE(tu)分别为tu时刻的不平衡功率、发电机机械功率、发电机电磁功率，Pa(t)、ωE(t)分别为t时刻的不平衡功率、发电机转速，t0为故障开始时刻，PmE(tr)、PeE(tr)、ωE(tr)分别为tr时刻发电机的机械功率、电磁功率、转速，tmin代表故障后不平衡功率达到最小值的时刻，Pa(tmin)代表故障后不平衡功率的最小值。
步骤2、计算轨迹灵敏度以及稳定裕度灵敏度
步骤201、针对上一步中不满足稳定裕度要求的故障，基于上述暂态时域仿真的数据，计算轨迹灵敏度初值。
当系统控制量发生变化时，会造成潮流分布的改变，所以首先对式(1)求偏导，得到静态潮流结果对于控制量的灵敏度然后对式(2)求偏导，得到暂态初值对于控制量的灵敏度具体公式如下：
由于状态变量不会突变，所以其灵敏度也不会突变，即而代数变量会发生突变，所以 可由式(4)计算得到。所以进行整理后，轨迹灵敏度的初值计算方程为：
步骤202、计算轨迹灵敏度和裕度灵敏度。
轨迹灵敏度方程可由式(3)和式(4)对控制量u求偏导得到：
式(14)中：Fx(t)、Fy(t)、Fu(t)、Gx(t)、Gy(t)和Gu(t)为随系统轨迹变化的时变偏导矩阵。
为了刻画控制变量对系统稳定裕度的影响，基于上述轨迹灵敏度，引入裕度灵敏度。对于不同的故障类型，有不同的裕度灵敏度表达式。
(1)对于极度失稳故障，其裕度灵敏度为：
(2)对于一般失稳故障，其裕度灵敏度为：
(3)对于稳定故障，当其稳定裕度较小而不满足稳定性要求时，仍需调整控制变量，所以其裕度灵敏度可表达为：
步骤203、基于上述各故障的裕度灵敏度，可构造暂态稳定性约束，可简写为如下形式：
式(18)中：η为各个故障稳定裕度构成的向量，为各故障裕度灵敏度构成的向量，u为系统控制变量构成的向量，ΔuT为系统控制变量在各故障下变化量构成的向量，ηmin为各故障稳定裕度最小值构成的向量，各故障稳定裕度构成的向量η包括：式(7)所示一般失稳故障的稳定裕度ηu、式(9)所示稳定故障的稳定裕度ηs、式(11)所示极度失稳故障的稳定裕度ηeu，ωE(tu)、δE(tu)分别为tu时刻发电机的转速、功角，tr为等值后系统首次达到FEP点的时刻，Pa(tr)、δE(tr)为分别为tr时刻的不平衡功率、发电机功角，tmin为不平衡功率达到最小值的时刻，Pa(tmin)为tmin时刻的不平衡功率值，PmE(t)、PeE(t)分别为等值后发电机在t时刻的机械功率、电磁功率，tcl为故障清除时刻。
一般失稳故障的稳定性约束目标是将一般失稳故障调整为稳定。由于极度失稳故障的裕度是根据故障后最小不平衡功率定义的，所以其稳定性约束的目标是将极度失稳调整为一般失稳故障，以后的迭代中再将一般失稳故障调整为稳定。在极度失稳故障和一般失稳故障的稳定性约束中，一般取ηmin＝0，在稳定故障中，根据实际电网需求，给定ηmin，使发生故障时系统稳定裕度不会过小，符合预防控制的要求。
步骤3、计算含简单暂态稳定约束的最优潮流
步骤301、由上述稳定性约束形成含简单暂态稳定约束的最优潮流模型：
Min f(u) (19)
式(19)为目标函数，如发电费用函数；式(21)为故障前稳态时需满足的运行约束方程，式(22)为稳定性约束方程，其中k＝1,…,nF代表第k个故障，nF为故障的总数目。
为提高计算准确性，对于较为严重的一般失稳故障，迭代时可在稳定性约束中增加一个校正因子λ＝0.8～0.9，构成校正的稳定性约束为：
该约束使迭代后系统为轻微失稳。而其余迭代中，校正因子λ＝1。采用轨迹灵敏度计算控制量时，会出现过度调整的现象。由于轨迹灵敏度是基于原轨迹的线性灵敏度，当控制量的变化较小时，应用灵敏度调整稳定裕度误差较小，精确性较高。通过这一处理方法，可使在预防控制过程中，系统从失稳调为稳定时，系统是由轻微失稳调整为稳定，控制量变化较小，调整精确性较高，过度调整不明显。
效果验证：
为验证本发明算法的有效性，针对IEEE10机39节点标准算例进行了测试，分别计算含单一故障和多故障的TSCOPF，并进行比较分析。发电机采用经典二阶模型，负荷采用恒阻抗模型，仿真时间为4s(考虑多摆失稳)，仿真步长为0.01s。flag是标识符，用于记录系统发生故障后的暂态稳定情况，0表示稳定，1表示一般失稳，2表示极度失稳。令符号L(i,j)表示线路i-j在母线i处0s时刻发生三相短路故障，tcl时刻切除故障线路。下面选取故障A1和故障A2，故障A1为L(29,28)，切除故障时间为0.32s，为极度失稳故障；故障A2为L(2,1)，切除故障时间为0.24s，为一般失稳故障。故障A1中38号机为临界群，其余为剩余群；故障A2中39号机为剩余群，其余为临界群，两者为不同失稳模式的故障。表1和表2分别为单独考虑故障A1和故障A2时计算TSCOPF的过程，表3为同时考虑故障A1和A2计算TSCOPF的过程，表中CCT表示临界切除时间。本算例分析中将发电机有功出力作为控制变量，取ηmin＝0.1。当一般失稳故障裕度小于-1时，取校正因子λ＝0.8。
表1故障A1的TSCOPF迭代过程
表1中，由原对偶内点法计算常规最优潮流，得到初始运行点。对故障A1仿真，判断为极度失稳故障，经过第1次迭代优化后，转为一般失稳故障，是严重的一般失稳故障。为了避免出现过度调整的情况，将校正的稳定性约束嵌入最优潮流计算，得到第2次迭代结果，此时仍然失稳，但裕度较大，经第3次迭代后，系统达到稳定。此时，系统仍不满足稳定性要求，经第4次迭代，达到稳定性要求。
表2故障A2的TSCOPF迭代过程
| 发电机 | 初始OPF | 第1次迭代 | 第2次迭代 | 第3次迭代 |
| 30(MW) | 243.96 | 242.94 | 241.45 | 250.89 |
| 31(MW) | 569.62 | 555.87 | 548.57 | 548.2 |
| 32(MW) | 646.06 | 648.39 | 647.35 | 647.03 |
| 33(MW) | 636.67 | 631.85 | 629.25 | 628.25 |
| 34(MW) | 513.11 | 509.66 | 507.5 | 506.94 |
| 35(MW) | 655.63 | 650.68 | 648.69 | 647.38 |
| 36(MW) | 563.31 | 558.3 | 556.29 | 555.09 |
| 37(MW) | 539.11 | 531.92 | 530.51 | 527.27 |
| 38(MW) | 837.3 | 812.29 | 813.47 | 809.78 |
| 39(MW) | 981.08 | 1041.86 | 1060.39 | 1062.3 |
| flag | 1 | 1 | 0 | 0 |
| 裕度(p.u.) | -1.4182 | -0.2537 | 0.0355 | 0.1321 |
| 费用($/h) | 61656.05 | 61689.91 | 61711.22 | 61717.5 |
| CCT(s) | 0.242 | 0.244 |
表3故障A1+A2的TSCOPF迭代过程
故障A2为多摆失稳故障，表2给出故障A2的TSCOPF计算过程，与故障A1类似，经3次迭代可达到稳定性要求。表3中，同时处理含故障A1和故障A2这两个不同失稳模式故障的暂态稳定最优潮流问题，经4次迭代可达到稳定性要求。对比表1、2和3可见，考虑多故障的发电费用略高于考虑单故障的发电费用。
由此可见，本发明提出的方法是合理有效的。
综上所述，本发明具有以下有益效果：
(1)采用EEAC法进行暂态稳定评估，考虑不平衡功率以及发电机相角构建了极度失稳、一般失稳和稳定情况下的暂态稳定判据，进而量化表示各故障的暂态稳定裕度，避免了只采用功角作为暂态稳定判据时无法量化暂态稳定性以及构建的暂态稳定约束不可靠的问题；
(2)针对各故障计算轨迹灵敏度，进而确定各故障的裕度灵敏度，再由轨迹灵敏度以及裕度灵敏度确定二次规划问题中暂态稳定约束的表达式，指导最优控制方向；
(3)本发明构建的将所述暂态稳定性约束嵌入常规最优潮流模型后得到含有简单暂态稳定约束的最优潮流模型，将一般失稳故障调整为稳定故障，将极度失稳故障调整为一般失稳故障，在以后的迭代中再将一般失稳故障调整为稳定故障，以满足给定的稳定性要求，实现系统由轻微失稳向稳定的调整，调整的过程中控制量变化较小，调整精确性较高且过度调整不明显；
(4)将暂态稳定约束最优潮流分解为暂态稳定评估、灵敏度分析、最优潮流三个子问题进行迭代求解，构建的各故障暂态稳定性约束嵌入最常规潮流计算后，其计算量与常规最优潮流计算几乎相同，减少了计算复杂性并解决了收敛困难的问题。
本发明所述方案不局限于具体实施例中涉及的常规潮流模型，本领域的技术人员在本发明技术方案的指导下，能够将本发明所述方案推广应用到其它常规潮流模型中，其它应用于任意常规潮流模型中的基于EEAC和轨迹灵敏度的暂态稳定约束最优潮流计算方法都落入本发明的保护范围。
Claims (5)
1.基于EEAC和轨迹灵敏度的暂态稳定约束最优潮流计算方法，其特征在于，包括如下步骤：
A、以常规最优潮流解为初始运行点计算暂态初值并设置预想故障集；
B、对预想故障集中的各故障进行暂态时域仿真、EEAC等值操作，对故障进行分类，将满足故障稳定裕度要求的故障从预想故障集中剔除，仅当预想故障集中所有故障都被剔除时直接输出最优解；
C、对于步骤B中不满足稳定裕度要求的故障：计算各故障的轨迹灵敏度，给出不同类型下各故障的裕度灵敏度，并根据各故障的裕度灵敏度构建暂态稳定性约束，
暂态稳定性约束为：η为各故障稳定裕度构成的向量，为各故障裕度灵敏度构成的向量，u为系统控制变量构成的向量，ΔuT为系统控制变量在各故障下的变化量构成的向量，ηmin为各故障稳定裕度最小值构成的向量，
各故障稳定裕度构成的向量η包括：一般失稳故障的稳定裕度ηu、稳定故障的稳定裕度ηs、极度失稳故障的稳定裕度ηeu：
ηeu＝-Pa(tmin)＝-min{PmE(t)-PeE(t),t＞tcl}，
其中：tu为等值后系统达到DSP点的时刻，ME为发电机的等值惯性时间常数，ωE(tu)、δE(tu)分别为tu时刻发电机的转速、功角，tr为等值后系统首次达到FEP点的时刻，Pa(tr)、δE(tr)为分别为tr时刻的不平衡功率、发电机功角，tmin为不平衡功率达到最小值的时刻，Pa(tmin)为tmin时刻的不平衡功率值，PmE(t)、PeE(t)分别为等值后发电机在t时刻的机械功率、电磁功率，tcl为故障清除时刻；
D、将所述暂态稳定性约束嵌入常规最优潮流模型后得到含有简单暂态稳定约束的最优潮流模型，重新计算常规潮流解，返回步骤A。
2.根据权利要求1所述的基于EEAC和轨迹灵敏度的暂态稳定约束最优潮流计算方法，其特征在于，步骤C中所述一般失稳故障的稳定裕度ηu，在其暂态稳定性约束中增加校正因子λ：
其中：λ＝0.8～0.9，为系统控制变量在一般失稳故障下变化量构成的向量，ηumin为一般失稳故障的稳定裕度最小值。
3.根据权利要求1所述的基于EEAC和轨迹灵敏度的暂态稳定约束最优潮流计算方法，其特征在于，步骤B中所述EEAC等值操作具体为：将对各故障进行暂态时域仿真得到的功角曲线进行EEAC等值，发电机等值后的模型为：
其中：δE、ωE、PmE、PeE分别为发电机的等值功角、等值转速、等值机械功率、等值电磁功率，Pa为等值后的不平衡功率。
4.根据权利要求3所述的基于EEAC和轨迹灵敏度的暂态稳定约束最优潮流计算方法，其特征在于，步骤B中所述对故障进行分类，采用如下判据：
一般失稳故障的判据：
稳定故障的判据：
极度失稳故障的判据：
其中：Pa(tu)、PmE(tu)、PeE(tu)分别为tu时刻的不平衡功率、发电机机械功率、发电机电磁功率，Pa(t)、ωE(t)分别为t时刻的不平衡功率、发电机转速，t0为故障开始时刻，PmE(tr)、PeE(tr)、ωE(tr)分别为tr时刻发电机的机械功率、电磁功率、转速，tend为暂态时域仿真结束时刻。
5.根据权利要求1至4中任意一项所述的基于EEAC和轨迹灵敏度的暂态稳定约束最优潮流计算方法，其特征在于，步骤D中采用非线性原对偶内点法求解含有简单暂态稳定约束的最优潮流模型。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510190459.7A CN104766142B (zh) | 2015-04-21 | 2015-04-21 | 基于eeac和轨迹灵敏度的暂态稳定约束最优潮流计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510190459.7A CN104766142B (zh) | 2015-04-21 | 2015-04-21 | 基于eeac和轨迹灵敏度的暂态稳定约束最优潮流计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104766142A CN104766142A (zh) | 2015-07-08 |
| CN104766142B true CN104766142B (zh) | 2019-03-05 |
Family
ID=53647956
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510190459.7A Active CN104766142B (zh) | 2015-04-21 | 2015-04-21 | 基于eeac和轨迹灵敏度的暂态稳定约束最优潮流计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104766142B (zh) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105184418B (zh) * | 2015-09-23 | 2019-04-02 | 重庆大学 | 基于灵敏度和约束等值的互联电网最优潮流计算方法 |
| CN106356885A (zh) * | 2016-09-19 | 2017-01-25 | 国家电网公司 | 一种风电并网系统暂态功角前二摆稳定性判断方法 |
| CN106972490B (zh) * | 2017-06-01 | 2019-10-22 | 清华-伯克利深圳学院筹备办公室 | 输电网预想事故分析的方法及装置 |
| CN108551177B (zh) * | 2018-04-02 | 2020-05-12 | 河海大学 | 基于灵敏度分析的直流受端系统暂态切负荷控制优化方法 |
| CN111125880B (zh) * | 2019-11-25 | 2022-07-22 | 国网四川省电力公司电力科学研究院 | 一种暂态稳定视角下电力系统仿真数据生成方法 |
| CN114021387B (zh) * | 2021-11-26 | 2023-05-02 | 国网四川省电力公司电力科学研究院 | 一种电力系统多场景智能仿真方法、系统及电子设备 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101969199A (zh) * | 2010-08-26 | 2011-02-09 | 天津大学 | 用于暂态功角稳定风险评估的故障损失估计方法 |
| CN102324743A (zh) * | 2011-09-21 | 2012-01-18 | 国网电力科学研究院 | 电力系统在线暂态安全稳定评估预想故障筛选方法 |
| CN102957166A (zh) * | 2012-09-13 | 2013-03-06 | 中国电力科学研究院 | 一种基于轨迹灵敏度的风电配置比例快速计算方法 |
| CN103457279A (zh) * | 2013-09-23 | 2013-12-18 | 广东电网公司电力调度控制中心 | 大规模电网的分区d-statcom综合优化配置方案的确定方法 |
-
2015
- 2015-04-21 CN CN201510190459.7A patent/CN104766142B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101969199A (zh) * | 2010-08-26 | 2011-02-09 | 天津大学 | 用于暂态功角稳定风险评估的故障损失估计方法 |
| CN102324743A (zh) * | 2011-09-21 | 2012-01-18 | 国网电力科学研究院 | 电力系统在线暂态安全稳定评估预想故障筛选方法 |
| CN102957166A (zh) * | 2012-09-13 | 2013-03-06 | 中国电力科学研究院 | 一种基于轨迹灵敏度的风电配置比例快速计算方法 |
| CN103457279A (zh) * | 2013-09-23 | 2013-12-18 | 广东电网公司电力调度控制中心 | 大规模电网的分区d-statcom综合优化配置方案的确定方法 |
Non-Patent Citations (3)
| Title |
|---|
| 基于二阶轨迹灵敏度的暂态稳定约束最优潮流计算;潮铸;《中国优秀硕士学位论文全文数据库》;20111215(第12期);正文第15-53页 |
| 暂态能量约束的最优潮流研究;夏世威;《中国博士学位论文全文数据库》;20150115(第1期);正文第63-74页 |
| 电力系统暂态稳定最优控制算法研究;黄玉龙;《中国博士学位论文全文数据库》;20121115(第11期);正文第16-21页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104766142A (zh) | 2015-07-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104766142B (zh) | 基于eeac和轨迹灵敏度的暂态稳定约束最优潮流计算方法 | |
| CN107093895B (zh) | 基于预想故障集自动筛选的在线暂态安全稳定评估方法 | |
| CN104156542B (zh) | 一种基于隐式投影的有源配电系统稳定性仿真方法 | |
| CN104716646B (zh) | 一种基于注入电流的节点耦合度分析方法 | |
| CN108183512A (zh) | 一种接入新能源的电力系统的可靠性评估方法 | |
| CN108667005B (zh) | 一种计及新能源影响的电网静动态结合脆弱性评估方法 | |
| CN102403720A (zh) | 一种基于暂态电压安全裕度的超实时重合时序整定方法 | |
| CN105610192A (zh) | 考虑大规模风电接入的在线风险评估方法 | |
| CN109995023A (zh) | 基于谱聚类区域划分的电力系统稳定性的判断方法和系统 | |
| CN103700036A (zh) | 一种适于电力系统多时间尺度的暂态稳定性投影积分方法 | |
| CN103915839B (zh) | 一种用于分析含风电电力系统随机稳定性的方法 | |
| Kontis et al. | Dynamic equivalencing of active distribution grids | |
| CN114970154A (zh) | 一种暂态电压支撑能力量化评估指标构建方法 | |
| CN106056305A (zh) | 一种基于状态聚类的发电系统可靠性快速评估方法 | |
| Zhang et al. | Impact of dynamic load models on transient stability-constrained optimal power flow | |
| Bento et al. | Analysis of the load growth direction variation in the dynamic security assessment | |
| CN110348598A (zh) | 一种电网网架规划方法 | |
| CN105140957A (zh) | 基于风电场和光伏电站聚合模型的机电振荡模式估算方法 | |
| CN102437643B (zh) | 确定多断面各输电断面可用输电能力的方法 | |
| CN113987848A (zh) | 一种负荷构成智能识别及精准负荷建模方法及系统 | |
| CN114298478A (zh) | 一种风电并网系统小扰动稳定辨识方法及系统 | |
| Ali et al. | Online voltage instability detection of distribution systems for smart-grid applications | |
| Zhao et al. | A novel SIME and sensitivity based model and algorithm of transient stability constrained optimal power flow | |
| Zhou et al. | Research on Key Technologies of Digital Twins System in Distribution Networks | |
| Bento et al. | Computing the worst case scenario for electric power system dynamic security assessment |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| EXSB | Decision made by sipo to initiate substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |