一种基于矢量匹配的电力系统次超同步谐波检测方法技术方案

本发明专利技术公开了一种基于矢量匹配的电力系统次超同步谐波检测方法，包括如下步骤：从电流序列x(n)中滤除工频分量得到序列y(n)；然后对序列y(n)的次同步和超同步频带范围内的频谱进行矢量匹配拟合，结合极值判别初步次超同步谐波；在预设的时间窗口内对初步次超同步谐波集合进行聚类，进而识别真实次超同步谐波。本发明专利技术充分借鉴电力系统间谐波分析范畴内的傅里叶理论和机器智能理论，实现在复杂时变运行环境下电网中次超同步谐波的准实时监测。在新能源并网后，大量电力电子控制装置并存的背景下优势更为突出。本发明专利技术对及时发现并切除电网中影响系统稳定性的谐波源具有重要指导意义。

【技术实现步骤摘要】
一种基于矢量匹配的电力系统次超同步谐波检测方法
本专利技术涉及电力系统动态振荡监测与分析领域，更具体地，涉及种基于矢量匹配的电力系统次超同步谐波检测方法。
技术介绍
随着风电、光伏等新能源的大规模开发和利用，以常规火电机组为主的电力系统正逐步转变为新能源电力系统。电力电子技术在增强新能源电力系统灵活性的同时，也给系统带来了新的扰动和隐患。与传统电网中同步、异步概念不同，系统中的电力电子装置诱发次超同步谐波后，可能仍会挂网运行并持续威胁机网，严重影响机网的安全运行。为了监测新能源电力系统的次超同步振荡现象，高精度量测装置(如次同步监测、稳控及PMU等)被部署于关键风电场和变电站内。它们通过测量线路的电压、电流和功率响应，并从滤波信号中获取次超同步信号，进而与整定值比较做出告警或动作保护。滤波信号所获得的次超同步信号是一个频带内的反映，并不能真实反映次超同步信号随时间变化的频率、幅值波动的现象。
技术实现思路
本专利技术克服了上述现有的关于次超同步谐振测量的技术缺陷，提供了一种新的基于矢量匹配的电力系统次超同步谐波检测方法，满足了电力系统上动态监测的准实时性要求。为解决上述技术问题，本专利技术的技术方案如下：一种基于矢量匹配的电力系统次超同步谐波检测方法，包括如下步骤：S1．从录波装置中获取电流采样序列x(n)，从x(n)中滤除工频分量得到序列y(n)，n为时域采样序列号；S2．在时间窗口TA内，将y(n)进行快速傅里叶变换FFT得到序列Y(k)，k为频域离散序列号；TA为预设值；S3．将次超同步谐波对应频段内的Y(k)进行矢量匹配，得到极点pl，通过pl中的共轭极点对得到对应的谐波频率f0；S4．若f0对应的点是Y(k)的极值点，且f0对应的点的幅值A1≥A0，则f0认定为初步次超同步谐波频率f1，执行S5，A1认定为初步次超同步谐波幅值A2；A0为预设值；若f0不符合要求，则执行S6；S5．将初步次超同步谐波频率f1和初步次超同步谐波幅值A2进行记录，得到初步次超同步谐波集合，所述的初步次超同步谐波集合中的每个元素包括三个维度：时间t，频率f1和幅值A2；S6．在时间窗口TB每相隔步长TC就移动时间窗口TA，若TA没有接触TB的结束时间，则执行S2；若TA已经接触TB的结束时间，则执行S7；所述的TB和TC为预设值；S7．对时间窗口TB内的初步次超同步谐波频率集合通过聚类算法得到聚类中心Cm(m=1，2，…，M，M为聚类中心的数目)，频率介于Cm-Δ和Cm+Δ之间的频率定义为真实次超同步谐波频率f3；Δ为预设值。在一种优选的方案中，所述的TC≤TA。在一种优选的方案中，所述的次超同步谐波对应频段分别是次同步谐波频段为10~40Hz，超同步谐波频段为60~90Hz。在一种优选的方案中，所述的矢量匹配是一种逼近s域频率响应函数f(s)的技术，其数学表达式为式中，L为矢量匹配拟合阶数；d、e、αl、pl是拟合参数，其中，一次项参数e和常数项参数d为可选项，对应非振荡模式；留数αl和极点pl为必选项，且两者存在两种情况：同时为实数或同时为虚数，当两者同时为虚数时，代表振荡模式。在一种优选的方案中，所述的f0通过以下公式进行计算：式中，imag（pl）表示取pl的虚部。在一种优选的方案中，所述的A0是x(n)在时间窗口TA中的工频分量的幅值的0.1‰。在一种优选的方案中，所述的TB是30s。在一种优选的方案中，所述的聚类算法是DBSCAN密度聚类算法。在一种优选的方案中，所述的Δ=0.5Hz。与现有技术相比，本专利技术技术方案的有益效果是：本专利技术能够实现在复杂时变运行环境下电网中次超同步谐波的准实时监测。在新能源并网后，大量电力电子控制装置并存的背景下优势更为突出。本专利技术对及时发现并切除电网中影响系统稳定性的谐波源具有重要指导意义。附图说明图1是本实施例流程图。图2是本实施例实测电流采样序列。图3是本实施例实测电流采样序列局部放大图。图4是本实施例初步次超同步谐波。图5是本实施例经过聚类识别后真实次超同步谐波。图6是本实施真实次超同步谐波的频率随时间变化图。图7是本实施真实次超同步谐波的幅值随时间变化图。具体实施方式附图仅用于示例性说明，不能理解为对本专利的限制；对于本领域技术人员来说，附图中某些公知结构及其说明可能省略是可以理解的。下面结合附图和实施例对本专利技术的技术方案做进一步的说明。如图1所示，一种基于矢量匹配的电力系统次超同步谐波检测方法，包括如下步骤：S1．从录波装置中获取电流采样序列x(n)，滤除工频分量得到序列y(n)：本实施例从风电汇集站PMU装置监测的线路电流进行采样，使用的带阻滤波器阶数为6，阻带截止频率为45Hz和55Hz，通带截止频率为40Hz和60Hz，通带内衰减<1dB,阻带内衰减>10dB，采样频率为1200Hz，采样时长30s。图2为采样信号图，由于采样点过多，图像过于密集，无法直观表达；图3为采样信号局部放大图。S2．设定时间窗口TA进行FFT得到Y(k)：在时间窗口TA内，将y(n)进行快速傅里叶变换FFT得到序列Y(k)，k为频域离散序列号；TA为预设值；本实施例中，TA=4s。S3．Y(k)进行矢量匹，得到谐波频率f0：将次超同步谐波对应频段内的Y(k)进行矢量匹配，得到极点pl，通过pl中的共轭极点对得到对应的谐波频率f0；本实施例中，次同步谐波频段为10~40Hz，超同步谐波频段为60~90Hz，矢量匹配的数学表达式是式中，L=6，可选项d和e未进行拟合，仅拟合必选项αl和pl；f0通过以下公式进行计算：式中，imag（pl）表示取pl的虚部。S4．f0对应的点是Y(k)的极值点，且f0对应的点的幅值A1≥A0：若f0对应的点是Y(k)的极值点，且f0对应的点的幅值A1≥A0，则f0认定为初步次超同步谐波频率f1，执行S5，A1认定为初步次超同步谐波幅值A2；A0为预设值；若f0不符合要求，则执行S6；本实施例中，A0是x(n)在时间窗口TA中的工频分量的幅值的0.1‰。S5．记录f0及其对应的A1，形成初步次超同步谐波集合：将初步次超同步谐波频率f1和初步次超同步谐波幅值A2进行记录，得到初步次超同步谐波集合，所述的初步次超同步谐波集合中的每个元素包括三个维度：时间t，频率f1和幅值A2；S6．设定时间窗口TB滑动时间窗口TA：在时间窗口TB每相隔步长TC就移动时间窗口TA，所述的TB和TC为预设值；本实施例中，TC=0.25s，TB=30s。图4为初步次超同步谐波集合，从中可以得知，部分谐波是随时间持续存在的，而部分谐波则是随时间偶然存在的。S7．TA没有接触TB的结束时间：若TA没有接触TB的结束时间，则执行S2；若TA已经接触TB的结束时间，则执行S8；S8．对初步次超同步谐波集合进行聚类，得到聚类中心：对时间窗口TB内的初步次超同步谐波频率集合通过聚类算法得到聚类中心Cm(m=1，2，…，M，M为聚类中心的数目)；S9．通过聚类中心识别真实次超同步谐波频率：频率介于Cm-Δ和Cm+Δ之间的频率定义为真实次超同步谐波频率f3；Δ为预设值。本实施例中，Δ=0.5Hz；使用了DBSCAN密度聚类算法从中检测出4种主要间谐波，其频率分别波动在75本文档来自技高网...

【技术保护点】
一种基于矢量匹配的电力系统次超同步谐波检测方法，其特征在于，包括如下步骤：S1.从录波装置中获取电流采样序列x(n)，从x(n)中滤除工频分量得到序列y(n)，n为时域采样序列号；S2.在时间窗口TA内，将y(n)进行快速傅里叶变换FFT得到序列Y(k)，k为频域离散序列号；TA为预设值；S3.将次超同步谐波对应频段内的Y(k)进行矢量匹配，得到极点pl，通过pl中的共轭极点对得到对应的谐波频率f0；S4.若f0对应的点是Y(k)的极值点，且f0对应的点的幅值A1≥A0，则f0认定为初步次超同步谐波频率f1，执行S5，A1认定为初步次超同步谐波幅值A2；A0为预设值；若f0不符合要求，则执行S6；S5.将初步次超同步谐波频率f1和初步次超同步谐波幅值A2进行记录，得到初步次超同步谐波集合，所述的初步次超同步谐波集合中的每个元素包括三个维度：时间t，频率f1和幅值A2；S6.在时间窗口TB每相隔步长TC就移动时间窗口TA，若TA没有接触TB的结束时间，则执行S2；若TA已经接触TB的结束时间，则执行S7；所述的TB和TC为预设值；S7.对时间窗口TB内的初步次超同步谐波频率集合通过聚类算法得到聚类中心Cm(m＝1,2,…,M，M为聚类中心的数目)，频率介于Cm‑Δ和Cm+Δ之间的频率定义为真实次超同步谐波频率f3；Δ为预设值。...

【技术特征摘要】
1.一种基于矢量匹配的电力系统次超同步谐波检测方法，其特征在于，包括如下步骤：S1.从录波装置中获取电流采样序列x(n)，从x(n)中滤除工频分量得到序列y(n)，n为时域采样序列号；S2.在时间窗口TA内，将y(n)进行快速傅里叶变换FFT得到序列Y(k)，k为频域离散序列号；TA为预设值；S3.将次超同步谐波对应频段内的Y(k)进行矢量匹配，得到极点pl，通过pl中的共轭极点对得到对应的谐波频率f0；S4.若f0对应的点是Y(k)的极值点，且f0对应的点的幅值A1≥A0，则f0认定为初步次超同步谐波频率f1，执行S5，A1认定为初步次超同步谐波幅值A2；A0为预设值；若f0不符合要求，则执行S6；S5.将初步次超同步谐波频率f1和初步次超同步谐波幅值A2进行记录，得到初步次超同步谐波集合，所述的初步次超同步谐波集合中的每个元素包括三个维度：时间t，频率f1和幅值A2；S6.在时间窗口TB每相隔步长TC就移动时间窗口TA，若TA没有接触TB的结束时间，则执行S2；若TA已经接触TB的结束时间，则执行S7；所述的TB和TC为预设值；S7.对时间窗口TB内的初步次超同步谐波频率集合通过聚类算法得到聚类中心Cm(m＝1,2,…,M，M为聚类中心的数目)，频率介于Cm-Δ和Cm+Δ之间的频率定义为真实次超同步谐波频率f3；Δ为预设值。2.根据权利要求1所述的电力系统次超同步谐波检测方法，其特征在于，所述的TC≤TA。3.根据权利要求1或2所述的电力系统次超同步谐波检测方法，其特征在于，所述的...

【专利技术属性】
技术研发人员：焦邵麟，曾耿晖，李一泉，刘玮，屠卿瑞，张智锐，
申请(专利权)人：广东电网有限责任公司电力调度控制中心，
类型：发明
国别省市：广东,44