计及保护与安自控制代价的在线预防控制风险决策方法与流程

技术特征：

1.计及保护与安自控制代价的在线预防控制风险决策方法，其特征在于，包括以下步骤：

1)基于电网当前的运行状态和设备模型及参数，生成反映电网当前运行状态、用于电网安全稳定在线分析的潮流数据和稳定数据，再结合电网预测信息和计划信息，生成预防控制候选措施集，分别作为最新生成的电网运行断面潮流数据、稳定数据和预防控制候选措施集，并将电网运行断面所对应的时刻记为t₀，进入步骤2)；

所述预防控制候选措施集是指当前运行方式下可选的预防控制设备及其控制变量、控制方向和可控量最大值的集合；

2)将预想故障发生的概率大于预先设定的门槛值p_cr的预想故障的集合记为预想故障集F₀，进入步骤3)；

所述预想故障是指电网在t₀至(t₀+Δt_f)的时段内可能发生的1个或多个设备故障的组合；

所述预想故障发生的概率是指电网在t₀至(t₀+Δt_f)的时段内该预想故障发生的次数；

其中，Δt_f为设定的预防控制防御时长；

3)基于最新生成的电网运行断面潮流数据和稳定数据，采用集群计算平台，分别针对F₀中每个预想故障，计及保护装置和安自装置的运行状态及其模型和参数，通过公式(1)计算出因预想故障直接触发的保护装置和安自装置动作所导致的风险代价，进入步骤4)；

L_rc.i＝(L_r.i+L_c.i)p_f.i i＝1,2,…,n (1)

式中，n为F₀中预想故障的数目，L_r.i为在最新生成的电网运行断面潮流数据所对应的电网运行状态下F₀中第i个预想故障发生后因预想故障直接触发的保护装置动作所导致的损失代价，L_c.i为在最新生成的电网运行断面潮流数据所对应的电网运行状态下F₀中第i个预想故障发生后因预想故障直接触发的安自装置动作所导致的损失代价，p_f.i为F₀中第i个预想故障发生的概率，L_rc.i为在最新生成的电网运行断面潮流数据所对应的电网运行状态下F₀中第i个预想故障发生后因预想故障直接触发的保护装置和安自装置动作所导致的风险代价；

4)按照预想故障发生的概率由大到小的顺序对F₀中的所有预想故障进行排序，得到预想故障排序表，其中对于发生概率相同的预想故障则按照步骤3)中计算得到的预想故障的风险代价由大到小的顺序确定这些发生概率相同的预想故障之间的顺序，进入步骤5)；

5)设置在线预防控制优化措施集C_p为空集，设置预想故障排序指针i_s为0，设置计算标志I_c为0，进入步骤6)；

所述在线预防控制优化措施集是指当前运行方式下选定的预防控制设备及其控制变量、控制方向和最终控制量的集合；

6)若当前时刻t与t₀之差小于设定的预防控制优化决策时长限值Δt_c，进入步骤7)，否则，将C_p作为最终的在线预防控制优化措施集提交给调度运行人员，结束本方法；

其中，Δt_c小于Δt_f；

7)若i_s小于n，则依据预想故障排序表中预想故障的排序，将排序自第(i_s+1)位开始直至第min{i_s+max[int(k_pn_p),1],n]}位的共{min{i_s+max[int(k_pn_p),1],n]}-i_s}个预想故障的集合记为F_c，并将i_s置为min{i_s+max[int(k_pn_p),1],n]}，进入步骤8)，否则，将C_p作为最终的在线预防控制优化措施集提交给调度运行人员，结束本方法；

其中，n_p为用于电网安全稳定在线预防控制优化决策计算的集群计算平台中所有CPU的总核数，k_p为设定的小于1的参数；

8)基于最新生成的电网运行断面潮流数据和稳定数据，采用集群计算平台，分别针对F_c中每个预想故障，计算出因预想故障直接触发的保护装置动作所导致的损失代价L_r，并计及保护装置和安自装置的运行状态及其模型和参数，进行预想故障下的电网安全稳定时域仿真计算，计算出因预想故障发生所导致的损失代价L_f，进入步骤9)；

9)若F_c中L_f大于L_r的预想故障数为0，则返回步骤6)；否则，将F_c中L_f大于L_r的预想故障的集合记为F_p0，将F_c中除F_p0中预想故障之外的预想故障的集合记为F_s0，通过公式(2)计算出F_p0中各个预想故障发生后所导致的风险代价总和SL_fp0，并设立集合F_p，将F_p的初始值设置为F_p0，将F_p中的各个预想故障发生后所导致的风险代价总和SL_fp的初始值设置为SL_fp0，进入步骤10)：

SL f p 0 = Σ k = 1 m p 0 ( L f p 0. k p f p 0. k ) - - - ( 2 ) ]]>

式中，m_p0为F_p0中预想故障的总数，L_fp0.k为步骤8)中计算得到的F_p0中第k个预想故障发生后所导致的损失代价，p_fp0.k为F_p0中第k个预想故障发生的概率；

10)将F_s0与F_p的并集记为F_s，采用集群计算平台，基于最新生成的潮流数据、稳定数据和预防控制候选措施集，计及F_s中各个预想故障的L_f与L_r分别相等的约束条件，进行电网安全稳定预防控制措施的优化决策计算；如能够得到满足F_s中各个预想故障安全稳定的预防控制优化措施增量集ΔC_p，则进入步骤12)，否则设置计算标志I_c为1，进入步骤11)；

所述预防控制优化措施增量集是指选定的预防控制设备及其控制变量、控制方向和控制量增量的集合；

11)若F_p中预想故障的总数m_p大于设定的预想故障数且SL_fp与SL_fp0的比值大于设定值λ，则首先按预想故障在预想故障排序表中的排序，从F_p中剔除排序在后的int(m_p/2)个预想故障，得到更新后的F_p和m_p，再通过公式(3)计算出SL_fp，返回步骤10)；否则将C_p作为最终的在线预防控制优化措施集提交给调度运行人员，结束本方法：

SL f p = Σ j = 1 m p ( L f p . j p f p . j ) - - - ( 3 ) ]]>

式中，L_fp.j为步骤8)中计算得到的F_p中第j个预想故障发生后所导致的损失代价，p_fp.j为F_p中第j个预想故障发生的概率；

12)基于最新生成的电网运行断面潮流数据，结合当前运行方式下保护装置的运行状态及其模型和参数，通过公式(4)计算出F₀中各个预想故障发生后因预想故障直接触发的保护装置动作所导致的风险代价总和SL_r；通过公式(5)计算出F_p中各个预想故障发生后除了因预想故障直接触发的保护装置动作所导致的风险代价之外所产生的风险代价总和SL_frp，进入步骤13)；

SL r = Σ i = 1 n ( L r . i p f . i ) - - - ( 4 ) ]]>

SL f r p = Σ j = 1 m p [ ( L f p . j - L r p . j ) p f p . j ] - - - ( 5 ) ]]>

式中，L_r.i为在最新生成的电网运行断面潮流数据所对应的电网运行状态下F₀中第i个预想故障发生后因预想故障直接触发的保护装置动作所导致的损失代价，p_f.i为F₀中第i个预想故障发生的概率；L_fp.j为步骤8)中计算得到的F_p中第j个预想故障发生后所导致的损失代价，L_rp.j为步骤8)中计算得到的F_p中第j个预想故障发生后因预想故障直接触发的保护装置动作所导致的损失代价，p_fp.j为F_p中第j个预想故障发生的概率；

13)在最新生成的电网运行断面潮流数据和稳定数据的基础上，生成实施ΔC_p之后的潮流数据和稳定数据，用于更新最新生成的电网运行断面潮流数据和稳定数据，进入步骤14)；

14)基于最新生成的电网运行断面潮流数据，结合当前运行方式下保护装置的运行状态及其模型和参数，通过公式(6)计算出步骤13)得到的最新生成的电网运行断面潮流数据所对应的电网运行状态下F₀中各个预想故障发生后因预想故障直接触发的保护装置动作所导致的风险代价总和SL'_r，进入步骤15)；

SL r , = Σ i = 1 n ( L r . i , p f . i ) - - - ( 6 ) ]]>

式中，L'_r.i为在步骤13)得到的最新生成的电网运行断面潮流数据所对应的电网运行状态下F₀中第i个预想故障发生后因预想故障直接触发的保护装置动作所导致的损失代价，p_f.i为F₀中第i个预想故障发生的概率；

15)计算在电网当前运行方式下实施C_p和ΔC_p所付出的调整代价与步骤14)中计算得到的SL'_r两者之和，若该和小于等于在电网当前运行方式下仅实施C_p所付出的调整代价与步骤12)计算得到的SL_r和SL_frp三者之和，则进入步骤16)，否则，返回步骤11)；

16)根据ΔC_p中选定的各个预防控制措施分别对C_p作如下处理后，进入步骤17)：

若ΔC_p中选定的预防控制措施的预防控制设备及其控制变量在C_p中已存在，则将与该选定的预防控制措施的预防控制设备及其控制变量相关的C_p中的最终控制量与ΔC_p中的控制量增量之和用于更新C_p中的最终控制量，否则，将该选定的预防控制措施的预防控制设备及其控制变量、控制方向和控制量增量分别作为新增的预防控制优化措施的预防控制设备及其控制变量、控制方向和最终控制量添加到C_p中；

17)若I_c为0，进入步骤18)，否则，将C_p作为最终的在线预防控制优化措施集提交给调度运行人员，结束本方法；

18)根据ΔC_p中选定的各个预防控制措施分别对最新生成的预防控制候选措施集作如下处理后，返回步骤6)：

第一步，若ΔC_p中选定的预防控制措施的预防控制设备及其控制变量是首次选定的，则从最新生成的预防控制候选措施集中剔除与该预防控制措施的预防控制设备及其控制变量相同但控制方向不同的候选措施；

第二步，从最新生成的预防控制候选措施集中搜索到与ΔC_p中选定的预防控制措施的预防控制设备及其控制变量相同的候选措施，若该候选措施的可控量最大值与ΔC_p中选定的预防控制措施的控制量增量相等，则从最新生成的预防控制候选措施集中剔除该候选措施，否则，采用该候选措施的可控量最大值与ΔC_p中选定的预防控制措施的控制量增量二者的差值更新该候选措施的可控量最大值。

2.根据权利要求1所述的计及保护与安自控制代价的在线预防控制风险决策方法，其特征在于，所述步骤2)中的门槛值p_cr取《电力系统安全稳定导则》中第二道安全防线所防御的所有故障中发生概率的最小值的十分之一。

3.根据权利要求1所述的计及保护与安自控制代价的在线预防控制风险决策方法，其特征在于，所述步骤2)中的Δt_f设置为15分钟。

4.根据权利要求1所述的计及保护与安自控制代价的在线预防控制风险决策方法，其特征在于，所述步骤6)中的Δt_c设置为5分钟。

5.根据权利要求1所述的计及保护与安自控制代价的在线预防控制风险决策方法，其特征在于，所述步骤7)中的k_p设置为0.9。

6.根据权利要求1所述的计及保护与安自控制代价的在线预防控制风险决策方法，其特征在于，所述步骤11)中的预想故障数为大于1的奇数。

7.根据权利要求6所述的计及保护与安自控制代价的在线预防控制风险决策方法，其特征在于，所述步骤11)中的预想故障数设置为5。

8.根据权利要求1所述的计及保护与安自控制代价的在线预防控制风险决策方法，其特征在于，所述步骤11)中的设定值λ设置为0.2。