电力系统失稳模式易变性的快速评估方法
【技术领域】
[0001] 本发明属于电力系统及其自动化技术领域，更准确地说，本发明涉及一种电力系 统失稳模式易变性的快速评估方法。
【背景技术】
[0002] 电力系统失稳模式的准确识别是进行电力系统暂态功角稳定量化分析和控制的 基础。现有的电力系统失稳模式的识别方法主要有逐步积分法及启发式方法。前者通过观 察由数值积分法得到的系统的运动轨线来判断系统的失稳模式；后者根据暂态过程特定时 刻的部分参量预估电力系统失稳模式。
[0003] 逐步积分法的优势在于系统轨线信息获得的准确性，然而它将消耗较大的计算 量；启发式方法的优势在于失稳模式识别速率的提升，不过由于未能计及全部暂态稳定信 息，识别可靠性并不理想。
[0004] 工程实用的失稳模式判别方法还要求同时满足判别准确性和快速性。前者要求样 本信息必须遍历研究算例暂态稳定全过程，后者则要求以尽量小的计算代价获得全部样本 信息。EEAC开发了识别失稳模式的有效方法，使得对于任何规模的工程系统都只需要对3 个或更少的候选失稳模式进行评估。
[0005] 然而，面对现代电力系统规模及模型的急剧增加，越来越多研究算例的失稳模式 呈现复杂特性，它受故障清除时间等因素的影响。对于研究算例失稳模式变化复杂特性 的详细分析要求实现所有感兴趣故障清除时间下全部失稳模式集的准确识别，现有的基于 EEAC理论失稳模式识别方法的计算量使得仿真分析仍然受到挑战。
[0006] 对于分析速度与精度近乎矛盾的要求使得失稳模式变化复杂特性的评估问题似 乎陷入死锁。而死锁的破解则在于针对具体研究算例，寻求适应的判别方法，从而在算例空 间中协调精度和速度。
[0007] 专利申请"电力系统暂态稳定评估预想故障快速强壮分类方法"（【申请号】 201410271454. 2)和"电力系统暂态稳定最优切机控制策略搜索的简化方法"（【申请号】 201410398838. 0)提供了较好的范例。它们利用简化程度不同的各算法分析结果间的相似 程度与算例时变程度的内在联系，为算例时变程度不同的各算例匹配恰当的分析算法，从 算例空间角度协调求解速度和精度。
[0008] 循着这一思路，本发明应用大步长泰勒级数展开获取研究算例各机转子角轨迹， 并综合恰当时间断面处各机转子角间隙信息快速评估研究算例失稳模式易变性，以期省却 对于失稳模式不变算例的逐个感兴趣故障清除时间下各失稳模式的详细识别。
【发明内容】
[0009] 本发明目的是：为协调电力系统失稳模式易变性评估的准确性和快速性，基于大 步长泰勒级数展开、综合由此获得的恰当时间断面处各机转子角间隙信息，提供一种电力 系统失稳模式易变性的快速评估方法。
[0010] 具体地说，本发明是采用以下技术方案实现的，包括以下步骤：
[0011] 1)选取预想故障全集中的一个算例，针对该算例，设置其故障清除时间为0.8s， 并由大步长泰勒级数展开获得该算例故障中、故障后各机转子角轨迹；
[0012] 2)选取η个时刻作为故障清除后各机组暂态稳定态势发展较为/完全明显的η个 时间断面，分别观察每一时间断面的各机组转子角，获取各时间断面最大、次最大转子角间 隙，并由此求得表征各时间断面各机组转子角间隙信息的指标h、k 2、…、kn;
[0013] 3)以各时间断面各机组转子角间隙信息指标kp k2、…、kn中的最大值K表征该 算例时变程度，K = max {kpk；；，···，!〇,并作如下判断：若K小于预先设定阈值ε i，则认为 该算例时变程度较弱、其失稳模式不随故障清除时间而变，且由识别得到的最大转子角间 隙将各机组分为领前、余下两群；否则，则认为该算例时变程度较强、其失稳模式随故障清 除时间而变，针对所有感兴趣的故障清除时间、由IEEAC算法对该算例依次进行详细仿真， 获得失稳模式随故障清除时间的演化情形；
[0014] 4)对该算例失稳模式易变性的评估结束，若预想故障全集中所有算例都已完成失 稳模式易变性评估，则结束本方法，否则取预想故障全集中的下一个算例执行步骤1)。
[0015] 上述技术方案的进一步特征在于，所述步骤1)中，分别以0. 4s、0.1 s为步长由泰 勒级数展开获得该算例故障中、故障后各机转子角轨迹。
[0016] 上述技术方案的进一步特征在于，所述步骤2)中，η = 5,并选取L= I. 2s、12 = I. 3s、t3= 1.4s、t4= 1.5s、t5= 1.6s五个时刻作为故障清除后各机组暂态稳定态势发展 较为/完全明显的五个时间断面，观察每一时间断面的各机组转子角。
[0017] 上述技术方案的进一步特征在于，所述步骤2)中，表征各时间断面各机转子角间 隙信息的指标，按以下公式计算：
[0018]
[0019] 上式中A1代表第i个时间断面处各机组转子角间隙信息指标，Gu代表第i个时 间断面处最大转子角间隙值，G il代表第i个时间断面处次最大转子角间隙值。
[0020] 上述技术方案的进一步特征在于，所述步骤3)中，阈值ε i设定为3. 4%。
[0021] 本发明的有益效果如下：本发明基于大步长泰勒级数展开技术，综合恰当的多个 时间断面处对应于较大故障清除时间的各机转角间隙信息，评估研究算例失稳模式随故障 清除时间的变化特性。本发明适用于所有工程实际系统，大量实验表明：应用本发明所述方 法可实现研究算例失稳模式随故障清除时间变化特性的全部可靠识别，且其计算代价小于 现有逐步积分方法计算量的10%，由此，本发明对于协调电力系统失稳模式易变性评估的 准确性和快速性，以深入理解暂态功角稳定本质、辅助控制策略制定具有重大的理论和工 程意义。
【附图说明】
[0022] 图1为本发明方法的流程图。
[0023] 图2为本发明实施例提供的预想故障全集中全体算例时变程度指标K按升序排列 情况示意图。
[0024] 图3为图2在K e [0·0,0·4]的局部示意图。
【具体实施方式】
[0025] 下面参照附图并结合实例对本发明作进一步详细描述。
[0026] 图1中步骤1描述的是研究算例各机转子角轨迹的获得。在预想故障全集中选取 一个算例，并针对该算例，设置其故障清除时间为〇. 8s，并由大步长泰勒级数展开获得该算 例故障中、故障后各机转子角轨迹（如可分别以〇. 4s、0.1 s为步长由泰勒级数展开获得该 算例故障中、故障后各机转子角轨迹）。
[0027] 图1中步骤2描述的是各机转子角间隙信息的提取。选取η个时刻作为故障清除 后各机组暂态稳定态势发展较为/完全明显的η个时间断面（如可选取t 1= I. 2s、12 = I. 3s、t3