针对智能电网中虚假数据注入攻击的数据修复方法技术

本发明专利技术涉及一种针对智能电网中虚假数据注入攻击的数据修复方法，从观测数据中分离出虚假数据注入攻击的修改，在分离过程中，尽可能保证真实测量数据中的合理的测量值扰动不被误判为攻击造成的修改。通过该方法可比较准确地还原被虚假数据注入攻击所篡改的智能电网观测数据。

【技术实现步骤摘要】

本专利技术涉及机器学习领域，特别是针对智能电网中虚假数据注入攻击的数据修复方法。
技术介绍
在智能电网中，管理者通过散布在整个智能电网网络中的传感器和计量表来获得对智能电网的观测数据，并根据这些观测数据所反映的智能电网状态采取对应策略，可有效地监控和管理电能生产、传输和使用等环节。虚假数据注入是一种新型的针对智能电网的攻击方式，意在篡改智能电网中部分被攻击的传感器和计量表的测量数据值，使得管理者所获得的部分观测数据为虚假数据，从而影响其对智能电网真实状态的判断及所采取的对应策略。采用预防的策略来防止虚假数据注入攻击往往具有一定缺陷。例如，由于智能电网中具有大量的传感器和计量表，因此对所有的传感器和计量表进行保护来阻止虚假数据注入攻击的代价较大；而仅选择其中重要传感器和计量表进行保护的方式操作起来具有一定难度，因为在一定网络拓扑结构中“重要传感器和计量表”的定义往往模糊，同时智能电网的网络拓扑结构也非一成不变。为了保证智能电网的安全性和鲁棒性，如何自动修复虚假数据注入攻击所篡改的数据、还原真实的测量结果变得尤为重要。同时，在修复数据时，还应尽量避免将观测数据真实扰动(源自智能电网状态的真实扰动和观测设备受到环境影响的合理测量误差)也看作虚假数据注入攻击的结果。否则，若将这些真实扰动当作虚假数据去除，则修复得到的观测值与真实值之间将存在偏差，从而影响数据修复结果的准确性。
技术实现思路
本专利技术提出一种针对智能电网中虚假数据注入攻击的数据修复方法，使电网观测数
据的真实扰动在数据修复时能得以更好保留，提高数据修复准确性。一般地，在一个遭受虚假数据注入攻击的电网中，某一个时间点k上的电网所有m个观测数据形成一个观测值向量zk＝(zk1,zk2,...,zkm)T，一个时间段内t个连续时间点上的观测值向量形成一个m行t列的观测值矩阵Za＝(z1,z2,...,zt)。令矩阵Zo表示若电网在该时间段上不受虚假数据注入攻击时的真实观测值矩阵，矩阵A表示在遭受攻击时各观测数据在各时间点上被注入的数据偏移量(即被注入的虚假数据)，则有Za＝Zo+A。本专利技术旨在准确地修复带虚假数据的观测值矩阵Za，即将注入的数据偏移量矩阵A由Za中分离开来，从而尽可能准确地还原真实观测值矩阵Zo。本专利技术所提出的数据修复方法源自对矩阵Zo和矩阵A相关特性的观察。具体来说，由于在一定时段内智能电网状态相对比较稳定，因此各真实测量数据在各时间点上变化不大，从而使得真实观测值矩阵Zo具有低秩的特点。此外，由于智能电网中的测量设备在物理距离上相距较远，攻击者在某一时间点上只可能对局部区域的少量传感器和计量表进行虚假数据注入攻击，且为了尽量避免被发现，攻击者一般不会长时间持续攻击，也即只有少量观测数据在部分时间节点上被注入了偏移量，且偏移量往往不会过于明显，从而使得矩阵A具有稀疏性的特点，即大部分矩阵元素为0，且非0元素所代表的偏移量一般不会过大。值得指出的是，虽然真实观测值矩阵Zo具有低秩的特点，但在修复数据时，一味最小化Zo的秩可能会导致Zo中观测数据的真实扰动也被当作被注入的偏移量而划归矩阵A中，从而降低数据修复的准确性。为避免这种情况，本专利技术提出在Zo的秩不大于给定上限r的情况下，通过寻找一个最优的Zo来最小化||A||1(即使A中所有元素绝对值之和最小)从而拟合矩阵A的稀疏性特点。为使矩阵Zo的秩不超过r，可将矩阵Zo分解为m行r列且各列两两正交的矩阵U和r行t列的矩阵W，即Zo＝UW。由于矩阵U、W的秩均不大于r，则还原得到的矩阵Zo的秩也将不大于r。通过以上矩阵分解，问题可转化为在限制条件UW+A-Za＝0下，寻找一个最优的列正交矩阵U及对应的矩阵
W以最小化||A||1，即求解minU,W||A||1，s.t.UW+A-Za＝0的最优化问题。本专利技术按下述方案解决该问题：步骤1，初始化矩阵U为一个随机的m行r列且各列两两正交的矩阵(r≤t)。步骤2，固定当前U的值不变，构造增广拉格朗日函数L(W,A,Y,ρ)，公式如下：其中，矩阵Y为对偶矩阵，参数ρ为惩罚参数，<·,·>表示两个矩阵的内积，||·||p表示矩阵的Lp-范式(即如果矩阵X＝[x1,x2,…,xn]T，则||X||p＝(|x1|p+|x2|p+…+|xn|p)1/p)。步骤3，初始化矩阵Y、矩阵A、以及惩罚因子ρ，令Y＝0,A＝0,ρ>0。步骤4，迭代更新W、A、Y、ρ直至W收敛，其中第i+1次迭代计算过程如下：ρ＝αρ 公式5其中，函数S1/ρ定义为S1/ρ(x)＝sgn(x)*max(|x|-1/ρ,0)(当x>0时，sgn(x)＝1；当x＝0时，sgn(x)＝0；当x<0时，sgn(x)＝-1)，0＜α＜1.8。步骤5，固定当前W、A、Y、ρ的值不变，采用梯度下降法更新矩阵U直至收敛。步骤6，如果第一次运行到本步骤，或者当前UW的值与上一次运行到本步骤时UW的值相比未能收敛，则返回步骤2；否则，返回真实观测值矩阵Zo＝UW。在步骤5中，采用梯度下降法更新矩阵U的过程如下：步骤5.1，对公式1中函数L求关于U的偏导LU：LU＝YW-2(Za-UW-A)WT 公式6步骤5.2，计算函数L在当前U上的梯度▽L：▽L＝LU-UUTLU 公式7步骤5.3，把U的值沿-▽L(梯度下降)方向从Ui更新为Ui+1：Ui+1＝UiVcos(Σηi)VT+Psin(Σηi)VT 公式8其中，PΣV是-▽L的一个奇异值分解；ηi是第i次更新U时梯度下降的步长，为保证梯度下降算法的快速收敛，该步长满足如下条件：且附图说明图1为本专利技术的总体框架图。具体实施方式下面通过实施例，并结合附图，对本专利技术的技术方案作进一步具体的说明。本专利技术提供的是一种针对智能电网中虚假数据注入攻击的数据修复方法，对给定一段时间内t个连续时间点m个观测项上的观测值向量形成的一个m行t列的观测值矩阵Za＝(z1,z2,...,zt)，求一个真实观测值矩阵Zo＝UW，使得Zo的秩不超过一个给定的上限r，UW+A-Za＝0，且||A||1取得最小值。请见图1，具体是采用包括以下步骤的方法：步骤1，初始化矩阵U为一个随机的m行r列且各列两两正交的矩阵(r≤t)。步骤2，固定当前U的值不变，构造增广拉格朗日函数L(W,A,Y,ρ)，公式如下：L(W,A,Y,ρ)=||A||1+<Y,UW+A-Za>+ρ2||UW+A-Za||22]]>其中，矩阵Y为对偶矩阵，参数ρ为惩罚参数，<·,·>表示两个矩阵的内积，||·||p表示矩阵的Lp-范式(即如果矩阵X＝[x1,x2,…,xn]T，则||X||p＝(|x1|p+|x2|p+…+|xn|p)1/p)。步骤3，初始化矩阵Y、矩阵A、以及惩罚因子ρ，令Y＝0,A＝0,ρ>0。步骤4，迭代更新W、A、Y、ρ直至W收敛，其中第i+1次迭代计算过程如下：Yi+1=argminYL(Wi+1,Ai+1,Yi)=Yi+ρ(UWi+1+Ai+1-Za)]]>ρ＝αρ其中，函数S1/ρ定义为S1/ρ(x)＝sgn(x)*max(|x|-1/ρ,0)(当x>0时，sgn本文档来自技高网...

【技术保护点】
针对智能电网中虚假数据注入攻击的数据修复方法，其特征在于，需要从对电网状态的观测值矩阵Za＝Zo+A中将虚假数据矩阵A分离出来以还原真实观测值矩阵Zo，其中，Za、Zo、A皆为m行t列矩阵；具体方法是：在保证Zo的秩不大于给定上限r的情况下，通过寻找一个最优的Zo来最小化||A||1以拟合矩阵虚假数据注入攻击的稀疏性特点；其中，为使矩阵Zo的秩不超过r，可将矩阵Zo分解为m行r列且各列两两正交的矩阵U和r行t列的矩阵W，即Zo＝UW，由于矩阵U、W的秩均不大于r，则还原得到的矩阵Zo的秩也将不大于r；也就是在限制条件UW+A‑Za＝0下，寻找一个最优的列正交矩阵U以及矩阵W以最小化||A||1，即获取最优化问题minU,W||A||1，s.t.UW+A‑Za＝0的最优解，具体求解步骤包括：步骤1，初始化矩阵U为一个随机的m行r列且各列两两正交的矩阵(r≤t)；步骤2，固定当前U的值不变，构造增广拉格朗日函数L(W,A,Y,ρ)，公式如下：L(W,A,Y,ρ)=||A||1+<Y,UW+A-Za>+ρ2||UW+A-Za||22]]>其中，矩阵Y为对偶矩阵，参数ρ为惩罚参数，<·,·>表示两个矩阵的内积，||·||p表示矩阵的Lp‑范式(即如果矩阵X＝[x1,x2,…,xn]T，则||X||p＝(|x1|p+|x2|p+…+|xn|p)1/p)；步骤3，初始化矩阵Y、矩阵A、以及惩罚因子ρ，令Y＝0,A＝0,ρ>0；步骤4，迭代更新W、A、Y、ρ直至W收敛，其中第i+1次迭代计算过程如下：Wi+1=argminLW(Wi,Ai,Ai)=(UTU)-1UT(Za-Ai-1ρAi)]]>Ai+1=argminAL(Wi+1,Ai,Yi)=S1/ρ(Za-UWi+1-1ρYi)]]>Yi+1=argminYL(Wi+1,Ai+1,Yi)=Yi+ρ(UWi+1+Ai+1-Za)]]>ρ＝αρ其中，函数S1/ρ定义为S1/ρ(x)＝sgn(x)*max(|x|‑1/ρ,0)(当x>0时，sgn(x)＝1；当x＝0时，sgn(x)＝0；当x<0时，sgn(x)＝‑1)，0＜α＜1.8；步骤5，固定当前W、A、Y、ρ的值不变，采用梯度下降法更新矩阵U直至收敛；步骤6，如果第一次运行到本步骤，或者当前UW的值与上一次运行到本步骤时UW的值相比未能收敛，则返回步骤2；否则，返回真实观测值矩阵Zo＝UW。...

【技术特征摘要】
1.针对智能电网中虚假数据注入攻击的数据修复方法，其特征在于，需要从对电网状态的观测值矩阵Za＝Zo+A中将虚假数据矩阵A分离出来以还原真实观测值矩阵Zo，其中，Za、Zo、A皆为m行t列矩阵；具体方法是：在保证Zo的秩不大于给定上限r的情况下，通过寻找一个最优的Zo来最小化||A||1以拟合矩阵虚假数据注入攻击的稀疏性特点；其中，为使矩阵Zo的秩不超过r，可将矩阵Zo分解为m行r列且各列两两正交的矩阵U和r行t列的矩阵W，即Zo＝UW，由于矩阵U、W的秩均不大于r，则还原得到的矩阵Zo的秩也将不大于r；也就是在限制条件UW+A-Za＝0下，寻找一个最优的列正交矩阵U以及矩阵W以最小化||A||1，即获取最优化问题minU,W||A||1，s.t.UW+A-Za＝0的最优解，具体求解步骤包括：步骤1，初始化矩阵U为一个随机的m行r列且各列两两正交的矩阵(r≤t)；步骤2，固定当前U的值不变，构造增广拉格朗日函数L(W,A,Y,ρ)，公式如下：L(W,A,Y,ρ)=||A||1+<Y,UW+A-Za>+ρ2||UW+A-Za||22]]>其中，矩阵Y为对偶矩阵，参数ρ为惩罚参数，<·,·>表示两个矩阵的内积，||·||p表示矩阵的Lp-范式(即如果矩阵X＝[x1,x2,…,xn]T，则||X||p＝(|x1|p+|x2|p+…+|xn|p)1/p)；步骤3，初始化矩阵Y、矩阵A、以及惩罚因子ρ，令Y＝0,A＝0,ρ>0；步骤4，迭代更新W、A、Y、ρ直至W收敛，其中第i...

【专利技术属性】
技术研发人员：黄浩，李宗鹏，王梁，王淞，
申请(专利权)人：武汉大学，
类型：发明
国别省市：湖北;42