考虑变电站网络攻击成功概率的电力系统风险评估方法

1.本发明涉及电力信息物理系统技术领域，尤其是涉及考虑变电站网络攻击成功概率的电力系统风险评估方法。


背景技术：

2.变电站中许多控制和通信设备的快速部署给电力系统带来巨大的网络风险，因此有必要通过考虑网络层变电站的脆弱性来扩展经典的n-k故障分析方法，以增强电力系统的安全性能；
3.双层优化作为一种流行的评估策略，用于识别将导致故障最严重负荷损失的关键线路，即传统的n-k模型，此外学者还提出了一种基于概率的n-k模型来识别具有最大风险的k条线路，其中假定每条线路停运概率都是独立的，在之前所有研究中，还探讨利用概率模型基于风险的安全评估有效性；
4.与可能直接摧毁电网输电线路的自然灾害或物理攻击不同，发生的网络攻击分为两个阶段，攻击者首先侵入变电站网络空间；然后发送跳闸命令，断开相关线路，电力系统造成严重影响，因此受到变电站网络攻击的线路停运之间存在明显的相关性，在网络攻击背景下，监控与数据采集(scada)系统和智能电子设备(ied)是网络攻击者最具吸引力的目标，在已有研究中，建立仅考虑scada系统上的网络攻击双层模型，但针对ied的网络攻击并没有反映在任何双层优化中；
5.在这种情况下，考虑网络攻击ied成功概率的双层优化模型，可以很好反映电力系统面临的风险指标，反映系统的各网络物理连接，以识别易受攻击的变电站及其相关线路。


技术实现要素：

6.本发明的目的在于提供考虑变电站网络攻击成功概率的电力系统风险评估方法。
7.本发明的目的可以通过以下技术方案实现：
8.考虑变电站网络攻击成功概率的电力系统风险评估方法，该方法具体包括下述步骤：
9.s1.变电站网络攻击成功概率的风险评估双层优化模型的建立：对电力信息物理融合系统进行攻击模拟，同时考虑信息网优化调度功能在潮流重新分配中的作用以及网络攻击变电站的成功概率，建立双层优化模型，计算电力系统中风险最大值；
10.s2.运用强对偶理论对双层优化模型进行转换：基于s1所建立的双层优化模型，根据强对偶理论，将双层优化模型进行转换成单层混合整数线性模型；
11.s3.根据转换后的单层混合整数线性模型，运用求解器进行求解，计算并筛选出风险数最大的变电站以及其连接线路。
12.进一步的，所述步骤s1的具体包括：
13.s11.变电站网络模型的建立：根据防御资源的分配和已识别的漏洞信息，运用贝叶斯网络模块的攻击图，全面完整的表示变电站的网络模型，其中每个网络组件的漏洞是
随机生成的。
14.s12.变电站成功入侵概率：使用贝叶斯网络计算变电站成功入侵概率，变电站的成功入侵概率随防御资源的分配和已识别的漏洞而变化；
15.s13.确定最大化网络风险值的关键变电站及其相关线路：基于电力信息物理交互作用下的故障模型，计算线路受到攻击后电力系统的潮流转移，最终得出物理传输系统上产生的负荷损失，通过乘以变电站攻击成功概率来计算网络风险值。
16.s14.调度层的负荷重新分配模拟：调度中心接收到根据经过状态估计处理后的受攻击量测数据进行电力系统优化调度，基于最优直流潮流模型对物理系统的发电机出力和负荷进行调整，尽可能地减少负荷损失，以模拟调度中心检测到攻击行为后的优化控制行为，系统到达新的运行稳态；
17.s15.基于风险指标的双层优化模型的建立：建立的双层优化模型分为上层模型和下层模型，上层模型的目标是确定关键变电站及其相关线路，以最大化网络风险，下层模型引入了直流最优潮流，其目标是通过重新分配发电量和甩负荷最小化电力系统总的负载损耗。
18.进一步的，步骤s15中的双层优化模型的建立，是通过将变电站的攻击成功概率与负荷损失量相乘，得到网络风险值，其模型的上层目标函数为：
[0019][0020]
式中，sipn、分别表示第n个变电站的成功入侵概率、电力系统总负载损耗量，νn为二进制变量，若νn＝1，则表示成功入侵变电站的概率为sipn，反之则νn＝0，具体的上层约束条件如下：
[0021][0022][0023][0024][0025][0026]
式中可见，当给定常数k1和k2时，通过上层约束可计算出最大化网络风险值相关的关键变电站及其相关线路，常数k1和k2分别表示在攻击资源有效的情况下，允许攻击变电站最大个数和允许断开线路的最大数值，约束(1b)-(1c)用于制定应用于变电站和输电线路的攻击预算，该预算可由系统操作员确定以进行模拟，约束(1c)的左侧确保至少k1条线路跳闸，因为至少一条连接到受攻击变电站的线路应跳闸，约束条件(1d)表明，只有在线路两端的变电站不受影响的情况下，线路才是安全的，约束(1f)-(1g)确保了攻击者是具有高智商的情况下，如果一个变电站受到他的攻击，至少一条连接到变电站的线路应该跳闸。
[0027]
进一步的，步骤s15中模型的下层目标函数为：
[0028][0029]
其代表总负载损耗量的最小值。具体的下层约束条件如下：
[0030][0031]
kl
·
f＝kp
·
p-kd
·
(d-δd)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1i)
[0032]
p
min
≤p≤p
max
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1j)
[0033]-f
max
≤f≤f
max
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1k)
[0034]
0≤δd≤d
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1l)
[0035]
式中可见，下层约束(1h)-(1l)引入了直流最优潮流，调度人员接收到被攻击后的状态估计之后，通过重新分配发电量和甩负荷，其目标是使总负载损耗最小化。
[0036]
进一步的，所述步骤s2是基于步骤s1所建立的双层优化模型，运用强对偶理论对其模型从双层混合整数模型转变为单层混合整数线性模型，有利于在求解器中进行求解。
[0037]
进一步的，所述步骤s3是基于步骤s2所得到的单层混合整数线性模型，求解其风险指标最大值，并得到与之结果相关的被成功入侵的变电站以及输电线路，实现了最大风险辨识以及攻击路径的预测。
[0038]
本发明的有益效果如下：
[0039]
1、本发明在电力信息物理交互融合的背景下，考虑网络攻击变电站对电力系统运行的影响，建立变电站攻击概率且基于风险性指标的双层优化模型；
[0040]
2、本发明考虑到变电站网络攻击的线路停运之间存在明显的相关性，所提指标能够更好的显示变电站及其输电线路的风险性，可以对攻击者使用不同攻击资源的情况下进行准确分析，能够从攻击者角度层面出发，对潜在的攻击行为进行预测；
[0041]
3、本发明不仅可以对拥有不同变电站攻击成功概率的现实事例进行建模，还能对传统的n-k故障进行更全面的分析；
[0042]
4、本发明有利于提升现代电力系统主动防御的能力，及时准确地对风险值较大的攻击路径进行预知，系统运营商可以进一步了解电力系统的脆弱性分析，能够指导电力系统管理者对关键变电站防护，并制定全面的保护策略。
附图说明
[0043]
下面结合附图对本发明作进一步的说明。
[0044]
图1为本实施例用于ieee 24节点电力信息物理系统的最大风险指标计算实现流程图。
[0045]
图2为本实施例在具体应用中的ieee 24节点输电线路拓扑图。
具体实施方式
[0046]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于
本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。
[0047]
请参阅图1所示，本发明为如图1所示，本实施例用于电力信息物理系统的攻击模拟方法的步骤包括：
[0048]
s1.变电站网络攻击成功概率的风险评估双层优化模型的建立：对电力信息物理融合系统进行攻击模拟，同时考虑信息网优化调度功能在潮流重新分配中的作用以及网络攻击变电站的成功概率，建立双层优化模型，计算电力系统中风险值最大值；
[0049]
s2.运用强对偶理论对双层优化模型进行转换：基于s1所建立的双层优化模型，根据强对偶理论，将双层优化模型进行转换成单层混合整数线性模型；
[0050]
s3.根据转换后的单层混合整数线性模型，运用求解器进行求解，计算并筛选出风险数最大的变电站以及其连接线路。
[0051]
本实施步骤s1所述的建立变电站网络攻击成功概率的风险评估双层优化模型具体包括：
[0052]
s11.变电站网络模型的建立：根据防御资源的分配和已识别的漏洞信息，运用贝叶斯网络模块的攻击图，全面完整的表示变电站的网络模型，其中每个网络组件的漏洞是随机生成的；
[0053]
s12.变电站成功入侵概率：使用贝叶斯网络计算变电站成功入侵概率，变电站的成功入侵概率随防御资源的分配和已识别的漏洞而变化；
[0054]
s13.确定最大化网络风险的关键变电站及其相关线路：基于电力信息物理交互作用下的故障模型，计算线路受到攻击后电力系统的潮流转移，最终得出物理传输系统上产生的负荷损失，通过乘以变电站攻击成功概率来计算网络风险；
[0055]
s14.调度层的负荷重新分配模拟：调度中心接收到根据经过状态估计处理后的受攻击量测数据进行电力系统优化调度，基于最优直流潮流模型对物理系统的发电机出力和负荷进行调整，尽可能地减少负荷损失，以模拟调度中心检测到攻击行为后的优化控制行为，系统到达新的运行稳态；
[0056]
s15.基于风险指标的双层优化模型的建立：建立的双层优化模型分为上层模型和下层模型，上层模型的目标是确定关键变电站及其相关线路，以最大化网络风险，下层模型引入了直流最优潮流，其目标是通过重新分配发电量和甩负荷最小化电力系统总的负载损耗；
[0057]
步骤s15中的双层优化模型的建立，是通过将变电站的攻击成功概率与负荷损失量相乘，得到网络风险值，其模型的上层目标函数为：
[0058][0059]
式中，sipn、分别表示第n个变电站的成功入侵概率、电力系统总负载损耗量，νn为二进制变量，若νn＝1，则表示成功入侵变电站的概率为sipn，反之则νn＝0，具体的上层约束条件如下：
[0060]
[0061][0062][0063][0064][0065]
式中可见，当给定常数k1和k2时，通过上层约束可计算出最大化网络风险值相关的关键变电站及其相关线路，常数k1和k2分别表示在攻击资源有效的情况下，允许攻击变电站最大个数和允许断开线路的最大数值，约束(1b)-(1c)用于制定应用于变电站和输电线路的攻击预算，该预算可由系统操作员确定以进行模拟，约束(1c)的左侧确保至少k1条线路跳闸，因为至少一条连接到受攻击变电站的线路应跳闸，约束条件(1d)表明，只有在线路两端的变电站不受影响的情况下，线路才是安全的，约束(1f)-(1g)确保了攻击者是具有高智商的情况下，如果一个变电站受到他的攻击，至少一条连接到变电站的线路应该跳闸。
[0066]
步骤s15中模型的下层目标函数为：
[0067][0068]
其代表总负载损耗量的最小值，具体的下层约束条件如下：
[0069][0070]
kl
·
f＝kp
·
p-kd
·
(d-δd)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1i)
[0071]
p
min
≤p≤p
max
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1j)
[0072]-f
max
≤f≤f
max
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1k)
[0073]
0≤δd≤d
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1l)
[0074]
式中，下层约束(1h)-(1l)引入了直流最优潮流，调度人员接收到被攻击后的状态估计之后，通过重新分配发电量和甩负荷，其目标是使总负载损耗最小化。
[0075]
本实施步骤s2是基于步骤s1所建立的双层优化模型，运用强对偶理论对其模型从双层混合整数模型转变为单层混合整数线性模型，有利于在求解器中进行求解。
[0076]
本实施步骤s3是基于步骤s2所得到的单层混合整数线性模型，求解其风险指标最大值，并得到与之结果相关的被成功入侵的变电站以及输电线路，实现了最大风险辨识以及攻击路径的预测；
[0077]
为验证本发明的有效性，以ieee24节点电力信息物理系统为例，采用本发明方法进行脆弱相关性的分析，具体该系统包含24个电力厂站节点，38条输电线路，17个带负荷节点，12个发电机，利用python3.7为仿真工具进行仿真，得到风险指标最大值，包括被攻击的变电站及其输电线路，为了便于分析与比较，线路攻击预算k2固定为5，我们依次将变电站攻击预算k1增加到5的计算模型，通过计算经典pn5模型进行比较，结果如表1所示：
[0078]
为了验证本发明还能对传统的n-k故障进行准确的分析，我们将每个变电站攻击
的成功概率sipn设为1，线路攻击预算k2固定为5，依次将变电站攻击预算k1增加到5的计算模型，通过计算经典n-k模型进行比较，结果如表2所示：
[0079]
表1 ieee 24节点系统的网络风险变化
[0080]
k1目标变电站目标线路网络风险11318,20,220.566029,108,9,10,12,130.500436,9,105,8,9,12,130.228044,6,9,105,4,9,12,130.093559,12,13,15,1615,17,18,23,270.0261pn5/4,8,23,27,290.1797
[0081]
表2 ieee 24节点系统n-k故障计算结果
[0082]
k1目标变电站目标线路负荷损失(p.u.)112,35.560021,22,3,4,59.10031,2,82,3,4,5,1112.60041,2,3,72,3,4,5,1112.60052,3,5,6,72,3,4,5,1112.600n-5/2,3,4,5,1112.600
[0083]
表1中，随着变电站攻击预算的增加，得出了ieee 24总线系统具有最大网络风险的攻击场景，由于目标函数中的联合概率减少超过了目标函数中的负荷损失，网络风险从0.566降至0.0261，这与入侵多个变电站更加困难的客观事实一致。
[0084]
计算pn5模型以筛选停运线路进行比较，结果列在表1的最后一行，使用线路列表，通过本发明计算网络风险值，攻击者将跳闸与pn5模型筛选出的线路相同的线路，pn5的相关风险低于当k1＝1,2,3的情况，这表明pn5模型在此情况下无法识别网络风险最大的场景；
[0085]
表2中，如果当k1小于3时，所提方法的计算结果与经典的n-5模型不一致，也即在此条件下，经典的n-5模型不能完全计算出所需断开的输电线路，当k1大于等于3时，所提方法与经典的n-5模型计算结果一致。
[0086]
以上内容仅仅是对本发明结构所作的举例和说明，所属本技术领域的技术人员对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代，只要不偏离发明的结构或者超越本权利要求书所定义的范围，均应属于本发明的保护范围。