CN103903196A - 一种考虑设备老化因素的电网调度运行实时风险评估方法 - Google Patents

Abstract

本发明公开了一种考虑设备老化因素的电网调度运行实时风险评估方法，步骤是：1)读取将要执行的调度操作指令序列，确定各步调度操作指令之间的调度运行状态，同时获取当前时刻的非电力信息，该非电力信息包括设备运行信息和社会状态信息；2)确定每个调度运行状态的预想故障场景集合；3)计算每个调度运行状态所对应的预想故障场景集合中每个预想故障场景的实时发生概率；4)确定使得预想故障场景从失效状态恢复到正常工作状态所需要最小控制费用时的负荷削减量；5)计算每个调度运行状态下的风险指标；6)根据实际情况制定风险控制措施，以保证调度操作指令序列的正常执行，有助于提高电网安全稳定运行能力。

Description

一种考虑设备老化因素的电网调度运行实时风险评估方法

技术领域

本发明涉及电网调度运行领域，特别是涉及一种考虑设备老化因素的电网调度运行实时风险评估方法。

背景技术

为了保证电网安全、经济和可靠运行，调度员必需采取合理的调度操作指令。目前调度员主要依据电网实时运行数据以及个人经验对电网进行操作。实际中，调度操作指令票通常为一个调度操作指令序列，包含若干条调度操作指令，只有前一条调度操作指令正常执行完成之后，才可以执行下一条指令。两条调度操作指令之间存在一调度运行状态，而调度操作指令执行过程中的恶劣天气和设备故障等因素可能造成调度运行状态风险水平过高，影响下一条调度操作指令的正常执行，威胁系统安全稳定运行。

当前已投入应用的电网运行安全实时风险评估方法只能够评估既定的电网运行方式，即某个调度指令票执行之前或完成后的电网稳定状态，还没有可以对调度操指令序列中的调度运行状态进行风险评估的方法。

发明内容

本发明的目的在于提出一种考虑设备老化因素的电网调度运行实时风险评估方法，能够在调度操作指令执行之前预先对每个调度运行状态进行风险评估，判断该调度运行状态风险水平是否可接受，是否影响下一条调度操作指令的正常执行，可以保证调度操作指令序列的正常执行，有助于提高电网安全稳定运行能力。

本发明一种考虑设备老化因素的电网调度运行实时风险评估方法，包括以下步骤：

步骤1)读取将要执行的调度操作指令序列，确定各步调度操作指令之间的调度运行状态，同时获取当前时刻的非电力信息，该非电力信息包括设备运行信息和社会状态信息；其中，所述设备运行信息是指设备的投运时间和距离当前时刻最近的一次检修时间，所述社会状态信息是指当前时期的保供电等级；

步骤2)根据系统当前设备投运情况，结合所述调度操作指令序列，确定每个调度运行状态对应的预想故障场景，从而确定每个调度运行状态的预想故障场景集合；

步骤3)根据电网设备历史数据、当前时刻的设备运行信息，确定系统设备的实时不可用率，并且应用状态枚举法计算每个调度运行状态所对应的预想故障场景集合中每个预想故障场景的实时发生概率；

步骤4)对每个预想故障场景进行潮流计算分析，如果出现设备过载或者节点电压越限，则判定为失效状态，对处于失效状态的预想故障场景进行基于最优潮流的最优控制，确定使得预想故障场景从失效状态恢复到正常工作状态所需要最小控制费用时的负荷削减量；

步骤5)根据步骤3)得到的每个预想故障场景的实时发生概率及步骤4)得到的负荷削减量，计算每个调度运行状态下的风险指标；

步骤6)分析整个调度操作指令序列的所有调度运行状态的风险指标数值，如果某调度运行状态风险指标数值超过可接受限值，需要发布风险预警信息，并根据实际情况制定风险控制措施，确保下一步调度操作指令正常执行。

其中，步骤1)中，

假定S₀表示系统当前状态，将要执行的调度操作指令序列为D，共包括m个调度操作指令，

D＝[d₁,d₂,...,d_m]

其中，d₁,d₂,…,d_m为m个调度操作指令，同时确定调度操作指令d₁,d₂,…,d_m的目的设备分别为e_k,e_k+1,…,e_k+m(k>1且k+m<n)；

调度操作指令序列D中共有m-1个调度运行状态，分别为：

S₁＝(S₀(±e_k))

S₂＝(S₁(±e_k+1))

............

S_m-1＝(S_m-2(±e_k+m-1))

其中，-e_k表示设备e_k退出运行，+e_k表示设备e_k投入运行，S₁＝(S₀(±e_k))调度运行状态S₁是在S₀的基础上设备e_k投入/退出运行，S₂＝(S₁(±e_k+1))表示调度运行状态S₂是在S₁的基础上设备e_k+1投入/退出运行，S_m-1＝(S_m-2(±e_k+m-1))调度运行状态S_m-1是在S_m-2基础上设备e_k+m-1投入/退出运行。

与现有技术相比，本发明的有益效果是：

本发明从电网调度运行出发，结合每个调度操作指令的目的设备，确定调度操作指令序列执行过程中的调度运行状态；应用电网当前普遍采用的N-1分析和N-2分析，确定每个调度运行状态的预想故障场景集合；根据系统当前实时数据，包括设备运行数据以及电网历史数据，计算出设备的实时不可用率，一方面可以反映出历史统计数据的可靠性，另一方面可以反映出实时影响因子的准确性，进而计算出每个预想故障场景的实时发生概率；选取预想故障场景进行潮流分析，对判定为失效状态的预想故障场景进行基于最优潮流的最优控制，得到使得预想故障场景从失效状态恢复到正常状态所需最小控制费用的负荷削减量；根据每个预想故障场景的实时发生概率和负荷削减量，计算本调度运行状态的期望缺供电量指标，如果本调度运行状态风险水平超出可接受限值，需要发布风险预警，并制定风险预控措施，以保证下一条调度操作指令正常执行。本发明能够充分考虑到设备运行因素和社会状态因素的影响，对调度运行状态进行风险评估，可以提前预测出风险水平超过可接受限值的调度运行状态，及时制定风险预控措施，保证整个调度操作指令序列的正常有序执行，有助于提高电网安全稳定运行能力。

附图说明

图1是本发明方法评估方法的流程图；

图2是本发明评估方法中设备老化因子示意图；

图3是本发明评估方法中设备检修因子示意图；

图4是本发明评估方法一具体实施例的系统接线图。

具体实施方式

下面结合附图和具体实施例对本发明技术方案作进一步详细描述。

本发明提出一种考虑设备老化因素的电网调度运行实时风险评估方法，该方法以电网历史数据及实时数据主要包括设备运行因素为基础，建立系统设备的时变概率模型，能够计算出系统设备的实时不可用率；同时应用基于最优潮流的最优控制方法对处于失效状态的预想故障场景进行负荷削减计算，能够提前对调度操作指令序列中的调度运行状态进行风险评估，预先对风险水平超过可接受限值的调度运行状态进行风险预警，以便及时制定风险控制措施，保证调度操作指令序列的正常有序执行，增强电网的安全稳定运行能力。

图1为本发明提出的一种设备老化因素的电网调度运行实时风险评估方法的具体流程图，其中，101-102表示根据读取的调度操作序列信息，确定调度运行状态以及每个调度运行状态所对应的预想故障场景，102-103表示对计算每个调度运行状态所对应预想故障场景的实时发生概率，102-104表示计算处于失效状态的预想故障场景的最小负荷削减量，103-105表示调度运行状态的风险指标计算需要103步骤的预想故障场景的实时发生概率，101-105表示调度运行状态的风险指标计算需要101步骤的设备运行信息与社会状态信息，104-105表示调度运行状态的风险指标计算需要104步骤的最小负荷削减量，105-106表示根据调度运行状态的风险指标确定是否是需要制定风险控制措施。为更好的理解本发明，下面以一具体实施案例进行说明，图4为IEEE RTS测试系统接线图；其中：编号1、2、7、13、14、15、16、18、21、22和23均为发动机节点，编号4、5、6、8、9、10、17、19和20均为负荷节点，编号11和12均为联络节点，#1B、#2B、#3B、#4B、#5B和#6B均为变压器；图中虚线框内为某变电站主接线图，站内采用双母并列运行方式(断路器D4处于闭合状态)，两台变压器分别由母线#1M和#2M供电，隔离开关G5、G6和G8处于拉开状态，隔离开关G1和G7处于闭合状态假定母线#2M上PT(电压互感器)一次侧发生故障，必须先将母线#2M退出运行，才能将PT退出检修。将母线#2M退出运行的调度操作方案为：先将变压器#2B停电倒至母线#1M上，再将母线#2M退出。设定调度运行状态可接受风险限值为1.5MW，当前时刻的保供电等级为一般时期，所有设备的设备老化因子和设备检修因子均为1.0，应用本发明对案例调度操作指令序列进行风险评估的具体步骤如下：

步骤1)读取将要执行的调度操作指令序列，以确定各步调度操作指令之间的调度运行状态，同时要获取当前时刻的非电力信息，包括设备运行信息和社会状态信息；

其中设备运行信息，主要指设备的投运时间和距离当前时刻最近的一次检修时间，社会状态信息，主要指当前时期的保供电等级；

假定S₀表示系统当前状态，将要执行的调度操作指令序列为D，共包括m个调度操作指令，

D＝[d₁,d₂,...,d_m]

其中，d₁,d₂,…,d_m为m个调度操作指令，同时确定调度操作指令d₁,d₂,…,d_m的目的设备分别为e_k,e_k+1,…,e_k+m(k>1且k+m<n)；

调度操作指令序列D中共有m-1个调度运行状态，分别为：

S₁＝(S₀(±e_k))

S₂＝(S₁(±e_k+1))

............

S_m-1＝(S_m-2(±e_k+m-1))

其中，-e_k表示设备e_k退出运行，+e_k表示设备e_k投入运行，S₁＝(S₀(±e_k))调度运行状态S₁是在S₀的基础上设备e_k投入/退出运行，S₂＝(S₁(±e_k+1))表示调度运行状态S₂是在S₁的基础上设备e_k+1投入/退出运行，S_m-1＝(S_m-2(±e_k+m-1))调度运行状态S_m-1是在S_m-2基础上设备e_k+m-1投入/退出运行。

调度操作指令序列及每条调度操作指令的目的设备，如表1所示；

表1调度操作指令序列

根据调度操作指令序列内容，可以确定本调度操作指令序列中共有4个调度运行状态，分别为第一步、第二步、第三步和第四步调度操作指令序列执行之后的系统状态，其中第五步调度操作指令执行之后的状态不属于调度运行状态，属于调度操作指令序列执行完毕状态，如下所示，

S₁＝(S₀(-D5))

S₂＝(S₁(-G7))

S₃＝(S₂(+G6))

S₄＝(S₃(+D5))

其中，S₀表示执行调度指令序列D之前的系统当前状态，如图4所示，所有变压器和线路都处于运行状态，站内采用双母并列运行方式(断路器D4处于闭合状态)，两台变压器分别由母线#1M和#2M供电，隔离开关G5、G6和G8处于拉开状态，隔离开关G1和G7处于闭合状态。

步骤2)根据系统当前设备投运情况，结合具体的调度操作指令序列，确定每个调度运行状态可能发生的预想故障场景，从而确定每个调度运行状态的预想故障场景集合；

调度运行状态的预想故障场景包括两种情况，第一种情况是单一设备退出运行(电网普遍采用的N-1预想事故分析)；第二种是两个设备同时退出运行；

假定调度运行状态S_k共有s个投运设备，分别为e₁,e₂,…,e_s，其中设备e₅和e₆为同塔双回线，设备e_s-1和e_s为同塔双回线，则调度运行状态S_k的单一设备退出运行预想故障场景集合为，

S_F1,k＝{S_F1,1,...,S_F1,4,S_F1,7,...,S_F1,s-2}

S_F1,j＝S_k(-e_j)(j＝1,2,...,4,7,8,...,s-2)

其中，S_F1,k表示调度运行状体S_k的单一设备退出运行预想故障场景集合，S_F1,k中共包括s-4个预想故障场景，S_F1,j＝S_k(-e_j)表示第j个单一设备退出运行预想故障场景是调度运行状态S_k基础上设备e_j退出运行；

调度运行状态S_k的两个设备同时退出运行预想故障场景集合为，

S_F2,k＝{S_F2,1,S_F2,2}

S_F2,1＝S_k(-e₅,-e₆)

S_F2,2＝S_k(-e_s-1,-e_s)

其中，S_F2,k表示调度运行状态S_k的两个设备同时退出运行预想故障场景集合，S_F2,k中共包括2个预想故障场景，S_F2,1＝S_k(-e₅,-e₆)表示第1个两个设备同时退出运行预想故障场景是调度运行状态S_k基础上设备e₅和e₆同时退出运行，S_F2,2＝S_k(-e_s-1,-e_s)表示第2个两个设备同时退出运行预想故障场景是调度运行状态S_k基础上设备e_s和e_s-1同时退出运行；

如图4所示，本例中以调度运行状态S₁为例，其第i个单一设备退出运行预想故障场景S_F1,i与第j个两个设备同时退出运行预想故障场景S_F2,j如下所示，

S_F1,i＝S₁(-L_17,18)

S_F2,j＝S₁(-L¹ _15,21,-L² _15,21)

其中，L_17,18表示图4中节点17与18之间的线路，L¹ _15,21表示图4中节点15与21之间的第一回线路，L² _15,21表示图4中节点15与21之间的第二回线路。

步骤3)根据电网设备历史数据、当前时刻的设备运行信息，确定系统设备的实时不可用率。并且应用状态枚举法计算每个调度运行状态所对应预想故障场景集合中每个预想故障场景的实时发生概率；

1)确定系统所有设备基于历史统计数据的平均不可用率，

$U_i\left(0\right)=\frac{{\mathrm{\&lambda;}}_i}{{\mathrm{\&lambda;}}_i+{\mathrm{\&mu;}}_i}$

其中，U_i(0)为系统第i个设备的平均不可用率，λ_i为系统第i个设备的失效率(失效次数/年)，μ_i为系统第i个设备的修复率(修复次数/年)，λ_i和μ_i可以通过从电网的历史数据库中得到；

2)根据设备运行信息计算系统所有设备的实时不可用率，

U_i(t)＝β_1,i(t)β_2,i(t)U_i(0)

其中，U_i(t)为第i个设备的实时不可用率，U_i(0)为系统第i个设备的平均不可用率，β_1,i(t)为第i个设备的设备老化因子，β_2,i(t)为第i个设备的设备检修因子；

β_1，i(t)为设备老化因子，其关于时间的变化曲线如图2所示，具体函数表达式如下，

$\left\{\begin{array}{cc}{\mathrm{\&beta;}}_{1,i}\left(t\right)=1+\frac{e^{-\frac{{(\ln t-\ln T_1)}^2}{2}}}{\sqrt{2\mathrm{\&pi;}}\&CenterDot;t[1-\mathrm{\&Phi;}(\ln t-\ln T_1)]},& t>T_1\\ {\mathrm{\&beta;}}_{1,i}\left(t\right)=1,& t\&le;T_1\end{array}\right.$

其中，T₁和t分别为第i个设备的有效寿命和服役年龄，Φ(t)为标准正态分布函数；

β_2,i(t)为设备检修因子，可用图3所示曲线近似表示，具体函数表达式如下，

$\left\{\begin{array}{cc}{\mathrm{\&beta;}}_{2,i}\left(t\right)=\frac{t}{T_2},& t>T_2\\ {\mathrm{\&beta;}}_{2,i}\left(t\right)=1,& t\&le;T_2\end{array}\right.$

其中，t为第i个设备距离上次检修的时间间隔，T₂为第i个设备的平均检修时间；

3)计算每个调度运行状态对应预想故障场景集合中每个预想故障场景的实时发生概率，对于单一设备退出运行预想故障场景，

$P_{F1}\left(t\right)=U\left(t\right)\underset{i=1}{\overset{n}{\mathrm{\&Pi;}}}(1-U_i\left(t\right))$

其中，P_F1(t)为单一设备退出运行预想故障场景的实时发生概率，U(t)为退出运行的单设备实时不可用率，n为处于正常工作状态设备数量，U_i(t)为第i个处于正常工作状态设备的实时不可用率；

对于两个设备同时退出运行预想故障场景，

$P_{F2}\left(t\right)=U^1\left(t\right)U^2\left(t\right)\underset{i=1}{\overset{n-1}{\mathrm{\&Pi;}}}(1-U_i\left(t\right))$

其中，P_F2(t)为两个设备同时退出运行预想故障场景的实时发生概率，U¹(t)为退出运行的第一个设备的实时不可用率，U²(t)为退出运行的第二个设备的实时不可用率，n-1为处于正常工作状态设备数量，U_i(t)为第i个处于正常工作设备的实时不可用率；

步骤4)对每个预想故障场景进行潮流计算分析，如果出现设备过载或者节点电压越限，则判定为失效状态，对处于失效状态的预想故障场景进行基于最优潮流的最优控制，确定使得预想故障场景从失效状态恢复到正常工作状态所需要最小控制费用时的负荷削减量；

1)对选取的预想故障场景进行潮流计算，得到每个节点的电压和每条支路的潮流，如果存在电压越限节点或潮流过载支路，则判定该预想故障场景为失效状态；否则判定该预想故障场景为正常状态，不会对系统产生不利影响；以下选取调度运行状态S₁的4个单一设备退出运行预想故障场景与4个两个设备同时退出运行预想故障场景，每个预想故障场景的判断结果如下表1所示，

表1预想故障场景判断结果

2)对判定为失效状态的预想故障场景进行基于最优潮流的最优控制，计算得到使得预想故障场景从失效状态恢复到正常状态所需最小控制费用时的负荷削减量，最优潮流模型如下，

$\min f({\mathrm{\&Delta;P}}_{G1},...,{\mathrm{\&Delta;P}}_{\mathrm{Gm}},...,{\mathrm{\&Delta;P}}_{L1},...,{\mathrm{\&Delta;P}}_{\mathrm{Ll}})=\underset{k=1}{\overset{m}{\mathrm{\&Sigma;}}}{C}_{\mathrm{Gk}}{\mathrm{\&Delta;P}}_{\mathrm{Gk}}+\underset{h=1}{\overset{l}{\mathrm{\&Sigma;}}}{C}_{\mathrm{Lh}}{\mathrm{\&Delta;P}}_{\mathrm{Lh}}$

$s.t.\underset{k=1}{\overset{m}{\mathrm{\&Sigma;}}}\frac{{\&PartialD;I}_{\mathrm{ij}}}{{\&PartialD;P}_{\mathrm{Gk}}}{\mathrm{\&Delta;P}}_{\mathrm{Gk}}+\underset{h=1}{\overset{l}{\mathrm{\&Sigma;}}}\frac{{\&PartialD;I}_{\mathrm{ij}}}{{\&PartialD;P}_{\mathrm{Lh}}}{\mathrm{\&Delta;P}}_{\mathrm{Lh}}\&le;I_{\mathrm{ij}}^{\max }-I_{\mathrm{ij}}^0$

$\underset{k=1}{\overset{m}{\mathrm{\&Sigma;}}}\frac{{\&PartialD;P}_s}{{\&PartialD;P}_{\mathrm{Gk}}}{\mathrm{\&Delta;P}}_{\mathrm{Gk}}+\underset{h=1}{\overset{l}{\mathrm{\&Sigma;}}}\frac{{\&PartialD;P}_s}{{\&PartialD;P}_{\mathrm{Lh}}}{\mathrm{\&Delta;P}}_{\mathrm{Lh}}=0$

$P_{\mathrm{Gk}}^{\min }-P_{\mathrm{Gk}}^0\&le;{\mathrm{\&Delta;P}}_{\mathrm{Gk}}\&le;P_{\mathrm{Gk}}^{\max }-P_{\mathrm{Gk}}^0,1\&le;k\&le;m$

$0\&le;{\mathrm{\&Delta;P}}_{\mathrm{Lh}}\&le;P_{\mathrm{Lh}}^0,1\&le;h\&le;l$

其中，m和l分别为发电机节点和负荷节点数量，△P_Gk为发电机节点k的有功出力变化量，△P_Lh负荷节点h的负荷削减量，C_Gk为发电机节点k的单位有功出力变化的控制费用，C_Lh为负荷节点h的单位负荷削减的控制费用，i和j为线路的首端节点和末端节点编号，Ii_j表示线路ij的电流，为线路ij允许通过电流的最大值，I_ij ⁰为线路ij进行最优控制之前的电流值，

和

表示线路ij通过电流对发电机节点k和负荷节点h的灵敏度系数，P_s表示平衡节点s的有功出力，

和

表示平衡节点s的有功出力对发电机节点k和负荷节点h的灵敏度系数，P_Gk ⁰为发电机节点k进行最优控制之前的有功出力，P_Gk ^max和P_Gk ^min分别为发电机节点k有功出力的上限和下限，P_Lh ⁰为负荷节点h进行最优控制之前的负荷值；

经过最优潮流分析可以确定哪些负荷节点发生负荷削减，从而得到预想故障场景的负荷削减总量，

${\mathrm{\&Delta;P}}_L=\underset{j=1}{\overset{x}{\mathrm{\&Sigma;}}}\mathrm{\&Delta;P}\left(j\right)$

其中，△P_L为预想故障场景的负荷削减总量，x表示出现负荷削减的节点数量，△P(j)为节点j的负荷削减量；

对表1中判定为失效状态的预想故障场景进行基于最优潮流的最优控制，计算结果如下表2所示，

表2预想故障场景基于最优潮流的最优控制结果

步骤5)根据每个预想故障场景的实时发生概率及负荷削减量等相关数据，计算得到每个调度运行状态的风险指标；

调度运行状态风险指标能够反映出系统的平均风险水平，风险指标计算公式为，

$R_L=\mathrm{\&psi;}\underset{S\&Element;A}{\mathrm{\&Sigma;}}P\left(S\right){\mathrm{\&Delta;P}}_L\left(S\right)$

其中，R_L表示期望缺供电量指标，A表示本调度运行状态所对应的出现负荷削减的预想故障场景集合，S表示出现负荷削减的预想故障场景，P(S)表示预想故障场景S的实时发生概率，△P_L(S)为预想故障场景S的负荷削减总量，ψ表示社会影响因子，可根据当前时刻的保供电等级确定，

表2社会影响因子

每个调度运行状态的期望缺供电量指标如表3所示。

表3期望缺供电量指标

步骤6)分析整个调度操作指令序列的所有调度运行状态的风险指标数值，如果某调度运行状态风险水平超过可接受限值，需要发布风险预警信息，并根据实际情况制定风险控制措施，确保下一步调度操作指令可以正常执行；

根据计算得到的整个调度操作指令序列所有调度运行状态的期望缺供电量指标，判断该风险水平是否超过可接受限值，如果某调度运行状态风险水平过高，需要制定风险预控措施，包括确定该调度运行状态所有预想故障场景出现的所有负荷削减节点，确保相应负荷节点的负荷切除装置能够正常快速动作，在发生故障时可以保证通过切除部分负荷维持系统稳定，同时执行下一步调度操作指令之前要加强对系统设备的巡检工作，降低系统设备发生故障概率；

如表3所示，执行第一步、第二步、第三步和第四步完毕之后的调度运行状态风险指标都小于可接受风险限值1.5MW，表明整个调度操作指令序列执行过程中，系统风险水平都处于可接受范围之内，不需要发布风险预警信息；但调度运行状态1和3的风险指标数值已非常靠近可接受风险限值，在第二步和第四步调度操作指令执行之前可制定风险控制措施，确保调度操作指令正常执行。

尽管上面结合附图对本发明进行了描述，但是本发明并不局限于上述的具体实施方式，上述的具体实施方式仅仅是示意性的，而不是限制性的，本领域的普通技术人员在本发明的启示下，在不脱离本发明宗旨的情况下，还可以做出很多变形，这些均属于本发明的保护之内。

Claims (2)

1.一种考虑设备老化因素的电网调度运行实时风险评估方法，其特征在于，包括以下步骤：

步骤1)读取将要执行的调度操作指令序列，确定各步调度操作指令之间的调度运行状态，同时获取当前时刻的非电力信息，该非电力信息包括设备运行信息和社会状态信息；其中，所述设备运行信息是指设备的投运时间和距离当前时刻最近的一次检修时间，所述社会状态信息是指当前时期的保供电等级；

步骤2)根据系统当前设备投运情况，结合所述调度操作指令序列，确定每个调度运行状态对应的预想故障场景，从而确定每个调度运行状态的预想故障场景集合；

步骤3)根据电网设备历史数据、当前时刻的设备运行信息，确定系统设备的实时不可用率，并且应用状态枚举法计算每个调度运行状态所对应的预想故障场景集合中每个预想故障场景的实时发生概率；

步骤4)对每个预想故障场景进行潮流计算分析，如果出现设备过载或者节点电压越限，则判定为失效状态，对处于失效状态的预想故障场景进行基于最优潮流的最优控制，确定使得预想故障场景从失效状态恢复到正常工作状态所需要最小控制费用时的负荷削减量；

步骤5)根据步骤3)得到的每个预想故障场景的实时发生概率及步骤4)得到的负荷削减量，计算每个调度运行状态下的风险指标；

步骤6)分析整个调度操作指令序列的所有调度运行状态的风险指标数值，如果某调度运行状态风险指标数值超过可接受限值，需要发布风险预警信息，并根据实际情况制定风险控制措施，确保下一步调度操作指令正常执行。

2.根据权利要求1所述考虑设备老化因素的电网调度运行实时风险评估方法，其中，步骤1)中，

假定S₀表示系统当前状态，将要执行的调度操作指令序列为D，共包括m个调度操作指令，

D＝[d₁,d₂,...,d_m]

其中，d₁,d₂,…,d_m为m个调度操作指令，同时确定调度操作指令d₁,d₂,…,d_m的目的设备分别为e_k,e_k+1,…,e_k+m(k>1且k+m<n)；

调度操作指令序列D中共有m-1个调度运行状态，分别为：

S₁＝(S₀(±e_k))

S₂＝(S₁(±e_k+1))

............

S_m-1＝(S_m-2(±e_k+m-1))

其中，-e_k表示设备e_k退出运行，+e_k表示设备e_k投入运行，S₁＝(S₀(±e_k))调度运行状态S₁是在S₀的基础上设备e_k投入/退出运行，S₂＝(S₁(±e_k+1))表示调度运行状态S₂是在S₁的基础上设备e_k+1投入/退出运行，S_m-1＝(S_m-2(±e_k+m-1))调度运行状态S_m-1是在S_m-2基础上设备e_k+m-1投入/退出运行。

Images