关于智能电网防御注入恶意攻击的疑问（更新）

flyrathabu

, T+ a5 E9 ]8 @; s$ `
背景交代：

智能电网SCADA以及AMI对于恶意攻击的识别和防御是近年来在能源领域的一个重要话题，智能电网在控制上涉及到的东西与其说是电力方向上的一个研究，更确切的说是一个关于控制概念在电力系统上的一个应用。
8 v! t2 G2 _! ?  R# s) ~8 k# F
一个电力网络被表示成测量量----通过AMI测量得到的电网当中的电压幅值（近年来关于电流以及电压幅角度也有谈及），有偿功率以及无偿功率数值。通过功率与电压的函数关系，来估计各个节点的电压幅值以及角度（状态量）。从而了解整个网络的状态并进行恰当的调节控制。它们的关系可以表示为
, W& l% a% W% p. Z( {+ w9 @
, x3 k, T( F, a0 t! O- KX：状态量
7 r+ [0 {9 l% q1 G! A& pY：测量量
U：干扰与错误
# t8 S2 R2 E4 ^/ k2 z+ r( A% JX（Ableitung）=AX+BU
Y= CX+DU。。。。。。。。。。。。。。。。。。。。（1）
0 h! E; V/ p* K' q这里，A,B,C,D是线性化后的矩阵，U在仿真模拟中被看成输入量
8 {$ \6 s, ~  E, ~+ ]& P- T4 H
或者可以表示为  Y = HX+U。。。。。。。。。。。。。。。。。。。（2）
其中H是这个不连续函数h（X）的雅可比矩阵，h（X）可以堪看成Y和X的关系函数，H则可以看成（1）的传递矩阵
7 s) C9 U* x) W4 [5 X4 e
在智能化的电网中，测量量的重要性显得越来越重要。因为状态估计的算法从某种意义上已经相对比较完善，所以保证测量量的正确性和以及识别错误的状态量和人为恶意攻击设置的虚假状态量就显得尤为重要。传统方法上，一般是利用残差的方式来识别不可利用的状态量，这类识别方法基本上都是基于方程（2）的延伸：
" ^2 g% _6 \" C5 B6 Q. H- Q||X（estimation估计量）-X（test测量量）||<t。。。。。。。。。。。。。。。（3）
0 \6 @# v$ |3 _) g当估计量和测量量的模制小于给定的t阀值的时候，可以认为这个数据是可信的（这里（3）主要涉及的是发现不良数据，检测出哪个是不良数据不展开说明）。

更为关键的地方在于如何识别人为的恶意攻击，即人为恶意篡改后的测量量
0 `7 B6 i' X6 g. S, w9 {这时候将X表述为
X（test）=X（real）+a。。。。。。。。。。。。。。。（4）
a是错误信息
  G4 [# s: d; g
3 X7 O2 _7 z5 P. @2 |& m' u, a一般而言，人为的恶意攻击也可以通过（3）的方法来识别，但是近年来有人通过数学的方式，发现了一种新型的恶意攻击模式，这种方式对于传统电网的危害比较有限，但是对于强烈依赖于AMI的智能电网则危害极大。这种恶意攻击的关键在于
# i4 X+ l; e' S; A2 T: l8 da=Hc。。。。。。。。。。。。。。（5）
, X! j9 C2 ?/ z  F  ]c是一个非零的向量
6 Z" g2 f, P: W0 t# }通过这种方式合成的恶意攻击无法通过传统基于（3）的识别系统识别，从而导致恶意数据进入到了控制中心环节。（在这里恶意攻击还有其他方式，这里主要考虑的是对于AMI端的攻击）
# z6 y" f$ x; Q% v" P
更多的论文中对于这类攻击有相当数量的探索，在发现恶意攻击上有了相对成熟的一些办法，但是在检测具体属于哪一个AMI这个环节一直缺乏一种简单有效的手段。


主题：
% H6 j1 Q6 G  s; q正如我在背景交代开始所说的，智能网络状态估计的问题与其说是能源领域的一个问题，不如说是一个控制系统上的问题。我的疑问是，整个系统可以表述成
Y=HX+Hc+U
2 E0 \* `$ g+ k3 Q4 T# B% H8 c5 ]Y是测量量（其中Y和X是有部分重合的，具体的说Y包括电压幅值，有偿功率和无偿功率，而X则是部分电压幅值和电压角度）
+ L0 D2 j3 M' e/ HX是待估计量
5 n, }; ^1 A, V5 T( A, nHc是恶意攻击，假设为定值
U是干扰，可以简单的极为一个高斯干扰，服从正态分布
（表述为A B C D方式的传递函数因为存在攻击模式的不同有不同的函数表述，所以暂时不展开）
' U( ?! U. U* s
补充：
或者从另一个角度去考虑
恶意攻击首先确切的来自于改变测量量
4 s( O3 I; f9 w2 Q5 q( T. n$ |那么可以叫Y表述成
/ `) O4 l# q" UY（test）=Y(real)+a
" \$ J# s4 {' S. w, y同时a=Hc
这时候传递函数就成为了
- f2 t/ n, U1 G. L( d: D7 u6 KX（Ableitung）=AX+BU
: K% M1 B+ D8 Q/ y5 o. d: D& a& TY=CX+DU-Hc
是否有办法从控制角度得到Y的值
' B, H# f1 c4 D  f8 |8 B6 ?1 R
是否可以从控制论的角度，通过反馈的方式来中和掉其中的Hc甚至是U？如果有思路，能否请您给与一个大致的解决方向的说明。有请各位大牛给与指点。

看客

不明觉厉

看客

不明觉厉

punishment

有点不明白式3。状态变量如果能测量，那测量量和状态变量之间的区别是什么？
6 r- \3 |! ]1 S; r  v  U3 Z, d

flyrathabu

punishment 发表于 2015-3-18 23:45
有点不明白式3。状态变量如果能测量，那测量量和状态变量之间的区别是什么？
& g0 h: V$ K$ h: l  a4 J- U ...

测量量是有误差的，而通过状态估计得到的状态量可以看成最接近于真实数据的电压幅值和角度。