一种考虑不等式约束的电力系统频率估计方法技术方案

本发明专利技术提出了一种考虑不等式约束的电力系统频率估计方法。该方法的步骤如下：步骤一：获取滤波目标电力信号；步骤二：确定滤波目标电力信号特征；步骤三：针对滤波目标信号，建立滤波模型，相应的状态方程、观测方程、状态不等式约束条件，以及离散滤波目标信号；步骤四：对滤波模型线性化后，设计等式约束条件下的扩展卡尔曼滤波器；步骤五：根据在线采样时刻的状态是否满足不等式条件，选择不同的扩展卡尔曼滤波器；步骤六：通过界面输出信号。频率估计对现代电力系统有着重要的意义，而其约束带来的一些数值计算错误会影响估计效果，因此在扩展卡尔曼滤波器设计中加入不等式约束条件，提高滤波的频率估计效率和精度。

【技术实现步骤摘要】
一种考虑不等式约束的电力系统频率估计方法
本专利技术属于电力系统领域，涉及考虑不等式约束的电力系统频率估计方法。
技术介绍
通常而言，在电力系统的参数估计中，频率估计是目前要解决的核心和难点，所以频率估计对于现代电力系统的运行有着重要意义。继电器采用基于频率的信号处理算法，从采样得到的电压信号或电流信号中，可估计系统状态信息。在电能质量监控领域，电力信号的电压幅值和相位检测皆是基于频率测量而进行的。在分布式发电系统并网过程中，根据频率测量值可分析供电设备与电网之间的动态能量平衡关系。但是电力系统中电力电子装置的广泛使用，可控硅转换装置和高频逆变器注入非正弦电流信号，导致电力信号受到谐波和噪声信号的破坏而失真。因此快速并准确对失真电压(电流)信号频率进行估计成为当前电力系统研究热点。卡尔曼滤波是一种高效率的递归滤波器，它能够从一系列的不完全包含噪声的测量中，估计动态系统的状态，采用时域状态空间的方法求解状态变量，是一种最优滤波器。广泛地应用于电力系统中信号处理、通讯、控制领域。在卡尔曼滤波的实际应用中，某些已知的状态信息经常被忽略，既容易造成数值计算错误，又极大影响估计的准确性，特别在电力系统频率估计时。于是将状态变量附加一定的等式或不等式约束条件在解决电力系统实际问题中具有重要的理论和实践意义。本专利技术提出了一种考虑不等式约束的电力系统频率估计方法，运用有效约束集合的思想，将不等式约束分类为有效约束和无效约束，指出了等式和不等式约束扩展卡尔曼滤波解性质的内在联系，提高了电力系统频率估计的计算效率和估计精度。
技术实现思路
本专利技术提出了一种考虑不等式约束的电力系统频率估计方法。此算法针对扩展卡尔曼滤波在电力系统实际应用中，某些信息经常被忽略，既容易造成数值计算错误，又极大影响估计得准确性等问题，提出了对状态变量附加一定的不等式约束条件的扩展卡尔曼滤波，并以此对电力系统的频率进行估计。为达到上述目的，本专利技术提供如下技术方案：一种考虑不等式约束的电力系统频率估计方法，包括以下步骤：步骤一：获取滤波目标电力信号；步骤二：确定滤波目标电力信号特征；步骤三：针对滤波目标信号，建立滤波模型，相应的状态方程、观测方程、状态不等式约束条件，以及离散滤波目标信号；步骤四：对滤波模型线性化后，设计等式条件下的扩展卡尔曼滤波器；步骤五：根据在线采样时刻的状态是否满足不等式条件，选择不同的扩展卡尔曼滤波器；步骤六：通过界面输出信号。进一步，步骤一具体包括，通过电压互感器获取电力系统中的失真信号，并将该失真信号调制为弱电信号制，从而得到滤波目标信号y。进一步，步骤二中，所述的滤波目标信号特征为，该信号为电力系统的所有正弦失真信号，包括三相或单相、电流或电压。进一步，步骤三中，具体包括如下步骤：首先，将连续非线性的电力系统失真信号转换为离散时间信号yk；之后，根据滤波目标信号的特征，采用连续三个采样时刻的电压等式关系构建滤波模型，包括状态方程和观测方程；在频率估计中，通过yk计算出然后计算出频率f：由数学知识可知|2cos2πfTs|≤2。然而在实际问题中，受噪声等因素影响，会出现状态向量中元素2cos2πfTs的绝对值大于2的现象(特别是在正弦波波峰和波谷)，f数值计算问题无解。因此状态向量一定要满足不等式其中F＝[100]，b＝2。需要注意的是，当谐波和衰减直流分量导致信号失真时，建模所采用的等式关系以及构建得到的模型都需改变。进一步，步骤四中，具体包括，对滤波模型线性化后，设计等式条件下的扩展卡尔曼滤波器；首先，进行一次无约束扩展卡尔滤波，得到状态预估然后结合状态向量不等式约束分类进行考虑。在电力系统中，状态变量之间的关系不能用确定的关系来表达，而是存在一定的不等式约束关系；大致可以将不等式约束分为两部分：有效约束和无效约束；当时，始终落在可行域中，去掉这类约束不会对解造成影响，称这类约束为无效约束；当时，那么不等约束问题随之变成一个有效等式约束的问题；于是我们得出结论，不等式约束卡尔曼滤波估值具有和等式约束卡尔曼滤波估值同样的统计性质；所以按照等式约束卡尔曼滤波的方法来处理，并设计等式扩展卡尔曼滤波器；当时，由于EKF实质是线性最优化估计，当超出线性不等式约束的边界时，视其为有效约束，即按上述第二种情况，等式约束卡尔曼滤波处理。进一步，步骤五中，对于无效约束，按照传统扩展卡尔曼滤波的公式，计算预估计误差协方差矩阵，计算卡尔曼增益矩阵，更新估计，计算更新后估计协方差矩阵；对于有效约束，采用时间更新，测量更新和约束更新这三步；接下来，两者都需要判断k与N的关系，当k小于N时(N是人为设定的扩展卡尔曼滤波结束时间)，重复以上的步骤，反之，则直接输出数据。进一步，步骤六中，将经过步骤1～5即不等式约束扩展卡尔曼滤波处理电力系统信号后，所得到的频率估计信号，通过用户界面反馈给电力系统工作人员。本专利技术提出了一种考虑不等式约束的电力系统频率估计方法。其优点在于：本专利技术方法有利于解决了不等式约束电力系统频率估计问题。在电力系统频率估计中，我们通过计算出然后计算出频率f。由数学知识知道|2cos2πfTs|≤2，然而在实际问题中，受噪声的因素的影响，得到的状态向量中元素大于2，频率估计时，f的数值计算出现无解的问题对于这个问题，目前学者们提出了两种方法。一种是忽略这一部分信息，缺点是某些信息被忽略，会造成的数值计算错误和影响估计准确性等问题。另一种是密集采样扩展卡尔曼滤波，虽然很大程度上减少这一部分信息的丢失，但是计算量过大，影响信号处理的实时性。本专利技术提出的一种考虑不等式约束的电力系统频率估计方法，完全避免了上述两种方法的缺点。本专利技术方法也具有传统卡尔曼滤波的优点，能从一系列的不完全包含噪声的测量中估计系统的状态，采用时域状态空间的方法求解状态变量。附图说明为了使本专利技术的目的、技术方案和有益效果更加清楚，本专利技术提供如下附图进行说明：图1为本专利技术所述的一种考虑不等式约束的电力系统频率估计方法流程示意图；图2为本专利技术所述的一种考虑不等式约束的电力系统频率估计方法算法流程图。具体实施方式下面将结合附图，对本专利技术的优选实施例进行详细的描述。图1为本专利技术的流程示意图，该电力系统失真信号参数估计方法包括五个步骤：一种考虑不等式约束的电力系统频率估计方法，包括以下步骤：步骤一：获取滤波目标电力信号；步骤二：确定滤波目标电力信号特征；步骤三：针对滤波目标信号，建立滤波模型，相应的状态方程、观测方程、状态不等式约束条件，以及离散滤波目标信号；步骤四：对滤波模型线性化后，设计等式条件下的扩展卡尔曼滤波器；步骤五：根据在线采样时刻的状态是否满足不等式条件，选择不同的扩展卡尔曼滤波器；步骤六：通过界面输出信号。1.由于本专利技术针对的是电力系统的失真信号，则步骤一中所述的获取滤波目标信号具体步骤为，通过电压互感器获取电力系统中的失真信号，并将该失真信号调制为弱电信号制，从而得到滤波目标信号yk。2.在本实施例中，步骤二中所述的滤波目标信号特征，首先需明确的一点便是该信号中混有附加噪声或谐波信号，本专利技术主要考虑电力系统信号混有零均值的高斯白噪声情况。其次，该滤波目标信号的特征还包括，该信号可以为单相或者三相非线性正弦波电压或电流信号。3.在本实施例中，步骤三采用的本文档来自技高网...

【技术保护点】
一种考虑不等式约束的电力系统频率估计方法，其特征在于，包括以下步骤：步骤一：获取滤波目标电力信号；步骤二：确定滤波目标电力信号特征；步骤三：针对滤波目标信号，建立滤波模型，相应的状态方程、观测方程、状态不等式约束条件，以及离散滤波目标信号；步骤四：对滤波模型线性化后，设计等式约束扩展卡尔曼滤波器；步骤五：根据在线采样时刻的状态是否满足不等式条件，选择不同的扩展卡尔曼滤波器；步骤六：通过界面输出信号。

【技术特征摘要】
1.一种考虑不等式约束的电力系统频率估计方法，其特征在于，包括以下步骤：步骤一：获取滤波目标电力信号；步骤二：确定滤波目标电力信号特征；步骤三：针对滤波目标信号，建立滤波模型，相应的状态方程、观测方程、状态不等式约束条件，以及离散滤波目标信号；步骤四：对滤波模型线性化后，设计等式约束扩展卡尔曼滤波器；步骤五：根据在线采样时刻的状态是否满足不等式条件，选择不同的扩展卡尔曼滤波器；步骤六：通过界面输出信号；步骤一，具体包括，通过电压互感器获取电力系统中的失真信号，并将该失真信号调制为弱电信号制，从而得到滤波目标信号y；步骤二中，所述的滤波目标信号特征为，该信号为电力系统的所有正弦失真信号，包括三相或单相、电流或电压的正弦失真信号；步骤三，具体包括如下步骤：首先，将连续非线性的电力系统失真信号转换为离散时间信号yk；之后，根据滤波目标信号的特征，采用连续三个采样时刻的电压等式关系构建滤波模型，包括状态方程和观测方程；在频率估计中，通过yk计算出然后计算出频率f：其中Ts为采样周期，和分别为第k-1时刻和k-2时刻的电压采样信号；由数学知识可知|2cos2πfTs|≤2；然而在实际问题中，受噪声等因素影响，会出现状态向量中元素2cos2πfTs的绝对值大于2的现象，f数值计算问题无解；因此状态向量一定要满足不等式其中F＝[100]，b＝2；需要注意的是，当谐波和衰减直流分量导致信号失真时，建模所采用的等式关系以及构建得到的模型都需改变；步骤四，具体包括，对滤波模...

【专利技术属性】
技术研发人员：魏善碧，柴毅，周展，陈淳，张利果，王诗年，
申请(专利权)人：重庆大学，
类型：发明
国别省市：重庆;85