一种融合递推状态估计的动态电力系统异常检测方法技术方案

本发明专利技术属于动态电力系统态势感知技术领域，公开了一种融合递推状态估计的动态电力系统异常检测方法，其包括以下步骤：首先根据电力系统动态特性建立简化的系统节点电压的非线性递推模型，合理表征电力系统负荷动态变化对系统节点电压的影响；然后基于递推状态估计滤波算法实现基于非线性递推模型的系统节点电压动态估计，在此基础上，进一步构造系统节点电压残差随机矩阵；最后构造基于特征谱均值‑方差统计分析的系统残差动态性能指标，以有效反映动态电力系统异常对系统节点电压残差矩阵特征值分布的影响，进而根据自适应统计阈值判定，最终实现对动态电力系统的有效状态评估和异常检测。

A Dynamic Power System Anomaly Detection Method Based on Recursive State Estimation

【技术实现步骤摘要】
一种融合递推状态估计的动态电力系统异常检测方法
本专利技术属于动态电力系统态势感知
，尤其涉及一种融合递推状态估计的动态电力系统异常检测方法。
技术介绍
作为一种新兴技术，动态电力系统态势感知可以有效获取、理解和预测导致系统状况变化的关键因素，为复杂电力系统智能运维和决策提供依据，在电力调度、输配电管理等领域具有广阔的应用前景。它可以促进电网自动化系统的集成，及时发现系统的弱点和威胁，并显着提高电力系统的智能化水平，从而保证电力系统安全稳定运行。此外，动态电力系统态势感知还可以与有效的系统决策和支持模块相结合，从而提高电力系统的可靠性。因此，动态电力系统态势感知是实现未来智能电网不可或缺的一部分。近年来，随着电力系统规模的不断增大和大数据、人工智能等先进数据分析技术的发展，对动态电力系统态势感知的研究受到广泛关注。电力系统态势感知的时间尺度可以从瞬态过程的几毫秒、长时状态评估的几分钟到电网规划的几个月甚至几年变化。在现代电力系统中，以实时向量测量单元为基础的监测数据集，其关联的数据管理至关重要。基于PMU数据流的数据驱动方法被认为非常适合于电力系统态势感知应用。其中，基于随机矩阵理论分析框架，提出了一种数据驱动分析架构来感知高维复杂电力网格的稳态运行状态，构造了一种基于圆环定理的平均谱半径统计指标进行系统运行异常事件检测。在此基础上，大量研究集中在基于RMT的电力系统稳态状态评估中，提出了多种基于线性特征值统计量的改进指标体系来评估电力系统稳定性，并分析了关联统计量的统计相关性。然而，线性特征值统计量指标的异常检测性能依赖于节点功率变化灵敏度。因此，必须根据经验设置其异常检测阈值，并且一些初期异常可能无法有效察觉。此外，在动态电力系统的长时状态监测中，由于电力负荷的正常变化，特别是当智能电网存在分布式可再生能源时，通常无法保证系统节点电压恒定不变。因此，上述电力系统态势感知性能指标是不合适的，无法实现对动态电力系统正常和异常运行状态的有效区分。值得注意的是，有效的系统状态观测器能够很好地估计动态电力系统的运行特征参数，有效地消除原始PMU数据流中存在的复杂时空关联，从而获得系统正常运行状态下具有显著随机分布特征的系统残差矩阵。
技术实现思路
本专利技术的目的在于提出一种融合递推状态估计的动态电力系统异常检测方法，以有效地实现动态电力系统长时状态评估和异常检测。本专利技术为了实现上述目的，采用如下技术方案：一种融合递推状态估计的动态电力系统异常检测方法，包括如下步骤：s1.对复杂电力系统进行多节点电气等效，进而对由电气等效后的多节点构成的复杂电力系统，根据电力系统长时动态特性建立简化的系统节点电压的非线性递推模型；s2.基于递推状态估计滤波算法实现基于非线性递推模型的系统节点电压动态递推估计，得到系统节点电压原始残差矩阵；在此基础上，进一步构造随机扰动注入的系统节点电压残差随机矩阵；s3.构造基于特征谱均值-方差统计量的系统残差动态性能指标，进而根据构造的自适应统计判定阈值，实现对动态电力系统的有效状态评估和异常检测。优选地，步骤s1中，经过n节点电气等效的复杂电力系统可表示为：xk+1＝f(xk)+qk(1)zk+1＝h(xk+1)+rk+1(2)其中，xk为(2n-1)×1维电压幅值状态向量，包括n个系统节点电压幅值与相角；zk为(2n+2m)×1维系统量测向量，包括n个系统节点电压幅值与相角，以及m条线路有功功率潮流与无功功率潮流；f(x)为等效非线性过程方程，用以表征电压幅值的状态向量；h(x)为等效非线性过程方程，用以表征系统量测向量的函数关系；qk～N(0,Qk)为n×1维系统状态扰动向量；rk+1～N(0,Rk+1)为m×1维系统量测扰动向量；Qk和Rk+1分别为系统状态扰动向量和系统量测扰动向量对应的协方差矩阵。优选地，非线性递推模型的表达式为：xk+1|k＝Sk+bk(3)Sk＝αHxk|k+(1-αH)xk|k-1(4)bk＝βH(Sk-Sk-1)+(1-βH)bk-1(5)其中，xk+1|k表示k时刻对k+1时刻系统状态的预测值；xk|k表示k时刻当前系统状态的量测值；xk|k-1表示k-1时刻对k时刻系统状态的预测值；Sk、Sk-1、bk以及bk-1均为中间过程变量；αH和βH分别为经验拟合系数。优选地，步骤s2中，递推状态估计滤波算法包括无迹卡尔曼滤波算法或粒子滤波算法。优选地，步骤s2中，系统节点电压原始残差矩阵为系统节点电压幅值与系统节点电压幅值递推估计值的差值矩阵表示为：其中，xi,k为系统节点电压幅值矩阵的第i行、第k列元素；为系统节点电压幅值递推估计矩阵的第i行、第k列元素；为系统节点电压原始残差矩阵的第i行、第k列元素；差值矩阵的列数为w，表示系统节点电压原始残差矩阵的采样窗口大小。优选地，步骤s2中，随机扰动注入的系统节点电压残差随机矩阵Xk为：其中，e为满足高斯分布的0均值随机注入扰动，其方差为优选地，步骤s3中基于特征谱均值-方差统计量的系统残差动态性能指标SAk的表达式为：其中，λj,k为k时刻随机扰动注入的系统节点电压残差随机矩阵的协方差矩阵特征值；和分别为系统正常状态下节点j对应特征值序列λj的均值和标准差，该系统正常状态节点特征值序列由满足高斯分布的0均值随机注入扰动e的均值和标准差近似替代。优选地，步骤s3中，自适应统计判定阈值δ为：其中，其中，f为施加于单一节点电压的异常分量值；λn为实非相关Wishart矩阵w-1EET的最小特征值；E为(2n-1)×w维矩阵，其中每一列分别对应随机注入扰动向量e，其方差为表示异常分量值f的平方的数学期望；表示实非相关Wishart矩阵w-1EET的最小特征值的数学期望；为实非相关Wishart矩阵w-1EET的第i特征值的方差。优选地，自适应统计判定阈值δ为单一异常幅值下异常分量值f检测的充分条件，当系统中多个节点叠加异常分量大于异常分量值f时，自适应统计判定阈值δ均能够可靠检测。本专利技术具有如下优点：如上所述，本专利技术融合递推状态估计理论，在对动态电力系统运行时空关联状态进行有效估计基础上，基于随机矩阵特征值统计分析方法，在设计有效的态势感知性能指标的同时，给出了明确的异常检测控制限，具有清晰严谨的数学物理意义，能更加有效的量化动态电力系统潜在异常特性，提高动态电力系统状态评估和异常检测的有效性和合理性。附图说明图1为本专利技术实施例中融合递推状态估计的动态电力系统异常检测方法实施流程图。图2为本专利技术实施例中系统节点电压原始数据与残差生成数据的正态性检验对比图。图3为本专利技术实施例中融合递推状态估计的动态电力系统异常检测方法实施效果图。具体实施方式下面结合附图以及具体实施方式对本专利技术作进一步详细说明：如图1所示，一种融合递推状态估计的动态电力系统异常检测方法，包括如下步骤：I.对复杂电力系统进行多节点电气等效，进而对由电气等效后的多节点构成的复杂电力系统，根据电力系统长时动态特性建立简化的系统节点电压的非线性递推模型。其中，经过n节点电气等效的复杂电力系统可表示为：xk+1＝f(xk)+qk(1)zk+1＝h(xk+1)+rk+1(2)其中，xk为(2n-1)×1维电压幅值状态向量，包括n个系统节点电压幅值与相角，本文档来自技高网...

【技术保护点】
1.一种融合递推状态估计的动态电力系统异常检测方法，其特征在于，包括如下步骤：s1.对复杂电力系统进行多节点电气等效，进而对由电气等效后的多节点构成的复杂电力系统，根据电力系统长时动态特性建立简化的系统节点电压的非线性递推模型；s2.基于递推状态估计滤波算法实现基于非线性递推模型的系统节点电压动态递推估计，得到系统节点电压原始残差矩阵；在此基础上，进一步构造随机扰动注入的系统节点电压残差随机矩阵；s3.构造基于特征谱均值‑方差统计量的系统残差动态性能指标，进而根据构造的自适应统计判定阈值，实现对动态电力系统的有效状态评估和异常检测。

【技术特征摘要】
1.一种融合递推状态估计的动态电力系统异常检测方法，其特征在于，包括如下步骤：s1.对复杂电力系统进行多节点电气等效，进而对由电气等效后的多节点构成的复杂电力系统，根据电力系统长时动态特性建立简化的系统节点电压的非线性递推模型；s2.基于递推状态估计滤波算法实现基于非线性递推模型的系统节点电压动态递推估计，得到系统节点电压原始残差矩阵；在此基础上，进一步构造随机扰动注入的系统节点电压残差随机矩阵；s3.构造基于特征谱均值-方差统计量的系统残差动态性能指标，进而根据构造的自适应统计判定阈值，实现对动态电力系统的有效状态评估和异常检测。2.根据权利要求1所述的动态电力系统异常检测方法，其特征在于，所述步骤s1中，经过n节点电气等效的复杂电力系统可表示为：xk+1＝f(xk)+qk(1)zk+1＝h(xk+1)+rk+1(2)其中，xk为(2n-1)×1维电压幅值状态向量，包括n个系统节点电压幅值与相角；zk为(2n+2m)×1维系统量测向量，包括n个系统节点电压幅值与相角，以及m条线路有功功率潮流与无功功率潮流；f(x)为等效非线性过程方程，用以表征系统节点电压幅值的状态向量；h(x)为等效非线性过程方程，用以表征系统量测向量的函数关系；qk～N(0,Qk)为n×1维系统状态扰动向量；rk+1～N(0,Rk+1)为m×1维系统量测扰动向量；Qk和Rk+1分别为系统状态扰动向量和系统量测扰动向量对应的协方差矩阵。3.根据权利要求1所述的动态电力系统异常检测方法，其特征在于，所述步骤s1中，非线性递推模型的表达式为：xk+1|k＝Sk+bk(3)Sk＝αHxk|k+(1-αH)xk|k-1(4)bk＝βH(Sk-Sk-1)+(1-βH)bk-1(5)其中，xk+1|k表示k时刻对k+1时刻系统状态的预测值；xk|k表示k时刻当前系统状态的量测值；xk|k-1表示k-1时刻对k时刻系统状态的预测值；Sk、Sk-1、bk以及bk-1均为中间过程变量；αH和βH分别为经验拟合系数。4.根据权利要求1所...

【专利技术属性】
技术研发人员：彭飞，李立伟，王凯，张明，
申请(专利权)人：青岛大学，
类型：发明
国别省市：山东,37