CN107132762B - 基于预想故障集自动筛选的在线静态安全评估方法 - Google Patents

Abstract

本发明公开了一种基于预想故障集自动筛选的在线静态安全评估方法，属于电力系统安全分析技术领域。本发明分别采用节点有功注入对支路有功的灵敏度和节点无功注入对节点电压的灵敏度，基于前一次预想故障下设备过载安全和静态电压安全的评估结果，以及当前运行状态下节点有功、无功注入相对于前一次静态安全评估所对应运行状态的变化量，估算当前运行状态下静态安全评估结果，只针对设备过载安全和静态电压安全裕度估算值较小的预想故障进行基于潮流计算的静态安全评估，并依据基于潮流计算的静态安全评估结果与估算值的差异，采用后验校正法对估算参数进行修正。本发明能够实现基于预想故障集自动筛选的在线静态安全快速评估。

Description

基于预想故障集自动筛选的在线静态安全评估方法

技术领域

本发明属于电力系统安全分析技术领域，更准确地说，本发明涉及一种适用于大电网预想故障集自动筛选的在线静态安全评估的方法。

背景技术

在线静态安全评估是EMS必备的一项应用功能，目前有两种技术方案：

第一种是由电网调度运行控制人员根据经验设置少量的预想故障，基于潮流计算进行静态安全评估。这种方案保证了静态安全评估的快速性，但当人工设置的预想故障集与电网实际运行状态的静态安全特性失配程度较大时，其静态安全评估结果不能反映电网的静态安全问题，给电网运行带来安全隐患。

第二种是基于电网的实际运行状态按照要求设置全部的预想故障，建立综合反映电网设备过载和电压安全的指标，通过评估指标对各支路导纳变化的灵敏度来反映相应支路故障对电网设备过载和电压安全的影响，实现对所有支路故障按严重程度排序，再依次针对排序在前位的严重故障进行基于潮流计算的静态安全评估，直至开断某条支路后设备过载和静态电压安全裕度大于设定门槛值为止。这种方案克服了预想故障集依赖人工设置、存在安全隐患的问题，但是由于不能预先确定需要进行基于潮流计算的静态安全评估的严重故障数目，需要通过设置静态安全裕度门槛值保证严重的预想故障都得到基于潮流计算的静态安全评估，因此，计算速度难以保障，稳定性也不够。此外，由于只对严重的预想故障进行了基于潮流计算的静态安全评估，其它预想故障没有静态安全评估结果，在电网存在静态安全问题需要进行预防控制决策计算时，难以避免预防控制决策实施后可能会导致没有进行静态安全评估的预想故障无法满足静态安全要求，给电网运行带来安全隐患。

发明内容

本发明目的是：针对上述现有技术中的不足，本发明提出一种不需要人工设置预想故障集且满足精度和快速性综合要求的数万个节点大电网在线静态安全评估方法。

本发明的基本原理在于：静态安全评估基于的是非线性潮流方程，若基于某个电网运行状态下节点有功注入对支路有功的灵敏度和节点无功注入对节点电压的灵敏度，直接计算相对于该电网运行状态变化后的支路有功变化量和节点电压变化量，其计算精度与反映该电网运行状态的潮流方程的非线性程度以及两个运行状态之间的差异程度二者紧密相关。若潮流方程的非线性程度不高，即使两个运行状态之间的差异程度较大，其计算精度也会比较高；若两个运行状态之间的差异程度不大，即使潮流方程的非线性程度较高，其计算精度仍会比较高。只有在潮流方程的非线性程度和两个运行状态之间的差异程度都较高的情况下，其计算精度才会比较差。为了应对这种情况，本发明引入估算参数，并采用后验校正法对估算参数进行动态修正以提高估算精度。

具体地说，本发明是采用以下技术方案实现的，包括以下步骤：

1)设基于预想故障集自动筛选的在线静态安全评估首次启动的电力系统运行时刻为t₀，t₀时刻电力系统运行状态为S₀，基于S₀，根据《电力系统安全稳定导则》中第一道安全防线静态安全评估预想故障集的设置要求，生成预想故障集，记为F₀，进入步骤2)；

2)针对S₀下F₀中所有预想故障分别进行基于潮流计算的静态安全量化评估，得到t₀时刻的在线静态安全评估结果，进入步骤3)；

3)设电力系统当前运行时刻为t₁，t₁时刻电力系统运行状态为S₁，基于S₁，根据《电力系统安全稳定导则》中第一道安全防线静态安全评估预想故障集的设置要求，生成预想故障集，记为F₁，进入步骤4)；

4)先将F₁中不属于F₀的预想故障以及F₁中引发S₁下电力系统解列的预想故障作为需要进行基于潮流计算的静态安全量化评估的预想故障集，记为F_c，再将F₁中不属于F_c的预想故障作为需要进行基于估算的静态安全量化评估的预想故障集，记为F_e，进入步骤5)；

5)分别针对F_e中的各个预想故障，基于其在S₀下的静态安全稳定量化评估结果以及S₀与S₁之间的差异，估算其在S₁下的静态安全稳定量化评估结果，将估算结果直接作为裕度估算值大于筛选门槛的预想故障在S₁下的静态安全稳定量化评估结果，将裕度估算值小于等于筛选门槛的预想故障添加到F_c中，并将添加到F_c中的预想故障组成的集合记为F_ce，进入步骤6)；

6)针对S₁下F_c中所有预想故障分别进行基于潮流计算的静态安全量化评估，得到静态安全量化评估结果，进入步骤7)；

7)将步骤5)得到的S₁下F₁中不属于F_c的静态安全量化评估结果和步骤6)得到的S₁下F_c的静态安全量化评估结果，作为t₁时刻的在线静态安全评估结果，并将S₁、F₁及其静态安全量化评估结果分别用于替代S₀、F₀及其静态安全量化评估结果，再根据F_ce中各个预想故障的静态安全稳定裕度估算值与潮流计算值之差，采用后验校正方法对估算参数进行修正，用于下一轮在线静态安全评估，返回步骤3)。

上述技术方案的进一步特征在于，，所述步骤2)中的t₀时刻的在线静态安全评估结果，包括与每个预想故障相对应的设备过载安全OLS薄弱支路及其裕度集和静态电压安全SVS薄弱节点及其裕度集；

所述OLS薄弱支路及其裕度集是指由与OLS主导支路裕度之差小于OLS薄弱支路筛选的裕度增量设定值的若干个OLS支路及其裕度组成的集合，用

来表示，其中，n_ol为OLS薄弱支路总数，NOL_i、η_ol.i分别为第i个OLS薄弱支路的名称和裕度，通过公式(1)计算η_ol.i；

式中，P_ol.lmt.i为S₀下预想故障的第i个OLS薄弱支路按电压和功率因数不变计算出的与其过载限值相对应的有功限值，P’_i为S₀下预想故障后其第i个OLS薄弱支路的稳态有功；

所述OLS主导支路是指OLS裕度最小的支路；

所述SVS薄弱节点及其裕度集是指由与SVS主导节点裕度之差小于SVS薄弱节点筛选的裕度增量设定值的若干个SVS节点及其裕度组成的集合，用

来表示，其中，n_sv为SVS薄弱节点总数，NSV_i、η_sv.i分别为第i个SVS薄弱节点的名称和裕度，通过公式(2)计算η_sv.i；

式中，V_sv.ulmt.i、V_sv.llmt.i分别为S₀下预想故障的第i个SVS薄弱节点的电压上、下限值，V’_i为S₀下预想故障后其第i个SVS薄弱节点的稳态电压；

所述SVS主导节点是指SVS裕度最小的节点。

上述技术方案的进一步特征在于，所述步骤5)中，针对F_e中的各个预想故障，按以下步骤进行处理：

5-1)记预想故障为f，首先将S₀下预想故障f的静态安全评估结果作为S₁下f的静态安全评估结果的初值；然后，根据S₀、S₁两个运行状态下设备的投/退变化，分别针对S₁下f的静态安全评估结果初值中各个OLS薄弱支路，从中剔除S₀下投运、S₁下退出的支路，针对S₁下f的静态安全评估结果初值中各个SVS薄弱节点，从中剔除S₀下投运、S₁下退出的节点；最后，分别通过公式(3)和公式(4)对S₁下f的静态安全评估结果初值中各个OLS薄弱支路的裕度和各个SVS薄弱节点的裕度进行估算；

式中，n_ol.f为S₁下f的OLS薄弱支路数，η_ol.f.i、η_ol.f.i.t分别为S₁下f的第i个OLS薄弱支路的裕度初值和估算值，J_gld为S₀、S₁下接入发电机、负荷和直流系统的节点总数，P_j、P_j.t分别为S₀和S₁下接入其中第j个节点的发电机、负荷和直流系统的有功注入之和，s_ol.i.j为S₀下第j个节点有功注入对S₀下f的第i个OLS薄弱支路有功的灵敏度，P_ol.lmt.i为S₀下f的第i个OLS薄弱支路按电压和功率因数不变计算出的与其过载限值相对应的有功限值，a为反映节点有功注入变化对OLS薄弱支路裕度的影响参数，初值取为1，L为S₀下投运、S₁下退出的线路和变压器集，若L中第l个线路或变压器是F₀中某个预想故障的设备且S₁下f的第i个OLS薄弱支路是S₀下该线路或变压器故障的OLS薄弱支路，则η_ol.l.i为S₀下L中第l个线路或变压器故障的OLS薄弱支路集中与S₁下f的第i个OLS薄弱支路相同的薄弱支路的裕度，P_l.i为该薄弱支路在S₀下的有功，否则，η_ol.l.i取为0，P_l.i取为P_ol.lmt.i；

n_sv.f为S₁下f的SVS薄弱节点数，η_sv.f.i、η_sv.f.i.t分别为S₁下f的第i个SVS薄弱节点的裕度初值和估算值，J_glx为S₀、S₁下接入发电机、负荷和无功补偿设备的节点总数，Q_k、Q_k.t分别为S₀和S₁下接入其中第k个节点的发电机、负荷和无功补偿设备的无功注入之和，s_sv.i.k为S₀下第k个节点无功注入对第i个SVS薄弱节点电压的灵敏度，V_sv.ulmt.i、V_sv.llmt.i分别为S₀下f的第i个SVS薄弱节点的电压上、下限值，b为反映节点无功注入变化对SVS薄弱节点裕度的影响参数，初值取为1，若L中第l个线路或变压器是F₀中某个预想故障的设备且S₁下f的第i个SVS薄弱节点是S₀下该线路或变压器故障的SVS薄弱节点，则η_sv.l.i为S₀下L中第l个线路或变压器故障的SVS薄弱节点集中与S₁下f的第i个SVS薄弱节点相同的薄弱节点的裕度，V_l.i为该薄弱节点在S₀下的电压，否则，η_sv.l.i取为0，V_l.i取为V_sv.llmt.i；

5-2)若S₁下f的所有OLS薄弱支路的裕度估算值和所有SVS薄弱节点的裕度估算值都分别大于设定的预想故障集自动筛选的OLS裕度门槛值和SVS裕度门槛值，则首先采用这些裕度估算值来替代S₁下f的静态安全评估结果中相应的裕度初值，然后将经过上述处理后的f在S₁下的静态安全评估结果初值直接作为其在S₁下的静态安全评估结果，不再对该预想故障进行基于潮流计算的静态安全量化评估，否则，将该预想故障添加到F_c中。

上述技术方案的进一步特征在于，所述步骤6)中的静态安全量化评估结果，包括与每个预想故障相对应的OLS薄弱支路及其裕度集和SVS薄弱节点及其裕度集。

上述技术方案的进一步特征在于，所述步骤7)中，通过以下步骤完成对估算参数的修正：

若S₁相对于S₀没有新投入线路或变压器，则分别通过求解公式(5)和公式(6)表示的优化函数，更新a和b，否则，不更新a和b；

式中，J_f为F_ce中预想故障数，η’_ol.d.j为步骤6)得到的S₁下F_ce中第j个预想故障的OLS主导支路的裕度，i_ol为步骤6)得到的S₁下F_ce中第j个预想故障的OLS主导支路在步骤5)得到的该预想故障OLS薄弱支路中的序号，

采用步骤5)中S₁下F_ce中第j个预想故障的第i_ol个OLS薄弱支路的裕度估算公式(3)等号的右端项来计算，a为优化变量；

η’_sv.d.j为步骤6)得到的S₁下F_ce中第j个预想故障的SVS主导节点的裕度，i_sv为步骤6)得到的S₁下F_ce中第j个预想故障的SVS主导节点在步骤5)得到的该预想故障SVS薄弱节点中的序号，

采用步骤5)中S₁下F_ce中第j个预想故障的第i_sv个SVS薄弱负荷的裕度估算公式(4)等号的右端项来计算，b为优化变量。

通过采用上述技术方案，本发明取得了下述技术效果：对于上万个节点的大电网，其N-1静态安全评估的预想故障集中故障数达到数万个，若不采用“先估算、再筛选”的计算策略，对所有预想故障分别进行基于潮流计算的静态安全评估，计算量很大，在目前调控中心的集群计算软硬配置条件下还不能满足电网运行控制对在线静态安全评估速度的要求。而本发明基于前一时间断面预想故障下静态安全的评估结果，以及当前运行状态下节点有功、无功注入相对于前一次静态安全评估所对应运行状态的变化量，直接估算当前运行状态下静态安全评估结果，只针对静态安全裕度估算值较小的预想故障进行基于潮流计算的静态安全评估，并依据基于潮流计算的静态安全评估结果与估算值的差异，采用后验校正法对估算参数进行动态修正，作为估算精度的保障措施。本发明既克服了现有技术中所述的第一种方案中预想故障集需要人工设置的不足，又克服了第二种技术方案中需要通过迭代才能确定严重故障、计算速度难以保障，以及非严重故障没有评估结果导致只针对严重故障的预防控制存在安全隐患的不足，满足了大电网静态安全评估精度和快速性的综合要求。采用本发明方法，通常情况下需要基于潮流计算的静态安全评估预想故障数不超过100个，基于目前调控中心的集群计算软硬件平台，只需要一次潮流计算的时间就能完成在线静态安全评估，可以满足电网运行控制对在线静态安全评估速度的要求。

附图说明

图1为本发明方法的流程图。

具体实施方式

下面参照附图对本发明作进一步详细描述。

图1中步骤1：设基于预想故障集自动筛选的在线静态安全评估首次启动的电力系统运行时刻为t₀，t₀时刻电力系统运行状态为S₀，基于S₀，根据《电力系统安全稳定导则》中第一道安全防线静态安全评估预想故障集的设置要求，生成预想故障集，记为F₀，进入步骤2)；

图1中步骤2：基于集群计算平台，针对S₀下F₀中所有预想故障分别进行基于潮流计算的静态安全量化评估，得到与每个预想故障相对应的设备过载安全OLS薄弱支路及其裕度集、静态电压安全SVS薄弱节点及其裕度集，作为t₀时刻的在线静态安全评估结果，进入步骤3)；

所述OLS薄弱支路及其裕度集是指由与OLS主导支路裕度之差小于OLS薄弱支路筛选的裕度增量设定值(通常设置为0.15)的若干个OLS支路及其裕度组成的集合，用

来表示，其中，n_ol为OLS薄弱支路总数，NOL_i、η_ol.i分别为第i个OLS薄弱支路的名称和裕度，通过公式(1)计算η_ol.i；

式中，P_ol.lmt.i为S₀下预想故障的第i个OLS薄弱支路按电压和功率因数不变计算出的与其过载限值相对应的有功限值，P’_i为S₀下预想故障后其第i个OLS薄弱支路的稳态有功；

所述OLS主导支路是指OLS裕度最小的支路；

所述SVS薄弱节点及其裕度集是指由与SVS主导节点裕度之差小于SVS薄弱节点筛选的裕度增量设定值(通常设置为0.15)的若干个SVS节点及其裕度组成的集合，用

来表示，其中，n_sv为SVS薄弱节点总数，NSV_i、η_sv.i分别为第i个SVS薄弱节点的名称和裕度，通过公式(2)计算η_sv.i；

式中，V_sv.ulmt.i、V_sv.llmt.i分别为S₀下预想故障的第i个SVS薄弱节点的电压上、下限值，V’_i为S₀下预想故障后其第i个SVS薄弱节点的稳态电压；

所述SVS主导节点是指SVS裕度最小的节点；

图1中步骤3：设电力系统当前运行时刻为t₁，t₁时刻电力系统运行状态为S₁，基于S₁，根据《电力系统安全稳定导则》中第一道安全防线静态安全评估预想故障集的设置要求，生成预想故障集，记为F₁，进入步骤4)；

图1中步骤4：先将F₁中不属于F₀的预想故障以及F₁中引发S₁下电力系统解列的预想故障作为需要进行基于潮流计算的静态安全量化评估的预想故障集，记为F_c，再将F₁中不属于F_c的预想故障作为需要进行基于估算的静态安全量化评估的预想故障集，记为F_e，进入步骤5)；

图1中步骤5：基于集群计算平台，针对F_e中的各个预想故障，分别按以下步骤处理后将添加到F_c中的预想故障组成的集合记为F_ce，进入步骤6)；

5-1)记预想故障为f，首先将S₀下预想故障f的静态安全评估结果作为S₁下f的静态安全评估结果的初值；然后，根据S₀、S₁两个运行状态下设备的投/退变化，分别针对S₁下f的静态安全评估结果初值中各个OLS薄弱支路，从中剔除S₀下投运、S₁下退出的支路，针对S₁下f的静态安全评估结果初值中各个SVS薄弱节点，从中剔除S₀下投运、S₁下退出的节点；最后，分别通过公式(3)和公式(4)对S₁下f的静态安全评估结果初值中各个OLS薄弱支路的裕度和各个SVS薄弱节点的裕度进行估算；

式中，n_ol.f为S₁下f的OLS薄弱支路数，η_ol.f.i、η_ol.f.i.t分别为S₁下f的第i个OLS薄弱支路的裕度初值和估算值，J_gld为S₀、S₁下接入发电机、负荷和直流系统的节点总数，P_j、P_j.t分别为S₀和S₁下接入其中第j个节点的发电机、负荷和直流系统的有功注入之和，s_ol.i.j为S₀下第j个节点有功注入对S₀下f的第i个OLS薄弱支路有功的灵敏度，P_ol.lmt.i为S₀下f的第i个OLS薄弱支路按电压和功率因数不变计算出的与其过载限值相对应的有功限值，a为反映节点有功注入变化对OLS薄弱支路裕度的影响参数，初值取为1，L为S₀下投运、S₁下退出的线路和变压器集，若L中第l个线路或变压器是F₀中某个预想故障的设备且S₁下f的第i个OLS薄弱支路是S₀下该线路或变压器故障的OLS薄弱支路，则η_ol.l.i为S₀下L中第l个线路或变压器故障的OLS薄弱支路集中与S₁下f的第i个OLS薄弱支路相同的薄弱支路的裕度，P_l.i为该薄弱支路在S₀下的有功，否则，η_ol.l.i取为0，P_l.i取为P_ol.lmt.i；

n_sv.f为S₁下f的SVS薄弱节点数，η_sv.f.i、η_sv.f.i.t分别为S₁下f的第i个SVS薄弱节点的裕度初值和估算值，J_glx为S₀、S₁下接入发电机、负荷和无功补偿设备的节点总数，Q_k、Q_k.t分别为S₀和S₁下接入其中第k个节点的发电机、负荷和无功补偿设备的无功注入之和，s_sv.i.k为S₀下第k个节点无功注入对第i个SVS薄弱节点电压的灵敏度，V_sv.ulmt.i、V_sv.llmt.i分别为S₀下f的第i个SVS薄弱节点的电压上、下限值，b为反映节点无功注入变化对SVS薄弱节点裕度的影响参数，初值取为1，若L中第l个线路或变压器是F₀中某个预想故障的设备且S₁下f的第i个SVS薄弱节点是S₀下该线路或变压器故障的SVS薄弱节点，则η_sv.l.i为S₀下L中第l个线路或变压器故障的SVS薄弱节点集中与S₁下f的第i个SVS薄弱节点相同的薄弱节点的裕度，V_l.i为该薄弱节点在S₀下的电压，否则，η_sv.l.i取为0，V_l.i取为V_sv.llmt.i；

5-2)若S₁下f的所有OLS薄弱支路的裕度估算值和所有SVS薄弱节点的裕度估算值都分别大于设定的预想故障集自动筛选的OLS裕度门槛值和SVS裕度门槛值(这两个门槛值通常都设置为0.2)，则首先采用这些裕度估算值来替代S₁下f的静态安全评估结果中相应的裕度初值，然后将经过上述处理后的f在S₁下的静态安全评估结果初值直接作为其在S₁下的静态安全评估结果，不再对该预想故障进行基于潮流计算的静态安全量化评估，否则，将该预想故障添加到F_c中；

图1中步骤6：基于集群计算平台，针对S₁下F_c中所有预想故障分别进行基于潮流计算的静态安全量化评估，得到与每个预想故障相对应的OLS薄弱支路及其裕度集和SVS薄弱节点及其裕度集，进入步骤7)；

图1中步骤7：将步骤5)得到的S₁下F₁中不属于F_c的静态安全量化评估结果和步骤6)得到的S₁下F_c的静态安全量化评估结果，作为t₁时刻的在线静态安全评估结果，并将S₁、F₁及其静态安全量化评估结果分别用于替代S₀、F₀及其静态安全量化评估结果，再根据F_ce中各个预想故障的静态安全稳定裕度估算值与潮流计算值之差，采用以下后验校正方法对估算参数进行修正，用于下一轮在线静态安全评估，返回步骤3)。

若S₁相对于S₀没有新投入线路或变压器，则分别通过求解公式(5)和公式(6)表示的优化函数，更新a和b，否则，不更新a和b；

式中，J_f为F_ce中预想故障数，η’_ol.d.j为步骤6)得到的S₁下F_ce中第j个预想故障的OLS主导支路的裕度，i_ol为步骤6)得到的S₁下F_ce中第j个预想故障的OLS主导支路在步骤5)得到的该预想故障OLS薄弱支路中的序号，

采用步骤5)中S₁下F_ce中第j个预想故障的第i_ol个OLS薄弱支路的裕度估算公式(3)等号的右端项来计算，a为优化变量；

η’_sv.d.j为步骤6)得到的S₁下F_ce中第j个预想故障的SVS主导节点的裕度，i_sv为步骤6)得到的S₁下F_ce中第j个预想故障的SVS主导节点在步骤5)得到的该预想故障SVS薄弱节点中的序号，

采用步骤5)中S₁下F_ce中第j个预想故障的第i_sv个SVS薄弱负荷的裕度估算公式(4)等号的右端项来计算，b为优化变量。

虽然本发明已以较佳实施例公开如上，但实施例并不是用来限定本发明的。在不脱离本发明之精神和范围内，所做的任何等效变化或润饰，同样属于本发明之保护范围。因此本发明的保护范围应当以本申请的权利要求所界定的内容为标准。

Claims (5)

1.基于预想故障集自动筛选的在线静态安全评估方法，其特征在于，包括以下步骤：

1)设基于预想故障集自动筛选的在线静态安全评估首次启动的电力系统运行时刻为t₀，t₀时刻电力系统运行状态为S₀，基于S₀，根据《电力系统安全稳定导则》中第一道安全防线静态安全评估预想故障集的设置要求，生成预想故障集，记为F₀，进入步骤2)；

2)针对S₀下F₀中所有预想故障分别进行基于潮流计算的静态安全量化评估，得到t₀时刻的在线静态安全评估结果，进入步骤3)；

3)设电力系统当前运行时刻为t₁，t₁时刻电力系统运行状态为S₁，基于S₁，根据《电力系统安全稳定导则》中第一道安全防线静态安全评估预想故障集的设置要求，生成预想故障集，记为F₁，进入步骤4)；

4)先将F₁中不属于F₀的预想故障以及F₁中引发S₁下电力系统解列的预想故障作为需要进行基于潮流计算的静态安全量化评估的预想故障集，记为F_c，再将F₁中不属于F_c的预想故障作为需要进行基于估算的静态安全量化评估的预想故障集，记为F_e，进入步骤5)；

5)分别针对F_e中的各个预想故障，基于其在S₀下的静态安全稳定量化评估结果以及S₀与S₁之间的差异，估算其在S₁下的静态安全稳定量化评估结果，将估算结果直接作为裕度估算值大于筛选门槛的预想故障在S₁下的静态安全稳定量化评估结果，将裕度估算值小于等于筛选门槛的预想故障添加到F_c中，并将添加到F_c中的预想故障组成的集合记为F_ce，进入步骤6)；

6)针对S₁下F_c中所有预想故障分别进行基于潮流计算的静态安全量化评估，得到静态安全量化评估结果，进入步骤7)；

7)将步骤5)得到的S₁下F₁中不属于F_c的静态安全量化评估结果和步骤6)得到的S₁下F_c的静态安全量化评估结果，作为t₁时刻的在线静态安全评估结果，并将S₁、F₁及其静态安全量化评估结果分别用于替代S₀、F₀及其静态安全量化评估结果，再根据F_ce中各个预想故障的静态安全稳定裕度估算值与潮流计算值之差，采用后验校正方法对估算参数进行修正，用于下一轮在线静态安全评估，返回步骤3)。

2.根据权利要求1所述的基于预想故障集自动筛选的在线静态安全评估方法，其特征在于，所述步骤2)中的t₀时刻的在线静态安全评估结果，包括与每个预想故障相对应的设备过载安全OLS薄弱支路及其裕度集和静态电压安全SVS薄弱节点及其裕度集；

所述OLS薄弱支路及其裕度集是指由与OLS主导支路裕度之差小于OLS薄弱支路筛选的裕度增量设定值的若干个OLS支路及其裕度组成的集合，用

来表示，其中，n_ol为OLS薄弱支路总数，NOL_i、η_ol.i分别为第i个OLS薄弱支路的名称和裕度，通过公式(1)计算η_ol.i；

式中，P_ol.lmt.i为S₀下预想故障的第i个OLS薄弱支路按电压和功率因数不变计算出的与其过载限值相对应的有功限值，P’_i为S₀下预想故障后其第i个OLS薄弱支路的稳态有功；

所述OLS主导支路是指OLS裕度最小的支路；

所述SVS薄弱节点及其裕度集是指由与SVS主导节点裕度之差小于SVS薄弱节点筛选的裕度增量设定值的若干个SVS节点及其裕度组成的集合，用

来表示，其中，n_sv为SVS薄弱节点总数，NSV_i、η_sv.i分别为第i个SVS薄弱节点的名称和裕度，通过公式(2)计算η_sv.i；

式中，V_sv.ulmt.i、V_sv.llmt.i分别为S₀下预想故障的第i个SVS薄弱节点的电压上、下限值，V’_i为S₀下预想故障后其第i个SVS薄弱节点的稳态电压；

所述SVS主导节点是指SVS裕度最小的节点。

3.根据权利要求2所述的基于预想故障集自动筛选的在线静态安全评估方法，其特征在于，所述步骤5)中，针对F_e中的各个预想故障，按以下步骤进行处理：

5-1)记预想故障为f，首先将S₀下预想故障f的静态安全评估结果作为S₁下f的静态安全评估结果的初值；然后，根据S₀、S₁两个运行状态下设备的投/退变化，分别针对S₁下f的静态安全评估结果初值中各个OLS薄弱支路，从中剔除S₀下投运、S₁下退出的支路，针对S₁下f的静态安全评估结果初值中各个SVS薄弱节点，从中剔除S₀下投运、S₁下退出的节点；最后，分别通过公式(3)和公式(4)对S₁下f的静态安全评估结果初值中各个OLS薄弱支路的裕度和各个SVS薄弱节点的裕度进行估算；

式中，n_ol.f为S₁下f的OLS薄弱支路数，η_ol.f.i、η_ol.f.i.t分别为S₁下f的第i个OLS薄弱支路的裕度初值和估算值，J_gld为S₀、S₁下接入发电机、负荷和直流系统的节点总数，P_j、P_j.t分别为S₀和S₁下接入其中第j个节点的发电机、负荷和直流系统的有功注入之和，s_ol.i.j为S₀下第j个节点有功注入对S₀下f的第i个OLS薄弱支路有功的灵敏度，P_ol.lmt.i为S₀下f的第i个OLS薄弱支路按电压和功率因数不变计算出的与其过载限值相对应的有功限值，a为反映节点有功注入变化对OLS薄弱支路裕度的影响参数，初值取为1，L为S₀下投运、S₁下退出的线路和变压器集，若L中第l个线路或变压器是F₀中某个预想故障的设备且S₁下f的第i个OLS薄弱支路是S₀下该线路或变压器故障的OLS薄弱支路，则η_ol.l.i为S₀下L中第l个线路或变压器故障的OLS薄弱支路集中与S₁下f的第i个OLS薄弱支路相同的薄弱支路的裕度，P_l.i为该薄弱支路在S₀下的有功，否则，η_ol.l.i取为0，P_l.i取为P_ol.lmt.i；

n_sv.f为S₁下f的SVS薄弱节点数，η_sv.f.i、η_sv.f.i.t分别为S₁下f的第i个SVS薄弱节点的裕度初值和估算值，J_glx为S₀、S₁下接入发电机、负荷和无功补偿设备的节点总数，Q_k、Q_k.t分别为S₀和S₁下接入其中第k个节点的发电机、负荷和无功补偿设备的无功注入之和，s_sv.i.k为S₀下第k个节点无功注入对第i个SVS薄弱节点电压的灵敏度，V_sv.ulmt.i、V_sv.llmt.i分别为S₀下f的第i个SVS薄弱节点的电压上、下限值，b为反映节点无功注入变化对SVS薄弱节点裕度的影响参数，初值取为1，若L中第l个线路或变压器是F₀中某个预想故障的设备且S₁下f的第i个SVS薄弱节点是S₀下该线路或变压器故障的SVS薄弱节点，则η_sv.l.i为S₀下L中第l个线路或变压器故障的SVS薄弱节点集中与S₁下f的第i个SVS薄弱节点相同的薄弱节点的裕度，V_l.i为该薄弱节点在S₀下的电压，否则，η_sv.l.i取为0，V_l.i取为V_sv.llmt.i；

5-2)若S₁下f的所有OLS薄弱支路的裕度估算值和所有SVS薄弱节点的裕度估算值都分别大于设定的预想故障集自动筛选的OLS裕度门槛值和SVS裕度门槛值，则首先采用这些裕度估算值来替代S₁下f的静态安全评估结果中相应的裕度初值，然后将经过上述处理后的f在S₁下的静态安全评估结果初值直接作为其在S₁下的静态安全评估结果，不再对该预想故障进行基于潮流计算的静态安全量化评估，否则，将该预想故障添加到F_c中。

4.根据权利要求3所述的基于预想故障集自动筛选的在线静态安全评估方法，其特征在于，所述步骤6)中的静态安全量化评估结果，包括与每个预想故障相对应的OLS薄弱支路及其裕度集和SVS薄弱节点及其裕度集。

5.根据权利要求4所述的基于预想故障集自动筛选的在线静态安全评估方法，其特征在于，所述步骤7)中，通过以下步骤完成对估算参数的修正：

若S₁相对于S₀没有新投入线路或变压器，则分别通过求解公式(5)和公式(6)表示的优化函数，更新a和b，否则，不更新a和b；

式中，J_f为F_ce中预想故障数，η'_ol.d.j为步骤6)得到的S₁下F_ce中第j个预想故障的OLS主导支路的裕度，i_ol为步骤6)得到的S₁下F_ce中第j个预想故障的OLS主导支路在步骤5)得到的该预想故障OLS薄弱支路中的序号，

采用步骤5)中S₁下F_ce中第j个预想故障的第i_ol个OLS薄弱支路的裕度估算公式(3)等号的右端项来计算，a为优化变量；

η'_sv.d.j为步骤6)得到的S₁下F_ce中第j个预想故障的SVS主导节点的裕度，i_sv为步骤6)得到的S₁下F_ce中第j个预想故障的SVS主导节点在步骤5)得到的该预想故障SVS薄弱节点中的序号，

采用步骤5)中S₁下F_ce中第j个预想故障的第i_sv个SVS薄弱负荷的裕度估算公式(4)等号的右端项来计算，b为优化变量。

Images