CN106099906A - 一种电力***安全稳定在线综合防御方法 - Google Patents

Abstract

本发明公开了一种电力***安全稳定综合防御方法，属于电力***安全稳定控制技术领域。本发明通过在电力***处于准稳态运行时分别计算针对预想故障集的安全稳定预防控制在线决策、针对各个预想故障的安全稳定紧急控制在线预决策和仅针对电力***实际运行状态的安全稳定紧急控制实时决策，采用这3类控制决策所针对的电力***运行时刻与当前时刻之间的时长是否超过相应的定值作为各类控制决策有效性的判据，结合安全稳定紧急控制实时决策优先、安全稳定紧急控制在线预决策次之、安全稳定预防控制在线决策最后的控制决策选择原则来确定用于执行的控制措施，为电力***安全稳定的在线综合防御提供了一种实用有效的方法。

Description

一种电力***安全稳定在线综合防御方法

技术领域

本发明属于电力***安全稳定控制技术领域，更准确地说，本发明涉及一种适用于综合应用在线预防控制、在线紧急控制和实时紧急控制多种控制手段防御电力***失去安全稳定性的方法。

背景技术

目前电力***安全稳定综合防御手段采用的主要是依据电力***安全稳定导则通过离线分析计算规划运行方式，辅助采用基于电力***实际运行断面数据的在线安全稳定分析技术来支撑电力***的安全稳定运行。随着间隙性新能源发电的大量接入、电力市场的逐步推进和未来能源互联网的实施，电力***的实际运行方式与常规离线策划的运行方式的偏差会越来越大。电力***的安全稳定综合防御将逐渐转变为离线以研究电力***运行特性和机理为主、安排运行方式为辅，在线以调度优化控制运行方式为主的电力***调度运行控制模式。

确保电力***安全稳定运行就是要能够应对电力***可能遭受的各种故障(或扰动)，实际上，电力***可能遭受的故障(或扰动)是无法穷尽的，受制于目前计算能力的限制，在线安全稳定分析也只能针对少量预想故障进行分析，因此，必然会有一些其它没有预想到的故障发生后会导致电力***失去安全稳定。除此之外，由于在线安全稳定分析总需要一定的计算时间，而在这段时间内电力***很可能已经发生了变化，即使采用基于未来运行方式的在线安全稳定分析，也会由于未来运行方式很可能与未来实际的运行方式的不一致，导致针对预想故障的在线安全稳定分析也难以保障这些预想故障真的发生后的安全稳定性。实际上，还有由于电力***模型和参数，特别是暂态和动态安全稳定计算紧密相关的模型和参数与实际的不一致，也会造成在线安全稳定分析结果与电力***实际的偏差。总之，不仅需要改进、完善现有的在线安全稳定分析技术，而且要设计综合应用针对预想故障集的安全稳定预防控制在线决策、针对各个预想故障的安全稳定紧急控制在线预决策和仅针对电力***实际运行状态的安全稳定紧急控制实时决策的防御方法，以提升电力***安全稳定综合防御的有效性和实用性。

发明内容

本发明目的是：弥补现有在线安全稳定分析中预防控制在线决策、紧急控制在线预决策与实时紧急控制决策三者彼此相互独立、没有统筹设计的不足，提出一种能够综合应用在线预防控制、在线紧急控制和实时紧急控制多种控制手段防御电力***失去安全稳定性的方法。

本发明的基本原理在于：对于各类对实时性要求高(秒级)的安全稳定问题只能通过在预想故障发生前对电力***采取预防控制或在预想故障发生后依靠安全自动装置对电力***采取紧急控制，对于各类对实时性要求不太高(分钟级)的安全稳定问题则可以通过及时有效的调度运行控制手段对电力***实施控制。因此，通过在线预防控制来防御包括安全自动装置防御的预想故障在内的尽可能多的预想故障发生后的各类对实时性要求高(秒级)的安全稳定问题，通过针对预想故障的在线紧急控制来防御预想故障后电力***过渡到准稳态时存在的各类对实时性要求不太高(分钟级)的安全稳定问题，通过针对实际运行状态的实时紧急控制来防御实际运行状态中存在的各类对实时性要求不太高(分钟级)的安全稳定问题，并采用这3类控制决策所针对的电力***运行时刻与当前时刻之间的时长是否超过相应的定值作为各类控制决策有效性的判据，结合安全稳定紧急控制实时决策优先、安全稳定紧急控制在线预决策次之、安全稳定预防控制在线决策最后的控制决策选择原则来确定用于执行的控制措施，可以有针对性地利用不同控制手段的优点，弥补不同控制手段的局限性，实现电力***安全稳定的在线综合防御。

具体地说，本发明具体采用以下技术方案：

一种电力***安全稳定在线综合防御方法包括电力***准稳态运行状态识别、电力***准稳态运行断面数据生成、安全稳定预防控制在线决策计算、安全稳定紧急控制在线预决策计算、安全稳定紧急控制实时决策计算和控制决策选择，其中：

所述电力***准稳态运行状态识别是指若在设定时长Δt内电力***中各母线频率的波动幅度都小于设定值Δf，则将当前运行状态识别为电力***处于准稳态运行状态，否则，将当前运行状态识别为电力***不处于准稳态运行状态；

所述电力***准稳态运行断面数据生成是指当电力***准稳态运行状态识别的结果为电力***处于准稳态运行状态时，首先，将当前时刻作为电力***准稳态运行断面数据所对应的时刻t_r，然后，对电力***当前运行状态进行网络拓扑分析，再针对各个计算节点，分别将各个计算节点在(t_r-Δt)时刻至t_r时刻之间的频率平均值作为该计算节点的稳态频率，并将所有计算节点的稳态频率的平均值作为电力***的稳态频率f_r，针对各个注入节点，分别将各个注入节点有功/无功在(t_r-Δt)时刻至t_r时刻之间的平均值作为该注入节点的稳态有功/无功，针对各个节点，分别将各个节点电压在(t_r-Δt)时刻至t_r时刻之间的平均值作为该节点的稳态电压，最后，根据网络拓扑及模型与参数、f_r、各个注入节点的稳态有功/无功和各节点的稳态电压生成电力***准稳态运行断面数据，并将包括等值发电机节点在内的发电机节点的发电机稳态有功之和作为电力***准稳态运行断面数据所对应的发电机有功出力之和P_G.sum.r；

所述注入节点包括发电机节点、等值发电机节点、负荷节点、等值负荷节点和无功设备接入节点；

所述安全稳定预防控制在线决策计算是指当电力***处于准稳态运行状态且最新生成的电力***准稳态运行断面数据所对应的t_r与当前时刻之间的时长小于Δt时，记所有需要防御的预想故障集为F，针对F，基于最新生成的电力***准稳态运行断面数据，考虑保护、安自装置和AGC/AVC的控制模型及参数，以满足F中各个预想故障发生后电力***都能保持暂态安全稳定、动态安全稳定和静态安全稳定进行优化决策计算，若在安全稳定预防控制在线决策计算过程中电力***不处于准稳态运行状态，则中止安全稳定预防控制在线决策计算；

所述安全稳定紧急控制在线预决策计算是指当电力***处于准稳态运行状态且最新生成的电力***准稳态运行断面数据所对应的t_r与当前时刻之间的时长小于Δt时，针对F中各个预想故障，分别基于最新生成的电力***准稳态运行断面数据，以满足相应的预想故障下暂态和动态过程结束后过渡到准稳态时电力***都能保持输电断面有功越限安全、设备过载安全、母线电压越限安全和电网频率越限安全进行优化决策计算，若在安全稳定紧急控制在线预决策计算过程中电力***不处于准稳态运行状态，则中止安全稳定紧急控制在线预决策计算；

所述安全稳定紧急控制实时决策计算是指当电力***处于准稳态运行状态且最新生成的电力***准稳态运行断面数据所对应的t_r与当前时刻之间的时长小于Δt时，基于最新生成的电力***准稳态运行断面数据，以满足电力***都能保持输电断面有功越限安全、设备过载安全、母线电压越限安全和电网频率越限安全进行优化决策计算，若在安全稳定紧急控制实时决策计算过程中电力***不处于准稳态运行状态，则中止安全稳定紧急控制实时决策计算；

所述控制决策选择包括以下步骤：

1)若安全稳定紧急控制实时决策计算的最新决策结果中有控制措施且电力***当前运行时刻与安全稳定紧急控制实时决策计算的最新决策结果所基于的电力***准稳态运行断面数据所对应的t_r之间的时长小于设定值Δt₁，则将安全稳定紧急控制实时决策计算的最新决策结果作为用于执行的控制决策，结束本次控制决策的选择，否则，进入步骤2)；

2)若F中某个预想故障的发生时刻t_f与安全稳定紧急控制在线预决策计算中该预想故障的最新决策结果所基于的电力***准稳态运行断面数据所对应的t_r之间的时长小于设定值Δt₂且t_r早于t_f，则将安全稳定紧急控制在线预决策计算中该预想故障的最新决策结果作为用于执行的控制决策，结束本次控制决策的选择，否则，进入步骤3)；

3)若电力***当前运行时刻与安全稳定预防控制在线决策计算的最新决策结果所基于的电力***准稳态运行断面数据所对应的t_r之间的时长小于设定值Δt₃，则将安全稳定预防控制在线决策计算的最新决策结果作为用于执行的控制决策，结束本次控制决策的选择，否则，认为不需要对电力***采取控制措施，结束本次控制决策的选择；

所述电力***准稳态运行状态识别、电力***准稳态运行断面数据生成、安全稳定预防控制在线决策计算、安全稳定紧急控制在线预决策计算、安全稳定紧急控制实时决策计算和控制决策选择在电力***的运行过程中分别按循环执行的模式独立运行。进一步，所述安全稳定紧急控制在线预决策计算，针对F中各个预想故障，分别进行如下步骤的计算：

A-1)基于最新生成的电力***准稳态运行断面数据，针对预想故障进行考虑保护、安自装置和AGC/AVC控制模型及参数的时域仿真，直至仿真到各个计算节点频率动态过程中相邻的极大值与极小值之差都小于设定值Δf时，终止时域仿真，得到各个计算节点的频率动态过程、各个注入节点的有功/无功动态过程、各个输电设备的有功动态过程和各个节点的电压动态过程，进入步骤A-2)；

所述输电设备包括线路和变压器；

A-2)针对各个计算节点，分别将计算各个节点频率动态过程结束前的极大值与极小值二者的平均值作为考虑预想故障后该计算节点的稳态频率，再将所有计算节点的稳态频率的平均值作为考虑预想故障后电力***的稳态频率f_e；针对各个注入节点，分别将各个注入节点有功/无功动态过程结束前的极大值与极小值二者的平均值作为考虑预想故障后该注入节点的稳态有功/无功，并将包括等值发电机节点在内的发电机节点的发电机稳态有功之和记为P_G.sum.e；针对各个输电设备，分别将各个输电设备有功动态过程结束前的极大值与极小值二者的平均值作为考虑预想故障后该输电设备的稳态有功，再针对各个用于安全稳定监视的输电通道，分别将各个输电通道中所有输电设备的稳态有功之和作为考虑预想故障后该输电通道的稳态有功；针对各个节点，分别将各个节点电压动态过程结束前的极大值与极小值二者的平均值作为考虑预想故障后该节点的稳态电压，进入步骤A-3)；

A-3)根据考虑预想故障后各个注入节点的稳态有功/无功和预想故障后电力***的网络拓扑及模型与参数，建立电力***的潮流方程，进入步骤A-4)；

A-4)通过求解公式(1)表示的优化函数，计算出针对最新生成的电力***准稳态运行断面数据所对应的t_r时刻该预想故障的安全稳定紧急控制在线预决策：

公式(1)中，I为包括等值发电机节点在内的发电机节点中预先设置的发电机有功可控节点的总数，C_G.i、ΔP_G.i分别为第i个发电机有功可控节点的发电机单位有功调整代价和发电机有功变化量，J为包括等值负荷节点在内的负荷节点中预先设置的负荷可控节点的总数，C_L.j、ΔP_L.j分别为第j个负荷可控节点的负荷单位有功调整代价和负荷有功变化量，P_G.i.0.e、P_G.i.max.e、P_G.i.min.e分别为预想故障后第i个发电机有功可控节点的发电机稳态有功、发电机稳态有功在设定的时长Δt_c内可控区间的上限和发电机稳态有功在Δt_c内可控区间的下限，P_L.i.0.e、P_L.i.min.e分别为预想故障后第j个负荷可控节点的负荷稳态有功和负荷稳态有功在Δt_c内可控区间的下限，K为电网中的输电设备数，P_k.0.e为预想故障后第k个输电设备的稳态有功，P_k.LMT为第k个输电设备的稳态有功限额，s_k.G.i.e为基于步骤A-3)潮流方程计算出的第i个发电机有功可控节点的发电机有功对第k个输电设备有功的灵敏度，s_k.L.j.e为基于步骤A-3)潮流方程计算出的第j个负荷可控节点的负荷有功对第k个输电设备有功的灵敏度，M为用于安全稳定监视的输电通道数，N_m为第m个用于安全稳定监视的输电通道中的输电设备数，P_m.0.e、P_m.LMT.e、P'_m.LMT.e分别为预想故障后第m个用于安全稳定监视的输电通道的稳态有功、正向安全稳定限额和反向安全稳定限额，s_m.n.G.i.e为基于步骤A-3)潮流方程计算出的第i个发电机有功可控节点的发电机有功对第m个用于安全稳定监视的输电通道中第n个输电设备有功的灵敏度，s_m.n.L.i.e为基于步骤A-3)潮流方程计算出的第j个负荷可控节点的负荷有功对第m个用于安全稳定监视的输电通道中第n个输电设备有功的灵敏度，P_G.sum.r、f_r分别为与最新生成的电力***准稳态运行断面数据所对应的P_G.sum.r和f_r，f_cr.low、f_cr.up分别为电网稳态频率越限安全的下限和上限，λ为考虑仿真与实际动态过程存在差异而设定的参数，X为预先设置的无功可控节点数，Q_x.0.e、ΔQ_x、Q_x.min.e、Q_x.max.e分别为预想故障后第x个无功可控节点的无功设备的稳态无功、稳态无功变化量、稳态无功最小值和稳态无功最大值，B为电网中的节点数，V_b.0.e、V_b.low.e、V_b.up.e分别为预想故障后第b个节点的稳态电压、稳态电压下限和稳态电压上限，s_vq.b.x.e为基于步骤A-3)潮流方程计算出的第x个无功可控节点的无功设备的无功对第b个节点电压的灵敏度；

所述输电设备的稳态有功限额是指根据输电设备的过载电流限额和额定电压，按指定功率因数计算出的稳态有功。

进一步，针对安全稳定紧急控制实时决策计算包括如下步骤：

B-1)根据最新生成的电力***准稳态运行断面数据，建立电力***的潮流方程，进入步骤B-2)；

2)通过求解公式(2)表示的优化函数，计算出针对最新生成的电力***准稳态运行断面数据所对应的t_r时刻电力***的安全稳定紧急控制实时决策：

公式(2)中，I为包括等值发电机节点在内的发电机节点中预先设置的发电机有功可控节点的总数，C_G.i、ΔP_G.i分别为第i个发电机有功可控节点的发电机单位有功调整代价和发电机有功变化量，J为包括等值负荷节点在内的负荷节点中预先设置的负荷可控节点的总数，C_L.j、ΔP_L.j分别为第j个负荷可控节点的负荷单位有功调整代价和负荷有功变化量，P_G.i.0、P_G.i.max、P_G.i.min分别为第i个发电机有功可控节点的发电机稳态有功、发电机稳态有功在Δt_c内可控区间的上限和发电机稳态有功在Δt_c内可控区间的下限，P_L.i.0、P_L.i.min分别为第j个负荷可控节点的负荷稳态有功和负荷稳态有功在Δt_c内可控区间的下限，K为电网中的输电设备数，P_k.0为第k个输电设备的稳态有功，P_k.LMT为第k个输电设备的稳态有功限额，s_k.G.i为基于步骤B-1)潮流方程计算出的第i个发电机有功可控节点的发电机有功对第k个输电设备有功的灵敏度，s_k.L.j为基于步骤B-1)潮流方程计算出的第j个负荷可控节点的负荷有功对第k个输电设备有功的灵敏度，M为用于安全稳定监视的输电通道数，N_m为第m个用于安全稳定监视的输电通道中的输电设备数，P_m.0、P_m.LMT、P'_m.LMT分别为第m个用于安全稳定监视的输电通道的稳态有功、正向安全稳定限额和反向安全稳定限额，s_m.n.G.i为基于步骤B-1)潮流方程计算出的第i个发电机有功可控节点的发电机有功对第m个用于安全稳定监视的输电通道中第n个输电设备有功的灵敏度，s_m.n.L.i为基于步骤B-1)潮流方程计算出的第j个负荷可控节点的负荷有功对第m个用于安全稳定监视的输电通道中第n个输电设备有功的灵敏度，P_G.sum.r、f_r分别为与最新生成的电力***准稳态运行断面数据所对应的P_G.sum.r和f_r，P_G.sum.r.0、f_r.0分别为与前一次生成的电力***准稳态运行断面数据所对应的P_G.sum.r和f_r，f_cr.low、f_cr.up分别为电网稳态频率越限安全的下限和上限，X为预先设置的无功可控节点数，Q_x.0、ΔQ_x、Q_x.min、Q_x.max分别为第x个无功可控节点的无功设备的稳态无功、稳态无功变化量、稳态无功最小值和稳态无功最大值，B为电网中的节点数，V_b.0、V_b.low、V_b.up分别为第b个节点的稳态电压、稳态电压下限和稳态电压上限，s_vq.b.x为基于步骤B-1)潮流方程计算出的第x个无功可控节点的无功设备的无功对第b个节点电压的灵敏度。

通过采用上述技术方案，本发明取得了下述技术效果：

在预防控制在线决策中考虑包括安全自动装置防御的预想故障在内的尽可能多的预想故障发生后的各类对实时性要求高(秒级)的安全稳定问题，通过在预想故障发生之前就采取控制，可以在较大程度上避免因安全自动装置的控制策略不合理而存在的各类对实时性要求高(秒级)的安全稳定隐患。通过针对预想故障的在线紧急控制来防御预想故障后电力***过渡到准稳态时存在的各类对实时性要求不太高(分钟级)的安全稳定问题，可以为预想故障发生后立即提供有针对性的控制决策支持，及时阻断安全稳定问题的扩大。通过针对实际运行状态的实时紧急控制来防御实际运行状态中存在的各类对实时性要求不太高(分钟级)的安全稳定问题，可以克服因针对预想故障的时域仿真与预想故障后电力***实际运行状态之间的偏差所带来的紧急控制在线预决策不合理，并弥补因实际发生的故障没有预先考虑到而缺失紧急控制在线预决策的不足。

采用这3类控制决策所针对的电力***运行时刻与当前时刻之间的时长是否超过相应的定值作为各类控制决策有效性的判据，结合安全稳定紧急控制实时决策优先、安全稳定紧急控制在线预决策次之、安全稳定预防控制在线决策最后的控制决策选择原则来确定用于执行的控制措施，可以有针对性地利用不同控制手段的优点，弥补不同控制手段的局限性，实现电力***安全稳定的在线综合防御。

附图说明

图1为本发明方法中针对某个预想故障的安全稳定紧急控制在线预决策计算的流程图。

图2为本发明方法中针对安全稳定紧急控制实时决策计算的流程图。

图3为本发明方法中控制决策选择的流程图。

具体实施方式

下面结合附图和实施例对本发明作进一步详细说明。

本实施例的一种电力***安全稳定在线综合防御方法，包括电力***准稳态运行状态识别、电力***准稳态运行断面数据生成、安全稳定预防控制在线决策计算、安全稳定紧急控制在线预决策计算、安全稳定紧急控制实时决策计算和控制决策选择，其中：

所述电力***准稳态运行状态识别是指若在设定时长Δt(通常设置为10s)内电力***中各母线频率的波动幅度都小于设定值Δf(通常设置为0.02Hz)，则将当前运行状态识别为电力***处于准稳态运行状态，否则，将当前运行状态识别为电力***不处于准稳态运行状态；

所述电力***准稳态运行断面数据生成是指当电力***准稳态运行状态识别结果为电力***处于准稳态运行状态时，首先，将当前时刻作为电力***准稳态运行断面数据所对应的时刻t_r，然后，对电力***当前运行状态进行网络拓扑分析，再针对各个计算节点，分别将各个计算节点在(t_r-Δt)时刻至t_r时刻之间的频率平均值作为该计算节点的稳态频率，并将所有计算节点的稳态频率的平均值作为电力***的稳态频率f_r，针对各个注入节点，分别将各个注入节点有功/无功在(t_r-Δt)时刻至t_r时刻之间的平均值作为该注入节点的稳态有功/无功，针对各个节点，分别将各个节点电压在(t_r-Δt)时刻至t_r时刻之间的平均值作为该节点的稳态电压，最后，根据网络拓扑及模型与参数、f_r、各个注入节点的稳态有功/无功和各节点的稳态电压生成电力***准稳态运行断面数据，并将包括等值发电机节点在内的发电机节点的发电机稳态有功之和作为电力***准稳态运行断面数据所对应的发电机有功出力之和P_G.sum.r；

所述注入节点包括发电机节点、等值发电机节点、负荷节点、等值负荷节点和无功设备接入节点；

所述安全稳定预防控制在线决策计算是指当电力***处于准稳态运行状态且最新生成的电力***准稳态运行断面数据所对应的t_r与当前时刻之间的时长小于Δt时，记所有需要防御的预想故障集为F，针对F，基于最新生成的电力***准稳态运行断面数据，考虑保护、安自装置和AGC/AVC的控制模型及参数，以满足F中各个预想故障发生后电力***都能保持暂态安全稳定、动态安全稳定和静态安全稳定进行优化决策计算，若在安全稳定预防控制在线决策计算过程中电力***不处于准稳态运行状态，则中止安全稳定预防控制在线决策计算；

其中，对于静态安全稳定约束可以通过设置不同安全稳定门槛值以区分对实时性要求高(秒级)静态安全稳定问题(例如严重过载)和对实时要求不太高(分钟级)的静态安全稳定问题(例如轻微过载)，安全稳定预防控制在线决策计算可以设置为只保障暂态安全稳定、动态安全稳定和对实时性要求高的静态安全稳定，也可设置为要保障暂态安全稳定、动态安全稳定和对实时性要求不太高的静态安全稳定；

所述安全稳定紧急控制在线预决策计算是指当电力***处于准稳态运行状态且最新生成的电力***准稳态运行断面数据所对应的t_r与当前时刻之间的时长小于Δt时，针对F中各个预想故障，分别基于最新生成的电力***准稳态运行断面数据，以满足相应的预想故障下暂态和动态过程结束后过渡到准稳态时电力***都能保持输电断面有功越限安全、设备过载安全、母线电压越限安全和电网频率越限安全进行优化决策计算，若在安全稳定紧急控制在线预决策计算过程中电力***不处于准稳态运行状态，则中止安全稳定紧急控制在线预决策计算；

针对F中某个预想故障的安全稳定紧急控制在线预决策计算流程如图1所示。

图1步骤1：基于最新生成的电力***准稳态运行断面数据，针对预想故障进行考虑保护、安自装置和AGC/AVC控制模型及参数的时域仿真，直至仿真到各个计算节点频率动态过程中相邻的极大值与极小值之差都小于设定值Δf时，终止时域仿真，得到各个计算节点的频率动态过程、各个注入节点的有功/无功动态过程、各个输电设备的有功动态过程和各个节点的电压动态过程，进入步骤2)；

所述输电设备包括线路和变压器；

图1步骤2：针对各个计算节点，分别将各个计算节点频率动态过程结束前的极大值与极小值二者的平均值作为考虑预想故障后该计算节点的稳态频率，再将所有计算节点的稳态频率的平均值作为考虑预想故障后电力***的稳态频率f_e；针对各个注入节点，分别将各个注入节点有功/无功动态过程结束前的极大值与极小值二者的平均值作为考虑预想故障后该注入节点的稳态有功/无功，并将包括等值发电机节点在内的发电机节点的发电机稳态有功之和记为P_G.sum.e；针对各个输电设备，分别将各个输电设备有功动态过程结束前的极大值与极小值二者的平均值作为考虑预想故障后该输电设备的稳态有功，再针对各个用于安全稳定监视的输电通道，分别将各个输电通道中所有输电设备的稳态有功之和作为考虑预想故障后该输电通道的稳态有功；针对各个节点，分别将各个节点电压动态过程结束前的极大值与极小值二者的平均值作为考虑预想故障后该节点的稳态电压，进入步骤3)；

图1步骤3：根据考虑预想故障后各个注入节点的稳态有功/无功和预想故障后电力***的网络拓扑及模型与参数，建立电力***的潮流方程，进入步骤4)；

图1步骤4：通过求解公式(1)表示的优化函数，计算出针对最新生成的电力***准稳态运行断面数据所对应的t_r时刻该预想故障的安全稳定紧急控制在线预决策。

公式(1)中，I为包括等值发电机节点在内的发电机节点中预先设置的发电机有功可控节点的总数，C_G.i、ΔP_G.i分别为第i个发电机有功可控节点的发电机单位有功调整代价和发电机有功变化量，J为包括等值负荷节点在内的负荷节点中预先设置的负荷可控节点的总数，C_L.j、ΔP_L.j分别为第j个负荷可控节点的负荷单位有功调整代价和负荷有功变化量，P_G.i.0.e、P_G.i.max.e、P_G.i.min.e分别为预想故障后第i个发电机有功可控节点的发电机稳态有功、发电机稳态有功在设定的时长Δt_c(通常设置为3分钟)内可控区间的上限和发电机稳态有功在Δt_c内可控区间的下限，P_L.i.0.e、P_L.i.min.e分别为预想故障后第j个负荷可控节点的负荷稳态有功和负荷稳态有功在Δt_c内可控区间的下限，K为电网中的输电设备数，P_k.0.e为预想故障后第k个输电设备的稳态有功，P_k.LMT为第k个输电设备的稳态有功限额，s_k.G.i.e为基于步骤3)潮流方程计算出的第i个发电机有功可控节点的发电机有功对第k个输电设备有功的灵敏度，s_k.L.j.e为基于步骤3)潮流方程计算出的第j个负荷可控节点的负荷有功对第k个输电设备有功的灵敏度，M为用于安全稳定监视的输电通道数，N_m为第m个用于安全稳定监视的输电通道中的输电设备数，P_m.0.e、P_m.LMT.e P'_m.LMT.e分别为预想故障后第m个用于安全稳定监视的输电通道的稳态有功、正向安全稳定限额和反向安全稳定限额，s_m.n.G.i.e为基于步骤3)潮流方程计算出的第i个发电机有功可控节点的发电机有功对第m个用于安全稳定监视的输电通道中第n个输电设备有功的灵敏度，s_m.n.L.i.e为基于步骤3)潮流方程计算出的第j个负荷可控节点的负荷有功对第m个用于安全稳定监视的输电通道中第n个输电设备有功的灵敏度，P_G.sum.r、f_r分别为与最新生成的电力***准稳态运行断面数据所对应的P_G.sum.r和f_r，f_cr.low、f_cr.up分别为电网稳态频率越限安全的下限和上限，λ为考虑仿真与实际动态过程存在差异而设定的参数(可以根据电网最近发生大功率缺失故障的分析计算得到，假设故障后过渡到准稳态时的电网实际频率与故障前实际频率的差值为a，针对该大功率缺失故障的仿真得到的电网准稳态频率与故障前实际频率的频差为b，则λ通常设置为a/b)，X为预先设置的无功可控节点数，Q_x.0.e、ΔQ_x、Q_x.min.e、Q_x.max.e分别为预想故障后第x个无功可控节点的无功设备的稳态无功、稳态无功变化量、稳态无功最小值和稳态无功最大值，B为电网中的节点数，V_b.0.e、V_b.low.e、V_b.up.e分别为预想故障后第b个节点的稳态电压、稳态电压下限和稳态电压上限，s_vq.b.x.e为基于步骤3)潮流方程计算出的第x个无功可控节点的无功设备的无功对第b个节点电压的灵敏度；

所述输电设备的稳态有功限额是指根据输电设备的过载电流限额和额定电压，按指定功率因数计算出的稳态有功。

所述安全稳定紧急控制实时决策计算是指当电力***处于准稳态运行状态且最新生成的电力***准稳态运行断面数据所对应的t_r与当前时刻之间的时长小于Δt时，基于最新生成的电力***准稳态运行断面数据，以满足电力***都能保持输电断面有功越限安全、设备过载安全、母线电压越限安全和电网频率越限安全而进行的优化决策计算，若在安全稳定紧急控制实时决策计算过程中电力***不处于准稳态运行状态，则中止安全稳定紧急控制实时决策计算；

针对安全稳定紧急控制实时决策计算流程如图2所示：

图2步骤1：根据最新生成的电力***准稳态运行断面数据，建立电力***的潮流方程，进入步骤2；

图2步骤2：通过求解公式(2)表示的优化函数，计算出针对最新生成的电力***准稳态运行断面数据所对应的t_r时刻电力***的安全稳定紧急控制实时决策。

公式(2)中，I为包括等值发电机节点在内的发电机节点中预先设置的发电机有功可控节点的总数，C_G.i、ΔP_G.i分别为第i个发电机有功可控节点的发电机单位有功调整代价和发电机有功变化量，J为包括等值负荷节点在内的负荷节点中预先设置的负荷可控节点的总数，C_L.j、ΔP_L.j分别为第j个负荷可控节点的负荷单位有功调整代价和负荷有功变化量，P_G.i.0、P_G.i.max、P_G.i.min分别为第i个发电机有功可控节点的发电机稳态有功、发电机稳态有功在Δt_c内可控区间的上限和发电机稳态有功在Δt_c内可控区间的下限，P_L.i.0、P_L.i.min分别为第j个负荷可控节点的负荷稳态有功和负荷稳态有功在Δt_c内可控区间的下限，K为电网中的输电设备数，P_k.0为第k个输电设备的稳态有功，P_k.LMT为第k个输电设备的稳态有功限额，s_k.G.i为基于步骤1)潮流方程计算出的第i个发电机有功可控节点的发电机有功对第k个输电设备有功的灵敏度，s_k.L.j为基于步骤1)潮流方程计算出的第j个负荷可控节点的负荷有功对第k个输电设备有功的灵敏度，M为用于安全稳定监视的输电通道数，N_m为第m个用于安全稳定监视的输电通道中的输电设备数，P_m.0、P_m.LMTP'_m.LMT分别为第m个用于安全稳定监视的输电通道的稳态有功、正向安全稳定限额和反向安全稳定限额，s_m.n.G.i为基于步骤1)潮流方程计算出的第i个发电机有功可控节点的发电机有功对第m个用于安全稳定监视的输电通道中第n个输电设备有功的灵敏度，s_m.n.L.i为基于步骤1)潮流方程计算出的第j个负荷可控节点的负荷有功对第m个用于安全稳定监视的输电通道中第n个输电设备有功的灵敏度，P_G.sum.r、f_r分别为与最新生成的电力***准稳态运行断面数据所对应的P_G.sum.r和f_r，P_G.sum.r.0、f_r.0分别为与前一次生成的电力***准稳态运行断面数据所对应的P_G.sum.r和f_r，f_cr.low、f_cr.up分别为电网稳态频率越限安全的下限和上限，X为预先设置的无功可控节点数，Q_x.0、ΔQ_x、Q_x.min、Q_x.max分别为第x个无功可控节点的无功设备的稳态无功、稳态无功变化量、稳态无功最小值和稳态无功最大值，B为电网中的节点数，V_b.0、V_b.low、V_b.up分别为第b个节点的稳态电压、稳态电压下限和稳态电压上限，s_vq.b.x为基于步骤1)潮流方程计算出的第x个无功可控节点的无功设备的无功对第b个节点电压的灵敏度。

所述控制决策选择的流程如图3所示。

图3中步骤1：若安全稳定紧急控制实时决策计算的最新决策结果中有控制措施且电力***当前运行时刻与安全稳定紧急控制实时决策计算的最新决策结果所基于的电力***准稳态运行断面数据所对应的t_r之间的时长小于设定值Δt₁(通常设置为30s)，则将安全稳定紧急控制实时决策计算的最新决策结果作为用于执行的控制决策，结束本次控制决策的选择，否则，进入步骤2；

图3中步骤2：若F中某个预想故障的发生时刻t_f与安全稳定紧急控制在线预决策计算中该预想故障的最新决策结果所基于的电力***准稳态运行断面数据所对应的t_r之间的时长小于设定值Δt₂(通常设置为2分钟)且t_r早于t_f，则将安全稳定紧急控制在线预决策计算中该预想故障的最新决策结果作为用于执行的控制决策，结束本次控制决策的选择，否则，进入步骤3；

图3中步骤3：若电力***当前运行时刻与安全稳定预防控制在线决策计算的最新决策结果所基于的电力***准稳态运行断面数据所对应的t_r之间的时长小于设定值Δt₃(通常设置为5分钟)，则将安全稳定预防控制在线决策计算的最新决策结果作为用于执行的控制决策，结束本次控制决策的选择，否则，认为不需要对电力***采取控制措施，结束本次控制决策的选择。

所述电力***准稳态运行状态识别、电力***准稳态运行断面数据生成、安全稳定预防控制在线决策计算、安全稳定紧急控制在线预决策计算、安全稳定紧急控制实时决策计算和控制决策选择在电力***的运行过程中分别按循环执行的模式独立运行。

需要进一步说明的是：对于异步联网的电力***，需要针对各个同步电网分别进行上述处理，并计及同步电网之间直流***换流站功率调整对送端同步电网和受端同步电网安全稳定的影响。

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

Claims (3)

1.一种电力***安全稳定在线综合防御方法，其特征在于，包括电力***准稳态运行状态识别、电力***准稳态运行断面数据生成、安全稳定预防控制在线决策计算、安全稳定紧急控制在线预决策计算、安全稳定紧急控制实时决策计算和控制决策选择，其中：

所述电力***准稳态运行状态识别是指若在设定时长Δt内电力***中各母线频率的波动幅度都小于设定值Δf，则将当前运行状态识别为电力***处于准稳态运行状态，否则，将当前运行状态识别为电力***不处于准稳态运行状态；

所述电力***准稳态运行断面数据生成是指当电力***准稳态运行状态识别的结果为电力***处于准稳态运行状态时，首先，将当前时刻作为电力***准稳态运行断面数据所对应的时刻t_r，然后，对电力***当前运行状态进行网络拓扑分析，再针对各个计算节点，分别将各个计算节点在(t_r-Δt)时刻至t_r时刻之间的频率平均值作为该计算节点的稳态频率，并将所有计算节点的稳态频率的平均值作为电力***的稳态频率f_r，针对各个注入节点，分别将各个注入节点有功/无功在(t_r-Δt)时刻至t_r时刻之间的平均值作为该注入节点的稳态有功/无功，针对各个节点，分别将各个节点电压在(t_r-Δt)时刻至t_r时刻之间的平均值作为该节点的稳态电压，最后，根据网络拓扑及模型与参数、f_r、各个注入节点的稳态有功/无功和各节点的稳态电压生成电力***准稳态运行断面数据，并将包括等值发电机节点在内的发电机节点的发电机稳态有功之和作为电力***准稳态运行断面数据所对应的发电机有功出力之和P_G.sum.r；

所述注入节点包括发电机节点、等值发电机节点、负荷节点、等值负荷节点和无功设备接入节点；

所述安全稳定预防控制在线决策计算是指当电力***处于准稳态运行状态且最新生成的电力***准稳态运行断面数据所对应的t_r与当前时刻之间的时长小于Δt时，记所有需要防御的预想故障集为F，针对F，基于最新生成的电力***准稳态运行断面数据，考虑保护、安自装置和AGC/AVC的控制模型及参数，以满足F中各个预想故障发生后电力***都能保持暂态安全稳定、动态安全稳定和静态安全稳定进行优化决策计算，若在安全稳定预防控制在线决策计算过程中电力***不处于准稳态运行状态，则中止安全稳定预防控制在线决策计算；

所述安全稳定紧急控制在线预决策计算是指当电力***处于准稳态运行状态且最新生成的电力***准稳态运行断面数据所对应的t_r与当前时刻之间的时长小于Δt时，针对F中各个预想故障，分别基于最新生成的电力***准稳态运行断面数据，以满足相应的预想故障下暂态和动态过程结束后过渡到准稳态时电力***都能保持输电断面有功越限安全、设备过载安全、母线电压越限安全和电网频率越限安全进行优化决策计算，若在安全稳定紧急控制在线预决策计算过程中电力***不处于准稳态运行状态，则中止安全稳定紧急控制在线预决策计算；

所述安全稳定紧急控制实时决策计算是指当电力***处于准稳态运行状态且最新生成的电力***准稳态运行断面数据所对应的t_r与当前时刻之间的时长小于Δt时，基于最新生成的电力***准稳态运行断面数据，以满足电力***都能保持输电断面有功越限安全、设备过载安全、母线电压越限安全和电网频率越限安全进行优化决策计算，若在安全稳定紧急控制实时决策计算过程中电力***不处于准稳态运行状态，则中止安全稳定紧急控制实时决策计算；

所述控制决策选择包括以下步骤：

1)若安全稳定紧急控制实时决策计算的最新决策结果中有控制措施且电力***当前运行时刻与安全稳定紧急控制实时决策计算的最新决策结果所基于的电力***准稳态运行断面数据所对应的t_r之间的时长小于设定值Δt₁，则将安全稳定紧急控制实时决策计算的最新决策结果作为用于执行的控制决策，结束本次控制决策的选择，否则，进入步骤2)；

2)若F中某个预想故障的发生时刻t_f与安全稳定紧急控制在线预决策计算中该预想故障的最新决策结果所基于的电力***准稳态运行断面数据所对应的t_r之间的时长小于设定值Δt₂且t_r早于t_f，则将安全稳定紧急控制在线预决策计算中该预想故障的最新决策结果作为用于执行的控制决策，结束本次控制决策的选择，否则，进入步骤3)；

3)若电力***当前运行时刻与安全稳定预防控制在线决策计算的最新决策结果所基于的电力***准稳态运行断面数据所对应的t_r之间的时长小于设定值Δt₃，则将安全稳定预防控制在线决策计算的最新决策结果作为用于执行的控制决策，结束本次控制决策的选择，否则，认为不需要对电力***采取控制措施，结束本次控制决策的选择；

所述电力***准稳态运行状态识别、电力***准稳态运行断面数据生成、安全稳定预防控制在线决策计算、安全稳定紧急控制在线预决策计算、安全稳定紧急控制实时决策计算和控制决策选择在电力***的运行过程中分别按循环执行的模式独立运行。

2.根据权利要求1所述的电力***安全稳定在线综合防御方法，其特征在于，所述安全稳定紧急控制在线预决策计算，针对F中各个预想故障，分别进行如下步骤的计算：

A-1)基于最新生成的电力***准稳态运行断面数据，针对预想故障进行考虑保护、安自装置和AGC/AVC控制模型及参数的时域仿真，直至仿真到各个计算节点频率动态过程中相邻的极大值与极小值之差都小于设定值Δf时，终止时域仿真，得到各个计算节点的频率动态过程、各个注入节点的有功/无功动态过程、各个输电设备的有功动态过程和各个节点的电压动态过程，进入步骤A-2)；

所述输电设备包括线路和变压器；

A-2)针对各个计算节点，分别将计算各个节点频率动态过程结束前的极大值与极小值二者的平均值作为考虑预想故障后该计算节点的稳态频率，再将所有计算节点的稳态频率的平均值作为考虑预想故障后电力***的稳态频率f_e；针对各个注入节点，分别将各个注入节点有功/无功动态过程结束前的极大值与极小值二者的平均值作为考虑预想故障后该注入节点的稳态有功/无功，并将包括等值发电机节点在内的发电机节点的发电机稳态有功之和记为P_G.sum.e；针对各个输电设备，分别将各个输电设备有功动态过程结束前的极大值与极小值二者的平均值作为考虑预想故障后该输电设备的稳态有功，再针对各个用于安全稳定监视的输电通道，分别将各个输电通道中所有输电设备的稳态有功之和作为考虑预想故障后该输电通道的稳态有功；针对各个节点，分别将各个节点电压动态过程结束前的极大值与极小值二者的平均值作为考虑预想故障后该节点的稳态电压，进入步骤A-3)；

A-3)根据考虑预想故障后各个注入节点的稳态有功/无功和预想故障后电力***的网络拓扑及模型与参数，建立电力***的潮流方程，进入步骤A-4)；

A-4)通过求解公式(1)表示的优化函数，计算出针对最新生成的电力***准稳态运行断面数据所对应的t_r时刻该预想故障的安全稳定紧急控制在线预决策：

$\left\{\begin{array}{c}\min \[\underset{i=1}{\overset{I}{\Σ}}\left(C_{G.i}\right|{\mathrm{\ΔP}}_{G.i}\left|\right)+\underset{j=1}{\overset{J}{\Σ}}\left(C_{L.j}\right|{\mathrm{\ΔP}}_{L.j}\left|\right)\]\\ s.t.\left\{\begin{array}{c}{P}_{G.i.\min .e}\≤P_{G.i\mathrm{.0.}e}+{\mathrm{\ΔP}}_{G.i}\≤P_{G.i.\max .e},i=1,2,\mathrm{...},I\\ P_{L.j.\min .e}\≤P_{L.j\mathrm{.0.}e}+{\mathrm{\ΔP}}_{L.j},j=1,2,\mathrm{...},J\\ P_{k\mathrm{.0.}e}+\underset{i=1}{\overset{I}{\Σ}}\left(s_{k.G.i.e}{\mathrm{\ΔP}}_{G.i}\right)+\underset{j=1}{\overset{J}{\Σ}}\left(s_{k.L.j.e}{\mathrm{\ΔP}}_{L.j}\right)\≤P_{k.LMT},k=1,2,\mathrm{...},K\\ P_{k\mathrm{.0.}e}+\underset{i=1}{\overset{I}{\Σ}}\left(s_{k.G.i.e}{\mathrm{\ΔP}}_{G.i}\right)+\underset{j=1}{\overset{J}{\Σ}}\left(s_{k.L.j.e}{\mathrm{\ΔP}}_{L.j}\right)\≥-P_{k.LMT},k=1,2,\mathrm{...},K\\ P_{m\mathrm{.0.}e}+\underset{n=1}{\overset{N_m}{\Σ}}\[\underset{i=1}{\overset{I}{\Σ}}\left(s_{m.n.G.i.e}{\mathrm{\ΔP}}_{G.i}\right)+\underset{j=1}{\overset{J}{\Σ}}\left(s_{m.n.L.i.e}{\mathrm{\ΔP}}_{L.j}\right)\]\≤P_{m.LMT.e},m=1,2,\mathrm{...},M\\ P_{m\mathrm{.0.}e}+\underset{n=1}{\overset{N_m}{\Σ}}\[\underset{i=1}{\overset{I}{\Σ}}\left(s_{m.n.G.i.e}{\mathrm{\ΔP}}_{G.i}\right)+\underset{j=1}{\overset{J}{\Σ}}\left(s_{m.n.L.i.e}{\mathrm{\ΔP}}_{L.j}\right)\]\≥-P_{m.n.LMT.e}^{,},m=1,2,\mathrm{...},M\\ \underset{i=1}{\overset{I}{\Σ}}{\mathrm{\ΔP}}_{G.i}-\underset{j=1}{\overset{J}{\Σ}}{\mathrm{\ΔP}}_{L.j}\≥\mathrm{\λ}\frac{P_{G.sum.r}-P_{G.sum.e}}{f_r-f_e}(f_{cr,low}-f_e)\\ \underset{i=1}{\overset{I}{\Σ}}{\mathrm{\ΔP}}_{G.i}-\underset{j=1}{\overset{J}{\Σ}}{\mathrm{\ΔP}}_{L.j}\≤\mathrm{\λ}\frac{P_{G.sum.r}-P_{G.sum.e}}{f_r-f_e}(f_{cr,up}-f_e)\\ Q_{x.\min .e}\≤Q_{x\mathrm{.0.}e}+{\mathrm{\ΔQ}}_x\≤Q_{x.\max .e},x=1,2,\mathrm{...},X\\ V_{b\mathrm{.0.}e}+\underset{x=1}{\overset{X}{\Σ}}\left(s_{vq.b.x.e}{\mathrm{\ΔQ}}_x\right)\≥V_{b.low.e},b=1,2,\mathrm{...},B\\ V_{b\mathrm{.0.}e}+\underset{x=1}{\overset{X}{\Σ}}\left(s_{vq.b.x.e}{\mathrm{\ΔQ}}_x\right)\≤V_{b.up.e},b=1,2,\mathrm{...},B\end{array}\right.\end{array}\right.---\left(1\right)$

公式(1)中，I为包括等值发电机节点在内的发电机节点中预先设置的发电机有功可控节点的总数，C_G.i、ΔP_G.i分别为第i个发电机有功可控节点的发电机单位有功调整代价和发电机有功变化量，J为包括等值负荷节点在内的负荷节点中预先设置的负荷可控节点的总数，C_L.j、ΔP_L.j分别为第j个负荷可控节点的负荷单位有功调整代价和负荷有功变化量，P_G.i.0.e、P_G.i.max.e、P_G.i.min.e分别为预想故障后第i个发电机有功可控节点的发电机稳态有功、发电机稳态有功在设定的时长Δt_c内可控区间的上限和发电机稳态有功在Δt_c内可控区间的下限，P_L.i.0.e、P_L.i.min.e分别为预想故障后第j个负荷可控节点的负荷稳态有功和负荷稳态有功在Δt_c内可控区间的下限，K为电网中的输电设备数，P_k.0.e为预想故障后第k个输电设备的稳态有功，P_k.LMT为第k个输电设备的稳态有功限额，s_k.G.i.e为基于步骤A-3)潮流方程计算出的第i个发电机有功可控节点的发电机有功对第k个输电设备有功的灵敏度，s_k.L.j.e为基于步骤A-3)潮流方程计算出的第j个负荷可控节点的负荷有功对第k个输电设备有功的灵敏度，M为用于安全稳定监视的输电通道数，N_m为第m个用于安全稳定监视的输电通道中的输电设备数，P_m.0.e、P_m.LMT.e、P’_m.LMT.e分别为预想故障后第m个用于安全稳定监视的输电通道的稳态有功、正向安全稳定限额和反向安全稳定限额，s_m.n.G.i.e为基于步骤A-3)潮流方程计算出的第i个发电机有功可控节点的发电机有功对第m个用于安全稳定监视的输电通道中第n个输电设备有功的灵敏度，s_m.n.L.i.e为基于步骤A-3)潮流方程计算出的第j个负荷可控节点的负荷有功对第m个用于安全稳定监视的输电通道中第n个输电设备有功的灵敏度，P_G.sum.r、f_r分别为与最新生成的电力***准稳态运行断面数据所对应的P_G.sum.r和f_r，f_cr.low、f_cr.up分别为电网稳态频率越限安全的下限和上限，λ为考虑仿真与实际动态过程存在差异而设定的参数，X为预先设置的无功可控节点数，Q_x.0.e、ΔQ_x、Q_x.min.e、Q_x.max.e分别为预想故障后第x个无功可控节点的无功设备的稳态无功、稳态无功变化量、稳态无功最小值和稳态无功最大值，B为电网中的节点数，V_b.0.e、V_b.low.e、V_b.up.e分别为预想故障后第b个节点的稳态电压、稳态电压下限和稳态电压上限，s_vq.b.x.e为基于步骤A-3)潮流方程计算出的第x个无功可控节点的无功设备的无功对第b个节点电压的灵敏度；

所述输电设备的稳态有功限额是指根据输电设备的过载电流限额和额定电压，按指定功率因数计算出的稳态有功。

3.根据权利要求1所述的电力***安全稳定在线综合防御方法，其特征在于，针对安全稳定紧急控制实时决策计算包括如下步骤：

B-1)根据最新生成的电力***准稳态运行断面数据，建立电力***的潮流方程，进入步骤B-2)；

2)通过求解公式(2)表示的优化函数，计算出针对最新生成的电力***准稳态运行断面数据所对应的t_r时刻电力***的安全稳定紧急控制实时决策：

$\left\{\begin{array}{c}\min \[\underset{i=1}{\overset{I}{\Σ}}\left(C_{G.i}\right|{\mathrm{\ΔP}}_{G.i}\left|\right)+\underset{j=1}{\overset{J}{\Σ}}\left(C_{L.j}\right|{\mathrm{\ΔP}}_{L.j}\left|\right)\]\\ s.t.\left\{\begin{array}{c}{P}_{G.i.\min }\≤P_{G.i.0}+{\mathrm{\ΔP}}_{G.i}\≤P_{G.i.\max },i=1,2,\mathrm{...},I\\ P_{L.j.\min }\≤P_{L.j.0}+{\mathrm{\ΔP}}_{L.j},j=1,2,\mathrm{...},J\\ P_{k.0}+\underset{i=1}{\overset{I}{\Σ}}\left(s_{k.G.i}{\mathrm{\ΔP}}_{G.i}\right)+\underset{j=1}{\overset{J}{\Σ}}\left(s_{k.L.j}{\mathrm{\ΔP}}_{L.j}\right)\≤P_{k.LMT},k=1,2,\mathrm{...},K\\ P_{k.0}+\underset{i=1}{\overset{I}{\Σ}}\left(s_{k.G.i}{\mathrm{\ΔP}}_{G.i}\right)+\underset{j=1}{\overset{J}{\Σ}}\left(s_{k.L.j}{\mathrm{\ΔP}}_{L.j}\right)\≥-P_{k.LMT},k=1,2,\mathrm{...},K\\ P_{m.0}+\underset{n=1}{\overset{N_m}{\Σ}}\[\underset{i=1}{\overset{I}{\Σ}}\left(s_{m.n.G.i}{\mathrm{\ΔP}}_{G.i}\right)+\underset{j=1}{\overset{J}{\Σ}}\left(s_{m.n.L.i}{\mathrm{\ΔP}}_{L.j}\right)\]\≤P_{m.LMT},m=1,2,\mathrm{...},M\\ P_{m.0}+\underset{n=1}{\overset{N_m}{\Σ}}\[\underset{i=1}{\overset{I}{\Σ}}\left(s_{m.n.G.i}{\mathrm{\ΔP}}_{G.i}\right)+\underset{j=1}{\overset{J}{\Σ}}\left(s_{m.n.L.i}{\mathrm{\ΔP}}_{L.j}\right)\]\≥-P_{m.n.LMT}^{,},m=1,2,\mathrm{...},M\\ \underset{i=1}{\overset{I}{\Σ}}{\mathrm{\ΔP}}_{G.i}-\underset{j=1}{\overset{J}{\Σ}}{\mathrm{\ΔP}}_{L.j}\≥\mathrm{\λ}\frac{P_{G.sum.r.0}-P_{G.sum.r}}{f_{r.0}-f_r}(f_{cr,low}-f_r)\\ \underset{i=1}{\overset{I}{\Σ}}{\mathrm{\ΔP}}_{G.i}-\underset{j=1}{\overset{J}{\Σ}}{\mathrm{\ΔP}}_{L.j}\≤\mathrm{\λ}\frac{P_{G.sum.r.0}-P_{G.sum.r}}{f_{r.0}-f_r}(f_{cr,up}-f_r)\\ Q_{x.\min }\≤Q_{x.0}+{\mathrm{\ΔQ}}_x\≤Q_{x.\max },x=1,2,\mathrm{...},X\\ V_{b.0}+\underset{x=1}{\overset{X}{\Σ}}\left(s_{vq.b.x}{\mathrm{\ΔQ}}_x\right)\≥V_{b.low},b=1,2,\mathrm{...},B\\ V_{b.0}+\underset{x=1}{\overset{X}{\Σ}}\left(s_{vq.b.x}{\mathrm{\ΔQ}}_x\right)\≤V_{b.up},b=1,2,\mathrm{...},B\end{array}\right.\end{array}\right.---\left(2\right)$

公式(2)中，I为包括等值发电机节点在内的发电机节点中预先设置的发电机有功可控节点的总数，C_G.i、ΔP_G.i分别为第i个发电机有功可控节点的发电机单位有功调整代价和发电机有功变化量，J为包括等值负荷节点在内的负荷节点中预先设置的负荷可控节点的总数，C_L.j、ΔP_L.j分别为第j个负荷可控节点的负荷单位有功调整代价和负荷有功变化量，P_G.i.0、P_G.i.max、P_G.i.min分别为第i个发电机有功可控节点的发电机稳态有功、发电机稳态有功在Δt_c内可控区间的上限和发电机稳态有功在Δt_c内可控区间的下限，P_L.i.0、P_L.i.min分别为第j个负荷可控节点的负荷稳态有功和负荷稳态有功在Δt_c内可控区间的下限，K为电网中的输电设备数，P_k.0为第k个输电设备的稳态有功，P_k.LMT为第k个输电设备的稳态有功限额，s_k.G.i为基于步骤B-1)潮流方程计算出的第i个发电机有功可控节点的发电机有功对第k个输电设备有功的灵敏度，s_k.L.j为基于步骤B-1)潮流方程计算出的第j个负荷可控节点的负荷有功对第k个输电设备有功的灵敏度，M为用于安全稳定监视的输电通道数，N_m为第m个用于安全稳定监视的输电通道中的输电设备数，P_m.0、P_m.LMT、P’_m.LMT分别为第m个用于安全稳定监视的输电通道的稳态有功、正向安全稳定限额和反向安全稳定限额，s_m.n.G.i为基于步骤B-1)潮流方程计算出的第i个发电机有功可控节点的发电机有功对第m个用于安全稳定监视的输电通道中第n个输电设备有功的灵敏度，s_m.n.L.i为基于步骤B-1)潮流方程计算出的第j个负荷可控节点的负荷有功对第m个用于安全稳定监视的输电通道中第n个输电设备有功的灵敏度，P_G.sum.r、f_r分别为与最新生成的电力***准稳态运行断面数据所对应的P_G.sum.r和f_r，P_G.sum.r.0、f_r.0分别为与前一次生成的电力***准稳态运行断面数据所对应的P_G.sum.r和f_r，f_cr.low、f_cr.up分别为电网稳态频率越限安全的下限和上限，X为预先设置的无功可控节点数，Q_x.0、ΔQ_x、Q_x.min、Q_x.max分别为第x个无功可控节点的无功设备的稳态无功、稳态无功变化量、稳态无功最小值和稳态无功最大值，B为电网中的节点数，V_b.0、V_b.low、V_b.up分别为第b个节点的稳态电压、稳态电压下限和稳态电压上限，s_vq.b.x为基于步骤B-1)潮流方程计算出的第x个无功可控节点的无功设备的无功对第b个节点电压的灵敏度。